Next Article in Journal
Ionic Liquid in Phosphoric Acid-Doped Polybenzimidazole (PA-PBI) as Electrolyte Membranes for PEM Fuel Cells: A Review
Next Article in Special Issue
Geometrical Influence on Particle Transport in Cross-Flow Ultrafiltration: Cylindrical and Flat Sheet Membranes
Previous Article in Journal
Normally-Off p-GaN Gated AlGaN/GaN MIS-HEMTs with ALD-Grown Al2O3/AlN Composite Gate Insulator
Previous Article in Special Issue
Fouling Mitigation via Chaotic Advection in a Flat Membrane Module with a Patterned Surface
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Numerical Procedure to Estimate the Electric Charge in the Pore from Filtration of Single-Salt Solutions

1
Institut de Science des Matériaux de Mulhouse (IS2M-UMR CNRS 7361), Université de Haute-Alsace, 3Bis Rue Alfred Werner, CEDEX, 68098 Mulhouse, France
2
Institut UTINAM (UMR CNRS 6213), Université de Bourgogne Franche-Comté, 16 route de Gray, CEDEX, 25030 Besançon, France
*
Author to whom correspondence should be addressed.
Membranes 2021, 11(10), 726; https://doi.org/10.3390/membranes11100726
Submission received: 10 September 2021 / Accepted: 19 September 2021 / Published: 23 September 2021
(This article belongs to the Special Issue Numerical Modelling in Membrane Processes)

Abstract

:
The assessment of physicochemical parameters governing the transport of ions through nanoporous membranes is a major challenge due to the difficulty in experimental estimation of the dielectric constant of the solution confined in nanopores and the volumetric membrane charge. Numerical identification by adjusting their values to fit experimental data is a potential solution, but this method is complicated for single-salt solutions due to the infinite number of couples that can describe a rejection curve. In this study, a novel procedure based on physical simplifications which allows the estimation of a range of values for these two parameters is proposed. It is shown here that the evolution of the interval of membrane charge with salt concentration can be described in all the experimental conditions by the Langmuir–Freundlich hybrid adsorption isotherm. Finally, it is highlighted that considering the mean dielectric constant and the adsorption isotherms assessed from a range of concentrations allowed a good prediction of rejection curves, irrespective of the salt and membrane considered.

1. Introduction

Pressure-driven membrane processes are often used to solve environmental issues such as desalination, wastewater treatment, drinking water production, etc. Among the potential applications, the removal of ionic contaminants is one of the main purposes for which membrane processes are particularly competitive due to the charged surface of commercial membranes. For decades, it has been well-know that the performances of nanoporous membranes (i.e., nanofiltration and low-molecular-weight cut-off ultrafiltration membranes) are governed by several physical complex mechanisms, and many research teams are still working on the development of accurate models which can describe and even predict the filtration performances of these membranes [1,2,3]. Many models are available in the literature, but a classical approach used by many researchers consists in considering an interfacial partitioning based on the Donnan equilibrium and steric exclusion, and describing the transport inside pores by the extended Nernst–Planck equation [4,5,6]. The vision of the transfer at the membrane/solution interfaces is usually improved by considering the dielectric exclusion due to a solvation energy barrier induced by changes in the dielectric constant of the solvent inside pores. This model, which was initially proposed by Bowen et al. [7], has proved a good ability to describe the experimental ion rejections obtained in many experimental conditions [8,9,10,11], provided that four parameters are known or adjusted. Two of these parameters are structural properties, whereas the two others are physicochemical properties of the membrane and the solution confined in nanopores. The two membrane structural parameters, namely intrinsic hydraulic permeability Lp and mean pore radius rp, are easy to assess and are usually estimated from pure water flux and filtration of a neutral solute, respectively [12]. The estimation of the volumetric membrane charge Xd and the dielectric constant of the solution inside pores εp is more problematic and the way to assess them represents the key for a better understanding of transport mechanisms and the development of a predictive numerical tool. Indeed, these parameters are very complicated to assess experimentally, even if recent studies have shown that it is not completely unrealistic to develop original techniques to do so. For instance, streaming potential or current measurements [13,14,15] have demonstrated a high potential to determine membrane charge. However, only the estimation of the charge at the membrane surface is possible for nanoporous membranes (tangential measurements), and the charge at the pore wall remains undeterminable. Moreover, the presence of conduction and streaming current potentially occurring in membrane porosity also complicate the signal analysis. Electrochemical impedance spectroscopy [16,17,18] or membrane potential [19,20] are potential techniques to assess the dielectric constant of the solution inside pores, but the difficulty in separating the contribution of the solution to the signal of the wet membrane is often problematic. In any case, these methods are not completely satisfying at the present time and a numerical way is often required [21].
It has already been highlighted that there exists an infinite number of couples (εp, Xd) which allow the description of a salt rejection curve [22]. εp and Xd, since they have different impacts on ions, can have their values assessed from ion mixtures, for which only one couple can describe the various curves and be extrapolated to single-salt solutions [23]. However, this procedure, which has been presented and validated during previous studies, presents the inconvenience of requiring many experiments at various proportions before being extrapolated for single salts.
The simpler approach proposed in this paper thus consists of squeezing the values by physical deduction to reduce the interval of possible values that correctly describe single-salt rejections. In this paper, this procedure is implemented to determine the interval of membrane charge which allows a correct description of experimental curves. These intervals of charge are then computed to estimate the best adsorption isotherms and investigate the predictive ability of the proposed procedure in several experimental conditions of salts and membranes.

2. Materials and Methods

Experimental filtrations were performed with a laboratory pilot plant provided by TIA (Techniques Industrielles Appliquées, Bollène, France). Three ceramic membranes were used in this work, namely a TiO2 membrane (1 kDa) provided by TAMI Industries (Nyons, France) and two home-made zeolite membranes (Na-mordenite and a bi-layer MFI/MFI). The membranes (25 cm length with an inner diameter of 7 mm) consist of a tubular support of porous alumina, on which the active layer (titania or zeolite) was deposited. The specific preparation of zeolite membranes and their characteristics are detailed in previous papers [11,24,25].
The membrane module is mounted on the lab scale setup in stainless steel. The solution is introduced into the feed tank (5 L) and the fluid flow inside the membrane module is provided by a volumetric pump (piston, 2.2 kW with frequency variation, Pmax = 3 MPa) with a flow rate of 700 L/h (i.e., fluid velocity > 5 m/s). The applied pressure ΔP is adjusted from 4 to 12 bar by closing a manual valve. Experimental tests were conducted at a constant temperature of 25 °C (controlled by a cooling unit) and with a constant feed concentration by recycling both feed and retentate solutions.
First of all, the membrane is conditioned with pure water by filtrating demineralized water (conductivity < 0.1 µS/cm) to reach steady hydraulic performances. The filtration of an uncharged solute (vitamin B12 from Alfa Aesar, purity 98%, 9 × 10−4 mol m−3) was performed to investigate steric effect and estimate the mean pore radius of the membrane. Retentate and permeate streams were sampled and analysed by UV–visible spectrometry at 362 nm with a Lambda 35 spectrophotometer from Elmer Instrument, Waltham, USA.
Filtrations of single-salt solutions containing NaF, NaCl, NaI and Na2SO4 (Sigma-Aldrich, St. Quentin Fallavier, France, purity > 99%) were performed at various concentrations from 0.5 to 100 mol m−3. For each operating condition, concentrations of retentate and permeate solutions (Cr and Cp, respectively) were measured by conductimetry (conductimeter PC 5000 L Phenomenal, VWR, Radnor, US). The rejection rate R is calculated as follows:
R = 1 C p C r
For each applied pressure, the permeation flux (Jv) was measured by weighing the permeate stream during a given time.
Between each filtration, the setup was washed, and pure water was filtered to estimate the intrinsic hydraulic membrane permeability Lp from the slope of the linear evolution of water flux with applied pressure Jw = f(ΔP).
J w = L p μ Δ P  
With μ as the dynamic viscosity of the fluid.

3. Numerical Modelling

3.1. Transport Model

The porous medium is assumed to be one-dimensional, constituted of cylindrical and uniform pores. The flow is assumed to be laminar without temporal variation (steady state). Variations according to radial direction (potential, concentrations) are neglected and constant values, defined as the radial average, are considered in this model.
The transport of ionic solutes in the porous media is modelled from the fact that the molar flux satisfies the continuity equation and divergence of the flux d i v ( J i ) is equal to zero. The transport of each solute (Ji) in the pores is considered as the sum of three contributions induced by diffusion due to the gradient of chemical potential ( μ i ) , electro-migration due to the electrical gradient Ψ and convection [4]. In steady state, this equation is:
d i v ( J i ) = d i v [ D i R T ( μ i + F z i Ψ ) + c i V ] = 0
where Di, zi and ci are the diffusivity, the valence and the local concentration inside pores of the solute i, respectively. R is the ideal gas constant, T is the temperature, F is the Faraday constant (96,487 C mol−1), Ψ is the gradient of electric potential in the pore and V is the fluid velocity in the pore.
The influence of pore shape and tortuosity on the diffusive and convective transports is taken into account by introducing hindered factors Ki,c and Ki,d [26]. Equation (3) can be rewritten and is known as the extended Nernst–Planck equation:
J i = K i , d D i d C i d x K i , d D i c i d l n γ i d x K i , d D i z i c i F R T d Ψ d x + K i , c c i V = V C i , p
where γi and Ci,p are the activity coefficient (calculated by using the extended Debye–Hückel equation) and the permeate concentration of the solute i, respectively.
Hindrance factors for convection and diffusion are calculated by equations proposed by Bowen et al. [27]:
K i , c = ( 2 ϕ ) ( 1.0 + 0.054 λ i 0.988 λ i 2 + 0.441 λ i 3 )
K i , d = 1.0 2.3 λ i + 1.154 λ i 2 + 0.224 λ i 3
where ϕ i = ( 1 λ i ) 2 is the steric partitioning coefficient, which depends on λi, which is defined as the ratio of the solute’s Stokes radius (ri,s) and pore radius (rp) λ i = r i , s r p [28].
The differential Equation (4) is solved for each solute i with two boundary conditions. The generalized chemical potentials of each solute at the interfaces between the bulk and the pore solution are equal at the pore inlet and outlet (both sides of the active layer x = 0, x = L). This equilibrium at the pore/solution interfaces can be described by the product of three contributions (steric, electric and dielectric), as proposed in Equation (5) [25].
The ratio between the concentration at the pore entry or outlet ci(x) and that of the feed or permeate solution Ci (depending on the interface considered) can be calculated by:
For   x = 0   or   L   c i ( x ) C i = ϕ i γ i , s γ i , p exp ( Δ W i ) exp ( z i F R T Δ Ψ D )
where Δ ψ D is the Donnan potential (difference in electrical potential at both sides of the interface, free solution/solution in the pore) and Δ W i is the solvation energy barrier due to the dielectric effect. This energy is computed in the model by considering an apparent dielectric permeability (εp) of solution in the pore in the Born model [29]:
Δ W i = ( z i e ) 8 π ε 0 r i , s k B T ( 1 ε p 1 ε b )
where e is the electronic charge (e = 1.602 10–19 C), kB is the Bolzmann constant (kB = 1.381 10−23 m2 kg s−2 K−1), ε0 is the permittivity of free space (ε0 = 8.85419 10−12 F m−1) and εb is the bulk dielectric constant (78.4 for water). It should be stressed that the impact of so-called “image forces” on dielectric exclusion, which is sometimes considered in Equation (7) [30,31], is not taken into account in this study.
For each ionic solute, Equations (4) and (7) are solved by respecting the electroneutrality condition in the feed solution (Equation (9)) and in the pore (Equation (10)):
i = 1 n z i C i = 0
for   0     x     L   i = 1 n z i c i ( x ) + X d = 0
where Xd is the volumetric membrane charge density in the pore.
The resolution of this set of equations has already been detailed in a previous paper [32].
The permeation flux is calculated from Equation (11) by taking into account the osmotic pressure difference ( Δ π ) due to different salt concentrations on both sides of the membrane.
J v = L p μ ( Δ P Δ π )
As the permeate concentration appears in the differential mass transfer equation (Equation (4)), the set of equations is iteratively solved by using an explicit Euler method up to convergence.
Four physical parameters are required (i.e., Lp, rp, Xd and εp) to implement simulations. First, the hydraulic permeability Lp was calculated from pure water filtration (Equation (2)). The mean pore radius is obtained from the filtration of a neutral solute (vitamin B12, ri,s = 0.72 nm and Di = 3.4 × 10−10 m2 s−1) by adjusting its value to fit experimental curve Ri = f(Jv). The mean pore radius is thus assessed by minimizing the quadratic error between experimental rejections and those calculated by a pore flow model based on the previously presented model, for which electro-migration, electric and dielectric contributions are neglected.
R i = 1 ϕ i K i , c 1 ( 1 ϕ i K i , c ) e x p ( K i , c r p 2 Δ P 8 μ K i , d D i )

3.2. Adsorption Models

After being assessed, the volumetric membrane charge densities are normalised by feed salt concentration Xd/Cf, and the evolution of this normalised charge is linked to concentration by various adsorption isotherms [33].
- Freundlich isotherm, which empirically describes multilayer adsorption without saturation:
X d C f = K F C f 1 n n
- Langmuir isotherm, which is a physical equation describing monolayer adsorption with surface saturation:
X d C f = Q m a x K L 1 + K L C f
- Langmuir–Freundlich isotherm (also called the Sips isotherm), which is a hybrid equation combining the two previous approaches:
X d C f = Q m a x K L F C f 1 n n 1 + K L F C f 1 n
In Equations (13)–(15), Cf represents the equilibrium feed concentration, 1/n the adsorption intensity and Qmax the maximum adsorption capacity. KL, KF and KLF are the Langmuir, Freundlich and Langmuir–Freundlich constants, which are all related to adsorption capacity.
It should be noted that the Langmuir–Freundlich model is usually more suitable for predicting adsorption on heterogeneous surfaces, avoiding the limitation of the two other isotherms. For this reason, such a hybrid model is more adequate for mineral membranes which exhibit various kind of porosities.
Hence, at a low adsorbate concentration this model can be equated to the Freundlich model, and at a high concentration, it can predict monolayer adsorption, similar to Langmuir model.

4. Results and Discussion

4.1. Preliminary Results

Before investigating electric and dielectric exclusion mechanisms, the intrinsic hydraulic permeability of the various studied membranes was assessed from water flux measurements. The hydraulic permeability was estimated between each experiment to ensure that the structural properties of the membranes do not vary during the experimental campaign. Additionally, filtration of the vitamin B12 solution was performed to estimate the sieving properties of the membrane via the calculation of the mean pore radius with Equation (11). The values obtained for each membrane are summarised in Table 1.
Structural properties given in Table 1 show that the three membranes are notably different in both terms of permeability and pore size, but with different behaviours. The TiO2 membrane appears to be very close to nanofiltration with a tight pore size, whereas the Na-mordenite exhibits the largest pore size, typical of an ultrafiltration membrane. The MFI–MFI membrane shows intermediate pore size close to low MWCO UF membranes. Concerning hydraulic permeability, the trend is completely different and the TiO2 membrane demonstrates the highest permeability compared to the two others. This behaviour is probably due to a thinner skin layer of the TiO2 membrane. The steric exclusion is thus clearly different between these three membranes.
It should be noted that the experimental error on permeability is very low, less than 4% of variation between all the experiments. The rp value is estimated by fitting a rejection curve with a pore flow model, so error is mainly due to model assumptions. The experimental variation in VB12 rejection is about 5%. The corresponding variation is less than 10% in pore radius except for the Na-mordenite membrane. However, steric effects are negligible with this membrane due to large pores, so uncertainties do not have any impact.

4.2. Numerical Procedure

Once these two parameters are known, the model only required two other parameters, namely the dielectric constant of the solution confined in nanopores (εp) and the volumetric membrane charge density (Xd), which can be adjusted to fit experimental curves. However, it has already been highlighted that there exists and infinite number of couples (εp, Xd) which allow the description of one rejection curve [22,34], as it is represented for two NaF solutions with different concentrations in Figure 1a,b.
Figure 1a shows experimental and simulated rejection curves obtained by filtration of two solutions containing 20 and 100 mol m−3 of NaF. The simulated curves were obtained by assessing the various couples (εp, Xd) and the corresponding couples are drawn in a Xd = f(εp) graph (Figure 1b). From such curves, it is impossible to obtain the relevant couple that describes the experimental rejection curve. For this reason, physical simplifications are required.
In the case of 1:1 monovalent salts, an increase in membrane charge, either in negative or positive values, leads to an increase in rejection, the curve depicting the possible couples is always symmetrical on the Xd = 0 axis.
First, zeta potential measurements were implemented on powders (synthesized in the same way than the active layers of the three membranes) during previous studies. The values do not have any importance for this study, but they can be found in references [24] for MFI, [25] for mordenite and [35] for TiO2. Indeed, only the sign of the membrane charge is required for physical simplification. In these references, it has been shown that the three membranes always exhibit a negative charge at natural pH, irrespective of the filtered solution considered. From this finding, only the negative part of the curve can be considered in Figure 1b.
The dielectric exclusion is induced by a change in the dielectric constant of the solution within membrane pores due to confinement. This decrease in the dielectric constant of the confined solution is mainly imputable to the orientation of solvent molecules due to the electric field induced by membrane charge and ions. The maximum dielectric exclusion thus necessarily corresponds to the highest concentration and the value of the dielectric constant estimated at the highest concentration is the lower limit, irrespective of the concentration of the solution. This deduction allows an estimation of the possible values for the dielectric constant of the solution. The dielectric constant cannot be lower than the minimum value obtained with the highest concentration. An example of the procedure is illustrated in Figure 1b. Considering the curves of potential (εp, Xd) obtained for 100 mol m−3 (curve in red), the minimum possible value of εp is 69, which means that the value cannot be lower irrespective of the concentration considered. Additionally, given the physical meaning of the dielectric constant, it is obvious that its value cannot be higher than that of unconfined water, i.e., 78.4.
It should be stressed that the minimum value of the εp range does not vary significantly between 75 (sometimes 50) and 100 mol m−3. Hence, there is no interest in increasing concentration above 100 mol m−3. For this reason, concentration was varied up to 100 mol m−3 and the minimum value of the dielectric constant is always identified at this concentration.
In this example, the apparent dielectric constant is thus in the range of 68 to 78.4 for all the NaF solutions filtered with the titania membrane. Using this range of εp, the corresponding values of volumetric membrane charge density estimated to fit the experimental curve obtained with the solution of NaF 20 mol m−3 are between −69 and −87 eq m−3 (surrounded by the blue line in Figure 1b).
Once the intervals of normalised membrane charge (Xd/Cf) are determined for all the concentrations, the latter are then linked to concentration through adsorption isotherms. In this study, three usual isotherms (Langmuir, Freundlich and Langmuir–Freundlich) have been investigated, but it was found that two-parameter models (Langmuir or Freundlich) do not allow a correct description in some cases. Oppositely, the hybrid three-parameter model always led to a correct description of membrane charge evolution. The parameters of the isotherms were assessed by minimizing the quadratic criteria between the predicted Xd/Cf-value and the minimum and maximum values of the Xd/Cf range.
An example of the description by the three adsorption models is provided in Figure 2. It can be seen that the best description is obtained by the Freundlich–Langmuir model and it was thus chosen to use only the Langmuir–Freundlich isotherm to describe the experimental trends presented in the next section.

4.3. Numerical Investigation

The procedure presented in the previous section was implemented with various salts containing halide and sulfate ions at various concentrations and for three ceramic membranes. The experimental rejection rates were numerically fitted to obtain the curves (Xd, εp) relative to each salt concentration and each membrane. By assuming the assumptions previously reported, the range of the electric charge in the pore is estimated for each filtration, leading to minimum and maximum values of Xd/Cf. Finally, the best-fitted Langmuir–Freundlich isotherm is determined, and two examples are provided in Figure 3 for the filtration of iodide sodium by TiO2 membrane and sodium fluoride by the Na-mordenite membrane.
From Figure 3a,b, it can be seen that the interval of Xd/Cf for each concentration can strongly vary depending on the salt and the membrane considered. For instance, in the case of iodide sodium by the TiO2 membrane, the proposed method is very precise, and the range of possible Xd/Cf is only a point. Oppositely, with fluoride and Na-mordenite, a larger interval of value was obtained. In all the cases, the experimental data are correctly described by the Langmuir–Freundlich isotherms and the possibility of using these isotherms instead of the adjustable parameter Xd for predictive purposes is investigated hereafter. It should be stressed that the dielectric constant considered for simulations was the mean value between the maximum and minimum values.
The fitted parameters of the Langmuir–Freundlich isotherms obtained with the various salts and membranes are given in Table 2.
The ranges of εp and Xd/Cf obtained with the various salts and membranes are summarised in Table 3 and examples are illustrated in Figure 4, Figure 5 and Figure 6.
A first example of experimental adsorption isotherm for Na2SO4 filtered by an MFI–MFI membrane is provided in Figure 4a and the simulations are compared with an experimental rejection curve obtained with 0.44 mol m−3 of Na2SO4 in Figure 4b.
It should be mentioned that the validation is carried out at concentrations used for the determination of adsorption isotherms. However, the parameters εp and Xd/Cf were not assessed at this concentration but by using adsorption isotherms estimated with the whole range of Xd/Cf (between min and max values) for all the experimental data and the mean value of the dielectric constant estimated from all the experimental curves. This allows for discussion of the reliability of the novel procedure based on an adsorption isotherm and a mean value of dielectric constants estimated from squeezed ranges of values.
From Figure 4, it is highlighted that the simulated curve is in good agreement with the experimental curve. It should be noted that the isotherm does not correctly describe the normalised membrane charge Xd/Cf at high concentration due to the possible positive value, which is not possible to describe with the adsorption isotherm.
Figure 5 shows another example of procedure relevance and it is highlighted from this figure that the performance obtained with the filtration of a 2 mol m−3 NaI solution is correctly predicted by the model and the Langmuir–Freundlich isotherm without adjusting parameters.
Finally, the quality of the procedure was studied on two concentrations of NaF using the MFI–MFI membrane. The adsorption isotherm correctly describes the evolution of Xd/Cf and its use in simulations allows a suitable description of the experimental curves obtained by filtering 2 and 50 mol m−3 NaF solutions, as can be seen in Figure 6.

5. Conclusions

It is well-known that the difficulty in modelling and predicting membrane filtration performances lies in the assessment of model parameters. The main objective of this study was to develop a numerical tool able to estimate the salt rejection from adsorption isotherms using a novel procedure based on physical simplifications. First, this novel numerical procedure was proposed to estimate a range of reasonable values for the two main parameters, namely the dielectric constant of the confined solution and the volumetric membrane charge. Various adsorption isotherms have been investigated in this study and it was found that only the hybrid Langmuir–Freundlich isotherm allowed a correct description of all the experimental curves for heterogenous mineral membranes. This procedure was applied to various sodium salts and three mineral membranes, and the results highlighted that the experimental performances are correctly predictable with the proposed approach without adjusting parameters on the specific experimental condition.

Author Contributions

Conceptualization, P.D. and T.J.D.; methodology, P.D. and S.D.; software, P.D. and S.D.; validation, P.D. and S.D.; formal analysis, P.D. and S.D.; investigation, P.D. and S.D.; resources, P.D.; data curation, P.D. and S.D.; writing—original draft preparation, P.D. and S.D.; project administration, P.D. and T.J.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available on request.

Conflicts of Interest

The authors declare no conflict of interest.

Glossary

ciConcentration of ion i within the pore (mol m−3)
Ci,pPermeate concentration of ion i (mol m−3)
Ci,fFeed concentration of ion i (mol m−3)
Di,Diffusion coefficient of ion i at infinite dilution (m2 s−1)
eElectronic charge (1.602 10−19 C)
FFaraday constant (96487 C mol−1)
jiFlux of ion i (mol m−2 s−1)
JvPermeation flux (m3 m−2 s−1)
JwPermeation flux of pure water (m3 m−2 s−1)
kBBoltzmann constant (1.381 10–23 m2 kg s−2 K−1)
Ki,cIonic hindrance factor for convection (dimensionless)
Ki,dIonic hindrance factor for diffusion (dimensionless)
LpHydraulic permeability (m3 m−2)
RGas constant (8.314 J mol−1 K−1)
ri,sStokes radius of ion i (m)
RiObserved rejection of ion i (dimensionless)
rpAverage pore radius (m)
TTemperature (K)
VSolvent velocity in the pore (m s−1)
xAxial position within the pore (m)
XdMembrane effective charge density in the pore (eq m−3)
ziValence of ion i (dimensionless)
Greek letters
γi,pActivity coefficient of ion i in the pore (dimensionless)
γi,sActivity coefficient of ion i in the solution side of the interface (dimensionless)
ΔPApplied pressure (Pa)
ΔWiDielectric exclusion energy (J)
ΔψDDonnan potential (V)
ΔπOsmotic pressure difference (Pa)
ε0Permittivity of free space (8.85419 × 10−12 F m−1)
εbBulk dielectric constant (dimensionless)
εpPore dielectric constant (dimensionless)
μDynamic viscosity (Pa s)
φiSteric partition coefficient (dimensionless)
ψElectrical potential within the pore (V)

References

  1. Fridman-Bishop, N.; Nir, O.; Lahav, O.; Freger, V. Predicting the Rejection of Major Seawater Ions by Spiral-Wound Nanofiltration Membranes. Environ. Sci. Technol. 2015, 49, 8631–8638. [Google Scholar] [CrossRef] [PubMed]
  2. Fernández de Labastida, M.; Yaroshchuk, A. Nanofiltration of multi-ion solutions: Quantitative control of concentration polarization and interpretation by solution-diffusion-electro-migration model. Membranes 2021, 11, 272. [Google Scholar] [CrossRef] [PubMed]
  3. Diaz, P.A.B.; de Araujo Kronemberger, F.; Habert, A.C. Predictive study of succinic acid transport via nanofiltration membranes using a Maxwell–Stefan approach. J. Chem. Technol. Biotechnol. 2021, 96, 785–800. [Google Scholar] [CrossRef]
  4. Tsuru, T.; Nakao, S.; Kimura, S. Calculation of ion rejection by extended Nernst-Planck equation with charged reverse osmosis membranes for single and mixed electrolyte solutions. J. Chem. Eng. Jpn. 1991, 24, 511–517. [Google Scholar] [CrossRef] [Green Version]
  5. Szymczyk, A.; Fievet, P. Investigating transport properties of nanofiltration membranes by means of a steric, electric and dielectric exclusion model. J. Membr. Sci. 2005, 252, 77–88. [Google Scholar] [CrossRef]
  6. Bandini, S.; Vezzani, D. Nanofiltration modeling: The role of dielectric exclusion in membrane characterization. Chem. Eng. Sci. 2003, 58, 3303–3326. [Google Scholar] [CrossRef]
  7. Bowen, W.R.; Mukhtar, H. Characterisation and prediction of separation performance of nanofiltration membranes. J. Membr. Sci. 1996, 112, 263–274. [Google Scholar] [CrossRef]
  8. Silva, V.; Geraldes, V.; Brites Alves, A.M.; Palacio, L.; Prádanos, P.; Hernández, A. Multi-ionic nanofiltration of highly concentrated salt mixtures in the seawater range. Desalination 2011, 277, 29–39. [Google Scholar] [CrossRef]
  9. Schaep, J.; Vandecasteele, C.; Mohammad, A.W.; Bowen, W.R. Modelling the retention of ionic components for different nanofiltration membranes. Sep. Purif. Technol. 2001, 22–23, 169. [Google Scholar] [CrossRef]
  10. Zerafat, M.M.; Shariati-Niassar, M.; Hashemi, S.J.; Sabbaghi, S.; Ismail, A.F.; Matsuura, T. Mathematical modeling of nanofiltration for concentrated electrolyte solutions. Desalination 2013, 320, 17–23. [Google Scholar] [CrossRef]
  11. Dutournié, P.; Limousy, L.; Blel, W.; Déon, S.; Fievet, P. Understanding of ion transport in a na-mordenite membrane: Use of numerical modeling to estimate surface-solute interactions in the pore. Ind. Eng. Chem. Res. 2014, 53, 8221–8227. [Google Scholar] [CrossRef]
  12. Bowen, W.R.; Mohammad, A.W.; Hilal, N. Characterisation of nanofiltration membranes for predictive purposes–Use of salts, uncharged solutes and atomic force microscopy. J. Membr. Sci. 1997, 126, 91–105. [Google Scholar] [CrossRef]
  13. Idil Mouhoumed, E.; Szymczyk, A.; Schäfer, A.; Paugam, L.; La, Y.H. Physico-chemical characterization of polyamide NF/RO membranes: Insight from streaming current measurements. J. Membr. Sci. 2014, 461, 130–138. [Google Scholar] [CrossRef] [Green Version]
  14. Datta, S.; Conlisk, A.T.; Kanani, D.M.; Zydney, A.L.; Fissell, W.H.; Roy, S. Characterizing the surface charge of synthetic nanomembranes by the streaming potential method. J. Colloid Interface Sci. 2010, 348, 85–95. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Déon, S.; Escoda, A.; Fievet, P. A transport model considering charge adsorption inside pores to describe salts rejection by nanofiltration membranes. Chem. Eng. Sci. 2011, 66, 2823–2832. [Google Scholar] [CrossRef]
  16. Efligenir, A.; Fievet, P.; Déon, S.; Salut, R. Characterization of the isolated active layer of a NF membrane by electrochemical impedance spectroscopy. J. Membr. Sci. 2015, 477, 172–182. [Google Scholar] [CrossRef]
  17. Montalvillo, M.; Silva, V.; Palacio, L.; Calvo, J.I.; Carmona, F.J.; Hernández, A.; Prádanos, P. Charge and dielectric characterization of nanofiltration membranes by impedance spectroscopy. J. Membr. Sci. 2014, 454, 163–173. [Google Scholar] [CrossRef]
  18. Bason, S.; Oren, Y.; Freger, V. Characterization of ion transport in thin films using electrochemical impedance spectroscopy: II: Examination of the polyamide layer of RO membranes. J. Membr. Sci. 2007, 302, 10–19. [Google Scholar] [CrossRef]
  19. Escoda, A.; Déon, S.; Fievet, P. Assessment of dielectric contribution in the modeling of multi-ionic transport through nanofiltration membranes. J. Membr. Sci. 2011, 378, 214–223. [Google Scholar] [CrossRef]
  20. Lanteri, Y.; Szymczyk, A.; Fievet, P. Influence of steric, electric, and dielectric effects on membrane potential. Langmuir 2008, 24, 7955–7962. [Google Scholar] [CrossRef]
  21. Dutournié, P.; Déon, S.; Limousy, L. Understanding the separation of anion mixtures by TiO2 membranes: Numerical investigation and effect of alkaline treatment on physicochemical properties. Chem. Eng. J. 2019, 363, 365–373. [Google Scholar] [CrossRef]
  22. Déon, S.; Dutournié, P.; Limousy, L.; Bourseau, P. Transport of salt mixtures through nanofiltration membranes: Numerical identification of electric and dielectric contributions. Sep. Purif. Technol. 2009, 69, 225–233. [Google Scholar] [CrossRef]
  23. Déon, S.; Escoda, A.; Fievet, P.; Dutournié, P.; Bourseau, P. How to use a multi-ionic transport model to fully predict rejection of mineral salts by nanofiltration membranes. Chem. Eng. J. 2012, 189–190, 24–31. [Google Scholar] [CrossRef]
  24. Said, A.; Limousy, L.; Nouali, H.; Michelin, L.; Halawani, J.; Toufaily, J.; Hamieh, T.; Dutournié, P.; Daou, T.J. Synthesis of mono- and bi-layer MFI zeolite films on macroporous alumina tubular supports: Application to nanofiltration. J. Cryst. Growth 2015, 428, 71–79. [Google Scholar] [CrossRef]
  25. Said, A.; Daou, T.J.; Limousy, L.; Bikai, J.; Halwani, J.; Toufaily, J.; Hamieh, T.; Dutournié, P. Surface energy modification of a Na-mordenite thin layer treated by an alkaline solution. Mater. Express 2015, 5, 451–456. [Google Scholar] [CrossRef]
  26. Bowen, W.R.; Welfoot, J.S. Modelling the performance of membrane nanofiltration-critical assessment and model development. Chem. Eng. Sci. 2002, 57, 1121–1137. [Google Scholar] [CrossRef]
  27. Bowen, W.R.; Mohammad, A.W. Diafiltration by nanofiltration: Prediction and optimization. AIChE J. 1998, 44, 1799–1812. [Google Scholar] [CrossRef]
  28. Ferry, J.D. Statistical evaluation of sieve constants in ultrafiltration. J. Gen. Physiol. 1935, 20, 95–104. [Google Scholar] [CrossRef] [Green Version]
  29. Born, M. Volumen and hydratationswärme der ionen. Z. Physik. Chem. 1920, 1, 45–48. [Google Scholar] [CrossRef]
  30. Vezzani, D.; Bandini, S. Donnan equilibrium and dielectric exclusion for characterization of nanofiltration membranes. Desalination 2002, 149, 477–483. [Google Scholar] [CrossRef]
  31. Yaroshchuk, A.E. Dielectric exclusion of ions from membranes. Adv. Colloid Interface Sci. 2000, 85, 193–230. [Google Scholar] [CrossRef]
  32. Déon, S.; Dutournié, P.; Bourseau, P. Modeling nanofiltration with Nernst-Planck approach and polarization layer. AIChE J. 2007, 53, 1952–1969. [Google Scholar] [CrossRef]
  33. Petkovska, M. Discrimination between adsorption isotherm models based on nonlinear frequency response results. Adsorption 2014, 20, 385–395. [Google Scholar] [CrossRef]
  34. Dutournié, P.; Limousy, L.; Zouaoui, N.; Mahzoul, H.; Chevereau, E. Facilitated transport of monovalent salt mixture through ultrafiltration Na-Mordenite membrane: Numerical investigations of electric and dielectric contributions. Desalination 2011, 280, 397–402. [Google Scholar] [CrossRef]
  35. Dutournié, P.; Limousy, L.; Anquetil, J.; Déon, S. Modification of the selectivity properties of tubular ceramic membranes after alkaline treatment. Membranes 2017, 7, 65. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. (a) Rejection curves of NaF obtained with TiO2 membrane for two feed concentrations (20 and 100 mol m−3). Symbols: experimental values and lines: simulated curves obtained with the various best-fitted parameters (Xd and εp) plotted in (b).
Figure 1. (a) Rejection curves of NaF obtained with TiO2 membrane for two feed concentrations (20 and 100 mol m−3). Symbols: experimental values and lines: simulated curves obtained with the various best-fitted parameters (Xd and εp) plotted in (b).
Membranes 11 00726 g001
Figure 2. Evolution of minimum and maximum normalised volumetric charge densities with concentration obtained by fitting rejection of NaCl by TiO2 membrane and the corresponding best-fitted adsorption isotherms.
Figure 2. Evolution of minimum and maximum normalised volumetric charge densities with concentration obtained by fitting rejection of NaCl by TiO2 membrane and the corresponding best-fitted adsorption isotherms.
Membranes 11 00726 g002
Figure 3. Evolution of minimum and maximum normalised volumetric charge density with concentration and the corresponding best-fitted Langmuir–Freundlich isotherms obtained by filtration of (a) NaI by the TiO2 membrane and (b) NaF by the Na-mordenite membrane.
Figure 3. Evolution of minimum and maximum normalised volumetric charge density with concentration and the corresponding best-fitted Langmuir–Freundlich isotherms obtained by filtration of (a) NaI by the TiO2 membrane and (b) NaF by the Na-mordenite membrane.
Membranes 11 00726 g003
Figure 4. (a) Evolution of minimum and maximum normalised volumetric charge density with concentration and (b) rejection curve obtained experimentally and simulated from the corresponding best-fitted Langmuir–Freundlich isotherms obtained by filtration of Na2SO4 0.44 mol m−3 by the MFI–MFI membrane.
Figure 4. (a) Evolution of minimum and maximum normalised volumetric charge density with concentration and (b) rejection curve obtained experimentally and simulated from the corresponding best-fitted Langmuir–Freundlich isotherms obtained by filtration of Na2SO4 0.44 mol m−3 by the MFI–MFI membrane.
Membranes 11 00726 g004
Figure 5. (a) Evolution of minimum and maximum normalised volumetric charge density with concentration and (b) rejection curve obtained experimentally and simulated from the corresponding best-fitted Langmuir–Freundlich isotherms obtained by filtration of NaI 2 mol m−3 by the MFI–MFI membrane.
Figure 5. (a) Evolution of minimum and maximum normalised volumetric charge density with concentration and (b) rejection curve obtained experimentally and simulated from the corresponding best-fitted Langmuir–Freundlich isotherms obtained by filtration of NaI 2 mol m−3 by the MFI–MFI membrane.
Membranes 11 00726 g005
Figure 6. Rejection curves obtained experimentally and simulated from the corresponding best-fitted Langmuir–Freundlich isotherms obtained by filtration by the MFI–MFI membrane of NaF (a) 50 mol m−3 and (b) 2 mol m−3.
Figure 6. Rejection curves obtained experimentally and simulated from the corresponding best-fitted Langmuir–Freundlich isotherms obtained by filtration by the MFI–MFI membrane of NaF (a) 50 mol m−3 and (b) 2 mol m−3.
Membranes 11 00726 g006
Table 1. Structural properties (Lp and rp) and maximum VB12 rejection of the studied membranes.
Table 1. Structural properties (Lp and rp) and maximum VB12 rejection of the studied membranes.
MembranesTiO2Na-MordeniteMFI–MFI
Lp × 1014 (m3 m−2)4.81.30.8
RVB12 max (%)62336
rp (nm)1.458.52.15
Table 2. Parameters of the Langmuir–Freundlich isotherms obtained with the various salts and membranes.
Table 2. Parameters of the Langmuir–Freundlich isotherms obtained with the various salts and membranes.
MembranesSaltsQmaxKLFn
TiO2NaCl3.371.200.022
NaI2.131.100.151
Na2SO41.121.050.057
Na-MORNaCl2.011.29−0.010
NaF3.191.150.051
MFI/MFINaI2.131.000.102
NaF1.261.100.050
Na2SO41.870.990.274
Table 3. Minimum and maximum values of εp and Xd/Cf estimated with the various salts and membranes.
Table 3. Minimum and maximum values of εp and Xd/Cf estimated with the various salts and membranes.
MembranesSaltsεp minεp maxXd/Cf minXd/Cf max
TiO2NaCl77.678.4−3.80
NaI78.478.4−2.00
Na2SO477.878.4−1.0+0.36
Na-MORNaCl77.678.4−2.30
NaF73.378.4−3.90
MFI/MFINaI77.678.4−2.30
NaF77.778.4−1.30
Na2SO477.978.4−1.9+0.25
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dutournié, P.; Daou, T.J.; Déon, S. A Novel Numerical Procedure to Estimate the Electric Charge in the Pore from Filtration of Single-Salt Solutions. Membranes 2021, 11, 726. https://doi.org/10.3390/membranes11100726

AMA Style

Dutournié P, Daou TJ, Déon S. A Novel Numerical Procedure to Estimate the Electric Charge in the Pore from Filtration of Single-Salt Solutions. Membranes. 2021; 11(10):726. https://doi.org/10.3390/membranes11100726

Chicago/Turabian Style

Dutournié, Patrick, T. Jean Daou, and Sébastien Déon. 2021. "A Novel Numerical Procedure to Estimate the Electric Charge in the Pore from Filtration of Single-Salt Solutions" Membranes 11, no. 10: 726. https://doi.org/10.3390/membranes11100726

APA Style

Dutournié, P., Daou, T. J., & Déon, S. (2021). A Novel Numerical Procedure to Estimate the Electric Charge in the Pore from Filtration of Single-Salt Solutions. Membranes, 11(10), 726. https://doi.org/10.3390/membranes11100726

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop