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

A newly proposed isotherm model to predict Cs exchange with crystalline silicotitanate in tank waste simulants

The Zheng Anthony Miller (ZAM) computer model, a multicomponent ion exchange model used to predict the exchange of Group I metals onto crystalline silicotitanate (CST), has historically been used to predict Cs distribution coefficients from Hanford and Savannah River Site (SRS) tank waste simulants. Comparison of experimentally determined Cs distribution coefficients from tank waste simulants with ZAM isotherm model predictions indicate overprediction of Cs and K distribution coefficients for simple and complex simulants with the engineered form of CST. Additionally, recent changes in chemical composition/manufacturing of IONSIV TM R9140-B have resulted in increased Cs capacity from high-salt, highly alkaline solutions. Here, this work served to assess different isotherm models and refine equilibrium parameters to develop a model that can be applied to Hanford and SRS tank waste Cs removal efforts. Toward this goal, the Campbell Westesen Peterson (CWP) model was developed. This model utilized the experimentally determined Cs capacity, and simplified ZAM equilibria expressions to include only the binary substitution of Cs + or K + on the Na + sites. Equilibrium constants for these equations were refined using experimentally determined distribution coefficients. Overall, the CWP model significantly improved our ability to predict both Cs and K loading capacity from complex matrices.

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

Internal Radiation Dose Evaluation For An Unruptured Post Release Tristructural Istropic Fuel Particle for Advanced and Micro-reactor Applications

There are unique benefits from advanced/micro reactor designs and fuel types that offer safety features in the case of an accident that may reduce environmental consequences compared to conventional reactors and fuels. Several micro reactor designs, for both civilian and military applications, are proposed using tristructural-isotropic (TRISO) fuel. TRISO particles are a robust advanced nuclear fuel type that leads to the unique question of how unruptured, activated TRISO particles will interact with humans. Human exposure to nuclear fuel, or radionuclides, includes all pathways of exposure. However, TRISO particles are 0.9 mm in size, and that particle size restricts internal dose assessment to the ingestion pathway; eliminating the inhalation and dermal absorption pathways. Furthermore, activity of the TRISO particle was established by High Temperature Engineering Test Reactor simulations assuming the TRISO particle was one day post-release from the reactor. International Commission of Radiological Protection publications 100 and 133 were used to develop the internal dosimetry model for an unruptured TRISO particle for the International Commission of Radiological Protection reference male and female dose assessment. The rectosigmoid colon had the highest internal organ dose for both reference male (2.1 Sv) and reference female (2.3 Sv). The internal dose from ingestion of the scenario-defined TRISO particle was 0.25 Sv for the reference male and 0.29 Sv for the refence female, which exceeds the annual occupational effective dose limit of 0.05 Sv in the code of federal regulations, 10.CFR.20 Subpart C.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

Mixed-precision iterative refinement using tensor cores on GPUs to accelerate solution of linear systems

Double-precision floating-point arithmetic (FP64) has been the de facto standard for engineering and scientific simulations for several decades. Problem complexity and the sheer volume of data coming from various instruments and sensors motivate researchers to mix and match various approaches to optimize compute resources, including different levels of floating-point precision. In recent years, machine learning has motivated hardware support for half-precision floating-point arithmetic. A primary challenge in high-performance computing is to leverage reduced-precision and mixed-precision hardware. We show how the FP16/FP32 Tensor Cores on NVIDIA GPUs can be exploited to accelerate the solution of linear systems of equations Ax = b without sacrificing numerical stability. The techniques we employ include multiprecision LU factorization, the preconditioned generalized minimal residual algorithm (GMRES), and scaling and auto-adaptive rounding to avoid overflow. We also show how to efficiently handle systems with multiple right-hand sides. On the NVIDIA Quadro GV100 (Volta) GPU, we achieve a 4×-5× performance increase and 5× better energy efficiency versus the standard FP64 implementation while maintaining an FP64 level of numerical stability.

GMRES↗
DOE OSTI2024

GPU Accelerated Sparse Cholesky Factorization

The solution of sparse symmetric positive definite linear systems is an important computational kernel in large-scale scientific and engineering modeling and simulation. We will solve the linear systems using a direct method, in which a Cholesky factorization of the coefficient matrix is performed using a right-looking approach and the resulting triangular factors are used to compute the solution. Sparse Cholesky factorization is compute intensive. In this work we investigate techniques for reducing the factorization time in sparse Cholesky factorization by offloading some of the dense matrix operations on a GPU. We will describe the techniques we have considered. We achieved up to 4x speedup compared to the CPU-only version.

Karsavuran, M Ozan↗
DOE OSTI2018

Geocellular model of St. Peter Sandstone for University of Illinois at Urbana-Champaign DDU Feasibility Study

The geocellular model of the St. Peter Sandstone was constructed for the University of Illinois at Urbana-Champaign DDU feasibility study. Starting with the initial area of review (18.0 km by 18.1 km [11.2 miles by 11.3 miles]) the boundaries of the model were trimmed down to 9.7 km by 9.7 km (6 miles by 6 miles) to ensure that the model enclosed a large enough volume so that the cones of depression of both the production and injection wells would not interact with each other, while at the same time minimizing the number of cells to model to reduce computational time. The grid-cell size was set to 61.0 m by 61.0 m (200 feet by 200 feet) for 160 nodes in the X and Y directions. The top surface of the St. Peter Sandstone was provided by geologists working on the project, and the average thickness of the formation was taken from the geologic prospectus they provided. An average thickness of 68.6 m (225 feet) was used for the St. Peter Sandstone, resulting in 45 layers for the model. Petrophysical data was taken from available rotary sidewall core data (Morrow et al., 2017). As geothermal properties (thermal conductivity, specific heat capacity) are closely related to mineralogy, specifically the percentage of quartz, available mineralogical data was assembled and used with published data of geothermal values to determine these properties (Waples and Waples, 2004; Robertson, 1988). The St. Peter Sandstone was divided into facies according to similar geothermal and petrophysical properties, and distributed according to available geophysical log data and prevailing interpretations of the depositional/diagenetic history (Will et al. 2014). Petrophysical and geothermal properties were distributed through geostatistical means according to the associated distributions for each lithofacies. The formation temperature was calculated, based on data from continuous temperature geophysical log from a deep well drilled into the Precambrian basement at the nearby Illinois Basin Decatur Project (IBDP) where CO2 is currently being sequestered (Schlumberger, 2012). Salinity values used in the model were taken from regional studies of brine chemistry in the St. Peter Sandstone, including for the IBDP (e.g., Panno et al. 2018). After being reviewed by the project's geologists, the model was then passed onto the geological engineers to begin simulations of the geothermal reservoir and wellbores.

15 GEOTHERMAL ENERGY↗
DOE OSTI2018

Geocellular Model of Mt. Simon Sandstone for University of Illinois at Urbana-Champaign DDU feasibility study

The geocellular model of the Mt. Simon Sandstone was constructed for the University of Illinois at Urbana-Champaign DDU feasibility study. Starting with the initial area of review (18.0 km by 18.1 km [11.2 miles by 11.3 miles]) the boundaries of the model were trimmed down to 9.7 km by 9.7 km (6 miles by 6 miles) to ensure that the model enclosed a large enough volume so that the cones of depression of both the production and injection wells would not interact with each other, while at the same time minimizing the number of cells to model to reduce computational time. The grid-cell size was set to 61.0 m by 61.0 m (200 feet by 200 feet) for 160 nodes in the X and Y directions. Within the model, 67 layers are represented that are parameterized with their sediment/rock properties and petrophysical data. The top surface of the Mt. Simon Sandstone was provided by geologists working on the project, and the average thickness of the formation was taken from the geologic prospectus they provided. An average thickness of 762 m (2500 feet) was used for the Mt. Simon Sandstone, resulting in 60 layers for the model. Petrophysical data was taken from available rotary sidewall core data (Morrow et al., 2017). As geothermal properties (thermal conductivity, specific heat capacity) are closely related to mineralogy, specifically the percentage of quartz, available mineralogical data was assembled and used with published data of geothermal values to determine these properties (Waples and Waples, 2004; Robertson, 1988). The Mt. Simon Sandstone was divided into three separate units (lower, middle, upper) according to similar geothermal and petrophysical properties, and distributed according to available geophysical log data and prevailing interpretations of the depositional/diagenetic history (Freiburg et al. 2016). Petrophysical and geothermal properties were distributed through geostatistical means according to the associated distributions for each lithofacies. The formation temperature was calculated, based on data from continuous temperature geophysical log from a deep well drilled into the Precambrian basement at the nearby Illinois Basin Decatur Project (IBDP) where CO2 is currently being sequestered (Schlumberger, 2012). Salinity values used in the model were taken from regional studies of brine chemistry in the Mt. Simon Sandstone, including for the IBDP (e.g., Panno et al. 2018). After being reviewed by the project's geologists, the model was then passed onto the geological engineers to begin simulations of the geothermal reservoir and wellbores.

15 GEOTHERMAL ENERGY↗
DOE OSTI2018

IC Annual Report for w17_solarcells

We used first principles simulations to guide engineering of nanoparticle solar cells by considering effects of temperature and nanoparticle interactions

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2018

Implementation of a Grid Connected Battery-Inverter Fleet Model

Batteries are designed to store electrical energy. The increasing variation in time value of energy has driven the use of batteries as controllable distributed energy resources (DER). This is enabled though low-cost power electronic inverters that are able to precisely control charge and discharge. This paper describes the software implementation of an open-source battery inverter fleet models in python. The Sandia BatterylnverterFleet class model can be used by scientists, researchers, and engineers to perform simulations of one or more fleets of similar battery-inverter systems connected to the grid. The program tracks the state- of-charge of the simulated batteries and ensures that they stay within their limits while responding to separately generated service requests to charge or discharge. This can be used to analyze control and coordination, placement and sizing, and many other problems associated with the integration of batteries on the power grid. The development of these models along with their python implementation was funded by the Grid Modernization Laboratory Consortium (GMLC) project 1.4.2. Definitions, Standards and Test Procedures for Grid Services from Devices.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Constrained Bayesian Optimization of Criticality Experiments at LLNL [Slides]

The design of criticality experiments is typically an iterative process that employs a Monte Carlo transport code. The goal is to find a design that optimizes some variable, like the sensitivity of a response to a cross section, while simultaneously ensuring criticality. The high fidelity of the Monte Carlo code is a great asset, but it makes exploring the design space computationally expensive. Herein, we present how a constrained Bayesian optimization algorithm can be used to efficiently design a criticality experiment. It uses Gaussian processes as a surrogate model to probe the design space and to reduce the number of code executions that are needed to find the optimum. We demonstrate constrained Bayesian optimization with a Pu-239/polyethylene solution system and a TEX experiment that is designed for critical ity safety validation of a nuclear waste model at the Hanford Site. For both systems, a global optimum was found within 75 Monte Carlo simulations.

42 ENGINEERING↗
DOE OSTI2025

Distribution Feeder Characteristics and Their Resiliency to Natural Hazards

This paper introduces a method and approach for initially screening actions to improve the reliability of distribution systems exposed to three natural hazards (wildlife, weather, and vegetation). The method is well-suited to needs and capabilities of smaller utilities because it relies on readily available information and does not depend on simulations involving detailed engineering models. The method involves first aggregating feeders into like types and then correlating them, separately by hazard, to standard measures of reliability, namely, System Average Interruption Frequency Index (SAIFI) and System Average Interruption Duration Index (SAIDI), and the individual constituents of these measures. The correlations provide both new insights into the efficacy of accepted reliability management practices while and re-confirm long-understood, accepted practices. Importantly, the method does not seek to replace the need for detailed engineering analysis. Instead, in view of the significant costs involved in conducting these more involved analyses, the method is intended to help smaller utilities prioritize their more limited resources to maximize the efficacy of actions they take to improve reliability.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Real-Time Characterization of Salt Aerosols Generated from Static and Sparged Molten Salt

The formation of radionuclide-bearing aerosols in the respirable size range has the potential to significantly influence offsite dose consequences and is, therefore, an important consideration in nuclear facility safety assessments. Molten salt reactor (MSR) developers will likely need to demonstrate an understanding of the conditions under which radionuclide-bearing aerosols may be generated from their reactor under normal operating and accident conditions, as well as the characteristics and transport behavior of these aerosols, to demonstrate to the U.S. Nuclear Regulatory Commission (NRC) that the facility can be operated safely. Recent reviews of the literature identified a lack of experimental data describing the mechanisms of formation and properties (size, concentration, and composition) of salt aerosol particles that are produced from molten salts. Experiments that identify the conditions that lead to radionuclide-bearing salt aerosol releases and quantify the characteristics of salt aerosols formed by different mechanisms are high-priority needs to support MSR licensing. This report describes tests that were conducted within the Argonne Salt Aerosol Test Stand (a sealed vessel and measurement system) to generate salt aerosols from static and sparged molten salts and measure their size and concentration in real-time. The results provide insight into salt aerosol formation by the vapor condensation and bubble bursting mechanisms and inform the potential radiological consequences of aerosol formation from molten fuel salt. Videos of the salt surface were taken during salt sparge tests to observe surface bubble behavior. The data in this report can be used to develop mechanistic source term and accident progression models for MSRs. The real-time salt aerosol characterization technique used in this study will be employed in future integral effects tests that are conducted at an engineering scale to simulate realistic MSR accidents and in future separate effects tests to address additional variables that may impact salt aerosol characteristics (e.g., presence of fission products in salt and humidity in atmosphere).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Connecting Risk and Resilience for a Power System Using the Portland Hills Fault Case Study

Active seismic faults in the Pacific Northwest area have encouraged electric utilities in the region to deeply contemplate and proactively intervene to support grid resilience. To further this effort this research introduces Monte Carlo (MC)-based power system modeling as a means to inform the Performance Based Earthquake Engineering method and simulates 100,000 sample earthquakes of a 6.8 magnitude (M6.8) Portland Hills Fault (PHF) scenario in the Portland General Electric (PGE) service territory as a proof of concept. This paper also proposes the resilience metric Seismic Load Recovery Factor (SLRF) to quantify the recovery of a downed power system and thus can be used to quantify earthquake economic risk. Using MC results, the SLRF was evaluated to be 19.7 h and the expected economic consequence cost of a M6.8 PHF event was found to be $\$180$ million with an annualized risk of $\$90,000$ given the event’s 1 in 2000 year probability of occurrence. The MC results also identified the eight most consequential substations in the PGE system—i.e., those that contributed to maximum load loss. This paper concludes that retrofitting these substations reduced the expected consequence cost of a M6.8 PHF event to $117 million.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

MACHINE LEARNING-ENABLED PREDICTION OF TRANSIENT INJECTION MAP IN AUTOMOTIVE INJECTORS WITH UNCERTAINTY QUANTIFICATION

Accurate prediction of injection profiles is a critical aspect of linking injector operation with engine performance and emissions. However, highly resolved injector simulations can take one to two weeks of wall-clock time, which is incompatible with engine design cycles with desired turnaround times of less than a day. Hence, it is important to reduce the time-to-solution of the internal flow simulations by several orders of magnitude to make it compatible with engine simulations. This work demonstrates a data-driven approach for tackling the computational overhead of injector simulations, whereby the transient injection profiles are emulated for a side-oriented, single-hole diesel injector using a Bayesian machine-learning framework. First, an interpretable Bayesian learning strategy was employed to understand the effect of design parameters on the total void fraction field. Then, autoencoders are utilized for efficient dimensionality reduction of the flowfields. Gaussian process models are finally used to predict the spatiotemporal void fraction field at the injector exit for unknown operating conditions. The Gaussian process models produce principled uncertainty estimates associated with the emulated flowfields, which provide the engine designer with valuable information of where the data-driven predictions can be trusted in the design space. The Bayesian flowfield predictions are compared with the corresponding predictions from a deep neural network, which has been transfer-learned from static needle simulations from a previous work by the authors. The emulation framework can predict the void fraction field at the exit of the orifice within a few seconds, thus achieving a speed-up factor of up to 38 x 10(6) over the traditional simulation-based approach of generating transient injection maps.

machine learning↗
DOE OSTI2024

teemi: An open-source literate programming approach for iterative design-build-test-learn cycles in bioengineering

Synthetic biology dictates the data-driven engineering of biocatalysis, cellular functions, and organism behavior. Integral to synthetic biology is the aspiration to efficiently find, access, interoperate, and reuse high-quality data on genotype-phenotype relationships of native and engineered biosystems under FAIR principles, and from this facilitate forward-engineering strategies. However, biology is complex at the regulatory level, and noisy at the operational level, thus necessitating systematic and diligent data handling at all levels of the design, build, and test phases in order to maximize learning in the iterative design-build-test-learn engineering cycle. To enable user-friendly simulation, organization, and guidance for the engineering of biosystems, we have developed an open-source python-based computer-aided design and analysis platform operating under a literate programming user-interface hosted on Github. The platform is called teemi and is fully compliant with FAIR principles. In this study we apply teemi for i) designing and simulating bioengineering, ii) integrating and analyzing multivariate datasets, and iii) machine-learning for predictive engineering of metabolic pathway designs for production of a key precursor to medicinal alkaloids in yeast. The teemi platform is publicly available at PyPi and GitHub.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Evaluation of Combustion Models for CFD Simulation of Pre-Chamber Ignition in a Natural Gas Engine

Natural gas is considered an attractive alternative fossil fuel as it produces 25% less carbon dioxide and emits fewer amounts of toxic pollutants compared to gasoline and diesel. A pre-chamber type ignition device in a natural gas engine is capable of extending the lean operating limit, and hence improving thermal efficiency, by providing spatially distributed ignition with multiple turbulent flame-jets and accelerating the combustion rate. However, the process of ignition in the main chamber triggered by turbulent jets emanating from a pre-chamber is very sensitive to the pre-chamber geometry and fueling strategy, and any effort to optimize these parameters requires fundamental understanding of the underlying physical and thermodynamic phenomena, which can be studied by a multi-dimensional numerical simulation approach. In this study, three turbulent combustion models widely adopted in the multi-dimensional engine modeling community were preliminary evaluated to simulate the combustion process in a natural gas engine with a pre-chamber spark-ignition system. One multi-zone well-stirred reactor model, and two flamelet-based combustion models, G-equation and Extended Coherent Flamelet Model, were employed for this assessment. Under a fuel-lean condition, the three combustion models were used to match experimental data and analyze the differences in terms of in-cylinder pressure and heat release rate.

combustion model evaluation↗
DOE OSTI2020

Mixing and detonation structure in a rotating detonation engine with an axial air inlet

High-fidelity simulations of an experimental rotating detonation engine with an axial air inlet were conducted. The system operated with hydrogen as fuel at globally stoichiometric conditions. Instantaneous data showed that the detonation front is highly corrugated, and is considerably weaker than an ideal Chapman–Jouguet wave. Regions of deflagration are present ahead of the wave, caused by mixing with product gases from the previous cycle, as well as the injector recovery process. It is found that as the post-detonation high pressure flow expands, the injectors recover unsteadily, leading to a transient mixing process ahead of the next cycle. The resulting flow structure not only promotes mixing between product and reactant gases, but also increases likelihood of autoignition. These results show that the detonation process is very sensitive to injector design and the transient behavior during the detonation cycle. Phase-averaged statistics and conditionally averaged data are used to understand the overall reaction structure. As a result, comparisons with available experimental data on this configuration show remarkable good agreement of the predicted reacting flow structure.

42 ENGINEERING↗
DOE OSTI2024

On the importance of species immiscibility in mixing-layer dynamics at supercritical pressures

The mixing dynamics of injected propellants is a key factor in determining the ignition performance and combustion-instability response of rocket engines, internal combustion engines, and gas turbines. A key source of uncertainty in the prediction of phase-exchange dynamics is the immiscibility of the injected propellants at supercritical pressure conditions, which induces liquid/vapor phase separation and surface-tension dynamics. While experimental observations indicate the presence of liquid/vapor interfacial structures, this interfacial dynamics is typically neglected in numerical analyses. To address this issue, the objective of the present study is to systematically evaluate the importance of species immiscibility on the phase-exchange dynamics of cryogenic LOX/GH2 mixing layers at typical rocket engine injection conditions. This is accomplished by comparing simulations of (i) a recently developed interface-capturing Regularized-Interface Method (RIM) formulation, and (ii) the commonly employed Diffuse-Interface Method (DIM) formulation. Analysis shows that the interfacial dynamics significantly impact the atomization and mixing of the propellants in the near-injector region, which is not captured by the DIM formulation. The findings of this study extend to other immiscible injection systems, such as LOX/kerosene and hydrocarbon-fuel/air, thereby highlighting the importance of resolving species immiscibility in simulating high-pressure combustion engines.

diffused-interface method↗