Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear data uncertainties”

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 1,135 records · Page 63

Advanced Simulation and Computing: ASC FY24 Implementation Plan

The DOE National Nuclear Security Administration (NNSA) Stockpile Stewardship Program (SSP) is an integrated technical program for maintaining the safety, security, and reliability of the U.S. nuclear stockpile. The SSP incorporates nuclear test data, computational modeling and simulation, and experimental facilities to advance understanding of nuclear weapons. The suite of data analyzed comes from activities including previous nuclear tests, stockpile surveillance, experimental research, and development and engineering programs. This integrated national program requires the continued use of experimental facilities and the computational capabilities to support the SSP missions. These component parts, in addition to an appropriately scaled production capability, enable NNSA to support stockpile requirements. The ultimate goal of the SSP, and thus of the Advanced Simulation and Computing (ASC) program, is to ensure that the U.S. maintains a safe, secure, and effective strategic deterrent. The ASC program is a cornerstone of the SSP, providing simulation capabilities and computational resources to support the annual stockpile assessment and certification process, study advanced nuclear weapons design and manufacturing processes, analyze accident scenarios and weapons aging, and provide the tools to enable stockpile Life Extension Programs (LEPs) and the resolution of Significant Finding Investigations (SFIs). This work requires a balance of resources, including technical staff, hardware, simulation software, and computer science solutions. The ASC program focuses on increasing the predictive capabilities in a three-dimensional (3D) simulation environment while maintaining support to the SSP. The Program continues to improve its unique tools for understanding and solving progressively more difficult stockpile problems (sufficient resolution, dimensionality, and scientific details), and quantifying critical margins and uncertainties. Resolving each issue requires increasingly difficult analyses because the aging process has progressively moved the stockpile further from the original test base. While the focus remains on the U.S. nuclear weapons program, where possible, the Program also enables the use of high-performance computing (HPC) and simulation tools to address broader national security needs, such as foreign nuclear weapon assessments and nuclear counterterrorism. The 2022 Nuclear Posture Review (NPR) calls for NNSA to “deliver a modern, adaptive nuclear security enterprise based on an integrated strategy for risk management, production-based resilience, science and technology innovation, and workforce initiatives.” Furthermore, “NNSA will establish a Science and Technology Innovation Initiative to accelerate the integration of science and technology (S&T) throughout its activities.” Executing this strategy necessitates the continued emphasis on developing and sustaining high-quality scientific and engineering staff, as well as supporting computational and experimental capabilities. These components constitute the foundation of the nuclear weapons program. The continued success of the SSP and LEPs is predicated upon the ability to credibly certify the stockpile, without a return to underground nuclear tests (UGTs). Shortly after the nuclear test moratorium entered into force in 1992, the Accelerated Strategic Computing Initiative (ASCI) was established to provide an extensive simulation capability to underpin stockpile certification. While computing and simulation have always been essential to the success of the nuclear weapons program, the program goal of ASCI was to execute NNSA’s vision of using these tools in support of the stockpile stewardship mission. The ASCI program was essential to the successful demonstration of the SSP, providing critical nuclear weapons simulation and modeling capabilities. ASCI officially evolved into the ASC program in fiscal year (FY) 2005, but the mission remains essentially the same: provide the simulation and computational capabilities that underpin the ability to maintain a safe, secure, effective nuclear weapon stockpile, without returning to underground nuclear testing. The capabilities that the ASC program provides at the national laboratories play a vital role in the nuclear security enterprise and are necessary for fulfilling the stockpile stewardship and life extension requirements outlined for NNSA. The Program develops modern simulation tools that provide insights into stockpile aging issues, provide the computational and simulation tools that enable designers and analysts to certify the current stockpile and life-extended nuclear weapons, and inform the decision-making process when any modifications in nuclear warheads or the associated manufacturing processes are deemed necessary. Furthermore, ASC is enhancing the predictive simulation capabilities that are essential to evaluate weapons effects, design experiments, and ensure test readiness. The ASC program continues to improve its unique tools to solve stockpile problems— with a focus on sufficient resolution, dimensionality, and scientific detail—to enable Quantification of Margins and Uncertainties (QMU) and to resolve the increasingly difficult analyses needed for stockpile stewardship. The needs of the Stockpile Management and Production Modernization programs (formerly Directed Stockpile Work) also drive the requirements for simulation and computational resources. These requirements include planned LEPs, stockpile support activities, and mitigation efforts against the potential for technical surprise. All of the weapons within the current stockpile are in some stage of the life extension process. The simulation and computational capabilities are crucial for successful execution of these life extensions and for ensuring NNSA can certify these life-extended weapons without conducting a UGT.

97 MATHEMATICS AND COMPUTING↗

Enhancing the Operational Resilience of Advanced Reactors with Digital Twins by Recurrent Neural Networks

Because of a lack of operation data during abnormal and accident scenarios, along with the existence of uncertainty in the evaluation model for transient and accident analysis, the established abnormal and emergency operating procedures can be biased in characterizing the reactor states and ensuring operational resilience. To improve state awareness and ensure operational flexibility for minimizing effects on the system due to anomaly, digital twin (DT) technology is suggested to support operator's decision-making by effectively extracting and using knowledge of the current and future plant states from the knowledge base. To demonstrate DT's capability for recovering the complete states of reactors and for predicting the future reactor behaviors, this paper develops and assesses both the diagnosis and prognosis DTs in a nearly autonomous management and control system for an Experimental Breeder Reactor-II simulator during different loss-of-flow scenarios.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Gamma-Ray Imaging to Monitor Accumulating Holdup

Over time, nuclear materials can accumulate in process equipment and keeping track of such material is of interest to many stakeholder communities including safeguards, process monitoring, criticality safety, and radiological safety. Passive gamma-ray imaging provides a visual and quantitative tool that can be applied to this problem, allowing one to determine the amount of material present from the images. In this work we emulate the accumulation of highly enriched uranium (HEU) in a large rectangular stainless-steel duct by stacking 5, 27 × 50 cm2 cards, each of which has approximately 11 g of 235U in oxide form distributed across its surface, behind each other. A series of images were collected using a coded-aperture gamma-ray imager, starting with a single card, adding a second card, etc. Long integrations were collected to allow evaluating the performance of the approach as a function of integration time. The images are analyzed using first principles to determine the amount of material in each image. The results are compared to the known amount of material in each configuration, and as a function of time. The analysis includes a first-principles error propagation and to further evaluate the uncertainty, a boot-strap sampling approach is also applied to the data. The work offers a proof of concept for unattended monitoring of uranium holdup with a quantitative gamma-ray imaging technique.

Ziock, Klaus-Peter↗

Ab-initio molecular dynamics study of eutectic chloride salt: MgCl2–NaCl–KCl

Ionic liquid materials are viable candidates as a heat transfer fluid (HTF) in a wide range of applications, notably within concentrated solar power (CSP) technology and molten salt reactors (MSRs). For next-generation CSP and MSR technologies that strive for higher power generation efficiency, a HTF with wide liquid phase range and energy storage capabilities is crucial. Studies have shown that eutectic chloride salts exhibit thermal stability at high temperatures, high heat storage capacity, and are less expensive than nitrate and carbonate salts. However, the experimental data needed to fully evaluate the potential of eutectic chloride salts as a HTF contender are scarce and entail large uncertainties. Considering the high cost and potential hazards associated with the experimental methods used to determine the properties of ionic liquids, molecular modeling can be used as a viable alternative resource. In this study, the eutectic ternary chloride salt MgCl 2 –NaCl–KCl is modeled using ab-initio molecular dynamics simulations (AIMDs) in the liquid phase. Using the simulated data, the thermophysical and transport properties of eutectic chloride salt can be calculated: density, viscosity, heat capacity, diffusion coefficient, and ionic conductivity. For an initial model validation, experimental pair-distribution function data were obtained from X-ray total scattering techniques and compared to the theoretical pair-distribution function. Additionally, theoretical viscosity values are compared to experimental viscosity values for a similar system. The results provide a starting foundation for a MgCl 2 –NaCl–KCl model that can be extended to predict other fundamental properties.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Million Worker Study

Radiation health issues have been an important aspect of DOE's Worker Safety and Health programs. The goal is to ensure that workers are adequately protected from the various radiological hazards associated with DOE sites and operations. Since the early 1940’s DOE has supported the conduct of epidemiologic studies of practically every DOE (AEC) facility and collected these data, now managed by the Oak Ridge Associated Universities (ORAU), in remarkable detail. In the early 2000s, the Office of Health, Safety and Security authorized access to specific DOE worker datasets. Shortly thereafter, the DOE Office of Science provided funds for a pilot study that confirmed the feasibility of the Million Worker Study (MWS), and then additional resources were provided (with other agencies) to extend the follow-up of many populations, not just DOE workers, but also atomic veterans, industrial radiographers, nuclear power plant workers and medical radiation workers (the Million Persons Study (MPS)). This proposal was a continuation of support provided by DOE, specifically to extend the follow-up of the worker populations at the Mallinckrodt Chemical Works (MCW) and Los Alamos National Laboratory (LANL). The work addresses DOE’s interest in clarifying the health risks of their workers as well as contributing knowledge on radiation risks that is relevant today with regard to compensation schemes. The findings are also important for the US public in light of the increased population exposure to medical imaging, environmental circumstances such as hydraulic fracturing, increased exposures during high altitude flights, and with regard to nuclear accidents such as Fukushima and possible terrorist events. Furthermore, it is important to consider reducing the uncertainty in current risk estimates by developing risk coefficients based on healthy American workers who are more representative of U.S. workers and the general public than 1945 Japanese survivors of the atomic bombs living in a war-torn country which experienced deprivation, malnourishment, and increased rates of infections and other diseases. But more importantly, it is important to learn whether radiation exposures received gradually over time (e.g., years) are more or less effective in causing health effects, cancer in particular, than if the radiation dose is received all at once in a fraction of a second as experienced by the Japanese atomic bomb survivors. Finally, the ability to evaluate and combine large populations with intakes of radionuclides such as uranium, radium, plutonium, americium and polonium will provide new quantitative knowledge on human health effects that hitherto has not been possible. This cost-efficient study has built on the investments made and foundations laid by investigators and government agencies, including DOE, over the past 30-40 years, which have established early worker cohorts that can now provide answers to questions on the lifetime human health risks associated with low-level radiation exposures. Collaborating institutions included: International Epidemiology Institute, Oak Ridge Associated Universities, Oak Ridge National Laboratory, Los Alamos National Laboratory, Landauer, Inc., and Vanderbilt University. Follow-up of these two populations and the integration of them with the many other cohorts (now a total of 31) in the MPS continues under a separate DOE grant (DE-AU0000046).

61 RADIATION PROTECTION AND DOSIMETRY↗

Conceptual Designs for Irradiation Creep Testing of SiC in HFIR

Understanding irradiation creep of nuclear fuel cladding is important to properly size the initial fuel-cladding gap and understand when pellet-cladding contact is expected to occur due to a combination of fuel swelling and cladding creep-down. Irradiation creep also plays a role in relaxing stresses that develop in-pile. Silicon carbide fiber–reinforced silicon carbide matrix (SiC/SiC) composites are the leading long-term accident-tolerant fuel cladding concept for light-water reactors (LWRs). Although some limited data are available regarding irradiation creep of the individual constituents (fibers, matrix), data regarding irradiation creep of SiC/SiC composites are currently insufficient. Additional data regarding irradiation creep compliance and the rupture lifetime (combination of creep and slow crack growth) are needed to understand material limitations. This work describes the design and development of two irradiation vehicles that are being pursued for testing SiC/SiC concepts in the High Flux Isotope Reactor (HFIR). The first is a passive experiment, referred to as the PRECISE experiment, that leverages the constant coolant pressure of HFIR to compress a metallic bellows and provide a well-characterized load to drive creep in a SiC/SiC dog bone specimen. The total creep strain would be quantified post-irradiation by measuring dimensional changes of the specimen length as well as local dimensional changes within the gauge region. Non-stressed specimens would also be irradiated under the same conditions to provide an indication of dimensional changes due to radiation-induced swelling in the absence of creep. A second, more complex experiment, referred to as the INSITE experiment, is being designed in parallel that would use pneumatics to pressurize a metal bellows and linear variable differential transformers (LVDTs) to measure the specimen displacement in situ during irradiation. Such an experiment would provide significantly more data regarding the evolution of the creep compliance as a function of dose and applied stress within a single experiment but would require significantly more development time and cost to execute. The primary concern with the INSITE experiment is the accuracy, reliability, and expected lifetime of the LVDTs during irradiation at elevated temperatures. Efforts are being made to adjust the experiment design and operating procedure to limit LVDT temperatures and mitigate or otherwise compensate for uncertainties due to factors such as temperature fluctuations, creep in the surrounding structural materials, and drift of the LVDTs. This work describes the experiment designs, thermal and structural analysis that were performed to ensure that the desired temperature and stress conditions can be achieved, some initial sensitivity analyses to predict the evolution of the radiation-induced specimen displacements, and potential sources of uncertainty in the measurements. Out-of-pile testing is being performed in parallel to confirm that the test trains achieve the expected stress states in the specimens and do not result in prohibitive stress concentrators (e.g., in the grip regions) that might risk pre-mature failure. The PRECISE experiments are proceeding toward fabrication and assembly with HFIR insertion planned during fiscal year 2026. The INSITE experiment is progressing toward out-of-pile demonstrations, which will provide more conclusive evidence regarding the feasibility of executing these tests in HFIR or whether alternative displacement monitoring techniques may need to be considered.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Probabilistic Groundwater Modeling of the RDX Plume at Los Alamos National Laboratory to Support Risk Assessment - 20359

A plume of the contaminant hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) with concentrations greater than the New Mexico tap water drinking standard (9.66 ppb) is present in the regional aquifer near the southwestern boundary of Los Alamos National Laboratory (LANL). A risk assessment is performed for exposure to regional aquifer groundwater, with long-term predictions of RDX concentrations provided by a calibrated, probabilistic, numerical fate and transport model. The structure of the model is hierarchical, with the RDX Regional Aquifer (RA) groundwater model acting as the primary tool for analysis of downgradient RDX concentrations. The RA model is deeply informed by the conceptual site model (CSM) and is calibrated using site RDX concentration data and hydraulic head measurements, along with other analyses. The model is calibrated using data through December 2019. Where data are scarce other lines of evidence are used to inform inputs, including the multiphase RDX Vadose Zone (VZ) model and Pipe and Disk analytical screening tool. Model inputs are described with informative prior distributions using a robust approach to development that incorporates all available lines of evidence for every parameter used as an input to the RA model. Calibration is performed using a classical nonlinear optimization routine, which is then used to initialize a Bayesian calibration. The Bayesian calibration constrains the uncertainty in the classical calibration, ultimately providing posterior distributions for all model parameters. The challenges of the calibration include high-dimensional parameter space, including spatially heterogeneous hydraulic conductivities, comparatively sparse data, and low RDX concentrations. Posterior parameter distributions developed in the calibration process are then used for stochastic predictive model runs into the future. The result of the forward runs is spatially and temporally explicit estimates of head and concentration with uncertainty at all points in space and time. The probabilistic modeling approach presented here includes innovative computational and statistical methods that leverage high-performance computing (HPC) resources. It makes use of prior modeling work performed at the chromium plume site in the central LANL area, with extensive updates. The risk assessment will ultimately be used to support decision-making at the site, using multiple appropriately weighted sources of information, as well as uncertainty, sensitivity, and value-of-information analyses. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

In-pile thermal conductivity of uranium dioxide at low burnup

In-pile thermal conductivity of uranium dioxide (UO 2 ) was investigated at low burnup levels (<0.2 kWd/kg-UO 2 ) to elucidate the concurrent effects of the fission damage and the thermal recovery. In-pile experiments were performed in the Halden reactor. The fuel centerline and the cladding temperatures were measured during the experiments, and the power level of specially designed test rodlets were also monitored. The uncertainty of the pellet-cladding gap's thermal resistivity was minimized with a liquid-metal bond, while comparable measurements were made in reference test rodlets without liquid metal bonding. The experimental data were analyzed using the inverse heat transfer approach. The heat conduction equation was solved by applying the measured temperature of the cladding as a boundary condition, and the calculated fuel centerline temperature was compared with the measured fuel centerline temperature to determine the relative matrix resistivity. Results showed that the in-pile matrix resistivity, the phonon-lattice interaction terms of the thermal conductivity, was ~1.5–2.5 times higher than its unirradiated value due to fission-induced damage. Finally, the in-pile annealing caused a significant recovery of the matrix resistivity, and the amount of recovery increased from ~20 to ~33%, while the annealing temperature increased from 700 to 1,000 °C.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Raman and infrared spectra of plutonium (IV) oxalate and its thermal degradation products

For over 80 years, plutonium dioxide has been routinely produced via thermal decomposition of hydrated plutonium(IV) oxalate. Despite the longstanding utility of this process, the chemical structures of starting materials and intermediates produced during this thermal conversion remain ill-defined. To help resolve this uncertainty, we measured high-resolution Raman and infrared spectra of Pu(C 2 O 4 ) 2 ·6H 2 O that was heated to 25, 100, 220, 250, 350, and 450 °C in air. Our measurements show that Pu(C 2 O 4 ) 2 ·6H 2 O has a rich vibrational spectrum with at least 15 Raman bands between 180 cm -1 and 1900 cm -1 and 9 infrared bands between 800 cm -1 and 4000 cm -1 . As Pu(C 2 O 4 ) 2 ·6H 2 O is heated, water is liberated, and the oxalate ligand decomposes to produce plutonium oxycarbide species. When heated to 350 °C or higher, vibrational spectra are consistent with PuO 2 with some residual carbon-containing species. Full vibrational spectra, powder X-ray diffraction, and scanning electron microscopy measurements of Pu(C 2 O 4 ) 2 ·6H 2 O and its thermal degradation products are presented herein along with approximate assignments for observed spectral bands. These data can be used to validate and potentially improve existing computational models that describe the chemical structure of compounds produced during thermal degradation of plutonium (IV) oxalate. Given the utility of plutonium (IV) oxalate in synthesizing plutonium dioxide, these results are expected to provide value in the fields of nuclear fuel processing, nuclear nonproliferation, and nuclear forensics.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

STAT7 v1.2 Verification and Validation Report

This report documents the software testing which has been performed for the STAT7 Version 1.2 (hereafter STAT7 v1.2) software, a code for the statistical propagation of uncertainties in steady-state thermal hydraulics (TH) analysis of plate-fueled research and test reactors. The key capabilities addressed in this testing were selected based on the use of the software for analysis work supporting the conversion to low-enriched uranium fuel of U.S. High-Performance Research Reactors (USHPRR) such as the Massachusetts Institute of Technology Research Reactor (MITR-II). This testing scope includes comparison of the code’s statistical processing capabilities using statistical tests, comparison of coolant property evaluations to NIST data tables, and systematic comparison of key TH solver capabilities to the results of hand calculations and to predictions from PLTEMP/ANL Version 4.4. The good agreement achieved in this testing campaign, within the acceptance criteria of a 1% relative difference, also demonstrates that these capabilities have been implemented in STAT7 v1.2 correctly.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Spectro-Timing Study of GX 339-4 in a Hard Intermediate State

We present an analysis of Nuclear Spectroscopic Telescope Array observations of a hard intermediate state of the transient black hole GX 339-4 taken in 2015 January. With the source softening significantly over the course of the 1.3 day long observation we split the data into 21 sub-sets and find that the spectrum of all of them can be well described by a power-law continuum with an additional relativistically blurred reflection component. The photon index increases from approx. 1.69 to approx. 1.77 over the course of the observation. The accretion disk is truncated at around nine gravitational radii in all spectra. We also perform timing analysis on the same 21 individual data sets, and find a strong type-C quasi-periodic oscillation (QPO), which increases in frequency from approx. 0.68 to approx. 1.05 Hz with time. The frequency change is well correlated with the softening of the spectrum. We discuss possible scenarios for the production of the QPO and calculate predicted inner radii in the relativistic precession model as well as the global disk mode oscillations model. We find discrepancies with respect to the observed values in both models unless we allow for a black hole mass of approx. 100 Mass compared to the Sun, which is highly unlikely. We discuss possible systematic uncertainties, in particular with the measurement of the inner accretion disk radius in the relativistic reflection model. We conclude that the combination of observed QPO frequencies and inner accretion disk radii, as obtained from spectral fitting, is difficult to reconcile with current models.

accretion disks↗

Benchmarking of Different Inverse Point Kinetics Implementations for an Autocorrected Reactimeter Algorithm

In November 2017, the Transient Reactor Test Facility returned to operation. Since that time, many transient test series have been completed; each has provided valuable data for materials performance, reactor safety that can be applied in future designs. During each experimental series, detector count rates provided important information on the core behavior during transients. However, a limitation of these data is that variations in the neutron distribution during experiments cause errors when attempting to infer reactivity evolution from detector signals. Neutron physics codes can be used to compute the flux shape variations. However, this is a poor solution when the experimental data is used to do verification, validation and uncertainty quantification (VVUQ) on codes. Indeed, if the output of the code is used both as a reference and to correct what the reference is compared to, the circular dependency limits the quality of the VVUQ approach. To overcome this problem, an Autocorrected Reactimeter Algorithm (ACRA) has been developed. This approach infers time-dependent reactivity evolution by testing different spatial corrections and selecting the one that minimizes reactivity variations when the core is in a frozen configuration (i.e. when there is no variation in parameters affecting reactivity). However, the scope of this method was limited to transients where there were negligible thermal feedback. Indeed, the core is never in a frozen configuration when the fuel temperature varies during the whole transient. This is our motivation for the development of an improved version of the ACRA which does not require frozen configurations

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Measurement of the inclusive isolated-photon production cross section in pp collisions at $\mathbf {\sqrt{\textit{s}}=13}$ TeV

The production cross section of inclusive isolated photons has been measured by the ALICE experiment at the CERN LHC in pp collisions at centre-of-momentum energy of $\mathbf {\sqrt{\textit{s}}=13}$ TeV collected during the LHC Run 2 data-taking period. The measurement is performed by combining the measurements of the electromagnetic calorimeter EMCal and the central tracking detectors ITS and TPC, covering a pseudorapidity range of and a transverse momentum range of |η γ | < 0.67 and a transverse momentum range 7 < p$^{γ}_{T}$ < 200 GeV/c. The result extends to lower p$^{γ}_{T}$ and x$^{γ}_{T}$ = 2p$^{γ}_{T}$/$\sqrt{s}$ ranges, the lowest x$^{γ}_{T}$ of any isolated photon measurements to date, extending significantly those measured by the ATLAS and CMS experiments towards lower p$^{γ}_{T}$ at the same collision energy with a small overlap between the measurements. The measurement is compared with next-to-leading order perturbative QCD calculations and the results from the ATLAS and CMS experiments as well as with measurements at other collision energies. The measurement and theory prediction are in agreement with each other within the experimental and theoretical uncertainties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Null Space Monte Carlo Evaluation of the Plateau to River Model

The Plateau to River Groundwater Model (P2R Model) is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by CH2M HILL Plateau Remediation Company at the Hanford Site in Washington State. Figure 1-1 illustrates the P2R Model extents, discretization, and boundary conditions. The P2R Model provides a computational framework to simulate the F&T of contaminants in groundwater associated with the 200-PO-1, 200-UP-1, 200-BP-5, and 200-ZP-1 Groundwater Operable Units (OUs) in the Hanford Site Central Plateau. In addition, the model includes adjacent areas and facilities (e.g., the State Approved Land Disposal Site). Intended and anticipated uses of the model include calculating water levels, hydraulic gradients, and groundwater flows throughout the model domain (encompassing the 200 West and 200 East Areas) for use in subsequent F&T calculations for contaminants of concern and developing scale-appropriate, telescopic-mesh refinement models for detailed evaluation of areas within the model domain where required. The overall objective of the modeling effort is to provide a basis for making informed remedial action decisions based on descriptions of current and expected future contaminant concentrations in groundwater at decision points within the OU boundaries. The objective for the model development phase is to create a common modeling platform that can be used for investigations of the Central Plateau groundwater OUs and areas downgradient toward the Columbia River. The P2R Model calibration to historical data observed at the Hanford Site is documented in CP-57037, Model Package Report for the Plateau to River Model Version 8.3. The purpose of this environmental calculation is to describe a null space Monte Carlo (NSMC) evaluation was conducted with the historic calibration of the P2R model. Use of numerical groundwater models is always accompanied with uncertainty in the results produced by a model because models are approximations of reality. Thus, by definition, lack the detail to fully represent observed behavior. Use of numerical techniques, such as a NSMC analysis, can help in identifying and quantifying the potential uncertainties associated with a numerical model such as the P2R Model. Use of the NSMC approach results in 100 groundwater flow models that are variants of the calibrated P2R Model. These variant models can be used to evaluate uncertainty in model predictions made by the calibrated P2R Model for other analyses. A secondary purpose of the environmental calculation is to establish these variant models for use with other applications.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Implementation of the D1S Methodology for Shutdown Dose Rate Calculations in the OpenMC Monte Carlo Particle Transport Code

We present an implementation of the direct one-step (D1S) methodology for shutdown dose rate (SDR) calculations in the OpenMC Monte Carlo particle transport code. In addition to being the first fully open-source D1S implementation, it is also the first to require no ad hoc source code or nuclear data library modifications. The code can seamlessly switch between production of prompt and decay photons based on a user input parameter, and the decay data needed for decay photon generation are made available through a depletion chain file, which is already used for OpenMC’s built-in depletion/activation solver. A set of Python functions significantly eases the burden of computing and applying time correction factors needed to properly account for the time dependence of radionuclide activity. To assess the accuracy of the D1S implementation, SDR calculations have been carried out for three problems: a prism of iron irradiated by 14-MeV neutrons, the ITER port plug computational benchmark, and the Frascati Neutron Generator (FNG) ITER dose rate benchmark problem from the Shielding INtegral Benchmark Archive and Database (SINBAD). For each of these problems, comparisons were made to calculations using the rigorous two-step (R2S) method. The results on the iron prism problem illustrate how the D1S method achieves superior spatial resolution compared to the R2S method without the need for spatial discretization of the activation regions. The D1S and R2S results for the ITER port plug benchmark agree well with previously reported results in the literature. While the D1S results are 10% to 15% lower than the R2S results, this may be due to stochastic uncertainty and/or spatial discretization in the R2S calculations. On the FNG dose rate benchmark problem, the D1S method produces dose rate estimates that are within 4% of the dose rates predicted using a cell-based R2S workflow. The D1S estimates of the SDR are also in reasonable agreement with the experimental measurements and show the same basic trends that have been observed in previous works. A qualitative analysis of the execution time and uncertainty for the R2S and D1S workflows suggests that the D1S method would attain a higher figure of merit.

D1S method↗

Computational Evaluation of Thermal Response of Open-Cell Foam With Circular Pore (Computional Evaluation of Thermal Response of Open-Cell Foam with Circular Pore)

The evaluation of effective material properties in heterogeneous materials (e.g., composites or multicomponent structures) critically relevant to a wide spectrum of applications, including nuclear power, electronic packaging, flame retardants, hypersonics, and gas turbine power. The work described in this paper is centered around the numerical assessment of the thermal behavior of porous materials obtained from finite element thermal modeling and simulation. Here, two-dimensional, steady state analyses were performed on unit cells with centered, circular pores using a second order accurate Galerkin finite element method (FEM). The effective thermal conductivities of the porous systems were examined, encompassing a range of porosities from 4.9% to 60.1%. The geometries of the models were generated based on ordered circular pores for each modeled porosity level. The system response quantity (SRQ) under investigation was the dimensionless effective thermal conductivity across the unit cell. The dimensionless effective thermal conductivity was compared across all simulated cases, producing a trend between porosity and effective thermal conductivity. In the presented investigation, the method of manufactured solutions (MMS) was used to perform code verification, and the grid convergence index (GCI) was employed to estimate discretization uncertainty as solution verification. Code verification concluded an approximately second order accurate Galerkin FEM solver. It was found that the introduction of porosity to the unit cell material structure reduces effective thermal conductivity, as anticipated. Numerical results obtained in this study are compared to an analytical solution and to a sample of empirical data.

36 MATERIALS SCIENCE↗

RAVEN Theory Manual

RAVEN is a software framework able to perform parametric and stochastic analysis based on the response of complex system codes. The initial development was aimed at providing dynamic risk analysis capabilities to the thermohydraulic code RELAP-7, currently under development at Idaho National Laboratory (INL). Although the initial goal has been fully accomplished, RAVEN is now a multi-purpose stochastic and uncertainty quantification platform, capable of communicating with any system code. In fact, the provided Application Programming Interfaces (APIs) allow RAVEN to interact with any code as long as all the parameters that need to be perturbed are accessible by input files or via python interfaces. RAVEN is capable of investigating system response and explore input space using various sampling schemes such as Monte Carlo, grid, or Latin hypercube. However, RAVEN strength lies in its system feature discovery capabilities such as: constructing limit surfaces, separating regions of the input space leading to system failure, and using dynamic supervised learning techniques. The development of RAVEN started in 2012 when, within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the need to provide a modern risk evaluation framework arose. RAVEN’s principal assignment is to provide the necessary software and algorithms in order to employ the concepts developed by the Risk Informed Safety Margin Characterization (RISMC) program. RISMC is one of the pathways defined within the Light Water Reactor Sustainability (LWRS) program. In the RISMC approach, the goal is not just to identify the frequency of an event potentially leading to a system failure, but the proximity (or lack thereof) to key safety-related events. Hence, the approach is interested in identifying and increasing the safety margins related to those events. A safety margin is a numerical value quantifying the probability that a safety metric (e.g. peak pressure in a pipe) is exceeded under certain conditions. Most of the capabilities, implemented having RELAP-7 as a principal focus, are easily deployable to other system codes. For this reason, several side activates have been employed (e.g. RELAP5-3D, any MOOSE-based App, etc.) or are currently ongoing for coupling RAVEN with several different software. The aim of this document is to provide a set of commented examples that can help the user to become familiar with the RAVEN code usage.

97 MATHEMATICS AND COMPUTING↗

Model Package Report: Central Plateau Vadose Zone Models

This model package report describes concisely the modeling objectives, conceptualization, implementation, uncertainty and sensitivity, configuration control, limitations, data needs, and recommendations for improvement of the Central Plateau Vadose Zone Models. This collection of models is developed to meet the vadose zone simulation needs of the Hanford Site Composite Analysis (CA), the Hanford Site Cumulative Impact Evaluation (CIE), and are expected to find other applications in Hanford Site remedial cleanup decision-making processes. This model package report describes and documents the development of the models themselves to fulfill technical approaches defined for the CA and the CIE. This report does not document any specific calculation using these models: instead, applications of individual vadose models to perform specific calculations will be documented in environmental calculation files, including inputs and results, as appropriate. Configuration management for the Central Plateau Vadose Zone Models is managed through assignment of unique sequential numerical version numbers and archival of the models in the Environmental Modeling Management Archive.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗