Search NASA⌕ Search

SEARCH · Search NASA

Results for “core 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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Multiphysics Demonstration of Temperature-Driven Assembly Bowing in SFRs using MOOSE-Based Codes

Core bowing is an important passive safety mechanism in liquid metal cooled fast reactors. When the core restraint system is properly designed, temperature and flux gradients influence assemblies in the core to bow into less reactive configurations during accident scenarios, resulting in negative reactivity feedback. Prediction of core bowing involves complex interplay of radiation transport, impacts of fluid flow and heat transfer on duct temperature, and mechanical responses to the induced temperature and flux gradients. Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program [1], an integrated multiphysics approach is being developed to model the core bowing phenomena in liquid metal-cooled fast reactors with the Multiphysics Object Oriented Simulation Environment (MOOSE) [2]. In this methodology, the MOOSE-based reactor physics code Griffin [3] will solve the neutron transport equation and determine the power distribution. With the detailed power distribution from Griffin, the subchannel analysis codes MOOSE-Subchannel [4] and Pronghorn [5] are utilized to calculate the assembly temperature distribution. MOOSE’s Solid Mechanics [6] and Contact [7] Modules are leveraged to calculate the thermal expansion and duct bowing displacement with the duct wall temperature from thermal hydraulics calculation. In this work, an initial one-way coupling demonstration of the integrated multiphysics approach has been performed on a seven-assembly problem based on the sodium-cooled fast reactor ABR-1000 design [8]. The neutronics calculation with Griffin is not yet involved in the current simulation. MOOSE-Subchannel and Pronghorn evaluate fluid and solid temperature based on a fixed power distribution. In addition, one-way coupling is utilized in this coupled calculation, via Pronghorn passing the duct temperature data to the MOOSE Solid Mechanics calculation. An assessment of the Solid Mechanics module was performed in parallel to verify duct bowing behavior with duct-to-duct contact phenomenon [9]. The displacement from MOOSE Solid Mechanics is not yet transferred back and utilized in the Pronghorn and MOOSE-Subchannel calculation. This model will be available on the National Reactor Innovation Center (NRIC) Virtual Test Bed (VTB) repository [10]. Future stages of this work will involve solving problems of increasing complexity as well as adding more physics (e.g. reactor physics) to the integrated workflow to reach the end goal of modeling the core bowing phenomenon with an integrated multiphysics workflow.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2014

UAS-Systems Integration, Validation, and Diagnostics Simulation Capability

As part of the Phase 1 efforts of NASA's UAS-in-the-NAS Project a task was initiated to explore the merits of developing a system simulation capability for UAS to address airworthiness certification requirements. The core of the capability would be a software representation of an unmanned vehicle, including all of the relevant avionics and flight control system components. The specific system elements could be replaced with hardware representations to provide Hardware-in-the-Loop (HWITL) test and evaluation capability. The UAS Systems Integration and Validation Laboratory (UAS-SIVL) was created to provide a UAS-systems integration, validation, and diagnostics hardware-in-the-loop simulation capability. This paper discusses how SIVL provides a robust and flexible simulation framework that permits the study of failure modes, effects, propagation paths, criticality, and mitigation strategies to help develop safety, reliability, and design data that can assist with the development of certification standards, means of compliance, and design best practices for civil UAS.

Buttrill, Catherine W.↗
DOE OSTI2025

Dislocation nano-hydrides in nickel: Nucleation, evolution and effects on dislocation behaviors

Nano-hydrides have been predicted to precipitate at the core of edge dislocations in the Ni-H system, a mechanism that may promote hydrogen embrittlement. However, nano-hydride nucleation, growth, and effects on dislocation behavior have seldom been explored. This work combines molecular dynamics grand canonical Monte Carlo (MD-GCMC) simulations and continuum modeling to uncover a wide range of phenomena linked to dislocation nano-hydrides. Simulations reveal that nano-hydrides can be stabilized at dislocation cores with all character angles, including screw segments, due to the hydrostatic stresses around the cores of the Shockley partials. Nano-hydride nucleation takes place in these regions, and growth is dictated by the character angle θ of the perfect dislocation. The equilibrium stacking fault width d eq varies dynamically to increase the local hydrostatic stress field and facilitate the formation of the nano-hydride, forming a constriction-like feature and leading to three distinct behaviors: d eq decreases for θ > 30°, d eq remains unchanged for θ = 30°, and d eq increases for θ < 30°. Remote hydrostatic and Escaig stresses are also shown to influence the nucleation stage, implying stress concentrations such as those ahead of crack tips may facilitate nano-hydride precipitation. Moreover, we identify a new hydrogen-induced 60° dislocation reaction that emits a Shockley partial on a conjugate plane, with potential implications for twin nucleation. Testable predictions from this study are then used to reinterpret previous results from the literature. In conclusion, these findings provide a comprehensive framework to assess nano-hydride formation and evolution at dislocations in nickel and other face-centered cubic metals, with important implications to hydrogen embrittlement.

Dislocations↗
DOE OSTI2026

From MAST-U to STEP: Power exhaust control challenges and opportunities

This work explores the challenges and opportunities for power exhaust control in STEP, informed by dedicated MAST-U experiments. The STEP system primarily relies on hydrogenic and argon gas injection into the divertors, compensating for transients originating in the plasma core. Fast transients (>10 Hz) that evolve too quickly for actuators to compensate must be absorbed passively; the enhanced transient buffering of the foreseen long-legged divertor provides a clear advantage over conventional geometries, as demonstrated in the MAST-U Super-X divertor. Experiments further indicate that STEP’s tight divertor baffling enables near-independent control of the upper and lower divertors. The most challenging transients are expected from power-sharing fluctuations, as MAST-U experiments observed extremely fast dynamics. Fluctuations from core pellet fuelling remain relatively benign due to the small pellet size according to simulations. Exhaust control in a reactor like STEP requires an integrated approach befitting its machine-critical nature. A predictive control element, integrated with core plasma control, allows pre-emptive preparation of the divertor for incoming transients. An observer-based diagnostic approach is foreseen to monitor the divertor in the challenging reactor environment, supported by an extended diagnostic set in the non-nuclear phase to validate the required dynamic models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Analysis of AP1000 Small-Break Loss-of-Coolant Accident Using Reactor Transient Simulator

The Westinghouse Electric Company’s Advanced Passive Reactor (AP1000) is characterized by the incorporation of passive safety systems (PSSs) designed to ensure core cooling during transient events. The assessment of PSSs requires evaluation of their performance through a combination of experiments and simulations employing various thermal-hydraulic codes. In addition, detailed evaluation of PSSs for a specific reactor system transient analysis such as loss-of-coolant-accident analysis supports understanding representative integral effects test facility development and the further evolution model development and assessment process. Developing a reactor system code is a complex and time-consuming process that requires significant engineering expertise and effort. It can take several months to even years to complete in the early stages of reactor system design and analysis. However, this process can be expedited through the use of transient simulator models for similar reactor systems, which can be used for lesson learning and training purposes. This study uses the Personal Computer Transient Analyzer (PCTRAN) code. The main advantage of PCTRAN is its ease of use and ability to run faster than real time. This study presents the results obtained for a small-break loss-of-coolant accident (SBLOCA) for two breaks using the full version (licensed) of PCTRAN. The purpose of this investigation is to evaluate the overall system behavior during the postulated SBLOCA event as well as assess the capability of the PCTRAN code to reproduce the system response during transient events. The obtained results were compared with the Westinghouse NOTRUMP system code. The PCTRAN code proved to be reliable in predicting the qualitative behavior of the system in both transient cases. As for the system response, it was found that it is contingent on the activation time of the PSSs. The differences in reactor coolant system pressure between the two codes were attributed to the critical flow model and simplification of mass and energy balance. Despite PCTRAN’s limitations, it can still provide a reasonable prediction of various reactor parameters such as pressure, mass flow rate, and void fraction during a SBLOCA scenario. It is worth noting that PCTRAN currently employs a bulk approach similar to that of the Modular Accident Analysis Program (MAAP) and MELCOR codes. However, the upcoming version of PCTRAN will include an artificial intelligence–based detection and accident prevention system, as well as different models for different reactor components. Consequently, PCTRAN has the potential to be upgraded to match the system thermal-hydraulic codes of the U.S. Nuclear Regulatory Commission and become more widely used in cybersecurity to safeguard nuclear power plants from cyberattacks.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2015

Implementation of an Open-Scenario, Long-Term Space Debris Simulation Approach

This paper provides a status update on the implementation of a flexible, long-term space debris simulation approach. The motivation is to build a tool that can assess the long-term impact of various options for debris-remediation, including the LightForce space debris collision avoidance concept that diverts objects using photon pressure [9]. State-of-the-art simulation approaches that assess the long-term development of the debris environment use either completely statistical approaches, or they rely on large time steps on the order of several days if they simulate the positions of single objects over time. They cannot be easily adapted to investigate the impact of specific collision avoidance schemes or de-orbit schemes, because the efficiency of a collision avoidance maneuver can depend on various input parameters, including ground station positions and orbital and physical parameters of the objects involved in close encounters (conjunctions). Furthermore, maneuvers take place on timescales much smaller than days. For example, LightForce only changes the orbit of a certain object (aiming to reduce the probability of collision), but it does not remove entire objects or groups of objects. In the same sense, it is also not straightforward to compare specific de-orbit methods in regard to potential collision risks during a de-orbit maneuver. To gain flexibility in assessing interactions with objects, we implement a simulation that includes every tracked space object in Low Earth Orbit (LEO) and propagates all objects with high precision and variable time-steps as small as one second. It allows the assessment of the (potential) impact of physical or orbital changes to any object. The final goal is to employ a Monte Carlo approach to assess the debris evolution during the simulation time-frame of 100 years and to compare a baseline scenario to debris remediation scenarios or other scenarios of interest. To populate the initial simulation, we use the entire space-track object catalog in LEO. We then use a high precision propagator to propagate all objects over the entire simulation duration. If collisions are detected, the appropriate number of debris objects are created and inserted into the simulation framework. Depending on the scenario, further objects, e.g. due to new launches, can be added. At the end of the simulation, the total number of objects above a cut-off size and the number of detected collisions provide benchmark parameters for the comparison between scenarios. The simulation approach is computationally intensive as it involves tens of thousands of objects; hence we use a highly parallel approach employing up to a thousand cores on the NASA Pleiades supercomputer for a single run. This paper describes our simulation approach, the status of its implementation, the approach to developing scenarios and examples of first test runs.

laser↗
DOE OSTI2026

Effects of density turbulence on helicon wave propagation in the core plasmas

Radio frequency wave propagation can be significantly affected by density irregularities, such as filaments in the scrape-off layer or instabilities in the core plasma. In this study, we examine the impact of edge turbulence on helicon wave propagation using the Petra-M simulation tool. To analyze the effect of edge turbulence, we utilize a realistic background plasma derived from XGC simulations, which includes spatial density fluctuations at the edge. This focus is particularly relevant for the DIII-D configuration, characterized by edge density turbulence of the core plasmas in a wide pedestal QH-mode. We focus on helicon wave propagation in the core plasma since the slow mode cannot propagate into the core plasma when the density is higher than the lower hybrid resonance, where a wide pedestal QH-mode occurs. The simulation results indicate that edge density fluctuations have a substantial impact on helicon wave coupling due to scattering in the core. We specifically demonstrate that the toroidal mode number and the level of density fluctuations are significant factors influencing scattering. While we have minimized the excitation of slow modes from the antenna, we still observe mode-converted slow modes resulting from the incoming helicon waves. The insights derived from these simulations will inform upcoming tokamak experiments regarding helicon antenna coupling in long pulse scenarios.

Kim, Eun-Hwa [Princeton Plasma Physics Laboratory ↗
DOE OSTI2026

Capturing Sub‐Kilometer Flood Inundation Dynamics During the California Rain‐on‐Snow Events of 2017

The severe impacts of rain‐on‐snow (ROS) extreme events have been widely recognized and studied. However, unlike hydrological processes, flood inundation dynamics and the relative contribution of rainfall and snowmelt during ROS events remain under‐investigated. We diagnosed and documented sub‐kilometer spatio‐temporal dynamics of flood inundation during the 2017 California ROS events, simulated by a 2‐dimensional hydrodynamic model, River Dynamical Core (RDycore). RDycore shows good performance in capturing fine‐scale flood inundation dynamics against gauge measurements (with a median correlation coefficient of 0.81) and satellite observations. On top of rainfall, snowmelt not only increases the mean maximum inundation depth (12.0%–25.1%), but also expands the total flooded area (19.9%–31.9%) and prolongs the mean flood duration (3.4%–7.1%) across the events. Our findings offer an explicit and accurate picture of when and where ROS flooding could occur and how snowmelt increases flood hazard, valuable for risk assessment and infrastructure planning.

Flood inundation↗
DOE OSTI2025

Efficient Neural Network Approaches for Conditional Optimal Transport with Applications in Bayesian Inference

In this work, we present two neural network approaches that approximate the solutions of static and dynamic conditional optimal transport (COT) problems. Both approaches enable conditional sampling and conditional density estimation, which are core tasks in Bayesian inference—particularly in the simulation-based (“likelihood-free”) setting. Our methods represent the target conditional distribution as a transformation of a tractable reference distribution. Obtaining such a transformation, chosen here to be an approximation of the COT map, is computationally challenging even in moderate dimensions. To improve scalability, our numerical algorithms use neural networks to parameterize candidate maps and further exploit the structure of the COT problem. Our static approach approximates the map as the gradient of a partially input convex neural network. It uses a novel numerical implementation to increase computational efficiency compared to state-of-the-art alternatives. Our dynamic approach approximates the conditional optimal transport via the flow map of a regularized neural ODE; compared to the static approach, it is slower to train but offers more modeling choices and can lead to faster sampling. We demonstrate both algorithms numerically, comparing them with competing state-of-the-art approaches, using benchmark datasets and simulation-based Bayesian inverse problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Mauka Energy FEVER Tool Dataset

Mauka Energy’s dataset, developed under the Forestry Electric Vehicle Energy Routing (FEVER) project and funded by the U.S. Department of Energy’s Small Business Innovation Research program, is a high-resolution geospatial resource designed to support energy modeling for electric log trucks in complex forestry environments. The dataset integrates detailed spatial and road network data to enable accurate simulation of vehicle performance across varied terrain. At its core, the dataset incorporates lidar-derived elevation models, road alignments, and surface classifications from Oregon State University’s McDonald-Dunn Research Forest. These data capture fine-scale variations in slope, curvature, and surface conditions across forest road systems, allowing for vehicle-level analysis of energy consumption and recovery. The dataset also includes data collected on the surrounding public and private road networks in Benton County, Oregon, used in real-world haul routes. These connecting segments provide critical context for modeling transitions between forest operations and regional transportation infrastructure, incorporating attributes such as grade profiles, elevation change, and speed constraints. This combined dataset underpins the development of Mauka Energy’s rolldown tool, which quantifies energy use and regenerative braking potential on downhill and variable-grade segments. By leveraging high-resolution terrain and road data, the FEVER project enables more accurate assessment of electric vehicle feasibility and performance in forestry applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Proton Beam Dynamics in Bare IOTA with Intense Space-Charge

We are commissioning a 2.5~MeV proton beam for the Integrable Optics Test Accelerator at Fermilab, allowing experiments in the strong space-charge regime with incoherent betatron tune shifts nearing 0.5. Accurate modelling of space-charge dynamics is vital for understanding planned experiments. We compare anticipated emittance growth and beam loss in the bare IOTA configuration using transverse space-charge models in Xsuite, PyORBIT, and MAD-X simulation codes. Our findings reveal agreement within 30\% in core phase-space density predictions up to 100~synchrotron periods at moderate beam currents, while tail distributions and beam loss show significant differences.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Center for Integrated Simulation of Fusion Relevant RF Actuators

This project was part of the “Center for Integrated Simulation of Fusion Relevant RF Actuators” SciDAC-4 project, led by Dr. Paul Bonoli (MIT). Rather than use an acronym (CISFRRFA), the project will be referred to in this document as the “RF-SciDAC4”. The larger SciDAC-4 project goals were to: 1. Develop an integrated simulation of the antenna + sheath + scrape-off-layer + core plasma system which fully utilizes leadership class computing resources to move towards a quantitative predictive capability for the response to RF power. 2. Work closely with the SciDAC-4 Whole Device Modeling (WDM) community to make both our new code development efforts, as well as the established hierarchy of RF tools, available within their environment, and to utilize WDM technologies to implement the couplings below. 3. Validate this predictive capability on appropriately diagnosed experiments including dedicated RF test stands, linear devices, and existing tokamaks. 4. Use these tools to inform design of robust, impurity-mitigating RF heating and current drive sources for future fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Uncertainty quantification and sensitivity analysis of a nuclear thermal propulsion reactor startup sequence

The research presented in this article describes progress in applying stochastic methods, uncertainty quantification, parametric studies, and variance-based sensitivity analysis (also known as Sobol sensitivity analysis) to a full-core model of a nuclear thermal propulsion (NTP) system simulated via the radiation transport code Griffin to simulate neutronics. Our goal is to develop a reduced-order (surrogate) model that can be rapidly sampled with perturbations to multiple input parameters. In this NTP system, reactivity and power feedback affect the rotation of control drums (CDs), which is itself controlled by a hybrid proportional-integral-derivative (PID) controller actuated by the power demand and reactivity feedback from the numerical model. This model uses reactor kinetic feedback (mean generation time [Λ] and effective delayed neutron fraction [ β eff ] from a transient Griffin simulation executed via Griffin’s improved quasi-static solver to provide the kinetic parameters) as inputs to functions that control the CD rotation angle. By investigating numerous stochastic approaches, we developed a dual-purpose surrogate model of the NTP system, using polynomial regression in the Multiphysics Object-Oriented Simulation Environment (MOOSE) Stochastic Tools Module (STM). The trained model can be rapidly sampled while simultaneously perturbing various input parameters, such as coefficients on the PID control or temperature (directly affecting the neutron cross section). The surrogate model delivers accurate (within 5%) results at speeds orders of magnitude faster (minutes, not days of computational time) than the base model. Once the surrogate model has been trained, distributions of the uncertain parameters can be changed at will to investigate the effects of perturbing multiple inputs as well as the effects of these inputs on the model output. For example, coefficients used in the PID control system may vary due to some type of physical interference, or uncertainty may exist in the temperature of the neutron cross sections in various regions of the reactor. A distribution can be placed on these parameters, and operational boundaries can be determined. The goal of this work is to support development of an advanced control system for operating CDs in a functioning NTP system. This work is a scoping study of the MOOSE STM.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
NASA NTRS1972

Acoustic fatigue and sound transmission characteristics of a ram composite panel design

An experimental study to determine the acoustic fatigue characteristics of a flat multi-layered structural panel is described. The test panel represented a proposed design for the outer skin of a research application module to be housed within the space shuttle orbiter vehicle. The test specimen was mounted in one wall of the Wyle 100,000 cu ft reverberation room and exposed to a broadband acoustic environment having an overall level of 145 db. The test panel was exposed to nine separate applications of the acoustic environment, each application consisting of 250 seconds duration. Upon completion of the ninth test run, the specimen was exposed to a simulated micrometeoroid impact near the panel center. One additional test run of 250 seconds duration was then performed to complete the overall simulation of 50 flight missions. The experimental results show that no significant fatigue damage occurred until the test specimen was exposed to a simulated micrometeoroid impact. The intermediate foam layer forming the core of the test specimen suffered considerable damage due to this impact, causing a marked variation in the dynamic characteristics of the overall test panel. During the final application of the acoustic environment, the strain and acceleration response spectra showed considerable variation from those spectra obtained prior to impact of the panel. Fatigue damage from acoustic loading however, was limited to partial de-bonding around the edges of the composite panel.

Cockburn, J. A.↗
NASA NTRS1975

N-body simulations of disks

The methods used in large-scale n-body simulations are discussed. However, the present review concentrates on results already obtained in n-body simulations using systems containing up to 200,000 simulation stars. Results are presented which show that the stability criterion developed for flattened systems applies only to truly axisymmetric instabilities. Purely stellar disks acquire rather large velocity dispersions, generally two or more times the velocity dispersion required by Toomre (1964) for axisymmetric stability. In computer simulations, the bar-forming instability can be prevented only by comparatively large velocity dispersions. However, simulations including the effects of the galactic halo and core as a fixed background field show that bar formation can be prevented for fixed halo components as large or larger than the self-consistent disk component. Experiments performed to determine the collisional relaxation time for large-scale gravitational n-body calculations show that these models are indeed 'collisionless'.

Hohl, F.↗
NASA NTRS1977

Liquid-circulating garment controls thermal balance

Experimental data and mathematical model of human thermoregulatory system have been used to investigate use of liquid-circulatory garment (LCG) to control thermal balance. Model proved useful as accurate simulator of such variables as sweat rate, skin temperature, core temperature, and radiative, evaporative, and LCG heat loss.

Kuznetz, L. H.↗
NASA NTRS1977

Generalized environmental control and life support system computer program (G189A) configuration control, phase 2

The G189A simulation of the Shuttle Orbiter ECLSS was upgraded. All simulation library versions and simulation models were converted from the EXEC2 to the EXEC8 computer system and a new program, G189PL, was added to the combination master program library. The program permits the post-plotting of up to 100 frames of plot data over any time interval of a G189 simulation run. The overlay structure of the G189A simulations were restructured for the purpose of conserving computer core requirements and minimizing run time requirements.

Mcenulty, R. E.↗
NASA NTRS1992

Solar models with enhanced energy transport in the core

The discrepancy between the observed and computed flux of solar neutrinos can be eliminated if the temperature gradient in the core of the model, and hence its central temperature, is reduced. This could in principle be accomplished by reducing the core opacity; alternatively, it has been proposed that a fraction of the energy transport in the solar core results from the motion of hypothetical particles (the so-called WIMPs). The resulting changes in solar structure have measurable effects on the frequencies of the solar 5-minute oscillations. Here, models are considered where the opacity in the solar core has been decreased in a manner which roughly simulates the effects of the WIMPs. The analysis of the models and their frequencies provides insight into the consequences of modifications to the physics of solar models. In particular, it is found that models with the observed neutrino flux are inconsistent with observations of low-degree solar oscillations.

Christensen-Dalsgaard, J.↗