Rayleigh–Taylor instability with gravity reversal
Not Available
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Effects of Atwood and Reynolds numbers on the evolution of buoyancy-driven homogeneous variable-density turbulence
The evolution of buoyancy-driven homogeneous variable-density turbulence (HVDT) at Atwood numbers up to 0.75 and large Reynolds numbers is studied by using high-resolution direct numerical simulations. To help understand the highly non-equilibrium nature of buoyancy-driven HVDT, the flow evolution is divided into four different regimes based on the behaviour of turbulent kinetic energy derivatives. The results show that each regime has a unique type of dependence on both Atwood and Reynolds numbers. It is found that the local statistics of the flow based on the flow composition are more sensitive to Atwood and Reynolds numbers compared to those based on the entire flow. It is also observed that, at higher Atwood numbers, different flow features reach their asymptotic Reynolds-number behaviour at different times. The energy spectrum defined based on the Favre fluctuations momentum has less large-scale contamination from viscous effects for variable-density flows with constant properties, compared to other forms used previously. The evolution of the energy spectrum highlights distinct dynamical features of the four flow regimes. Thus, the slope of the energy spectrum at intermediate to large scales evolves from -7/3 to -1 , as a function of the production-to-dissipation ratio. The classical Kolmogorov spectrum emerges at intermediate to high scales at the highest Reynolds numbers examined, after the turbulence starts to decay. Finally, the similarities and differences between buoyancy-driven HVDT and the more conventional stationary turbulence are discussed and new strategies and tools for analysis are proposed.
42 ENGINEERING↗
Impact of Biofuel Blends on Black Carbon Emissions from a Gas Turbine Engine
With a view towards renewable feedstock while decreasing the CO2 footprint associated with jet travel, biofuels are receiving increasing interest. Presented here is an overview of non-volatile particulate matter (nvPM) emissions, i.e. “soot” as gauged by varied TEM analyses for a J-85 turbojet fueled with Jet-A as well as with blends of Jet-A using Camelina biofuel. A unifying explanation is provided encompassing fuels and spanning length scales across three orders of magnitude. The variation of primary particle size, aggregate size and nanostructure are analyzed as a function of biofuel blend across range of engine thrust levels. The postulate is based on where fuels start along the soot formation pathway. Increasing biofuel content lowers aromatic concentration while placing increasing dependence upon fuel pyrolysis reactions to form the requisite concentration of aromatics for particle inception and growth. The required “kinetic” time for pyrolysis reactions to produce benzene and multi-ring PAHs allows increased fuel-air mixing by turbulence, diluting the fuel-rich soot-forming regions, effectively lowering their equivalence ratio. With a lower precursor concentration, particle inception is slowed, the resulting concentration of primary particles is lowered and smaller aggregates. The lower equivalence ratio also results in smaller primary particles given the lower concentration of growth species.
Kumal, Raju R.↗
Demonstration of an Automated System for Vertical Profiles of Volatile Organic Compounds
Volatile organic compounds (VOCs) play important roles throughout the atmosphere, many of which are altitude dependent. This highlights the need for easily deployable devices to sample VOCs across different atmospheric layers. To address this, we present the design and initial application of a Time Resolved Automated Volatile organIc compounds Sampling system (TRAVIS). VOCs are collected on sorbent tubes, which are subsequently analyzed by a thermal desorption gas chromatography mass spectrometry pipeline. TRAVIS leverages a piezoelectric pump with an integrated pressure sensor for precise (0.1% flow rate relative standard deviation) and accurate (−3 ± 2% error in VOC quantitation) measurements. Via deployment on a tethered balloon system over an agricultural area, TRAVIS is used to show consistent vertically resolved VOC concentrations in a well-mixed (i.e., turbulent) atmosphere (e.g., 5% relative standard deviation for isoprene) and vertically dependent concentrations for a stratified atmosphere (e.g., prior to boundary layer development). Furthermore, we also show VOC information from an intermittent plume via both targeted and untargeted analysis, highlighting future applications for spurious events in agriculture, air quality monitoring, and environmental impact. Overall, the development of TRAVIS represents a lightweight, accurate, sensitive, and precise VOC sampling module for the scientific community.
Aerosols↗
Improving a Multilevel Turbulence Closure Model for a Shallow Lake in Comparison With Other 1-D Models
Lakes differ from lands in water availability, heat capacity, albedo, and roughness, which affect local surface-atmospheric interactions. This study modified a multilevel upper ocean model (UOM) for lake applications and evaluated its performance in Lake Taihu (China) with comprehensive measurements against three popular one-dimensional (1-D) lake models. These models were based on different concepts, including the self-similarity (FLake), the wind-driven eddy diffusion (LISSS), the k-ε turbulence closure (SIMSTRAT), and a simplified turbulence closure (UOM). The surface flux scheme in these models was unified to exclude the discrepancies in representing air-lake exchanges. All models in their default formulations presented obvious cold water temperature biases and largely underestimated the lake surface temperature (LST) diurnal range. For each model, these deficiencies were significantly reduced by incorporating new physics schemes or calibrated tunable parameters based on systematic sensitivity tests. The primary modifications for UOM included (1) a new scheme of decreased surface roughness lengths to better characterize the shallow lake, (2) a solar radiation penetration scheme with increased light extinction coefficient and surface absorption fraction to account for the high water turbidity, and (3) turbulent Prandtl number increased by a factor of 20 to reduce the turbulent vertical mixing. All other models were improved in these three aspects (roughness, extinction, and mixing) within their original formulations. Given these improvements, UOM showed superior performance to other models in capturing LST diurnal cycle and daily to seasonal variations, as well as summer-autumn vertical stratification changes. The new UOM is well suited for application in shallow lakes.
54 ENVIRONMENTAL SCIENCES↗
Kelvin–Helmholtz instability under stabilizing parallel magnetic field in nonhomogeneous compressible MHD flows
We study the Kelvin–Helmholtz instability (KHI) for the general case of a compressible, nonhomogeneous, magnetized plasma flow. The study is limited to a vortex sheet interface with an imposed parallel magnetic field. We introduce a new formalism based on a convective Mach number M c , a convective Alfvénic Mach number M Ac , and a total convective Mach number that combines the two. We derive an analytic expression of the KHI growth rate for a homogeneous flow (i.e., zero Atwood number, A=0) that converges toward both the expression for unmagnetized compressible flow and Chandrasekhar's expression for magnetized incompressible flow. Otherwise, the dispersion relation is solved numerically and allows deriving general stability diagrams of magnetized KHI for the triplet (A, M c , β −plasma) parameters. We show these parameters uniquely define all configurations for a parallel magnetic field. We also construct diagrams with respect to the convective Alfvénic Mach number, the β − plasma parameter, or the magnetic field showing which magnetic field strength is required for stabilizing a given shear flow. The theoretical growth rates are compared with 18 simulations made with the GAMERA code, currently used for 3D magnetospheric simulations. Finally, we apply our results to the analysis of a past KHI experiment performed at the OMEGA laser facility, showing linear theory succeeds to provide accurate estimates of the growth rate at early times. We further discuss how our results can inform future experiments in the high-Mach magnetized regime at the National Ignition Facility. Possible limitations of the study due to resistive, mixing, or turbulence effects are discussed.
compressible flows↗
Finding evidence for inflation and the origin of galactic magnetic fields with CMB surveys
The origin of the μG magnetic fields observed in galaxies is unknown. One promising scenario is that magnetic fields generated during inflation, larger than 0.1 nG on Mpc scales, were adiabatically compressed to μG strengths in galaxies during structure formation. Thus, detecting a scale-invariant primordial magnetic field (PMF) above 0.1 nG on Mpc scales just after recombination would indicate an inflationary origin of galactic magnetic fields. Furthermore, this would also provide compelling evidence that inflation occurred since only an inflationary mechanism could generate such a strong, scale-invariant magnetic field on Mpc scales. In contrast, constraining the scale-invariant PMF strength to be below 0.1 nG would imply an inflationary scenario is not the primary origin, since such weak PMFs cannot be amplified enough via adiabatic compression to produce the strength of the galactic fields we observe today. We find that measurements of anisotropic birefringence by future CMB surveys will be able to improve the sensitivity to Mpc-scale inflationary PMFs by an order of magnitude, and, in particular, that CMB-HD would lower the upper bound to 0.072 nG at the 95% C.L., which is below the critical 0.1 nG threshold for ruling out a purely inflationary origin. If inflationary PMFs exist, we find that a CMB-HD survey would be able to detect them with about 3σ significance or higher, providing evidence for inflation itself.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Lagrangian analysis for turbulent transport in variable-density turbulence
In this report Lagrangian analysis of materially conserved scalars is applied to the problem of turbulent transport in variable-density flows. The consequences of an additional material conserved quantity, the density, is generally not acknowledged and leads to significant and meaningfully different expressions for turbulent transport in the moment equations. The formal Lagrangian analysis produces gradient transport expressions substantially different from those obtained by the physically intuitive “argument by analogy” method used in computational models. Various intuitive arguments, in Favre and Reynolds averaged settings, are contrasted to the formal Lagrangian results. Using expressions from the formal analysis, we derive consistent gradient transport closures for the turbulent transport terms that appear in the first- and second-order Favre moment equations. Results for coupled multispecies turbulent transport are given. The analysis is limited to variable-density turbulence in which the dilatation of the fluctuating velocity is small. The results are applicable to turbulent combustion and to stellar convection problems in which the density fluctuations are on the order of the mean density.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Mass and momentum transport in the Tilted Rocket Rig experiment
Here, high-fidelity large eddy simulations of the inclined Rayleigh-Taylor instability in the Tilted Rocket Rig experimental configuration [Smeeton and Youngs, AWE Report No. O 35/87 (1987)] are performed using a tenth-order compact finite difference code. These simulations are analyzed for spatial distributions of turbulent kinetic energy, turbulent mass flux velocity, species mass fraction flux, species mass fraction variance, and Favre-averaged Reynolds stresses at two time instances, t = 37 ms and t = 52 ms. Additionally, the vertical distribution of the components of the unclosed budget equations over the center of the domain for these quantities are examined. The dominant terms of these budget equations are further decomposed to examine the principal contributions to these terms in each axis. Notably, the principal contribution to the horizontal turbulent mass flux velocity budget is found to be from a term which is commonly neglected in many Reynolds-averaged Navier-Stokes models.
42 ENGINEERING↗
Scalar Power Spectra and Scalar Structure Function Evolution in the Richtmyer–Meshkov Instability Upon Reshock
The Richtmyer–Meshkov instability of a twice-shocked gas interface is studied using high-speed planar laser-induced fluorescence in the Wisconsin Shock Tube Laboratory's vertical shock tube. The initial condition is a shear layer with broadband diffuse perturbations at the interface between a helium–acetone mixture and argon. Furthermore, this initial condition is accelerated by a shock of nominal strength M = 1.9, and then accelerated again by the transmitted shock that reflects off the end wall of the tube. Three individual experiments are analyzed, the energy spectrum and the structure functions of the light gas mole fraction field are calculated and compared.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Droplet Growth or Evaporation Does Not Buffer the Variability in Supersaturation in Clean Clouds
Abstract Water vapor supersaturation in clouds is a random variable that drives activation and growth of cloud droplets. The Pi Convection–Cloud Chamber generates a turbulent cloud with a microphysical steady state that can be varied from clean to polluted by adjusting the aerosol injection rate. The supersaturation distribution and its moments, e.g., mean and variance, are investigated for varying cloud microphysical conditions. High-speed and collocated Eulerian measurements of temperature and water vapor concentration are combined to obtain the temporally resolved supersaturation distribution. This allows quantification of the contributions of variances and covariances between water vapor and temperature. Results are consistent with expectations for a convection chamber, with strong correlation between water vapor and temperature; departures from ideal behavior can be explained as resulting from dry regions on the warm boundary, analogous to entrainment. The saturation ratio distribution is measured under conditions that show monotonic increase of liquid water content and decrease of mean droplet diameter with increasing aerosol injection rate. The change in liquid water content is proportional to the change in water vapor concentration between no-cloud and cloudy conditions. Variability in the supersaturation remains even after cloud droplets are formed, and no significant buffering is observed. Results are interpreted in terms of a cloud microphysical Damköhler number (Da), under conditions corresponding to , i.e., the slow-microphysics regime. This implies that clouds with very clean regions, such that is satisfied, will experience supersaturation fluctuations without them being buffered by cloud droplet growth. Significance Statement The saturation ratio (humidity) in clouds controls the growth rate and formation of cloud droplets. When air in a turbulent cloud mixes, the humidity varies in space and time throughout the cloud. This is important because it means cloud droplets experience different growth histories, thereby resulting in broader size distributions. It is often assumed that growth and evaporation of cloud droplets buffers out some of the humidity variations. Measuring these variations has been difficult, especially in the field. The purpose of this study is to measure the saturation ratio distribution in clouds with a range of conditions. We measure the in-cloud saturation ratio using a convection cloud chamber with clean to polluted cloud properties. We found in clouds with low concentrations of droplets that the variations in the saturation ratio are not suppressed.
54 ENVIRONMENTAL SCIENCES↗
Coupling coarse-mesh CFD with fine-mesh CFD for modeling molten-salt reactors in the Virtual Test Bed (VTB)
The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims at developing a simulation tool kit to accelerate the development and deployment of nuclear power technologies. NEAMS multiphysics tools have been designed to provide numerical simulation support for the design and licensing of GEN IV reactors. Pronghorn is NEAMS's coarse-mesh computational fluid dynamics (CFD) tool, which is designed to run 3D core transients in GEN IV reactors at a reduced computational cost. To increase their accuracy, coarse-mesh CFD simulations require calibrated closure coefficients. One way of computing these coefficients is via the Nek5000, NEAMS's high-fidelity CFD tool. This article discusses our current research lines in informing Pronghorn closure coefficients via Nek5000 to enable multiphysics simulations of the core cavity of the molten-salt fast reactor. We present an application in which Pronghorn mixing length turbulent viscosity has been calibrated from Nek5000 simulations. The resulting Pronghorn thermal-hydraulics model is then coupled to Griffin, the NEAMS neutron transport solver, to solve for the steady-state configuration of the molten-salt fast reactor. (authors)
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Experimental Validation of Dense Plasma Transport Models using the Z-Machine
Mixing of cold, higher-Z elements into the fuel region of an inertial confinement fusion target spoils the fusion burn efficiency. This mixing process is driven by both "turbulent" and "atomic" mixing processes, the latter being modeled through transport corrections to the basic hydrodynamic models. Recently, there has been a surge in the development of dense plasma transport modeling and the associated transport coefficients; however, experimental validation remains in its infancy.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Advancement of hybrid fluid-kinetic modeling for HEDP and ICF science
We report on the development progress of a hybrid fluid-kinetic code for simulating fluids and plasmas in a wide range of environments, such as laser–matter interactions, inertial confinement fusion, magnetic confinement fusion, and pulsed power. The suite of numerical tools under development utilizes heterogeneous computer architectures and leverages the benefits of particle–based simulation techniques. By working to combine the kinetic particle-in-cell (PIC) model with a particle-based fluid simulation technique, such as smoothed particle hydrodynamics, we are developing a flexible framework capable of accurately modeling complex flows within and between kinetic and fluid regimes. The TriForce code is under development as a C++ framework for parallel, 3D, particle-based, hybrid fluid-kinetic plasma simulations. The fluid half of TriForce will be based upon the meshless smoothed-particle-hydrodynamics (SPH) approach, well-suited for shear, mixing, and turbulence, whereas the kinetic half resembles a traditional particle-in-cell (PIC) code; other particle-based approaches to fluid modeling that do use a mesh are also possible to use and are under investigation. Maxwell’s electromagnetic field equations are solved either via explicit or implicit algorithms or approximated via resistive magnetohydrodynamics (MHD) using an Ohm’s law and resulting induction equation (extended MHD is under development). A primary goal of enabling direct comparisons, from the same code, between results from the variants of MHD and implicit electromagnetic solutions is to improve our fundamental understanding of systems with magnetic fields. The code is under development to recover results from both radiation-MHD and fully kinetic codes in those limits, and is continuing to be developed from other follow-on grants to operate in between where both descriptions may co-exist and interact. For certain applications, it is desired for a simulation to contain fluid ions and electrons as well as kinetic ions and electrons. Typically, it is too computationally intensive to model a full-scale ICF or HEDP experiment fully kinetically since many cycles are expended with very small time steps on modeling the fluid part of a material that is well treated by the fluid approximation. In this case, many traditional PIC particles can be replaced with a single fluid particle representing the thermal part of the distribution function, and there are fewer needed kinetic particles, which describe the non-thermal part and can be sub-cycled relative to the fluid particle advance. Furthermore, a pure fluid code may, depending on the problem, simply lack many physically important details that are beyond the scope of its reduced approximations and assumptions. In this report, we summarize the objectives achieved in the development of the collisional and kinetic half of the code, and the physics problems to which the code has been applied in the areas of advanced and innovative fusion concepts, pulsed power, and magneto-inertial fusion.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Chaos, Percolation and the Coronavirus Spread
The dynamics of the spreading of the COVID-19 virus has similar features to the turbulent flow, chaotic maps and other non-linear systems: a small seed grows exponentially and eventually saturates. Like in the percolation model, the virus is mostly dangerous if the probability of transmission (or the bond probability p in the percolation model) is high. This suggests a relation with the population density, ρ s , which must be higher than a certain value (ρ s > 1000 persons/km 2 ). A ‘seed’ implanted in such towns grows vigorously and affects nearby places at distance x. Thus, the spreading is governed by the ratio ρ = ρ s /x. Assuming a power law dependence τ of the number of positives to the virus N + from ρ, we find τ = 0.55, 0.75 and 0.96 for S.Korea, Italy and China respectively.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Atmospheric effects on sonic-boom signatures
Atmospheric refraction, grazing incidence wave impingement, and turbulence interaction effects on sonic boom ground pressure signatures
TURBULENT MIXING↗
Spectrum of stationary homogeneous magnetohydrodynamic turbulence.
Spectral distribution of kinetic and magnetic energies from stationary homogeneous MHD TURBULENCE
MAGNETOHYDRODYNAMIC TURBULENCE↗
Thermal and chemi-ionization processes in afterburning rocket exhausts.
Mathematical model for thermal and chemi- ionization processes in turbulent nonequilibrium afterburning rocket exhausts plume, studying neutral and charged species distribution
DENSITY DISTRIBUTION↗