Search NASA⌕ Search

SEARCH · Search NASA

Results for “Heterogeneous media”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Finite Element Modeling of Diffusion in Fractured Porous Media by Using Hierarchical Material Properties

Abstract Fractured porous media challenge modeling approaches due to high computational costs and excessive mesh refinement imposed by the extreme scale variability of fractures and the heterogeneity of the surrounding porous rock. To overcome such difficulties, we utilize the hierarchical finite element method ( Hi ‐FEM) that has been developed previously to simulate the electrical potential distribution in complex geologic environments. The method employs the hierarchical basis functions in classical finite element analysis to enable representation of material properties on each dimensional component of a given 3D unstructured finite element, thereby inherently allowing for interactions at the boundary between fracture and a host rock. In this study, we extend its application to transient fluid flow and heat conduction in the Laplace domain. Time‐domain flow solutions are obtained by numerical inverse Laplace transform. We evaluate the accuracy of the method using different flow models and demonstrate its robustness for large‐scale, rock mass models featuring complex fracture networks. Moreover, for the computation of nodal Darcian velocity fields in fractured porous media where the fractures are represented as 2D features, a new approach that employs the Yeh's Galerkin model for both volume and facet elements is proposed. Results show that Hi ‐FEM can produce accurate flow solutions for fractured porous media without any need of coupling or transfer mechanism while still being computationally economical and numerically robust, even for large‐scale simulations.

Beskardes, G. D.↗
DOE OSTI2025

Permeability Prediction Using Vision Transformers

Accurate permeability predictions remain pivotal for understanding fluid flow in porous media, influencing crucial operations across petroleum engineering, hydrogeology, and related fields. Traditional approaches, while robust, often grapple with the inherent heterogeneity of reservoir rocks. With the advent of deep learning, convolutional neural networks (CNNs) have emerged as potent tools in image-based permeability estimation, capitalizing on micro-CT scans and digital rock imagery. This paper introduces a novel paradigm, employing vision transformers (ViTs)—a recent advancement in computer vision—for this crucial task. ViTs, which segment images into fixed-sized patches and process them through transformer architectures, present a promising alternative to CNNs. We present a methodology for implementing ViTs for permeability prediction, its results on diverse rock samples, and a comparison against conventional CNNs. The prediction results suggest that, with adequate training data, ViTs can match or surpass the predictive accuracy of CNNs, especially in rocks exhibiting significant heterogeneity. This study underscores the potential of ViTs as an innovative tool in permeability prediction, paving the way for further research and integration into mainstream reservoir characterization workflows.

58 GEOSCIENCES↗
DOE OSTI2026

Spatiotemporal and Statistical Mapping of Transition Metal Equilibria in Alkaline Media

Transition metal dissolution and redeposition (D/R) kinetics in alkaline media play a critical role in various chemical and electrochemical processes. Competitive reaction kinetics between different transition metals can modulate individual metal behavior in these processes. To date, these phenomena have remained largely unmeasured, and even when captured, they are difficult to statistically characterize due to their dynamic nature, simultaneous occurrence, and spatially heterogeneous nature. Here, in this study, we develop a statistical analysis framework based on in situ and operando X-ray fluorescence microscopy (XFM) to investigate the relative D/R kinetics of multiple transition metals in alkaline media. By employing statistical analysis, we quantify the spatial distribution of D/R species and assess the rate at which the system reaches equilibrium under varying reaction conditions. We show that pH does not simply change the rate of dissolution and redeposition, but reorganizes the cross-element kinetic correlations among Ni, Fe, and Mn and accelerates the spatial equilibration of D/R events, as quantified through correlation analysis, reaction-rate estimation, probability function distributions, and texture-based monitoring statistics. Additionally, we demonstrate how modifying the solvent environment can influence D/R kinetics, providing a pathway for tuning materials synthesis and process optimization. Our study offers valuable insights into the complex interplay between different transition metals and provides a reliable statistical framework for spatial analysis of diverse imaging data sets, enabling deeper extraction of latent information across multiple modalities.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Hidden order across online extremist movements can be disrupted by nudging collective chemistry

Disrupting the emergence and evolution of potentially violent online extremist movements is a crucial challenge. Extremism research has analyzed such movements in detail, focusing on individual- and movement-level characteristics. But are there system-level commonalities in the ways these movements emerge and grow? Here we compare the growth of the Boogaloos, a new and increasingly prominent U.S. extremist movement, to the growth of online support for ISIS, a militant, terrorist organization based in the Middle East that follows a radical version of Islam. We show that the early dynamics of these two online movements follow the same mathematical order despite their stark ideological, geographical, and cultural differences. The evolution of both movements, across scales, follows a single shockwave equation that accounts for heterogeneity in online interactions. These scientific properties suggest specific policies to address online extremism and radicalization. We show how actions by social media platforms could disrupt the onset and ‘flatten the curve’ of such online extremism by nudging its collective chemistry. Our results provide a system-level understanding of the emergence of extremist movements that yields fresh insight into their evolution and possible interventions to limit their growth.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

Earth Sciences Applications for Energy and Global Security: A Few Vignettes [Slides]

Understanding the physics of flow and transport in the subsurface environment is crucial to several National Security applications. In the Energy Security realm, some examples include characterizing radioactive waste disposal, exploring natural gas extraction and sequestering CO2. An example in the Global Security arena is to be able to distinguish between anthropogenic and naturally occurring seismic disturbances, the former being of particular concern as underground nuclear testing activity has seen an increase in recent years. A common theme across these applications is the ability to accurately model fluid flow in porous and/or fractured subsurface media. For decades, researchers have grappled with how to include the effects of uncertainties in these models, including but not limited to heterogeneity, unknown initial and boundary conditions and issues with scalability. More recently advances in computing and machine learning methodologies have allowed a broader exploration of the uncertainty space as we strive for real time decision making. I will present an overview of the field and discuss a few examples of how my own research in these topics has evolved over nearly two decades.

58 GEOSCIENCES↗
DOE OSTI2020

Grain detachment and transport clogging during mineral dissolution in carbonate rocks with permeable grain boundaries

In this paper, we explore the effect of permeable grain boundaries on carbonate rock dissolution in order to improve our understanding of grain detachment and migration during reactive flow. To do so we investigated the effects of grain size and enhanced reactivity along grain boundaries on global dissolution kinetics. Variations of permeability and porosity were calculated, and exponential relationships were observed in fractured rocks. Our model employed a reactive transport framework based on the Darcy-Brinkman equation to simulate calcite dissolution in carbonate rocks composed of microporous grains and, included fluid transport along grain boundaries. The model includes fluid flow, solute transport by advection–diffusion, heterogeneous reaction between fluid and minerals and grain detachment with subsequent grain transport in macropores. The migration of solid particles due to dissolution was based on cluster analysis and local movement, and the results show that grain detachment can lead to significant decreases in permeability due to the clogging of transport pathways. Microporous media with smaller grain sizes (fine grains in this study) showed a higher average reaction rate than those with coarser grains. In addition, as the overall rates of geochemical processes are commonly affected by the presence of texture heterogeneities such as fractures, a single fracture (macropore) introduced into a microporous matrix composed of permeable grains and grain boundaries was modelled as a large channel connecting the inlet and outlet of the simulation domain. It was found that macropore clogging by grain detachment temporary decreased permeability, lowering the long-term global reaction rates. Lastly, we observed that increases in flow rate can reduce detachment and fracture clogging by reducing local dissolution along grain boundaries.

58 GEOSCIENCES↗
DOE OSTI2022

Dynamic and single cell characterization of a CRISPR-interference toolset in Pseudomonas putida KT2440 for β-ketoadipate production from p -coumarate

We report Pseudomonas putida KT2440 is a well-studied bacterium for the conversion of lignin-derived aromatic compounds to bioproducts. The development of advanced genetic tools in P. putida has reduced the turnaround time for hypothesis testing and enabled the construction of strains capable of producing various products of interest. Here, we evaluate an inducible CRISPR-interference (CRISPRi) toolset on fluorescent, essential, and metabolic targets. Nuclease-deficient Cas9 (dCas9) expressed with the arabinose (8K)-inducible promoter was shown to be tightly regulated across various media conditions and when targeting essential genes. In addition to bulk growth data, single cell time lapse microscopy was conducted, which revealed intrinsic heterogeneity in knockdown rate within an isoclonal population. The dynamics of knockdown were studied across genomic targets in exponentially-growing cells, revealing a universal 1.75 ± 0.38 hour quiescent phase after induction where 1.5 ± 0.35 doublings occur before a phenotypic response is observed. To demonstrate application of this CRISPRi toolset, β-ketoadipate, a monomer for performance-advantaged nylon, was produced at a 4.39 ± 0.5 g/L and yield of 0.76 ± 0.10 mol/mol from p-coumarate, a hydroxycinnamic acid that can be derived from grasses. These cultivation metrics were achieved by using the higher strength IPTG (1K)-inducible promoter to knockdown the pcaIJ operon in the βKA pathway during early exponential phase. This allowed the majority of the carbon to be shunted into the desired product while eliminating the need for a supplemental carbon and energy source to support growth and maintenance.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Scalable simulation of coupled adsorption and transport of methane in confined complex porous media with density preconditioning

The growing significance of shales and tight formations in the transition to less carbon-intensive and clean energy drives the research endeavor to understand the physics of gas flow within these systems. However, shales are composed of massively heterogeneous physical and chemical features. Most nano-sized pores connect to millimeter-scale fractures, leading to multiscale transport. These nano-scale pore throats demonstrate non-classical flow behavior, such as non-negligible slip velocities and adsorbed gas layers at the boundary. As a result, classical computational fluid dynamics models do not capture the physics. In this work, we develop a coupling scheme for the multiple-relaxation-time (MRT) lattice Boltzmann (LB) method that integrates the Peng-Robinson equation of state into a pseudo-potential interaction model to capture the physics of methane flow in irregular networks of channels that represent nano-scale porous media. We use atomistic simulations to calibrate and validate our model in slit nano-channels. We propose a preconditioning scheme to initialize the coupled transport and adsorption simulation of methane in complex porous media. The results of this implementation of LB agree with Direct Simulation Monte Carlo (DSMC) and Molecular Dynamics (MD) simulations. We then scale up the LB implementation through vectorization and indirect addressing. We parallelize it using Message Passing Interface (MPI) and OpenMP frameworks to simulate transport and adsorption in complex media with a million lattices. Additionally, we analyze the differences between coupled and transport-only simulations in two case studies and show that considering phase behavior, i.e., adsorption, can significantly change the flow behavior. This work constitutes an important step towards bridging the gap between molecular flow and system-scale behavior of complex disordered porous media.

42 ENGINEERING↗
DOE OSTI2021

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗
DOE OSTI2020

Three-Dimensional Imaging and Quantification of Gas Storativity in Nanoporous Media via X-rays Computed Tomography

This study provides the engineering science underpinnings for improved characterization and quantification of the interplay of gases with kerogen and minerals in shale. Natural nanoporous media such as shale (i.e., mudstone) often present with low permeability and dual porosity, making them difficult to characterize given the complex structural and chemical features across multiple scales. These structures give nanoporous solids a large surface area for gas to sorb. In oil and gas applications, full understanding of these media and their sorption characteristics are critical for evaluating gas reserves, flow, and storage for enhanced recovery and CO 2 sequestration potential. Other applications include CO 2 capture from industrial plants, hydrogen storage on sorbent surfaces, and heterogeneous catalysis in ammonia synthesis. Therefore, high-resolution experimental procedures are demanded to better understand the gas–solid behavior. In this study, CT imaging was applied on the sub-millimeter scale to shale samples (Eagle Ford and Wolfcamp) to improve quantitative agreement between CT-derived and pulse decay (mass balance) derived results. Improved CT imaging formulations are presented that better match mass balance results, highlighting the significance of gas sorption in complex nanoporous media. The proposed CT routine implemented on the Eagle Ford sample demonstrated a 17% error reduction (22% to 5%) when compared to the conventional CT procedure. These observations are consistent in the Wolfcamp sample, emphasizing the reliability of this technique for broader implementation of digital adsorption studies in nanoporous geomaterials.

42 ENGINEERING↗
DOE OSTI2024

Quantifying Dissolution Dynamics in Porous Media Using a Spatial Flow Focusing Profile

Abstract The diverse range of patterns in porous media formed by dissolution processes depends on the relative magnitude of flow, transport, and chemical reactions at pore surfaces. However, distinguishing between regimes often relies solely on qualitative, visual comparisons of emergent structures. Here, we propose a quantitative measure capable of identifying different regimes using the concept of the spatial flow focusing profile, which segments the medium into cross sections along the flow direction to calculate the flow focusing index for each section. We employ this measure in numerical simulations of a dissolving porous medium using a pore network model. We obtain a morphological phase diagram of dissolution patterns, which we characterize using the flow focusing profile. In particular, we demonstrate that analyzing the temporal changes in the profile allows one to quantitatively distinguish between wormholing and channeling. The transition between them is shown to be affected by the heterogeneity of the system.

58 GEOSCIENCES↗
DOE OSTI2024

EFRC-MUSE: Multi-Scale Fluid-Solid Interactions in Architected and Natural Materials

Phase interactions and fluid properties in geological and other environments are critical in applications ranging from hydrogen production and geologic storage and recovery, carbon dioxide storage and sequestration, and the sustainable use of water resources. The four goals of EFRC-MUSE: Multi-Scale Fluid-Solid Interactions in Architected and Natural Materials were based on the priority directions articulated in the Basic Research Needs documents, and the scientific needs in nanoscience. 1. Develop a fundamental understanding of confinement and surface interactions in mesoscale media with nanometer-sized pores on the phase behavior, thermodynamic and multiphase flow properties of multicomponent fluid mixtures. 2. Examine the impact of mineralogy and material heterogeneity on mechanical properties to better understand chemo-mechanical interactions in material failure. 3. Determine in-operando cross-scale structural and nanostructural material properties with fluids in confinement and flow under realistic condi

58 GEOSCIENCES↗
DOE OSTI2025

Hydrodynamic characterization of the redox chemistry of crown-encapsulated uranyl complexes

The redox properties of actinide-containing species strongly influence their reactivity, speciation, and interfacial behavior, but the experimental quantification of the electrochemical characteristics of molecular actinide complexes in nonaqueous media has not received the attention it deserves. Here, results from hydrodynamic methods and electrochemical simulations of U(VI)/U(V) redox are reported, including quantification of heterogeneous electron-transfer kinetics and estimation of chemical reversibility of U(VI)/U(V) interconversion at electrodes in acetonitrile-based electrolyte. The complexes investigated are recently reported U(VI) and U(V) complexes in which the uranyl ion (UO 2 n+ ) is encapsulated in a macrocyclic 18-crown-6-like moiety templated by a Pt(II) center. These complexes feature the most positive value U VI /U V reduction potential yet reported and are thus particularly relevant to study of facile U(V) generation from U(VI) precursors as well as uranium electroanalysis. Rotating disk electrode (RDE) studies have been used to quantify the diffusion coefficients of the U(VI) and U(V) complexes, and standard heterogeneous electron transfer rate constants ( k 0 ) for the redox have been determined using a conventional Koutecký-Levich analysis. Rotating ring-disk electrode (RRDE) studies have been used to directly interrogate the chemical reversibility of U(VI)-U(V) interconversion, confirming that reduction of the U(VI) complex at an Au disk is associated with formation of the U(V) analogue that can be readily re-oxidized at a Pt ring under hydrodynamic (rotating) conditions. Because measurements of the type reported here are generally associated with current flows that are larger than those found in corresponding quiescent (unstirred) conditions, our findings suggest that hydrodynamic methods could be advantageous for design of electroanalytical approaches to detection of actinide species and study of their redox properties.

actinides↗
DOE OSTI2021

Computationally Efficient Multiscale Neural Networks Applied to Fluid Flow in Complex 3D Porous Media

Abstract The permeability of complex porous materials is of interest to many engineering disciplines. This quantity can be obtained via direct flow simulation, which provides the most accurate results, but is very computationally expensive. In particular, the simulation convergence time scales poorly as the simulation domains become less porous or more heterogeneous. Semi-analytical models that rely on averaged structural properties (i.e., porosity and tortuosity) have been proposed, but these features only partly summarize the domain, resulting in limited applicability. On the other hand, data-driven machine learning approaches have shown great promise for building more general models by virtue of accounting for the spatial arrangement of the domains’ solid boundaries. However, prior approaches building on the convolutional neural network (ConvNet) literature concerning 2D image recognition problems do not scale well to the large 3D domains required to obtain a representative elementary volume (REV). As such, most prior work focused on homogeneous samples, where a small REV entails that the global nature of fluid flow could be mostly neglected, and accordingly, the memory bottleneck of addressing 3D domains with ConvNets was side-stepped. Therefore, important geometries such as fractures and vuggy domains could not be modeled properly. In this work, we address this limitation with a general multiscale deep learning model that is able to learn from porous media simulation data. By using a coupled set of neural networks that view the domain on different scales, we enable the evaluation of large ( $$>512^3$$ > 512 3 ) images in approximately one second on a single graphics processing unit. This model architecture opens up the possibility of modeling domain sizes that would not be feasible using traditional direct simulation tools on a desktop computer. We validate our method with a laminar fluid flow case using vuggy samples and fractures. As a result of viewing the entire domain at once, our model is able to perform accurate prediction on domains exhibiting a large degree of heterogeneity. We expect the methodology to be applicable to many other transport problems where complex geometries play a central role.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Transport Upscaling under Flow Heterogeneity and Matrix-Diffusion in Three-Dimensional Discrete Fracture Networks

For this work, we investigate the combined effects of network scale flow variability and retention due to matrix-diffusion on the scaling behavior of transport through fractured media. Two of the principal mechanisms controlling the transport of solutes through fractured low-permeability media are broad distributions of flow velocities and retention times in the solid matrix. We study the relative impact of these two processes under different initial conditions using a set of three-dimensional discrete fracture network simulations. We use these simulations to develop and calibrate an upscaled continuous time random walk (CTRW) approach for advective transport based on an Ornstein-Uhlenbeck model for the particle velocities that accounts for the fracture-matrix coupling using a compound Poisson process. This CTRW model can be conditioned on the initial solute distribution and allows to observe late-time scaling behavior at distances beyond what is feasible to observe using high-fidelity direct numerical simulations. We determine that the initial distribution of particles leads to marked differences in the persistent long-term scale behavior in the solute travel time distributions, even those undergoing retention due to matrix diffusion through implementation and analysis of the model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Spectroscopic Insight on Neodymium Solvation in Lithium Borohydride-Supported Electrolyte

Borohydride-based electrolytes have recently emerged as promising media for the electrodeposition of electropositive metals, including rare earth (RE) elements. While the presence of supporting alkali metal cations and RE counteranions provides essential electrochemical conductivity for achieving fast metal electrodeposition, interactions between the host ligand and solvated neodymium (Nd) complexes remain unclear. This study provides insights into the coordination structure of a concentrated and directly solvated Nd salt in a lithium borohydride-supported electrolyte. Our spectroscopic results indicate that the RE coordination environment is significantly influenced by the solvation mechanism, which can vary between metathesis and complexation pathways, primarily dictated by stoichiometric factors. Under dilute conditions, nearly complete metathesis of anions leads to a high coordination number for the host ligand (borohydride), consistent with the previously reported solvated Nd speciation in chlorine-free electrolytes. In contrast, concentrated dissolution of the Nd salt in the supported electrolyte is dominated by a complexation pathway featuring a Li-ion-paired complex with a low coordination number of the host ligand. Density functional theory (DFT) calculations indicated that the observed blue shift in the borohydride vibration was the result of an increase in electron density drawn into the terminal B–H interbond region from the hydride as the coordination changed from Li to Nd. In conjunction with DFT results, vibrational analyses allowed correlation of the experimental shifts associated with changes in Nd ligation and coordination spheres, further consolidating the prevalence of highly chloride-coordinated species under concentrated conditions. In conclusion, the outcomes of this work illuminate the distinctive and heterogeneous coordination structures that the electroactive RE species can adopt at high concentrations in lithium borohydride-supported electrolytes, as a key step to comprehend the reported metal electrodeposition performance in these media.

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

A data-driven framework for permeability prediction of natural porous rocks via microstructural characterization and pore-scale simulation

Understanding the microstructure–property relationships of porous media is of great practical significance, based on which macroscopic physical properties can be directly derived from measurable microstructural informatics. However, establishing reliable microstructure–property mappings in an explicit manner is difficult, due to the intricacy, stochasticity, and heterogeneity of porous microstructures. In this paper, a data-driven computational framework is presented to investigate the inherent microstructure–permeability linkage for natural porous rocks, where multiple techniques are integrated together, including microscopy imaging, stochastic reconstruction, microstructural characterization, pore-scale simulation, feature selection, and data-driven modeling. A large number of 3D digital rocks with a wide porosity range are acquired from microscopy imaging and stochastic reconstruction techniques. A broad variety of morphological descriptors are used to quantitatively characterize pore microstructures from different perspectives, and they compose the raw feature pool for feature selection. Here high-fidelity lattice Boltzmann simulations are conducted to resolve fluid flow passing through porous media, from which reliable permeability references are obtained. The optimal feature set that best represents permeability is identified through a performance-oriented feature selection process, upon which a cost-effective surrogate model is rapidly fitted to approximate the microstructure-permeability mapping via data-driven modeling. This surrogate model exhibits great advantages over empirical/analytical formulas in terms of prediction accuracy and generalization capacity, which can predict reliable permeability values spanning four orders of magnitude. Besides, feature selection also greatly enhances the interpretability of the data-driven prediction model, from which new insights into the mechanism of how microstructural characteristics determine intrinsic permeability are obtained.

58 GEOSCIENCES↗
DOE OSTI2025

Pore Resolved Simulations of Joule Heating in Fibrous Media using an Embedded Boundary Method

Joule heating has been regarded as an energy-efficient and sustainable method for heating materials and gases at large scales. The modeling of local temperature effects at pore-resolved scales for such systems, however, has been difficult to achieve due to challenges in coupling thermo-chemical processes in complex porous media and in large representative volume elements (RVEs). To this end, we developed an electro-thermal model at the pore scale to study Joule heating effects in large heterogeneous systems with different microstructures. This was achieved using the level set method to implicitly delineate distinct regions within the domain, and an embedded boundary method to facilitate heat exchange across the fluid-solid interface. Moreover, we applied this method to investigate unsteady non-linear electro-thermal effects in non-woven fibrous graphite conductors for RVEs with characteristic lengths of 2 mm, with different fiber orientations, porosity (80% – 90%) and fiber diameters (10 – 20µm). The coupled equations were solved numerically and they produced peak temperatures greater than 2000 K resulting in heating rates as high as 80,000 K/s. Moreover, the results depended strongly on the microstructure of the fiber skeleton and current density. Geometries with large fibers (∼ 20µm) had the highest average and peak temperatures with the mean temperature increasing by 3.9 % while the peak temperature increased by 9.9 %. Anisotropic domains on the other hand had the lowest mean and peak temperatures with peak and mean temperatures of 2293 K and 1437.7K respectively representing a corresponding 12.1% and 5.1% drop in the temperatures. An increase in porosity from 80% to 90%, however, led to an increase in the peak temperature by 5.1%.

Joule heating↗