Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “initial state”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Deciphering competing elementary steps to correlate electrocatalyst chemical state with activity

The overpotential in multielectron transfer heterogeneous electrocatalysis fundamentally arises from thermodynamic and kinetic disparities among elementary steps; however, deciphering coupled and competing steps has long remained a challenge. Here, we establish an electrochemical deconvolution paradigm based on key processes in electrocatalytic reactions, such as charge accumulation, electron/proton transfer, and intermediate evolution, to resolve competing elementary steps. Taking the oxygen evolution reaction as a prototypical reaction, we design a model catalyst featuring a precise isolated cation-anion vacancy pair and track the previously elusive electrochemical behavior of lattice oxygen by disentangling interference from adsorbed oxygen intermediates. Mechanistically, the lattice oxygen oxidation pathway originates from the spontaneous, nonelectrochemical deprotonation of replenished water molecules coordinated to unsaturated cation sites. Alternating current techniques further reveal that although lattice oxygen oxidation requires a higher potential than metal oxidation, it exhibits faster kinetics, providing insight into its superior catalytic activity. These findings establish a direct experimental correlation between the initial chemical state and the catalytic activity and prove surface-confined lattice oxygen cycling. Furthermore, expanding conventional potential-current analysis into a multidimensional framework enables disentanglement of thermodynamic and kinetic contributions of key elementary steps, thereby guiding the rational optimization of various complex multielectron transfer reactions.

OER↗
DOE OSTI2025

Operando APXPS for direct probing of Li ion battery LCO electrode/electrolyte interface chemistry during lithiation/delithiation

The real-time interface chemistry between the lithium cobalt oxide (LCO) working electrode and the LiClO 4 /propylene carbonate (PC) electrolyte is investigated during lithiation/delithiation using dip-and-pull ambient pressure photoelectron spectroscopy (APXPS). The APXPS results appear to exhibit the seldom discussed Co 2+ state in the LCO structure, where the operando measurements indicate electron transfer among Co 2+ , Co 3+ , and Co 4+ states. Specifically, the lithiation of LCO reduces the Co 4+ state to both Co 3+ and Co 2+ states, where, as a function of voltage, reduction to the Co 2+ state is initially more pronounced followed by Co 3+ formation. In addition, a delay in surface delithiation is observed during the reverse potential steps. This is discussed in terms of overpotential at the interface measurement position as a consequence of the dip-and-pull setup for this experiment. Finally, the shifts in the apparent binding energies of the spectral features corresponding to the electrolyte and LCO at their interface show that the electrochemical potentials at delithiation voltage steps are different from the lithiation steps at the same applied voltages. This further explains the non-responsive delithiation. The BE shift observed from the LCO surface is argued to be dominantly due to the semi-conductive nature of the sample. Overall, this article shows the importance of operando APXPS for probing non-equilibrium states in battery electrodes for understanding electron transfer in the reactions.

Liu, Qianhui [Uppsala Univ. (Sweden)] (ORCID:00000↗
DOE OSTI2024

Search for Θ+ in KLp→K+n reaction in KLF at JLab

The possibility of the existence of multiquark hadrons made of 4-quarks for mesons and 5-quarks for baryons was predicted by Gell-Mann [Phys. Lett. 8, 214 (1964)]. The renewed interest for the search for exotic pentaquark states was initiated by the paper by Diakonov, Petrov, and Polyakov [Z. Phys. A 359, 305 (1997)]. The 2003 experimental reports on the observation of [Formula: see text] pentaquark with [Formula: see text] quark content created a great excitement and many following experiments have reported its observation [K. H. Hicks, Eur. Phys. J. H 37, 1 (2012)]. After high-statistics experiments at JLab, which did not confirm previous claims by the CLAS collaboration, the community concluded that the [Formula: see text] pentaquark either does not exist at all or has an extremely small cross-section, making it currently unobserved. There were different review papers on this subject, questioning the existence of [Formula: see text] or trying to explain the reasons why reaching a conclusion based on production experiments is challenging [M. Amaryan, Eur. Phys. J. Plus 137, 684 (2022)]. To address the challenge of minimal 3-body final states, a formation experiment with a projectile-kaon beam is proposed. In the following, we discuss how the [Formula: see text] could be observed in the [Formula: see text] reaction in the KLF experiment at JLab [M. Amaryan et al. [KLF Collaboration], arXiv:2008.08215 [nucl-ex]].

Astronomy & Astrophysics↗
DOE OSTI2021

Influence of ethylenediaminetetraacetic acid on the long-term oxidation state distribution of plutonium

Here, spectrophotometry was used to study the effect of EDTA on plutonium oxidation state distribution as a function of time, pH, and ligand-to-metal ratio (L/M) under anoxic conditions. Novel Pu(V)-EDTA absorption bands were identified at 571, 993, 1105, and 1150 nm with molar absorption coefficients of 15 ± 1, 6 ± 1, 10 ± 1, and 10 ± 1 cm –1 M –1 , respectively. Pu(V)-EDTA spectral changes occurred at L/M < 1, indicating only Pu V O 2 (EDTA) 3- formed with logK = 3.6 ± 0.3. Time-resolved experiments showed EDTA drastically increased the Pu(V/VI) reduction rate, which we propose is driven by amine lone-pair electron donation and the oxidative decarboxylation of EDTA. Oxidation of Pu(III)-EDTA to Pu(IV)-EDTA occurred on a slower time scale (110–237 days) than previously reported (<15 min) and is hypothesized to be radiolysis driven. Pu(V/VI)-EDTA and Pu(III)-EDTA both approached Pu(IV)-EDTA stabilization over time, yet Pu(V/VI)-EDTA solubility data was ≥ 1.0 log 10 units higher than predicted by Pu(IV)-EDTA solubility models, indicating that current thermodynamic models are incomplete. Ultimately, the data show EDTA preferentially stabilizes Pu(IV) over time regardless of initial oxidation state, but Pu(V)-EDTA can persist under environmentally-relevant conditions, emphasizing the need to continue investigating redox reactions, speciation, and behavior of these complexes to support the transuranic waste disposal and surface remediation/containment efforts.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2026

Xanthos-Lake Dataset

The Xanthos-Lake v1.0 dataset provides the input data, trained machine-learning models, and simulation outputs needed to characterize lake water balance, snow and ice conditions, and mixing-layer temperature within the Xanthos global hydrological modeling framework. The dataset supports lake representation across a wide range of lake sizes and hydroclimatic conditions by combining xLSIM, a basin-specific machine-learning emulator of lake snow, ice, ice-cover fraction, and mixing-layer temperature, with the Xanthos-Lake water-balance model. The archive contains NetCDF datasets used to train and evaluate xLSIM, trained model weights, processed meteorological and lake-property inputs, and basin- and lake-category-specific simulation outputs. These materials are organized into four primary data groups, described below. Snowice_model_inputs: Contains the NetCDF input data used to train xLSIM. The xLSIM machine-learning framework uses three lake-based datasets. The meteorological forcing dataset provides monthly relative humidity, specific humidity, surface wind speed, maximum and minimum air temperature, downward longwave and shortwave radiation, snowfall, surface air pressure, and total precipitation. Lake surface area is included as an additional static predictor. The target-state dataset provides lake ice thickness, snow depth, snow cover, and lake mixing-layer temperature, while a companion lake-surface dataset provides the lake ice-cover fraction. Before training, ice thickness and snow depth are converted from meters to centimeters, mixing-layer temperature is converted from kelvin to degrees Celsius and constrained to nonnegative values, and ice-cover fraction is converted from a fraction to a percentage. The predictor variables are normalized using statistics calculated across the selected lakes and time steps. Snowice_model_outputs: Contains the NetCDF outputs generated by xLSIM. For each basin, xLSIM produces a file containing observed and predicted lake-state variables for the training, validation, and testing periods. The modeled variables include lake ice thickness, snow depth, snow cover, mixing-layer temperature, and lake ice-cover fraction. For basins without a sufficiently persistent snow-and-ice signal, the emulator predicts only mixing-layer temperature. The outputs also include training and validation loss histories, the selected model configuration, identifiers of the lakes used in training, and SHAP-based feature-importance information at the global, lake, and seasonal-regime levels. The trained machine-learning model weights are provided separately within the dataset archive. Together, these files support model evaluation and subsequent coupling with the Xanthos-Lake water-balance framework. XanthosLAKES: Contains the NetCDF input data used by the Xanthos-Lake framework. Monthly meteorological inputs include relative and specific humidity, downward shortwave and longwave radiation, mean, maximum, and minimum air temperature, wind speed, precipitation, snowfall, and surface air pressure. Static lake-property datasets provide lake identifiers, geographic locations, surface area, volume, mean depth, elevation, drainage area, fetch, outlet-routing information, and associated Xanthos grid-cell attributes. Separate bathymetric datasets provide the coefficients of the area–depth and volume–depth relationships for each aggregated lake unit. GLEV-based records provide observed lake surface area and evaporation data used to initialize lake states, define reference conditions, and calibrate and evaluate the model. Xanthos-Lake Outputs: Contains the basin- and lake-category-specific NetCDF outputs generated by Xanthos-Lake. Monthly variables include lake surface area, storage volume, outlet discharge, evaporation rate, evaporation volume, lake–groundwater exchange, lake inflow, ice thickness, snow depth, snow-cover fraction, ice-cover fraction, and mixing-layer temperature. The files also contain lake-specific calibration and validation statistics, including normalized root-mean-square error, mean absolute error, Nash–Sutcliffe efficiency, Kling–Gupta efficiency, and percent bias. Stored calibrated and derived parameters include the weir discharge coefficient, fractional freeboard, groundwater exchange coefficient, reference water level, corresponding reference surface area and storage volume, weir-width adjustment factor, and the fraction of routed inflow entering the lake. Basin identifiers, lake category, simulation period, calibration and validation periods, and parameter-schema information are retained as NetCDF metadata.

Abeshu, Guta [Pacific Northwest National Laborator↗
DOE OSTI2025

Structural and Mechanistic Advances in the Chemistry of Methyl-Coenzyme M Reductase (MCR)

Methane represents 34% of U.S. energy consumption and is a major greenhouse gas related to the global carbon cycle and energy production. However, current industrial practices significantly increase atmospheric methane levels, necessitating a deeper understanding of its biosynthesis and oxidation. Methyl-coenzyme M reductase (MCR) is central to biological methane metabolism, catalyzing the final step of methanogenesis and the first step in anaerobic methane oxidation. It is also a key target for strategies to capture and transform methane into value-added chemicals.The active site of MCR is a buried Ni-based cofactor only accessible by the substrates via a 50 Å long tunnel. Although the Ni(I) state is required to initiate catalysis, capturing this state remains a challenge for the current structural techniques. Recent advances in structural biology using X-ray Free-Electron Laser serial crystallography have provided insights into MCR's inactive Ni(II) state at room temperature and show promise for capturing its active Ni(I) form.Our team has established several critical aspects of the MCR mechanism using a combination of experimental and computational studies. MCR uses CH 3 -SCoM and CoBSH as substrates, producing methane and a disulfide product CoMSSCoB. Kinetic analysis showed that productive substrate binding requires CH 3 -SCoM to bind first, inducing conformational changes that optimize the active site for subsequent CoBSH binding. Following substrate binding, four proposed methane production/oxidation mechanisms were examined, establishing whether the reaction proceeds through an organometallic methyl-nickel(III), methyl anion ion, or methyl radical intermediate. Experimental measurements using CoBSH analogs successfully slowed the reaction, allowing for mechanistic insight that demonstrated the methyl radical pathway, where the initial interactions involve homolytic cleavage of the methyl-sulfur bond, generating a methyl radical that quickly abstracts the thiol hydrogen atom of CoBSH to form methane. Computational studies further confirmed that, compared to other mechanisms, the methyl radical mechanism is thermodynamically more favorable and accessible under physiological conditions.Spectroscopic and computational studies challenged the conventional understanding of substrate binding in MCR by proposing an alternative positioning of CH 3 -SCoM and CoMSSCoB in the active site pocket. The research suggested that CH 3 -SCoM (substrate) and CoMSSCoB (product) bind via their sulfonate groups to the Ni(I) center of cofactor F 430 . This binding allows for the reaction without substrate reorganization in the pocket but would require a long-range electron transfer. Altogether, the work summarized in this review reflects our current understanding of the enzyme's catalytic mechanism and structural dynamics. This is essential for developing efficient methane conversion technologies that could mitigate its environmental impact while harnessing energy-storage properties.

03 NATURAL GAS↗
DOE OSTI2023

Field‐Induced Ferroelectric Phase Evolution During Polarization “Wake‐Up” in Hf 0.5 Zr 0.5 O 2 Thin Film Capacitors

Abstract As an emerging nonvolatile memory technology, HfO 2 ‐based ferroelectrics exhibit excellent compatibility with silicon CMOS process flows; however, the reliability of polarization switching in these materials remains a major challenge. During repeated field programming and erase of the polarization state of initially pristine HfO 2 ‐based ferroelectric capacitors, the magnitude of the measured polarization increases, a phenomenon known as “wake‐up”. In this study, the authors attempt to understand what causes the wake‐up effect in Hf 0.5 Zr 0.5 O 2 (HZO) capacitors using nondestructive methods that probe statistically significant sample volumes. Synchrotron X‐ray diffraction reveals a concerted shift in HZO Bragg peak position as a function of polarization switching cycle number in films prepared under conditions such that they exhibit extremely large (≈3000%) wake‐up. In contrast, a control sample with insignificant wake‐up shows no such peak shift. Capacitance – voltage measurements show evolution in the capacitance loop with switching cycle number for the wake‐up sample and no change for the control sample. Piezoresponse force microscopy measurements are utilized to visualize the domain switching with wake‐up. The combination of these observations clearly demonstrates that wake‐up is caused by a field‐driven phase transformation of the tetragonal phase to the metastable ferroelectric orthorhombic phase during polarization switching of HZO capacitors.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Linker Redox Mediated Control of Morphology and Properties in Semiconducting Iron‐Semiquinoid Coordination Polymers**

Abstract The emergence of conductive 2D and less commonly 3D coordination polymers (CPs) and metal–organic frameworks (MOFs) promises novel applications in many fields. However, the synthetic parameters for these electronically complex materials are not thoroughly understood. Here we report a new 3D semiconducting CP Fe 5 (C 6 O 6 ) 3 , which is a fusion of 2D Fe‐semiquinoid materials and 3D cubic Fe x (C 6 O 6 ) y materials, by using a different initial redox‐state of the C 6 O 6 linker. The material displays high electrical conductivity (0.02 S cm −1 ), broad electronic transitions, promising thermoelectric behavior ( S 2 σ=7.0×10 −9 W m −1 K −2 ), and strong antiferromagnetic interactions at room temperature. This material illustrates how controlling the oxidation states of redox‐active components in conducting CPs/MOFs can be a “pre‐synthetic” strategy to carefully tune material topologies and properties in contrast to more commonly encountered post‐synthetic modifications.

Wang, Lei↗
DOE OSTI2020

Electronic couplings for singlet fission: Orbital choice and extrapolation to the complete basis set limit

Abstract For the search for promising singlet fission candidates, the calculation of the effective electronic coupling, which is required to estimate the singlet fission rate between the initially excited state (S 0 S 1 ) and the multiexcitonic state ( 1 TT, two triplets on neighboring molecules, coupled into a singlet), should be sufficiently reliable and fast enough to explore the configuration space. We propose here to modify the calculation of the effective electronic coupling using a nonorthogonal configuration interaction approach by: (a) using only one set of orbitals, optimized for the triplet state of the molecules, to describe all molecular electronic states, and (b) only taking the leading configurations into consideration. Furthermore, we also studied the basis set convergence of the electronic coupling, and we found, by comparison to the complete basis set limit obtained using the cc‐pV n Z series of basis sets, that both the aug‐cc‐pVDZ and 6–311++G** basis sets are a good compromise between accuracy and computational feasibility. The proposed approach enables future work on larger clusters of molecules than dimers.

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

AIMD‐Based Protocols for Modeling Exciplex Fluorescence Spectra and Inter‐System Crossing in Photocatalytic Chromophores

ABSTRACT This study introduces a computational protocol for modeling the emission spectra of exciplexes using excited‐state ab initio molecular dynamics (AIMD) simulations. The protocol is applied to a model exciplex formed by oligo‐p‐phenylenes (OPPs) and triethylamine (TEA), which is of interest in the context of photocatalytic reduction of . AIMD facilitates efficient sampling of the conformational space of OPP3 and OPP4 exciplexes with TEA, offering a dynamic alternative to previously employed static methods. The AIMD‐based protocol successfully reproduces experimental emission spectra for OPP‐TEA exciplexes, agreeing with previous computational and experimental findings. The results show that AIMD simulations provide an efficient means of sampling the conformational space of these exciplexes, requiring less user input and, in some instances, fewer computational resources than multiple excited‐state optimizations initiated from user‐specified initial structures. The study also evaluates the yield of intersystem crossing (ISC) using AIMD and Landau‐Zener probability. The results suggest that ISC is a minor decay channel for OPP3 and OPP4. This work provides new insights into the structural flexibility and emission characteristics of OPP‐TEA photoredox catalyst systems, potentially contributing to improved design strategies for organic chromophores in reduction applications.

Giudetti, Goran [Department of Chemistry Universit↗
DOE OSTI2026

On the Trotter Error in Many-body Quantum Dynamics with Coulomb Potentials

Efficient simulation of many-body quantum systems is central to advances in physics, chemistry, and quantum computing, with a key question being whether the simulation cost scales polynomially with the system size. Here, in this work, we analyze many-body quantum systems with Coulomb interactions, which are fundamental to electronic and molecular systems. We prove that Trotterization for such unbounded Hamiltonians achieves a 1/4-order convergence rate, with explicit polynomial dependence on the number of particles. The result holds for all initial wavefunctions in the domain of the Hamiltonian, and the 1/4-order convergence rate is optimal, as previous work has numerically demonstrated that it can be saturated by a specific initial ground state. The main challenges arise from the many-body structure and the singular nature of the Coulomb potential. Our proof strategy differs from prior state-of-the-art Trotter analyses, addressing both difficulties in a unified framework. Our analysis treats the Coulomb potential as an unbounded operator without modification or regularization, and does not rely on spatial discretization, making it compatible with both first- and second-quantized circuit constructions.

Fang, Di [Duke Univ., Durham, NC (United States)]↗
DOE OSTI2023

Attaining freshwater and estuarine-water soil saturation in an ecosystem-scale coastal flooding experiment

Abstract Coastal upland forests are facing widespread mortality as sea-level rise accelerates and precipitation and storm regimes change. The loss of coastal forests has significant implications for the coastal carbon cycle; yet, predicting mortality likelihood is difficult due to our limited understanding of disturbance impacts on coastal forests. The manipulative, ecosystem-scale Terrestrial Ecosystem Manipulation to Probe the Effects of Storm Treatments (TEMPEST) experiment addresses the potential for freshwater and estuarine-water disturbance events to alter tree function, species composition, and ecosystem processes in a deciduous coastal forest in MD, USA. The experiment uses a large-unit (2000 m 2 ), un-replicated experimental design, with three 50 m × 40 m plots serving as control, freshwater, and estuarine-water treatments. Transient saturation (5 h) of the entire soil rooting zone (0–30 cm) across a 2000 m 2 coastal forest was attained by delivering 300 m 3 of water through a spatially distributed irrigation network at a rate just above the soil infiltration rate. Our water delivery approach also elevated the water table (typically ~ 2 m belowground) and achieved extensive, low-level inundation (~ 8 cm standing water). A TEMPEST simulation approximated a 15-cm rainfall event and based on historic records, was of comparable intensity to a 10-year storm for the area. This characterization was supported by showing that Hurricane Ida’s (~ 5 cm rainfall) hydrologic impacts were shorter (40% lower duration) and less expansive (80% less coverage) than those generated through experimental manipulation. Future work will apply TEMPEST treatments to evaluate coastal forest resilience to changing hydrologic disturbance regimes and identify conditions that initiate ecosystem state transitions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Influence of native oxide film on corrosion behavior of additively manufactured stainless steel 316L

The influence of the native oxide film on passive film properties and localized corrosion of additively manufactured SS 316L was studied in 1 wt% HCl by XPS characterization, electrochemical polarization curves, and post-test morphology analysis by SEM. Increased Cr oxidation kinetics was observed in the as-polished sample with the native oxide film resulting in formation of an overall more protective and compact film compared to the cathodically-activated sample. Electrochemical analysis showed that corrosive attack varied between dislocation cell boundaries to cell interiors depending on the initial surface state and polarization conditions. In conclusion, a corrosion mechanism is proposed to explain this variation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A thermodynamic analysis of hypervelocity impacts on metals

Hypervelocity impacts on solid surfaces generate and propagate shock waves and create extreme density, pressure, and temperature conditions. Understanding hypervelocity impacts requires a multi-physics approach at the intersection of fluid dynamics, solid mechanics, and quantum theory, which has the capability of spanning multiple different dominating physics regimes. These strong shock conditions found in hypervelocity impacts are studied in the context of high-speed collisions between spacecraft and meteoroids and debris. We study the problem of hypervelocity impacts on metals with the lens of understanding the processes that lead to the creation of an impact generated plasma. In order to characterize the impacting particle as well as its threat to spacecraft from the subsequent radiating plasma an analysis was performed that highlights how the kinetic energy of the impactor is deposited into a target, how the shock wave and shocked material in the target evolve, and how these effects can lead towards the creation of an impact plasma through heating. The analysis suggests that at low velocities the majority of the energy from a hypervelocity impact is stored in reversible elastic energy and does not contribute to strong heating. At higher impact velocities, nucleus and electron thermal components dominate, increasing the energy available for ionization. Overall, this work provides the complete set of state variable initial conditions for any study of hypervelocity impact plasma expansion and its effects.

42 ENGINEERING↗
DOE OSTI2020

Ultracentral Collisions of Small and Deformed Systems at RHIC

We study a range of collision systems involving deformed ions and compare the elliptic and triangular flow harmonics produced in a hydrodynamics scenario versus a color glass condensate (CGC) scenario. For the hydrodynamics scenario, we generate initial conditions using TRENTO and work within a linear response approximation to obtain the final flow harmonics. For the CGC scenario, we use the explicit calculation of two-gluon correlations taken in the high-p T “(semi)dilute-(semi)dilute” regime to express the flow harmonics in terms of the density profile of the collision. Here, we consider ultracentral collisions of deformed ions as a testbed for these comparisons because the difference between tip-on-tip and side-on-side collisions modifies the multiplicity dependence in both scenarios, even at zero impact parameter. We find significant qualitative differences in the multiplicity dependence obtained in the initial conditions+hydrodynamics scenario and the CGC scenario, allowing these collisions of deformed ions to be used as a powerful discriminator between models. We also find that sub-nucleonic fluctuations have a systematic effect on the elliptic and triangular flow harmonics which are most discriminating in 0–1% ultracentral symmetric collisions of small deformed ions and in 0–10% d197Au collisions. The collision systems we consider are 238 U 238 U, d 197 Au, 9 Be 197 Au, 9 Be 9 Be, 3 He 3 He, and 3 He 197 Au.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

A passive model for the evolution of subgrid-scale instabilities in turbulent flow regimes

In this work, a zero-dimensional model for the development of turbulent mixing layers driven by Rayleigh–Taylor (RT) or Richtmyer–Meshkov (RM) instabilities is presented, for application in Reynolds-Averaged Navier–Stokes (RANS) simulations. The model approximates the behavior of the turbulent mixing layer at a passively advected tracer particle that sits on the material interface of interest in the simulation, and is used to provide initial conditions for the RANS model once the mixing layer width has developed to a size that may be resolved on the computational mesh. The model relies on approximations to the Besnard–Harlow–Rauenzahn (BHR) RANS closure equations and is thus consistent the turbulent development of the mixing layer, unlike previous models of this type which have focused on the early-time laminar development of RT or RM instabilities. These instabilities will often grow through an initially laminar state before developing into turbulence, and criteria for switching from an existing laminar model to the proposed turbulent model are considered. The proposed model reproduces the mixing layer development of RT and RM instabilities seen in mixed width resolving BHR RANS simulations to a reasonable degree of accuracy, and relaxes the need for the RANS simulation to resolve the mixed width as soon as the mixing layer becomes turbulent, which may be impractical in cases where the mixing layer grows through a wide range of scales.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗