Bulk Flow Parameterization of Temporally Evolving Turbulent Stratified Open Channel Flow Subject to Radiative Heating

The transient response of an initially neutral turbulent open channel subject to sudden stable stratification through radiative surface heating is investigated through direct numerical simulations. We find the convergence of the vertical buoyancy and momentum fluxes toward their respective analytical equilibrium values to be a global monotonic process such that locally, buoyancy and momentum flux equilibrium is obtained simultaneously at all locations within the channel. We present scaling arguments to show that the evolution of the flux convergence ratios toward equilibrium scale directly with Riτ−1/2(t/Tτ) and reach equilibrium at Riτ−1/2(t/Tτ)≈2, consistent with results from the literature, where Riτ is the friction Richardson number, Tτ is the friction time scale, and t is the measured time from the initial isothermal state. We define the upper laminar layer thickness δi and present scaling arguments to show that the laminar layer thickness is well-parameterized by a hybrid bulk parameter such that δi/δ≈0.45(CR1/2PrRiτ−1Reτ)−1/2 across the full external parameter set. Here δ is the channel height, CR is the equilibrium convergence ratio, Pr is the molecular Prandtl number, and Reτ is the friction Reynolds number.