Solution of 2D Boussinesq systems with freefem++: the flat bottom case

Abstract

Abstract -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 and which was studied in [7,9,10].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.