I'm not familiar with IH2VOF, but in most environmental-scale models turbulent motions are not resolved, and the diffusive effect of turbulence is parameterized by a closure (such as k-epsilon). If you turn the turbulence closure off, you could still (very roughly) capture turbulent mixing by...