SEARCH · Engineering Papers

Results for “Turbulent Mixing”

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 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2024

How turbulence spreading improves power handling in quiescent high confinement fusion plasmas

Viable magnetic fusion devices necessitate combining good confinement with effective power flux handling. A major concern for ITER, and devices beyond, is the divertor heat load width, which sets peak boundary heat loads on the plasma-facing materials. Current estimates of the heat flux width are narrow for future reactors. Here, we demonstrate how pedestal turbulence can expand into, or entrain, the stable scrape-off-layer and so broaden the heat flux width beyond these neoclassical predictions. Employing combined theoretical, computational, and experimental approaches, we focus on quiescent high confinement discharges on the DIII-D tokamak, but the results are of broader significance. Our findings uncover common trends in the edge turbulence intensity flux, the pressure perturbation skewness, and the turbulence mixing length, which together determine the heat flux width. This research demonstrates the physics of scrape-off-layer broadening by turbulence and highlights the promise of a turbulent pedestal for successful core-edge integration in ITER and future fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Analytical Investigation of the Role of Lateral Mixing in the Evolution of Nonprecipitating Cumulus. Part II: Dissolving Stage

A minimalistic analytical model allowing analysis of the dissolving stage of nonprecipitating convective clouds is proposed. The model takes into account two mechanisms: turbulent mixing with a dry environment and cloud volume settling. The temporal changes in the spatial structure of a cloud and in its immediate environment in the course of cloud dissolving are analyzed. The comparison of the effects of a temperature increase in the course of cloud descent and mixing with dry surrounding air shows that the descent is a dominating factor determining a decrease in the liquid water content (LWC), while mixing has a stronger effect on the cloud shape. Narrowing/broadening of clouds due to lateral mixing with dry air during cloud dissolving is determined by the potential evaporation parameter proportional to the ratio of the saturation deficit in the cloud environment to LWC. An equation for cloud dissolving time is obtained. After a cloud totally dissolves, it leaves behind an area with humidity exceeding that of the environment. The main parameter determining the dissolving time is the downdraft velocity. It should exceed 50 cm s -1 in order to provide reasonable dissolving time. The turbulent intensity, LWC, and humidity of the environment air also have an impact on dissolving time: the lower the LWC and the humidity of environment air, the faster cloud dissolving is. The simple solution presented in this paper can be useful for evaluation of cloud characteristics at the dissolving stage and can be included in procedures of parameterization of cloud cover formed by nonprecipitating or slightly precipitating cumulus clouds (Cu). Values of the environment humidity and temperature, LWC at cloud top, cloud width, vertical velocity of downdraft, and the turbulent coefficient should be parameters of this parameterization.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Improved Representation of Horizontal Variability and Turbulence in Mesoscale Simulations of an Extended Cold-Air Pool Event

Abstract Cold-air pools (CAPs), or stable atmospheric boundary layers that form within topographic basins, are associated with poor air quality, hazardous weather, and low wind energy output. Accurate prediction of CAP dynamics presents a challenge for mesoscale forecast models in part because CAPs occur in regions of complex terrain, where traditional turbulence parameterizations may not be appropriate. This study examines the effects of the planetary boundary layer (PBL) scheme and horizontal diffusion treatment on CAP prediction in the Weather Research and Forecasting (WRF) Model. Model runs with a one-dimensional (1D) PBL scheme and Smagorinsky-like horizontal diffusion are compared with runs that use a new three-dimensional (3D) PBL scheme to calculate turbulent fluxes. Simulations are completed in a nested configuration with 3-km/750-m horizontal grid spacing over a 10-day case study in the Columbia River basin, and results are compared with observations from the Second Wind Forecast Improvement Project. Using event-averaged error metrics, potential temperature and wind speed errors are shown to decrease both with increased horizontal grid resolution and with improved treatment of horizontal diffusion over steep terrain. The 3D PBL scheme further reduces errors relative to a standard 1D PBL approach. Error reduction is accentuated during CAP erosion, when turbulent mixing plays a more dominant role in the dynamics. Last, the 3D PBL scheme is shown to reduce near-surface overestimates of turbulence kinetic energy during the CAP event. The sensitivity of turbulence predictions to the master length-scale formulation in the 3D PBL parameterization is also explored. Significance Statement In this article, we demonstrate how a new framework for modeling atmospheric turbulence improves cold pool predictions, using a case study from January 2017 in the Columbia River basin (U.S. Pacific Northwest). Cold pools are regions of cold, stagnant air that form within valleys or basins, and improved forecasts could help to mitigate the risks they pose to air quality, transportation, and wind energy production. For the chosen case study, our tests show a reduction in temperature and wind speed errors by up to a factor of 2–3 relative to standard model options. These results strongly motivate continued development of the framework as well as its application to other complex weather events.

17 WIND ENERGY↗

DOE OSTI2026

MRV Challenge 3: velocity and passive scalar comparison in a 3D turbulent flow

The third iteration of a magnetic resonance velocimetry (MRV) challenge comparison activity consisting of 3D flow measurements for a turbulent flow in a water channel moving past centrally positioned cubic obstacles is reported. In this challenge iteration, MRV measurements are coupled with either temperature or concentration measurements to extend the diagnostic utility. Five research teams from around the world conducted the measurements on a single shared apparatus. The water channel included partial elements along the channel sidewalls that varied in height and precluded easy optical accessibility. In addition to a turbulent mainstream flow, a secondary flow entered the channel between the second and third cubic elements from a square injector hole at the channel bottom wall. The injector jet interacts with the mainstream flow and mixes turbulently as it advects downstream. For the selected flow regime and water solutions used as working fluids, the mixing of the higher temperature or concentration through the secondary flow with the mainstream flow satisfies the same dynamics so that both temperature and concentration fields behave as passive scalars and can be directly compared. The measurements are explained in detail by each participating team, and the results are interpolated onto a common grid and compared using line profiles, contour plots, and isosurfaces.

Benson, Mike [ORNL] (ORCID:000000023210116X)↗

DOE OSTI2023

Enhancements in Cloud Condensation Nuclei Activity From Turbulent Fluctuations in Supersaturation

Abstract The effect of aerosols on the properties of clouds is a large source of uncertainty in predictions of weather and climate. These aerosol‐cloud interactions depend critically on the ability of aerosol particles to form cloud droplets. A challenge in modeling aerosol‐cloud interactions is the representation of interactions between turbulence and cloud microphysics. Turbulent mixing leads to small‐scale fluctuations in water vapor and temperature that are unresolved in large‐scale atmospheric models. To quantify the impact of turbulent fluctuations on cloud condensation nuclei (CCN) activation, we used a high‐resolution Large Eddy Simulation of a convective cloud chamber to drive particle‐based cloud microphysics simulations. We show small‐scale fluctuations strongly impact CCN activity. Once activated, the relatively long timescales of evaporation compared to fluctuations causes droplets to persist in subsaturated regions, which further increases droplet concentrations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Impact of Turbulence Representation on the Relationship Between Cloud Feedback and Aerosol‐Cloud Interaction in an E3SMv2 Perturbed Parameter Ensemble

Recent studies reveal an anti‐correlation between global cloud feedback (CF) and effective radiative forcing due to aerosol‐cloud interaction (ERFaci) in Earth system models, but the physical mechanisms underlying it remain uncertain. Here we investigate how different turbulence representations contribute to this relationship over the global ocean using an ensemble of Energy Exascale Earth System Model version 2 simulations with perturbed turbulence parameters. The anti‐correlation appears only in the tropical ascent regime. In the Northern Hemisphere midlatitude and high latitude regimes, there is no significant correlation, and in the tropical marine low cloud and Southern Ocean regimes, the correlation is positive. These opposite correlations are primarily driven by opposing CF responses to perturbed parameters. We find that the mean‐state turbulent mixing strength affects both CF and ERFaci, enabling strong correlations in certain regimes. This study highlights the complex linkages between CF and ERFaci through turbulent processes across diverse cloud regimes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Upper ocean mixing, surface heat fluxes, and heat content variability in the upper 150 m during Hurricane Laura (2020)

Observations (from thermistor sensors and a hurricane glider) and model-based temperature were used to answer questions including the following: 1) What physical and environmental conditions explain the rapid intensification of Hurricane Laura? 2) How did the pre-existing warm mixed layer at Stone mooring (StM) influence the degree of cooling in the mixed layer? 3) How did vertical mixing, surface heat fluxes, and advective processes collectively shape the mixed layer heat evolution at StM? Observations measured by thermistor sensors showed that Hurricane Laura induced a cooling of -1.2oC in the mixed layer when Hurricane Laura’s wind speed increased to 47 m/s on 26 August. This contrasts with the -1.04oC mixed layer temperature change estimated from model-based temperature. At StM, the presence of a 31oC warm mixed layer and elevated heat content (60-80kJ/cm2) effectively preconditioned the upper ocean ahead of Hurricane Laura’s passage. After Hurricane Laura passed, turbulent mixing associated with the storm transported the heat anomalies downward into subsurface layers. In addition, this study shows that a surface heat flux of 5.04 kJ/cm2 supported the intensity of Hurricane Laura as it traversed over the StM. Using a mixed layer heat balance model, this study shows that entrainment flux and surface flux mostly contributed to the observed change in the mixed layer temperature.

Tsai, Senam↗

DOE OSTI2021

Rayleigh–Taylor and Richtmyer–Meshkov instabilities: A journey through scales

Hydrodynamic instabilities such as Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities usually appear in conjunction with the Kelvin–Helmholtz (KH) instability and are found in many natural phenomena and engineering applications. They frequently result in turbulent mixing, which has a major impact on the overall flow development and other effective material properties. This can either be a desired outcome, an unwelcome side effect, or just an unavoidable consequence, but must in all cases be characterized in any model. The RT instability occurs at an interface between different fluids, when the light fluid is accelerated into the heavy. The RM instability may be considered a special case of the RT instability, when the acceleration provided is impulsive in nature such as that resulting from a shock wave. Here in this pedagogical review, we provide an extensive survey of the applications and examples where such instabilities play a central role. First, fundamental aspects of the instabilities are reviewed including the underlying flow physics at different stages of development, followed by an overview of analytical models describing the linear, nonlinear and fully turbulent stages. RT and RM instabilities pose special challenges to numerical modeling, due to the requirement that the sharp interface separating the fluids be captured with fidelity. These challenges are discussed at length here, followed by a summary of the significant progress in recent years in addressing them. Examples of the pivotal roles played by the instabilities in applications are given in the context of solar prominences, ionospheric flows in space, supernovae, inertial fusion and pulsed-power experiments, pulsed detonation engines and Scramjets. Progress in our understanding of special cases of RT/RM instabilities is reviewed, including the effects of material strength, chemical reactions, magnetic fields, as well as the roles the instabilities play in ejecta formation and transport, and explosively expanding flows. The article is addressed to a broad audience, but with particular attention to graduate students and researchers who are interested in the state-of-the-art in our understanding of the instabilities and the unique issues they present in the applications in which they are prominent.

42 ENGINEERING↗

DOE OSTI2020

Turbulence/flame/wall interactions in non-premixed inclined slot-jet flames impinging at a wall using direct numerical simulation

In this paper, three-dimensional turbulent non-premixed oblique slot-jet flames impinging at a wall were investigated using direct numerical simulation (DNS). Two cases are considered with the Damköhler number (Da) of case A being twice that of case B. A 17 species and 73-step mechanism for methane combustion was employed in the simulations. It was found that flame extinction in case B is more prominent compared to case A. Reignition in the lower branch of combustion for case A occurs when the scalar dissipation rate relaxes, while no reignition occurs in the lower branch for case B due to excessive scalar dissipation rate. A method was proposed to identify the flame quenching edges of turbulent non-premixed flames in wall-bounded flows based on the intersections of mixture fraction and OH mass fraction iso-surfaces. The flame/wall interactions were examined in terms of the quenching distance and the wall heat flux along the quenching edges. There is essentially no flame/wall interaction in case B due to the extinction caused by excessive turbulent mixing. In contrast, significant interactions between flames and the wall are observed in case A. The quenching distance is found to be negatively correlated with wall heat flux as previously reported in turbulent premixed flames. The influence of chemical reactions and wall on flow topologies was identified. The FS/U and FC/U topologies are found near flame edges, and the NNN/U topology appears when reignition occurs. The vortex-dominant topologies, FC/U and FS/S, play an increasingly important role as the jet turbulence develops.

42 ENGINEERING↗

DOE OSTI2024

Direct numerical simulation of a high-pressure hydrogen micromix combustor: Flame structure and stabilisation mechanism

A high-pressure hydrogen micromix combustor has been investigated using direct numerical simulation with detailed chemistry to examine the flame structure and stabilisation mechanism. The configuration of the combustor was based on the design by Schefer et al., using numerical periodicity to mimic a large square array. A precursor simulation of an opposed jet-in-crossflow was first conducted to generate appropriate partially-premixed inflow boundary conditions for the subsequent reacting simulation. The resulting flame can be described as an predominantly-lean inhomogeneously-premixed lifted jet flame. Five main zones were identified: a jet mixing region, a core flame, a peripheral flame, a recirculation zone, and combustion products. The core flame, situated over the jet mixing region, was found to burn as a thin reaction front, responsible for over 85% of the total fuel consumption. The peripheral flame shrouded the core flame, had low mean flow with high turbulence, and burned at very lean conditions (in the distributed burning regime). It was shown that turbulent premixed flame propagation was an order-of-magnitude too slow to stabilise the flame at these conditions. Stabilisation was identified to be due to ignition events resulting from turbulent mixing of fuel from the jet into mean recirculation of very lean hot products. Ignition events were found to correlate with shear-driven Kelvin-Helmholtz vortices, and increased in likelihood with streamwise distance. At the flame base, isolated events were observed, which developed into rapidly burning flame kernels that were blown downstream. Further downstream, near-simultaneous spatially-distributed ignition events were observed, which appeared more like ignition sheets. The paper concludes with a broader discussion that considers generalising from the conditions considered here.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Instability and turbulent relaxation in a stochastic magnetic field

An analysis of instability dynamics in a stochastic magnetic field is presented for the tractable case of the resistive interchange. Externally prescribed static magnetic perturbations convert the eigenmode problem to a stochastic differential equation, which is solved by the method of averaging. The dynamics are rendered multi-scale, due to the size disparity between the test mode and magnetic perturbations. Maintaining quasi-neutrality at all orders requires that small-scale convective cell turbulence be driven by disparate scale interaction. Here, the cells in turn produce turbulent mixing of vorticity and pressure, which is calculated by fluctuation-dissipation type analyses, and are relevant to pump-out phenomena. The development of correlation between the ambient magnetic perturbations and the cells is demonstrated, showing that turbulence will 'lock on' to ambient stochasticity. Magnetic perturbations are shown to produce a magnetic braking effect on vorticity generation at large scale. Detailed testable predictions are presented. The relations of these findings to the results of available simulations and recent experiments are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Simulation and flow physics of a shocked and reshocked high-energy-density mixing layer

This work describes a computational investigation of multimode instability growth and multimaterial mixing induced by multiple shock waves in a high-energy-density (HED) environment, where pressures exceed 1 Mbar. The simulations are based on a series of experiments performed at the National Ignition Facility (NIF) and designed as an HED analogue of non-HED shock-tube studies of the Richtmyer–Meshkov instability and turbulent mixing. A three-dimensional computational modelling framework is presented. It treats many complications absent from canonical non-HED shock-tube flows, including distinct ion and free-electron internal energies, non-ideal equations of state, radiation transport and plasma-state mass diffusivities, viscosities and thermal conductivities. The simulations are tuned to the available NIF data, and traditional statistical quantities of turbulence are analysed. Integrated measures of turbulent kinetic energy and enstrophy both increase by over an order of magnitude due to reshock. Large contributions to enstrophy production during reshock are seen from both the baroclinic source and enstrophy–dilatation terms, highlighting the significance of fluid compressibility in the HED regime. Dimensional analysis reveals that Reynolds numbers and diffusive Péclet numbers in the HED flow are similar to those in a canonical non-HED analogue, but conductive Péclet numbers are much smaller in the HED flow due to efficient thermal conduction by free electrons. It is shown that the mechanism of electron thermal conduction significantly softens local spanwise gradients of both temperature and density, which causes a minor but non-negligible decrease in enstrophy production and small-scale mixing relative to a flow without this mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Features of mid- and high-latitude low-level clouds and their relation to strong aerosol effects in the Energy Exascale Earth System Model version 2 (E3SMv2)

The E3SMv2 model, like various other global models that include representations of aerosol–cloud interactions, uses an empirically chosen lower bound on the simulated in-cloud cloud droplet number concentration (CDNC) to help constrain the effective radiative forcing of anthropogenic aerosols, ERFaer. This study identifies where ultra-low CDNCs (i.e., concentrations lower than 10 cm−3) occur in the stratiform and shallow convective clouds simulated by E3SMv2 and which of the occurrences have the strongest impact on ERF aer . Process-level analyses are presented to reveal characteristics of the cloud droplet formation and removal processes associated with impactful ultra-low CDNCs. Simulations performed with present-day emissions show that ultra-low CDNCs are most frequently found over the mid- and high-latitude oceans in both hemispheres, while the occurrences are also frequent in polluted continental regions despite the high aerosol concentrations. Ultra-low CDNCs with the largest impacts on the simulated regional and global mean ERF aer are found in the lower troposphere in the Northern Hemisphere middle and high latitudes. These cases are typically associated with large cloud fractions, strong water vapor condensation, weak turbulence, and lack of cloud droplet nucleation from aerosol activation. Under such atmospheric conditions, boosting aerosol activation and enhancing turbulent mixing of cloud droplet number can increase the simulated CDNCs, although the magnitude of the global mean ERF aer increases undesirably. The reason for this model behavior is discussed. Overall, our study suggests that mid- and high-latitude low-level stratus occurring under weak turbulence is a cloud regime worth further investigating for the purpose of identifying and addressing the root causes of ultra-low CDNCs and strong ERF aer in E3SM.

Aerosol-Cloud Interaction↗

DOE OSTI2022

Three-dimensional simulations of reshocked inclined Richtmyer-Meshkov instability: Effects of initial perturbations

The effect of initial perturbations on the evolution of the inclined Richtmyer-Meshkov turbulent mixing layer before and after reshock initiated by a shock wave with Mach number 1.55 is investigated through three-dimensional (3D) simulations using the flash code. The 3D simulations aim to reproduce both predominantly single-mode and multimode interfaces between light and heavy gases (N 2 -CO 2 , Atwood number, A≈0.22; amplitude to wavelength ratio of 0.088) which were created in an inclined shock tube facility to analyze the effects of initial conditions on mixing development in the entire flow field. The two-dimensional center slices of 3D simulations are compared with the experimental results to validate the computational code. Mixing width, mixed mass, mixed-mass thickness, and circulation in addition to concentration fields are shown to be in good agreement with the experimental data. The three-dimensional density and vorticity fields are first presented to qualitatively describe the flow behavior before and after reshock. Several measured density/velocity-related quantities indicate that the growth of the mixing material is strongly dependent on initial conditions. Before reshock and at early times after reshock, flow is clearly maintaining the memory of initial perturbations. However, at late time after reshock, although the large wavelength feature still dominates the flow motion, and the morphology of the two different interfaces indicates several differences, by breakdown of large-scale coherent structures to much finer scales, the memory of small scales of the multimode initial perturbation is not as clear as pre-reshock. Regarding three-dimensionality of the flow, before reshock in the multimode case, the baroclinic vorticity production, circulation, turbulent kinetic energy, and turbulent mass flux suggest that the small-scale roll-up features along the large inclined wavelength quickly evolves in all three dimensions. The coherent vortex tubes break down to smaller wormlike vortex structures, and turbulent fluctuations in the out-of-plane dimension are comparable to the spanwise direction. After reshock, this three-dimensionality of mixing growth was observed in the flow for both initial conditions. The results of this work represent a significant extension of previous computational studies performed on this specific topic. A different code with a different numerical method is validated through comparison with the experimental data. The initial perturbations are directly measured from the experimental results. Moreover, the entire three-dimensional experimental shock tube domain is simulated, and more quantities are investigated to understand the mixing mechanism and instability evolution in all three dimensions.

42 ENGINEERING↗

DOE OSTI2020

Temperature statistics in a radiatively heated particle-laden turbulent square duct flow

Radiation absorption by preferentially concentrated particles in a turbulent square duct flow is studied experimentally. The particle-laden flow is exposed to near-infrared radiation, and the gas phase temperature statistics are measured along the wall bisector of the duct. It is found that the instantaneous temperature fluctuations are comparable to the overall mean temperature rise. The temperature statistics at the duct centerline and near the wall are qualitatively different. The former reflects preferential concentration in isotropic flows while the latter displays evidence of particle clustering into streamwise elongated streaks. Comparison of the experimental data to a simplified heat transfer model suggests that the Lagrangian evolution of particle clusters and voids, and turbulent mixing in the vicinity of particle clusters, are important. This work was motivated by particle solar receiver technology, but the findings are also relevant to systems where there is localized heat release or mass transfer from disperse particles or droplets. It shows that obtaining Lagrangian histories of particle trajectories is an important next step towards understanding thermal transport phenomena in particle-laden turbulent flows.

42 ENGINEERING↗

DOE OSTI2020

Dynamics of Subsiding Shells in Actively Growing Clouds with Vertical Updrafts

The dynamics of a subsiding shell at the edges of actively growing shallow cumulus clouds with updrafts is analyzed using direct numerical simulation. The actively growing clouds have a fixed in-cloud buoyancy and velocity. Turbulent mixing and evaporative cooling at the cloud edges generate a subsiding shell that grows with time. A self-similar regime is observed for first- and second-order moments when normalized with respective maximum values. Internal scales derived from integral properties of the flow problem are identified. A self-similarity analysis using these scales reveals that contrary to classical self-similar flows, the turbulent kinetic energy budget terms and velocity moments scale according to the buoyancy and not with the mean velocity. The shell thickness is observed to increase linearly with time. The buoyancy scale remains time invariant and is set by the initial cloud–environment thermodynamics. The shell accelerates ballistically with a magnitude set by the saturation value of the buoyancy of the cloud–environment mixture. In this regime, the shell is buoyancy driven and independent of the in-cloud velocity. Relations are obtained for predicting the shell thickness and minimum velocities by linking the internal scales with external flow parameters. The values thus calculated are consistent with the thickness and velocities observed in typical shallow cumulus clouds. The entrainment coefficient is a function of the initial state of the cloud and the environment, and is shown to be on the same order of magnitude as fractional entrainment rates calculated for large-scale models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Coherent structures in stably stratified wall-bounded turbulent flows

To date, a growing body of literature has documented the existence and impacts of coherent structures known as large- and very-large-scale motions within wall-bounded turbulent flows under neutral and unstable thermal stratification. These coherent structures can account for a considerable fraction of the overall turbulent transport and have been found to modulate small-scale turbulent fluctuations near the wall. In the context of stably stratified flows, however, the examination of such coherent structures has garnered relatively little attention. Stable stratification limits vertical transport and turbulent mixing within flows, which makes it unclear the extent to which previous findings on coherent structures under unstable and neutral stratification are applicable to stably stratified flows. In this study, we investigate the existence and characteristics of coherent structures under stable stratification with a wide range of statistical and spectral analyses. Outer peaks in premultiplied spectrograms under weak stability indicate the presence of large-scale motions, but these peaks become weaker and eventually vanish with increasing stability. Quadrant analysis of turbulent transport efficiencies (the ratio of net fluxes to their respective downgradient components) demonstrates dependencies on both stability and height above ground, which is evidence of morphological differences in the coherent structures under increasing stability. Amplitude modulation by large-scale streamwise velocity was found to decrease with increasing gradient Richardson number, whereas modulation by large-scale vertical velocity was approximately zero across all stability ranges. For sufficiently stable stratification, large eddies are suppressed enough to limit any inner–outer scale interactions.

Greene, Brian R. (ORCID:0000000343766818)↗

DOE OSTI2023

A Scaling Analysis of Post-Detonation Mixing With Detailed Chemical Kinetics

This report documents a scaling study of hemispherical pentaerythritol tetranitrate explosives detonated in ambient air. The focus is on the mixing layer and evolution of the chemical species. The HyBurn code developed at the University of Florida is used to simulate center initiated 12, 38, 120, and 1200 mm diameter charges. The results show blast structures commonly found in the literature. Azimuthal averaging is used to radially average the flow variables to get mean values. The charges exhibit self-similar behavior at early times in a primarily gas dynamic regime, but discrepancies between the scales form at later times due to the finite rate chemistry. The turbulent mixing layer is calculated and used to define a region to conditionally average by mixture fraction. The results of the conditional averages imply a transition from a kinetic-limited regime to a mixing-limited regime with increasing scale. The scaling of the total heat release rate confirms an early time gas dynamic regime, and a later time regime where mixing and chemistry becomes impactful.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗