Search NASASearch

SEARCH · Search NASA

Results for “Criticality Calculations”

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

Calculating Critical Inertia of a Power System

The increasing integration of renewable energy sources in modern power systems has led to a decline in system inertia, raising concerns about frequency stability following large disturbances. Determining critical inertia is essential to prevent excessive frequency decline and ensure grid stability. This paper evaluates four different methods for calculating critical inertia, using the Electric Reliability Council of Texas (ERCOT) system as a test case. The results highlight the importance of employing multiple methodologies to capture the full spectrum of inertia requirements.

Hakim sneha, Fariha [University of Tennessee, Knox

Updating Critical Temperature Calculations for Several Secondary High Explosives

In designing experiments where high explosives (HEs) are heated, it is important to have an understanding of where thermal runaway may occur. This determination is often done by using the Frank-Kamenetskii (FK) equation. For several highly-studied HEs, the parameters necessary for FK calculations are typically referenced from decades-old literature, and more recent experimental data have shown that these values require adjustment. For example, some of the legacy parameters dangerously overpredict the critical temperature by tens of degrees relative to more recent observed values. Herein, we first summarize historical results and highlight the importance of insulation effects when estimating critical temperatures for small samples, which is especially relevant to small-scale aging, compatibility, and characterization experiments. Here, we then present updated FK parameters for several common secondary CHNO HEs like RDX, HMX, PETN, and TATB, as well as the first reported values for 2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105). Our updated parameters produce critical-temperature curves that are consistent with both mm-scale differentical scanning calorimetry (DSC) measurements and cm-scale data from the Lawrence Livermore National Laboratory one-dimensional time-to-explosion (ODTX) experiment. In analyzing the critical temperature versus sample size curves derived from the transcendental FK equation, we discovered that these curves are remarkably well described by a simple power-law function with a universal power of 0.142 and a HE-specific scaling factor.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF

SCALE Non-LWR Models for NRC Volume 5

This dataset contains input and result files of computational simulations with the SCALE code system. The simulations cover radionuclide inventory generation, criticality calculations, and dose rate/shielding analyses of various advanced reactors. Users wanting to reproduce results from this dataset are required to obtain a license to the SCALE code system for which details on the distribution can be found here: https://www.ornl.gov/scale/releases

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS

Modeling a generic TRISO-fueled heat pipe microreactor using SCALE: Depletion, transportation criticality, and shielding

This paper demonstrates the applicability of the SCALE code system to tristructural-isotropic (TRISO)-fueled heat pipe microreactors through depletion, transportation criticality, and shielding analyses of a generic reference design. The study conducted supports US Nuclear Regulatory Commission code readiness efforts for advanced non–light-water reactor technologies and is intended as a code capability demonstration rather than as an optimization of a specific microreactor design. The modeled reactor employs high-assay low-enriched uranium (HALEU) uranium oxycarbide (UCO) TRISO fuel and beryllium oxide (BeO) reflectors and operates at 7.5 MWth with a nominal lifetime of about 3 effective full power years. Representative cases for fresh and irradiated cores were selected to exercise SCALE methods relevant to reactor operation and post-irradiation transport. The discharged-core decay heat is approximately 6% of operating power immediately after shutdown. Transportation criticality calculations show that internal water ingress is the dominant reactivity effect, with fully flooded fresh core and irradiated core configurations remain above the subcriticality criterion, even with the available control mechanisms. Shielding calculations for a simplified transportation package indicate that normal-condition dose rates are governed mainly by shielding thickness and cooling time, whereas the breached hypothetical accident case is governed primarily by cooling time. Overall, the study shows that SCALE supports depletion, transportation criticality, and shielding evaluations efficiently for TRISO-fueled heat pipe microreactors within a single code system.

Criticality

Multiphysics Analyses of the Protected and Unprotected Loss of Forced Cooling Accidents in the HTR-PM

Here, we present the multiphysics simulation results for the protected pressurized and depressurized loss of forced cooling (PLOFC and DLOFC) events in the High-Temperature gas-cooled Reactor--Pebble-bed Module (HTR-PM) equilibrium core using the Griffin-Pronghorn coupled code system. Additionally, this paper discusses the strategy for estimating the spontaneous fission neutron source needed for unprotected events and re-criticality calculations. The solutions of the protected PLOFC and DLOFC events were verified against similar solutions obtained for temperature evolutions from the open literature. Both the average and maximum pebble surface temperatures behaved as expected during the DLOFC and remained below 1800 K. The PLOFC results are highly dependent on the ability to resolve the natural circulation in the core, which is impacted by the mesh resolution in Pronghorn. Furthermore, we present the results of the unprotected DLOFC transient to predict the timing of the re-criticality event, which occurred 47 hours after the onset of the transient, and the new steady-state power of 1.3 MW.

42 - ENGINEERING

Developing a Nuclear Quality Assurance Compliant Design Methodology for Neutronic Analysis of Xe-100 Design

The primary objective of this work is to develop a design methodology compliant with nuclear quality assurance standards for the Xe-100 neutronic design verification studies. To achieve this, a Monte Carlo model of the Xe-100 reactor was constructed using the exclusion principle, transformation technique, and universe-based level specification following Idaho National Laboratory (INL) NQA level-1 compliant standards and an NQA-1 compliant version of MCNP6. The model encompasses the entire reactor core structures, including the upper plenum, core region, and lower plenum sections, along with all sub-components. The active core section was represented using the spectral regions, each comprising a particular fuel composition and temperature averaged over the considered zone, calculated by X-energy using Very Superior Old Programs (VSOP). Additionally, a component-wise temperature map was implemented into the model, not only for the core region but also for the structural components. Temperature-dependent cross-section libraries, along with thermal scattering law libraries, generated using INL NQA-1 compliant version of NJOY21, were utilized for each isotope in the burnt fuel and the structural materials. Furthermore, the volume of each modeled component was estimated using a stochastic approach with the ray tracing method in MCNP and criticality calculations were performed.

22 GENERAL STUDIES OF NUCLEAR REACTORS

Strengthening mechanisms for microstructures containing unimodal and bimodal γ' precipitates in ATI 718Plus

Here, the influence of γ' precipitate size distribution on the deformation mechanisms under tensile loading in ATI 718Plus was studied. A set of aging treatments within the temperature range of 720 °C–900 °C was performed on solution-treated samples to obtain various γ' precipitate size distributions. Unimodal and bimodal γ' precipitate size distributions were achieved through single-step and two-step aging sequences, respectively, and such microstructures were tensile tested to failure to assess their yield strength, ultimate tensile strength, and elongation-to-failure. Some of the tensile samples were interrupted after achieving 3–4 % plastic strain, and the deformed microstructures were examined using transmission electron microscopy to investigate the γ' precipitate-dislocation interactions. For the unimodal γ' precipitate size distribution samples with the smaller γ' precipitates (radius ~ 7 nm), dislocations sheared through the precipitates. Both dislocation loops and paired dislocations were observed for the microstructures containing larger γ' precipitates (radius ~ 24 nm). The microstructure containing a bimodal γ' precipitate size distribution, which included average γ' precipitate radii of ~6 nm and ~28 nm, exhibited shearing as the dominant deformation mechanism, and this microstructure exhibited the highest strength values. The experimental observations were rationalized based on the theoretically-calculated critical resolved shear stress values for shearing and looping and a modified model for predicting the yield strength for bimodal microstructures was introduced.

36 MATERIALS SCIENCE

An Analytic Benchmark for Neutron Boltzmann Transport with Downscattering—Part IV: PFNS and $\bar{ν}$ Uncertainty Propagation

An analytic benchmark with continuous-energy cross sections was previously derived to validate criticality calculations. Here, to extend the utility of the analytic benchmark to verify the implementation of $\bar{ν}$ and prompt fission neutron spectrum (PFNS) uncertainty propagation methods, new simplified forms that are dependent on the incident (fission-causing) neutron energy, as well as the outgoing neutron energy for the PFNS, are introduced in this work. The analytical forms for the flux and adjoint flux are derived for the extended benchmark and used to determine the 𝑘-eigenvalue sensitivity to $\bar{ν}$ and PFNS. The 𝑘-eigenvalue uncertainty due to $\bar{ν}$ and PFNS is calculated for the analytic benchmark using simplified$\bar{ν}$ and PFNS representations based on the ENDF-B/VIII.0 239 Pu evaluation. Because of the low sensitivity of the analytic benchmark to the physical PFNS, a nonphysical high-sensitivity PFNS is also presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS

COG Software Architecture Design Description Document

This COG Software Architecture Design Description Document describes the organization and functionality of the COG Multiparticle Monte Carlo Transport Code for radiation shielding and criticality calculations, at a level of detail suitable for guiding a new code developer in the maintenance and enhancement of COG. The intended audience also includes managers and scientists and engineers who wish to have a general knowledge of how the code works. This Document is not intended for end-users.

61 RADIATION PROTECTION AND DOSIMETRY

Preliminary Investigation of Target keff Search for Fuel Feed in MSRs

Fluid fuel molten salt reactor (MSR) concepts may achieve continuous-at-power operation with continuous or batch fresh fuel feed to counter reactivity loss from fuel burnup and fission product buildup. Accurate modeling and simulation of MSRs requires consideration of online fuel feed during depletion as well as a capability to determine the required feed amount to maintain a target multiplication factor (keff). As a first step supporting a full implementation of a new search capability into the SCALE code system’s TRITON reactor physics and depletion sequence, a preliminary investigation was performed to identify the best approach for a critical feed search capability considering computation time and accuracy. Using a simple iterative approach in which initial feed amounts were guessed based on previous depletion steps, the main drivers of computation time were found as the Monte Carlo neutron transport settings and the applied tolerance for the target keff. In the applied model, with tight convergence criteria (keff convergence to 8 pcm and a tolerance of ±15 pcm), the critical search calculation required twice the number of neutron transport calculations compared to a depletion calculation without critical search. In contrast, with a relaxed tolerance of ±50 pcm, the overall computation time for the critical search for merely increased by 30-40%. Further studies will be performed before implementing the critical search capability into TRITON.

Faulkner, Jonathon [ORNL]

COG11.3 Abstract

COG is a high-resolution code for the Monte Carlo simulation of coupled particle transport in arbitrary3-D geometry. COG will transport neutrons, protons, deuterons, alpha particles with energies up to hundreds of GeV, and photons with energy ranges limited by the available cross section sets and physics models. Electrons can be transported via the EGS5 electron transport kernel, electrons can also be transported. The COG code is a significant upgrade from earlier Monte Carlo transport codes and has been written specifically to makes it more versatile, accurate, and easy to use. COG has provisions for calculating deep penetration (shielding) problems, criticality problems, and neutron activation problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS

Availability of Critical Benchmark Experiments for the Pebble Tanker Transportation Model for Nuclear Criticality Safety Validation of TRISO Pebbles

This study addresses the need for comprehensive investigations into TRi-structural ISOtropic (TRISO) fuel pebble transportation validation. In this work, an exploratory model, the pebble tanker(PT), was developed with the aim of facilitating the validation of nuclear criticality safety calculations in the context of industrial-scale transportation of TRISO fuel. The PT model was designed to investigate the availability and applicability of critical benchmark experiments crucial for assessing the transportation of these pebbles. This work incorporated sensitivity/uncertainty (S/U) similarity studies to quantify the applicability of critical benchmark experiments and to address nuclear data uncertainties in the context of TRISO transportation. Two container models were investigated: one for the Hermes-type pebble and one for the Pebble Bed Modular Reactor (PBMR)–type pebble. The models were simplified, considering fuel, containment, and either water or air, to enable a focus on the underlying physics of applications involving TRISO fuel pebbles using the PT model. A crucial aspect under consideration was the capacity of the transport package to hold pebbles while ensuring subcriticality in the flooded state. An approach in the criticality validation process involves assessing the similarity between systems through an integral index parameter evaluation. This involves calculating a correlation coefficient (referred to as c k ) based on shared nuclear data–induced uncertainty between a benchmark experiment and the application of the PT model. To facilitate this analysis, the SCALE tools, particularly the CSAS6-Shift, TSUNAMI-3D-Shift, and TSUNAMI-IP sequences, were employed for comprehensive studies in neutronics and S/U analysis. Our findings showed that there are sufficient critical experimental benchmarks to perform this validation of the PT model in the most reactive state, i.e. when the tanker is flooded. This paper provides valuable insights into validating a transport package for Generation IV TRISO fuel pebbles.

22 GENERAL STUDIES OF NUCLEAR REACTORS

A Python Tool for Aqueous Plutonium Nitrate Density Law Input Preprocessing in MCNP6

Here, this work develops a predictive density tool in Python, named Plutonium Nitrate Solutions (PuNS), to reduce bias and uncertainty in nuclear criticality safety calculations for plutonium nitrate systems. The Pitzer method and an empirical method were implemented into the PuNS tool to generate atom densities for use in MCNP6 material cards. These material cards are directly prepared into an MCNP6 input text file and are calculated based on customizable user inputs of plutonium content, nitric acid content, temperature, and plutonium isotope weight percentages. The PuNS tool is validated and verified against the International Criticality Safety Benchmark Evaluation Project Handbook experiments and is observed to predict densities within a root mean square error of 0.89% for the Pitzer method and 1.82% for the empirical method. These errors in density lead to up to 1569 pcm difference in MCNP6 calculated k eff for the Pitzer method and up to a 1751 pcm difference for the empirical method when compared to experimental benchmarks. Simultaneous work is also being performed at Los Alamos National Laboratory and the University of New Mexico to create a similar tool for plutonium chloride solutions, named Plutonium Chloride Solution, which aims to provide the accreditation of the chlorine absorption. These capabilities will not only provide more accurate models but also facilitate an improved understanding of solution systems and a potential relaxation in the conservatism of current aqueous plutonium processing criticality safety limits.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA

Defluorination Mechanisms and Real-Time Dynamics of Per- and Polyfluoroalkyl Substances on Electrified Surfaces

Per- and polyfluoroalkyl substances (PFAS) are persistent environmental contaminants found in groundwater sources and a wide variety of consumer products. In recent years, electrochemical approaches for the degradation of these harmful contaminants have garnered a significant amount of attention due to their efficiency and chemical-free modular nature. However, these electrochemical processes occur in open, highly non-equilibrium systems, and a detailed understanding of PFAS degradation mechanisms in these promising technologies is still in its infancy. To shed mechanistic insight into these complex processes, we present the first constant-electrode potential (CEP) quantum calculations of PFAS degradation on electrified surfaces. These advanced CEP calculations provide new mechanistic details about the intricate electronic processes that occur during PFAS degradation in the presence of an electrochemical bias, which cannot be gleaned from conventional density functional theory calculations. We complement our CEP calculations with large-scale ab initio molecular dynamics simulations in the presence of an electrochemical bias to provide time scales for PFAS degradation on electrified surfaces. Taken together, our CEP-based quantum calculations provide critical reaction mechanisms for PFAS degradation in open electrochemical systems, which can be used to prescreen candidate material surfaces and optimal electrochemical conditions for remediating PFAS and other environmental contaminants.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Criticality Safety SU based USL Calculation for UCl3-NaCl Fuel Salt Operations

The Molten Chloride Reactor Experiment (MCRE) is a fast spectrum, molten salt fueled reactor that is planned to be constructed and operated in the Laboratory for Operation and Testing in the U.S. (LOTUS) testbed, formerly known as the ZPPR cell, at Idaho National Laboratory (INL). MCRE will provide valuable data to support design, licensing, and operation of full scale commercialized molten salt reactor designs.

99 - GENERAL AND MISCELLANEOUS

Criticality Safety S/U based USL Calculation for UCl3-NaCl Fuel Salt Operations

Recent experimental data has shown inconsistencies with the 35Cl(n,p) cross-section. The cross-section uncertainty data does not account for the recent data. The 1s relative uncertainty is assumed to be 100% for the 35Cl(n,p) cross section over all energies above 0.017 MeV TerraPower and LANL have recently done cross-section measurements and developed new cross-sections for 35Cl.

99 - GENERAL AND MISCELLANEOUS

Guiding Principles for Geochemical/Thermodynamic Model Development and Validation in Nuclear Waste Disposal: A Close Examination of Recent Thermodynamic Models for H + —Nd 3+ —NO 3 - (—Oxalate) Systems

Development of a defensible source-term model (STM), usually a thermodynamical model for radionuclide solubility calculations, is critical to a performance assessment (PA) of a geologic repository for nuclear waste disposal. Such a model is generally subjected to rigorous regulatory scrutiny. In this article, we highlight key guiding principles for STM model development and validation in nuclear waste management. We illustrate these principles by closely examining three recently developed thermodynamic models with the Pitzer formulism for aqueous H + —Nd 3+ —NO 3 - (—oxalate) systems in a reverse alphabetical order of the authors: the XW model developed by Xiong and Wang, the OWC model developed by Oakes et al., and the GLC model developed by Guignot et al., among which the XW model deals with trace activity coefficients for Nd(III), while the OWC and GLC models are for concentrated Nd(NO 3 ) 3 electrolyte solutions. The principles highlighted include the following: (1) Principle 1. Validation against independent experimental data: A model should be validated against experimental data or field observations that have not been used in the original model parameterization. We tested the XW model against multiple independent experimental data sets including electromotive force (EMF), solubility, water vapor, and water activity measurements. The results show that the XW model is accurate and valid for its intended use for predicting trace activity coefficients and therefore Nd solubility in repository environments. (2) Principle 2. Testing for relevant and sensitive variables: Solution pH is such a variable for an STM and easily acquirable. All three models are checked for their ability to predict pH conditions in Nd(NO 3 ) 3 electrolyte solutions. The OWC model fails to provide a reasonable estimate for solution pH conditions, thus casting serious doubt on its validity for a source-term calculation. In contrast, both the XW and GLC models predict close-to-neutral pH values, in agreement with experimental measurements. (3) Principle 3. Honoring physical constraints: Upon close examination, it is found that the Nd(III)-NO 3 association schema in the OWC model suffers from two shortcomings. Firstly, its second stepwise stability constant for Nd(NO 3 ) 2+ (log K 2 ) is much higher than the first stepwise stability constant for NdNO 3 2+ (log K 1 ), thus violating the general rule of (log K 2 –log K 1 ) < 0, or $\frac{K1}{K2}$>1. Secondly, the OWC model predicts abnormally high activity coefficients for Nd(NO 3 ) 2 + (up to ~900) as the concentration increases. (4) Principle 4. Minimizing degrees of freedom for model fitting: The OWC model with nine fitted parameters is compared with the GLC model with five fitted parameters, as both models apply to the concentrated region for Nd(NO 3 ) 3 electrolyte solutions. The latter appears superior to the former because the latter can fit osmotic coefficient data equally well with fewer model parameters. The work presented here thus illustrates the salient points of geochemical model development, selection, and validation in nuclear waste management.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W

Automated Direct Perturbation Calculations with SCALE TSUNAMI [Abstract]

In nuclear criticality safety analysis, the sensitivity of the eigenvalue keff to uncertainties in nuclear data and its evaluation are crucial. The TSUNAMI sequences within the SCALE code system offer users various options with both multigroup (MG) and continuous-energy (CE) 3D Monte Carlo (MC) transport capabilities for calculating keff sensitivity coefficients and storing them in a sensitivity data file (SDF). Each methodology available in TSUNAMI offers distinct advantages and limitations, and its effectiveness can vary based on the specific problem being solved. As a best practice, practitioners typically use the direct perturbation (DP) method as a confirmatory step alongside their sensitivity calculations to verify the accuracy of the sensitivity data generated. In this process, DP calculations are usually performed on select nuclides, those considered most important for validating their total sensitivities. However, because of code limitations, analysts use a workaround method when conducting DP calculations for a single nuclide: rather than perturbing the nuclide's microscopic cross section, an equivalent number density for this nuclide is calculated to reflect the effect of a change in the macroscopic cross section due to a perturbation in the microscopic cross section. The current approach requires rerunning the CSAS criticality calculation several times with model changes. Although this method can yield results with acceptable accuracy, it is labor-intensive and prone to errors.

AZURE: SAMMY