Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Advanced reactor modeling and 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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Bayesian Calibration of Nuclear Graphite Property Models Accounting for Model Inadequacy and Impacts on Component Performance

Nuclear-grade structural graphite is extensively utilized in the core designs of various advanced nuclear reactors. In the reactor environment, graphite is subjected to prolonged exposure to extreme conditions, including high temperatures, radiation, and potentially molten salt and oxygen. Such exposure can induce several degradation mechanisms in graphite, including nonuniform volumetric strains caused by irradiation and thermal expansion, leading to stresses that may compromise the performance of graphite components. Assessing component integrity requires accurate models of graphite's thermomechanical response. This report documents the Bayesian calibration of thermomechanical properties for nuclear-grade graphite and their application to graphite component modeling and simulation using the Grizzly code. As part of this work, uncertainty-quantified models were developed for the elastic modulus, coefficient of thermal expansion, irradiation-induced dimensional change, and irradiation-induced creep for graphite grades IG-110, NBG-18, NBG-17, PCEA, and 2114. Using a hierarchical Bayesian approach, multiple experimental data sources were combined to develop Gaussian process models for the properties. Using the Kennedy O'Hagan framework, the uncertainties due to inadequacies in the model and the inherent spread in the experimental data were quantified for three different models. These uncertainty-quantified models, with a model-form correction, were subsequently applied to a coupled-physics simulation of representative graphite components, revealing that the uncertainties have a large impact on the components' deformation.

36 - MATERIALS SCIENCE↗
DOE OSTI2022

Accelerated statistical failure analysis of multifidelity TRISO fuel models

Statistical nuclear fuel failure analysis is critical for the design and development of advanced reactor technologies. Although Monte Carlo Sampling (MCS) is a standard method of statistical failure analysis for fuels, the low failure probabilities of some advanced fuel forms and the correspondingly large number of required model evaluations limit its application to low-fidelity (e.g., 1-D) fuel models. In this paper, we present four other statistical methods for fuel failure analysis in Bison, considering tri-structural isotropic (TRISO)-coated particle fuel as a case study. The statistical methods considered are Latin hypercube sampling (LHS), adaptive importance sampling (AIS), subset simulation (SS), and the Weibull theory. Using these methods, we analyzed both 1-D and 2-D representations of TRISO models to compute failure probabilities and the distributions of fuel properties that result in failures. The results of these methods compare well across all TRISO models considered. Overall, SS and the Weibull theory were deemed the most efficient, and can be applied to both 1-D and 2-D TRISO models to compute failure probabilities. Moreover, since SS also characterizes the distribution of parameters that cause TRISO failures, and can consider failure modes not described by the Weibull criterion, it may be preferred over the other methods. Finally, a discussion on the efficacy of different statistical methods of assessing nuclear fuel safety is provided.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Thermal Fluid Modeling Approaches in SAM for High Temperature Gas-Cooled Reactor Applications

This paper provides an overview of the unique modeling and simulation challenges and needs in the system and safety analysis of two common types of High Temperature Gas-cooled Reactor (HTGR) designs. The challenges are associated with the complex core geometry configurations and the change of the dominant heat removal mechanisms between normal operating conditions and decay heat removal transients. This requires that simulation tools utilize models with reasonable length scales and capture both the localized heat transfer and the core-scale heat transfer at the same time. To meet such analysis needs, different methodologies were developed. In conclusion, this includes adopting existing methods and further developing/improving them, as well as proposing innovative methods that leverage advanced computational frameworks such as the Multiphysics Object Oriented Simulation Environment (MOOSE).

HTGR↗
DOE OSTI2021

The advanced tokamak path to a compact net electric fusion pilot plant

Abstract Physics-based simulations project a compact net electric fusion pilot plant with a nuclear testing mission is possible at modest scale based on the advanced tokamak concept, and identify key parameters for its optimization. These utilize a new integrated 1.5D core-edge approach for whole device modeling to predict performance by self-consistently applying transport, pedestal and current drive models to converge fully non-inductive stationary solutions, predicting profiles and energy confinement for a given density. This physics-based approach leads to new insights and understanding of reactor optimization. In particular, the levering role of high plasma density is identified, which raises fusion performance and self-driven ‘bootstrap currents’, to reduce current drive demands and enable high pressure with net electricity at a compact scale. Solutions at 6–7 T, ∼4 m radius and 200 MW net electricity are identified with margins and trade-offs possible between parameters. Current drive comes from neutral beam and ultra-high harmonic (helicon) fast wave, though other advanced approaches are not ruled out. The resulting low recirculating power in a double null configuration leads to a divertor heat flux challenge that is comparable to ITER, though reactor solutions may require more dissipation. Strong H-mode access (x2 margin over L–H transition scalings) and ITER-like heat fluxes are maintained with ∼20%–60% core radiation, though effects on confinement need further analysis. Neutron wall loadings appear tolerable. The approach would benefit from high temperature superconductors, as higher fields would increase performance margins while potential for demountability may facilitate nuclear testing. However, solutions are possible with conventional superconductors. An advanced load sharing and reactive bucking approach in the device centerpost region provides improved mechanical stress handling. The prospect of an affordable test device which could close the loop on net-electric production and conduct essential nuclear materials and breeding research is compelling, motivating research to validate the techniques and models employed here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

On FIRE mode in KSTAR

We report on the status of Fast Ion Regulated Enhancement (FIRE) mode experiments in the Korea Superconducting Tokamak Advanced Research. This regime is being developed for high-performance, steady-state operation which features a stationary ion internal transport barrier, enabling a central ion temperature approaching 10 keV to be sustained for up to 50 s, without the need for delicate profile control and with no significant impurity accumulation. As its key novelty lies in the significant contribution of fast ions that stabilize core turbulence, the regime has been named FIRE mode. To achieve this regime, neutral beam injection is applied at moderate power levels near the L–H power threshold, while maintaining low plasma density to avoid the L–H transition. The scenario is typically established in diverted magnetic configurations. The core features of FIRE mode were investigated through power balance analysis and fluctuation measurements, revealing a clear transport bifurcation in the ion channel. At the plasma edge, FIRE mode occasionally exhibits I-mode characteristics, particularly in unfavorable magnetic null configurations with q 95 ∼ 4, including the presence of weakly coherent modes. In terms of MHD activity, sawtooth oscillations are observed but appear to be stabilized during the high-performance phase. Fast ion-driven Alfvénic eigenmodes (AEs), indicated by strong frequency chirping near 200 kHz, are also observed. Additionally, lower-frequency MHD activities, distinct from the fast ion-driven AEs, are present and have some influence on plasma performance. The enhancement of core confinement is primarily attributed to fast ion effects, which were evaluated with respect to dilution, alpha stabilization, and resonant interactions with turbulence using gyrokinetic analyses. Among these, the dilution effect was found to be the most dominant. The characteristics of FIRE mode were compared with those of other hot ion plasma scenarios, such as supershot and hot ion mode. While they share many similarities, FIRE mode is distinguished by the accessibility to conditions with T i ≈ T e , presence of I-mode edge features, and its long-duration sustainment. Predictive simulations of FIRE mode were performed using integrated transport modeling with TRIASSIC, employing the TGLF anomalous transport model. These simulations confirmed the critical role of fast ions in achieving this regime. The future prospect of FIRE mode for application in fusion reactors is discussed, with an emphasis on possibly extending the regime to higher density operation.

FIRE mode↗
DOE OSTI2023

Characterization of the acoustic cavitation in ionic liquids in a horn-type ultrasound reactor

Most ultrasound-based processes root in empirical approaches. Because nearly all advances have been conducted in aqueous systems, there exists a paucity of information on sonoprocessing in other solvents, particularly ionic liquids (ILs). In this work, we modelled an ultrasonic horn-type sonoreactor and investigated the effects of ultrasound power, sonotrode immersion depth, and solvent’s thermodynamic properties on acoustic cavitation in nine imidazolium-based and three pyrrolidinium-based ILs. The model accounts for bubbles, acoustic impedance mismatch at interfaces, and treats the ILs as incompressible, Newtonian, and saturated with argon. Following a statistical analysis of the simulation results, we determined that viscosity and ultrasound input power are the most significant variables affecting the intensity of the acoustic pressure field (P), the volume of cavitation zones (V), and the magnitude of the maximum acoustic streaming surface velocity (u). V and u increase with the increase of ultrasound input power and the decrease in viscosity, whereas the magnitude of negative P decreases as ultrasound power and viscosity increase. Probe immersion depth positively correlates with V, but its impact on P and u is insignificant. 1-alkyl-3-methylimidazolium-based ILs yielded the largest V and the fastest acoustic jets – 0.77cm 3 and 24.4ms -1 for 1-ethyl-3-methylimidazolium chloride at 60W. 1-methyl-3-(3-sulfopropyl)-imidazolium-based ILs generated the smallest V and lowest u – 0.17cm 3 and 1.7ms -1 for 1-methyl-3-(3-sulfopropyl)-imidazolium p-toluene sulfonate at 20W. Sonochemiluminescence experiments validated the model.

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

High-Fidelity Energy Deposition Ignition Model Coupled with Flame Propagation Models at Engine-like Flow Conditions

With the heightened pressure on car manufacturers to increase the efficiency and reduce the carbon emissions of their fleets, more challenging engine operation has become a viable option. Highly dilute, boosted, and stratified charge, among others, promise engine efficiency gains and emissions reductions. At such demanding engine conditions, the spark-ignition process is a key factor for the flame initiation propagation and the combustion event. From a computational standpoint, there exist multiple spark-ignition models that perform well under conventional conditions but are not truly predictive under strenuous engine operation modes, where the underlying physics needs to be expanded. In this paper, a hybrid Lagrangian-Eulerian spark-ignition (LESI) model is coupled with different turbulence models, grid sizes, and combustion models. The ignition model, previously developed, relies on coupling Eulerian energy deposition with a Lagrangian particle evolution of the spark channel, at every time-step. The spark channel is attached to the electrodes and allowed to elongate at a speed derived from the flow velocity. The LESI model is used to simulate spark ignition in a nonquiescent crossflow environment at engine-like conditions, using converge commercial computational fluid dynamics (CFD) solver. The results highlight the consistency, robustness, and versatility of the model in a range of engine-like setups, from typical with Reynolds-averaged Navier-Stokes (RANS) and a larger grid size to high fidelity with large-eddy simulation (LES) and a finer grid size. The flame kernel growth is then evaluated against Schlieren images from an optical constant volume ignition chamber with a focus on the performance of flame propagation models, such as G-equation and thickened flame model, versus the baseline well-stirred reactor model. Finally, future development details are discussed.

Advanced ignition modeling↗
DOE OSTI2020

Multiscale-Informed Modeling of High Temperature Component Response with Uncertainty Quantification

This report summarizes a joint effort between Argonne National Laboratory, Idaho National Laboratory, and Los Alamos National Laboratory to develop and deploy constitutive models targeted at predicting the life of Grade 91 alloy components subjected to high temperature environments typical of those that structural components in advanced nuclear reactors would experience. Two distinct, but complementary constitutive modeling approaches have been taken here. The first employs a phenomenological viscoplastic model for which parameters have been calibrated based on experimental data for a wide range of Grade 91 alloy that has undergone a variety of processing. A Bayesian approach was used to derive distributions of uncertain parameters for this model based on this data set. The second approach is a reduced order model suitable for engineering-scale analysis that is based on the results of a large set of mesoscale simulations. Mesoscale models allow for the microstructure and composition of a particular alloy to be directly taken into account in the computation of the viscoplastic response, but are computationally expensive, which makes it impractical to directly call those models for the material constitutive response in an engineering-scale simulation. The reduced-order representation of the response of the underlying model used here allows for an engineering-scale model to take into account the characteristics of the underlying microstructure, while only incurring a reasonable computational expense. Both of these approaches have different strengths, and are applicable for different parts of the design/analysis process. The phenomenological models can be readily parameterized based on a set of experimental data for a given class of materials and used for scoping calculations. Once a specific material is chosen and adequately characterized, the reduced order models can accurately predict the response of that specific alloy, and because the models are based on predictive models of the underlying microstructure, they can be used to more confidently predict the response under conditions in regions where there is limited experimental data. Both of these models have been integrated in the Grizzly code, which is used here to perform proof-of-concept uncertainty quantification analyses of a simple component under prototypical conditions. The built- in stochastic analysis capabilities in the MOOSE framework that Grizzly is built on are used here to run large sets of simulations for this uncertainty quantification analysis. As would be expected, because the reduced order models are developed for a much more tightly defined alloy, they predict tighter distributions of the time to failure than the phenomenological models, which are calibrated to a broader set of data. Also important is that these simulations demonstrate that a reduced order modeling approach can be successfully deployed to propagate uncertainties from the material scale to practical engineering-scale component simulations.

42 ENGINEERING↗
DOE OSTI2024

Utilizing Physics-Informed Synthetic Data to Train a Digital Twin for Predicting Reactor Operations

Understanding techniques to strengthen the nuclear nonproliferation regime is crucial in reducing the creation of nuclear weaponry on the basis of advancements in the nuclear energy industry. Prior to construction of a nuclear power plant, it is necessary to understand the proliferation potential of the plant’s reactor. Digital twins serve as a unique solution to recognizing reactor behavior indicative of nuclear proliferation. The following research conducted serves as a validation to training a digital twin on synthetic data fabricated via means of reactor physics simulations based on parameters of Idaho State University’s AGN-201 reactor. The synthetic data is utilized to train a long short-term memory (LSTM) recurrent neural network model. The accuracy of the predicted data is measured against real operational data to verify the reliability of the synthetic data creation methods and if these methods should be used in the future to inform inspectors of a reactor’s proliferation capabilities.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2025

Ceramography and Thermal Modelling of Irradiated TRISO Particles using BISON

This study presents a comprehensive analysis of irradiated TRISO fuel particles from the Advanced Gas Reactor (AGR) Fuel Development and Qualification Program conducted at the Advanced Test Reactor (ATR) as part of my internship. Failure Analysis using high-resolution microscopy was performed to characterize mechanisms in TRISO particle layers, including the kernel, buffer, inner pyrolytic carbon (IPyC), silicon carbide (SiC), and outer pyrolytic carbon (OPyC) coatings. Failure modes were systematically classified. Thermal conductivity properties of each TRISO layer were evaluated through experimental measurements and computational models for use in BISON finite element code. Experimental thermal conductivity values ranged from 0.5 W/m·K for the buffer layer to 168 W/m·K for SiC, while BISON simulations provided averaged values of 0.57 W/m·K (buffer), 4.0 W/m·K (PyC layers), and 15.95 W/m·K (SiC). Temperature-dependent correlations from PARFUME were implemented for kernel thermal properties.

BISON↗
DOE OSTI2021

Validation of SAS4A/SASSYS-1 for predicting steady-state single-phase natural circulation

The validation of system analysis codes for nuclear reactor systems is required for the development and application of these computational tools. Designed as a comprehensive system analysis code for advanced nuclear reactors, SAS4A/SASSYS-1 requires validation of its physics model for capturing single-phase natural circulation behavior. To support the validation of SAS4A/SASSYS-1, high-precision experiments are performed capturing steady-state single-phase natural circulation on a scaled facility with comprehensive instrumentation. In this work, dedicated tests are performed quantifying the critical modeling parameters, and a single-phase natural circulation benchmark dataset is obtained with well-documented uncertainty and comprehensive facility description. The validation is then performed against the dataset examining the capability of SAS4A/SASSYS-1 in simulating steady-state single-phase natural circulation. The experimental facility is modeled in the candidate code. Solution verification is performed using Richardson-extrapolation-based estimators which quantify and restrict numerical errors from discretization. Input uncertainty provided by the benchmark dataset is forward propagated through the candidate code, quantifying the output uncertainty in a Monte Carlo approach. The composition of the output uncertainty is also quantified through a variance-based sensitivity analysis. With the uncertainty quantified for each individual condition, a detailed comparison between the simulation results and experimental data is performed covering the whole dataset. The results show consistent agreement for all primary parameters. The current validation activity provides a valuable benchmark dataset for the validation of system analysis codes in capturing single-phase natural circulation and demonstrates satisfactory prediction capability of SAS4A/SASSYS1 for steady-state single-phase natural circulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

A quinary WTaCrVHf nanocrystalline refractory high-entropy alloy withholding extreme irradiation environments

In the quest of new materials that can withstand severe irradiation and mechanical extremes for advanced applications (e.g. fission & fusion reactors, space applications, etc.), design, prediction and control of advanced materials beyond current material designs become paramount. Here, through a combined experimental and simulation methodology, we design a nanocrystalline refractory high entropy alloy (RHEA) system. Compositions assessed under extreme environments and in situ electron-microscopy reveal both high thermal stability and radiation resistance. We observe grain refinement under heavy ion irradiation and resistance to dual-beam irradiation and helium implantation in the form of low defect generation and evolution, as well as no detectable grain growth. The experimental and modeling results—showing a good agreement—can be applied to design and rapidly assess other alloys subjected to extreme environmental conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Advances in prediction of tokamak experiments with theory-based models

The successful validation of theory-based models of transport, MHD stability, heating and current drive, with tokamak measurements over the last 20 years, has laid the foundation for a new era where these models can be routinely used in a "predict first" approach to design and predict the outcomes of experiments on tokamaks today. The capability to predict the plasma confinement and core profiles with a quantified uncertainty, based on a multi-machine, international, database of experience, will provide confidence that a proposed discharge will remain within the operational limits of the tokamak. Developing this predictive capability for the first generation of burning plasma devices, beginning with ITER, and progressing to tokamak demonstration reactors, is a critical mission of fusion energy research. Major advances have been made implementing this predict first methodology on today's tokamaks. In this work, an overview of several of these recent advances will be presented, providing the integrated modeling foundations of the experimental successes. The first steps to include boundary plasmas, and tokamak control systems, have been made. A commitment to predicting experiments as part of the planning process is needed in order to collect predictive accuracy data and evolve the models and software into a robust whole discharge pulse design simulator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A Metallic Fuel Performance Benchmark Problem Based on the IFR-1 Experiment

Metallic nuclear fuels are an active area of research and development for use in advanced reactors. Robust and accurate metallic fuel performance models are necessary for the design, safety analysis, and licensing of such reactors. However, metallic fuel performance models are not as mature as ceramic fuel performance models and require additional development. In order to better gauge the accuracy of existing models, identify high-priority models for development, and quantify the improvements made by future model development, ORNL has developed a benchmark problem based on the IFR-1 experiment. This work gathers all relevant information on the IFR-1 experiment and uses it to develop the benchmark problem. It simulates the problem using the fuel performance code BISON and compare the results to post-irradiation examination data from the experiment. Then a sensitivity analysis is performed on the BISON model to determine the benchmark problem’s sensitivity to uncertainty in the input parameters. The results suggest that BISON’s mechanical models require additional development. Plastic deformation is under-predicted in the cladding and axial swelling is over-predicted in the fuel. These problems may be related. Additionally, there is a bias in the temperature which may be a result of uncertainty in the input parameters rather than an issue with the fuel performance models.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

A Data-driven approach to Core Power distribution reconstruction in a Nuclear Reactor

This report presents the initial development of a data-driven approach for reconstructing the core power distribution in a nuclear reactor (power shape synthesis) using ex-core sensors. Traditional techniques rely on deploying a large number of detectors throughout the reactor core. However, this approach is not feasible for innovative reactor concepts like Advanced Reactors and Microreactors. First, the tight lattice pitch, designed to maximize power density, limits the space available for sensors. Secondly, the harsh operating conditions are not compatible with commercially available detectors. The method proposed in this work integrates high-fidelity modeling with data-driven techniques to accurately reconstruct power distribution across various reactor types, thereby reducing the reliance on in-core sensors. Purdue University Reactor One (PUR-1) was selected as the test case. The CAD model representing the latest configuration of the PUR-1 core was imported into the OpenMC simulation framework, and the model was built. Additionally, the previously developed MCNP6 model was updated. The two models were assessed against the data collected during an experimental campaign conducted in July 2024. Thirty gold foils were placed in three Irradiation Assemblies in PUR-1 core. Using the measured activity of the irradiated foils, the neutron flux at different core locations was estimated.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Numerical Simulation of a Natural Convection–Driven Air-Cooled Reactor Cavity Cooling System Experiment

Ensuring the efficient removal of decay heat from the reactor vessel is essential for the safety of advanced reactor technologies. Several Generation-IV concepts incorporate variations in the reactor vessel cooling systems to achieve this objective. High-temperature gas-cooled reactors utilize a reactor cavity cooling system (RCCS), a passive ex-vessel system designed to operate without active components or external power during accident conditions. The RCCS removes decay heat primarily through radiative and convective heat transfer mechanisms. Here, this study presents a comprehensive validation of a computational fluid dynamics Reynolds-averaged Navier-Stokes model for the University of Wisconsin-Madison air-cooled RCCS facility. Validation was conducted for both high- and low-power natural convection cases under a uniform heating profile. Near-wall resolution was found to be critical for accurately modeling natural convection in the RCCS; employing an all-𝑦 + wall treatment resulted in wall temperature discrepancies exceeding 50 °⁢𝐶 compared to a wall-resolved mesh. Thermal-hydraulic behaviors under natural and forced convection conditions were compared within the heated cavity and RCCS. A turbulence model sensitivity analysis indicated that low-Reynolds number k-ɛ, k-ω shear stress transport (SST), and Reynolds stress transport models produce similar wall temperature predictions. A buoyancy modeling sensitivity study revealed that the Boussinesq approximation significantly underpredicted thermal-hydraulic behavior in the RCCS. Based on these findings, modeling recommendations are provided. The validated data set along with identified sensitivities refine the modeling of natural convection in the RCCS. The information produced by this study supports RCCS design, optimization, and safety evaluations, enabling the calibration and verification of reduced-order thermal-hydraulic models.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Performance Demonstration of Self-Powered Neutron Detectors for Steady-State Reactor Operations

The irradiation testing of sensors in reactors is a crucial step towards calibrating and qualifying sensors prior to their deployment in experiments. This report details the process toward qualifying and calibrating custom-designed rhodium-based self-power neutron detectors (Rh-SPNDs) for steady-state reactor irradiations. This process serves to both demonstrate the performance capabilities of Rh-SPNDs as well as to provide experimental data for development of a sensor sensitivity model. Two designs of Rh-SPNDs were tested in various reactors to demonstrate: detection resolution in a low neutron flux environment, a delayed-response compensation technique, output linearity in a large range of neutron flux, and measurement accuracy verified with dosimetry. The detection resolution and compensation technique was demonstrated in the AGN-201m reactor at Idaho State University. The irradiation confirmed the sensors’ capability to perform steady-state operations in a low neutron flux of ~2E8 n/cm 2 -sec. Sensor output linearity coupled with the delayed-response compensation was investigated at the neutron radiography reactor at Idaho National Laboratory. A Rh-SPND was irradiated to neutron fluxes ranging from 2E8 to 2E13 n/cm2-sec range. The measured data demonstrated a wide and linear range of operation with a measured linear sensitivity of 1.0129 ×10 -13 A/W with a correlation-squared value of r 2 =0.9927. The measurement accuracy was investigated at the Advanced Test Reactor Critical reactor. The SPNDs were inserted into a test vehicle with collocated flux wires. Two irradiations with different flux levels were performed, and the SPNDs relative measurement between the two irradiations was calculated to be 1.2613 ± 0.0153 for the small SPND design and 1.1809 ± 0.0108 for the large SPND design. Both SPND measurements fell between the co-axial dosimetry result, which reported 1.218 ± 0.047. Additionally, the preliminary MCNP model for calculating SPND sensitivity was developed in parallel to this work. Modeled neutron spectrum with measured magnitude was used for inputs to determine the simulated SPND output. The results showed an overestimation of signal strength by a factor of 5, which was expected because of model simplification. This leads to future modeling work to account for signal losses from additional physical properties, including high temperature environments for FY-21.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Modeling a Thermal-Spectrum LEU-fueled Molten Salt Reactor Co-fueled with Thorium in SCALE 6.3.1 and Serpent-2

Interest in the development of advanced nuclear fission reactors for commercial electricity production has risen in recent years, with Generation IV reactor designs offering numerous advantages in safety, efficiency, fuel cycle sustainability, and waste management. In particular, thorium-fueled molten salt reactors (MSRs) are considerably promising for enhancing fuel cycle sustainability in their ability to breed fissile 233U fuel from thorium, a presently untapped and widely abundant resource. A novel MSR fuel cycle concept, the ”Sourdough” refueling and waste management strategy, has previously been demonstrated with a traditional uranium-based fuel cycle in a thermal-spectrum MSR operating with low enriched uranium (LEU) fuel with favorable neutronic performance. However, the ability to use this unique fuel cycle approach with thorium-based molten salt fuels has not yet been studied. In this work, the Sourdough fuel cycle was implemented in a small, thermal-spectrum MSR fueled with high assay low enriched uranium (HALEU) and fertile 232Th for breeding 233U fuel. Relevant neutronic data, including fuel and isothermal temperature feedback behavior, was studied using the SCALE 6.3.1 and Serpent-2 code systems, and keff data was measured during simulated depletion at 400 MWth. The Sourdough fuel cycle concept is shown to perform favorably with a thorium-fueled MSR model, thus warranting further study into its use in other MSR designs.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗