Search NASA⌕ Search

SEARCH · Search NASA

Results for “Permeability”

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

Effectiveness of nature-based solutions to reduce flooding in Quad Cities Metro Area (QCMA) using SWMM-HEC based flood model

Nature-based solutions (NbS) have gained significant attention as strategies for addressing urban environmental challenges, particularly since the establishment of the UN Sustainable Development Goals (SDGs) for 2030. However, the current research on NbS for urban flood management lacks comprehensive methodological approaches for identifying suitable areas and evaluating their effectiveness across different urban settings. Here, this study attempts to fill this gap by proposing a methodological framework integrating multi-criteria analysis with a SWMM-HEC-based hydrologic and hydraulic (HH) model to assess the suitability of NbS for the Quad Cities Metro Area (QCMA), consisting of Davenport, Bettendorf, Moline, and Rock Island. Eight NbS options-green roofs, rain gardens, infiltration trenches, permeable pavements, vegetative swales, dry detention basins, retention ponds, and rain barrels/cisterns - were considered based on volumetric efficiency and runoff reduction efficiency. The study reveals that implementing the proposed NbS could have substantially reduced flood depths in key historical flood events by 21% in 1993, 15% in 2008, 16% in 2011, 23% in 2014, 40% in 2019, and 10% in 2023. The findings highlight a critical trade-off between peak runoff and NbS implementation: while NbS effectively reduce flood impacts, they also enhance volumetric efficiency by approximately 43%. In high-density areas of the QCMA, flood depth reductions of around 20% suggest that NbS are a viable solution for dense urban environments with limited space. This shows the potential for integrating NbS into existing infrastructure, offering a promising approach for cities facing increasing flooding risks. The proposed methodology provides a practical framework for incorporating NbS into urban stormwater management, addressing gaps in optimizing NbS performance, and offering a pathway to scale their application in other urban areas with various environmental and social contexts.

CMIP6↗
DOE OSTI2025

Extending TOUGH + HYDRATE with a parallel particle transport simulator: numerical investigation of sand production during gas production from hydrate deposits

A new parallel code for simulating particle transport in porous media is integrated with the TOUGH + HYDRATE simulator to investigate sand production associated with gas production from unconsolidated gas hydrate-bearing sediments (HBS). Here, the parallel coupled simulator is named THMPT and uses the integral finite difference method to describe the Darcian and non-Darcian flow of fluids and heat transport, the finite element method to describe the associated geomechanical changes, and the discrete element method to track the trajectory of individual sand particles within the HBS. The THMPT simulator is written in Fortran, incorporates multiple optimized algorithms, and can comprehensively address the coupled flow, thermal, chemical, geomechanical, and particle transport processes that characterize the system behaviors during gas production from HBS. The simulator can capture all processes involved in sand particle transport in porous media, including sand detachment, collision, clogging (i.e., bridging), and migration. A benchmark case study of sand production in the course of depressurization-induced gas production from a representative HBS reveals various distinct microscopic particle migration mechanisms and the adverse impact of sand particle detachment, transport, and clogging. The numerical investigation also examines the effect of bottomhole pressure on mitigating sand production. The simulation results indicate that sand clogging near the wellbore significantly reduces permeability, decreasing gas production by at least 50%. Lastly, the efficiency of gravel packing in mitigating sand production is numerically evaluated, revealing that the structure of the porous media appears to profoundly influence the macroscopic motion behavior of sand particles and sand clogging characteristics.

discrete element method↗
DOE OSTI2024

Splitting tensile strength of shale cores: intact versus fractured and sealed with ureolysis-induced calcium carbonate precipitation (UICP)

Ureolysis-induced calcium carbonate precipitation (UICP) is a biomineral solution where the urease enzyme converts urea and calcium into calcium carbonate. The resulting biomineral can bridge gaps in fractured shale, reduce undesired fluid flow, limit fracture propagation, better store carbon dioxide, and potentially enhance well efficiency. The mechanical properties of shale cores were investigated using a modified Brazilian indirect tensile strength test. An investigation of intact shale using Eagle Ford and Wolfcamp cores was conducted at varying temperatures. Results show no significant difference between shale types (average tensile strength = 6.19 MPa). Eagle Ford displayed higher strength at elevated temperature, but temperature did not influence Wolfcamp. Comparatively, cores with a single, lengthwise heterogeneous fracture were sealed with UICP and further tested for tensile strength. UICP was delivered via a flow-through method which injected 20–30 sequential patterns of ureolytic microorganisms and UICP-promoting fluids into the fracture until permeability reduced by three orders of magnitude or with an immersion method which placed cores treated with guar gum and UICP-promoting fluids into a batch reactor, demonstrating that guar gum is a suitable inclusion and may reduce the number of flow-through injections required. Tensile results for both delivery methods were variable (0.15–8 MPa), and in some cores the biomineralized fracture split apart, possibly due to insufficient sealing and/or heterogeneity in the composite UICP-shale cores. Notably in other cores the biomineralized fracture remained intact, demonstrating more cohesion than the surrounding shale, indicating that UICP may produce a strong seal for subsurface application.

58 GEOSCIENCES↗
DOE OSTI2024

Phase-field modeling of thermally-grown oxide and damage evolution in environmental barrier coatings

Silicon carbide-based ceramic matrix composites protected by environmental barrier coatings (EBCs) present a promising materials solution for next-generation gas turbines. Developing more robust and efficient EBCs is therefore of significant technological importance. During the service in high-temperature oxidative environments, there is a thermally grown oxide (TGO) layer, spontaneously formed in the EBC system. TGO is recognized as a critical factor for the degradation and failure of EBCs, yet the detailed mechanisms of TGO growth and its effect on EBC failure remain unclear. In this study we develop a comprehensive chemo-mechano-phase-field model to simulate growth of the TGO in EBCs, factoring in creep and deformation, and especially the cracking behaviors. The volume expansion due to TGO growth and the resulting large inelastic deformation are addressed by using our recently developed, so-called incremental realization of inelastic deformation (IRID) algorithm, in combination with an adapted Hu-Chen spectral solver for elasticity. Simulations of TGO growth are performed considering different growth modes of TGOs determined mainly by the ratio of oxidant permeability in the topcoat to that in the TGO itself. Large-scale three-dimensional (3D) simulations are performed to model the formation of interconnecting vertical/channel cracks (often called ‘mud cracks’). The simulated crack morphology are in excellent agreement with the experimental observations from the literature. The simulations also provide insights into the cracking of EBCs and its dependence on the structure and constituent properties of the coating system. Furthermore, these results demonstrate the developed damage model can be a useful tool for design of more durable EBCs.

Cracking↗
DOE OSTI2026

Stochastic fracture generation and thermo-hydro-mechanical modeling in an equivalent continuum framework for enhanced geothermal systems

Enhanced geothermal systems (EGS) involve fracturing low permeability material to establish well connectivity and then injecting and circulating fluid into the fractured subsurface for geothermal power production. Changes in fracture aperture from contraction of the cooling matrix rock may alter network connectivity and risk thermal short-circuiting. Thermo-hydro-mechanical (THM) models are a useful tool to study these processes. However, as fracture networks are complex, and data may be limited, fracture networks in THM models are often stochastically generated. Given reliance on stochastic fracture networks and THM modeling to represent the subsurface and assess productivity of EGS, increased understanding of the influence of such statistically derived fracture networks on flow and heat transport in THM models is needed. Here, a new fracture process model is developed in the reactive transport code PFLOTRAN to stochastically generate fracture families and simulate changes in fracture aperture over time due to temperature changes of the rock matrix. Sixty-four different fracture networks ranging from well to poorly-connected, are modeled in PFLOTRAN with and without mechanical processes (THM vs TH). Results indicate that for well-connected fracture networks, thermal short-circuiting is less of a concern due to the abundance of available alternative flowpaths. For poorly-connected fracture networks, inclusion of mechanical processes showed steep thermal drawdown coincident with increase in fracture aperture along developing colder flowpaths, demonstrating the risk of thermal short-circuiting. Simulations with additional, larger fractures engineered to establish connectivity in a poorly-fractured subsurface, indicate that while stochastic variation of fracture orientation of the background network had limited influence, such variation in the engineered fractures significantly affected flow and heat transport.

Discrete fracture networks (DFN)↗
DOE OSTI2024

Scalable and compact magnetocaloric heat pump technology

Magnetocaloric heat pumping (MCHP) promises to be more efficient than traditional vapor compression while also eliminating the deleterious effects of gaseous refrigerants. While MCHP devices have shown the temperature spans and efficiencies needed for different heating and cooling applications, they struggle to become commercially viable due to their large size and mass, and resultant high cost. This paper evaluates a baseline MCHP device and explores methods to boost its system power density (SPD). The key components of the baseline system are the gadolinium packed-particle bed active magnetic regenerator (AMR) and a magnetic source composed of permanent magnets and high permeability magnetic steel. To enhance the SPD, the paper evaluates maximizing the AMR volume, opting for first-order magnetocaloric materials, optimizing the magnet and AMR geometry, and reducing the size of magnets and magnetic steel parts. At larger thermal powers, increasing the AMR diameter and the number of magnetic poles were evaluated. Using finite element models, solid models, and estimates of magnetocaloric material performance, thermal powers ranging from 37 W to 44 kW at a nominal 10 K temperature span were projected, and SPD was estimated to improve from 6 W/kg to 81 W/kg. Neglecting end effects, an upper limit of 114 W/g is estimated. Compared to SPD of off-the-shelf compressors with similar environment temperatures, MCHP power density using gadolinium is competitive up to roughly 200 W of cooling power. This is extended to 1 kW when using LaFeSi alloys and up to 3 kW in the limiting case. In conclusion, these results indicate that the performance and mass of MCHP can match that of compressors, which is a critical step toward cost-competitive magnetocaloric technology.

42 ENGINEERING↗
DOE OSTI2024

Rational design of semi-interpenetrating network based on hyperbranched polyglycerol grafted MXene/polyurethane–epoxy for compressed hydrogen storage

Here, hydrogen (H 2 ) energy has emerged as a promising alternative in the automotive sector to replace fossil fuel; however, its storage with volumetric efficiency remains a challenge. The latest type IV storage vessel featuring a polymeric liner with composite overwrap incurs hydrogen saturation, which eventually causes failure. To prevent the H 2 dissolution in the liner, we have developed a barrier coating comprised of polyurethane/epoxy semi-interpenetrating network (S-IPN) combined with hyperbranched polyglycerol grafted MXene (h-MXene). The hyperbranched structure facilitated the dispersion of MXene in the coating solution by enhancing the exfoliation and improved the polymer filler interaction by forming covalent and H-bonding through end-terminal hydroxyl groups. Leveraging the dense network of S-IPN combined with the dispersed layer structured h-MXene significantly increased the tortuosity of the spray-coated barrier film applied to the nylon 6 liner. The coating exhibited a 90 % reduction in the H 2 gas permeability at only 2 wt% h-MXene concentration, which further improved to above 98 % at 10 wt% loading. Additionally, the h-MXene considerably improved the adhesion with the liner even in a highly stretched condition, signifying the durability of the coating under cyclic pressurization and depressurization of the storage vessel.

08 HYDROGEN↗
DOE OSTI2024

Dielectric and magnetic properties of microwave-absorbing FeAl x O y catalysts fabricated via solution combustion synthesis

Iron-based alumina (FeAl x O y ) nanocomposites are microwave-absorbers and catalysts, which makes them promising for emerging microwave-assisted thermocatalytic technologies. Solution combustion synthesis (SCS) has been used to synthesize FeAl x O y powders, and prior work has demonstrated that adjusting SCS parameters significantly changes phase composition and specific surface area of the products. However, it is unclear how synthesis parameters affect their microwave-absorbing properties, which are essential for optimizing microwave-assisted technologies. To address this challenge, in the present work, twelve different FeAl x O y products were synthesized at different combinations of the SCS parameters such as two fuels (citric acid and glycine), two heating modes (hotplate and muffle furnace), and three Fe:Al molar ratios (2:1, 1:1, 1:2). Dielectric and magnetic properties of the products were characterized using a network analyzer and a vibrating sample magnetometer. Based on the measured permittivity and permeability, penetration depth and reflection loss were calculated as a function of frequency and bed thickness. The products were heated by microwaves at 2.45 GHz and then examined with X-ray diffraction (XRD) analysis. For all products, the magnetic saturation was lower than for bulk iron oxides because of the small crystallite size and aluminum substitution. The use of glycine induced high dielectric losses and enabled fast microwave-heating rates compared to citric acid. Higher Fe:Al ratio also led to higher dielectric and magnetic losses. With glycine fuel, SCS in a furnace induced larger penetration depth and lower microwave absorption than SCS on a hotplate. The minimization of reflected power was more sensitive to the thickness of the product bed than to the frequency of the electromagnetic field. Post-heating XRD analysis revealed different phase transformations in the FeAl x O y powders depending on the SCS parameters. As a result, an FeAl x O y material, synthesized via incipient wetness impregnation, lacked magnetic losses and did not heat well as compared to the SCS products.

Combustion synthesis↗
DOE OSTI2025

Interaction between dissolution and precipitation during olivine carbonation: Implications for CO 2 mineralization

Large-scale carbonation of olivine is considered a promising approach for in situ mineral carbonation, offering a permanent and stable method for CO 2 storage. A critical aspect of this process is understanding how dissolution and precipitation interact, as this could drive fracturing and enhance further reactions. In this study, we conducted carbonation experiments on olivine using CO 2 -saturated aqueous solutions of NaHCO 3 and NaCl. Two experimental setups were used: one representing an open geochemical system and the other a closed system, corresponding to reaction-limited and flow-limited scenarios, respectively. Further, post-reaction textural analysis using scanning electron microscopy (SEM) revealed surface coatings of reaction products in the closed system, while etch pits and etch channels were prevalent in the open system. Although no direct evidence of reaction-driven fracturing was observed, etch pits and etch channels may serve as initiation points for subcritical crack formation and growth, potentially maintaining permeability and exposing new unreacted surfaces. Using linear elastic fracture mechanics (LEFM) model, we estimate that microcracks could propagate under a pressure of 0.1 GPa if reaction products accumulate within the etch pits. Our findings offer new insights into the mechanisms governing olivine carbonation.

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

A numerical study of process complexity in permafrost dominated regions

Numerical modeling of permafrost dynamics requires adequate representation of atmospheric and surface processes, a reasonable parameter estimation strategy, and site-specific model development. The three main research objectives of the study are: (i) to propose a novel methodology that determines the required level of surface process complexity of permafrost models by conducting parameter sensitivity and calibration, (ii) to design and compare three numerical models of increasing surface process complexity, and (iii) to calibrate and validate the numerical models at the Yakou catchment on the Qinghai-Tibet Plateau as an exemplary study site. The calibration was carried out by coupling the Advanced Terrestrial Simulator (numerical model) and PEST (calibration tool). Simulation results showed that (i) A simple numerical model that considers only subsurface processes can simulate active layer development with the same accuracy as other more complex models that include surface processes. (ii) Peat and mineral soil layer permeability, Van Genuchten alpha, and porosity are highly sensitive. (iii) Liquid precipitation aids in increasing the rate of permafrost degradation. (iv) Deposition of snow insulated the subsurface during the thaw initiation period. We have developed and released an integrated code that couples the numerical software ATS to the calibration software PEST. The numerical model can be further used to determine the impacts of climate change on permafrost degradation.

Calibration↗
DOE OSTI2026

Effects of oxide surface chemistry on diffusioosmosis

Diffusioosmosis is the movement of fluid induced by gradients in solute concentration. Recent studies suggest that in low permeability rocks, it may be the dominant mode of reactant transport and thus control rates of diagenesis, which cannot be adequately explained by pressure-driven flow alone. In this paper we investigate how the equilibrium between an oxide mineral and the surrounding fluid phase influences the diffusioosmotic velocity. We have developed a theory for the case of a thin double layer, where the Debye length is smaller than the characteristic pore size. Several factors contribute to the total fluid velocity: the chemical structure of the mineral surface, the electrolyte type and concentration gradient, and the solution pH. Individual factors can act in concert or in opposition, leading to widely varying magnitudes and directions of the velocity. The numerical results are within the range of the limited experimental data. Our results highlight how surface charging and surface complexation impact the flow, and how they depend on pH.

Diagenesis↗
DOE OSTI2025

Applying a Compact Porous Media Model to Numerically Derive Resistance Coefficients for Lattice Structures

Additive Manufacturing allows for exploring various geometries to achieve specific engineering criteria. Lattices are one geometry with unique properties, including being periodically repeating structures which allow flow through them to be represented as a porous media according to Darcy-Forchheimer equations. These equation’s coefficients are generally experimentally derived, but this work demonstrates the ability to numerically derive them with CFD. Simulations were performed using three-dimensional stead state Reynolds-averaged Navier-Stokes with a k-ω Shear Stress Transport turbulence model using Ansys Fluent. Three lattice geometries were investigated and drag coefficients were derived. The method was validated against externally published data for similar geometries demonstrating strong agreement, and grid convergence for all simulations was calculated with a Grid Convergence Index method. Wall roughness is demonstrated to have a non-negligible impact on results and roughness values are considered for the primary focus Octahedral geometry where both smooth wall and rough wall coefficients were derived. The porosity coefficients for the Octahedral geometry at 1.0 [m/s] were found to be 2.89×10 6 and 2.90×10 6 [1/(Pa*m*s)] for the permeability coefficients, 6.37×10 1 and 5.44×10 1 [m 2 /kg] for the inertial resistance coefficients, and with a max pressure drop of 5116.7 [Pa] and 4429.5 [Pa] for the smooth walls and rough walls, respectively. The derived numerical method enables rapid exploration and optimization of new lattice designs for diverse engineering applications.

42 ENGINEERING↗
DOE OSTI2024

Crack opening calculation in phase-field modeling of fluid-filled fracture: A robust and efficient strain-based method

The phase-field method has become popular for the numerical modeling of fluid-filled fractures, thanks to its ability to represent complex fracture geometry without algorithms. However, the algorithm-free representation of fracture geometry poses a significant challenge in calculating the crack opening (aperture) of phase-field fracture, which governs the fracture permeability and hence the overall hydromechanical behavior. Although several approaches have been devised to compute the crack opening of phase-field fracture, they require a sophisticated algorithm for post-processing the phase-field values or an additional parameter sensitive to the element size and alignment. Here, we develop a novel method for calculating the crack opening of fluid-filled phase-field fracture, which enables one to obtain the crack opening without additional algorithms or parameters. Here we transform the displacement-jump-based kinematics of a fracture into a continuous strain-based version, insert it into a force balance equation on the fracture, and apply the phase-field approximation. Through this procedure, we obtain a simple equation for the crack opening which can be calculated with quantities at individual material points. We verify the proposed method with analytical and numerical solutions obtained based on discrete representations of fractures, demonstrating its capability to calculate the crack opening regardless of the element size or alignment.

58 GEOSCIENCES↗
DOE OSTI2026

Peridynamic modeling of cementitious materials for nuclear waste management

Radioactive and hazardous waste generated from fuel processing plants, nuclear reactors, and hospitals, requires effective management strategies. Cementitious materials are widely applied for these needs, serving as structural materials, reactive barriers, or waste forms. In these applications cracking poses a significant risk to performance, driven by inconsistent shrinkage of the components and the varying strength and permeability of their interfaces. Traditional modeling approaches face challenges in representing the complex fracture behavior of cementitious materials due to the reliance on spatial derivatives and difficulties with mesh generation. Here, to overcome these limitations, we employ peridynamics, a novel continuum mechanics formulation that uses integrals to describe mechanical equilibrium, avoiding discontinuities associated with traditional methods. Through incorporation of a bi-linear softening model and quasistatics, an experimentally validated model for Portland cement concrete samples was created. Mechanical parameters, including compressive strength and elastic modulus were validated and variation due to aggregate packing was evaluated. Additionally, sensitivity analysis of the peridynamic parameters to the bulk material properties was established. The results lay the groundwork for evaluating the impact of unique conditions of cementitious waste forms that can be assessed to improve the reliability of waste management strategies.

Aggregates↗
DOE OSTI2025

The impact of capillary heterogeneity on CO 2 flow and trapping across scales

Capillary heterogeneity has been identified over the last decade as a key control on subsurface CO 2 flow behavior during geological CO 2 sequestration. These heterogeneities can be formed in all sedimentary rocks, ranging from slight variations in the sand grain sizes to extensive sequences of interbedded sands, shales, and limestones. Capillary heterogeneity has been largely, although not entirely, overlooked in subsurface flow modeling because it is assumed to only directly influence fluid redistribution over scales of centimeters to meters. However, even small-scale fluid movements can result in dramatic impacts on the mobility and trapping of the CO 2 over kilometers. Therefore, neglecting capillary heterogeneity at multiple scales could potentially lead to errors in modeling and predicting field-scale plume migration. In this review paper, we aim to provide a consistent overview to (1) establish that capillary heterogeneity can have a major impact on CO 2 plume migration, (2) establish the respective length scales at which capillary heterogeneity matters, and (3) provide guidance for numerical modeling. This review covers pertinent literature and extracts key observations from the core to the field scales. Experimental studies have shown that millimeter-decimeter scale capillary heterogeneity can cause the so-called capillary heterogeneity trapping in addition to pore-scale residual trapping. Even at such a small scale, capillary heterogeneity can already lead to complex upscaled constitutive relationships, such as flow-rate dependent and anisotropic relative permeability, which affects field-scale CO 2 migration even when field-scale heterogeneities are present. Under gravity-dominated flow regimes, centimeter-meter scale capillary heterogeneity can entrap a significant amount of CO 2 at field scale, not just after imbibition but also during drainage. In certain cases, the presence of capillary heterogeneity can even completely stop the vertical movement of the CO 2 plume, hence greatly reducing leakage risks. At meter-kilometer scale, the influence of capillary heterogeneity is more pronounced and can hinder or redirect CO 2 migration in both lateral and vertical directions. The impact of capillary heterogeneity across multiple spatial scales poses a great challenge in modeling CO 2 migration at field scale, because it is practically impossible to build a field-scale earth model with grid blocks at millimeter scale. We recommend a hierarchical modeling approach to address this challenge. At field scale, earth models are built to capture geological features and heterogeneities in high but still practical grid resolutions. For each facies or rock type of the field-scale model, high- resolution meter-scale “conceptual” models are built with millimeter-scale grid blocks to capture representative fine-scale bedding geometries and heterogeneities in various environments of deposition, bridging the gap from subcore scale to the size of a field-scale simulation grid block. Upscaling is then used to preserve the smaller-scale flow dynamics of various rock types in field-scale simulations. Here, future work is needed to (1) refine, improve, and validate the hierarchical modeling approach; (2) build libraries of fine-scale bedding models for facies in various environments of deposition; (3) quantify multiscale capillary heterogeneity effects under subsurface uncertainties; (4) gain learning from different storage formations; and (5) establish best practices that balance accuracy and computational speed.

Capillary heterogeneity↗
DOE OSTI2025

Geologic hydrogen: From natural occurrences to anthropogenic generation – A review of fundamentals, potential, challenges and prospects

Growing demand for hydrogen is exposing the environmental and economic limits of reforming-based and carbon-managed supply chains, while the scale-up of electrolytic capacity remains capital-constrained. Geologic hydrogen, defined as molecular H₂ generated and stored within the Earth's crust offers a complementary, potentially lower-cost resource, yet exploration is still ad hoc. This review (1) revisits a global inventory of confirmed hydrogen seeps and subsurface occurrences; (2) analyzes the controlling reactions, migration pathways, and trapping conditions governing these occurrences; (3) proposes a process-based geologic hydrogen system concept analogous to, yet distinct from, the petroleum system; and (4) evaluates potential geologic hydrogen systems within the United States as a representative case study. Here, we contrast natural systems powered by serpentinization, mantle degassing or radiolysis with anthropogenic systems that stimulate the same reactions or convert in-situ hydrocarbons. Stable hydrogen accumulations require generation rates that outpace combined physical, chemical and microbial losses; the Bourakébougou field (Mali) exemplifies a self-recharging, free-gas reservoir sustained by meteoric-water serpentinization beneath an efficient caprock. Prospective geologic hydrogen resources are likely to occur in regions where iron-rich lithologies, deep-seated faults, and low-permeability sealing formations coexist. Applying this principle, we highlight three promising hydrogen play types in U.S. geological terrains: ophiolite belts (Appalachian and Californian regions), the Midcontinent Rift and the Lake Superior banded‑iron formations. Multiphysics numerical models and positive-unlabeled machine-learning workflows help to accelerate play screening and de-risk future production; yet, reaction kinetics, stimulation strategies, and full techno-economic and life-cycle assessments remain pivotal knowledge gaps.

Anthropogenic hydrogen generation↗
DOE OSTI2025

Bridging interfacial properties and cell performance: A multiscale model for proton-exchange-membrane fuel cells

Here, to elucidate the impact of local interfaces on mass-transport resistance and overall cell performance of low-loaded proton-exchange-membrane fuel cells (PEMFCs), we present a multiscale modeling framework incorporating a novel modified agglomerate model. The model considers three distinct Pt-electrolyte interfaces: Pt on the carbon surface covered by either ionomer or water film and Pt inside carbon nanopores. Detailed mass-transport voltage-loss breakdowns reveal that coupled agglomerate-interface-scale mass transport dominates the mass-transport loss. The ionomer poisons the exterior-Pt surface through suppressing O 2 adsorption and intrinsic ORR activity, leading to low current-density performance. Conversely, interior-Pt interface enhances the kinetic performance but limits high current-density performance due to its low interfacial permeability. The exterior-Pt/water interface demonstrates superior kinetic performance and mass transport, though its practical implementation requires ensuring proton transport. By coupling the multiscale CL properties with ink parameters, the model identifies an optimal I to C ratio of approximately 0.5, a moderate value where the ionomer content is sufficient to guarantee proton transport without fully covering the Pt surface and forming large agglomeration, thus allowing the utilization of the Pt-water interface and avoiding high mass-transport loss. Overall, the model helps unravel limiting phenomena across different operating regimes and provides routes for optimizing performance.

Cell diagnostic↗
DOE OSTI2026

A collective trigger for widespread planetesimal formation revealed by accretion ages

The formation of planetesimals was an integral part of the cascading series of processes that built the terrestrial planets. To illuminate planetesimal formation, here we develop a refined thermal evolution model to calculate the formation ages of meteorite parent planetesimals. This model includes chemical reactions and phase changes during heating, as well as natural variations in the proportions of the constituent phases of these planetesimals. We find that the parent bodies of non-carbonaceous (NC) and carbonaceous (CC) iron meteorites start forming at very similar times ( ∼ 0.95 Myr after calcium-aluminium-rich inclusion [CAI] formation) and occupy overlapping time windows. NC and CC chondrite parent bodies formed later during non-overlapping periods. We combine these ages with proportions of isotopic end-members we recover from mixing models to construct records of motion throughout the protoplanetary disk. These records argue that NC and CC material traversed the barrier in the disk after ∼ 0.95 Myr after CAI formation. The onset of this motion coincided with planetesimal formation, indicating that the phenomenon that drove motion also triggered planetesimal formation. We argue that this feature also served as the semi-permeable barrier in the disk. Although its identity is uncertain, the effects this phenomenon had on the timing of planetesimal formation and motion through the disk can now serve as constraints on models of disk evolution. Models that reproduce these effects would elucidate the nature and implications of this phenomenon, which is key to unlocking a holistic model of terrestrial planet building.

Astronomy and AstroPhysics↗