Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “magnitude”

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

Analysis of curtailment at The Geysers geothermal Field, California

Geothermal energy has traditionally been viewed as a baseload energy source, but the rapid growth of intermittent renewable energy has led to a need for more flexibility in power generation to avoid mandatory curtailment imposed by grid operators. This study of curtailment at The Geysers provides insights into the magnitude, duration, frequency, temporal and spatial distribution, and potential causes of curtailment events between 2013 and 2018. Annual levels of curtailment range during this period from 9 to 47 GW h, representing 0.15 to 0.81 % of the net generation. Most curtailments occurred at the power plants connected to a lower capacity transmission line and may result from transmission constriction. There is a clear link between negative pricing and economic curtailment, especially when solar production is higher. Economic curtailment events tend to be only a few hours and vary in magnitude up to almost 300 MW, whereas transmission-related curtailment events can be up to several weeks in duration. It is likely that curtailment of geothermal power will be an increasing concern, and could be mitigated by flexible generation strategies and increases in energy storage. It is critical to know the nature of curtailment events so that flexible generation options can be assessed properly.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Induced seismicity red-light thresholds for enhanced geothermal prospects in the Netherlands

Prospects for geothermal energy in the Netherlands have renewed concerns around induced earthquakes. Risks from induced earthquakes are managed by traffic light protocols (TLPs), where the red-light is chosen as the stop-point before exceeding a tolerance to risk. Here, we simulate post-shut-in earthquake scenarios based on realistic information for the Netherlands. We focus on three risk metrics: aggregates like nuisance and damage impacts and also local personal risk (LPR) – a likelihood of building collapse fatality for an individual. Our results show that the severity of these risks varies spatially by orders of magnitude. Prior induced seismicity (e.g., the 2012 Huizinge event) provides a reference baseline to calibrate the Dutch earthquake risk tolerances. We find that these calibrated risk tolerances are similar to those observed in North America, suggesting an underlying sociological ‘license to operate.’ Furthermore, the use of calibrated risk tolerances results in nuisance concerns completely eclipsing the other two metrics. We compare our results to a hypothetical Groningen geothermal operation and find that our approach sets red-light thresholds approximately one magnitude unit below the ML 3.6 Huizinge event. Overall, our results provide a first-order recommendation for red-light thresholds and proactive management of Dutch enhanced geothermal induced seismicity.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Geological study of an outburst flood event in the upper Yangtze River and risk of similar extreme events

Here, outburst flood sediments are studied on the upper course of the Yangtze River (known as the Jinsha River) where is a key region for the development of hydroelectric power. As the history and magnitude of outburst flooding in this stretch of the Jinsha River during historic times has remained largely undocumented, we have endeavored as part of this study to conduct field investigations into, and provide numerical dating (using OSL and 14 C dating techniques) of outburst flood evidence found in the First Bend area of the Yangtze River. Carbonized nutshells and charcoal from the ancient settlement yielded 14 C ages that precisely pinpointed an outburst flood event ~1200–1300 cal yr BP, during China's Tang Dynasty (618–907 CE). However, OSL ages using the minimum age model for all outburst flood samples ranged from 9.1 to 1.8 ka, which was ~0.5–7.8 ka older than the 14 C-dated outburst flood. This discrepancy suggests a potential overestimation of quartz age when dating outburst flood sediments from the Holocene period. Furthermore, using 2D hydraulic modeling, we estimated the magnitude of the aforementioned outburst flood, revealing a reconstructed peak flow discharge of 51,000–55,300 m 3 /s. This discharge was three times greater than the 10,000-year flood for the hydroelectric power station nearest to the study area. Our study underscores the significant role of combining archaeological and geological evidence in enhancing paleoflood hydrology research. Through a multidisciplinary approach, our findings emphasize the importance of paleoflood hydrology studies in comprehending both the latent and catastrophic flood risks in alpine mountain areas amid the context of global warming.

14C dating↗
DOE OSTI2025

Characterizing the pore-scale impacts of enhanced-solubilization solutions on organic-liquid morphology and fluid-fluid interfacial area with microtomography

Synchrotron X-ray microtomography was used to directly quantify the impact of enhanced-solubilization solutions on mobilization of a non-wetting fluid in a natural porous medium. Polyoxyethylene sorbitan monooleate (tween 80), hydroxypropyl-β-cyclodextrin (HPCD), sodium dodecyl sulfate (SDS), and ethanol were used as the solubilization reagents. Tetrachloroethene (PCE) was used as the entrapped organic immiscible liquid. Microtomography images were collected prior to and after successive floods with each reagent solution. The total trapping number was varied by more than three orders of magnitude by changing the velocity and viscosity of the aqueous phase as well as the interfacial tension between fluid phases. The results were compared to those obtained from equivalent experiments conducted with water flooding. Inspection of images indicated that mobilization caused a change in PCE morphology and distribution for all reagents with the exception of HPCD, for which a lower trapping number was determined compared to the other experiments. Mobilization of organic liquid was observed at total trapping numbers that were consistent with prior work. The enhanced-solubilization solutions had minimal impact on the configuration or magnitude of fluid-fluid interfacial area.

Interfacial area↗
DOE OSTI2022

Numerical simulation of falling film sensible heat transfer over round horizontal tubes

Numerical simulations are performed to explore the heat transfer characteristics of falling films over horizontal round tubes with uniform heat flux imposed for a range of Reynolds numbers spanning the droplet, jet, and sheet regimes. Simulations results agree well with available experimental measurements. The study analyzes the local as well as average heat transfer behavior under the different flow modes. The numerical results show that the local Nusselt number (Nu) distribution depends on the flow features in each mode and varies substantially in all directions for the respective mode. In the droplet mode, the Nu value varies significantly as the droplet impinges and the remnant liquid-bridge retracts (peak instantaneous Nu near 6), followed by wave propagation over the tube surface with peak Nu around 0.25. For the jet modes, the local maximum in the heat transfer occurs off-center to the impingement location with magnitudes of peak Nu = 3.1 for the inline jet mode and Nu = 2.7 for the staggered jet mode, while on the rest of the tube surface, it has an inverse relation with the liquid film thickness. Substantial variations in the heat transfer value are also recorded in the middle of the two impinging jets with Nu = 0.95 in the inline jet mode where the neighboring jets do not interact, and Nu = 0.60 in the crest region of the staggered jet mode where the neighboring jets interact with each other. In the sheet mode, the Nu was seen to depend on the thickness of the liquid waves traversing over the tube surface. Lower Nu values were recorded beneath the crest location of the liquid waves, which increases (1.4–11.6% depending on circumferential location) abruptly in magnitude at the advancing fronts of the waves. The temperature distribution in the liquid film in each of the modes was examined to evaluate the mechanism of heat transfer process. This study also compares the local heat transfer coefficient distribution with the analytical heat transfer models derived to predict heat transfer performance over horizontal tube surfaces.

42 ENGINEERING↗
DOE OSTI2022

High temperature stability and transport characteristics of hydrogen in alumina via multiscale computation

Here, the impact of hydrogen charge states on the stability and transport characteristics of hydrogen interstitials in alumina polymorphs is evaluated by multiscale computational methods including density functional theory (DFT), ab initio molecular dynamics (AIMD) and machine learned force fields. Thermodynamic calculations show that the protonic H i +1 interstitial is the most stable defect species for most values of the electronic bandgap in both and amorphous alumina (Al 2 O 3 ). Further, active learned Gaussian approximation potentials (GAP) were developed using AIMD data to study temperature dependent long time proton diffusion in alumina. Diffusivity calculations from GAP-MD simulations are found to be comparable with of the AIMD data, while being ~340 times faster and scalable to larger systems. Comparisons with diffusivity values for other interstitial charge states (H i 0 and H i -1 ) and published experimental literature indicate that H i +1 diffusion is the likely mechanism of hydrogen transport. A good agreement is obtained between H i +1 diffusivity calculated in α-Al 2 O 3 from DFT: 5.05 10 -3 exp(-0.81 eV/k B /T) cm 2 /s and reported experiment: 9.7X10 -4 exp(-0.83 eV/k B /T) cm 2 /s. Computationally and experimentally calculated energy barriers (0.81 and 0.83 eV respectively) only differ by 2.5%. Similarly, the pre-exponential diffusion coefficients only differ by 0.5 orders of magnitude. Moreover, the diffusivity of H i +1 in amorphous Al 2 O 3 in the 1000–2000 K range is calculated to be 2.53X10 -2 exp(-0.89 eV/k B /T), just one order of magnitude higher than the corresponding value in α-Al 2 O 3 . This suggests that local structural disorder does not significantly affect the energy landscape and diffusion behavior of H i +1 in Al 2 O 3 . Overall, these results show promise for the application of alumina polymorphs as hydrogen permeation barriers.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

A generalized power-law criterion for rocks based on Mohr failure theory

Nonlinearity of Mohr envelopes confirmed by past laboratory test data can be of considerable importance in theoretical interpretation of rock strength behaviors and many practical applications. Current nonlinear failure criteria are introduced with parameters to empirically fit test data, and the physical and mechanical meanings of these fitting parameters are often not known. Predicted strengths by theoretical criteria (e.g., Griffith criterion) generally do not agree well with test data. Considering the nonlinear properties of a Mohr envelope, a nonlinear failure criterion in a generalized power-law form is derived in this paper based on strength magnitude coefficient (a) and failure envelope curvature coefficient (n) which are expressed in terms of uniaxial tensile strength, uniaxial compressive strength, and angle of internal friction. Test data on granite, marble, dolomite, and andesite from published literature were chosen to test the performance of the proposed criterion along with Griffith and Hoek-Brown criteria. Statistical analysis was performed to evaluate the application of different criteria in predicting the strength of rocks. In addition, some Griffith-typed criteria are mathematically shown to be the special cases of the proposed criterion. The proposed and Hoek-Brown criteria have similar power-law forms with a significant difference in the coefficient functions that govern the shear strength magnitude of rock. The different tensile and compressive failure mechanisms are also investigated from the perspective of the relationship between the uniaxial tension Mohr circle and the proposed failure envelope.

58 GEOSCIENCES↗
DOE OSTI2021

A new coupled multiphase flow–finite strain deformation–fault slip framework for induced seismicity

Production of hydrocarbons and water from subsurface reservoirs are known to cause permanent deformation of the reservoir and seismicity along faults both of which are detrimental to sustainable development of natural resources. Most of the prior studies on understanding fluid flow-induced plasticity and seismicity have focused on one or the other phenomenon due to the numerical difficulty associated with simultaneous modeling of the two failure phenomena because stress and deformation evolve non-linearly in both plasticity and seismicity. However, in reservoirs undergoing long-term production, plastic failure can alter the stress paths of points on a fault such that the onset, location, and magnitude of actual seismic events can no longer be predicted by a poroelastic simulation due to inaccurate stress and deformation history. We present a computational framework for coupled multiphase flow, finite strain poroplastic deformation, and dynamic fault slip and use it to understand the impact of plastic deformation on the onset, location, and magnitude of induced fault slip events. We evaluate the impact of plasticity on reservoir pressure, deformation, induced stress, and fault slip by comparing infinitesimal strain elastic and finite strain poro-elastoplastic models. For real-world applications, we consider different scenarios where the reservoir is either mechanically weaker or stronger than the caprock. We analyze the stress evolution as a function of the change in reservoir pressure to understand the role of contrast in reservoir and caprock elastic moduli on geomechanical stability. The results show that the poroplastic reservoir exhibits larger vertical deformation and delayed slip than the poroelastic reservoir after the same amount of oil production. For the same amount of pressure drop, a reservoir with a smaller modulus than the caprock displays a larger vertical displacement and an earlier onset of both plastic failure and fault slip. For a reservoir with a larger modulus than the caprock, vertical displacement is larger on the reservoir top boundary and smaller on the ground surface, and a higher pressure drop is needed to induce plastic failure and fault slip.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Integrating a ponderomotive guiding center algorithm into a quasi-static particle-in-cell code based on azimuthal mode decomposition

High fidelity modeling of plasma based acceleration (PBA) requires the use of three dimensional, fully nonlinear, and kinetic descriptions based on the particle-in-cell (PIC) method. In PBA an intense particle beam or laser (driver) propagates through a tenuous plasma whereby it excites a plasma wave wake. Three-dimensional PIC algorithms based on the quasi-static approximation (QSA) have been successfully applied to efficiently model the interaction between relativistic charged particle beams and plasma. In a QSA PIC algorithm, the plasma response to a charged particle beam or laser driver is calculated based on forces from the driver and self-consistent forces from the QSA form of Maxwell's equations. These fields are then used to advance the charged particle beam or laser forward by a large time step. Since the time step is not limited by the regular Courant-Friedrichs-Lewy (CFL) condition that constrains a standard 3D fully electromagnetic PIC code, a 3D QSA PIC code can achieve orders of magnitude speedup in performance. Recently, a new hybrid QSA PIC algorithm that combines another speedup technique known as an azimuthal Fourier decomposition has been proposed and implemented. This hybrid algorithm decomposes the electromagnetic fields, charge and current density into azimuthal harmonics and only the Fourier coefficients need to be updated, which can reduce the algorithmic complexity of a 3D code to that of a 2D code. Modeling the laser-plasma interaction in a full 3D electromagnetic PIC algorithm is very computationally expensive due the enormous disparity of physical scales to be resolved. In the QSA the laser is modeled using the ponderomotive guiding center (PGC) approach. We describe how to implement a PGC algorithm compatible for the QSA PIC algorithms based on the azimuthal mode expansion. Here this algorithm permits time steps orders of magnitude larger than the cell size and it can be asynchronously parallelized. Details on how this is implemented into the QSA PIC code that utilizes an azimuthal mode expansion, QPAD, are also described. Benchmarks and comparisons between a fully 3D explicit PIC code (OSIRIS), as well as a few examples related to laser wakefield acceleration, are presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Efficient Bayesian inference with latent Hamiltonian neural networks in No-U-Turn Sampling

When sampling for Bayesian inference, one popular approach in the computational field is to use Hamiltonian Monte Carlo (HMC) and specifically the No-U-Turn Sampler (NUTS), which automatically decides the end time of the Hamiltonian trajectory. However, HMC and NUTS can require numerous numerical gradients of the target density and can prove slow in practice when relying on computationally expensive forward models. We propose Latent Hamiltonian neural networks (L-HNNs) with HMC and NUTS for solving Bayesian inference problems. Once trained, L-HNNs do not require numerical gradients of the target density during sampling, and hence numerous evaluations of the forward computational model. Moreover, L-HNNs satisfy important properties such as perfect time reversibility and Hamiltonian conservation, making them well-suited for use within HMC and NUTS because stationarity can be shown. We also propose the integration of L-HNNs in an online error monitoring scheme, in which numerical gradients of the target density are used for a few samples whenever the L-HNNs prediction errors are large. This online error monitor scheme prevents sample degeneracy in regions of low probability density and ensures robust uncertainty quantification. We demonstrate L-HNNs in NUTS with online error monitoring on several analytical examples involving complex, heavy-tailed, and high-local-curvature probability densities. We then demonstrate the applicability of L-HNNs in NUTS to two computational case studies, namely the Allen-Cahn stochastic partial differential equation and an elliptic partial differential equation with 25 and 50 inference parameters, respectively. Overall, the L-HNNs in NUTS with online error monitoring satisfactorily inferred these probability densities. In conclusion, compared to traditional NUTS, L-HNNs in NUTS with online error monitoring required 1–2 orders of magnitude fewer numerical gradients of the target density and improved the effective sample size (ESS) per gradient (which is a measure of both the sampling quality and the computational expense) by an order of magnitude.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Comparison of source-location algorithms for atmospheric samplers

Numerous algorithms have been developed to determine the source characteristics for an atmospheric release, e.g., (Bieringer et al., 2017). This study compares three models that have been applied to the data collected by the International Monitoring System operated by the Comprehensive Nuclear-Test-Ban Treaty Preparatory Commission to estimate source event parameters. Each model uses a different approach to estimate the parameters. A deterministic model uses a possible source region (PSR) approach (Ringbom et al., 2014) that is based on the correlation between predicted and measured sample values. An unnamed model developed at Pacific Northwest National Laboratory (herein called BAYEST) uses a Bayesian formulation (Eslinger et al., 2019; Eslinger et al., 2020; Eslinger and Schrom, 2016). The FREAR model uses a different Bayesian formulation (De Meutter and Hoffman, 2020; De Meutter et al., 2021; Pieter De Meutter, 2021). The performance of the three source-location models is evaluated with 100 synthetic release cases for the single xenon isotope, 133 Xe. The release cases resulted in detections in a fictitious network with 120 noble gas samplers. All three source-location models use the same sampling data. The two Bayesian models yield more accurate location estimates than the deterministic PSR model, with FREAR having slightly better location performance than BAYEST. Samplers with collection periods of 3, 6, 8, 12, and 24-h were used. Results from BAYEST show that location accuracy improves with each reduction in sample collection length. The BAYEST model is slightly better for estimating the start time of the release. The PSR model has about the same spread in start times as the FREAR model, but the PSR results have a better average start time. The Bayesian source-location algorithms give more accurate results than the PSR approach, and provide release magnitude estimates, while the base PSR model does not estimate the release magnitude. This investigation demonstrates that a reasonably dense sampling grid will sometimes yield poor location and time estimates regardless of the model. The poor estimates generally coincide with cases where there is a much larger distance between the release point and the first detecting sampler than the average sampler spacing.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2021

Competitive exchange between divalent metal ions [Cu(II), Zn(II), Ca(II)] and Hg(II) bound to thiols and natural organic matter

Mercuric Hg(II) ion forms exceptionally strong complexes with various organic ligands, particularly thiols and dissolved organic matter (DOM) in natural water. Few studies, however, have experimentally determined whether or not the presence of base cations and transition metal ions, such as Ca(II), Cu(II), and Zn(II), would compete with Hg(II) bound to these ligands, as concentrations of these metal ions are usually orders of magnitude higher than Hg(II) in aquatic systems. Different from previous model predictions, a significant fraction of Hg(II) bound to cysteine (CYS), glutathione (GSH), or DOM was found to be competitively exchanged by Cu(II), but not by Zn(II) or Ca(II). About 20–75% of CYS-bound-Hg(II) [at 2:1 CYS:Hg(II)] and 14–40% of GSH-bound-Hg(II) [at 1:1 GSH:Hg(II)] were exchanged by Cu(II) at concentrations 1–3 orders of magnitude greater than Hg(II). Competitive exchange was also observed between Cu(II) and Hg(II) bound to DOM, albeit to a lower extent, depending on relative abundances of thiol and carboxylate functional groups on DOM and their equilibrium time with Hg(II). When complexed with ethylenediaminetetraacetate (EDTA), most Hg(II) could be exchanged by Cu(II) and Zn(II), as well as Ca(II) at increasing concentrations. Furthermore, these results shed additional light on competitive exchange reactions between Hg(II) and coexisting metal ions and have important implications in Hg(II) chemical speciation and biogeochemical transformation, particularly in contaminated environments containing relatively high concentrations of Hg(II) and metal ions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Epsilon Metal: A Waste Form for Noble Metals from Used Nuclear Fuel

Epsilon metal (e-metal) is the metallic phase that forms as inclusions at the grain boundaries in the UO2 fuel during reactor operation. This metal is composed of Pd, Mo, Rh, Ru, and Tc. These metallic inclusions are insoluble in strong acid and remnants of these metallic inclusions have been found in the UO2 matrix that remains from the natural reactors in Gabon that were active 1.8 billion years ago, therefore e-metal should be an excellent waste form for the immobilization of the long-lived isotopes 107Pd (6.5 ? 106 a) and 99Tc (2.13 ? 105 a), with 99Tc being the isotope of interest for repository performance. Therefore, the chemical durability of this potential waste form is assessed in this study. Typically, corrosion rates for metallic materials are measured electrochemically because they are quick, inexpensive, and can reveal the mechanism by which a metal corrodes, at least initially. However, in a repository the waste form would be subjected to slowly flowing water without an applied electrical potential over long time periods. Therefore, the corrosion rates of e-metal specimens were measured with both electrochemical tests and the single-pass flow-through test (SPFT). Potentiodynamic and potentiostatic polarization results suggest that a thin passive film exists on the alloy surface, which seems to be responsible for its high corrosion resistance. Additionally, X-Ray photoelectron spectroscopic results suggest that Pd oxides are significantly enriched in the passive film Results from the SPFT show that the dissolution rates were weakly dependent on pH. Only Mo and Re were found in solution and were used for the calculation of the dissolution rates. In general, the electrochemically determined corrosion rates agree reasonably well with the initial dissolution rate measured with the SPFT test, but they are about one or two orders of magnitude higher than the steady state rates. The causes for this discrepancy are discussed.

epsilon-metal, nuclear waste form, corrosion, sing↗
DOE OSTI2022

The response of accident tolerant fuel cladding to LOCA burst testing: A comparative study of leading concepts

Accident tolerant fuel claddings seek to improve safety margins during loss-of-coolant accident (LOCA) scenarios by reducing cladding steam oxidation rates, hydrogen production, delaying or mitigating cladding burst, and reducing cladding deformation. After rapid and extensive progress in recent years, Cr coated Zr-based alloys, ferritic FeCrAl, and SiC/SiC ceramic matrix composites (CMCs) have emerged as the most promising cladding concepts, but there are few studies comparing the accident tolerance of all three concepts simultaneously. Here, in this study, LOCA burst testing has been conducted on Cr high impulse power magnetron sputtered Zry-4 (Cr/Zry-4), C26M (FeCrAl), and SiC/SiC SiGA (CMC). Experimental observations indicate Cr/Zry-4 burst at higher temperatures than Zry-4 however the effectiveness of the coating on mitigating ZrO 2 formation via Cr 2 O 3 formation depended on physical proximity to the burst opening. No dramatic improvement in diametrical strain nor opening geometry was observed. C26M burst at higher temperatures than Zry-4 at lower pressures, exhibiting small amounts of strain and smaller openings, and formed slow-growing Al 2 O 3 regardless of location. The CMCs did not burst or show any macroscopic signs of deformation, all the way up to an internal cladding pressure of 15.2 MPa. High temperature expanding Nb plug testing indicated that CMC room temperature strength may persist up to 1650°C, with an approximately 40% loss of strength at 1900°C. CMC steam oxidation was also compared to reference chemically vapor deposited SiC up to 1700°C, and CMCs were found to oxidize at slightly accelerated rates that were still orders of magnitude lower than bare Zry-4. Projected H 2 generation during LOCA burst testing showed that when accounting for cladding burst exposing inner diameters to steam, Cr coated Zry-4 reduced H 2 generation by a factor of two compared with bare Zry-4, while both C26M and SiC/SiC reduced H 2 generation by orders of magnitude relative to bare Zry-4.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Microstructure-process relationships in monolithic U-10Mo fuel foil single-pass rolling: A parametric simulation study

In this work, a previously validated coupling of Kinetic Monte Carlo (KMC) Potts Model and finite element method (FEM) simulations was implemented to investigate the effects of microstructural features in as-cast and homogenized monolithic U-10Mo foils on the emergent microstructure after rolling and reheating. Parameters that could potentially affect recrystallization behavior of the rolled U-10Mo foils were considered: grain size distribution, uranium carbide (UC) size distribution, UC volume fraction, spatial distribution of UC, and rolling reduction magnitude. Grain structure and the magnitude of rolling reduction have the strongest influence on recrystallization kinetics and the fabricated grain size distribution. The UC distribution had only a weak effect on the recrystallization kinetics and final microstructures. While particle-stimulated nucleation (PSN) occurred in simulation more frequently as grain size increased, its incidence did not appear to considerably affect the recrystallization kinetics or grain size distribution.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Cluster dynamics simulations of tritium and helium diffusion in lithium ceramics

Tritium (T) and He diffusion in LiAlO 2 and LiAl 5 O 8 phases influences the performance of tritium producing burnable absorber rods (TPBARs) by affecting the gas release, swelling and thermal conductivity of Li-bearing ceramic pellets. Frenkel pair defects and clusters created by irradiation can attract T and He interstitials and form clusters of the type He i x Li, He i x Al, He i x O, T i x Li, T i x Al, and T i x O, 1 ≤ x ≤ 4 in a Li, Al or O vacancy site (notation denotes x He or T atoms in a 1 Li, 1 Al or 1 O vacant site). The concentration and mobility of each of these clusters collectively contribute to the diffusion of the He and T gases in LiAlO 2 and LiAl 5 O 8 . Here, in this work, free energy cluster dynamics simulations implemented in the Centipede code, are used to obtain the concentration and diffusivities of these clusters which are then used to calculate the total diffusivity of T and He gases in LiAlO 2 and LiAl 5 O 8 . The results show that diffusivity of T is at least one order of magnitude higher in LiAlO 2 as compared to that in LiAl 5 O 8 whereas He diffusion is 2–13 orders of magnitude higher in LiAlO 2 as compared to that in LiAl 5 O 8 . There is a higher concentration of highly diffusive species (T interstitials and T i 03 Li for the case of tritium and He i 01 Li, He i 02 Li, and He i 03 Li for the case of He) in LiAlO 2 than in LiAl 5 O 8 which increase the total diffusion of T and He in LiAlO 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2024

Radiation damage study of POCO ZXF-5Q graphite for neutrino production targets using 4.5 MeV helium ions

To address the challenges of increased beam power and target survivability associated with next-generation particle production beam lines, high dose, high-energy proton beam conditions are simulated using irradiation from low-energy ion beams. A low-energy ion irradiation study of POCO ZXF-5Q graphite under conditions similar to those of the NuMI NT-02 neutrino production target at the Fermi National Accelerator Laboratory is reported. Helium ion irradiation was performed at 100 ∘ C to a maximum damage level of 0.9 displacements per atom (DPA). Irradiation induced hardening, swelling of the irradiated region, inter-plane lattice expansion, and intraplane lattice contraction with increasing ion fluence was observed using micromechanical (nanoindentation, atomic force microscopy) and electron microscopy (high-resolution imaging, selected area diffraction) characterization. Similar changes were also observed in post irradiation examination of the NT-02 target indicating that ion irradiation can be a valuable tool for estimating radiation damage in proton beam targets. Caution must be exercised though, because the hardening, lattice alteration, and swelling occur to different magnitudes for a given damage level. The observed hardening and embrittlement were greater for ion irradiated graphite. For He ion irradiated samples the lattice spacing changes were smaller at low damage levels (78% less expansion and 71% less contraction at 0.1 DPA) and larger at high damage levels (38% more expansion and 5% more contraction at 0.9 DPA) relative to that observed in the NT-02 target. The magnitude of swelling was 8.5× greater under ion irradiation which is influenced by the differing damage gradients and inclusion of implanted He ions in the region of interest.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Building and grid system benefits of demand flexibility and energy efficiency

The global electricity system is undergoing a dramatic transformation. Technology advancements in clean energy generation (e.g., solar photovoltaics and wind) and electrified demand (e.g., transportation and building heating) are providing a cost-competitive pathway to mitigate the urgent and severe consequences associated with a changing climate. In addition, as extreme events continue to increase in frequency and magnitude, solutions are needed for adapting to the current climate. Through energy efficiency and demand flexibility, buildings play a prominent role in meeting these challenges; recently in Joule, a paper by Langevin et al. provides clear insight into the magnitude of buildings’ technical potential.

24 POWER TRANSMISSION AND DISTRIBUTION↗