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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

RX-PSA

This code is built on top of the ML-PSA utility and implements the ability to implement a multi-fidelity parallel simulated annealing optimization for a range of engineering problems. RX-PSA provides specializations to interact with modern nuclear reactor codes such as VERA to perform assembly and core optimization in a robust way. RX-PSA also provides LWR specific objective functions and constraints to provide a straightforward user interface that reactor designers are familiar with.

Collins, BenjaminS.↗
DOE OSTI2022

DPC Direct Disposal Postclosure Thermal Modeling

Performance of geologic radioactive waste repositories depends on near-field and far-field processes, including km-scale flow and transport in engineered and natural barriers, that may require simulations of up to 1 M years of regulatory period. For a relatively short time span (less than 1000 years), the thermohydro-mechanical-chemical (THMC) coupled processes caused by heat from the waste package will influence near-field multiphase flow, chemical/reactive transport, and mechanical behaviors in the repository system. This study integrates the heat-driven perturbations in thermo-hydro-mechanical characteristics into thermo-hydro-chemical simulations using PFLOTRAN to reduce dimensionality and improve computational efficiency by implementing functions of stress-dependent permeability and saturation-temperature-dependent thermal conductivity. These process couplings are developed for spent nuclear fuel in dual-purpose canisters in two different hypothetical repositories: a shale repository and a salt repository.

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

SISGR: Chemomechanics of Far-From-Equilibrium Interfaces (COFFEI)

Portable, reliable, and deployable devices for energy storage and conversion require fundamental changes in design of solid-state composites comprising ceramics and metals. These materials comprise the electrodes and electrolytes of next-generation solid-oxide fuel cells and all solid-state batteries, forming solid-state functional composites. The advent of solid-state batteries – which replace liquid electrolytes with solid electrolytes capable of lithium ion transport for reliable energy storage in portable batteries – and the increased demand space for all solid-state fuel cells capable of oxygen reduction at intermediate temperatures remain important challenges for improved material stability and decreased system cost. However, little is understood about three fundamental facets of materials that enable such solid-state energy applications. First, how do such materials deform, fracture, or delaminate under operando conditions? Second, how does such mechanical deformation limit or facilitate electronic and ionic transport within and across such material interfaces? Third, how we can predictably design interface-rich composites to engineer both structural and electrochemical stability? This COFFEI Group comprised expertise from Materials Science & Engineering and Nuclear Science & Engineering to integrate unique in situ experiments, simulations, and fabricated interfaces that address these fundamental questions in solid-state interfaces of nanoscale composites that will guide solid-state electrochemistry, transport kinetics, and mechanical deformation for nonstoichiometric materials that enable such applications. In particular, we built on COFFEI’s understanding of chemomechanical coupling among defect concentrations, ionic transport, electron transport, and stored elastic energy that is particularly acute in the far-from-equilibrium conditions typical of energy device applications. By tailoring our focus to solid-state interfaces, we addressed these important issues by (a) developing and applying advanced in situ and ex situ characterization tools to characterize model materials and interfaces synthesized with molecular-level control, under both laboratory-controlled and extreme environments representative of energy device operation; and (b) employing computational modeling and simulation frameworks to predict transport mechanisms, reactivity and stability of these model materials and interfaces under significant chemical strains typical of energy device operation. Recent progress provided insights to additional materials systems and electrochemomechanical fatigue and fracture that were not fully envisioned when the program was initiated. Specifically, in the final three years of COFFEI we pursued two integrated thrusts, with complementary focus. Thrust I focused on failure-resistant electrochemomechanical composites, while Thrust II focused on strain-modulated conductivity and reactivity across interfaces. In contrast to our initial COFFEI focus, these thrusts concentrated wholly on solid-state material interfacial interactions and included greater integration of multiscale visualization including in situ electron microscopy of strained structures/interactions and mesoscale simulations. Successful development of functionally superior and long-lived battery and fuel cell systems and stress adaptable oxides requires a deeper, fundamental understanding of the coupling among the historically important subfields of solid-state electrochemistry, transport kinetics, and mechanical deformation for nonstoichiometric metal oxide electrodes. In this program, the understanding and the application of chemomechanical coupling of defect concentrations, ionic transport, electro-catalytic activity and stored elastic energy, particularly acute in the far-from-equilibrium conditions typical of energy device applications, are being refined and implications for device operation clarified, including for miniaturized solid-state batteries and fuel cells.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Simulating target tests at AP0 facility for Mu2e

The predicted production of muons and energy deposition for the production target of the Mu2e experiment are dependent on the documentation of cross-sections values, the properties of the target’s material, etc. Hence, making real tests on similar targets before its use on Mu2e will be a good indicator of the accuracy of the simulations done for the experiment. The main goal of the project is to simulate a target test at the AP0 Target Facility, previously used to create the muon beam for the Muon g-2 experiment. Results from these simulations will be used as input to engineering models of the target to produce predicted stress and surface temperature distributions, parameters that will later be measured in future beam runs to validate underlying cross section models.

Delgado, Amii Matamoros↗
DOE OSTI2026

Waveform Simulation Framework: User Manual with Tutorials

This manuscript describes the Waveform Simulation Framework (WSF), a Python-based framework that provides a unified, programmable interface for generating synthetic seismograms for applications such as seismic array design, method development, and special event analysis. WSF standardizes how users define sources, receivers, and velocity models while abstracting simulator-specific configuration details, enabling workflows that are largely independent of the underlying numerical engine. The document provides installation guidance and tutorial-driven examples for three WSF simulator wrappers—WSF PyFK, WSF SW4, and WSF SPECFEM2D—illustrating end-to-end workflows from forward waveform simulation to common post-processing tasks (e.g., visualization and backprojection) using consistent data products (e.g., ObsPy Stream objects and SAC files).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Fabricated devices for performing bacterial-fungal interaction experiments across scales

Diverse and complex microbiomes are found in virtually every environment on Earth. Bacteria and fungi often co-dominate environmental microbiomes, and there is growing recognition that bacterial-fungal interactions (BFI) have significant impacts on the functioning of their associated microbiomes, environments, and hosts. Investigating BFI in vitro remains a challenge, particularly when attempting to examine interactions at multiple scales of system complexity. Fabricated devices can provide control over both biotic composition and abiotic factors within an experiment to enable the characterization of diverse BFI phenotypes such as modulation of growth rate, production of biomolecules, and alterations to physical movements. Engineered devices ranging from microfluidic chips to simulated rhizosphere systems have been and will continue to be invaluable to BFI research, and it is anticipated that such devices will continue to be developed for diverse applications in the field. This will allow researchers to address specific questions regarding the nature of BFI and how they impact larger microbiome and environmental processes such as biogeochemical cycles, plant productivity, and overall ecosystem resilience. Devices that are currently used for experimental investigations of bacteria, fungi, and BFI are discussed herein along with some of the associated challenges and several recommendations for future device design and applications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Numerical Modeling of Supersonic Combustion in a Non- Premixed Rotating Detonation Engine

A computational fluid dynamics (CFD) model was developed to simulate a representative non-premixed RDE combustor. Unsteady Reynolds-Averaged Navier-Stokes (RANS) simulations were performed for the full combustor geometry (including the separate fuel and air injection ports), with hydrogen as the fuel and air as the oxidizer. Adaptive mesh refinement (AMR) was incorporated to achieve a good trade-off between model accuracy and computational expense. A finite-rate chemistry model was employed along with a 10-species detailed kinetic mechanism to describe the H2-Air combustion chemistry. Two different operating conditions were simulated, corresponding to the same global equivalence ratio of unity but different fuel and air mass flow rates. Numerical results were compared against available experimental data on detonation wave frequency/height, reactant fill height, oblique shock angle and axial pressure distribution in the channel. The CFD model predicted the sensitivity of these wave characteristics to operating conditions with good accuracy, both qualitatively and quantitatively. The present CFD model offers a potential capability to perform rapid design space exploration and/or performance optimization studies for realistic full-scale RDE configurations.

Pal, Pinaki↗
DOE OSTI2018

Auto-ignitive deflagration of CH4 blended DME/Air mixtures with composition stratification

Front propagation speeds from one dimensional simulations of DME/Air and DME/CH4/Air mixtures under engine relevant conditions are presented with reduced kinetics and transport. Different time scales of monochromatic inhomogeneities in DME concentration and varying DME-CH4 blending ratios are explored to understand its subsequent effect on the propagation speed as well as the ensuing transition from deflagration to spontaneous propagation. It is found that for a given timescale, the instantaneous peak propagation speed achieved by a front as well as its overall variation is significantly affected by the level of CH4 concentration in the mixture. Specifically, at relatively smaller time scales, the magnitude of variation in the instantaneous propagation speed is found to subside as the level of CH4 concentration in the mixture is increased. An opposite trend is observed at comparatively larger time scales, wherein, a rapid change in the instantaneous propagation speed is observed with an increase in the level of CH4 concentration. This study provides evidence as to how the propagation speed of a single fuel flame varies compared to a dual fuel flame under the influence of composition stratification.

Desai, Swapnil S.↗
DOE OSTI2021

Thermal and Mechanical Energy Performance Analysis of Closed-loop Systems in Hot-Dry-Rock and Hot-Wet-Rock Reservoirs

To understand the potential and limitations for recovering thermal and mechanical energy from closed-loop geothermal systems a collaborative study is underway that will investigate an array of system configurations, working fluids, geothermal reservoir characteristics, operational periods, and heat transfer enhancements. Closed-loop geothermal systems are distinguished from hydrothermal or enhanced geothermal systems (EGS) in that the working fluid only circulates through drilled boreholes. Principal objectives of this study are to determine upper limits for thermal and mechanical energy recovery and optimal operational and configuration parameters for each scenario. Teams of scientists and engineers are applying a suite of numerical simulation and analytical tools to model the heat recovery from closed-loop geothermal systems, and then optimizing operational and configuration parameters to maximize the thermal and mechanical energy recovery. Results from the suite of numerical simulators and analytical tools, such as outlet and inlet states and temperature profiles in the geothermal reservoir over time are intercompared to increase confidence in the analysis. This paper documents the study findings for closed-loop systems in hot-dry-rock and hot-wet-rock reservoirs, where water is the working fluid. The characteristics of the hot-dry-rock reservoir were based on the U.S. Department of Energy's Utah Frontier Observatory for Research in Geothermal Energy (FORGE) site, near Milford Utah. Two objective functions are defined to optimize the operational and configuration parameters of the system, one each for the recovery of mechanical and thermal energy over the period of operation. For both objective functions, a surface plant thermal to mechanical energy conversion factor and an energy drilling cost is required. In keeping with the study objectives the surface plant conversion factor is determined from a second-law of thermodynamics analysis of a generic binary plant, and drilling costs are based on those from the Utah FORGE site and current national electrical costs.

closed-loop geothermal systems↗
DOE OSTI2021

Thermal and Mechanical Energy Performance Analysis of Closed-loop Systems in Hot-Dry-Rock and Hot-Wet-Rock Reservoirs

To understand the potential and limitations for recovering thermal and mechanical energy from closed-loop geothermal systems a collaborative study is underway that will investigate an array of system configurations, working fluids, geothermal reservoir characteristics, operational periods, and heat transfer enhancements. Closed-loop geothermal systems are distinguished from hydrothermal or enhanced geothermal systems (EGS) in that the working fluid only circulates through drilled boreholes. Principal objectives of this study are to determine upper limits for thermal and mechanical energy recovery and optimal operational and configuration parameters for each scenario. Teams of scientists and engineers are applying a suite of numerical simulation and analytical tools to model the heat recovery from closed-loop geothermal systems, and then optimizing operational and configuration parameters to maximize the thermal and mechanical energy recovery. Results from the suite of numerical simulators and analytical tools, such as outlet and inlet states and temperature profiles in the geothermal reservoir over time are intercompared to increase confidence in the analysis. This paper documents the study findings for closed-loop systems in hot-dry-rock and hot-wet-rock reservoirs, where water is the working fluid. The characteristics of the hot-dry-rock reservoir were based on the U.S. Department of Energy’s Utah Frontier Observatory for Research in Geothermal Energy (FORGE) site, near Milford Utah. Two objective functions are defined to optimize the operational and configuration parameters of the system, one each for the recovery of mechanical and thermal energy over the period of operation. For both objective functions, a surface plant thermal to mechanical energy conversion factor and an energy drilling cost is required. In keeping with the study objectives the surface plant conversion factor is determined from a second-law of thermodynamics analysis of a generic binary plant, and drilling costs are based on those from the Utah FORGE site and current national electrical costs.

Closed-loop geothermal systems, hot-dry-rock, hot-↗
DOE OSTI2022

Novel Solver Algorithms for Nearly Singular Linear Systems Arising in Combustion Modelling

Direct Numerical Simulations of realistic combustion devices are extremely challenging due to the wide separation of scales in the simulation, for example an internal combustion (IC) engine chamber, and the flame thickness of a high-pressure flame. The PeleLMeX solver uses adaptive mesh refinement (AMR) to evolve multi-species reacting flows in the low Mach number limit at the Exascale and relies on an embedded boundary (EB) approach to represent complex geometries. In that framework, the EB geometries often give rise to very small cut-cells along the boundary, which translate into extreme ill-conditioning of the pressure-projection, with eigenvalues that span 15-16 orders of magnitude. In this talk, we focus on the case of a typical IC piston bowl geometry for which we present on a novel approach towards solving these nearly singular linear systems with ILU-based, C-AMG smoothers on massively parallel architectures. In particular, we use scaling and equilibration algorithms to handle the non-normality of the upper triangular factors. This enables us to approximate the highly sequential triangular solve algorithm, embedded in the AMG smoothing-solve phase, with Jacobi iterations. This approximation can be written as a convergent Neumann series whose terms are composed of highly parallel sparse matrix vector multiplications. The result is an algorithm that substantially decreases setup and solve time, compared to state-of-the-art, for these challenging linear systems.

combustion modelling↗
DOE OSTI2018

Numerical Simulation of Recirculating Flow and Physical Model of Slag–Metal Behavior in an RH Reactor: Application to Desulfurization

Computational fluid dynamics (CFD) techniques and a 1:7.5 physical model of an RH degasser have been used to evaluate the flow of gas and metal inside an RH reactor for vacuum degassing of liquid steel. The effect of gas injection on the gas spatial distribution, steel circulation rate and flow field inside the ladle, snorkels and vacuum chamber have been assessed. N-pentane oil was employed to evaluate the average residence time as well as the slag droplet size distribution. The predicted radial gas distribution and liquid circulation rate have been validated against experimental data from a physical model. The results with incorporation of the virtual mass force coefficient of 0.25 and the turbulence dispersion force showed better predictions of gas distribution in the up-snorkel as well as circulation rate. Full-scale simulations were performed, and the predicted circulation rate was significantly affected by the argon bubble expansion. Data from these simulations were used to analyze the degree of desulfurization performed by the addition of desulfurizing agents inside the vacuum chamber. A model of the kinetics of desulfurization based on the results from the physical model and CFD simulation and on slag dispersion inside liquid steel yields degrees of desulfurization similar to the industrial trials reported in the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Investigation of Cycle-to-Cycle Variations in Internal Combustion Engine Using Proper Orthogonal Decomposition

The understanding, modeling and control of the cycle-to-cycle variation (CCV) in the modern internal combustion engine (ICE) is a key scientific challenge to achieve stable engine operation. High CCV in the engine combustion chamber may contribute to partial burn, misfire and knock, which adversely affects the engine performance and may potentially damage the engine. The objective of the current study is to leverage high-fidelity numerical simulations to improve the understanding of the causes of CCV. Using the massively parallel code, Nek5000, multi-cycle, wall-resolved large-eddy simulations (LES) were performed for the General Motors (GM), Transparent Combustion Chamber (TCC-III) optical engine under motored operating conditions. Further, the large-scale structures of the in-cylinder flow were investigated using a triple proper orthogonal decomposition (POD) technique to explore the characteristics of different parts of the flow and their contributions to CCV. The kinetic energy of the subset of flow structures were determined and correlated between the intake and compression strokes. The insights from the analysis of the large-scale flow structures will be used to assist the development of improved engine designs with reduced CCV and enhance the engine performance.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Model Reduction by Generalized Falk Method for Efficient Field-Circuit Simulations

The Generalized Falk Method (GFM) for coordinate transformation, together with two model-reduction strategies based on this method, are presented for efficient coupled field-circuit simulations. Each model-reduction strategy is based on a decision to retain specific linearly-independent vectors, called trial vectors, to construct a vector basis for coordinate transformation. The reduced-order models are guaranteed to be stable and passive since the GFM is a congruence transformation of originally symmetric positive definite systems. We also show that, unlike the Pad´e-via-Lanczos (PVL) method, the GFM does not generate unstable positive poles while reducing the order of circuit problems. Further, the proposed GFM is also faster when compared to methods of the type Lanczos (or Krylov) that are already widely used in circuit simulations for electrothermal and electromagnetic problems. The concept of response participation factors is introduced for the selection of the trial vectors in the proposed model-reduction methods. Further, we present methods to develop simple equivalent circuit networks for the field component of the overall field-circuit system. The implementation of these equivalent circuit networks in circuit simulators is discussed. With the proposed model-reduction strategies, significant improvement on the efficiency of the generalized Falk method is illustrated for coupled field-circuit problems.

42 ENGINEERING↗
DOE OSTI2018

Predicting the 3D fatigue crack growth rate of small cracks using multimodal data via Bayesian networks: In-situ experiments and crystal plasticity simulations

Small crack propagation accounts for most of the fatigue life of engineering structures subject to high cycle fatigue loading conditions. Determining the fatigue crack growth rate of small cracks propagating into polycrystalline engineering alloys is critical to improving fatigue life predictions, thus lowering cost and increasing safety. In this work, cycle-by-cycle data of a small crack propagating in a beta metastable titanium alloy is available via phase and diffraction contrast tomography. Crystal plasticity simulations are used to supplement experimental data regarding the micromechanical fields ahead of the crack tip. Experimental and numerical results are combined into a multimodal dataset and sampled utilizing a non-local data mining procedure. Furthermore, to capture the propensity of body-centered cubic metals to deform according to the pencil-glide model, a non-local driving force is postulated. The proposed driving force serves as the basis to construct a data-driven probabilistic crack propagation framework using Bayesian networks as building blocks. The spatial correlation between the postulated driving force and experimental observations is obtained by analyzing the results of the proposed framework. Results show that the above correlation increases proportionally to the distance from the crack front until the edge of the plastic zone. Moreover, the predictions of the propagation framework show good agreement with experimental observations. Finally, we studied the interaction of a small crack with grain boundaries (GBs) utilizing various slip transmission criteria, revealing the tendency of a crack to cross a GB by propagating along the slip directions minimizing the residual Burgers vector within the GB.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Full-System Linearization for Floating Offshore Wind Turbines in OpenFAST

The wind engineering community relies on multiphysics engineering software to run nonlinear time-domain simulations (e.g., for design-standards-based loads analysis). Although most physics involved in wind energy are nonlinear, linearization of the underlying nonlinear system equations is often advantageous to understand the system properties and exploit well-established methods and tools for analyzing linear systems. This paper presents the development of the new linearization functionality of the open-source engineering tool OpenFAST for floating offshore wind turbines, as well as the concepts and mathematical background needed to understand and apply it.

17 WIND ENERGY↗
DOE OSTI2023

Preliminary Nuclear Containment Vessel Modeling for Multi-Hazard Probabilistic Risk Assessment under Seismic Hazards and Concrete Degradation

The current practice for natural phenomena hazards (NPH) risk assessment of nuclear facilities is to compute the risk for each hazard independently and then compound the total risk as a combination of single hazard risks. This state of practice does not consider correlations between hazards and the cascading impacts to structures, systems, and components (SSCs), and could thus underestimate the NPH risk or overestimate the nuclear facility safety. Events such as the Fukushima Daiichi accident have highlighted the importance of multi-hazard risk considerations to nuclear power plants (NPPs) that quantify the cascading damage effects to SSCs in the risk models. Moreover, the current fleet of NPPs in the United States is aging; these NPPs are now expected to operate well beyond their initially planned design life. Aging-related deterioration can potentially decrease the capacity of critical structures such as containment vessels to withstand NPH. Such aging considerations may not be adequately accounted for by the current NPH risk assessment guidelines. This paper presents a preliminary modeling and simulation of a representative reinforced concrete containment vessel subjected to seismic mainshock and aftershock considering concrete degradation due to alkali silica reaction. The broader aim is to develop multi-hazard time-dependent fragility functions that could be subsequently used in the probabilistic risk assessment (PRA) model. The multi-hazard component comes into play due to the consideration of damage to the containment vessel under seismic loads and concrete degradation. Consideration of concrete degradation also brings into play the time-dependent nature of the containment vessel response. The response of the containment vessel under varying degrees of concrete degradation to seismic loads is investigated. The results presented are simulated using the Multi-hazard Analysis for STOchastic time-DOmaiN phenomena (MASTODON) software for seismic analysis and the Blackbear software for concrete degradation and damage modeling. Both software are open source and developed within the Multiphysics Object-Oriented Simulation Environment (MOOSE).

42 ENGINEERING↗
DOE OSTI2018

Full-System Linearization for Floating Offshore Wind Turbines in OpenFAST: Preprint

The wind engineering community relies on multiphysics engineering software to run nonlinear time-domain simulations, e.g. for design-standards-based loads analysis. Although most physics involved in wind energy are nonlinear, linearization of the underlying nonlinear system equations is often advantageous to understand the system properties and exploit well-established methods and tools for analyzing linear systems. This paper presents the development of the new linearization functionality of the open-source engineering tool OpenFAST for floating offshore wind turbines (FOWTs), as well as the concepts and mathematical background needed to understand and apply it correctly.

17 WIND ENERGY↗