A Novel Immersed Boundary Approach for Irregular Topography with Acoustic Wave Equations

Edward Caunt,Rhodri Nelson,Fabio Luporini,Gerard Gorman

2023-09-07

Abstract:Irregular terrain has a pronounced effect on the propagation of seismic and acoustic wavefields but is not straightforwardly reconciled with structured finite-difference (FD) methods used to model such phenomena. Methods currently detailed in the literature are generally limited in scope application-wise or non-trivial to apply to real-world geometries. With this in mind, a general immersed boundary treatment capable of imposing a range of boundary conditions in a relatively equation-agnostic manner has been developed, alongside a framework implementing this approach, intending to complement emerging code-generation paradigms. The approach is distinguished by the use of N-dimensional Taylor-series extrapolants constrained by boundary conditions imposed at some suitably-distributed set of surface points. The extrapolation process is encapsulated in modified derivative stencils applied in the vicinity of the boundary, utilizing hyperspherical support regions. This method ensures boundary representation is consistent with the FD discretization: both must be considered in tandem. Furthermore, high-dimensional and vector boundary conditions can be applied without approximation prior to discretization. A consistent methodology can thus be applied across free and rigid surfaces with the first and second-order acoustic wave equation formulations. Application to both equations is demonstrated, and numerical examples based on analytic and real-world topography implementing free and rigid surfaces in 2D and 3D are presented.

Numerical Analysis,Mathematical Software

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to use the acoustic wave equation to more accurately simulate the propagation of seismic wave fields under irregular terrain conditions. Specifically, the paper focuses on how to effectively represent irregular terrain boundary conditions in the structured finite - difference method, because these conditions have a significant impact on the propagation of seismic wave fields. Traditional treatment methods are either of limited application scope or difficult to be applied to actual complex geometries. For this reason, the author has developed a general immersed boundary method, which can impose various boundary conditions in a relatively equation - independent manner and is compatible with existing code - generation paradigms.

### Main problems
1. **Influence of irregular terrain**: Irregular terrain has a significant impact on the propagation of seismic wave fields, which will lead to complex path effects, making originally clear reflections become overlapping arrival waves. The interaction between waves and terrain will cause phenomena such as diffraction and scattering, focusing and defocusing, etc., and these must be considered when simulating wave - field behavior.
2. **Limitations of existing methods**: The methods described in the current literature are usually of limited application scope or difficult to be applied to actual complex geometries. For example, some methods approximate the free surface by adding a low - impedance layer on the surface, but this method may produce large errors under complex geometries.
3. **Challenges of the structured finite - difference method**: The structured finite - difference method faces challenges when dealing with irregular terrain because it is difficult to accurately represent sharp and inhomogeneous material discontinuities on a regular grid. Ideally, it is hoped that complex curved terrains can be accurately represented while retaining the advantages of the structured grid.

### Solutions
1. **Immersed boundary method**: The paper proposes a new immersed boundary method, which imposes boundary conditions by using the high - dimensional Taylor series extrapolation method near the boundary. This method utilizes a set of points distributed near the boundary and modifies the derivative template to achieve the imposition of boundary conditions.
2. **Consistent methodology**: This method can be consistently applied to free surfaces and rigid surfaces, and is suitable for the formulations of first - order and second - order acoustic wave equations. This makes the method highly flexible and accurate when dealing with different types of boundary conditions.
3. **Numerical verification**: The paper verifies the effectiveness of this method through numerical examples, and shows the two - dimensional and three - dimensional numerical experiment results based on analytical and actual terrains.

### Formula examples
- **Taylor series expansion**:
\[
f(x)=\sum_{m = 0}^{M}\frac{(x - x_0)^m}{m!}\frac{\partial^m f(x_0)}{\partial x^m}+O((x - x_0)^{M + 1})
\]
This formula is used to perform Taylor series expansion of the function \(f(x)\) at a certain point \(x_0\).

- **Multidimensional Taylor series**:
\[
f(x)=\sum_{m = 0}^{\infty}\sum_{n = 0}^{\infty}\sum_{l = 0}^{\infty}\frac{(x - x_0)^m(y - y_0)^n(z - z_0)^l}{m!n!l!}\frac{\partial^{m + n + l}f(x_0)}{\partial x^m\partial y^n\partial z^l}
\]
This formula is used to perform Taylor series expansion of the function \(f(x,y,z)\) in three - dimensional space.

- **Linear system**:
\[
A\delta=\begin{pmatrix}
f(x_{-M/2})\\
\vdots\\
f(x_{M/2})
\end{pmatrix}
\]
This linear system is used to solve derivatives through discrete point values.

### Conclusion
The paper proposes a new immersed boundary method, which can effectively handle the simulation problem of the acoustic wave equation under irregular terrain. This method not only has high flexibility and accuracy in theory, but also shows good performance in practical applications. Verified by numerical experiments, this method has significant advantages when dealing with complex terrains.