Robust untangling of curvilinear meshes

JOURNAL OF COMPUTATIONAL PHYSICS

卷:254

Article

Toulorge, Thomas¹  Geuzaine, Christophe²  Remacle, Jean-Francois¹  Lambrechts, Jonathan^1,3 

[1] Catholic Univ Louvain, Inst Mech Mat & Civil Engn iMMC, B-1348 Louvain, Belgium

[2] Univ Liege, Dept Elect Engn & Comp Sci, Montefiore Inst B28, B-4000 Liege, Belgium

[3] Fonds Natl Rech Sci, B-1050 Brussels, Belgium

关键词: High order methods;    Finite elements;    Mesh generation;   

DOI  :  10.1016/j.jcp.2013.07.022

来源: Elsevier

【 摘 要 】

This paper presents a technique that allows to untangle high-order/curvilinear meshes. The technique makes use of unconstrained optimization where element Jacobians are constrained to lie in a prescribed range through moving log-barriers. The untangling procedure starts from a possibly invalid curvilinear mesh and moves mesh vertices with the objective of producing elements that all have bounded Jacobians. Bounds on Jacobians are computed using the results of Johnen et al. (2012, 2013) [1,2]. The technique is applicable to any kind of polynomial element, for surface, volume, hybrid or boundary layer meshes. A series of examples demonstrate both the robustness and the efficiency of the technique. The final example, involving a time explicit computation, shows that it is possible to control the stable time step of the computation for curvilinear meshes through an alternative element deformation measure. (C) 2013 Elsevier Inc. All rights reserved.

【 授权许可】

Free