Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Neutron Multiplication”

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 775 records · Page 43

Measurement of the sodium and iodine scintillation quenching factors across multiple NaI(Tl) detectors to identify systematics

The amount of light produced by nuclear recoils in scintillating targets is strongly quenched compared to that produced by electrons. A precise understanding of the quenching factor is particularly interesting for weakly interacting massive particles (WIMP) searches and coherent elastic neutrino-nucleus scattering ( CE ν NS ) measurements since both rely on nuclear recoils, whereas energy calibrations are more readily accessible from electron recoils. There is a wide variation among the current measurements of the quenching factor in sodium iodide (NaI) crystals, especially below 10 keV, the energy region of interest for dark matter and CE ν NS studies. A better understanding of the quenching factor in NaI(Tl) is of particular interest for resolving the decades-old puzzle in the field of dark matter between the null results of most WIMP searches and the claim for dark matter detection by the DAMA/LIBRA collaboration. In this work, we measured sodium and iodine quenching factors for five small NaI(Tl) crystals grown with similar thallium concentrations and growth procedures. Unlike previous experiments, multiple crystals were tested, with measurements made in the same experimental setup to control systematic effects. The quenching factors agree in all crystals we investigated, and both sodium and iodine quenching factors are smaller than those reported by DAMA/LIBRA. The dominant systematic effect was due to the electron equivalent energy calibration originating from the nonproportional behavior of the NaI(Tl) light yield at lower energies, potentially the cause for the discrepancies among the previous measurements. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Exploring entanglement and spectral split correlations in three-flavor collective neutrino oscillations

In environments with prodigious numbers of neutrinos, such as core-collapse supernovae, neutron star mergers, or the early Universe, neutrino-neutrino interactions are dynamically significant. They can dominate neutrino flavor evolution and force it to be nonlinear, causing collective neutrino oscillations. Such collective oscillations have been studied numerically, for systems with up to millions of neutrinos, using mean-field or one-particle effective approximations. However, such a system of interacting neutrinos is a quantum many-body system, wherein quantum correlations could play a significant role in the flavor evolution—thereby motivating the exploration of many-body treatments that follow the time evolution of these correlations. In many-body flavor evolution calculations with two neutrino flavors, the emergence of spectral splits in the neutrino energy distributions has been found to be correlated with the degree of quantum entanglement across the spectrum. In this work, for the first time, we investigate the emergence of spectral splits in the three-flavor many-body collective neutrino oscillations. We find that the emergence of spectral splits resembles the number and location found in the mean-field approximation but not in the width. Moreover, unlike in the two-flavor many-body calculations, we find that additional degrees of freedom make it more difficult to establish a correlation between the location of the spectral splits and the degree of quantum entanglement across the neutrino energy spectrum. The observation from the two-flavor case, that neutrinos nearest to the spectral split frequency exhibit the highest level of entanglement, is more difficult to ascertain in the three-flavor case because of the presence of multiple spectral splits across different pairwise combinations of flavor and/or mass states. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Ternary molecular switching in a single-crystal optical actuator with correlated crystal strain

A growing portfolio of single-crystal optical actuators is forging a new class of photonic materials that hold prospects for quantum technologies. Ruthenium-based complexes that exhibit this phenomenon via SO 2 -linkage photoisomerisation are of particular interest since they display multiple metastable states, once induced by green light; yet, complete photoconversion into each SO 2 -isomeric state is rarely achieved. We discover a new complex, trans-[Ru(SO 2 )(NH 3 ) 4 (4-bromopyridine)]tosylate 2 , that produces 100% photoconverted η 1 -OSO isomeric crystal structures at 90 K, which fully transition into η 2 -(OS)O photoisomers upon warming to 100 K, while the dark-state η 1 -SO 2 structure is wholly recovered by heating the crystal to room temperature. Crystal structures and optical-absorption profiles of each state are captured via in-situ light-induced single-crystal X-ray diffraction and optical-absorption spectroscopy. Results show that both photoisomeric species behave as optical switches, but with distinct optical properties. The photoisomerisation process causes thermally-reversible micro- and nanoscopic crystal strain, as characterised by optical microscopy and in-situ light-induced atomic-force microscopy.

actuators↗

Integrated Steady-State System Package for Nuclear Thermal Propulsion Analysis Using Multi-Dimensional Thermal Hydraulics and Dimensionless Turbopump Treatment

Nuclear thermal propulsion is an evolving technology that can be utilized for long-distance space travel. This technology yields the advantage of a high thrust and specific impulse, but requires an examination of the potential design adjustments necessary to enhance its feasibility. The development of nuclear thermal propulsion requires a comprehensive understanding of the system-level behavior during transient and steady-state operation. This paper extends our previous research by including the proper handling of turbomachinery with multi-channel thermal hydraulic simulations only for steady-state solutions. The system-level approach presented here enables the treatment of the turbopump components through non-dimensional analysis that eliminates the assumption of constant efficiencies. All the other components within the system (e.g., reflector and core) can be discretized to multiple channels and layers, in which the full thermal hydraulic solution is established. The approach chosen here enables the realistic modeling of the propellant flow within the expander cycle by capturing the pressure losses, mass flow rate splits, and enthalpy gain for various operational conditions. The verification of the package is completed through point comparisons of previous investigations into similar system designs. Furthermore, sensitivity studies are used to benchmark the capabilities of the package and investigate solution variations due to the perturbation of operational conditions and regimes. The sensitivity studies performed here are important to capture variation in flow characteristics (e.g., temperature, pressure, mass flow rates) for different design objectives such as the thrust and specific impulse. This work demonstrates that system-level simulations lacking multi-channel capability and proper turbomachinery treatment may yield higher uncertainties in understanding the engine’s response and characteristics to changing various requirements. This is extremely important when screening the design space of such propulsion systems and when transient simulations are required.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NucMesh: nuclear reactor geometry creation and mesh generation module in NEMoSys

NucMesh is a parameterized geometry and mesh generator for nuclear reactors developed within the Nuclear Energy Modeling System NEMoSys at Illinois Rocstar. NEMoSys is a platform developed for mesh generation, adaptive refinement, and solution verification. NucMesh is implemented to be generalized and extensible with a robust computer aided design engine and multiple mesh generation algorithms for unstructured triangular, quad-dominant, and structured quadrilateral meshing. In this paper, we present the geometric and meshing features of NucMesh. Geometrically objects are constructed bottom-to-top and overlaps are addressed automatically. A sophisticated object tracking algorithm prevents data from being lost for segmented objects. We discuss the primitive objects of circle and polygons that constitute the module and show how they are used with example inputs. Arrays of primitives and arrays of arrays are utilized to build large assemblies of objects. The concept of saved objects is discussed to demonstrate how repetitive objects can be reused easily and augmented in place. Three dimensional meshes can be obtained through mesh extrusion where all materials and side sets are extended to three dimensions. We show that side sets can be defined nearly anywhere within the geometry and can then be applied to the mesh. Finally, example reactor meshes are demonstrated for the Idaho National Laboratory Advanced Test Reactor and Los Alamos National Lab Empire reactor, both of which use control drums that NucMesh handles easily. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Estimation of stagnation performance metrics in magnetized liner inertial fusion experiments using Bayesian data assimilation

Here we present a new analysis methodology that allows for the self-consistent integration of multiple diagnostics including nuclear measurements, x-ray imaging, and x-ray power detectors to determine the primary stagnation parameters, such as temperature, pressure, stagnation volume, and mix fraction in magnetized liner inertial fusion (MagLIF) experiments. The analysis uses a simplified model of the stagnation plasma in conjunction with a Bayesian inference framework to determine the most probable configuration that describes the experimental observations while simultaneously revealing the principal uncertainties in the analysis. We validate the approach by using a range of tests including analytic and three-dimensional MHD models. An ensemble of MagLIF experiments is analyzed, and the generalized Lawson criterion χ is estimated for all experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

High-temperature steam oxidation study of irradiated FeCrAl defueled specimens

Post irradiation examinations (PIE) were performed on irradiated iron-chromium-aluminum (FeCrAl) specimens. These FeCrAl specimens were fabricated at the US Department of Energy's Oak Ridge National Laboratory (ORNL). The experimental setup involved subjecting FeCrAl cladding, along with UO 2 pellets, to irradiation in the Idaho National Laboratory Advanced Test Reactor (ATR). In parallel, the FeCrAl alloy tubing without UO 2 pellets was irradiated at ORNL's High Flux Isotope Reactor (HFIR). After irradiation, the ATR-irradiated rodlet was transported to an ORNL hot cell, where it was sectioned into multiple samples for the PIE and severe-accident testing. The sectioning process revealed that the fuel was not bonded to the cladding and could be easily detached from sectioned cladding slices. Microstructural analysis of the fuel cross sections demonstrated no significant interaction between the fuel and the cladding. Additionally, high-temperature steam oxidation tests on defueled cladding segments showed minimal oxygen uptake even at 1200 °C. Here, the ATR-irradiated specimens began to exhibit signs of enhanced oxidation upon reaching a temperature of 1300 °C. Furthermore, enhanced oxidation was observed on the inner surface of the ATR-irradiated FeCrAl specimen, which had been subjected to 1300 °C for a duration of 1 min. By contrast, high-temperature steam oxidation experiments indicated that the HFIR-irradiated FeCrAl cladding provided good thermal stability when exposed to 1300 °C for up to 4 h. Comparative analysis encompassing the oxidation behavior of the ATR-irradiated fueled FeCrAl, HFIR-irradiated unfueled FeCrAl, and unirradiated FeCrAl suggests that the fuel–cladding interaction, although not visible via standard microscale electron microscopy measurements, may accelerate the deterioration of FeCrAl cladding in beyond-design-basis accident scenarios.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Direct observation of C 3 S particle dissolution using fast nano X-ray computed tomography

Tricalcium silicate (C 3 S) occupies 50 % to 70 % of ordinary portland cement (OPC) by mass and it is an important component affecting the hydration of OPC [1], [2], [3], [4], [5], [6], [7]. Generally, the hydration of C 3 S is described by two processes: the dissolution of C 3 S particles and the precipitation of hydration products. While it is understood that the dissolution rates of C 3 S vary with time, more precise measurements are needed to understand this process. Many mechanisms have been proposed to explain the time-evolving dissolution rates of C 3 S [2]. The metastable barrier hypothesis suggests that a thin metastable layer of hydrates forms around the C 3 S particle surface and prohibits the access of grains to the aqueous solution [8], [9], [10], [11], [12]. The slow dissolution step hypothesis suggests that the increased ion concentration from the initial reaction delays the C 3 S dissolution [2], [13], [14], [15], [16], [17]. More recent publications suggest that C 3 S may react differently depending on the existence of crystallographic defects [18], [19], [20]. Etch pits are thought to open on the particle surface during the initial reaction and this contributes to the C 3 S dissolution [21], [22]. As hydrates precipitate and cover these highly reactive surfaces, hydration slows down and the induction period starts [18], [23], [24], [25], [26]. Many experiments have been conducted to test the aforementioned mechanisms. Some hydration studies utilize bulk measurements, such as isothermal calorimetry [27], [28], [29], pore solution analysis [30], quasi-elastic neutron scattering [31], and nuclear magnetic resonance spectroscopy [32], [33]. One limitation of these measurements is that they do not provide direct and detailed information on the individual C 3 S particles. Some other studies utilize imaging techniques, such as scanning electron microscopy (SEM) [34], [35], [36] and transmission electron microscopy (TEM) [37]. However, SEM/TEM cannot track the evolution of individual particles throughout hydration [34], [35], [38], [39], [40] and they do not give insights into the microstructure of materials before hydration [34], [35], [39]. This makes it challenging to draw strong conclusions from only SEM or TEM observations. Synchrotron X-ray tomography techniques have been used more broadly in recent years to study cement hydration. They are not only non-destructive but also able to image a sample in full 3D with resolutions that can reach from micron to nanoscale. Nano computed tomography (nCT) is one technique that has been applied to study cement hydration at the nanoscale [26], [41]. A typical nCT can reach a pixel size from 15 to 65 nm, providing enough detail for observing features <1 μm. However, nCT often takes >0.5 h to finish one scan. This makes the application of this technique on continuous scans for in-situ observations challenging. Fast X-ray computerized tomography (fCT) is another technique that has shown success in studying the time-evolving cement microstructures [20], [42], [43], [44], [45], [46], [47], [48], [49]. Due to the high flux of the X-ray beam from the synchrotron ring, fCT allows a scan to be captured within 1 min at a pixel size of 1 μm. This allows a paste sample to be continuously scanned during the hydration process. However, the micron-sized resolutions limit does not provide detailed insights for particles <5 μm [20], [49]. Fortunately, the combination of nCT and fCT has allowed the development of fast nano X-ray computed tomography (fnCT). fnCT can capture a 3D data set in <2 min at a pixel size of 50 nm. This makes this procedure an exciting method to evaluate hydrating pastes. fnCT collects multiple X-ray radiographs at various rotation angles and generates a 3D model of the scanned sample, which is also referred to as a 3D tomography [50], [51]. In one tomography, the X-ray absorptions of different components (e.g., C 3 S and hydrates) differ as functions of density and chemistry [52], [53]. These X-ray absorption contrasts can be used to extract detailed information about the 3D microstructure [26], [54], [55]. In this paper, fnCT is used to collect time-lapse tomographs of hydrating C 3 S paste from 18 min after mixing to 7 h of hydration. The bulk measurements of anhydrous C 3 S, as well as the microstructural changes of individual C 3 S particles, are directly observed, quantified, and discussed. The dissolution behavior of C 3 S particles at various size scales is systematically analyzed and compared. This work aims to find the relationship between the size of C 3 S particle sizes and their dissolution rates. This provides significant insights into the early-age hydration of C 3 S on length and time scales not previously possible. Because of the magnitude of the data and the substantial amount of observations, this work will solely focus on the change in the anhydrous particles. Changes in the hydration products will be reported in future work.

42 ENGINEERING↗

Low-mode nonuniformity in direct-drive ICF implosions due to laser smoothing techniques employed on OMEGA

For successful laser-direct-drive inertial confinement fusion implosions, the laser irradiation must be highly uniform over the target surface. On OMEGA, multiple laser beams are used to illuminate targets quasi-uniformly. High-mode-number nonuniformities due to laser speckle on each individual beam are reduced by splitting each beam into two orthogonal polarizations (i.e., polarization smoothing, or PS) and a range of wavelengths (i.e., smoothing by spectral dispersion) that are dispersed at the target plane. However, cross-beam energy transfer (CBET) is sensitive to both the polarizations and wavelengths of the interacting beams, so the interplay between CBET and the laser-smoothing schemes results in unique intensity variation across each beam profile, which is a systematic source of low-mode drive nonuniformity on OMEGA. Here, we model these effects and find that the predicted ℓ = 1 mode in the laser-absorption distribution is consistent with the systematic core-flow direction that has been determined from the OMEGA implosion database. We also observe good agreement with the measured core-flow directions for two specific sets of implosions (one with PS, the other without PS) when we also account for the measured beam mispointing and the beam power imbalance.

Crossed beam scattering↗

Micro- and Nano-techniques for the Study of Fission Product Precipitation in SiC Layer

The complexity of fission product distribution and composition within the SiC layer of tristructural isotropic (TRISO) coated particles, irradiation effects on the SiC structure as well as the variable nature of metallic fission product release, has led to the exploration of multiple micro- and nano-characterization techniques. The fine scale of the fission product precipitates necessitates unique learnings and application of electron microscopic techniques with irradiated fuel during the Advanced Gas Reactor (AGR)-1 experiment analysis. A summarized discussion on the specific advanced techniques with associated method development for TRISO coated particles is provided followed by the down selected techniques currently considered to provide the highest impact. Techniques considered are electron probe micro-analyzer (EPMA), scanning transmission electron microscopy (STEM), and Precession electron diffraction (PED). SiC grain boundary characteristics are evaluated due to the mobility of specific fission products through grain boundaries. PED was performed for understanding the role of grain boundary character on fission product transport. One set of comparative analyses between AGR-1 and AGR-2 particles leads to a finding that the AGR-1, high-Ag-retention particle, had statistically more coincidence site lattice (CSL)-related grain boundaries (but fewer low-angle grain boundaries), compared to the AGR-2 particle. This implies that CSL-related grain boundaries may have a direct influence on Ag retention. However, precipitates were found on a relatively small fraction of CSL-related grain boundaries in both particles. In another set of comparison between AGR-1 and AGR-2 particles that were subjected to safety testing, there is significant differences in average SiC grain boundary distributions, but the fission product distributions on these grain boundaries are very similar.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Charge collection and electrical characterization of neutron irradiated silicon pad detectors for the CMS High Granularity Calorimeter

The replacement of the existing endcap calorimeter in the Compact Muon Solenoid (CMS) detector for the high-luminosity LHC (HL-LHC), scheduled for 2027, will be a high granularity calorimeter. It will provide detailed position, energy, and timing information on electromagnetic and hadronic showers in the immense pileup of the HL-LHC. The High Granularity Calorimeter (HGCAL) will use 120-, 200-, and 300-$\mu\textrm{m}$ thick silicon (Si) pad sensors as the main active material and will sustain 1-MeV neutron equivalent fluences up to about $10^{16}~\textrm{n}_\textrm{eq}\textrm{cm}^{-2}$. In order to address the performance degradation of the Si detectors caused by the intense radiation environment, irradiation campaigns of test diode samples from 8-inch and 6-inch wafers were performed in two reactors. Characterization of the electrical and charge collection properties after irradiation involved both bulk polarities for the three sensor thicknesses. Since the Si sensors will be operated at -30 $^\circ$C to reduce increasing bulk leakage current with fluence, the charge collection investigation of 30 irradiated samples was carried out with the infrared-TCT setup at -30 $^\circ$C. TCAD simulation results at the lower fluences are in close agreement with the experimental results and provide predictions of sensor performance for the lower fluence regions not covered by the experimental study. All investigated sensors display 60$\%$ or higher charge collection efficiency at their respective highest lifetime fluences when operated at 800 V, and display above 90$\%$ at the lowest fluence, at 600 V. The collected charge close to the fluence of $10^{16}~\textrm{n}_\textrm{eq}\textrm{cm}^{-2}$ exceeds 1 fC at voltages beyond 800 V.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The impact of detection rate changes and correlations on random-coincidence background measurements

Coincidence detection of multiple particles emitted during an experiment can yield a new depth of understanding of the underlying process under study. However, the probability of detecting particles that are generated from the same physical event within a given coincidence time window is generally much lower than that of detecting particles that appear in the same coincidence time window, but were not created from the same physical event, and are therefore detected randomly in coincidence with each other. Thus, accurate and precise methods of measuring this random-coincidence background are essential for a wide variety of fields of science. A method to determine this background directly using the data themselves without any additional experimental run time or fake signals introduced in the data was recently established (O’Donnell, 2016). This method yields a statistical uncertainty on the random-coincidence background that is orders of magnitude smaller than that of the true coincidence data, though the potential for systematic errors of backgrounds from this method was never explored. In this work, we discuss common varieties of correlated and uncorrelated changes in the detection rates of each particle detected in an experiment. Here we demonstrate here that a correlation between particle detection rates from, for example, an incident particle beam that initiates a physical process of interest, creates systematic errors in the random-coincidence background measurement. We also discuss the impact of a variety of other realistic scenarios for rate changes in experiments. Lastly, a method is introduced to correct for errors in the random-coincidence background from any source, yielding an optimization between statistical precision and eliminating potential lingering systematic errors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Performance of neutron and proton irradiated AC-LGAD sensors

Characterization of strip and pixel AC-LGAD devices with both laser TCT and probe station (IV/CV) will be shown on AC-LGADs irradiated with 1 MeV reactor neutrons at JSI/Ljubljana and with 400 MeV protons at FNAL ITA to fluences from 1e13~$n_{eq}/cm^2$ to a few times 1e15~$n_{eq}/cm^2$. This study was conducted within the scope of the ePIC detector time of flight (TOF) layer R&D program at the EIC, which will feature AC-LGADs with strip and pixel geometry. Sensors in the TOF layer will receive up to 1e13 $n_{eq}/cm^2$ fluence over the lifetime of the experiment.

AC-LGAD strips↗

Improved Cross Section Generation Capability of Griffin in FY22

The Griffin code is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. The code includes a variety of deterministic steady-state transport solvers for fixed source, k-eigenvalue, adjoint, and subcritical multiplication as well as transient solvers for spatial dynamics with the improved quasi-static method. Griffin uses cross section data in the ISOXML format generated from external deterministic or Monte Carlo cross section generation codes. In recent years, the MC 2 -3 modules have been added to Griffin for fast reactor cross section generation, and the self-shielding application programming interface (SSAPI) was implemented in the ISOXML module for thermal reactor cross section generation. The on-the-fly slowing down method and double-heterogeneity treatment have been implemented to SSAPI and verified against particulate fuel-bearing graphite-moderated thermal reactor problems with high accuracy. This year, work has been focused on improving the cross section generation capability of ISOXML and streamlining the cross section generation procedures. In addition, the form function data were added to ISOXML in order to support the pin power reconstruction capability that was newly implemented in Griffin in this fiscal year. To facilitate the cross section generation using MC 2 -3 and SSAPI in Griffin, the cross section generation workflows have been set up for both fast and thermal spectrum reactors. The MOOSE action system tool was devised for fast spectrum problems, and the MOOSE stochastic tool was adopted to the branch calculation procedure for thermal spectrum problems. Meanwhile, to ensure the accuracy of group-constants, the thermal up-scattering kernel calculator accounting for resonance scattering was implemented in ISOXML, demonstrating the accurate computation of a Doppler-broadened scattering kernel of any Legendre order within a reasonable timescale. Other aspects of ISOXML, such as deletion solver and data, documentation, ISOXML file management, and the interface for Mixture, were improved as well.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

AGR-3/4 TRISO Fuel Compact Ceramography

The combined third and fourth irradiation in the Advanced Gas Reactor (AGR) program (AGR-3/4) contained tristructural isotropic (TRISO)-coated particle fuel and designed-to-fail (DTF) fuel particles. The DTF particles were only coated with highly anisotropic pyrocarbon (PyC) so they would purposely fail during the AGR-3/4 irradiation and provide a source of fission products for measurement. To observe the post-irradiation morphology of these DTF particles and the TRISO-coated particles, three AGR-3/4 compacts were mounted in epoxy, sectioned above their centerlines, and ground/polished. Three rounds of grinding/polishing and optical microscopy were performed so that the particles could be observed at multiple planes. Each compact contained approximately 1,872 TRISO particles and exactly 20 DTF particles. The three AGR-3/4 compacts examined covered a wide temperature range and featured both the hottest average irradiation temperature (1375°C) and the coldest/lowest burnup (872°C and 5.5% fissions per initial metal atom) of any compacts to undergo post-irradiation examination (PIE) in the AGR program to date. A total of 29 DTF particles were located and observed via microscopy across the three compacts. All the DTF particles observed via microscopy had completely failed PyC coatings. Some different DTF kernel and PyC morphologies were observed that have been attributed to the differences in irradiation temperature and/or burnup. At low and medium irradiation temperatures (roughly 850 to 1050°C), it appears irradiation-induced dimensional changes in the DTF PyC caused it to fracture and fold in on itself, and the kernel deformed to accommodate this PyC deformation. The extent of DTF kernel deformation to accommodate DTF PyC buckling tended to be greater for the medium-temperature fuel compared to the cold fuel. In the high-irradiation-temperature compact (1375°C), the DTF PyC layer appears to have completely reacted chemically with the kernel material. The only material surrounding the DTF particles in this hot compact resembles the compact matrix material. Observations of the TRISO-coated driver fuel particles and fuel compact graphitic matrix were also made. The TRISO particle kernels in the medium-temperature compact (1047°C), and many in the high-temperature compact (1375°C), showed morphologies consistent with what was seen in AGR-1 and AGR-2. In other TRISO particles from the high-temperature compact, a spatial gradient in the kernel appearance was evident. Here the cool side of the kernel (away from the center of the compact) had a much darker appearance (like that of the buffer and pyrocarbon). This is the first observation of this kind of kernel spatial gradient in AGR fuel. It is believed that the high irradiation temperature of that compact (an average of 1375°C) activates or accelerates an unidentified chemical reaction, and the temperature gradient within the fuel causes a sharp spatial gradient. The low-temperature/low-burnup compact also showed unique kernel morphologies where the TRISO-coated particles had clearly distinguishable oxide and carbide phases at the center of the kernel, an oxide rind surrounding this, and remnants of a carbide skin surrounding the oxide rind. These are features observed in the as-fabricated fuel that are no longer present in higher burnup fuel. Occasional gaps between the outer pyrolytic carbon (OPyC) and graphitic matrix material were observed in all compacts. Gaps were most often found in the small, matrix-filled spaces between adjacent particles. Finally, buffer morphologies in AGR-3/4 TRISO particles were like those observed in AGR-1 and AGR-2, and the buffer fracture frequencies in AGR-3/4 and AGR-2 TRISO particles were plotted versus fast neutron fluence and irradiation temperature. In AGR-3/4 the buffer fracture rate was highest (23%) for an irradiation temperature of 1047°C and a fast fluence of 5.18E25 n/m 2 . Increasing the irradiation temperature to 1375°C reduced the buffer fracture rate to 14%. Reducing the irradiation temperature to 872°C and the fast fluence by a factor of 3 reduced the buffer fracture rate to 7%. The temperature and fluence dependencies observed for AGR-3/4 buffer fracture were consistent with those observed in AGR-2. Lower fluences and/or higher temperatures significantly reduced the buffer fracture frequencies. Lower buffer fracture frequencies were also observed at lower temperatures as long as the fluence was also significantly reduced. This is believed to be due to lower fluence resulting in less irradiation-induced dimensional change (shrinkage) of the buffer, and a higher temperature promoting enhanced creep relaxation of stresses within the buffer, leading to less buffer fracture.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Different to the core: The pre-supernova structures of massive single and binary-stripped stars

The majority of massive stars live in binary or multiple systems and will interact with a companion during their lifetimes, which helps to explain the observed diversity of core-collapse supernovae. Donor stars in binary systems can lose most of their hydrogen-rich envelopes through mass transfer. As a result, not only are the surface properties affected, but so is the core structure. However, most calculations of the core-collapse properties of massive stars rely on single-star models. We present a systematic study of the difference between the pre-supernova structures of single stars and stars of the same initial mass (11–21 M ⊙ ) that have been stripped due to stable post-main-sequence mass transfer at solar metallicity. We present the pre-supernova core composition with novel diagrams that give an intuitive representation of the isotope distribution. As shown in previous studies, at the edge of the carbon-oxygen core, the binary-stripped star models contain an extended gradient of carbon, oxygen, and neon. This layer remains until core collapse and is more extended in mass for higher initial stellar masses. It originates from the receding of the convective helium core during core helium burning in binary-stripped stars, which does not occur in single-star models. We find that this same evolutionary phase leads to systematic differences in the final density and nuclear energy generation profiles. Binary-stripped star models have systematically higher total masses of carbon at the moment of core collapse compared to single-star models, which likely results in systematically different supernova yields. In about half of our models, the silicon-burning and oxygen-rich layers merge after core silicon burning. We discuss the implications of our findings for the “explodability”, supernova observations, and nucleosynthesis of these stars. Our models are publicly available and can be readily used as input for detailed supernova simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Performance Improvements for the Griffin Transport Solvers

Griffin is a Multiphysics Object-Oriented Simulation Environment based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. Griffin includes a variety of deterministic radiation transport solvers for fixed source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. A code assessment performed in FY-20 identified two significant issues with the transport solvers in Griffin: first, the primary heterogeneous SN (discrete ordinates) transport solver based on continuous finite element methods required significant mesh refinement and higher memory usage compared to solvers based on the method of characteristic for equivalent accuracy. Second, the homogeneous PN (spherical harmonics expansion) transport solver did not adequately support polynomial refinement, which is a feature usually required for problems with spatial homogenization and pronounced streaming, typical in fast or gas-cooled reactor systems. To address the first issue, the development effort focused on the more promising discontinuous finite element method (DFEM)-based SN transport solver in Griffin. The addition of an asynchronous parallel transport sweeper and coarse mesh finite difference (CMFD) acceleration have rendered a superior heterogeneous SN transport capability for multiphysics problems that requires far less computing resources in terms of both CPU time and memory usage. This is demonstrated with typical thermal- and fast-spectrum reactor benchmark problems, including 2D Transient Reactor Test, 3D Advanced Burner Test Reactor (ABTR), and 2D and 3D Empire microreactor. For the second issue, the development effort focused on a new transport solver based on the hybrid finite element PN method (HFEM-PN), equivalent to the variational nodal method, as well as a new diffusion solver based on HFEM-Diffusion. This solver is intended for homogenized domains with multiphysics coupling (i.e., supports mesh displacement, seamless temperature feedback, etc.). Initial calculations with the HFEM-Diffusion implementation show very good parallel efficiency for the residual evaluations with the 2D ABTR benchmark. A future development effort will be centered on further improvements to the CMFD, HFEM-PN, and DFEM diffusion solvers to ensure Griffin meets performance and software quality assurance requirements for advanced reactor design and analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗