Difference between revisions of "Anisotropic Mesh Generation"

From crtc.cs.odu.edu
Jump to: navigation, search
(30K Complexity)
 
(16 intermediate revisions by the same user not shown)
Line 1: Line 1:
==AFLR==
+
== Analytic Metric Functions ==
==CDT3D==
+
===Linear metric===
 
 
=== Analytic Metric Functions ===
 
* Linear metric given by  [[File:Linear metric.png| 300px]]
 
  
 +
Metric given by  [[File:Linear metric.png| 300px|link=]]
 +
where h<sub>x</sub> =0.1, h<sub>y</sub>=0.1, h<sub>0</sub> = 0.001, h<sub>z</sub> = h<sub>0</sub> + 2(0.1-h<sub>0</sub>)|z-0.6|
  
 
{|
 
{|
Line 14: Line 13:
 
!colspan="2"| Surface Mesh, cross-cut and detail of the generated volume mesh.
 
!colspan="2"| Surface Mesh, cross-cut and detail of the generated volume mesh.
 
|}
 
|}
 +
 +
{|
 +
|-
 +
|[[File:Cube_qcircle_edges.png|350px]]
 +
|[[File:Cube_qcircle_qual.png|350px]]
 +
|-
 +
!colspan="2"| Quality statistics in the metric space for generated volume mesh.94 % of the edges has length between [0.5,1.5].
 +
For the mean ratio 96%,of the elements have quality higher that 0.7.
 +
|}
 +
 +
=== Quarter of a Cylinder ===
 +
Metric given by [[File:Qcircle_metric.png| 700px]]
 +
 +
where h<sub>x</sub> = min(0.002  5^a, h<sub>max</sub>),
 +
h<sub>y</sub> = min(0.05 &middot; 2<sup>a</sup>, h<sub>max</sub>),
 +
h<sub>z</sub>= h<sub>max</sub>,
 +
h<sub>max</sub> = 0.1,
 +
θ = arctan(x,y),
 +
a = 10 &middot; | 0.75 -
 +
<!-- code for square root -->
 +
&radic;<span style="text-decoration:overline">
 +
(<i>x</i><span style="font-size: 10px;vertical-align:+25%;">2</span>+
 +
<i>y</i><span style="font-size: 10px;vertical-align:+25%;">2</span>)&nbsp;</span>
 +
<!-- code for square root -->
 +
|
 +
{|
 +
|-
 +
|[[File:Cube_qcircle1.png|300px|link=]]
 +
|[[File:Cube_qcircle2.png|300px|link=]]
 +
|[[File:Cube_qcircle_detail.png|300px|link=]]
 +
|-
 +
!colspan="2"| Surface Mesh, cross-cut and detail of the generated volume mesh.
 +
|}
 +
 +
{|
 +
|-
 +
|[[File:Cube_linear_edges.png|350px|link=]]
 +
|[[File:Cube_linear_qual.png|350px|link=]]
 +
|-
 +
!colspan="2"| Quality statistics in the metric space for generated volume mesh.
 +
91 % of the edges has length between [0.5,1.5]. For the mean ratio 79%
 +
of the elements have quality higher that 0.7.
 +
|}
 +
 +
 +
== Solution Based Adaptation ==
 +
 +
=== Onera M6 ===
 +
==== 30K Complexity ====
 +
Surface and background meshes were acquired from [https://github.com/UGAWG UGAWG repository]
 +
{|
 +
|-
 +
|[[File:Onera_initial mesh.png|500px]]
 +
|[[File:RefinedSurface_onera03.png|500px]]
 +
|-
 +
| '''Initial surface mesh acquired from [https://github.com/UGAWG/solution-adapt-cases/blob/master/onera-m6/initial-grid/onera_curvAdapt.meshb]'''
 +
| '''Refined surface mesh using as a 30K complexity background mesh acquired from [https://github.com/UGAWG/solution-adapt-results/blob/master/onera-m6/sa-m084-a306-fun3d/onera03-metric.solb]'''
 +
|-
 +
| [[File:Onera03_side.png|500px]]
 +
| [[File:Onera03_cross_section.png|500px]]
 +
|-
 +
| colspan="2" style="text-align: center;" | ''' cross cuts of the final mesh'''
 +
|}
 +
 +
 +
The final mesh (available [http://cs.odu.edu/crtc/CDT3D/Anisotropic/onera03_CDT3D.7z here] ) was compared against the results from the 2018 Scitech paper ''' Unstructured Grid Adaptation and Solver Technology for Turbulent Flows'''
 +
available at https://arc.aiaa.org/doi/abs/10.2514/6.2018-1103 .
 +
 +
{|
 +
|-
 +
|[[File:Histogram ratio onera03.png|500px]]
 +
|[[File:Histogram qual onera03.png|500px]]
 +
|-
 +
|}
 +
 +
==== 50K Complexity ====
 +
 +
Same input surface but with different metric complexity :
 +
 +
{|
 +
| [[File:Onera05_side.png|500px]]
 +
| [[File:Onera05_cross_section.png|500px]]
 +
|-
 +
| colspan="2" style="text-align: center;" | ''' cross cuts of the final mesh'''
 +
|}
 +
 +
The background mesh is available [https://github.com/UGAWG/solution-adapt-results/blob/master/onera-m6/euler-m084-a306-fun3d/onera05-metric.solb here ]
 +
 +
The final mesh (available [http://cs.odu.edu/crtc/CDT3D/Anisotropic/onera05_CDT3D.7z here] )

Latest revision as of 14:56, 13 April 2018

Analytic Metric Functions

Linear metric

Metric given by Linear metric.png where hx =0.1, hy=0.1, h0 = 0.001, hz = h0 + 2(0.1-h0)|z-0.6|

Cube linear.png Cube linear2.png Cube linear detail.png
Surface Mesh, cross-cut and detail of the generated volume mesh.
Cube qcircle edges.png Cube qcircle qual.png
Quality statistics in the metric space for generated volume mesh.94 % of the edges has length between [0.5,1.5].

For the mean ratio 96%,of the elements have quality higher that 0.7.

Quarter of a Cylinder

Metric given by Qcircle metric.png

where hx = min(0.002 5^a, hmax), hy = min(0.05 · 2a, hmax), hz= hmax, hmax = 0.1, θ = arctan(x,y), a = 10 · | 0.75 - √ (x2+ y2 |

Cube qcircle1.png Cube qcircle2.png Cube qcircle detail.png
Surface Mesh, cross-cut and detail of the generated volume mesh.
Cube linear edges.png Cube linear qual.png
Quality statistics in the metric space for generated volume mesh.

91 % of the edges has length between [0.5,1.5]. For the mean ratio 79% of the elements have quality higher that 0.7.


Solution Based Adaptation

Onera M6

30K Complexity

Surface and background meshes were acquired from UGAWG repository

Onera initial mesh.png RefinedSurface onera03.png
Initial surface mesh acquired from [1] Refined surface mesh using as a 30K complexity background mesh acquired from [2]
Onera03 side.png Onera03 cross section.png
cross cuts of the final mesh


The final mesh (available here ) was compared against the results from the 2018 Scitech paper Unstructured Grid Adaptation and Solver Technology for Turbulent Flows available at https://arc.aiaa.org/doi/abs/10.2514/6.2018-1103 .

Histogram ratio onera03.png Histogram qual onera03.png

50K Complexity

Same input surface but with different metric complexity :

Onera05 side.png Onera05 cross section.png
cross cuts of the final mesh

The background mesh is available here

The final mesh (available here )