SEARCH · Search NASA

Results for “Engine Simulation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 469 records · Page 26

All sources NASA NTRS DOE OSTI

DOE OSTI2025

PUFFIn Software Modeling for Quality Management

PUFFIn (PENELOPE User Friendly Fast Interface) was designed as a fast and simple Monte Carlo simulation tool for the transport of photons and electrons, with a primary purpose as a learning and education tool for a broad range of static configurations in the radiation processing industry. Development of the PUFFIn software is funded by the Office of Radiological Security (ORS) within the United States National Nuclear Security Administration (NNSA). PUFFIn helps fill the education and knowledge gaps in the industry, as identified in reports by Fermilab (2017) and the IAEA (2020). PUFFIn uses the PENELOPE (NEA-2023) physics engine to perform simulations on static configurations. PUFFIn has support for multiple geometry types from simple, single material simulations to full 3D configurations created from CAD input files or images from X-Ray Tomography scans. PUFFin was designed to be easy for the novice user, it will generate the input and geometry files required by PENLOPE and will display the output plots within the PUFFIn interface. PUFFin is distributed for free but requires a free workshop so users can be adequately trained in its use. Workshops have been presented in the past at Texas A&M university, the Aerial-CRT facility in Strasbourg France and Jakarta Indonesia. PUFFIn simulations have been validated by 10 MeV ebeam experiments done at Aerial-CRT in France (Radiation Physics and Chemistry 222 (2024) 111774). Further user experimental comparisons were made at the medical product hands on workshop at Texas A&M in October 2024.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

A multiphase flow model of water droplets dielectrophoretic-induced air dehumidification phenomena

Air humidity in indoor spaces plays a critical role in human comfort and health. Dehumidification systems are used for building humidity controls, but they can take significant energy consumption, especially in geographic locations with high outdoor humidity and warm climates. Consequently, there is a growing demand for innovative dehumidification processes that consume minimal energy. Dielectrophoretic air dehumidification represents one such promising approach. However, it has not garnered significant attention due to the absence of engineering models and simulation tools capable of evaluating its performance and limitations at large-scale airflows. A new numerical multiphase CFD model, which is also experimentally validated, is developed in a customized Reacting Foam solver based on OpenFOAM® version 9. The newly developed model seeks to decrease substantial energy consumption and lower costs by leveraging the dielectrophoretic phenomenon to regulate moisture levels in the air. The solver integrates a hybrid Eulerian-Lagrangian framework to track the droplet's trajectory and growth rate while solving the continuum equations for the moist air. An electrospray produces electrically charged droplets, which grow during their in-flight trajectories as water vapor condenses onto their surfaces. The role of electrostatic forces in promoting vapor condensation within a high-gradient electrical field is investigated, and the dielectrophoretic vapor nucleation process on charged water droplets is discussed. The CFD model was validated against results from the literature and from proof-of-concept experiments conducted by the authors, which showed a 2 % air dehumidification with a single electrospray and airflow rate of 5 cubic feet per minute. The simulation results indicated that augmenting the number of electrically charged spray droplets increased the dehumidification of the air to 25 %. The initial mean droplet diameter, the orientation of the injector and relative humidity significantly influence the assessment of dehumidification. As a result, scaling up this approach to larger airflow volumes is identified as a potential future research direction.

42 ENGINEERING↗

DOE OSTI2026

Accelerating kinetic plasma simulations with machine-learning-generated initial conditions

Computational models of plasma technologies often solve for the system operating conditions by time-stepping an initial value problem to a quasi-steady solution. However, the strongly nonlinear and multi-timescale nature of plasma dynamics often necessitate millions, or even hundreds of millions, of steps to reach convergence, reducing the effectiveness of these simulations for computer-aided engineering. We consider acceleration of kinetic plasma simulations via data-driven machine-learning-generated initial conditions, which initialize the simulations close to their final quasi-steady-state, thereby reducing the number of steps to reach convergence. Three machine-learning models are developed to predict the density and ion kinetic profiles of capacitively coupled plasma discharges relevant to the microelectronics industry. The models are trained on kinetic simulations over a range of device operating frequencies and pressures. Best performance was observed when simulations were initialized with ion kinetic profiles generated by a convolutional neural network, reducing the mean number of steps to reach convergence by 17.1× when compared to initialization with a zero-dimensional global model. We also outline a workflow for continuous data-driven model improvement and simulation speedup, with the aim of generating sufficient data for full device digital twins.

Artificial neural networks↗

DOE OSTI2022

Caldera_Grid

Caldera Grid is part of Caldera software platform, a suite of collective, open-source tools that was developed to improve the state of the art in modeling the impacts of Electric Vehicle (EV) charging on the grid. Caldera Grid is a co-simulation framework implemented using Hierarchical Engine for Large scale Infrastructure Co-Simulation (HELICS). The framework facilitates the co-simulation of EV charging and Smart Charge Management (SCM) strategies in Caldera Infrastructure Charge Model (ICM) with distribution level grid models in OpenDSS. Vehicle energy needs and charge session requirements generated using Caldera Charge Decision Module (CDM) are fed in as input to Caldera Grid. The EV charging models in Caldera ICM simulates both the uncontrolled charging loads as well as loads modified by the SCM strategies to develop distributed load profiles for each grid node hosting an Electric Vehicle Supply Equipment (EVSE) – also known as chargers. These loads were simulated in OpenDSS alongside existing distribution feeder loads. Each of these models are co-simulated in a HELICS federate. The HELICS co-simulation framework facilitates communication and synchronization between the federates. By co-simulating EV loads and distribution feeder loads, Caldera Grid can assess the potential grid impacts of EV charging on distribution feeders under various grid conditions. A control strategy federate is implemented with an interface where control strategies based on feedback from EV charging status and grid conditions can be developed and implemented. The platform can also support multiple control strategies in a single co-simulation.

Sundarrajan, ManojKumar Cebol↗

DOE OSTI2015

Upscaling ab-initio chemistry models to non-equilibrium flow simulations

Thermochemical nonequilibrium is generated by hypersonic cruise vehicles, high speed combustion and detonations. Nonequilibrium interactions affect reaction rates and energy transfer that are critical for many engineering applications. Ab-initio potential energy surfaces can be used in quasi-classical trajectory (QCT) calculations to calculate state-specific dissociation rates and cross sections, as well as state-to-state energy exchange rates and cross sections. These rates and cross sections are needed for simulating such nonequilibrium flows. However, there are two main challenges in extending QCT calculations to modeling all chemical processes found in typical engineering flowfields. First, ab-initio surfaces for many collision partners are not available and when available their fidelity is not quantified. Second, state-to-state cross section datasets are typically on the order of 1000 4, , and are too large to be efficiently implemented in engineering simulations. This work focuses on O 2 +O interactions because of prevalence of oxygen in many engineering application; however, all of the models developed here may be applied to other diatom-atom interactions and can be extended to other diatom-diatom combinations. It is shown in this work that at temperatures between 2,500 and 20,000K, state-specific cross sections are not highly sensitive to potentials and a Morse additive pairwise (MAP) potential is proposed for modeling collisions for which ab-initio potentials are not available. Next, this work suggests an approach to model state-to-state cross sections based on maximum entropy considerations. The calculated cross sections are fit to a new model with just 8 parameters, which can be easily scaled up to flow simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Application of physics-informed neural networks (PINNs) solution to coupled thermal and hydraulic processes in silty sands

Abstract The accurate modeling of water and heat transport in soils is crucial for both geo-environmental and geothermal engineering. Traditional modeling methods are problematic because they require well-defined boundaries and initial conditions. Recently, physics-informed neural networks (PINNs), which incorporate partial differential equations (PDEs) to solve forward and inverse problems, have attracted increasing attention in machine learning research. In this study, we applied PINNs to tackle hydraulic and thermal transport coupling forward problems in silty sands. A fully connected deep neural network was utilized for training. This neural network model leverages automatic differentiation to apply the governing equations as constraints, based on the mathematical approximations established by the neural network itself. We conducted forward problems and compared the solutions derived from PINNs with those from Finite Element Method (FEM) simulations. The forward problem results demonstrate the PINNs model’s capability in predicting hydraulic transport, heat transport, and thermal–hydraulic coupling in silty sands under various boundary conditions. The PINNs exhibited great performance in simulating the thermal–hydraulic coupling problem. The accuracy of the PINNs solutions shows its potential for simulation in geotechnical engineering.

Feng, Yuan↗

DOE OSTI2021

Coupling of regional geophysics and local soil-structure models in the EQSIM fault-to-structure earthquake simulation framework

Accurate understanding and quantification of the risk to critical infrastructure posed by future large earthquakes continues to be a very challenging problem. Earthquake phenomena are quite complex and traditional approaches to predicting ground motions for future earthquake events have historically been empirically based whereby measured ground motion data from historical earthquakes are homogenized into a common data set and the ground motions for future postulated earthquakes are probabilistically derived based on the historical observations. This procedure has recognized significant limitations, principally due to the fact that earthquake ground motions tend to be dictated by the particular earthquake fault rupture and geologic conditions at a given site and are thus very site-specific. Historical earthquakes recorded at different locations are often only marginally representative. There has been strong and increasing interest in utilizing large-scale, physics-based regional simulations to advance the ability to accurately predict ground motions and associated infrastructure response. However, the computational requirements for simulations at frequencies of engineering interest have proven a major barrier to employing regional scale simulations. In a U.S. Department of Energy Exascale Computing Initiative project, the EQSIM application development is underway to create a framework for fault-to-structure simulations. This framework is being prepared to exploit emerging exascale platforms in order to overcome computational limitations. This article presents the essential methodology and computational workflow employed in EQSIM to couple regional-scale geophysics models with local soil-structure models to achieve a fully integrated, complete fault-to-structure simulation framework. Here, the computational workflow, accuracy and performance of the coupling methodology are illustrated through example fault-to-structure simulations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2019

On-line Generation and Error Handling for Surrogate Models within Multifidelity Uncertainty Quantification

Uncertainty quantification is recognized as a fundamental task to obtain predictive numerical simulations. However, many realistic engineering applications require complex and computationally expensive high-fidelity numerical simulations for the accurate characterization of the system responses. Moreover, complex physical models and extreme operative conditions can easily lead to hundreds of uncertain parameters that need to be propagated through high-fidelity codes. Under these circumstances, a single fidelity approach, i.e. a workflow that only uses high-fidelity simulations to perform the uncertainty quantification task, is unfeasible due to the prohibitive overall computational cost. In recent years, multifidelity strategies have been introduced to overcome this issue. The core idea of this family of methods is to combine simulations with varying levels of fidelity/accuracy in order to obtain the multifidelity estimators or surrogates with the same accuracy of their single fidelity counterparts at a much lower computational cost. This goal is usually accomplished by defining a prioria sequence of discretization levels or physical modeling assumptions that can be used to decrease the complexity of a numerical realization and thus its computational cost. However ,less attention has been dedicated to low-fidelity models that can be built directly from the small number of high-fidelity simulations available. In this work we focus our attention on Reduced-Order Models that can be considered a particular class of data-driven approaches. Our main goal is to explore the combination of multifidelity uncertainty quantification and reduced-order models to obtain an efficient framework for propagating uncertainties through expensive numerical codes.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

TAILORED BIOBLENDSTOCKS WITH LOW ENVIRONMENTAL IMPACT TO OPTIMIZE MCCI ENGINES

The overall objective of the project is to develop and demonstrate a microalgae bio-blendstock with greater than 60% greenhouse gas reduction potential relative to petroleum diesel, that can reduce sooting propensity, increase cetane number and improve engine thermal efficiency relative to a baseline diesel engine operating on conventional fuel. Specific objectives include: (1) development of a new framework for both LCA and TEA that explicitly considers temporal variation in productivity and the frequency of crop loss; (2) determination of how fuel compounds that can be produced from the algal biomass can be ‘bio-tailored’ based on the species composition and biological production process, and subsequent processing via HTL to biocrude, and upgrading of the biocrude; (3) execution of a feedback loop (algae production  biocrude refining  combustion optimization  feedback to refining stage), for optimization of fuels for MCCI combustion; (4) optimization of MCCI combustion and emissions performance, accounting first for the biological processes that dictate the chemical composition of biocrude oil, and second for the subsequent chemical processes that comprise mixing controlled compression ignition combustion; (6) simulation of MCCI engine combustion processes to demonstrate the incorporation of relevant fuel chemistry that captures the specific impacts of optimized algal fuels.

09 BIOMASS FUELS↗

DOE OSTI2020

Tailored Bioblendstocks with Low Environmental Impact to Optimize MCCI Engines

09 BIOMASS FUELS↗

DOE OSTI2023

Phase transition in bilayer MoS 2 under tensile loading: a molecular dynamics study

Abstract Molybdenum disulfide (MoS 2 ), especially single-layer MoS 2 , has been experimentally and computationally discovered to exist in several different polymorphs exhibiting various electronic and mechanical properties. The morphology of MoS 2 can be tuned through strain engineering. Molecular dynamics simulations are conducted to systematically study the phase transition of single-layer MoS 2 and bilayer MoS 2 under the uniaxial tensile condition at room temperature. The roles of edge and S-line vacancy are investigated. Phase transitions are always triggered near the edge and vacancy sites. The initiation of the metastable T″ phase can release the tensile stress in the lattice, followed by I4/mmm phase initiation, regardless of the edge conditions. The growth of the I4/mmm phase can cause the local buckling of the MoS 2 plane. With a tilted S-line vacancy, I4/mmm phase is first initiated to reduce the local shear stress accumulated near the vacancy line. Overall, the phase transition mechanism of single layer and bilayer MoS 2 under the uniaxial tensile loading is provided, which guides the future strain engineering of MoS 2 in nanoelectronics applications.

Materials Science↗

DOE OSTI2025

CFD modeling of non-catalytic, partial-oxidation engine reformer for flare mitigation

Flaring associated natural gas is commonly employed in the oil and gas industry to reduce methane (CH 4 ) emissions but generates carbon dioxide (CO 2 ) and harmful pollutants, significantly contributing to air pollution and posing risks to public health. To mitigate this impact, M2X Energy Inc. has developed a small-scale, modular gas-to-methanol system. This system features an engine reformer that performs fuel-rich partial oxidation of wellhead gas to produce syngas—a mixture of carbon monoxide (CO) and hydrogen (H 2 )—followed by a downstream reactor for methanol synthesis. This study focused on computational fluid dynamics (CFD) modeling of the engine reformer to simulate partial oxidation chemistry, predict the rich-burn operating limit, and assess syngas quality, ultimately aiding in design and operational optimization. The CFD model, developed within a Reynolds-Averaged Navier-Stokes (RANS) turbulence framework, incorporated sub-models for turbulent combustion, a chemical mechanism with polycyclic aromatic hydrocarbon (PAH) pathways, and soot emissions to accurately capture the fuel-rich, turbulent jet ignition and combustion processes. Model validation against experimental data showed good agreement across pre- and main-chamber pressures, apparent heat release rates, and exhaust gas concentrations of key species (H 2 , CO, CO 2 , CH 4 ) for varying intake equivalence ratios. Here, the model identified a rich-burn operating limit near a fuel-air equivalence ratio of 2.35, consistent with experimental observations. Furthermore, syngas quality analysis revealed that extending the rich-burn limit through engine reformer optimization could enhance syngas production, contributing to higher methanol synthesis efficiency.

Computational Fluid Dynamics↗

DOE OSTI2020

A phenomenological rate of injection model for predicting fuel injection with application to mixture formation in light-duty diesel engines

Fuel injection rate laws are one of the most important pieces of information needed when modeling engine combustion with computational fluid dynamics. In this study, a simple phenomenological model of a common-rail injector was developed and calibrated for the Bosch CRI2.2 platform. The model requires three tunable parameter fits, making it relatively easy to calibrate and suitable for injector modeling when high-fidelity information about the internal injector’s geometry and electrical circuit details are not available. Each injection pulse is modeled as a sequence of up to four stages: an injection needle mechanical opening transient, a full-lift viscous flow inertial transient, a Bernoulli steady-state stage, and a needle descent transient. Parameters for each stage are obtained as polynomial fits from measured injection rate properties. The model enforces total injected mass, and the intermediate stages are only introduced if the injection pulse duration is long enough. Experimental rates of injection from two separate campaigns on the same injector were used to calibrate the model. The model was first validated against measured injection rate laws featuring pilot injections, short partially premixed combustion pulses, and conventional diesel combustion injection strategies. Then, it was employed as an input to engine computational fluid dynamics simulations, which were run to simulate experiments of mixture formation in an optically accessible light-duty diesel engine. It was found that, though simple, this model is capable of predicting both pilot and main injection pulse mass flow rates well: the simulations yielded accurate predictions of in-cylinder equivalence ratio distributions from injection strategies for both partially premixed combustion and pilot injections. Also, once calibrated, the model produced appropriate results for a wide range of injected mass and rail pressure values. Finally, it was observed that usage of such a relatively simple model can be a good choice when high-fidelity injection rate input and highly detailed information of the injector’s geometry and operation are not available, particularly as noticeable discrepancies can be present also among different experimental campaigns on similar hardware.

42 ENGINEERING↗

DOE OSTI2020

SENSITIVITY ANALYSIS OF FUEL PHYSICAL PROPERTY EFFECTS ON SPARK IGNITION ENGINE PERFORMANCE

Alternative fuels are of interest to automakers and regulators due to their potential to reduce net greenhouse gas emissions from transportation sources. Alternative fuels also have fuel properties which may enable advanced combustion modes with higher engine thermal efficiencies. There has been previous work to identify the relationship between various fuel properties and engine performance, but most of this work has been experiments or simulations where the change in properties was obtained through changing the fuel composition, making isolating the effects of individual fuel properties difficult. In this study, numerical simulations have been used to investigate the effects of individual fuel physical properties such as viscosity or heat of vaporization (HoV) on engine performance. Simulations have been performed of two different engine platforms, the first an optical, single-cylinder research engine and the second a multi-cylinder production engine. Both engines are direct-injection spark-ignition engines with pent-roof heads and are designed for automotive applications. Each engine was run at a different operating condition, one stable and one knock-limited. Different base fuels provided a variety of simulated conditions. Up to six different fuel properties were varied as part of Global Sensitivity Analyses performed for each of the engines with multiple performance targets including thermal efficiency, combustion efficiency and combustion phasing. Results show trends that are largely consistent with previous experimental findings using multiple fuels. The engine thermal efficiency was primarily sensitive to the fuel's HoV, with other fuel physical properties having smaller effects. For optical engine results, the magnitude of the effect was greater in this study than expected based on previous experimental results were many fuel physical and chemical properties were varied simultaneously. However, for the multi-cylinder production engine, the relationship between thermal efficiency and HoV was slightly smaller.

engine simulation↗

DOE OSTI2025

Breaking the passivation barrier via d-p orbital optimization for stable hydrogen production and sulfion upgrading

The development of energy-efficient hydrogen production technologies represents a critical pathway toward achieving global carbon neutrality objectives. This work provides fundamental insights into overcoming catalyst passivation challenges in sulfide oxidation reaction (SOR)-coupled hydrogen evolution reaction (HER) systems through precise orbital hybridization engineering. Our theoretical simulations reveal that sulfur-passivated ruthenium surfaces can effectively modulate d-p orbital hybridization, significantly reduce d-electron activity while stabilizing long-chain S 8 species and decreasing intermediate adsorption energies. Furthermore, metal carbides/ruthenium heterostructure (MC/Ru, M = V, Mo, W) was designed to achieve simultaneous optimization of both HER (∆G H* = −0.11 eV) and SOR (∆G RDS = 1.51 eV) via work-function-mediated interfacial electron transfer, which effectively tailors surface electronic states. Guided by theoretical predictions, we successfully synthesized a series of metal carbides/ruthenium/nitrogen-doped carbon catalysts based on a solid-phase reaction and designated as MC/Ru@NC (M = V, Mo, W). The optimized VC/Ru@NC catalyst exhibits exceptional performance in a membrane-free two-electrode system, achieving an ultralow cell voltage of 0.76 V at 200 mA cm −2 with outstanding stability over 1600 h, while maintaining 97.5 % Faradaic efficiency for hydrogen production and 73.8 % sulfur recovery efficiency.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

Numerical study of a micro flow reactor at engine pressures: Flames with repetitive extinction and ignition and simulations with a reduced chemical model

Combustion in a narrow channel with an imposed temperature gradient is studied numerically at elevated pressure with engine-relevant fuels. Here, the focus is placed on unsteady flames with repetitive extinction and ignition (FREI) to determine the potential of this regime for fuel testing and calibration of reduced chemical mechanisms. First, it is shown that the FREI regime does occur at elevated pressures for sufficiently small tube diameters. The sensitivity of the extinction and ignition temperatures to low-temperature chemistry is found to be significantly enhanced at 25 bar compared to atmospheric conditions. The ignition and extinction temperatures differ by up to 100 K between PRF mixtures with varying octane numbers. Ternary mixtures of iso-octane/n-heptane/toluene and iso-octane/n-heptane/ethanol at similar research and motor octane numbers are also studied. Second, the potential of using data from the micro flow reactor to infer reaction rates is assessed. A reduced chemical mechanism combining a small fuel-dependent submechanism with a detailed fuel-independent submechanism for the core species chemistry is used for that purpose. Only the most sensitive fuel-dependent reactions are inferred. the calibrated reduced model is compared to a detailed model and good agreement in ignition delay times and laminar flame speeds is observed. Finally, this illustrates the potential of a micro flow reactor with a controlled temperature profile to test fuels and infer kinetic data.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2026

Parametric amplification of spin-motion coupling in three-dimensional trapped-ion crystals

Three-dimensional (3D) crystals offer a route to scaling up trapped-ion systems for quantum sensing and quantum simulation applications; however, engineering coherent spin-motion couplings and effective spin-spin interactions in large crystals poses technical challenges associated with decoherence and prolonged timescales to generate appreciable entanglement. Here, we explore the possibility of speeding up these interactions in 3D crystals via parametric amplification. For this purpose, we derive a general Hamiltonian for the parametric amplification of spin-motion coupling that is broadly applicable to normal modes with motion transverse to or along the spatial extent of the crystal. Unlike in lower-dimensional crystals, we find that the ability to faithfully (uniformly) amplify the spin-spin interactions in 3D crystals depends on the physical implementation of the spin-motion coupling. We consider the light-shift gate, and the so-called phase-insensitive and phase-sensitive Mølmer-Sørensen (MS) gates, and we find that only the phase-sensitive MS gate can be faithfully amplified in general 3D crystals. We discuss a situation where nonuniform amplification can be advantageous. We also reconsider the effect of counter-rotating terms on parametric amplification and find that they are not as detrimental as previous studies suggest.

Optical parametric oscillators & amplifiers↗

DOE OSTI2024

Machine Learning Aided Modeling of Granular Materials: A Review

Artificial intelligence (AI) has become a buzzy word since Google’s AlphaGo beat a world champion in 2017. In the past five years, machine learning as a subset of the broader category of AI has obtained considerable attention in the research community of granular materials. This work offers a detailed review of the recent advances in machine learning-aided studies of granular materials from the particle-particle interaction at the grain level to the macroscopic simulations of granular flow. This work will start with the application of machine learning in the microscopic particle-particle interaction and associated contact models. Then, different neural networks for learning the constitutive behaviour of granular materials will be reviewed and compared. Finally, the macroscopic simulations of practical engineering or boundary value problems based on the combination of neural networks and numerical methods are discussed. We hope readers will have a clear idea of the development of machine learning-aided modelling of granular materials via this comprehensive review work.

42 ENGINEERING↗