A high-order elliptic PDE based level set reinitialisation method using a discontinuous Galerkin discretisation

JOURNAL OF COMPUTATIONAL PHYSICS

卷:379

Article

Adams, Thomas¹  Giani, Stefano¹  Coombs, William M.¹ 

[1] Univ Durham, Dept Engn, South Rd, Durham DH1 3LE, England

关键词: Level set;    Reinitialisation;    Discontinuous Galerkin;   

DOI  :  10.1016/j.jcp.2018.12.003

来源: Elsevier

【 摘 要 】

In this paper, an efficient, high-order accurate, level set reinitialisation method is proposed, based on the elliptic reinitialisation method (Basting and Kuzmin, 2013 [1]), which is discretised spatially using the discontinuous Galerkin (DG) symmetric interior penalty method (SIPG). In order to achieve this a number of improvements have been made to the elliptic reinitialisation method including; reformulation of the underlying minimisation problem driving the solution; adoption of a Lagrange multiplier approach for enforcing a Dirichlet boundary condition on the implicit level set interface; and adoption of a narrow band approach. Numerical examples confirm the high-order accuracy of the resultant method by demonstrating experimental orders of convergence congruent with optimal convergence rates for the SIPG method, that is h(p+1) and h(p) in the L-2 and DG norms respectively. Furthermore, the degree to which the level set function satisfies the Eikonal equation improves proportionally to h(p), and the often ignored homogeneous Dirichlet boundary condition on the interface is shown to be satisfied accurately with a rate of convergence of at least h 2 for all polynomial orders. (C) 2019 The Authors. Published by Elsevier Inc.

【 授权许可】

Free