Solution of 2D Boussinesq systems with FreeFem++: The flat bottom case
This work provides a numerical tool for simulating coastal wave dynamics, but it is an incremental implementation of existing methods.
The authors developed a FreeFem++ code to solve 2D Boussinesq systems for long surface waves, using P1 finite elements and a second-order Runge-Kutta scheme with mesh adaptation. Their results match those in the literature.
We consider here different family of Boussinesq systems in two space dimensions. These systems approximate the three-dimensional Euler equations and consist of three coupled nonlinear dispersive wave equations that describe propagation of long surface waves of small amplitude in ideal fluids over a horizontal bottom. We present here a FreeFem++ code aimed at solving numerically these systems where a discretization using P1 finite element for these systems was taken in space and a second order Runge-Kutta scheme in time. We give the detail of our code where we use a mesh adaptation technique. An optimization of the used algorithm is done and a comparison of the solution for different Boussinesq family is done too. The results we obtained agree with those of the literature.