Search NASA⌕ Search

SEARCH · Search NASA

Results for “Source term calculation”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Notice of Submittal – 2023 Radionuclide Air Emissions Report for Los Alamos National Laboratory

This report describes the emissions of airborne radionuclides from operations at Los Alamos National Laboratory (LANL) for calendar year 2023 and the resulting off-site dose from these emissions. This document fulfills the requirements established by the National Emissions Standards for Hazardous Air Pollutants in 40 CFR 61, Subpart H – Emissions of Radionuclides other than Radon from Department of Energy Facilities, commonly referred to as the Radionuclide NESHAP or Rad-NESHAP. Compliance with this regulation and preparation of this document is the responsibility of LANL’s Rad NESHAP compliance program, which is part of the Environmental Protection and Compliance (EPC) Division. The information in this report is required under the Clean Air Act and is being submitted to the U.S. Environmental Protection Agency (EPA) Headquarters and EPA Region 6. The highest effective dose equivalent (EDE) to an off-site member of the public was calculated using procedures specified by the EPA and described in this report. LANL’s EDE was 0.43 for 2023. The annual limit is 10 millirem per year, established by the EPA in 40 CFR 61 Subpart H. All measured air emissions are modeled to a single location, known as the Maximally Exposed Individual (MEI). During calendar year 2023, LANL continuously monitored radionuclide emissions at 28 “major” release points, or stacks. The Laboratory estimates emissions from an additional 59 “minor” release points using radionuclide usage source terms in lieu of stack monitoring. Also, LANL uses an EPA approved network of air samplers around the Laboratory perimeter to monitor ambient airborne levels of radionuclides. To provide data for dispersion modeling and dose assessment, LANL maintains and operates several meteorological monitoring towers. From these various systems, a comprehensive evaluation is conducted to calculate the MEI dose for the Laboratory. The MEI can be any member of the public at any off-site location where there is a residence, school, business, or office. In 2023, this MEI location was a business at 129 New Mexico State Road 4 (NM-4), located in the northern end of White Rock. The primary contributors to the off-site dose at this location are the ambient air data at that location combined with the collected potential emissions from unmonitored (minor) sources. Overall, the MEI dose in 2023 is similar to that which has been observed in recent years, and it remains well below the EPA’s 10 millirem per year limit. Doses reported to the EPA for the past 10 years are shown in Table E1.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

2022 LANL Radionuclide Air Emissions Report (Rev. 2)

This report describes the emissions of airborne radionuclides from operations at Los Alamos National Laboratory (LANL) for calendar year 2022 and the resulting off-site dose from these emissions. This document fulfills the requirements established by the National Emissions Standards for Hazardous Air Pollutants in 40 CFR 61, Subpart H – Emissions of Radionuclides other than Radon from Department of Energy Facilities, commonly referred to as the Radionuclide NESHAP or Rad-NESHAP. Compliance with this regulation and preparation of this document is the responsibility of LANL’s Rad NESHAP compliance program, which is part of the Environmental Protection and Compliance (EPC) Division. The information in this report is required under the Clean Air Act and is being submitted to the U.S. Environmental Protection Agency (EPA) Headquarters and EPA Region 6. The highest effective dose equivalent (EDE) to an off-site member of the public was calculated using procedures specified by the EPA and described in this report. LANL’s EDE was 0.45 for 2022. The annual limit is 10 millirem per year, established by the EPA in 40 CFR 61 Subpart H. All measured air emissions are modeled to a single location, known as the Maximally Exposed Individual (MEI). During calendar year 2022, LANL continuously monitored radionuclide emissions at 27 “major” release points, or stacks. The Laboratory estimates emissions from an additional 34 “minor” release points using radionuclide usage source terms in lieu of stack monitoring. Also, LANL uses an EPA approved network of air samplers around the Laboratory perimeter to monitor ambient airborne levels of radionuclides. To provide data for dispersion modeling and dose assessment, LANL maintains and operates several meteorological monitoring towers. From these various systems, a comprehensive evaluation is conducted to calculate the MEI dose for the Laboratory. The MEI can be any member of the public at any off-site location where there is a residence, school, business, or office. In 2022, this MEI location was a business at 95 Entrada Drive, located in the eastern end of Los Alamos town site. The primary contributors to the off-site dose at this location are the ambient air data at that location combined with radioactive gas emissions from the LANSCE facility and the collected potential emissions from unmonitored (minor) sources. Overall, the MEI dose in 2022 is similar to that which has been observed in recent years, and it remains well below the EPA’s 10 millirem per year limit. Doses reported to the EPA for the past 10 years are shown in Table E1.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Depletion-driven thermochemistry of molten salt reactors: review, method, and analysis

Molten salt reactors (MSRs) are innovative advanced nuclear reactors that utilize nuclear fuel by dissolving it in a high-temperature liquid salt. This unique feature differentiates MSRs from other types of reactors and allows for enhanced safety and economic performance. The liquid fuel also entails several multiphysics effects that can complicate reactor design and operation. One primary effect termed here as depletion-driven thermochemistry is a driving force in altering the multiphysics behavior of the reactor. Essentially, depletion-driven thermochemistry is the effect that fuel depletion has on changing the chemical redox potential of the fuel salt over time. As the fuel is consumed, the redox potential shifts toward a more oxidizing state. Without active control, the changing chemistry due to depletion increases corrosion thereby limiting reactor component lifetimes. Additionally, the changing redox potential of the fuel salt alters the vapor pressures of chemical species dissolved in the fuel salt. Changing vapor pressures of species in the fuel salt is an important parameter to understand when off-gassing volatile species during normal reactor operation, and for source term characterization during accident scenario transients. The present work represents a fundamental step toward modeling and coupling the driving physics (i.e., neutronics and chemistry) involved in altering the redox potential in an MSR. Here, the neutronic code Griffin models the depletion of the fuel-salt system, while the chemical equilibrium code Thermochimica calculates the thermochemical state of the isotopic inventory, using the Molten Salt Thermodynamic Database - Thermochemical (MSTDB-TC). These two codes are tightly coupled to predict the impact of fuel depletion in altering the chemistry in MSR systems. Redox potential control methods are discussed and can be modeled using this multiphysics approach. The vapor pressures of chemical species that could be extracted to an off-gas system, as determined by the reactor’s thermochemical state, are examined. The neutronics-chemistry coupling developed in this work is expected to have potential application for analyzing corrosion, source term evolution, and material safeguards in MSR systems. Lastly, suggestions for areas of further improvements of the models to expand these capabilities by incorporating other coupled physics effects is provided.

Walker, Samuel A.↗
DOE OSTI2022

Thermodynamic assessment of lithium halide reciprocal salt systems for energy applications

Lithium halides are of interest because of their unique applications as fuel or coolants in molten salt reactors (MSRs) and electrolytes in thermally activated batteries. Because iodine is an important fission product for MSRs, the Molten Salt Thermal Properties Database-Thermochemical (MSTDB-TC) is being expanded to include relevant iodide reciprocal salt systems. In this work, we assess the thermochemical properties of lithium halide pseudobinary systems based on reported experimental phase diagrams, eutectic temperatures, and enthalpies of mixing. Within the framework of the modified quasi-chemical model in the quadruplet approximation (MQMQA), extrapolations to the related pseudoternary and pseudoquaternary systems have been performed. The pseudoquaternary representations of phase equilibria, enthalpy, entropy, and heat capacity are represented in the resultant models. Of interest for generating source terms for accident analysis is the ability to compute the vapor pressures of the dominant iodide species for the systems at any composition. Here, an example of such calculations is provided for the LiF-LiI system at x LiI = 0.01 and x LiI = 0.70, in this case indicating that the LiI vapor species are predominant.

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

Modeling Radionuclide Vaporization from Sodium Pools for SFR Mechanistic Source Term Analysis

To assist both the advanced reactor industry and U.S. Nuclear Regulatory Commission (NRC) in the pursuit of reactor design and licensing, the U.S. Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has established a mechanistic source term (MST) research project under the Multiphysics Applications technical area. As part of this effort, an MST mod/sim development pathway was developed in FY21 under the NEAMS program, which outlines the high-level objectives and near-term tasks necessary to achieve the project objectives. One such recommended task was the assessment of a previously developed thermodynamic database to model radionuclide vaporization from liquid sodium pools (specifically in the context of MST analyses), which is the focus of this topic report. As such, a thermodynamic database was created in FactSage for liquid sodium pools containing oxygen and relevant elements which represent radionuclides pertinent to MST analyses. The database was demonstrated in calculating vapor fraction curves for each element in the system at concentrations relevant to SFR MST.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

The Release of the Simplified Radionuclide Transport (SRT) Code (V.2.1)

Reactor licensing centers on the protection of the public and environment from the inadvertent release of radioactive material. Therefore, mechanistic source term analysis, or the realistic evaluation of radionuclide transport from the source to the environment for specific transient scenarios, is vital to reactor licensing efforts. Developed to resolve a gap in mechanistic source term modeling capabilities for sodium fast reactors (SFRs), the Simplified Radionuclide Transport (SRT) code created by Argonne National Laboratory (Argonne) is now utilized by advanced reactor vendors, universities, and research institutions to support a multitude of design, licensing, and research efforts. Recently, SRT version 2.1 was released, which includes improvements to code models and verification and validation (V&V) suite to support the SRT user community. The following work provides an overview of the improvements made as part of the release of SRT version 2.1. This effort is supported by the U.S. Department of Energy Office of Nuclear Energy (DOE:NE) Advanced Reactor Technologies (ART) Fast Reactor Program (FRP), as part of the program’s support of national laboratory design and safety analysis computer codes utilized by the fast reactor industry. The expansion of SRT code capabilities and improvements to code V&V associated with version 2.1 are in response to user requests and lessons learned from recent source term analyses performed by Argonne and advanced reactor vendors. They also align with the evolving role of SRT, from research and development tool to software utilized for reactor licensing calculations. The report is structured in alignment with the code improvements, as outlined in Figure 1-1. Section 2 provides background information on SRT, including its history, capabilities, and utilization. Section 3 details new code capabilities as part of version 2.1, while Section 4 focuses on the expansion of the code’s V&V suite. Lastly, Section 5 provides a summary and discussion of next steps.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

The Legendre Polynomial Axial Expansion Method

This work presents a new formulation of the axial expansion transport method explicitly using Legendre polynomials for arbitrarily high-order expansions. This new formulation also features an alternative method of axial leakage calculation to allow for nonextruded flat source region meshes. This alternative axial leakage is introduced alongside a balance equation requirement to ensure that neutron balance is preserved in the coarse mesh for a given axial leakage formulation, which allows for effective coarse mesh finite difference acceleration. A matrix exponential table method is derived to allow for fast computations of arbitrarily high-order matrix exponentials for this work and precludes the need for further research into matrix exponential calculations for this method. Numerical results are presented that demonstrate the stability of the axial expansion method in systems with voidlike regions, showcase the speedup from matrix exponential tables, and investigate the axial convergence of the method in terms of both expansion order and mesh size.

Herring, Nicholas↗
DOE OSTI2026

Experimental observation and integrated modelling of proton-beryllium fusion in He and D plasmas at JET

Validated integrated modelling of JET ITER-like wall experiments in which fusion performance is driven by reactions between fast ions and intrinsically present metal wall impurities is presented. A steady-state L-mode plasma with dominant proton-beryllium fusion and neutron yields of up to ≈ 6·10 13 s -1 is developed in He and D, via radiofrequency heating of a H minority. The fusion drive is unambiguously confirmed by the neutral particle analyser, fast ion loss detector, and γ-ray diagnostics. Experiments are analysed via an integrated modelling framework, developed to model the two-stage proton beryllium-fusion chain and produce high-fidelity fusion product source terms. The modelling chain comprises TRANSP and JETTO for plasma core modelling, LOCUST for full orbit product tracking and collisional slowing-down, DRESS to resolve two- and three-body fusion kinematics, and MCNP for neutron transport calculations. Modelling shows that the primary 9 Be(p,n) 9 B reaction is the dominant neutron emitter at naturally present concentrations of beryllium in these experiments. The yield contribution of secondary reactions between fusion products and beryllium, 9 Be(d,n) 10 B and 9 Be(α,n) 12 C, is found to be negligible. The proton-deuteron knock-on effect in D plasmas is modelled, which is calculated to contribute ≈ 25% to the total neutron yield. For both He and D discharges the total computed neutron rates match fission chamber (FC) measurements within the combined experimental and computational uncertainty, with an average discrepancy of ≈ ± 20%. Realistic proton-beryllium neutron sources are propagated through JET’s MCNP neutron transport model which shows that 235 U FCs’ response is sensitive to p–Be source changes, with up to ≈ 10% variation compared to a D–D neutron source. We show that the high-energy tail of the fast proton minority can be studied with multi-foil neutron activation. The framework is also applied to the study of interactions between fast protons and boron impurities, of relevance to ITER. We calculate that in JET conditions a significant alpha source with DT-like energies could be generated through 11 B(p, α)2α fusion, and detected via γ-emission in secondary interactions between fast alphas and boron. The work represents an important step towards validating predictive integrated modelling capabilities for non-standard fusion reactions.

JET↗
DOE OSTI2023

Meteorological, reference evapotranspiration and estimated transpiration data, July 2020-Dec 2021, Slate River Floodplain, Crested Butte, Colorado

This data package includes a time series of meteorological data (air temperature, relative humidity, barometric pressure, precipitation, wind speed, soil heat flux, incoming shortwave radiation, incoming longwave radiation, outgoing shortwave radiation and outgoing longwave radiation) from the Slate River floodplain outside Crested Butte, Colorado, a focus field site for the SLAC Floodplain Hydro-Biogeochemistry SFA. The data was generated as part of the work targeting the overarching research question for the SLAC SFA: How do ubiquitous subsurface interfaces mediate molecular-scale biogeochemical processes and groundwater quality in floodplains and watersheds? The package also includes reference evapotranspiration calculated via the Penman-Monteith equation and an estimate of willow (Salix spp.) evapotranspiration using the method of Irmak et al., J. Hydrology, 2013. Due to two malfunctioning sensors (incoming shortwave radiation from Oct 2020-May 2021 and incoming longwave radiation from May 2021-March 2022), observed net radiation and calculated evapotranspiration are only included from July-October 2020. However, users may be able to extend calculated evapotranspiration using estimates of these two parameters from other sources. Both sensors have been replaced and data beyond 2021 be uploaded to this repository when it is available.Within the data package, "FLMD.csv" describes file-level metadata and "dd.csv" defines column headers and universal terms across the dataset. The data package includes 3 "*data.csv" files, two for meteorological data for 2020 and 2021 and one for reference evapotranspiration and transpiration data from 2020 ("SLT_OBJ2_PenmanMonteith_2020_data.csv"). Each "*data.csv" file of meteorological data has a corresponding "*_InstallationMethods.csv" file that describes the location, sensor model, sensor serial number and other metadata corresponding to each measured parameter. We also include a Jupyter notebook ("PenmanMonteithET_Calculations.ipynb") that details the assumptions and equations used to calculate daily and hourly reference evapotranspiration and daily willow transpiration from raw meteorological data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The importance of delayed neutron precursors in gamma dose calculations for activated primary heat exchanger components in the Molten Salt Breeder Reactor

The Molten Salt Reactor (MSR) Multiphysics Applications technical area in the Nuclear Energy Advanced Modeling and Simulation program has supported the development of 3D Monte Carlo models of the Molten Salt Breeder Reactor (MSBR) over the last couple of years. This MSBR model was previously run with the Shift Monte Carlo code to perform radiation shielding calculations in the reactor cell area. The MSBR is a 2250 MWth (1000 MWe) liquid-fueled molten salt reactor design developed at Oak Ridge National Laboratory in the 1970s. Determining the source terms from activated primary heat exchanger (HX) components is important because delayed and prompt neutron fluxes incident on these components affect the dose rate in the primary HX maintenance areas. This information can be used in the development of remote handling procedures required during shutdown for maintenance. A methodology has been developed and is proposed in this paper to quantify the activated source term from the primary HX components as a result of the movement of the delayed neutron precursors in flowing primary fuel salt through the primary HXs in the MSBR. The goal of this research is to evaluate the gamma dose rates in the maintenance hatches above the primary HX using the activated HX source terms. The study showed that the gamma dose rates are approximately two orders of magnitudes higher when accounting for the neutron activation from the movement of delayed neutron precursors through the HXs than when flowing fuel is not considered. Thus, the movement of delayed neutron precursors must be taken into account for accurately predicting the neutron activation of primary loop components.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Development, Implementation, and Verification of Partially-Ionized Collisional Multifluid Plasma Models in Drekar

We consider the development of multifluid models for partially ionized multispecies plasmas. The models are composed of a standard set of five-moment fluid equations for each species plus a description of electromagnetics. The most general model considered utilizes a full set of fluid equations for each charge state of each atomic species, plus a set of fluid equations for electrons. The fluid equations are coupled through source terms describing electromagnetic coupling, ionization, recombination, charge exchange, and elastic scattering collisions in the low-density coronal limit. The form of each of these source terms is described in detail, and references for required rate coefficients are identified for a diverse range of atomic species. Initial efforts have been made to extend these models to incorporate some higher-density collisional effects, including ionization potential depression and three- body recombination. Some reductions of the general multifluid model are considered. First, a reduced multifluid model is derived which averages over all of the charge states (including neutrals) of each atomic species in the general multifluid model. The resulting model maintains full consistency with the general multifluid model from which it is derived by leveraging a quasi-steady-state collisional ionization equilibrium assumption to recover the ionization fractions required to make use of the general collision models. Further reductions are briefly considered to derive certain components of a single-fluid magnetohydrodynamics (MHD) model. In this case, a generalized Ohm's law is obtained, and the standard MHD resistivity is expressed in terms of the collisional models used in the general multifluid model. A number of numerical considerations required to obtain robust implementations of these multifluid models are discussed. First, an algebraic flux correction (AFC) stabilization approach for a continuous Galerkin finite element discretization of the multifluid system is described in which the characteristic speeds used in the stabilization of the fluid systems are synchronized across all species in the model. It is demonstrated that this synchronization is crucial in order to obtain a robust discretization of the multifluid system. Additionally, several different formulations are considered for describing the electromagnetics portion of the multifluid system using nodal continuous Galerkin finite element discretizations. The formulations considered include a parabolic divergence cleaning method and an implicit projection method for the traditional curl formulation of Maxwell's equations, a purely- hyperbolic potential-based formulation of Maxwell's equations, and a mixed hyperbolic-elliptic potential-based formulation of Maxwell's equations. Some advantages and disadvantages of each formulation are explored to compare solution robustness and the ease of use of each formulation. Numerical results are presented to demonstrate the accuracy and robustness of various components of our implementation. Analytic solutions for a spatially homogeneous damped plasma oscillation are derived in order to verify the implementation of the source terms for electromagnetic coupling and elastic collisions between fluid species. Ionization balance as a function of electron temperature is evaluated for several atomic species of interest by comparing to steady-state calculations using various sets of ionization and recombination rate coefficients. Several test problems in one and two spatial dimensions are used to demonstrate the accuracy and robustness of the discretization and stabilization approach for the fluid components of the multifluid system. This includes standard test problems for electrostatic and electromagnetic shock tubes in the two-fluid and ideal shock-MHD limits, a cylindrical diocotron instability, and the GEM challenge magnetic reconnection problem. A one-dimensional simplified prototype of an argon gas puff configuration as deployed on Sandia's Z-machine is used as a demonstration to exercise the full range of capabilities associated with the general multifluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Kinetic simulations of collision-less plasmas in open magnetic geometries

Laboratory plasmas in open magnetic geometries can be found in many different applications such as (a) scrape-of-layer (SOL) and divertor regions in toroidal confinement fusion devices, (b) linear divertor simulators, (c) plasma-based thrusters and (d) magnetic mirrors etc. A common feature of these plasma systems is the need to resolve, in addition to velocity space, at least one physical dimension (e.g. along flux lines) to capture the relevant physics. In general, this requires a kinetic treatment. Fully kinetic particle-in-cell (PIC) simulations can be applied but at the expense of large computational effort. A common way to resolve this is to use a hybrid approach: kinetic ions and fluid electrons. In the present work, the development of a hybrid PIC computational tool suitable for open magnetic geometries is described which includes (a) the effect of non-uniform magnetic fields, (b) finite fully-absorbing boundaries for the particles and (c) volumetric particle sources. Analytical expressions for the momentum transport in the paraxial limit are presented with their underlying assumptions and are used to validate the results from the PIC simulations. A general method is described to construct discrete particle distribution functions in a state of mirror-equilibrium. This method is used to obtain the initial state for the PIC simulation. Collisionless simulations in a mirror geometry are performed. The results show that the effect of magnetic compression is correctly described and momentum is conserved. The self-consistent electric field is calculated and is shown to modify the ion velocity distribution function in a manner consistent with analytic theory. Based on this analysis, the ion distribution function is understood in terms of a loss-cone distribution and an isotropic Maxwell-Boltzmann distribution driven by a volumetric plasma source. Finally, the inclusion of a Monte Carlo based Fokker-Planck collision operator is discussed in the context of future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Comparison of Radiation Models for a Turbulent Piloted Methane/Air Jet Flame: A Frozen-Field Study

Numerical modeling of radiative transfer in nongray reacting media is a challenging problem in computational science and engineering. The choice of radiation models is important for accurate and efficient high-fidelity combustion simulations. Different applications usually involve different degrees of complexity, so there is yet no consensus in the community. In this paper, the performance of different radiative transfer equation (RTE) solvers and spectral models for a turbulent piloted methane/air jet flame are studied. The flame is scaled from the Sandia Flame D with a Reynolds number of 22,400. Three classes of RTE solvers, namely the discrete ordinates method, spherical harmonics method, and Monte Carlo method, are examined. The spectral models include the Planck-mean model, the full-spectrum k-distribution (FSK) method, and the line-by-line (LBL) calculation. The performances of different radiation models in terms of accuracy and computational cost are benchmarked. The results have shown that both RTE solvers and spectral models are critical in the prediction of radiative heat source terms for this jet flame. The trade-offs between the accuracy, the computational cost, and the implementation difficulty are discussed in detail. The results can be used as a reference for radiation model selection in combustor simulations.

David, Chloe↗
DOE OSTI2022

Neural Networks for Supercapacitor Modeling

This software implements a one-dimensional model of a supercapacitor and calibrates it to experimental data. A system of coupled partial differential equation models the salt concentration and the ionic and electronic potentials. The term in the equation modeling the capacitance contains a neural network, whose parameters are optimized to fit the provided experimental data. The software is written in python using the open-source library JAX to calculate the gradients needed for the optimization.

Salazar de Troya, Miquel↗
DOE OSTI2020

AniMACCS User Guide

AniMACCS is a utility code in the MELCOR Accident Consequence Code System (MACCS) software suite that allows for certain MACCS output information to be visually displayed and overlaid onto a geospatial map background. AniMACCS was developed by Sandia National Laboratories for the U.S. Nuclear Regulatory Commission. MACCS is designed to calculate health and economic consequences following a release of radioactive material in the atmosphere. MACCS accomplishes this by modeling the atmospheric dispersion, deposition, and consequences of the release, which depend on several factors including the source term, weather, population, economic, and land-use characteristics of the impacted geographical area. From these inputs, MACCS determines the characteristics of the plume, as well as ground and air concentrations as a function of time and radionuclide.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2022

Analytic error analysis of cross section interpolation methods in nodal diffusion codes - II: Numerical results

This paper is the second part of a two-part paper that documents the numerical results for the partial derivatives model presented in part I. In this paper, we derive the error bounds for the analytical point-wise error expression and verify our bounds with numerical experiments. The point-wise error expressions make available, and bound, the sources that contribute to the total error of the interpolated cross section in terms of the Lagrange interpolation errors and the model form error. MPACT is used to generate two-group homogenized cross sections for Westinghouse's AP1000 Region 4 lattice to evaluate the accuracy of the bounds. Error bounds calculated over a grid are compared to numerical data for uni-variate and multi-variate interpolation. The point-wise error bounds of a typical case matrix - two branches in each state variable - are displayed for bi-variate interpolation in the state variables: moderator density, fuel temperature, and boron concentration. The error bounds are shown to be highly accurate compared to numerical results, and in accordance with the underlying physics. We then discuss and show how the sources of error contribute to the total error, and consider the improvement of each error source. Finally, we mention future work such as propagating our cross section error bounds through a reactivity calculation. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Update of the Nuclear Criticality Slide Rule Calculations: Plutonium systems – Delayed Fission Gamma

IRSN (France), LLNL (USA) and ORNL (USA) began a long-term collaboration effort in 2015 to update the nuclear criticality Slide Rule for the emergency response to a nuclear criticality accident. The Slide Rule permits the estimation of neutron and gamma dose rates and integrated doses based upon estimated fission yields, as a function of distance from the fission source, and time after criticality accidents for different critical systems. This paper presents results from the fourth phase of the current update of the Slide Rule project, in which delayed fission-product gamma (DFG) dose rates of unreflected plutonium critical systems were compared by several modern 3D radiation transport codes (MCNP, COG, SCALE), using updated flux-to-dose conversion factors. Dose rates are calculated for fissile material at five moderation ratios (H/Pu), at 1 m above the ground as a function of distance (between 30 cm and 1,200 m) from the external surface of the source to the center of the detector, and for periods between 1 s and 1,000 min after the critical instantaneous event. Further efforts have been devoted to the determination of the delayed gamma source, by comparing the time-dependent energy spectra obtained from several methods. Overall, DFG dose rates calculated by each participant led to consistent results. Extra effort is under way to identify the cause of the remaining differences, by comparing precisely the gamma source, and particularly nuclides inventories.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

"Source Term Modeling for Advanced Gas Micro-Reactors"

Maintaining the safety of the public, environment, and operating personnel is the most important factor in designing, operating, maintaining, and decommissioning nuclear reactors. In recent years, there has been a growing interest in the development of micro-reactors employing TRi-structural ISOtropic (TRISO)-coated particle fuel. In gas reactors, TRISO fuel plays an important role in the safety case for high temperature reactors because of the fission product retention properties of the fuel. This ability enables the use of a functional containment strategy for the reactor where multiple barriers are used to prevent fission product release to the environment. Part of the safety analysis of these advanced reactors is the assessment of radionuclide releases under normal and accident conditions through the multiple credited safety barriers. Using conservative assumptions, a mechanistic analysis can be performed to quantify these releases that combines the probabilistic assessment of failure with analytic solutions to radionuclide transport equations. Source term modeling for TRISO fuel has been performed for previous reactor designs; however, these models are outdated, in many cases proprietary, and need updates to be applied to the current state of TRISO fuel technology and alternative gas reactor core configurations [1]. Currently, the only publicly available source term assessment for gas reactors is an expert-based Monte Carlo simulation based on the effectiveness of the fuel kernel, coating layers, and graphite block in a modular high temperature gas reactor [2]. Thus, there is a need to develop a simple, versatile, and mechanistic model of fission product release and transport in gas reactor cores that could be applied to a variety of reactors through user inputs and reactor-specific radionuclide inventories. The release is calculated by the diffusion of the key safety important fission products through the kernel, silicon carbide (SiC), graphite for both intact and defective TRISO particles based on fuel and graphite temperatures in the reactor under normal operation. These releases from the fuel enter the coolant where they can plate-out on cooler surfaces. A clean-up model is included for designs with a coolant purification system to remove fission gases. This initial distribution of fission products in the reactor serves as an initial condition for potential releases under postulated accident conditions. The model then can calculate the fission product release for any transient temperature profile and fission product releases can then be used to assess radiological dose to the workers and the public using conventional dose tools. Data on the diffusion of fission products is based on historic German TRISO experiments and the more current Department of Energy (DOE) Advanced Gas Reactor (AGR) TRISO fuel development program. The model is coded in python with inputs and outputs in excel spreadsheets, as well as python plotting utilities to aid in the interpretation of the results. References: [1] INL, NGNP Mechanistic Source Term White Paper, INL-10-17997, July 2010. [2] David A. Petti, Richard R. Hobbins, Peter Lowry, Hans Gougar, “Representative Source Terms and The Influence of Reactor Attributes on Functional Containment in Modular High Temperature Gas-cooled Reactors,” Nuclear Technology, Vol. 184, p. 181-197, Nov. 2013.

07 ISOTOPE AND RADIATION SOURCES↗