No Arabic abstract
We present a new physically-motivated parameterization, based on the ratio of Thorpe and Ozmidov scales, for the irreversible turbulent flux coefficient $Gamma_{mathcal M}= {mathcal M}/epsilon$, i.e. the ratio of the irreversible rate ${mathcal M}$ at which the background potential energy increases in a stratified flow due to macroscopic motions to the dissipation rate of turbulent kinetic energy. Our parameterization covers all three key phases (crucially, in time) of a shear-induced stratified turbulence life cycle: the initial, `hot growing phase, the intermediate energetically forced phase, and the final `cold fossilization decaying phase. Covering all three phases allows us to highlight the importance of the intermediate one, to which we refer as the `Goldilocks phase due to its apparently optimal (and so neither too hot nor too cold, but just right) balance, in which energy transfer from background shear to the turbulent mixing is most efficient. $Gamma_{mathcal M}$ is close to 1/3 during this phase, which we demonstrate appears to be related to an adjustment towards a critical or marginal Richardson number for sustained turbulence $sim 0.2-0.25$.
We present experimental measurements of a wall-bounded gravity current, motivated by characterizing natural gravity currents such as oceanic overflows. We use particle image velocimetry and planar laser-induced fluorescence to simultaneously measure the velocity and density fields as they evolve downstream of the initial injection from a turbulent channel flow onto a plane inclined at 10$^circ$ with respect to horizontal. The turbulence level of the input flow is controlled by injecting velocity fluctuations upstream of the output nozzle. The initial Reynolds number based on Taylor microscale of the flow, R$_lambda$, is varied between 40 and 120, and the effects of the initial turbulence level are assessed. The bulk Richardson number $Ri$ for the flow is about 0.3 whereas the gradient Richardson number $Ri_g$ varies between 0.04 and 0.25, indicating that shear dominates the stabilizing effect of stratification. Kelvin-Helmholtz instability results in vigorous vertical transport of mass and momentum. We present baseline characterization of standard turbulence quantities and calculate, in several different ways, the fluid entrainment coefficient $E$, a quantity of considerable interest in mixing parameterization for ocean circulation models. We also determine properties of mixing as represented by the flux Richardson number $Ri_f$ as a function of $Ri_g$ and diapycnal mixing parameter $K_rho$ versus buoyancy Reynolds number $Re_b$. We find reasonable agreement with results from natural flows.
In linearly stable shear flows at moderate Re, turbulence spontaneously decays despite the existence of a codimension-one manifold, termed the edge of chaos, which separates decaying perturbations from those triggering turbulence. We statistically analyse the decay in plane Couette flow, quantify the breaking of self-sustaining feedback loops and demonstrate the existence of a whole continuum of possible decay paths. Drawing parallels with low-dimensional models and monitoring the location of the edge relative to decaying trajectories we provide evidence, that the edge of chaos separates state space not globally. It is instead wrapped around the turbulence generating structures and not an independent dynamical structure but part of the chaotic saddle. Thereby, decaying trajectories need not cross the edge, but circumnavigate it while unwrapping from the turbulent saddle.
Turbulent-laminar intermittency, typically in the form of bands and spots, is a ubiquitous feature of the route to turbulence in wall-bounded shear flows. Here we study the idealised shear between stress-free boundaries driven by a sinusoidal body force and demonstrate quantitative agreement between turbulence in this flow and that found in the interior of plane Couette flow -- the region excluding the boundary layers. Exploiting the absence of boundary layers, we construct a model flow that uses only four Fourier modes in the shear direction and yet robustly captures the range of spatiotemporal phenomena observed in transition, from spot growth to turbulent bands and uniform turbulence. The model substantially reduces the cost of simulating intermittent turbulent structures while maintaining the essential physics and a direct connection to the Navier-Stokes equations. We demonstrate the generic nature of this process by introducing stress-free equivalent flows for plane Poiseuille and pipe flows which again capture the turbulent-laminar structures seen in transition.
The Lagrangian (LA) and Eulerian Acceleration (EA) properties of fluid particles in homogeneous turbulence with uniform shear and uniform stable stratification are studied using direct numerical simulations. The Richardson number is varied from $Ri=0$, corresponding to unstratified shear flow, to $Ri=1$, corresponding to strongly stratified shear flow. The probability density functions (pdfs) of both LA and EA have a stretched-exponential shape and they show a strong and similar influence on the Richardson number. The extreme values of the EA are stronger than those observed for the LA. Geometrical statistics explain that the magnitude of the EA is larger than its Lagrangian counterpart due to the mutual cancellation of the Eulerian and convective acceleration, as both vectors statistically show an anti-parallel preference. A wavelet-based scale-dependent decomposition of the LA and EA is performed. The tails of the acceleration pdfs grow heavier for smaller scales of turbulent motion. Hence the flatness increases with decreasing scale, indicating stronger intermittency at smaller scales. The joint pdfs of the LA and EA indicate a trend to stronger correlations with increasing Richardson number and at larger scales of the turbulent motion. A consideration of the terms in the Navier--Stokes equation shows that the LA is mainly determined by the pressure-gradient term, while the EA is dominated by the nonlinear convection term.
This paper presents a method for calculating the wall shear rate in pipe turbulent flow. It collapses adequately the data measured in laminar flow and turbulent flow into a single flow curve and gives the basis for the design of turbulent flow viscometers. Key words: non-Newtonian, wall shear rate, turbulent, rheometer