Project A2 • Numerical methods for wave problems with nontrivial boundary conditions and nonlocal material laws

Principal investigators

  Prof. Dr. Marlis Hochbruck (7/2015 - )
  Prof. Dr. Christian Lubich (7/2015 - )

Project summary

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 several types of nontrivial boundary conditions:

  • linear and nonlinear dynamic boundary conditions for the treatment of wave phenomena in bounded domains
  • generalized impedance boundary conditions in the context of wave scattering
  • transparent boundary conditions that allow one to truncate an unbounded domain to a bounded non-convex domain

Linear, semilinear, and nonlinear wave equations with 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.

The animations show the results of three numerical simulations with different boundary conditions. On the left, we have homogeneous Neumann boundary conditions. In the middle and on the right, the animations show simulations with dynamic boundary conditions. Note that for kinetic boundary conditions (middle) there are small waves that stick at the boundary. Acoustic boundary conditions (right) introduce a second variable at the boundary which is visualized by black arrows.

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.

The picture show two meshs of the blue disc. On the left, we have a triangulation that approximates the disc very crudely. The right picture shows the corresponding isoparametric triangulation with quadratic elements that cover the disc much better.

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 to prove error bounds in the energy norm for non-conforming space discretizations of linear and semi-linear wave equations, which combines a variational approach with semigroup theory. This systematic analysis provides a priori error bounds in terms of interpolation, data and conformity errors of the method. These bounds can be used to derive convergence rates for a large class of equations and their space discretizations in a simple, systematic and modular way ([HHS19, Hip17, Lei17]). Our unified error analysis is also used in projects A3, A4 and A7, e.g., to analyze the error of a heterogeneous multiscale method for Maxwell’s equation [Fre22, HMS19] or for quasilinear wave equations in [HM22, Mai22]. In addition, the unified error analysis serves as a first step 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 [Lei21] and show that the full discretization error is bounded by the sum of the (stiff) time integration error and the space discretization error.
  • $L^2$ error estimates for space discretizations of wave equations require different techniques than bounds in the energy norm. In [HK21], 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 aim at designing time integration schemes that exploit the structure of the corresponding problems. These methods improve the efficiency of standard methods such as, e.g., the implicit Crank-Nicolson method or the explicit leapfrog scheme. For the efficient time integration of problems with locally Lipschitz continuous nonlinearities we proposed a second-order implicit-explicit (IMEX) scheme that treats the linear part of the equation implicitly and the nonlinear part explicitly [HL21]. We proved that the time step size is not restricted by the spatial mesh width but only by the nonlinear part. In order to ensure the stability of the numerical scheme for wave equations with fully nonlinear dynamic boundary conditions, we discretized the nonlinearities in such a way that the quasi-monotonicity of the continuous operators is preserved. Furthermore, we managed to bound the discretization errors of the nonlinearities. While both are straightforward for conforming discretizations, it turns out to be much more involved in the non-conforming case.
  • Energy norm error estimates for wave equations with kinetic or acoustic boundary conditions. Since these problems are stated on smooth domains, we consider isoparametric finite element approximations for the space discretization. Using the unified error analysis, we prove convergence rates in an energy norm for the wave equation with acoustic boundary conditions in [HHS19] and kinetic boundary conditions in [Hip17]. Moreover, we provide a rigorous error analysis of the wave equation with semi-linear kinetic boundary conditions [HL20]. We proved novel optimal full discretization error bounds in the energy norm for wave equations with nonlinear damped kinetic boundary conditions in [Lei21] and with acoustic boundary conditions with nonlinear coupling in [Lei22].

Generalized impedance boundary conditions

Generalized impedance boundary conditions are effective, approximate boundary conditions that describe the scattering of waves in situations where the wave interaction with the material involves a small scale along the boundary of the scatterer.

These three simulations show numerical results in 3D modeling acoustic scattering from a thin-layer boundary condition, strongly absorbing materials, and acoustic boundary conditions. As the incoming wave, a planar acoustic wave has been chosen, which interacts with a halfpipe shape. The plots depict the resulting scattered waves, on the plane which intersects the halfpipe through the middle. Details on the simulation are found in the final section of [BLN22].

The main motivation for this subject arises from materials with a thin coating (with the thickness of the coating $\varepsilon$ as the small scale), which are modeled with the effective boundary condition \begin{align*} \partial_{\nu} u &=\varepsilon(\partial_t^2u -\Delta_{\Gamma} u) \quad\text{ on }\Gamma. \end{align*}

Another effective model of interest has been derived in the study of strongly absorbing materials, originally formulated by Leontovich in the 1930s to study the interaction of the earth's surface with radio waves, reads in the time-domain \begin{align*} \partial_{\nu} u &=\dfrac{1}{\varepsilon}\partial_t^{1/2} u \quad\text{ on }\Gamma. \end{align*} Further, we studied acoustic boundary conditions. All of these boundary conditions were considered in the context of wave scattering, namely on an exterior (unbounded) domain, which is the complement of a bounded scatterer. An incoming wave, initially away from the boundary, illuminates a medium at rest and interacts with the material.

To give a complete numerical analysis of problems of this type, the following steps are taken:

  • stable time-dependent boundary integral equations (BIEs) are derived and analyzed
  • the convolution quadrature method is employed to discretize the BIEs in time
  • the boundary element method is employed to discretize the BIEs in space
  • an error analysis is conducted, either based on Laplace domain techniques or on energy techniques
Research and results
  • In [BLN22], the boundary conditions described above have been treated in the context of acoustic scattering. The boundary integral equations are discretized by polynomial boundary elements and the convolution quadrature method based on multistep methods. Error bounds are proved and illustrated by numerical experiments.
  • The electromagnetic counterparts of the described boundary conditions appear in many applications. In [NKL22], we transfer and expand the results of the acoustic setting to the electromagnetic setting. Along the way, several new fundamental time-harmonic bounds are shown. Applying the Runge–Kutta based convolution quadrature method combined with a Raviart–Thomas boundary element discretization yields stable and convergent schemes for the corresponding boundary conditions. Numerical experiments demonstrate the use of the proposed method.
  • Finally, in [Nic22] we considered a class of nonlinear boundary conditions. Whereas the boundary integral equation and the corresponding full discretization follows from the previous paper [NKL22], the (nonlinear) error analysis is fundamentally different. Based on energy techniques relying on time-discrete transmission problems error bounds are shown, by leveraging several new innovations in the error analysis. Numerical experiments illustrate the theoretical results and show the effectiveness of the method.

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.

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.

Publications

  1. , , and . Optimal $W^{1,\infty}$-estimates for an isoparametric finite element discretization of elliptic boundary value problems. Electron. Trans. Numer. Anal., 58:1–21, January . URL https://doi.org/10.1553/etna_vol58s1. [preprint] [bibtex]

  2. . A unified error analysis for nonlinear wave-type equations with application to acoustic boundary conditions. Numer. Math., 152:907–936, October . URL https://doi.org/10.1007/s00211-022-01326-8. [preprint] [files] [bibtex]

  3. and . Full discretization error analysis of exponential integrators for semilinear wave equations. Math. Comp., 91(336):1687–1709, July . URL https://doi.org/10.1090/mcom/3736. [preprint] [files] [bibtex]

  4. , , and . Time-dependent electromagnetic scattering from thin layers. Numer. Math., 150(4):1123–1164, April . URL https://doi.org/10.1007/s00211-022-01277-0. [preprint] [files] [bibtex]

  5. , , and . Time-dependent acoustic scattering from generalized impedance boundary conditions via boundary elements and convolution quadrature. IMA J. Numer. Anal., 42(1):1–26, January . URL https://doi.org/10.1093/imanum/draa091. [preprint] [bibtex]

  6. and . An implicit-explicit time discretization scheme for second-order semilinear wave equations with application to dynamic boundary conditions. Numer. Math., 147(4):869–899, April . URL https://doi.org/10.1007/s00211-021-01184-w. [preprint] [files] [bibtex]

  7. and . Finite element error analysis of wave equations with dynamic boundary conditions: $L^2$ estimates. IMA J. Numer. Anal., 1:638–728, January . URL https://doi.org/10.1093/imanum/drz073. [preprint] [bibtex]

  8. , , and . Correction to: Stable and convergent fully discrete interior-exterior coupling of Maxwell's equations. Numer. Math., 147(4):997–1000, 4 . URL https://doi.org/10.1007/s00211-021-01196-6. [bibtex]

  9. and . Finite element discretization of semilinear acoustic wave equations with kinetic boundary conditions. Electron. Trans. Numer. Anal., 53:522–540, August . URL https://doi.org/10.1553/etna_vol53s522. [preprint] [bibtex]

  10. , , and . A convergent evolving finite element algorithm for mean curvature flow of closed surfaces. Numer. Math., 143(4):797–853, December . URL https://doi.org/10.1007/s00211-019-01074-2. [bibtex]

  11. and . Runge–Kutta convolution coercivity and its use for time-dependent boundary integral equations. IMA J. Numer. Anal., 39(3):1134–1157, July . URL https://doi.org/10.1093/imanum/dry033. [preprint] [bibtex]

  12. , , and . Unified error analysis for nonconforming space discretizations of wave-type equations. IMA J. Numer. Anal., 39(3):1206–1245, July . URL https://doi.org/10.1093/imanum/dry036. [preprint] [bibtex]

  13. , , , and . Analytical and numerical analysis of linear and nonlinear properties of an rf-SQUID based metasurface. Phys. Rev. B, 99(7):075401, February . URL https://doi.org/10.1103/PhysRevB.99.075401. [preprint] [bibtex]

  14. and . Stability and convergence of time discretizations of quasi-linear evolution equations of Kato type. Numer. Math., 138(2):365–388, February . URL https://doi.org/10.1007/s00211-017-0909-3. [preprint] [bibtex]

  15. and . Stable and convergent fully discrete interior-exterior coupling of Maxwell's equations. Numer. Math., 137(1):91–117, September . URL https://doi.org/10.1007/s00211-017-0868-8. [preprint] [bibtex]

  16. and . Numerical analysis of parabolic problems with dynamic boundary conditions. IMA J. Numer. Anal., 37(1):1–39, January . URL https://doi.org/10.1093/imanum/drw015. [preprint] [bibtex]

Preprints

  1. , , and . Time-dependent electromagnetic scattering from dispersive materials. CRC 1173 Preprint 2023/22, Karlsruhe Institute of Technology, November . [files] [bibtex]

  2. . Numerical analysis for electromagnetic scattering from nonlinear boundary conditions. CRC 1173 Preprint 2022/37, Karlsruhe Institute of Technology, August . [files] [bibtex]

Theses

  1. . Wave scattering from nontrivial boundary conditions. PhD thesis, Karlsruhe Institute of Technology (KIT), February . [bibtex]

  2. . A unified error analysis for the numerical solution of nonlinear wave-type equations with application to kinetic boundary conditions. PhD thesis, Karlsruhe Institute of Technology (KIT), February . [bibtex]

  3. . Semilineare Wellengleichungen mit dynamischen Randbedingungen. Master's thesis, Karlsruhe Institute of Technology (KIT), September . [bibtex]

  4. . A unified error analysis for spatial discretizations of wave-type equations with applications to dynamic boundary conditions. PhD thesis, Karlsruhe Institute of Technology (KIT), June . [bibtex]

Other references

  1. and . Error analysis for space discretizations of quasilinear wave-type equations. IMA J. Numer. Anal., 42(3):1963–1990, July . URL https://doi.org/10.1093/imanum/drab073. [preprint] [files] [bibtex]

  2. . The heterogeneous multiscale method for dispersive Maxwell systems. Multiscale Model. Simul., 20(2):769–797, June . URL https://doi.org/10.1137/21M1443960. [preprint] [bibtex]

  3. . Error analysis for full discretizations of quasilinear wave-type equations with two variants of the implicit midpoint rule. IMA J. Numer. Anal., drac010, May . URL https://doi.org/10.1093/imanum/drac010. [preprint] [bibtex]

  4. . An implementation and numerical experiments of the FEM‐BEM coupling for the elastodynamic wave equation in 3D. ZAMM Z. Angew. Math. Mech., 99(12):e201900050, December . URL https://doi.org/10.1002/zamm.201900050. [bibtex]

  5. , , and . Heterogeneous multiscale method for Maxwell's equations. Multiscale Model. Simul., 17(4):1147–1171, October . URL https://doi.org/10.1137/18M1234072. [preprint] [bibtex]

  6. and . Upwind discontinuous Galerkin space discretization and locally implicit time integration for linear Maxwell's equations. Math. Comp., 88(317):1121–1153, May . URL https://doi.org/10.1090/mcom/3365. [preprint] [bibtex]

  7. . The elastic wave equation and the stable numerical coupling of its interior and exterior problems. ZAMM Z. Angew. Math. Mech., 98(7):1261–1283, July . URL https://doi.org/10.1002/zamm.201600236. [bibtex]

Former staff
Name Title Function
Dr. Postdoctoral researcher
M.Sc. Doctoral researcher
PD Dr. Member
Dr. Postdoctoral and doctoral researcher
Dr. Postdoctoral researcher
Dr. Postdoctoral and doctoral researcher
Dr. Postdoctoral and doctoral researcher