Pore-scale reactive transport and flow behavior involving dissolution and precipitation in heterogeneous porous media
Explore the source record for details and available documents.
Yoon, Hongkyu↗
SEARCH · Engineering Papers
Results for “Porous flow”
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.
Microcapsule Transport in Porous Media for Treatment of Non-Ideal Flow Zones
Explore the source record for details and available documents.
Kibikas, William Martin [Sandia National Laborator↗
Microcapsule Transport in Porous Media for Treatment of Non-Ideal Flow Zones
Explore the source record for details and available documents.
Kibikas, William Martin [Sandia National Laborator↗
A Dynamic Pore Network Model for Imbibition Simulation Considering Corner Film Flow
Wetting films can develop in the corners of angular pores under strong wetting conditions. Modeling the dynamics of corner film remains elusive using direct numerical simulations because of the significant scale difference between main meniscus and corner film flow. In this paper, the modified interacting capillary bundle model (ICB), developed in our previous work to describe accurately corner film dynamics in a single square tube, is incorporated into a single-pressure dynamic pore network model (DPNM) to simulate imbibition in strongly wetting porous media with corner film flow. The traditional pore network is decomposed into several layers of interacting subpore networks where the 0th layer of subpore network simulates the main meniscus flow and higher layers the corner film flow. The fluid flow between different layers is captured by interlayer throats. In addition, the snap-off mechanism caused by the thickening of wetting corner film is considered. The accuracy of the developed model is validated for four cases: spontaneous imbibition in a single square tube, wetting fluid redistribution through corner films under a capillary pressure difference, snap off in a narrow throat connecting two large pores, and imbibition dynamics in a real microfluidic porous geometry. The validated model is then used to simulate both spontaneous and controlled imbibition in a pore network with random pore size distribution. Finally, the interaction between corner film and main meniscus flow in porous media is analyzed from a pore-scale perspective.
58 GEOSCIENCES↗
3D Printed Nickel–Molybdenum-Based Electrocatalysts for Hydrogen Evolution at Low Overpotentials in a Flow-Through Configuration
Three-dimensional (3D) printed, hierarchically porous nickel molybdenum (NiMo) electrocatalysts were synthesized and evaluated in a flow-through configuration for the hydrogen evolution reaction (HER) in 1.0 M KOH(aq) in a simple electrochemical H-cell. 3D NiMo electrodes possess hierarchically porous structures because of the resol-based aerogel precursor, which generates superporous carbon aerogel as a catalyst support. Relative to a traditional planar electrode configuration, the flow-through configuration allowed efficient removal of the hydrogen bubbles from the catalyst surface, especially at high operating current densities, and significantly decreased the overpotentials required for HER. An analytical model that accounted for the electrokinetics of HER as well as the mass transport with or without the flow-through configuration was developed to quantitatively evaluate voltage losses associated with kinetic overpotentials and ohmic resistance due to bubble formation in the porous electrodes. Furthermore, the chemical composition, electrochemical surface area (ECSA), and roughness factor (RF) were also systematically studied to assess the electrocatalytic performance of the 3D printed, hierarchically porous NiMo electrodes. An ECSA of 25163 cm2 was obtained with the highly porous structures, and an average overpotential of 45 mV at 10 mA cm–2 was achieved over 24 h by using the flow-through configuration. The flow-through configuration evaluated in the simple H-cell achieved high electrochemical accessible surface areas for electrochemical reactions and provided useful information for adaption of the porous electrodes in flow cells.
08 HYDROGEN↗
Enforcing global constraints for the dispersion closure problem: τ 2 -SIMPLE algorithm
Permeability and effective dispersion tensors are critical parameters to characterize flow and transport in porous media at the continuum scale. Homogenization theory defines a framework in which such effective properties are first computed from solving a closure problem in a repeating unit cell of the periodic microstructure and then used in a macroscopic formulation for efficient computation. The closure problem is formulated as a local boundary value problem subjected to global constraints, which guarantee the uniqueness of the solution and can be difficult to satisfy for complex geometries and at high flow conditions. These constraints also ensure that pore-scale pressure, velocity, and concentration fields can be accurately reconstructed from the closure variable. Building on a previous work, here we present a framework that allows to satisfy global constraints associated to both the permeability and the dispersion closure problems by introducing two artificial time scales. The algorithm, called τ 2 -SIMPLE, computes both permeability and effective dispersion given an arbitrarily complex geometry and flow condition. Furthermore, this algorithm is demonstrated to be accurate for both 2D and 3D geometries across varying flow conditions, and thus it can be used to quickly characterize effective properties from porous media images in many applications.
97 MATHEMATICS AND COMPUTING↗
Non-destructive evaluation of specific surface area of porous powder systems using experimental and virtual air permeametry
The correlation of flow-permeable surface area (FSSA) of porous initiator powder to its performance is well-established. More specifically, with age, such powders tend to coarsen (i.e., lose surface area), with an accompanying deterioration in performance. Commercially available permeametry apparatuses, such as that from Fisher Scientific or a more modern version from Micromeritics suffers from a few drawbacks, including: (1) the need to work with a fixed mass equivalent of 1 cm 3 of solid material; (2) the need to re-compact the material to monitor aging effects (which can introduce unwanted errors in the derived results); (3) a superfluous design that involves a reference plug, and for the older Fisher Scientific model, the need to eyeball water meniscus height that can lead to inaccuracies; (4) calibrated and tested only for samples packed in cylindrical geometries.
36 MATERIALS SCIENCE↗
Improving the modeling of near-wall interphase heat transfer in porous media models of Pebble Bed Reactors
Here, this work aims to improve capabilities for modeling localized effects in porous media models of Pebble Bed Reactors. The wall-channeling effect is the primary local phenomenon of interest in a PBR, where the presence of the reflector wall disrupts the pebble packing, causing the pebbles near the wall to pack less efficiently and creating large void regions. Accurate modeling of the near-wall region is important as it will affect core bypass flow and temperature predictions. Porous media models are commonly used for design scoping and plant-level simulations of PBRs. Although these models have some capabilities to model the near-wall region, the correlations that are available in porous media codes are often inaccurate when a multi-region model is used to discretize the near-wall region. This work employs a high-to-low analysis to study the accuracy of available interphase heat transfer closures. NekRS, a spectral element computational fluid dynamics code, is used to perform Large Eddy Simulations. These LES simulation results are compared to porous media model results from the Pronghorn porous media code. The friction term of the KTA drag closure is first improved, reducing the error in the prediction of the near-wall velocity from over 50% to less than 5%. This is combined with improvements to the form term from previous works to produce a drag closure that is capable of accurately modeling the wall-channeling effect across a variety of flow conditions. The Nusselt number predictions of several heat transfer correlations are compared to the high-fidelity results where it is found that the KTA heat transfer correlation is capable of accurately predicting the local Nusselt numbers that were determined in the high-fidelity simulation. Comparison of the radial solid temperature profiles, however, reveal discrepancies between NekRS and Pronghorn. It is discovered that the implementation of the interphase heat transfer coefficient that exists in many current porous media codes is not valid when local porosities are modeled. Instead, it is suggested that the interphase heat transfer coefficient should be dependent on the local porosity, the Nusselt number, and the local solid surface-to-volume ratio. Implementation of this change produces improvement in the agreement between the results obtained by NekRS and Pronghorn while using the KTA heat transfer correlation.
interphase heat transfer↗
Colloidal Three-Dimensional Covalent Organic Frameworks and Their Application as Porous Liquids
Developing solid porous materials into free-flowing liquids with permanent porosity is a promising strategy for overcoming certain limitations of conventional sorbent materials employed in gas storage and separation applications. The ability to control the pore size and chemical functionalities of organic frameworks gives these particular nanoporous materials distinct advantages over other small cage-like molecules or hollow particles that have been developed into porous liquids. Here, we describe the synthesis of a 3D imine-linked colloidal covalent organic framework (COF)-based porous liquid, designed for efficient size-exclusion of solvent, as well as for long-term stability. By tethering ionic liquids to the colloid surface, the colloids can be dried, purified, and resuspended in a variety of solvents without irreversible aggregation typically observed of COF colloids. Colloid size could be controlled between 50 and 400 nm, with surface areas as high as 800 m2/g. The 3D intertwining morphology of the colloids had pore sizes ranging from 5 to 14 Å, allowing them to efficiently sizeexclude bulky ionic liquids. The COF colloids were stable towards flocculation in an ionic liquid for > 1 year. Permanent porosity was confirmed with a combination of 19F NMR measurements and gas sorption techniques. CO2 and CH4 uptake in these porous liquids increased more than 10 and 20-fold, respectively, over non-porous, neat ionic liquid control samples. The work not only advances the state of COF-based colloid science but also represents a practical advance towards developing more robust, tunable framework-based porous liquid materials for a host of gas storage and separation applications.
Mow, Rachel E.↗
Fast and Accurate Core Analysis by the Full-Immersion Pressure-Pulse Decay: Part 1—Theory
Summary Core-scale measurements are considered the ground truth that oil and gas or electric utility companies use to predict the migration of fluids such as oil, natural gas, carbon dioxide, or brine deep underground during their extraction or injection operations. To provide a greater understanding of petrophysical properties of low-permeability geologic formations such as shales and tight gas sandstones, this study introduces a novel core-analysis procedure. The technique follows conventional pressure-pulse-decay permeametry, where the pressure in an inlet chamber adjacent to a cylindrical core plug undergoes a rapid pressurization, the system is shut in, and the pressure reaches a new equilibrium. However, unlike a standard unidirectional pressure-pulse decay, the full-immersion pressure-pulse decay (Hannon 2019) applies a pressure disturbance to the entire outer surface area of the sample. This article covers the numerical simulator designed to model flow through an anisotropic porous sample in this scenario. The model assumes distinct but uniform permeabilities along the radial and axial directions of the cylindrical plug sample. When extracting a plug vertically (or perpendicular to bedding), the permeability along the radial direction associates with the horizontal permeability (i.e., parallel to bedding), whereas in the axial direction, flow occurs perpendicular to bedding (a vertical permeability). Investigations of these model outputs demonstrate an approximately 20-fold decrease in time to complete a full-immersion experiment compared with conventional pressure-pulse decay. Furthermore, the pressure-decay curves resulting from the full-immersion method have slightly different shapes than those resulting from other unidimensional transient methods. These differences begin to demonstrate that, under achievable experimental conditions, the analysis of pressure data from one full-immersion test could enable the simultaneous estimation of the apparent permeabilities parallel and perpendicular to bedding of a cylindrical sample in addition to its porosity. A follow-up article finalizes the proof of this capability with a parameter-estimation procedure and presents experimental verification through a proof-of-concept study. Described in greater detail in that article, the parameter estimator requires multiple forward simulations to analyze the data, which behooves minimizing the compute time for each simulation. By using an alternating direction implicit time-marching scheme and a structured but variably spaced grid, the numerical simulator built for this purpose provides a forward model output with suitable accuracy in approximately 0.5 seconds.
Energy & Fuels↗
Study of Cathode Gas Diffusion Architecture for Improved Oxygen Transport in Hydroxide Exchange Membrane Fuel Cells
The high pH environment in hydroxide exchange membrane fuel cells (HEMFCs) has the potential to reach lower costs than the current proton exchange membrane fuel cells (PEMFCs), the incumbent technology. A significant difference between HEMFCs and PEMFCs is the location of water production within the cell. In PEMFCs, the water is produced on the cathode, limiting oxygen transport. In HEMFCs, the water is produced on the anode where the fuel is pure hydrogen. This allows the cathode to be optimized for oxygen transport without the presence of excess liquid water. Limiting current analysis, a technique previously used in PEMFCs, is adopted in HEMFCs to evaluate the oxygen mass transport resistances for different sections of the cathode. Through elimination of the microporous layer (MPL), gas diffusion layer (GDL), and traditional flow field and using porous nickel foam for gas distribution, the transport resistance at an operating condition of 150 kPa(g) and with the cell temperature at 80 °C was decreased from 112 s m −1 to 48 s m −1 , effectively halved. The optimal configuration for performance was found with Ni foam and a GDL, eliminating the MPL and traditional flow field, which vastly improved oxygen transport while maintaining adequate electrical contact with the cathode catalyst layer.
25 ENERGY STORAGE↗
Constraint energy minimizing generalized multiscale finite element method for multi-continuum Richards equations
In fluid flow simulation, the multi-continuum model is a useful strategy. When the heterogeneity and contrast of coefficients are high, the system becomes multiscale, and some kinds of reduced order methods are demanded. Combining these techniques with nonlinearity, we will consider in this paper a dual-continuum model which is generalized as a multi-continuum model for a coupled system of nonlinear Richards equations as unsaturated flows, in complex heterogeneous fractured porous media; and we will solve it by a novel multiscale approach utilizing the constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM). In particular, such a nonlinear system will be discretized in time and then linearized by Picard iteration (whose global convergence is proved theoretically). Subsequently, we tackle the resulting linearized equations by the CEM-GMsFEM and obtain proper offline multiscale basis functions to span the multiscale space (which contains the pressure solution). More specifically, we first introduce two new sources of samples, and the GMsFEM is used over each coarse block to build local auxiliary multiscale basis functions via solving local spectral problems, that are crucial for detecting high-contrast channels. Second, per oversampled coarse region, local multiscale basis functions are created through the CEM as constrainedly minimizing an energy functional. Various numerical tests for our approach reveal that the error converges with the coarse-grid size and that only few oversampling layers as well as basis functions are needed.
97 MATHEMATICS AND COMPUTING↗
Electrochemical Residence Time Distribution as a Diagnostic Tool for Redox Flow Batteries
The fluid dynamic and electrochemical performance of redox flow batteries (RFBs) stems from the relationship between the flow field and the porous electrode, whose interplay determines how active species move and react during device operation. While characterization techniques, such as residence time distribution, offer insights into species mobility within a reactive volume for a traditional chemical reactor, electrochemical reactors also enable simultaneous measurement of the redox reactions, unlocking another dimension of analysis. Herein, we demonstrate how potentiodynamic measurements, using injections of electrolyte examined through moment analysis, can provide electrode-specific performance scaling relationships across a matrix of carbon paper and cloth electrodes with flow through and interdigitated flow fields. We further combine experimental campaigns with multiphysics simulations to demonstrate how electrode surface area can be estimated with this technique, which we then validate with activated and unactivated commercial carbon cloth electrodes. These studies reveal the multiscale observations that potentiodynamic measurements afford, augmenting existing electrochemical techniques for holistic electrochemical reactor diagnostics.
Electrochemistry↗
June 2022 Permeability Measurements in P-Tunnel, Nevada National Security Site
This document describes gas permeability measurements made in June 2022 in P-Tunnel at the Nevada National Security Site (NNSS). Four boreholes (AC-1, GI-5, GI-6, and HF-1) were tested with a dual packer system. The pressure versus flow rate data that were collected are analyzed using a numerical simulation method that captures compressible gas flow in partially saturated porous media using a close approximation to the geometry of the test setup. The first part of this document begins with data collection and a description of the packer assembly followed by details on data acquisition, measurement instrumentation, and equipment followed by a description of pre-test system verification. Next, the implementation of testing in the field is described, including how specific locations were chosen for testing. The second part of this document describes the numerical simulations used to back-out permeability from the pressure versus flow rate data. First, the numerical model is described, including assumptions. Next the details of the mesh that represents the experiment are described. A method for generating type curves in pressure versus flow rate space is presented and finally the data are mapped to permeability. We conclude with a discussion of the results.
54 ENVIRONMENTAL SCIENCES↗
Data-Driven Reduced-Order Modeling of Convective Heat Transfer in Porous Media
This work presents a data-driven Reduced-Order Model (ROM) for parametric convective heat transfer problems in porous media. The intrusive Proper Orthogonal Decomposition aided Reduced-Basis (POD-RB) technique is employed to reduce the porous medium formulation of the incompressible Reynolds-Averaged Navier–Stokes (RANS) equations coupled with heat transfer. Instead of resolving the exact flow configuration with high fidelity, the porous medium formulation solves a homogenized flow in which the fluid-structure interactions are captured via volumetric flow resistances with nonlinear, semi-empirical friction correlations. A supremizer approach is implemented for the stabilization of the reduced fluid dynamics equations. The reduced nonlinear flow resistances are treated using the Discrete Empirical Interpolation Method (DEIM), while the turbulent eddy viscosity and diffusivity are approximated by adopting a Radial Basis Function (RBF) interpolation-based approach. The proposed method is tested using a 2D numerical model of the Molten Salt Fast Reactor (MSFR), which involves the simulation of both clean and porous medium regions in the same domain. For the steady-state example, five model parameters are considered to be uncertain: the magnitude of the pumping force, the external coolant temperature, the heat transfer coefficient, the thermal expansion coefficient, and the Prandtl number. For transient scenarios, on the other hand, the coastdown-time of the pump is the only uncertain parameter. The results indicate that the POD-RB-ROMs are suitable for the reduction of similar problems. The relative L2 errors are below 3.34% for every field of interest for all cases analyzed, while the speedup factors vary between 54 (transient) and 40,000 (steady-state).
German, Péter (ORCID:0000000307285283)↗
Quantifying local and global mass balance errors in physics-informed neural networks
Physics-informed neural networks (PINN) have recently become attractive for solving partial differential equations (PDEs) that describe physics laws. By including PDE-based loss functions, physics laws such as mass balance are enforced softly in PINN. This paper investigates how mass balance constraints are satisfied when PINN is used to solve the resulting PDEs. We investigate PINN’s ability to solve the 1D saturated groundwater flow equations (diffusion equations) for homogeneous and heterogeneous media and evaluate the local and global mass balance errors. We compare the obtained PINN’s solution and associated mass balance errors against a two-point finite volume numerical method and the corresponding analytical solution. We also evaluate the accuracy of PINN in solving the 1D saturated groundwater flow equation with and without incorporating hydraulic heads as training data. We demonstrate that PINN’s local and global mass balance errors are significant compared to the finite volume approach. Tuning the PINN’s hyperparameters, such as the number of collocation points, training data, hidden layers, nodes, epochs, and learning rate, did not improve the solution accuracy or the mass balance errors compared to the finite volume solution. Mass balance errors could considerably challenge the utility of PINN in applications where ensuring compliance with physical and mathematical properties is crucial.
54 ENVIRONMENTAL SCIENCES↗
PFLOTRAN-SIP: A PFLOTRAN Module for Simulating Spectral-Induced Polarization of Electrical Impedance Data
Spectral induced polarization (SIP) is a non-intrusive geophysical method that collects chargeability information (the ability of a material to retain charge) in the time domain or its phase shift in the frequency domain. Although SIP is a temporal method, it cannot measure the dynamics of flow and solute/species transport in the subsurface over long times (i.e., 10–100 s of years). Data collected with the SIP technique need to be coupled with fluid flow and reactive-transport models in order to capture long-term dynamics. To address this challenge, PFLOTRAN-SIP was built to couple SIP data to fluid flow and solute transport processes. Specifically, this framework couples the subsurface flow and transport simulator PFLOTRAN and geoelectrical simulator E4D without sacrificing computational performance. PFLOTRAN solves the coupled flow and solute-transport process models in order to estimate solute concentrations, which were used in Archie’s model to compute bulk electrical conductivities at near-zero frequency. These bulk electrical conductivities were modified while using the Cole–Cole model to account for frequency dependence. Using the estimated frequency-dependent bulk conductivities, E4D simulated the real and complex electrical potential signals for selected frequencies for SIP. These frequency-dependent bulk conductivities contain information that is relevant to geochemical changes in the system. This study demonstrated that the PFLOTRAN-SIP framework is able to detect the presence of a tracer in the subsurface. SIP offers a significant benefit over ERT in the form of greater information content. It provided multiple datasets at different frequencies that better constrained the tracer distribution in the subsurface. Consequently, this framework allows for practitioners of environmental hydrogeophysics and biogeophysics to monitor the subsurface with improved resolution.
54 ENVIRONMENTAL SCIENCES↗
Geometrical-Based Generative Adversarial Network to Enhance Digital Rock Image Quality
X-ray microcomputed tomography (micro-CT) is a common tool for the study of porous media structures and properties. High-quality micro-CT data are required to accurately capture pore structures. Acquiring high-quality micro-CT data, however, is not always possible, owing to application limitations and experimental constraints. Therefore, we propose a geometrical-based generative adversarial network (GAN) to rapidly restore noisy micro-CT images to their clean counterparts. The training data and related ground-truth (GT) data are scanned for 7 min and 9.5 h, respectively. To evaluate the performance of the geometrical-based GAN, a 6003 voxel image that has never been used for training is reconstructed and compared with the corresponding GT image. Histogram matching and linear normalization are implemented to adjust the histogram of the reconstructed image to that of the GT image. A watershed-based segmentation method is then applied to delineate pore and solid phases. Lastly, we measure the Minkowski functionals and petrophysical properties, including absolute permeability, pore size distribution, drainage capillary pressure-saturation curve, and imbibition relative permeability, to estimate the physical accuracy of the denoised image. The results show that the proposed geometrical-based GAN can accurately restore noisy micro-CT data. By reducing the scanning time from 9.5 h to 7 min, the expenditure of collecting micro-CT can be decreased significantly. This is particularly important for applications where time-lapse images of a dynamic process are required, high-throughput imaging is necessary for real-time data analysis, or where the quantification of large sample volumes is required.
58 GEOSCIENCES↗