An official website of the United States government
Ion transport within saturated porous media is an intricate process in which efficient ion delivery is desired in many engineering problems. However, controlling the behavior of ion transport proves challenging, as ion transport is influenced by a variety of driving mechanisms, which requires a systematic understanding. Herein, we study a coupled advection–diffusion–electromigration system for controlled ion transport within porous media using the scaling analysis. Using the Lattice–Boltzmann–Poisson method, we establish a transport regime classification based on an Advection Diffusion Index (ADI) and a novel Electrodiffusivity Index (EDI) for a two-dimensional (2D) microchannel model under various electric potentials, pressure gradients, and concentration conditions. The resulting transport regimes can be well controlled by changing the applied electric potential, the pressure field, and the injected ions concentration. Furthermore, we conduct numerical simulations in a synthetic 2D porous media and an x-ray microcomputed tomography sandstone image to validate the prevailing transport regime. The simulation results highlight that the defined transport regime observed in our simple micromodel domain is also observed in the synthetic two- and three-dimensional domains, but the boundary between each transport regime differs depending on the variation of the pore size within a given domain. Consequently, the proposed ADI and EDI emerge as dimensionless indicators for controlled ion transport. Overall, our proof-of-concept for ion transport control in porous media is demonstrated under advection–diffusion–electromigration transport, demonstrating the richness of transport regimes that can develop and provide future research directions for subsurface engineering applications.
more »
« less
A Pore-Scale Model for Electrokinetic In situ Recovery of Copper: The Influence of Mineral Occurrence, Zeta Potential, and Electric Potential
AbstractElectrokinetic in-situ recovery is an alternative to conventional mining, relying on the application of an electric potential to enhance the subsurface flow of ions. Understanding the pore-scale flow and ion transport under electric potential is essential for petrophysical properties estimation and flow behavior characterization. The governing physics of electrokinetic transport is electromigration and electroosmotic flow, which depend on the electric potential gradient, mineral occurrence, domain morphology (tortuosity and porosity, grain size and distribution, etc.), and electrolyte properties (local pH distribution and lixiviant type and concentration, etc.). Herein, mineral occurrence and its associated zeta potential are investigated for EK transport. The new Ek model which is designed to solve the EK flow in complex porous media in a highly parallelizable manner includes three coupled equations: (1) Poisson equation, (2) Nernst–Planck equation, and (3) Navier–Stokes equation. These equations were solved using the lattice Boltzmann method within X-ray computed microtomography images. The proposed model is validated against COMSOL multiphysics in a two-dimensional microchannel in terms of fluid flow behavior when the electrical double layer is both resolvable and unresolvable. A more complex chalcopyrite-silica system is then obtained by micro-CT scanning to evaluate the model performance. The effects of mineral occurrence, zeta potential, and electric potential on the three-dimensional chalcopyrite-silica system were evaluated. Although the positive zeta potential of chalcopyrite can induce a flow of ferric ion counter to the direction of electromigration, the net effect is dependent on the occurrence of chalcopyrite. However, the ion flux induced by electromigration was the dominant transport mechanism, whereas advection induced by electroosmosis made a lower contribution. Overall, a pore-scale EK model is proposed for direct simulation on pore-scale images. The proposed model can be coupled with other geochemical models for full physicochemical transport simulations. Meanwhile, electrokinetic transport shows promise as a human-controllable technique because the electromigration of ions and the applied electric potential can be easily controlled externally. Graphical abstract
more »
« less
- Award ID(s):
- 2324787
- PAR ID:
- 10536456
- Publisher / Repository:
- Springer
- Date Published:
- Journal Name:
- Transport in Porous Media
- Volume:
- 150
- Issue:
- 3
- ISSN:
- 0169-3913
- Page Range / eLocation ID:
- 601 to 626
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract This paper reports the development of a numerical solver aimed to simulate the interaction between the space charge (i.e. ions) distribution and the electric field in liquid argon time projection chamber (LArTPC) detectors. The ion transport equation is solved by a time-accurate, cell-centered finite volume method and the electric potential equation by a continuous finite element method. The electric potential equation updates the electric field which provides the drift velocity to the ion transport equation. The ion transport equation updates the space charge density distribution which appears as the source term in the electric potential equation. The interaction between the space charge distribution and the electric field is numerically simulated within each physical time step. The convective velocity in the ion transport equation can include the background flow velocity in addition to the electric drift velocity. The numerical solver has been parallelized using the Message Passing Interface (MPI) library. Numerical tests show and verify the capability and accuracy of the current numerical solver. It is planned that the developed numerical solver, together with a Computational Fluid Dynamics (CFD) package which provides the flow velocity field, can be used to investigate the space charge effect on the electric field in large-scale particle detectors.more » « less
-
A model is constructed to describe the flow field and arbitrary deformation of a drop or vesicle that contains and is embedded in an electrolyte solution, where the flow and deformation are caused by an applied electric field. The applied field produces an electrokinetic flow, which is set up on the charge-up time scale $$\tau _{*c}=\lambda _{*} a_{*}/D_{*}$$ , where $$\lambda _{*}$$ is the Debye screening length, $$a_{*}$$ is the inclusion length scale and $$D_{*}$$ is an ion diffusion constant. The model is based on the Poisson–Nernst–Planck and Stokes equations. These are reduced or simplified by forming the limit of strong electrolytes, for which dissolved salts are completely ionised in solution, together with the limit of thin Debye layers. Debye layers of opposite polarity form on either side of the drop interface or vesicle membrane, together forming an electrical double layer. Two formulations of the model are given. One utilises an integral equation for the velocity field on the interface or membrane surface together with a pair of integral equations for the electrostatic potential on the outer faces of the double layer. The other utilises a form of the stress-balance boundary condition that incorporates the double layer structure into relations between the dependent variables on the layers’ outer faces. This constitutes an interfacial boundary condition that drives an otherwise unforced Stokes flow outside the double layer. For both formulations relations derived from the transport of ions in each Debye layer give additional boundary conditions for the potential and ion concentrations outside the double layer.more » « less
-
The accurate measurement of wall zeta potentials and solute–surface interaction length scales for electrolyte and non-electrolyte solutes, respectively, is critical to the design of many biomedical and microfluidic applications. We present a novel microfluidic approach using diffusioosmosis for measuring either the zeta potentials or the characteristic interaction length scales for surfaces exposed to, respectively, electrolyte or non-electrolyte solutes. When flows containing different solute concentrations merge in a junction, local solute concentration gradients can drive diffusioosmotic flow due to electrokinetic, steric, and other interactions between the solute molecules and solid surfaces. We demonstrate a microfluidic system consisting of a long, narrow pore connecting two large side channels in which solute concentration gradients drive diffusioosmosis within the pore, resulting in predictable fluid velocity/pressure and solute profiles. Furthermore, we present analytical results and a methodology to determine the zeta potential or interaction length scale for the pore surfaces based on the solute concentrations in the main side channels, the flow rate in the pore, and the pressure drop across the pore. We apply this method to the experimental data of Lee et al. to predict the zeta potentials of their system, and we use 3D numerical simulations to validate the theory and show that end effects caused by the junctions are negligible for a wide range of parameters. Because the dynamics in the proposed system are driven by diffusioosmosis, this technique does not suffer from certain disadvantages associated with the use of sensitive electronics in traditional zeta potential measurement approaches such as streaming potential, streaming current, or electroosmosis. To the best of our knowledge this is the first flow-based approach to characterize surface/solute interactions with non-electrolyte solutes.more » « less
-
Ion concentration polarization (ICP) accomplishes preconcentration for bioanalysis by localized depletion of electrolyte ions, thereby generating a gradient in electric field strength that facilitates electrokinetic focusing of charged analytes by their electromigration against opposing fluid flow. Such ICP focusing has been shown to accomplish up to a million-fold enrichment of nucleic acids and proteins in single-stage preconcentrators. However, the rate at which the sample volume is swept is limited, requiring several hours to achieve these high enrichment factors. This limitation is caused by two factors. First, an ion depleted zone (IDZ) formed at a planar membrane or electrode may not extend across the full channel cross section under the flow rate employed for focusing, thereby allowing the analyte to “leak” past the IDZ. Second, within the IDZ, large fluid vortices lead to mixing, which decreases the efficiency of analyte enrichment and worsens with increased channel dimensions. Here, we address these challenges with faradaic ICP (fICP) at a three-dimensional (3D) electrode comprising metallic microbeads. This 3D-electrode distributes the IDZ, and therefore, the electric field gradient utilized for counter-flow focusing across the full height of the fluidic channel, and its large area, microstructured surface supports smaller vortices. An additional bed of insulating microbeads restricts flow patterns and supplies a large area for surface conduction of ions through the IDZ. Finally, the resistance of this secondary bed enhances focusing by locally strengthening sequestering forces. This easy-to-build platform lays a foundation for the integration of enrichment with user-defined packed bed and electrode materials.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1007/s11242-023-02023-2
