A Space-Time Cut Finite Element Method with quadrature in time

We consider convection-diffusion problems in time-dependent domains and present a space-time finite element method based on quadrature in time which is simple to implement and avoids remeshing procedures as the domain is moving. The evolving domain is embedded in a domain with fixed mesh and a cut finite element method with continuous elements in space and discontinuous elements in time is proposed. The method allows the evolving geometry to cut through the fixed background mesh arbitrarily and thus avoids remeshing procedures. However, the arbitrary cuts may lead to ill-conditioned algebraic systems. A stabilization term is added to the weak form which guarantees well-conditioned linear systems independently of the position of the geometry relative to the fixed mesh and in addition makes it possible to use quadrature rules in time to approximate the space-time integrals. We review here the space-time cut finite element method presented in Hansbo et al. (Comput. Methods Appl. Mech. Eng. 307: 96–116, 2016) where linear elements are used in both space and time and extend the method to higher order elements for problems on evolving surfaces (or interfaces). We present a new stabilization term which also when higher order elements are used controls the condition number of the linear systems from cut finite element methods on evolving surfaces. The new stabilization combines the consistent ghost penalty stabilization (Burman, C. R. Acad. Sci. Paris, Ser. I 348(21–22):1217–1220, 2010) with a term controlling normal derivatives at the interface.

[1]  Arnold Reusken,et al.  Analysis of trace finite element methods for surface partial differential equations , 2015 .

[2]  Christoph Lehrenfeld,et al.  The Nitsche XFEM-DG Space-Time Method and its Implementation in Three Space Dimensions , 2014, SIAM J. Sci. Comput..

[3]  Peter Hansbo,et al.  Cut finite element methods for partial differential equations on embedded manifolds of arbitrary codimensions , 2016, ESAIM: Mathematical Modelling and Numerical Analysis.

[4]  Shinichi Kawahara Galerkin-type approximations which are discontinuous in time for parabolic equations in a variable domain , 1977 .

[5]  Christoph Lehrenfeld,et al.  High order unfitted finite element methods on level set domains using isoparametric mappings , 2015, ArXiv.

[6]  Mark Sussman,et al.  An Efficient, Interface-Preserving Level Set Redistancing Algorithm and Its Application to Interfacial Incompressible Fluid Flow , 1999, SIAM J. Sci. Comput..

[7]  Peter Hansbo,et al.  Cut finite element methods for coupled bulk–surface problems , 2014, Numerische Mathematik.

[8]  Kuan-Yu Chen,et al.  A conservative scheme for solving coupled surface-bulk convection-diffusion equations with an application to interfacial flows with soluble surfactant , 2014, J. Comput. Phys..

[9]  R. I. Saye,et al.  High-Order Quadrature Methods for Implicitly Defined Surfaces and Volumes in Hyperrectangles , 2015, SIAM J. Sci. Comput..

[10]  Peter Hansbo,et al.  Characteristic cut finite element methods for convection–diffusion problems on time dependent surfaces , 2015 .

[11]  Peter Hansbo,et al.  Stabilized CutFEM for the convection problem on surfaces , 2015, Numerische Mathematik.

[12]  Peter Hansbo,et al.  A stabilized cut finite element method for partial differential equations on surfaces: The Laplace–Beltrami operator , 2013, 1312.1097.

[13]  Mats G. Larson,et al.  Stabilization of high order cut finite element methods on surfaces , 2017, IMA Journal of Numerical Analysis.

[14]  Peter Hansbo,et al.  A cut finite element method for coupled bulk-surface problems on time-dependent domains , 2015, 1502.07142.

[15]  V. Thomée Galerkin Finite Element Methods for Parabolic Problems (Springer Series in Computational Mathematics) , 2010 .

[16]  L Liggieri,et al.  Adsorption and partitioning of surfactants in liquid-liquid systems. , 2000, Advances in colloid and interface science.

[17]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[18]  Maxim A. Olshanskii,et al.  A Finite Element Method for Elliptic Equations on Surfaces , 2009, SIAM J. Numer. Anal..

[19]  Sashikumaar Ganesan,et al.  Arbitrary Lagrangian-Eulerian finite-element method for computation of two-phase flows with soluble surfactants , 2012, J. Comput. Phys..

[20]  Charles M. Elliott,et al.  Unfitted Finite Element Methods Using Bulk Meshes for Surface Partial Differential Equations , 2013, SIAM J. Numer. Anal..

[21]  C. Peskin Numerical analysis of blood flow in the heart , 1977 .

[22]  Christoph Lehrenfeld,et al.  Analysis of a High-Order Trace Finite Element Method for PDEs on Level Set Surfaces , 2016, SIAM J. Numer. Anal..

[23]  Maxim A. Olshanskii,et al.  An Eulerian Space-Time Finite Element Method for Diffusion Problems on Evolving Surfaces , 2013, SIAM J. Numer. Anal..

[24]  P. Hansbo,et al.  Fictitious domain finite element methods using cut elements , 2012 .

[25]  T. Hou,et al.  Removing the stiffness from interfacial flows with surface tension , 1994 .

[26]  Maxim A. Olshanskii,et al.  Error Analysis of a Space-Time Finite Element Method for Solving PDEs on Evolving Surfaces , 2014, SIAM J. Numer. Anal..

[27]  T. Fries Towards higher‐order XFEM for interfacial flows , 2015 .

[28]  P. Hansbo,et al.  A cut finite element method for a Stokes interface problem , 2012, 1205.5684.

[29]  Jörg Grande,et al.  Eulerian Finite Element Methods for Parabolic Equations on Moving Surfaces , 2014, SIAM J. Sci. Comput..