Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mixture modeling”

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 1,279 records · Page 71

IC Report for project “Developing nonlinear laser-plasma instability models for high-fidelity, multi-physics simulation capability for ICF/HED” [Slides]

We’ve achieved our technical goals: (1) Perform large-scale kinetic simulations of SRS and CBET a range of plasma and laser conditions; (2) Develop a new $\delta f$-Gaussian-mixture algorithm to represent non-Maxwellian distribution functions from particle trapping and time-dependent plasma response; (3) Develop physics-based preliminary nonlinear LPI models for coupling to LRT/ Mazinisin and AMP; and (4) Begin implementation of LPI sub-grid models into Mazinisin.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High altitude chemically reacting gas particle mixtures. Volume 1: A theoretical analysis and development of the numerical solution

The overall contractual effort and the theory and numerical solution for the Reacting and Multi-Phase (RAMP2) computer code are described. The code can be used to model the dominant phenomena which affect the prediction of liquid and solid rocket nozzle and orbital plume flow fields. Fundamental equations for steady flow of reacting gas-particle mixtures, method of characteristics, mesh point construction, and numerical integration of the conservation equations are considered herein.

Smith, S. D.↗

On the use of a global hydrogen-air combustion model in the calculation of turbulent reacting flows

A two-step global model describing the combustion of hydrogen in air at one atmosphere pressure is developed by comparing the temperature histories obtained from a 28-reaction H-O mechanism. Using criteria discussed herein, good agreement is obtained in the range of initial mixture temperatures from 1000K to 2000K and in the range of equivalence ratios from 0.2 to 2.0. The two-step global model is compared with the results obtained using an eight-step, H-O reaction mechanism in a computer program describing the turbulent diffusion of hydrogen in supersonic, axisymmetric and two-dimensional reacting flows. Comparisons of profiles of temperature, fuel mass flux, and pitot pressure are presented. The global model is judged to be adequate in flows which are not dominated by long ignition delay times.

Rogers, R. C.↗

Directional and polarimetric effects in the optical domain

Papers focused on land surface, atmospheric, and ocean properties are reported. Specific comments pertaining to polarization, models and inversion, and measurements, are given. Recommendations are: continued research into the application potential of the BRDF (Bidirectional Reflectance Distribution Function) and polarization properties of ground surface and atmospheric targets; three dimensional models, which account for the statistical behavior of remotely sensed data, should be extended and inverted in order to support analysis of data potentially covering rolling terrain such that pixels represent heterogeneous mixtures of surface cover types and project ground footprints with sizes between 10 to 6 km, the ground pixel sizes of planned future sensors; available reflectance models should be further validated by means of multi dimensional (directional, spectral, temporal) field data and existing models should be intercompared in more depth to evaluate their performance and limitations; existing methods for model inversion should be validated in more depth in order to quantify the practical limitations and the expected accuracy of the parameters retrieved and new approaches should be developed based upon apriori knowledge of plant canopy development and spectral BRDF properties; there is a need to establish a protocol of validation and intercomparison of the indices and compositing techniques which have been proposed during these last years.

Vanderbilt, V. C.↗

Impact of recent ENDF nuclear data update, high initial enrichment and high burnup fuel on critical experiments applicability determination via the integral index c k for burnup credit validation

In 2012, NUREG/CR-7109 reported on the validation of burnup credit calculations involving major and minor actinides and major fission products which was investigated for pressurized and boiling water reactor (PWR and BWR) fuel enrichments up to 5 wt% 235 U and assembly-average burnups up to 60 GWd/MTU. Recently, there has been interest in increasing the maximum enrichment used in PWR fuel as high as 8 wt% 235 U and correspondingly increasing the maximum assembly-average burnups to approximately 75 GWd/MTU. These proposed increases in enrichment and burnup necessitate reinvestigation of the validation basis for k eff calculations for this expanded application space. Additionally, the 2012 study was performed by using the Evaluated Nuclear Data File (ENDF)/B-VII.0 nuclear data with the SCALE 6 covariance library, and the effects of using the newly released ENDF/B-VII.1 and ENDF/B-VIII.0 nuclear data and covariance libraries should be evaluated. In this work, published in NUREG/CR-7309 in 2025, the validation assessment was performed consistently with NUREG/CR-7109: modeling irradiated fuel assemblies in the Generic Burnup Credit (GBC)-32 cask defined in NUREG/CR-6747. The TSUNAMI-3D sequence was used to generate sensitivity data for the application model, and the data were compared with sensitivity data from select benchmark models. The integral parameter c k is the metric of similarity used in this study and is consistent with NUREG/CR-7109, where a c k value in excess of 0.8 indicates sufficient similarity for use in validation. A new set of benchmark experiments with sensitivity data has been assembled for this effort. The number of experiments with available sensitivity data is now 2,104, compared to 474 in NUREG/CR-7109. This increase was facilitated by the efforts of the Nuclear Energy Agency to generate sensitivity data for a majority of the experiments in the International Criticality Safety Benchmark Evaluation Project (ICSBEP) Handbook to supplement the data available in the Oak Ridge National Laboratory (ORNL) Verified, Archived Library of Inputs and Data (VALID). The complete set of benchmarks considered here includes experiments for low-enriched uranium (LEU), intermediate enriched uranium (IEU), and a mixture of uranium and plutonium (MIX) from the ICSBEP Handbook and VALID, as well as ORNL models of the Haut Taux de Combustion (HTC) experiments and other potentially relevant models not included in VALID. The updated similarity study shows that none of the extended burnup and higher enrichment combinations considered show a significant decrease in the number of potentially applicable experiments, meaning sufficient critical experiments exist for the validation of BUC criticality safety calculations, with initial enrichments up to 8 wt% 235 U and burnups up to 80 GWd/MTU. Additionally, both the ENDF/B-VII.1 and ENDF/B-VIII.0 nuclear data libraries can be used for validation since the number of critical experiments applicable for validation increases for most cases with the most recent nuclear data compared to the previous one. As in previous BUC validation studies, the French HTC experiments are the most similar in a majority of the application cases studied, especially from representative discharge burnups ranging from 40 to 80 GWd/MTU. In conclusion, these results match the conclusions presented in NUREG/CR-7109 regarding validation of the primary actinides in BUC analyses.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Tar Production from Biomass Pyrolysis in a Fluidized Bed Reactor: A Novel Turbulent Multiphase Flow Formulation

A novel multiphase flow model is presented for describing the pyrolysis of biomass in a 'bubbling' fluidized bed reactor. The mixture of biomass and sand in a gaseous flow is conceptualized as a particulate phase composed of two classes interacting with the carrier gaseous flow. The solid biomass is composed of three initial species: cellulose, hemicellulose and lignin. From each of these initial species, two new solid species originate during pyrolysis: an 'active' species and a char, thus totaling seven solid-biomass species. The gas phase is composed of the original carrier gas (steam), tar and gas; the last two species originate from the volumetric pyrolysis reaction. The conservation equations are derived from the Boltzmann equations through ensemble averaging. Stresses in the gaseous phase are the sum of the Newtonian and Reynolds (turbulent) contributions. The particulate phase stresses are the sum of collisional and Reynolds contributions. Heat transfer between phases, and heat transfer between classes in the particulate phase is modeled, the last resulting from collisions between sand and biomass. Closure of the equations must be performed by modeling the Reynolds stresses for both phases. The results of a simplified version (first step) of the model are presented.

Bellan, J.↗

Development of a Machine-Learned Cruise Guide Indicator for Rotorcraft

This paper presents a machine-learned virtual cruise guide indicator (vCGI) for Chinook helicopters. Two temporal neural networks were trained and evaluated on measured data from 55 flight tests, one for the fore rotor and another for the aft rotor, to predict a vCGI value, which protects 23 components from fatigue damage during steady-state conditions. Three different classes of machine learning architectures were evaluated for prediction of the vCGI from time sequences: a temporal convolutional neural network with 1D dilated causal convolutions, a long short-term memory recurrent neural network, and an attention-based transformer architecture. The final average model accuracy on unseen flight data is currently greater than 93% for CGI values which could result in fatigue damage and 90% for normal operation CGI values. Model accuracy was improved through a series of advancements in:(1) selection of optimal training data using temporal collective variables and unsupervised learning, (2) dataset augmentation with maximum-entropy temporal collective variables, and (3) implementation of a mixture-of-experts classification- regression approach using an adversarial classification approach to assign maneuver labels. The results are presented for each advancement in model development along with lessons learned in training machine learning models on real- world, time-dependent rotorcraft data.

Boyer, Mathew↗

Control of core–shell nanoparticles properties through plasma synthesis: a computational study

The improved properties of core–shell nanoparticles (CSNPs) over homogeneous nanoparticles (NPs) have expanded and diversified the applications of these nanomaterials. However, controlling the properties of CSNPs can be a challenging task. Low temperature plasmas have proven to be an effective method of producing NPs with uniform size and morphology, and high yield. That said, NP transport and growth dynamics are sensitive to LTP properties. We report on a computational investigation of the evolution of Ge–Si CSNP properties as a function of operating conditions through the modeling of a flowing, two-zone inductively coupled plasma (ICP) reactor. Ar/GeH 4 and Ar/SiH 4 gas mixtures were supplied to separate plasma zones at a pressure of 1 Torr to promote growth of Ge cores and Si shells. The negatively charged CSNPs are trapped electrostatically in the vicinity of the antennas where the plasma is generated and where the majority of particle growth occurs. Particles that grow to a critical size are then de-trapped by fluid drag due to neutral gas flow. A two-dimensional hybrid plasma model coupled with a three-dimensional kinetic NP transport model were utilized to resolve plasma chemistry and NP growth processes that take place on distinct timescales. The trends in CSNP properties and trapping mechanisms associated with flow rate, applied ICP power and inlet precursor fraction are discussed. While the spatial distribution of plasma produced radical species can have significant impact on the NP growth process, the NP transport dynamics are what ultimately dictates the growth environment that is unique to each particle and so determines their final dimension and composition. The key to optimizing reactor conditions involves controlling the spatial density of growth species and plasma profile as a means to tailor particle trapping dynamics suitable to produce CSNPs for a specific application.

36 MATERIALS SCIENCE↗

Continued Testing and Modeling of the AirSquared All-Metal Tritium Scroll Pump

The implementation of fusion energy requires the processing of the deuterium-tritium mixture used to fuel the reaction. Many gas processing operations require either a gas-motive force or a pressure gradient across a system; both are achieved using vacuum pumps. However, because of the material challenges of working with tritium, specialized pumps are required that use only metal on wetted surfaces. This prevents tritium retention and degradation of materials such as oils and polymers and reduces the generation of tritiated waste. Typical gas handling systems for tritium operations involve metal bellows pumps, all-metal scroll pumps, and metal turbomolecular pumps, depending on the pressure requirements. With the industry standard Normetex Type 15 m 3 /h (typically referred to as the Normetex 15) pump no longer available, the fusion market requires viable alternatives. In this study, we test an American-made pump, the AirSquared V16H034A-C01, and compare its performance to the Normetex. A performance model of the AirSquared pump is also produced, which can be applied to general process models to generate an accurate description of pressures and flows through various unit operations.

08 HYDROGEN↗

Fluid hydrogen at high density - Pressure dissociation

A model for the Helmholtz free energy of fluid hydrogen at high density and high temperature is developed. This model aims at describing both pressure and temperature dissociation and ionization and bears directly on equations of state of partially ionized plasmas, as encountered in astrophysical situations and high-pressure experiments. This paper focuses on a mixture of hydrogen atoms and molecules and is devoted to the study of the phenomenon of pressure dissociation at finite temperatures. In the present model, the strong interactions are described with realistic potentials and are computed with a modified Weeks-Chandler-Andersen fluid perturbation theory that reproduces Monte Carlo simulations to better than 3 percent. Theoretical Hugoniot curves derived from the model are in excellent agreement with experimental data.

Saumon, Didier↗

Micro Acoustic Resonant Chambers for Heating/Agitating/Mixing (MARCHAM)

A variety of applications require the mixing and/or heating of a slurry made from a powder/fluid mixture. One of these applications, Sub Critical Water Extraction (SCWE), is a process where water and an environmental powder sample (sieved soil, drill cuttings, etc.) are heated in a sealed chamber to temperatures greater than 200 degrees Celsius by allowing the pressure to increase, but without reaching the critical point of water. At these temperatures, the ability of water to extract organics from solid particulate increases drastically. This paper describes the modeling and experimentation on the use of an acoustic resonant chamber which is part of an amino acid detection instrument called Astrobionibbler [Noell et al. 2014, 2015]. In this instrument we use acoustics to excite a fluid- solid fines mixture in different frequency/Amplitude regimes to accomplish a variety of sample processing tasks. Driving the acoustic resonant chamber at lower frequencies can create circulation patterns in the fluid and mixes the liquid and fines, while driving the chamber at higher frequencies one can agitate the fluid and powder and create a suspension. If one then drives the chamber at high amplitude at resonance heating of the slurry occurs. In the mixing and agitating cell the particle levitation force depends on the relative densities and compressibility’s of the particulate and fluid and on the kinetic and potential energy densities associated with the velocity and pressure fields [Glynne-Jones, Boltryk and Hill 2012] in the cell. When heating, the piezoelectric transducer and chamber is driven at high power in resonance where the solid/fines region is modelled as acoustic transmission line with a large loss component. In this regime, heat is pumped into the solution/fines mixture and rapidly heats the sample. We have modeled the piezoelectric transducer/chamber/ sample using Mason’s equivalent circuit. In order to assess the validity of the model we have built and tested a variety of chambers. This paper describes the experimental results which are in general agreement with theory within the limitations of the modeling.

fluid particle interaction↗

Tri-Gas Pressurization System Testing and Modeling for Cryogenic Applications

The use of Tri‐gas in rocket propulsion systems is somewhat of a new technology. This paper defines Tri‐gas as a mixture of gases composed largely of helium with a small percentage of a stoichiometric mixture of hydrogen and oxygen. When exposed to a catalyst the hydrogen and oxygen in the mixture combusts, significantly raising the temperature of the mixture. The increase in enthalpy resulting from the combustion process significantly decreases the required quantity of gas needed to pressurize the ullage of the vehicle propellant tanks. The objective of this effort was to better understand the operating characteristics of Tri‐gas in a pressurization system with low temperature applications. In conjunction with ongoing programs at NASA Marshall Space Flight Center, an effort has been undertaken to evaluate the operating characteristics of Tri‐gas through modeling and bench testing. Through improved understanding of the operating characteristics, the risk of using this new technology in a launch vehicle propulsion system was reduced. Bench testing of Tri‐gas was a multistep process that targeted gas characteristics and performance aspects that pose a risk to application in a pressurization system. Pressurization systems are vital to propulsion system performance. Keeping a target ullage pressure in propulsions tanks is necessary to supply propellant at the conditions and flow rates required to maintain desired engine functionality. The first component of testing consisted of sampling Tri‐gas sources that had been stagnant for various lengths of time in order to determine the rate at which stratification takes place. Second, a bench test was set up in which Tri‐gas was sent through a catalyst bed. This test was designed to evaluate the performance characteristics of Tri‐gas, under low temperature inlet temperatures, in a flight‐like catalyst bed reactor. The third, most complex, test examined the performance characteristics of Tri‐gas at low temperature temperatures in a test configuration built to more closely resemble a vehicle pressurization system. The results of these bench tests address various risks that all relate to underperformance of Tri‐gas pressurization systems.

Taylor, B.↗

Large Primtive Asteroids: Thermal and Dynamical Context

Primitiveasteroids, most of which are located in the outer belt and Jupiter’s Trojan clouds, provide information related to the origin and evolution of the solar system and the conditions in which the solar nebula was formed. These asteroids are widely thought to be the origin of the least-altered car-bonaceous chondrite that allow us to put crucial con-straints on the current dynamical and thermal theories of the formation and evolution of the early solar sys-tem. The nature of surface composition of large and low-albedo asteroids, like (1) Ceres, (10) Hygiea and (52) Europa, is still under intense debate and different interpretations have beenput forth to explain the ab-sorption features in these objects (e.g.,[1, 2, 3,4,5,6]). Laboratory and spectroscopic experiments on me-teorites that represent all nine carbonaceous chondrite types also found no spectral matches for these large asteroids (e.g.,[7]). Previous studies of asteroid Ceres (the largest primitive asteroid in the solar system) have been conducted to constrain and estimate its surface composition (e.g., [2, 5]). Using linear mixing,[2]found hydroxide brucite, serpentines, and carbonates, to be consistent with Ceres’ ground-based spectra.[5]estimated the surface composition of Ceres and found evidence of widespread NH3-phyllosilicates across its surface using best-fit solutions to Dawn’s NIR spectra. The presence of NH3-phyllosilicates implies that mate-rial from the outer solar system was incorporated into large primitive asteroids, either during their formation at great heliocentric distance or by incorporation of material transported into the Main Belt region. Here we presentnew largeand primitive asteroidsthat share the same spectral similaritieswith the largest asteroid in the solar system, Ceres.We also present the context of these new observations in terms oftheirthermal and dynamicalevolution.Large Primitive Asteroids: Dynamical and Thermal Context:Orbital Distribution of LargePrimitiveAsteroids: Constraining the mineralogy and surface composition of large primitive asteroids will place crucial con-straints on existing dynamical and thermal theories of the formation of the early solar system. We identified several additional asteroids in the Europa-like group in addition to asteroids Europa, Euphrosyne, and Patien-tia, that wereidentified in [8]using NASA IRTFtele-scope (e.g.,Figure 1). These new asteroids along with the already-observed Ceres-and Europa-like group members are localized in the 2.8 < a < 3.4 AU region and characterized by larger sizes, showing an interest-ing orbital distribution (Figure 2). Ceres-likeand Eu-ropa-like groups, which include the largest asteroids in the solar system, show an interesting orbital distribu-tion. These groups are located in the 2.6 < a < 3.6 AU region that contains the snow-line. The snow-line’s location may have been driftedinwards due to the disk’s cooling and evolution[9,10]. Recent dynamical models[11,12]suggested that a substantial fraction of primitive asteroids originated between or beyond the giant planets (a > 5 AU), where water ice would have been stable, and then implanted in the outer Main Belt region because of the giant planets’ growth.Figure 1. Two asteroids, 94 Aurora and 423 Diotima, showing spectra similar to asteroids (1)Ceres and (52)Europa.Thermal modeling and evolution of primitive aster-oids:Primitive water-rich asteroids are thought to be originally composed of mixtures of anhydrous materi-als and water ice that waslater melted by heating sources such as the decay of 26Al, reacting with anhy-drous materials to form H2O/OH-rich minerals. Calcu-lations of the evolution of the temperature and struc-ture of icy planetesimals were performedusing a 1D finite differences thermal evolution model[13, 14] for 26Al-heated planetesimals. In particular, thermally activated compaction due to hot pressing of bodies with an initially unconsolidated porous structure is included. An ice-rich initial composition that leads to a material dominated by phyllosilicates upon aqueous alteration (with 25 vol% H2O and a rock fraction that contains 85 vol% phyllosilicates and 15 vol% olivine upon aque-ous alteration, similar to CIand CM chondrites) was assumed. Atypical initial porosity of 40%[15] is re-duced following the change of the strain rate that is calculated as Voigt approximation from the strain rates of components[16]. Material properties (thermal con-ductivity, density, heat capacity, etc.) correspond to the composition assumed and are adjusted with tempera-ture and porosity. Melting of the water ice as well as water-rock separation are included[14]. Both short-and long-lived radionuclides are consideredas heat sources. Figure 3 shows the maximum temperature calculated as a function of radius and accretion time. A variety of internal structures is obtained, ranging from primordial (no melting of water ice) over partially melted or partially differentiated (melting of water ice, hydration, formation of a rocky core and water ocean below an undifferentiated layer) to completely differ-entiated ones (rocky core, water mantle, Enceladus-like case). The heating and differentiation of planetesi-mals is determined by the availability of 26Al, i.e., by the accretion time t0 relative to the formation of the calcium-aluminum-rich inclusions (CAIs), such that maximum temperatures and structures vary stronglyfor t0< 6 Myr rel. to CAIs. However, for a later accre-tion only the size of the body determines its maximum tempera-ture and structure due to the nearly constant heating by long-lived radionuclides. Average densities of Ceres-and Europa-like group membersimply high-ly porous interiors and, consequently, relatively late accretion at t0> 3 Myr rel. to CAIs with a maximum temperature of < 600 K (Figure3)

D Takir↗

An Experimental Kinetics Study of Isopropanol Pyrolysis and Oxidation behind Reflected Shock Waves

Isopropanol has potential as a future bio-derived fuel and is a promising substitute for ethanol in gasoline blends. Even so, little has been done in terms of high-temperature chemical kinetic speciation studies of this molecule. To this end, experiments were conducted in a shock tube using simultaneous CO and H2O laser absorption measurements. Water and CO formation during isopropanol pyrolysis was also examined at temperatures between 1127 and 2162 K at an average pressure of 1.42 atm. Species profiles were collected at temperatures between 1332 and 1728 K and at an average pressure of 1.26 atm for equivalence ratios of 0.5, 1.0, and 2.0 in highly diluted mixtures of 20% helium and 79.5% argon. Species profiles were also compared to four modern C3 alcohol mechanisms, including the impact of recent rate constant measurements. The Li et al. (2019) and Saggese et al. (2021) models both best predict CO and water production under pyrolysis conditions, while the AramcoMech 3.0 and Capriolo and Konnov models better predict the oxidation experimental profiles. Additionally, previous studies have collected ignition delay time (τign) data for isopropanol but are limited to low pressures in highly dilute mixtures. Therefore, real fuel–air experiments were conducted in a heated shock tube with isopropanol for stoichiometric and lean conditions at 10 and 25 atm between 942 and 1428 K. Comparisons to previous experimental results highlight the need for real fuel–air experiments and proper interpretation of shock-tube data. The AramcoMech 3.0 model over predicts τign values, while the Li et al. model severely under predicts τign. The models by Capriolo and Konnov and Saggese et al. show good agreement with experimental τign values. A sensitivity analysis using these two models highlights the underlying chemistry for isopropanol combustion at 25 atm. Additionally, modifying the Li et al. model with a recently measured reaction rate shows improvement in the model’s ability to predict CO and water profiles during dilute oxidation. Finally, a regression analysis was performed to quantify τign results from this study.

09 BIOMASS FUELS↗

Mixed Material Scintillator Systems Position Reconstruction Modelling Report

A Mixed-Material Scintillator System (MMSS) is a radiation detector using a scintillator made of a heterogeneous, structured mixture of two or more materials. MMSSs are designed so that the structured mixing of materials encodes properties of the radiation detected in the scintillator. This allows for new radiation detectors with advantages over detectors using traditional homogeneous scintillators.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Subpixel Snow-Covered-Area and Snow Grain Size From Mixture Analysis with AVIRIS Data

Snow-covered-area (SCA) and snow grain size are crucial inputs to hydrologic and climatologic modeling of alpine and other seasonally snow-covered regions. SCA is necessary to parameterize energy budget calculations in climate models, to determine in which regions point snowmelt models are to be run for distributed snowmelt modeling efforts and to provide a basis from which estimates of snow water equivalent (SWE) may be made. Snow grain size, SWE and snow impurities determine the spectral albedo of snow, which controls the net solar flux at the snowpack surface. Snow albedo is of the utmost importance in snowmelt modeling, yet the difficulty with which grain size, SWE, and impurities are mapped has left the spatial distribution of snow albedo in alpine catchments poorly understood. The Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) has been used to estimate sub-pixel snow-covered-area and snow grain size independently. In this paper we present a technique which improves estimates of both snow parameters by treating their mapping simultaneously.

Green Robert O.↗

Subpixel Snow-Covered-Area and Snow Grain Size from Mixture Analysis with Aviris Data

Snow-covered-area (SCA) and snow grain size are crucial inputs to hydrologic and climatologic modeling of alpine and other seasonally snow-covered regions. SCA is necessary to parameterize energy budget calculations in climate models, to determine in which regions point snowmelt models are to be run for distributed snowmelt modeling efforts and to provide a basis from which estimates of snow water equivalent (SWE) may be made.

hydrologic↗

Numerical modeling of hydrogen mixing in a direct-injection engine fueled with gaseous hydrogen

Hydrogen is considered as one of the most promising options to achieve effective decarbonization of the energy and transportation sectors. As such, it has recently been receiving increasing attention because of its promising potential as an energy carrier for advanced energy and propulsion systems. With a focus on internal combustion engines, direct injection (DI) of gaseous hydrogen during the compression stroke offers great potential for high engine efficiency and specific power while reducing the risk of backfiring and pre-ignition issues. Therefore, many experimental and numerical efforts have recently been dedicated to understanding the physical and chemical behaviors of hydrogen in engine during mixing and combustion. This study focuses on computational fluid dynamics (CFD) modeling of the hydrogen DI process in a hydrogen optical research engine. Under the conditions studied, gaseous hydrogen is injected into the combustion chamber via a centrally located single-hole injector at a pressure of 100 bar. Two configurations, namely low-and high-tumble, are investigated to understand the impact of different in-cylinder flow patterns on the fuel-air mixture preparation. Simulations are carried out using the commercial CFD software CONVERGE. Here, the in-cylinder turbulence is modeled with an unsteady Reynolds-averaged Navier-Stokes (URANS) formulation closed by the renormalization group (RNG) k-ε model. Several numerical methods and model constants, including but not limited to turbulent Schmidt number, are evaluated. The numerical results are systematically compared against experimental measurements of velocity and hydrogen concentration fields on the vertical center plane to assess the performance of the CFD model, unveil the physics of hydrogen mixing, and establish best practices for modeling hydrogen DI under relatively high injection pressure conditions.

33 ADVANCED PROPULSION SYSTEMS↗