High-order 2D mesh curving methods with a piecewise linear target and application to Helmholtz problems

Abstract High-order simulation techniques typically require high-quality curvilinear meshes. In most cases, mesh curving methods assume that the exact geometry is known. However, in some situations only a fine linear FEM mesh is available and the connection to the CAD geometry is lost. In other applications, the geometry may be represented as a set of scanned points. In this paper, two curving methods are described that take a piecewise fine linear mesh as input: a least squares approach and a continuous optimization in the H 1 -seminorm. Hierarchic, modal shape functions are used as basis for the geometric approximation. This approach allows to create very high-order curvilinear meshes efficiently ( q > 4 ) without having to optimize the location of non-vertex nodes. The methods are compared on two test geometries and then used to solve a Helmholtz problem at various input frequencies. Finally, the main steps for the extension to 3D are outlined.

[1]  Xevi Roca,et al.  Defining an 2-disparity Measure to Check and Improve the Geometric Accuracy of Non-interpolating Curved High-order Meshes , 2015 .

[2]  Amaury Johnen,et al.  The Generation of Valid Curvilinear Meshes , 2015 .

[3]  H. Bériot,et al.  Plane wave basis in Galerkin BEM for bidimensional wave scattering , 2010 .

[4]  Mark S. Shephard,et al.  Tracking Adaptive Moving Mesh Refinements in 3D Curved Domains for Large-Scale Higher Order Finite Element Simulations , 2008, IMR.

[5]  Alfredo Bermúdez,et al.  An optimal perfectly matched layer with unbounded absorbing function for time-harmonic acoustic scattering problems , 2007, J. Comput. Phys..

[6]  Krzysztof J. Fidkowski,et al.  Improving High-Order Finite Element Approximation Through Geometrical Warping , 2015 .

[7]  M. Shephard,et al.  Geometry representation issues associated with p-version finite element computations , 1997 .

[8]  M. Floater.  Mean value coordinates , 2003, Computer Aided Geometric Design.

[9]  I. Doležel,et al.  Higher-Order Finite Element Methods , 2003 .

[10]  Helmut Alt,et al.  Approximate matching of polygonal shapes , 1995, SCG '91.

[11]  J. Remacle,et al.  The Influence of Geometric Approximation on the Accuracy of High Order Methods , 2002 .

[12]  K. Morgan,et al.  The generation of arbitrary order curved meshes for 3D finite element analysis , 2013 .

[13]  Alan M. McIvor,et al.  A comparison of local surface geometry estimation methods , 1997, Machine Vision and Applications.

[14]  Jörg Stiller,et al.  Generation of High-Order Polynomial Patches from Scattered Data , 2014 .

[15]  Xiangmin Jiao,et al.  Identification of C1 and C2 discontinuities for surface meshes in CAD , 2008, Comput. Aided Des..

[16]  David Moxey,et al.  23rd International Meshing Roundtable (IMR23) A thermo-elastic analogy for high-order curvilinear meshing with control of mesh validity and quality , 2014 .

[17]  Hadrien Bériot,et al.  A comparison of high-order polynomial and wave-based methods for Helmholtz problems , 2016, J. Comput. Phys..

[18]  Nonisoparametric formulations for the three-dimensional boundary element method , 1988 .

[19]  S. Rebay,et al.  High-Order Accurate Discontinuous Finite Element Solution of the 2D Euler Equations , 1997 .

[20]  Onur Atak,et al.  HIGH-ORDER CURVED MESH GENERATION BY USING A FINE LINEAR TARGET MESH , 2016 .

[21]  Xevi Roca,et al.  Defining Quality Measures for Validation and Generation of High-Order Tetrahedral Meshes , 2013, IMR.

[22]  R. Cárdenas.  NURBS-Enhanced Finite Element Method (NEFEM) , 2009 .

[23]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[24]  Xevi Roca,et al.  Optimization of a regularized distortion measure to generate curved high‐order unstructured tetrahedral meshes , 2015 .

[25]  Onur Atak,et al.  Comparison of 2D boundary curving methods with modal shape functions and a piecewise linear target mesh , 2017 .

[26]  Hadrien Beriot,et al.  Efficient implementation of high‐order finite elements for Helmholtz problems , 2016 .

[27]  Tomás Vejchodský,et al.  Imposing orthogonality to hierarchic higher-order finite elements , 2007, Math. Comput. Simul..

[28]  Christophe Geuzaine,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

[29]  J. Hesthaven,et al.  Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications , 2007 .

[30]  Christophe Geuzaine,et al.  Robust untangling of curvilinear meshes , 2013, J. Comput. Phys..

[32]  Ludwig Fahrmeir,et al.  Regression: Models, Methods and Applications , 2013 .

[33]  Wim Desmet,et al.  Curved Boundary Treatments for the Discontinuous Galerkin Method Applied to Aeroacoustic Propagation , 2009 .