SEARCH · Engineering Papers

Results for “Turbulence Stabilization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Advancing the physics basis for quasi-helically symmetric stellarators

A new optimized quasi-helically symmetric configuration is described that has the desirable properties of improved energetic particle confinement, reduced turbulent transport by three-dimensional shaping and non-resonant divertor capabilities. The configuration presented in this paper is explicitly optimized for quasi-helical symmetry, energetic particle confinement, neoclassical confinement and stability near the axis. Post optimization, the configuration was evaluated for its performance with regard to energetic particle transport, ideal magnetohydrodynamic stability at various values of plasma pressure and ion temperature gradient instability induced turbulent transport. The effects of discrete coils on various confinement figures of merit, including energetic particle confinement, are determined by generating single-filament coils for the configuration. Preliminary divertor analysis shows that coils can be created that do not interfere with expansion of the vessel volume near the regions of outgoing heat flux, thus demonstrating the possibility of operating a non-resonant divertor.

Physics↗

DOE OSTI2026

Delving into the depths of NGC 3783 with XRISM V. Broad-band modelling of ionized outflows

The Seyfert 1 galaxy NGC 3783 hosts a multi-phase warm absorber (WA) that has been extensively studied in the X-ray band. High-resolution spectra from 2000−2001 revealed a complex outflow with multiple ionization and velocity components. Two decades later, new XMM-Newton and XRISM observations allow us to investigate the long-term evolution of these outflows. We performed joint spectral modelling of the XMM-Newton /RGS and XRISM/Resolve time-averaged spectra using the pion photoionization code within SPEX. We derived the ionization parameter, column density, turbulent velocity, and outflow velocity for each absorption component, and investigated their thermal stability and absorption measure distribution (AMD) to characterize the physical and dynamical properties of the WA in NGC 3783 in 2024. We compared these results with the 2000−2001 epoch to assess long-term variability, stability, and possible changes in the absorber population. We identify eight WA components spanning log ξ = 1.08 − 3.38 and outflow velocities of 480−1230 km s −1 . The ranges of column densities and turbulent velocities remain broadly consistent with the WAs from 2000−2001, but the earlier data contained more low-ionization, high-velocity components. The total column density of all the outflows in 2024 is 1.5 times larger than in 2000−2001, which means that it has been replenished by fresh material. The dominant unresolved transition array (UTA) absorber (component B3) has increased its column density by a factor of three while maintaining a similar ionization parameter. The WAs in NGC 3783 have undergone significant structural and dynamical evolution over the past 24 years.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

PeleMP: The Multiphysics Solver for the Combustion Pele Adaptive Mesh Refinement Code Suite

Combustion encompasses multiscale, multiphase reacting flow physics spanning a wide range of scales from the molecular scales, where chemical reactions occur, to the device scales, where the turbulent flow is affected by the geometry of the combustor. This scale disparity and the limited measurement capabilities from experiments make modeling combustion a significant challenge. Recent advancements in high-performance computing (HPC), particularly with the Department of Energy's Exascale Computing Project (ECP), have enabled high-fidelity simulations of practical applications to be performed. The major physics submodels, including chemical reactions, turbulence, sprays, soot, and thermal radiation, exhibit distinctive computational characteristics that need to be examined separately to ensure efficient utilization of computational resources. This paper presents the multiphysics solver for the Pele code suite, called PeleMP, which consists of models for spray, soot, and thermal radiation. Here, the mathematical and algorithmic aspects of the model implementations are described in detail as well as the verification process. The computational performance of these models is benchmarked on multiple supercomputers, including Frontier, an exascale machine. Results are presented from production simulations of a turbulent sooting ethylene flame and a bluff-body swirl stabilized spray flame with sustainable aviation fuels to demonstrate the capability of the Pele codes for modeling practical combustion problems with multiphysics. This work is an important step toward the exascale computing era for high-fidelity combustion simulations providing physical insights and data for predictive modeling of real-world devices.

42 ENGINEERING↗

DOE OSTI2020

Efficacy of the Cell Perturbation Method in Large-Eddy Simulations of Boundary Layer Flow over Complex Terrain

A challenge to simulating turbulent flow in multiscale atmospheric applications is the efficient generation of resolved turbulence motions over an area of interest. One approach is to apply small perturbations to flow variables near the inflow planes of turbulence-resolving simulation domains nested within larger mesoscale domains. While this approach has been examined in numerous idealized and simple terrain cases, its efficacy in complex terrain environments has not yet been fully explored. Here, we examine the benefits of the stochastic cell perturbation method (CPM) over real complex terrain using data from the 2017 Perdigão field campaign, conducted in an approximately 2-km wide valley situated between two nearly parallel ridges. Following a typical configuration for multiscale simulation using nested domains within the Weather Research and Forecasting (WRF) model to downscale from the mesoscale to a large-eddy simulation (LES), we apply the CPM on a domain with horizontal grid spacing of 150 m. At this resolution, spurious coherent structures are often observed under unstable atmospheric conditions with moderate mean wind speeds. Results from such an intermediate resolution grid are often nested down for finer, more detailed LES, where these spurious structures adversely affect the development of turbulence on the subsequent finer grid nest. We therefore examine the impacts of the CPM on the representation of turbulence within the nested LES domain under moderate mean flow conditions in three different stability regimes: weakly convective, strongly convective, and weakly stable. In addition, two different resolutions of the underlying terrain are used to explore the role of the complex topography itself in generating turbulent structures. We demonstrate that the CPM improves the representation of turbulence within the LES domain, relative to the use of high-resolution complex terrain alone. During the convective conditions, the CPM improves the rate at which smaller-scales of turbulence form, while also accelerating the attenuation of the spurious numerically generated roll structures near the inflow boundary. During stable conditions, the coarse mesh spacing of the intermediate LES domain used herein was insufficient to maintain resolved turbulence using CPM as the flow develops downstream, highlighting the need for yet higher resolution under even weakly stable conditions, and the importance of accurate representation of flow on intermediate LES grids.

17 WIND ENERGY↗

DOE OSTI2020

Boundary-consistent B-spline filtering schemes and application to high-fidelity simulations of turbulence

A filtering operation, based on B-spline discretizations, is introduced to target weakly growing mesh-scale oscillations that can arise in high-fidelity turbulence simulations. This is a spectral regularization that can be described using the singular values of a banded matrix operator, with the filtering strength set by a scalar- or vector-valued penalty parameter. The penalty parameter can be specified though it can also be advantageously selected to minimize the generalized cross validation (GCV) measure of distance between the pre- and post-filtered solutions. Efficient algorithms are developed to compute both the scalar and vector penalty parameters. The B-spline filter has a sharper localization to high-wavenumber than compact or explicit filters of the same stencil width and is demonstrated for solutions of the Burgers' equation, decaying Burgers' turbulence, and compressible Navier–Stokes turbulent channel flow. Furthermore, these simulations confirm the scheme's numerical stability and ability to narrowly target the high wavenumber components of numerical solutions. An advantage over finite-difference filters is that these B-spline filters are stable on bounded domains and even preserve formal order of accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Modulation of North Atlantic Polar Low Activity and Associated Flow Patterns by Sudden Stratospheric Warmings

Abstract Polar lows (PLs) are intense maritime mesoscale cyclones that can pose hazards to coastal communities and marine operation in the Arctic. This study examines the impacts of sudden stratospheric warmings (SSWs) on PL activity in the subarctic North Atlantic. The 20 days following SSWs are characterized by tropospheric circulation anomalies resembling the negative phase of the North Atlantic Oscillation. PL activity decreases significantly over the Labrador Sea, which can be attributed to the infrequent occurrence of low static stability and strong environmental baroclinicity, as well as reduced surface turbulent heat fluxes. These results suggest that a skillful prediction of SSWs can improve the extended-range forecast of PL activity over the Labrador Sea. For the Nordic seas, the results imply that the spatial structure of an SSW event is important for the PL modulation through different tropospheric circulation patterns. Situations with increased PL frequency in the Nordic seas are characterized by SSWs centered close to northern Greenland occurring over a smaller area, and a tropospheric response featuring enhanced cold-air outbreaks over the Norwegian Sea. Conversely, PL activity is suppressed over the Nordic seas when the SSW favors the formation of a tropospheric anticyclone above Greenland and Scandinavia. Significance Statement This study investigates the relationships between polar lows (PLs) and sudden stratospheric warmings (SSWs) over the subarctic North Atlantic. A better understanding of the effect of SSWs on PL development has the potential to improve extended-range forecasts of PLs. It is shown that SSWs are responsible for the significantly suppressed regional PL activity over the Labrador Sea, suggesting that SSWs can serve as a predictor for the extended-range forecast of PLs over this region. Following SSW events, the thermodynamic state of atmosphere becomes more stable over the Labrador Sea and hinders the convective development of PLs. For the northern Nordic seas, the impacts of SSWs on PL activity are sensitive to the spatial structure of stratospheric warming.

Chang, Chuan-Chieh↗

DOE OSTI2025

A High-Order Discontinuous Galerkin Spectral Element Method for Compressible Reacting Flows

High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large eddy simulations due to their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reactive Euler equations encountered in high-speed combustion. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of DG approach. Thus, the framework is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. The numerical method is implemented within the spectral element solver Nek5000. Validation cases are conducted for both non-reactive and reactive discontinuous flows to demonstrate the solver capability. In particular, canonical one-dimensional and two-dimensional detonation simulations are performed and the high-order numerical results are validated against available literature data.

computational fluid dynamics (CFD)↗

DOE OSTI2022

The Two-Energies Turbulence Scheme Coupled to the Assumed PDF Method

An update of the two-energy turbulence scheme is presented, the 2TE + APDF scheme. The original version of the two-energy scheme is able to successfully model shallow convection without the need of an additional parameterization for non-local fluxes. However, the performance of the two-energy scheme is worse in stratocumulus cases, where it tends to overestimate the erosion of the stable layers. We have identified the causes: the non-local stability parameter does not consider local stratification, the scheme lacks an internal parameter that could distinguish between a shallow convection regime and a stratocumulus regime, and it uses an inflexible turbulence length scale formulation. To alleviate this problem, we propose several modifications: an update of the stability parameter, a modified computation of the turbulence length scale, and the introduction of the entropy potential temperature to distinguish between a shallow convection and a stratocumulus regime. In addition, the two-energy scheme is coupled to a simplified assumed probability density function method in order to achieve a more universal representation of the cloudy regimes. The updated turbulence scheme is evaluated for several idealized cases and one selected real case in the ICOsahedral Nonhydrostatic (ICON) modeling framework. The results show that the updated scheme corrects the overmixing problem in the stratocumulus cases. The performance of the updated scheme is comparable to the operational setup, and can be thus used instead of the operational turbulence and shallow convection scheme in ICON. Additionally, the updated scheme improves the coupling with dynamics, which is beneficial for the modeling of coherent flow structures in the atmospheric boundary layer.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Stability and transport of gyrokinetic critical pedestals

Abstract A gyrokinetic threshold model for pedestal width–height scaling prediction is applied to multiple devices. A shaping and aspect ratio scan is performed on National Spherical Torus Experiment (NSTX) equilibria, finding Δ ped = 0.92 A 1.04 κ − 1.24 0.38 δ β θ , ped 1.05 for the wide-pedestal branch with pedestal width Δ ped , aspect ratio A , elongation κ , triangularity δ , and normalized pedestal height β θ , ped . The width–transport scaling is found to vary significantly if the pedestal height is varied either with a fixed density or fixed temperature, showing how fueling and heating sources affect the pedestal density and temperature profiles for the kinetic-ballooning-mode (KBM) limited profiles. For an NSTX equilibrium, at fixed density, the wide branch is Δ ped = 0.028 ( q e / Γ e − 1.7 ) 1.5 ∼ η e 1.5 and at fixed temperature Δ ped = 0.31 ( q e / Γ e − 4.7 ) 0.85 ∼ η e 0.85 , where q e and Γ e are turbulent electron heat and particle fluxes and η e = ∇ ln ⁡ T e / ∇ ln ⁡ n e for an electron temperature T e and density n e . Pedestals close to the KBM limit are shown to have modified turbulent transport coefficients compared to the strongly driven KBMs. The role of flow shear is studied as a width–height scaling constraint and pedestal saturation mechanism for a standard and lithiated wide pedestal discharge. Finally, the stability, transport, and flow shear constraints are combined and examined for an NSTX experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Floquet stability analyses of stratified oscillating boundary layers on adiabatic slopes

The presence of a no-slip, impermeable, adiabatic, sloped boundary in an otherwise quiescent, stably stratified, Boussinesq flow generates baroclinic vorticity within a diffusive boundary layer. Such conditions are typical of the oscillating boundary layers on adiabatic abyssal slopes, sloped lake bathymetry and sloped coastal bathymetry in the absence of high-wavenumber internal waves, mean flows, far-field turbulence on larger scales, and resonant tidal–bathymetric interaction. We investigate the linear stability of the oscillating flow within non-dimensional parameter space typical of the $M_2$ tide and hydraulically smooth, middle-latitude abyssal slopes through Floquet linear stability analysis. The flow dynamics depends on three non-dimensional variables: the Reynolds number for Stokes’ second problem, the Prandtl number, and a frequency ratio that accounts for the resonance conditions ($C$, criticality) of the buoyant restoring force and the tidal forcing. The Floquet analysis results suggest that oscillating laminar boundary layers on adiabatic abyssal slopes are increasingly unstable as Reynolds number, criticality parameter and/or spanwise disturbance wavenumber are increased. We also show that the two-dimensional Floquet linear instability necessarily generates three-dimensional baroclinic vorticity, which suggests that the evolution of the gravitational instabilities may be nonlinear as $t→∞$.

42 ENGINEERING↗

DOE OSTI2025

Towards stabilization of pressurized oxy-coal diffusion flames

Designing an effective burner is vital for the development of pressurized oxy-fuel combustion technologies. Turbulent jet diffusion burners are adopted for the pressurized oxy-combustor in this research and a bluff-body is employed to create a recirculation zone to stabilize the flame in pressurized oxy-combustor. The objective of this numerical work is to perform a systematic analysis of the characteristics of this pressurized non-premixed flame. Specifically, a 15-bar POFC combustor of power 100 kWth is modeled employing the Ansys FLUENT commercial platform, using Reynolds-averaged Navier-Stokes (RANS) modeling. The present computational work focuses on identifying the aerothermodynamic features of the isothermal and reacting flow with a disk-shape bluff-body. Results show that fuel-oxidizer stream momentum ratio has great impacts on the temperature profile of the down-fired, co-axial, pressurized oxy-coal, diffusion flame. Parametric study of the width of the burner tip identifies an optimal range for this pressurized burner.

Li, Lei↗

DOE OSTI2022

Large-eddy simulation of non-vaporizing sprays using the spectral-element method

Predictive simulations of high-pressure sprays require accurate representation of the turbulent gaseous flow field generated by liquid jet. Typically, the accuracy that can be obtained with low-order numerical methods (e.g. finite volume, finite element) is limited by stability issues in fine grids and the order of convergence of the method. In this work, we resolve the turbulent flow field in an Eulerian manner using the high-order spectral element method, coupled with a Lagrangian parcels approach to model the atomizing liquid jet. Large eddy simulations of single-hole sprays under non-evaporative conditions were conducted and compared against experimental data from Margot et al. (2008) and Spray A data from the Engine Combustion Network. The sensitivity of liquid penetration and droplet sizes to different breakup model parameters was studied. The effect of different numerical parameters, such as polynomial order of the solution (grid resolution), on liquid penetration was also analyzed. The method achieved grid-independent results using p-refinement, achieving finer resolution (by a factor of x 1.7 - x 3.5) in the gas-phase solution than in state-of-the-art simulations using the finite-volume method. Here the results showed good agreement with experimental data, demonstrating the ability of the current method to accurately capture liquid penetration and the shape of the spray.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Reynolds-Averaged Turbulence Modeling Using Deep Learning with Local Flow Features: An Empirical Approach

Reynolds-Averaged Navier-Stoke (RANS) models offer an alternative avenue in predicting flow characteristics when the corresponding experiments are difficult to achieve due to geometry complexity, limited budget, or knowledge. RANS models require the knowledge of subgrid scale physics to solve conservation equations for mass, energy, and momentum. Mechanistic turbulence models, such as k-epsilon, are generally evaluated and calibrated for specific flow conditions with various degrees of uncertainty. These models have limited capability to assimilate a substantial amount of data due to model form constraints. Meanwhile, deep learning (DL) has been proven to be universal approximators with the potential to assimilate available, relevant, and adequately evaluated data. Moreover, deep neural networks (DNNs) can create surrogate models without knowing function forms. Such a data-driven approach can be used in updating fluid models based on observations as opposed to hard-wiring models with precalibrated correlations. The paper presents progress in applying DNNs to model Reynolds stress using two machine learning (ML) frameworks. A novel flow feature coverage mapping is proposed to quantify the physics coverage of DL-based closures. It can be used to examine the sufficiency of training data and input flow features for data-driven turbulence models. The case of a backward-facing step is formulated to demonstrate that not only can DNNs discover underlying correlation behind fluid data but also they can be implemented in RANS to predict flow characteristics without numerical stability issues. Finally, the presented research is a crucial stepping-stone toward the data-driven turbulence modeling, which potentially benefits the design of data-driven experiments that can be used to validate fluid models with ML-based fluid closures.

42 ENGINEERING↗

DOE OSTI2021

Gravitomagnetic vorticity generation in black hole accretion discs: a potential spatial constraint on plasma flow stability

We calculate the vorticity generation rate in the accretion disc near a slowly rotating black hole in the low-velocity, weak-field limit of general relativity. Specifically, we find that the frame-dragging effect due to the black hole’s rotation – manifested through the gravitomagnetic field – can generate vorticity in a moving plasma in the accretion disc. The mechanism remains operational as long as the accretion disc has non-negligible vertical height and is independent of the exact thermodynamical profile of the disc. The enstrophy density generation rate, as a measure of turbulence and dissipation, is presented, which indicates that the frame-dragging effect can disrupt the stability of the disc away from the z = 0 plane.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Estimating and Evaluating Roughness Length and Displacement Height in Heterogeneous Urban Environments

The roughness length (z 0 ) and displacement height (z d ) are essential surface-layer parameters in numerical models (e.g., weather, climate, wall-modeled LES, etc.). This work evaluates the consistency of z 0 and z d estimates from morphometric and anemometric methods using data from two eddy-covariance flux towers (AmeriFlux US-INg and US-INc) in Indianapolis, IN. Results show inconsistencies in estimated z 0 and z d values depending on the chosen method. The two evaluated anemometric methods estimate non-physical values of z d when compared to roughness elements surrounding both towers. Additionally, predictions of mean wind speed using surface-layer similarity theory with morphometric estimates exhibit a bias during near-neutral and stable conditions relative to observations. The overestimation of mean wind speed by surface layer similarity theory is consistent with previous observational and modeling studies in urban areas, suggesting that the application of similarity theories to urban environments may have limitations. Differentiation of vegetation from built structures appears to impact morphometric z 0 and z d estimates, particularly where vegetation is abundant; however, it has little impact on correcting biases in the similarity theory. Specifically, we find that existing similarity theories using morphometric estimates underestimate integral velocity and length scales, and the degree of underestimation depends on the stability conditions. Accounting for the degree of anisotropy in surface-layer turbulence helps reduce the biases between similarity theories and observations during unstable conditions, but not in near-neutral cases. Future work is needed to identify the cause of such biases for near-neutral conditions.

Aerodynamic roughness length↗

DOE OSTI2025

Exploring the Effects of Varying Pre-Chamber Geometry in a Heavy-Duty Natural Gas Optical Engine under Dilution Conditions

Pre-chamber combustion is an advanced ignition strategy that has been shown to enhance spark ignition (SI) combustion stability in natural gas (NG) engines by providing distributed ignition sites from turbulent jets and enhancing main-chamber turbulence. Pre-chamber combustion has been proven especially advantageous compared to SI in ultra-lean and dilute operating conditions. This work involves experimental investigation of the effects of varying passive pre-chamber nozzle configuration on pre-chamber and main chamber combustion under simulated exhaust gas recirculation (EGR) dilution (0 and 20%) conditions in a heavy-duty, single-cylinder, optically accessible NG engine at stoichiometric fuel-air ratio. Pre-chamber nozzle configurations include four pre-chambers with constant nozzle area to pre-chamber volume ratio (A/V) with different nozzle sizes and orientations and one configuration with larger nozzles. The optical engine is operated in a skip-fire sequence consisting of 18 motored cycles followed by two consecutive fired cycles to elucidate the effect of combustion residuals (internal EGR) on combustion evolution. Pressure-based diagnostics are used to monitor pre-chamber and subsequent main chamber combustion, and optical diagnostics include high-speed OH* chemiluminescence to visualize the development of pre-chamber jets and the resulting ignition of the main chamber charge. Heat release analysis of the in-cylinder pressure data indicates that a faster pre-chamber pressure rise does not always translate into faster main-chamber combustion. The pre-chamber with the smallest nozzle diameter produced the highest pre-chamber pressure rise and fastest combustion under non-diluted conditions. However, dilution delays the main chamber ignition for smaller nozzles despite a comparable rise in pre-chamber pressure compared to configurations with larger nozzles. This effect is more pronounced for cycles with in-cylinder combustion residuals in addition to external dilution. Additionally, it was observed that pre-chambers with swirling nozzles have a faster pressure rise in the pre-chamber and main chamber under dilute conditions. Optical diagnostics suggest that the main reason for the delay between the pre-chamber pressure rise and main-chamber combustion lies in jet quenching and delayed re-ignition, which can even lead to misfire if jets emitted from small nozzles combined with dilution fail to re-ignite.

Dhotre, Akash [University Of Minnesota-Twin Cities↗

DOE OSTI2021

Stability Functions in the Stable Surface Layer Derived from the Mellor–Yamada–Nakanishi–Niino (MYNN) Scheme

It is desirable that a surface layer scheme in an atmospheric numerical model is consistent with an atmospheric boundary layer scheme incorporated in the same model. In this study, stability functions based on Monin–Obukhov similarity theory for momentum and heat, φ m and φ h , in the stable surface layer are derived from the Mellor–Yamada–Nakanishi–Niino (MYNN) scheme. The resulting stability functions are approximated by φ m = 1 + 4.8z/L and φ h = 0.74 + 6.0z/L, which can be analytically integrated with respect to height z to obtain momentum and heat fluxes, where L is the Obukhov length. The fluxes thus, obtained are compared with those acquired from stability functions in four previous studies: they are nearly the same as those from two of them, and show better agreement with observational data of the Surface Heat Budget of the Arctic Ocean experiment (SHEBA) over sea ice than those from the other two studies. Detailed comparisons of the results of the MYNN scheme with the SHEBA data suggest that significant variations of the fluxes observed during “winter” when the ice was covered with dry snow may have been caused by those of the surface roughness around the observational site. The stability functions obtained from the MYNN scheme predict that the bulk and flux Richardson numbers approach critical values of 0.26 and 0.21, respectively, in the limit of z/L → ∞. These critical values result from the Kolmogorov hypothesis applied to the turbulent dissipation in the MYNN scheme and are considered to correspond to a transition from Kolmogorov to non-Kolmogorov turbulence.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Uniform momentum and temperature zones in unstably stratified turbulent flows

Wall-bounded turbulent flows exhibit a zonal arrangement, in which streamwise velocity organizes into uniform momentum zones (UMZs), separated by thin layers of elevated interfacial shear. While significant research efforts have focused on these structural features in neutrally stratified flows, the effects of unstable thermal stratification on UMZs and on analogous uniform temperature zones (UTZs) have not been considered previously. In this article, statistical properties of UMZs and UTZs are investigated using a suite of large eddy simulations of unstably stratified turbulent channel flow spanning weakly to highly convective conditions. When normalized by the friction velocity and stability-dependent mixing length, the mean velocity gradient based on UMZ interfacial velocity jumps and the vorticity thickness exhibits good collapse for all stabilities, establishing a link between UMZ properties and scaling predictions from Monin–Obukhov similarity theory. A similar relationship is found between UTZ properties and surface-layer scaling of the mean temperature gradient. In the mixed layer, mean UMZ depth is quasi-constant with wall-normal distance, while the deepest UTZs are found in the centre of the boundary layer. These instantaneous structures are found to be linked to the well-mixed velocity and temperature profiles in the convective mixed layer. Conditional averaging indicates that both UMZ and UTZ interfaces are associated with ejections of momentum and warm updrafts below the interface and sweeps of momentum and cool downdrafts above the interface. These results demonstrate a tangible connection between instantaneous structural features, mean properties and scaling laws in unstably stratified flows.

Mechanics↗