No Arabic abstract
Three-dimensional laminar flow structures with mixing, chemical reaction, normal strain, and shear strain qualitatively representative of turbulent combustion at the small scales are analyzed. A mixing layer is subjected to counterflow in the transverse y- and z-directions. Both non-reactive and reactive flows are examined. Reduction of the three-dimensional boundary-layer equations to a one-dimensional similar form is obtained allowing for heat and mass diffusion with variations in density and properties. In steady configurations, a set of ODEs governs the three velocity components as well as the scalar-field variables. A flamelet model for individual diffusion flames with combined shear and normal strain is developed. Another model with solution in similar form is obtained for a configuration with a dominant diffusion flame and a weaker fuel-rich premixed flame. Results for the velocity and scalar fields are found for ranges of Damkohler number Da, normal strain rate due to the counterflow, streamwise-velocity ratio across the mixing layer, Prandtl number, and Mach number. For the flamelet model, a conserved scalar is cast as the independent variable to give an alternative description of the results. The imposed normal strain decreases mixing-layer thickness and increases scalar gradients and transport rates. There is indication of diffusion control for partially premixed flames in the multi-branched flame situation. The enhancement of the mixing and combustion rates by imposed normal strain on a shear layer can be very substantial. Also, the imposition of shear strain and thereby vorticity on the counterflow can be substantial indicating the need for flamelet models with both shear strain and normal strain.
Previous works for a liquid suddenly contacting a gas at a supercritical pressure show the coexistence of both phases and the generation of diffusion layers on both sides of the liquid-gas interface due to thermodynamic phase equilibrium. A related numerical study of a laminar mixing layer between the liquid and gas streams in the near field of the splitter plate suggests that mass, momentum and thermal diffusion layers evolve in a self-similar manner at very high pressures. In this paper, the high-pressure, two-phase, laminar mixing-layer equations are recast in terms of a similarity variable. A liquid hydrocarbon and gaseous oxygen are considered. Freestream conditions and proper matching conditions at the liquid-gas interface are applied. To solve the system of equations, a real-fluid thermodynamic model based on the Soave-Redlich-Kwong equation of state is selected. A comparison with results obtained by directly solving the laminar mixing-layer equations shows the validity of the similarity approach applied to non-ideal two-phase flows. Even when the gas is hotter than the liquid, condensation can occur at high pressures while heat conducts into the liquid. Finally, a generalized correlation is proposed to represent the evolution of the mixing layer thickness for different problem setups.
Numerical analysis of a shear layer between a cool liquid n-decane hydrocarbon and a hot oxygen gas at supercritical pressures shows that a well-defined phase equilibrium can be established. Variable properties are considered with the product of density and viscosity in the gas phase showing a nearly constant result within the laminar flow region with no instabilities. Sufficiently thick diffusion layers form around the liquid-gas interface to support the case of continuum theory and phase equilibrium. While molecules are exchanged for both species at all pressures, net mass flux across the interface shifts as pressure is increased. Net vaporization occurs for low pressures while net condensation occurs at higher pressures. For a mixture of n-decane and oxygen, the transition occurs around 50 bar. The equilibrium values at the interface quickly reach their downstream asymptotes. For all cases, profiles of diffusing-advecting quantities collapse to a similar solution (i.e., function of one independent variable). Validity of the boundary layer approximation and similarity are shown in both phases for Reynolds numbers greater than 239 at 150 bar. Results for other pressures are also taken at high Reynolds numbers. Thereby, the validity of the boundary layer approximation and similarity are expected. However, at very high pressures, the similar one-dimensional profiles vary for different problem constraints.
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.
On its way to turbulence, plane Couette flow - the flow between counter-translating parallel plates - displays a puzzling steady oblique laminar-turbulent pattern. We approach this problem via Galerkin modelling of the Navier-Stokes equations. The wall-normal dependence of the hydrodynamic field is treated by means of expansions on functional bases fitting the boundary conditions exactly. This yields a set of partial differential equations for the spatiotemporal dynamics in the plane of the flow. Truncating this set beyond lowest nontrivial order is numerically shown to produce the expected pattern, therefore improving over what was obtained at cruder effective wall-normal resolution. Perspectives opened by the approach are discussed.
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$.