skip to main content


Title: Pore‐scale Study of Ion Transport Mechanisms in Inhomogeneously Charged Nanoporous Rocks: Impacts of Interface Properties on Macroscopic Transport
Abstract

The electrokinetic transport mechanisms of multispecies ions through 3‐D nanoporous rocks with chemical reaction at the solid‐aqueous solution interfaces are investigated. We systematically study the multiphysics transport phenomena by considering either inhomogeneous (local surface charge based on local pH and ion concentrations) or prescribed homogeneous surface charge at solid‐aqueous solution interface while the pores are screened via electric double layers. We develop a lattice Boltzmann numerical framework to solve the set of governing equations (Poisson‐Nernst‐Planck plus Navier–Stokes). Our modeling results reveal that the averaged local electric potential of the nanoporous rock is significantly underestimated (about 83%) when a homogeneous surface charge is prescribed based on the bulk solution properties. It is shown that increasing the porosity of the nanoporous media considerably increases the absolute values and inhomogeneity of the surface charge, which means that while the electric double layers screened the pores, increasing the porosity enhances the ion selectivity of the porous medium. When the scenario with inhomogeneous charge is taken into account, the predicted electroosmotic permeability and tortuosity are higher in comparison with the prescribed homogeneous case. Moreover, we have studied the electrostatic tortuosity, coupling coefficient, and the effective excess charge density of the nanoporous rocks. The results demonstrate that ignoring the inhomogeneity of surface charges may cause erroneous prediction of the ion transport through porous rocks with chemically active surfaces.

 
more » « less
NSF-PAR ID:
10457245
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
124
Issue:
6
ISSN:
2169-9313
Page Range / eLocation ID:
p. 5387-5407
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Advances in solution-phase graphene patterning has provided a facile route for rapid, low-cost and scalable manufacturing of electrochemical devices, even on flexible substrates. While graphene possesses advantageous electrochemical properties of high surface area and fast heterogenous charge transport, these properties are attributed to the edge planes and defect sites, not the basal plane. Herein, we demonstrate enhancement of the electroactive nature of patterned solution-phase graphene by increasing the porosity and edge planes through the construction of a multidimensional architecture via salt impregnated inkjet maskless lithography (SIIML) and CO 2 laser annealing. Various sized macroscale pores (<25 to ∼250 μm) are patterned directly in the graphene surface by incorporating porogens ( i.e. , salt crystals) in the graphene ink which act as hard templates for pore formation and are later dissolved in water. Subsequently, microsized pores (∼100 nm to 2 μm in width) with edge plane defects are etched in the graphene lattice structure by laser annealing with a CO 2 laser, simultaneously improving electrical conductivity by nearly three orders of magnitude (sheet resistance decreases from >10 000 to ∼50 Ω sq −1 ). We demonstrate that this multidimensional porous graphene fabrication method can improve electrochemical device performance through design and manufacture of an electrochemical organophosphate biosensor that uses the enzyme acetylcholinesterase for detection. This pesticide biosensor exhibits enhanced sensitivity to acetylthiocholine compared to graphene without macropores (28.3 μA nM −1 to 13.3 μA nM −1 ) and when inhibited by organophosphate pesticides (paraoxon) has a wide linear range (10 nM to 500 nM), low limit of detection (0.6 nM), and high sensitivity (12.4 nA nM −1 ). Moreover, this fabrication method is capable of patterning complex geometries [ i.e. interdigitated electrodes (IDEs)] even on flexible surfaces as demonstrated by an IDE supercapacitor made of SIIML graphene on a heat sensitive polymer substrate. The supercapacitor demonstrates a high energy density of 0.25 mW h cm −3 at a power density of 0.3 W cm −3 . These electrochemical devices demonstrate the benefit of using SIIML and CO 2 laser annealing for patterning graphene electrodes with a multidimensional porous surface even on flexible substrates and is therefore a platform technology which could be applied to a variety of different biosensors and other electrochemical devices. 
    more » « less
  2. null (Ed.)
    In this paper, we consider an important problem for modeling complex coupled phenomena in porous media at multiple scales. In particular, we consider flow and transport in the void space between the pores when the pore space is altered by new solid obstructions formed by microbial growth or reactive transport, and we are mostly interested in pore-coating and pore-filling type obstructions, observed in applications to biofilm in porous media and hydrate crystal formation, respectively. We consider the impact of these obstructions on the macroscopic properties of the porous medium, such as porosity, permeability and tortuosity, for which we build an experimental probability distribution with reduced models, which involves three steps: (1) generation of independent realizations of obstructions, followed by, (2) flow and transport simulations at pore-scale, and (3) upscaling. For the first step, we consider three approaches: (1A) direct numerical simulations (DNS) of the PDE model of the actual physical process called BN which forms the obstructions, and two non-DNS methods, which we call (1B) CLPS and (1C) LP. LP is a lattice Ising-type model, and CLPS is a constrained version of an Allen–Cahn model for phase separation with a localization term. Both LP and CLPS are model approximations of BN, and they seek local minima of some nonconvex energy functional, which provide plausible realizations of the obstructed geometry and are tuned heuristically to deliver either pore-coating or pore-filling obstructions. Our methods work with rock-void geometries obtained by imaging, but bypass the need for imaging in real-time, are fairly inexpensive, and can be tailored to other applications. The reduced models LP and CLPS are less computationally expensive than DNS, and can be tuned to the desired fidelity of the probability distributions of upscaled quantities. 
    more » « less
  3. null (Ed.)
    The focus of this paper is a numerical simulation study of the flow dynamics in a periodic porous medium to analyse the physics of a symmetry-breaking phenomenon, which causes a deviation in the direction of the macroscale flow from that of the applied pressure gradient. The phenomenon is prominent in the range of porosity from 0.43 to 0.72 for circular solid obstacles. It is the result of the flow instabilities formed when the surface forces on the solid obstacles compete with the inertial force of the fluid flow in the turbulent regime. We report the origin and mechanism of the symmetry-breaking phenomenon in periodic porous media. Large-eddy simulation (LES) is used to simulate turbulent flow in a homogeneous porous medium consisting of a periodic, square lattice arrangement of cylindrical solid obstacles. Direct numerical simulation is used to simulate the transient stages during symmetry breakdown and also to validate the LES method. Quantitative and qualitative observations are made from the following approaches: (1) macroscale momentum budget and (2) two- and three-dimensional flow visualization. The phenomenon draws its roots from the amplification of a flow instability that emerges from the vortex shedding process. The symmetry-breaking phenomenon is a pitchfork bifurcation that can exhibit multiple modes depending on the local vortex shedding process. The phenomenon is observed to be sensitive to the porosity, solid obstacle shape and Reynolds number. It is a source of macroscale turbulence anisotropy in porous media for symmetric solid-obstacle geometries. In the macroscale, the principal axis of the Reynolds stress tensor is not aligned with any of the geometric axes of symmetry, nor with the direction of flow. Thus, symmetry breaking in porous media involves unresolved flow physics that should be taken into consideration while modelling flow inhomogeneity in the macroscale. 
    more » « less
  4. null (Ed.)
    Highly efficient capacitive deionization (CDI) relies on unimpeded transport of salt ions to the electrode surface. Graphene is an ideal candidate to provide superb conditions for ion adsorption as it possesses high theoretical surface area and electrical conductivity. When ions are stored solely within the electric double layers (EDLs), a hydrophilic graphene surface with hierarchical pores can maximize the accessible surface area and promote the ion transport. In the case of synergistic ion storage via electrostatic adsorption and faradaic redox reaction, graphene can act as both the electron highway and the reciprocal spacer to provide surface-confined effects. Substantially, structural and chemical engineering towards graphene can enhance the ion removal capacity and rate, and improve the charge efficiency and ion selectivity. In this review, we keep pace with the in-depth studies of CDI technologies and recent progress on graphene-based materials for CDI. Major challenges in the rational assembly of the desired material functionalities in terms of surface area, pore structure, and hydrophilicity are addressed. As electrode materials develop, the ultimate goal is to achieve highly efficient, energy-saving, and environment-friendly CDI. 
    more » « less
  5. Electro-osmotic flow (EOF) is a phenomenon where fluid motion occurs in porous materials or micro/nano-channels when an external electric field is applied. In the particular example of single-molecule electrophoresis using single nanopores, the role of EOF on the translocation velocity of the analyte molecule through the nanopore is not fully understood. The complexity arises from a combination of effects from hydrodynamics in restricted environments, electrostatics emanating from charge decorations and geometry of the pores. We address this fundamental issue using the Poisson–Nernst–Planck and Navier–Stokes (PNP–NS) equations for cylindrical solid-state nanopores and three representative protein nanopores (α-hemolysin, MspA, and CsgG). We present the velocity profiles inside the nanopores as a function of charge decoration and geometry of the pore and applied electric field. We report several unexpected results: (a) The apparent charges of the protein nanopores are different from their net charge and the surface charge of the whole protein geometry, and the net charge of inner surface is consistent with the apparent charge. (b) The fluid velocity depends non-monotonically on voltage. The three protein nanopores exhibit unique EOF and velocity–voltage relations, which cannot be simply deduced from their net charge. Furthermore, effective point mutations can significantly change both the direction and the magnitude of EOF. The present computational analysis offers an opportunity to further understand the origins of the speed of transport of charged macromolecules in restricted space and to design desirable nanopores for tuning the speed of macromolecules through nanopores. 
    more » « less