Pore network model of electrokinetic transport through charged porous media

Influence of Mesoscale Interactions on Proton, Water, and Electrokinetic Transport in Solvent-Filled Membranes: Theory and Simulation

Transport of protons and water through water-filled, phase-separated cation-exchange membranes occurs through a network of interconnected nanoscale hydrophilic aqueous domains. This paper uses numerical simulations and theory to explore the role of the mesoscale network on water, proton, and electrokinetic transport in perfluorinated sulfonic acid (PFSA) membranes, pertinent to electrochemical energy-conversion devices. Concentrated-solution theory describes microscale transport. Network simulations model mesoscale effects and ascertain macroscopic properties. An experimentally consistent 3D Voronoi-network topology characterizes the interconnected channels in the membrane. Measured water, proton, and electrokinetic transport properties from literature validate calculations of macroscopic properties from network simulations and from effective-medium theory. The results demonstrate that the hydrophilic domain size affects the various microscale, domain-level transport modes dissimilarly, resulting in different distributions of microscale coefficients for each mode of transport. As a result, the network mediates the transport of species nonuniformly with dissimilar calculated tortuosities for water, proton, and electrokinetic transport coefficients (i.e., 4.7, 3.0, and 6.1, respectively, at a water content of 8 H₂O molecules per polymer charge equivalent). The dominant water-transport pathways across the membrane are different than those taken by the proton cation. Finally, the distribution of transport properties across the network induces local electrokinetic flows that couple water and proton transport; specifically, local electrokinetic transport induces water chemical-potential gradients that decrease macroscopic conductivity by up to a factor of 3. Macroscopic proton, water, and electrokinetic transport coefficients depend on the collective microscale transport properties of all modes of transport and their distribution across the hydrophilic domain network.

Microscopic Simulations of Electrochemical Double-Layer Capacitors

Electrochemical double-layer capacitors (EDLCs) are devices allowing the storage or production of electricity. They function through the adsorption of ions from an electrolyte on high-surface-area electrodes and are characterized by short charging/discharging times and long cycle-life compared to batteries. Microscopic simulations are now widely used to characterize the structural, dynamical, and adsorption properties of these devices, complementing electrochemical experiments and in situ spectroscopic analyses. In this review, we discuss the main families of simulation methods that have been developed and their application to the main family of EDLCs, which include nanoporous carbon electrodes. We focus on the adsorption of organic ions for electricity storage applications as well as aqueous systems in the context of blue energy harvesting and desalination. We finally provide perspectives for further improvement of the predictive power of simulations, in particular for future devices with complex electrode compositions.

Nonuniform Collective Dissolution of Bubbles in Regular Pore Networks

Understanding the evolution of solute concentration gradients underpins the prediction of porous media processes limited by mass transfer. Here, we present the development of a mathematical model that describes the dissolution of spherical bubbles in two-dimensional regular pore networks. The model is solved numerically for lattices with up to 169 bubbles by evaluating the role of pore network connectivity, vacant lattice sites and the initial bubble size distribution. In dense lattices, diffusive shielding prolongs the average dissolution time of the lattice, and the strength of the phenomenon depends on the network connectivity. The extension of the final dissolution time relative to the unbounded (bulk) case follows the power-law function, $${B^k/\ell }$$ B k / ℓ , where the constant $$\ell$$ ℓ is the inter-bubble spacing, B is the number of bubbles, and the exponent k depends on the network connectivity. The solute concentration field is both the consequence and a factor affecting bubble dissolution or growth. The geometry of the pore network perturbs the inward propagation of the dissolution front and can generate vacant sites within the bubble lattice. This effect is enhanced by increasing the lattice size and decreasing the network connectivity, yielding strongly nonuniform solute concentration fields. Sparse bubble lattices experience decreased collective effects, but they feature a more complex evolution, because the solute concentration field is nonuniform from the outset.

Random pore-network model for polymer electrolyte membranes

A random pore-network model for polymer electrolyte membranes (PEM) is presented that couples the flow of protons and water through cylindrical channels to the swelling of the membrane. While the flows are determined by closed-form solutions of the Poisson-Nernst-Planck-Stokes equations, the fluid-structure interaction is described by a pressure balance at the channel walls. Macroscopic membrane properties, such as the conductivity, permeability and electro-osmotic coefficient, are computed and compared to experimental data in the literature. In light of the model simplifications, the results compare favourably to data but they also point to the importance of describing proton diffusion in PEM nanopores accurately.

Tuning the selective permeability of polydisperse polymer networks

We study the permeability and selectivity ('permselectivity') of model membranes made of polydisperse polymer networks for molecular penetrant transport, using coarse-grained, implicit-solvent computer simulations. In our work, permeability P is determined on the linear-response level using the solution-diffusion model, P = KDin, i.e., by calculating the equilibrium penetrant partition ratio K and penetrant diffusivity Din inside the membrane. We vary two key parameters, namely the network-network interaction, which controls the degree of swelling and collapse of the network, and the network-penetrant interaction, which tunes the selective penetrant uptake and microscopic energy landscape for diffusive transport. We find that the partitioning K covers four orders of magnitude and is a non-monotonic function of the parameters, well interpreted by a second-order virial expansion of the free energy of transferring one penetrant from a reservoir into the membrane. Moreover, we find that the penetrant diffusivity Din in the polydisperse networks, in contrast to highly ordered membrane structures, exhibits relatively simple exponential decays. We propose a semi-empirical scaling law for the penetrant diffusion that describes the simulation data for a wide range of densities and interaction parameters. The resulting permeability P turns out to follow the qualitative behavior (including maximization and minimization) of partitioning. However, partitioning and diffusion are typically anti-correlated, yielding large quantitative cancellations, controlled and fine-tuned by the network density and interactions, as rationalized by our scaling laws. We finally demonstrate that even small changes of network-penetrant interactions, e.g., by half a kBT, modify the permselectivity by almost one order of magnitude.

Blessing and Curse: How a Supercapacitor's Large Capacitance Causes its Slow Charging

The development of novel electrolytes and electrodes for supercapacitors is hindered by a gap of several orders of magnitude between experimentally measured and theoretically predicted charging time scales. Here, we propose an electrode model, containing many parallel stacked electrodes, that explains the slow charging dynamics of supercapacitors. At low applied potentials, the charging behavior of this model is described well by an equivalent circuit model. Conversely, at high potentials, charging dynamics slow down and evolve on two relaxation time scales: a generalized RC time and a diffusion time, which, interestingly, become similar for porous electrodes. The charging behavior of the stack-electrode model presented here helps to understand the charging dynamics of porous electrodes and qualitatively agrees with experimental time scales measured with porous electrodes.

Tuning the Permeability of Dense Membranes by Shaping Nanoscale Potentials

Permeability is one of the most fundamental transport properties in soft matter physics, material engineering, and nanofluidics. Here, we report by means of Langevin simulations of ideal penetrants in a nanoscale membrane made of a fixed lattice of attractive interaction sites, how the permeability can be massively tuned, even minimized or maximized, by tailoring the potential energy landscape for the diffusing penetrants, depending on the membrane attraction, topology, and density. Supported by limiting scaling theories we demonstrate that the observed nonmonotonic behavior and the occurrence of extreme values of the permeability is far from trivial and triggered by a strong anticorrelation and substantial (orders of magnitude) cancellation between penetrant partitioning and diffusivity, especially within dense and highly attractive membranes.

Reverse electrodialysis through nanochannels with inhomogeneously charged surfaces and overlapped electric double layers

Modeling of electro-chemo-mechanical transport phenomena in simple (nanochannel) or complex (nanoporous media) geometries with inhomogeneous surface charge and overlapped electric double layers remains challenging. This bottleneck originates from lack of a comprehensive model to predict the local surface charge density based on the variable local solution properties. This work aims to propose a model, so-called representative bulk layer (RBL), which makes the chemically non-isolated solid-liquid interfaces (due to the electric double layers interaction) as isolated interfaces by introducing a local effective bulk ion concentration. Using RBL together with the electrical triple layer model to provide boundary conditions for the multi-physio-chemical transport equations (PNP + NS), we investigate the reverse electrodialysis (RED) in nanochannels. Our modeling results indicate that the length of an ion-selective membrane not only influences the ionic current but also the logarithm of the slope of current-voltage curve increases linearly with the ratio of nanochannel length to height. This interesting finding inspires us to propose a dimensionless relation for the current-voltage curve that is independent of the nanochannel dimensions. The present contribution numerical framework could shed light on the electro-chemo-mechanical transport mechanism through nanofluidic devices and membranes.

Impact of Nanoporosity on Hydrocarbon Transport in Shales' Organic Matter

In a context of growing attention for shale gas, the precise impact of organic matter (kerogen) on hydrocarbon recovery from unconventional reservoirs still has to be assessed. Kerogen's microstructure is characterized by a very disordered pore network that greatly affects hydrocarbon transport. The specific structure and texture of this organic matter at the nanoscale is highly dependent on its origin. In this study, by the use of statistical physics and molecular dynamics, we shed some new lights on hydrocarbon transport through realistic molecular models of kerogen at different level of maturity [ Bousige et al. Nat. Mater. 2016 , 15 , 576 ]. Despite the apparent complexity, severe confinement effects controlled by the porosity of the various kerogens allow linear alkanes (from methane to dodecane) transport to be studied only via the self-diffusion coefficients of the species. The decrease of the transport coefficients with the amount of adsorbed fluid can be described by a free volume theory. Ultimately, the transport coefficients of hydrocarbons can be expressed simply as a function of the porosity (volume fraction of void) of the microstructure, thus paving the way for shale gas recovery predictions.

Upscaling scheme for long-term ion diffusion in charged porous media

Description of long-term (over years) ion diffusion at the pore scale is a huge challenge since the characteristic time of diffusion in a typical representative elementary volume is around microseconds, generally ten orders of magnitude lower than the time we were concerned with. This paper presents a numerical upscaling scheme for ion diffusion with electrical double-layer effects (electrodiffusion) considered in charged porous media. After a scaling analysis for the nondimensional governing equations of ion transport at the pore scale, we identify the conditions for decoupling of electrical effect and diffusion, and therefore are able to choose apposite temporal and spatial scales for corresponding directions of the electrodiffusion process. The upscaling scheme is therefore proposed based on a numerical framework for governing equations using a lattice Boltzmann method. The electrical potential or concentration profiles from steady- or unsteady-state electrodiffusion in the long, straight channel, calculated by this upscaling scheme, are compared with the well-meshed full-sized simulations with good agreement. Furthermore, this scheme is used to predict tracer-ion throughdiffusion and outdiffusion in hardened cement pastes. All numerical results show good agreement with the full-sized simulations or experiment data without any fitting parameters. This upscaling scheme bridges the ion diffusion behaviors in different time scales, and may help to improve the understanding of long-term ion transport mechanisms in charged porous media.

Multiscale Model for Electrokinetic Transport in Networks of Pores, Part I: Model Derivation

We present an efficient and robust numerical model for the simulation of electrokinetic phenomena in porous media and microstructure networks considering a wide range of applications including energy conversion, deionization, and microfluidic-based lab-on-a-chip systems. Coupling between fluid flow and ion transport in these networks is governed by the Poisson-Nernst-Planck-Stokes equations. These equations describe a wide range of phenomena that can interact in a complex fashion when coupled in networks involving multiple pores with variable properties. Capturing these phenomena by direct simulation of the governing equations in multidimensions is prohibitively expensive. We present here a reduced-order model that treats a network of many pores via solutions to 1D equations. Assuming that each pore in the network is long and thin, we derive a 1D model describing the transport in the pore's longitudinal direction. We take into account the cross-sectional nonuniformity of potential and ion concentration fields in the form of area-averaged coefficients in different flux terms representing fluid flow, electric current, and ion fluxes. These coefficients are obtained from the solutions to the Poisson-Boltzmann equation and are tabulated against dimensionless surface charge and dimensionless thickness of the electric double layer (EDL). Although similar models have been attempted in the past, distinct advantages of the present framework include a fully conservative discretization with zero numerical leakage, fully bounded area-averaged coefficients without any singularity in the limit of infinitely thick EDLs, a flux discretization that exactly preserves equilibrium conditions, and extension to a general network of pores with multiple intersections. In part II of this two-article series, we present a numerical implementation of this model and demonstrate its applications in predicting a wide range of electrokinetic phenomena in microstructures.

Free Volume Theory of Hydrocarbon Mixture Transport in Nanoporous Materials

Despite recent focus on shale gas, hydrocarbon recovery from the ultraconfining and disordered porosity of organic matter in shales (kerogen) remains poorly understood. Key aspects such as the breakdown of hydrodynamics at the nanoscale and strong adsorption effects lead to unexplained non-Darcy behaviors. Here, molecular dynamics and statistical mechanics are used to elucidate hydrocarbon mixture transport through a realistic molecular model of kerogen [ Bousige, C.; et al. Nat. Mater. 2016 , 15 , 576 ]. Owing to strong adsorption effects, velocity cross-correlations between the mixture components and between molecules of the same species are shown to be negligible. This allows estimation of each component permeance from its self-diffusivity, which can be obtained from single-component data. These permeances are found to scale with the reciprocal of the alkane length and decrease with the number of adsorbed molecules following a simple free volume theory, therefore allowing mixture transport prediction as a function of the amount of trapped fluid.

Electro-osmosis in inhomogeneously charged microporous media by pore-scale modeling

Surface charge at solid-electrolyte interface is generally coupled with the local electrolyte properties (ionic concentration, pH, etc.), and therefore not as assumed homogeneous on the solid surfaces in the previous studies. The inhomogeneous charge brings huge challenges in predictions of electro-osmotic transport and has never been well studied. In this work, we first propose a classification of electro-osmosis based on a dimensionless number which is the ratio of the Debye length to the characteristic pore size. In the limit of thin electrical double layer, we establish a pore-scale numerical model for inhomogeneously charged electro-osmosis including four ions: Na⁺,Cl^-,H⁺ and OH^-. Based on reconstructed porous media, we simulate the electro-osmosis with inhomogeneous charge using lattice Boltzmann method. The nonlinear response of electro-osmotic velocity to applied electrical field and the reverse flow have been observed and analyzed.

Revisiting Morrison and Osterle 1965: the efficiency of membrane-based electrokinetic energy conversion

We revisit Morrison and Osterle (1965) who derived a phenomenological expression for the 'figure-of-merit' [Formula: see text] of the electrokinetic energy conversion (EKEC) of a pressure difference into electric energy (and vice versa) using charged nanotubes, nanopores or ion-exchange membranes. We show the equivalence with Morrison and Osterle of a novel expression of [Formula: see text] derived by Bentien et al (2013). We analyze two physical models for ionic and solvent flow which directly relate [Formula: see text] to nanopore characteristics such as pore size and wall charge density. For the uniform potential model, we derive an analytical expression as a function of pore size, viscosity, ion diffusion coefficients and membrane charge density, and compare results with the full space-charge model by Osterle and co-workers as a function of pore size and ion diffusion coefficient. We present a novel expression for [Formula: see text] for salt solutions with ions with unequal diffusion coefficients (mobilities) and show that to increase [Formula: see text] the counterion mobility must be low and the coion mobility high.

Lattice Boltzmann model for Coulomb-driven flows in dielectric liquids

In this paper, we developed a unified lattice Boltzmann model (LBM) to simulate electroconvection in a dielectric liquid induced by unipolar charge injection. Instead of solving the complex set of coupled Navier-Stokes equations, the charge conservation equation, and the Poisson equation of electric potential, three consistent lattice Boltzmann equations are formulated. Numerical results are presented for both strong and weak injection regimes, and different scenarios for the onset and evolution of instability, bifurcation, and chaos are tracked. All LBM results are found to be highly consistent with the analytical solutions and other numerical work.

Multiscale adsorption and transport in hierarchical porous materials

This review presents the state-of-the-art of multiscale adsorption and transport in hierarchical porous materials.