Animations comparing surface wave simulations with standard and "stabilized" two-equation turbulence models

It is emphasized that all of the simulations shown in these animations (as well as all results shown in Larsen and Fuhrman, 2018, as indicated on their p. 434) have the buoyancy production term active in the k-equation, as suggested recently in this context by e.g. Devolder et al. (2017). This only creates a local sink in the turbulent kinetic energy near the air-water interface, but clearly does not stabilize the models. This is because the problem of over-production of turbulence in two-equation models is not confined to the near-surface region; Rather the entirety of the nearly-potential flow region beneath the surface waves is unstable, as again proved unconditionally by Larsen and Fuhrman (2018). It is the simple modification of the eddy viscosity, coupled with its elimination from either the \omega or \epsilon production term, which formally stabilizes the models, as shown therein.