A code is developed to simulate incompressible free surface flows using the Roe's fluxdifference splitting scheme. An interface of two fluids is considered as a moving contact discontinuity. The continuities of pressure and normal velocity across the interface are enforced by the conservation law in the integral sense. The fluxes are computed using the Roe's flux-difference splitting scheme for two incompressible fluids. The interface can be identified based on the computed density distribution. However, no additional treatment is required along the interface during the whole computations. Complicated time evolution of the interface including topological change can be captured without any difficulties. The developed code is applied to simulate the Rayleigh-Taylor instability of two incompressible fluids in the density ratio of 7.2:1 and the broken dam problem of water-air. The present results are compared with other available results and good agreements are achieved for the both cases.

SUMMARYIn this study, a method is developed to simulate the interaction between free surface flows and moving or deforming boundaries using the flux-difference splitting scheme on the hybrid Cartesian/immersed boundary method. At each physical time step, the boundary is defined by an unstructured triangular surface grid. Immersed boundary (IB) nodes are distributed inside an instantaneous fluid domain based on edges crossing the boundary. At an IB node, dependent variables are reconstructed along the local normal line to the boundary. Inviscid fluxes are computed using Roe's flux-difference splitting scheme for immiscible and incompressible fluids. The free surface is considered as a contact discontinuity in the density field. The motion of free surface is captured without any additional treatment along the fluid interface. The developed code is validated by comparisons with other experimental and computational results for a piston-type wave maker, impulsive motion of a submerged circular cylinder, flow around a submerged hydrofoil, and Rayleigh-Taylor instability. The developed code is applied to simulate wave generation due to a continuously deforming bed beneath the free surface. The violent motion of a free surface caused by sloshing in a spherical tank is simulated. In this case, the free surface undergoes breakup and reconnection.