Space-time CutFEM on overlapping meshes II: simple discontinuous mesh evolution

Mats G. Larson,Carl Lundholm

DOI: https://doi.org/10.1007/s00211-024-01413-y

IF: 2.5

2024-05-29

Numerische Mathematik

Abstract:We present a cut finite element method for the heat equation on two overlapping meshes: a stationary background mesh and an overlapping mesh that evolves inside/"on top" of it. Here the overlapping mesh is prescribed by a simple discontinuous evolution, meaning that its location, size, and shape as functions of time are discontinuous and piecewise constant . For the discrete function space, we use continuous Galerkin in space and discontinuous Galerkin in time, with the addition of a discontinuity on the boundary between the two meshes. The finite element formulation is based on Nitsche's method. The simple discontinuous mesh evolution results in a space-time discretization with a slabwise product structure between space and time which allows for existing analysis methodologies to be applied with only minor modifications. We follow the analysis methodology presented by Eriksson and Johnson (SIAM J Numer Anal 28(1):43–77, 1991; SIAM J Numer Anal 32(3):706–740, 1995). The greatest modification is the introduction of a Ritz-like "shift operator" that is used to obtain the discrete strong stability needed for the error analysis. The shift operator generalizes the original analysis to some methods for which the discrete subspace at one time does not lie in the space of the stiffness form at the subsequent time. The error analysis consists of an a priori error estimate that is of optimal order with respect to both time step and mesh size. We also present numerical results for a problem in one spatial dimension that verify the analytic error convergence orders.

mathematics, applied

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to handle the simple discontinuous motion between a fixed background mesh and an overlapping mesh moving on it when using the Cut Finite Element Method (CutFEM) to solve the heat equation on two overlapping meshes. Specifically, the paper focuses on the case when the position, size and shape of the overlapping mesh are discontinuous and piecewise - constant as a function of time. This method results in a space - time discretization with a piecewise product structure between space and time, which allows the existing analysis methods to still be applicable after some minor modifications. The main contribution of the paper lies in introducing a "shift operator" similar to the Ritz projection to obtain discrete strong stability, which is crucial for error analysis. Through this method, the author is able to provide a priori error estimates of the optimal order of time and mesh size and verify these theoretical results through numerical results.