Next Article in Journal
Conductive Nanofilms with Oppositely Charged Reduced Graphene Oxides as a Base for Electroactive Coatings and Sensors
Next Article in Special Issue
Driven Transport of Dilute Polymer Solutions through Porous Media Comprising Interconnected Cavities
Previous Article in Journal
Recent Advances in Nanomaterials for Dermal and Transdermal Applications
Previous Article in Special Issue
Numerical Simulation of the Distribution of In Situ Cigarette Combustion-Generated Particulate Matter
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Electroosmosis and Electric Conduction of Electrolyte Solutions in Charge-Regulating Fibrous Media

Department of Chemical Engineering, National Taiwan University, Taipei 10617, Taiwan
*
Author to whom correspondence should be addressed.
Colloids Interfaces 2021, 5(1), 19; https://0-doi-org.brum.beds.ac.uk/10.3390/colloids5010019
Submission received: 15 February 2021 / Revised: 14 March 2021 / Accepted: 15 March 2021 / Published: 19 March 2021
(This article belongs to the Special Issue Locomotion of Colloidal Particles)

Abstract

:
An analytical study of the electroosmosis and electric conduction of electrolyte solutions in a fibrous medium composed of parallel charge-regulating cylinders with arbitrary electric double layer thickness is presented. A linearized charge regulation model was adopted for the association and dissociation reactions occurring at the amphoteric functional groups over the surfaces of the cylinders, and a unit cell model was employed to accommodate interactions among the cylinders. The electrokinetic equations governing the ionic concentration, electric potential, and liquid flow fields were solved at low zeta potential for the cylinders. Explicit formulas for the electroosmotic mobility and effective electric conductivity in the fiber matrix were obtained. The results indicate that the charge regulation characteristics, such as the equilibrium constants of the reactions occurring at the cylinders’ surfaces and the bulk concentration of the charge-determining ions, influence the surface charge density and potential, electroosmotic mobility, and effective electric conductivity substantially.

1. Introduction

When a charged narrow channel or porous medium filled with an electrolyte solution is subject to an external electric field, an electroosmotic flow of the fluid occurs. Electroosmotic flows are of fundamental and practical interest in a variety of scientific and technological areas. The basic formulas of electroosmosis have been derived mainly for single microchannels or nanochannels [1,2,3,4,5].
In practical applications of electroosmosis, porous media are usually encountered. To avoid the difficulty of multiple-entity boundary value problems, the unit cell model was frequently used to predict the effects of interactions among the small solid entities constructing the porous medium on the electroosmotic velocity and effective electric conductivity. Although different shapes of the cell may be employed, a spherical [6,7,8,9,10,11,12,13,14] or circular cylindrical [15,16,17,18,19,20] one for the fictitious envelope of fluid surrounding each spherical or circular cylindrical entity is most convenient.
The previous studies of the electroosmotic velocity and electric conductivity in porous media, it was often assumed that the charge density or electric potential at the solid surface is constant. However, the charge density and electric potential of most biological and some inorganic surfaces are governed by some association and dissociation reactions, and the number density of their ionogenic functional groups. The extent of these reactions depends on the local concentration of the charge-determining or potential-determining ions, and other relevant parameters [21,22,23,24,25,26,27]. When such a charge-regulating surface is subject to an external electric field, the electric potential and ionic concentration distributions around will be influenced by the charge regulation characteristics and the fluid adjacent to the surface undergoes electroosmosis accordingly.
The electroosmosis and electric conduction of electrolyte solutions in long capillary tubes or slits whose inside walls were covered by a layer of adsorbed charge-regulating polymer of arbitrary thickness were analytically studied for the case of an arbitrary electric double layer, and closed-form formulas for the local and average electroosmotic velocities and electric current density were obtained for a linearized charge regulation model [28]. Recently, the electroosmotic mobility and ionic conductance in a long charge-regulating nanochannel of arbitrary cross-section but low surface potential have been semi-analytically examined, and the geometrical configuration of the nanochannel was shown to play a key role in the determination of these transport properties [29]. On the other hand, the electroosmosis/electrophoresis and electric conduction in a suspension or porous system of an assemblage of charge-regulating spherical entities have been investigated using the cell model, and explicit formulas for the electroosmotic/electrophoretic mobility and effective electric conductivity were derived as functions of the porosity and other relevant parameters [30].
In this paper, the unit cell model is used to analyze the electroosmotic mobility and effective electric conductivity in a fiber matrix constructed by a uniform array of parallel charge-regulating circular cylinders surrounded by arbitrary double layers. The closed-form results enable these transport properties to be easily predicted for various values of the charge regulation characteristics and porosity of the fibrous system.

2. Analysis

Consider an arbitrary fluid solution having M ionic species filled in a fiber matrix constructed by a homogeneous array of parallel, circular, dielectric cylinders. Each cylinder is long and has a charge-regulating surface whose ionogenic functional groups are in chemical equilibrium with the ambient solution. When the fiber matrix is subject to a transversely imposed electric field E e x , the fluid undergoes electroosmosis with a superficial velocity U e x , where e x is the unit vector in the x direction. As illustrated in Figure 1, we use a unit cell model in which each cylinder with radius a is enveloped by a coaxial annulus of the fluid with outer radius b , and the fluid-to-cell volume ratio equals the porosity 1 ϕ of the fibrous system; i.e., ϕ = ( a / b ) 2 . The origin of the polar coordinate system ( ρ ,   φ ) is set at the cylinder axis with x = ρ cos φ .

2.1. Governing Equations

The magnitude of the fluid velocity is taken to be small, such that the electric double layer encompassing the cylinder is just distorted slightly from equilibrium. Thus, the electric potential distribution ψ ( ρ , φ ) , the concentration distribution n m ( ρ , φ ) of the ionic species m , and the pressure distribution p ( ρ , φ ) may be written as in Equation (1):
ψ = ψ ( eq ) + δ ψ         ,   n m = n m ( eq ) + δ n m         ,   p = p ( eq ) + δ p ,
where ψ ( eq ) ( ρ ) , n m ( eq ) ( ρ ) , and p ( eq ) ( ρ ) represent the equilibrium profiles of the electric potential, concentration of the species m , and pressure, respectively (with n m ( eq ) relating to ψ ( eq ) by the Boltzmann equation), and δ ψ ( ρ ,   φ ) , δ n m ( ρ ,   φ ) , and δ p ( ρ ,   φ ) are the small perturbations to equilibrium. Note that δ ψ ( ρ ,   φ ) contains the electric potentials contributed by both the applied electric field itself and the charge regulation effect at the cylinder’s surface induced by the applied electric field.
The perturbed quantities δ ψ , δ n m , and δ p along with the fluid velocity distribution u ( ρ , φ ) satisfy the continuity equation u = 0 and the following linearized electrokinetic equations [10]:
2 δ ψ = e ε k T m = 1 M z m n m exp ( z m e ψ ( eq ) k T ) ( δ μ m z m e δ ψ ) ,
2 δ μ m = z m e k T [ ψ ( eq ) δ μ m k T D m ψ ( eq ) u ] ,   m = 1 ,   2 ,   ,   M ,
η 2 u = δ p ε [ 2 ψ ( eq ) δ ψ + 2 δ ψ   ψ ( eq ) ] .
In Equations (2) and (3), δ μ m is the perturbed ionic electrochemical potential energy defined by
δ μ m = z m e δ ψ + k T n m ( eq ) δ n m ,
n m , D m , and z m are the bulk concentration, diffusion coefficient, and valence, respectively, of the species m ; η and ε are the viscosity and dielectric permittivity, respectively, of the fluid; e is the charge of a proton; k is the Boltzmann constant; T is the absolute temperature. Note that the first parts of Equation (1) and Equation (2) may not be needed to solve for the linear electroosmotic velocity (in terms of the applied electric field strength) of the fluid if the charge regulation effect is absent at the cylinder surface.

2.2. Boundary Conditions

The boundary conditions for δ μ m and u at the surface of the cylinder are
ρ = a : δ μ m ρ = 0 ,   u = 0 ,
which are obtained from the assumptions that the shear plane coincides with the cylinder surface and no ions can penetrate into the cylinder. To obtain the boundary condition for the perturbed quantity δ ψ at the charge-regulating surface, we adopt the linearized regulation model proposed by Carnie and Chan [22] and express the surface charge density σ as a linear function of the surface potential ψ S ,
σ = σ ( eq ) + ( d σ d ψ S ) ψ S = ζ   δ ψ S ,  
where the constants σ ( eq ) and ζ are the values of σ and ψ S , respectively, at equilibrium. Note that ζ does not coincide with the actual interfacial electric potential, which is not assumed to be constant when the external electric field is present. The substitution of Equations (1) and (7) into the Gauss condition at the cylinder surface,
ρ = a : ψ ρ = σ ε ,  
results in
ρ = a : δ ψ ρ = L δ ψ = 0 ,  
where the charge regulation coefficient L is defined and related to measurable quantities by Equation (A6) in Appendix A for the mechanism of charge regulation. The limits L and L = 0 reduce to the simple cases of constant surface potential and constant surface charge density, respectively.
The boundary conditions at the outer (virtual) surface of the unit cell may be
ρ = b : δ ψ ρ = E cos Φ ,
δ μ m ρ = z m e E cos Φ ,
u ρ =   U cos Φ ,
τ ρ φ = η [ ρ ρ ( u φ ρ ) + 1 ρ u ρ Φ ] = 0
for the Happel model;
( × u ) z = 1 ρ ρ ( ρ u φ ) 1 ρ u ρ Φ = 0
for the Kuwabara model, where u ρ and u Φ are the ρ and Φ components, respectively, of u and U is the superficial electroosmotic velocity of the fluid (may also be taken as the electrophoretic velocity of the circular cylinder in the opposite direction in a quiescent fluid) to be determined. In Equations (10a) and (10b), the local electric field consistent with the imposed electric field E e x is set according to the Levine-Neale electrostatic model [6].
Alternatively, the electric potential at the outer surface of the cell may be set as the profile resulting in the imposed electric field E e x , known as the Zharkikh–Shilov electrostatic model [7], and the boundary conditions in Equations (10a) and (10b) become
ρ = b :   δ ψ = E ρ cos Φ ,   δ μ m = z m e E ρ cos Φ .
Note that the electric field strength at infinity E in the previous equation might be replaced more appropriately by the volume-mean electric field strength [12].

2.3. Equilibrium Electric Potential

To solve for the electroosmotic velocity in a unit cell, we need the equilibrium electric potential distribution ψ ( eq ) satisfying the Poisson–Boltzmann equation and the boundary conditions
ρ = a :   ψ ( eq ) = ζ ,
ρ = b :   d ψ ( eq ) d ρ = 0 .
For the case of small equilibrium surface (zeta) potential ζ of the cylinder, it can be shown that
ψ ( eq ) = ψ eq 1 ( ρ ) ζ ¯
to the order of ζ ¯ , where ζ ¯ = ζ e / k T ,
ψ eq 1 ( ρ ) = k T e I 1 ( κ b ) K 0 ( κ ρ ) + I 0 ( κ ρ ) K 1 ( κ b ) I 1 ( κ b ) K 0 ( κ a ) + I 0 ( κ a ) K 1 ( κ b ) ,
κ is the Debye screening parameter equal to e ( m = 1 M z m 2 n m / ε   k T   ) 1 / 2 , and I n and K n are the modified Bessel functions of order n of the first and second kinds, respectively.
By substituting Equations (14) and (15) into the Gauss condition, one obtains the following formula relating the zeta potential and surface charge density of the cylinder at equilibrium:
σ ( eq ) = ε κ ζ I 1 ( κ b ) K 1 ( κ a ) I 1 ( κ a ) K 1 ( κ b ) I 1 ( κ b ) K 0 ( κ a ) + I 0 ( κ a ) K 1 ( κ b ) .
The zeta potential ζ of the charge-regulating cylinder can be evaluated in terms of measurable quantities by using Equations (16) and (A3) (with δ μ s = 0 at equilibrium). For the limiting case of ϕ = ( a / b ) 2 = 0 , Equation (16) reduces to the simple relation σ ( eq ) = ε κ ζ K 1 ( κ a ) / K 0 ( κ a ) for an isolated charged cylinder.

2.4. Perturbed Quantities and Electroosmotic Velocity

The small perturbed quantities δ ψ , δ μ m , v r , v θ , and δ p satisfying Equations (2)–(11) and the continuity equation v = 0 can be solved in terms of ψ eq 1 in Equation (15), and the results to the sufficient (zeroth or first) orders of ζ ¯ for determining the electroosmotic velocity in the fibrous medium correct to the order of ζ ¯ are
δ ψ = E F ψ ( ρ ) cos ϕ ,
δ μ m = E [ F m 0 + ζ ¯ F m 1 ( ρ ) ] cos ϕ ,
u ρ = E ζ ¯ F ρ ( ρ ) cos ϕ ,
u φ = E ζ ¯ F φ ( ρ ) sin ϕ ,
δ p = E ζ ¯ [ η a F p ( ρ ) + ε κ 2 F ψ ( ρ ) ψ eq 1 ( ρ ) ] cos ϕ ,
For these solutions, the functions F m 1 ( ρ ) , F ρ ( ρ ) , F ϕ ( ρ ) , and F p ( ρ ) were given in [18] but with ψ eq 1 ( ρ ) expressed by Equation (15):
F m 0 = z m e ρ χ ( 1 + a 2 ρ 2 ) ,
F ψ ( ρ ) = I 1 ( κ ρ ) { A K ( ρ ) 2 h K + κ g K [ h K A I ( a ) + h I A K ( a ) ] κ [ g I h K h I g K ] } K 1 ( κ ρ ) { A I ( ρ )
+ 2 h I + κ g I [ h K A I ( a ) + h I A K ( a ) ] κ [ g I h K h I g K ] }
if Equation (10a) is used,
F ψ ( ρ ) = I 1 ( κ ρ ) { A K ( ρ ) b h K K 1 ( κ b ) [ h K A I ( a ) + h I A K ( a ) ] I 1 ( κ b ) h K + h I K 1 ( κ b ) } K 1 ( κ ρ ) { A I ( ρ )
+ b h I + I 1 ( κ b ) [ h K A I ( a ) + h I A K ( a ) ] I 1 ( κ b ) h K + h I K 1 ( κ b ) }
if Equation (11) is used, where
A I ( ρ ) = 1 χ ρ b ( κ ρ ) 2 I 1 ( κ ρ ) [ 1 + ( a ρ ) 2 ] d ρ ,
A K ( ρ ) = 1 χ ρ b ( κ ρ ) 2 K 1 ( κ ρ ) [ 1 + ( a ρ ) 2 ] d ρ ,
g I = I 0 ( κ b ) + I 2 ( κ b ) ,       g K = K 0 ( κ b ) + K 2 ( κ b ) ,
h I = κ a I 0 ( κ a ) ( 1 + a L ) I 1 ( κ a ) ,       h K = κ a K 0 ( κ a ) + ( 1 + a L ) K 1 ( κ a ) ,
and the coefficient χ equals 1 ϕ for the Levine–Neale model in Equations (10a) and (10b) and 1 + ϕ for the Zharkikh–Shilov model in Equation (11). Note that the perturbed electric potential δ ψ given by Equations (17) and (21) corrected by the charge regulation effect at the cylinder surface induced by the applied electric field does not affect the functions F m 0 ( ρ ) , F m l ( ρ ) , F ρ ( ρ ) , F ϕ ( ρ ) , and F ρ ( ρ ) for the perturbed ionic electrochemical potential energy δ μ m and the fluid flow, but it does contribute to the perturbed pressure δ p via Equation (19c) and influence the force balance on the unit cell. The electric potential distribution along the cylinder surface ψ s may be determined from Equations (1), (14), (15), (17), and (21) by taking ρ = a , and then the surface charge density distribution σ is obtained from Equation (A3).
The hydrodynamic drag force exerted on the unit cell at its virtual surface per unit length is
F h = 2 b 0 π { δ p e ρ + η [ u + ( u ) T ] e ρ } ρ = b d φ .
Since the whole cell is electrically neutral, both the electrostatic force and the net (hydrodynamic) force exerted on it must be zero. By applying this constraint to Equation (25) (with the substitution of Equation (19)), we obtain the electroosmotic velocity of the fluid correct to the order of ζ ¯ as
U = ε E ( κ a ) 2 ζ ¯ 8 η ( 1 + ϕ 2 ) χ { χ F ψ ( b ) ψ eq 1 ( b ) a ϕ 1 / 2 ω   a   b ( 1 + a 2 ρ 2 ) [ 1 ρ 2 a 2 ( 1 2 ln ρ a )
+ ϕ 2 ρ 2 a 2 ( 1 + 2 ln ρ a ρ 2 a 2 ) ] d ψ eq 1 d ρ d ρ }
for the Happel hydrodynamic model, and
U = ε E ( κ a ) 2 ζ ¯ 8 η χ { 3 χ F ψ ( b ) ψ eq 1 ( b ) 2 a ϕ 1 / 2 ω   a   b ( 1 + a 2 ρ 2 ) [ 1 ρ 2 a 2 ( 1 2 ln ρ a )
1 2 ϕ   ( 1 2 ρ 2 a 2 + ρ 4 a 4 ) ] d ψ eq 1 d ρ d ρ }
for the Kuwabara hydrodynamic model, where
ω = ( 1 ϕ 2 + ln ϕ + ϕ 2 ln ϕ ) 1 ,       ω = ( 1 4 3 ϕ + 1 3 ϕ 2 + 2 3 ln ϕ ) 1 .
When Equation (11) for the Zharkikh–Shilov electrostatic model is used, Equation (26) (with χ = 1 + ϕ ) for U will not depend on the charge regulation coefficient L . Equation (26) can also be used to express the electrophoretic velocity of a suspension of parallel circular cylindrical particles in the direction normal to their axes.

2.5. Electric Conductivity

For an array of parallel charge-regulating circular cylinders subject to a transversely imposed electric field, the effective electric conductivity correct in the order of ζ ¯ can be obtained as
Λ = Λ + e π b k T m = 1 M z m n m D m   0 2 π { F m 0 ρ F m 0 ρ ) + ζ ¯   [ ρ ( z m e k T ψ eq 1 F m 0 ρ F m 1 ρ )
+ F m 1 z m e k T ψ eq 1 ( F m 0 z m e F ψ ) ] } ρ = b d ϕ ,
Where Λ = e 2 m = 1 M z m 2 n m D m / k T is the electric conductivity of the ionic fluid in the absence of the cylinders.
By substituting Equations (15), (20), and (21) together with F m 1 ( ρ ) into Equation (28) and making relevant calculations, we can obtain the effective electric conductivity in the fibrous medium as
Λ = Λ ( H + H 2 I β ζ ¯ )
for the Levine–Neale model, and
Λ = Λ ( H + I β ζ ¯ )
for the Zharkikh–Shilov model. In these expressions,
β = m = 1 M z m 3 n m D m m = 1 M z m 2 n m D m ,
H = 1 ϕ 1 + ϕ ,
I =   e k T [ φ 1 / 2 a ψ eq 1 ( b ) F ψ ( b ) + 2 φ χ ψ eq 1 ( b ) φ χ 2   a   b ( a 2 ρ 2 ρ 2 a 2 ) d ψ eq 1 d ρ d ρ ] ,
Note that the dimensionless coefficient I is irrelevant to the hydrodynamic boundary condition at the outer edge of the unit cell since the effect of electroosmosis on Λ does not appear to the order of ζ ¯ . When Equation (11) for the Zharkikh-Shilov electrostatic model is used, the result for the coefficient I [given by Equation (32) with χ = 1 + ϕ ] will not depend on the charge regulation coefficient L .

3. Results and Discussion

3.1. Charge Regulation Characteristics

Before calculating the electroosmotic velocity and electric conductivity in the array of charge-regulating circular cylinders from Equations (26) and (29)–(32), we present typical results of the dimensionless equilibrium surface charge density e σ ( eq ) / ε κ k T , zeta potential ζ ¯ , and charge regulation parameter L a as functions the bulk concentration n 1 of the charge-determining ions, equilibrium constants K + and K of the surface reactions in Equations (A1) and (A2), respectively, and porosity 1 ϕ of the fibrous medium in Figure 2, Figure 3 and Figure 4 using Equation (16), (A3), and (A6) (independent of the electrokinetic conditions at the outer boundary of the unit cell). The aqueous 1-1 electrolyte solution at T = 298   K (with dielectric constant 78.54) and the cylinders of radius a = 0.1   μ m with amphoteric functional group density N = 1 × 10 16   site / m 2 are chosen. When the bulk concentration of the charge-determining ions is smaller than that of the isoelectric point [ n 1 < ( K + K ) 1 / 2 ], the values of σ ( eq ) and ζ are negative as expected; the magnitude of ζ ¯ increases monotonically with a decrease in n 1 keeping the other parameters unchanged, whereas the magnitude of e σ ( eq ) / ε κ k T also increases monotonically with a decrease in n 1 for the case of ϕ = 0 but may have a maximal magnitude at a finite value of n 1 for the case of nonzero ϕ . When n 1 > ( K + K ) 1 / 2 , the values of e σ ( eq ) / ε κ k T and ζ ¯ are positive and may not be monotonic functions of n 1 .
The magnitudes of e σ ( eq ) / ε κ k T and ζ ¯ decrease as K / K + decreases (implying that the concentration of unionized surface groups increases). The cylinders are less positively charged or more negatively charged as K + K increases (implying that the concentration of negatively charged surface groups increases or that of positively charged surface groups decreases). The magnitude of e σ ( eq ) / ε κ k T decreases but that of ζ ¯ increases as ϕ increases. The charge regulation parameter L a is not a monotonic function of n 1 , K + K , K / K + , and ϕ . In general, the dimensionless quantities e σ ( eq ) / ε κ k T , ζ ¯ , and L a vary sensitively with changes in n 1 , K + , and K , but their dependencies on ϕ are negligible when n 1 is greater than about 0.1 mM.

3.2. Electroosmotic Velocity

Since the electroosmotic velocity of an arbitrary electrolyte solution in the charge-regulating fiber matrix predicted from the Zharkikh-Shilov electrostatic model is independent of the charge regulation coefficient L , we plot the normalized electroosmotic velocity η U / ε ζ E calculated from Equation (26) with χ = 1 ϕ (the Levine-Neale electrostatic model) for various values of the electrokinetic radius κ a , charge regulation parameter L a , and porosity 1 ϕ in Figure 5 and Figure 6 for the Happel and Kuwabara hydrodynamic models, respectively. The result is shown up to ϕ = 0.9 , corresponding to the minimum attainable porosity for an array of identical parallel circular cylinders [31]. The velocity η U / ε ζ E increases moderately and monotonically with a decrease in L a for fixed values of ϕ and κ a , but becomes independent of L a in the limiting cases of ϕ = 0 , κ a = 0 , and κ a . For specified values of L a and ϕ , this normalized velocity increases monotonically with an increase in κ a . When κ a = 0 , U = ε ζ E / 2 η as ϕ = 0 and U = 0 for any finite value of ϕ .
For the Kuwabara model, η U / ε ζ E decreases monotonically and significantly with an increase in ϕ for given values of L a and κ a and always equals unity as κ a . For the Happel model, however, η U / ε ζ E is no longer a monotonic function of ϕ but has a maximum (whose value can even be larger than unity) at some moderate value of ϕ for a constant value of κ a larger than about 0.7. The location of this maximum shifts to smaller ϕ as κ a decreases. For a fixed value of κ a smaller than about 0.6, similarly to the Kuwabara model, η U / ε ζ E becomes a monotonic decreasing function of ϕ . The Happel model always predicts a greater value of η U / ε ζ E than the Kuwabara model does for any combination of κ a , L a , and ϕ . Our results in the limit L a = 0 (solid curves in Figure 5 and Figure 6 for the case of constant surface charge density) agree with those obtained previously [18]. Note that, although the normalized electroosmotic velocity η U / ε ζ E just varies moderately (for the Levine-Neale model) or does not vary at all (for the Zharkikh-Shilov model) with the change in the charge regulation parameter L a , the dimensional electroosmotic velocity U still varies substantially with the change in the charge regulation characteristics such as n 1 , K + , and K since the zeta potential ζ varies strongly with these characteristics as shown in Figure 3.

3.3. Electric Conductivity

The dimensionless coefficient I in Equation (29) for the effective electric conductivity in the charge-regulating fibrous medium is independent of the hydrodynamic boundary condition prescribed at the virtual surface of the unit cell and that predicted from the Zharkikh-Shilov electrostatic model [given by Equation (32) with χ = 1 + ϕ ] is independent of the charge regulation coefficient L . In Figure 7, the coefficient H 2 I in Equation (29) calculated from Equations (31) and (32) with χ = 1 ϕ (the Levine-Neale electrostatic model) is plotted for various values of the electrokinetic radius κ a , charge regulation parameter L a , and porosity 1 ϕ . For the case of L a = 0 (constant surface charge density), H 2 I is always negative and thus the presence of the surface charge decreases and increases the effective electric conductivity for any porosity of the fiber matrix if β ζ > 0 and β ζ < 0 , respectively. The coefficient H 2 I decreases with a decrease in L a for specified values of κ a and ϕ , but this dependence becomes weak as κ a is large or ϕ is small. Interestingly, for the combination of small κ a , large ϕ , and large L a , the value of H 2 I may become positive and thus the presence of the surface charge increases and decreases the effective conductivity if β ζ > 0 and β ζ < 0 , respectively.
The coefficient H 2 I has a minimum at a finite value of ϕ and vanishes as ϕ = 0 for fixed values of κ a and L a , and a maximum of H 2 I may also exist at a larger value of ϕ when L a is large. The locations of these extrema shift to smaller ϕ as κ a decreases or L a increases. The coefficient H 2 I increases with an increase in κ a if either ϕ or L a is small, is not a monotonic function of κ a if both ϕ and L a are large, and vanishes as κ a irrespective of ϕ and L a . Our results in the limit L a = 0 (solid curves in Figure 7) agree with those obtained previously [18]. Although H 2 I may vary moderately (for the Levine-Neale model as κ a is large or ϕ is small) or does not vary at all (for the Zharkikh-Shilov model) with the change in the charge regulation parameter L a , the effective electric conductivity Λ in Equation (26) still varies substantially with the change in the charge regulation characteristics such as n 1 , K + , and K since ζ varies strongly with these characteristics.

4. Concluding Remarks

In this work, the electroosmosis and electric conduction of an arbitrary electrolyte solution in a fiber matrix constructed by a homogeneous array of parallel charge-regulating circular cylinders with an arbitrary electric double layer thickness were analytically investigated. A linearized charge regulation model was employed for the association and dissociation reactions occurring at the ionogenic functional groups on the surfaces of the cylinders. By solving the linearized electrokinetic equations applicable to a unit cell, we have determined the electroosmotic mobility of the fluid and effective electric conductivity in the fibrous medium. Our results show that the charge regulation characteristics such as the equilibrium constants of the association and dissociation reactions occurring at the cylinder surfaces and the bulk concentration of the potential-determining ions, and the porosity of the fiber matrix significantly, affect the surface charge density and potential, electroosmotic mobility, and effective electric conductivity.
Equations (26) and (29) were obtained on the basis of a small electric potential around the charge-regulating cylinder in a unit cell. Similar formulas for the electrophoretic mobility and effective electric conductivity of suspensions of charged spherical particles with a low zeta potential were shown to give good approximations for the situation of moderate zeta potentials (with errors less than 5% for ζ ¯ 2 in KCl solutions) [30,32]. Thus, our results might be used tentatively for in cases of reasonably high zeta potentials.
When the Levine–Neale electrokinetic cell model with Equations (10a) and (10b) is used, Equation (26) for the electroosmotic mobility and Equation (32) for the effective electric conductivity of an arbitrary electrolyte solution in the charge-regulating fibrous medium will be functions of the charge regulation parameter L a , and these transport properties will not depend on L a when the Zharkikh–Shilov electrokinetic cell model (Equation (11)) is used. Since both the electroosmotic flow and the electric conduction are considered within the linear approach in terms of the external electric field, and the electroosmotic mobility and effective electric conductivity are expressed by the only surface parameter ζ , in which the charge regulation characteristics such as n 1 , K + , and K are embedded, this dependency might be the same for any type of adsorption equilibrium. Accordingly, the predictions containing implicit dependencies on the charge regulation parameter L a by the Levine–Neale model seem to be incorrect. This outcome probably comes from the insufficiency of this model that the tangential (angular) component of the external electric field at the virtual boundary ( ρ = b ) of the unit cell is not specified. Indeed, in a cell-model analysis of the electrophoresis in a suspension of spherical particles, the Zharkikh–Shilov model predicts the electrophoretic mobility consistently with the ensemble-averaged result (obtained by using the concept of statistical mechanics), whereas the Levine–Neale model does not [33].

Author Contributions

Conceptualization, H.J.K.; methodology, H.J.K. and W.L.C.; investigation, H.J.K. and W.L.C.; writing—original draft preparation, H.J.K.; writing—review and editing, H.J.K.; supervision, H.J.K.; funding acquisition, H.J.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Science and Technology, Taiwan (Republic of China), grant number MOST 106-2221-E-002-167-MY3.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

A general model for the charge-regulating surface in contact with an electrolyte solution to result in fixed charge includes the following association and dissociation reactions occurring at the amphoteric functional groups AB distributed uniformly at the surface:
AB 2 Z + AB + B Z + ,
AB A Z + B Z + ,
where B Z + denotes the potential- or charge-determining ion with the valence Z (e.g., H + with Z = 1 ). By utilizing the Boltzmann distribution for the equilibrium concentration of B Z + , one can express the surface charge density as
σ = Z e N [ AB 2 Z + ] [ A Z ] [ AB ] + [ AB 2 Z + ] + [ A Z ] = Z e N Δ sinh { [ Z e ( ψ N ψ s ) + δ μ s ] / k T } 1 + Δ cosh { [ Z e ( ψ N ψ s ) + δ μ s ] / k T } ,
where
Δ = 2 ( K K + ) 1 / 2 ,
ψ N = k T Z e ln n 1 ( K + K ) 1 / 2 ,
N is the number of surface groups AB per unit area, ψ s is the surface potential, δ μ s is the perturbed electrochemical potential energy (defined by Equation (5)) of B Z + at the surface from equilibrium, n 1 is the bulk concentration of B Z + , K + and K are the equilibrium constants of the surface reactions (A1) and (A2), respectively, and ψ N is the Nernst potential, which vanishes at the point of zero charge (with n 1 = ( K + K ) 1 / 2 ). Equation (A3) leads to the charge regulation coefficient:
L = 1 ε ( d σ d ψ s ) ψ s = ζ = Z 2 e 2 N Δ { Δ + cosh [ Z e ( ψ N ( eq ) ζ ) / k T ] } ε k T { 1 + Δ cosh [ Z e ( ψ N ( eq ) ζ ) / k T ] } 2 ,
whose reciprocal is a length characterizing the charge regulation at the surface, and ζ is the value of ψ s at equilibrium.

References

  1. Rice, C.L.; Whitehead, R. Electrokinetic Flow in a Narrow Cylindrical Capillary. J. Phys. Chem. 1965, 69, 4017–4024. [Google Scholar] [CrossRef]
  2. Keh, H.J.; Liu, Y.C. Electrokinetic Flow in a Circular Capillary with a Surface Charge Layer. J. Colloid Interface Sci. 1995, 172, 222–229. [Google Scholar] [CrossRef]
  3. Keh, H.J.; Tseng, H.C. Transient Electrokinetic Flow in Fine Capillaries. J. Colloid Interface Sci. 2001, 242, 450–459. [Google Scholar] [CrossRef] [Green Version]
  4. Chang, C.C.; Wang, C.-Y. Electro-osmotic flow in a sector microchannel. Phys. Fluids 2009, 21, 42002. [Google Scholar] [CrossRef]
  5. Luo, R.H.; Keh, H.J. Electrokinetic flow and electric conduction of salt-free solutions in a capillary. Electrophor. 2020, 41, 1503–1508. [Google Scholar] [CrossRef] [PubMed]
  6. Levine, S.; Neale, G.; Epstein, N. The prediction of electrokinetic phenomena within multiparticle systems. J. Colloid Interface Sci. 1976, 57, 424–437. [Google Scholar] [CrossRef]
  7. Zharkikh, N.I.; Shilov, V.N. Theory of collective electrophoresis of spherical particles in the Henry approximation. Colloid J. USSR Engl. Transl. 1982, 43, 865–870. [Google Scholar]
  8. Kozak, M.W.; Davis, E. Electrokinetics of concentrated suspensions and porous media. J. Colloid Interface Sci. 1989, 127, 497–510. [Google Scholar] [CrossRef]
  9. Ohshima, H. Electrical Conductivity of a Concentrated Suspension of Spherical Colloidal Particles. J. Colloid Interface Sci. 1999, 212, 443–448. [Google Scholar] [CrossRef]
  10. Ding, J.M.; Keh, H.J. The Electrophoretic Mobility and Electric Conductivity of a Concentrated Suspension of Colloidal Spheres with Arbitrary Double-Layer Thickness. J. Colloid Interface Sci. 2001, 236, 180–193. [Google Scholar] [CrossRef] [Green Version]
  11. Carrique, F.; Cuquejo, J.; Arroyo, F.J.; Jimenez, M.L.; Delgado, A.V. Influence of cell-model boundary conditions on the con-ductivity and electrophoretic mobility of concentrated suspensions. Adv. Colloid Interface Sci. 2005, 118, 43–50. [Google Scholar] [CrossRef]
  12. Zholkovskij, E.K.; Masliyah, J.H.; Shilov, V.N.; Bhattacharjee, S. Electrokinetic Phenomena in concentrated disperse systems: General problem formulation and Spherical Cell Approach. Adv. Colloid Interface Sci. 2007, 134-135, 279–321. [Google Scholar] [CrossRef]
  13. Liu, H.C.; Keh, H.J. Electrophoresis and electric conduction in a suspension of charged soft particles. Colloid Polym. Sci. 2016, 294, 1129–1141. [Google Scholar] [CrossRef]
  14. Lin, W.C.; Keh, H.J. Diffusiophoresis in Suspensions of Charged Soft Particles. Colloids Interfaces 2020, 4, 30. [Google Scholar] [CrossRef]
  15. Kozak, M.W.; Davis, E. Electrokinetic phenomena in fibrous porous media. J. Colloid Interface Sci. 1986, 112, 403–411. [Google Scholar] [CrossRef]
  16. Ohshima, H. Electroosmotic Velocity in Fibrous Porous Media. J. Colloid Interface Sci. 1999, 210, 397–399. [Google Scholar] [CrossRef] [PubMed]
  17. Hsu, W.T.; Keh, H.J. Electric conductivity in a fibrous porous medium with thin but polarized double layers. Colloid Polym. Sci. 2003, 282, 985–992. [Google Scholar] [CrossRef]
  18. Keh, H.J.; Wu, Y.Y. Electroosmotic Velocity and Electric Conductivity in a Fibrous Porous Medium in the Transverse Direction. J. Phys. Chem. B 2011, 115, 9168–9178. [Google Scholar] [CrossRef]
  19. Chiang, C.C.; Keh, H.J. Transient electroosmosis in the transverse direction of a fibrous porous medium. Colloids Surfaces A: Physicochem. Eng. Asp. 2015, 481, 577–582. [Google Scholar] [CrossRef]
  20. Su, Y.W.; Keh, H.J. Electrokinetic Flow of Salt-Free Solutions in a Fibrous Porous Medium. J. Phys. Chem. B 2019, 123, 9724–9730. [Google Scholar] [CrossRef] [PubMed]
  21. Ninham, B.W.; Parsegian, V. Electrostatic potential between surfaces bearing ionizable groups in ionic equilibrium with physiologic saline solution. J. Theor. Biol. 1971, 31, 405–428. [Google Scholar] [CrossRef]
  22. Carnie, S.L.; Chan, D.Y. Interaction Free Energy between Plates with Charge Regulation: A Linearized Model. J. Colloid Interface Sci. 1993, 161, 260–264. [Google Scholar] [CrossRef]
  23. Ding, J.M.; Keh, H.-J. Electrophoretic Mobility and Electric Conductivity in Dilute Suspensions of Charge-Regulating Composite Spheres. Langmuir 2003, 19, 7226–7239. [Google Scholar] [CrossRef]
  24. Hill, R.J. On the electrophoretic mobility of succinoglycan modelled as a spherical polyelectrolyte: From Hermans-Fujita theory to charge regulation in multi-component electrolytes. J. Colloid Interface Sci. 2016, 482, 131–134. [Google Scholar] [CrossRef] [Green Version]
  25. Philipse, A.; Tuinier, R.; Kuipers, B.; Vrij, A.; Vis, M. On the Repulsive Interaction Between Strongly Overlapping Double Layers of Charge-regulated Surfaces. Colloid Interface Sci. Commun. 2017, 21, 10–14. [Google Scholar] [CrossRef] [Green Version]
  26. Lin, C.Y.; Keh, H.J. Sedimentation Velocity and Potential in Dilute Suspensions of Charge-Regulating Porous Spheres. J. Phys. Chem. B 2019, 123, 3002–3009. [Google Scholar] [CrossRef]
  27. Bharti; Gopmandal, P.P.; Sinha, R.K.; Ohshima, H. Electrophoresis of pH-regulated zwitterionic soft particle: A semi-analytical study. Colloid Polym. Sci. 2020, 298, 79–89. [Google Scholar] [CrossRef]
  28. Keh, H.J.; Ding, J.M. Electrokinetic flow in a capillary with a charge-regulating surface polymer layer. J. Colloid Interface Sci. 2003, 263, 645–660. [Google Scholar] [CrossRef]
  29. Sadeghi, M.; Saidi, M.H.; Moosavi, A.; Sadeghi, A. Geometry effect on electrokinetic flow and ionic conductance in pH-regulated nanochannels. Phys. Fluids 2017, 29, 122006. [Google Scholar] [CrossRef]
  30. Keh, H.-J.; Ding, J.M. Electrophoretic Mobility and Electric Conductivity of Suspensions of Charge-Regulating Colloidal Spheres. Langmuir 2002, 18, 4572–4583. [Google Scholar] [CrossRef]
  31. Berryman, J.G. Random close packing of hard spheres and disks. Phys. Rev. A 1983, 27, 1053–1061. [Google Scholar] [CrossRef]
  32. O’Brien, R. The electrical conductivity of a dilute suspension of charged particles. J. Colloid Interface Sci. 1981, 81, 234–248. [Google Scholar] [CrossRef]
  33. Wei, Y.K.; Keh, H.-J. Diffusiophoresis and Electrophoresis in Concentrated Suspensions of Charged Colloidal Spheres. Langmuir 2001, 17, 1437–1447. [Google Scholar] [CrossRef]
Figure 1. Geometrical sketch of the unit cell model for an array of long parallel circular cylinders.
Figure 1. Geometrical sketch of the unit cell model for an array of long parallel circular cylinders.
Colloids 05 00019 g001
Figure 2. The dimensionless equilibrium surface charge density e σ ( eq ) / ε κ k T versus the bulk concentration n 1 of an aqueous 1-1 electrolyte solution in a fibrous medium with a = 100 nm and N = 1 × 10 16 site/m2: (a) K + K = 10 6   M 2 ; (b) K / K + = 10 2 . The solid and dashed curves represent the cases ϕ = 0.3 and ϕ = 0 , respectively.
Figure 2. The dimensionless equilibrium surface charge density e σ ( eq ) / ε κ k T versus the bulk concentration n 1 of an aqueous 1-1 electrolyte solution in a fibrous medium with a = 100 nm and N = 1 × 10 16 site/m2: (a) K + K = 10 6   M 2 ; (b) K / K + = 10 2 . The solid and dashed curves represent the cases ϕ = 0.3 and ϕ = 0 , respectively.
Colloids 05 00019 g002
Figure 3. The dimensionless zeta potential ζ ¯ versus the bulk concentration n 1 of an aqueous 1-1 electrolyte solution in a fibrous medium with a = 100 nm and N = 1 × 10 16 site/m2: (a) K + K = 10 6   M 2 ; (b) K / K + = 10 2 . The solid and dashed curves represent the cases ϕ = 0.3 and ϕ = 0 , respectively.
Figure 3. The dimensionless zeta potential ζ ¯ versus the bulk concentration n 1 of an aqueous 1-1 electrolyte solution in a fibrous medium with a = 100 nm and N = 1 × 10 16 site/m2: (a) K + K = 10 6   M 2 ; (b) K / K + = 10 2 . The solid and dashed curves represent the cases ϕ = 0.3 and ϕ = 0 , respectively.
Colloids 05 00019 g003
Figure 4. The dimensionless charge regulation parameter L a versus the bulk concentration n 1 of an aqueous 1-1 electrolyte solution in a fibrous medium with a = 100 nm and N = 1 × 10 16 site/m2: (a) K + K = 10 6   M 2 ; (b) K / K + = 10 2 . The solid and dashed curves represent the cases ϕ = 0.3 and ϕ = 0 , respectively.
Figure 4. The dimensionless charge regulation parameter L a versus the bulk concentration n 1 of an aqueous 1-1 electrolyte solution in a fibrous medium with a = 100 nm and N = 1 × 10 16 site/m2: (a) K + K = 10 6   M 2 ; (b) K / K + = 10 2 . The solid and dashed curves represent the cases ϕ = 0.3 and ϕ = 0 , respectively.
Colloids 05 00019 g004
Figure 5. The normalized electroosmotic velocity η U / ε ζ E of an electrolyte solution through a fibrous medium calculated from Equation (26a) for the Happel/Levine-Neale cell model versus the parameters κ a (a) and ϕ (b). The solid, dot-dash, and dashed curves represent the cases of the charge regulation parameter L a equal to 0, 5, and , respectively.
Figure 5. The normalized electroosmotic velocity η U / ε ζ E of an electrolyte solution through a fibrous medium calculated from Equation (26a) for the Happel/Levine-Neale cell model versus the parameters κ a (a) and ϕ (b). The solid, dot-dash, and dashed curves represent the cases of the charge regulation parameter L a equal to 0, 5, and , respectively.
Colloids 05 00019 g005
Figure 6. The normalized electroosmotic velocity η U / ε ζ E of an electrolyte solution through a fibrous medium calculated from Equation (26b) for the Kuwabara/Levine-Neale cell model versus the parameters κ a (a) and ϕ (b). The solid, dot-dash, and dashed curves represent the cases of the charge regulation parameter L a equal to 0, 5, and , respectively.
Figure 6. The normalized electroosmotic velocity η U / ε ζ E of an electrolyte solution through a fibrous medium calculated from Equation (26b) for the Kuwabara/Levine-Neale cell model versus the parameters κ a (a) and ϕ (b). The solid, dot-dash, and dashed curves represent the cases of the charge regulation parameter L a equal to 0, 5, and , respectively.
Colloids 05 00019 g006
Figure 7. The dimensionless coefficient H 2 I in Equation (29a) for the Levine-Neale cell model for the effective electric conductivity in a fibrous medium calculated from Equations (31) and (32) versus the parameters κ a (a) and ϕ (b). The solid, dot-dash, and dashed curves represent the cases of the charge regulation parameter L a equal to 0, 5, and , respectively.
Figure 7. The dimensionless coefficient H 2 I in Equation (29a) for the Levine-Neale cell model for the effective electric conductivity in a fibrous medium calculated from Equations (31) and (32) versus the parameters κ a (a) and ϕ (b). The solid, dot-dash, and dashed curves represent the cases of the charge regulation parameter L a equal to 0, 5, and , respectively.
Colloids 05 00019 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, W.L.; Keh, H.J. Electroosmosis and Electric Conduction of Electrolyte Solutions in Charge-Regulating Fibrous Media. Colloids Interfaces 2021, 5, 19. https://0-doi-org.brum.beds.ac.uk/10.3390/colloids5010019

AMA Style

Chen WL, Keh HJ. Electroosmosis and Electric Conduction of Electrolyte Solutions in Charge-Regulating Fibrous Media. Colloids and Interfaces. 2021; 5(1):19. https://0-doi-org.brum.beds.ac.uk/10.3390/colloids5010019

Chicago/Turabian Style

Chen, Wei L., and Huan J. Keh. 2021. "Electroosmosis and Electric Conduction of Electrolyte Solutions in Charge-Regulating Fibrous Media" Colloids and Interfaces 5, no. 1: 19. https://0-doi-org.brum.beds.ac.uk/10.3390/colloids5010019

Article Metrics

Back to TopTop