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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Cost‐Effective Layered Oxide – Olivine Blend Cathodes for High‐Rate Pulse Power Lithium‐Ion Batteries

Abstract Sustainability and supply‐chain concerns require lithium‐ion batteries (LIBs) free from critical minerals, such as nickel and cobalt. While recent advances provide encouraging signs that cobalt can be removed, the question remains how much Ni can be removed from Co‐free layered oxide cathodes before sacrificing critical performance metrics. This study highlights the effect of reducing Ni by benchmarking several Co‐free cathodes with decreasing Ni content. Keeping the energy density the same by increasing the charge voltage, cathodes below 80% Ni content exhibit worsened capacity fade due to increasing oxygen release and electrolyte decomposition. Charge transfer and diffusion kinetics are also hindered with increasing Mn content and exacerbated by resistive surface phases formed at high voltages, rendering lower‐Ni, Co‐free cathodes less competitive than high‐Ni cathodes for high energy and power applications. It is demonstrated blending layered oxide with olivine as an effective alternative to deliver energy density and cycling stability comparable to lower‐Ni cathodes with moderate charging voltages. Blending with 30 wt% olivine LiMn 0.5 Fe 0.5 PO 4 (LMFP) virtually eliminates the diffusion limitation of layered oxides at low state‐of‐charge, with enhanced pulse power characteristics rivaling the high‐Ni counterparts. Cathode blending can further reduce the overall Ni content and cost without the performance limitations of lower‐Ni, Co‐free cathodes.

Lee, Steven↗

DOE OSTI2025

Understanding the impact of the gas diffusion layer structure on catalyst utilization in the PEM water electrolyzer

A multiphysics half-cell model of a polymer electrolyte membrane water electrolyzer (PEMWE) was developed to probe impacts of the detailed 3-dimensional pore structure of the gas diffusion layer (GDL) on performance characteristics. We show that pores in the titanium GDL mesh led to significant underutilization of the catalyst layer (CL), with only 45% of the catalyst effectively utilized. This contradicts the assumption of uniform electron flow across the CL, as shown in graphical abstract (a), as near-zero current was observed near GDL pore regions and the current distribution in CL was influenced by GDL structure, as shown in graphical abstract (b). Instead, oxygen generation was primarily concentrated under the solid titanium regions, diffusing out around the pore walls. High current density peaks were also noted at the GDL-catalyst contact, correlating with degradation hotspots that were directly observed in companion experiments, as shown in graphical abstract (c). Collectively, these findings point to the critical importance of the heterogeneous GDL porous architecture not only for PEMWE efficiency but also for uneven degradation of the CL.

COMSOL↗

DOE OSTI2025

Co‐Doping Approach for Enhanced Electron Extraction to TiO 2 for Stable Inorganic Perovskite Solar Cells

Inorganic perovskite CsPbI 3 solar cells hold great potential for improving the operational stability of perovskite photovoltaics. However, electron extraction is limited by the low conductivity of TiO 2 , representing a bottleneck for achieving stable performance. In this study, a co‐doping strategy for TiO 2 using Nb(V) and Sn(IV), which reduces the material's work function by 80 meV compared to Nb(V) mono‐doped TiO 2 , is introduced. To gain fundamental understanding of the processes at the interfaces between the perovskite and charge‐selective layer, transient surface photovoltage measurements are applied, revealing the beneficial effect of the energetic and structural modification on electron extraction across the CsPbI 3 /TiO 2 interface. Using 2D drift‐diffusion simulations, it is found that co‐doping reduces the interface hole recombination velocity by two orders of magnitude, increasing the concentration of extracted electrons by 20%. When integrated into n–i–p solar cells, co‐doped TiO 2 enhances the projected T S80 lifetimes under continuous AM1.5G illumination by a factor of 25 compared to mono‐doped TiO 2 . This study provides fundamental insights into interfacial charge extraction and its correlation with operational stability of perovskite solar cells, offering potential applications for other charge‐selective contacts.

14 SOLAR ENERGY↗

DOE OSTI2025

Mapping of fracture and ionic conductivity changes in ion implanted solid electrolytes: Insights from molecular dynamics

Ion implantation emerges as a promising technique to address the persistent challenge of lithium (Li) filament growth in solid-state electrolytes as it can induce compressive stresses inhibiting crack growth and deflect dendrites, de facto mitigating early electrolyte failure. In this study, we examine the potential paradox of ion implantation: while aiming to enhance electrolyte performance, the radiation damage associated with implantation might inadvertently compromise both the ionic conductivity and the intrinsic fracture toughness of the material, rendering the material unsuitable for battery applications. Specifically, we employed molecular dynamics simulations to examine the scope of the downsides of ion implantation, specifically: (i) reduced ionic conductivity (due to radiation-induced amorphization) and (ii) mechanical stability (due to radiation-induced embrittlement) in ion-implanted Li 7 La 3 Zr 2 O 12 (LLZO) solid-state electrolytes. We explore how radiation damage impacts LLZO’s crystalline structure, Li-ion diffusion, and fracture properties at various temperatures and radiation damage levels. The study aims to provide insights into the competing effects of ion implantation and suggest potential engineering strategies for developing more robust solid-state electrolytes with improved conductivity and dendrite resistance.

Conductivity↗

DOE OSTI2025

Bayesian learning with Gaussian processes for low-dimensional representations of time-dependent nonlinear systems

This work presents a data-driven method for learning low-dimensional time-dependent physics-based surrogate models whose predictions are endowed with uncertainty estimates. We use the operator inference approach to model reduction that poses the problem of learning low-dimensional model terms as a regression of state space data and corresponding time derivatives by minimizing the residual of reduced system equations. Standard operator inference models perform well with accurate training data that are dense in time, but producing stable and accurate models when the state data are noisy and/or sparse in time remains a challenge. Another challenge is the lack of uncertainty estimation for the predictions from the operator inference models. Our approach addresses these challenges by incorporating Gaussian process surrogates into the operator inference framework to (1) probabilistically describe uncertainties in the state predictions and (2) procure analytical time derivative estimates with quantified uncertainties. The formulation leads to a generalized least-squares regression and, ultimately, reduced-order models that are described probabilistically with a closed-form expression for the posterior distribution of the operators. The resulting probabilistic surrogate model propagates uncertainties from the observed state data to reduced-order predictions. Furthermore, we demonstrate the method is effective for constructing low-dimensional models of two nonlinear partial differential equations representing a compressible flow and a nonlinear diffusion–reaction process, as well as for estimating the parameters of a low-dimensional system of nonlinear ordinary differential equations representing compartmental models in epidemiology.

Data-driven model reduction↗

DOE OSTI2025

Thermal Conductivity of Graphene Coated Copper Under Uniaxial Tensile Mechanical Strain

Graphene continues to demonstrate promise as a highly effective barrier coating, even at only one atom thick. The thermal properties of this coating are also promising to allow diffusion of heat across the surface, as the isolated graphene is an intrinsically good thermal conductor, yet this and its behavior under mechanical deformation has been less extensively studied. This report demonstrates that the in-plane thermal conductivity and interfacial thermal conductance of graphene coatings on copper are affected by mechanical strain. By inducing a strain in the copper substrate, the Raman-active 2D peak exhibits a change in position as well as a change in laser power dependence as the copper substrate is uniaxially elongated to a maximum of 0.5%. Non-linear trends in thermal conductivity are observed with tensile strain in samples with differing strain transfer rates from the substrate, indicating the close correlation between intrinsic thermal conduction and interfacial properties in atomically thin coatings transferred onto metals.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Universality in diffusion-controlled nucleation and growth

Nucleation and growth are studied in a system that undergoes diffusion-controlled condensation under gradual changes in parameters, such as cooling. Here, it is demonstrated that when the Gibbs–Thompson effect becomes negligible, the system falls into a universal regime. i.e., the final droplet size distribution remains invariant under certain rescaling of system parameters. An approximate yet very accurate analytic form is obtained for the final droplet size distribution in this regime.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Galactic echoes

Gaia has revealed a variety of substructures in the phase space of stars in the Solar neighbourhood, including the vertical ‘Snail’ $(z, v_z)$ in space. Such substructures are often interpreted as the incompletely phase-mixed response of the disc stars to a single perturbation, such as an impulsive encounter with a satellite galaxy. In this paper, we consider the possibility that such structures contain manifestations of phase-space echoes. First established in plasma physics in the 1960s, echoes arise when a collisionless system is perturbed twice: the macroscopic responses to both perturbations mix to small scales in phase space, whereupon they couple non-linearly, producing a third macroscopic ‘echo’ response without the need for a third perturbation. We derive the galactic analogue of the plasma echo theory using angle-action variables and apply it to a one-dimensional model of vertical motion in the Milky Way. We verify the predicted echo behaviour using idealized test particle simulations, both with and without the inclusion of diffusion through orbital scattering off molecular clouds. While we conclude that the Gaia Snail itself is unlikely a (pure) echo effect, the basic physics we uncover is sufficiently generic that we expect phase-space echoes to be common in disc galaxies.

Galaxy: evolution↗

DOE OSTI2024

Slant, Fan, and Narrow: The Response of Stellar Streams to a Tilting Galactic Disk

Stellar streams are sensitive tracers of the gravitational potential, which is typically assumed to be static in the inner Galaxy. However, massive mergers like Gaia-Sausage-Enceladus can impart torques on the stellar disk of the Milky Way that result in the disk tilting at rates of up to 10°–20° Gyr –1 . Here, we demonstrate the effects of disk tilting on the morphology and kinematics of stellar streams. Through a series of numerical experiments, we find that streams with nearby apocenters (r apo ≲ 20 kpc) are sensitive to disk tilting, with the primary effect being changes to the stream's on-sky track and width. Interestingly, disk tilting can produce both more diffuse streams and more narrow streams, depending on the orbital inclination of the progenitor and the direction in which the disk is tilting. Our model of Pal 5's tidal tails for a tilting rate of 15° Gyr –1 is in excellent agreement with the observed stream's track and width, and reproduces the extreme narrowing of the trailing tail. We also find that failure to account for a tilting disk can bias constraints on shape parameters of the Milky Way's local dark matter distribution at the level of 5%–10%, with the direction of the bias changing for different streams. Disk tilting could therefore explain discrepancies in the Milky Way's dark matter halo shape inferred using different streams.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Impact of a Bluff Body on a Pressurized Oxy-fuel Flame

Designing an effective burner is vital for the development of pressurized oxy-fuel combustion technologies. In the present work, turbulent jet diffusion burners are adopted for a pressurized oxy-combustor, with a bluff-body employed to create a recirculation zone, thereby stabilizing the flame in such a combustor. The objective of this numerical study is to perform a systematic analysis of the characteristics of such a pressurized non-premixed flame. Specifically, a 15-bar pressurized oxy-fuel combustor of power 100 kWth is modeled by means of the Ansys FLUENT commercial platform, using the Reynolds-averaged Navier-Stokes (RANS) approach. The present work focuses on identifying the aerothermodynamic features of the pressurized oxy-fuel burner with a disk-shape bluff-body. It is shown that the fuel-to-oxidizer stream momentum ratio has a great impact on the temperature profile of the down-fired, co-axial, pressurized oxy-coal, diffusion flame. A parametric study of the blockage ratio of the burner identifies an optimal range for this pressurized burner.

Li, Lei↗

DOE OSTI2025

The Crucial Role of Vacancy Concentration in Enabling Superatomic Diffusion in Lithium Intermetallics

Anode-free solid-state Li batteries promise significant increases in energy densities compared to current commercial batteries that rely on liquid electrolytes. Major challenges persist in controlling morphological evolution during the plating and stripping of lithium metal at the anode current collector. Elemental additives that alloy with lithium have been found to modify the plating and stripping behavior of lithium. Many alloying elements form intermetallics with lithium and the mobility of Li through these intermetallics is believed to have an important effect on morphological evolution. This study shows that Li transport coefficients through intermetallics span a wide range in values, with the B32 LiAl intermetallic predicted to have a Li tracer diffusion coefficient as high as 10 –6 cm 2 /s at room temperature, which is 8 orders of magnitude larger than that of isostructural B32 LiZn. This work demonstrates the crucial role of vacancy concentration in controlling the mobility of Li atoms through intermetallics. While the migration barriers for Li-vacancy exchanges in both LiAl and LiZn are remarkably low, the superatomic conductivity in LiAl is shown to arise from the unique electronic structure of the B32 LiAl compound, which favors high concentrations of vacancies.

25 ENERGY STORAGE↗

DOE OSTI2024

Development of a Griffin model of the advanced test reactor

In the pursuit of a higher fidelity deterministic simulation capability of the Advanced Test Reactor, it is important to have a fast yet accurate deterministic neutronics model. Here, to achieve this, we employed an advanced two-step method. The first step involves generating homogenized cross sections using OpenMC, a cutting-edge Monte Carlo neutron transport code. OpenMC offers excellent modular capabilities, allowing for easy component integration and flexibility in incorporating new designs into the model. The second step involves deterministic transport calculations, which are performed using Griffin, a reactor physics application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE). To ensure the accurate spatial resolution and assignment of material cross sections, a Cubit-generated mesh for the Advanced Test Reactor is utilized as an intermediate step between the OpenMC and Griffin models; Griffin utilizes the mesh for its finite element solution, while OpenMC material identifications are written to the mesh file to be used in Griffin material assignments. Additionally, a Python-based script converts the cross sections generated by OpenMC into the ISOXML format required by Griffin. Initial comparisons using the Griffin diffusion solver indicated good agreement between the neutron multiplication factors obtained from the standalone OpenMC model and the Griffin model, with differences of less than 10 pcm in the 2D geometry configuration; it was later determined that this agreement was likely due to compensating effect and was more likely on the order of –700 pcm relative to the OpenMC solution. However, in three-dimensional calculations, an unacceptably large error (almost 8,000 pcm) was found in the Griffin solution with the diffusion solver. Subsequent calculations using Griffin’s discrete ordinates solver demonstrated substantially improved agreement, within 116 pcm of the OpenMC solution used to generate the cross sections for Griffin. Building on this capability, future work will seek to perform more detailed validation calculations. The ultimate goal is to evaluate both transient and multiphysics simulations of the reactor.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2024

Dynamic density functional theory of polymers with salt in electric fields

Here we present a dynamic density functional theory for modeling the effects of applied electric fields on the local structure of polymers with added salt (polymer electrolytes). Time-dependent equations for the local electrostatic potential and volume fractions of polymer, cation, and anion of added salt are developed using the principles of linear irreversible thermodynamics. For such a development, a field theoretic description of the free energy of polymer melts doped with salts is used, which captures the effects of local variations in the dielectric function. Connections of the dynamic density functional theory with experiments are established by relating the three phenomenological Onsager’s transport coefficients of the theory to the mutual diffusion of electrolyte, ionic conductivity, and transference number of one of the ions. The theory is connected with a statistical mechanical model developed by Bearman and Kirkwood [J. Chem. Phys. 28, 136 (1958)] after relating the three transport coefficients to friction coefficients. The steady-state limit of the dynamic density functional theory is used to understand the effects of dielectric inhomogeneity on the phase separation in polymer electrolytes. The theory developed here provides not only a way to connect with experiments but also to develop multi-scale models for studying connections between local structure and ion transport in polymer electrolytes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Intrinsic kinetics of water-inhibited ultra-lean methane oxidation over PtPd-Mg/ θ -Al 2 O 3 catalyst

Here, this study develops intrinsic methane oxidation kinetics for ultra-lean methane conditions in the presence of water over a highly active and stable PtPd–Mg/θ-Al 2 O 3 catalyst. Comprehensive laboratory experiments were conducted over a wide range of methane concentrations (150–1200 ppm CH 4 ), water contents (1–5% H 2 O), and industrially relevant space velocities (80,000 ≤ GHSV ≤ 160,000 h -1 ). These systematic experiments informed a two-dimensional, axisymmetric, multiscale reactor model that was used to develop and validate methane oxidation kinetics under practically relevant conditions, including non-isothermal operation and high conversion regimes. Combined experimental and modeling results revealed significant intraparticle diffusion resistance and transport-induced reaction exotherm at elevated temperatures, which limited catalyst utilization despite high intrinsic activity. These transport effects were explicitly incorporated into the reactor model, enabling accurate estimation of intrinsic kinetic parameters without reliance on conventional effectiveness-factor corrections. The resulting kinetic model successfully captured both kinetically controlled and mass-transfer-limited regimes and reliably predicted CH 4 conversion across broad ranges of temperature, methane concentration, and water content

Heat and mass transfer limitations↗

DOE OSTI2024

Influence of misfit dislocations on ionic conductivity at oxide interfaces

Mismatched complex oxide thin films and heterostructures have gained significant traction for use as electrolytes in intermediate temperature solid oxide fuel cells, wherein interfaces exhibit variation in ionic conductivity as compared to the bulk. Although misfit dislocations present at interfaces in these structures impact ionic conductivity, the fundamental mechanisms responsible for this effect are not well understood. To this end, a kinetic lattice Monte Carlo (KLMC) model was developed to trace oxygen vacancy diffusion at misfit dislocations in SrTiO 3 /BaZrO 3 heterostructures and elucidate the atomistic mechanisms governing ionic diffusion at oxide interfaces. The KLMC model utilized oxygen vacancy migration energy barriers computed using molecular statics. While some interfaces promote oxygen vacancy diffusion, others impede their transport. Fundamental factors such as interface layer chemistry, misfit dislocation structure, and starting and ending sites of migrating ions play a crucial role in oxygen diffusivity. Molecular dynamics (MD) simulations were further performed to support qualitative trends for oxygen vacancy diffusion. Overall, the agreement between KLMC and MD is quite good, though MD tends to predict slightly higher conductivities, perhaps a reflection of nuanced structural relaxations that are not captured by KLMC. The current framework comprising KLMC modeling integrated with molecular statics offers a powerful tool to perform mechanistic studies focused on ionic transport in thin film oxide electrolytes and facilitate their rational design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Puzzling 𝐵⁡(𝐸⁢2;0 + → 2 + ) strength in the proton dripline nucleus 36 Ca

Recent measurements of the 𝐸⁢2 transition rate from the ground state to the first 2 + excited state of the proton dripline nucleus 36 Ca show an unusual pattern when compared to its isotopic neighbor 38 Ca: despite having a higher 𝐸 𝑥 ⁡(2$^{+}_{1}$) excitation energy, the 𝐵⁡(𝐸⁢2;0$^{+}_{1}$ → 2$^{+}_{1}$) rate in 36 Ca is larger. The question that naturally arises is to what extent this observation can be attributed to the unbound character of the 2$^{+}_{1}$ state. To understand the influence of the continuum space on the low-energy properties of 36 Ca, we carried out Gamow shell model calculations that can account for the continuum coupling effects associated with the occupation of unbound 𝑓⁡𝑝 shells. We found that in the threshold 2 + state, 36 Ca is spatially diffused, which impacts the observed 𝐵⁡(𝐸⁢2) trend.

20 ≤ A ≤ 38↗

DOE OSTI2024

Tusqh

SAND2025-00675O Tusqh is a software tool that generates cubical meshes in 2D and 3D and computes the homology of these meshes using persistent homology. It includes a grid cell in the output if its volume-fraction is above a selectable threshold, estimated by sampling points within the cell. Tusqh incorporates anti-aliasing algorithms to mitigate grid orientation and scale effects. It is designed for creating finite element meshes for simulations and can be used in various applications such as heat diffusion, mechanical simulations, and computer graphics rendering. The software outputs meshes in an open format compatible with downstream software. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗

DOE OSTI2026

Formation of Non-Doped Cubic Lithium Lanthanum Zirconium Oxide Nanofibers: Insights from In Situ Synchrotron X-Ray Scattering

This study investigates the formation mechanism of non-doped cubic lithium lanthanum zirconium oxide (c-LLZO) nanofibers using in situ synchrotron X-ray scattering techniques. Electrospun polymer precursor nanofibers were annealed at temperatures up to 800 °C, enabling real-time tracking of phase transitions via simultaneous small-angle X-ray scattering (SAXS), wide-angle X-ray scattering (WAXS), and evolved CO 2 gas analysis. The results reveal a three-step transformation pathway: polymer decomposition, formation of La 2 Zr 2 O 7 (LZO), and direct conversion of LZO to c-LLZO without intermediate tetragonal phases detected within the sensitivity of our in situ WAXS measurement. Cryo-electron energy loss spectroscopy (EELS) further elucidates the role of lithium diffusion, showing Li enrichment at fiber surfaces and Li deficiency in the interior, which stabilizes the cubic phase. This Li segregation effect in nanostructured LLZO materials extends beyond the previously reported size effect. This work advances the understanding of c-LLZO formation mechanisms and provides practical insights for optimizing synthesis routes to achieve phase-pure c-LLZO for solid-state battery applications.

LLZO phase stability↗