A weak coupling algorithm for seabed–wave interaction analysis Original Research Article

This paper proposes a weak coupling algorithm for the seabed–wave interaction analysis. This algorithm is based on the error integration over a time interval and space along fluid–seabed interfaces. With this weak coupling algorithm, seabed–wave interaction analysis can be easily carried out with a fluid solver in fluid domain and a meshless solver in seabed domain independently. The fluid solver is based on a two-step projection method for Navier–Stokes equation. The free surface of seawaters during wave propagation is traced through the volume of fluid (VOF). The meshless solver is based on the Biotʹs consolidation theory of saturated poroelasticity which considers the interaction of soil deformation and pore water pressure in seabed. Data exchange at the interface between fluid and porous medium domains is carried out through a point interpolation method. Interface condition is approximately satisfied in a variational sense through an iteration scheme between two solvers. This algorithm does allow not only non-matching mesh along the interface but also non-matching time step in different subdomains. Numerical examples show that this weak coupling algorithm is able to capture the deformation of interface and improve the computation accuracy and efficiency.