|Prof. Dr. Marlis Hochbruck||(7/2015 - )|
|Prof. Dr. Christian Lubich||(7/2015 - )|
Boundary conditions are an essential part of mathematical models of wave phenomena that are confined in a bounded spatial domain and describe the behavior of the waves at the boundary.
During the rise of computing power and demand for simulation, nontrivial boundary conditions became increasingly important, for example to effectively model complex surface phenomena or as artificial boundary conditions imposed on domains truncated for computational purposes. Often, such nontrivial boundary conditions are not covered by standard numerical analysis.
The aim of this project is the construction of stable numerical schemes and to provide a rigorous error analysis for wave equations with two types of nontrivial boundary conditions: First, transparent boundary conditions that allow to truncate an unbounded domain to a bounded non-convex domain suitable for numerical approximation. Second, dynamic boundary conditions which are used in sophisticated models for surface phenomena.
Transparent boundary conditions
Transparent boundary conditions are often used in the numerical approximation of models for wave-type phenomena on unbounded domains, where the initial data and the source is contained in a bounded domain. Many efficient techniques (perfectly matched layers, fast multipole methods, etc.) are known for convex domains. This subproject is concerned with the case of non-convex domains. A particular difficulty is the fact that waves can leave and reenter the domain.
We devised stable and convergent numerical methods which only involve the interior of this domain and its boundary. The transparent boundary conditions are imposed as a time-dependent boundary integral equation, in which the appearing Calderon operator is non-local in both space and time. Continue reading. Collapse content.
Particular difficulties of the analysis lay in finding the suitable functional analytic setting which permits the right definition of boundary integral operators, proving the coercivity of the Calderon operator, and ensuring stability of the numerical methods.
Research and results
- Maxwell's equations [KL17]: The suitable abstract framework, allowing a boundary integral formulation, is developed. We show the coercivity of the Calderon operator, upon which the stability of the scheme relies. We prove stability and convergence of a numerical method, using a discontinuous Galerkin method (cf. [HS19] from A4) and leapfrog scheme in the interior coupled to boundary elements and (BDF2) convolution quadrature on the boundary.
- Elastodynamics [Ebe18]: Again the key issue is finding the coercivity of the Calderon operator encoding the transparent boundary conditions, and the abstract framework, which allows to prove its coercitvity. Stability and convergence of a numerical scheme, coupling finite elements and leapfrog time-stepping in the interior with boundary elements and (BDF2) convolution quadratures on the boundary, is proved. Numerical experiments and the handling of implementation issues involving the hyper-singular operator are presented in [Ebe19].
- Runge—Kutta coercivity [BL19]: The coercivity property plays a central role in showing stability of the interor—exteror coupling. Convolution quadratures based on A-stable multistep methods are known to preserve this property. In this paper we study the question: “which Runge—Kutta-based convolution quadratures inherit this property?” It is shown that this holds without any restriction for the third-order Radau IIA method, and on permitting a shift in the Laplace domain variable, this holds for all algebraically stable Runge–Kutta methods and hence for methods of arbitrary order. The results are illustrated with various numerical experiments. Continue reading. Collapse content.
Dynamic boundary conditions
In contrast to Dirichlet or Neumann boundary conditions, dynamic boundary conditions do not neglect the momentum of the wave on the boundary. They are particularly suited to model flexible boundaries like an eardrum and often appear as effective models for wave propagation in in thin layers. Hence, the simulation of these kind of boundary conditions is important in many applications.
A particular difficulty for the space discretization are domains with curved boundaries which naturally appear for certain dynamic boundary conditions. Only few numerical methods can treat such domains exactly. Instead, most numerical methods are applied on an approximated domain. This makes the implementation as well as the error analysis more involved. Continue reading. Collapse content.
Despite the fact that the error analysis of such problems mostly use the same ideas, it is usually carried out for a particular combination of wave equation and numerical method. To avoid repetitive work and ease the derivation of new results, we developed the unified error analysis which applies to many different kinds of wave equations and numerical methods.
Research and results
- Unified error analysis of non-conforming space discretizations for wave-type equations: We propose an abstract framework for non-conforming space discretizations of linear and semi-linear wave equations, which combines a variational approach with semigroup theory. The error estimates we derive in this framework allow applicants to easily evaluate the accuracy of a method for their concrete problem by inserting basic information about the numerical method. [DHS19, Hip17, Lei17]. Our unified error analysis is also used in other projects, e.g. to analyze the error of a heterogeneous multiscale method for Maxwell’s equation [Bernhards MA]. Our future aim is to further extend this ansatz to cover nonlinear cases.
- Energy norm error estimates for space discretizations: We consider isoparamteric finite element approximations. Using the unified error analysis, we prove convergence rates in an energy norm for the wave equation with acoustic boundary conditions in [DHS19] and kinetic boundary conditions in [Hip17]. Moreover, we provide a rigorous error analysis of the wave equation with semi-linear kinetic boundary conditions. [Lei17]
- \(L^2\) error estimates for space discretizations of wave equations require more direct techniques then the ones available via the unified error analysis (where we derive estimates in the energy norm). Hence the final convergence results themselves are less 'unified' [HK18]. We show semi-discrete convergence estimates in \(L^2\) norm, using non-conforming finite elements for four representative examples of wave equations with dynamic boundary conditions, working in the unified setting as long as the analysis allows. Fully discrete \(L^2\) norm convergence results, using algebraically stable and coercive Runge-Kutta methods, are also shown.
- Efficient time integration: We are particularly interested in time integration schemes that exploit the structure of the corresponding problems. Such methods can be much more efficient than traditional implicit methods as e.g. the Crank-Nicolson method while having superior stability properties compared to explicit methods as the Leapfrog scheme. We have first results for semi-linear problems which show the second-order convergence of a linear implicit method, which does not require the solution of nonlinear systems of equations [in preparation].
- Full discretization: The framework of the unified error analysis also allows to analyze fully discrete approximations obtained by the method of lines. As examples, we consider the Crank-Nicolson method [Hip17] and algebraically stable and coercive Runge-Kutta methods [BL19] and show that the full discretization error is bounded by the sum of the (stiff) time integration error and the space discretization error.