Turbulence modeling by use of the renormalization group （RNG） κ-ε model for Reynolds-stress closure is carried out to reveal the evolution dynamics for lock release gravity currents with the so-called slumping, inviscid and viscous phases. Field evolution of the turbulent current is investigated, and time transition of global energy balance is presented between the terms of potential energy, averaged kinetic energy, turbulent kinetic energy, turbulent dissipation and viscous dissipation. It is well illustrated that turbulent dissipation and viscous force are respectively dominant in the inviscid and viscous phases, while inertia effect accounts for the slumping.
The motion of a lock-release oil slick as an immiscible two-fluid gravity current is numerically studied by a finite dif ference algorithm based on the volume of fluid （VOF） method for the basic formulation and a rigid cover approximation for the open free surface. Detailed numerical simulation with careful model validation reveals the existence of turbulence and the adaptability of the renormalization group （RNG） k - ε model for the Reynolds-stress closure in the case of the oil slick. The time evolution and spatial distribution of the mean velocity, turbulence kinetic energy and turbulent viscosity are characterized. The mechanism for the transition from an initial gravity-inertial phase to a second gravity-vinous phase is shown to be the relaminarization effect of the initially highly turbulent slick. Compared well with known theoretical analyses and experimental observations, the turbulence modeling results in self-similar spreading laws in terms of the fact that the oil slick passes through the initial gravity-inertial phase with the front speed decreasing as t ^-1/3（where t is the time measured from lock release） and the second gravity-viscous phase with the front speed decreasing as t^-5/8.
