Search NASA⌕ Search

SEARCH · Search NASA

Results for “Decay rate”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Water-Triggered Photoinduced Electron Transfer in Acetonitrile–Water Binary Solvent. Solvent Microstructure-Tuned Reactivity of Hydrophobic Solutes

The solvent-composition dependence of quenching triplet states of benzophenone (3 BP) by anisole in acetonitrile-water (ACN-H2O) mixtures was investigated by laser flash photolysis over the water mole fraction (xw) increasing from 0 to 0.92. Single exponential decay of 3 BP was observed over the whole composition range. The quenching rate constant consistently increased with the water content but increased far more rapidly with xw > 0.7. The water-triggered electron-transfer (ET) mechanism was confirmed by a steeply growing quantum yield of the benzophenone ketyl radical anion, escaping back-ET when the partial water volume exceeded the acetonitrile one. The water-content influence on the 3 BP quenching rate was described by a kinetic model accounting for the microheterogeneous structure of the ACN-H2O mixtures and the very different solubility of the reactants in the solvent components. According to the model, the ET mechanism occurs at a rate constant of 1.46 × 109 M-1 s -1 and is presumably assisted by the ACN-H2O hydrogen-bonding interaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Antineutrino energy spectrum unfolding based on the Daya Bay measurement and its applications

The prediction of reactor antineutrino spectra will play a crucial role as reactor experiments enter the precision era. The positron energy spectrum of 3.5 million antineutrino inverse beta decay reactions observed by the Daya Bay experiment, in combination with the fission rates of fissile isotopes in the reactor, is used to extract the positron energy spectra resulting from the fission of specific isotopes. This information can be used to produce a precise, data-based prediction of the antineutrino energy spectrum in other reactor antineutrino experiments with different fission fractions than Daya Bay. The positron energy spectra are unfolded to obtain the antineutrino energy spectra by removing the contribution from detector response with the Wiener-SVD unfolding method. Consistent results are obtained with other unfolding methods. A technique to construct a data-based prediction of the reactor antineutrino energy spectrum is proposed and investigated. Given the reactor fission fractions, the technique can predict the energy spectrum to a 2% precision. In addition, we illustrate how to perform a rigorous comparison between the unfolded antineutrino spectrum and a theoretical model prediction that avoids the input model bias of the unfolding method.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Lagrangian statistics of a shock-driven turbulent dynamo in decaying turbulence

Small-scale fluctuating magnetic fields of order nG are observed in supernova shocks and galaxy clusters, where its amplification is likely caused by the Biermann battery mechanism. However, these fields cannot be amplified further without the turbulent dynamo, which generates magnetic energy through the stretch-twist-fold (STF) mechanism. Thus, we present here novel 3D magnetohydrodynamic (MHD) simulations of a laser-driven shock propagating into a stratified, multiphase medium, to investigate the post-shock turbulent magnetic field amplification via the turbulent dynamo. The configuration used here is currently being tested in the shock tunnel at the National Ignition Facility (NIF). In order to probe the statistical properties of the post-shock turbulent region, we use 384 × 512 × 384 tracers to track its evolution through the Lagrangian framework, thus providing a high-fidelity analysis of the shocked medium. Our simulations indicate that the growth of the magnetic field, which accompanies the near-Saffman kinetic energy decay (E kin ∝ t -1.15 ) without turbulence driving, exhibits slightly different characteristics as compared to periodic box simulations. Seemingly no distinct phases exist in its evolution because the shock passage and time to observe the magnetic field amplification during the turbulence decay are very short (~0.3 of a turbulent turnover time). Yet, the growth rate is still consistent with those expected for compressive (curl-free) turbulence driving in subsonic, compressible turbulence. Phenomenological understanding of the dynamics of the magnetic and velocity fields are also elucidated via Lagrangian frequency spectra, which are consistent with the expected inertial range scalings in the Eulerian–Lagrangian bridge.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Classification of Dissolution Events Using Fusion of Effluents Measurements and Classifiers

Classifiers for dissolution events at a radiochemical processing facility are studied using gamma spectra measurements of effluents collected by a high purity germanium detector located at its off-gas stack. Data sets collected at the Oak Ridge National Laboratory’s Radiochemical Engineering Development Center under a Pu dissolution campaign spanning a three months period are utilized. Features corresponding to the activity levels of 15 radionuclides, including isotopes of iodine, krypton, and xenon, that are indicated by the target decay chains, are computed from the spectra at 1 hour intervals. A conceptualization diagram is developed to reflect the steps from the source to measurement to feature computation that depend on fission products indicated by decay chains, chemical processing, and effluents transport to the off-gas stack. A diverse set of eight classifiers based on different design principles are trained using the ground truth data for this campaign, and the outputs of top three classifiers, namely, classification trees, Ensemble of Trees (EOT), and k-nearest neighbor, are combined using EOT classifier-fuser. Our results show that for 5-fold cross validation, features associated with isotopes of xenon provide the lowest classification error among the different elements across the classifiers; the classification error is furthered improved when all 15 isotope features are used by each classifier, and it is again improved by the fusion of three classifiers. Further reduction in classification error is achieved by using a measurement window of 1-2 days which is identified based on half-life time estimates of the isotopes; it is long enough for the stabilization of feature estimates while being short enough not to be affected by the follow on dissolution events. As a net result of feature and classifier fusion, combined with the incorporation of decay chain and isotope half-life information, this approach achieves 98% detection rate while maintaining a false alarm rate under 2% for this data set.

Rao, Nageswara↗
DOE OSTI2022

SCALE Analysis of a Fluoride Salt-Cooled High-Temperature Reactor in Support of Severe Accident Analysis

As part of a US Nuclear Regulatory Commission–sponsored project to assess the modeling and simulation capabilities for accident progression, source term, and consequence analysis for advanced reactor technologies with SCALE and MELCOR, SCALE was used for the modeling and simulation of a fluoride salt-cooled high-temperature reactor (FHR). Based on the preconceptual design for a small modular 236 MWth FHR developed by the University of California, Berkeley (PB-FHR-Mk1), a SCALE model of the PB-FHR-Mk1 reactor core was developed. The reactor was modeled at equilibrium state with different fuel compositions in different regions of the reactor. An iterative approach was used to interpolate and mix the burnup-dependent fuel compositions obtained through the depletion calculation of a core slice model. After demonstrating the applicability of SCALE’s multigroup (MG) approach for the simulation of the PB-FHR-Mk1, the resulting equilibrium core was studied in terms of the power profile, the flux profile, temperature reactivity coefficients, and the xenon reactivity. Furthermore, the tritium production rate in the salt coolant was determined, and the dependence of the one-group cross sections on the burnup and location in the reactor core was studied. The results obtained with SCALE were post-processed to provide the MELCOR team with the core inventory and decay heat of the equilibrium core, a zone-wise power profile, temperature feedback coefficients, the tritium production rate, and the xenon worth.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Measurement of 𝐵 meson production fraction ratios in proton-proton collisions at √𝑠 = 13 TeV using open-charm and charmonium decays

Production fraction ratios of 𝐵 + , 𝐵 0 , and 𝐵$^{0}_{𝑠}$ mesons are measured in proton-proton collisions at √𝑠 = 13 TeV using a special dataset recorded in 2018 with high-rate triggers designed to collect an unbiased sample of 1⁢010 𝑏 hadrons with the CMS experiment at the LHC. These data allow the study of the open-charm decays of 𝐵 mesons (𝐵 (𝑠) → 𝜋⁢𝐷 (𝑠) ) where the 𝐷 meson decays into fully hadronic final states. By utilizing known branching fractions and precise theoretical calculations, production fraction ratios as functions of 𝐵 meson transverse momentum (𝑝 T ) and rapidity (𝑦) are measured using open-charm decays in the kinematic range of 8 < 𝑝 T < 60 GeV and |𝑦| <2.25. In addition, the same dataset is used to measure the relative production fraction ratios with the charmonium decay channels [𝐵 (𝑠) →𝑋⁢𝐽/𝜓 with 𝑋 indicating a 𝐾 + , 𝐾*(892) 0 , or 𝜙⁡(1020) meson] where the 𝐽/𝜓 meson decays into a pair of muons. The open-charm results are used to normalize the relative production fraction ratios obtained from the charmonium samples. Measurements of the ratios of branching fractions of 𝐵 meson decays to charmonium and open-charm final states are also reported, which will improve the world-average values of these ratios. Finally, we test isospin invariance in 𝐵 meson production in proton-proton collisions and observe that it holds within the experimental precision.

Bottom mesons↗
DOE OSTI2024

Radiative decay branching ratio of the Hoyle state

The triple-alpha process is a vital reaction in nuclear astrophysics, characterized by two consecutive reactions [2α ⇌ 8 Be (α, γ) 12 C] that drive carbon formation. The second reaction occurs through the Hoyle state, a 7.65 MeV excited state in 12 C with J π = 0 + . The rate of the process depends on the radiative width, which can be determined by measuring the branching ratio for electromagnetic decay. Recent measurements by Kibédi et al. conflicted with the adopted value and resulted in a significant increase of nearly 50% in this branching ratio, directly affecting the triple-alpha reaction. Here, this work aims to utilize charged-particle spectroscopy with magnetic selection as a means to accurately measure the total radiative branching ratio (Γ rad /Γ) of the Hoyle state in 12 C.

6 ≤ A ≤ 19↗
DOE OSTI2026

Axionlike particles and sterile neutrinos solve the 𝐵 → 𝐾⁢$𝜈⁢\bar{𝜈}$ and 𝐵 → 𝜋⁢𝐾 puzzles

The recent measurement of the branching ratio of 𝐵 + → 𝐾 + +inv (where “inv” denotes invisible states) by the Belle II Collaboration is enhanced relative to the standard model expectation by 2.7⁢𝜎. An older puzzle persists in measurements of the branching ratios and 𝐶⁢𝑃 asymmetries of 𝐵 → 𝜋⁢𝐾 decays. We address these two anomalies in flavor-changing neutral current 𝐵 decays, with a short-lived axionlike particle (ALP) with mass close to that of the 𝜋 0 . In the model with the minimum number of new couplings, the ALP has couplings to the photon, top quark and a heavy sterile neutrino. The ALP contributes to the 𝐵 → 𝜋 0 ⁢𝐾 decays by mixing with the 𝜋 0 . It contributes to 𝐵 + → 𝐾 + + inv by its off-shell coupling to sterile neutrino pairs. The model can explain the excess in the total rate, but not the observed distribution of signal events. We make predictions for all 𝐵 → 𝐾 (*) + inv modes and for the rare kaon decays, 𝐾 + → 𝜋 + + inv and 𝐾 𝐿 → 𝜋 0 + inv. We find an appreciable contribution to the magnetic moment of the muon, and negligible contributions to the magnetic moment of the electron and 𝑏 → 𝑠⁢𝑒 + ⁢𝑒 − .

axion-like particles↗
DOE OSTI2023

Observation of the rare decay of the $\eta$ meson to four muons

A search for the rare $\eta$$\to$$\mu^+\mu^-\mu^+\mu^-$ double-Dalitz decay is performed using a sample of proton-proton collisions, collected by the CMS experiment at the CERN LHC with high-rate muon triggers in 2017-2018 and corresponding to an integrated luminosity of 101 fb$^{-1}$. A signal having a statistical significance well in excess of 5 standard deviations is observed. Using the \emm decay as normalization, the branching fraction $\mathcal{B}($η$ \to \mu^+\mu^-\mu^+\mu^-)$ = ( 5.0 $\pm$ 0.8 (stat) $\pm$ 0.7 (syst) $\pm$ 0.7 $\mathcal{B}_{2\mu}$ ) $\times$ 10$^{-9}$ is measured, where the last term is the uncertainty in the normalization channel branching fraction. This is the first measurement of this branching fraction and is found to be in agreement with theoretical predictions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Novel estimation of aerosol processes with particle size distribution measurements: a case study with the TOMAS algorithm v1.0.0

Atmospheric aerosol microphysical processes are a significant source of uncertainty in predicting climate change. Specifically, aerosol nucleation, emissions, and growth rates, which are simulated in chemical transport models to predict the particle size distribution, are not understood well. However, long-term size distribution measurements made at several ground-based sites across Europe implicitly contain information about the processes that created those size distributions. This work aims to extract that information by developing and applying an inverse technique to constrain aerosol emissions as well as nucleation and growth rates based on hourly size distribution measurements. We developed an inverse method based upon process control theory into an online estimation technique to scale aerosol nucleation, emissions, and growth so that the model–measurement bias in three measured aerosol properties exponentially decays. The properties, which are calculated from the measured and predicted size distributions, used to constrain aerosol nucleation, emission, and growth rates are the number of particles with a diameter between 3 and 6 nm, the number with a diameter greater than 10 nm, and the total dry volume of aerosol ( N 3–6 , N 10 , V dry ), respectively. In this paper, we focus on developing and applying the estimation methodology in a zero-dimensional “box” model as a proof of concept before applying it to a three-dimensional simulation in subsequent work. The methodology is first tested on a dataset of synthetic and perfect measurements that span diverse environments in which the true particle emissions, growth, and nucleation rates are known. The inverse technique accurately estimates the aerosol microphysical process rates with an average and maximum error of 2 % and 13 %, respectively. Next, we investigate the effect that measurement noise has on the estimated rates. The method is robust to typical instrument noise in the aerosol properties as there is a negligible increase in the bias of the estimated process rates. Finally, the methodology is applied to long-term datasets of in situ size distribution measurements in western Europe from May 2006 through June 2007. At Melpitz, Germany, and Hyytiälä, Finland, the average diurnal profiles of estimated 3 nm particle formation rates are reasonable, having peaks near noon local time with average peak values of 1 and 0.15 cm –3 s –1 , respectively. The normalized absolute error in estimated N 3–6 , N 10 , and V dry at three European measurement sites is less than 15 %, showing that the estimation framework developed here has potential to decrease model–measurement bias while constraining uncertain aerosol microphysical processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Differential Organic Carbon Mineralization Responses to Soil Moisture in Three Different Soil Orders Under Mixed Forested System: Supporting Data

This data contains data from 90-day long incubation study which aimed to look at the soil moisture-texture relationship on soil organic carbon (SOC) cycling. Soils were collected from three distinct soil textures from mixed forests in 2017: sandy (Georgia, 2017-05-01), loamy (Missouri, 2017-06-14) and clayey (Texas, December 2017) were incubated at different soil moisture levels (air-dried, 25% water holding capacity (WHC), 50% WHC, 100% WHC and 175% WHC) at room temperature for a period of 90 days. Files contain microbial respiration, active and slow SOC pools, and their respective mineralization rates, extractable organic carbon (C), and C-acquiring extracellular enzymes. Findings from these data were used in Singh et al. (2021). This study aimed to examine the interactive effect of soil moisture and texture on SOC mineralization. Soil samples of three distinct textures (sandy, loamy, and clayey) were collected from mixed forests of Georgia, Missouri, and Texas, respectively. Soil cores of 5 cm diameter were collected from numerous random locations at each site from 0-15 cm depth after scraping the litter layer and mixed thoroughly to obtain a composite sample per site. Three additional soil cores were collected to determine the WHC using pressure plate extractors. Soil samples were composited, and triplicate soil samples were incubated in mason jars for a period of 90 days at room temperature under different moisture regimes: air dried, 25% WHC, 50% WHC, at WHC and 100% saturation. Soil respiration was measured weekly, and destructive sampling was conducted at 1, 15, 60, and 90 days to determine extractable organic C, C acquiring enzyme activity, and active and slow SOC pools with their respective mineralization rates. The C acquiring enzyme activity was the total activity of α-glucosidase, β-glucosidase, cellobiohydrolase, and β-xylosidase enzymes. Gas samples for microbial respiration measurements were collected from headspace of incubation jars through the sampling ports on the lids and then analyzed using a Shimadzu Gas Chromatograph (GC-2014). Prior to sampling, the vials were evacuated. Blank correction was also done by collecting gas samples from empty incubation jars. Double pool exponential decay model was used in SigmaPlot to determine the active and slow SOC pools and their mineralization rates (Farrar et al., 2012; Jagadamma et al., 2014). The C-acquiring extracellular enzymes were measured using the microplate method by German et al., (2011). Microbial community structure was determined using the phospholipid fatty acid (PLFA) and neutral lipid fatty acid (NLFA) analyses (Buyer and Sasser, 2012). This dataset has seven data files provided in comma-separate (*.csv) format. Additional metadata are provided: seven data dictionaries and a file-level metadata file in comma separate (*.csv) format and a user guide in PDF (*.pdf) format.

Carbon acquiring enzyme activity↗
DOE OSTI2023

Stability study of aqueous foams under high-temperature and high-pressure conditions relevant to Enhanced Geothermal Systems (EGS)

In Enhanced Geothermal System (EGS) hydraulic fracturing is carried out by injecting cold water into deep Hot Dry Rocks (HDR) under carefully controlled conditions to create new or reopen existing fractures. Water-based fracturing fluids demonstrate some challenges including immense quantity of water usage, water sensitivity of the formations, water blocking, and lack of proppant carrying capacity and transportation. Thus, an alternative is to use foam-based fracturing fluid which offers potential advantage over conventional water-based fracturing fluid such as minimum water usage, reduced wellbore damage, high proppant carrying capacity, and less environmental damage. However, foams are complex mixture of gaseous phase and liquid phase which are thermodynamically unstable at downhole conditions, and their stability can decrease over time due to liquid drainage, bubble coarsening, and coalescence. This paper shows laboratory experiments executed to study foam stability at high temperature (200°C) and high pressure (6.9MPa) conditions which simulates the geothermal environment. Foam stability was characterized by half-life of foam, which is defined as the time taken by the foam to decreases by 50% of its original height due to drainage. In this paper, two types of gaseous phases, nitrogen (N 2 ) and carbon dioxide (CO 2 ) were investigated. Also, based on successful practice of foam-based fracturing fluid in oil and gas industries, four surfactants, including Alpha olefin sulfonate (AOS), Sodium dodecyl sulfonate (SDS), Tergitol™ (NP – 40), and Cetyltrimethylammonium chloride (CTAC) at optimum concentration of 1wt.% were tested for best stability performance. In addition, different stabilizing agents including guar gum, bentonite clay, crosslinker, silicon dioxide nanoparticles (SiO 2 ), graphene oxide (GO) were also studied. Experimental results showed that N 2 foams were more stable than CO 2 foams. It was observed that foam half-life decreased with the increase in temperature. Among all the surfactants, AOS foams showed the most promising thermal stability at high temperatures. Moreover, with the addition of stabilizing agents, foam's half-life was enhanced. Stabilizing agents such as crosslinker and GO dispersion showed the most stable foams with half-life recorded at 20min and 17min, respectively, at 200°C and 6.9MPa. Finally, pressure also showed a positive effect on foam stability; with increased pressure, foam half-life was increased. Based on the experimental data, analytical models for the effect of temperature and pressure were developed, considering foam degradation is a first-order kinetic reaction that linearly depends on the foam drainage mechanism. The effect of temperature on foam half-life was studied as an exponential decay model. In this model, foam half-life is a function of drainage rate constant (D A ) and activation energy (E a ) of the foam system. The effect of pressure on foam half-life was found to obey a power-law model where an increase in pressure showed an increase in foam half-life. Furthermore, a linear relation was studied for the effect of pressure on foam activation energy and drainage rate. Then the combined effects of temperature and pressure were studied, which yielded an analytical model to predict the foam stabilities in terms of half-life for different foam compositions. In conclusion, this research indicates that with an appropriate selection of surfactants and stabilizing agents, it is possible to obtain stable foams, which could replace conventional water fracturing fluid under EGS conditions.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Mismatching integration-enabled strains and defects engineering in LDH microstructure for high-rate and long-life charge storage

Layered double hydroxides (LDH) have been extensively investigated for charge storage, however, their development is hampered by the sluggish reaction dynamics. Herein, triggered by mismatching integration of Mn sites, we configured wrinkled Mn/NiCo-LDH with strains and defects, where promoted mass & charge transport behaviors were realized. The well-tailored Mn/NiCo-LDH displays a capacity up to 518 C g -1 (1 A g -1 ), a remarkable rate performance (78%@100 A g -1 ) and a long cycle life (without capacity decay after 10,000 cycles). We clarified that the moderate electron transfer between the released Mn species and Co 2+ serves as the pre-step, while the compressive strain induces structural deformation with promoted reaction dynamics. Theoretical and operando investigations further demonstrate that the Mn sites boost ion adsorption/transport and electron transfer, and the Mn-induced effect remains active after multiple charge/discharge processes. This contribution provides some insights for controllable structure design and modulation toward high-efficient energy storage.

25 ENERGY STORAGE↗
DOE OSTI2024

Growth of Native Oxides of Niobium Thin Films

Both XPS and ToF-SIMS data and analysis suggest that Nb oxides grow rapidly and overwhelm Nb metal within the first 24 hours following a logarithmic growth. Nb2O5 growth is dominant out of all types of Nb oxide with reaching 80% At conc within 80 days. NbO2 is quickly saturated, maintaining 10% At conc over time with a logarithmic growth. NbO on the opposite follows an exponential decay. Oxide thickness growth also implies sample preparation techniques impact oxide saturation and growth rate.

Lei, Zhicheng↗
DOE OSTI2025

Radiative processes on a quantum computer

Radiative processes, where a photon/neutrino is emitted as a result of a collision or decay of a particle, play a central role in atomic, nuclear, and particle physics. Their rate is determined by certain off-diagonal matrix elements between different initial and final states. We propose a method to compute them using quantum computers. It relies on a single extra qubit that, in a certain sense, represents the photon/neutrino. The generic formula relating this matrix element to the amplitude and frequency of oscillations of the extra qubit is derived for the near-resonance case. Furthermore, we demonstrate the feasibility of the method by using it in actual quantum computations and simulations of simple systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Observation of the Rare Decay of the 𝜂 Meson to Four Muons

A search for the rare η → μ + μ − μ + μ − double-Dalitz decay is performed using a sample of proton-proton collisions, collected by the CMS experiment at the CERN LHC with high-rate muon triggers in 2017-2018 and corresponding to an integrated luminosity of 101 fb −1 . A signal having a statistical significance well in excess of 5 standard deviations is observed. Using the η → μ + μ − decay as normalization, the branching fraction B(η → μ + μ − μ + μ − ) = [5.0 ± 0.8 (stat) ± 0.7 (syst) ± 0.7 (B 2μ )] × 10 −9 is measured, where the last term is the uncertainty in the normalization channel branching fraction. This work achieves an improved precision of over five orders of magnitude compared to previous results, leading to the first measurement of this branching fraction, which is found to agree with theoretical predictions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Machine-Learning Analysis of Radiative Decays to Dark Matter at the LHC

The search for weakly interacting matter particles (WIMPs) is one of the main objectives of the High Luminosity Large Hadron Collider (HL-LHC). In this work we use Machine-Learning (ML) techniques to explore WIMP radiative decays into a Dark Matter (DM) candidate in a supersymmetric framework. The minimal supersymmetric WIMP sector includes the lightest neutralino that can provide the observed DM relic density through its co-annihilation with the second lightest neutralino and lightest chargino. Moreover, the direct DM detection cross section rates fulfill current experimental bounds and provide discovery targets for the same region of model parameters in which the radiative decay of the second lightest neutralino into a photon and the lightest neutralino is enhanced. This strongly motivates the search for radiatively decaying neutralinos which, however, suffers from strong backgrounds. We investigate the LHC reach in the search for these radiatively decaying particles by means of cut-based and ML methods and estimate its discovery potential in this well-motivated, new physics scenario. We demonstrate that using ML techniques would enable access to most of the parameter space unexplored by other searches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Estimate of Gamma Dose Rates from Arrays of Fermi-1 Blanket Elements During the MEDE Process

The Enrico Fermi Atomic Power Plant Unit 1 (Fermi-1) was a sodium-cooled fast breeder reactor located in Monroe County, Michigan. The reactor was powered by a core of enriched uranium metal alloy driver fuel, which was enveloped by an axial and radial blanket material consisting of depleted uranium metal alloyed with 3 wt.% molybdenum. There are 406 axial and 559 radial irradiated sodium-bonded Fermi-1 blanket assemblies in storage at INL, totaling 34 metric tons of heavy metal. Disposal of the Fermi-1 blanket material directly into a geological repository is prohibited due to the reactive characteristic of its bond sodium. A melt drain evaporate (MEDE) process can effectively remove bond sodium from Fermi-1 blanket material. Consequently, planning is underway to apply a MEDE process to treat the 34 metric tons of heavy metal of irradiated Fermi-1 blanket material. Given the irradiated Fermi-1 blanket material’s relatively low power history and 50+ years of decay, it can be treated in a shielded glovebox. To assess the requisite shielding, the dose rates of the Fermi-1 blanket elements in various process configurations are needed. The method to perform the dose rate calculations is to first generate an average source term and associated photon source spectra for the Fermi-1 blanket material in SCALE and then use the associated photon spectra to calculate dose rates using MCNP6.2 in seven representative geometries. Dose rates in rem/h were calculated on contact (1 cm from outer geometry surface), 30 cm away from outer geometry surface, and 1 m from outer geometry surface at axial heights spanning the length of the blanket material within the geometry. The maximum average dose rate for the single Fermi-1 radial blanket element is ~0.60 rem/h, for the Fermi-1 radial blanket assembly ~1.10 rem/h, for the single Fermi-1 axial blanket element ~0.5 rem/h, for the Fermi-1 axial blanket assembly ~1.01 rem/h, for the MEDE can ~ 0.98 rem/h, for the Fermi Storage Canister ~ 0.83 rem/h, and for the MEDE cans in the DOE Standard Canister ~0.51 rem/h. Based on these maximum average values, the bounding dose rates are assumed to be 3× the average dose rates calculated for the axial blanket material and 7.6× the average dose rates calculated for the radial blanket material, based on the distribution of Fermi-1 blanket material assembly burnup.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗