Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic simulations”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Computational modeling of charge hopping dynamics along a disordered one-dimensional wire with energy gradients in quantum environments

This computational study investigates the effects of energy gradients on charge hopping dynamics along a one-dimensional chain of discrete sites coupled to quantum bath, which is modeled at the level of Pauli master equation (PME). This study also assesses the performance of different approximations for the hopping rates. Three different methods for solving the PME, a fourth order Runge–Kutta method, numerical diagonalization of the rate matrix followed by analytic propagation, and kinetic Monte Carlo simulation method, are tested and confirmed to produce virtually identical values of time dependent mean square displacement, diffusion constant, and mobility. Five different rate expressions, exact numerical evaluation of Fermi’s Golden Rule (FGR) rate, stationary phase interpolation (SPI) approximation, semiclassical approximation, classical Marcus rate, and Miller–Abrahams rate, are tested to help understand the effects of approximations in representing quantum environments in the presence of energy gradients. The results based on direct numerical evaluation of FGR rate exhibit transition from diffusive to non-diffusive behavior with the increase in the gradient and show that the charge transport in the quantum bath is more sensitive to the magnitude of the gradient and the disorder than in the classical bath. Among all the four approximations for the hopping rates, the SPI approximation is confirmed to work best overall. A comparison of two different methods to calculate the mobility identifies drift motion of the population distribution as the major source of non-diffusive behavior and provides more reliable information on the contribution of quantum bath.

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

Overestimated Eddy Kinetic Energy in the Eddy‐Rich Regions Simulated by Eddy‐Resolving Global Ocean–Sea Ice Models

Abstract The performance of eddy‐resolving global ocean–sea ice models in simulating mesoscale eddies is evaluated using six eddy‐resolving experiments forced by different atmospheric reanalysis products. Interestingly, eddy‐resolving ocean general circulation models (OGCMs) tend to simulate more (less) energetic eddy‐rich (eddy‐poor) regions with a smaller (larger) spatial extent than satellite observation, which finally shows that larger (smaller) mesoscale energy intensity (EI) is simulated in the eddy‐rich (eddy‐poor) regions. Quantitatively, there is an approximately 27%–60% overestimation of EI in the eddy‐rich regions, which are mainly located in the Kuroshio–Oyashio Extension, the Gulf Stream, and the Antarctic Circumpolar Currents regions, although the global mean EI is underestimated by 25%–45%. Apparently, the eddy kinetic energy in the eddy‐poor region is underestimated. Further analyses based on coherent mesoscale eddy properties show that the overestimation in the eddy‐rich regions is mainly attributed to mesoscale eddies’ intensity and is more prominent when mesoscale eddies are in their growth stage.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Electroconvection near an ion-selective surface with Butler–Volmer kinetics

We study the effects of interfacial kinetics on the electro-hydrodynamics of ion transport near an ion-selective surface using a combination of linear stability analysis and numerical simulation. The finite kinetics of the electrolyte–electrode interface affects the ion transfer and electroconvection in many ways. On a surface of fixed topography, such as a metal surface of slow and stable ion deposition or covered by a polymer membrane, the finite kinetics reduces the current in one-dimensional ion diffusion/migration, increases the critical voltage for the onset of the electroconvective instability, changes the dynamics of the electroconvection and the overlimiting current, and enhances the lateral ion diffusion within the interfacial layer. The first three effects are indirectly caused by the reaction kinetics and can be characterized by an effective voltage difference across the liquid electrolyte. In comparison, the last effect is controlled by a direct interplay between kinetics and nonlinear electroconvection. Scaling laws for ion transport and features of electroconvection are proposed. We also analyse the linear stability of a surface which evolves under ion deposition and find that the finite kinetics decreases the growth rate of both electroconvective and morphological instabilities and therefore modifies the wavenumber of the most unstable mode.

Mechanics↗
DOE OSTI2020

Magnetopause Reconnection as Influenced by the Dipole Tilt Under Southward IMF Conditions: Hybrid Simulation and MMS Observation

Abstract Using a three‐dimensional (3‐D) global‐scale hybrid code, the Magnetospheric Multiscale (MMS) reconnection event around 02:13 UT on 18 November 2015, highlighted in the Geospace Environment Modeling (GEM) Dayside Kinetic Challenge, is simulated, in which the interplanetary magnetic field (IMF) points southward and the geomagnetic field has a −27° dipole tilt angle. Strong southward plasma jets are found near the magnetopause as a result of the dayside reconnection. Our results indicate that the subsolar magnetopause reconnection X line shifts from the subsolar point toward the Northern Hemisphere due to the effect of the tilted geomagnetic dipole angle, consistent with the MMS observation. Subsequently, the reconnection X lines or sites and reconnection flux ropes above the equator propagate northward along the magnetopause. The formation and global distribution of the X lines and the structure of the magnetopause reconnection are investigated in detail with the simulation. Mirror mode waves are also found in the middle of the magnetosheath downstream of the quasi‐perpendicular shock where the plasma properties are consistent with the mirror instability condition. As a special outcome of the GEM challenge event, the spatial and temporal variations in reconnection, the electromagnetic power spectra, and the associated D‐shaped ion velocity distributions in the simulated reconnection event are compared with the MMS observation.

Guo, Zhifang↗
DOE OSTI2023

Catalyst: Fast and flexible modeling of reaction networks

We introduce Catalyst.jl, a flexible and feature-filled Julia library for modeling and high-performance simulation of chemical reaction networks (CRNs). Catalyst supports simulating stochastic chemical kinetics (jump process), chemical Langevin equation (stochastic differential equation), and reaction rate equation (ordinary differential equation) representations for CRNs. Through comprehensive benchmarks, we demonstrate that Catalyst simulation runtimes are often one to two orders of magnitude faster than other popular tools. More broadly, Catalyst acts as both a domain-specific language and an intermediate representation for symbolically encoding CRN models as Julia-native objects. This enables a pipeline of symbolically specifying, analyzing, and modifying CRNs; converting Catalyst models to symbolic representations of concrete mathematical models; and generating compiled code for numerical solvers. Leveraging ModelingToolkit.jl and Symbolics.jl, Catalyst models can be analyzed, simplified, and compiled into optimized representations for use in numerical solvers. Finally, we demonstrate Catalyst’s broad extensibility and composability by highlighting how it can compose with a variety of Julia libraries, and how existing open-source biological modeling projects have extended its intermediate representation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

A Novel use of Direct Simulation Monte-Carlo to Model Dynamics of COVID-19 Pandemic Spread

In this report, we evaluate a novel method for modeling the spread of COVID-19 pandemic. In this new approach we leverage methods and algorithms developed for fully-kinetic plasma physics simulations using Particle-In-Cell (PIC) Direct Simulation Monte-Carlo (DSMC) models. This approach then leverages Sandia-unique simulation capabilities, and High-Performance Computer (HPC) resources and expertise in particle-particle interactions using stochastic processes. Our hypothesis is that this approach would provide a more efficient platform with assumptions based on physical data that would then enable the user to assess the impact of mitigation strategies and forecast different phases of infection. This work addresses key scientific questions related to the assumptions this new approach must make to model the interactions of people using algorithms typically used for modeling particle interactions in physics codes (kinetic plasma, gas dynamics). The model developed uses rational/physical inputs while also providing critical insight; the results could serve as inputs to, or alternatives for, existing models. The model work presented was developed over a four-week time frame, thus far showing promising results and many ways in which this model/approach could be improved. This work is aimed at providing a proof-of-concept for this new pandemic modeling approach, which could have an immediate impact on the COVID-19 pandemic modeling, while laying a basis to model future pandemic scenarios in a manner that is timely and efficient. Additionally, this new approach provides new visualization tools to help epidemiologists comprehend and articulate the spread of this and other pandemics as well as a more general tool to determine key parameters needed in order to better predict pandemic modeling in the future. In the report we describe our model for pandemic modeling, apply this model to COVID-19 data for New York City (NYC), assess model sensitivities to different inputs and parameters and , finally, propagate the model forward under different conditions to assess the effects of mitigation and associated timing. Finally, our approach will help understand the role of asymptomatic cases, and could be extended to elucidate the role of recovered individuals in the second round of the infection, which is currently being ignored.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Challenges in integrating dissolved organic matter chemodiversity into kinetic models of soil respiration

The chemodiversity of dissolved organic matter (DOM) in soil has been proposed to influence the microbial metabolism and fate of belowground organic carbon (C). However, integrating DOM chemistry into soil C cycle models to improve predictions of C stocks and fluxes—beyond simply considering DOM pool size—remains a challenge. While recent research suggests that incorporating DOM chemodiversity into models can improve predictions of microbial respiration, there is still a lack of mechanistic understanding describing how DOM chemodiversity affects microbial metabolism and soil respiration. Here, we evaluated whether DOM chemodiversity was a determinant of soil respiration using paired measurements of high-resolution DOM chemistry, obtained from Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS), and potential soil respiration rates from across the United States (U.S.), all data provided by the Molecular Observation Network. Our objectives were to (1) assess statistical relationships between DOM chemodiversity and microbial respiration, and (2) evaluate the ability of kinetic models to leverage DOM chemistry to explain empirical relationships found in statistical models. Statistical regressions revealed that DOM chemodiversity (alpha diversity) was nonlinearly related to potential soil respiration rates, both independently and through its interactions with DOM and total C concentrations. In soils with relatively high DOM but low total C concentrations, potential soil respiration rates were negatively correlated with DOM alpha diversity, whereas in soils with relatively low DOM and high total C concentrations showed the opposite trend. However, when metabolic transition theory kinetic models were modified to include chemodiversity, their performance was comparable to traditional Monod kinetics approaches, which simulate respiration rates as a function of DOM concentration. The inability to account for nonlinearities in DOM chemodiversity–respiration relationships highlight an opportunity to advance substrate uptake kinetics by establishing causal links between DOM chemodiversity, microbial metabolism trade-offs, and potential interactions under varied environmental conditions.

Bioenergetic model↗
DOE OSTI2023

Modeling–Experiment–Theory Analysis of Reactions Initiated from Cl + Methyl Formate

Methyl formate (MF; CH 3 OCHO) is the smallest representative of esters, which are common components of biodiesel. The present study characterizes the thermal dissociation kinetics of the radicals formed by H atom abstraction from MF—CH 3 OCO and CH 2 OCHO—through a combination of modeling, experiment, and theory. For the experimental effort, excimer laser photolysis of Cl 2 was used as a source of Cl atoms to initiate reactions with MF in the gas phase. Time-resolved species profiles of MF, Cl 2 , HCl, CO 2 , CH 3 , CH 3 Cl, CH 2 O, and CH 2 ClOCHO were measured and quantified using photoionization mass spectrometry at temperatures of 400–750 K and 10 Torr. The experimental data were simulated using a kinetic model, which was informed by ab initio-based theoretical kinetics calculations and included chlorine chemistry and secondary reactions of radical decomposition products. Here, we calculated the rate coefficients for the H-abstraction reactions Cl + MF → HCl + CH 3 OCO (R1a) and Cl + MF → HCl + CH 2 OCHO (R1b): k 1a,theory = 6.71 × 10 –15 ·T 1.14 ·exp(—606/T) cm 3 /molecule·s; k 1b,theory = 4.67 × 10 –18 ·T 2.21 ·exp(—245/T) cm 3 /molecule·s over T = 200–2000 K. Electronic structure calculations indicate that the barriers to CH 3 OCO and CH 2 OCHO dissociation are 13.7 and 31.6 kcal/mol and lead to CH 3 + CO 2 (R3) and CH 2 O + HCO (R5), respectively. The master equation-based theoretical rate coefficients are k 3,theory (P = ∞) = 2.94 × 10 9 ·T 1.21 ·exp(—6209/T) s –1 and k 5,theory (P = ∞) = 8.45 × 10 8 ·T 1.39 ·exp(—15132/T) s –1 over T = 300–1500 K. The calculated branching fractions into R1a and R1b and the rate coefficient for R5 were validated by modeling of the experimental species time profiles and found to be in excellent agreement with theory. Additionally, we found that the bimolecular reactions CH 2 OCHO + Cl, CH 2 OCHO + Cl 2 , and CH 3 + Cl 2 were critical to accurately model the experimental data and constrain the kinetics of MF-radicals. Inclusion of the kinetic parameters determined in this study showed a significant impact on combustion simulations of larger methyl esters, which are considered as biodiesel surrogates.

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

Isotope effects on transport in LHD

Abstract Isotope effects are one of the most important issues for predicting future reactor operations. Large helical device (LHD) is the presently working largest stellarator/helical device using super conducting helical coils. In LHD, deuterium experiments started in 2017. Extensive studies regarding isotope effects on transport have been carried out. In this paper, the results of isotope effect studies in LHD are reported. The systematic studies were performed adjusting operational parameters and nondimensional parameters. In L mode like normal confinement plasma, where internal and edge transport barriers are not formed, the scaling of global energy confinement time ( τ E ) with operational parameters shows positive mass dependence ( M 0.27 ; where M is effective ion mass) in electron cyclotron heating plasma and no mass dependence ( M 0.0 ) in neutral beam injection heating plasma. The non-negative ion mass dependence is anti-gyro-Bohm scaling. The role of the turbulence in isotope effects was also found by turbulence measurements and gyrokinetic simulation. Better accessibility to electron and ion internal transport barrier (ITB) plasma is found in deuterium (D) plasma than in hydrogen (H). Gyro kinetic non-linear simulation shows reduced ion heat flux due to the larger generation of zonal flow in deuterium plasma. Peaked carbon density profile plays a prominent role in reducing ion energy transport in ITB plasma. This is evident only in plasma with deuterium ions. New findings on the mixing and non-mixing states of D and H particle transports are reported. In the mixing state, ion particle diffusivities are higher than electron particle diffusivities and D and H ion density profiles are almost identical. In the non-mixing state, ion particle diffusivity is much lower than electron diffusivity. Deuterium and hydrogen ion profiles are clearly different. Different turbulence structures were found in the mixing and non-mixing states suggesting different turbulence modes play a role.

Physics↗
DOE OSTI2024

Phase field modeling of hydrogen release in nitrile-butadiene rubber composites after high-pressure hydrogen exposure

Experimental data showed that a simple diffusion model doesn’t suffice to explain the hydrogen (H 2 ) release kinetics in nitrile butadiene rubber (NBR) composites (one with no filler or plasticizer, one with plasticizer only, one with fillers and plasticizer, and one with fillers only) after exposed to hydrogen gas at a pressure of 27.6 MPa. In this work, we developed a phase field model that considers the effect of H 2 gas bubble evolution on H 2 release kinetics. With the model, we simulated the effect of thermodynamic and kinetics properties on the nucleation, growth, and shrinkage of gas bubbles and H 2 release kinetics during and after decompression. Here, the results demonstrated that 1) the model built upon certain thermodynamic and kinetics properties of a given material can well describe its H 2 release kinetics measured in experiments, and 2) the model can predict the gas bubble evolution, which is associated with the material property degradation and failures, with accurate thermodynamic and kinetics properties or H 2 release data from experiments.

08 HYDROGEN↗
DOE OSTI2023

Molecular Dynamics Simulation of Pore-Size Effects on Gas Adsorption Kinetics in Zeolites

Strong gas-mineral interactions or slow adsorption kinetics require a molecular-level understanding of both adsorption and diffusion for these interactions to be properly described in transport models. In this combined molecular simulation and experimental study, noble gas adsorption and mobility is investigated in two naturally abundant zeolites whose pores are similar in size (clinoptilolite) and greater than (mordenite) the gas diameters. Simulated adsorption isotherms obtained from grand canonical Monte Carlo simulations indicate that both zeolites can accommodate even the largest gas (Rn). However, gas mobility in clinoptilolite is significantly hindered at pore-limiting window sites, as seen from molecular dynamics simulations in both bulk and slab zeolite models. Further, experimental gas adsorption isotherms for clinoptilolite confirm the presence of a kinetic barrier to Xe uptake, resulting in the unusual property of reverse Kr/Xe selectivity. Finally, a kinetic model is used to fit the simulated gas loading profiles, allowing a comparison of trends in gas diffusivity in the zeolite pores.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Air Classification of Forestry Residues for Fast Pyrolysis

Understanding critical biomass attributes through efficient fractionation is crucial for advancing sustainable pyrolysis for renewable energy and chemical production. This study investigates the intricate relationship between biomass preprocessing and pyrolysis product yields, employing the air classification technique for the treatment of loblolly pine residues with varying moisture content. A comprehensive exploration of the physicochemical properties of air-classified loblolly pine informs a sophisticated pyrolysis simulation model. Given the complex and multifaceted nature of biomass pyrolysis, operating across diverse temporal and spatial scales, a pyrolysis kinetics-based CFD–DEM simulation method is employed to predict product yields. Results showed that the elevated moisture content amplifies particle adhesiveness, necessitating augmented air velocities for effective separation, thereby influencing the efficiency of the separation process. While carbon and hydrogen contents exhibit relative stability across diverse moisture contents and blower frequencies, the oxygen content undergoes noticeable changes. For example, the oxygen contents were measured as 29.2 and 38.6 wt% in the light fraction of 30% moisture content sample at blower frequencies of 10 and 20 Hz, respectively. An intriguing finding emerges from pyrolysis simulation, indicating that a lower blower frequency in air classification moderately enhances bio-oil yield and significantly improves its quality, particularly in terms of water content. For instance, the water content in the bio-oil was about 1.5% and 10% in the heavy and light fractions, respectively from 10% moisture sample under 15 Hz blower frequency.

09 - BIOMASS FUELS↗
DOE OSTI2026

Catalytic resonance theory for the kinetics of photon-promoted catalysis

The illumination of catalytic surfaces with a continuous or pulsed stream of photons dynamically modulates surface chemistry for faster rates, non-equilibrium conversion, or product selectivity control. To establish fundamental principles of dynamic photon-modulated catalysis, the photocatalytic conversion of a generic surface reaction was simulated using the kinetic Monte Carlo method to understand the kinetic implications of an independent stream of photons that promotes surface product desorption. The time-averaged photocatalytic rate at differential conditions for varying photon flux and temperatures indicated three kinetic regimes described by product thermal desorption control, surface reaction control, and an intermediate kinetic regime with a zero slope Arrhenius plot, consistent with a degree of rate control dominated by the photon arrival frequency (i.e., per-site photon flux). Here, the maximum photocatalytic rate occurred orders of magnitude above the Sabatier limit at the resonance frequency, identified as the photon arrival frequency matching the surface reaction rate constant.

Canavan, Jesse R. [University of Minnesota, Minnea↗
DOE OSTI2021

Computational Investigation of the Role of Active Site Heterogeneity for aSupported Organovanadium(III) Hydrogenation Catalyst

A crucial consideration for supported heterogeneous catalysts is the nonuniformity of the active sites, particularly for supported organometallic catalysts. Standard spectroscopic techniques, such as X-ray absorption spectroscopy, reflect the nature of the most populated sites, which are often intrinsically structurally distinct from the most active catalytic sites. Additionally, with computational models, often, only a few representative structures are used to depict catalytic active sites on a surface, even though there are numerous observable factors of surface heterogeneity that contribute to the kinetically favorable active species. A previously reported study on the mechanism of a surface organovanadium-(III) catalyst [(SiO 2 )V III (Mes)(THF)] for styrene hydrogenation yielded two possible mechanisms: heterolytic cleavage and redox cycling. These two mechanistic scenarios are challenging to differentiate experimentally since the kinetic readouts of the catalyst are identical. To showcase the importance of modeling surface heterogeneity and its effect on catalytic activity, density functional theory (DFT) computational models of a series of potential active sites of [(SiO 2 )V III (Mes)(THF)] for the reaction pathways are applied in combination with kinetic Monte Carlo (kMC) simulations. Computed results were then compared to the previously reported experimental kinetic study: (1) DFT free-energy reaction pathways indicated the likely active site and pathway for styrene hydrogenation, a heterolytic cleavage pathway requiring a bare tripodal vanadium site. (2) From the kMC simulations, a mixture of different bond lengths from the support oxygen to the metal center was required to qualitatively describe the experimentally observed kinetic aspects of a supported organovanadium(III) catalyst for olefin hydrogenation. This work underscores the importance of modeling surface heterogeneity in computational catalysis.

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

Modeling discontinuous dynamic recrystallization containing second phase particles in magnesium alloys utilizing phase field method

Pre-existing second phase particles can induce particle stimulated nucleation (PSN) and pinning effect on grain boundaries during DRX process of magnesium alloys. The interaction among pinning effect, PSN mechanism and strain-induced grain boundary migration (SIBM) nucleation mechanism imposes challenge in quantitative study of microstructure evolution. A phase field model that describes discontinuous dynamic recrystallization process considering second phase particles (PF-DDRX-SP) which mainly distribute along grain boundaries is developed. The second phase order parameter with diffusion interface is introduced to the free energy function. With considering the pinning effect of second phase particles on grain boundaries, the grain boundary energy couples the interaction between second phase and grains. The effect of second phase on dislocation density evolution is considered in the model, and the PSN and SIBM mechanism are both executed by limiting the nucleation position to the grain boundary and the second phase interface. For the validation, the initial topology consistent with initial microstructure is applied to the simulation of DRX kinetics and flow behavior of AZ80 magnesium alloy, and great reliability and accuracy of developed model is indicated. Furthermore, through the model, the microstructure evolution with different particle size and with different volume fraction of round particle during DRX process were simulated. Finally, the results show that particles with higher dispersion reduce the number of DRX nuclei, which indicate that second phase particles suppress the nucleation at grain boundaries because of the pinning effect. More importantly, the proposed model could also quantitatively predict the relationships between the parameters of second phase particles and DRX behaviors, and enable to optimize the initial second phase structure in a uniform grain structure during thermomechanical process.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Comparing ab initio and quantum-kinetic approaches to electron transport in warm dense matter

Accurate knowledge of the electronic transport properties of warm dense matter is one of the main concerns of research in high-energy-density physics. Three modern approaches with vastly different levels of fidelity are reviewed and compared: the Kubo–Greenwood (KG) approach based on density-functional-theory molecular dynamics simulations (QMD), quantum kinetic theory based on average-atom models, and time-dependent density functional theory. Throughout, emphasis is placed on the connection between static properties of the electrons (e.g., density of states) and transport properties. Overall, it is found that whenever the conduction electrons can be modeled as being nearly free, fair to excellent agreement is found between QMD and kinetic theory approaches. Such a circumstance is required for modeling warm dense matter as a plasma of ions and free electrons, which is assumed in most kinetic theory approaches. The sensitivity of transport properties to the electronic structure is further highlighted by comparing different exchange–correlation approximations in QMD and KG calculations. It is found that the inclusion of exact exchange via thermal hybrid functionals can make a pronounced impact on electrical and thermal conduction in warm dense matter. We also investigate dynamic screening physics via kinetic theory and time-dependent density functional theory calculations of the mean free path of an electron in a hot dense plasma. In sum, we identify three axes along which to make progress in predicting electron transport in warm dense matter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Improved Chemical Kinetics and Algorithms for More Accurate, Faster Simulations

Internal combustion engine design is increasingly driven by computational models used to predict change in performance as a result of design changes, which previously would have been made by limited intuition or expensive and time-consuming physical testing. Improved model capabilities shorten design cycles and enable the production of cleaner and more efficient engines. This project focuses on advancing the state-of-the-art in internal combustion engine simulations. The overarching goal is to enable predictive models and reduced time to solution for simulations that impact combustion engine design.

02 PETROLEUM↗
DOE OSTI2021

Improved Chemical Kinetics and Algorithms for More Accurate, Faster Simulations

Internal combustion engine design is increasingly driven by computational models used to predict change in performance due to change of design. The design process previously depended on limited intuition and expensive and time-consuming physical testing. Improved model capabilities shorten design cycles and enable the production of cleaner and more efficient engines. This project focuses on advancing the state of the art in internal combustion engine simulations. The overarching goal is to enable predictive models and reduced time to solution for simulations that impact combustion engine design.

02 PETROLEUM↗