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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Mixing length in low Reynolds number compressible turbulent boundary layers

The paper studies the effect of low Reynolds number in high-speed turbulent boundary layers on variations of mixing length. Boundary layers downstream of natural transition on plates, cones and cylinders, and boundary layers on nozzle walls without laminarization-retransition are considered. The problem of whether low Reynolds number amplification of shear stress is a result of transitional flow structure is considered. It is concluded that a knowledge of low Reynolds number boundary layer transition may be relevant to the design of high-speed vehicles.

Bushnell, D. M.↗
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↗
NASA NTRS2019

CFD Jet Mixing Model Validation Against Zero-Boil-Off Tank (ZBOT) Microgravity Experiment

The Zero-Boil-Off Tank (ZBOT) experiments are a series of small scale two-phase fluids physics experiments conducted aboard the International Space Station (ISS) in 2017. A goal of these experiments is to study various fundamental fluid flow, heat and mass transport, and phase change phenomena that control storage tank pressurization and pressure control in microgravity. In this paper, results of a state-of-the-art two-phase CFD (Computational Fluid Dynamics) storage tank model are presented and validated against jet mixing data provided by microgravity ZBOT experiments. Jet-ullage interaction was studied in the experiment for different jet speeds representing different flow regimes, and tank fill levels. Effects of various physical parameters, such as, jet angle, gravitational acceleration and surface tension value were investigated. Computational parameters, such as turbulence and surface tension models and interface capturing scheme, affecting jet-ullage interaction were also studied. Computational results were compared to experimental results for model validation.

Turbulence↗
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↗
NASA NTRS1979

On turbulence modeling for unsteady transonic flows

A detailed examination of the turbulent field in an unsteady transonic flow undergoing shock-induced separation is conducted. Ensemble-averaged mean and fluctuating velocities, obtained from conditionally sampled laser velocimeter data, are described and analyzed to assess the applicability of modeling concepts usually employed in steady-flow problems. Some comparisons with computations employing the Reynolds averaged Navier-Stokes equations with a mixing length turbulence model are then presented to illustrate the status of current predictive capabilities. The results appear to imply that turbulence models developed for steady flows apply and that the model need not reflect all the fine details of the turbulent structure but rather account in an approximate way for the production and destruction of the turbulence.

Marvin, J. G.↗
NASA NTRS1984

Dynamic response of shock waves in transonic diffuser and supersonic inlet - An analysis with the Navier-Stokes equations and adaptive grid

An existing method which solves the multi-dimensional ensemble-averaged compressible time-dependent Navier-Stokes equations in conjunction with mixing length turbulence model and shock capturing technique has been extended to include the shock-tracking adaptive grid systems. The numerical scheme for solving the governing equations is based on a linearized block implicit approach. The effects of grid-motion and grid-distribution on the calculated flow solutions have been studied in relative detail and this is carried out in the context of physically steady, shocked flows computed with non-stationary grids. Subsequently, the unsteady dynamics of the flows occurring in a supercritically operated transonic diffuser and a mixed compression supersonic inlet have been investigated with the adaptive grid systems by solving the Navier-Stokes equations.

Liu, N.-S.↗
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↗
NASA NTRSDate not supplied

Jet Noise Modeling Using an Acoustic Analogy: From Theoretical Formulation to Practical Noise Predictions

In this presentation an overview of the development of an acoustic analogy formulation, and its implementation into a reduced-order prediction method for turbulent jet mixing noise, will be given. The definition and rationale of an acoustic analogy approach to the prediction of flow-generated noise will be provided. A high-level description of the mathematical processes involved will be outlined in the context of a particular formulation developed at NASA Glenn Research Center. The modelling approximations needed to implement the general theoretical formulation into a practical noise prediction tool will be discussed. Sample results obtained using a noise prediction code based on the NASA Glenn formulation, and comparisons with experimental data, will be presented for a series of turbulent jet test cases. A view on the prospects for continued use of acoustic analogy-based methods for the prediction of flow-generated noise and the assessment of noise-reduction technologies will conclude the presentation.

Jet noise↗
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↗
NASA NTRS1974

Calculation of compressible turbulent free shear layers

Rudy and Bushnell's (1972) mixing length theory for the prediction of mean velocity fields in turbulent shear layers is extended to include supersonic free turbulent mixing. The normal momentum equation is coupled with the conventional equations of motion and is solved iteratively. Results show that transverse static pressure variation has very little effect on the mean flow variables in compressible free turbulent shear layers.

Oh, Y. H.↗
NASA NTRS1990

Application of mixing-controlled combustion models to gas turbine combustors

Gas emissions were studied from a staged Rich Burn/Quick-Quench Mix/Lean Burn combustor were studied under test conditions encountered in High Speed Research engines. The combustor was modeled at conditions corresponding to different engine power settings, and the effect of primary dilution airflow split on emissions, flow field, flame size and shape, and combustion intensity, as well as mixing, was investigated. A mathematical model was developed from a two-equation model of turbulence, a quasi-global kinetics mechanism for the oxidation of propane, and the Zeldovich mechanism for nitric oxide formation. A mixing-controlled combustion model was used to account for turbulent mixing effects on the chemical reaction rate. This model assumes that the chemical reaction rate is much faster than the turbulent mixing rate.

Nguyen, Hung Lee↗
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↗
NASA NTRS2023

Variably Premixed Rotating Detonation Engine for Evaluation of Detonation Cycle Dynamics

A variably premixed rotating detonation engine using gaseous hydrogen and air reactants is introduced to enable investigation of key cycle processes while varying the homogeneity of the reactant inlet conditions. Two chamber configurations are investigated, the first with reactants filling the entire span of a straight annular channel and the second with a slightly larger channel and a backward-facing step. The first configuration permits both premixed and non-premixed fuel injection, enabling mixing quality modulation. The second configuration is operated only at fully premixed conditions. Operating modes and detonation wave speeds are characterized using exhaust-plume imaging, while the chamber heat release field is captured by transverse imaging through a transparent outer body. Tests using the first configuration were characterized by unstable counterpropagating modes with low detonation wave speeds regardless of the state of premixing, while the second configuration rendered single-wave behavior with wave speeds up to 86% of the Chapman–Jouguet velocity. Comparisons with a simple computational fluid dynamics model of the second configuration indicate that reactant preheating significantly influences the detonation wave topology, highlighting the potential utility of the test platform for isolating key physics associated with the effects of reactant premixing, preheating, and chamber geometry on rotating detonation engine operation.

Detonation↗
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↗