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,153 records · Page 64

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↗

Estimating Cosmological Constraints from Galaxy Cluster Abundance using Simulation-Based Inference

Inferring the values and uncertainties of cosmological parameters in a cosmology model is of paramount importance for modern cosmic observations. In this paper, we use the simulation-based inference (SBI) approach to estimate cosmological constraints from a simplified galaxy cluster observation analysis. Using data generated from the Quijote simulation suite and analytical models, we train a machine learning algorithm to learn the probability function between cosmological parameters and the possible galaxy cluster observables. The posterior distribution of the cosmological parameters at a given observation is then obtained by sampling the predictions from the trained algorithm. Our results show that the SBI method can successfully recover the truth values of the cosmological parameters within the 2σ limit for this simplified galaxy cluster analysis, and acquires similar posterior constraints obtained with a likelihood-based Markov Chain Monte Carlo method, the current state-of the-art method used in similar cosmological studies.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Application of On-Line Monitoring and Real-Time Characterization of Low Level Waste Samples from Hanford Tanks - 20380

Waste and treatment stream characterization is an essential part of safe, efficient, and cost-effective processing of Hanford materials. The demand for materials characterization, including chemical, radionuclide, and physical attributes during waste management, transfer, and staging operations in the Hanford tank farms, is expected to significantly increase as the U.S. Department of Energy Office of River Protection (ORP) progresses toward Hanford waste treatment and immobilization. The application of on-line monitoring can greatly improve operation time frames and reduce costs by providing real-time data and supporting the reduction of necessary grab sample collection. Recently, our team applied Raman online monitoring and chemometric modeling to the supernate of Hanford tank 241-AP-105, and demonstrated the quantitative measurement of nine analytes within this waste. Application of this method has been expanded to additional Hanford tank wastes, AW-102 and AP-107. The performance of the process monitoring instruments and the accuracy of the chemometric models will be assessed. On-line monitoring results will be compared to anion results determined from the current standard, ion-chromatographic and inductively coupled plasma mass spectrometry methods. Comparisons of both determined values and uncertainties of results will be presented. Results will be shared with personnel working on flowsheet interfaces for the Direct-Feed Low Activity Waste mission at Hanford. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Central speed of sound, the trace anomaly, and observables of neutron stars from a perturbative analysis of scaled Tolman-Oppenheimer-Volkoff equations

The central speed of sound (SS) measures the stiffness of the equation of state (EOS) of superdense neutron star (NS) matter. Its variations with density and radial coordinate in NSs in conventional analyses often suffer from uncertainties of the specific nuclear EOS used. Using the central SS and NS mass/radius scaling obtained from solving perturbatively the scaled Tolman-Oppenheimer-Volkoff (TOV) equations, we study the variations of SS, trace anomaly and several closely related properties of NSs in an EOSmodel- independent manner. We find that the SS increases with the reduced central pressure Pˆ c ≡ P c =ε c (scaled by the central energy density εc), and the conformal bound for SS tends to break down for NSs with masses higher than about 1.9M ⊙ . The ratio P=ε is upper bounded as P=ε ≲ 0.374 around the centers of stable NSs. We demonstrate that it is an intrinsic property of strong field gravity and is more relevant than the perturbative QCD bound on it. While a sharp phase transition at high densities characterized by a sudden vanishing of SS in cores of massive NSs are basically excluded, the probability for a continuous crossover signaled by a peaked radial profile of SS is found to be enhanced as Pˆ c decreases, implying it likely happens near the centers of massive NSs. Moreover, a new and more stringent causality boundary as R max =km ≳ 4.73$M$$^{max}_{NS}$ =M ⊙ þ 1.14 for the NS mass-radius curve is found to be excellently consistent with observational data on NS masses and radii. Here, new constraints on the ultimate energy density and pressure allowed in NSs before collapsing into black holes are obtained and compared with earlier predictions in the literature

79 ASTRONOMY AND ASTROPHYSICS↗

Covariance Testing at BNL [Slides]

This presentation discusses testing implemented in FUDGE & ADVANCE in regard to proper ENDF formatting: STAN, STANEF, CHECKR, FUDGE, NJOY, convert to BOXR format, plotting, endfcovreport.py which is an overview of ENDF library’s covariance contents, and covariance Matrix class implements that are the core of all GNDS covariance. Several checks are performed, including symmetric, real, positive definite, condition # of Eigenvalues, and too big uncertainties. The presentation states that it can fix too big and/or small uncertainties and that it can plot matrices. It states that further testing is recommended by CSEWG but has not yet been implemented and suggests methods such as documentation which needs a human to execute, completeness which is partially complete, and reasonableness which has not been started. The presentation concludes by elaborating on reasonableness and states that only cross section is documented, that similar limits on CP reactions cross sections are needed, that maximum limits are needed too which is tricky since ENDF assumes Normal PDFs, and that there may be other observables.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measurement of ( n, γ ) cross section of 90 Zr [Abstract]

The isotopes of Zr with A = [90, 91, 92, 94] make up more than 97% of naturally occurring Zr and are important to many nuclear applications such as nuclear reactors. One of the attractive qualities of naturally occurring Zr isotopes is that they have a low σγ/σt ratio at most neutron energies, such that they improve the neutron economy in reactors by preferentially scattering neutrons rather than absorbing them. This same quality also presents a challenge to measuring the capture cross section, σγ, of Zr isotopes. The ENDF/B-VIII.0 library has a relative uncertainty of approximately 10-20% for incident neutron energies < 0.1 MeV, and uncertainty greater than 20% for energies > 0.1 MeV for the majority of natural Zr isotopes. This motivated the Nuclear Criticality Safety Program (NCSP) to embark on a campaign to accurately measure and evaluate these isotopes of Zr. Here we will show energy-dependent neutron capture cross section measurements for the first enriched sample to be measured: 90 Zr. The measurements of isotopically enriched samples are being carried out at the Geel Electron Linear Accelerator (GELINA) facility of the Joint Research Center - Geel (JRC-Geel) of the European Union. As isotopic enrichment is a costly process we are careful not to activate any of the samples, as this may hinder future radiation-sensitive measurements. The activation analysis is presented in a report by Brown et al.. Once we were satisfied that the Zr samples would not be activated by the measurements, the 90 Zr sample was fabricated at Oak Ridge National Laboratory (ORNL) and shipped to GELINA. The dimensions of the cylindrical sample are approximately 0.12 cm thick and a radius of 2.5 cm. Since 90 Zr is not chemically reactive to air, bare metallic samples were employed. The sample was measured at a flight path (FP) length of 60 m, using four C 6 D 6 detectors on FP14. The final paper will include experimental details and measured cross section data for 90 Zr compared to current evaluated nuclear data libraries.

07 ISOTOPE AND RADIATION SOURCES↗

Evaluation and Uncertainty of Radioxenon Transport with a Mesoscale Model after the February 2013 Underground Test in North Korea

Abstract The transport of radioxenon released from the February 2013 underground nuclear weapons test in North Korea was analyzed at two receptors—one at the Comprehensive Nuclear Test Ban Treaty site Rn58 in Russia (400 km downwind) and a second at Rn38 in Japan (1000 km downwind). Transport was modeled with two ensembles of mesoscale simulations, one generated with varying initial and lateral boundary conditions taken from the Global Forecasting System uncertainty ensemble, and a second created from different parameterizations and surface conditions. The wind variability was similar for the two ensembles and consistent with observations at 925 mb (1 mb = 1 hPa) but not at the surface. Biases in calculated surface winds and the radioxenon concentration in the ensembles were attributed mainly to poor simulation of the sea breeze at both locations and mountain lee affects at Rn38 in Japan. These wind regimes affected the timing of the surface radioxenon plume at Rn58 and its duration at Rn38. Surface wind variability induced by terrain and land–sea contrast (the sea breeze) also had a significant effect on the surface winds and plume dynamics, including blocking of flow approaching elevated terrain near Vladivostok and the west side of Japan. Increased plume uncertainty was seen at night because of surface wind variability. Measured surface chemical variability was larger than found in the first European Tracer Experiment in central Europe. The study found that horizontal model resolution contributes to uncertainty but not as much as vertical resolution, boundary layer parameterizations, and assimilation of surface meteorological data near the receptor.

54 ENVIRONMENTAL SCIENCES↗

Influence of Shear Strength Assumptions on BISON Debonding Simulations

Accurately predicting the thermomechanical response of buffer–IPyC debonding in TRISO fuel particles requires reliable mechanical property inputs for each coating layer, particularly the normal and shear strengths that influence interlayer delamination and stress concentrations. Micro tensile testing of AGR-2 fuel particles provided experimentally measured normal strengths for the buffer, IPyC, and buffer–IPyC interface; however, shear strength was not measured. As a result, BISON simulations of interface debonding must rely on assumed shear strength values, typically estimated as 20–40% of the measured ultimate tensile strength. This study evaluates how these assumed shear strength values influence cohesive zone model (CZM) predictions of buffer–IPyC separation in AGR 2 TRISO particles. Using micro tensile data from three AGR 2 compacts (2 1 3, 5 1 3, and 6 3 3), BISON simulations were performed with multiple shear strength assumptions to quantify their effect on radial and tangential stress evolution, debonding, and gap propagation. The results show that shear strength is a high sensitivity parameter: increasing the assumed shear strength significantly alters the stress distribution at the buffer–IPyC junction, shifts the predicted debonding location, and changes the extent of partial gap formation. While normal strength controls the initiation of interface separation, shear strength strongly influences the mode mixity of the failure process and the resulting stress concentrations transmitted to the IPyC and SiC layers. These findings highlight a critical gap in current TRISO mechanical characterization. Without experimentally measured shear strength, BISON simulations must rely on approximations that introduce uncertainty into predictions of coating layer integrity and fission product barrier performance. Future fuel qualification campaigns should therefore consider measurement of shear strength at the interlayer interfaces to reduce model uncertainty and improve the fidelity of TRISO fuel performance simulations.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A new HALEU isotopic certified reference material

Candidate certified reference material (CRM) U196 was prepared by gravimetrically mixing and dissolving NNSA NRMP CRM 112-A (natural uranium metal) and CRM 116-A (highly enriched 93.2% 235 U metal) to achieve a high-assay low-enriched uranium (HALEU) isotopic abundance reference material. CRM U196 will be certified for the 233 U/ 238 U, 234 U/ 238 U, 235 U/ 238 U, and 236 U/ 238 U ratios. Gravimetrically determined uranium isotope amount ratios and values derived from the isotope data were calculated from the buoyancy corrected mass data and from certified values for the metal CRMs used to create the solution. The 235 U isotope amount fraction (·100) and expanded uncertainty is 19.5846 ± 0.0028. For attribute verification, total evaporation and modified total evaporation isotope ratio measurements were performed. In conclusion, the reference material is to be distributed as a unit containing 5 mg uranium as dry uranyl nitrate.

High-assay low-enriched uranium↗

Modeling Plutonium Decorporation in a Female Nuclear Worker Treated with Ca-DTPA after Inhalation Intake

The present work models plutonium (Pu) biokinetics in a female former nuclear worker. Her bioassay measurements are available at the US Transuranium and Uranium Registries. The worker was internally exposed to a plutonium-americium mixture via acute inhalation at a nuclear weapons facility. She was medically treated with injections of 1 g Ca-DTPA on days 0, 5, and 14 after the intake. Between days 0 and 20, fecal and urine samples were collected and analyzed for 239 Pu and 241 Am. Subsequently, she was followed up for bioassay monitoring over 14 y, with additional post-treatment urine samples collected and analyzed for 239 Pu. The uniqueness of this dataset is due to the availability of: (1) both early and long-term bioassay data from a female with plutonium intake; (2) data on chelation therapy for a female; and (3) fecal measurement results. Chelation therapy with Ca- and/or Zn-salts of DTPA is known to aid in reducing the internal radiation dose by enhancing the excretion of plutonium and americium from the body. Such enhancement affects plutonium biokinetics in the human body, posing a challenge to the internal dose assessment. The current radiation dose assessment practice is to exclude the data affected by Ca-DTPA from the analysis. The present analysis is the first to explicitly model the chelation-affected bioassay data in a female by using a newly developed chelation model. Thus, the bioassay data collected during and after the Ca-DTPA administrations were used for biokinetic modeling and dose assessment. The Markov Chain Monte Carlo method was used to investigate model parameter uncertainty, based on the bioassay data and assumed prior probability distributions. A χ 2 /nData (number of data points) ≈ 1 was observed in this study, which indicates self-consistency of the data with the model. Results of this study show that the worker’s 239 Pu intake was 12 Bq, with a committed effective dose to the whole-body of 1.2 mSv and a committed equivalent dose to the bone surfaces, liver, and lungs of 37.8, 9.1, and 0.8 mSv, respectively. This study also discusses the worker’s dose reduction due to chelation treatment.

61 RADIATION PROTECTION AND DOSIMETRY↗

Demonstration of sub-micron UCN position resolution using room-temperature CMOS sensor

High spatial resolution of ultracold neutron (UCN) measurement is of growing interest to UCN experiments such as UCN spectrometers, UCN polarimeters, quantum physics of UCNs, and quantum gravity. Here we utilize physics informed deep learning to enhance the experimental position resolution and to demonstrate sub-micron spatial resolutions for UCN position measurements obtained using a room-temperature CMOS sensor, extending our previous work that demonstrated a position uncertainty of 1.5 microns. We explore the use of the open-source software Allpix Squared to generate experiment-like synthetic hit images with ground-truth position labels. We use physics-informed deep learning by training a fully connected neural network (FCNN) to learn a mapping from input hit images to output hit position. The automated analysis for sub-micron position resolution in UCN detection combined with the fast data rates of current and next generation UCN sources will enable improved precision for future UCN research and applications.

10B nanometer thin film↗

Picosecond measurements of plastic scintillator pulse shapes from gamma-ray interactions

An understanding of the pulse shape of organic scintillators can provide insight into scintillation mechanisms and inform the selection of the optimal detector material for a given application. Although the timing properties of organic scintillators have been extensively studied, significant discrepancies persist in reported rise and decay times. New plastic scintillating media have also been developed in recent years for which no literature data exist. The goal of this work is to provide high-precision measurements of the pulse shape of a suite of fast plastic organic scintillators from Eljen Technology—EJ-200, EJ-204, EJ-208, EJ-230, EJ-232, and EJ-232Q (with 0.5% benzophenone)—under excitation from $γ$-ray sources. The contributors to the system temporal resolution were quantified, and the dominant source of uncertainty was identified as the determination of the start time of a scintillation event. A pulse shape model was applied to the reconstructed temporal distributions, and rise and decay times were extracted and compared with manufacturer specifications. This work outlines a measurement approach to obtain the pulse shape of organic scintillators in response to $γ$-rays with temporal resolution on the order of tens of picoseconds, yielding data required for modeling plastic scintillator based detection systems for a range of applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Final Design for Thermal/Epithermal eXperiments (TEX) with Chloride Absorbers to Provide Validation Benchmarks for Y-12 Electrorefining Facility

Uranium electrorefining operations at Y-12 require validation for chlorine absorption and reflection. Plutonium chloride solution operations at Los Alamos National Laboratory require chlorine absorption validation. The growing need for chlorine validation is evident to the wider criticality safety community, with multiple attendees at the recent TEX 2.0 meeting at Lawrence Livermore National Laboratory (May 2023) requesting validation for chlorine (Idaho National Laboratory (INL), Institute de radioprotection et de surete nucleaire (IRSN), Savannah River Nuclear Site (SRNS), LANL, and Y-12). Los Alamos National Laboratory has recently performed, and benchmarked, an experiment titled Chlorine Worth Study (CWS) for the internal operations at the lab, but due to the difficulties in precisely characterizing the material compositions of the chlorine absorbers it is advantageous to perform a complimentary study with a different chlorine-based absorber material. Furthermore, having a uranium-based vs plutonium-based experiment provides a separate and important validation basis for criticality safety and nuclear data evaluation. The original final design report for the study of chlorine absorption using the TEX-HEU experimental base was presented in 2022, but used the same chlorine-bearing materials that were found to be difficult to characterize in the LANL benchmark. A complete redesign of the experiment has been performed looking at alternative absorber materials in various forms to produce an experiment that is fully characterizable. This report presents five novel chlorine experiments using the TEX-HEU test bed which provides direct comparison to the Y-12 and INL/Terrapower application needs, utilizing sodium chloride (NaCl) absorber plates. The absorber plates will consist of granulated NaCl (≥99.5% pure), which will be fully encapsulated in aluminum tins, providing a simple but effective chlorine-based absorber material that can be completely characterized. Of the five configurations presented in this report, it is expected that two or three configurations will be down selected for the actual experiment, with the other configurations being alternates. Three of the configurations are in the standard configuration, where the absorber is placed directly on the HEU fuel plates, and two in the sandwich configuration, where the absorber is surrounded by polyethylene moderators to force additional neutron thermalization prior to reaching the absorber. There are two thicknesses of NaCl absorber plates: 3/16” and 1/4” active thicknesses (i.e. not including the encapsulation). Both variations of the absorbers have an active absorber radius of 6” and a total radius of 7.5” to match the diameter of the HEU plates, with the outer 1.5” being aluminum encapsulation. The high-density polyethylene (HDPE) moderators are of the thicknesses: 27/16”, 7/4”, 1/8”, 11/16”, and 3/4”. The final configurations have six (one sandwich and one standard configuration), eight (one sandwich and one standard configuration), and 18 total fuel layers (standard configuration). The standard and sandwich configurations were designed such that the differences in moderator thicknesses are 1” HDPE, which were already procured for the original CED-2. The proposed configurations were precisely tuned to closely match the sensitivity profiles and neutron spectra of the Y-12 upset cases and were also compared to the INL/Terrapower upset cases. The assessment of experimental uncertainties of the non-absorber components was predicted to be 0.00114 Δk eff . The assessment of uncertainties resulting from the absorbers was predicted to be 0.00029 Δk eff . This results in a total uncertainty of 0.00118 Δk eff . Many of the largest uncertainties, namely the moderator densities, may be reduced with precision dimensional inspection of the components. The 1” HDPE moderators as well as the HDPE reflectors from the original CED-2 were incorporated in the final designs presented here. Additional HDPE moderator plates must be fabricated to complete the configurations. NaCl absorber plates will by fabricated at LLNL. The total additional cost is expected to be $\$$54,250 for the remaining components. Precise inspection, including dimensional, mass, density, and impurity, is recommended for all components. LLNL estimates that these costs are around $\$$12,000. It is expected, based on previous TEX-HEU experiments, that three weeks of experimental facility time is needed to complete the experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

THERMAL MODELING OF HANFORD LEAD CANISTER’S HEATER BENCH TESTS

A computational fluid dynamics (CFD) model was built to simulate planned testing of heater assemblies for the Hanford Lead Canister (HLC) project. The HLC is a canister storage system that will contain heaters to simulate the decay heat of nuclear material and provide the canister storage system with environmental conditions equivalent to the operating conditions on a dry storage pad. The HLC will be equipped with long-term data collection and monitoring systems to provide an early warning of corrosion, pitting, cracking, or other signs of canister degradation that might threaten the integrity of the containment boundary over the potentially long term of dry storage. An important part of the HLC development is to confirm the function and ability of the electric heater assemblies that were specially designed to provide heating similar to the decay heat of nuclear material contained within the canister storage system. Heater bench testing is planned for early 2022 in a test configuration that does not include the canister. The goal of the bench testing is to verify that the heaters can replicate the decay heat of a canister with nuclear material and to validate the thermal models, which are critical to understanding the HLC’s thermal environment, including the local air flow within the canister storage system. Testing of the heater assemblies inside the canister system are planned in the future to validate canister level thermal models, and rigorous pre-deployment testing of the complete HLC cask and canister system is intended to be completed before the HLC is deployed in the 2025-2026 timeframe. This study presents the pre-test temperature predictions of the bench testing. A description of the heater assemblies and planned bench testing is presented. The model was developed with the commercial CFD code STAR-CCM+. An uncertainty analysis was run with the CFD model to determine the uncertainty in the temperature predictions and provide a range over which the predicted temperatures are expected to vary. The uncertainty analysis was performed by coupling STAR-CCM+ with the software Dakota, which provides advanced parametric analyses, including quantification of margins and uncertainty with computational models. This work is expected to provide insight into SNF canister behavior.

Suffield, Sarah R.↗

Reaction Rate Distribution Measurements in ZPPR-15

Small metal foils were irradiated in ZPPR-15 to measure distributions for 239 Pu(n,f), 235 U(n,f), 238 U(n,f), 238 U(n,γ), 115 In(n,n’) and 115 In(n,γ) reaction rates in ZPPR-15. The measured data were reported as basic reaction rate distributions for individual foils, cell-averaged reaction rate distributions and spectral indices or reaction rate ratios. All of the measured data were analyzed and tabulated in such a manner that the data can be readily compared with computed values. The published uncertainties for ZPPR-15 the reaction rate measurements include the statistical uncertainties in the measurements, reproducibility uncertainty and small uncertainties related to corrections for items such as foil compositions and post-irradiation decay between reactor shutdown and foil measurement. There are additional uncertainty components related to detector calibration that are not included in the published uncertainties. A full uncertainty analysis was performed for each foil measurement and for quantities derived from the basic foil measurements. The significant uncertainties were quantified, and a total uncertainty was derived for each measured value. The ZPPR-15 experimental records were used to create detailed as-built Monte Carlo models for the foil reaction rate measurements in ZPPR-15 loadings 123, 134 and 203 and for the fission chamber measurements of 237 Np and 242 Pu fission rates in ZPPR-15 loading 141. Because there were many foils distributed in the ZPPR-15 assembly and because an individual foil occupies a negligible volume in the ZPPR-15 assembly, direct Monte Carlo calculations of the measured foil reaction rates would require an impractical computational effort. The as-built Monte Carlo models and related files were prepared for conversion to deterministic models that can be calculated with realistic computational effort.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Improvement of the $\mathrm{BISON U_3Si_2}$ modeling capabilities based on multiscale developments to modeling fission gas behavior

Uranium silicide (U 3 Si 2 ) is a concept explored as a potential alternative to UO 2 fuel used in light water reactors (LWRs) since it may improve accident tolerance and economics due to its higher thermal conductivity and increased uranium density. U 3 Si 2 has been previously used in research reactors in the form of dispersion fuel, but operated at lower temperatures than commercial LWRs. The research reactor data illustrated that significant gaseous swelling occurs as the fuel burnup increases. Therefore, it is imperative to understand the fission gas behavior of U 3 Si 2 under higher temperature LWR operating conditions. In this work, molecular dynamics and phase-field modeling techniques are used to reduce the uncertainty in select modeling assumptions made in developing the fission gas behavior model for U 3 Si 2 in the BISON fuel performance code. These lower length scale informed models are then utilized in the validation of BISON U 3 Si 2 modeling capabilities to simulate the ATF-1 experiments irradiated in the Advanced Test Reactor (ATR). Sensitivity analysis (SA) and uncertainty quantification (UQ) are included as part of the validation process to identify where further experiments and lower length scale modeling would be beneficial. Here, the multiscale modeling approach utilized in this work can be applied to new fuel concepts being explored for both LWRs and advanced reactors (e.g., uranium nitride, uranium carbide).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

An updated estimate of the Mu2e experiment sensitivity

The Mu2e experiment at Fermilab will search for the conversion of a negative muon into an electron inside the field of a nucleus. This process does not conserve charged-lepton flavour and is heavily suppressed in the Standard Model (SM), with a branching ratio < 10-50. Any evidence of it would be an undeniable evidence of new physics beyond the SM. The project sets out to achieve a single event sensitivity of $\sim 3 × 10^{-17}$ on the ratio between the probability for a conversion of a negative muon into an electron and the one for a muon capture by the nucleus. Such a sensitivity would represent a 4 orders of magnitude improvement on the previous upper limit for the process, making possible to test predictions of different extensions of the SM. Mu2e uses three superconducting solenoids to produce and measure the muon conversions. In the first solenoid, the Production Solenoid, pions and kaons are produced, together with other secondary products, by 8 GeV kinetic energy proton interactions in a tungsten target. A gradient magnetic field is specifically designed to direct low momentum particles into the Transport Solenoid, an S-shaped magnet that filters out particles with unwanted charge and momentum. Muons, produced by pion and kaon decays, finally reach the aluminum Stopping Target in the Detector Solenoid, where they eventually stop and convert. The result of the conversion process is a monochromatic electron of ~105 MeV/c momentum. The Detector Solenoid also hosts the two main detectors: a straw tube tracker and two CsI calorimeter disks, both providing measurement of event kinematics. A germanium detector and a LaBr crystal are located downstream of the Detector Solenoid to measure the X and gamma rays produced by the muon captures in the Stopping Target. A veto system of scintillators covering the Detector Solenoid and half of the Transport Solenoid is used to identify and reject cosmic rays interactions. With respect to the initial project, the Mu2e running plan has evolved to a staged configuration with 2 years at reduced intensity before the 2025 accelerator shutdown for the neutrino beam upgrade and 2 or 3 years at full intensity after that. This, together with geometry changes and a better knowledge of detector performances obtained by the first slice tests, has required a full revision of the signal over background selection that is the subject of this thesis. In order to achieve a new estimate of Mu2e sensitivity, the simulation of the data corresponding to the first 2 years of data acquisition has been performed. This includes both conversion electrons (CE) and the main sources of background: cosmics, decay-in orbits (DIO), radiative pion captures (RPC) and antiprotons. The characteristics of the signal and of the main backgrounds have been studied to define the best selection variables for CE. A special effort has been devoted to the evaluation of the antiproton background. The lack of experimental data for antiproton production cross section makes the systematic uncertainty on this background significant. A new parameterization of the cross section has been developed to fit the existing data and to provide a more reliable estimate of the systematic uncertainty by comparing the results of the old and the new model. A special effort has been devoted to the optimization of the antiprotons Monte Carlo generator. This study has also revealed that the dominant component within this background is represented by antiprotons produced in the opposite direction with respect to the Transport Solenoid entrance and then redirected to it by back-scattering processes in the Production Target. This ultimately highlights the sensitivity of the background estimate to the G4 handling of antiproton interactions in the 1 to 3 GeV/c momentum range. Finally, the final experimental sensitivity has been studied. The momentum and time selection have been optimized to obtain the 5-sigma discovery reach or, in case of no signal, the upper limit on conversion probability. The results confirm that, in the firs t two years of data taking, Mu2e will be able to improve the current experimental sensitivity for muon-to-electron conversion in an atomic field by more than 3 order of magnitudes.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Advanced Reactor Cyber Analysis and Development Environment (ARCADE) for System-Level Design Analysis

Cybersecurity is a persistent concern to the safety and security of Nuclear Power Plants (NPPs), but has lacked data-driven, evidence-based research. Rigorous cybersecurity analysis is critical for the licensing of advanced reactors using a performance-based approach. One tool that enables cybersecurity analysis is modeling and simulation. The nuclear industry makes extensive use of modeling and simulation throughout the decision process but lacks a method to incorporate cybersecurity analysis with existing models. To meet this need, the Advanced Reactor Cyber Analysis and Development Environment (ARCADE) was developed. ARCADE is a suite of publicly available tools that can be used to develop emulations of industrial control system devices and networks and integrate those emulations with physics simulators. This integration of cyber emulations and physics models enables rigorous cyber-physical analysis of cyber-attacks on NPP systems. This report provides an overview of key considerations for using ARCADE with existing physics models and demonstrates ARCADE’s capabilities for cybersecurity analysis. Using a model of the Small Modular Advanced High Temperature Reactor (SmAHTR), ARCADE was able to determine the sensitivity of the primary heat exchangers (PHX) to coordinated cyber-attacks. The analysis determined that while the PHX’s failures cause disruption to the reactor, they did not cause any safety limits to be exceeded because of the plant design, including passive safety features. Further development of ARCADE will enable rigorous, repeatable, and automated cyber-physical analysis of advanced reactor control systems. These efforts will also help reduce regulatory uncertainty by presenting similar types of cybersecurity analyses in a common format, driving standard approaches and reporting.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗