Search NASA⌕ Search

SEARCH · Search NASA

Results for “time-dependent data”

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

Ultra-sensitive radon assay using an electrostatic chamber in a recirculating system

Rare event searches such as neutrinoless double beta decay and Weakly Interacting Massive Particle detection require ultra-low background detectors. Radon contamination is a significant challenge for these experiments, which employ highly sensitive radon assay techniques to identify and select low-emission materials. This work presents the development of ultra-sensitive electrostatic chamber (ESC) instruments designed to measure radon emanation in a recirculating gas loop, for future lower background experiments. Unlike traditional methods that separate emanation and detection steps, this system allows continuous radon transport and detection. This is made possible with a custom-built recirculation pump. A Python-based analysis framework, PyDAn, was developed to process and fit time-dependent radon decay data. Radon emanation rates are given for various materials measured with this instrument. A radon source of known activity provides an absolute calibration, enabling statistically-limited minimal detectable activities of 20 µBq. These devices are powerful tools for screening materials in the development of low-background particle physics experiments.

47 OTHER INSTRUMENTATION↗

Integration of Online Cross-Section Generation Capability with Depletion and Transient Solvers in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE)-based reactor multiphysics analysis application jointly developed by Argonne and Idaho National Laboratories under the DOENE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In FY25, an online crosssection generation capability based on the Self-Shielding Application Programming Interface (SSAPI) was demonstrated for TRISO-fueled reactor problems under steady-state conditions. This fiscal year, that capability was extended to support depletion and transient multiphysics calculations, enabling high-fidelity analyses that generate self-shielded cross sections on the fly from the actual evolving composition and temperature states rather than from pre-tabulated libraries. For depletion, a two-way coupling was established in which SSAPI computes compact-averaged self-shielded cross sections that the depletion solver then uses to advance the Bateman equations, with the updated compositions returned to SSAPI at each step; the depletion module was refactored to support both library-based and SSAPI-based cross sections, and additional logic was added to track daughter isotopes and to exclude minor isotopes for efficiency. For transient analysis, the SSAPI multigroup library was extended with the kinetics data required for time-dependent calculations, the Improved Quasi-Static (IQS) scheme was coupled with SSAPI, and several supporting capabilities were implemented, including a self-shielding treatment that lets control rods and drums move within a self-shielded model, which had previously been impossible and had ruled out rod- and drum-movement transients with on-the-fly cross sections altogether, a new mixing scheme for delayed-neutron precursor decay constants, a checkpoint-based restart workflow, and performance improvements such as pointwise cross-section interpolation and the bypassing of unnecessary Dancoff factor calculations. The implemented capabilities were verified against Serpent Monte Carlo solutions. For depletion, a prismatic pin-cell problem based on a Next Generation Nuclear Plant (NGNP) Very High Temperature Reactor benchmark showed excellent agreement, with eigenvalue differences within 200 pcm over the entire burnup range (up to 140 MWD/kgU) and fission-product and actinide inventories agreeing to within 0.8% and 2.5%, respectively; a heat-pipe microreactor assembly problem with a much higher fuel loading confirmed the same behavior and quantified the bias introduced when the multigroup equivalence effect is neglected. For transient analysis, a pin-cell problem with a step reactivity insertion and temperature feedback reproduced the analytically expected asymptotic power and showed close agreement between the direct and IQS solutions, and a two-dimensional microreactor core problem with control-drum rotation exercised the new moving-drum self-shielding treatment and demonstrated successful coupling of the online crosssection generation with both the direct and IQS transient methods. The capability was further exercised on a full-core pebble-bed problem, in which Griffin was coupled with the System Analysis Module (SAM) to simulate load-following operation of the gPBR with the Doppler feedback resolved at the TRISO fuel kernel temperature. These developments in Griffin provide a convenient, high-fidelity approach to cross-section generation for advanced thermal reactors with geometrically complex and highly heterogeneous configurations, including TRISO-fueled prismatic and pebble-bed systems, and support steady-state, depletion, and transient multiphysics calculations. They also enable self-shielded cross sections to be evaluated directly at the actual coupled state of the system, thereby establishing a foundation for high-fidelity, fully coupled multiphysics analysis of advanced reactors

Park, H.↗

The Role of Nuclear Data Sensitivities in Prompt α-Eigenvalue Predictions of Delayed Critical Benchmarks

Alpha (α) eigenvalues, which describe the logarithmic time derivative of the neutron population in a multiplying system, are integral to time-dependent behavior and diagnostic applications. However, uncertainties in the evaluated nuclear data can significantly impact the accuracy of transport simulations for such quantities. This work explores the use of machine learning models to predict two key outputs, α-eigenvalues and keff bias, using input features derived from α-eigenvalue sensitivities to nuclear data. The criticality safety benchmark models used in this study come from the International Handbook of Evaluated Criticality Safety Benchmark Experiments. Three models, random forest, XGBoost, and NGBoost, are trained on both energy-resolved and energy-summed α sensitivities. For the α-eigenvalue bias prediction, NGBoost achieved the highest R 2 (0.9476) using energy-resolved features, while XGBoost performed best using summed sensitivities. In contrast, when predicting the keff bias, all the models showed moderate predictive capability (best R 2 ≈ 0.72), as the mapping from the static α-sensitivities to the static keff bias was less direct. SHAP (SHapley Additive exPlanations) analysis was used to interpret the model predictions. Across both prediction tasks, the features associated with neutron capture [H-1 (n, γ)], uranium scattering reactions (such as 235 U elastic/inelastic), and actinide capture/fission reactions (such as 239 Pu and 234 U) were consistently identified as the most impactful. This highlights the key role of specific nuclear reactions and energy ranges in shaping both time-dependent and steady-state criticality behavior. These results demonstrated that α-sensitivities, despite being computed for time-dependent metrics, can provide valuable insights for predicting both α-eigenvalues and the keff bias. Moreover, machine learning models offer a promising pathway for uncovering important nuclear data dependencies and guiding future data evaluation efforts.

Nuclear data↗

Electrochemical Modeling of PID Leakage Current of PV Modules: Steady-State Current, Transient Current, and RC-Equivalent Circuit

Potential-induced degradation (PID) remains a significant reliability concern for photovoltaic (PV) modules, arising when a voltage difference between the module frame and the solar cells drives unintended leakage current through the glass-encapsulant stack. Although PID ultimately manifests as PID-s, PID-p, or PID-c, the underlying behavior of the leakage current-its magnitude and time dependence-requires clearer electrochemical interpretation. Traditional explanations attribute the initial transient current to bulk capacitive elements of the glass, encapsulant, and antireflection coatings, and the steady-state current according to their effective ohmic resistance. More recent studies, however, indicate that electrochemical charge-transfer processes at the encapsulant-metallization interface can play a dominant role in defining the leakage-current path. This paper develops a unified electrochemical framework for modeling PID leakage current. First, an RC-equivalent circuit is formulated by combining conventional RC elements with a Randles-type interface to capture transient leakage current through double-layer capacitance and faradaic processes at ionic-electronic boundaries. Second, the steady-state current-voltage behavior is explained using a linearized Butler-Volmer relationship, showing that the measured ohmic response corresponds to the low-overpotential limit of charge-transfer kinetics. Analytical results demonstrate that, for typical module materials-3.2-mm soda-lime glass and 0.45-mm encapsulant-the dominant modulators to PID leakage current are the glass surface resistance (under dry-surface conditions), the glass bulk capacitance, and the encapsulant resistance (under wet-surface conditions), with soda lime glass surface and EVA/POE encapsulant resistances primarily governing steady-state current. The proposed electrochemical model is validated against measured leakage-current data, showing good agreement in both the magnitude and the time-dependent evolution of PID leakage current.

14 SOLAR ENERGY↗

Autoregressive long-horizon prediction of plasma edge dynamics *

Accurate modeling of scrape-off layer (SOL) and divertor-edge dynamics is vital for designing plasma-facing components in fusion devices. High-fidelity edge fluid/neutral codes such as SOLPS-ITER capture SOL physics with high accuracy, but their computational cost limits broad parameter scans and long transient studies. We present transformer-based, autoregressive surrogates for efficient prediction of 2D, time-dependent plasma edge state fields. Trained on SOLPS-ITER spatiotemporal data for the KSTAR tokamak, the surrogates forecast electron temperature, electron density, and radiated power over extended horizons. We evaluate model variants trained with increasing autoregressive horizons (1–100 steps) on short- and long-horizon prediction tasks. Longer-horizon training systematically improves rollout stability and mitigates error accumulation, enabling stable predictions over hundreds to thousands of steps and reproducing key dynamical features such as the motion of high-radiation regions. Measured end-to-end wall-clock times show the surrogate is orders of magnitude faster than SOLPS-ITER, enabling rapid parameter exploration. Prediction accuracy degrades when the surrogate enters physical regimes not represented in the training dataset, motivating future work on data enrichment and physics-informed constraints. Overall, this approach provides a fast, accurate surrogate for computationally intensive plasma edge simulations, supporting rapid scenario exploration, control-oriented studies, and progress toward real-time applications in fusion devices.

autoregressive deep learning↗

Systems Engineering and Analysis in Support of a US Federal Staging Facility for UNF

The US Department of Energy Office of Nuclear Energy (DOE-NE) Office of Spent Fuel and High-Level Waste Disposition is examining a set of system options and conducting supporting analyses to inform the development of an integrated waste management system, which may include one or more federal staging facilities (FSFs) for used nuclear fuel (UNF ) sited using a collaborative siting process. This paper focuses on the ongoing activities in two systems engineering and analysis work areas: (1) data and tools development, validation, and maintenance and (2) systems engineering execution. Within the first work area, the STANDARDS 5.0 UNF data and analysis tool, formerly known as UNF-ST&DARDS, is being developed as a foundational resource to assist in the management of UNF data. It has the key capability to model UNF throughout the entire back end of the fuel cycle. STANDARDS also includes several compatible analysis tools for the time-dependent characterization of UNF and related systems by interfacing with the SCALE code system for nuclear analysis and COBRA-SFS for thermal analysis. Also, within the data and tools area is the Next Generation System Analysis Model (NGSAM), which is an agent-based simulation software tool expressly designed to be capable of modeling the waste management system, including the transportation of UNF to and from a FSF. NGSAM has been developed to enable informed decision-making by providing the capability to analyze various potential system options for the management of UNF and high-level radioactive waste. Finally, in the systems engineering execution area, the team has begun to apply a disciplined systems engineering approach at the system level along with supporting analysis to guide the development of the FSF project requirements (including associated transportation infrastructure). Systems engineering principles and practices and their adaptation/application to design and development activities will ensure that the waste management system is effectively implemented as work proceeds. Other activities include investigating the implications of changes in various assumptions and parameters related to waste management systems, such as UNF acceptance rates, receipt logic, facility capacities and capabilities, use of standardized canisters, and different assumed facility operation start dates. Keywords: federal staging facility (FSF), used nuclear fuel (UNF), integrated waste management (IWM) system, Next Generation System Analysis Model (NGSAM), STANDARDS, systems engineering

Joseph, Robert↗

Synthetic modeling of soft x-ray emissivity for magnetic island analysis in LTX-β

We present a synthetic soft x-ray (SXR) forward-modeling framework that characterizes the emissivity structure of rotating magnetic islands in the Lithium Tokamak eXperiment-β across space and time. Magnetic islands associated with tearing modes produce modulations in line-integrated SXR brightness. Traditional tomographic methods struggle to resolve these structures in devices with limited sightlines, resulting in an under-determined inversion problem, or otherwise necessitate equilibrium reconstruction data, preventing their use in active control. Our approach removes this challenge by combining a fully three-dimensional ray-tracing model with a time-dependent emissivity prescription of a $m/n$ = 2/1 magnetic island. The model is validated through helical island geometry, rotation measured by magnetic diagnostics, and equilibrium constraints from PSI-Tri reconstructions. Synthetic brightness signals are generated for each photodiode sightline from a modeled emissivity profile and directly compared with experimental data from a tangential midplane SXR array. By fitting the synthetic diagnostic output to the observed brightness across time, we simultaneously infer key geometric parameters-such as island radial location, island width, and rotation frequency-without requiring full tomographic reconstruction, all with an average absolute deviation of less than 3%. This work demonstrates that forward modeling with a single tangential array can extract key island parameters in a spherical tokamak, providing a computationally efficient alternative to conventional SXR tomography and providing a pathway toward real-time magnetic-island characterization in future devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Total absorption spectroscopy of two isomers in 70 Cu influencing nucleosynthesis signatures

Isomers have long been known to be important for astrophysical nucleosynthesis processes, yet they are often neglected in network calculations due to computational limitations or lack of data. "Astromers" are astrophysically metastable nuclear states that can greatly impact nucleosynthesis pathways. In this work we show that astromers further impact the time-dependent electromagnetic signal during heavy element nucleosynthesis. In an experiment performed at the National Superconducting Cyclotron Laboratory, three 𝛽-decaying states of 70 Cu (6 − ground state, and two isomeric states: 3 − and 1 + ) were produced. 𝛽-feeding values were extracted from experimental spectra and compared to shell-model and QRPA+PVC calculations. Here, average 𝛾-ray energies from the 𝛽-decay events were incorporated into simulations of heavy element nucleosynthesis and were found to exhibit different energy release profiles over time, which may impact, in aggregate, time-dependent observable signals.

59 ≤ A ≤ 89↗

IER-620 CED-3b: Experiment Execution Summary for the Pulsed-Neutron Die-Away Experiments (PNDA) with Propylene Glycol and Mobilmet 423

There is a strong need for new benchmarks to validate neutron thermal scattering laws (TSLs). Lawrence Livermore National Laboratory (LLNL) has designed a Pulsed-Neutron Die Away (PNDA) testbed for this purpose. The experiment has a deuterium-tritium (D-T) neutron generator that impinges a 10 -4 s, mono-energetic pulse of 14.1 MeV neutrons on a target sample. After the pulse, the neutron population moderates and establishes a thermal equilibrium within the sample, with a fundamental spatial mode and characteristic decay-time eigenvalue, ⍺. The ⍺ eigenvalue can be extracted from the experimental measurements of the time-dependent neutron flux coming off the surface of the sample and can then be used as an integral parameter (similar to k eff ) to validate nuclear data involved with neutron migration, thermalization, and absorption. For moderating materials and geometric configurations, the ⍺ eigenvalue is heavily dependent on thermal neutron scattering of the target material. For that reason, a PNDA experiment can have a higher sensitivity to TSLs than is commonly available with the k eff parameter in critical experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Comprehensive analytical model of the dynamic 𝑍 pinch

In this work we present an analytical 1D axisymmetric model describing the evolution of the dynamic 𝑍 pinch. This model is capable of predicting the trajectories of the imploding sheath's magnetic piston and preceding shock front, along with the velocity, pressure, density, and magnetic field profiles, for any time-dependent current, spatially varying initial density profile, and weak initial axial field. The implosion is divided into stages, with each stage described by a set of coupled ordinary differential equations derived from the ideal MHD equations. Comparisons with experimental data from the COBRA pulsed-power facility are quite promising and imply this model could prove useful in designing and analyzing future pulsed-power experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Constraining the Galactic bar using the M92 stellar stream

Stellar streams are excellent probes of the gravitational potential in which they evolve. In the Milky Way (MW), globular cluster (GC) streams are routinely used to infer properties about time-dependent perturbations of the underlying potential. This implies that streams with Galactocentric radii small enough to be perturbed by the MW bar should offer constraints on it, such as its pattern speed, which currently has a wide range of values reported in the literature and is important when studying stellar kinematics. The GC M92 has a small pericentre and should be affected by the bar. It has a diffuse stellar stream, but confirming stream members has previously been hindered by a lack of spectroscopic data. In this paper, we use Dark Energy Spectroscopic Instrument (DESI) observations together with photometric and astrometric data to obtain spectroscopic members of the M92 stream for the first time. We identify a clear spatial distribution and gradients in distance moduli, proper motions, and radial velocities that confirm the stream's existence. We compare the observed stream to mock streams generated in different barred potentials and estimate the MW bar's pattern speed $Ω= 29.1^{+0.7}_{-0.4}$ km s$^{-1}$ kpc$^{-1}$ and $\dot Ω= 0.7^{+3.5}_{-2.3}$ km s$^{-1}$ kpc$^{-1}$ Gyr$^{-1}$. This is the first time a stellar stream is used to probabilistically infer these bar properties, and it opens up an exciting realm of inner Galactic potential characterisation using stellar streams.

Byström, Amanda [Edinburgh U., Inst. Astron.] (ORC↗

Part-scale evolution of fine-scale microstructural heterogeneity in solid-state additive manufacturing

Current solid-state additive manufacturing methods, refined through costly and time-consuming trial and error, have spurred interest in computational models that replicate material behavior under typical thermomechanical conditions (e.g., strain-rate ~ 102 s−1). These models, however, struggle to capture time-dependent microstructural evolution. In this work, Additive Friction-Stir Deposition (AFSD) is used as a representative case study for part-scale quantification of microstructure evolution at a fine spatial resolution (200 μm) by examining a liquid-nitrogen-cooled stop-action build via energy-dispersive X-ray diffraction coupled with a multi-channel detector. These results inform modeling efforts by linking process asymmetry to stored plastic strain, residual elastic strain, and texture development, and unlike current state-of-the-art characterization methods (e.g., EBSD or neutron diffraction), this approach provides both the spatial resolution and collection efficiency necessary to quantify fine-scale microstructural heterogeneity over large component volumes. As such, this technique provides essential validation data for computational models, e.g., crystal plasticity, enabling future prediction of heterogeneous behavior in AFSD and other additive manufacturing processes.

Franz, Cole [ORNL] (ORCID:0000000213465881)↗

FAST-1.2.2: A Computer Code for Thermal-Mechanical Nuclear Fuel Analysis under Steady-state and Transients

Fuel Analysis under Steady-state and Transients (FAST) is the U.S. Nuclear Regulatory Commission (NRC)’s computer code that calculates the steady-state and transient response of nuclear reactor fuel rods during long-term in-reactor burnup, anticipated operational occurrences (AOOs), design basis accidents (DBAs), and dry storage conditions. The code calculates the temperature, pressure, and deformation of a fuel rod as functions of time-dependent fuel rod power and coolant boundary conditions. The phenomena modeled by the code include heat conduction through the fuel and other materials, heat transfer from the cladding-to-coolant, cladding elastic and plastic deformation (including creep), fuel-cladding mechanical interaction, fission gas release from the fuel, rod internal pressure, void volume, and cladding oxidation. The code contains necessary material and coolant properties, as well as clad-to-coolant heat transfer correlations, for normal operation through postulated accidents and AOOs for today’s U.S.-based light water reactor (LWR) fuel designs. FAST-1.2.2 also contains preliminary materials and models for new LWR fuel concepts, such as accident tolerant fuel (ATF), and non-LWR fuel concepts such as metallic fuels for sodium fast reactors (SFRs). FAST has been developed for use on Windows and Linux operating systems. This document describes FAST-1.2.2 and is one of a series of documents on the code; the other documents detail the material properties used by FAST as well as its integral assessment to experiments and commercial data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Parallel-in-Time Solution of Allen-Cahn Equations by Integrating Operator Learning into the Parareal Method

While recent advances in deep learning have shown promising efficiency gains in solving time-dependent partial differential equations (PDEs), matching the accuracy of conventional numerical solvers still remains a challenge. One strategy to improve the accuracy of deep learning-based solutions for time-dependent PDEs is to use the learned model as the coarse propagator in the Parareal method and a traditional numerical method as the fine solver. However, successful integration of deep learning into the Parareal method requires consistency between the coarse and fine solvers, particularly for PDEs exhibiting rapid changes such as sharp transitions. Here, to ensure this consistency, we propose using convolutional neural networks (CNNs) to learn the fully discrete time-stepping operator defined by the same numerical scheme employed as the fine solver. We demonstrate the effectiveness of the proposed method in solving the classical and mass-conservative Allen–Cahn (AC) equations. Through iterative updates in the Parareal algorithm, our approach achieves a significant computational speedup compared to traditional fine solvers while converging to high-accuracy solutions. Our results highlight that the proposed hybrid Parareal algorithm effectively accelerates simulations, particularly when implemented on multiple GPUs, and converges to the desired accuracy in only a few iterations. Another advantage of our method is that the CNN model is trained on trajectory-based data generated from random initial conditions, such that the trained model can be used to solve the AC equations with various initial conditions without retraining. This work demonstrates the potential of integrating neural network methods into parallel-in-time frameworks for efficient and accurate simulations of time-dependent PDEs.

97 MATHEMATICS AND COMPUTING↗

Egg yolk as a model for gelation: From rheometry to flow physics

Egg yolks are an excellent model for studying sol-gel transitions, particularly the power law viscoelasticity that defines the critical point of gelation. However, prior studies lack comprehensive datasets and fail to visualize flow behavior linked to temperature and time-dependent linear and nonlinear rheology. Here, we present a detailed dataset characterizing egg yolk viscoelasticity across temperature, time, and forcing amplitude using oscillatory shear, step strain, step stress, and constant high strain rate. Novel protorheology visualizations link rheological properties with observable flow behavior. Our findings highlight the nuanced determination of the critical gel point, emphasizing observation timescale dependencies. We compare methods to identify critical temperatures for gelation, including power law viscoelasticity, moduli crossover, diverging zero-shear viscosity, and emerging equilibrium elastic modulus, while visualizing flow consequences near these transitions. Egg yolk is an accessible, realistic, and nontoxic material relevant to the physicist and the chef alike, making it ideal for understanding the rheology of critical gels. By integrating protorheology photos and videos with rigorous rheometric data, we deepen the understanding of critical gels, with broader impacts for studying other materials with sol-gel transitions.

Marsh, Maxwell C. [Department of Mechanical Scienc↗

DEPRECATED AI-Batt-OS (Autonomous Identification of Battery Life Models - Open Source) [SWR 21-17]

DEPRECATED. This repository was archived by the owner on Jun 30, 2026. It is now read-only. Open source implementation of some of the methods utilized by AI-Batt, a battery lifetime modeling and analysis toolkit provided by the National Laboratory of the Rockies (NLR). This software demonstrates the use of bi-level optimization and symbolic regression techniques to semi-autonomously identify algebraic models predicting the capacity fade of lithium-ion batteries during calendar aging. Modeling the degradation of batteries is a complex task, due to the difficulty in separating the time-dependent and time-independent factors impacting cell level degradation, across multiple data series with different numbers of measurements and/or data quality. Bi-level optimization enables model parameters to be optimized to either the entire data set or to individual data series, allowing statistical disambiguation of global behaviors (data series independent) and local behaviors (data series dependent). Symbolic regression is used to automatically search for optimal low-dimesional models predicting the variation of locally optimized parameters versus time-independent experimental variables from millions of possible models, resulting in a more accurate and repeatable model identification process than is possible by a manual search. The provided tools also implement cross-validation and bootstrap resampling schemes, empowering statistical model comparison/selection and quantification of model uncertainties. An example script replicates the results from the manuscript "Challenging Practices of Algebraic Battery Life Models through Statistical Validation and Model Identification via Machine-Learning", submitted to ECS. All code is written in MATLAB. Requires the Statistics and Machine Learning Toolbox. Contact Dr. Paul Gasper at Paul.Gasper@nlr.gov for any questions.

Gasper, Paul↗

Pulsed-Neutron Experiments at the Inherently Safe Subcritical Assembly

The pulsed neutron technique is a powerful, dynamic method to assay the reactivity of a multiplying system. This work presents the novel application of the pulsed neutron technique to the Inherently Safe Subcritical Assembly, an experimental configuration accepted by the International Criticality Safety Benchmark Evaluation Project Handbook. The experiments were replicated with COG11.3, a continuous-energy Monte Carlo code. The pulsed neutron data were analyzed using the Sjöstrand and Gozani area-ratio methods and by extracting the prompt neutron decay constant. Subsequent static k-eigenvalue and 𝛼-eigenvalue simulations were also performed for the same configurations. Neutron detector dead-time effects from the experiments were corrected using the Backwards Extrapolation Method and shown to have a negligible impact on the estimated reactivities. The results highlight that capturing time-dependent effects like delayed neutron precursor buildup are essential to accurately reproduce experimental results. They also highlight the importance of shielding the detectors from generator source neutrons in deeply subcritical configurations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

High-Fidelity Modeling Approach for Coupling Thermal, Fluid, and Neutronic Analysis in an NTP Engine

Nuclear thermal propulsion (NTP), with its increased efficiency compared to traditional rocket engines, is a potentially enabling technology for space missions of higher complexity, such as delivering heavy machinery or crew to the moon or other planets. Thus, exploring modeling approaches that can effectively couple the thermal, fluid, and neutronic behaviors of NTP reactors is essential to mitigate cost and schedule challenges with hardware development. This paper details one such approach which utilizes Cardinal to unify reactor physics analysis performed using OpenMC, heat transfer analysis performed using the MOOSE framework, and thermal hydraulic analysis performed using the MOOSE Thermal Hydraulic Module. This approach is demonstrated with a proof of concept that focuses on a legacy reactor geometry from the historic NERVA program. This preliminary model found an overprediction in, when compared to the reference data, axial power peaking and the temperatures of both solids and fluids by roughly 20%. It is clear, however, that the model conserved a number of broad trends, and a few model refinements should reduce the observed error. This approach has definite merit and should enable the production of informative results for more modern NTP systems. This work is a first step in creating a transient reactor model that can be updated across a sequence of time steps based on continually fluctuating, time-dependent boundary conditions.

Nuclear Thermal Propulsion↗