Search NASA⌕ Search

SEARCH · Search NASA

Results for “pore-scale”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Pore-Scale Study on the Positive Feedback Between Stress and Porosity Caused by Pressure Solution in Porous Media

Pressure solution is an important process in the evolution of sedimentary rocks, which provide storage space for most of our petroleum resources. It directly influences the generation, migration, and storage of petroleum fluids in subsurface sedimentary rocks. Here, in this paper, we develop a pore-scale, mechanochemical model to demonstrate a possible positive feedback between the local porosity and pore surface stress, in which a higher local porosity causes a higher local pore surface stress, thus enhancing pressure solution and consequently further increasing the local porosity. Pore surface stress represents stress on a solid grain adjacent to a pore. Specifically, the pore-scale, mechanochemical model directly simulates the stress distribution over solid and pore surfaces using a finite element model. The dissolution of solids at the solid-pore interfaces under a far-from-equilibrium condition is simulated using a first-order kinetics model that accounts for the local stress distribution. The updated pore geometry, caused by pore surface dissolution, is then used in the stress simulation in the next numerical iteration. Two types of porous media, the Oriskany sandstone and an artificial porous medium with spherical pores, were tested in the mechanochemical simulation. The positive stress-porosity feedback during pressure solution was observed in both samples. In addition, the model quantitatively illustrated the distribution of local mineral dissolution rates on all pore surfaces, as well as its relation to the effective mineral dissolution rate of the entire sample. Based on the comparison between the two porous media, the local mineral dissolution was regulated by pore space distribution, geometry, and coalescence during pressure solution. This work is the first that uses direct, pore-scale numerical simulation to demonstrate the positive stress-porosity feedback during pressure solution, which has the potential to advance the understanding of the mechanical-chemical (MC) coupling in many geological processes that are relevant to subsurface energy systems, such as the recovery of petroleum hydrocarbons and geothermal energy.

CT scanning↗
DOE OSTI2015

Pore-scale investigation on stress-dependent characteristics of granular packs and the impact of pore deformation on fluid distribution

Understanding the effect of changing stress conditions on multiphase flow in porous media is of fundamental importance for many subsurface activities including enhanced oil recovery, water drawdown from aquifers, soil confinement, and geologic carbon storage. Geomechanical properties of complex porous systems are dynamically linked to flow conditions, but their feedback relationship is often oversimplified due to the difficulty of representing pore-scale stress deformation and multiphase flow characteristics in high fidelity. In this work, we performed pore-scale experiments of single- and multiphase flow through bead packs at different confining pressure conditions to elucidate compaction-dependent characteristics of granular packs and their impact on fluid flow. A series of drainage and imbibition cycles were conducted on a water-wet, soda-lime glass bead pack under varying confining stress conditions. Simultaneously, X-ray micro-CT was used to visualize and quantify the degree of deformation and fluid distribution corresponding with each stress condition and injection cycle. Micro-CT images were segmented using a gradient-based method to identify fluids (e.g., oil and water), and solid phase redistribution throughout the different experimental stages. Changes in porosity, tortuosity, and specific surface area were quantified as a function of applied confining pressure. Results demonstrate varying degrees of sensitivity of these properties to confining pressure, which suggests that caution must be taken when considering scalability of these properties for practical modeling purposes. Changes in capillary number with confining pressure are attributed to the increase in pore velocity as a result of pore contraction. Furthermore, this increase in pore velocity was found to have a marginal impact on average phase trapping at different confining pressures.

58 GEOSCIENCES↗
DOE OSTI2026

Pore-Scale Transport Effects in Electrochemical CO 2 Reduction on Gold via Coupled Microkinetic-Transport Modeling

A pore-resolved modeling framework is developed to quantify how pore-scale transport affects the intrinsic microkinetics of CO 2 -to-CO on Au. A DFT-informed microkinetic model is coupled self-consistently to a Generalized-Modified Poisson–Nernst–Planck (GMPNP) transport description in a single, electrolyte-filled cylindrical pore, allowing local concentrations and potential to feed back into site-specific reaction rates. FIB-SEM is used to determine pore sizes within realistic electrode materials. Across pore diameters, d p = 10–6000 nm, the surface-averaged CO 2 reduction rate is systematically reduced relative to the ideal microkinetic baseline where mass transport is not accounted for; the effectiveness factor 𝜂 𝑠,CO 2 , which quantifies this ratio, decreases rapidly at more negative potentials and is about 1% near −1.0 V vs SHE due to reactant depletion. Spatial maps reveal pore-bulk alkalization that emerges at higher cathodic bias, with a small, near-wall pH dip due to electrostatic repulsion of hydroxide at the cathode interface. For a fixed aspect ratio L p /d p , narrower pores exhibit larger 𝜂 𝑠,CO 2 by shortening diffusion paths, whereas variations in the aspect ratio L p /d p play a secondary role. A dimensionless analysis (surface/bulk Damköhler numbers) delineates operating regimes. In conclusion, this work offers a concept for incorporating microkinetic models into homogenized porous-electrode models through effectiveness factors and pore-size distribution.

Au-catalyst↗
DOE OSTI2024

Accessing Fungal Contributions to the Birch Effect: Real-Time Respiration from Pore-Scale Microfluidics

Drying and rewetting of soil stimulates soil carbon emission. The Birch effect, driven by these cycles, leads to CO2 efflux, which can be monitored using real-time mass spectrometry (RTMS). Although soil fungi retain water during droughts, their contribution to CO2 release during drying–rewetting cycles remains unclear. In this study, we present the first demonstration of integrating micromodels with RTMS to monitor the Birch effect by simulating drought and rewetting. Micromodels were inoculated with axenic fungal culture and dried to assess moisture retention. After drying, RTMS quantified CO2 release upon rewetting with H218O mixtures. Our results showed that soil fungi released CO2 upon rehydration and immediately utilized the external water source at the pore scale by generating subsequent 46CO2. This work is the first to integrate RTMS with microsystems to investigate pore-scale biogeochemistry and the involvement of fungi in the Birch effect.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Characterizing the pore-scale impacts of enhanced-solubilization solutions on organic-liquid morphology and fluid-fluid interfacial area with microtomography

Synchrotron X-ray microtomography was used to directly quantify the impact of enhanced-solubilization solutions on mobilization of a non-wetting fluid in a natural porous medium. Polyoxyethylene sorbitan monooleate (tween 80), hydroxypropyl-β-cyclodextrin (HPCD), sodium dodecyl sulfate (SDS), and ethanol were used as the solubilization reagents. Tetrachloroethene (PCE) was used as the entrapped organic immiscible liquid. Microtomography images were collected prior to and after successive floods with each reagent solution. The total trapping number was varied by more than three orders of magnitude by changing the velocity and viscosity of the aqueous phase as well as the interfacial tension between fluid phases. The results were compared to those obtained from equivalent experiments conducted with water flooding. Inspection of images indicated that mobilization caused a change in PCE morphology and distribution for all reagents with the exception of HPCD, for which a lower trapping number was determined compared to the other experiments. Mobilization of organic liquid was observed at total trapping numbers that were consistent with prior work. The enhanced-solubilization solutions had minimal impact on the configuration or magnitude of fluid-fluid interfacial area.

Interfacial area↗
DOE OSTI2026

Enhanced Manganese Oxidation at the Biofilm–Fluid Interface Drives Pore-Scale Patterns in Mineral Precipitation

Microbial oxidation of manganese (Mn) from aqueous Mn(II) to solid-phase Mn(III, IV) minerals catalyzes Mn(II) removal in natural and engineered porous systems. However, little is known about the spatiotemporal evolution of Mn biomineralization in confined spaces that experience simultaneous Mn(II) delivery and Mn oxide precipitation. Here, we combine time-lapse microscopy, image analysis, and mass spectrometry to quantify the extent and rate of Mn biomineralization by Pseudomonas putida GB-1 in an optically transparent two-dimensional porous medium. We found that Mn(II) oxidation initially occurred within biofilms but shifted over time toward the edges of biofilms in contact with pore fluid. Minerals precipitated outside of the initial biofilm footprint likely due to surface-mediated oxidation of Mn(II) by nascent biogenic Mn oxides, reinforcing a gradient in mineral accumulation from the Mn(II) source near the reactor inlet to the outlet. The rate of mineral precipitation outside the biofilm footprint surpassed the rate of mineral accumulation inside biofilms within 6 h and accounted for two-thirds of the total Mn oxide mass in the pore space at the end of the experiment. This work advances a mechanistic understanding of coupled biotic and abiotic Mn oxidation in porous environments while providing a novel platform to quantify microbe-mineral-fluid interactions.

Biofilms↗
DOE OSTI2024

Pore-scale visualization of natural hydrate-bearing sediments

Accurate modeling of gas hydrate reservoir productivity and geomechanical risks associated with subsurface dissociation of natural gas hydrates (NGH) requires the determination of model parameters through physical testing on natural hydrate-bearing sediments (HBS). This involves investigating the hydro-mechanical behavior of undisturbed hydrate samples from nature under in situ conditions using pressure core characterization and analysis, which provides a unique opportunity for research. By employing state-of-the-art micro computed tomography imagery on cryogenically preserved, hydrate-bearing sediment samples, we can determine hydrate saturation as well as permeability with and without the presence of hydrates in the sediment. Furthermore, utilizing a machine learning based image segmentation technique, it is possible to extract pore space and grain information. Subsections of the entire image volume were used to determine anisotropic permeabilities using a finite-difference method Stokes solver (FDMSS). Additionally, permeability measurements on whole pressure and temperature preserved hydrate-bearing core were analyzed by utilizing the National Energy Technology Laboratory’s (NETL) Pressure Core Characterization and X-ray CT Visualization Tool (PCXT) to manipulate, cut, and analyze pressure preserved sediment. Permeabilities were measured under a broad range of vertical stress states to simulate expected pressure changes during production scenarios, and the results show that permeabilities derived from images are in agreement with those from traditional core derived experiments. The collected stress-dependent permeability, permeability anisotropy, and corresponding gas hydrate saturations provide valuable input into numerical simulations of reservoir productivity. These properties have been proven to be key parameters determining a long-term reservoir response under depressurization.

Liu, Mengwei [Oak Ridge Institute for Science and ↗
DOE OSTI2026

Iron as a driver of organic carbon fate in permafrost regions

Iron (Fe) and its biogeochemical interactions with organic carbon (OC) exert critical controls over carbon storage and water quality across Arctic and sub-Arctic landscapes. Permafrost thaw affects Fe–organic carbon interactions by exposing thawed material and driving shifts in hydrologic regimes. Here we propose a conceptual framework describing how Fe–OC interactions respond to permafrost thaw and the implications for global carbon cycling. The framework builds on current understanding of hydrologically driven processes that are altered by permafrost thaw, including microbial Fe reduction that influences OC decomposition and carbon emissions; Fe cycling through redox transformations and sulfide oxidation; dissolution and reprecipitation of Fe–OC assemblages; and Fe export to river networks. We put forward four testable hypotheses—on topics from Fe supply to its role in soil OC storage and transfer to rivers—to advance understanding of the coupled environmental controls on pore-scale reactions (that is, microscale reactions at mineral–water–microorganism interfaces) that are central to global Fe–OC cycling. Our framework highlights how soil-profile- and landscape-scale processes, such as hydrologic connectivity and redox dynamics, regulate pore-scale Fe–OC reactions.

Herndon, Elizabeth [ORNL] (ORCID:0000000291945493)↗
DOE OSTI2026

Geologic stress modulates fluid mixing at fracture intersections

Fracture intersections are critical links that enable flow and transport in subsurface fracture networks, and their behavior strongly influences fluid mixing in a network. Although all subsurface fractures are subjected to geological stress, we lack a fundamental understanding of how fracture intersection geometry evolves under stress and how these changes influence fluid mixing. Here, we combine 3D printing, 3D X-ray tomographic imaging, and 3D pore-scale numerical simulations to reveal stress-induced changes in intersection geometry and their impact on mixing. Mixing is found to be strongly affected by partial closure of an intersection under stress. As an intersection closes, the void area for fluid flow and diffusion decreases leading to substantial deviations between conventional mixing models and full pore-scale modeling. To address this, we propose a modified mixing model that accounts for intersection deformation, which is essential for accurate modeling of solute transport and mixing through fracture networks.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

How are Heterogeneous Nucleation Rate Observations Influenced by Instrument Resolution?

Experimental measurements of the heterogeneous nucleation rate rely on counting the number of nuclei with time. However, the size of a thermodynamically stable nucleus is often a few nanometers in diameter and is below the resolution of most (in situ) measurement techniques that provide a statistically valid sample. Due to the finite resolution of the instruments and analysis methods, it is challenging to capture the incipient nuclei and the subsequent evolution of nuclei density over time. In this work, we demonstrate the impact of instrument resolution on observed nuclei densities by comparing numerical modeling with experimental results. Further, to achieve this, we implemented heterogeneous nucleation within the pore-scale reactive transport modeling framework using classical nucleation theory (CNT). We compared the modeling results with nucleation rates measured using X-ray nanotomography (XnT) and evaluated how these impact the apparent values of the prefactor and interfacial energy based on CNT and the crystal growth rate. Specifically, we applied a resolution threshold (artificial resolution limit) in the model during nuclei counting to resemble an experimental resolution, ranging from 15 to 500 nm. The findings reveal that the instrument resolution significantly impacts the apparent prefactor and interfacial energy. Both apparent prefactor and interfacial energy decrease with a decrease in the instrument resolution. While deviation in the prefactor due to resolution is anticipated, those in the interfacial energy are unexpected. The approach described here allows one to correct apparent nucleation rates that depend on the instrument’s resolution to derive “intrinsic” CNT parameters for the prefactor and interfacial energy.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

The rigorous upscaling of advection-dominated transport in heterogeneous porous media via the Method of Finite Averages

Systems involving advection-dominated transport through heterogeneous porous and fractured media are ubiquitous in subsurface engineering applications. However, upscaling such systems continues to challenge rigorous modeling efforts, particularly when advection is stronger than diffusion at fine spatial scales (i.e., when the Péclet number is greater than one at length scales that characterize a system’s unit-cells, representative elementary volumes, or averaging regions). Here, in this work, we propose and validate a strategy for extending the Method of Finite Averages (MoFA), a rigorous upscaling methodology for heterogeneous porous media, to upscale transport systems experiencing stronger advection than diffusion at fine scales (i.e., fine-scale Péclet numbers greater than one). We detail the strategy, the physical conditions under which it can be applied while retaining a priori modeling error guarantees, and implement the strategy to obtain a MoFA model for advective-diffusive transport that accommodates advective physics at fine spatial scales. We then perform two numerical experiments considering systems with system-scale Péclet numbers of 300 and 1000 — which correspond to fine-scale Péclet numbers of 30 and 100, respectively — to verify that the error guarantees are met under the strategy. After, we conduct a numerical study to demonstrate the strategy’s advantages over the original MoFA methodology. The results suggest that rigorously-upscaled transport models for heterogeneous porous media experiencing advective physics at finer spatial scales can be derived through MoFA and resolved orders of magnitude faster than their pore-scale counterparts. The results also suggest that the presented strategy is limited to modeling shallow concentration gradients when there are large differences between the time scales related to advection and a system’s temporally-varying boundary conditions. This limitation hinders the strategy’s practicality in modeling more advective systems, and as such, opportunity exists for developing additional strategies that accommodate rapidly-varying boundary conditions — and consequentially, steeper concentration gradients — while modeling advective systems with MoFA.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A statistical and simulation-informed model for estimating permeability from pore size distribution in saturated geomaterials

Accurate permeability estimation is essential across subsurface engineering applications but remains challenging due to the complex pore structures of natural geomaterials. Traditional empirical methods and simplified theoretical models often inadequately capture the role of pore size distribution and connectivity. Here, this study develops a statistical and simulation-informed permeability model that collapses pore-scale complexity into a compact scaling of the form k = αϕμ d 2 , where ϕ is porosity, μ d is mean pore size, and α is a weakly varying coefficient. By combining pore network simulations with statistical analysis of unimodal and bimodal pore size distributions, we identify three key findings: (i) permeability is much more sensitive to mean pore size than to porosity; (ii) across extensive datasets, the ratio σ d /μ d (standard deviation to mean) clusters around a characteristic value ∼0.4, allowing the effects of the full pore size distribution to be represented by μ d and a narrowly varying α ≈ 0.05; and (iii) for bimodal systems, there exists a critical fraction of small pores ∼0.78 above which flow becomes small-pore dominated, enabling the definition of an effective flow-controlling pore population and facilitating simplified permeability estimation for such systems. The resulting model, which requires only porosity and a representative mean pore size as inputs, is validated against comprehensive experimental datasets (>1700 samples) spanning diverse soils and rocks and achieves good predictive accuracy. Overall, this work provides a physically grounded yet practically simple permeability estimator suitable for subsurface engineering, environmental protection, and resource management applications.

Permeability↗
DOE OSTI2024

PYSIMFRAC: A Python library for synthetic fracture generation and analysis

In this paper, we introduce PYSIMFRAC, an open-source python library for generating 3-D synthetic fracture realizations, integrating with fluid simulators, and performing analysis. PYSIMFRAC allows the user to specify one of three fracture generation techniques (Box, Gaussian, or Spectral) and perform statistical analysis including the autocorrelation, moments, and probability density functions of the fracture surfaces and aperture. This analysis and accessibility of a python library allows the user to create realistic fracture realizations and vary properties of interest. In addition, PYSIMFRAC includes integration examples to two different pore-scale simulators and the discrete fracture network simulator, dfnWorks. The capabilities developed in this work provides opportunity for quick and smooth adoption and implementation by the wider scientific community for accurate characterization of fluid transport in geologic media. We present PYSIMFRAC along with integration examples and discuss the ability to extend PYSIMFRAC from a single complex fracture to complex fracture networks.

58 GEOSCIENCES↗
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

Controlling matric potential in microfluidics to examine microbial dynamics in unsaturated porous media

The use of microfluidics for the study of soil microbial ecology is an emerging field. Most microfluidic studies of biological systems, however, have been performed under fully saturated conditions that are not representative of natural soil. Therefore, while microfluidics offer many unique capabilities that other methodologies cannot, they are not currently suited to address the effects of matric potential, an important variable defining the microbial moisture niche. Here, a methodology is presented that allows the user to control the aqueous conditions within microfluidic networks by manipulating matric potential using a hanging water column. The method relies on hydrophilic surface treatment of the microfluidic device using polyvinyl alcohol (PVA) and incorporating a bed of small pores at the network boundaries, which serve as a porous ceramic plate analogue (PPA). The method was validated on a simple capillary bundle and then on a more complex pore network. A water retention curve, exhibiting hysteresis, was generated for the pore network over a narrow matric potential range of 0 to – 5 kPa. Both the drainage and wetting curves were reproducible, as were the spatial configuration and the number of fragmented moisture niches in the pore network, particularly on the drainage curve. In contrast, the wetting curve exhibited greater variability in spatial configuration due to the “ink bottle effect,” where capillarity was interrupted by wider pore bodies. Ultimately, the methodology provides realistic pore-scale moisture conditions that can be easily manipulated and maintained, enabling new opportunities to explore soil biophysics and microbial biogeography in unsaturated porous media. As a brief example, images showing the localization of fluorescently tagged Pantoea sp. YR343 at −4.3 kPa are presented, highlighting bacterial distributions in water films and air-water interfaces.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Impact of Rock Wettability and Mineralization on CO 2 Storage Efficiency in Basalt Reservoirs

Basalt formations have emerged as highly promising targets for geological CO 2 storage due to their abundance of reactive silicate minerals that can rapidly convert dissolved CO 2 into stable carbonate minerals. The efficiency of CO 2 trapping in basalts, however, depends not only on their geochemical reactivity but also on reservoir wettability, which governs how CO 2 partitions among structural, residual, solubility, and mineral trapping mechanisms. However, wettability and mineral kinetics have largely been examined separately, wettability for plume migration and kinetics in terms of mineralization, leaving their coupled impact unquantified. Here, we systematically integrate th ese processes in a reservoir-scale reactive transport model to quantify the impact of three distinct wettability conditions (water-wet, mixed-wet, and CO 2 -wet) on plume migration, mineral precipitation, pore-scale blockage, and long-term immobilization. Results show that under ideal conditions (no mineralization), water-wet systems exhibited plume aspect ratios up to 28% higher than CO 2 -wet systems and promoted at least 60% greater residual trapping. When mineralization was incorporated, mineral trapping accounted for 57%, 53%, and 46% of the total immobilized CO 2 in water-wet, mixed-wet, and CO 2 -wet systems, respectively. Dissolved CO 2 contributed 43-46%, while residual trapping remained below 5%, and mobile CO 2 was substantially reduced due to mineralization, demonstrating that as the system becomes more water-wet, CO 2 is more effectively converted into stable carbonate minerals, although extensive mineral precipitation may lead to pore clogging. A global Sobol-Morris sensitivity analysis further reveals that the kinetic parameters of reactive minerals are the most influential factor, followed by wettability, while salinity has only a minor influence. Furthermore, these findings highlight that wettability is not merely a flow property but a fundamental control on CO 2 mineralization in basalts, with direct implications for site selection for storage projects.

Basalts↗
DOE OSTI2024

Quantitative local state of charge mapping by operando electrochemical fluorescence microscopy in porous electrodes

We introduce operando quantitative electrochemical fluorescence state of charge mapping (QEFSM), a non-invasive technique to study operating electrochemical systems along with a new design of optically transparent microfluidic redox flow cells compatible with the most demanding optical requirements. QEFSM allows quantitative mappings of the concentration of a particular oxidation state of a redox-active species within a porous electrode during its operation. In this study, we used confocal microscopy to map the fluorescence signal of the reduced form of 2,7-anthraquinone disulfonate (AQDS) in a set of multistep-chronoamperometry experiments. Calibrating these images and incorporating an analytical model of quinhydrone heterodimer formation with no free parameters, and accounting for the emission of each species involved, we determined the local molecular concentration and the state of charge (SOC) fields within a commercial porous electrode during operation. With this method, electrochemical conversion and species advection, reaction and diffusion can be monitored at heretofore unprecedented transverse and axial resolution (1 μm and 25 μm, respectively) at frame rates of 0.5 Hz, opening new routes to understanding local electrochemical processes in porous electrodes. Here, we observed pore-scale SOC inhomogeneities appearing when the fraction of electroactive species converted in a single pass through the electrode becomes large.

42 ENGINEERING↗
DOE OSTI2024

Multiscale Dynamics of Reactive Fronts in the Subsurface

Understanding and predicting flow and reactive transport in rocks (i.e. geologic porous media) is critical to many technologies at the heart of the energy transition, including CO 2 sequestration and H2 storage. However, accurate modeling and prediction of these systems is very complex because physico-chemical processes that occur at very small spatial scales, i.e. in the pores of the rocks, can dramatically control the system performance at the field scale (km). For example, precipitation reactions at the pore-scale can lead to large permeability changes at the field scale and dramatically alter the migration of stored gases. Properly accounting for multi-scale coupling effects is critical to achieve predictivity and confidence in model outputs, which then can guide design and optimization at the system-scale. This can be achieved through the development of rigorous mathematical models that can appropriately account for fine-scale effects at the large scale. The final report of the Early Career award DE-SC0019075 “Multiscale dynamics of reactive fronts in the subsurface” summarizes the mathematical, numerical and experimental advancements to study and predict reactive transport in geologic porous media across scales.

58 GEOSCIENCES↗