Modification of the k-epsilon scheme and its application for describing turbulence in inland water bodies

<p>A parameterization of the Prandtl number as a function of the gradient Richardson number is proposed in order to correctly take into account stratification when calculating the thermohydrodynamic regime of inland water bodies. This parameterization allows the existence of turbulence at any values ​​of the Richardson number.</p><p>The proposed function is used to calculate the turbulent thermal conductivity coefficient in a k-epsilon mixing scheme. Modification is implemented in the three-dimensional hydrostatic model developed at the Research Computing Center of Moscow State University.</p><p>It is demonstrated that the proposed modification (in contrast to the standard scheme with a constant Prandtl number) leads to smoothing all sharp changes in vertical distributions of turbulent mixing parameters (turbulent kinetic energy, temperature and thickness of the shock layer) and imposes a Richardson number-dependent relation on the empirical constants of k-epsilon turbulent mixing scheme.</p><p>The work was supported by grants of the RF President’s Grant for Young Scientists (MK-1867.2020.5) and by the RFBR (19-05-00249, 20-05-00776). </p>