Multiphase turbulence modeling of the flow in the wake of a transom stern

The objective of this effort is to develop and assess multiphase turbulence closure models for incompressible highly variable density turbulent (IHVDT) flows such as the two-phase flow in the wake of a transom stern. These flows, which have an Atwood number $At=(\rho_{2} - \rho_{1})/(\rho_{2} + \rho_{1})\approx 1$, are characterized by significant turbulent mass flux for which there is little guidance in turbulence closure modeling for both the momentum and the scalar transport. In this work, high-resolution numerical simulations are performed on the wake of a canonical transom stern at large scales using conservative Volume-of-Fluid (cVOF) and implicit Large Eddy Simulation (iLES). Boundary Data Immersion Method (BDIM) is used to simulate the dry transom stern wake region at three different Froude numbers and two different effective viscosities. Analysis of the simulation results for the turbulent anisotropy, turbulent kinetic energy and turbulent mass flux budget, as well as \emph{a priori} closure model testing will be presented.