H. Karjoun,A. Beljadid
Abstract:In this study, a numerical model preserving a class of nontrivial steady-state solutions is proposed to predict waves propagation and waves run-up on coastal zones. The numerical model is based on the Saint-Venant system with source terms due to variable bottom topography and bed friction effects. The resulting nonlinear system is solved using a Godunov-type finite volume method on unstructured triangular grids. A special piecewise linear reconstruction of the solution is implemented with a correction technique to ensure the accuracy of the method and the positivity of the computed water depth. Efficient semi-implicit techniques for the friction terms and a well-balanced formulation for the bottom topography are used to exactly preserve stationary steady-state s solutions. Moreover, we prove that the numerical scheme preserves a class of nontrivial steady-state solutions. To validate the proposed numerical model against experiments, we first demonstrate its ability to preserve nontrivial steady-state solutions and then we model several laboratory experiments for the prediction of waves run-up on sloping beaches. The numerical simulations are in good agreement with laboratory experiments which confirms the robustness and accuracy of the proposed numerical model in predicting waves propagation on coastal areas.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to develop a numerical model capable of accurately predicting the waves run - up phenomenon in coastal areas. Specifically, the model aims to:
1. **Retain non - trivial steady - state solutions**: Ensure that the numerical model can precisely maintain a class of non - trivial steady - state solutions, which is crucial for simulating complex coastal topography and hydrodynamic flow.
2. **Predict wave propagation and run - up**: Predict the wave propagation and run - up phenomena in coastal areas by solving the Saint - Venant system of equations that takes into account variable seabed topography and bed friction effects.
### Key points of the numerical model
- **Mathematical model basis**: Based on the Saint - Venant system of equations (including source terms) that takes into account variable seabed topography and bed friction effects. These equations are nonlinear and are solved on unstructured triangular meshes using the Godunov - type finite - volume method.
The Saint - Venant system of equations is:
\[
\begin{cases}
\frac{\partial h}{\partial t}+\frac{\partial(hu)}{\partial x}+\frac{\partial(hv)}{\partial y} = 0,\\
\frac{\partial(hu)}{\partial t}+\frac{\partial\left(hu^{2}+\frac{g}{2}h^{2}\right)}{\partial x}+\frac{\partial(huv)}{\partial y}=-gh\frac{\partial B}{\partial x}-C_{f}u\sqrt{u^{2}+v^{2}},\\
\frac{\partial(hv)}{\partial t}+\frac{\partial(huv)}{\partial x}+\frac{\partial\left(hv^{2}+\frac{g}{2}h^{2}\right)}{\partial y}=-gh\frac{\partial B}{\partial y}-C_{f}v\sqrt{u^{2}+v^{2}},
\end{cases}
\]
where \(h\) is the water depth, \(u\) and \(v\) are the velocity components along the \(x\) and \(y\) directions respectively, \(g\) is the gravitational acceleration, \(B(x,y)\) is the seabed topography, and \(C_{f}\) is the bed friction coefficient.
- **Numerical method**: Adopt the central - upwind finite - volume method, combined with an efficient semi - implicit technique to handle the friction terms and a balanced discretization to handle the bottom - topography source terms, in order to ensure the stability and accuracy of the model.
- **Retention of steady - state solutions**: A piecewise - linear reconstruction technique is specially designed to ensure that the calculated water depth is non - negative, and it is proved that this numerical scheme can maintain a class of non - trivial steady - state solutions, especially the steady - state flow on an inclined surface.
### Verification and application
To verify the proposed numerical model, the author first demonstrates its ability to maintain non - trivial steady - state solutions, and then models the waves run - up through multiple laboratory experiments. The numerical simulation results match the experimental data well, verifying the robustness and accuracy of this numerical model in predicting wave propagation in coastal areas.
In conclusion, this research proposes a numerical model that can effectively predict the wave propagation and run - up phenomena on complex coastal topographies, which has important practical application value.