Engineering PapersSearch

SEARCH · Engineering Papers

Results for “lateral 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 19 records

Observations of Subduction, Downward Heat Flux and Dense Filament Collapse in the Northern Gulf of Mexico

Submesoscale processes are important contributors to the global heat budget and generally support upward heat transport through restratification. However, in salinity‐stratified regions, such as the northern Gulf of Mexico with its influx of freshwater from the Mississippi‐Atchafalaya river system, temperature can act like a passive tracer and submesoscale processes can contribute to downward heat transport. Oceanic heat content is a factor in many environmental risks the region faces, for example, hurricane intensification, and marine heatwaves. During the 2022 field campaign of the Submesoscales Under Near‐Resonant Inertial Shear Experiment, a sampling plan was developed to study such submesoscale processes in high resolution. Over 31 hr, four assets (two research ships and two remotely controlled boats) drove in parallel across a dense filament, capturing its evolution in time and space. The observations show that surface waters, warmed by daytime solar radiation, were subducted and that the associated overturning circulation transported heat below the surface layer where it was later irreversibly mixed away. The estimated downward heat flux was as strong as the concurrent net air‐sea heat flux into the ocean. The filament was then observed to rapidly collapse which we attribute to boundary layer turbulence and the breakdown of geostrophic balance. The collapsing fronts display behaviors indicative of gravity currents. These observations highlight how in salinity‐stratified regions, frontal dynamics can be associated with downward heat flux and how the submesoscale can play an important role in the oceanic heat budget.

58 GEOSCIENCES

Extended Heat Transfer Model Dataset

Data repository for the paper: Molerus and Wirth's Heat Transfer Model for Bubbling Fluidized Beds: Proposal for an Extended Model Including Immersed Tube Banks and Particle Cross-Flow Powder Technology 2025

contact time

Size-resolved Eddy-Covariance Particle Flux Measurement during the TRACER Campaign (Final Report)

The main goal of the TRacking Aerosol Convection interactions ExpeRiment (TRACER) campaign was to study aerosol–cloud interactions during deep convection over the Houston area. This project deployed a suite of instrumentation with the aim to (1) quantify turbulent vertical particle fluxes during at DOE-ARM sites, including TRACER, (2) assess hygroscopic growth factors and hygroscopicity parameters of the material driving modal aerosol growth during new particle formation and growth events, (3) derive turbulent aerosol mass fluxes using co-located Doppler LIDAR measurements, and (4) create quality-controlled PI data products to support future research utilizing data collected during the TRACER campaign. This report summarized the main findings from the deployments at two DOE-ARM sites. Briefly, we found that new particle formation may occur aloft, in a residual layer, near the top of the boundary layer. Small grown particles appear later due to downward mixing with daytime turbulence. The species that are responsible for aerosol modal growth had hygroscopicity parameters varying between 0.05 and 0.34. These values systematically depended on the wind sector, suggesting that the chemical composition of the precursors differed. This work demonstrated that lidar retrievals of the elastic backscatter and Doppler velocity can be used to obtain surface number emissions of particles with a diameter greater than 0.53 µm. During TRACER, emission particle number fluxes peaked near ∼ 100 cm−2 s−1. Multiple quality-controlled PI data products that will support future TRACER related science were generated and made publically available.

54 ENVIRONMENTAL SCIENCES

Five Year Comparison of Mixing Height Determinations at the Savannah River Site

Air quality dispersion modeling is performed for the Savannah River Site (SRS) to demonstrate compliance with applicable regulations. The AMS/EPA Regulatory Model (AERMOD) modeling system is an EPA recommended model for air quality applications with a data preprocessor (AERMET) to incorporate meteorological data collected on site. AERMET parameterizes or calculates meteorological variables that are not directly measured onsite. One of the parameters estimated by AERMET is the atmospheric mixing height. While the mixing height is not currently a measurement input into AERMET, SRS has the capability to measure the local mixing height. The Savannah River National Laboratory (SRNL) operates a Vaisala CL31 Lidar Ceilometer which estimates mixing height from aerosol backscatter. This study compares the parameterized mixing height from AERMET to the ceilometer estimated mixing height for the current regulatory period at SRS incorporating data from 2015-2019. Results from this study showed the average daily minimum values (morning) from AERMET were an order of magnitude lower than the commonly used Holzworth (1972) method and the ceilometer estimated mixing heights. Additionally, on average, the ceilometer exhibited a daily maximum mixing height value that occurred 1-3 hours later than the AERMET estimated maximum. This difference is likely due to the nighttime atmospheric mixing height assumptions and calculations used by AERMET. The AERMET algorithm cuts off mixing height growth at sunset while the ceilometer data show ongoing evening convection typical of the southeastern United States. These results suggest that the AERMET parametrization scheme assumptions may not be representative of a forested landscape and evening convection which could account for more mixing overnight. The results obtained in this study are significant for air dispersion modeling applications for regulatory purposes and worker safety. Mixing height can impact model estimated pollutant concentrations. A greater mixing height will provide more volume for pollutant dispersion. This report documents efforts to quantify the dependence of mixing height inputs toward a conservative estimated pollutant concentration.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W

Burning conditions and transportation pathways determine biomass-burning aerosol properties in the Ascension Island marine boundary layer

Abstract. African biomass-burning aerosol (BBA) in the southeast Atlantic Ocean (SEA) marine boundary layer (MBL) is an important contributor to Earth's radiation budget, yet its representation remains poorly constrained in regional and global climate models. Data from the Layered Atlantic Smoke Interactions with Clouds (LASIC) field campaign on Ascension Island (7.95° S, 14.36° W) provide insight into how burning conditions, fuel type, transport pathways, and atmospheric processing affect the chemical, microphysical, and optical properties of BBA between June and September 2017. A total of 10 individual plume events characterize the seasonal evolution of the BBA properties. Early-season inefficient fires, determined by low refractory black carbon to above-background carbon monoxide mixing ratios (rBC : ΔCO), led to enhanced concentrations of organic- and sulfate-rich aerosols. Mid-season efficient fires, determined by higher rBC : ΔCO values, led to rBC-enriched BBA. A mix of efficient and inefficient fires later in the season resulted in conflicting BBA properties. Prolonged transport (∼ 10 d) through the MBL and lower free troposphere (FT) facilitated chemical and aqueous-phase processing, which led to a reduction in organic aerosol mass concentrations. This resulted in lower organic aerosol (OA) to rBC (OA : rBC) mass ratios (2–5) in the MBL compared to higher values (5–15) in the nearby FT. These atmospheric and cloud oxidation processes yield more light-absorbing BBA and explain the notably low single-scattering albedo at 530 nm (SSA530) values (< 0.80) observed in the MBL. This study establishes a robust correlation between SSA530 and OA : rBC across the MBL and FT, underscoring the dependency of optical properties on chemical composition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Shock-driven three-fluid mixing with various chevron interface configurations

When a shock wave crosses a density interface, the Richtmyer–Meshkov instability causes perturbations to grow. Richtmyer–Meshkov instabilities arise from the deposition of vorticity from the misaligned density and pressure gradients at the shock front. In many engineering applications, microscopic surface roughness will grow into multi-mode perturbations, inducing mixing between the fluid on either side of an initial interface. Applications often have multiple interfaces, some of which are close enough to interact in the later stages of instability growth. In this study, we numerically investigate the mixing of a three-layer system with periodic zigzag (or chevron) interfaces, calculating the dependence of the width and mass of mixed material on properties such as the shock timing, chevron amplitude, multi-mode perturbation spectrum, density ratio, and shock mach number. The multi-mode case is also compared with a single-mode perturbation. The Flash hydrodynamic code is used to solve the Euler equations in three dimensions with adaptive grid refinement. Key results include a significant increase in mixed mass when changing from a single-mode to a multi-mode perturbation on one of the interfaces. The mixed width is mainly sensitive to the density ratio and chevron amplitude, whereas the mixed mass also depends on the multi-mode spectrum. In conclusion, steeper initial perturbation spectra have lower mixed mass at early times but a greater mixed mass after the reflected shock transits back across the layer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC

Detect the Unobservable: Abnormality Detection in mixed Autonomy for Lane Change Maneuver with Following Vehicles’ Trajectories Only

Highly Automated Vehicles (HAVs) and Advanced Driver-Assistance Systems (ADAS) are transforming modern transportation with enhanced mobility, safety, and efficiency. Despite their advantages, cybersecurity vulnerabilities in these systems can lead to abnormal behavior, posing significant risks to surrounding human-driven vehicles (HDVs) in mixed traffic environments. Here, this article addresses the challenge of detecting abnormal lateral movements of HAVs/ADAS vehicles using only trajectory profiles of following HDVs. Specifically, we propose a novel modeling approach that captures both normal and abnormal lateral behaviors through vehicle kinematics, integrated decision-making processes, vehicle control using symbolic regression for lane change vehicles. Additionally, we introduce an abnormality detection framework that relies on observable HDV data, even in occlusion scenarios. The framework evaluates the sensitivity of various car-following models to detect abnormal behaviors, providing insights into the interaction between HAVs/ADAS and HDVs in mixed autonomy systems.

Connected and Automated vehicles

Boomerang mechanism explaining the excess radio background

We propose a boomerang mechanism for the explanation of the excess radio background detected by ARCADE 2. In an early stage of the Universe, at a temperature 𝑇 in the range ∼ 0.1 keV−⁢1 MeV, a fraction of relic neutrinos is resonantly converted into dark neutrinos by mixing induced by a preexisting lepton asymmetry. Dark neutrinos decay much later into a dark-standard photon state and a dark fermion, with a lifetime longer than the age of the Universe, as required by a solution to the excess radio background. This scenario circumvents the upper bound on the neutrino magnetic moment but still implies a testable lower bound.

Dev, P. S. Bhupal [Washington Univ., St. Louis, MO

Aerosol Particle Mixing State and Composition at AMF3 (Interim Field Campaign Report)

The campaign is designed to support deployment of aerosol sampling equipment at the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s third ARM Mobile Facility (AMF3) to enable both singe-particle and bulk characterization of the aerosol particles at the field site in Bankhead National Forest, northern Alabama. We are proposing to do this work as part of our proposal “Generalizing aerosol mixing state: synthesis from observations and connection to models,” which has the pre-proposal tracking number PRE-0000035878. The full proposal has Dr. Rachel O’Brien as the Principal Investigator (PI) and Dr. Andy Ault and Dr. Nicole Riemer as Co-PIs. We sought to deploy instruments in winter 2025 and summer 2025 to obtain samples in different seasons, but were delayed by AMF3 set-up and availability of students and PIs. We just completed our first deployment this summer (2025) and we are interested in obtaining samples during this upcoming fall and winter to see the seasonal trends in aerosol mixing state at the site. This will provide a comparison of the impact of biogenic emissions on the mixing state and the chemical composition of the aerosol particles. Overall, the proposed study will provide key information regarding the mixing state of aerosol particles at AMF3 with direct ties to impacts of the particles on cloud formation and light scattering/absorption. From our deployment that ended las week we have some visible observations from our samples. We are observing a regular amount of black carbon on stage 6 of our Micro-Orifice Uniform Deposit Impactor (MOUDI) and we have also observed some short, large increases in aerosol loading. We hope to continue observations of these two features to determine the sources/aging for the black carbon and the physical mixing state as well as the sources for the short (a few minutes) increases in the aerosol loading and determine if the source for these is dust or sea salt or another source. We are also interested in comparing the black carbon we see here to samples collected during the agriculture burns in January/February and in comparing the larger particles to bioaerosols in the spring. To obtain these comparisons, we will need to collect more samples in the other seasons. Finally, we have some size-resolved particle collectors that can be deployed at the same time as the tethered balloon launches. We have discussed with Swarup China about trying to get the timing to work in August and we are interested in coordinating again later, if possible. This will give us ground-level and vertical aerosol mixing state sample sets that will be very exciting to analyze.

54 ENVIRONMENTAL SCIENCES

Reduced order modeling of a fluidized bed particle receiver for concentrating solar power with thermal energy storage

Oxide particles can serve as both the heat transfer and thermal energy storage (TES) media for next-generation concentrating solar power (CSP) plants where high-temperature TES enables dispatchable electricity from efficient power cycles with firing temperatures above 600 °C. Transferring heat to flowing particles at such high temperatures in a MW-scale central tower receiver remains a challenge for the CSP community. For indirect receivers with external walls to contain the particles, maintaining wall temperatures below the limits of structural metal alloys requires high heat transfer coefficients between the wall and the moving particle stream. Bubbling fluidization of downward-flowing particles can sustain high bed-wall heat transfer coefficients (> 1000 W m -2 K -1 ). Using experimentally calibrated correlations for bed-wall heat transfer and vertical particle dispersion, this study implements an axially discretized zonal model of a counterflow fluidized bed receiver to explore how bubbling fluidization may enable indirect cavity particle receivers. High bed-wall heat transfer coefficients support solar fluxes on angled cavity walls > 200 kW m -2 at peak aperture fluxes of 980 kW m -2 while maintaining external wall temperatures < 950 °C. Lateral particle dispersion enables hotter particles near the receiver leading edge to mix with cooler particles further from the leading edge to lower maximum external wall temperatures. Parametric studies identify how mass fluxes, particle dispersion, and solar concentrations impact indirect receiver thermal efficiency and uniformity for a CSP plant. These studies provide a basis for the design of indirect fluidized-bed cavity receivers that can maintain particle outlet temperatures for TES above 750 °C.

14 SOLAR ENERGY

Multi‐Material Gradient Printing Using Meniscus‐enabled Projection Stereolithography (MAPS)

Light‐based additive manufacturing methods are widely used to print high‐resolution 3D structures for applications in tissue engineering, soft robotics, photonics, and microfluidics, among others. Despite this progress, multi‐material printing with these methods remains challenging due to constraints associated with hardware modifications, control systems, cross‐contamination, waste, and resin properties. Here, a new printing platform coined Meniscus‐enabled Projection Stereolithography (MAPS) is reported, a vat‐free method that relies on generating and maintaining a resin meniscus between a crosslinked structure and bottom window to print lateral, vertical, discrete, or gradient multi‐material 3D structures with no waste and user‐defined mixing between layers. MAPS is compatible with a wide range of resins shown and can print complex multi‐material 3D structures without requiring specialized hardware, software, or complex washing protocols. MAPS's ability to print structures with microscale variations in mechanical stiffness, opacity, surface energy, cell densities, and magnetic properties provides a generic method to make advanced materials for a broad range of applications.

bioprinting

Effects of 5-ion 6-beam sequential irradiation in the presence and absence of hindlimb or control hindlimb unloading on behavioral performances and plasma metabolic pathways of Fischer 344 rats

Introduction Effects and interactions between different spaceflight stressors are expected to be experienced by crew on missions when exposed to microgravity and galactic cosmic rays (GCRs). One of the limitations of previous studies on simulated weightlessness using hindlimb unloading (HU) is that a control HU condition was not included. Methods We characterized the behavioral performance of male Fischer rats 2 months after sham or total body irradiation with a simplified 5-ion 6-mixed-beam exposure representative of GCRs in the absence or presence of HU. Six months later, the plasma, hippocampus, and cortex were processed to determine whether the behavioral effects were associated with long-term alterations in the metabolic pathways. Results In the open field without and with objects, interactions were observed for radiation × HU. In the plasma of animals that were not under the HU or control HU condition, the riboflavin metabolic pathway was affected most for sham irradiation vs. 0.75 Gy exposure. Analysis of the effects of control HU on plasma in the sham-irradiated animals showed that the alanine, aspartate, glutamate, riboflavin, and glutamine metabolisms as well as arginine biosynthesis were affected. The effects of control HU on the hippocampus in the sham-irradiated animals showed that the phenylalanine, tyrosine, and tryptophan pathway was affected the most. Analysis of effects of 0.75 Gy irradiation on the cortex of control HU animals showed that the glutamine and glutamate metabolic pathway was affected similar to the hippocampus, while the riboflavin pathway was affected in animals that were not under the control HU condition. The effects of control HU on the cortex in sham-irradiated animals showed that the riboflavin metabolic pathway was affected. Animals receiving 0.75 Gy of irradiation showed impaired glutamine and glutamate metabolic pathway, whereas animals receiving 1.5 Gy of irradiation showed impaired riboflavin metabolic pathways. A total of 21 plasma metabolites were correlated with the behavioral measures, indicating that plasma and brain biomarkers associated with behavioral performance are dependent on the environmental conditions experienced. Discussion Phenylalanine, tyrosine, and tryptophan metabolism as well as phenylalanine and tryptophan as plasma metabolites are biomarkers that can be considered for spaceflight as they were revealed in both Fischer and WAG/Rij rats exposed to simGCRsim and/or HU.

Physiology

Effects of external magnetic fields on the multi-mode Rayleigh–Taylor instability

The magneto-Rayleigh–Taylor instability (mRTI) is a key process in inertial confinement fusion and is thought to be widespread in the interstellar medium, where it can concentrate plasma into discrete structures. We present resistive MHD simulations of the nonlinear evolution of multi-mode mRTI in both two and three dimensions, examining the effects of uniform external magnetic fields oriented either parallel or perpendicular to the initial interface. In both 2D and 3D, weak parallel fields enhance mixing-zone growth, whereas stronger fields suppress it. For perpendicular fields, growth is initially inhibited but becomes enhanced at later times. These behaviors arise from magnetic tension, which modifies flow anisotropy, buoyancy, drag, and vortex dynamics. The interplay of these mechanisms governs the distinct ways in which magnetic fields influence mRTI evolution.

Flow instabilities

Effect of layer bending on montmorillonite hydration and structure from molecular simulation

Conceptual models of smectite hydration include planar (flat) clay layers that undergo stepwise expansion as successive monolayers of water molecules fill the interlayer regions. However, X-ray diffraction (XRD) studies indicate the presence of interstratified hydration states, suggesting non-uniform interlayer hydration in smectites. Additionally, recent theoretical studies have shown that clay layers can adopt bent configurations over nanometer-scale lateral dimensions with minimal effect on mechanical properties. Therefore, in this study we used molecular simulations to evaluate structural properties and water adsorption isotherms for montmorillonite models composed of bent clay layers in mixed hydration states. Results are compared with models consisting of planar clay layers with interstratified hydration states (e.g. 1W–2W). The small degree of bending in these models (up to 1.5 Å of vertical displacement over a 1.3 nm lateral dimension) had little or no effect on bond lengths and angle distributions within the clay layers. Except for models that included dry states, porosities and simulated water adsorption isotherms were nearly identical for bent or flat clay layers with the same averaged layer spacing. Similar agreement was seen with Na- and Ca-exchanged clays. In conclusion, while the small bent models did not retain their configurations during unconstrained molecular dynamics simulation with flexible clay layers, we show that bent structures are stable at much larger length scales by simulating a 41.6×7.1 nm 2 system that included dehydrated and hydrated regions in the same interlayer.

58 GEOSCIENCES

A collective trigger for widespread planetesimal formation revealed by accretion ages

The formation of planetesimals was an integral part of the cascading series of processes that built the terrestrial planets. To illuminate planetesimal formation, here we develop a refined thermal evolution model to calculate the formation ages of meteorite parent planetesimals. This model includes chemical reactions and phase changes during heating, as well as natural variations in the proportions of the constituent phases of these planetesimals. We find that the parent bodies of non-carbonaceous (NC) and carbonaceous (CC) iron meteorites start forming at very similar times ( ∼ 0.95 Myr after calcium-aluminium-rich inclusion [CAI] formation) and occupy overlapping time windows. NC and CC chondrite parent bodies formed later during non-overlapping periods. We combine these ages with proportions of isotopic end-members we recover from mixing models to construct records of motion throughout the protoplanetary disk. These records argue that NC and CC material traversed the barrier in the disk after ∼ 0.95 Myr after CAI formation. The onset of this motion coincided with planetesimal formation, indicating that the phenomenon that drove motion also triggered planetesimal formation. We argue that this feature also served as the semi-permeable barrier in the disk. Although its identity is uncertain, the effects this phenomenon had on the timing of planetesimal formation and motion through the disk can now serve as constraints on models of disk evolution. Models that reproduce these effects would elucidate the nature and implications of this phenomenon, which is key to unlocking a holistic model of terrestrial planet building.

Astronomy and AstroPhysics

Tuning Anisotropic Optical Properties of Inorganic and Hybrid Organic–Inorganic MXenes via Topochemical Surface Modification

Surface groups are central to the properties of MXenes, yet their role in optical anisotropy remains largely unexplored. Here, we use a topochemical route to synthesize single crystals of stacked Ti 3 C 2 Cl 2 and hybrid organic–inorganic MXenes (h-MXenes) with lateral sizes of 38–75 μm, rotational registry, and tunable interlayer spacing. Solid-state NMR spectroscopy shows that topochemical substitution generates mixed amido, imido, and hydride surface motifs, which modify the electronic structure of the Ti 3 C 2 inorganic core. Imaging spectroscopic ellipsometry with micron-scale spatial resolution enables reconstruction of the complex dielectric tensor of individual multilayer crystals. Ti 3 C 2 Cl 2 exhibits a type-II hyperbolicity above 930 nm, whereas h-MXenes do not display hyperbolicity within the measured 300–1700 nm window, instead showing reduced in-plane conductivity, suppressed out-of-plane light absorption, and a chain-length-dependent blue shift of a near-infrared absorption feature. These results demonstrate topochemical surface modification as a direct handle for engineering MXenes as surface-programmable optical media.

Hybrid materials

The new ITER baseline, research plan and open R&D issues

A new baseline (NB) has been proposed by the ITER Project to ensure a robust achievement of the Projects’ goals, in view of past challenges including delays incurred due to the Covid-19 pandemic, technical challenges in completing first-of-a-kind components and in nuclear licensing. The NB includes modifications to the configuration of the ITER device and its ancillaries (e.g. change from beryllium to tungsten as first wall material, modification of the heating and current drive mix, etc.) as well as additional testing of components (e.g. toroidal field coils) or phased installation (start with inertially cooled first wall before later installation of the final actively water-cooled components) to minimise operational risks. In the NB, the ITER research plan (IRP) will be divided into three main phases: (a) start of research operation, with 40 MW of ECH and 10 MW of ICH, which will focus on the demonstration of 15 MA operation in L-mode, commissioning of all required systems, including disruption mitigation, and the demonstration of H-mode plasma operation in deuterium; (b) DT-1, with 60–67 MW of ECH, 33 MW of neutral beam injection (NBI) and 10–20 MW of ICH, which will demonstrate robust operation in high confinement H-mode plasmas in DT up to Q ⩾ 10 and for burn durations of 300–500 s within an accumulated neutron fluence of ∼1% of the ITER machine’s lifetime total, and; (c) DT-2, with up to 67 MW of ECH, up to 49.5 MW of NBI and up to 20 MW of ICH, with the ITER tokamak and ancillaries in their final configuration to demonstrate routine operation in DT plasmas at high Q and the Q ⩾ 5 long-pulse and steady-state scenarios to the final neutron fluence and to perform R&D on nuclear fusion reactor issues. The logic, physics basis, modelling and experimental evaluations carried out to support the NB and the associated IRP are described. These include the impact of the tungsten wall on plasma scenarios and associated risk mitigation measures, as well as the optimisation of the tokamak components and ancillaries to minimise Project risks. Open R&D issues related to these evaluations and mitigation measures are also described together with experimental, modelling and validation activities required to address them.

ITER

Multidimensional Modeling of Mixture Formation in a Hydrogen-Fueled Heavy-Duty Optical Engine With Direct Injection

Hydrogen (H 2 ), as a carbon-free fuel, is considered as one of the most promising solutions to reduce the carbon footprint of hard-to-decarbonize energy and transportation sectors. As such, hydrogen-fueled internal combustion engines (H 2 ICEs) have recently been receiving increasing attention, particularly in applications such as on-road/off-road heavy-duty transport and combined heat and power. The direct injection (DI) of gaseous hydrogen into the combustion chamber offers great potential for achieving high power density and high engine efficiency, while mitigating the risk of backfire and reducing pre-ignition. However, the numerical simulation of H 2 DI system remains a formidable challenge associated with the high computational cost of reproducing compressible supersonic flow and shocks in narrow injector passages and in near-nozzle regions. In general, there is a lack of well-established and validated practices for the modeling of high-pressure H 2 DI in large-bore engines. Here, to this end, this study focuses on computational fluid dynamics (CFD) modeling of the mixture formation process in a heavy-duty optical engine employing a medium-pressure H 2 DI system. Both large eddy simulations (LES) and Reynolds Averaged Navier–Stokes (RANS) simulations are performed and evaluated against optical data. Gaseous hydrogen is injected into the combustion chamber via a centrally located outward opening hollow-cone injector at a pressure of 40 bar. Simulations are carried out for two injection timings, namely, −120 and −60 °CA. The numerical predictions for H 2 distribution in different horizontal and vertical planes during the compression stroke are systematically compared against optical data obtained through planar laser-induced fluorescence (PLIF) measurements. Overall, the LES approach using the Dynamic Structure model is found to have good predictive capabilities for the early jet penetration in terms of length and shape, as well as the later H 2 distributions. However, the unsteady RANS approach with the renormalization group $k - ϵ$ model, which is widely used by industry to model heavy-duty ICEs, significantly underpredicts the H 2 mixing, even at similar mesh resolution to that used in LES. These results indicate that there is a need for the improvement of mixing submodels within the RANS approach when applied to H 2 DI simulations.

LES