Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EXPECTATION VALUE”

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,405 records · Page 78

Recent increases in annual, seasonal, and extreme methane fluxes driven by changes in climate and vegetation in boreal and temperate wetland ecosystems

Climate warming is expected to increase global methane (CH 4 ) emissions from wetland ecosystems. Although in situ eddy covariance (EC) measurements at ecosystem scales can potentially detect CH 4 flux changes, most EC systems have only a few years of data collected, so temporal trends in CH 4 remain uncertain. Here, we use established drivers to hindcast changes in CH 4 fluxes (FCH 4 ) since the early 1980s. We trained a machine learning (ML) model on CH 4 flux measurements from 22 [methane-producing sites] in wetland, upland, and lake sites of the FLUXNET-CH 4 database with at least two full years of measurements across temperate and boreal biomes. The gradient boosting decision tree ML model then hindcasted daily FCH 4 over 1981-2018 using meteorological reanalysis data. We found that, mainly driven by rising temperature, half of the sites (n = 11) showed significant increases in annual, seasonal, and extreme FCH 4 , with increases in FCH 4 of ca. 10% or higher found in the fall from 1981–1989 to 2010–2018. The annual trends were driven by increases during summer and fall, particularly at high-CH 4 -emitting fen sites dominated by aerenchymatous plants. We also found that the distribution of days of extremely high FCH 4 (defined according to the 95th percentile of the daily FCH 4 values over a reference period) have become more frequent during the last four decades and currently account for 10–40% of the total seasonal fluxes. The share of extreme FCH 4 days in the total seasonal fluxes was greatest in winter for boreal/taiga sites and in spring for temperate sites, which highlights the increasing importance of the non-growing seasons in annual budgets. Our results shed light on the effects of climate warming on wetlands, which appears to be extending the CH 4 emission seasons and boosting extreme emissions.

54 ENVIRONMENTAL SCIENCES↗

Higher than expected N 2 O emissions from soybean crops in the Pampas Region of Argentina: Estimates from DayCent simulations and field measurements

In developing countries, agriculture generally represents a large fraction of GHG emissions reported in National Inventories, and emissions are typically estimated using Tier 1 IPCC guidelines. However, field data and locally adapted simulation models can improve the accuracy of IPCC estimations. In this report we aimed to quantify anthropogenic N2O emissions from croplands of Argentina through field measurements, model simulations and IPCC guidelines. Here we measured N 2 O emissions and their controlling factors in 62 plots of the Pampas Region with corn, soybean and wheat/soybean crops and in unmanaged grasslands. We accounted for gross emissions from crops and background emissions from unmanaged grasslands to calculate net anthropogenic emissions from crops as the difference between them. We calibrated and evaluated the DayCent model and then simulated different weather and management scenarios. Finally, we applied IPCC guidelines to estimate anthropogenic N 2 O emissions at the same plots. The DayCent model accurately simulated annual N 2 O emission for all crops as compared to measured data (RMSE = 1.4 g N ha -1 day -1 ). Measured and simulated emissions in soybean crops were higher than in corn and wheat/soybean crops. Gross N 2 O emissions ranged from 1.4 to 5.1 kg N ha -1 yr -1 for current environmental (soil and weather) and management (crops and fertilizer doses) conditions. Background emissions ranged between 1.1 and 1.3 kg N ha -1 yr -1 , and therefore net anthropogenic emissions ranged from 0.3 to 4.0 kg N ha-1 yr -1 . IPCC Tier 1 emission factors underestimated N2O releases from soybean, that were on average 4.87 times greater when estimated with DayCent and observations (0.53 vs 2.47 and 2.69 kg N ha -1 yr -1 , respectively). On the contrary, IPCC estimates for corn and wheat/soybean crops were similar to modeled and measured values. Our results suggest that N 2 O emissions from the vast 15 million ha of soybean croplands in the Pampas Region may be substantially underestimated.

59 BASIC BIOLOGICAL SCIENCES↗

Isotropic nanocrystalline Sm(Fe,Co) 11.3 Ti 0.7 magnets modified with B and Zr

Rare-earth-lean Sm(Fe,Co,Ti) 12 alloys with the ThMn 12 crystal structure and less than one Ti atom per formula unit have the potential of exceptionally powerful permanent magnets, but all prior attempts to develop high coercivity in bulk alloys, especially coercivity combined with crystallographic texture, have fallen short of the expectations. This study was aimed at improvement of the currently best Sm(Fe,Co,Ti) 12 magnets prepared through melt-spinning which are inherently isotropic. Modifications of the alloys with B and Zr, already demonstrated in earlier studies to be effective separately, have been implemented simultaneously. Here, a systematic study of Sm 1.1-x (Fe,Co) 11.3-y Ti 0.7 B y alloys melt-spun at a tangential speed of 50 m/s and annealed at 600–950 °C allowed for monitoring the continuous evolution of the two consecutive crystal structures, those of the TbCu 7 and ThMn 12 types. Zirconium was found to facilitate the formation of the 1:12 structure at the expense of the 1:7, whereas boron has the opposite effect, at certain concentrations completely suppressing the 1:12. When the two alloying elements are introduced simultaneously, they inhibit growth of the 1:12 crystallites at annealing temperatures higher than 800 °C, thus allowing for the development of a higher coercivity. Because of instrumental limitations, bulk magnets were prepared through a two-step process – compaction of the melt-spun ribbons at 650 °C and additional treatment at a higher temperature – and they were characterized by a reduced, 90–93%, density. Nevertheless, an isotropic Sm 0.9 Zr 0.2 (Fe,Co) 10.8 Ti 0.7 B 0.5 magnet exhibited fair values of the remanence (7.4 kG), maximum energy product (8.5 MGOe) and coercivity (5.4 kOe), as well as high Curie temperature of 525 °C and remarkably small temperature coefficient of the coercivity, -0.25%/°C.

36 MATERIALS SCIENCE↗

Phase Equilibria and Thermochemistry of Advanced Fuels: Modeling Burnup Behavior (Final Report)

Achieving the goal of developing advanced fuel concepts that meet the DOE objectives of being robust, demonstrating high performance, and are more tolerant of accident conditions than current fuel systems will require a thorough understanding of the thermophysical and thermochemical properties of the constituent materials. Non-oxide fuel systems are being explored under the Advanced Fuels Program that hold significant promise for improved performance and accident tolerance, including the uranium silicide-based system considered in the current work. Prospective cladding materials currently considered that contribute to improved accident tolerance include silicon carbide composites and ferritic alloys (Fe-Cr-Al base compositions). Thus, the effort developed thermochemical models and values, supported with targeted experiments, to evaluate the ferritic alloy and silicon carbide composite cladding systems in contrast to current zirconium alloy cladding. The developed detailed understanding will serve to aid in relatively early screening of candidate systems to avoid wasted effort, guide development of new fuel forms, and to provide a basis for predicting and modeling fuel performance. Major deliverables for the project included: Thermochemical assessment and models of phases in the U-Si and U-Si-N systems; thermochemical evaluation supported by experimental measurements of fuel-cladding interactions of silicide fuel with baseline zirconium, ceramic composite, and ferritic alloy cladding; thermochemical assessment and models of phases supported by experimental measurements for silicide fuel with key fission products provided in a dataset and reported in refereed publications. Within the project a significantly refined U-Si phase diagram was developed and reported that now includes homogeneity ranges for key phases, such at the U 3 Si 2 proposed fuel phase, and settles issues with regard to uncertainty in the stability of some phases. Computational efforts together with key experiments has determined phase formation in interactions between U 3 Si 2 and Zircaloy-4 cladding material, a ferritic FeCrAlY alloy of interest as an advanced cladding material, and silicon carbide, also of interest as a fiber-reinforced composite cladding. As expected, very significant reactions occur between U3Si2 and Zircaloy-4, with much less interaction at higher temperatures for the ferritic alloy, and finally interactions with SiC only in the region of contact. A potentially major issue is the stability of U 3 Si 2 fuel that has undergone significant burnup. The result is the loss of the uranium metal, liberating silicon, and the formation of concomitant fission product elements that either dissolve in the U 3 Si 2 phase or form independent, and possibly silicide phases. A combination of experimental determinations of phase formation of U 3 Si 2 reacted with representative fission products yttrium, gadolinium, cerium, zirconium, and molybdenum and first principles calculations has helped understand the fuel chemistry. The behavior of the U 3 Si 2 phase and the partitioning of silicon to possible fission product phases with burnup was thus determined, with significant dissolution in U 3 Si 2 of cerium, gadolinium, zirconium, and plutonium predicted along with independent phase formation of a U-Mo-Si ternary phase, yttrium silicide, and elemental selenium. It can be concluded that at significant burnup there will be a very minor amount of the U 3 Si 2 fuel phase that will decompose to a lower silicide or a uranium alloy as silicon preferentially forms a secondary phase.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Multi-channel, multi-template event reconstruction for SuperCDMS data using machine learning

SuperCDMS SNOLAB uses kilogram-scale germanium and silicon detectors to search for dark matter. Each detector has Transition Edge Sensors (TESs) patterned on the top and bottom faces of a large crystal substrate, with the TESs electrically grouped into six phonon readout channels per face. Noise correlations are expected among a detector's readout channels, in part because the channels and their readout electronics are located in close proximity to one another. Moreover, owing to the large size of the detectors, energy deposits can produce vastly different phonon propagation patterns depending on their location in the substrate, resulting in a strong position dependence in the readout-channel pulse shapes. Both of these effects can degrade the energy resolution and consequently diminish the dark matter search sensitivity of the experiment if not accounted for properly. We present a new algorithm for pulse reconstruction, mathematically formulated to take into account correlated noise and pulse shape variations. This new algorithm fits N readout channels with a superposition of M pulse templates simultaneously - hence termed the N$\times$M filter. We describe a method to derive the pulse templates using principal component analysis (PCA) and to extract energy and position information using a gradient boosted decision tree (GBDT). We show that these new N$\times$M and GBDT analysis tools can reduce the impact from correlated noise sources while improving the reconstructed energy resolution for simulated mono-energetic events by more than a factor of three and for the 71Ge K-shell electron-capture peak recoils measured in a previous version of SuperCDMS called CDMSlite to $<$ 50 eV from the previously published value of $\sim$100 eV. These results lay the groundwork for position reconstruction in SuperCDMS with the N$\times$M outputs.

Albakry, M. F. [British Columbia U.; TRIUMF]↗

Property Measurements of the NaCl-PuCl 3 Molten Salt System

Thermal properties of several compositions of the binary NaCl-PuCl 3 salt were measured to confirm the eutectic composition and provide high quality data for use by MSR developers. The PuCl 3 was generated by reducing PuO 2 to plutonium metal in the presence of calcium metal and then chlorinating the metallic plutonium with NH 4 Cl. The resulting PuCl 3 was used to make seven binary NaCl-PuCl 3 compositions with between 59.9 mol % PuCl 3 (Salt 1) and 20.0 mol % PuCl 3 (Salt 7). The eutectic temperature was measured to be 457 ± 4 °C by onset determination in DSC analyses of the seven mixtures, which is consistent with the values in the literature. Additional transitions were observed at approximately 332 and 363 °C in analyses of the PuCl 3 -rich compositions (Salts 1-3) at about 376 and 439 °C in analyses the Na-rich compositions (Salts 4-7). The solid state heat capacity decreased with increasing PuCl 3 content. The liquid state heat capacity was measured for Salt 3 (37.4 mol % PuCl 3 ) and Salt 4 (38.3 mol % PuCl 3 ). A higher heat capacity was measured for Salt 4 and measurements with both salts showed a positive correlation with temperature between 520 and 730 °C. Cells fabricated from nickel and molybdenum for use in DSC measurements at high temperatures could not be adequately sealed because these materials are not sufficiently malleable, even after high temperature annealing. It is recommended that future development focus on corrosion-resistant materials that have mechanical and thermal properties similar to the commercially available gold cells, such as Pt-Rh (80–20) and pure platinum. Those materials are chemically inert, soft, and malleable like gold, but have melting temperatures higher than 1000 °C (Rakhtsaum, 2013). Use of the Pt-Rh alloy should be evaluated first based on superior machinability for making precision parts and its current use in commercial DSC high temperature crucibles. The commercial Pt-Rh cells are not hermetically sealable, but it is expected that cells made from thinner stock can be sealed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A Refined BISON-Griffin Workflow For Designing Targeted TRISO Transient Experiments In The Transient Reactor Test Facility

U.S. industry is expected to deploy tristructural isotropic (TRISO) particle fuel technologies for commercial reactors within the next decade. In our previous work, we defined a preliminary transient design space for TRISO fuels, identified potential gaps in the available data, and began to develop multiphysics modeling tools that could be applied to design targeted Transient Reactor Test Facility (TREAT) experiments to fill these gaps. This work builds on that foundation by (1) updating BISON fuel performance and Griffin reactor physics models to reflect the current TREAT experiment tube and capsule designs,(2) coupling the codes to improve the accuracy and usability of the transient design analyses, and (3) demonstrating their use over an expanded design space that includes fuel burnup. The simulated mechanical responses of the TRISO particles were complex functions of fission product accumulation, fission gas release, and irradiation-induced dimensional change in the pyrolytic carbon layers. Predicted tangential stresses in particle silicon carbide layers were least compressive for preheated tests involving fresh fuels but remained compressive throughout the ranges of temperature, heat rate, and burnup considered in this work. Finally, comparisons between the potential TREAT transients and historical test reactor irradiations showed that the TREAT tests would produce significantly lower average energy deposition rates, yielding less severe transients with greater relevance to near-term commercial applications. Use of these predictive capabilities has the potential to increase the value of each test, improving the overall efficiency and cost effectiveness of transient testing for TRISO and other advanced fuels.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Monte Carlo Simulation and Reconstruction: Assessment of Myocardial Perfusion Imaging of Tracer Dynamics With Cardiac Motion Due to Deformation and Respiration Using Gamma Camera With Continuous Acquisition

Purpose: Myocardial perfusion imaging (MPI) with single photon emission computed tomography (SPECT) is routinely used for stress testing in nuclear medicine. Recently, our group extended its potential going from 3D visual qualitative image analysis to 4D spatiotemporal reconstruction of dynamically acquired data to capture the time variation of the radiotracer concentration and the estimated myocardial blood flow (MBF) and coronary flow reserve (CFR). However, the quality of reconstructed image is compromised due to cardiac deformation and respiration. The work presented here develops an algorithm that reconstructs the dynamic sequence of separate respiratory and cardiac phases and evaluates the algorithm with data simulated with a Monte Carlo simulation for the continuous image acquisition and processing with a slowly rotating SPECT camera. Methods: A clinically realistic Monte Carlo (MC) simulation is developed using the 4D Extended Cardiac Torso (XCAT) digital phantom with respiratory and cardiac motion to model continuous data acquisition of dynamic cardiac SPECT with slowly rotating gamma cameras by incorporating deformation and displacement of the myocardium due to cardiac and respiratory motion. We extended our previously developed 4D maximum-likelihood expectation-maximization (MLEM) reconstruction algorithm for a data set binned from a continuous list mode (LM) simulation with cardiac and respiratory information. Our spatiotemporal image reconstruction uses splines to explicitly model the temporal change of the tracer for each cardiac and respiratory gate that delineates the myocardial spatial position as the tracer washes in and out. Unlike in a fully list-mode data acquisition and reconstruction the accumulated photons are binned over a specific but very short time interval corresponding to each cardiac and respiratory gate. Reconstruction results are presented showing the dynamics of the tracer in the myocardium as it continuously deforms. These results are then compared with the conventional 4D spatiotemporal reconstruction method that models only the temporal changes of the tracer activity. Mean Stabilized Activity (MSA), signal to noise ratio (SNR) and Bias for the myocardium activities for three different target-to-background ratios (TBRs) are evaluated. Dynamic quantitative indices such as wash-in (K1) and wash-out (k2) rates at each gate were also estimated. Results: The MSA and SNR are higher with higher TBRs while biases were improved with higher TBRs to less than 10%. The correlation between exhalation-inhalation sequence with the ground truth during respiratory cycle was excellent. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Estimated values of K1 and k2 were also consistent with the ground truth. Conclusion: The continuous image acquisition for dynamic scan using conventional two-head gamma cameras can provide valuable information for MPI. Our study demonstrated the viability of using a continuous image acquisition method on a widely used clinical two-head SPECT system. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Precise implementation of reconstruction algorithms, better segmentation techniques by generating images of different tissue types and background activity would improve the feasibility of the method in real clinical environment.

60 APPLIED LIFE SCIENCES↗

Refinement and demonstration of a coupled BISON-Griffin workflow for designing targeted TRISO transient experiments in TREAT

The U.S. nuclear industry is expected to deploy tristructural isotropic (TRISO) particle fuel technologies for commercial reactors within the next decade. In previous work, we defined a preliminary transient design space for TRISO fuels, identified potential gaps in the available data, and began to develop multiphysics modeling tools that could be applied to design targeted Transient Reactor Test Facility (TREAT) experiments to fill these gaps. Here, this work builds on that foundation by (1) updating BISON fuel performance and Griffin reactor physics models to reflect the current TREAT experiment tube and capsule designs, (2) coupling the codes to improve the accuracy and usability of the transient design analyses, and (3) demonstrating their use over an expanded design space that includes fuel burnup. The simulated mechanical responses of the TRISO particles were complex functions of fission product accumulation, fission gas release, and irradiation-induced dimensional change in the pyrolytic carbon layers. The predicted tangential stresses in the particles' silicon carbide layers were least compressive for preheated tests involving fresh fuels but remained compressive throughout the ranges of temperature, heat rate, and burnup considered in this work. Finally, comparisons between the potential TREAT transients and historical test reactor irradiations showed that the TREAT tests would produce significantly lower average energy deposition rates, yielding less severe transients with greater relevance to near-term commercial applications. The use of these predictive capabilities has the potential to increase the value of each test, improving the overall efficiency and cost-effectiveness of transient testing for TRISO and other advanced fuels.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

H-Canyon Jumper Gasket Sealing Evaluation

H-Canyon Documented Safety Analysis approved December 2019 credits the leak integrity of H-Canyon (HCA) process tanks, piping, and jumpers in preventing the release of radioactive material in a post seismic environment. Demonstration of leak tightness capability for this equipment was achieved through substantial analytical effort. However, one potential weak point that affects leak integrity exists at the gasketed joint between Hanford Connectors (HC) and nozzles. Review of the HCA Crane Log over a three-year period indicates the HCs periodically leak. Savannah River National Laboratory (SRNL) was requested to perform gasket compression tests that evaluate gasket relaxation characteristics, i.e., creep, and provide improvements to the HC installation process that can reduce or eliminate the frequency of leaks at HCs. A HC torque value of 312 ft-lbf is the target value used because it is the minimum torque value, with 95% confidence, when an unlubricated HC is impacted for standard sealing period of 9 seconds. In the HCA pipe-jumper connections are installed remotely by using HCs. HCs inherently add to the uncertainty of a gasketed joint because a calibrated torque value cannot be applied. The gasket material is a special material made of asbestos-fiber impregnated with Teflon. This material has been used for decades and generally performs well but does leak intermittently as gaskets creep and preloaded torque is lost. Gasket types are either snap-ring (SR) or lampshade (LS). Both gasket types cover the same seating area. The SR gasket is designed to remain with the jumper when the jumper is removed leaving the nozzle faces clean of any debris. SR gaskets are initially installed when the jumper is fabricated. Changing gaskets is time consuming and the radioactive environment presents risks to personnel. When jumpers are removed LS gaskets could be used to better seal the connection to reduce down time and ALARA concerns. That is, an LS gasket would be remotely positioned over the nozzle and the HC reinstalled. Testing also included two situations not typical of plant operation. One situation was the use of the non-standard double-layered gaskets (DG), and the other situation was the use of a LS without removal of the installed SR gasket. The DG was used when the supply of LS gaskets ended, so that testing could continue. It was of interest to demonstrate this gasket’s creep response because the DG, which was two gaskets glued together, was thicker than the LS gasket and the thicker material was expected to elicit a greater creep characteristics. The second situation of using two gaskets together by installing an LS gasket without removing the installed gasket is not allowed by the ASME Codes. However, by demonstrating the creep history of such a gasket arrangement it could be considered as an option that could possibly seal a connection while significantly reducing operational down time and personnel risk. Forty tests were conducted, and the overall results demonstrate the initial applied force and how the gasket changes with time. The change in sealing force, in terms of torque, was then measured to obtain a force history as the gaskets creep. This report documents the results of those tests as well as the method used to perform the tests.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The Contribution of Proton-Donor pKa on Reactivity Profiles of [FeFe]-hydrogenases

The [FeFe]-hydrogenases are enzymes that catalyze the reversible activation of H 2 coupled to the reduction–oxidation of electron carriers. Members of the different taxonomic groups of [FeFe]-hydrogenases display a wide range of preference, or bias, for H 2 oxidation or H 2 production reactions, despite sharing a common catalytic cofactor, or H-cluster. Identifying the properties that control reactivity remains an active area of investigation, and models have emerged that include diversity in the catalytic site coordination environments and compositions of electron transfer chains. The kinetics of proton-coupled electron transfer at the H-cluster might be expected to be a point of control of reactivity. To test this hypothesis, systematic changes were made to the conserved cysteine residue that functions in proton exchange with the H-cluster in the three model enzymes: CaI, CpII, and CrHydA1. CaI and CpII both employ electron transfer accessory clusters but differ in bias, whereas CrHydA1 lacks accessory clusters having only the H-cluster. Changing from cysteine to either serine (more basic) or aspartate (more acidic) modifies the sidechain pKa and thus the barrier for the proton exchange step. The reaction rates for H 2 oxidation or H 2 evolution were surveyed and measured for model [FeFe]-hydrogenases, and the results show that the initial proton-transfer step in [FeFe]-hydrogenase is tightly coupled to the control of reactivity; a change from cysteine to more basic serine favored H 2 oxidation in all enzymes, whereas a change to more acidic aspartate caused a shift in preference toward H 2 evolution. Overall, the changes in reactivity profiles were profound, spanning 10 5 in ratio of the H 2 oxidation-to-H 2 evolution rates. The fact that the change in reactivity follows a common trend implies that the effect of changing the proton-transfer residue pKa may also be framed as an effect on the scaling relationship between the H-cluster di(thiolmethyl)amine (DTMA) ligand pKa and E m values of the H-cluster. Experimental observations that support this relationship, and how it relates to catalytic function in [FeFe]-hydrogenases, are discussed.

59 BASIC BIOLOGICAL SCIENCES↗

Passive heat removal in horizontally oriented micro-HTGRs

There are novel, horizontally oriented high-temperature gas-cooled micro-reactor (HTGR) designs under consideration because of their portability, simplicity, and reliable operation. These features of the micro-HTGR design can meet the fission battery attributes, such as economic, standardize, installed, unattended and reliable. Although there have been HTGR related safety studies, micro-HTGRs when horizontally oriented are expected to exhibit significantly different thermal physics. During the unavailability of helium circulation, the internal reactor core is designed to cool by block-to-block conduction and radiation, and the reactor vessel surface is cooled by the ambient air. This scenario is anticipated during the transport of the micro-HTGR in a shipping container. The conduction and radiation between the prismatic micro-HTGR blocks in the core can be influenced by variances in the thermal contacts. This work investigated the conduction within a simulated horizontal HTGR core. An experimental setup was used to validate a numerical model. The experimental setup consisted of a hexagonal assembly with scaled prismatic blocks placed within a high-temperature vacuum environment. The gaps between the blocks were well controlled and monitored. The experimental setup was designed in such a way that the temperature variation in the axial direction was minimal, such that the experiment could be observed as a 2D (r, ) heat transfer problem. The pressurized tube was equipped with Infra-Red (IR) transparent windows, and an IR camera was used to capture the spatio-temporal temperature evolution in the hexagonal block assembly during the cooldown experiments. The experimental scenario was computationally modeled with a finite element analysis (FEA) program. Once validated, the computational model was used to investigate the impact of gap conductance on overall decay heat removal. Using a conservative estimate for gap conductance value (100 W/m 2 -K) between the prismatic blocks, there is a negligible increase in temperature observed during decay heat generation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Phase Field Dislocation Dynamics (PFDD) version 2.x

This disclosure is for version 2.x of a mesoscale model called Phase Field Dislocation Dynamics (PFDD). PFDD is used for investigating deformation in nanoscale (grain sizes of ~300 nm and less) materials, such as metals and alloys. This approach models the motion and interaction of individual defects, namely dislocations, in the material using scalar-valued phase field variables, also called order parameters. The system is evolved through energy minimization thus the model calculates the total energy density in terms of the phase field variables. The energy minimization is completed using the Ginzburg-Landau equation, and is implemented with explicit time integration. The total system energy can be comprised of several terms, including the strain energy (which describes dislocation-dislocation interactions), the energy due to an applied stress (dislocation interactions with the applied stress), and a core/lattice (perfect dislocations) or generalized stacking fault (partial dislocations) energy (described the dislocation core structure). The latter term in particular may vary based on the crystal structure being modeled and is typically informed using lower length scale (e.g., atomistic) approaches, although no such (atomistic) calculations are completed within the PFDD algorithm. This basic formulation was previously reviewed by Los Alamos National Laboratory and released under license number C17113. This previously reviewed version we will henceforth refer to as PFDD v1.0. PFDD v1.0 consisted of 2 codes (one parallel and one serial) plus input files, all written in the C language. This new disclosure is addressing the next versions of the PFDD, versions 2.x. There have been several enhancements of PFDD v1.0, which are described here and included in the attached code, which we will refer to as PFDD v2.0. There are also several new features described here that are either planned or already in process and are expected to be subsequent releases, i.e., v2.1, v2.2, ...v2.x.

Hunter, Abigail↗

Tensile Property and Lifetime Prediction for Low-Temperature Aged Uranium-Niobium Alloys

Thermal aging models and lifetime predictions for uranium-niobium (U-Nb) alloys were created using an approach similar to those previously employed. Lifetime estimates for generic U-6Nb components were thus updated; the reported value being 800 years. This update represents a small change in lifetime vs. that of the 2012 assessment (540 years). This lifetime estimate emerged from consideration of several model fits specific to the aging datasets and properties chosen. Aging was quantified using quasi-static tensile properties measured on specimens artificially aged for up to 10 years. The major change relative to the most recent 2012 LANL assessment was that a more comprehensive body of U-Nb literature data was mined, in addition to being augmented by the latest LANL and UK AWE data. The tensile data compilation was published separately (LANL report LA-14493, December 2016). Recognizing the chemical banding of industrially produced U-6Nb, models were developed for the mid-range (6 wt.%) and extrema (4 and 8 wt.%) compositions. Lifetime estimates were calculated for all three alloy classes (4, 6, 8 wt.% nominal) and two measures of total tensile elongation (TE) to failure, namely TE-ext. — extensometer method, and TE-NCD — normalized crosshead displacement method. The conservative assumption was made that whichever composition (4, 6, or 8 wt.% Nb) and property (TE-ext or TE-NCD) was the first to cross the ductility failure threshold would limit the lifetime of the entire component. Tensile strength properties did not figure into the lifetime predictions, but could be useful as age-sensitive diagnostics and were also modeled. Of these, only first yield strength is expected to show a change at 40°C aging vs. time = 0 over the ~100-year timespan of engineering interest. Second yield strength evolves more slowly, and ultimate tensile strength slower still. Among all the models, the apparent activation energies for aging were mostly in the narrow 29– 37 kcal/mol range, which is close to that for diffusion of Nb in gamma-uranium. This agreement may be coincidental. The data from recent long-term aging studies substantially improved the model fit quality and robustness of the lifetimes. Appendices document sensitivity studies of the model fits and lifetimes with respect to using more limited datasets. These results highlight the limitations of relying solely on data from scattered literature studies and smaller datasets more generally.

36 MATERIALS SCIENCE↗

Markets and Economic Requirements for Fission Batteries and Other Nuclear Systems

Fission Batteries (FBs) are nuclear reactors defined by five characteristics which enable large-scale deployment: cost competitive, standardized sizes for economic mass production, easy installation and removal, secure and safe unattended operation with high reliability. FBs are not defined by technology or power level. Technical and market considerations suggest that most FBs will produce 20 to 30 MWt. This proceedings reports on the outcomes of two workshops that were held in January 2021 to better define markets and economic challenges for FBs. Three major markets were identified. The largest market is the industrial and commercial heat market. There are about 4000 industrial users (excluding utilities) that require more than one megawatt of heat. The number of customers versus size of heat demand was determined. In a low-carbon world there is the potential for many additional customers—including expanded biofuels production and district heat. The second market is for non-grid electricity. This includes co-generation plants that produce heat and electricity for a single customer. The third market is the maritime market with ~100,000 ships worldwide. In the United States, natural gas is the low-cost energy option today and will remain so unless constraints or taxes impact its use. If restrictions on greenhouse gas emissions, the FB competition includes natural gas with carbon capture, biofuels, hydrogen and grid electricity. Natural gas with carbon capture is not economically viable on a small scale. Biofuels may be expensive but may be the economically preferred option for locations with small energy demands of a few megawatts. Hydrogen is a potential competitor with many of the characteristics of natural gas. Grid electricity is not a competitive source of heat. For FBs to be economically competitive, the price of delivered heat must be $20-50/MWh ($6-15/million BTU). The economically competitive range for non-grid electricity is estimated at $70-100/MWh. These electricity prices are competitive with the retail prices of electricity in many parts of the United States for the customer. FBs are not expected to be competitive selling wholesale electricity to the grid. To achieve the aforementioned cost targets for heat and electricity markets, FB designers must (1) maximize the power output within the constraints of a FB (e.g., truck transportability, passive decay heat removal), (2) drastically reduce the size of onsite staff, (3) adopt core designs with low fuel costs (enrichment and fabrication), and (4) develop a system design that is efficiently manufactured in factories. The business case depends upon more than being just a replacement for natural gas. The largest incentives for adoption of FBs is where they create new markets and new sources of revenue. An example is the paper and pulp industry that burns biomass wastes to provide heat and electricity to make paper. An external heat source could meet the demand for heat and electricity by the paper process and enable converting waste biomass into liquid biofuels rather than burning to provide heat. Other markets, such as data centers, are driven by special energy requirements such as extreme reliability. Most customers are not in the energy business but need heat and electricity to produce a product—a manufactured good, education, retail sales (shopping malls), marine transport or some other product. As a consequence, there will be large incentives to lease rather than own FBs. Leasing avoids the regulatory challenges that remain with the owner of the FB. Leasing creates large incentives for FP standardization of sizes and transportability to maintain the value of the FB at the end of the lease—similar to the leasing of jet engines and aircraft. The economic constraints combined with technical constraints suggest competitive FBs will likely have outputs exceeding 10 MWt. There appear to be little incentives for very long-lived reactor cores because such machines require much larger inventories of fuel. Maintenance requirements and the options to provide technology updates may favor shorter lifetimes (~5 years). The assessment is that there is the potential for FBs to be economically viable and play a major role in global decarbonization in three markets: heat, non-grid electricity and maritime applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Superconductivity at carrier density 10 17 cm – 3 in quasi-one-dimensional Li 0.9 Mo 6 O 17

Superconductivity is a fascinating phenomenon that involves an attractive interaction by which electrons are paired and can move without electrical resistance below a critical temperature T c . A paradigm for understanding superconductive pairing, the successful “BCS” theory, developed more than 65 years ago, is challenged by certain classes of materials for which parameters (e.g. the value of T c , the density of electrons, etc.) fall outside the range for which BCS theory is applicable. The study of such materials is important because new physics may be revealed and because superconductivity has many potential applications in sensing and electronics. Superconductivity at very low carrier density is a rare and fascinating phenomenon that is of considerable importance and interest in condensed matter physics. The conditions of high density of states and effective screening of Coulomb repulsion, both key ingredients of stable Cooper pairing in the BCS theory, may not be met at low carrier density. Thus, such systems may also harbor unconventional electron pairing (e.g. non-phononic) or pairing without phase coherence as a precursor to Bose-Einstein condensation. An intriguing route toward extreme low-density (and possibly unconventional) superconductivity may be found in quasi-one-dimensional metals for which density-wave instabilities compete for the ground state and can lead to gapping of substantial portions of the Fermi surface (FS), leaving a residual (reconstructed) FS comprising a very small carrier density. Our work highlights a novel material in this unusual category: Li 0.9 Mo 6 O 17 (T c =2 K) (also known as "lithium purple bronze" or "LiPB"). We present both electrical and thermoelectric magnetotransport measurements that reveal superconductivity to occur at carrier density 2 x10 17 cm -3 , among the lowest known for any superconductor, and four orders of magnitude lower than is expected from its chemical valence and band structure. With its exceptionally large ratio of superconducting to Fermi temperature, T c /TF ≳ 0.1, LiPB is properly placed at the top of a very small group of superconductors (including cuprates, heavy-fermion and Fe-based superconductors) for which this ratio falls in the BCS-BEC crossover regime and for which pairing is believed to be spin mediated.

36 MATERIALS SCIENCE↗

Mechanical Loads on Spent Nuclear Fuel in the General 30 cm Package Drop Scenario

The U.S. Department of Energy Office of Nuclear Energy (DOE-NE) Spent Fuel and Waste Science & Technology (SFWST) research program is guided by the high-level goal of closing prioritized knowledge gaps related to spent nuclear fuel (SNF) storage and transportation, which are summarized by Saltzstein et al. (2020). One of the high-priority knowledge gaps is the identification and quantification of mechanical loads that are expected to affect SNF during normal conditions of transportation and storage to inform the range of physical SNF test programs. This report uses modeling and analysis methods to estimate the mechanical loads on spent nuclear fuel (SNF) in the general 30 cm package drop scenario. The drop scenario assumes impact limiters are in place in the transportation configuration and the impact surface is perfectly rigid. The goal of this analysis is to consider the universe of potential mechanical loading conditions that can happen to SNF and present the results in a manner that is useful for materials testing, decision making, and regulatory rule making purposes. This study uses validated finite element models and methods to perform a broad parametric study of key variables that can affect the mechanical loads on SNF during a hypothetical package free drop scenario. Physical drop test data from a cask and fuel assembly drop test campaign is the basis for model validation. Additionally, the results of the parametric study are used to inform a damage model, which uses multiple nonlinear regression to estimate the relationships between input variables and output response. The parametric finite element analyses consider thousands of input variable combinations, while the damage model estimates millions of combinations. The breadth of this study provides confidence that the potential range of mechanical loads that SNF might experience during the general 30 cm package drop scenario are characterized well enough to consider this knowledge gap closed. While this report documents the overall peak values calculated in this study, the 95 th percentile values, the histograms, and the observed trends are equally important. This study covered a large range of SNF temperatures, room temperature to 300°C, and burnups, 10 GWd/MTU to 62 GWd/MTU. Each temperature and burnup combination has a different cladding yield strain, so it is more meaningful to summarize the calculated cladding strain response as its factor of safety, which is defined relative to the yield strain. The factor of safety is calculated as the yield strain divided by the peak cladding strain. A factor of safety greater than unity indicates that the cladding remains below yield, whereas a value less than unity is indicative of plastic deformation. In all cases of this study a safety factor over 1.0 was calculated, although in the most limiting case at 300°C the safety factor was only 1.01, which suggests that yielding could occur when additional loads like rod internal pressure are included. When the temperature is restricted to 200°C the limiting safety factor increases to 1.28, which has significant margin to accommodate internal pressure and potential local cladding defects that could cause a local stress concentration. An important trend in the calculated fuel rod mechanical loads is that the 2 nd highest loaded fuel rod in an assembly tends to be significantly lower than the highest loaded rod. The implication is that even if one rod in an assembly experiences a failure the loads would have to be significantly higher to cause two or more rods to fail.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Is there warming in the pipeline? A multi-model analysis of the Zero Emissions Commitment from CO 2

Abstract. The Zero Emissions Commitment (ZEC) is the change in global mean temperature expected to occur following the cessation of net CO 2 emissions and as such is a critical parameter for calculating the remaining carbon budget. The Zero Emissions Commitment Model Intercomparison Project (ZECMIP) was established to gain a better understanding of the potential magnitude and sign of ZEC, in addition to the processes that underlie this metric. A total of 18 Earth system models of both full and intermediate complexity participated in ZECMIP. All models conducted an experiment where atmospheric CO 2 concentration increases exponentially until 1000 PgC has been emitted. Thereafter emissions are set to zero and models are configured to allow free evolution of atmospheric CO 2 concentration. Many models conducted additional second-priority simulations with different cumulative emission totals and an alternative idealized emissions pathway with a gradual transition to zero emissions. The inter-model range of ZEC 50 years after emissions cease for the 1000 PgC experiment is -0.36 to 0.29°C, with a model ensemble mean of -0.07 °C, median of -0.05°C, and standard deviation of 0.19°C. Models exhibit a wide variety of behaviours after emissions cease, with some models continuing to warm for decades to millennia and others cooling substantially. Analysis shows that both the carbon uptake by the ocean and the terrestrial biosphere are important for counteracting the warming effect from the reduction in ocean heat uptake in the decades after emissions cease. This warming effect is difficult to constrain due to high uncertainty in the efficacy of ocean heat uptake. Overall, the most likely value of ZEC on multi-decadal timescales is close to zero, consistent with previous model experiments and simple theory.

MacDougall, Andrew H.↗