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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Simplified calculation of the area specific impedance for solid-state battery design

Simplified algebraic area specific impedance (ASI) correlations have been developed for solid-state composite battery electrodes made of a single ion conducting electrolyte, conductive additive, and intercalation active material. Two ASI expressions were developed, one for short times ( i.e. , pulsed power operation) and another for the pseudo steady state operation (i.e., sustained discharge for energy estimation). A full electrochemical model based on porous electrode theory was developed to examine the accuracy of the simplified ASI expressions. The simplified expressions agree favorably with full model results over a wide range of parameters (i.e. , electrode thicknesses, electrolyte conductivities, solid-state diffusion coefficients, specific surface areas, etc.) and conditions (i.e. , C-rates, states of charge, and pulse times). Under most conditions, the error between the full model and the correlations is well below 7 %. Higher errors were observed for the pseudo steady state expression at high/ low states of charge where the assumption of uniform reaction distributions loses validity. Here, the short time ASI has higher error at low states of charge due to the nonlinearity of the open circuit voltage equation, which is assumed linear in the formulation of the simplified algebraic expression.

25 ENERGY STORAGE↗
DOE OSTI2023

A conformal heat-drying direct ink writing 3D printing for high-performance lithium-ion batteries

High areal capacity electrodes hold great potential for high-energy density lithium-ion batteries (LIBs), but their poor electrochemical kinetics limit their power density. Here, in this study, high areal capacity 3D-structured LiNi 0.8 Mn 0.1 Co 0.1 O 2 cathodes (4.3 mAh cm -2 ) are prepared via 3D printing with a manner of direct ink writing. The electrodes had an enlarged electrode–electrolyte contact area, shortened diffusion pathway, and reduced intercalation-induced stress, thereby delivering enhanced rate capability and cyclability in LIBs, which is 143.6 mAh g -1 at 3C and a 60.2 % capacity retention over 800 cycles at 1C. Moreover, at electrode level, the 3D-NMC exhibits an energy and power densities of 313.1 Wh kg -1 and 657.9 W kg -1 , respectively. Furthermore, the theoretical calculation suggests that reducing the gap width will be highly beneficial to the energy and power densities. This work establishes a milestone in understanding the cycling effect on the electrode local structure, including the void area and the LiNi 0.8 Mn 0.1 Co 0.1 O 2 region, which confirms the effectiveness of 3D printing for electrode preparation.

25 ENERGY STORAGE↗
DOE OSTI2020

Predicting cell-to-cell failure propagation and limits of propagation in lithium-ion cell stacks

Thermal runaway of lithium-ion batteries is a risk that is magnified when stacks of lithium-ion cells are used for large scale energy storage. When limits of propagation can be identified so that systems can be designed to prevent large scale cascading failure even if a failure does occur, these systems will be safer. This work addresses the prediction of cell-to-cell failure propagation and the propagation limits in lithium-ion cell stacks to better understand and identify safe designs. A thermal-runaway model is presented based on recent developments in thermochemical source terms. It is noted that propagating failure is characterized by temperatures above which calorimetry data is available. Results show high temperature propagating failure predictions are too rapid unless an intra-particle diffusion limit is included, introducing a Damköhler number limiter into the rate expression. This new model form is evaluated against cell-to-cell failure propagation where the end cell of a stack is forced into thermal runaway through a nail-induced short circuit. Limits of propagation for this configuration are identified. Results show cell-to-cell propagation predictions are consistent with measurements over a range of cell states of charge and with the introduction of metal plates between cells to add system heat capacity representative of structural members. This consistency extends from scenarios where propagation occurs through scenarios where propagation is prevented.

25 ENERGY STORAGE↗
DOE OSTI2024

Freeze Tape Casting Electrode with Bilayered Architecture for High-Performance Lithium-Ion Batteries

A facile freeze tape casting (FTC) strategy is utilized to prepare bilayered 4 mAh cm –2 high-loading LiNi 0.6 Co 0.2 Mn 0.2 O 2 cathodes. The bottom layer is a conventional nonaqueous electrode, which has a dense structure for high-energy purposes. Further, the top layer is prepared by the proposed FTC, exhibiting a porous feature for high-power requirement. With the assistance of FTC, the bilayered electrodes successfully deliver enhanced rate and cyclic performance due to the improved lithium-ion diffusion kinetics and pathways. Therefore, the proposed FTC strategy and its delivered electrodes are promising for energy- and power-density lithium-ion batteries, potentially enlightening the research and development of lithium-ion battery manufacturing.

25 ENERGY STORAGE↗
DOE OSTI2024

Ce- and Ni-Codoped Double PrBaMn 2 O 5 Perovskite as a Ceramic SOFC Anode

This study explores the efficacy of cerium introduction, both on its own and in combination with nickel, into PrBaMn 2 O 5+δ (PBM) structures to enhance solid oxide fuel cell (SOFC) anodes. We synthesized Pr 1–x Ce x BaMn 2 O 5+δ compositions for x values of 0.05 (PrCe5) and 0.1 (PrCe10), as well as a nickel-doped variant (PrCe5Ni), assessing their performance under H 2 – 3% H 2 O reducing conditions pertinent to SOFC anode operations. Our findings reveal that the PrCe5 composition exhibits a thermal expansion coefficient (TEC) that not only improves upon that of the Ce-free counterpart but also aligns closely with the TEC standards of prevalent SOFC electrolytes. A notable advancement was achieved with the application of a 3.5 μm gadoliniadoped ceria (GDC) buffer layer through physical vapor deposition, effectively mitigating chemical interactions between PBM-based anodes and yttria-stabilized zirconia (YSZ) electrolytes, a concern highlighted by in situ neutron diffraction analyses. Electrochemical impedance spectroscopy, conducted over 220 h in a H 2 –3% H 2 O atmosphere at 750 °C, demonstrated that the optimal 3.5 μm thickness of the GDC buffer layer significantly minimizes the area-specific resistance (ASR) degradation rate to 0.002 Ω cm 2 /h, markedly outperforming both thinner (1 μm) and thicker (8 μm) GDC layers, which showed higher degradation rates of 0.15–0.2 Ω cm 2 /h due to the diffusion of Ba ions or delamination. Moreover, cerium doping fosters superior microstructural stability and obviates barium diffusion, thereby suggesting an enhanced durability of the doped anodes over their lifespan. Integrating nickel into the PrCe5 structure halved the ASR to 0.5 Ω cm 2 at 750 °C, situating it well within the ideal performance range for SOFC anodes. Furthermore, the enhancement brought about by simultaneously doping PBM with cerium and nickel, which fundamentally relies on the critical contributions of defect chemistry and crystal structure, highlights the significance of these fields in creating sophisticated materials for energy related applications.

25 ENERGY STORAGE↗
DOE OSTI2020

Soil Respiration Response to Rainfall Modulated by Plant Phenology in a Montane Meadow, East River, Colorado, USA

Soil respiration is a primary component of the terrestrial carbon cycle. However, predicting the response of soil respiration to climate change remains a challenge due to the complex interactions between environmental drivers, especially plant phenology, temperature, and soil moisture. In this study, we use a 1-D diffusion-reaction model to calculate depth-resolved CO 2 production rates from soil CO 2 concentrations and surface efflux observations in a subalpine meadow in the East River watershed, CO. Modeled rates are compared to in situ soil temperature and moisture conditions and MODIS satellite enhanced vegetation index (EVI) representing plant phenology across three hydrologically distinct growing seasons from 2016-2018. While soil respiration correlated with temperature on diel timescales (p<0.05), seasonal variability was dominated by soil moisture and plant phenology (p<0.05). We observed significant respiration increases in response to precipitation events; however, magnitude and duration were significantly higher in 2017 than 2016 despite similar wetting characteristics. Based on MODIS EVI, we suggest the respiration response to rainfall is controlled by plant phenology, which in turn reflects the capacity of plants to respond to precipitation via increased photosynthesis and autotrophic respiration, behavior that is not captured in typical soil respiration pulse models. Projected changes in montane climate such as earlier snowmelt and prolonged fore-summer drought may decrease soil respiration fluxes by decreasing the overlap between peak productivity and the summer monsoon. Lastly, we observed significant late-season CO 2 fluxes from the deep sub-soil (>165 cm) that support growing evidence for the importance of sub-soil processes in driving integrated respiration fluxes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Charge carrier extraction and recombination effects in GaInAs/GaAsP multi-quantum well solar cells

The carrier extraction and transport mechanisms as well as the relative contributions of radiative and non-radiative recombination processes are investigated in high-quality strain-balanced GaInAs/GaAsP multi-quantum well solar cells recently implemented in record efficiency multijunction solar cells. A comprehensive suite of complementary characterization techniques including temperature- and suns-dependent photoluminescence and photovoltaic measurements are employed to analyze thermal escape and tunneling rates, which demonstrate the need to move beyond simple drift-diffusion models of p–n junctions. This study examines the processes that best characterize the operation of these devices across varying temperatures using a simple two-diode model, incorporating multiple transport protocols, and provides insights into the performance-limiting processes and pathways for their optimization.

14 SOLAR ENERGY↗
DOE OSTI2026

Investigating shock-induced chemical reactions in Ni+Al multilayers: A continuum-based mesoscale approach with Arrhenius kinetics and artificial thermal conduction

A new continuum-based mesoscale modeling approach for shock-induced chemical reactions (SICRs) in Ni+Al multilayers is demonstrated in Sandia’s shock physics hydrocode, CTH. The approach utilizes Arrhenius-type kinetics and artificial thermal conduction. Our work builds upon previous efforts to parameterize equations of state for Ni x Al y [J. Appl. Phys. 137, 075102 (2025)], as well as simulations of inert shocks in realistic 2D microstructures [J. Appl. Phys. 137, 225301 (2025)]. To calibrate the reaction kinetics, pairs of the reaction coordinate, R′, vs time are extracted from the molecular dynamics (MD) literature. Here, the MD-informed kinetics are used to simulate the dynamic evolution of pressure and temperature in 2D mesoscale simulations. Overall, the MD-informed kinetics obtained for planar interfaces are too slow, as initial reaction is not observed on a nanosecond time scale. Even with quasi-periodic shock focusing leading to the formation of so-called “hot-spots,” the hot spots are unable to grow and coalesce using the fitted Arrhenius rate constants. However, by increasing the rate constants by two orders of magnitude, SICRs are observed at a shock pressure near 30 GPa, which is supported by experiments. Consequently, these mesoscale simulations suggest that unresolved shear-based mechanical mixing might possibly account for the discrepancies in kinetic rates, with shock-generated intense perturbations, interfacial vortical flows, and elevated temperatures serving as favorable reaction conditions. Future work will calibrate a shear-dependent reaction rate from the MD simulations with realistic interfaces that are wavy, diffuse, and disordered.

Kittell, David E. [Sandia National Laboratories (S↗
DOE OSTI2022

Pseudocapacitive Charge Storage in Electrochromic Transition-Metal Oxide Thin Films

Electrochromic pseudocapacitive transition-metal oxide materials, such as tungsten oxide, which combine fast response, high energy density, and optical effects, can play a significant role as energy storage materials. Here we investigate the electrochemical kinetics of thin films of tungsten oxide, which turn transparent to sky-blue color in the lithiated state due to the reduction of W 6+ to W 5+ . In this study, we investigated the charge density, charge transfer, ion diffusion, and interfacial behavior upon Li + insertion/de-insertion in WO 3 . The electrochromic thin film's pseudocapacitive and electrical double layer mechanism was differentiated based on the power-law. Faradaic diffusion-controlled process dominates over the surface capacitive behavior at scan rates below 40 mV s –1 . These films exhibit an areal charge density of around 100 mC cm –2 and a capacitance of 80 mF cm –2 , superior to most comparable electrochromic materials and supercapacitors. This work combines electrochromics and energy storage properties and provides a fundamental understanding of pseudocapacitive and electrochromic mechanisms in WO 3 .

25 ENERGY STORAGE↗
DOE OSTI2024

Electrochemical Reactivity of Hydrogen Storage Materials: Exploring Borohydride and Hydrazinium Salt Mixtures

Electrochemical characterization of hydrogen storage materials was conducted in a non-aqueous environment to investigate the direct electrochemical release and consumption of hydrogen and the potential for regeneration. We first address the challenge of minimal solubility of the synthetic precursors, sodium borohydride (NaBH 4 ) and hydrazinium bromide (N 2 H 5 Br), in both organic and inorganic solvents. We next determine and calibrate a reference electrode formulation compatible with our non-aqueous media and analytes that demonstrates a stable reference potential. We employ cyclic voltammetry (CV) to characterize the precursors and mixtures thereof. Each CV peak is assigned to a corresponding electrochemical reaction. Using the rate-dependent CV method and Randles–Ševčík equation, we calculate the diffusion coefficient of each chemical (NaBH 4 and N 2 H 5 Br). Analysis of the CVs, coupled with 11B NMR analysis, reveals a room temperature chemical transformation of NaBH 4 and N 2 H 5 Br mixtures into hydrazine borane (N 2 H 4 BH 3 ). These results are particularly significant, considering the limited information available on the electrochemical characterization of metal borohydride and hydrazinium salt in non-aqueous media. This work establishes a foundation for adapting a non-aqueous electrochemical system to further study the borohydride family of chemistries and to design and develop electrochemical devices for direct electrical and chemical energy interconversion with hydrogen storage materials.

08 HYDROGEN↗
DOE OSTI2021

Adaptive Mesh Refinement for Parallel in Time Methods

The project applied the multigrid-reduction-in-time (MGRIT) algorithm to an existing sub-cycled adaptive mesh refinement (AMR) code to investigate the performance of flows dominated by inertial physics. Previous work demonstrated good performance from MGRIT+AMR applied to flows dominated by diffusive physics. Consistent with previous experience, inertial physics negatively affected convergence rates and performance. Efforts to circumvent this issue by appealing to the physics of turbulence were investigated. It has been demonstrated that scales can be effectively transferred between multigrid levels for a turbulent flow resulting in a) partial convergence observed and b) nearly identical results to sequential time-stepping. Performance improvements have not yet been demonstrated - attempts at coarsening the grid on coarser MG levels compromises the solution quality and leads to divergence. This report summarize the accomplishments for the time-frame from 10/5/2020 to 12/31/2020 with an informal no-cost-extension to 05/20/2021

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Laboratory Measurements of Cloud Scavenging of Interstitial Aerosol in a Turbulent Environment (Final Technical Report)

The link between aerosol particles and clouds is well established. Every cloud droplet forms on a preexisting particle; clouds are the primary removal mechanism for aerosol particles with diameters in the range of 0.1 to 1 micrometer. Traditionally, those connections have been considered in the absence of the chaotic air motions (turbulence) that are ubiquitous in Earth’s atmosphere. Recent work at Michigan Tech indicates that turbulent fluctuations in temperature and water vapor, which couple aerosols and cloud droplets, play an important role in whether particles grow to the size of cloud droplets and scavenging by diffusion. In this project, we measured, in a turbulent environment, the rate at which aerosol particles are incorporated into cloud droplets through both these processes, which Michigan Tech's turbulent cloud chamber enables because we can sustain a mixing cloud almost indefinitely.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Simple, Tough, and Effective Porous Coating for Anti-Corrosion Application

This presentation discusses how a change in polymeric binder, the material used for assembling hydrophobic nanoparticles into a porous structure, led to better mechanical toughness and a simpler fabrication process based on 100% commercially available components. Most importantly, the optimized porous coating can reduce the corrosion rate by two orders of magnitude (0.003 mm/year) due to higher diffusion barrier based on a more resistant tortuous diffusion path. The combination of simplicity, toughness, effectiveness, and scalability greatly enhanced the prospect of commercializing this technology as an internal anti-corrosion coating for natural gas pipelines.

internal corrosion↗
DOE OSTI2025

Transition to Petschek Reconnection in Subrelativistic Pair Plasmas: Implications for Particle Acceleration

While relativistic magnetic reconnection in pair plasmas has emerged in recent years as a candidate for the origin of radiation from extreme astrophysical environments, the corresponding subrelativistic pair-plasma regime has remained less explored, leaving open the question of how relativistic physics affects reconnection. In this paper, we investigate the differences between these regimes by contrasting two-dimensional particle-in-cell simulations of reconnection in pair plasmas with relativistic magnetization (σ ≫ 1) and subrelativistic magnetization (σ < 1). By utilizing unprecedentedly large domain sizes and outflow boundary conditions, we demonstrate that lowering the magnetization results in a change in the reconnection geometry from a plasmoid chain to a Petschek geometry, where laminar exhausts bounded by slow-mode shocks emanate from a single diffusion region. We attribute this change to the reduced plasmoid production rate in the low-σ case: When the secondary tearing rate is sufficiently low, plasmoids are too few in number to prevent the system from relaxing into a stable Petschek configuration. This geometric change also affects particle energization: We show that while high-σ plasmoid chains generate power-law energy spectra, low-σ Petschek exhausts merely heat incoming plasma and yield negligible nonthermal acceleration. These results have implications for predicting the global current sheet geometry and the resulting energy spectra in a variety of systems.

Plasma astrophysics↗
DOE OSTI2023

CompLaB v1.0: a scalable pore-scale model for flow, biogeochemistry, microbial metabolism, and biofilm dynamics

Abstract. Microbial activity and chemical reactions in porous media depend on the local conditions at the pore scale and can involve complex feedback with fluid flow and mass transport. We present a modeling framework that quantitatively accounts for the interactions between the bio(geo)chemical and physical processes and that can integrate genome-scale microbial metabolic information into a dynamically changing, spatially explicit representation of environmental conditions. The model couples a lattice Boltzmann implementation of Navier–Stokes (flow) and advection–diffusion-reaction (mass conservation) equations. Reaction formulations can include both kinetic rate expressions and flux balance analysis, thereby integrating reactive transport modeling and systems biology. We also show that the use of surrogate models such as neural network representations of in silico cell models can speed up computations significantly, facilitating applications to complex environmental systems. Parallelization enables simulations that resolve heterogeneity at multiple scales, and a cellular automaton module provides additional capabilities to simulate biofilm dynamics. The code thus constitutes a platform suitable for a range of environmental, engineering and – potentially – medical applications, in particular ones that involve the simulation of microbial dynamics.

58 GEOSCIENCES↗
NASA NTRS1973

Magnetic field reconnection and particle acceleration

After a brief review of steady state reconnection models, the processes are discussed which place an upper limit on the reconnection rate. In many applications the requirement that the width of the diffusion region exceed the ion gyroradius leads to a severe restriction on this limit. Unless the conductivity is much less than the value given by Spitzer's formula, solar flare times may result which are too long by a factor of 1000 or more. Turbulent and gyro effects which may lead to a reduced conductivity are shown to be operative, but only for sufficiently small ratios of plasma to magnetic pressure in the inflow. Particle energization processes associated with magnetic field reconnection and particle reflection by fast shocks are considered. These processes are effective as a first stage in the acceleration of charged particles, but additional multistep processes are needed to reach proton energies above a few MeV.

Sonnerup, B. U. O.↗
NASA NTRS1978

Energy spectrum of medium energy gamma-rays from the galactic center region

A balloon-borne magnetic core digitized spark chamber with two assemblies of spark-chambers above and below the scintillation counters was used to measure the medium energy gamma ray flux from the galactic center region. Gamma ray calculations are based on the multiple scattering of the pair electrons in 15 aluminum plates interleaved in the spark chamber modules. Counting rates determined during ascent and at ceiling indicate the presence of diffuse component in this energy range. Preliminary results give an integral flux between 15 and 70 MeV compared to the differential points in other results.

Palmeira, R. A. R.↗