Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “percolation”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Non-equilibrium simulations of hydraulic permeation: Role of mechanical boundary conditions in dense membranes

Understanding the mechanisms of water transport in reverse osmosis membranes is critical for improving membrane performance and guiding material design. While classical models describe transport as either solution-diffusion (SD)—involving concentration-driven diffusion through a homogeneous medium—or pore-flow (PF)—involving pressure-driven convection through percolated water channels—their applicability to crosslinked polyamide membranes remains debated. Here, using non-equilibrium molecular dynamics simulations, we investigate the impact of mechanical support conditions on pressure-driven water transport in polyamide membranes across varying crosslink densities and pressure differentials (1000–5000 bar). Two support conditions are considered: graphene-restrained, representing experimentally relevant supported membranes, and freeze-restrained, mimicking a self-supported structure. In graphene-restrained systems, water concentration gradients and constant pressure profiles emerge, consistent with SD theory and incompatible with PF assumptions due to the absence of percolated pores and sub-nanometer voids. In contrast, freeze-restrained systems display uniform water concentration and linearly decreasing pressure at 1000 bar, and exhibit compressibility-induced water gradients and partial percolation at 5000 bar, resembling PF-like behavior. However, the underlying assumptions of PF theory—continuous solvent pathways and pressure transmission through water-filled pores—are not met under most conditions. Our results demonstrate that accurate modeling of reverse osmosis membranes must incorporate realistic mechanical boundary conditions to distinguish between transport mechanisms. For dense polyamide membranes supported by porous substrates, graphene-restrained simulations best reflect experimental setups and support the SD model as the dominant mechanism of water permeation.

molecular dynamics↗
DOE OSTI2026

Correlated Ion Transport Governed by Dynamic Local Structure in High Concentration and Localized High Concentration Electrolytes

Understanding the dynamics of cluster formation and network percolation provides the mechanistic link between microscopic solvation structure and transport in concentrated electrolytes, including localized high-concentration electrolytes (LHCEs). Although recent studies have shown LHCEs to form micelle-like aggregates at specific compositions, a quantitative understanding of how solvation structures and transport properties depend on salt–solvent–diluent ratios remains limited. Here, we integrate molecular dynamics with Onsager transport analyses to chart the evolution of solvation microstructure and associated ionic transport in LiFSI/DMC/TTE, achieving good agreement with experimental conductivities across composition. We show that micelle formation arises from a dynamic instability of the cation–anion network, and that network percolation is the primary determinant of conductivity. Two compositional thresholds emerge: A critical network concentration (CNC) and a critical micelle concentration (CMC) that delineate transitions from extended percolating networks to micelle-like clusters and then to fragments. These structural transitions rationalize the nonintuitive conductivity decrease at intermediate dilution despite monotonically decreasing viscosity and provide composition-level design rules for LHCEs.

Mohanakrishnan, Rohith Srinivaas [University of Ca↗
NASA NTRS1993

Unusual radar echoes from the Greenland ice sheet

In June 1991, the NASA/Jet Propulsion Laboratory airborne synthetic-aperture radar (AIRSAR) instrument collected the first calibrated data set of multifrequency, polarimetric, radar observations of the Greenland ice sheet. At the time of the AIRSAR overflight, ground teams recorded the snow and firn (old snow) stratigraphy, grain size, density, and temperature at ice camps in three of the four snow zones identified by glaciologists to characterize four different degrees of summer melting of the Greenland ice sheet. The four snow zones are: (1) the dry-snow zone, at high elevation, where melting rarely occurs; (2) the percolation zone, where summer melting generates water that percolates down through the cold, porous, dry snow and then refreezes in place to form massive layers and pipes of solid ice; (3) the soaked-snow zone where melting saturates the snow with liquid water and forms standing lakes; and (4) the ablation zone, at the lowest elevations, where melting is vigorous enough to remove the seasonal snow cover and ablate the glacier ice. There is interest in mapping the spatial extent and temporal variability of these different snow zones repeatedly by using remote sensing techniques. The objectives of the 1991 experiment were to study changes in radar scattering properties across the different melting zones of the Greenland ice sheet, and relate the radar properties of the ice sheet to the snow and firn physical properties via relevant scattering mechanisms. Here, we present an analysis of the unusual radar echoes measured from the percolation zone.

Rignot, E. J.↗
DOE OSTI2026

Exploring sodium dynamics in the dilute oxygen regime of mixed oxy-sulfide NaPSiSO glasses

All-solid-state‑sodium batteries have great potential to lower the cost of grid-scale energy storage systems. However, they require a low cost, high conductivity solid electrolyte to become a reality. The electrochemical properties of mixed oxy-sulfide glasses NaPSiSO exceed those of either the pure sulfide or the pure oxide glass alone. In this work, we have tested models of sodium ion dynamics in the NaPSiSO series: zNa 2 S+(1-z)[(1-y)[(1-x)SiS 2 + xPS 5/2 ] + yNaPO 3 ] where 0.57 ≤ z ≤ 0.588, 0 ≤ x ≤ 0.35, and 0.15 ≤ y ≤ 0.35 using 23 Na nuclear magnetic resonance to directly probe the local fields and their fluctuations at the nuclear sites. We determined the distribution of activation energies in these materials arising from distributions of local environments and provided quantitative information about the connectivity of the ionic pathways through which the sodium ions move. We found that the sodium ion conduction dynamics were described as an ionic percolation through a Gaussian distribution of energy barriers and that the composition dependence of the ionic conductivity in the dilute oxygen regime is well-described by the Nernst-Einstein relation modified for percolation.

25 ENERGY STORAGE↗
DOE OSTI2025

Tuning Cation (Dis)Order in Cr‐Based Li‐Excess Oxide Cathode Materials to Improve Li+ Transport Properties

Abstract Li‐excess disordered rocksalts (DRXs) hold promise as next‐generation cathodes for Li‐ion batteries due to their high capacity and energy density, along with the potential to eliminate the need for Co and Ni. However, due to their low Li + diffusivity, DRXs need to be pulverized into nanoparticles to achieve high performance. Herein, a new strategy for overcoming this limitation is demonstrated, involving the design of as‐synthesized partially disordered oxides with a structure that lies between ordered layered and fully disordered, exhibiting a varying degree of local (dis)order. This unique structure activates new Li + diffusion channels, improving percolation and transport properties. This strategy allows a large content of Li + to be accessed in material with large, micron‐sized particles through highly reversible Cr 3+/6+ and O redox, yielding a first discharge capacity of 286 mAh g −1 (881 Wh kg −1 ). The Li + percolation network is further improved by substituting Ti with a mixture of multiple metals, which appears to locally decrease the migration barrier through lattice distortion. Proper tuning of the chemical composition, especially the content of metals with empty d orbitals, is established as a crucial factor for controlling the degree of disorder and mitigating voltage fade and hysteresis growth upon cycling.

Moździerz, Maciej↗
DOE OSTI2025

Origin of enhanced performance when Mn-rich rocksalt cathodes transform to δ -DRX

Most Mn-rich cathodes are known to undergo phase transformation into structures resembling spinel-like ordering upon electrochemical cycling. Recently, the irreversible transformation of Ti-containing Mn-rich disordered rock-salt cathodes into a phase — named δ — with nanoscale spinel-like domains has been shown to increase energy density, capacity retention, and rate capability. However, the nature of the boundaries between domains and their relationship with composition and electrochemistry are not well understood. In this work, we discuss how the transformation into the multi-domain structure results in eight variants of Spinel domains, which is crucial for explaining the nanoscale domain formation in the δ -phase. We study the energetics of crystallographically unique boundaries and the possibility of Li-percolation across them with a fine-tuned CHGNet machine learning interatomic potential. Energetics of 16 d vacancies reveal a strong affinity to segregate to the boundaries, thereby opening Li-pathways at the boundary to enhance long-range Li-percolation in the δ structure. Defect calculations of the relatively low-mobility Ti show how it can influence the extent of Spinel ordering, domain morphology and size significantly; leading to guidelines for engineering electrochemical performance through changes in composition.

Anand, Shashwat↗
DOE OSTI2025

Unraveling the pathway towards superionic transport in polymer electrolytes

Ionic transport in polymers is critical for Li-ion batteries, fuel cells, flow batteries and many other energy storage and conversion technologies. A significant enhancement of ion conductivity in polymers may be achieved through an increase in the polarity of side chains and their self-organization into specific morphologies, which can potentially act as percolated ionic structures. However, higher polarity increases attractive interactions within a polymer matrix and slows down its segmental dynamics, which conversely hinders ionic transport. To overcome this tradeoff, we designed the functionalization of a Li salt-doped polymer matrix by tailored amounts of zwitterionic (ZI) groups. Our results suggest the emergence of a self-assembled percolation conductivity regime above a specific ZI concentration, in which ion hopping decouples from segmental dynamics by up to ten orders of magnitude. Consequently, in the highly concentrated ZI regime, our polymeric materials exhibit in their glassy state energy barriers for ion hopping similar to, or even smaller than, those reported for superionic ceramics. Our study also reveals that ion dynamics in the poly(zwitterion) with all monomers carrying ZI groups is significantly faster than that of a monomeric ZI compound, although the latter has much faster structural relaxation. Furthermore, this result highlights the crucial role played by the local morphology on the ion transport of polymer electrolytes and opens a new pathway for the design of superionic polymers, significantly expanding the current limited portfolio of solid-state electrolytes for energy applications.

Ion conductivity↗
DOE OSTI2025

Larger than uncorrelated vacancy diffusion contributions in chemically disordered crystalline materials

Arising from a variational approach to compute diffusion coefficients, we introduce “superkinetic kinosons”, which are contributions to the diffusion flux of mobile defects and atomic species by single jump mechanisms that can exceed their uncorrelated values. In vacancy-mediated diffusion in crystalline random and high entropy materials, these contributions primarily arise from intermittent jumps of slow moving atoms that remove the vacancy out of correlation traps and enable a quantification of such effects on the diffusion coefficients. They can be significant contributors to diffusion, even when faster moving atoms have infinite percolation networks. Furthermore, their existence and importance provide fundamental insights about underlying aspects of vacancy diffusion such as formation of localized correlation traps and deviation of vacancy diffusion from percolation behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Kinetic Monte Carlo Framework for Coupled Degradation and Dehydration of Anion Exchange Membranes

Kinetic Monte Carlo (kMC) simulations, augmented with temporal-acceleration schemes, can efficiently handle stiff reaction-transport networks when fast processes rapidly relax to quasi-equilibrium on a fixed lattice. However, in glassy anion-exchange membranes (AEM), rare and irreversible chemical degradation events continuously reshape the nanoscale morphology, and the associated hydration and transport degrees of freedom remain far from a well-defined local equilibrium. This combination of evolving state space and nonequilibrated fast dynamics lies outside the scope of existing kMC acceleration frameworks. Here, to address this challenge, we introduce an auxiliary-particle kinetic Monte Carlo (AP-kMC) scheme. In AP-kMC, short-lived mobile particles spawned at degradation sites execute hop, water-elimination, and decay moves, enforcing rapid local relaxation of the hydration structure while preserving the stochastic rules of kMC. Parameterized with molecular-dynamics morphologies and experimental solution degradation kinetics, AP-kMC reproduces the evolution of ion-exchange capacity, water uptake, and conductivity, and reveals a feedback loop in which poorly hydrated sites degrade first and each degradation event induces further local dehydration. The resulting thinning and fragmentation of water channels cause loss of hydrophilic percolation and abrupt conductivity collapse well before complete charge loss. AP-kMC thus reframes AEM durability as a coupled degradation–drying–percolation problem and provides a transferable strategy to simulate reactive, out-of-equilibrium polymer electrolytes where local solvation controls reactivity.

organic↗
DOE OSTI2024

Vacancy-Mediated Increases in Brine–Salt Surface Energies

Salt formations have been explored for the permanent isolation of spent nuclear fuel based on their high thermal conductivity, self-healing nature, and low hydraulic permeability to brine flow. Vacancy defect concentrations in salt complicate fracture mechanics not driven by dislocation dynamics and can influence the resulting surface structure. Classical molecular dynamic simulations were used to simulate tensile testing of salt crystals (halite) with vacancy defect concentrations of up to 0.5 defects/nm 3 . Increasing defect concentrations resulted in a decrease in ultimate tensile strength and fracture surface energies, driven by increased surface roughness rather than changes in the amount of surface area. Brine–salt surface energies of the fractured surfaces were 0.22 to 0.26 J/m 2 , significantly higher than values reported for atomically flat (100) surfaces at the same brine composition. This change in surface energy increased the brine–salt dihedral angle by ~27°. The dihedral angle threshold for percolation in salt is 60°, and a 27° increase due to rough fracture surfaces identifies a reduction in porosity percolation and a decrease in salt permeability. Therefore, bedded salt and salt domes may be even more stable than those previously predicted from dihedral angle calculations.

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

Investigating Morphology and Diffusion in Simulations of Precise Anion-Conducting Polymers

Using atomistic molecular dynamics simulations, we investigate the morphology and transport properties of a new class of polymers which are functionalized with quaternary ammonium groups for use as anion exchange membranes. The polymers are precision polyolefins with either a trimethylammonium (p5CNMe3) or a dimethyl-hexyl ammonium (p5CNMe2Hx) pendant group at every fifth carbon along a polyethylene backbone. Simulations are performed at hydration levels of 5, 10, 15, and 20 water molecules per ammonium group. The hydrated polymers form nanoscale, percolated hydrophilic domains (water channels) in the hydrophobic polymer matrix that become wider with increasing water content. Water and hydroxide anion diffusion coefficients also increase with increasing water content. The morphology of the water domains is similar in both polymers, while the diffusion coefficients are somewhat lower in p5CNMe2Hx at fixed water content. Furthermore, the diffusion coefficients in both polymers fall on the same curve as a function of the fractal dimension of the percolated water channels, which appears to be a useful scalar measure of the effects of the nanoscale morphology on water and hydroxide anion transport.

Anion exchange membrane↗
DOE OSTI2025

Defect-Mediated Diffusion Pathways in Spodumene Accelerate Lithium Transport

Lithium extraction from naturally occurring α-spodumene is hindered by poor lithium diffusivity, necessitating high-temperature phase transformation to a low-density β polymorph. Although β spodumene exhibits up to 5 orders of magnitude higher lithium-ion diffusivity, both phases have diffusion activation energies between 0.8 and 1 eV, indicating that polymorph density is not the controlling factor over diffusivity. We show that aluminum vacancies facilitate lithium-ion diffusion in α-spodumene by reducing the migration barrier from 2.4 to 0.9 eV. Bond valence site energy and nudged elastic band calculations show a new lithium local minimum site which promotes a one-dimensional percolation network by reducing the lithium intersite distance from 4.5 Å to 2.9 Å. However, aluminum vacancies are energetically unfavorable to percolate through the whole structure, resulting in very low net lithium diffusivity and highlighting the critical role of nonstoichiometric defects in facilitating lithium transport in rigid aluminosilicate structures.

Chemical structure↗
DOE OSTI2025

Tetrahedral Lithium Stuffing in Disordered Rocksalt Cathodes for High-Power-Density and Energy-Density Batteries

Li-rich cation-disordered rocksalt (DRX) materials introduce new paradigms in the design of high-capacity Li-ion battery cathode materials. However, DRX materials show strikingly sluggish kinetics due to random Li percolation with poor rate performance. Here, in this study, we demonstrate that Li stuffing into the tetrahedral sites of the Mn-based rocksalt skeleton injects a novel tetrahedron-octahedron-tetrahedron diffusion path, which acts as a low-energy-barrier hub to facilitate high-speed Li transport. Moreover, the enhanced stability of lattice oxygen and the suppression of transition metal migration preserve the efficacy of the Li percolation network during cycling. Overall, the tetrahedral Li stuffing DRX material exhibits high energy density (311 mAh g -1 , 923 Wh kg -1 ) and high power density (251 mAh g -1 , 697 Wh kg -1 at 1000 mA g -1 ). Our results highlight the potential to develop high-performance and earth-abundant cathode materials within the extensive range of rocksalt compounds.

Disordered Rocksalt Cathodes↗
DOE OSTI2025

Tape-in-tape-out (TITO): a new approach for in-situ contact resistance measurements in high temperature superconducting CORC® cables

One of the ongoing development challenges with ReBCO high-temperature superconducting (HTS) cables is normal zone initiation and local heating, which is associated with over-critical current flowing through and around local performance reductions in individual tapes. Although inter-tape contact resistances are well-reported for individual tapes and HTS cables, these measurements are a challenge in CORC® cables as current percolates through much of the cable. In this work, developments in tape-in-tape-out (TITO) automated individual tape powering experiments are presented, and a simplified modeling approach for current percolation in CORC® cables is developed. An optimization is formulated to fit the model parameters to a large set of TITO experiments on a single cable, allowing the layer-dependent inter-tape contact resistance to be extracted. Measurements are presented and discussed for a straight CORC® cable and a cable bent to a 152 mm and 76 mm radius. The approach provides both quantitative insight and facilitates better understanding of current sharing phenomena in HTS cables.

Teyber, Reed↗
DOE OSTI2025

Interplay of superconducting, metallic, and crystalline states of composite fermions at 𝜈 = $\frac{1}{6}$ in wide quantum wells

Evidence for developing fractional quantum Hall effect (FQHE) at filling fraction 𝜈 = 1/6 and 1/8 was recently reported in wide GaAs quantum wells [Wang et al., Phys. Rev. Lett. 134, 046502 (2025)]. In this article, we theoretically investigate the nature of the state at 𝜈 = 1/6 as a function of the quantum well width and the density by considering composite-fermion (CF) crystals, CF Fermi sea, and various kinds of paired CF states. The 𝑓-wave paired state has the lowest energy among the paired CF states. However, for parameters of interest, the energies of the CF crystal, the CF Fermi liquid, and the 𝑓-wave paired CF state are too close to distinguish. We, therefore, predict that 𝑖𝑓 the FQHE at 𝜈 = 1/6 is experimentally confirmed, this state would be an 𝑓-wave paired state of CFs, which can be verified by measurement of its thermal Hall conductance. Exact diagonalization studies on clean systems with up to eight electrons show that the ground states at 𝜈 = 𝑛/(6⁢𝑛 ± 1) are incompressible for all widths and densities we have considered, and are well described by the corresponding Laughlin and Jain states. We propose a phase diagram for large quantum well widths and densities in which at zero disorder, incompressible FQHE states are stabilized at 𝜈 = 𝑛/(6⁢𝑛 ± 1) and 𝜈 = 1/6, but in between these fillings the CF crystal is stabilized. We also present a qualitative discussion on the effects of disorder and propose a schematic phase diagram based on it. With disorder, which creates a spatial variation in the filling factor, two regimes are identified: (i) for small disorder, when the incompressible states percolate at the special fillings, FQHE with quantized Hall plateaus and vanishing longitudinal resistance should occur; and (ii) for larger disorder, when the CF crystal percolates, the longitudinal resistance rises with decreasing temperature but the domains of FQHE liquid produce minima at the special filling factors. Here, experiments are consistent with the latter scenario. We also mention a possible connection of the phase diagram presented here to a puzzling behavior observed for the fractional quantum anomalous Hall effect in pentalayer graphene.

Composite fermions↗
DOE OSTI2025

SOC Microstructural Property Estimator

This pre-trained ML model is a tool that uses basic compositional parameters for porous solid oxide cell (SOC) electrodes - the phase fractions and mean particle/pore diameters – as inputs and uses them to estimate additional electrochemical performance parameters: active (i.e., connected) TPB density, all tortuosity factors, and phase pair specific interfacial areas. The electrode is assumed to be composed of two solid phases and a pore phase. The property calculations are performed using neural network regression models trained on a large bank of synthetic electrode microstructural data that NETL has generated using the program DREAM3D (that bank is also hosted on EDX: https://edx.netl.doe.gov/dataset/soc-synthetic-microstructure-bank). This means the generated parameters are based on training from actual measured properties from 3D microstructures, not estimated from geometric simplifications. This tool was developed and is intended to replace percolation theory calculations in models that use hypothetical electrode properties. An example use case would be running SOC performance simulations across a parametric sweep of electrode designs (e.g., varying phase fractions and particle sizes) and assessing how it impacts the electrochemical performance of the SOC. Within the parameter space of the training data (statistics of that parameter space is provided in the readme file), this model achieves sub-5% mean absolute percent errors, an order of magnitude less error than percolation theory across the same parameter space. However, be aware that this tool was developed with parametric simulations in mind, and users are encouraged to assess accuracy for their own specific use case rather than taking accuracy metrics at face value. More info, including a usage guide, is in the included readme file. This tool should be cited with the DOI number provided.

Electrode Microstructure↗
NASA NTRS1988

Multiphase flow in porous media

A development history and current status evaluation are presented for the theory of permeability and percolation. The microscale phenomena treated in this field have proven difficult to analyze due both to their tortuous geometry and the influence of capilarity. Capilary effects may be not only important but predominant, and are differentiated into those at the fluid-fluid interface, and those involving the existence of a contact line between the solid substrate and this interface. Percolation theory has been borrowed from physics and adapted to the two-phase engineering context.

Adler, Pierre M.↗
NASA NTRS1991

Measures of large-scale structure in the CfA redshift survey slices

Variations of the counts-in-cells with cell size are used here to define two statistical measures of large-scale clustering in three 6 deg slices of the CfA redshift survey. A percolation criterion is used to estimate the filling factor which measures the fraction of the total volume in the survey occupied by the large-scale structures. For the full 18 deg slice of the CfA redshift survey, f is about 0.25 + or - 0.05. After removing groups with more than five members from two of the slices, variations of the counts in occupied cells with cell size have a power-law behavior with a slope beta about 2.2 on scales from 1-10/h Mpc. Application of both this statistic and the percolation analysis to simulations suggests that a network of two-dimensional structures is a better description of the geometry of the clustering in the CfA slices than a network of one-dimensional structures. Counts-in-cells are also used to estimate at 0.3 galaxy h-squared/Mpc the average galaxy surface density in sheets like the Great Wall.

De Lapparent, Valerie↗