Generation and propagation of a Tsunami wave : a new mesh adaptation technique

Georges Sadaka

DOI: https://doi.org/10.48550/arXiv.1210.4260

2018-01-22

Abstract:The importance of the study of the propagation of a Tsunami wave came from the complex phenomenon and its natural disasters which represents a major risk for populations. To model this phenomena, we will consider a simplified Boussinesq system of Benjamin, Bona and Mahony type derived by D. Mitsotakis over a flat bottom then over a variable bottom in space and in time and apply this system, first, using a mesh generated using a photo of the Mediterranean sea, second, using a mesh generated using an imported xyz bathymetry for the sea near Java island and then we will consider a realistic example of the Tsunami wave near Java island which happened in <a class="link-external link-http" href="http://2006.We" rel="external noopener nofollow">this http URL</a> choose here to use FreeFem++ software which simplifies the construction of the domain, in particular, one of the advantage of FreeFem++ is that we can build a mesh using a photo and we can easily export bathymetric data in order to consider more realistic simulations where a special adapt mesh technique applied for these two methods is detailed in the sequel.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is related to the simulation of tsunami wave generation and propagation, especially for complex natural phenomena and the natural disaster risks they bring. Specifically, the article focuses on the following points:

1. **Model selection**: Use the simplified Boussinesq system (sBBM), which was proposed by D. Mitsotakis, to describe tsunami waves in cases of flat seabed and seabed varying with space and time.

2. **Mesh generation**: Two different methods are proposed to generate meshes for the Mediterranean Sea area and the sea area near Java Island for simulation:
- Generate a mesh using a photo of the Mediterranean Sea.
- Generate a mesh near Java Island using bathymetric data in XYZ format downloaded from the NOAA website.

3. **Adaptive mesh technique**: In order to simulate tsunami waves more accurately, especially around the initial data, a special adaptive mesh technique is applied to improve the accuracy and efficiency of the simulation.

4. **Actual case study**: The effectiveness of the proposed method is verified by simulating the actual tsunami event that occurred near Java Island in 2006.

5. **Numerical simulation**: Multiple numerical simulation experiments are carried out, including checking the convergence of the code, simulating wave propagation on the flat seabed in the Mediterranean Sea, and simulating wave propagation on the variable - bottom topography near Java Island.

### Formula summary

The main formula involved in the paper is the simplified Boussinesq system (sBBM):

\[
\begin{aligned}
&\eta_t+\nabla\cdot((D + \eta+\zeta)V)+\zeta_t+\tilde{A}\nabla\cdot(D^2\nabla\zeta_t)\\
&+\nabla\cdot\left(AD^2[\nabla(\nabla D\cdot V)+\nabla D\nabla\cdot V]-bD^2\nabla\eta_t\right)=0,\\
&V_t + g\nabla\eta+\frac{1}{2}\nabla|V|^2 + B_gD[\nabla(\nabla D\cdot\nabla\eta)+\nabla D\Delta\eta]-dD^2\Delta V_t - B_D\nabla\zeta_{tt}=0,
\end{aligned}
\]

where:
- \( \eta \) represents the deviation of the free surface,
- \( V=(u, v)^T \) represents the horizontal velocity of the fluid at a certain height,
- \( D(x, y) \) represents the seabed depth,
- \( \zeta(x, y, t) \) represents the seabed deformation,
- \( g = 9.81 \) is the acceleration due to gravity.

### Conclusion

This paper aims to more accurately simulate the generation and propagation processes of tsunami waves through improved numerical methods and adaptive mesh techniques, thereby providing a scientific basis for disaster prediction and prevention and mitigation.