Engineering Papers⌕ Search

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.

At least 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effects of variable deceleration periods on Rayleigh-Taylor instability with acceleration reversals

The dynamics of an interfacial flow that is initially Rayleigh-Taylor unstable but becomes statically stable for some intermediate period due to the reversal of the externally imposed acceleration field is studied. We discuss scenarios that consider both single and double-acceleration reversals. The accel-decel (AD) case consists of a single reversal imposed at an instant after the constant acceleration instability has entered a self-similar regime. The layer of mixed fluid ceases to grow upon acceleration reversal, and the dominant mechanics are due to internal wave oscillations. Variation of mass flux and the Reynolds stress anisotropy is observed due to the action of the internal waves. Here, a second reversal of the AD case that is termed as accel-decel-accel, ADA is then explored; the response of the mixing layer is shown to depend strongly on the duration and the periodicity of the Reynolds stress anisotropy of the mixing layer during the deceleration period. We explore the effect of this variable deceleration period after the second acceleration reversal where the flow once again becomes Rayleigh-Taylor unstable based on metrics that include the integral mixing-layer width, bubble and spike amplitudes, mass flux, Reynolds stress anisotropy tensor, and the molecular mixing parameter.

42 ENGINEERING↗
DOE OSTI2020

Ducted fuel injection

Various technologies presented herein relate to heating a fuel and enhancing fuel and charge-gas mixing inside a combustion chamber to enable minimal, or no, generation of soot and/or other undesired emissions during ignition and subsequent combustion of the locally premixed mixtures. To enable sufficient mixing of the fuel and charge-gas, a jet of fuel can be directed to pass through a bore of a duct causing charge-gas to be drawn into the bore creating turbulence to mix the fuel and the drawn charge-gas. The duct can be heated to provision heating of the fuel and charge-gas mixture. The duct can be located proximate to an opening in a tip of a fuel injector. An ignition assist component can be located downstream of the duct to facilitate ignition of the fuel/charge-gas mixture.

42 ENGINEERING↗
DOE OSTI2021

Scalar mixing in a Kelvin-Helmholtz shear layer and implications for Reynolds-averaged Navier-Stokes modeling of mixing layers

Large-eddy simulation of a temporally evolving Kelvin-Helmholtz (KH) mixing layer is performed with the tenth-order compact difference code miranda to examine the steady-state behavior of a passive scalar in a shear-driven mixing layer. It is shown that the integral behavior of scalar variance in a KH mixing layer behaves similarly to the integral behavior of scalar variance in a Rayleigh-Taylor (RT) mixing layer, and mixedness of the simulated KH shear layer tends towards a value of about 0.8. It is further shown that if the k-L-a-V Reynolds-averaged Navier-Stokes (RANS) model [B. E. Morgan et al., Phys. Rev. E 98, 033111 (2018)], calibrated to reproduce steady-state mixing in an RT layer, is applied to simulate a KH mixing layer, the RANS model will significantly overpredict the magnitude of scalar variance in the KH layer. A straightforward addition to the k-L-a-V model is then suggested, and self-similarity analysis is applied to determine constraints on model coefficients. Furthermore, it is shown that with the addition of a buoyancy production term in the model equation for scalar variance, it becomes possible to eliminate the model deficiency and match steady-state mixedness in simulations of both RT and KH mixing layers with a single model calibration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

QES-Plume v1.0: a Lagrangian dispersion model

Low-cost simulations providing accurate predictions of transport of airborne material in urban areas, vegetative canopies, and complex terrain are demanding because of the small-scale heterogeneity of the features influencing the mean flow and turbulence fields. Common models used to predict turbulent transport of passive scalars are based on the Lagrangian stochastic dispersion model. The Quick Environmental Simulation (QES) tool is a low-computational-cost framework developed to provide high-resolution wind and concentration fields in a variety of complex atmospheric-boundary-layer environments. Part of the framework, QES-Plume, is a Lagrangian dispersion code that uses a time-implicit integration scheme to solve the generalized Langevin equations which require mean flow and turbulence fields. Here, QES-Plume is driven by QES-Winds, a 3D fast-response model that computes mass-consistent wind fields around buildings, vegetation, and hills using empirical parameterizations, and QES-Turb, a local-mixing-length turbulence model. In this paper, the particle dispersion model is presented and validated against analytical solutions to examine QES-Plume’s performance under idealized conditions. In particular, QES-Plume is evaluated against a classical Gaussian plume model for an elevated continuous point-source release in uniform flow, the Lagrangian scaling of dispersion in isotropic turbulence, and a non-Gaussian plume model for an elevated continuous point-source release in a power-law boundary-layer flow. In these cases, QES-Plume yields a maximum relative error below 6 % when compared with analytical solutions. In addition, the model is tested against wind-tunnel data for a uniform array of cubical buildings. QES-Plume exhibits good agreement with the experiment with 99 % of matched zeros and 59 % of the predicted concentrations falling within a factor of 2 of the experimental concentrations. Furthermore, results also emphasize the importance of using high-quality turbulence models for particle dispersion in complex environments. Finally, QES-Plume demonstrates excellent computational performance.

58 GEOSCIENCES↗
DOE OSTI2023

Scaling of Turbulence and Microphysics in a Convection–Cloud Chamber of Varying Height

Abstract The convection–cloud chamber enables measurement of aerosol and cloud microphysics, as well as their interactions, within a turbulent environment under steady‐state conditions. Increasing the size of a convection–cloud chamber, while holding the imposed temperature difference constant, leads to increased Rayleigh, Reynolds and Nusselt numbers. Large–eddy simulation coupled with a bin microphysics model allows the influence of increased velocity, time, and spatial scales on cloud microphysical properties to be explored. Simulations of a convection–cloud chamber, with fixed aspect ratio and increasing heights of H = 1, 2, 4, and (for dry conditions only) 8 m are performed. The key findings are: Velocity fluctuations scale as H 1/3 , consistent with the Deardorff expression for convective velocity, and implying that the turbulence correlation time scales as H 2/3 . Temperature and other scalar fluctuations scale as H −3/7 . Droplet size distributions from chambers of different sizes can be matched by adjusting the total aerosol injection rate as the horizontal cross‐sectional area (i.e., as H 2 for constant aspect ratio). Injection of aerosols at a point versus distributed throughout the volume makes no difference for polluted conditions, but can lead to cloud droplet size distribution broadening in clean conditions. Cloud droplet growth by collision and coalescence leads to a broader right tail of the distribution compared to condensation growth alone, and this tail increases in magnitude and extent monotonically as the increase of chamber height. These results also have implications for scaling within turbulent, cloudy mixed‐layers in the atmosphere, such as fog layers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Two self-similar Reynolds-stress transport models with anisotropic eddy viscosity

Two Reynolds-averaged Navier-Stokes models with full Reynolds-stress transport (RST) and tensor eddy viscosity are presented. These new models represent RST extensions of the $k−2L−a−\mathscr{C}$ and $k−ϕ−L−a−\mathscr{C}$ models by Morgan. Self-similarity analysis is used to derive constraints on model coefficients required to reproduce expected growth parameters for a variety of canonical flows, including Rayleigh-Taylor (RT) and Kelvin-Helmholtz (KH) mixing layers. Both models are then applied in one-dimensional simulation of RT and KH mixing layers, and the expected self-similar growth rates and anisotropy are obtained. Next, models are applied in two-dimensional simulation of the so-called “tilted rocket rig” inclined RT experiment and in simulation of a shock-accelerated localized patch of turbulence. Here it is found that RST is required to capture the qualitative growth of the shock-accelerated patch, and an anisotropic eddy viscosity provides substantial improvement over a Boussinesq treatment for the tilted rocket rig problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Capturing Secondary Kinetic Instabilities in Three‐Dimensional Dayside Reconnection Using an Improved Gradient‐Based Closure

Magnetic reconnection is a highly dynamic process that excites a wide variety of kinetic waves and instabilities. Transverse current sheet instabilities such as the lower-hybrid drift and secondary drift-kink instabilities in particular have been shown by kinetic simulations to modify the reconnection and introduce significant turbulence and mixing to the reconnection layer. Past studies using the ten-moment fluid model to capture important kinetic physics such as the electron inertia and full representation of the pressure tensor proved advantageous to a two-fluid representation of reconnection, but the model struggled when using a local relaxation closure for the heat flux to replicate the current sheet instabilities and subsequent mixing seen in kinetic simulations. This work uses the Gkeyll software framework to perform simulations of asymmetric reconnection based on the 16 October 2015 MMS crossing of a diffusion region, the Burch event. An improved gradient-based heat flux closure is implemented, showing significant improvement in secondary kinetic instabilities that grow in the current sheet. These instabilities generate turbulence which leads to growth of secondary magnetic islands and flux ropes.

Bradshaw, K. [Princeton University, NJ (United Sta↗
DOE OSTI2020

Rayleigh–Taylor Instability With Varying Periods of Zero Acceleration

We present our findings from a numerical investigation of the acceleration-driven Rayleigh–Taylor Instability, modulated by varying periods without an applied acceleration field. It is well known from studies on shock-driven Richtmyer–Meshkov instability that mixing without external forcing grows with a scaling exponent as ≈ t 0.20-0.28 When the Rayleigh–Taylor Instability is subjected to varying periods of “zero” acceleration, the structural changes to the mixing layer remain remarkably small. After the acceleration is re-applied, the mixing layer quickly resumes the profile of development it would have had if there had been no intermission. As a result, this behavior contrasts in particular with the strong sensitivity that is found to other variable acceleration profiles examined previously in the literature.

42 ENGINEERING↗
DOE OSTI2025

Evaluation of probability density function descriptions for three-component Rayleigh–Taylor mixing

Results from simulations of a three-component Rayleigh–Taylor (RT) mixing problem are presented. These simulations are conducted in heavy–light–heavy and heavy–intermediate–light configurations, and each of these configurations are further considered in high- and low-Reynolds-number regimes. This results in RT-unstable flow with one or both interfaces initially unstable, permitting the influence of problem configuration on the statistical description of three-component RT-driven mixing to be considered. Mass fraction covariances are observed to undergo a sign change through the mixing layer in all four configurations considered. This appears to be unique to the multi-component case and represents another way in which multi-component RT mixing differs from the two-component case. Qualitative and quantitative comparisons of joint and marginal probability density function (PDF) descriptions of species concentration are made. Three-, five-, and six-parameter model PDFs are compared against simulation data to assess how accurately they describe the mixing, and it is found that three-component mixing requires at least a five-parameter model PDF to accurately describe the mixing. Notably, the marginal distributions of three-component mixing do not appear to conform to a beta distribution, representing a departure from the classical two-component RT case. In conclusion, statistical neutrality also appears to influence the optimal choice of model PDF, which is found to be a function of problem configuration.

Large-eddy simulation↗
DOE OSTI2024

Transition in ICF Capsule Implosions

Longstanding design and reproducibility challenges in inertial confinement fusion (ICF) capsule implosion experiments involve recognizing the need for appropriately characterized and modeled three-dimensional initial conditions and high-fidelity simulation capabilities to predict transitional flow approaching turbulence, material mixing characteristics, and late-time quantities of interest—e.g., fusion yield. We build on previous coarse graining simulations of the indirect-drive national ignition facility (NIF) cryogenic capsule N170601 experiment-a precursor of N221205 which resulted in net energy gain. We apply effectively combined initialization aspects and multiphysics coupling in conjunction with newly available hydrodynamics simulation methods, including directional unsplit algorithms and low Mach-number correction-key advances enabling high fidelity coarse grained simulations of radiation-hydrodynamics driven transition. Our presentation includes discussion of the capsule initialization and implosion dynamics, analysis of the vorticity production budget, transition signatures, quantities of interest—late-time ion temperature and fusion-neutron yield, numerical uncertainty quantification, and comparisons with NIF data.

42 ENGINEERING↗
DOE OSTI2024

RANSBox: A zero-dimensional modular software package for Reynolds-averaged Navier-Stokes modeling

RANSBox is a zero-dimensional software package for Reynolds-averaged Navier-Stokes (RANS) modeling developed at Lawrence Livermore National Laboratory to support common implementation of RANS models across a variety of host codes with different numerical schemes and code bases. Herein this work describes the key features of RANSBox including “model-agnostic integration,” which allows new models to be implemented in RANSBox and quickly deployed to host codes without additional changes to the host code base. Three one-dimensional test problems with analytical asymptotic properties are described which can be used to verify correctness of RANSBox integration. Results for these problems are then compared across different RANS models in a single code and across different codes with a common RANS model. By using a common model implementation in RANSBox, host codes with different numerical schemes and formal orders of accuracy are demonstrated to predict the expected behavior for the three test problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Demonstrating non-isothermal vacuum membrane air dehumidification for efficient next-generation air conditioning

Vacuum membrane air dehumidification has gained significant interest in recent years as a highly efficient means of air dehumidification. Prior theoretical modeling work by the authors introduced the Active Membrane Energy Exchanger (AMX) concept, which combines active heat exchange and vacuum membrane dehumidification into one non-isothermal system, and found that it could outperform other air conditioning technologies under many conditions. However, no experimental literature exists on combining active heat exchangers and vacuum membrane dehumidification. Here, the goal of the present study is to evaluate the dehumidification (mass transfer) performance of the AMX concept relative to isothermal membrane dehumidification through three main methodologies: (1) membrane material-level characterization, (2) experimental prototype development and testing, and (3) computational fluid dynamics (CFD) simulations. The dense membranes used in this work showed up to a 40% increase in water vapor permeance at cooler temperatures, and the prototype system showed up to a 6% increase in humidity removal when the air was simultaneously cooled. Furthermore, the membrane module-integrated heat exchange tubes provided additional mixing and turbulence, leading to a 4–8% increase in humidity removal. The upper limit coefficient of performance was equal to approximately 2.2, necessitating advanced system designs to improve efficiency. This study is the first to highlight that combining the cooling and vacuum dehumidification processes can improve dehumidification performance.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2023

Advanced two-phase subchannel method via non-linear iteration

A fast-running, robust two-phase flow, sub-channel model is presented based on non-linear solution of the steady-state subchannel fluid flow equations. The drift-flux model solves for conservation of liquid and vapor mass, mixture energy, and axial and transverse mixture momentum as part of an efficient planar marching scheme and nonlinear, nested outer and inner iteration. Here, models based on mechanistic subcooled boiling, two-phase turbulent void mixing, and drift are included. Solution verification and mesh convergence studies were performed for modern GE 10 × 10 fuel geometry and are shown to have excellent convergence behavior. Run time performance for a 50 axial mesh model showed 2.2 seconds on a single CPU core to tightly converge all 3D distributions (flow, void, pressure) for the GE 10 × 10 fuel geometry, supporting its efficient use within the Virtual Environment for Reactor Applications boiling water reactor framework.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

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↗
DOE OSTI2020

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.↗
DOE OSTI2025

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↗
DOE OSTI2020

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↗