Search NASA⌕ Search

SEARCH · Search NASA

Results for “effective diffusivity”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Ni coarsening under humid atmosphere in the electrode of solid oxide cells: A combined study of density-functional theory and phase-field modeling

Ni coarsening in solid oxide cell (SOC) electrodes is known to be significantly faster under a humid atmosphere. The underlying mechanisms, though, have not been fully understood. In this work, we examine the surface diffusion of Ni(OH) x by a combination of density-functional theory and phase-field modeling. Here, density-functional theory is used to evaluate the adsorption and diffusion of the Ni(OH) x species on Ni (111) surface, and the results are used to obtain the effective surface diffusivity of Ni versus steam and hydrogen gas partial pressures. This diffusivity is used as an input for a phase-field model to investigate the Ni coarsening in SOC electrodes. It is found that Ni(OH) x formed through chemical reaction cannot accelerate coarsening unless an extremely high steam to hydrogen ratio is reached. However, if Ni(OH) x is formed through electrochemical reactions near triple phase boundaries (TPBs), surface diffusion of Ni(OH) x may cause faster Ni coarsening in fuel cell mode under a large overpotential. Specifically, surface diffusion of Ni(OH) may be comparable to or faster than that of Ni under an overpotential that is large but still possible under fuel cell operating conditions.

25 ENERGY STORAGE↗
DOE OSTI2023

Methane mass transfer in mesoporous silica saturated with liquid hydrocarbons

Mass transfer across gas/liquid interfaces plays a central role in many industrial applications. In particular, gas dissolution and diffusion in liquid hydrocarbon mixtures, confined in nanometer-sized pores, is an essential mechanism during enhanced oil recovery (EOR) from unconventional formations. In this work, we have measured methane (C 1 ) diffusion in n-decane (C 10 ), n-hexadecane (C 16 ), and mixtures of C 10 + C 16 in a mesoporous material with an average pore size of 4 nm at 50 °C and ~ 8 MPa of gas pressure. A key conclusion of this work is that the diffusivities, measured in the bulk phase, for the relevant binary systems, are sufficient to predict the diffusion behavior of the corresponding multicomponent systems in the porous medium by using Wilke’s equation for the evaluation of effective component diffusivities, combined with an accurate equation of state (EOS) representation of the (bulk) phase behavior. This observation can facilitate accurate prediction of recovery processes in unconventional formations and, potentially, guide other applications that entail gas-liquid interface mass transfer in mesoporous materials.

42 ENGINEERING↗
DOE OSTI2023

Evolving radon diffusion through earthen barriers at uranium waste disposal sites

Field measurements of Rn-222 fluxes from the tops and bottoms of compacted clay radon barriers were used to calculate effective Rn diffusion coefficients (D Rn ) at four uranium waste disposal sites in the western United States to assess cover performance after more than 20 years of service. Values of D Rn ranged from 7.4 × 10 -7 to 6.0 × 10 -9 m 2 /s, averaging 1.42 × 10 -7 . Water saturation (S W ) from soil cores indicated that there was relatively little control of D Rn by S W , especially at higher moisture levels, in contrast to estimates from most steady-state diffusion models. Further, this is attributed to preferential pathways intrinsic to construction of the barriers or to natural process that have developed over time including desiccation cracks, root channels, and insect burrows in the engineered earthen barriers. A modification to some models in which fast and slow pathway D Rn values are partitioned appears to give a good representation of the data; 4% of the fast pathway was needed to fit the data regression. For locations with high S w and highest D Rn (and fluxes) at each site, the proportion of fast pathway ranged from 1.7% to 34%, but for many locations with lower fluxes, little if any fast pathway was needed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Curvature effects in the reaction-diffusion system governing frontal polymerization

Frontal polymerization is a fast and energy-efficient curing process of thermoset polymers involving the self-sustained propagation of a sharp front that transforms a liquid monomer or gel into a polymer. The temperature, speed, and geometry of the polymerization front are closely intertwined quantities that can affect the curing process and the mechanical properties of the polymer. This work investigates the effect of the front curvature on the front propagation speed based on a reaction-diffusion model. The proposed analytical relations between front speed and curvature are verified by finite element simulations and validated through experiments. The models show that a concave front has a higher front temperature and speed than its convex counterpart. When the curvature of a convex front reaches a critical value, the reaction can no longer be sustained, and the front is quenched. The relation between the speed of a steadily propagating front and its geometry is derived, showing good agreement with experimental observations. This study also provides analytical estimates of the thickness of the polymerization front and insight into the relationship between the front curvature and the kinetics of the exothermic reaction.

42 ENGINEERING↗
DOE OSTI2022

Flexible machine-learning interatomic potential for simulating structural disordering behavior of Li 7 La 3 Zr 2 O 12 solid electrolytes

Batteries based on solid-state electrolytes, including Li 7 La 3 Zr 2 O 12 (LLZO), promise improved safety and increased energy density; however, atomic disorder at grain boundaries and phase boundaries can severely deteriorate their performance. Machine-learning (ML) interatomic potentials offer a uniquely compelling solution for simulating chemical processes, rare events, and phase transitions associated with these complex interfaces by mixing high scalability with quantum-level accuracy, provided that they can be trained to properly address atomic disorder. To this end, we report the construction and validation of an ML potential that is specifically designed to simulate crystalline, disordered, and amorphous LLZO systems across a wide range of conditions. The ML model is based on a neural network algorithm and is trained using ab initio data. Performance tests prove that the developed ML potential can predict accurate structural and vibrational characteristics, elastic properties, and Li diffusivity of LLZO comparable to ab initio simulations. As a demonstration of its applicability to larger systems, we show that the potential can correctly capture grain boundary effects on diffusivity, as well as the thermal transition behavior of LLZO. Here these examples show that the ML potential enables simulations of transitions between well-defined and disordered structures with quantum-level accuracy at speeds thousands of times faster than ab initio methods.

25 ENERGY STORAGE↗
DOE OSTI2024

Characterization of electromigration-induced short-range stress development in Al(0.25 at. % Cu) conductor line

Scanning x-ray microbeam topography and fluorescence experiments were conducted in situ to study the electromigration behavior of a 0.5 μm thick, 10 μm wide, and 200 μm long Al(0.25 at. % Cu) conductor line with 1.5 μm-thick SiO2 passivation on a single crystal Si substrate. The strain sensitivity of x-ray topography measurement allowed detailed examination of the electromigration-induced stress distribution and evolution in the conductor line in response to the depletion of Cu solute early in the electromigration process. Upon electromigration at 0.4 MA/cm2 and 303 °C, a short-range stress gradient was quickly induced by Al migration in the Cu-depleted cathode region to counteract further Al flow. The stress gradient was fully developed during the 5.3 h incubation time, extending over the critical Blech length of about 66 μm from the cathode end. Plastic deformation then occurred at the downstream end of the Cu-depleted region. The preferential electromigration of Cu did not cause detectable stress change outside the Cu-depleted region, except for the significant stress development from the Al2Cu precipitation at the anode end which appeared to initiate the fracture in the passivation. Preliminary finite difference modeling was undertaken to simulate the experimental observations, from which important parameters dictating electromigration in Al(Cu) line were extracted: an apparent effective valence of −5.6 and −1.9 for Cu and Al in Al(Cu), respectively, and a critical Cu concentration of 0.16 at. % above which Al grain boundary diffusion is effectively blocked.

Physics↗
DOE OSTI2024

Onset of kinetic effects on Rayleigh–Taylor instability: Advective–diffusive asymmetry

In nature and engineering applications, the Rayleigh–Taylor instability (RTI) occurs over a wide range of Atwood, Reynolds, Mach, and Knudsen numbers. At low Atwood, Mach, and Knudsen numbers, the classic advective instability causes quasi-symmetric bubble and spike growth on the two sides of the interface. However, recent findings suggest that at high degrees of rarefaction, advective effects are suppressed and molecular diffusion leads to planar growth of the density fronts on either side of the interface. This study aims to investigate the flow physics of the transition from advective to diffusive behavior, focusing on the onset of the kinetic effects. Using the gas kinetic methodology, RTI is simulated over a range of Knudsen and Mach numbers in the transition regime. The simulation results reveal the various stages of transformation from advective instability to diffusive transport. For the first time, the study demonstrates the existence of a Knudsen–Mach parameter regime where the bubble side exhibits advective instability, while the other side shows a planar density front due to molecular diffusion, rather than the canonical advective spike shape. The dominance of different mechanisms on the two sides of the interface leads to the advective–diffusive asymmetry. In conclusion, the findings of this study can lead to a more comprehensive understanding of RTI over a wide range of Mach and Knudsen numbers.

42 ENGINEERING↗
DOE OSTI2020

Pattern formation in a coupled membrane-bulk reaction-diffusion model for intracellular polarization and oscillations

Reaction-diffusion systems have been widely used to study spatio-temporal phenomena in cell biology, such as cell polarization. Coupled bulk-surface models naturally include compartmentalization of cytosolic and membrane-bound polarity molecules. Here we study the distribution of the polarity protein Cdc42 in a mass-conserved membrane-bulk model, and explore the effects of diffusion and spatial dimensionality on spatio-temporal pattern formation. We first analyze a one-dimensional (1-D) model for Cdc42 oscillations in fission yeast, consisting of two diffusion equations in the bulk domain coupled to nonlinear ODEs for binding kinetics at each end of the cell. In 1-D, our analysis reveals the existence of symmetric and asymmetric steady states, as well as anti-phase relaxation oscillations typical of slow-fast systems. We then extend our analysis to a two-dimensional (2-D) model with circular bulk geometry, for which species can either diffuse inside the cell or become bound to the membrane and undergo a nonlinear reaction-diffusion process. We also consider a nonlocal system of PDEs approximating the dynamics of the 2-D membrane-bulk model in the limit of fast bulk diffusion. In all three model variants we find that mass conservation selects perturbations of spatial modes that simply redistribute mass. In 1-D, only anti-phase oscillations between the two ends of the cell can occur, and in-phase oscillations are excluded. In higher dimensions, no radially symmetric oscillations are observed. Instead, the only instabilities are symmetry-breaking, either corresponding to stationary Turing instabilities, leading to the formation of stationary patterns, or to oscillatory Turing instabilities, leading to traveling and standing waves. Codimension-two Bogdanov—Takens bifurcations occur when the two distinct instabilities coincide, causing traveling waves to slow down and to eventually become stationary patterns. Our work clarifies the effect of geometry and dimensionality on behaviors observed in ma.ss-conserved cell polarity models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The response of the Amazon ecosystem to the photosynthetically active radiation fields: integrating impacts of biomass burning aerosol and clouds in the NASA GEOS Earth system model

The Amazon experiences fires every year, and the resulting biomass burning aerosols, together with cloud particles, influence the penetration of sunlight through the atmosphere, increasing the ratio of diffuse to direct photosynthetically active radiation (PAR) reaching the vegetation canopy and thereby potentially increasing ecosystem productivity. In this study, we use the NASA Goddard Earth Observing System (GEOS) model with coupled aerosol, cloud, radiation, and ecosystem modules to investigate the impact of Amazon biomass burning aerosols on ecosystem productivity, as well as the role of the Amazon's clouds in tempering this impact. The study focuses on a 7-year period (2010–2016) during which the Amazon experienced a variety of dynamic environments (e.g., La Niña, normal years, and El Niño). The direct radiative impact of biomass burning aerosols on ecosystem productivity – called here the aerosol diffuse radiation fertilization effect – is found to increase Amazonian gross primary production (GPP) by 2.6 % via a 3.8 % increase in diffuse PAR (DFPAR) despite a 5.4 % decrease in direct PAR (DRPAR) on multiyear average during burning seasons. On a monthly basis, this increase in GPP can be as large as 9.9 % (occurring in August 2010). Consequently, the net primary production (NPP) in the Amazon is increased by 1.5 %, or ~92 Tg C yr -1 – equivalent to ~37 % of the average carbon lost due to Amazon fires over the 7 years considered. Clouds, however, strongly regulate the effectiveness of the aerosol diffuse radiation fertilization effect. The efficiency of this fertilization effect is the highest in cloud-free conditions and linearly decreases with increasing cloud amount until the cloud fraction reaches ~0.8, at which point the aerosol-influenced light changes from being a stimulator to an inhibitor of plant growth. Nevertheless, interannual changes in the overall strength of the aerosol diffuse radiation fertilization effect are primarily controlled by the large interannual changes in biomass burning aerosols rather than by changes in cloudiness during the studied period.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Polyelectrolyte Surface Diffusion in a Nanoslit Geometry

The surface diffusion of poly-l-lysine (PLL) in a planar nanoslit was studied using convex lens-induced confinement (CLiC) single-molecule tracking microscopy. Here, three surface chemistries were employed to understand the interplay of electrostatic and short-range interactions: an amine-functionalized silica surface, an oligo(ethylene oxide) (OEG)-modified surface, and a 1:1 mixture of the two ligands. Effective surface diffusion coefficients increased rapidly with slit height until saturating for slit heights <30 nm. While diffusion at a semi-infinite interface was significantly faster for OEG surfaces, the diffusion coefficient increased most rapidly with slit height for amine-functionalized surfaces, resulting in surface diffusion within very thin slits being nearly independent of surface chemistry. Intermittent random walks were simulated within a planar slit geometry, using experimentally measured parameters obtained from diffusion at a single interface to account for the characteristic short-range interactions between PLL and each surface chemistry, and were in good agreement with experimental measurements.

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

K-space algorithmic reconstruction (KAREN): a robust statistical methodology to separate Bragg and diffuse scattering

Diffuse scattering occurring in the Bragg diffraction pattern of a long-range-ordered structure represents local deviation from the governing regular lattice. However, interpreting the real-space structure from the diffraction pattern presents a significant challenge because of the dramatic difference in intensity between the Bragg and diffuse components of the total scattering function. In contrast to the sharp Bragg diffraction, the diffuse signal has generally been considered to be a weak expansive or continuous background signal. In this paper, using 1D and 2D models, it is demonstrated that diffuse scattering in fact consists of a complex array of high-frequency features that must not be averaged into a low-frequency background signal. To evaluate the actual diffuse scattering effectively, an algorithm has been developed that uses robust statistics and traditional signal processing techniques to identify Bragg peaks as signal outliers which can be removed from the overall scattering data and then replaced by statistically valid fill values. This method, described as a 'K-space algorithmic reconstruction' (KAREN), can identify Bragg reflections independent of prior knowledge of a system's unit cell. KAREN does not alter any data other than that in the immediate vicinity of the Bragg reflections, and reconstructs the diffuse component surrounding the Bragg peaks without introducing discontinuities which induce Fourier ripples or artifacts from underfilling `punched' voids. The KAREN algorithm for reconstructing diffuse scattering provides demonstrably better resolution than can be obtained from previously described punch-and-fill methods. The superior structural resolution obtained using the KAREN method is demonstrated by evaluating the complex ordered diffuse scattering observed from the neutron diffraction of a single plastic crystal of CBr 4 using pair distribution function analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2021

EGR Dilution and Fuel Property Effects on High-Efficiency Spark-Ignition Flames

Modern spark ignition internal combustion engines rely on fast combustion rates and high dilution to achieve high brake thermal efficiencies. To accomplish this, new engine designs have moved towards increased tumble ratios and stroke-to-bore ratios. Increased tumble ratios correlate positively with increases in turbulent kinetic energy and improved fuel and residual gas mixing, all of which favor faster and more efficient combustion. Longer stroke-to-bore ratios allow higher geometric compression ratios and use of late intake valve closing to control peak compression pressures and temperatures. The addition of dilution to improve efficiency is limited by the resulting increase in combustion instabilities manifested by cycle-to-cycle variability. A number of effects - preferential diffusion, turbulence-combustion interactions, stochastic flow patterns, laminar-turbulent flame kernel transitions, and relative length and velocity scales between flame and turbulence - are believed to be responsible for the increase in cycle-to-cycle variations, where their contributions are likely interlinked. Several studies have shown the influence of stochastic flow characteristics on the nature of combustion instabilities, such as velocity patterns on flame kernel formation and cycle-to-cycle variations in residual gas. However, few have focused on the specific effects of fuel properties. The objective of this work is to contrast the effects of dilution on propane stoichiometric combustion against gasoline. Dilution tolerance experiments were conducted in a purpose-built high stroke-to-bore ratio single cylinder engine with both gasoline and LPG. Three-dimensional full cycle computational fluid dynamics (CFD) simulations employing a level-set combustion approach and Reynolds averaged Navier-Stokes (RANS) turbulence modeling was used to qualitatively assess the changes in length and velocity scales for turbulence and the flame. The experimental results showed that LPG can tolerate higher exhaust gas recirculation (EGR) dilution under a variety of conditions. Analysis of CFD simulations showed that propane flames are likely less sensitive to influences from the flow field due to less thickening of the flame and higher effective flame speeds.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Charging of radioactive and environmental airborne particles

The charging of various airborne particles was investigated using single-particle levitation and charge-balance equations. Though radioactive decay and triboelectrification can induce charging, it is typically assumed that the aerosols in a radioactive plume will not carry significant charge at steady state since atmospheric particles can have their charge neutralized through the capture of adjacent counter-ions (i.e., diffusion charging). Here to assess this assumption, we directly measured the surface charge and charge density of various triboelectrically charged aerosols including radioactive uranium oxide (<1 μm), urban dust, Arizona desert dust, hydrophilic and hydrophobic silica nanoparticles, and graphene oxide powders using an electric field-assisted particle levitator in air. Of these particles, uranium oxide aerosols exhibited the highest surface charge density. Charge balance equations were employed to predict the average charge gained from radioactive decay as a function of time and to evaluate the effects of diffusion charging on triboelectrically charged radioactive and non-radioactive particles in the atmosphere. Simulation results show that particles, initially charged through triboelectrification, can be quickly discharged by diffusion charging in the absence of radioactive decay. Nevertheless, simulation results also indicate that particles can be strongly charged when they carry radionuclides. These experimental and simulation results suggest that radioactive decay can induce strong particle charging that may potentially affect atmospheric transport of airborne radionuclides.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Analysis of Silver Release from Furnace-Tested TRISO Particles

High temperature testing of intact TRISO particles previously irradiated in the AGR-5/6/7 experiment was performed in the Furnace for Irradiated TRISO Testing (FITT) to directly confirm diffusive release of silver and europium from intact TRISO particles. Testing was conducted from 1,100–1,600°C for exposure times up to 100 h to directly confirm silver through layer release below safety testing temperatures. The FITT analysis showed highest levels of silver release at 1,300–1,400°C which confirms athermal release behaviors previous observed in step-wise and transient safety tests. Additionally, release was non-uniform with some particles releasing a majority of their inventory while others did not appear to release silver under identical testing conditions, which was consistent with historic observations. An assessment of the effective silver diffusion coefficient in the SiC layer was conducted and indicated maximum values in the 1,300–1,400°C range. The magnitude of the calculated diffusion coefficients also exceeded currently accepted diffusion coefficients. The release behavior of europium from intact particles was also analyzed in FITT for at 1,450–1,550°C for 500 h to 750 h. Direct confirmation of europium release below safety testing temperatures was confirmed absent contributions from matrix release. The analysis indicated europium release follows a general Arrhenius behavior and suggests general uniform release and indicates irradiation conditions influence observed release response. Calculated diffusion kinetics agree well with historic experiments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Local-scale heterogeneity of soil thermal dynamics and controlling factors in a discontinuous permafrost region

In permafrost regions, the strong spatial and temporal variability in soil temperature cannot be explained by the weather forcing only. Understanding the local heterogeneity of soil thermal dynamics and their controls is essential to understand how permafrost systems respond to climate change and to develop process-based models or remote sensing products for predicting soil temperature. In this study, we analyzed soil temperature dynamics and their controls in a discontinuous permafrost region on the Seward Peninsula, Alaska. We acquired one-year temperature time series at multiple depths (at 5 or 10 cm intervals up to 85 cm depth) at 45 discrete locations across a 2.3 km 2 watershed. We observed a larger spatial variability in winter temperatures than that in summer temperatures at all depths, with the former controlling most of the spatial variability in mean annual temperatures. We also observed a strong correlation between mean annual ground temperature at a depth of 85 cm and mean annual or winter season ground surface temperature across the 45 locations. We demonstrate that soils classified as cold, intermediate, or warm using hierarchical clustering of full-year temperature data closely match their co-located vegetation (graminoid tundra, dwarf shrub tundra, and tall shrub tundra, respectively). We show that the spatial heterogeneity in soil temperature is primarily driven by spatial heterogeneity in snow cover, which induces variable winter insulation and soil thermal diffusivity. These effects further extend to the subsequent summer by causing variable latent heat exchanges. Finally, we discuss the challenges of predicting soil temperatures from snow depth and vegetation height alone by considering the complexity observed in the field data and reproduced in a model sensitivity analysis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Grain Boundary Effects on Radiation Damage and Tritium Diffusion in Li–Al–O Ceramics from Molecular Dynamics and Experiments

Understanding irradiation damage and tritium transport in LiAlO2 ceramics is essential for their deployment in tritium-producing burnable absorber rods (TPBARs). Grain boundaries (GBs) play an important role in governing radiation response and tritium transport in LiAlO2 and its secondary phase LiAl5O8. Using molecular dynamics simulations, we investigated defect evolution and tritium diffusion during displacement cascades in single-crystal and bicrystal LiAlO2 and LiAl5O8 at 600 K. The results reveal that GBs suppress the damage in the LiAlO2 bicrystals as compared to single crystals, via interstitial emission. The effect is more pronounced for Li defects, than for Al and O. In LiAl5O8 the damage reduction due to introduction of GB is less pronounced overall. Tritium diffusion coefficients at GBs are higher by factors of 2–10 relative to diffusion inside grain bulk. The effect is particularly pronounced in LiAlO2, where S5, S17, and S25 GBs promote rapid tritium transport, whereas LiAl5O8 exhibits slower diffusion due to reduced free volume at its GBs. These findings evince a trade-off that while GBs mitigate radiation damage by absorbing interstitials, they simultaneously provide fast pathways for tritium migration which is not desirable for TPBAR application. The mechanistic insights gained here establish a foundation for microstructural design strategies to balance radiation tolerance and tritium retention in ceramic breeder materials.

Roy, Ankit↗
DOE OSTI2022

Influence of the complete basis set approximation, tight weighted-core, and diffuse functions on the DLPNO-CCSD(T1) atomization energies of neutral H,C,O-compounds

The impact of complete basis set extrapolation schemes (CBS), diffuse functions, and tight weighted core functions on enthalpies of formation predicted via the DLPNO-CCSD(T1) reduced Feller Peterson-Dixon approach has been examined for neutral H,C,O-compounds. All tested three-point (TZ/QZ/5Z) extrapolation schemes result in mean unsigned deviation (MUD) below 2 kJ mol -1 relative to the experiment. Here, the two-point QZ/5Z and TZ/QZ CBS 1/$l^3_{max}$ extrapolation schemes are inferior to their inverse power counterpart (1/($l_{max}$ + 1/2) 4 ) by 1.3 and 4.3 kJ mol -1 . The CBS extrapolated frozen core atomization energies are insensitive (within 1 kJ mol -1 ) to augmentation of the basis set with tight weighted core functions. The core-valence correlation effects converge already at triple-ζ, although double- ζ/triple- ζ CBS extrapolation performs better and is recommended. The effect of diffuse function augmentation converges slowly, and cannot be reproduced with double- ζ or triple- ζ calculations as these are plagued with basis set superposition and incompleteness errors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗