Using the medial axis to represent flow features for feature-aligned unstructured quad-dominant mesh generation

Abstract A method to represent complex flow structures as geometric entities, using the medial axis for the purpose of generating flow feature-aligned meshes, is presented. The medial curves representing the flow features are embedded in the domain to guide mesh alignment and location of background sources, which influence element spacing and orientation when the mesh is regenerated. The resulting meshes possess regions of high-quality anisotropic quadrilateral elements locally aligned with the flow features. To further optimise the mesh based on the solution, a nodal movement adaptation scheme is applied which allows full alignment of the cell interfaces with the flow features, making the mesh more suitable for the one-dimensional approximate Riemann solvers and higher-order reconstruction schemes. The technique and its potential for improving the numerical simulation is presented for a range of test cases for the solution of the Euler and Navier–Stokes equations for compressible flows, which include shock waves, boundary layers and steady trailing wakes as the dominant flow features.

[1]  M. Fortin,et al.  Anisotropic mesh adaptation: towards user‐independent, mesh‐independent and solver‐independent CFD. Part II. Structured grids , 2002 .

[2]  Meng-Sing Liou,et al.  A Flux Splitting Scheme with High-Resolution and Robustness for Discontinuities(Proceedings of the 12th NAL Symposium on Aircraft Computational Aerodynamics) , 1994 .

[3]  D. Mavriplis.  UNSTRUCTURED GRID TECHNIQUES , 1997 .

[4]  Ted D. Blacker,et al.  Paving: A new approach to automated quadrilateral mesh generation , 1991 .

[5]  M. Fortin,et al.  Anisotropic mesh adaptation: towards user‐independent, mesh‐independent and solver‐independent CFD. Part I: general principles , 2000 .

[6]  Ko-Foa Tchon,et al.  Automated refinement of conformal quadrilateral and hexahedral meshes , 2004 .

[7]  Pascal Frey,et al.  Anisotropic mesh adaptation for CFD computations , 2005 .

[8]  W. K. Anderson,et al.  An implicit upwind algorithm for computing turbulent flows on unstructured grids , 1994 .

[9]  Frédéric Alauzet,et al.  Multi-Dimensional Continuous Metric for Mesh Adaptation , 2006, IMR.

[10]  Dimitri J. Mavriplis,et al.  Unstructured-Mesh Discretizations and Solvers for Computational Aerodynamics , 2007 .

[11]  M. Carpenter,et al.  Accuracy of Shock Capturing in Two Spatial Dimensions , 1999 .

[12]  Michel Fortin,et al.  A directionally adaptive methodology using an edge-based error estimate on quadrilateral grids , 1996 .

[13]  Ning Qin,et al.  Flow feature aligned grid adaptation , 2006 .

[14]  T. Barth.  A 3-D upwind Euler solver for unstructured meshes , 1991 .

[15]  Aramais R. Zakharian,et al.  Two-dimensional Riemann solver for Euler equations of gas dynamics , 2001 .

[16]  Bharat K. Soni,et al.  Handbook of Grid Generation , 1998 .

[17]  P. George,et al.  Mesh Generation: Application to Finite Elements , 2007 .

[18]  V. Venkatakrishnan.  On the accuracy of limiters and convergence to steady state solutions , 1993 .

[19]  D. Holmes,et al.  Solution of the 2D Navier-Stokes equations on unstructured adaptive grids , 1989 .

[20]  J. Peiro,et al.  Adaptive remeshing for three-dimensional compressible flow computations , 1992 .

[21]  M. Leatham,et al.  Automatic mesh generation for rapid-response Navier-Stokes calculations , 2000 .

[22]  Nigel P. Weatherill,et al.  A procedure for efficient generation of solution adapted unstructured grids , 1995 .

[23]  D. Ait-Ali-Yahia,et al.  Anisotropic mesh adaptation for 3D flows on structured and unstructured grids , 2000 .

[24]  Alan M. Shih,et al.  Solution Adaptive Mesh Generation Using Feature-Aligned Embedded Surface Meshes , 2007 .

[25]  T. Gerhold,et al.  Calculation of Complex Three-Dimensional Configurations Employing the DLR-tau-Code , 1997 .

[26]  D. Scott McRae,et al.  r-Refinement grid adaptation algorithms and issues , 2000 .

[27]  H. Borouchaki,et al.  Adaptive triangular–quadrilateral mesh generation , 1998 .

[28]  Robert Haimes,et al.  Shock detection from computational fluid dynamics results , 1999 .

[29]  Alan M. Shih,et al.  Solution Adaptive Mesh Generation Using Feature-Aligned Embedded Surface Meshes , 2009 .

[30]  Yasushi Ito,et al.  A Solution-Based Adaptive Redistribution Method for Unstructured Meshes , 2006, IMR.

[31]  M. F. Peeters,et al.  A 3-D adaptive anisotropic method for external and internal flows , 1998 .

[32]  Herbert Edelsbrunner,et al.  Three-dimensional alpha shapes , 1994, ACM Trans. Graph..