Search NASA⌕ Search

SEARCH · Search NASA

Results for “TAYLOR INSTABILITY”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Effects of a refractory metal protective barrier on uranium carbonitride ceramic–metallic composites exposed to a high temperature hydrogen environment

Uranium nitride has been under consideration at NASA and the Department of Energy for nuclear thermal propulsion applications intended for deep space exploration vehicles. It is critical to the success of the technology development that the behavior of the fuel system is properly characterized in operational environments. This body of work fabricated ceramic metal composites consisting of uranium ceramic fuel in a molybdenum and molybdenum–tungsten metal matrix, both with and without an exterior molybdenum barrier. Specimens were introduced to a high temperature hydrogen environment and characterized through x-ray diffraction, scanning electron microscopy, and energy dispersive spectroscopy to identify characteristic changes. A portion of the specimens was fabricated with a 1.27 mm thick molybdenum barrier in order to isolate the material from the hydrogen environment. Instability of the uranium nitride occurred as low as 2073 K, leading to aggressive instability approaching 2658 K. Carbon reactions with the metal matrix were observed as well as alloying between the uranium and metal matrix constituents. The liquid uranium resulting from the instability of the uranium nitride was observed to attack the grain boundaries of the refractory metal, resulting in the formation of a liquid uranium molybdenum eutectic melt. Specimens encased in a barrier experienced similar results. The authors conclude this is a significant challenge to applications of this cermet at these temperatures. Alternate configurations with resilient hermetic seals or lower temperature applications such as power reactors may have more success.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Causes of fuel–ablator mix inferred from modeling of monochromatic time-gated radiography of OMEGA cryogenic implosions

Here, we present evidence, in the context of OMEGA cryogenic target implosions, that laser imprint, known to be capable of degrading laser-direct-drive target performance, plays a major role in generating fuel–ablator mix. OMEGA cryogenic target implosions show a performance boundary correlated with acceleration-phase shell stability; for sufficiently low adiabats (where the adiabat is the ratio of the pressure to the Fermi pressure) and high in-flight aspect ratios (IFAR's), the neutron-weighted shell areal density and neutron yield relative to the clean simulated values sharply decline. Direct evidence of Rayleigh–Taylor fuel–ablator mixing was previously obtained using a Si Heα backlighter driven by an ~20-ps short pulse generated by OMEGA EP. The shadow cast by the shell shortly prior to stagnation, as diagnosed using backlit radiographs, shows a softening near the limb, which is evidence of an ablator–fuel mix region for a low-adiabat implosion (α ~ 1.9, IFAR = 14) but not for a moderate adiabat implosion (α ~ 2.5, IFAR = 10). We find good agreement between experimental and synthetic radiographs in simulations that model laser imprint and account for uncertainty in the initial ablator thickness. Furthermore, we further explore the role of other mechanisms such as classical instability growth at the fuel–ablator interface, species concentration diffusion, and long-wavelength drive and target asymmetries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Three-dimensional Core-collapse Supernova Simulations with 160 Isotopic Species Evolved to Shock Breakout

We present three-dimensional simulations of core-collapse supernovae using the FLASH code that follow the progression of the explosion to the stellar surface, starting from neutrino radiation hydrodynamic simulations of the neutrino-driven phase performed with the Chimera code. We consider a 9.6 M ⊙ zero-metallicity progenitor starting from both 2D and 3D Chimera models and a 10 M ⊙ solar-metallicity progenitor starting from a 2D Chimera model, all simulated until shock breakout in 3D while tracking 160 nuclear species. The relative velocity difference between the supernova shock and the metal-rich Rayleigh–Taylor (R-T) "bullets" determines how the metal-rich ejecta evolves as it propagates through the density profile of the progenitor and dictates the final morphology of the explosion. We find maximum 56 Ni velocities of ~1950 and ~1750 km s -1 at shock breakout from 2D and 3D 9.6 M ⊙ Chimera models, respectively, due to the bullets' ability to penetrate the He/H shell. When mapping from 2D, we find that the development of higher-velocity structures is suppressed when the 2D Chimera model and 3D FLASH model meshes are aligned. The development of faster-growing spherical-bubble structures, as opposed to the slower-growing toroidal structure imposed by axisymmetry, allows for interaction of the bullets with the shock and seeds further R-T instabilities at the He/H interface. We see similar effects in the 10 M ⊙ model, which achieves maximum 56 Ni velocities of ~2500 km s -1 at shock breakout.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A Preliminary Feasibility Study on Multi-Cavity Cryomodule Integration for the Electron Ion Collider Energy Recover Linac Cooler

The Electron-Ion Collider (EIC) is a state-of-the-art accelerator designed for collisions between highly polarized electrons and ions. To achieve optimal luminosity, the ion beam is cooled using an electron beam sourced from an energy recovery linac (ERL). In the current ERL baseline design, one BNL type RF cavity is used per cryomodule, leading to spatial and cost challenges. This study examines the feasibility of using more compact PERLE (Powerful Energy Recovery Linac for Experiments) type cavities, which can house multiple cavities within a single cryomodule, by evaluating their Beam Breakup (BBU) instability performance. Higher Order Modes (HOM) parameters were obtained through frequency scaling, assuming constant quality factor Q_L, R/Q, and linearly scaling HOM frequencies. To predict the BBU threshold current and ensure accuracy, two different BBU tracking simulations are used for cross-verification. Although the reduced footprint of PERLE-type cavities is advantageous, maintaining sufficient damping of HOMs remains crucial. Finally, we compare the HOM damping efficiency of both cavity types and suggest a pathway forward.

Deitrick, K.↗
DOE OSTI2020

Beyond Strain: Controlling the Surface Chemistry of CsPbI 3 Nanocrystal Films for Improved Stability against Ambient Reactive Oxygen Species

Colloidal halide perovskite nanocrystals (NCs) have the possibility of easy scale-up due to their batch synthesis and have demonstrated excellent optoelectronic properties. In particular, perovskite NCs have remarkably high photoluminescence quantum yields in solution and as thin films and impressive open circuit voltages in photovoltaic devices. Despite these promising results, little work has been done to understand the stability of CsPbI 3 NCs for optoelectronic device applications. It has been previously shown that the ligands impart tensile surface strain, which stabilizes the black three-dimensional (3D) perovskite phase against phase degradation, making CsPbI 3 NCs some of the most structurally robust inorganic halide perovskites to date. However, understanding exactly how CsPbI 3 NCs degrade under ambient conditions is critical. Additionally, we demonstrate that the degradation mechanism of NCs is unique from, and 2 orders of magnitude slower than, their polycrystalline thin-film counterparts. Under specific conditions, CsPbI 3 NC films show a compositional instability instead of the phase instability seen in large grain CsPbI 3 . This is mediated through reactions with superoxide and other reactive oxygen species, which are initiated from surface defect states, O 2 and light. We then use this mechanistic insight to identify multiple strategies to prolong the lifetimes of CsPbI 3 NC films, by going beyond surface strain to mitigate key surface chemistries. We demonstrate that (1) minimizing the number of surface defects (2) using an alkylammonium bromide ligand surface treatment and (3) encapsulation with an oxygen scavenging layer all increase NC film lifetimes by inhibiting various steps in the photo-oxidation degradation reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

High-temperature carbon dioxide capture in a porous material with terminal zinc hydride sites

Carbon capture can mitigate point-source carbon dioxide (CO 2 ) emissions, but hurdles remain that impede the widespread adoption of amine-based technologies. Capturing CO 2 at temperatures closer to those of many industrial exhaust streams (>200°C) is of interest, although metal oxide absorbents that operate at these temperatures typically exhibit sluggish CO 2 absorption kinetics and instability to cycling. Here, in this study, we report a porous metal–organic framework featuring terminal zinc hydride sites that reversibly bind CO 2 at temperatures above 200°C—conditions that are unprecedented for intrinsically porous materials. Gas adsorption, structural, spectroscopic, and computational analyses elucidate the rapid, reversible nature of this transformation. Extended cycling and breakthrough analyses reveal that the material is capable of deep carbon capture at low CO 2 concentrations and high temperatures relevant to postcombustion capture.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Magnetized liner inertial fusion platform development to assess performance scaling with drive parameters

Magnetized liner inertial fusion (MagLIF) experiments have demonstrated fusion-relevant ion temperatures up to 3.1 keV and thermonuclear production of up to 1.1 × 1013 deuterium–deuterium neutrons. This performance was enabled through platform development that provided increases in applied magnetic field, coupled preheat energy, and drive current. Advanced coil designs with internal reinforcement enabled an increase from 10 to 20 T. An improved laser pulse shape, beam smoothing, and thinner laser entrance foils increased preheat energy coupling from less than 1 to 2.3 kJ. A redesign of the final transmission line and load region increased peak load current from 16 to 20 MA. The wider range of input parameters was leveraged to study target performance trends with preheat energy, applied magnetic field, and peak load current. Ion temperature and neutron yield generally followed trends in two-dimensional clean Lasnex calculations. Stagnation performance improved with peak load current when other input parameters were also increased such that convergence was maintained. This dataset suggests that reducing convergence to less than 30 would improve predictability of target performance. Lasnex was used to identify a simulation-optimized scaling path, which suggests 10+ kJ of fusion yield is possible on the Z facility with achievable input parameters. This path also indicates >10 MJ could be generated through volume burn on a future facility with a path to high yield (>200 MJ) using cryogenic dense fuel layers. The newly developed MagLIF platform enables exploration of both this simulation optimized scaling path and a recently developed similarity-scaling path.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Determination of Phenolic Groups in Bio-Oils Using Revised Folin-Ciocalteu Methods: Single Cuvette and Plate Reader. Laboratory Analytical Procedure (LAP), Issue Date: May 12, 2022

Phenol components are ubiquitous in wood-derived bio-oils and biocrudes. Their reaction with other functional groups (e.g., aldehydes) may contribute to the formation of carbonaceous species and the expected thermal instability of pyrolysis oils. During hydrotreating, phenols can be recalcitrant species, requiring higher reaction temperatures than other oxygen-containing functional groups. Phenols are also present in upgraded products and have been shown to lead to catalyst deactivation during hydrotreating. They also are the first oxygenated functional group to re-appear in the upgraded product, signaling catalyst deactivation. This procedure covers the determination of phenolic compounds in fast pyrolysis oils. This Laboratory Analytical Procedure (LAP) includes two methods, the first allowing for shorter analysis time at increased reaction temperature, while the second employs a longer analysis time but at room temperature. Additionally, the use of both a single cuvette and of a plate reader are also presented.

09 BIOMASS FUELS↗
DOE OSTI2026

Analytical model for balling defects in laser melting using rivulet theory and solidification

In laser welding and additive manufacturing communities, the balling (humping) defect is primarily attributed to the Plateau–Rayleigh fluid instability (PRI) with a few authors suggesting fluid jetting and volume conservation as alternative mechanisms. As analytical descriptions of these mechanisms are unavailable, combining them into a single formalism is unfeasible. We present a new model of PRI with higher accuracy, accounting for competition with solidification, to compare the expected behavior with known experimental trends when fluid jetting is neglected. We adapt a rivulet instability model from fluid physics to account for the stabilizing effects of the substrate, which the traditional cylindrical-jet geometry does not account for and estimate the instability growth rate. Our model yields a continuous transition from non-balling to balling hitherto lacking in current literature and predicts instability growth at higher wavelengths with strong sensitivity to the solidification front curvature. While the fluid surface is most unstable for shallow melt pools, the absolute magnitude of balling relevant to printing defects scales with melt pool depth and has a maximum for a given melt pool geometry. Synchrotron-based x-ray radiography of thin samples indicates that PRI growth rates and solidification can be comparable in magnitude and thus compete, as we find in this work. We predict that deviations between model predictions and our experimental results demonstrate the importance of fluid flows and heat transport in the balling process. Our experiments further demonstrate at least one mechanism by which the melt pool length and the balling wavelength are not equivalent, as commonly claimed.

Taylor, Z. [Stanford Univ., CA (United States); SL↗
DOE OSTI2021

Scaling of magnetic dissipation and particle acceleration in ABC fields

Using particle-in-cell numerical simulations with electron–positron pair plasma, we study how the efficiencies of magnetic dissipation and particle acceleration scale with the initial coherence length λ 0 in relation to the system size L of the two-dimensional ‘Arnold–Beltrami–Childress’ (ABC) magnetic field configurations. Topological constraints on the distribution of magnetic helicity in two-dimensional systems, identified earlier in relativistic force-free simulations, that prevent the high-(L/λ 0 ) configurations from reaching the Taylor state, limit the magnetic dissipation efficiency to about ϵ diss ≃60% . Here we find that the peak growth time scale of the electric energy τ E,peak scales with the characteristic value of initial Alfvén velocity β A,ini like τ E,peak ∝(λ 0 /L)β$^{–3}_{A,ini}$ . The particle energy change is decomposed into non-thermal and thermal parts, with non-thermal energy gain dominant only for high initial magnetisation. The most robust description of the non-thermal high-energy part of the particle distribution is that the power-law index is a linear function of the initial magnetic energy fraction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A 2–4 keV multilayer mirrored channel for the NIF Dante system

During inertial confinement fusion experiments at the National Ignition Facility (NIF), a capsule filled with deuterium and tritium (DT) gas, surrounded by a DT ice layer and a high-density carbon ablator, is driven to the temperature and densities required to initiate fusion. In the indirect method, 2 MJ of NIF laser light heats the inside of a gold hohlraum to a radiation temperature of 300 eV; thermal x rays from the hohlraum interior couple to the capsule and create a central hotspot at tens of millions degrees Kelvin and a density of 100–200 g/cm3. During the laser interaction with the gold wall, m-band x rays are produced at ∼2.5 keV; these can penetrate into the capsule and preheat the ablator and DT fuel. Preheat can impact instability growth rates in the ablation front and at the fuel–ablator interface. Monitoring the hohlraum x-ray spectrum throughout the implosion is, therefore, critical; for this purpose, a Multilayer Mirror (MLM) with flat response in the 2–4 keV range has been installed in the NIF 37° Dante calorimeter. Precision engineering and x-ray calibration of components mean the channel will report 2–4 keV spectral power with an uncertainty of ±8.7%.

Rubery, M. S. (ORCID:0000000344152818)↗
DOE OSTI2025

Amplification of laser imprint in the presence of strong, externally imposed, target-normal magnetic fields

Here, an experiment was performed on the OMEGA EP Laser System to investigate the effects of magnetic fields on laser imprint. A 30 µ⁢m thick CH target was driven by a single beam delivering 3.0 kJ of UV energy in a 5 ns square pulse without smoothing by spectral dispersion. The Rayleigh-Taylor amplification of this beam's imprint on the target surface was monitored using face-on gated x-ray radiography from a Gd backlighter. Magnetic fields of up to 45 T were applied normal to the target surface. Analysis of the resulting radiographs shows a 60 ± 13% increase in the spectrally resolved surface perturbation amplitudes, consistent at all times and for all unsaturated frequencies. This consistency indicates that the increase in perturbation amplitudes was due to a change in the initial amplitudes rather than Rayleigh-Taylor growth. This is supported by the trajectory of individual modes and the inferred time-averaged perturbation growth rates. Laser imprint is therefore inferred to have increased due to strong magnetic fields that remain in the conduction zone of the target, suppressing off-axis electron motion and limiting the effects of thermal smoothing.

direct drive↗
DOE OSTI2024

How numerical simulations helped to achieve breakeven on the NIF

The inertial confinement fusion program relies upon detailed simulations with inertial confinement fusion (ICF) codes to design targets and to interpret the experimental results. These simulations treat as much physics from essential principles as is practical, including laser deposition, cross beam energy transfer, x-ray production and transport, nonlocal thermal equilibrium kinetics, thermal transport, hydrodynamic instabilities, thermonuclear burn, and transport of reaction products. Improvements in radiation hydrodynamic code capabilities and vast increases in computing power have enabled more realistic, accurate 3D simulations that treat all known asymmetry sources. We describe how numerical simulations helped to guide the program, assess the impediments to breakeven, and optimize every aspect of target design. A preshot simulation of the first National Ignition Facility experiment that surpassed breakeven predicted an increased yield that matches the experimental result, within the preshot predicted uncertainty, with a target gain of 1.5. We will cover the key developments in Lawrence Livermore National Laboratory ICF codes that enabled these simulations and give specific examples of how they helped to guide the program.

Marinak, M. M. (ORCID:0009000331127696)↗
DOE OSTI2021

Entropy generation from hydrodynamic mixing in inertial confinement fusion indirect-drive targets

The increase in entropy from the physical mixing of two adjacent materials in inertial confinement fusion (ICF) implosions and gas-filled hohlraums is analytically assessed. An idealized model of entropy generation from the mixing of identical ideal-gas particles across a material interface in the presence of pressure and temperature gradients is applied. Physically, mix-driven entropy generation refers to the work done by the gases in expanding into a larger common volume from atomic mixing under the condition of no internal energy change, or work needed to restore the initial unmixed state. Furthermore, the effect of a mix-generated entropy increase is analytically shown to lead to less compression of the composite ICF fluid under adiabatic conditions. The amount of entropy generation is estimated to be ~10 J for a Rayleigh–Taylor-induced micrometer-scale annular mixing layer between the solid deuterium–tritium fuel and (undoped) high-density carbon pusher of an imploding capsule at the National Ignition Facility (NIF). This level of entropy generation is consistent with lower-than-expected fuel compressions measured on the NIF [Hurricane et al., Phys. Plasmas 26, 052704 (2019)]. The degree of entropy increase from mixing of high-Z hohlraum wall material and low-Z, moderate- to high-density gas fills is estimated to lead to ~100 kJ of heat generation for NIF-scale experiments [Moody et al., Phys. Plasmas 21, 056317 (2014)]. This value represents a significant fraction of the inferred missing x-ray drive energy based on observed delays in capsule implosion times compared with mainline simulations [Jones et al., Phys. Plasmas 19, 056315 (2012)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Evolution of the electrothermal instability from thick rod z pinches subject to dynamically and statically applied axial magnetic field

LDRD Project 229427 aimed to determine how electrothermal instability (ETI) driven heating on a z-pinch rod pulsed with intense current evolves under mixed magnetic field (azimuthal + axial) conditions, which is pertinent to pulsed-power-driven magnetically-insulated transmission lines and physics targets. Experiments focused on diagnosing ETI-driven heating from deliberately-machined and well-characterized micron-scale surface defects (referred to as engineered defects or ED). Prior to the start of this project, understanding of how unmagnetized (B z =0) ED evolve had been obtained—simulations largely reproduce the experimentally observed high temperature spots which develop at the poles of bare/uncoated ED. Project 229427 extended the Mykonos Facility ED experimental platform to include axial field. In the first class of experiments, axial field was provided “dynamically” via a helical return can (HRC). In this case, B z and B θ rise at the same rate. Generally, the HRC generated magnetic field at a fixed polarization angle Φ B =arctan(B z /B θ )=15° on the rod's surface. In the second class of experiments, axial field was provided “statically” via a slow-rising (millisecond) external Helmholtz coil pair. In this case, B z was effectively constant/static throughout the 100 ns rise of the Mykonos current. For either case, a primary goal was to determine whether ETI provides a helical seed perturbation for the subsequent growth of the helical magneto Rayleigh-Taylor modes observed in MagLIF (static B z ) and dynamic screw pinch (DSP, dynamic B z ) experiments. When dynamic field was applied using an HRC, emissions from individual ED aligned toward Φ B , while emissions from ED within pairs elongated and preferentially merged along Φ B . These data strongly support that for a randomized defect distribution, heating from nearby current-density perturbations will favorably merge about Φ B to generate an extended seed perturbation that aligns toward the surface-field polarization, and this may impact the orientation of subsequent MRT growth on imploding liners. The results from the static field experiments were largely inconclusive, as any ETI heating rotation, if present, was obscured/overwhelmed by local/random heating from ED rim imperfections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Hydrostatic and Non‐Hydrostatic Baroclinic Instability in the Dynamical Core of the DOE Global Climate Model

Abstract The dynamical core of the Department of Energy global climate model is used to understand the role of non‐hydrostatic dynamics in the simulation of dry and moist baroclinic waves. To reduce computational cost, the Diabatic Acceleration and REscaling approach is adopted. Scale analysis and numerical simulations suggest that the model solution is not distorted by the change of spherical metric terms due to the change of Earth radius, neither are the inter‐scale interactions strongly altered as indicated by the analysis of spectral flux. Compared with hydrostatic simulations, the onset of baroclinic instability is delayed under non‐hydrostatic dynamics as the associated weaker vertical motions tend to increase the critical wavelength and narrow the range of waves that can be baroclinically unstable. During the development of baroclinic waves, non‐hydrostatic dynamics tends to induce vertical motions in the upper troposphere, accelerating the eastward propagation of upper‐level ridges through their impact on local vorticity tendency. These processes reduce the westward tilt of the vertical ridge axes and suppress the conversion of mean flow available potential energy to eddy kinetic energy, leading to weaker baroclinic eddies than in hydrostatic simulations. The contrasts between hydrostatic and non‐hydrostatic settings hold in supplemental experiments that vary the background flow Rossby number, the amount of water vapor content and vertical resolution. We also find that mesoscale and smaller‐scale activities can be considerably under‐represented when the vertical resolution is limited to that typically used in global climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Evaluation of a nonlinear variational multiscale method for fluid transport problems

Diverse transport problems, especially those based on fluid flow models, are intrinsically multiscale and nonlinear, characteristics that often lead to intricate dynamics such as the development of instabilities and turbulence. Computational simulations that resolve all scales in these problems are often unfeasible, prompting to coarse-grained simulation strategies in which small-scale features are modeled instead of resolved. Variational Multiscale (VMS) methods, and particularly residual-based Large-Eddy Simulation (LES) approaches, have proven effective and robust for the coarse-grained simulation of complex transport problems. VMS methods avoid the assumption of separable nonlinearity and the reliance on empirical small-scale models by using a variational decomposition of scales together with a residual-based approximation of the small-scales. Evaluation of a nonlinear VMS approach, denoted as VMS n , is presented for the coarse-grained simulation of transient-advective-diffusive-reactive (TADR) transport problems arising from fluid flow models. In contrast to classical VMS approaches that neglect the effect of the small scales on the transport operator, VMS n treats the inter-dependence between large- and small-scales upfront. The treatment of inter-scale coupling involves the solution of a local algebraic nonlinear system describing the evolution of the small-scales. The VMS n approach is complemented with two algebraic approximations of the small-scales: one based on the main diagonal of the transport matrices and another that preserves transport fluxes and is suitable for generic TADR systems. The suitability of the VMS n approach for handling general TADR problems and regimes is evaluated with benchmark incompressible, compressible, and magnetohydrodynamic laminar flow problems, the incompressible Taylor-Green vortex flow, the turbulent free jet, and the two-temperature arc in crossflow. Furthermore, simulation results show that VMS n leads to minor improvements in accuracy with respect to the classical VMS for the laminar flow problems, but to significantly greater accuracy for the turbulent flows and the unsteady plasma flow problems, while using the same cohesive numerical formulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗