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 1,711 records · Page 95

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Dose Measurement and Monte Carlo FLUKA simulations for the FE and FOE of the 27-ID HEX Beamline

The 27-ID High Energy Engineering X-ray Scattering (HEX) beamline, commissioned in 2024, operates within the Hard X-ray Scattering and Spectroscopy program at Brookhaven Lab’s National Synchrotron Light Source II (NSLS-II). The high energy spectrum at the HEX beamline extends up to several hundred keV, making it well suited for applications such as imaging large engineering samples, including fully assembled batteries, fuel cells, etc.

43 PARTICLE ACCELERATORS↗

DOE OSTI2019

High fidelity simulation and validation of crossflow through a tube bundle and the onset of vibration

Flow induced vibrations can be detrimental to all manner of engineering applications with fluid flow. Heat exchangers with bundles of thin tubes experiencing a crossflow are especially susceptible to the unexpected onset of vibration. Heat-exchanger design necessarily involves extensive modeling and testing to ensure that significant vibrations cannot develop for any expected flow conditions. A properly validated numerical simulation can work in conjunction with physical experiments to identify problematic vibrations and allow for rapid iteration of design at relatively low expense. A high-fidelity fluid–structure interaction code has been developed by fully coupling CFD LES/DNS code Nek5000 and CSM code Diablo. The coupled code is used to simulate crossflow through a tube bundle in a geometry recreated after a physical experiment. Validation against this experiment involves comparing the amplitude and frequency spectrum for three different flow velocities. The velocities compared straddle the onset of large magnitude vibrations. The simulation accurately captures the onset of the vibrations with a marginally lower predicted frequency of vibration. When the onset velocity is compensated for the difference in natural frequency, the simulation results closely match the experiment. Finally, the continuous nature of the simulation measurements helps illustrate the fluid mechanics behind the pin motion and reveals the front to back propagation of the vibrations as the fluid velocity increases.

42 ENGINEERING↗

DOE OSTI2022

Less can be more: Insights on the role of electrode microstructure in redox flow batteries from two-dimensional direct numerical simulations

Understanding how to structure a porous electrode to facilitate fluid, mass, and charge transport is key to enhancing the performance of electrochemical devices, such as fuel cells, electrolyzers, and redox flow batteries (RFBs). Here, using a parallel computational framework, direct numerical simulations are carried out on idealized porous electrode microstructures for RFBs. Strategies to improve an electrode design starting from a regular lattice are explored. By introducing vacancies in the ordered arrangement, it is possible to achieve higher voltage efficiency at a given current density, thanks to improved mixing of reactive species, despite reducing the total reactive surface. Careful engineering of the location of vacancies, resulting in a density gradient, outperforms disordered configurations. Our simulation framework is a new tool to explore transport phenomena in RFBs, and our findings suggest new ways to design performant electrodes.

25 ENERGY STORAGE↗

DOE OSTI2017

SimApp 1.8.x

System for converting DESIRE (Direct Executing Simulation in Real Time) simulation models to executable code models (C# / C++). These converted models are used as the basis for 2 products, one of which is the focus of this review: an engineering design interface for Fissile Solution Systems. ModelConverter.exe (C#) converts DESIRE model inputs to a C++ plug-in subclass. SimApp.exe (C# / C++) has a C++ dll that is a numerical integration engine that the Converted models plug into. The engineering design interface is a C# GUI that calls on the C++ dll to execute a particular converted model. The GUI allows the user to specify the length of a simulation, run the simulation, view data plots and data tables. Stability plots may also be produced and manipulated. The engineering design interface will be released as source code and executable; the model converter source is also intended for release.

Determan, JohnC.↗

DOE OSTI2017

SimApp 1.8.x

Determan, John C.↗

DOE OSTI2016

Co-Optimization of CO 2 -EOR and Storage Processes in Mature Oil Reservoirs

This article presents an optimization methodology for CO 2 enhanced oil recovery in partially depleted reservoirs. A field-scale compositional reservoir flow model was developed for assessing the performance history of an active CO 2 flood and for optimizing both oil production and CO 2 storage in the Farnsworth Unit (FWU), Ochiltree County, Texas. A geological framework model constructed from geophysical, geological, and engineering data acquired from the FWU was the basis for all reservoir simulations and the optimization method. An equation of state was calibrated with laboratory fluid analyses and subsequently used to predict the thermodynamic minimum miscible pressure (MMP). Initial history calibrations of primary, secondary and tertiary recovery were conducted as the basis for the study. After a good match was achieved, an optimization approach consisting of a proxy or surrogate model was constructed with a polynomial response surface method (PRSM). The PRSM utilized an objective function that maximized both oil recovery and CO 2 storage. Experimental design was used to link uncertain parameters to the objective function. Control variables considered in this study included: water alternating gas cycle and ratio, production rates and bottom-hole pressure of injectors and producers. Other key parameters considered in the modeling process were CO 2 purchase, gas recycle and addition of infill wells and/or patterns. The PRSM proxy model was ‘trained’ or calibrated with a series of training simulations. This involved an iterative process until the surrogate model reached a specific validation criterion. A sensitivity analysis was first conducted to ascertain which of these control variables to retain in the surrogate model. A genetic algorithm with a mixed-integer capability optimization approach was employed to determine the optimum developmental strategy to maximize both oil recovery and CO 2 storage. The proxy model reduced the computational cost significantly. The validation criteria of the reduced order model ensured accuracy in the dynamic modeling results. The prediction outcome suggested robustness and reliability of the genetic algorithm for optimizing both oil recovery and CO 2 storage. The reservoir modeling approach used in this study illustrates an improved approach to optimizing oil production and CO 2 storage within partially depleted oil reservoirs such as FWU. Lastly, this study may serve as a benchmark for potential CO 2 –EOR projects in the Anadarko basin and/or geologically similar basins throughout the world.

02 PETROLEUM↗

DOE OSTI2024

FLOWERS AEP: An Analytical Model for Wind Farm Layout Optimization

Annual energy production (AEP) is commonly used in objective functions for wind farm layout optimization. AEP is proportional to wind farm power production integrated over an annual distribution of free-stream wind conditions. Physics-based estimates of wind farm power production typically rely on low-fidelity engineering wake models that approximate the steady-state wind farm flow field. AEP estimates are then obtained by performing independent simulations for discrete wind conditions and using rectangular quadrature to account for each condition's expected frequency of occurrence. Depending on the number of simulated discrete wind conditions, this numerical integral could be hampered by poor accuracy or high computational costs. The FLOWERS AEP model instead poses an analytical integral of the engineering wake model over the variable wind conditions, yielding a closed-form, analytical function for wind farm AEP. This paper derives the analytical functions for FLOWERS AEP and its derivatives with respect to turbine position, which are useful for gradient-based wind farm layout optimization, in nondimensional form. We then analyze the benefits of the FLOWERS AEP model over conventional reference models, focusing on its low cost, adequate wake loss predictions, and smooth design space. Although the FLOWERS approach is found to predict the exact value of AEP with some error relative to the reference model (within 14% on average), it dramatically reduces computation time by an order of magnitude, produces a qualitatively similar design space at relatively low resolution, and yields comparable optimal layouts. This significant speed improvement is critical in layout optimization applications, where determining an optimal layout in an efficient manner is more important than precise AEP prediction.

17 WIND ENERGY↗

DOE OSTI2018

CFD Modeling of Boron Removal from Liquid Silicon with Cold Gases and Plasma

The present study focuses on a specific step of the metallurgical path of purification to provide solar-grade silicon: the removal of boron through the injection of H{sub 2}O(g)-H{sub 2}(g)-Ar(g) (cold gas process) or of Ar-H{sub 2}-O{sub 2} plasma (plasma process) on stirred liquid silicon. We propose a way to predict silicon and boron flows from the liquid silicon surface by using a CFD model (©Ansys Fluent) combined with some results on one-dimensional diffusive-reactive models to consider the formation of silica aerosols in a layer above the liquid silicon. The comparison of the model with experimental results on cold gas processes provided satisfying results for cases with low and high concentrations of oxidants. This confirms that the choices of thermodynamic data of HBO(g) and the activity coefficient of boron in liquid silicon are suitable and that the hypotheses regarding similar diffusion mechanisms at the surface for HBO(g) and SiO(g) are appropriate. The reasons for similar diffusion mechanisms need further enquiry. We also studied the effect of pressure and geometric variations in the cold gas process. For some cases with high injection flows, the model slightly overestimates the boron extraction rate, and the overestimation increases with increasing injection flow. A single plasma experiment from SIMaP (France) was modeled, and the model results fit the experimental data on purification if we suppose that aerosols form, but it is not enough to draw conclusions about the formation of aerosols for plasma experiments.

36 MATERIALS SCIENCE↗

DOE OSTI2019

A Comparative Study of the Efficacy of Local/Global and Parametric/Nonparametric Machine Learning Methods for Establishing Structure–Property Linkages in High-Contrast 3D Elastic Composites

Reduced-order structure–property (S-P) linkages play a pivotal role in the tailored design of materials for advanced engineering components. There is a critical need to distill these from the simulation datasets aggregated using sophisticated, computationally expensive, physics-based numerical tools (e.g., finite element methods). The recent emergence of materials data science approaches has opened new avenues for addressing this challenge. In this paper, we critically compare the relative merits of the application of four distinct machine learning approaches for their efficacy in extracting microstructure-property linkages from the finite element simulation data aggregated on high-contrast elastic composites with different microstructures. The machine learning approaches selected for the study have included different combinations of local/global and parametric/nonparametric approaches. Furthermore, the nonparametric approaches selected for this study are based on Gaussian Process (GP) models that allow for a formal treatment of uncertainty quantification in the predicted values. The predictive performances of these different approaches have been compared against each other using rigorous cross-validation error metrics. Furthermore, their sensitivity to both the dataset size and dimensionality has been investigated.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Sustainable hydrogen manufacturing via renewable-integrated intensified process for refueling stations

The widescale consumer adoption of hydrogen fuel cell electric vehicles (HFCEVs) is currently hindered by the high cost of small-scale hydrogen generation and the lack of extensive hydrogen refueling infrastructure. Natural gas-based hydrogen is cheaper when produced in large volumes but is also associated with high CO 2 emissions. To counter these challenges, we propose a hybrid approach where both natural gas and renewables are integrated in a synergistic manner using a dynamic process intensification technology that can be deployed on-site for meeting local demands of refueling stations. The technology is based on sorption enhanced steam methane reforming (SE-SMR) that utilizes a combination of reaction with in-situ CO 2 adsorption for enhancing process modularity, productivity and efficiency thereby outperforming conventional SMR at small scale. We develop a mixed integer linear programming (MILP)-based optimization framework for simultaneous design and scheduling of the SE-SMR process. The simultaneous optimization provides a synergistic combination whereby the renewables allow sustainable hydrogen manufacturing and the dynamic SE-SMR allows optimal use of the intermittency of the renewables. The U.S. nationwide analysis indicates that for futuristic renewable prices and a hydrogen production capacity of 2 ton/day, hydrogen can be produced at 50% less cost compared to the current cost of small-scale hydrogen generation. Finally, the city-wise analysis with varying hydrogen demand shows that even with just 5% HFCEV market penetration level, hydrogen production cost less than $3/kg can be obtained at small scales across the United States with even cheaper hydrogen for large cities.

08 HYDROGEN↗

DOE OSTI2021

Assessment of model parameters in MFiX particle-in-cell approach

The limitations in numerical treatment of solids-phase in conventional methods like Discrete Element Model and Two-Fluid Model have facilitated the development of alternative techniques such as Particle-In-Cell (PIC). However, a number of parameters are involved in PIC due to its empiricism. In this work, global sensitivity analysis of PIC model parameters is performed under three distinct operating regimes common in chemical engineering applications, viz. settling bed, bubbling fluidized bed and circulating fluidized bed. Simulations were performed using the PIC method in Multiphase Flow with Interphase eXchanges (MFiX) developed by National Energy Technology Laboratory (NETL). A non-intrusive uncertainty quantification (UQ) based approach is applied using Nodeworks to first construct an adequate surrogate model and then identify the most influential parameters in each case. This knowledge will aid in developing an effective design of experiments and determine optimal parameters through techniques such as deterministic or statistical calibration.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2019

Regimes of premixed turbulent spontaneous ignition and deflagration under gas-turbine reheat combustion conditions

Here we address the problem of turbulent spontaneous ignition vs. deflagration propagating against a mean residence time gradient, particularly relevant to the stabilization of quasi-steady diesel flames and gas-turbine reheat combustion, by performing a large series of statistically one-dimensional premixed turbulent flame DNS. An autoignitve mixture of hydrogen-vitiated air representative of gas-turbine reheat combustion is considered. First, we demonstrate that turbulence can trigger the transition from spontaneous ignition to a deflagration. We then construct a spontaneous ignition vs. deflagration regime diagram based on inflow velocity, turbulent velocity fluctuations and integral length scale and identify the mechanism that governs the transition from one propagation mode to the other. Finally, we propose a theoretical model that captures accurately this transition and assess the validity of a proxy based on chemical explosive mode analysis to determine the local propagation mode in complex configurations.

42 ENGINEERING↗

DOE OSTI2022

In-situ imaging of flexure-induced fracture in tape-laminate composites using high-resolution X-ray computed tomography

Here, this study presents development of an in-situ flexure test for imaging progressive inter- and intralaminar fracture in tape-laminate composites using X-ray computed tomography (CT). The intent of this test is to provide detailed experimental observations of ply-level damage that can be used to validate existing, and develop new, progressive damage analysis (PDA) tools. The test consists of a vertically mounted specimen which is flexed using two eccentric compressive loads using an in-situ uniaxial load stage. The flexure specimen contains a starter notch mid-span which promotes initiation of composite failure within the X-ray field of view. Specimens with two different laminate stacking sequences were tested, imaged, analyzed. For a quasi-isotropic laminate with large angle changes between adjacent plies, there was near simultaneous growth of transverse cracks and delaminations below the midplane of the laminate. For a laminate with small angles between adjacent plies, there was extensive formation of transverse crack networks which penetrated the laminate thickness without delamination growth. In addition to imaging fracture, the X-ray CT data from both specimen types were used to quantify the variability in thickness and analyze the local orientation of individual plies. Overall, the proposed test and the image-data analysis methodology provided an important insight into the fracture processes in tape laminates and highlights the inherent ply-level geometrical variabilities that should be accounted for in PDA simulations.

42 ENGINEERING↗

DOE OSTI2026

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗

DOE OSTI2021

Joint impedance and facies inversion of time-lapse seismic data for improving monitoring of CO 2 incidentally stored from CO 2 EOR

Time-lapse seismic monitoring is an effective and proven technology for mapping the distribution of CO 2 in a subsurface reservoir. When injected CO 2 displaces other reservoir fluids, porous-medium properties are changed and thus the seismic impedance changes, causing time-lapse seismic amplitude differences in the injection zones. The analysis and interpretation of images created from these amplitude differences can provide information about reservoir architecture and the CO 2 migration within the reservoir. Incorporating seismic inversion and rock physics into the interpretation of time-lapse seismic data can considerably improve the modeling and monitoring to detect and assess the location of CO 2 over time. The joint inversion method presented in this paper has an integral representation of the geology in the inversion algorithm using elastic facies, which provides information about the spatial distribution of the geologic heterogeneities controlling the movement of fluids in the reservoir. The method was successfully applied to time-lapse seismic data from a mature oil field undergoing CO 2 enhanced oil recovery. The estimated seismic acoustic impedances and facies reflect the characteristics of individual geologic facies and fluid conditions of the reservoir subject to CO 2 injection. Finally, the probabilities estimated by the joint impedance and facies inversion for the reservoir's litho-fluid facies can be used for forecasting CO 2 saturation and pressure changes within the target reservoir.

4D seismic↗

DOE OSTI2019

Self-induced temperature correction for inter-phase heat transfer in Euler-Lagrange point-particle simulation

In a two-way coupled Euler-Lagrange simulation, particles are approximated as point sources and their momentum and energy exchange with the surrounding flow are modeled through point-particle force and heat transfer models. As the particle size increases and approaches the Eulerian grid, the feedback force and heat transfer become large and strongly affect the local flow at the particle location. The fluid velocity and temperature computed in the EL simulation when interpolated to the particle location will be substantially different from the corresponding undisturbed values. In this work we will follow the approach pursued in Balachandar et al. (2018) [1] for self-induced velocity correction to develop an analogous correction procedure for self-induced temperature. We obtain analytical solutions for the self-induced thermal perturbation for a steady uniform cross flow past a Gaussian filtered heat source under both steady and unsteady conditions. The analytical results are extended to finite Peclet number using numerical simulations. The resulting quasi-steady and unsteady models are tested with a simple problem of exponential thermal evolution against the corresponding EL simulation results. We also provide a simple criterion for when the self-induced temperature correction is needed in an Euler-Lagrange simulation.

42 ENGINEERING↗

DOE OSTI2020

Effect of isotropic and anisotropic porous microstructure on electrochemical performance of Li ion battery cathodes: An experimental and computational study

Liquid electrolyte mass transport is a major limitation affecting high-power Li ion batteries. Fast discharging causes Li salt depletion in the current collector region of the cathode which produces overpotential in the electrolyte and consequently a drop of cell voltage to below the cut-off voltage, especially at higher electrode thickness and discharge rate. In this study, through experiment and simulation, we have investigated the effect of electrode thickness, mass loading, discharge rate and tortuosity on electrolyte mass transport and final derived areal capacity and specific energy for electrodes having isotropic (normal tape casting) and anisotropic (freeze tape casting) porous microstructure. Furthermore, the macroporous channels in freeze tape cast electrodes facilitate Li salt transport and reduce the Li salt mass transport limitations even at high electrode thickness and discharge rates, and high electrode tortuosity. In our results, computer simulations show that freeze tape cast electrodes may be fully discharged up to 750 μm thickness at 1 C rate compared to 300 μm for normal tape cast electrodes with the same mass loading. Freeze tape cast electrodes also show stable maximum areal capacity for C rates about double the maximum C rates of their normal tape cast electrode counterparts with the same mass loading.

25 ENERGY STORAGE↗

DOE OSTI2021

Benchmarking the Performance of Neuromorphic and Spiking Neural Network Simulators

Software simulators play a critical role in the development of new algorithms and system architectures in any field of engineering. Neuromorphic computing, which has shown potential in building brain-inspired energy-efficient hardware, suffers a slow-down in the development cycle due to a lack of flexible and easy-to-use simulators of either neuromorphic hardware itself or of spiking neural networks (SNNs), the type of neural network computation executed on most neuromorphic systems. While there are several openly available neuromorphic or SNN simulation packages developed by a variety of research groups, they have mostly targeted computational neuroscience simulations, and only a few have targeted small-scale machine learning tasks with SNNs. Evaluations or comparisons of these simulators have often targeted computational neuroscience-style workloads. In this work, we seek to evaluate the performance of several publicly available SNN simulators with respect to non-computational neuroscience workloads, in terms of speed, flexibility, and scalability. We evaluate the performance of the NEST, Brian2, Brian2GeNN, BindsNET and Nengo packages under a common front-end neuromorphic framework. Our evaluation tasks include a variety of different network architectures and workload types to mimic the computation common in different algorithms, including feed-forward network inference, genetic algorithms, and reservoir computing. We also study the scalability of each of these simulators when running on different computing hardware, from single core CPU workstations to multi-node supercomputers. Our results show that the BindsNET simulator has the best speed and scalability for most of the SNN workloads (sparse, dense, and layered SNN architectures) on a single core CPU. However, when comparing the simulators leveraging the GPU capabilities, Brian2GeNN outperforms the others for these workloads in terms of scalability. NEST performs the best for small sparse networks and is also the most flexible simulator in terms of reconfiguration capability NEST shows a speedup of at least 2x compared to the other packages when running evolutionary algorithms for SNNs. The multi-node and multi-thread capabilities of NEST show at least 2x speedup compared to the rest of the simulators (single core CPU or GPU based simulators) for large and sparse networks. We conclude our work by providing a set of recommendations on the suitability of employing these simulators for different tasks and scales of operations. We also present the characteristics for a future generic ideal SNN simulator for different neuromorphic computing workloads.

97 MATHEMATICS AND COMPUTING↗