Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion 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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Estimation of 14 MeV neutron rate from triton burn-up in future W7-X deuterium plasma campaigns

Fast ion confinement is of major importance for the ignition of a burning fusion plasma. In future deuterium plasma campaigns of the Wendelstein 7-X stellarator, W7-X, the amount of triton burn-up is one possible measure for fast ion confinement. A well-established technique to observe triton burn-up is the 14 MeV neutron rate. In this paper, it is estimated whether an existing scintillating fibre neutron detector is also suited to measure triton burn-up in W7-X with sufficient accuracy. In this work, an estimation is presented, which can be applied to any tokamak or stellarator design and is one-dimensional in the minor radius. The inputs are profiles of density, temperature, and differential volume element as well as the triton slowing-down time. The estimation calculates the thermal deuteron fusion rate and the associated deuteron-triton fusion rate; thus, the triton burn-up generated 14 MeV neutron rate. It neither takes triton diffusion nor explicit losses into account. This thermally generated fusion rate is compared to the neutral beam injection heating induced beam-plasma fusion rate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Degradation of Poly- and Perfluoroalkyl Substances (PFAS) in Water via High Power, Energy-Efficient Electron Beam Accelerator

The goal of the 2-year workplan was to see if electron beam (EB) could be used to break down a sub-set of the larger chemical family of per and polyfluoroalkylated substances (PFAS) in an energy efficient and economical manner when compared to conventional water treatment technologies. Year one (Y1) work focused on sample EB treatment work in the Fermi National Accelerator Laboratory’s (FNALs) Accelerator Applications Demonstration and Development (A2D2) EB accelerator. While there are reportedly thousands of types of PFAS, for the point of most of the work herein, a small subset was examined, typically perfluorooctane sulfonate (PFOS) and perfluorooctanoate (PFOA). PFOA and PFOS are two of the most well studied PFAS and are studied for baseline evaluations and are considered most useful. The work from Y1 provided information about the optimal operating parameters and additives to use when treating PFOS and PFOA via EB. The data were then used to see where in a water treatment system an EB accelerator would be best suited to treat PFAS. A conventional water treatment technology, GAC, was then compared to e-beam treatment technology with respect to energy and costs for treatment. In year two (Y2), several conventional e-beam accelerator designs, and FNAL’s developmental compact SRF accelerator design, were evaluated for their suitability in PFAS treatment, from an energy efficiency and cost standpoint. Several EB parameters were evaluated and optimized for the removal of PFOA and PFOS from water at normal pressure and temperature, measured as total PFAS removal. Under the optimized test conditions both PFOA showed complete destruction to inorganic fluoride, and PFOS to inorganic fluoride and sulfate, with mass balance. The effect on PFAS removal relative to solution pH, total EB dose, EB dose rate, dissolved oxygen concentration (DO), temperature, and initial PFAS concentration were evaluated. In general, PFOA was easier to destroy than PFOS. Degradation products, typically observed under less-than-optimal EB conditions, provided insight to degradation mechanisms. Products were identified to rule out possible deleterious biproduct formation. The water radiolysis radical reaction kinetics with PFOS and PFOA were not dependent on the initial concentration over 5-orders of magnitude from 2 μg/L to 20 mg/L. This is thought to be because there was an overabundance of the reactive water radiolysis radicals relative to PFAS molecules and largely attributed to aqueous electrons. The reaction rates appeared to be diffusion limited. Testing at higher concentrations (100-200 mg/L) showed a decrease in removal efficiency, suggesting alternative kinetics, possibly second order rates, at higher concentrations. In all, we successfully defined a set of optimal EB parameters to treat PFOA and PFOS at concentrations of 20 mg/L in water with destruction efficiencies near 100%. We further tested the optimized EB parameters with other types of PFAS, including shorter and longer fluorocarbon chain homologs of PFOA and PFOS, and PFAS with alternative functional groups such as sulfonamides. Based on our results EB can be optimized as an effective destructive technology for removing PFAS from water. The conditions optimized for PFOA and PFOS were less effective with ultra-short fluorocarbon compounds like TFMS, PFES, PFPS and PFBS, and likely require re-optimization of parameters to them. In all, it was determined that from a cost and energy efficiency standpoint, EB would be best applied to waste streams with relatively high concentrations of PFOS and PFOA and is not as cost effective as GAC treatment for removing low concentrations of PFAS from water. Higher concentrations of PFAS can be found in the wastewater of conventional treatment processes such as RO and IE and therefore EB may be used to supplement such treatment technologies. Some real-world IE regeneration wash water and RO reject water containing higher concentrations of PFAS and obtained from pilot scale industrial wastewater treatment system at a fluorochemical manufacturing facility, showed that EB could remove PFAS from such types of wastewaters. The IE regenerant wash water appeared to be the most efficient of the two types of wastewaters tested. However, some further optimization of the EB parameters for the specific PFAS types present in those wastewaters may be required. Also, the effects of co-present TOC and mineral salts should be considered during such optimization efforts. From the experimental Y1 results it was seen that the aqueous electron drives degradation of the PFAS. In a hypothetical water treatment skid using EB for PFAS destruction the parameters of the system should be optimized to promote aqueous electron production. Before EB treatment, the PFAS should be preconcentrated when possible, the pH should be raised to pH 10 or higher to enhance aqueous electron production, and the water should be nitrogen purged to remove dissolved oxygen to minimize aqueous electron scavenging. An excel spreadsheet was created that calculates optimal conditions based on inlet PFAS concentration and desired outlet concentration, by optimizing the accelerator power, dose rate, water treatment rate, pH and dissolved oxygen levels to reach the desired endpoint. Given this information on accelerator operating conditions five different EB accelerator systems were compared. One EB system was a continuous-wave, linear superconducting accelerator being designed at Fermilab. Three other EB systems (IMPELA at 5% and 25% duty factor and the ILU-14) were normal conducting pulsed linear accelerators. The fifth system was an IBA Rhodotron which is a normal conducting, circular, continuous-wave accelerator. The accelerator efficiency (% of the incoming power that is used in water treatment) was the dominating factor in accelerator choice. The radio frequency (RF) power supply and the accelerator design (superconducting versus warm technology) drive the accelerator efficiency. The IBA Rhodotron was seen to be the most energy efficient commercially available technology with a wall-plug (total) power efficiency of 43% at 400 kW. The Fermilab design, with a prototype for a different application currently being fabricated, was the most energy efficient at 55% when driven by a Klystron RF power supply and as high as 77% when powered by a magnetron. As the Fermilab design was the most energy efficient by approximately 10-30%, further design work was done on the accelerator and beam delivery system specific to the destruction of PFAS in water. The Fermilab design is unique from industrial accelerators in that is superconducting. Superconducting technology allows for the acceleration of electrons without losses. The accelerator must be cooled to below the point where it is superconducting and is operated around 4 degrees Kelvin. The bulk of the design work for the accelerator is on making the accelerator as energy efficient as possible so that it does not require liquid helium and can be cooled with conduction cooling via cryocoolers. Final design work resulted in an EB accelerator that would operate at minimally 200 kW and 10 MeV. Prototype construction would cost $\$ $7.8 million dollars when driven by a Klystron power supply. A second version of the same accelerator would cost $\$ $5.5 million dollars when driven by a magnetron that is still under development. The commercially available 300 kW IBA Rhodotron cost was estimated at approximately $\$ $9 million. While it is hard to directly compare, an operational GAC system used by 3M for groundwater treatment capital cost (2022 dollars) was estimated to cost $\$ $3.3 million. While the capital expense of the EB accelerator systems was higher than GAC, the accelerator EB treatment would result in destruction of the PFAS and not just sequestration of PFAS to form a new waste stream that requires further treatment or disposal. The operating cost to destroy the PFAS via 400 kw EB system was less than $\$ $1000/kg of PFAS destroyed when treating at a 20 mg/L PFAS concentration, compared to GAC with operating costs that calculated at $\$ $27,530 per kg of PFAS sequestered when treating 100 μg/L PFOA and PFOS combined concentration.

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

Surface diffusion manifestation in electrodeposition of metal anodes

Metal anode-based battery systems have been deemed indispensable towards energy storage renaissance engendering extensive research into strategies countering dendritic growth of metal electrodeposition. Fundamentally, the morphological evolution of a material is uniquely characterized by the heights of its self-diffusion barrier across multiple pathways. Herein, based on a coarse-grained kinetic Monte Carlo method, we derive insights into the nucleation and growth of metallic electrodeposits in liquid electrolytes, governed by surface self-diffusion characteristics cognizant of the diverse diffusion routes including terrace, away from step and interlayer pathways. We deconvolve the roles played by each of these surface diffusion mechanisms in conjunction with the electrochemical reaction rate on the deposition morphology regime (film vs. mossy vs. fractal). We identify interlayer diffusion as the predominant morphology-determining mechanism; dendrite-free deposition even at moderate current rates constrains this diffusion barrier to an upper limit. Additionally, we highlight subtle features amidst the realm of the morphological growth assortment that connect to the cell's electrochemical performance. Finally, we delineate morphological features of Li, Na, Mg and Al based on their respective surface diffusion barriers and applied overpotentials, and provide a baseline for the interpretation of experimental observations. This fundamental study sheds light on the mesoscale underpinnings of morphological variances in mono-valent and multi-valent metal electrodeposition.

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

Structure and Electrochemical Properties of Bronze Phase Materials Containing Two Transition Metals

Bronze phase transition-metal oxides have recently attracted attention as high-rate lithium-ion battery anode materials. Their crystal structures are distinguished by large tunnels and an open framework, facilitating lithium-ion diffusion and high-rate charge–discharge properties. The presence of two transition metals also offers a route to achieve high energy density from multielectron redox. In this paper, we report the chemistry, structure, and electrochemical properties of two different bronze phase compositions having the same stoichiometry: W 3 Nb 2 O 14 and Mo 3 Nb 2 O 14 . These materials provide insight into how the transition metals affect the electrochemical behavior and structural stability of bronze phase materials. Mo 3 Nb 2 O 14 exhibits greater than 1 electron redox per transition metal leading to lithium capacities above 200 mAh g –1 at C/2 but is unable to maintain this high capacity at high rates due to incomplete Mo redox reactions. In contrast, W 3 Nb 2 O14 exhibits reversible redox reactions and retains its open structure on cycling. This study highlights the potential of bronze phase materials containing two transition metals to exhibit fast charging properties with a high energy density.

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

Long-range exciton transport and slow annihilation in two-dimensional hybrid perovskites

Two-dimensional hybrid organic-inorganic perovskites with strongly bound excitons and tunable structures are desirable for optoelectronic applications. Exciton transport and annihilation are two key processes in determining device efficiencies; however, a thorough understanding of these processes is hindered by that annihilation rates are often convoluted with exciton diffusion constants. Here we employ transient absorption microscopy to disentangle quantum-well-thickness-dependent exciton diffusion and annihilation in twodimensional perovskites, unraveling the key role of electron-hole interactions and dielectric screening. The exciton diffusion constant is found to increase with quantum-well thickness, ranging from 0.06 ± 0.03 to 0.34 ± 0.03 cm 2 s –1 , which leads to long-range exciton diffusion over hundreds of nanometers. The exciton annihilation rates are more than one order of magnitude lower than those found in the monolayers of transition metal dichalcogenides. The combination of long-range exciton transport and slow annihilation highlights the unique attributes of two-dimensional perovskites as an exciting class of optoelectronic materials.

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

Nuclear Magnetic Resonance Dynamics of LiTFSI–Pyrazole Eutectic Solvents

Deep Eutectic Solvents (DESs) have emerged as promising candidates to replace conventional organic solvents in various technological applications due to their low vapor pressure, non-flammability, and ease of preparation at low costs. In particular, Type IV DESs, which are composed of metal salts and hydrogen bond donors, are possible replacements for lithium-ion battery electrolytes. In this study, we investigate the molecular dynamics of solvents of lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) and pyrazole (PYR) at varying LiTFSI:PYR molar ratios (1:2, 1:3, 1:4, 1:5) using Nuclear Magnetic Resonance Dispersion (NMRD) and Pulsed Field Gradient (PFG) Nuclear Magnetic Resonance (NMR). PFG NMR reveals composition-dependent diffusion trends, while NMRD provides molecular-level insights into the longitudinal relaxation rate (R 1 = 1/T 1 ). Notably, the LiTFSI:PYR (1:2) sample shows distinct behavior across both techniques, exhibiting enhanced relaxation rates and lower self-diffusion for 1 H compared to the other nuclei ( 19 F and 7 Li), suggestive of stronger and more efficient Li + –pyrazole interactions, as confirmed by the modeling of the relaxation profiles. Our study advances understanding of ion dynamics in azole-based eutectic solvents, supporting their potential use in safer battery electrolytes.

FFC-NMR↗
DOE OSTI2021

Unlocking bulk and surface oxygen transport properties of mixed oxide-ion and electron conducting membranes with combined oxygen permeation cell and oxygen probe method

Surface exchange kinetics and bulk diffusion of oxygen are of paramount importance to the activity of oxygen electrocatalysis and performance of electrochemical devices such as fuel cell, metal-air batteries, and oxygen separation membranes. Conventional approaches to obtaining these transport properties are often limited to single property under a specific non-operation related condition. Here we use a combined oxygen permeation cell and oxygen probe methodology to simultaneously attain rates of oxygen surface exchange and bulk conductivity/chemical diffusivity of three representative mixed oxide-ion and electron conductors, namely SrCo 0.9 Ta 0.1 O 3-δ (SCT), La 0.6 Sr 0.4 CoO 3-δ (LSC) and La 0.6 Sr 0.4 FeO 3-δ (LSF), operated under a steady-state oxygen flux. The results explicitly show that SCT exhibit the highest oxide-ion conductivity/chemical diffusivity, fastest rates of surface oxygen exchange kinetics, thus promising to be the best oxygen electrocatalyst. Here, we have also mapped out the distribution of oxygen chemical potential gradient across the membranes and applied B-transport number concept to illustrate the rate-limiting steps in the overall oxygen permeation process.

42 ENGINEERING↗
DOE OSTI2020

Computational modeling of charge hopping dynamics along a disordered one-dimensional wire with energy gradients in quantum environments

This computational study investigates the effects of energy gradients on charge hopping dynamics along a one-dimensional chain of discrete sites coupled to quantum bath, which is modeled at the level of Pauli master equation (PME). This study also assesses the performance of different approximations for the hopping rates. Three different methods for solving the PME, a fourth order Runge–Kutta method, numerical diagonalization of the rate matrix followed by analytic propagation, and kinetic Monte Carlo simulation method, are tested and confirmed to produce virtually identical values of time dependent mean square displacement, diffusion constant, and mobility. Five different rate expressions, exact numerical evaluation of Fermi’s Golden Rule (FGR) rate, stationary phase interpolation (SPI) approximation, semiclassical approximation, classical Marcus rate, and Miller–Abrahams rate, are tested to help understand the effects of approximations in representing quantum environments in the presence of energy gradients. The results based on direct numerical evaluation of FGR rate exhibit transition from diffusive to non-diffusive behavior with the increase in the gradient and show that the charge transport in the quantum bath is more sensitive to the magnitude of the gradient and the disorder than in the classical bath. Among all the four approximations for the hopping rates, the SPI approximation is confirmed to work best overall. A comparison of two different methods to calculate the mobility identifies drift motion of the population distribution as the major source of non-diffusive behavior and provides more reliable information on the contribution of quantum bath.

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

Transport Processes in a Li-ion Cell during an Internal Short-Circuit

Internal short-circuit in a lithium-ion cell causes an abrupt increase in cell temperature and triggers subsequent thermal runaway. In this work, we present a detailed electrochemical-thermal model to investigate the physical behavior during an internal short-circuit. Simulations at wide range of heat transfer coefficients and short-circuit resistances are conducted to illustrate electrochemical and thermal behavior under a wide range of conditions. The Joule heating at the shorted region promotes electrochemical reactions nearby, causing in-plane non-uniformity of electrolyte and active material transport. Furthermore, it is found that diffusion in solid active materials plays a significant role at very high shorting currents (~20 C), as electrochemical reactions rate are being controlled progressively more by availability of Li + at the interface, due to limitations in diffusion through the active material with increasing discharge rates. This diffusion limitation causes a drop in available energy, and subsequently a decrease in cell equilibrium potential and the heat generation rate at the location of the short. On the other hand, rapid depletion of lithium concentration in the electrolyte and accumulation on the electrode surface results in highly non-uniform transport properties resulting in higher heat generation rates. Hence, the heating regime shifts from "local heating" to "global heating". Based on the findings, important design parameters for battery safety are discussed.

25 ENERGY STORAGE↗
DOE OSTI2021

Analytical modeling for redox flow battery design

Deeper market penetration of redox flow batteries requires optimization of the cell performance. Though important for optimization, detailed analytical solutions have not been developed for electrolyte flow, mass and charge transport, and reaction kinetics within redox flow batteries. To this end, here we present analytical solutions to active species concentration and over-potential based on advection-diffusion transport for ions and Bulter-Volmer model for interface reaction kinetics. The solutions were validated with results from a finite element model. These solutions were then applied to investigate the relationship between over-potential and state of charge, current density, reaction rate constant, flow velocity, diffusivity, total active species concentration, and electrode structure. Explicit formulas were identified for minimum activation over-potential and limiting current density as well as their dependence on electrolyte properties, operation conditions, and electrode structure. With our new mathematical formulas, this work provides a theoretical framework for flow battery design.

25 ENERGY STORAGE↗
DOE OSTI2023

Mass transfer in catalytic depolymerization: External effectiveness factors and serendipitous processivity in stagnant and stirred melts

Several heterogeneous catalysts are being developed to recycle plastics. Most operate in viscous polymer melts, where external mass transfer effects could limit the supply of co-reactants to active sites. External mass transfer can also impede the diffusion of long chain products away from the catalyst after each cut. Product egress limitations could potentially confer unintentional processivity to catalyst operation, i.e. a tendency for the catalyst to repeatedly cut the same chain after an initial encounter. We formulate reaction–diffusion equations to quantify mass transfer effects on the co-reactant transport to the catalyst and the degree of serendipitous processivity. Results are developed for catalysts in stagnant or stirred melts, with simple expressions involving Damkohler, Peclet, and Sherwood numbers, i.e. dimensionless combinations of rate constants, catalyst particle size, polymer diffusivities, and shear rates (where applicable). In conclusion, we estimate the impact of these effects for a spherical core–shell catalyst.

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

Visualization and analysis of coupling between plasmas self-organization and plasma-induced fluid circulation in 1 atm DC glows with liquid anode

The physical processes prevailing when plasma contacts liquid water are poorly understood; however, it is this very interaction that is the basis of new, exciting technologies for water purification, the treatment of cancer and disease and the production of new, high value chemical products. A chief impediment to widespread application of plasmas for the aforementioned applications is scale up and cost. Here we investigate plasma induced flows which facilitate transport of reactivity from the plasma to liquid water at rates much faster than diffusion. We find that this plasma induced motion varies depending on whether or not the plasma is self organized. The insight from this work paves the way to a better understanding of transport and ultimately informing plasma treatment technologies on how to best optimize the rate at which processes such as plasma based water purification occurs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Modeling mesoscale fission gas behavior in UO2 by directly coupling the phase field method to spatially resolved cluster dynamics

Abstract Fission gas release within uranium dioxide nuclear fuel occurs as gas atoms diffuse through grains and arrive at grain boundary (GB) bubbles; these GB bubbles grow and interconnect with grain edge bubbles; and grain edge tunnels grow and connect to free surfaces. In this study, a hybrid multi-scale/multi-physics simulation approach is presented to investigate these mechanisms of fission gas release at the mesoscale. In this approach, fission gas production, diffusion, clustering to form intragranular bubbles, and re-solution within grains are included using spatially resolved cluster dynamics in the Xolotl code. GB migration and intergranular bubble growth and coalescence are included using the phase field method in the MARMOT code. This hybrid model couples Xolotl to MARMOT using the MultiApp and Transfer systems in the MOOSE framework, with Xolotl passing the arrival rate of gas atoms at GBs and intergranular bubble surfaces to MARMOT and MARMOT passing evolved GBs and bubble surface positions to Xolotl. The coupled approach performs well on the two-dimensional simulations performed in this work, producing similar results to the standard phase field model when Xolotl does not include fission gas clustering or re-solution. The hybrid model performs well computationally, with a negligible cost of coupling Xolotl and MARMOT and good parallel scalability. The hybrid model predicts that intragranular fission gas clustering and bubble formation results in up to 70% of the fission gas being trapped within grains, causing the increase in the intergranular bubble fraction to slow by a factor of six. Re-solution has a small impact on the fission gas behavior at 1800 K but it has a much larger impact at 1000 K, resulting in a twenty-times increase in the concentration of single gas atoms within grains. Due to the low diffusion rate, this increase in mobile gas atoms only results in a small acceleration in the growth of the intergranular bubble fraction. Finally, the hybrid model accounts for migrating GBs sweeping up gas atoms. This results in faster intergranular bubble growth with smaller initial grain sizes, since the additional GB migration results in more immobile gas clusters reaching GBs.

Kim, Dong-Uk↗
DOE OSTI2022

Mechanistic calculation of the effective silver diffusion coefficient in polycrystalline silicon carbide: Application to silver release in AGR-1 TRISO particles

The silicon carbide (SiC) layer in tristructural isotropic (TRISO) fuel particles serves as a barrier to prevent the escape of fission products produced and not retained in the fuel kernel. The release of silver (Ag) is a concern due to the long half-life of the 110m Ag isotope. However, accurately determining the fission gas release rate requires knowing the diffusion coefficient through the SiC layer. In this study, we leverage atomistic calculations of Ag diffusivity in SiC bulk and grain boundaries (GBs) to develop a mesoscale effective Ag diffusion coefficient (D eff ) in SiC. Since GBs serve as pathways for Ag diffusion, D eff is defined as a function of temperature and microstructure variables. In particular, the size of SiC grains in the direction perpendicular to diffusion is shown to significantly affect Ag diffusion. The prediction of the mechanistic, mesoscale approach falls within one order of magnitude of empirical values. The temperature and microstructure-dependent effective Ag diffusivity in SiC is implemented in the fuel performance code Bison with a correction factor to predict Ag release from AGR-1 TRISO fuel particles. We hereby quantify the impact of SiC grain size on Ag release and improve Bison’s predictions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Complete quasilinear model for the acceleration-driven lower hybrid drift instability and a computational assessment of its validity

A complete quasilinear model is derived for the electrostatic acceleration-driven lower hybrid drift instability in a uniform two-species low-beta plasma in which current is perpendicular to the background magnetic field. The model consists of coupled nonlinear velocity space diffusion equations for the volume-averaged ion and electron distribution functions. Each species' diffusion coefficient depends on a time-evolving spectral density of the electric-field energy per unit volume and a time-evolving dispersion relation. The dispersion relation is expressed analytically in integral form without the use of asymptotic limits and applies to arbitrary distribution functions, so long as they can be expressed as a function of one velocity coordinate, e.g., f⁡(vy) or f⁡(v⊥). The quasilinear model conserves energy and is complete in that it fully describes the evolution of the distribution functions, including resonant and nonresonant particle-wave interactions, while accounting for distribution-function-dependent mixed-complex frequencies. Further, the quasilinear diffusion model is solved numerically and self-consistently using a Crank-Nicolson temporal discretization and a second-order finite-volume velocity-space discretization. Numerical solutions are compared to nonlinear fourth-order accurate continuum kinetic Vlasov-Poisson simulations. Evolution of electric-field energy, growth rates, distribution functions, and diffusion coefficients are shown to be in agreement with Vlasov simulations. The quasilinear model is shown to predict anomalous transport terms, like resistivity and heating, to within a factor of order unity. Discrepancies between the quasilinear model and Vlasov simulations are assessed and attributed primarily to lack of damping in the quasilinear description and to the use of unperturbed-orbit susceptibilities in the linear theory dispersion relation. The results illuminate the predictive accuracy of the quasilinear model, place approximate bounds on its validity, and provide much needed vetting of quasilinear theory's ability to predict the nonlinear state of a microturbulent plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Hybrid quantum-classical simulations of magic angle spinning dynamic nuclear polarization in very large spin systems

Solid-state nuclear magnetic resonance can be enhanced using unpaired electron spins with a method known as dynamic nuclear polarization (DNP). Fundamentally, DNP involves ensembles of thousands of spins, a scale that is difficult to match computationally. This scale prevents us from gaining a complete understanding of the spin dynamics and applying simulations to design sample formulations. We recently developed an ab initio model capable of calculating DNP enhancements in systems of up to ~1000 nuclei; however, this scale is insufficient to accurately simulate the dependence of DNP enhancements on radical concentration or magic angle spinning (MAS) frequency. We build on this work by using ab initio simulations to train a hybrid model that makes use of a rate matrix to treat nuclear spin diffusion. We show that this model can reproduce the MAS rate and concentration dependence of DNP enhancements and build-up time constants. Finally, we then apply it to predict the DNP enhancements in core–shell metal-organic-framework nanoparticles and reveal new insights into the composition of the particles’ shells.

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

Evaluation of Oxygen Consumption as a Sensitive Measure of Electrical Cable Polymer Insulation Degradation

Polymers used as dielectric materials in nuclear grade electrical cables, such as cross-linked polyethylene (XLPE) and ethylene-propylene rubber (EPR), are long-lived and durable. Estimating the decades-long useful service life of polymers experimentally in a feasible timeframe of weeks or months requires greatly accelerated aging in the laboratory at elevated temperatures and/or gamma dose rates. Many limitations of accelerated aging, such as due to artefacts of the process including diffusion limited oxidation and dose rate effects, have been identified that can make accurate lifetime prediction at service conditions challenging. This can be due to differences in degradation mechanisms between rapid and extended aging or to the presence of induction periods during which early aging is undetectable using common metrics of polymer insulation characterization such as tensile elongation at break. Conceptually, an ideal measure of aging would track clearly and predictably with material changes throughout the lifecycle of the evolving material state from early aging through end of useful life. With such a metric accelerated aging could be performed at much milder conditions, conditions closer to those found in long term service, and end of life calculated from the early data. Oxygen consumption has previously been proposed as a sensitive measure of insulation polymer oxidation that might meet such desirable characteristics. In this work, nuclear cable insulation samples of the two most common types, XLPE and EPR, from two of the most sourced manufacturers, RSCC and The Okonite Company, were subjected to gamma irradiation at dose rates of 100, 200, and 300 Gy/hr for up to 42 days at 26°C. The aged samples were used to evaluate oxygen consumption as a sensitive measure of aging for these materials and these conditions. For comparison and validation, the aged samples were also characterized using the more conventional tests of mass change, density, carbonyl index, and yellowness index. For both materials, oxygen consumption during irradiation was found to track consistently with applied dose—oxygen consumption rate was linear with time. Mass change and carbonyl index were also clearly proportional to exposure dose. Yellowness index, a surface sensitive measure, was observed to increase linearly with dose at low exposure levels before plateauing. Sample density was not observed to change under the conditions explored. Oxygen consumption was found to be a particularly sensitive metric to measure the degradation of polymeric cable insulation exposed to radiation. Oxygen consumption was compared to Arrhenius and time-temperature-superposition methods for activation energy calculation using literature date and found to be a useful and viable alternative for lifetime prediction. The techniques developed in this work are anticipated to be valuable in the estimation of insulation aging rates under mild conditions approaching service conditions. Direction comparison of useful life prediction by oxygen consumption with prediction via Arrhenius methodology for thermal aging and existing methods for radiation and combined thermal/radiation aging would serve to further validate this promising approach and refine our understanding of material-specific aging rates.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Antiphase boundary migration as a diffusion mechanism in a P3 sodium layered oxide

The electrode materials of Na-ion batteries must exhibit fast ion transport, a property that is often negatively impacted when mobile cations order. Many layered Na (and K) intercalation compounds adopt the P3 crystal structure, which tends to stabilize cation-vacancy orderings over a wide range of Na concentrations. The ordered phases have been predicted to consist of domains that are periodically separated by antiphase boundaries (APBs) to accommodate variations in Na concentration. The mechanisms and rates with which Na diffuses within such ordered phases of P3 are not understood. We have conducted a first-principles study of Na diffusion in the P3 structure, using Na x CoO 2 as a model system. We find that most simple nearest-neighbor hops in Na-vacancy ordered phases in P3 result in unstable configurations and are therefore forbidden. We identify an alternative mechanism of Na transport that is mediated by APB migration via the formation and expansion of kinks along domain boundaries. The limiting kinetic barriers for this mechanism are quite low (0.03 to 0.30 eV). Vacancy defects are not required, and furthermore have formation energies that are higher than those of the APB kink defects. Furthermore, our results suggest that APB migration is a fundamental diffusion mechanism in the P3 structure.

36 MATERIALS SCIENCE↗