Anisotropic Mesh Generation

From crtc.cs.odu.edu
Revision as of 14:56, 13 April 2018 by Ctsolakis (talk | contribs) (Solution Based Adaptation)
Jump to: navigation, search

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 [[File:RefinedSurface_onera03.png|500px]
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 )