Search NASA⌕ Search

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.

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

At least 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A real-time distributed solid oxide electrolysis cell (SOEC) model for cyber-physical simulation

System integration and dynamic operability between SOEC and balance-of-plant (BoP) components are major technical challenges before realizing rapid load following of SOEC systems. Cyber-physical simulation (CPS) is a leading-edge digital engineering approach and is regarded as the next step beyond Digital Twins. CPS approach can be used to research SOEC system integration and develop dynamic controls prior to actual pilot testing without using a real SOEC. To seamlessly couple with BoP hardware and access non-observable operational parameters (e.g., local temperature gradient) during transients, a distributed one-dimensional (1D) real-time SOEC model was developed. Its real-time execution was demonstrated for 20 to 640 nodes at the fixed time step of 5 ms. A higher excess air ratio enabled smaller local temperature gradients on SOEC solid materials and faster transients upon current density step change from 0.15 to 0.55 A cm -2 . During the transients, the magnitude of the peak temperature gradient nearly doubled in 10 s from -3.5 to -5.9 °C cm -1 . This represents a significant operating risk that can impact the dynamic operability of SOEC systems. In addition, the local temperature gradient was found to change directions on all nodes in SOEC solid materials, with the greatest impact on the upstream nodes. The SOEC model was also tested at the thermal neutral voltage using actual process air flow parameters as variable model inputs. Variable process air temperatures were found to induce alternating local temperature gradients on SOEC solid materials. These are new operational mechanisms for SOEC degradation relevant for load following operational modes yet distinct from previous reports. To mitigate these unfavorable features, the SOEC can be operated at voltages that are slightly (±20 mV) deviated from the thermal neutral voltage. Here, the corresponding net thermal energy change was less than 1.6% of the electric power consumption. This 1D real-time SOEC model established the basis of cyber-physical simulation of SOEC hybrid systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Next Gen High Efficiency Boosted Engine Development

This work represents an advanced engineering research project partially funded by the U.S. Department of Energy (DOE). Ford Motor Company, FEV North America, and Oak Ridge National Laboratory collaborated to develop a next generation boosted spark ignited engine concept. The project goals, specified by the DOE, were 23% improved fuel economy and 15% reduced weight relative to a 2015 or newer light-duty vehicle. The fuel economy goal was achieved by designing an engine incorporating high geometric compression ratio, high dilution tolerance, low pumping work, and low friction. The increased tendency for knock with high compression ratio was addressed using early intake valve closing (EIVC), cooled exhaust gas recirculation (EGR), an active pre-chamber ignition system, and careful management of the fresh charge temperature. Engine weight reduction measures were implemented throughout the engine system making use of composite materials, advanced manufacturing techniques, and architectural choices. This report outlines the analytical, design, fabrication, and test work conducted for the duration of the project. The combustion system stability, EGR tolerance, and knock resistance were validated on a single cylinder engine. An inline six-cylinder engine was then designed targeting application in the Ford F150. Multi-cylinder engines were produced and tested achieving the target vehicle fuel economy improvement of 23% assessed using measured engine fuel consumption combined with a vehicle drive cycle simulation. Actions were identified and designs were demonstrated to achieve the 15% weight reduction target. This project included items covering a range of technology readiness levels. Some of the technologies explored are production ready, while others were investigated to understand the limitations for what can be achieved in a stoichiometric, gasoline-fueled, spark-ignited internal combustion engine.

42 ENGINEERING↗
DOE OSTI2017

Thermal-Hydraulic Results for the Aboveground Configuration of a Dry Cask Simulator

The thermal performance of a commercial nuclear spent fuel dry storage cask is evaluated through detailed analytical modeling. These modeling efforts are performed by the vendor to demonstrate the performance as well as regulatory compliance and are independently verified by the Nuclear Regulatory Commission (NRC). The majority of commercial dry storage casks in use today are aboveground. Cooling of the assemblies located inside the sealed canister is achieved by naturally-convective air flowing in the bottom of the cask and exiting out the top. Carefully measured data sets generated from testing of full sized casks or smaller cask analogs are widely recognized as vital for validating design and performance models. Numerous studies have been previously conducted. Recent advances in dry storage cask designs have significantly increased the maximum thermal load allowed in a cask in part by increasing the efficiency of internal conduction pathways and by increasing the internal convection through greater canister helium pressure. These vertical, canistered cask systems rely on ventilation between the canister and the overpack to convect heat away from the canister to the environment for both above and belowground configurations. While several testing programs have been previously conducted, these earlier validation efforts did not capture the effects of elevated helium pressures or accurately portray the external convection of aboveground and belowground canistered dry cask systems. Thus the enhanced performance of modern dry storage casks cannot be fully validated using previous studies.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Machine learning assisted modeling of mixing timescale for LES/PDF of high-Karlovitz turbulent premixed combustion

Accurate modeling of mixing in the transported probability density function (PDF) method remains a great challenge, especially for turbulent premixed combustion under extreme conditions such as high Karlovitz number Ka. Recently, a power-law based mixing timescale model was developed for the large-eddy simulations (LES)/PDF modeling of high-Ka number turbulent premixed flames. It is found in this work that the power-law mixing timescale model is highly sensitive to the model parameters. It is thus critically needed to develop accurate calibration of these model parameters. The empirical specification of the model parameters developed in Zhang et. al. is found to be inadequate for accurate modeling of the mixing timescale. Here, machine learning is introduced as an attractive alternative in this work for the specification of the model parameters. A high-Ka number DNS jet flame is used as the training and validation of the machine learning models. The choices of the input parameters are discussed and compared for the machine learning models. The effect of differential molecular diffusion on mixing is examined by including the effect of the Lewis number in the training of the machine learning models. The performance of different machine learning algorithms is compared for the specification of the mixing model parameters. Overall, excellent performance of the machine learning models is observed for assisting the mixing modeling. The feasibility, interpretability, applicability, generality, and portability of using machine learning are discussed in general to provide a perspective on applying data-driven machine learning for turbulent combustion modeling studies.

42 ENGINEERING↗
DOE OSTI2018

Study of a Single-Power Two-Circuit ESR Process with Current-Carrying Mold: Mathematical Simulation of the Process and Experimental Verification

A novel single-power two-circuit ESR process (ESR-STCCM) with current-carrying mold has been investigated via numerical simulation and experimental research in this paper. A 2D quasi-steady-state mathematical model is developed to describe ESR-STCCM. The electromagnetic field, flow field, slag pool temperature distribution, and the shape of a molten steel pool in ESR-STCCM have been investigated by FLUENT software as well as user-defined functions (UDF). The results indicate that ESR-STCCM is different from the conventional ESR process. The maximum electromagnetic force, current density, Joule heat, and slag pool flow velocity are located in the lower part of the conductor in the ESR-STCCM process. The direction of the maximum electromagnetic force inclines upward. There are two distinct vortices in the slag pool. The larger swirl rotates counterclockwise near the conductor, with a value of 0.0263 m s{sup −1} due to the interaction of the electromagnetic force and gravity. The maximum temperature of the slag pool is 2070 K (1797 °C) and is located in the center of the swirl with a filling ratio of 0.6 and a 20 mm electrode immersion depth. The depth of a molten steel pool is shallower, which is conducive to improving solidification quality. In addition, the filling ratio of 0.6 is conducive to controlling steel solidification quality. Some experiments have been done, and the numerical model is confirmed by experimental results.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Nanoparticle Injection Technology for Remediating Leaks of CO 2 Storage Formation (Final Report)

Leakage of a well can occur through the wellbore, the annulus between well tubing and steel casing, or on the outside of the casing. The overall goal of the project is to develop a new technology that can repair or reduce the leakage of wellbore. The new technology is based on simultaneously injecting nanoparticles electrochemically and removing some of the harmful ions in the system so the distressed cement materials can be densified and the corrosion risk of steel casing can be reduced. There were two objectives. One is the development of the injection technology for leakage repair and the other is the development of new numerical simulation models that can simulate and predict the performance of the new repair technology. An electro-migration testing unit and a in-house wellbore testing system were developed; healing agents with nanoparticles were selected, numerical models were developed to simulate the electro-migration technology, and the effectiveness of the new technology were evaluated.

42 ENGINEERING↗
DOE OSTI2025

Nanoscale engineering and dynamic stabilization of mesoscopic spin textures

Thermalization, while ubiquitous in physics, has traditionally been viewed as an obstacle to be mitigated. In contrast, we demonstrate here the use of thermalization in the generation, control, and readout of “shell-like” spin textures with interacting 13 C nuclear spins in diamond, wherein spins are polarized oppositely on either side of a critical radius. The textures span several nanometers and encompass many hundred spins; they are created and interrogated without manipulating the nuclear spins individually. Long-time stabilization is achieved via prethermalization to a Floquet-engineered Hamiltonian under the electronic gradient field: The texture is therefore metastable and robust against spin diffusion. This enables the state to endure over multiple minutes before it decays. Our work on spin-state engineering paves the way for applications in quantum simulation and nanoscale imaging.

42 ENGINEERING↗
DOE OSTI2023

Direct Numerical Simulation of Partial Fuel Stratification Assisted Lean Premixed Combustion for Assessment of Hybrid G-Equation/Well-Stirred Reactor Model

Partial fuel stratification (PFS) is a promising fuel injection strategy to stabilize lean premixed combustion in spark-ignition (SI) engines. PFS creates a locally stratified mixture by injecting a fraction of the fuel, just before spark timing, into the engine cylinder containing homogeneous lean fuel/air mixture. Further, this locally stratified mixture, when ignited, results in complex flame structure and propagation modes similar to partially premixed flames and allows for faster and more stable flame propagation than a homogeneous lean mixture. This study focuses on understanding the detailed flame structures associated with PFS-assisted lean premixed combustion. First, a two-dimensional direct numerical simulation (DNS) is performed using detailed fuel chemistry, experimental pressure trace, and realistic initial conditions mapped from a prior engine large-eddy simulation (LES), replicating practical lean SI operating conditions. DNS results suggest that the conventional triple flame structure is prevalent during the initial stage of flame kernel growth. Both premixed and nonpremixed combustion modes are present with the premixed mode contributing dominantly to the total heat release. Detailed analysis further reveals the effects of flame stretch and fuel pyrolysis on flame displacement speed. Based on the DNS findings, the accuracy of a hybrid G-equation/well-stirred reactor (WSR) combustion model is assessed for the PFS-assisted lean operation in the LES context. It is found that the G-equation model qualitatively captures the premixed branches of the triple flame, while the WSR model predicts the nonpremixed branch of the triple flame. Finally, potential needs for improvements to the hybrid G-equation/WSR modeling approach are discussed.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2017

Evaluation of Rapid for a UNF cask benchmark problem

To ensure criticality safety and maintain a high level of security of Used Nuclear Fuel (UNF), it is necessary to perform detailed Particle Transport simulations. The commonly used technique is the 3-D, continuous energy Monte Carlo method that is time-consuming and, for eigenvalue calculations, may have difficulties with the source convergence. Wenner and Haghighat and Dufek demonstrated that in criticality problems with high dominance ratio and/or under-sampling possibilities, the fission matrix (FM) methodology should be used in order to be able to obtain unbiased solution. The RAPID (Real-time Analysis for Particle transport and In-situ Detection) code system is based on the Multi-stage Response-function Transport (MRT) methodology, and uses the FM formulation to determine eigenvalue, subcritical multiplication and fission density distribution of nuclear systems. This paper uses a benchmark problem similar to the GBC-32 Cask to evaluate the accuracy and performance of RAPID against the results of a reference MCNP solution. We performed a detailed sensitivity analysis on a single-assembly model in order to determine the MCNP eigenvalue parameters necessary to obtain a reliable reference solution. These results provided a guide for the selection of an appropriate set of parameters for the eigenvalue calculation for the whole cask. The paper is organized as follows: initially, the RAPID code system is briefly discussed; then, the benchmark problem is introduced. The methodology used to analyze results is proposed and some of the results of different cases are analyzed to identify the best set of eigenvalue parameters for a reference MCNP calculation. The selected case is then compared with the RAPID calculation for both a single assembly and a full cask model. In the end some conclusions and future work perspectives are drawn.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2018

Direct numerical simulation of flame stabilization assisted by autoignition in a reheat gas turbine combustor

A three-dimensional direct numerical simulation (DNS) is performed for a turbulent hydrogen-air flame, represented with detailed chemistry, stabilized in a model gas-turbine combustor. The combustor geometry consists of a mixing duct followed by a sudden expansion and a combustion chamber, which represents a geometrically simplified version of Ansaldo Energia’s GT26/GT36 sequential combustor design. In this configuration, a very lean blend of hydrogen and vitiated air is prepared in the mixing duct and convected into the combustion chamber, where the residence time from the inlet of the mixing duct to the combustion chamber is designed to coincide with the ignition delay time of the mixture. The results show that when the flame is stabilized at its design position, combustion occurs due to both autoignition and flame propagation (deflagration) modes at different locations within the combustion chamber. A chemical explosive mode analysis (CEMA) reveals that most of the fuel is consumed due to autoignition in the bulk-flow along the centerline of the combustor, and lower amounts of fuel are consumed by flame propagation near the corners of the sudden expansion, where the unburnt temperature is reduced by the thermal wall boundary layers. An unstable operating condition is also identified, wherein periodic auto-ignition events occur within the mixing duct. These events appear upstream of the intended stabilization position, due to positive temperature fluctuations induced by pressure waves originating from within the combustion chamber. This DNS investigation represents the initial step of a comprehensive research effort aimed at gaining detailed physical insight into the rate-limiting processes that govern the sequential combustor behavior and avoid the insurgence of the off-design auto-ignition events.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2017

Magnetar-powered Supernovae in Two Dimensions. II. Broad-line Supernovae Ic

Nascent neutron stars (NSs) with millisecond periods and magnetic fields in excess of 10{sup 16} Gauss can drive highly energetic and asymmetric explosions known as magnetar-powered supernovae. These exotic explosions are one theoretical interpretation for supernovae Ic-BL, which are sometimes associated with long gamma-ray bursts. Twisted magnetic field lines extract the rotational energy of the NS and release it as a disk wind or a jet with energies greater than 10{sup 52} erg over ∼20 s. What fraction of the energy of the central engine go into the wind and the jet remain unclear. We have performed two-dimensional hydrodynamical simulations of magnetar-powered supernovae (SNe) driven by disk winds and jets with the CASTRO code to investigate the effect of the central engine on nucleosynthetic yields, mixing, and light curves. We find that these explosions synthesize less than 0.05 M{sub ⊙} of {sup 56}Ni and that this mass is not very sensitive to central engine type. The morphology of the explosion can provide a powerful diagnostic of the properties of the central engine. In the absence of a circumstellar medium, these events are not very luminous, with peak bolometric magnitudes of M{sub b}∼−16.5 due to low {sup 56}Ni production.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Replicating HCCI-like autoignition behavior: What gasoline surrogate fidelity is needed?

This work seeks to characterize the fidelity needed in a gasoline surrogate with the intent to replicate the complex autoignition behavior exhibited within advanced combustion engines, and specifically Homogeneous Charge Compression Ignition (HCCI). A low-temperature gasoline combustion (LGTC) engine operating in HCCI mode and a rapid compression machine (RCM) are utilized to experimentally quantify fuel reactivity, through autoignition and preliminary heat release characteristics. Fuels considered include a research grade E10 U.S. gasoline (RD5-87), three multi-component surrogates (PACE-1, PACE-8, PACE-20), and a binary surrogate (PRF88.4). Each fuel was studied at lean/HCCI-like conditions covering a wide range of temperatures and pressures that are representative of naturally aspirated to high boost engine operation. Detailed chemical kinetic modeling is also undertaken using a recently updated gasoline surrogate kinetic model to simulate the RCM experiments and to provide chemical insight into surrogate-to-surrogate differences. The LGTC engine experiments demonstrate nearly identical reactivity between PACE-20 and RD5-87 across studied conditions, while faster phasing is seen for both PACE-1 and PACE-8 due to their stronger intermediate- and low-temperature heat release (ITHR/LTHR) at naturally aspirated and boosted conditions, respectively. The RCM experiments reveal typical low-temperature, negative temperature coefficient (NTC) and intermediate-temperature autoignition behaviors at all pressure conditions for RD5-87, which are qualitatively reproduced by all surrogates. Quantitative discrepancies in both autoignition and preliminary heat release are observed for all surrogates, while their ability to replicate RD5-87 autoignition behavior follows the order of PACE-20 > PACE-1 > PACE-8 > PRF88.4. Excellent mapping is obtained between the LGTC engine and the RCM, where the engine pressure-time trajectories can be characterized by the regimes represented by the RCM autoignition isopleths. The kinetic model performs commendably when simulating both autoignition and preliminary heat release of PACE-20, while typically overpredicting ignition delay times for PACE-1, PACE-8 and PRF88.4 at high-pressure and low-temperature/NTC conditions. Sensitivity and rate of production (ROP) analyses highlight surrogate-to-surrogate differences in the governing chemical kinetics where n-pentane initiates rapid OH branching at a faster rate and an earlier timing for PACE-20 than iso-pentane does for PACE-1 and PACE-8, making it computationally more reactive than the other surrogates. The current study highlights the need to include non-standardized properties, such as the lean/HCCI-like autoignition characteristics, in addition to ASTM properties (e.g., RON, MON) as metrics of fuel reactivity and targets to be matched when formulating high-fidelity surrogates that fully capture gasoline advanced combustion behavior such as HCCI-like autoignition.

42 ENGINEERING↗
DOE OSTI2020

Driven-dissipative quantum mechanics on a lattice: Simulating a fermionic reservoir on a quantum computer

The driven-dissipative many-body problem remains one of the most challenging unsolved problems in quantum mechanics. The advent of quantum computers may provide a unique platform for efficiently simulating such driven-dissipative systems. But, there are many choices for how one can engineer the reservoir. One can simply employ ancilla qubits to act as a reservoir and then digitally simulate them via algorithmic cooling. A more attractive approach, which allows one to simulate an infinite reservoir, is to integrate out the bath degrees of freedom and describe the driven-dissipative system via a master equation, that can also be simulated on a quantum computer. In this work, we consider the particular case of noninteracting electrons on a lattice driven by an electric field and coupled to a fermionic thermostat. Then, we provide two different quantum circuits: the first one reconstructs the full dynamics of the system using Trotter steps, while the second one dissipatively prepares the final nonequilibrium steady state in a single step. We run both circuits on the IBM quantum experience. For circuit (i), we achieved up to five Trotter steps. When partial resets become available on quantum computers, we expect that the maximum simulation time can be significantly increased. Lastly, the methods developed here suggest generalizations that can be applied to simulating interacting driven-dissipative systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Determination of the Most Likely Wear Zones of Pipe Bends in the Pulsed Pneumatic Conveyor at the RT-1 Plant

The results of an experimental investigation and numerical modeling of the most likely high-wear zones in 90° pipeline bends in the pulsed pneumatic conveyor system at the RT-1 plant of the Mayak Industrial Association are presented. The experimental studies were conducted on a model of a pulsed pneumatic conveyor system; the numerical modeling of the gas-dynamic flow in the pipeline channel was performed in a model of a two-dimensional non-stationary motion of a compressible medium. It was found on the basis of the experimental results obtained using the model and numerical modeling in the homogeneous approximation that local regions with high wear of walls in the 90° pipe bends in the pulsed pneumatic conveyor do indeed exist and their experimentally found positions agree with the calculations.

42 ENGINEERING↗
DOE OSTI2019

Modeling of the Heat-Affected and Thermomechanically Affected Zones in a Ti-6Al-4V Inertia Friction Weld

Inertia friction welding has been used across the aerospace, automotive, and power-generation industries for the fabrication of complex axisymmetric components for over forty years. The process involves one axisymmetric piece being held stationary and another piece being brought into contact set to rotate about its axis of symmetry by a flywheel with the system under an applied load across the joint. Plasticization at the joint interface through the frictional heating sees the two pieces bond together. The titanium alloy Ti-6Al-4V has been widely studied for inertia welding applications. A successful selection of processing parameters (flywheel energy and mass, applied load) allows an inertia welding process which produces a very high-integrity weld, with a minimal heat-affected zone (HAZ) and thermomechanically affected zone (TMAZ), formed as a narrow band at the interface and extending further into the material. The width of this narrow band of heated material is dependent upon the process parameters used. A series of experimental inertia friction welds were performed using Ti-6Al-4V, and a finite element (FE) modeling framework was developed using the FE code Deform in order to predict the widths of the HAZ and TMAZ at the weld interface. The experimentally observed HAZ boundaries were correlated with the thermal fields from the FE model, while TMAZ boundaries were correlated with the Von Mises plastic strain fields.

42 ENGINEERING↗
DOE OSTI2018

Blowing polar skyrmion bubbles in oxide superlattices

Particle-like topological structures such as skyrmions and vortices have garnered ever-increasing interests due to their rich physical insights and potential broad applications in spintronics. Here, in this work, we discover the reversible switching between polar skyrmion bubbles and ordered vortex arrays in ferroelectric superlattices under an electric field, reminiscent of the Plateau-Raleigh instability in fluid mechanics. An electric field phase diagram is constructed, showing a wide stability window for the observed polar skyrmions. A “volcano”-like pontryagin density distribution is formed, indicating the formation of a smooth circular skyrmion. The topological charge Q at different applied field is calculated, verifying the field-driven topological transition between Q = 0 and Q = ±1 states. This study is a demonstration for the computational design of field-induced topological phase transitions, giving promise for the design of next-generation nanoelectronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Evaluation of enzymatic depolymerization of PET, PTT, and PBT polyesters

Millions of tons of waste polyester plastics, including polyethylene terephthalate (PET), polytrimethylene terephthalate (PTT), and polybutylene terephthalate (PBT), end up in the environment as soil and water contaminants. Recent advances in enzymatic polyester degradation have motivated researchers toward biorecycling of plastic wastes. Leaf-branch compost cutinase (LCC) enzymes have been proven to be effective in the biodegradation of PET. This study focuses on enzymatic depolymerization of waste PET, PTT, and PBT materials by using an ICCG variant of LCC (LCC ICCG ) produced from Escherichia coli BL21(DE3). The degradation efficiency of the polyesters was determined by the monomer terephthalic acid (TPA) released from the depolymerization reaction. It was found that the most efficient depolymerization was achieved for PET, followed by PTT and PBT. A kinetic model based on Langmuir adsorption and the Michaelis-Menten equation was developed to describe the enzymatic depolymerization of PET, PTT, and PBT with various enzyme and substrate loadings. The model simulation results revealed that the LCCICCG enzyme loading should be linearly increased as the work capacity of the polyester substrate increases. A specific enzyme loading of 0.91 mg/g PET is suggested to achieve 90% depolymerization of PET within three days. Furthermore, the experimental data and model simulation results can be used to help further engineer the enzyme and process to achieve a complete biodegradation of polyester wastes at large scale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗