A Full Eulerian Wave-Body Interaction Solver for Incompressible Two-Phase Flows with High Density and Viscosity Ratios

A force-free immersed boundary method is developed based on a primitive variable formulation of the Navier-Stokes equations, using Chorin's projection method. In the proposed method, the predetermined motion of a rigid structure is considered, and the rigid domain velocity field is imposed on the formulation by patching the rigid-body motion onto the fluid domain velocity field in the predictor step of the formulation. The capability of the method in handing two-phase
fluid interaction with a structure is examined. For a two-phase uid with large density and viscosity ratios, the Laplacian operator becomes discontinuous. The interface evolution under the effect of an oscillating heaving cylinder based on the setup presented in [1] is considered. Before examining the method's accuracy in solving fluid-structure interaction (FSI) problems comprising two uids with high density and viscosity ratio, the classical problem of an oscillating cylinder in a quiescent fluid is considered, and the results are compared with [2] and [3].