Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “subcritical experiments”

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 181 records · Page 10

Evaluation of Neutron Detection Technology and Analysis Techniques to Support Fukushima Daiichi Decommissioning [Poster]

Decommissioning efforts and associated fuel debris retrieval from the damaged Units 1-3 require robust methods to monitor and assure subcriticality. These include adequate detection technologies capable to withstand high gamma radiation conditions within the units; as well as analysis methods capable to provide robust signatures to assure agile criticality monitoring in situation where material distribution is changing. This poster presents overview of experimental and analytical efforts focused on evaluation of neutronbased technologies and development of real-time analysis using state-of-the art capabilities available within the 1F Fuel Retrieval and Monitoring Experiments (1FRAME) project.

1FRAME↗

Chemical controls on the propagation and healing of subcritical fractures

Human activities involving subsurface reservoirs—resource extraction, carbon and nuclear waste storage—alter thermal, mechanical, and chemical steady-state conditions in these systems. Because these systems exist at lithostatic pressures, even minor chemical changes can cause chemically assisted deformation. Therefore, understanding how chemical effects control geomechanical properties is critical to optimizing engineering activities. The grand challenge in predicting the effect of chemical processes on mechanical properties lays in the fact that these phenomena take place at molecular scales, while they manifest all the way to reservoir scales. To address this fundamental challenge, we investigated chemical effects on deformation in model and real systems spanning molecular- to centimeter scales. We used theory, experiment, molecular dynamics simulation, and statistical analysis to (1) identify the effect of simple reactions, such as hydrolysis, on molecular structures in interfacial regions of stressed geomaterials; (2) quantify chemical effects on the bulk mechanical properties, fracture and displacement for granular rocks and single crystals; (3) develop initial understanding of universal scaling for individual displacement events in layered geomaterials; and (4) develop analytic approximations for the single-chain mechanical response utilizing asymptotically correct statistical thermodynamic theory. Taken together, these findings advance the challenging field of chemo-mechanics.

36 MATERIALS SCIENCE↗

Sensitivity/Uncertainty Comparison Study Involving IRSN, LANL, and ORNL Tools to Support Validation

Under a DOE Nuclear Criticality Safety Program (NCSP) task involving Analytical Methods, three Laboratories collaborated in a comparison of results obtained from Sensitivity/Uncertainty (S/U) packages relevant to validation of transport codes. The task involves Institut de Radioprotection et de Sûreté Nucléaire (IRSN), Los Alamos National Laboratory (LANL), and Oak Ridge National Laboratory (ORNL) comparing results of MORET 5/MACSENS V3.0, MCNP6.2/Whisper-1.1, and SCALE 6.2.3/TSUNAMI/USLSTATS respectively. All Monte Carlo transport code results utilize nuclear data from ENDF/B-VII.1 evaluation. This study examines five cases from the International Handbook of Evaluated Criticality Safety Benchmark Experiments (ICSBEP Handbook) selected as application models: IEU-MET- FAST-002-001, LEU-COMP-THERM-001-001, LEU-SOL-THERM-004-001, MIX-COMP- THERM-001-001, and U233-SOL-THERM-001-001. This is a continuation of a previous study to examine Pu and HEU cases: HEU-MET-FAST-013-001, HEU-SOL-THERM-001-008, PU-MET- FAST-022-001, and PU-SOL-THERM-001-001. Ultimately, comparison is made between Upper Subcritical Limits (USLs) obtained using each code package for each application case. Since differences exist in whether packages take into account margin of subcriticality (MOS), the USL is computed using only bias and bias uncertainty, also known as the calculational margin (CM) in ANSI/ANS-8.24. Results comparison appears to show that benchmark selection has a greater influence on the USL than the method used for calculation of bias and bias uncertainty.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Current Progress of the Final Design of a Subcritical Assembly at the Oak Ridge National Laboratory

During the past decade, the US Department of Energy (DOE)/National Nuclear Security Administration (NNSA) Nuclear Criticality Safety Program (NCSP) has been coordinating nuclear criticality safety (NCS) courses to provide training and qualification to new staff, supervisors, managers, regulators, and other professionals involved with NCS. To address the growing need for more NCS training possibilities, a subcritical assembly was designed at Oak Ridge National Laboratory (ORNL). The Oak Ridge Subcritical Assembly (ORSA), operated at ORNL, could be used in a new course for fissile material handlers/operators or as a backup training facility for currently existing NCSP hands-on NCS training. This paper describes some of the design steps necessary for obtaining a final ORSA design, beginning with a successful feasibility study performed in 2020. Numerous parameters within each main element of the assembly were considered, such as the reflector used, the type of fuel coating or plating used, and the number of fuel plates used. The design choices related to each parametric study described in this paper are detailed and justified at the end of each subsection. An overview of the current final design characteristics of ORSA and the experiments’ possibilities is also provided. ORSA will be built with fuel plates from legacy Aerojet General Nucleonics (AGN)-201M research reactors currently stored at the Y-12 National Security Complex (Y-12). The cylindrical fuel plates comprise UO 2 particles embedded in polyethylene and have a nominal enrichment of 19.5 ± 0.5 wt % 235 U. The basis for the final design was the feasibility study from 2020. In the feasibility study, ORSA achieved k eff = 0.95, corresponding to a multiplication factor of 20, with approximately 620 g of 235 U (four 4 cm plates, three 2 cm plates, and two 1 cm plates) and an 8.5 cm thick graphite reflector radially and on top of the assembly. The bottom of the assembly has a graphite thickness of 5 cm, and the chosen graphite density was 2.3 g/cm 3 . All design calculations were performed using SCALE 6.2.4 KENO-VI and with the ENDF/B VII.1 Continuous Energy cross-section library.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

IER-519 CED-2: Final Design for Thermal/Epithermal eXperiments (TEX) with Absorbers to Provide Validation Benchmarks for Hanford Tank Farms

The Hanford tank farms contain 56 million gallons of waste across 177 tanks. The primary criticality safety concern for the waste is the plutonium inventory in waste solids – approximately 670 kg in total. Criticality safety analysis credits the absorption and dilution properties of the large quantities of other elements (aluminum, chromium, iron, manganese, nickel, silicon, sodium, and zirconium) present in the waste. Of these, iron and manganese are by far the most significant neutron absorbers, particularly for the waste compositions of highest criticality safety concern. The criticality safety analyses at the Hanford Waste Treatment and Vitrification Plant (WTP) and the Savannah River tank farms also credit iron and manganese as the primary neutron absorbers to demonstrate subcriticality.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Overview of the Tolerance Limit Calculations with Application to TSURFER

To establish confidence in the results of computerized physics models, a key regulatory requirement is to develop a scientifically defendable process. The methods employed for confidence, characterization, and consolidation, or C3, are statistically involved and are often accessible only to avid statisticians. This manuscript serves as a pedagogical presentation of the C3 process to all stakeholders—including researchers, industrial practitioners, and regulators—to impart an intuitive understanding of the key concepts and mathematical methods entailed by C3. The primary focus is on calculation of tolerance limits, which is the overall goal of the C3 process. Tolerance limits encode the confidence in the calculation results as communicated to the regulator. Understanding the C3 process is especially critical today, as the nuclear industry is considering more innovative ways to assess new technologies, including new reactor and fuel concepts, via an integrated approach that optimally combines modeling and simulation and minimal targeted validation experiments. This manuscript employs intuitive, analytical, numerical, and visual representations to explain how tolerance limits may be calculated for a wide range of configurations, and it also describes how their values may be interpreted. Various verification tests have been developed to test the calculated tolerance limits and to help delineate their values. The manuscript demonstrates the calculation of tolerance limits for TSURFER, a computer code developed by the Oak Ridge National Laboratory for criticality safety applications. The goal is to evaluate the tolerance limit for TSURFER-determined criticality biases to support the determination of upper, subcritical limits for regulatory purposes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Interaction between dissolution and precipitation during olivine carbonation: Implications for CO 2 mineralization

Large-scale carbonation of olivine is considered a promising approach for in situ mineral carbonation, offering a permanent and stable method for CO 2 storage. A critical aspect of this process is understanding how dissolution and precipitation interact, as this could drive fracturing and enhance further reactions. In this study, we conducted carbonation experiments on olivine using CO 2 -saturated aqueous solutions of NaHCO 3 and NaCl. Two experimental setups were used: one representing an open geochemical system and the other a closed system, corresponding to reaction-limited and flow-limited scenarios, respectively. Further, post-reaction textural analysis using scanning electron microscopy (SEM) revealed surface coatings of reaction products in the closed system, while etch pits and etch channels were prevalent in the open system. Although no direct evidence of reaction-driven fracturing was observed, etch pits and etch channels may serve as initiation points for subcritical crack formation and growth, potentially maintaining permeability and exposing new unreacted surfaces. Using linear elastic fracture mechanics (LEFM) model, we estimate that microcracks could propagate under a pressure of 0.1 GPa if reaction products accumulate within the etch pits. Our findings offer new insights into the mechanisms governing olivine carbonation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

TEX-HEU: Integral Experiment Execution with Polyethylene at Very Low Temperatures

The low-temperature variant of the TEX HEU (called Low-Temperature TEX or sometimes LT TEX) campaign is a highly anticipated and necessary experimental series by the greater nuclear science community. Fundamentally, the need for low-temperature integral experiments is required to perform validation of cross sections below room temperature. There has been substantial international interest in low-temperature benchmarks to validate below room temperature cross sections, namely talks given at the 2019 International Conference on Nuclear Criticality (ICNC): UK (Watson, 2019), France (Milin, 2019), and UK (Gan & Wilson, 2019). Additionally, NCSP funded thermal scattering laws (TSLs) were produced by North Carolina State University and require low-temperature benchmarks to validate them. Validation of low-temperature cross sections is also necessary for criticality safety applications. One particularly important application is to ensure that during transportation, fissile materials must remain subcritical under normal ambient conditions which is defined as temperatures down to -40°C/°F by the United States 10 CRF 71 as well as a regulation put forward by the International Atomic Energy Agency (IAEA).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Review of Neutron Yield from (α, n) Reactions: Data, Methods, and Prospects

Understanding the radiogenic neutron production rate through the (α, n) reaction is crucial in many areas of physics, including dark matter searches, neutrino studies, and nuclear astrophysics. In addition to its relevance for fundamental research, the (α, n) reaction also plays a significant role in nuclear energy technologies, for example by contributing to neutron production in subcritical systems using UO2, as well as in applications such as medical physics. This review examines the current state of (α, n) yield calculations and neutron spectra, describes the computational tools used for their estimation, and discusses the available cross-section data. We discuss the uncertainties affecting (α, n) yield estimations and propose a strategy to enhance their accuracy. Furthermore, this paper discusses and emphasizes the need for new measurements of (α, n) cross-sections for a variety of relevant materials. Such measurements are essential for improving neutron flux predictions, which are crucial for reducing uncertainties in sensitivity estimates for next-generation physics experiments operating in the keV—MeV range.

(α, n)reactions↗

A closer look at turbulence spreading: How bistability admits intermittent, propagating turbulence fronts

In magnetic fusion plasmas, mounting evidence indicates the possibility of sustained turbulence below the linear stability threshold or more generally global turbulence bistability. The usual reduced models for turbulence spreading are unistable/supercritical and incompatible with this result. The older models further cannot realistically support fronts connecting laminar and turbulent domains. In this work, a minimal model for “subcritical” turbulence spreading is introduced and analyzed. The model may be viewed as phenomenological or derived directly by considering the effect of profile corrugations in an E×B staircase. The model, which is related to the FitzHugh–Nagumo system, supports the robust coexistence of multiple turbulence levels via bistability. We show that this model predicts stronger penetration of turbulence into a linearly stable region as well as the formation of intermittent turbulence fronts that resemble avalanches. We derive the critical size that a localized slug of turbulence must exceed in order to spread. Lastly, we make a prediction of global hysteretic behavior associated with the bistability, which should be testable via experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Insights into ENSO Diversity from an Intermediate Coupled Model. Part II: Role of Nonlinear Dynamics and Stochastic Forcing

Abstract In this study, we investigate how a single leading linear El Niño–Southern Oscillation (ENSO) mode, as studied in Part I, leads to the irregular coexistence of central Pacific (CP) and eastern Pacific (EP) ENSO, a phenomenon known as ENSO spatiotemporal diversity. This diversity is fundamentally generated by deterministic nonlinear pathways to chaos via the period-doubling route and, more prevailingly, the subharmonic resonance route with the presence of a seasonally varying basic state. When residing in the weakly nonlinear regime, the coupled system sustains a weak periodic oscillation with a mixed CP/EP pattern as captured by the linear ENSO mode. With a stronger nonlinearity effect, the ENSO behavior experiences a period-doubling bifurcation. The single ENSO orbit splits into coexisting CP-like and EP-like ENSO orbits. A sequence of period-doubling bifurcation results in an aperiodic oscillation featuring irregular CP and EP ENSO occurrences. The overlapping of subharmonic resonances between ENSO and the seasonal cycle allows this ENSO irregularity and diversity to be more readily excited. In the strongly nonlinear regime, the coupled system is dominated by regular EP ENSO. The deterministic ENSO spatiotemporal diversity is thus confined to a relatively narrow range corresponding to a moderately unstable ENSO mode. Stochastic forcing broadens this range and allows ENSO diversity to occur when the ENSO mode is weakly subcritical. A close relationship among a weakened mean zonal temperature gradient, stronger ENSO activity, and more (fewer) occurrences of EP (CP) ENSO is noted, indicating that ENSO–mean state interaction may yield ENSO regime modulations on the multidecadal time scale.

Geng, Licheng↗

Laboratory Experiments Examining the Effect of Thermal and Mechanical Processes on Hydraulic Transmissivity Evolution

Using laboratory slide-hold-slide experiments, at temperatures from 22 to 200 degrees C, to examine effects of fracture reactivation and quasi-static loading on the evolution of fluid transport properties of simulated fractures in Westerly granite. At all temperatures, the in-plane hydraulic transmissivity consistently decays during hold periods resulting in an overall reduction in transmissivity. During the first three to fifteen hours of an experiment, transmissivity decreases rapidly due to the generation of wear products, development of a sliding surface, and compaction of the resulting gouge. Once the sliding surface has developed, the long-term transmissivity decay rate at 22 and 100 degrees C is significantly lower than the transmissivity decay rate during the initial 3-15 hours of the experiment. However, at 200 degrees C, the decay of hydraulic transmissivity remains high throughout the experiment. The long-term decay of hydraulic transmissivity can be fitted with a power law model with more rapid reduction of hydraulic transmissivity at higher temperature. Periods of sliding on the fracture surface result in transient increases in the transmissivity, due to shear dilation, as is expected for Coulomb materials. These transients are superimposed on the long-term decay. When sliding ceases and a new hold period commences, there is a rapid reduction in transmissivity and return to the long-term rate of transmissivity decay. The rate of decay of the transmissivity transients is inversely correlated with temperature, in contrast to the long-term decay and the expected behavior for processes like subcritical crack growth and indentation creep. The higher decay rates that are observed during the initial 3-15 hours of the tests and following sliding, are associated with times that the porosity of the gouge is expected to be high. The difference in decay rates suggests that when the gouge is driven far from equilibrium by active shearing, densification may be dominated by a different mechanism from long-term compaction.

15 GEOTHERMAL ENERGY↗

Electrochemical techniques and mechanisms for the corrosion of metals and alloys in sub- and supercritical aqueous systems

Corrosion has become a decisive factor hindering supercritical water technologies' long-cycle, safe implementation. Electrochemical experiments performed in supercritical aqueous systems can capture essential, real-time information about the micro/nanoscale processes of corrosion and their kinetics. As a starting point, this paper first summarizes the critical chemical properties of supercritical aqueous systems and analyzes their role in determining the duality of corrosion mechanisms (whether electrochemical or chemical micro-processes predominate). Next, it reviews the available electrochemical techniques such as pH sensors, electrochemical noise measurements, and related applications for quantitatively elaborating the corrosion issues in supercritical aqueous systems. With thick and heavy colors, a series of mechanistic models, such as the mixed potential model and the coupled environment models were introduced from the fundamental theories to the related applications, which have been applied widely in ambient and subcritical aqueous systems and exhibit great values and potential to be employed in high-density supercritical water environments. The passivity, depassivation, and passivity breakdown of metals and alloys in terms of the Point Defect Model (PDM) were also successfully implemented and discussed. Further, the newly built SCW_PDM describes the film growth kinetics at the atomic level of metals and alloys in supercritical water and is being developed to quantitatively analyze and resolve a series of corrosion issues such as the passivity breakdown and the stress corrosion cracking that occur in supercritical aqueous systems. This review is helpful by promoting the development of the electrochemical corrosion theory in supercritical aqueous systems and by identifying corrosion prevention and control methods for the safe implementation of relevant supercritical water technologies.

36 MATERIALS SCIENCE↗

Development of an accelerator-based neutron source to prototype Mo-99 production, Part II: A liquid LBE loop under a high vacuum

To provide US domestic supply of Mo-99 without using high-enriched uranium (HEU), a subcritical uranium target assembly (UTA) is irradiated by an accelerator-based neutron source to create Mo-99 through fission. Part I of this work discusses the design of a liquid lead–bismuth eutectic (LBE) windowless target for an accelerator-based neutron source development. Part II discusses how to couple this windowless target to an accelerator operating at an ultra-high vacuum and the subcritical UTA cooled by water at room temperature. Due to the windowless design of the target, the liquid LBE flow shares an ultra-high vacuum (<1.3 × 10 -7 Pa or 10 -9 Torr) space with the accelerator. As a result of this shared vacuum space, the LBE system must operate at a high vacuum (10 -3 ~10 -6 Pa or 10 -5 ~10 -8 Torr). A magnetic rotary motion feedthrough unit utilizes magnetic fluid to allow rotation of the pump while maintaining a high vacuum environment. Prior to testing the LBE system under vacuum, a pump curve measurement is performed to estimate flowrate in the system. This measurement also generates data on orifice loss coefficients, which are compared to correlations in literature. The second experiment investigates vacuum level in the LBE system during operation. High vacuum is maintained (10 -3 ~10 -5 Pa or 10 -5 ~10 -7 Torr) during system operation, and a residual gas analyzer (RGA) scan shows that partial pressures of residual gases in the LBE system lower over the duration of LBE system operations. The third experiment investigates the gravity driven liquid LBE flowing out of the target chamber in the return line, which is partially full. If the liquid LBE is not drained quickly enough, flooding in the target chamber could occur. The coefficient n in the Manning equation is found to be around 0.008 s/m 1/3 . The last experiment performed is a demonstration that a vacuum jacket could provide sufficient thermal insulation to allow coupling between 300 °C LBE loop and a water tank at room temperature. In conclusion, the results from these experiments have influenced the development of the neutron source for the future commercial scale Mo-99 production system.

43 PARTICLE ACCELERATORS↗

TRIPOLI-4 neutronics calculations for IAEA-CRP benchmark of CEFR start-up tests using new libraries JEFF-3.3 and ENDF/B-VIII

Sodium-cooled Fast Reactors (SFR) are one of the promising Generation IV fast reactors. The China Experimental Fast Reactor (CEFR) is a 65 MWth pool-type SFR with a high neutron leakage core using high enriched uranium oxide fuel. The CEFR start-up tests in 2010 consist of series of neutronics experiments. Essential experimental data are available from the 2018 IAEA-CEFR-CRP benchmark document and useful for the validation of neutron transport codes and nuclear data libraries. The TRIPOLI-4 Monte Carlo transport code is a general-purpose neutronics code using continuous-energy nuclear data libraries. It has a rich validation database covering different computational and experimental benchmark data sets to assure the accuracy and credibility of numerical studies. Previous TRIPOLI-4 SFR core physics calculations used mainly MOX fuels. SFR control rod worth studies with TRIPOLI-4 were mainly on big-size cores. To test the modeling capability of TRIPOLI-4 on fuel loading patterns and on variable control rods positions for SFRs and to check recent nuclear data libraries, it is interesting to investigate the CEFR start-up tests in this work, including core states from subcritical to supercritical, control rods and rod group worth, point kinetics parameters, and radial foil activation measurements. TRIPOLI-4 calculation results using new data libraries JEFF-3.3, ENDF/B-VIII, and those of ENDF/B-VII.1 for different CEFR core states are reported here. By means of different options of TRIPOLI-4, control-rod SAs reactivity worth, point kinetics parameters, and radial {sup 237}Np(n,f) fission rate distribution were successfully evaluated. Both ENDF/B-VIII and JEFF-3.3 nuclear data libraries provided reliable results with a difference of 170 +/- 11 pcm in k{sub eff}.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Development of an accelerator-based neutron source to prototype Mo-99 production, part I: A liquid LBE windowless target

In this article, Molybdenum-99 (Mo-99)’s decay product, technetium-99 (Tc-99 m), is one of the most critical isotopes for medical diagnostics. To provide U.S. domestic supply of Mo-99 without using high-enriched uranium (HEU), a subcritical uranium target assembly (UTA) is irradiated by an accelerator-based neutron source to create Mo-99 through fission. This study discusses the development of the accelerator-based neutron source. The high-energy electrons from the accelerator irradiate a liquid lead-bismuth eutectic (LBE) target to produce neutrons. Part I of this work focuses on numerical and experimental analysis towards the development of a liquid LBE windowless target. Unlike the existing windowless targets in literature, the current design creates a vertical free surface for a beam to irradiate. First, a hydrodynamic analysis of the LBE windowless target is performed. Simplified analytical calculations are assisted by 2D computational fluid dynamics (CFD) simulations to design the target, with the focus on eliminating recirculation zones and avoiding cavitation. With the optimized geometry, the experimental study is performed to investigate the flow hydrodynamics using liquid LBE. The experiments (1) compare pressure drop in the system to correlation predictions; (2) visualize the free surface liquid LBE flow from the beam view; (3) validate the LBE flow profile using temperature sensitive paint from the side view; and (4) validate the liquid LBE film thickness using gamma densitometer measurements. Second, the power handling capability of the designed windowless target is investigated. The divider plate in the current design is susceptible to overheating due to the thin LBE film in front. As LBE erosion and corrosion is likely to occur at an LBE velocity of 2.0 m/s and temperature above 500 °C, a power limit of 10 kW of beam power was established to prevent this corrosion from occurring, which is calculated by a Nusselt number correlation. The divider plate surface temperature at 10 kW agrees well with the 3D CFD simulation results. Part I demonstrates the fundamental physics in liquid LBE windowless target design and associated testing. A companion paper, Part II will demonstrate how to couple this windowless target into the Mo-99 production system, including an accelerator system operating under an ultra-high vacuum and the UTA cooled by water at room temperature.

43 PARTICLE ACCELERATORS↗

EUCLID Sensitivity Database

This report documents the EUCLID sensitivity database along with its several use-cases. EUCLID computed sensitivities for the following integral responses: Criticality of ICSBEP critical assemblies, LLNL pulsed-sphere neutron-leakage spectra, reaction rates in selected ICBSEP critical assemblies, delayed neutron factions of selected ICSBEP critical assemblies, reactivity coefficients in two ICSBEP critical assemblies, sub-critical assembly responses and Rossi-alpha of selected critical assemblies. It is described for each response what the reported observable constitutes, the method we used to obtain the sensitivities, and which integral experiments were studied. It is also documented briefly in what format these sensitivities are stored. These sensitivities were used for many aspects of the EUCLID project, like ML-supported large-scale nuclear-data validation, or optimization of integral experiments. But these sensitivities can also be applied for more established processes in the nuclear-data application field such as adjustment or assessing the upper sub-critical limit.

Delayed Neutron Fraction↗

Effects of Chlorine Capture and a Proposed Density Law on the Reactivity of Plutonium Solution Systems

During fissionable material processing, all normal and credible abnormal conditions must remain safely subcritical. Nuclear Criticality Safety (NCS) uses a number of methods to determine subcriticality, one of which is the use of neutron transport codes such as MCNP6. In order to create models for use with MCNP6, both the geometry and materials in fissionable material processes must be known, or assumptions must be made and quantified for the impact to bias. One of the systems with a significant amount of bias due to material modeling assumptions is in the area of aqueous plutonium processing. These solutions are typically plutonium nitrate solutions or plutonium chloride solutions, which are modeled as fictitious plutonium metal-water mixtures because little is known about the actual density of the solution and there is no current predictive capability approved for use at Los Alamos National Laboratory (LANL) for modeling them. This research is currently underway to fill the gap and develop an algorithm for use with MCNP6 to model the density of plutonium chloride solutions. The method is to be validated with experimental data for density, and also validated with critical experiments using MCNP6. Note that the Chlorine Worth Study (CWS) was performed in December 2021 to help bridge the gap in chlorine data for critical experiments, and is currently awaiting International Criticality Safety Benchmark Evaluation Project (ICSBEP) review. This study was performed by LANL at the National Criticality Experiments Research Center (NCERC) at the Nevada National Security Site (NNSS). Additional information regarding this experiment may currently be found in LA-UR- 22-29180. Additionally, the Chemistry-Actinide Analytical Chemistry (C-AAC) at LANL has performed a number of solution density measurements for PuCl 3 -HC 1 -H 2 O, allowing for such data be used to create a semi-empirical density via the Pitzer method. The published dataset for the measurements is documented in LA-UR-22-25454. This method has already been tested successfully for aqueous plutonium nitrate solutions in SCALE. Current solution density measurements exist of plutonium concentrations of 0-~142g/L, all at 2M HC1, at temperatures 20-40°C. Additional data was taken for HC1-corrected density values, which essentially mimics the data for a pure PuCl x -water solution. The calculations in this report aim to support the current research by demonstrating the difference in system reactivity for the current modeling method when compared to the new proposed modeling with a density law implementation, which is being written as a Python tool to be used with MCNP6.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗