Modeling and parallel computation of the non-linear interaction of rigid bodies with incompressible multi-phase flow

A computational tool is developed to capture the interaction of solid object with two-phase flow. The full two-dimensional Navier–Stokes equations are solved on a regular structured grid to resolve the flow field. The level set and the immersed boundary methods are used to capture the free surface of a fluid and a solid object, respectively. A two-step projection method along with Multi-Processing (OpenMP) is employed to solve the flow equations. The computational tool is verified based on numerical and experimental data with three scenarios: a cylinder falling into a rectangular domain due to gravity, transient vertical oscillation of a cylinder by releasing above its equilibrium position, and a dam breaking in the presence of a fixed obstacle. In the first two validation simulations, the accuracy of the immersed boundary method is verified. However the accuracy of the level set method while the computational tool can model the high density ratio is confirmed in the dam breaking simulation. The results obtained from the current method are in good agreement with experimental data and other numerical studies. The applicability of the current computational tool for the interaction of a buoy in a water wave tank with two types of waves; symmetrical and asymmetrical waves; has also been studied.