Investigation on the Use of a Multiphase Eulerian CFD Solver to Simulate Breaking Waves

The main challenge in CFD multiphase simulations of breaking waves is the wide range of interfacial length scales occurring in the flow: from the free surface measurable in meters down to the entrapped air bubbles with size of a fraction of a millimeter. This paper presents a preliminary investigation on a CFD model capable of handling this problem. The model is based on a solver, available in the open-source CFD toolkit OpenFOAM, which combines the Eulerian multi-fluid approach for dispersed flows with a numerical interface sharpening method. The solver, enhanced with additional formulations for mass and momentum transfer among phases, was satisfactorily tested against an experimental bubble column flow. The model was then used to simulate the propagation of a laboratory solitary breaking wave. The motion of the free surface was successfully reproduced up to the breaking point. Further implementations are needed to simulate the air entrainment phenomenon.