Premium
High‐order curvilinear mesh generation technique based on an improved radius basic function approach
Author(s) -
Zhao Zhong,
Li Ming,
He Lei,
Shao Shuai,
Zhang Laiping
Publication year - 2019
Publication title -
international journal for numerical methods in fluids
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.938
H-Index - 112
eISSN - 1097-0363
pISSN - 0271-2091
DOI - 10.1002/fld.4741
Subject(s) - mesh generation , solver , curvilinear coordinates , interpolation (computer graphics) , boundary (topology) , boundary layer , triangle mesh , grid , signed distance function , computational fluid dynamics , algorithm , geometry , computer science , topology (electrical circuits) , mathematics , finite element method , mathematical optimization , mathematical analysis , polygon mesh , engineering , mechanics , computer graphics (images) , structural engineering , physics , animation , combinatorics
Summary A high‐order curvilinear hybrid mesh generation technique is developed for high‐order numerical method (eg, discontinuous Galerkin method) applications to improve the accuracy for problems with curve boundary. The grid generation technique is based on an improved radius basic function (RBF) approach by which the straight‐edge mesh is converted into high‐order curve mesh. Firstly, an initial straight‐edge mesh is prepared by traditional grid generation software. Then, high‐order interpolation points are inserted into the mesh entities such as edges, faces, and cells according to the final demand of mesh order. To preserve the original geometry, the inserted points on solid wall are then projected onto the CAD model using an open source tool “Open Cascade.” Finally, other inserted points in the field near the solid wall are moved to appropriate positions by the improved RBF approach to avoid tangled cells. If we use the original RBF approach, then the inserted points on the edge and face entities normal to the solid boundary in the region of boundary layer will move to improper positions. To overcome this problem, a weighting based on the local grid aspect ratio between normal direction and tangential direction is introduced into the baseline RBF approach. Three typical configurations are tested to validate the mesh generator. Meanwhile, a third‐order solution of subsonic flow over an analytical 3D body of revolution in the second International Workshop on High‐Order CFD Methods is supplied by a discontinuous Galerkin solver. These numerical tests demonstrate the potential capability of present technique for high‐order simulations of complex geometries.