A three-dimensional parallel CIP-based model is presented in this paper, which is accelerated by parallel computation based on an open Message Passing Interface (OpenMPI) algorithm and open Multi-Processing (OpenMP). The Navier-Stokes equations are solved by the projection method. A constrained interpolation profile (CIP) method is applied for the convection term. The HYPRE library is adopted for the Poisson equation. The tangent of hyperbola for interface capturing with slope weighting (THINC/SW) method is used to simulate the variation of the free surface. The benchmarks are the flow over the sphere and the interaction between the dam-break wave and a stationary square cylinder. The drag coefficient, free-surface elevations and the force on the cylinder can be accurately predicted. A convincing agreement between the numerical predictions and the published data is obtained. Numerical investigation of three-dimensional problems based on a reliable model can help us better understand the complex physical phenomena and provide a reference for practical engineering.


Wave-structure interaction is an important topic in many scientific and engineering applications, e.g. offshore platforms, marine renewable energy and disaster prevention and mitigation engineering. The evolution process of dam-break will induce shock waves. It can be used to simulate the interaction between disaster waves and structures. The flow physics of wave-structure interaction can guide the safety design of ocean engineering. However, there are still many challenges to be addressed. Firstly, the theoretical method simplifies the interaction process between wave and structure. However, the process is strongly nonlinear. Secondly, the experimental method has a longer period to obtain the whole flow field and the slamming information of the structure. Finally, slamming and air entrainment are still a challenge to the stability and accuracy of the numerical model.

Numerical simulation of the interaction between dam-break waves and structures is studied by several scholars. De et al. (2013) used the SPH method to simulate two-dimensional nonlinear wave-structure interaction problems. Dam-break flood waves in a dry channel with a hump were investigated by Ozmen-Cagatay et al. (2014) and only twodimensional models are used. The SPH method was adopted for the sequential dam-break wave problem by Kocaman and Dal (2020). Five successive stages were observed in the downstream reservoirs: dambreak wave propagation, overtopping, reflection wave, run-up, and oscillations. These scholars have done a lot of research. However, it is still not easy to study three-dimensional problems. In fact, the hydrodynamic characteristics of three dimensions are different from those of two dimensions. Moreover, the three-dimensional characteristics of wave structure interaction are obvious. Therefore, it is necessary to develop a three-dimensional numerical model.

This content is only available via PDF.
You can access this article if you purchase or spend a download.