Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “production simulation codes”

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 523 records · Page 29

Monte Carlo Methodology Serves Up a Software Success

Widely used for the modeling of gas flows through the computation of the motion and collisions of representative molecules, the Direct Simulation Monte Carlo method has become the gold standard for producing research and engineering predictions in the field of rarefied gas dynamics. Direct Simulation Monte Carlo was first introduced in the early 1960s by Dr. Graeme Bird, a professor at the University of Sydney, Australia. It has since proved to be a valuable tool to the aerospace and defense industries in providing design and operational support data, as well as flight data analysis. In 2002, NASA brought to the forefront a software product that maintains the same basic physics formulation of Dr. Bird's method, but provides effective modeling of complex, three-dimensional, real vehicle simulations and parallel processing capabilities to handle additional computational requirements, especially in areas where computational fluid dynamics (CFD) is not applicable. NASA's Direct Simulation Monte Carlo Analysis Code (DAC) software package is now considered the Agency s premier high-fidelity simulation tool for predicting vehicle aerodynamics and aerothermodynamic environments in rarified, or low-density, gas flows.

Source record↗

Optimizing Simulation Fidelity in Direct-Drive Inertial Confinement Fusion with Cassio

Recently, the National Ignition Facility (NIF) demonstrated that inertial confinement fusion (ICF) is capable to achieve thermonuclear (TN) ignition in the laboratory, making it a crucial method on the path to replicate the Sun’s power production mechanism on Earth. However, the physics governing the high-energy density environments is very complex and remains a challenge to fully understand and model. For example, dopants in the TN fuel are important diagnostic tools to extract the thermodynamic conditions of the plasma. However, if their concentration is chosen too high, they can significantly degrade the performance of an ICF capsule. In this study, we use the Los Alamos National Laboratory radiation-hydrodynamics code Cassio to model ICF implosions of capsules which contain deuterium fuel with high-Z dopants like Krypton and Argon from the high-Z campaign conducted 15 years ago. We focus on how chosen computational and physics parameters influence the implosion outcomes. By systematically changing the resolution of the computational mesh and the photon energies as well as modifying settings for the laser drive and TN fuel pre-heat effects, we assess the impact on experimentally measured performance metrics like neutron production from TN burn and x-ray emission during the implosion. Our results will help to improve the fidelity of simulations and guide future numerical studies and experimental designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

AN OPEN-SOURCE NEUTRINO RADIATION HYDRODYNAMICS CODE FOR CORE-COLLAPSE SUPERNOVAE

We present an open-source update to the spherically symmetric, general-relativistic hydrodynamics, core-collapse supernova (CCSN) code GR1D. The source code is available at http://www.GR1Dcode.org. We extend its capabilities to include a general-relativistic treatment of neutrino transport based on the moment formalisms of Shibata et al. and Cardall et al. We pay special attention to implementing and testing numerical methods and approximations that lessen the computational demand of the transport scheme by removing the need to invert large matrices. This is especially important for the implementation and development of moment-like transport methods in two and three dimensions. A critical component of neutrino transport calculations is the neutrino–matter interaction coefficients that describe the production, absorption, scattering, and annihilation of neutrinos. In this article we also describe our open-source neutrino interaction library NuLib (available at http://www.nulib.org). We believe that an open-source approach to describing these interactions is one of the major steps needed to progress toward robust models of CCSNe and robust predictions of the neutrino signal. We show, via comparisons to full Boltzmann neutrino-transport simulations of CCSNe, that our neutrino transport code performs remarkably well. Furthermore, we show that the methods and approximations we employ to increase efficiency do not decrease the fidelity of our results. We also test the ability of our general-relativistic transport code to model failed CCSNe by evolving a 40-solar-mass progenitor to the onset of collapse to a black hole.

79 ASTRONOMY AND ASTROPHYSICS↗

MEXIICO: a new equipment to study the pressure impact on the irradiated fuel behavior

The new MEXIICO experimental loop has recently been started in the LECA-STAR hot-cell facility of the CEA Cadarache. It is dedicated to the understanding of the basic mechanisms occurring in nuclear fuels during nuclear reactor power transients and accidental situations. Under such conditions the nuclear fuels are submitted to thermo-mechanical stresses. These stresses play a major role on fission product release and fuel fragmentation. Evaluating such fuel behavior represents a key issue for reactor safety analyses. In the first part of the paper, we describe the complete MEXIICO device, which main component is a hot isostatic press allowing to perform analytical experiments in the [20 -1600 C. degrees / 10 MPa - 160 MPa] region. As shown by simulations (ALCYONE fuel performance code), the upper value of the pressure (160 MPa) is consistent with the one experienced by fuels under nuclear power transients. We also present the experimental loop starting from the argon pressurized feed-line, the furnace and its associated instrumentation and the exhaust line equipped with specific valves and flowmeters controlling the pressure and flowrate. We then focus on the gamma spectrometry measurement which allows recording the fission gas release kinetics. In the second part of the paper, we detail the capability of the MEXIICO device to provide new experimental data regarding the influence of pressure on the fission gas release and fuel fragmentation, and to determine a pressure threshold above which fission gas release and fuel fragmentation may be inhibited. This device also allows investigating the impact of pressure on the microstructure evolution, coupling it with post-test microanalysis examinations: this procedure will allow characterizing the effect of temperature and pressure on the bubble evolution and interconnection. Evolutions of the fuel modeling, needed to simulate correctly these new experimental data, will be also presented. In the last section of the paper, we present the commissioning tests of the MEXIICO device showing the good reproducibility of the tests performed first in conventional laboratory and then in hot-cell laboratory. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Monte Carlo electron-photon transport using GPUs as an accelerator: Results for a water-aluminum-water phantom

An electron-photon coupled Monte Carlo code ARCHER - Accelerated Radiation-transport Computations in Heterogeneous Environments - is being developed at Rensselaer Polytechnic Institute as a software test bed for emerging heterogeneous high performance computers that utilize accelerators such as GPUs. In this paper, the preliminary results of code development and testing are presented. The electron transport in media was modeled using the class-II condensed history method. The electron energy considered ranges from a few hundred keV to 30 MeV. Moller scattering and bremsstrahlung processes above a preset energy were explicitly modeled. Energy loss below that threshold was accounted for using the Continuously Slowing Down Approximation (CSDA). Photon transport was dealt with using the delta tracking method. Photoelectric effect, Compton scattering and pair production were modeled. Voxelised geometry was supported. A serial ARHCHER-CPU was first written in C++. The code was then ported to the GPU platform using CUDA C. The hardware involved a desktop PC with an Intel Xeon X5660 CPU and six NVIDIA Tesla M2090 GPUs. ARHCHER was tested for a case of 20 MeV electron beam incident perpendicularly on a water-aluminum-water phantom. The depth and lateral dose profiles were found to agree with results obtained from well tested MC codes. Using six GPU cards, 6x10{sup 6} histories of electrons were simulated within 2 seconds. In comparison, the same case running the EGSnrc and MCNPX codes required 1645 seconds and 9213 seconds, respectively, on a CPU with a single core used. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Manuscript model outputs, model source code, and figure scripts: The role of geomorphology in mediating biomass allocation impacts on salt-marsh resilience and carbon accumulation

This data package provides model source code (C++), model outputs, and figure generation scripts (R) needed to reproduce the following manuscript: Bruns, Nicholas E., Genevieve L. Noyce, and Matthew L. Kirwan. "The Role of Geomorphology in Mediating Biomass Allocation Impacts on Salt-Marsh Resilience and Carbon Accumulation." Estuarine, Coastal and Shelf Science 327 (December 2025): 109549. https://doi.org/10.1016/j.ecss.2025.109549. This manuscript investigates how geomorphology mediates the impact of biomass allocation shifts on salt marsh persistence and carbon (C) sequestration under sea level rise. We use a 1-D soil-column model (Kirwan and Mudd 2012) to perform experiments across a range of static root:shoot ratios (RSR = 1-4) spanning observed values. The model explicitly simulates interactions between tidal inundation, productivity, inorganic sediment deposition, and organic matter accumulation. Experiments determine whether geomorphic feedbacks amplify, dampen, or leave unchanged the ecosystem response to biomass allocation shifts. A first experiment uses constant sea level rise (2.5 mm/yr) to examine equilibrium dynamics and their influence on carbon accumulation rates across different suspended sediment concentrations (SSC). A second set of experiments calculates threshold sea level rise rates for marsh drowning across RSR and SSC combinations. Final experiments apply accelerating sea level rise scenarios (2000-2200) derived from NOAA projections (Sweet et al. 2022) to generate an envelope of expected responses, quantifying the importance of biomass allocation shifts on marsh carbon accumulation and persistence. All experiments are repeated across SSC ranging from 5-50 mg/L to investigate how these interactions vary in micro-tidal marshes with different sediment supplies. Package contents: * README.txt with detailed description of package contents and instructions for reproducing manuscript figures and model outputs * R scripts for generating all manuscript figures * C++ baseline model code from Kirwan and Mudd (2012) * C++ source code for 4 experimental model variants used in the manuscript, extending above baseline code * Model input files (.csv, .txt) including sea level rise scenarios * Model output files (.txt) used in manuscript analyses Temporal coverage: Model simulations span years 1900-2200, with accelerating sea level rise scenarios for 2000-2200. Key variables: Root:shoot ratio, suspended sediment concentration, carbon accumulation rate, vertical accretion rate, marsh elevation, inundation depth, threshold sea level rise rate.

54 ENVIRONMENTAL SCIENCES↗

A Fluka study of underground cosmogenic neutron production

Neutrons produced by cosmic muon interactions are important contributors to backgrounds in underground detectors when searching for rare events. Typically such neutrons can dominate the background, as they are particularly difficult to shield and detect. Since actual data is sparse and not well documented, simulation studies must be used to design shields and predict background rates. Thus validation of any simulation code is necessary to assure reliable results. This work compares in detail predictions of the FLUKA simulation code to existing data, and uses this code to report a simulation of cosmogenic backgrounds for typical detectors embedded in a water tank with liquid scintillator shielding.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Module-based Hybrid Uncertainty Quantification for Multi-physics Applications: Theory and Software

In this project we proposed to develop an innovative uncertainty quantification methodology that captures the best of the two competing approaches in UQ, namely, intrusive and non-intrusive approaches. The idea is to develop the mathematics and the associated computational framework and algorithms to facilitate the use of intrusive or non-intrusive UQ methods in different modules of a multi-physics multi-module simulation model in a way that physics code developers for different modules are shielded (as much as possible) from the chores of accounting for the uncertain ties introduced by the other modules. As the result of our research and development, we have produced a number of publications, conference presentations, and a software product.

97 MATHEMATICS AND COMPUTING↗

The Role of Dredge-up in Double White Dwarf Mergers

We present the results of an investigation of the dredge-up and mixing during the merger of two white dwarfs (WDs) with different chemical compositions by conducting hydrodynamic simulations of binary mergers for three representative mass ratios. In all the simulations, the total mass of the two WDs is $\lesssim$1.0 M ⊙ . Mergers involving a CO and a He WD have been suggested as a possible formation channel for R Coronae Borealis (RCB)–type stars, and we are interested in testing if such mergers lead to conditions and outcomes in agreement with observations. Even if the conditions during the merger and subsequent nucleosynthesis favor the production of 18 O, the merger must avoid dredging up large amounts of 16 O, or else it will be difficult to produce sufficient 18 O to explain the oxygen ratio observed to be of order unity. Here, we performed a total of nine simulations using two different grid-based hydrodynamics codes using fixed and adaptive meshes and one smooth particle hydrodynamics (SPH) code. We find that in most of the simulations, >10 -2 M ⊙ of 16 O is indeed dredged up during the merger. However, in SPH simulations where the accretor is a hybrid He/CO WD with a ~0.1 M ⊙ layer of helium on top, we find that no 16 O is being dredged up, while in the q = 0.8 simulation <10 -4 M ⊙ of 16 O has been brought up, making a WD binary consisting of a hybrid CO/He WD and a companion He WD an excellent candidate for the progenitor of RCB stars.

79 ASTRONOMY AND ASTROPHYSICS↗

Axisymmetric Radiative Transfer Models of Kilonovae

The detailed observations of GW170817 proved for the first time directly that neutron star mergers are a major production site of heavy elements. The observations could be fit by a number of simulations that qualitatively agree, but can quantitatively differ (e.g., in total r-process mass) by an order of magnitude. We categorize kilonova ejecta into several typical morphologies motivated by numerical simulations, and apply a radiative transfer Monte Carlo code to study how the geometric distribution of the ejecta shapes the emitted radiation. We find major impacts on both spectra and light curves. The peak bolometric luminosity can vary by two orders of magnitude and the timing of its peak by a factor of five. These findings provide the crucial implication that the ejecta masses inferred from observations around the peak brightness are uncertain by at least an order of magnitude. Mixed two-component models with lanthanide-rich ejecta are particularly sensitive to geometric distribution. A subset of mixed models shows very strong viewing angle dependence due to lanthanide "curtaining," which persists even if the relative mass of lanthanide-rich component is small. The angular dependence is weak in the rest of our models, but different geometric combinations of the two components lead to a highly diverse set of light curves. We identify geometry-dependent P Cygni features in late spectra that directly map out strong lines in the simulated opacity of neodymium, which can help to constrain the ejecta geometry and to directly probe the r-process abundances.

79 ASTRONOMY AND ASTROPHYSICS↗

“Max-tech FHR”: What is the maximum technically achievable water delivery capacity from a single 120 V plug?

Electrifying water heating will require affordable, direct drop-in replacements for existing gas water heaters that meet consumer expectations for hot water delivery. For many replacement scenarios, affordability requires a replacement option with like-for-like physical product dimensions, and existing consumer expectations presume a large delivery capacity from a small water heater. This work explores the technical potential for hot water delivery capacity for a heat pump water heater (HPWH) powered from a common 120 V plug. Delivery capability in this work is quantified by two metrics: the First Hour Rating (FHR) defined in the Code of Federal Regulations, and a real-world draw profile that resulted in hot water runout during field studies of 120 V HPWHs. The work focuses on the common “tall-and-slim” gas water heater product category, with external product dimensions of 20 inch diameter and 60 inch height. Using a simulation model in DOE/ORNL Heat Pump Design Model, validated by experimental baseline testing of a commercially available 120 V HPWH model, several measures were evaluated to maximize the delivery capacity of a 120 V HPWH. The 120 V HPWH was constrained by 1) external product dimensions of 20 in diameter and 60 inch height, 2) electrical power of 1.5 kW (85% of the available power from a dedicated 120 V, 15 A circuit), and 3) commercially-available components. The measures implemented included a large compressor, a pumped loop with plate heat exchanger to maximize tank stratification, enlarged heat exchangers, and relatively thin insulation to maximize the water storage volume.

Gluesenkamp, Kyle↗

Performance Analysis of Speculative Parallel Adaptive Local Timestepping for Conservation Laws

Stable simulation of conservation laws, such as those used to model fluid dynamics and plasma physics applications, requires the satisfaction of the so-called Courant-Friedrichs-Lewy condition. By allowing regions of the mesh to advance with different timesteps that locally satisfy this stability constraint, significant work reduction can be attained when compared to a time integration scheme using a single timestep size. However, parallelizing this algorithm presents considerable difficulty. Since the stability condition depends on the state of the system, dependencies become dynamic and potentially non-local. In this article, we present an adaptive local timestepping algorithm using an optimistic (Timewarp-based) parallel discrete event simulation. We introduce waiting heuristics to limit misspeculation and a semi-static load balancing scheme to eliminate load imbalance as parts of the mesh require finer or coarser timesteps. Last, we outline an interface for separating the physics of the specific conservation law from the temporal integration allowing for productive adoption of our proposed algorithm. We present a misspeculation study for three conservation laws, demonstrating both the productivity of the local timestepping API, for which 74% of the lines of code are reused across different conservation laws, and the robustness of the waiting heuristics—at most 1.5% of element updates are rolled back. Our performance studies demonstrate up to a 2.8× speedup versus a baseline unoptimized local timestepping approach, a 4x improvement in per-node throughput compared to an MPI parallelization of synchronous timestepping, and scalability up to 3,072 cores on NERSC’s Cori Haswell partition.

97 MATHEMATICS AND COMPUTING↗

NuGrid stellar data set – III. Updated low-mass AGB models and s-process nucleosynthesis with metallicities Z= 0.01, Z = 0.02, and Z = 0.03

ABSTRACT The production of the neutron-capture isotopes beyond iron that we observe today in the Solar system is the result of the combined contribution of the r-process, the s-process, and possibly the i-process. Low-mass asymptotic giant branch (AGB) (1.5 < M/M⊙ < 3) and massive (M > 10 M⊙) stars have been identified as the main site of the s-process. In this work we consider the evolution and nucleosynthesis of low-mass AGB stars. We provide an update of the NuGrid Set models, adopting the same general physics assumptions but using an updated convective-boundary-mixing model accounting for the contribution from internal gravity waves. The combined data set includes the initial masses MZAMS/M⊙ = 2, 3 for Z = 0.03, 0.02, 0.01. These new models are computed with the mesa stellar code and the evolution is followed up to the end of the AGB phase. The nucleosynthesis was calculated for all isotopes in post-processing with the NuGrid mppnp code. The convective-boundary-mixing model leads to the formation of a 13C-pocket three times wider compared to the one obtained in the previous set of models, bringing the simulation results now in closer agreement with observations. Using these new models, we discuss the potential impact of other processes inducing mixing, like rotation, adopting parametric models compatible with theory and observations. Complete yield data tables, derived data products, and online analytic data access are provided.

79 ASTRONOMY AND ASTROPHYSICS↗

Simulated Fission Gas Behavior in Silicide Fuel at LWR Conditions

As a promising candidate for the accident tolerant fuel (ATF) used in light water reactors (LWRs), the fuel performance of uranium silicide (U 3 Si 2 ) at LWR conditions needs to be well-understood. However, existing experimental post-irradiation examination (PIE) data are limited to the research reactor conditions, which involve lower fuel temperature compared to LWR conditions. This lack of appropriate experimental data significantly affects the development of fuel performance codes that can precisely predict the microstructure evolution and property degradation at LWR conditions, and therefore evaluate the qualification of U 3 Si 2 as an AFT for LWRs. Considering the high cost, long timescale, and restrictive access of the in-pile irradiation experiments, this study aims to utilize ion irradiation to simulate the inpile behavior of the U 3 Si 2 fuel. Both in situ TEM ion irradiation and ex situ high-energy ATLAS ion irradiation experiments were employed to simulate different types of microstructure modifications in U 3 Si 2 . Multiple PIE techniques were used or will be used to quantitatively analyze the microstructure evolution induced by ion irradiation so as to provide valuable reference for the development of fuel performance code prior to the availability of the in-pile irradiation data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Production and characterization of a nitrogen-implanted Fe standard to calibrate PIGE measurements

Three calibration standard was produced by ion implantation of nitrogen in samples of Armco iron (99.7% iron). The samples was irradiated with nitrogen ion beams at several different energies (between 4 keV and 40 keV), and the ion doses were adjusted to obtain an uniform depth profile, using simulations with SRIM code. Two standards, one thick and other a foil (1.62mg/cm{sup 2}), was irradiated at same time with total nominal dose of 6.6×10{sup −16} atoms/cm{sup 2} distributed in a region of 100 nm in depth, with an average concentration of 9.0% nitrogen in iron. The third sample uses the same profile, but with a small dose, 1.1×10{sup −16} atoms/cm{sup 2} and average concentration of 1.5% nitrogen. The characterization of the implanted samples was done using RBS and NRA techniques to quantification of nitrogen.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

New risk metrics and mathematical tools for risk analysis: Current and future challenges

The current status of the food safety supply world wide, has led Food and Agriculture Organization (FAO) and World Health Organization (WHO) to establishing Risk Analysis as the single framework for building food safety control programs. A series of guidelines and reports that detail out the various steps in Risk Analysis, namely Risk Management, Risk Assessment and Risk Communication is available. The Risk Analysis approach enables integration between operational food management systems, such as Hazard Analysis Critical Control Points, public health and governmental decisions. To do that, a series of new Risk Metrics has been established as follows: i) the Appropriate Level of Protection (ALOP), which indicates the maximum numbers of illnesses in a population per annum, defined by quantitative risk assessments, and used to establish; ii) Food Safety Objective (FSO), which sets the maximum frequency and/or concentration of a hazard in a food at the time of consumption that provides or contributes to the ALOP. Given that ALOP is rather a metric of the public health tolerable burden (it addresses the total ‘failure’ that may be handled at a national level), it is difficult to be interpreted into control measures applied at the manufacturing level. Thus, a series of specific objectives and criteria for performance of individual processes and products have been established, all of them assisting in the achievement of FSO and hence, ALOP. In order to achieve FSO, tools quantifying the effect of processes and intrinsic properties of foods on survival and growth of pathogens are essential. In this context, predictive microbiology and risk assessment have offered an important assistance to Food Safety Management. Predictive modelling is the basis of exposure assessment and the development of stochastic and kinetic models, which are also available in the form of Web-based applications, e.g., COMBASE and Microbial Responses Viewer), or introduced into user-friendly softwares, (e.g., Seafood Spoilage Predictor) have evolved the use of information systems in the food safety management. Such tools are updateable with new food-pathogen specific models containing cardinal parameters and multiple dependent variables, including plate counts, concentration of metabolic products, or even expression levels of certain genes. Then, these tools may further serve as decision-support tools which may assist in product logistics, based on their scientifically-based and “momentary” expressed spoilage and safety level.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Measuring key Sm isotope ratios in irradiated UO 2 for use in plutonium discrimination nuclear forensics

As part of an experimental validation of a nuclear forensics methodology for Pu source reactor-type discrimination, destructive analysis has been performed here on two irradiated UO 2 pellets with different irradiation histories. Analysis has focused on measuring key Sm fission product isotope ratios used in a previously published maximum likelihood methodology to determine the most likely irradiation history of the pellets. A total of 21 Sm isotope ratios were measured within the irradiated pellets, and generally agreed within 20% of the irradiations as simulated using the Monte Carlo Radiation Transport and material depletion code, MCNP6. Results indicate the chosen approach can accurately measure the isotope ratios within 5% experimental error.

07 ISOTOPE AND RADIATION SOURCES↗

LDES-Sizing-and-siting-model (Power grid test cases for long-duration energy storage (LDES) siting) [SWR-25-75]

This repository contains code and data for performing a siting analysis of long-duration energy storage (LDES) in the 5-bus and RTS systems using the Sienna suite developed by the National Renewable Energy Laboratory for production cost modeling (PCM). This analysis involves moving the LDES component to different buses in the system, running a simulation, and considering the production cost of the simulation. Different system configurations are also analyzed with these scripts (such as moving load or renewable dispatch generators to different buses) to observe the impacts the system configuration has on optimal siting. This repository contains three sub directories discussed below. Both the 5-bus and RTS systems have two different initial configurations for the renewable energy components in them, one that is predominantly PV-driven and one that is predominantly wind-driven. Sienna suite can be found here: https://github.com/NREL-Sienna

Cole, David [University of Wisconsin]↗