Transition to stably stratified states in open channel flow with radiative surface heating

Abstract

Direct numerical simulations of turbulent stratified flow in an open channel with an internal heat source following the Beer-Lambert law from the surface are used to investigate the transition from neutral to strongly stable flow. Our buoyancy bulk parameter is defined through the ratio of the domain height δ to L, a bulk Obukhov length scale for the flow. We cover the range λ = δ/L = 0−2.0, from neutral conditions to the onset of the stable regime, with the Reynolds number range Reτ = 200 − 800, at a Prandtl number of 0.71. The result is a boundary layer flow where the effects of stratification are weak in the wall region but progressively stronger in the outer layer up to the free surface. At λ 1 the turbulent kinetic energy budget is in local equilibrium over a region extending from the near wall region to a free surface affected region a distance lν from the surface with lν/δ ∼ Re−1/2. In this equilibrium region the flow can be characterised by the flux Richardson number Rf and the local Obukhov length scale Λ. At higher λ local mixing limit conditions are observed over an extended region. At λ = 2 the flux Richardson number approaches critical limit values of Rf,c 0.18 and gradient Richardson number Ric 0.2. At high λ, we obtain a flow field where buoyancy interacts with the smallest scales of motion the turbulent shear stress and buoyancy flux are suppressed to molecular levels. We find this regime can be identified in terms of the parameter ReL,c = Luτ /ν 200 − 400 which is related to the L∗ parameter of Flores & Riley (2011) and buoyancy Reynolds number R. With energetic equilibrium attained, the local buoyancy Reynolds number ReΛ = Λ u w 1/2/ν, is directly related to separation of the Ozmidov (lO) and Kolmogorov (η) length scales in the outer boundary layer by ReΛ R ≡ (lO/η)4/3. The inner wall region has the behaviour R ∼ ReL Reτ in contrast to stratified boundary layer flows where the buoyancy flux is non-zero at the wall and R ∼ ReL .