Abstract

We develop a family of cut finite element methods of different orders based on the discontinuous Galerkin framework, for hyperbolic conservation laws with stationary interfaces in both one and two space dimensions, and for moving interfaces in one space dimension. Interface conditions are imposed weakly and so that both conservation and stability are ensured. A CutFEM with discontinuous elements in space is developed and coupled to standard explicit time stepping schemes for linear advection problems and the acoustic wave problem with stationary interfaces. In the case of moving interfaces, we propose a space-time CutFEM based on discontinuous elements both in space and time for linear advection problems. We show that the proposed CutFEM are conservative and energy stable. For the stationary interface case an a priori error estimate is proven. Numerical computations in both one and two space dimensions support the analysis, and in addition demonstrate that the proposed methods have the expected accuracy.

Highlights

  • A finite element method (FEM) that uses discontinuous piecewise polynomial spaces as trial and test spaces is commonly called a discontinuous Galerkin (DG) method

  • The implementation of the space-time unfitted finite element method we propose is straightforward and simple starting from an implementation of Cut Finite Element Methods (CutFEM) for a stationary interface

  • We have presented two high order CutFEM based on the DG framework, applicable to conservation laws with discontinuous coefficients in the flux across stationary and moving interfaces, respectively

Read more

Summary

Introduction

A finite element method (FEM) that uses discontinuous piecewise polynomial spaces as trial and test spaces is commonly called a discontinuous Galerkin (DG) method. A common technique in connection with Cut Finite Element Methods (CutFEM) is to add ghost penalty stabilization terms in the weak form [3,4]. CutFEM based on discontinuous piecewise polynomial spaces and ghost penalty stabilization has been developed, e.g. see [12] where a time independent linear advection-reaction problem is considered and see [10] for time dependent nonlinear conservation laws. 3, we consider a stationary interface, propose a discontinuous cut finite element discretization in space and perform a stability analysis, and an a priori error estimate is given for the scalar problems. The stability of the semi-discrete scheme is analyzed and we present some examples to show that the method can simulate the moving interface problem with expected accuracy and with conservation.

Model Problem
84 Page 4 of 39
Mesh and Spaces
Weak Formulation
84 Page 8 of 39
The Acoustic System
84 Page 10 of 39
Numerical Results
84 Page 12 of 39
Scalar Problem
84 Page 14 of 39
84 Page 16 of 39
Moving Interface
84 Page 20 of 39
Stability Analysis of the Semi-discrete Scheme
84 Page 22 of 39
Quadrature in Time
Scalar Problem with a Moving Interface
84 Page 24 of 39
A Locally Implicit Method
84 Page 26 of 39
Numerical Examples
Extension to Two Space Dimensions
The Finite Element Method
Numerical Example
Convergence Study
Conservation Study
84 Page 30 of 39
Conclusion
A Positivity of the S-Matrix
B Proof of Theorem 2
84 Page 36 of 39
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call