hp-version time domain boundary elements for the wave equation on quasi-uniform meshes

Heiko Gimperlein,Ceyhun Oezdemir,David Stark,Ernst P. Stephan

DOI: https://doi.org/10.1016/j.cma.2019.07.018

2019-07-08

Abstract:Solutions to the wave equation in the exterior of a polyhedral domain or a screen in $\mathbb{R}^3$ exhibit singular behavior from the edges and corners. We present quasi-optimal $hp$-explicit estimates for the approximation of the Dirichlet and Neumann traces of these solutions for uniform time steps and (globally) quasi-uniform meshes on the boundary. The results are applied to an $hp$-version of the time domain boundary element method. Numerical examples confirm the theoretical results for the Dirichlet problem both for screens and polyhedral domains.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: in three - dimensional space, the singular behavior on the boundary of the solution of the wave equation for polyhedral domains or outside the screen. Specifically, the paper focuses on how to use the hp - version of the time - domain boundary element method (BEM) on quasi - uniform meshes to approximate the Dirichlet and Neumann traces of the wave equation and obtain the optimal convergence rate.

### Specific description of the problem

1. **Boundary value problem of the wave equation**
The paper considers the initial - boundary value problem of the wave equation outside the polyhedral domain or screen \(\Gamma\):
\[
c^{-2} \frac{\partial^2 u}{\partial t^2}(t,x) - \Delta_x u(t,x) = 0 \quad \text{in } \mathbb{R}_t^+ \times \Omega_x,
\]
where \(\Omega\subset\mathbb{R}^3\) is the exterior of the polyhedral domain or screen \(\Gamma\), and the boundary conditions are the given Dirichlet \(u|_\Gamma = g\) or Neumann \(\frac{\partial u}{\partial n}|_\Gamma = g\) conditions.

2. **Singular behavior**
The solution \(u\) exhibits singular behavior at non - smooth boundary points (such as edges and corners). These singular behaviors pose challenges to numerical solutions, especially in the time domain.

3. **hp - version of the time - domain boundary element method**
The paper introduces a space - time hp - version of the time - domain boundary element method for solving the wave equation with inhomogeneous Dirichlet or Neumann boundary conditions. This method is based on tensor - product - form space and time piecewise polynomial functions and is discretized on quasi - uniform meshes.

### Main objectives

- **Approximation properties**
Study the approximation properties of the Dirichlet and Neumann traces of the solution of the wave equation on quasi - uniform meshes, especially improving the approximation accuracy by increasing the polynomial order \(p\) and reducing the mesh size \(h\).

- **Convergence rate**
Obtain quasi - optimal convergence rates, especially the convergence analysis in space - time anisotropic Sobolev spaces.

### Methods and results

- **Asymptotic expansion**
Using the work of Plamenevskii et al., give the asymptotic expansion of the solution of the wave equation, including explicit singular functions and smoother terms.

- **Numerical experiments**
Verify the theoretical results through numerical experiments, especially for the wave equations outside the circular screen and the icosahedron, and confirm the predicted convergence rates.

### Conclusion
The paper successfully demonstrates how to use the hp - version of the time - domain boundary element method on quasi - uniform meshes to efficiently approximate the solution of the wave equation and obtain quasi - optimal convergence rates. This method shows significant advantages in dealing with the wave equation outside polyhedral domains or screens.