Search NASA⌕ Search

SEARCH · Search NASA

Results for “predictive 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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Accuracy Enhancement of Nuclear Power Plant Simulators Utilizing High Accuracy Simulation Predictions

More recently, reactor core simulators for core designs associated with commercial nuclear power plants that utilize what is believed to be higher fidelity models have been developed. Features such as neutronics models that utilize transport equation solvers with fine spatial meshes and many energy-groups, thermal-hydraulic models that utilize sub-channel solvers with fine spatial mesh and capable of treating a wide range of fluid conditions, and fuel-coolant chemistry interaction models capable of treating CRUD deposition are to be found in these higher fidelity core simulators. These reactor core simulators require access to higher performance computers, characterized by many processors, cores and large memory. So associated with utilization of these simulators is access to high performance computers and ability to accommodate in one’s workflow longer execution times. By contrast, currently used core simulators by the nuclear industry can execute on engineering workstations and have execution times of seconds to minutes. The desirability for having short execution times is not only desired for support of time critical tasks but supports the mental process of decision making by engineers. The goal of the work reported upon here has the objective of retaining the fidelity of higher fidelity models while retaining the ability to utilize engineering workstations. Beyond the core simulator goal, additional goals of this work include incorporating the just described core simulator capability into a Nuclear Steam Supply System (NSSS) simulator, and to incorporate the resulting capability into an environment supportive of design and operational decision making associated with nuclear power stations. The model selected for the core neutronics model is the NESTLE code, for the core thermal-hydraulic model is the CTF code utilizing coarse mesh, and for the NSSS model is the RELAP5-3D code. WSC’s proprietary 3KEYMASTERTM platform is being used to provide software coupling, user interface, visualization, and reporting. The NESTLE core neutronics simulator was first integrated with the CTF core thermal-hydraulic simulator using CTF developed communication commands which are also used for CTF to communicate with RELAP5-3D under WSC’s proprietary 3KEYMASTERTM platform. To assure NESTLE prediction consistency with higher fidelity core neutronic simulators, buffer codes have been created to automatically generate from output files written by the VERA core simulator the NESTLE nodal neutronic parameter’ library, geometry, and pin-power reconstruction input files, thereby avoiding a number of challenges associated with utilizing lattice physics codes and providing consistency with VERA predictions. To treat absorber rod effects a multi-set library is utilized, where a set refers to a specific absorber rod fully inserted pattern. A coarse spatial mesh CTF model was developed with features added that support using CTF as envisioned in the engineering quality simulator. A hybrid meshing approach was implemented to allow for automated construction of models with mixed levels of refinement. Specifically, a core model could resolve some assemblies at a nodal level (4 subchannels per assembly) and others at a pin-resolution (one subchannel per coolant subchannel in the assembly). The intention is that this will allow for better resolution of limiting conditions such as DNBR and PCT, which are based on local rod and subchannel conditions. Further development was done of features that enhance the capabilities for the envisioned engineering quality simulator that has been developed, but now for RELAP-3D. The RELAP5-3D code development includes ability to model more than 999 components and the addition of the cross-channels turbulence mixing model and the void drift model that are implemented in CTF, aiming to achieve closer prediction agreement of the two codes for transient simulations, specifically, more accurate matches of the overall mass, momentum, and energy exchanges of both the liquid and gas phases between the neighboring core assemblies. Graphics were also developed for the Instructor Station for this project under WSC’s proprietary 3KEYMASTERTM platform to facilitate design and operational decision making.

42 ENGINEERING↗
DOE OSTI2026

Quantitative measurements of dislocations in metals for advancing predictive simulations

LLNL applications require scientists to predict how materials evolve under various thermomechanical conditions. While this is achieved through physics-based simulations, uncertainty in the predictions of mechanical properties remains a serious challenge that limits the predictive capabilities of models because we lack methods to compare predictions of atomic-scale defects (dislocations) with experimental measurements. High energy X-ray diffraction (HEXRD) is the most relevant technique that can provide the necessary statistical information on dislocations. However, this technique is not yet quantitative because we lack a precise understanding of the relationship between X-ray diffraction patterns and the underlying material dislocation content and arrangements. To address this need, we used our novel computational X-ray diffraction method to simulate the effect of dislocations on the diffraction patterns. We compared virtual and experimental diffraction patterns. Results allowed us to clearly establish the relationship between X-ray diffraction patterns and the underlying dislocation structures, proving that it is feasible to quantitatively measure dislocation statistics with HEXRD. This project delivered a method that can provide the missing piece to LLNL’s mechanical property simulations in advanced metals by obtaining experimentally long-needed quantitative dislocation data, which could fully enable predictive capabilities.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Toward an improved nonlocal thermodynamic equilibrium model for more predictive simulations of ignition scale hohlraums

Recently, nonlocal thermodynamic equilibrium (NLTE) modeling has been identified as the primary reason for discrepant predictions of the peak neutron production time in indirectly driven inertial confinement fusion (ICF) platforms. It has also been observed that predictions of collisional excitation rates differ by as much as 50% from measurements. Theoretical uncertainties in dielectronic recombination rates have also been posited as possibly contributing to errors in NLTE predictions. This work examines the impact of multipliers on collisional excitation and dielectronic recombination rates on simulations of a directly driven gold sphere and an indirect drive ICF implosion. It is found that multipliers on the collisional excitation rates have a strong impact on radiant intensity and electron temperature and a weaker impact on ionization state, whereas multipliers on dielectronic recombinations rates strongly impact ionization state with a smaller impact on radiant intensity and electron temperature. A self-consistent NLTE model which places multipliers on differing transitions, as motivated by experimental measurements and more detailed atomic physics predictions, improves agreement but does not completely eliminate discrepancies with measurements of the radiant intensity within the 2–4 keV spectral range.

Farmer, W. A. [Lawrence Livermore National Laborat↗
DOE OSTI2025

Plateau to River Model Predictive Simulations for All Ensemble Realizations to Support Modeling Work in Fiscal Year 2025

The purpose of this environmental calculation file (ECF) is to document predictions of flow and hydraulic head on the Central Plateau of the Hanford Site using the Plateau-to-River (P2R) Model (CP-57037, Model Package Report for the Plateau-to-River Model: Version 9.1). This calculation documents the simulation of the groundwater for the parent model domain of the P2R Model as a basis for use in other applications of the P2R Model. This application is unique from the standpoint that it will simulate all ensemble member models of the P2R Model whereas other applications may only utilize specific ensemble members. These simulations provide results that can be used in the process of selecting an appropriate subset of ensemble members for other applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Simulative Prediction of Solar Illuminance and Application of the Du-Sharples Model in Estimating Adapted Daylighting Metrics for an Urban Environment

The practice of daylighting in indoor spaces can significantly reduce electricity consumption and carbon emissions, improve human productivity, and enhance mood and cognitive perception. This work discussed the recent developments in daylighting science and practice, computed the periodic variations in average diurnal daylight levels for each month, quantified in terms of global horizontal illuminance and diffuse horizontal illuminance, for Kolkata, India, a city with tropical wet and dry climate, with two empirical luminous efficacy models of estimating solar illuminance, and assessed daylighting metrics with the Du-Sharples model. A program was formulated that could compute and generate daylight data with monthly-hourly solar irradiation data and the Du-Sharples model was utilized to predict dirt-corrected daylighting metrics for three glazing transmittance values and five elemental carbon deposition levels on glazing material. The highest monthly average global horizontal illuminance is recorded in April (64.05 klx for Littlefair model and 66.82 klx for Muneer-Kinghorn model) and the highest monthly average diffuse horizontal illuminance is recorded in July (33.23 klx for Littlefair model and 30.63 klx for Muneer-Kinghorn model). Further, the computed yearly average global and diffuse horizontal illuminance levels agree well with a previous study that applied the Perez model. Yearly average horizontal work surface illuminance level remained >1.5 klx for window-towall area ratio >30 %. The approach adopted in this work and the temporal variation charts of computed exterior daylight level data may assist building service engineers, architects, and indoor lighting practitioners in making informed policy decisions at different stages of building planning.

Engineering↗
DOE OSTI2025

High Energy Density Physics of Inertial Confinement Fusion Ablator Materials

The historic December 5, 2022 experiment at Lawrence Livermore National Lab’s (LLNL) National Ignition Facility (NIF) reached fusion energy ignition for the first time. This is the most important scientific breakthrough of the 21st century paves the way to future clean inertial fusion energy (IFE). The diamond (high density carbon (HDC)) ablator material used in this experiment displays detrimental effects due to the development of hydrodynamic instabilities at the diamond/fuel interface under shock compression. New alternatives to diamond ablators are required to step up the energy yield in ICF experiments. The unique combination of mechanical strength (approaching that of diamond), the ability to accommodate high-Z dopants (in contrast to diamond), and the tunability of the properties (through synthesis material with varying sp 3 content) make amorphous carbon (a-C) a promising material for next-generation IFE ablative capsules. However, despite its critical importance to the IFE program, the behavior of a-C carbon at extreme temperatures and pressures remains largely unexplored. The primary goals of this project were to perform groundbreaking dynamic compression experiments and predictive simulations to uncover the fundamental high-energy-density physics of amorphous carbon. Our goals were (1) to uncover the metastability range of amorphous carbon and probe phase transitions to diamond or metastable supercooled liquid carbon; (2) to acquire high-quality equation of state (EOS) data and develop an experimentally validated EOS from machine-learning MD simulations of the complex states of carbon; and (3) to uncover the complex behavior of carbon liquid in both thermodynamically stable and metastable supercooled states by accessing large areas of carbon phase diagram with amorphous samples with variable sp 3 content. Our proposed experimental program included measurements of equation of state and diffraction measurements using the Omega EP laser at the Laboratory of Laser Energetics at the University of Rochester. The theoretical/simulation program involved the development of machine-learning models of the complex response of amorphous carbon under dynamic compression by performing molecular dynamics simulations at experimental time and length scales using leadership class DOE supercomputers. Simulations guided experiments to observe predicted phenomena and acquire critical experimental data in specific pressure-temperature domains to validate theoretical models. This research delivered fundamental properties of novel amorphous carbon IFE ablator material, including phase diagram and EOS. These results will aid in IFE target design and implosion experiments. A unique combination of predictive simulations and dynamic and static experiments provided a highly inspirational intellectual environment for graduate students and postdocs involved in this project.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A Cellular Automaton Simulation for Predicting Phase Evolution in Solid-State Reactions

New computational tools for solid-state synthesis recipe design are needed in order to accelerate the experimental realization of novel functional materials proposed by high-throughput materials discovery workflows. This work contributes a cellular automaton simulation framework for predicting the time-dependent evolution of intermediate and product phases during solid-state reactions as a function of precursor choice and amount, reaction atmosphere, and heating profile. The simulation captures the effects of reactant particle spatial distribution, particle melting, and reaction atmosphere. Reaction rates based on rudimentary kinetics are estimated using density functional theory data from the Materials Project and machine learning estimators for the melting point and the vibrational entropy component of the Gibbs free energy. The resulting simulation framework allows for the prediction of the likely outcome of a reaction recipe before any experiments are performed. We analyze five experimental solid-state recipes for BaTiO 3 , CaZrN 2 , and YMnO 3 found in the literature to illustrate the performance of the model in capturing reaction selectivity and reaction pathways as a function of temperature and precursor choice. This simulation framework offers an easier way to optimize existing recipes, aid in the identification of intermediates, and design effective recipes for yet unrealized inorganic solids in silico .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Performance prediction applying different reduced turbulence models to the SMART tokamak

The SMall Aspect Ratio Tokamak (SMART) is currently being commissioned at the University of Seville and will be able to compare the performance of positive and negative triangularity plasmas at low aspect ratio. Predictive simulations have been performed for different machine scenarios and heating schemes using the TRANSP code. The objectives of these simulations are to predict the parameters expected in positive triangularity plasmas, to guide diagnostic development, and to validate transport models. Several reduced turbulence models have been used to predict electron and ion temperatures for the operational phase 2. All models provide similar results from approximately mid-radius to the separatrix but important discrepancies are found in the core region. These positive triangularity results are compared with experiments from a similar size machine like GLOBUS-M2. The multi-mode model (MMM) shows the best agreement. Simulations with different boundary conditions have been performed and no strong differences have been observed between them. The impact of neutral beam injection (NBI) on the predicted profiles has also been addressed. Rotation reduces turbulence levels so higher temperatures are achieved when included in the simulations. Studying the different contributions to the thermal diffusivities, it is observed that electron temperature gradient (ETG) turbulence dominates at the plasma core while micro-tearing modes (MTM) dominate at the edge in the electron channel. In the ion channel, the neoclassical contribution is dominant at the core and at the very edge while the Weiland component, which includes ion temperature gradient mode (ITG), trapped electron mode (TEM), kinetic ballooning mode (KBM), peeling mode (PM) and collisionless and collision dominated magnetohydrodynamic (MHD) modes governs the mid-radius region. For phase 3, two plasmas with different electron densities have been studied. The case with lower density matches well a specific discharge of GLOBUS-M2. The higher density plasma shows high performance with β N ≈ 3.8.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Modeling Dust Production, Growth, and Destruction in Reionization-era Galaxies with the CROC Simulations. II. Predicting the Dust Content of High-redshift Galaxies

Abstract We model the interstellar dust content of the reionization era with a suite of cosmological, fluid-dynamical simulations of galaxies with stellar masses ranging from ∼10 5 to 10 9 M ⊙ in the first 1.2 Gyr of the Universe. We use a post-processing method that accounts for dust creation and destruction processes, allowing us to systematically vary the parameters of these processes to test whether dust-dependent observable quantities of galaxies at these epochs could be useful for placing constraints on dust physics. We then forward model observable properties of these galaxies to compare to existing data. We find that we are unable to simultaneously match existing observational constraints with any one set of model parameters. Specifically, the models that predict the largest dust masses D / Z ≳ 0.1 at z = 5—because of high assumed production yields and/or efficient growth via accretion in the interstellar medium—are preferred by constraints on total dust mass and infrared (IR) luminosities, but these models produce far too much attenuation in the ultraviolet (UV), preventing them from matching observations of β UV . To investigate this discrepancy, we analyze the relative spatial distribution of stars and dust as probed by IR and UV emission, which appear to exhibit overly symmetric morphologies compared to existing data, likely due to the limitations of the stellar feedback model used in the simulations. Our results indicate that the observable properties of the dust distribution in high redshift galaxies are a particularly strong test of stellar feedback.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Predictive ALEGRA Simulations of the Mach Wave Concept on the Z Machine

This high level presentation introduces the Sandia Z machine, dynamic materials properties experiments including the Mach Wave concept, and Alegra magnetohydrodynamic modeling for a technical audience. Alegra models are shown to be highly capable, with predictive capabilities limited only by the quality and accuracy of material models.

Boerner, Jeremiah James [Sandia National Laborator↗
DOE OSTI2025

Meshfree simulation and prediction of recrystallized grain size in friction stir processed 316L stainless steel

Friction stir processing (FSP) is a promising solid-phase microstructural modification technique that can repair and enhance damaged stainless steel surfaces exposed to harsh environments. The quality of the repaired material is closely correlated to the recrystallized grain size in the stir zone (SZ), which is influenced by the thermomechanical conditions dictated by FSP process parameters. Thus, establishing a reliable relationship between these parameters and recrystallized grain size in the SZ is crucial for optimizing repair quality. However, existing experimental approaches often rely on indirect temperatures measured far from the SZ, along with rough strain rate estimations, which are imprecise and time-consuming. Meanwhile, existing mesh-based modeling methods usually face numerical challenges when dealing with the large material deformations inherent in FSP. Here, to address these issues, this study introduces a meshfree process model for FSP based on the smoothed particle hydrodynamics (SPH) method, aimed at predicting process conditions under different parameters. The model is validated using experimental data from 11 combinations of tool traverse and rotation speeds on 316 L stainless steel. Correlations between process parameters, material flow, temperature, strain, strain rate, and recrystallized grain size are revealed through SPH simulations and electron backscatter diffraction (EBSD) imaging. The results show that in situ SZ temperatures range from 1071 to 1322°C, which exceed the tool temperature by over 300°C. Furthermore, SZ temperature, strain rate, and grain size increase monotonically with higher tool temperature and faster traverse speed. A relationship is then established between the model-predicted Zener-Hollomon parameter and the recrystallized grain size based on EBSD data, expressed as ln(d) = -0.364 ln(Z) + 14.673. Finally, this relationship exhibits satisfactory accuracy with errors of less than 26.9% in predicting grain sizes at various SZ locations, which offers valuable insights for optimizing FSP repair processes for 316 L stainless steel.

316L stainless steel↗
DOE OSTI2025

Uncertainties in the effects of organic aerosol coatings on polycyclic aromatic hydrocarbon concentrations and their estimated health effects

We used the CAM5 model to examine how different particle-bound polycyclic aromatic hydrocarbon (PAH) degradation approaches affect the spatial distribution of benzo(a)pyrene (BaP). Three approaches were evaluated: NOA (no effect of OA coatings state on BaP), shielded (viscous OA coatings shield BaP from oxidation under cool and dry conditions) and ROI-T (viscous OA coatings slow BaP oxidation in response to temperature and humidity). Results show that BaP concentrations vary seasonally, influenced by emissions, deposition, transport and degradation approach, all of which are influenced by meteorological conditions. All simulations predict higher population-weighted global average (PWGA) fresh BaP concentrations during December–January–February (DJF) compared to June–July–August (JJA), due to increased emissions from household activities and reduced removal processes during colder months. The shielded and ROI-T approaches, which account for OA coatings, result in 2–6 times higher BaP concentrations in DJF compared to NOA. The shielded simulation predicts the highest PWGA fresh BaP concentration (1.3 ng m −3 ), with 90 % of BaP protected from oxidation. In contrast, the ROI-T approach forecasts lower concentrations in middle to low latitudes, as it assumes less effective OA coatings under warmer, more humid conditions. Evaluations against observed BaP concentrations show the shielded approach performs best, with a normalized mean bias (NMB) within ± 20 %. The combined incremental lifetime cancer risk (ILCR) for both fresh and oxidized PAHs is similar across simulations, emphasizing the importance of considering both forms in health risk assessments. This study highlights the critical role of accurate degradation approaches in PAH modeling.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2025

Magnetohydrodynamic simulations preliminarily predict the habitability and radio emission of TRAPPIST-1e

TRAPPIST-1e, an Earth-sized exoplanet in the habitable zone of the nearby M dwarf TRAPPIST-1, may experience magnetospheric responses that vary with stellar space weather, which could potentially influence both its habitability and radio emissions. Our objective is to investigate how different Earth-like magnetospheric configurations of TRAPPIST-1e – specifically variations in dipolar magnetic field strength and axial tilt – respond to diverse stellar space weather conditions, including events analogous to coronal mass ejections (CMEs), and to assess their implications for potential habitability and expected radio emissions. We conducted 3D magnetohydrodynamic simulations of the TRAPPIST-1e system using the PLUTO code in spherical coordinates. The planetary magnetic field was modelled as dipolar, with equatorial strengths from Earth-like to several times stronger. The dipole axis spans a representative range of axial tilts. We investigate four stellar wind environments, from sub-Alfvénic flow to CME-like disturbances. Planetary shielding was quantified based on the magnetopause standoff distance, and radio powers were estimated via empirical scaling laws. Our simulations show that both shielding and radio power depend strongly on the magnetic configuration. Stronger fields increase protection, while larger tilts reduce it. Radio power increases with both field strength and tilt across all wind regimes. An Earth-like magnetic field can provide effective shielding even under intense CMEs, whereas high tilts require stronger fields. Predicted radio powers reach ~10 20 erg s −1 during CMEs, making bursts from close-in, magnetised planets more detectable. However, for TRAPPIST-1e, the maximum cyclotron frequency lies below the Earth’s ionospheric cutoff (~10 MHz), making ground-based detection currently infeasible.

Sun: magnetic fields↗
DOE OSTI2024

Orbital-Free Quantum Simulation Methods for Application to Warm Dense Matter (Final Technical Report)

Predictive simulations for prediction of condensed system behavior in state conditions far from ambient is increasingly crucial to DOE priorities. Warm dense matter (WDM) is the paradigm: temperature T > 1-15 eV, pressures P to 1 Mbar or greater. Experiments under such state conditions are difficult and costly. We summarize work driven by the need and opportunity to make free-energy density functional theory (DFT) as powerful a tool for ab initio simulation of matter under such extreme conditions as ground state DFT is for ordinary matter Advancing orbital-free DFT (OF-DFT) to eliminate the Kohn-Sham (KS) scaling bottleneck in such simulations is the other priority. The concurrent challenge for both goals is the intrinsic complexity of WDM. We summarize 15 years of successes and major progress on (1) free energy exchange-correlation functionals; (2) non-interacting free energy functionals (counterpart to T=0 Kohn-Sham kinetic energy density functionals); (3) rigorous results and constraints for free-energy DFT; (4) software for free energy DFT calculations in both conventional Kohn-Sham and OF-DFT form; (5) de-orbitalization of advanced orbital-dependent ground state functionals for use in OF-DFT; (6) demonstration calculations; (7) ancillary achievements (e.g. major review articles, secondary explorations motivated by primary goals).

36 MATERIALS SCIENCE↗
DOE OSTI2024

Versatile stochastic model for predictive KMC simulation of fcc metal nanostructure evolution with realistic kinetics

Stochastic lattice-gas models provide the natural framework for analysis of the surface diffusion-mediated evolution of crystalline metal nanostructures on the appropriate time scale (often 10 1 –10 4 s) and length scale. Model behavior can be precisely assessed by kinetic Monte Carlo simulation, typically incorporating a rejection-free algorithm to efficiently handle the broad range of Arrhenius rates for hopping of surface atoms. The model should realistically prescribe these rates, or the associated barriers, for a diversity of local surface environments. However, commonly used generic choices for barriers fail, even qualitatively, to simultaneously describe diffusion for different low-index facets, for terrace vs step edge diffusion, etc. We introduce an alternative Unconventional Interaction–Conventional Interaction formalism to prescribe these barriers, which, even with few parameters, can realistically capture most aspects of behavior. Here, the model is illustrated for single-component fcc metal systems, mainly for the case of Ag. It is quite versatile and can be applied to describe both the post-deposition evolution of 2D nanostructures in homoepitaxial thin films (e.g., reshaping and coalescence of 2D islands) and the post-synthesis evolution of 3D nanocrystals (e.g., reshaping of nanocrystals synthesized with various faceted non-equilibrium shapes back to 3D equilibrium Wulff shapes).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Simulation of Deuterium Oxidation in a Room with a Fire

A series of experiments where deuterium was released in trace amounts into a room with a fire were performed at Sandia National Laboratories. This report describes the corresponding effort to model the test series using SIERRA/FUEGO. The objective of this modeling effort was to produce a simulation test matrix that can be utilized to help interpret the corresponding experiments and be used to assess the credibility of using the SIERRA/FUEGO simulations as a surrogate for real tritium reaction data for mock fire and release studies. The simulations assessed multiple configurations of fire size and location and overall found that there was low conversion of deuterium to oxidized deuterium. There are some notable differences between the simulation predictions and the experimental results, such as the simulations generally predicting higher amounts of deuterium conversion compared to the experiment. A strategy to address modeling uncertainties and unknowns is also given in this report.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗