Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nuclear Reactor 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 379 records · Page 21

2025 Intern Poster

The Hot Fuel Examination Facility (HFEF) at the Materials and Fuels Complex (MFC) houses the largest U.S. inert atmosphere hot cell for nuclear material research. Key features include the precision gamma scanning (PGS), Fuel Accident Condition Simulator (FACS), Neutron Radiography Reactor (NRAD), and the focus of this project, the Metallograph Loading Cell (MET Cell). The MET Cell performs tests on spent nuclear fuel, such as microhardness testing, microscopy, and neutron radiography. However, the MET Cell’s existing pressure and lighting control systems are outdated and inefficient, with inadequate documentation for system changes over time. This project aims to design a new automated control system for the MET Cell, ensuring longevity (minimum ten years), ease of troubleshooting/repair, and integration into the building monitoring system. The design process addressed challenges such as space restrictions, varied voltages within enclosures, sourcing new components, and security limitations. Compliance with NFPA 70, UL508A, MFC Physical Security, and INL Engineering standards was essential. The project involves repurposing an existing PLC to manage lighting and pressure control using digital and analog signals, simplifying wiring, and ensuring thorough documentation for future reference.

42 - ENGINEERING↗

Stereolithographic geometry model of the IBR-2M experimental facility

The IBR-2M is a fast research reactor that operates in supercritical condition for ∼ 800 μs every 200 ms. Two reflector parts in nickel rotate in opposite directions generating 1.8 GWth peak power when they align with the fuel zone changing the reactor status from deep subcritical to supercritical. The reactor core uses high-enriched plutonium fuel and is cooled by sodium. This reactor has been modeled by MCNP and SERPENT computer programs. The MCNP model uses combinatorial geometry, whereas the SERPENT model employs the Stereolithographic (STL) geometry representation that can be used by 3D printers. The STL geometry was constructed using the CUBIT computer program. The CUBIT program was also used for a three-dimensional visualization of the Monte Carlo models. SERPENT and MCNP models use the same geometry, material specifications, and nuclear data. The latter are based on the ENDF/B-7.0 library. SERPENT and MCNP using same geometry and same material specifications produce similar k{sub eff} values within 120 pcm.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Initial study on cross-section generation requirements for a PBR equilibrium core

A Serpent model of the equilibrium core HTR-PM small modular nuclear reactor in China, was developed for use in cross-section preparation studies in order to guide methods development for the Griffin reactor multiphysics application. The model includes detailed isotopics for 10 distinct pebble burnup groups in 126 core zones with unique fuel and moderator temperatures obtained from a coupled neutronics-thermal-fluids equilibrium core calculation using Griffin-Pronghorn. A sensitivity study of the fuel and moderator temperatures for various core regions was performed with the MOOSE stochastic tools. The results show that the uncertainties are, not unexpectedly, dominated by the value of the fluid temperature and that the power level, heat transfer coefficient and effective conduction to neighboring pebbles and fluid constitute, at best, second order effects. The temperature uncertainty range varies from 28 K to 57 K between the core entry and exit planes, respectively, but these values are probably higher. We still have to quantify the significance of these uncertainties in the preparation of cross-sections in future work. In addition, we verified that the effective pebble approximation used in the PEBBED and V.S.O.P. computer codes works well for the preparation of region averaged cross-sections. Nevertheless, there are some discrepancies in the cross-sections when compared to the multi-pebble model, which could affect the prediction of peak values and the depletion calculation. We conclude that is highly desirable for future studies with Griffin to be able to handle both the 'effective' pebble approximation and the multi-pebble approach for various pebble burnup groups. This enables Griffin users with the flexibility to perform higher-fidelity studies. Finally, we initiated the preparation of cross-sections for various core regions from the full core Serpent reference model. We quantified the differences in 26 group cross-sections from infinite domain models versus the full core approach. These reference cross-sections will serve to verify the double heterogeneity, self-shielding, and spectrum-correction methods in Griffin. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Deployment of the HFIRCON transport and depletion tool for plutonium-238 production studies

Irradiation of {sup 237}Np-bearing targets in Oak Ridge National Laboratory's (ORNL) High Flux Isotope Reactor (HFIR) results in the efficient production of {sup 238}Pu, which, in the form of heat source PuO{sub 2}, is used as a reliable power source for deep-space and planetary NASA missions. A technology demonstration subproject was initiated at ORNL in 2011 to develop and implement the technology required to establish a {sup 238}Pu supply chain. A systematic progression of NpO{sub 2}/Al cermet (20 vol.% NpO{sub 2}) activities to date has successfully demonstrated target fabrication, irradiation, and chemical recovery processes. Recent program tasks have included the development of the HFIRCON transport and depletion tool for efficient reactor physics analyses and the evaluation of increased NpO{sub 2} loadings (i.e., beyond 20 vol.%) and NpN-based targets. This paper documents the deployment of the HFIRCON code to assess various Np concentrations in NpO{sub 2}- and NpN-based targets in HFIR's inner small vertical experiment facilities. Results indicate that {sup 238}Pu production and quality can be enhanced with increased Np loadings; however, target conversion rates are reduced. The results recorded in this paper, thermal and material balance evaluations, and testing requirement planning will be used to determine whether increased NpO{sub 2} loadings or NpN-based targets will be further considered. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Coupled Monte Carlo and thermal-hydraulics modeling of a prismatic gas reactor fuel assembly using Cardinal

Cardinal is a MOOSE application that couples OpenMC Monte Carlo transport and NekRS computational fluid dynamics to the MOOSE framework, closing the neutronics and thermal-fluid gaps in conducting tightly-coupled, high-resolution multiscale and multiphysics analyses. By leveraging MOOSE's interfaces for wrapping external codes, Cardinal overcomes many challenges encountered in earlier multiphysics coupling works, such as file-based I/O or overly-restrictive geometry mapping requirements. In this work, we leverage a subset of the multiphysics interfaces in Cardinal to perform coupling of OpenMC neutron transport, MOOSE heat conduction, and THM thermal-fluids for steady-state modeling of a prismatic gas reactor fuel assembly.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

An IIPG-based finite element framework in MOOSE for modeling fiber reinforced composite failure governed by extrinsic cohesive laws

Fiber reinforced composites (FRC) are widely applied to many important engineering structures such as silicon carbide ceramic matrix composite (SiC CMC) of nuclear fuel claddings. Many material and structure failures are initiated at the interface between two distinctive materials, such as the debonding of interfaces of matrix/fiber of composites or bulk matrix/inclusion in metals at micro-scales. Accurate and robust modeling of these interface failures is critical to the development of novel material and structure with improved functionality. This paper presents a discontinuous Galerkin (DG) finite element method based on the Incomplete Interior Penalty scheme and its implementation in MOOSE for modeling material interface behaviors. The proposed method naturally accommodates important extrinsic cohesive zone material models and is capable of modeling the failures of material interfaces across which material properties are highly distinctive. We demonstrate a good performance of the proposed method using single fiber and multiple fiber debonding problems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Development of a point kinetics subroutine for Molten Salt reactors in RELAP5-3D

RELAP5-3D has historically focused on accurate modeling of design-basis and beyond-design-basis accidents in solid fuel light-water reactors. The assumption used in solid fuel reactors that delayed neutron precursors remain at the location in which they were produced is inherently incorrect for a molten salt reactor with flowing fuel. Therefore, new methods for modeling MSR kinetics were required. Subroutine used in RELAP5-3D were updated to analyze point kinetics in a molten salt reactor during both steady state and transient analyses. The phenomena related to MSR kinetics that were added to RELAP5-3D include arrays to store previous values for independent and dependent variables, arrays to add terms to existing kinetics equations, and a new reactivity bias term. A custom MATLAB assessment code was used to establish a suite of verification tests. This capability is undergoing final documentation and testing for public release with a future version of RELAP5-3D.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Supo Thermal Model Development IV

This report describes the continuation of the Computational Fluid Dynamics (CFD) model of the Supo cooling system described in the reports, Supo Thermal Model Development, II, and III. In the previous reports, the solution natural convection of Supo was simulated using a 2D, axisymmetric model that altered the actual geometry. A 3D geometry of Supo will more accurately estimate the natural convection heat transfer coefficient (HTC) of the system since Supo contains several non-symmetric internal components that can only be modeled in 3D. This, in turn, will improve the natural convection HTC results for comparison purposes with past operational data and previous 2D CFD results. Unfortunately, the complexities of performing multiphase 3D CFD simulations has proven too large of a task for our current capabilities. Instead, this report will focus on the 3D model development for the 2.39 kW/L (31.3 kW) case and limitations and challenges associated with the 3D simulations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Performance of a thermionic converter module utilizing emitter and collector heat pipes

A thermionic converter module simulating a configuration for an out-of-core thermionic nuclear reactor was designed, fabricated, and tested. The module consists of three cylindrical thermionic converters. The tungsten emitter of the converter is heated by a tungsten, lithium heat pipe. The emitter heat pipes are immersed in a furnace, insulated by MULTI-FOIL thermal insulation, and heated by tungsten radiation filaments. The performance of each thermionic converter was characterized before assembly into the module. Dynamic voltage, current curves were taken using a 60 Hz sweep and computerized data acquisition over a range of emitter, collector, and cesium-reservoir temperatures. An output power of 215 W was observed at an emitter temperature of 1750 K and a collector temperature of 855 K for a two diode module. With a three diode module, an output power of 270 W was observed at an average emitter temperature of 1800 K and a Collector temperature of 875 K.

Kroeger, E. W.↗

Development of a Molecular Tagging Velocimetry Technique for Non-Intrusive Velocity Measurements in Low-Speed Gas Flows

N2O molecular tagging velocimetry (N2O-MTV) is developed for use in very-high-temperature reactor environments. Tests were carried out to determine the optimum excitation wavelength, tracer concentration, and timing parameters for the laser system. Using NO tracers obtained from photo-dissociation of N2O, velocity profiles are successfully obtained in air, nitrogen, and helium for a large range of parameters: temperature from 295 to 781 K, pressure from 1 to 3 bars, with a velocity precision of 0.01 m/s. Furthermore, by using two read pulses at adjustable time delays, the velocity dynamic range can be increased. An unprecedented dynamic range of 5,000 has been obtained to successfully resolve the flow during a helium blowdown from 1000 m/s down to 0.2 m/s. This technique is also applied to the high-temperature test facility (HTTF) at Oregon State University (OSU) during a depressurized condition cooldown (DCC) event. Details of these measurements are presented in a companion paper. This technique shows a strong potential for fundamental understanding of gas flows in nuclear reactors and to provide benchmark experimental data to validate numerical simulations.

Andre, M. A.↗

Coupled Reactor Multiphysics and Mass Scalability Assessment for Crewed Megawatt-Class NEP System Architectures

Nuclear Electric Propulsion (NEP) is an in-space propulsion technology capable of enabling opposition and conjunction class crewed Mars missions. NEP subsystems include the reactor for heat generation, a power conversion system (PCS), power management and distribution (PMAD), electric propulsion subsystem (EPS), and a primary heat rejection system. Specific mass, or αe (kg / kWe), is a key performance parameter (KPP) of the propulsion system which is directly scalable with the performance and mass estimates for individual components. To inform technology maturation planning, full system and component level parametric modeling is ongoing to explore the design trade space and illustrate the effect of subsystem design choices on the system KPPs. In this study, scaling of high-assay, low-enriched uranium (HALEU) reactor designs is assessed through coupled reactor physics and thermal hydraulics analyses. Scaling analyses evaluate the impact of system performance parameters (power level, interface temperatures) on mass for direct gas cooled, pumped liquid metal, and passively cooled heat pipe reactor concepts. Each concept requires specific geometries and working fluids to reach the performance goals of PCS interface conditions (temperature, pressure, flow rate) and system mass. The reactor assembly includes the active core (fuel, moderator, cladding, working fluid), axial and radial neutron reflectors, control drums, structural support / pressure vessel, and external radiation shielding. Each of these components are parametrically sized based on performance parameters for a megawatt-class power cycle. Results of this scaling analysis increase NEP propulsion system modeling fidelity and ultimately aim to support technology down-selection along with related technology development planning. The reactor and shield αe are a function of several PCS design choices, and reactor scaling with these parameters must be considered to enable an informed decision on reactor geometry and working fluid combination.

Nuclear Electric Propulsion↗

Coupled Reactor Multiphysics and Mass Scalability Assessment for Crewed Megawatt-Class NEP System Architectures

Nuclear Electric Propulsion (NEP) is an in-space propulsion technology capable of enabling opposition and conjunction class crewed Mars missions. NEP subsystems include the reactor for heat generation, a power conversion system (PCS), power management and distribution, electric propulsion system, and heat rejection system. Specific mass, or α (kg/kWe), is a key performance parameter (KPP) of the propulsion system which is directly scalable with the performance and mass predictions for each individual component. To inform technology maturation planning activities, full system and component level parametric modeling is ongoing to explore the design trade space and illustrate the effect of subsystem design choices on the system KPPs. In this study, scaling of high-assay, low-enriched uranium reactor designs is assessed through coupled reactor physics and thermal hydraulics analyses. Scaling analyses evaluate the impact of system performance parameters (power level, interface temperatures) on mass for direct gas cooled, pumped liquid metal, and passively-cooled heat pipe reactor concepts. Each concept requires specific geometries, fluids, and power conversion interface conditions (temperature, pressure, flow rate) to meet desired performance and mass. The reactor assembly includes the active core (fuel, moderator, cladding, working fluid), axial and radial neutron reflectors, control drums, structural support / pressure vessel, and external radiation shielding. Each of these components are parametrically sized based on performance parameters for a megawatt-class power cycle. Results of this scaling analysis increase NEP propulsion system modeling fidelity and ultimately aim to support concept down-selection along with related technology development planning. The reactor and shield α are a function of several PCS and heat rejection system design choices, and reactor scaling with these parameters must be considered to enable an informed decision on an optimal reactor geometry and working fluid combination.

Nuclear Electric Propulsion↗

New modelling capabilities in IDT

This work concerns the enhanced modelling capabilities of the discrete ordinates transport solver IDT. The novelties introduced allow for modelling unstructured geometries composed by a collection of X/Y segments and circles, and the use of reciprocity and conservation relations reduce the memory imprint as well as the computational cost of the method. IDT decomposes geometries in modular Cartesian patterns, which are the so-called Heterogeneous Cartesian Cells (HCCs), containing a chunk of the original unstructured geometries. Each HCC can be then discretized by superimposing a XY grid to refine locally the HCC. Unlike the most popular MOC, IDT performs the spatial sweeping by directional collision probabilities instead of trajectories. The sources and interface angular fluxes are expanded up to linear order. The accuracy of ray-tracing, the memory imprints together with the novel mesh refinement capabilities have been verified. A first set of preliminary results on PWR lattice problems will be presented in this paper.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A preconceptual design of an inverted Stable-Salt Reactor

This paper presents a preconceptual design of an inverted stable-salt reactor (SSR), with the unique feature that the static fuel salt is accommodated in one single bulk fuel container, while the coolant salt is circulating through the core via flow tubes penetrating the fuel container.The inverted SSR is expected to substantially decrease the capital cost compared to the state-of-the-art SSR such as the SSR-W designed by Moltex Energy, because one single gas vent at the top of the fuel container is required in the inverted SSR to prevent the leakage of the radioactive fission product gases instead of one gas vent in each of the hundreds of thousands of fuel tubes in SSR-W. In addition, the operation of the inverted SSR is largely simplified. However, numerous challenges associated with the design of the inverted SSR still remain to be overcome. Preliminary neutronics and thermal-hydraulics performance characteristics of the inverted SSR are assessed in this study. A viable core configuration of the inverted SSR is identified, and the superior characteristics of the inverted SSR, including a better neutron economy and a lower maximum fuel temperature, are demonstrated in this paper. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Transient MOC with frequency transform and DSA on unstructured mesh

We present an implementation of the transient method of characteristics (MOC) with isotropic time derivatives, accelerated by diffusion synthetic acceleration (DSA). The fully implicit frequency transform method is used to solve the transient problem with analytic precursor integration. The code works on meshes composed of almost any of the commonly used non-curvilinear finite element types, and can handle the deformation of geometry in time-dependent transport calculations. We present results of a continuous Fourier analysis for the transient multigroup DSA problem, and representative benchmarking results are presented for the C5G7-TD benchmark in 2D showing reasonable performance and agreement compared to other codes. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

SERPENT 2 and TRIPOLI-4 neutronics benchmarking of the Jules Horowitz reactor

Jules Horowitz Reactor (JHR) is a future center for nuclear material testing in Europe, currently under construction at CEA Cadarache in southern France. It has been modelled extensively with CEA's HORUS3D/N neutron calculation package, that includes a Monte-Carlo neutron transport code TRIPOLI-4. In order to support the future computational analyses on the reactor core physics, VTT's Monte Carlo neutronics code SERPENT 2 has been used to create an alternative model of the JHR, taking advantage of the pre-existing Constructive Solid Geometry (CSG) structures on which the TRIPOLI-4 model was built. This work describes a benchmark, in which two JHR core configurations were analyzed using both of these models. The first configuration is a fresh fuel core at the beginning of life, in which aluminum based dummy-devices are installed in the reactor test positions. The second configuration represents the end-of-cycle equilibrium core with a heterogeneous burnup distribution, in which CALIPSO devices containing cladding material samples are installed at the in-core experimental positions, while the reflector positions are occupied by DEN-REP devices containing UO{sub 2} fuel. Neutron flux spectra and heat deposition are the parameters of interest in this study, due to their high importance in creating the desired experimental irradiation environment. The impact of delayed gamma energy deposition will be studied in detail, since the two models account for it with different methods. This work aims to establish a SERPENT 2 model of JHR, understanding its differences and agreements with TRIPOLI-4, and discussing its potential future applications. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Investigation of methods for targeted search of dominant higher modes in subcritical systems

The reactor physics modeling of subcritical systems requires a substantially different approach from that of a critical reactor due to the presence of the so-called higher modes. Earlier investigations showed that the dominant higher modes lie in the inner part of the eigenvalue spectrum, making the conventional eigenvalue searches unfeasible and calling for a targeted search. This paper investigates the possibility of targeted and multitudinous eigenvalue calculation by testing the Krylov-Schur and the Dynamic Modes Decomposition (DMD) method for simple, analytically tractable problems. While the DMD method clearly showed its efficiency in a one-group homogeneous problem, modal analysis of a three-dimensional reflected reactor in two-groups approximation demonstrated that the increasing complexity poses challenges for both methods. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Neutronic analysis of a PWR-type SMR core using duplex ThO{sub 2}-UO{sub 2} in TRISO fuel particles

Currently, studies on small modular reactors (SMRs) present an important development due to the potential they represent in terms of safety, operational flexibility, economy, and non-energy applications. Furthermore, there is increasing interest in studying the use of thorium as fuel, as an effective way to solve problems such as the shortage of uranium reserves, reduction of nuclear waste and nuclear proliferation. Also, the use of thorium combined with highly enriched uranium TRISO particles has been studied, showing proper performance. In this work, the concept of ThO{sub 2}-UO{sub 2} duplex fuel is used, for the core configuration of a PWR type SMR that uses TRISO fuel, designed to achieve extended fuel cycles. Three distribution cases of ThO{sub 2} and UO{sub 2} in TRISO particles inside the fuel rods are compared. First, particles composed by a homogeneous mixture of ThO{sub 2} and UO{sub 2} are distributed inside of fuel elements. Second, the fuel zone of the fuel elements is divided into two radially, an internal one where the TRISO particles composed of ThO{sub 2} are distributed and the external one where the TRISO particles composed of UO{sub 2} are distributed. The third case is like the previous one, except that the particles that contain UO{sub 2} are distributed in the inner zone and those that contain ThO{sub 2} in the outer zone of the fuel elements. The comparison of the cases is carried out in terms of cycle main isotopes' mass transmutation, power distribution and temperature reactivity coefficients. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗