Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nuclear structure functions”

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,171 records · Page 65

Use of Unmanned Aerial Systems for Post-closure Waste Site Maintenance and Surveillance - 20143

The P- and R-Area Reactor buildings, located at SRS near Aiken, SC, were in-situ decommissioned by grouting below-grade portions of the buildings and demolishing some above-grade structures to grade level. Other reactor building structures were left above-grade and sealed to prevent human or animal access. Because the above grade-structures are expected to continue in their present state for hundreds of years, the condition of the building roofs is critical to mitigate rainwater intrusion. For this reason, building roof areas were strengthened with high strength concrete to ensure the long-term integrity of the roof structure. Periodic inspections of the roof structures are required to ensure that the roofs are functioning properly and to identify damage areas. The traditional method of inspecting closed reactor buildings requires the use of a helicopter and photographer to generate photographs and videos for review. This method proved sufficient but lacked the resolution and clarity that is required for a thorough inspection of the building structure. Fortunately, the Savannah River National Laboratory (SRNL) has established a Small Unmanned Aircraft System (sUAS) program using commercially available and custom-built remote-controlled aircraft. EC and ACP contacted the SRNL program to evaluate whether the sUAS technology could be employed for periodic reactor building inspections. The sUAS can fly within a few meters of the reactor buildings and hover, allowing for a more thorough aerial inspection. The initial sUAS inspection of the P-Area Reactor building was completed in February 2018 and at the R-Area Reactor building in August 2018 and were successful in providing higher resolution photos and videos. The sUAS inspections also revealed that vegetation had begun to grow on the roofs that could potentially damage the structure integrity. SRNL partnered with Virginia Polytechnic Institute and State University to build a custom heavy lift sUAS that could be equipped with herbicides to remotely treat the vegetative growth. Herbicides were successfully applied to the R-Area Reactor building roof using the custom sUAS in September 2018. The use of unmanned aircraft systems at SRS to perform aerial inspections and herbicide treatment in otherwise inaccessible areas has proven to be an efficient and cost-effective technology that provides high value, increases knowledge of facility conditions, and provides for early detection of damage. Periodic inspections of the reactor building roofs using sUAS technology are performed safely, efficiently, and at a significant cost and schedule savings while reducing emissions, noise, and fossil fuel use. The use of sUAS equipment for building inspections and herbicide application at SRS is unique within the DOE complex. The goal of SRS is to apply this remote technology to other waste unit operations and maintenance activities in the future. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Exclusive π 0 p electroproduction off protons in the resonance region at photon virtualities 0.4 GeV 2 ≤ Q 2 ≤ 1 GeV 2

The limited electroproduction process e p → e ' p ' π 0 was measured in the range of photon virtualities Q 2 = 0.4 – 1.0 GeV 2 and the invariant mass range of the p π 0 system of W = 1.1 – 1.8 GeV. These kinematics are covered in exclusive π 0 electroproduction off the proton with nearly complete angular coverage in the p π 0 center-of-mass system and with high statistical accuracy. Nearly 36 000 cross-section points were measured, and the structure functions σ T + ε σ L , σ L T , and σ T T , were extracted via fitting the φ π 0 dependence of the cross section. A Legendre polynomial expansion analysis demonstrates the sensitivity of our data to high-lying N * and Δ * resonances with M > 1.6 GeV. As part of a broad effort to determine the electrocouplings of the N * and Δ * resonances using both single- and double-pion electroproduction, this dataset is key for the reliable extraction of the high-lying resonance electrocouplings from the combined isospin analysis of the N π and π + π - p channels.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Examination of decay heat measurements and their relevance for understanding the origin of the reactor antineutrino anomaly

Measurements of the decay energy released as a function of time following the thermal neutron induced fission of 235 U and 239,241 Pu were performed in the 1970s at Oak Ridge National Laboratory with the purpose of quantifying possible loss of coolant accident scenarios. The derivative of this decay energy with respect to time, known in technical parlance as decay heat, is mainly composed of two terms, that of the electrons produced together with antineutrinos in the $β$-minus decay of the neutron-rich fission products, and that of the $γ$ rays produced in the subsequent decay of excited nuclear levels. In this work we study if this extensive set of decay energy measurements can be used to assess the reactor antineutrino anomaly, that is, the approximately 5% deficit of electron antineutrinos produced by nuclear reactors, first deduced by Mention and collaborators in 2011, and observed by the major reactor antineutrino experiments since. Furthermore, with the assistance of nuclear databases, we are able to obtain the ratio of electron spectra under equilibrium conditions for 235 U to 239 Pu, in better agreement with the lower trend recently reported by Kopeikin and collaborators, as well as those for 235 U to 241 Pu and 241 Pu to 239 Pu, which do not agree well with those measured at the Institut Laue-Langevin in the 1980s. We conclude that a new experimental campaign is needed to measure the electron spectra utilizing a high resolution and signal-to-noise-ratio electron spectrometer and a highly precise fission normalization procedure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modernization of Technical Requirements for Licensing of Advanced Non-Light Water Reactors: Risk-Informed and Performance-Based Evaluation of Defense-in-Depth Adequacy

This document supports the work contained in Nuclear Energy Institute (NEI) 18-04 “Risk-Informed Performance-Based Technology Inclusive Guidance for Advanced Reactor Licensing Basis Development” Revision 0. NEI 18-04 presents a modern, technology-inclusive, risk-informed, and performance-based (TI-RIPB) process for selection of Licensing Basis Events (LBEs); safety classification of structures, systems, and components (SSCs) and associated risk-informed special treatments; and determination of defense-in-depth (DID) adequacy for non-LWRs. The NEI guidance document provides one acceptable means for addressing the aforementioned topics as part of demonstrating a specific design provides reasonable assurance of adequate radiological protection. This report provides the framework and associated methodology guidelines and discussion for establishing, then evaluating, confirming, and documenting the adequacy of defense-in-depth (DID) for advanced non-light-water reactor technologies. It was developed as part of the Licensing Modernization Project led by Southern Company and cost-shared by the United States Department of Energy and has benefited from considerable NRC formal reviews and public workshops. The methodology converts the DID philosophy into a structured process that is implementable, embraces existing United States and international definitions and philosophies of DID that set the foundation for the process. It builds on the DID framework developed in the Department of Energy Next Generation Nuclear Plant Project and earlier works on this subject. The approach to establishing DID adequacy involves the incorporation of DID attributes into the plant capabilities and programmatic elements of DID. The integrated evaluation of DID adequacy includes both quantitative elements to incorporate risk-informed and performance-based (RIPB) considerations and qualitative elements that address uncertainties and limitations in the quantitative models and supporting data. Demonstration of DID adequacy ensures that there are multiple layers of defense for risk-significant challenges to the design and that the plant capabilities and programs that support each layer are provided in a manner that minimizes dependencies among these layers. The focus of this report is assurance of DID adequacy with respect to protection of the public from radiological exposures resulting from accidental releases of radioactive material. While other hazards are not specifically addressed, this methodology is expected to be beneficial for determining DID adequacy for them as well. Risk-informed evaluation of DID considers the integrated performance of all plant SSCs and associated programs to manage daily operational activities, transients, and accidents, including the evaluation of strategies for accident prevention and mitigation. The RIPB LBE scenario methodology used in this evaluation defines the challenges to the plant safety features included in the plant design basis and beyond, and the scope of all deterministic and probabilistic safety evaluations. By examining event sequences across the whole spectrum of LBEs, a systematic assessment of DID can be accomplished.

42 ENGINEERING↗

Development of Analysis Methods that Integrate Numeric and Textual Equipment Reliability Data

Within the Light Water Reactor Sustainability (LWRS) program, the Risk-Informed Systems Analysis (RISA) Pathway is performing collaborative research on the development and deployment of technologies designed to assist operating nuclear power plants (NPPs) to reduce operating costs improve plant reliability and availability. One of the RISA research areas is focusing on the development of methods and tools designed to optimize plant operations (e.g., maintenance/replacement schedules, optimal maintenance postures for plant structures, systems, and components [SSCs]) in a manner that is more cost effective than current approaches and makes better use of available SSC health data. The Risk-Informed Asset Management (RIAM) project targets this research area by creating a direct bridge between component equipment reliability (ER) data and system engineer decision making regarding maintenance activity scheduling and component aging management. In this respect, one challenge that NPP system engineers are facing is that the amount of ER data being continuously generated is not only extremely large in size, but it comes in different forms: textual (e.g., condition or maintenance reports) and numeric (e.g., generated by monitoring systems). All these data elements provide them with valuable insights and information regarding: 1) the discovery of anomalous behaviors or degradation trends, 2) the identification of the possible causes behind such behaviors/trends, and 3) the prediction of their direct consequences. However, several challenges have proved to be roadblocks to this process. While some of these challenges are technical in nature (i.e., data are often distributed over several physical servers/databases), others are conceptual in nature: data elements come in different formats (e.g., numeric or textual), and measured values have different scales (e.g., vibration spectra and oil temperature). The activities performed by the RIAM project during FY23 directly tackles the need to simultaneously integrate the analysis of ER data in all its forms, numeric and textual. Note that such task has never been performed before due to the complexity of the systems under consideration but, most importantly, because of the technical challenges behind the harmonization of ER data formats and the lack of adequate computational methods to analyze them. Our approach borrows ideas and concepts from the medical field where integration of several data sources is vital to assist medical practitioners to perform correct diagnosis and indicate optimal treatments. In our view a NPP asset is equivalent to a patient in a medical context. The main difference is the complexity of a human body is a magnitude more complex when compared to typical assets commonly present in NPPs (e.g., centrifugal pumps, or motor operated valves). This simplifies our first requirement when analyzing heterogenous ER data formats: to put data into “context”. Context is here intended as the additional piece of information that is needed by ER data analysis tools to understand what these data elements are referring to, i.e., which king of knowledge they are generating. In our context, this knowledge can be translated into models that capture the form and functional architecture of assets/systems, their dependencies, and how they interact. These models actually emulate the knowledge that that NPP system engineers possess about assets and systems; this is their key of success when analyzing ER data, their challenge is ability to handle large amount of data. Here, we employ model-based system engineering (MBSE) models of systems and assets to represent and capture their architecture and functional, i.e. cause-effect, relations. Then, ER data elements are processed by identifying first of all which elements of the developed MBSE elements they are referring to. For numeric ER data this task is fairly easy since it is possible to precisely pinpoint what MBSE elements the corresponding sensor are observing (e.g., bearing temperature of a centrifugal pump). Task is much harder for textual data since the information contained in issue or maintenance reports needs to “be understood” by a computational tool. Here we called this process as “knowledge extraction”. Once again, we borrow the experience in the medical field where methods to extract knowledge from textual data have been developed in the past decade. The missing element for us is the availability of a complete dictionary of NPP related concepts (in addition to the MBSE models presented earlier) that can put “text into context”. In FY23, such dictionary has been developed along with all the computational elements required for knowledge extraction. Lastly, once numeric and textual ER data elements have been processed and “understood”, then the last step is the discovery of possible cause-effect relations among them. This is performed by observing if a logical connection through the MBSE models exists, and if the

97 MATHEMATICS AND COMPUTING↗

Study of γ γ → γ ψ ( 2 S ) at Belle

Using 980 fb -1 of data at and around the Υ(nS)(n = 1, 2, 3, 4, 5) resonances collected with the Belle detector at the KEKB asymmetric-energy e+e- collider, the two-photon process γγ → γψ(2S) is studied from the threshold to 4.2 GeV for the first time. Two structures are seen in the invariant mass distribution of γψ(2S): one at M R 1 = 3922.4 ± 6.5 ± 2.0 MeV/c 2 with a width of Γ R 1 = 22 ± 17 ± 4 MeV, and another at M R 2 = 4014.3 ± 4.0 ± 1.5 MeV/c 2 with a width of Γ R 2 = 4 ± 11 ± 6 MeV; the signals are parametrized with the incoherent sum of two Breit- Wigner functions. The first structure is consistent with the X(3915) or the χ c2 (3930), and the local statistical significance is determined to be 3.1σ with the systematic uncertainties included. The second matches none of the known charmonium or charmonium like states, and its global significance is determined to be 2.8σ including the look-elsewhere effect. The production rates are Γ γγ B(R 1 → γψ(2S)) = 9.8 ± 3.6 ± 1.2 eV assuming (J PC , |λ|) = (0 ++ , 0) or 2.0 ± 0.7 ± 0.2 eV with (2 ++ , 2) for the first structure and Γ γγ B(R 2 → γψ(2S)) = 6.2 ± 2.2 ± 0.8 eV with (0 ++ , 0) or 1.2 ± 0.4 ± 0.2 eV with (2 ++ , 2) for the second one. Furthermore, the first errors are statistical and the second systematic, and λ is the helicity.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Czochralski Growth and Characterization of a Compositionally Complex Rare Earth Aluminum Garnet Scintillator: (Gd 1/4 Y 1/4 Tb 1/4 Lu 1/4 ) 3 Al 5 O 12 :Ce

Compositionally complex oxides have garnered increasing interest for their enhanced phase stability and tunable functional properties, yet their development as bulk single crystal scintillators remains limited. Herein, we report the Czochralski growth and characterization of (Gd 1/4 Y 1/4 Tb 1/4 Lu 1/4 ) 3 Al 5 O 12 :Ce (GYTLAG), a compositionally complex garnet incorporating four dodecahedrally coordinated principal rare earth elements. The garnet phase was confirmed by powder and single crystal X-ray diffraction, and macroscopic defects are described. X-ray absorption near-edge structure measurements confirm the 3+ oxidation state of all rare earths and support their occupation of the same crystallographic site; white line intensity variations correlate with the anticipated segregation behavior. Elemental segregation is quantified by SEM/EDS and ICP-OES, and a linear trend was established between the segregation coefficient and the difference between each rare earth’s ionic radius (r) and the average ionic radius (AIR) of the dodecahedral site. This trend offers a predictive framework for compositional control in future REAG crystals grown by the Czochralski method. Photoluminescence and radioluminescence measurements reveal both Ce 3+ and Tb 3+ emission. Scintillation pulses exhibit four-component decay with dominant ~230 µs and ~1.2 ms components, and the light yield is estimated to be up to 43,000 ph/MeV under 137Cs γ-ray excitation. GYTLAG also demonstrates a strong radioluminescence efficiency and 50% lower afterglow at 20 ms compared to a LuAG:Ce reference, underscoring its promise for scintillator applications.

Compositionally complex oxide, high entropy oxide,↗

Using DAPPER to extract the photon strength function of 58 Fe using the inverse Oslo and shape methods

The photon strength function of 58 Fe has been extracted using both the Oslo and Shape methods from particle–γ coincidence data measured using the Detector Array for Photons, Protons, and Exotic Residues, which probes nuclei utilizing (d,p) reactions in inverse kinematics. Four particle–γ coincidence matrices, each constructed with different treatments of the γ–ray energies, are explored in order to observe the impact on the resulting nuclear level density and photon strength. The final photon strength function reported is found to agree well with previous Oslo measurements of other iron isotopes. Systematic uncertainties are included, using different model parameters and their reported errors to perform the Oslo method normalization. The model-independent Shape method is explored and the functional form of the photon strength function obtained is in agreement with the Oslo method results. A low-energy enhancement is not reported for 58 Fe in this work given possible subtraction issues originating from strongly populated states.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Coulomb excitation of 102 Ru with 12 C and 16 O

Here, the Coulomb excitation of 102 Ru was performed with beams of 12 C and 16 O ions. The beam particles scattered at forward angles were momentum analyzed with a magnetic spectrograph. The resolution achieved enabled the populations of the $2^{+}_{1}$ state, the unresolved $2^{+}_{2}$/ $4^{+}_{1}$, and $2^{+}_{4}$/$3^{-}_{1}$, doublets of states, and the $3^{-}_{2}$ state to be determined as a function of the scattering angle. These populations are compared with gosia calculations, yielding B(E2; $2^{+}_{1}$→$0^{+}_{1}$) = 41.5 ± 2.3 W.u., B(E2; $2^{+}_{2}$→$0^{+}_{1}$) = 1.75 ± 0.11 W.u., B(E3; $3^{-}_{1}$→$0^{+}_{1}$) = 31.5 ± 3.5 W.u., and B(E3; $3^{-}_{2}$→$0^{+}_{1}$) = 6.8 ± 0.5 W.u. The B(E3; $3^{-}_{1}$→$0^{+}_{1}$) value is significantly larger than previously measured. The weakly populated $2^{+}_{3}$ state, presumed to be a member of the band built on the $0^{+}_{2}$ state, was observed clearly for a single angle only, and a fit to its population results in B(E2; $2^{+}_{3}$→$0^{+}_{1}$) = 0.053 ± 0.011 W.u. Using the known γ -ray branching ratios for the $2^{+}_{3}$ level, the B(E2; $2^{+}_{3}$→$0^{+}_{2}$) value is calculated to be 18 ± 4 W.u., substantially less than the B(E2; $2^{+}_{1}$→$0^{+}_{1}$) . This suggests that the deformation of the $0^{+}_{2}$ state is lower than that of the $0^{+}_{1}$ state. The results are compared with beyond-mean-field calculations with the Gogny-D1S interaction using the symmetry-conserving configuration-mixing method.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

PORFLOW modeling of Vadose zone flow and transport for the E-area intermediate level vault

In support of the E-Area Performance Assessment, a two-dimensional model of water flow and radionuclide transport through the E-Area Intermediate Level Vault (ILV) and local vadose zone has been developed using the PORFLOWTM software. The purpose of the model is to calculate flux to the water table for radionuclides eluted from the ILV during its operational life, the period of institutional control, and times following site closure. Results of model calculations will be used by a three-dimensional PORFLOW model of transport through the aquifer to determine radionuclide concentrations at a hypothetical 100 meter well and at the site boundary where contaminated groundwater could be assessed by members of the public following site closure. While newly developed, the model structure closely follows that used in the 2008 PA while incorporating a refined computational mesh, updated material properties for the vadose zone soil and vault concrete, revised infiltration rates, and new solid-liquid distribution coefficients. The model also addresses degradation of concrete hydraulic properties by blending soil and concrete water retention curves over a 500 year period. That is, the hydraulic properties of ILV concrete start out as fresh concrete at the time of site closure and degrade to soil properties over a 500 year time period. This approach has not been used previously in E-Area PA’s. The model has been used to calculate water flow through the ILV vadose zone and, on a trial basis, transport of a limited number of radionuclides from the waste region inside the ILV to the water table. Figure ES-1 provides a view of the model and an example of water flow after the concrete has substantially degraded. Figure dimensions are in centimeters. The red area in the left hand side figure is the waste disposal region within the ILV. Colored bands in the figures indicate different soil and concrete regions. Figure ES-2 shows a calculation of flux to the water table (mol/yr) for U-234 and its one-year half-life decay chain (U-234 → Th-230 → Ra-226 → Pb-210). Pb-210 flux is too small to show on the same scale used for Th-230 and Ra-226 in Figure ES-2. The model has also been used to make a preliminary evaluation of the release of tritium from TPBAR disposal containers that will be place in the ILV. A separately calculated tritium release rate as a function of time was used as a source in the PORFLOW model. The TPBAR disposal schedule and tritium release rate differ from the values used in the 2008 PA. Results from these initial trial runs are compared to calculations made for the 2008 PA ILV analysis. In brief, a new PORFLOW model of flow and transport through the ILV vadose zone has been developed and tested. The model is ready to be used for PA calculations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Measurement of the central exclusive production of charged particle pairs in proton-proton collisions at $\sqrt{s}$ = 200 GeV with the STAR detector at RHIC

We report on the measurement of the Central Exclusive Production of charged particle pairs h+h- (h = π, K, p) with the STAR detector at RHIC in proton-proton collisions at $\sqrt{s}$ = 200 GeV. The charged particle pairs produced in the reaction pp → p' + h + h - + p' are reconstructed from the tracks in the central detector and identified using the specific energy loss and the time of flight method, while the forward-scattered protons are measured in the Roman Pot system. Exclusivity of the event is guaranteed by requiring the transverse momentum balance of all four final-state particles. Differential cross sections are measured as functions of observables related to the central hadronic final state and to the forward-scattered protons. They are measured in a fiducial region corresponding to the acceptance of the STAR detector and determined by the central particles’ transverse momenta and pseudorapidities as well as by the forward-scattered protons’ momenta. This fiducial region roughly corresponds to the square of the four-momentum transfers at the proton vertices in the range 0.04 GeV 2 < -t 1 , -t 2 < 0.2 GeV 2 , invariant masses of the charged particle pairs up to a few GeV and pseudorapidities of the centrally-produced hadrons in the range |η| < 0.7. The measured cross sections are compared to phenomenological predictions based on the Double Pomeron Exchange (DPE) model. Structures observed in the mass spectra of π + π - and K + K - pairs are consistent with the DPE model, while angular distributions of pions suggest a dominant spin-0 contribution to π + π - production. For π + π - production, the fiducial cross section is extrapolated to the Lorentz-invariant region, which allows decomposition of the invariant mass spectrum into continuum and resonant contributions. The extrapolated cross section is well described by the continuum production and at least three resonances, the f 0 (980), f 2 (1270) and f 0 (1500), with a possible small contribution from the f 0 (1370). Fits to the extrapolated differential cross section as a function of t 1 and t 2 enable extraction of the exponential slope parameters in several bins of the invariant mass of π + π - pairs. These parameters are sensitive to the size of the interaction region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Tuning Anisotropic Optical Properties of Inorganic and Hybrid Organic–Inorganic MXenes via Topochemical Surface Modification

Surface groups are central to the properties of MXenes, yet their role in optical anisotropy remains largely unexplored. Here, we use a topochemical route to synthesize single crystals of stacked Ti 3 C 2 Cl 2 and hybrid organic–inorganic MXenes (h-MXenes) with lateral sizes of 38–75 μm, rotational registry, and tunable interlayer spacing. Solid-state NMR spectroscopy shows that topochemical substitution generates mixed amido, imido, and hydride surface motifs, which modify the electronic structure of the Ti 3 C 2 inorganic core. Imaging spectroscopic ellipsometry with micron-scale spatial resolution enables reconstruction of the complex dielectric tensor of individual multilayer crystals. Ti 3 C 2 Cl 2 exhibits a type-II hyperbolicity above 930 nm, whereas h-MXenes do not display hyperbolicity within the measured 300–1700 nm window, instead showing reduced in-plane conductivity, suppressed out-of-plane light absorption, and a chain-length-dependent blue shift of a near-infrared absorption feature. These results demonstrate topochemical surface modification as a direct handle for engineering MXenes as surface-programmable optical media.

Hybrid materials↗

Alpha-decay width of a near proton-threshold resonance in the 7 Li( α, α ) channel

We investigate the alpha-decay width of a near proton-threshold resonance in 11 B by the excitation functions of 7 Li(α, α) and 7 Li(α, α’) reactions. This resonance is an example of loosely bound atomic nuclei, understood as small open quantum systems, that exhibit properties significantly influenced by their coupling to the continuum. The present experiment focuses on the formation of a specific state in 11 B at excitation energy in the region near 11.4 MeV by bombarding a 7 Li target with a 4 He beam at energies ranging from 3.92 to 4.56 MeV in the laboratory frame. The study provides data complementary to previous observations of the proton emission in the β-decay of the neutron-rich halo nucleus 11 Be. An R-matrix fit to the data provides resonance energy and partial widths consistent with J π =1/2 + for a narrow near-threshold state in 11 B at 11.400(20) MeV, that is 171(20) keV above the proton emission threshold, with partial α-decay width Γ α =2.5$^{+3.5}_{–1.5}$ keV in elastic scattering and Γ α' =15.8$^{+1.9}_{–0.4}$ keV in inelastic scattering. These findings contribute to the understanding of the structure of this open quantum system resonance.

11B↗

Advanced Reactor Designs Security Analysis, Risk, and Recommendations: Risks, Consequences, and Possible by-Design Mitigation Approaches Associated with Select Advanced Reactors

Next-generation advanced reactors (ARs) incorporate enhanced safety systems, have smaller source terms, and feature compact modular designs, which should lessen their collective risk profiles. However, to fully evaluate risk, security needs to be a part of the equation. Without taking security into consideration, safety systems and components in the new ARs may be vulnerable to sabotage. These base attributes, coupled with enhanced security features specific to AR design through sound engineering and security-by-design (SeBD), should provide developers and operators with lower inherent security risk profiles. Building security early into the AR design may remove or passively secure potential critical targets from an adversary’s reach , thereby increasing overall safety and security. An integrated approach and diverse design team that includes engineering, operations, and security experts are fundamental to building security into the design without sacrificing fundamental operational efficiencies and principles. The objective of this project was to evaluate the security and safety interfaces for five classes of reactors, identify potential security vulnerabilities of structures, systems, and components (SSC), and underscore the need to consider security alongside safety in the design o f these concepts. The five reactor classes evaluated in this project and presented in this report are molten-salt reactors (MSR), high temperature gas reactors (HTGR), sodium-fast reactors (SFR), advanced light-water reactors (ALWR), and microreactors. These designs were selected because they reflect the concepts that are closest to market deployment and have received significant resource investments from the public and private sector. This project assesses the inherent security risks posed by common classes of ARs, provides a methodology and framework to assess security along with safety, and offers an analysis of potential mitigation strategies that could be incorporated. For each AR technology, the SSCs that relate to radionuclide source safety functions are discussed to understand the SSC contribution to safety and relative importance in the protective strategy for the design. The assumptions that went into evaluating each reactor concept originated from generic publicly available nonproprietary information and should not directly be used to qualify an absolute risk profile nor to rank specific AR designs. Instead, the purpose of the analysis is to understand and compare the generic inherent security risks of different AR technologies.

98 - NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL↗

Molten Salt Reactor Experiment Simulation using Shift/Griffin

The Department of Energy (DOE)’s NEAMS focuses its efforts on the development of advanced modeling and simulation (M&S) tools for light-water reactors (LWRs) and non–LWRs (i.e., molten salt reactors, high-temperature gas reactors, microreactors, and fast reactors). In the previous fiscal year, the Multiphysics Applications Driver Technical Area funded molten salt reactor (MSR) M&S at Oak Ridge National Laboratory (ORNL) to generate multigroup macroscopic cross sections with Shift for a MSRE 2D lattice model in Griffin. In addition, Shift’s capability to calculate gamma dose rates from activated components in the primary exchangers in a molten salt breeder reactor was also demonstrated. In fiscal year 2023, ORNL generated multigroup macroscopic cross sections using Shift for a 3D MSRE core model. MSRE depletion calculations using Griffin were also demonstrated in this fiscal year. For the depletion calculation, one-group microscopic cross sections for the 3D MSRE core were generated using Shift, and the decay transmutation library from ORIGEN was converted to an ISOXML file, which is required as input in Griffin. Several Monte Carlo codes, such as OpenMC and Serpent, were also used to benchmark and supplement multigroup cross sections generated by Shift. Multigroup libraries were generated with 8 and 20 group structures, and the study found the 8-group structure to be more accurate when comparing Griffin results to continuous energy (CE) Monte Carlo results. The average flux from CE Shift calculations is up to ~6% higher than the CE Serpent calculations because of different values applied for the energy released per fission (κ values). The average flux in the fuel salt calculated by Griffin using cross sections generated with Shift agrees well with the reference CE Shift solution; the same is valid for the corresponding Serpent results. The maximum relative error is ~6% and ~2% compared to the CE Shift and Serpent reference solutions, respectively. Meanwhile, the average flux calculated by Griffin in the graphite moderator shows a higher difference in the thermal range when compared to both reference Monte Carlo solutions; this result suggests a need for improvement in cross section generation for the graphite moderator in the thermal range in both Monte Carlo codes. Griffin depletion calculations using cross sections from Shift and Serpent were performed and compared against ORIGEN calculations, and the nuclide densities calculated by Griffin were found to be generally in agreement with those of ORIGEN. Because a different approach was taken to calculate the energy released per fission ( κ values) in Shift and Serpent, a difference in nuclide densities from differences in the average flux was observed between Griffin using Serpent and Shift cross sections. Griffin calculations with Shift cross sections produced higher average flux in the salt than with Serpent cross sections, leading to higher consumption of 235 U and higher production of 135 Xe. For time-dependent depletion calculations, cross sections were generated with Serpent, and Griffin’s results using these cross sections were compared to CE Serpent depletion results, demonstrating good agreement. The average difference in keff between Serpent and Griffin as a function of burnup is about 155 pcm. Similarly, good agreement with small differences up to ~0.5% was also noticed in the nuclide density of 235 U and 135 Xe. More details regarding the methodologies invoked to generate the cross section to make code-to-code comparisons are discussed further in this report. User feedback on Griffin and Shift capabilities that will enhance these calculations is provided in this report for future consideration. The work performed this fiscal year can be extended further for multiphysics coupling of Griffin-Pronghorn/SAM with Mole to study precursor flow and salt chemistry.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DIF3D-VARIANT 12.0: Updates and New Features

The DIF3D code has been a workhorse of fast reactor analysis work at Argonne National Laboratory for over 40 years. In 1995, a transport option called VARIANT was added to DIF3D to improve the flux solutions for fast reactor problems which we term DIF3D-VARIANT today. DIF3D-VARIANT performs nodal neutron transport calculations using P N or SP N theory in Cartesian and hexagonal two- and three-dimensional geometries. The limited computing capabilities of the time restricted DIF3D-VARIANT to use at most a 6 th order spatial approximation combined with a P3 flux approximation and P1 scattering kernel for a 33 group structure on most studied reactor problems. Computer capabilities have increased steadily since 1995 and today much larger space-angle-energy approximations are possible. This manuscript serves as an update to the theory section of the original DIF3D-VARIANT manual and details more than twenty years of changes made to DIF3D to make version 12 which was released on November 1 st , 2024. The primary focus of the initial work was to extend the space-angle approximations available in DIF3D-VARIANT such that the error due to transport approximations could be better understood. This work was started and completed in 2002 and marked the official version 10. Unfortunately, those higher order approximations could not be used at that time due to the memory constraints of the BPOINTER part of DIF3D (limited to 2 GB). In version 11, completed in 2012, BPOINTER was circumvented in DIF3D-VARIANT for the largest arrays by introducing a Fortran 90 module called LMA (Large Memory Array). This seamlessly replaces all of the functionality of the BPOINTER concept, but it allows 64 bit addressing for every array such that they can be larger than 2 GB. It is now common for DIF3D-VARIANT jobs to consume 50 GB of memory on modern workstations when using high order space-angle approximations and a large number of groups. Many improvements were made to version 11 from 2012 to 2022 when work to create version 12 started. For version 12, several parts of DIF3D were updated to improve performance and thread parallelism was introduced to further reduce the runtime. Numerous minor bugs were discovered in DIF3D-VARIANT as part of the process of creating the perturbation and sensitivity code PERSENT. All of these algorithmic problems were identified in the transition from version 10 to version 11 which prevented DIF3D-VARIANT from running efficiently and reliably. Firstly, the coarse mesh rebalance scheme would routinely diverge and a study detailed in this report demonstrates how it was also typically not effective. This is not a failure of the coarse mesh rebalance methodology, but a failure of its implementation in DIF3D-VARIANT for hexagonal geometries. The fission source extrapolation algorithm was also found to be unreliable on larger group structure problems, leading to divergence in some cases and a negligible improvement in performance overall. Finally, the “Omega” acceleration applied to the partial current solver routine of DIF3D-VARIANT was found to cause DIF3D-VARIANT to converge to the wrong answer. To resolve these issues, both the coarse mesh rebalance and fission source extrapolation were permanently disabled in version 11. The Tchebychev acceleration was put in as a temporary reliable alternative but it is generally inferior to coarse mesh rebalance or coarse mesh finite difference. For the Omega acceleration, the factor was restricted to guarantee that it would not cause follow-on errors in PERSENT. Due to limited funding to support maintenance and development of DIF3D in the last 10 years, no effort was spent since to resolve the outer iteration acceleration. Except for the threading work, all of the changes discussed in this manuscript refer to changes made between version 10 and version 11. Performance comparisons are done to demonstrate the improvements from version 9 to version 12. As will be demonstrated, the updated versi

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Projections of two-particle correlations onto transverse rapidity in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV at STAR

Two-particle correlation measurements projected onto two-dimensional, transverse rapidity coordinates (y T1 ,y T2 ), allow access to dynamical properties of the QCD medium produced in relativistic heavy-ion collisions that angular correlation measurements are not sensitive to. Here, we report non-identified charged-particle correlations for Au + Au minimum-bias collisions at $\sqrt{s_{NN}}$ = 200 GeV taken by the STAR experiment at the Relativistic Heavy-Ion Collider (RHIC). Correlations are presented as 2D functions of transverse rapidity for like-sign, unlike-sign and all charged-particle pairs, as well as for particle pairs whose relative azimuthal angles lie on the near-side, the away-side, or at all relative azimuth. The correlations are constructed using charged particles with transverse momentum p T ≥ 0.15 GeV/c, pseudorapidity from –1 to 1, and azimuthal angles from –π to π. The significant correlation structures that are observed evolve smoothly with collision centrality. The major correlation features include a saddle shape plus a broad peak with maximum near y T ≈ 3, corresponding to p T ≈ 1.5 GeV/c. The broad peak is observed in both like- and unlike-sign charge combinations and in near- and away-side relative azimuthal angles. The all-charge, all-azimuth correlation measurements are compared with the theoretical predictions of HIJING and EPOS. The results indicate that the correlations for peripheral to mid-central collisions can be approximately described as a superposition of nucleon + nucleon collisions with minimal effects from the QCD medium. Strong medium effects are indicated in mid- to most-central collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Predicting Saturated Hydraulic Conductivity over Time for Degrading Saltstone Vault Concrete - 20376

At the Savannah River Site (SRS), low-level liquid tank waste is mixed with cementitious materials to produce a solid, monolithic waste form referred to as saltstone. The liquid saltstone mixture is pumped into massive concrete structures referred to as saltstone disposal units (SDUs) for curing and long-term storage. SDUs are vaults made of high performance concrete forming a barrier between saltstone and the surrounding environment. Predicting the long-term performance of SDU vault concrete, as well as the saltstone material itself, is critical to understanding the potential releases of radionuclides to the surrounding environment as described in the Saltstone Disposal Facility (SDF) performance assessment (PA). The hydraulic conductivity profile for SDU vault concrete was conservatively assumed to degrade linearly from the initial saturated hydraulic conductivity of the intact concrete over a long period of time. The initial hydraulic properties of saturated SDU vault concrete have been measured, and the final hydraulic properties of degraded vault concrete are assumed to be the same as the porous backfill material that will surround the vault after closure. The assumed final condition of the vault concrete is also conservative in that it provides a greater estimate of water penetration through the vault than would generally be expected. Therefore, insights that provide more realistic but conservative estimates as to how degradation processes may change the hydraulic properties of the SDU vault concrete, the time dependent relationship of the vault concrete hydraulic properties to degradation (and thus time), and the final hydraulic properties of the vault concrete provides useful additions to the SDF PA. This paper focuses on the assumption that the saltstone vault concrete saturated hydraulic conductivity would decrease linearly as a function of degradation time; this assumption resulted from a worst-case, weighted (by relative layer thickness), arithmetic average of intact and backfill material hydraulic conductivities representing flow parallel to a layered system. This research indicates that the effective saturated hydraulic conductivity of the degraded SDU vault concrete should be based on the relative thickness-weighted geometric mean of the intact and degraded layers hydraulic conductivities. This alternative, which is supported extensively in the relevant literature and evaluated in this paper, leads to a more realistic and defensible estimate of the effective saturated hydraulic conductivity in degraded vault concrete while remaining conservative. The impact of using the more realistic estimate of the effective saturated hydraulic conductivity for the degraded concrete decreases the effective saturated hydraulic conductivity by up to several orders of magnitude for early times (i.e., when the SDU vault concrete is assumed degraded via sulfate attack and carbonation). This revised estimate of the effective saturated hydraulic conductivity has been implemented in the recently revised SDF PA. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗