An unfitted finite element method using level set functions for extrapolation into deformable diffuse interfaces

• Published in: Journal of computational physics 2022-07, Vol.461, p.111218, Article 111218
• Main Authors: Kuzmin, Dmitri, Bäcker, Jan-Phillip
• Format: Article
• Language: English
• Subjects: Boundary conditions; Computational physics; Conservation laws; Delta function; Diffuse interface; Domains; Extension velocity; Extrapolation; Finite element analysis; Finite element method; Formability; Ghost penalty; Integrals; Interfaces; Interpolation; Level set algorithm; Mathematical analysis; Quadratures; Regularized delta function; Unfitted finite element method
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2022.111218

•Weak imposition of interface conditions in fictitious domain methods.•Efficient closest-point search using an approximate distance function.•Simple reconstruction and extrapolation of data defined on level sets.•Damping functions that restrict computations to a narrow band.•New 2D test problems with fixed and moving embedded boundaries. We explore a new way to handle flux boundary conditions imposed on level sets. The proposed approach is a diffuse interface version of the shifted boundary method (SBM) for continuous Galerkin discretizations of conservation laws in embedded domains. We impose the interface conditions weakly and approximate surface integrals by volume integrals. The discretized weak form of the governing equation has the structure of an immersed boundary finite element method. That is, integration is performed over a fixed fictitious domain. Source terms are included to account for interface conditions and extend the boundary data into the complement of the embedded domain. The calculation of these extra terms requires (i) construction of an approximate delta function and (ii) extrapolation of embedded boundary data into quadrature points. We accomplish these tasks using a level set function, which is given analytically or evolved numerically. A globally defined averaged gradient of this approximate signed distance function is used to construct a simple map to the closest point on the interface. The normal and tangential derivatives of the numerical solution at that point are calculated using the interface conditions and/or interpolation on uniform stencils. Similarly to SBM, extrapolation of data back to the quadrature points is performed using Taylor expansions. Computations that require extrapolation are restricted to a narrow band around the interface. Numerical results are presented for elliptic, parabolic, and hyperbolic test problems, which are specifically designed to assess the error caused by the numerical treatment of interface conditions on fixed and moving boundaries in 2D.