Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation results”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Benchmarking CO₂ storage simulations: Results from the 11 th Society of Petroleum Engineers Comparative Solution Project

The 11 th Society of Petroleum Engineers Comparative Solution Project (shortened SPE11 herein) benchmarked simulation tools for geological carbon dioxide (CO 2 ) storage. A total of 45 groups from leading research institutions and industry across the globe signed up to participate, with 18 ultimately contributing valid results that were included in the comparative study reported here. This paper summarizes the SPE11 results. A comprehensive introduction and qualitative discussion of the submitted data are provided, together with an overview of online resources for accessing the full depth of data. A global metric for analyzing the relative distance between submissions is proposed and used to conduct a quantitative analysis of the submissions. This analysis attempts to statistically resolve the key aspects influencing the variability between submissions. The study shows that the major qualitative variation between the submitted results is related to thermal effects, dissolution-driven convective mixing, and resolution of facies discontinuities. Moreover, a strong dependence on grid resolution is observed across all three versions of the SPE11. However, our quantitative analysis suggests that the observed variations are predominantly influenced by factors not documented in the technical responses provided by the participants. We therefore identify that unreported variations due to human choices within the process of setting up, conducting, and reporting on the simulations underlying each SPE11 submission are at least as impactful as the computational choices reported.

Nordbotten, Jan M. [Univ. of Bergen (Norway); Norw↗
DOE OSTI2025

Simulation Results of a Thermal Power Dispatch System from a Generic Pressurized Water Reactor in Normal and Abnormal Operating Conditions

Amid economic pressures in the U.S. electricity market, nuclear utilities are exploring new revenue streams, including hydrogen production. A generic pressurized water reactor simulator was modified to incorporate a novel design for a TPD system coupled to a hydrogen production plant. Standard malfunctions were included in the simulation design, including steam line breaks at various system locations and flow interruptions in the hydrogen plant due to multiple faults, reflecting anticipated operational challenges. It is imperative that the TPD system operation has a minimal effect on the reactor power, primary coolant system, and turbine system operation and performance. Due to the specific design and application of this TPD system, with the proposed turbine control system changes, the overall impact on the existing plant systems is low. Normal TPD operating scenarios resulted in minor effects on the existing plant systems: reactor power changes by at most 0.2%, and gross generator output changes by 20.5 MWe from 100 MWt of TPD. The most severe malfunction analyzed in this work is a full TPD steam line break downstream of the extraction location, which results in an increase in reactor power of about 0.5%. The gross generator output decreases by 36 MWe, a total decrease of 60 MWe from the full power steady state (FPSS) condition. These results indicate that an industrial hydrogen production plant could be coupled thermally to a nuclear power plant with limited effects on the existing system operation and safety.

08 HYDROGEN↗
DOE OSTI2023

Simulation results for a low energy nuclear recoil yields measurement in liquid xenon using the MiX detector

Measuring the scintillation and ionization yields of liquid xenon in response to ultra-low energy nuclear recoil events is necessary to increase the sensitivity of liquid xenon experiments to light dark matter. Neutron capture on xenon can be used to produce nuclear recoil events with energies below 0.3 keV NR via the asymmetric emission of γ rays during nuclear de-excitation. The feasibility of an ultra-low energy nuclear recoil measurement using neutron capture was investigated for the Michigan Xenon (MiX) detector, a small dual-phase xenon time projection chamber that is optimized for a high scintillation gain. Simulations of the MiX detector, a partial neutron moderator, and a pulsed neutron generator indicate that a population of neutron capture events can be isolated from neutron scattering events. Additionally, the rate of neutron captures in the MiX detector was optimized by varying the thickness of the partial neutron moderator, neutron pulse width, and neutron pulse frequency.

detector calibration↗
DOE OSTI2025

Comparison of laser-produced plasma spatio-temporal electron density evolution measured using interferometry with simulation results

Due to the difficulties associated with experimental measurements of laser-produced plasma (LPP) properties during the earliest stages of plasma evolution, radiation hydrodynamic codes are often used. However, although these codes have been extensively validated in the higher intensity regimes, validation at low to moderate intensities has been limited. In this study, the spatio-temporal electron density evolution in an LPP generated at moderate laser intensities and at various laser wavelengths was validated against the FLASH code for times up to 20 ns. The LPP was generated by focusing the fundamental and various harmonics radiation (1064, 532, and 266 nm) from a 6 ns full width half maximum Nd:YAG laser, at a laser intensity of 10 GW cm -2 , onto a copper target. The spatio-temporal density evolution of the expanding plasma was analyzed using Nomarski interferometry. Experimental measurements were found to be consistent with FLASH simulations, and the dependence of electron densities on wavelength was found to be in agreement with analytical models, varying as n e ∝ λ -0.7 . However, slight differences were noted in the widths and shape functions of the experimental and simulated electron density profiles.

2D fast imaging↗
DOE OSTI2021

Cryogenic Vibrationally Resolved Photoelectron Spectroscopy of OH–(H2O): Confirmation of Multidimensional Franck-Condon Simulation Results for the Transition State of the OH + H2O Reaction

We present a transition state spectroscopic study of the OH + H2O reaction using the experimental technique of cryogenic Negative Ion Photoelectron Spectroscopy (NIPES). The recorded NIPE spectrum at 193 nm exhibits multiple vibrational progressions that include excitations to the shared H atom antisymmetric stretching mode with an interval of 0.32 eV as well as other progressions, mainly involving the H bending and O…O symmetric stretching modes. The Vertical Detachment Energy (VDE) was measured at 3.53 eV, whereas an upper limit for the Adiabatic Detachment Energy (ADE) was estimated at 2.90 eV. These values are in excellent agreement with the theoretically computed values of 3.51 eV and 2.83 eV respectively, obtained at the CCSD(T)/aug-cc-pV5Z level of theory. The recorded NIPE spectrum is in very good agreement when compared to the one recently reported from four-dimensional Franck-Condon simulations, in which a similar spectral profile was predicted. Besides observing the ground state, a charge-transfer excited state in the form of [OH–(H2O)+] is identified with a relative energy of 1.39 eV, well matching the previous prediction of 1.36 eV. This work was supported by U.S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences, and performed using EMSL, a national scientific user facility sponsored by DOE’s Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory, which is operated by Battelle Memorial Institute for the DOE. The theoretical calculations were conducted on EMSL’s “Cascade” Supercomputer. This research also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231.

Cao, Wenjin↗
DOE OSTI2020

BEPAM Model Code and CABBI Simulation Results for "Assessing the Additional Carbon Savings with Biofuel"

This dataset consists of various input data that are used in the GAMS model. All the data are in the format of .inc which can be read within GAMS or Notepad. Main data sources include: acreage data (acre), crop budget data ($/acre), crop yield data (e.g. bushel/acre), Soil carbon sequestration data (KgCO2/ha/yr). Model details can be found in the "Assessing the Additional Carbon Savings with Biofuel" and GAMS model package. ## File Description (1) GAMS Model.zip: This includes all the input files and scripts for running the model (2) Table*.csv: These files include the data from the tables in the manuscript (3) Figure2_3_4.csv: This contains the data used to create the figures in the manuscript (4) BaselineResults.csv: This includes a summary of the model results. (5) SensitivityResults_*.csv: Model results from the various sensitivity analyses performed (6) LUC_emission.csv: land use change emissions by crop reporting district for changes of pasturelands to annual crops.

Anticipated baseline approach↗
DOE OSTI2020

Compartmentalized Storage Reservoir Model and Simulation Results

A hypothetical closed (or compartmentalized) reservoir was created based on the numerical mesh used for Wainwright et al. (2013). The reservoir extent is smaller than the original model and limited between the two faults. The caprock properties are assumed in the rest of the storage formation.

Carbon Storage↗
DOE OSTI2026

Plasma processing of SRF cavities at Jefferson Lab: Experiment results and simulation insight

Plasma processing of superconducting radio frequency (SRF) cavities has been an active research effort at Jefferson Lab (JLab) since 2019, aimed at enhancing cavity performance by removing hydrocarbon contaminants and reducing field emission. In this experiment, processing using argon-oxygen and helium-oxygen gas mixtures to find minimum ignition power at different cavity pressure was investigated. Ongoing simulations are contributing to a better understanding of the plasma surface interactions and the fundamental physics behind the process. These simulations, combined with experimental studies, guide the optimization of key parameters such as gas type, RF power, and pressure to ignite plasma using selected higher-order mode (HOM) frequencies. This paper presents experimental data from argon-oxygen and helium-oxygen gas mixture C75 and C100 cavity plasma ignition studies, as well as simulation results for the C100-type cavity based on the COMSOL model previously applied to the C75 cavity.

Accelerator Physics↗
DOE OSTI2026

Effect of macroscopic surface defects on dynamic damage: An experimental and numerical study

This study examines the impact of macroscopic surface defects on the dynamic ductile damage behavior of polycrystalline metals using plate-impact experiments. Defects of various shapes (flat, round, and point) were manufactured on the free or impact surfaces of annealed copper specimens. The experiments were diagnosed with photon Doppler velocimetry measurements and soft recovery techniques. The experimental results revealed that defect shape and location significantly affect velocity–time profiles and void distribution. In order to understand the dynamics of shock propagation and corresponding ductile damage evolution within the target specimens, the impact experiments were modeled with a calibrated Tepla model, a dynamic ductile damage model for polycrystalline metals [Nguyen et al., Int. J. Solids Struct. 329, 113833 (2026)]. Overall, our resulting simulated velocity showed good agreement with measured velocity, and our simulated porosity distributions qualitatively matched experimental data. Based on our simulation results, defects on the free surface were found to distort rarefaction waves and therefore the corresponding spall planes. On the other hand, defects on the impact surface generated a delayed shock when the flyer plate and the defected target area were in contact, leading to distortion of the spall plane. The distortion of the spall plane resulted in a non-uniform distribution of voids within the defected specimens. Lastly, we examined the role of local work hardening due to the defect manufacturing process on the velocity and porosity distribution, using Tepla simulations with a simplified representation of local hardening. This investigation highlights the importance of defect geometry, location, and local hardening associated with defect generation in dynamic ductile damage processes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Dark Energy Survey Year 3 results: Simulation-based cosmological inference with wavelet harmonics, scattering transforms, and moments of weak lensing mass maps. Validation on simulations

Beyond-two-point statistics contain additional information on cosmological as well as astrophysical and observational (systematics) parameters. In this methodology paper we provide an end-to-end simulationbased analysis of a set of Gaussian and non-Gaussian weak lensing statistics using detailed mock catalogs of the Dark Energy Survey (DES). Here, we implement: 1) second and third moments; 2) wavelet phase harmonics (WPH); 3) the scattering transform (ST). Our analysis is fully based on simulations, it spans a space of seven $νw$CDM cosmological parameters, and it forward models the most relevant sources of systematics of the data (masks, noise variations, clustering of the sources, intrinsic alignments, and shear and redshift calibration). We implement a neural network compression of the summary statistics, and we estimate the parameter posteriors using a likelihood-free-inference approach. We validate the pipeline extensively, and we find that WPH exhibits the strongest performance when combined with second moments, followed by ST, and then by third moments. The combination of all the different statistics further enhances constraints with respect to second moments, up to 25 percent, 15 percent, and 90 percent for S 8 , Ω m , and the figure-of-merit FoM S8;Ωm , respectively. We further find that non-Gaussian statistics improve constraints on w and on the amplitude of intrinsic alignment with respect to second moments constraints. The methodological advances presented here are suitable for application to Stage IV surveys from Euclid, Rubin-LSST, and Roman with additional validation on mock catalogs for each survey. In a companion paper we present an application to DES Year 3 data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Dark Energy Survey Year 3 results: Simulation-based cosmological inference with wavelet harmonics, scattering transforms, and moments of weak lensing mass maps. II. cosmological results

Here, we present a simulation-based cosmological analysis using a combination of Gaussian and non-Gaussian statistics of the weak lensing mass (convergence) maps from the first three years of the Dark Energy Survey. We implement the following: (1) second and third moments; (2) wavelet phase harmonics; (3) the scattering transform. Our analysis is fully based on simulations, spans a space of seven 𝑤 Cold Dark Matter (𝑤⁢ CDM) cosmological parameters, and forward models the most relevant sources of systematics inherent in the data: masks, noise variations, clustering of the sources, intrinsic alignments, and shear and redshift calibration. We implement a neural network compression of the summary statistics, and we estimate the parameter posteriors using a simulation-based inference approach. Including and combining different non-Gaussian statistics is a powerful tool that strongly improves constraints over Gaussian statistics (in our case, the second moments); in particular, the figure of merit (𝑆 8 , Ω m ) is improved by 70% (Λ ⁢CDM) and 90% (𝑤 ⁢CDM). When all the summary statistics are combined, we achieve a 2% constraint on the amplitude of fluctuations parameter 𝑆 8 ≡ 𝜎 8 ⁢(Ω m /0.3) 0.5 , obtaining 𝑆 8 = 0.794 ±0.017 (Λ⁢ CDM) and 𝑆 8 = 0.817 ±0.021 (𝑤 ⁢CDM), and a ∼10% constraint on Ω m , obtaining Ω m =0.259 ±0.025 (Λ ⁢CDM) and Ω m = 0.273 ±0.029 (𝑤⁢ CDM). In the context of the 𝑤⁢ CDM scenario, these statistics also strengthen the constraints on the parameter 𝑤, obtaining 𝑤 <−0.72. The constraints from different statistics are shown to be internally consistent (with a 𝑝-value>0.1 for all combinations of statistics examined). We compare our results to other weak lensing results from the first three years of the Dark Energy Survey data, finding good consistency; we also compare with results from external datasets, such as planck constraints from the cosmic microwave background, finding statistical agreement, with discrepancies no greater than <2.2⁢𝜎.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Dimensional stability of low-cost thermoplastic composite molds

This report addresses a key question for using additively manufactured tooling from carbon fiber reinforced thermoplastic composites. How does the distortion of the molded part change with the print orientation used in the manufacture of tooling? The objective of this research was to discuss the distortion of parts molded with such tooling, in relation to the deflections in the tooling during the compression molding process. This was achieved through both model simulations and experimental molding runs. A vertical orientation was used in printing the mold with a distinctly non-planar surface: this orientation allowed for higher performance along the press closing direction -- higher thermal conductivity, higher stiffness, and lower coefficient of thermal expansion. Epoxy-carbon fiber twill weave fabric prepreg was then compression molded in this mold. The thermomechanical properties of the tooling material and the molded part were measured and used in simulation of mold deformation as well as part distortion. The thermomechanical anisotropy of the mold is quite different from that of the molded part because in the mold, the stiffest direction is the z-axis, while in the molded thermoset part, the z-axis is the weakest direction. The connection between mold deformation during the compression molding process and the final part distortion can be seen from the simulation results. The simulation results for the case where the mold was four times as stiff as the part along the press closing direction compared well with experiment. When the mold stiffness was lowered in relation to the part stiffness, the mold deformation during the compression molding increased; but this led to a smaller extent of part distortion.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Impact of undoped substrates on high performance silicon solar cells

Today's highest-efficiency silicon solar cells typically operate near the threshold between low-level and high-level injection. It is not well understood if pushing further into a regime in which the cell operating point is solidly in high-level injection at all times of the day has further benefits for the solar cell performance. From a reliability perspective, cells fabricated on lower doped silicon have a larger breakdown voltage. This advantage can affect the design of modules allowing higher voltages and a relaxation of the number of bypass diodes needed. In this project, we present a comprehensive assessment, both experimental and using simulation, of how bulk resistivity, light intensity, and operation temperature impact the performance of the silicon solar cell. This work incorporates a comprehensive device physics analysis assisted by numerical simulation. The source code is now available under General Public License (GPL3), and to further leverage the findings of this project and outreach people outside of the scientific community, we are working with www.pveducation.org (which receives >1 million visitors a year) to create interactive content using our simulation results. The simulation results indicate that high bulk resistivity wafers (>>10 Ωcm) require bulk Shockley-Read-Hall (SRH) lifetimes in the millisecond range to outperform wafers with standard bulk resistivities (<10 Ωcm). Additionally, above bulk resistivities of 10 Ωcm (the exact value depends on the bulk characteristics of the wafer), the cell efficiency is weakly dependent on the bulk resistivity. As a result, ingot manufactures may have an opportunity to further reduce wafer cost by growing higher resistivity ingots that are more tolerant to resistivity variations. This project is particularly relevant today, as solar cell architectures with improved surface passivation and milliseconds lifetimes wafers are commercially available, leveraging potential benefits of using higher bulk resistivities. Outside of interdigitated back contact (IBC) cell, reported studies on high resistivity silicon (>100 Ωcm) are very limited. To the best of our knowledge, this project provides for the first-time experimental insight on solar cells fabricated on wafers with bulk resistivities up to several thousand Ωcm, delivering a comprehensive vision of their performance under real-world temperature and light intensity operation conditions. We manufactured and characterized solar cells with bulk resistivities in the range of 1 Ωcm to >15k Ωcm. Under standard testing conditions (STC), we measured solar cells efficiencies over 20% over the entire range of bulk resistivities, using our baseline cell processing. To evaluate the cell performance in real-world operation conditions, the solar cells were measured at different temperatures (25-80°C) and at different light intensities (0.1-1 suns). The measurements show that the bulk resistivity does not impact the solar cell response to temperature and light intensity. Similar thermal coefficients (TC) were measured for standard and high bulk resistivities, and they are comparable with the TC values reported in the literature for standard bulk resistivities <10 Ωcm. After light soaking, the solar cell didn’t show signs of light-induced degradation (LID). This result was expected since n-type float zone (FZ) wafers were used in this work, i.e. low traces of boron and low concentration of oxygen (oxygen is typically found in the seed end of Czochralski (CZ) ingots). We measured for high bulk resistivities (>10 Ωcm) extremely high breakdown voltages (>1000V). In conclusion, the insight provided by this project can positively impact the levelized cost of energy (LCOE) of the photovoltaic systems through its effect on cell and ingot manufacturing yield, silicon cell power output, and module reliability.

14 SOLAR ENERGY↗