Next Article in Journal
Effects of Different TiO2/CNT Coatings of PVDF Membranes on the Filtration of Oil-Contaminated Wastewaters
Previous Article in Journal
Engineering of Multifunctional Nanocomposite Membranes for Wastewater Treatment: Oil/Water Separation and Dye Degradation
Previous Article in Special Issue
Modelling the Performance of Electrically Conductive Nanofiltration Membranes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Ionic Conductivity of Cation Exchange Membranes in Mixed Sulfates Using the Two-Phase Model

Ministry of Education’s Research Center for Comprehensive Utilization and Clean Process Engineering of Phosphorous Resources, School of Chemical Engineering, Sichuan University, Chengdu 610065, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Membranes 2023, 13(10), 811; https://doi.org/10.3390/membranes13100811
Submission received: 10 August 2023 / Revised: 15 September 2023 / Accepted: 19 September 2023 / Published: 26 September 2023

Abstract

:
The concentration dependence of the conductivity of ion exchange membranes (IEMs), as well as other transport properties, has been well explained by the contemporary two-phase model (Zabolotsky et al., 1993) considering a gel phase and an inter-gel phase filled with electroneutral solution. Here, this two-phase model has been adopted and first applied in electrolytes containing mixed counter-ions to investigate the correlation between the membrane ionic conductivity and its microstructure. For three representative commercial cation exchange membranes (CEMs), the total membrane conductivity ( κ T ) when in equilibrium with mixed MgSO4 + Na2SO4 and H2SO4 + Na2SO4 electrolytes could be well predicted with the experimental composition of counter-ions in the gel and inter-gel phase, as well as the counter-ion mobility in the gel phase when the membrane is in a single electrolyte. It is found that the volume fraction of the inter-gel phase ( f 2 ) has little impact on the predicted results. The accuracy of the model can be largely improved by calculating the inter-gel phase conductivity ( κ i n ) with the ionic mobility being the same as that in the external solution (obtained via simulation in the OLI Studio), rather than simply as equivalent to the conductivity of the external solution ( κ s ). Moreover, a nonlinear correlation between the CEMs’ conductivities and the counter-ion composition in the gel phase is observed in the mixed MgSO4 + Na2SO4 solution, as well as for the Nafion117 membrane in the presence of sulfuric acid. For CEMs in mixed MgSO4 + Na2SO4 electrolytes, the calculated conductivity values considering the interaction parameter σ , similar to the Kohlrausch’s law, are closer to the experimental ones. Overall, this work provides new insights into membrane conductivity with mixed counter-ions and testifies to the applicability of the contemporary two-phase model.

1. Introduction

Ion exchange membranes (IEMs) are one of the most advanced membranes and now are widely applied in several separation technologies due to their high permselectivity and energy-efficient performance [1,2,3]. It is known that the IEMs, including the so-called homogeneous ones, are actually spatially nonuniform [4]. This nonuniformity on the sub-microscopical scale will eventually influence the properties of IEMs. Therefore, many researchers in the IEMs field have been exploring the correlation between the fundamental membrane properties and their microstructures by introducing static and dynamic parameters [5,6,7,8].
In early studies, the IEMs are considered as a single homogeneous mixture of matrix polymer chains, ion-exchanged and mobile ions, and water [9,10]. However, the oversimplification of this mathematical model makes it difficult to predict the phenomenological coefficients. Considering the membrane microstructural inhomogeneity, the IEMs could be regarded as a system of two or several microphases [11,12,13]. Then the membrane properties could be correlated with the microphases within the membrane. Ionic conductivity is one of the most essential properties of IEMs [14,15], which has been found as a function of the corresponding microphase properties [16]. According to the different ionic states (ions due to ion exchange with fixed-charge groups and due to the Donnan sorption) within IEMs, A. Yamauchi et al. divided the total conductivity of membranes into two parts: the ion-exchanged conductivity and the Donnan electrolytes conductivity [17]. Based on the Nernst–Planck equation [18], V. I. Zabolotsky et al. [7] developed the contemporary two-phase model for the transport properties, including conductivity, by introducing the structural parameters into the theoretical framework. The model was initially formulated and tested for 1:1 strong electrolytes. N.P. Gnusin et al. [19] successfully applied this contemporary two-phase model in describing the membrane electro-transport properties and in solving the dynamic problems of electrodialysis. Y. Sedkaoui et al. [20] determined the volume fraction of the inter-gel phase of membranes in NaCl and KCl solutions and provided a literature review of these values in NaCl solution. Indeed, there are experimental results reported by V. V. Sarapulova et al. for the membrane conductivity and diffusion permeability treated using the two-phase model for the cases where the solution contains NaCl, CaCl2, or Na2SO4 [21,22,23]. However, this model has only been verified in IEMs in equilibrium with a single electrolyte. When such a model is applied to weak electrolytes, or to acids and alkalis instead of salts, the researchers might produce results not anticipated by the model since the model does not take into account the additional chemical interactions. Except for the widely accepted two-phase model, there are also papers developing similar multiphase models [24]. T. W. Xu et al. [6] also proposed a simple evaluation method based on conductivity for evaluating the membrane microstructure and transport parameters. However, due to a series of assumptions made in this framework, this method has been less considered. Nevertheless, the membrane behavior in the mixture of these electrolytes was not studied from this theoretical point of view based on the two-phase model.
Recent studies concentrate more on membrane conductivity measurement techniques and the concentration dependence related to the membrane microstructure [20,25,26,27,28,29]. Theoretical considerations about the correlation between the membrane structure and its ionic conductivity are relatively less reported [30]. In this work, the contemporary two-phase model has been adopted and applied to three representative commercial cation exchange membranes (CEMs) in two different mixed electrolytes (MgSO4/Na2SO4 and H2SO4/Na2SO4) under a total sulfate concentration range of 0.1~1.0 M. It is known that the proton exhibits an unusual transport mechanism: except for the diffusion, migration, and convection in a typical electro-membrane process [9], there also exists the so-called structural diffusion (Grotthuss) mechanism [31]. Therefore, the contemporary two-phase model is also applied for the mixed electrolytes containing sulfuric acid, which has been rarely reported in literature. Overall, the structural parameters in the contemporary two-phase model will be systematically investigated in the mixed electrolytes, and the correlation between the CEM microstructure and the ionic conductivity will also be discussed here.

2. Theoretical Background

2.1. Contemporary Two-Phase Model

The importance of studying the effect of IEM structural inhomogeneity has drawn much attention [13]. Based on the thermodynamic knowledge of IEMs [32], the membrane is generally considered as being composed of two microphases, according to the contemporary two-phase model: the gel phase composed of functional charged groups and hydrophilic matrix polymer chain, and the inter-gel phase filled with the Donnan electrolytes [7].
According to the work of V. I. Zabolotsky et al. [7], as shown in Figure 1, the membrane total conductivity ( κ T ) is the sum of corresponding conductivities from two parts: the gel phase conductivity ( κ g ) and the inter-gel phase conductivity ( κ i n ). Here, κ i n is implicitly taken as the conductivity of the external solution ( κ s ). Thus, the total conductivity is further expressed by introducing the structural parameters (Equation (1)).
κ T = ( f 1 ( κ g ) α + f 2 ( κ s ) α ) 1 / α
where f 1 , f 2 represents the volume ratio of the gel and the inter-gel phase, respectively, and f 1 + f 2 = 1 . α is the structural parameter, which reflects the reciprocal arrangement of the microphase elements within the membrane. Two extreme conditions, α = 1 and −1, correspond to the situation when all the gel and the inter-gel microphase elements are arranged in parallel and in series, respectively, to the flow direction of the current. Actually, Equation (1) is generalized from these two extreme conditions to reflect the actual arrangement situation of the microphase elements in the membrane.
The previous study, ref. [7] found that the total conductivity depends slightly on α . When α 1 , Equation (1) becomes:
κ T = ( κ g ) f 1 ( κ s ) f 2
Then a linear relationship can be anticipated in a ln κ T ln κ s correlation:
l n κ T = f 1 l n κ g + f 2 l n κ s
However, the correlation above (Equation (1)) has only been verified in a single-electrolyte solution with one monovalent counter-ion. As far as the authors would know, this work is the first to investigate its applicability in a mixed electrolytes system.
In equilibrium with a single electrolyte, the membrane’s gel phase conductivity is dominated by the counter-ion’s concentration and its mobility in the gel phase:
κ i g = z i F c i g u i g = F c g u i g
where z i is the valence of counter-ion i, F is the Faraday constant, c i g is the molarity of counter-ion i in the gel phase (mol·m−3). u i g is the mobility of counter-ion i and c g is the equivalent counter-ion concentration in the gel phase (equ·m−3).
In our recent work [33], for mixed electrolytes containing two different counter-ions, the conductivity of the gel phase is determined by the mobilities of the two counter-ions and the interaction between them. Then, two situations exist:
(I) If the interaction between different counter-ions is negligible, the conductivity of the gel phase is expressed as:
κ i , j g = z i F c i g u i g + z j F c j g u j g
Then, the total conductivity of the IEM in equilibrium with mixed electrolytes could be expressed by combining Equations (2a) and (4).
κ T = ( z i F c i g u i g + z j F c j g u j g ) f 1 ( κ s ) f 2
Here, the mobility of counter-ion i ( u i g ) is obtained through Equation (3) when the IEM is in a single electrolyte. The concentration of counter-ion ( c i g ) is determined experimentally.
(II) If the interaction between counter-ions is non negligible, the Kohlrausch’s law is then employed for the description of the counter-ion mobility variations [34].
κ i , j g = x i m κ i g σ x j m + x j m κ j g
σ is the interaction parameter, which is determined by the experimental conductivity data. x i m is equivalent molar ratio of counter-ion i in the gel phase and is expressed as x i m = z i c i g / ( z i c i g + z j c j g ) . Then, Equation (2a) is further expressed as:
κ T = ( x i m κ i g σ x j m + x j m κ j g ) f 1 ( κ s ) f 2

2.2. The Inter-Gel Phase Conductivity

In Section 2.1, the conductivity of the inter-gel phase is simply considered equivalent to the conductivity of the external solution, which is implicitly indicated in the work of V.I. Zabalotsky et al. [7]. However, the value of the inter-gel phase conductivity has a direct influence on the structural parameter fitting results in the contemporary two-phase model. Here, the inter-gel phase conductivity is also calculated via two different approaches.
In the first approach, the ionic mobilities in the infinitely diluted solution are applied in calculating the inter-gel phase conductivity ( κ i n ):
κ i n = F ( z + c + D u + d i l u t e + z c D u d i l u t e )
where c + D , c D are the concentrations of cations and anions, respectively, due to the Donnan sorption. u + d i l u t e , u d i l u t e are the mobility of corresponding ions in the infinitely diluted solution. Table 1 lists the values of ionic mobilities employed in this work [3].
In the second approach, the conductivity of the inter-gel phase is also calculated by assuming that the mobilities of ions comprising the Donnan electrolytes are the same as the corresponding ions’ mobilities in the external solution. In a single-electrolyte solution, its ionic conductivity is expressed as:
κ s = F ( z + c + s u + s + z c s u s )
c + s , c s are the concentrations of cations and anions in the external solution; u + s , u s are the ionic mobilities in the external solution, respectively. According to the electro-neutrality of the external solution, z + c + s = z c s = c e q u s , Equation (9) becomes:
κ s = F c e q u s ( u + s + u s )
where c e q u s is the equivalent electrolyte concentration in the external solution. Similarly, the conductivity of the inter-gel phase is expressed as:
κ i n = F c e q u D ( u + D + u D )
where c e q u D is the equivalent Donnan electrolyte concentration in the inter-gel phase. u + D , u D are the ionic mobilities within the Donnan electrolytes, respectively.
Considering that the ionic mobilities are consistent in both the Donnan electrolytes and the external solution ( u + s = u + D , u s = u D ), the value of the inter-gel phase conductivity could be calculated from the conductivity of the external solution, as:
κ i n = κ s c e q u D c e q u s
Then the total conductivity of IEMs could be predicted by the fitting results of the κ T κ i n correlation as well as by replacing the κ s in Equations (5) and (7) with κ i n .

3. Materials and Methods

3.1. Cation Exchange Membranes

Three representative commercial CEMs with different microstructures are investigated in this work. The CSE membrane (Astom Co., Tokyo, Japan) is a standard homogeneous CEM used for desalination. Nafion117 is a representative perfluorosulfonic acid membrane from Chemours, America. A heterogeneous membrane, 3361BW (Shanghai Shanghua Water Treatment Material Co., Shanghai, China), is also investigated. The fundamental properties of the membranes have been determined and reported in our previous report [33], and are also provided here in Table 2.

3.2. Membrane Ionic Composition

The composition of both the Donnan electrolytes within the inter-gel phase and counter-ions in the gel phase has been investigated under a total sulfate concentration range of 0.1~1.0 M. Figure 2 schematically shows the experimental procedures for determining the CEM ionic composition. Before the partitioning experiments, all the membrane samples are cut into 4 × 4 cm and immersed in 0.5 M Na2SO4 solution for 48 h to convert the membrane into Na+ form. Then, all the samples are washed several times with deionized water and transferred to the mixed-electrolyte solutions. After equilibrium in a mixed-electrolytes solution for 48 h, the membrane surface moisture is quickly removed with filter papers. Next, the membranes are immersed in excessive deionized water for 24 h to completely elute the Donnan electrolytes within the inter-gel phase, and all the elution is collected to give a final volume of 250 mL. Then, the wet membranes are prepared for another step: the determination of the counter-ions’ composition in the gel phase. For the mixed MgSO4/Na2SO4 electrolytes system, a membrane sample is ion-exchanged with 100 mL 0.5 M MgSO4 solution twice, while the membrane sample in the mixed H2SO4/Na2SO4 electrolytes system is ion-exchanged with 100 mL 0.5 M H2SO4 twice. All the elution solutions are then gathered to give a total volume of 250 mL. The concentration of cations in the solution is analyzed by ICP-OES (inductively coupled plasma, optically emission spectrometer, PE Avio 200, Waltham, MA, USA) or an automatic potentiometric titrator (916, Metrohm, Herisau, Switzerland).

3.3. Conductivity Measurement

Based on electrochemical impedance spectroscopy (EIS), a direct-contact method is applied to determine the through-plane ionic conductivity of CEMs. The conductivity cell installation (Figure 3) is the same as that used in our previous work [33]. After equilibrium in the external solution, the investigated membrane is removed from the solution and directly sandwiched between two electrodes. Then, the EIS spectra are acquired from an electrochemical workstation (Zennium IM6, Zahner, Gundelsdorf, Germany), together with a buffer unit supplying an AC current of 5 mV amplitude in the frequency range of 1~300 kHz. All the measurements are conducted at room temperature (25 °C). More details can be found in our previous report [33].

4. Results and Discussions

4.1. Determination of the Membrane Structural Parameters

Before applying the contemporary two-phase model in mixed electrolytes, the membrane structural parameters, f 2 and α , fitted from either the nonlinear (Equation (1)) or the linear (Equation (2b)) correlation, need to be determined. The total conductivity ( κ T ) of CEMs in a single electrolyte is measured under a sulfate concentration range of 0.1 ~ 1.0 M. The actual conductivity of the external solution ( κ s ) is obtained from the OLI Studio simulation results. In this work, the membrane structural parameters α and the volume fraction of two membrane microphases f i are obtained through two different approaches.

4.1.1. Fitted from the External Solution Conductivity κ s

As illustrated in Figure 4a, a good linear relationship can be observed in the ln κ T ln κ s correlation for the heterogeneous membrane 3361BW under different single-electrolyte solutions. The slop represents the volume fraction of the inter-gel phase ( f 2 ). It is clear that f 2 changes with the type of the external electrolyte. The value of f 2 for the 3361BW membrane reaches its maximum in the H2SO4 solution (~0.48). However, for the homogeneous CEMs (CSE and Nafion117, Figure 4b,c), generally poor linear fitting results are observed.
Considering that the linear relationship is not so obvious for all three investigated CEMs under different electrolytes, the nonlinear correlation method (Equation (1)) has also been applied in acquiring the membrane structural parameters. Table 3 shows the nonlinear fitting results for the CEMs under different electrolytes. Though the variance of the nonlinear fitting is overall larger than the linear fitting results, the value of f 2 is too large to be real for some CEMs under a certain electrolyte. For example, for the heterogenous membrane 3361BW in sulfuric acid, f 2 = 0.818, for the CSE membrane under MgSO4, f 2 = 0.998, and under H2SO4, f 2 = 0.999. These unrealistic values of f 2 would be very likely numerical artifacts of the nonlinear fitting. In the report of V. I. Zabolotsky et al. [7], the value of α is in the range of 0.1~0.3 for most CEMs. It is observed here that the fitting values of α for CEMs could be negative, especially in sulfuric acid solution. These negative values can be anticipated and have been reported in the work of N.D. Pismenskaya et al. [30], and they are discussed in the review of V. Nikonenko et al. [35].

4.1.2. Fitted from the Inter-Gel Phase Conductivity κ i n

In Section 4.1.1, the conductivity of the inter-gel phase ( κ i n ) is simply taken as equivalent to the external solution conductivity κ s . According to our previous research [33], the composition of the Donnan electrolytes within the membrane is in fact not the same as the external solution. Therefore, in order to better describe the relationship between the CEM conductivity and the membrane microstructure, κ i n is also acquired by assuming that the ionic mobility of the Donnan electrolytes is the same as that in the external electrolyte.
The membrane structural parameter results calculated via two different approaches (Equations (8) and (12)) are shown in Table 4 and Table 5, respectively. Different from the results fitted with the external solution conductivity κ s (Table 3), the f 2 value here is in the range of 0 ~ 0.3, which is consistent with the literature results [4,7,36]. Interestingly, the values of α for CEMs in H2SO4 electrolyte solution here are all greater than 1, which is not reported in previous research. Numerical artifact alone would not explain these unphysical values that contradict the model. At the moment, we do not have a clear understanding of why α would be larger than 1. Moreover, there are some cases in which the nonlinear fitting process does not converge (the missing data in Table 5). Therefore, the membrane structural parameters fitted from the linear method are applied in predicting the total conductivity of CEMs under mixed electrolytes in subsequent sections.

4.2. Prediction of Membrane Conductivity in Mixed Electrolytes

4.2.1. κ T Prediction with the External Solution Conductivity κ S

In Section 4.1, the membrane structural parameters fitted via different approaches with the membrane conductivity data in the single-electrolyte solution have been systematically discussed. Then the conductivity of CEMs could be predicted through Equation (5) or Equation (7).
Figure 5 and Figure 6 show the variation of the total membrane conductivities with the cation composition of the external mixed-sulfate solutions. Here, the total membrane conductivities are obtained both experimentally (via the direct-contact method) and theoretically (Equation (5)). It is indicated from the results that the calculated values provide a relatively accurate prediction of the trend in the total conductivity of CEMs under mixed electrolytes (the red solid line in the graph represents the trend of the predicted values). Interestingly, the predicted values are overall larger than the experimental ones, whether in the mixed MgSO4 + Na2SO4 or the H2SO4 + Na2SO4 electrolytes.
As can be seen in Section 4.1 (Table 4 and Table 5), the value of f 2 will change with the type of external electrolyte. Therefore, the fitted values of f 2 from different electrolytes are adopted here in predicting the total conductivities of CEMs under mixed electrolytes. The orange, blue, and red symbols in Figure 5 and Figure 6 represent the total membrane conductivities calculated with different f 2 values. It is clear that the variation of f 2 has little impact on the predicted membrane ionic conductivity, especially for the perfluorosulfonic acid membrane Nafion117. Overall, the difference between the predicted and experimental values is relevant to the membrane’s microscopic channels and the water content. For the heterogeneous membrane 3361BW, with a larger inter-gel phase (Table 4 and Table 5) and higher water content (Table 2), the accuracy of the contemporary two-phase model is obviously lower than for the homogeneous CEMs.

4.2.2. κ T Prediction with the Inter-Gel Phase Conductivity κ i n

In Section 4.2.1, the predicted total conductivities of CEMs are obtained by simply incorporating the conductivity of the external solution ( κ s ) into the contemporary two-phase model. Here, the membrane total conductivity is also predicted by assuming that the ionic mobility in the inter-gel phase is the same as that in the external solution.
As mentioned in Section 4.1.2, there are two different approaches to calculating the ionic mobility in the inter-gel phase. In the first approach, the ionic mobility under infinite dilution is used for calculating the κ i n (Equation (8)). As shown in Figure 7a, the results predicted via this approach with the contemporary two-phase model are contrary to the experimental ones for the CSE membrane under mixed MgSO4 + Na2SO4 solutions. This discrepancy can be attributed to the over-simplification of the ionic mobility in the inter-gel phase and that under infinite dilution (Table 1). The results for other CEMs (3361BW and Nafion117) are provided in the Supporting Information (Figures S1 and S2).
In the second approach, the mobilities of ions in the inter-gel phase are represented with the actual mobilities of ions in the corresponding external solutions (Equations (10)–(12)). Then, κ i n is obtained from the simulated conductivity data of the external solution ( κ s ) analyzed in the OLI Studio, together with the experimental composition of the Donnan electrolyte within the membrane ( c e q u D , Equation (12)). It is obvious from the revised results (Figure 7b) that the accuracy of the contemporary two-phase model can be largely improved by employing the actual ionic mobility in the external solution, especially for the CSE membrane in the di-/monovalent cations (Mg2+/Na+) system. Thus, the second approach is reasonable and suggested.
Furthermore, the numerical results of the cations mobility in the gel phase calculated from two different approaches (Table 6) can better explain the discrepancy observed in Figure 7a. First, the fitting results of the ln κ T ln κ i n correlation (Equation (2b)) for a single-electrolyte solution give the conductivity of the membrane gel phase ( κ g ). Then, values of the cation mobility in the membrane gel phase ( u i g ) can be calculated via Equation (3). The two indications in Table 6, “infinite dilution” and “OLI simulation”, represent the approaches when κ i n is calculated with the ionic mobility in the inter-gel phase as the value under infinite dilution and as the value the same as that in the external solution simulated in the OLI Studio. It can be concluded that the main reason for the conductivity discrepancy of the CSE membrane is the different Mg2+ mobility in the membrane phase.
Figure 8 and Figure 9 show the total conductivity results of the CEMs under two different kinds of mixed electrolytes. Compared to the method of simply incorporating the solution conductivity ( κ s ) into the contemporary two-phase model (Figure 5 and Figure 6), the results obtained with the ionic mobility in the inter-gel phase as identical to that in the external solution are closer to the experimental values, especially for the mixed-electrolyte solutions containing sulfuric acid. However, for the heterogenous membrane 3361BW under the mixed MgSO4 + Na2SO4 electrolytes, the accuracy of this approach is lower than that for the homogenous CEMs.
The numerical data of predicted CEMs’ total conductivity in mixed-sulfate solutions of different total sulfate concentration are provided in the Supporting Information (Tables S1–S6).

4.3. Correlation between Conductivity and Counter-Ions Composition in the Gel Phase

In the discussions above, the interaction between different counter-ions in the membrane gel phase is ignored. However, a nonlinear correlation between the conductivity and the counter-ions’ composition in the gel phase of IEMs has been proved both in our previous work [33,37] and in the literature [34,38]. In this work, Kohlrausch’s law (Equation (6)) is introduced to describe this nonlinear relationship. It is shown in Figure 10 that the nonlinearity between the membrane gel phase conductivity and the cation composition is non-negligible for three CEMs under mixed MgSO4 + Na2SO4 electrolytes, especially when the total sulfate concentration is low (0.1 M).
Interestingly, a perfect linear relationship is observed for the 3361BW and CSE membranes in mixed electrolytes containing sulfuric acid (Figure 11a,b). The gel phase conductivity hardly changes with the total sulfate concentration in the external solution. However, for Nafion117 (Figure 11c), nonlinear correlations can clearly be observed for different total sulfate concentrations, which is consistent with the results in the literature [33,34]. The numerical fitting results are provided in the Supporting Information (Tables S7 and S8).
Therefore, the contemporary two-phase model revised with Kohlrausch`s interaction parameter σ (Equations (6) and (7)) has also been applied in calculating the CEMs’ total conductivities under mixed electrolytes of different compositions. Figure 12 shows the calculated total membrane conductivities. Compared with the linear method (Equation (5) and Figure 8), the total membrane conductivity results calculated with the κ g values obtained from this nonlinear correlation are closer to experimental ones for CEMs under mixed MgSO4 + Na2SO4 electrolytes. The total conductivity values of CEMs in mixed-electrolyte solutions containing sulfuric acid are provided in the Supporting Information (Figure S3). It is observed from these results that the interaction of counter-ions in the membrane gel phase for the Na+/Mg2+ pair is more obvious than that for the Na+/H+ pair.

5. Conclusions

In this work, the contemporary two-phase model of ion exchange membranes (IEMs) has been systematically investigated for three representative cation exchange membranes (CEMs) in equilibrium with mixed-sulfate electrolytes of two different counter-ions (Na+/Mg2+ and Na+/H+) using a total sulfate concentration of 0.1~1.0 M. After experimentally obtaining the membrane ionic composition, the total conductivity of CEMs ( κ T ) in mixed electrolytes could be well predicted via the two-phase model. Though the membrane structural parameter f in this model changes with the nature of the external electrolyte, the volume fraction of the inter-gel phase f 2 has little impact on the calculated membrane total conductivity. The accuracy of the two-phase model could be largely improved by replacing the external solution conductivity ( κ s ) with the inter-gel phase conductivity ( κ i n ). A further improvement can be achieved by calculating the κ i n as the Donnan electrolytes’ conductivity, with the ionic mobility in the inter-gel phase being the same as that in the external solution, rather than simply with the ionic mobilities under infinite dilution, especially for the CSE membrane under mixed MgSO4 + Na2SO4 electrolytes. Moreover, the calculated total conductivity revised with the interaction parameter ( σ ) derived from Kohlrausch’s law is closer to the experimental values in mixed MgSO4 + Na2SO4 electrolytes, while the linear correlation method is more suitable for mixed-sulfate electrolytes in the presence of sulfuric acid. Overall, the contemporary two-phase model could be successfully applied in mixed-sulfate electrolytes for CEMs, and the accuracy of this model is better for mixed mono-/monovalent counter-ion pairs than for di-/monovalent counter-ions pairs.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/membranes13100811/s1, supplementary file.

Author Contributions

Conceptualization, L.W. and T.L.; Data curation, L.W. and H.J.; Formal analysis, L.W. and H.J.; Funding acquisition, T.L.; Investigation, L.W. and H.J.; Methodology, L.W., H.J. and T.L.; Project administration, T.L. and X.W.; Resources, X.W.; Supervision, T.L. and X.W.; Validation, L.W. and H.J.; Visualization, L.W., H.J. and T.L.; Writing—original draft, L.W. and H.J.; Writing—review and editing, T.L. All authors have read and agreed to the published version of the manuscript.

Funding

The authors acknowledge the financial support from the National Natural Science Foundation of China (22108176), and partial support from the Research Funds of Sichuan Province (2022NSFSC0219).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Data will be available upon request.

Acknowledgments

The authors would like to acknowledge the critical suggestions from two reviewers to help improve the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jiang, S.; Sun, H.; Wang, H.; Ladewig, B.P.; Yao, Z. A comprehensive review on the synthesis and applications of ion exchange membranes. Chemosphere 2021, 282, 130817. [Google Scholar] [CrossRef] [PubMed]
  2. Shehzad, M.A.; Yasmin, A.; Ge, X.; Wu, L.; Xu, T. A Review of Nanostructured Ion-Exchange Membranes. Adv. Mater. Technol. 2021, 6, 2001171. [Google Scholar] [CrossRef]
  3. Luo, T.; Abdu, S.; Wessling, M. Selectivity of ion exchange membranes: A review. J. Membr. Sci. 2018, 555, 429–454. [Google Scholar] [CrossRef]
  4. Bryk, M.; Zabolotsky, V.; Atamanenko, D.; Dvorkina, G. Structural inhomogeneity of ion-exchange membranes in swelling state and methods of its investigations. Khim. Tekhnol. Vody 1989, 11, 491–498. [Google Scholar]
  5. Auclair, B.; Nikonenko, V.; Larchet, C.; Metayer, M.; Dammak, L. Correlation between transport parameters of ion-exchange membranes. J. Membr. Sci. 2002, 195, 89–102. [Google Scholar] [CrossRef]
  6. Xu, T.W.; Li, Y.; Wu, L.; Yang, W.H. A simple evaluation of microstructure and transport parameters of ion-exchange membranes from conductivity measurements. Sep. Purif. Technol. 2008, 60, 73–80. [Google Scholar] [CrossRef]
  7. Zabolotsky, V.I.; Nikonenko, V.V. Effect of Structural Membrane Inhomogeneity on Transport-Properties. J. Membr. Sci. 1993, 79, 181–198. [Google Scholar] [CrossRef]
  8. Pourcelly, G.R. Conductivity and selectivity of ion exchange membranes: Structure-correlations. Desalination 2003, 147, 359–361. [Google Scholar] [CrossRef]
  9. Strathmann, H. Ion Exchange Membrane Separation Processes; Elsevier: Amsterdam, The Netherlands, 2004. [Google Scholar]
  10. Lauger, P. Transport phenomena in membranes. Angew. Chem. Int. Ed. Engl. 1969, 8, 42–54. [Google Scholar] [CrossRef]
  11. Narębska, A.; Wódzki, R. Diffusion of electrolytes across inhomogeneous permselective membranes. Angew. Makromol. Chem. 1979, 80, 105–118. [Google Scholar] [CrossRef]
  12. Sata, T. Ion Exchange Membranes: Preparation, Characterization, Modification and Application; The Royal Society of Chemistry: Cambridge, UK, 2004. [Google Scholar]
  13. Petropoulos, J.H. Membrane-Transport Properties in Relation to Microscopic and Macroscopic Structural Inhomogeneity. J. Membr. Sci. 1990, 52, 305–323. [Google Scholar] [CrossRef]
  14. Steymans, C.L. Kurzvorträge zum Hauptthema. Untersuchungen über Leitfähigkeitsverhalten von Austauschermembranen. I. Mitteilung: Messung der Membranleitfähigkeit. Berichte der Bunsen-Gesellschaft fur Physikalische Chemie 1967, 71, 818–822. [Google Scholar]
  15. Koter, S.; Narebska, A. Conductivity of ion-exchange membranes-ii.Mobilities of ions and water. Electrochim. Acta 1987, 32, 455–458. [Google Scholar] [CrossRef]
  16. Narcbska, R.W.A.; Koter, S. Composition and structure of the cation permselective membranes I. Evaluation of the electrochemical models. Angew. Makromol. Chem. 1980, 86, 157–170. [Google Scholar] [CrossRef]
  17. Yamauchi, Y.H.A.A. Additivity between donnan salt and ion-exchanged salt in the specific conductance of membranes. J. Membr. Sci. 1990, 48, 25–31. [Google Scholar]
  18. Buck, R.P. Kinetics of bulk and interfacial ionic motion: Microscopic bases and limits for the Nernst—Planck equation applied to membrane systems. J. Membr. Sci. 1984, 17, 1–62. [Google Scholar] [CrossRef]
  19. Gnusin, N.P.; Berezina, N.P.; Kononenko, N.A.; Dyomina, O.A. Transport structural parameters to characterize ion exchange membranes. J. Membr. Sci. 2004, 243, 301–310. [Google Scholar] [CrossRef]
  20. Sedkaoui, Y.; Szymczyk, A.; Lounici, H.; Arous, O. A new lateral method for characterizing the electrical conductivity of ion-exchange membranes. J. Membr. Sci. 2016, 507, 34–42. [Google Scholar] [CrossRef]
  21. Sarapulova, V.; Pismenskaya, N.; Butylskii, D.; Titorova, V.; Wang, Y.; Xu, T.; Zhang, Y.; Nikonenko, V. Transport and Electrochemical Characteristics of CJMCED Homogeneous Cation Exchange Membranes in Sodium Chloride, Calcium Chloride, and Sodium Sulfate Solutions. Membranes 2020, 10, 165. [Google Scholar] [CrossRef]
  22. Sarapulova, V.; Pismenskaya, N.; Titorova, V.; Sharafan, M.; Wang, Y.; Xu, T.; Zhang, Y.; Nikonenko, V. Transport Characteristics of CJMAED™ Homogeneous Anion Exchange Membranes in Sodium Chloride and Sodium Sulfate Solutions. Int. J. Mol. Sci. 2021, 22, 1415. [Google Scholar] [CrossRef]
  23. Sarapulova, V.V.; Titorova, V.D.; Nikonenko, V.V.; Pismenskaya, N.D. Transport Characteristics of Homogeneous and Heterogeneous Ion-Exchange Membranes in Sodium Chloride, Calcium Chloride, and Sodium Sulfate Solutions. Membr. Membr. Technol. 2019, 1, 168–182. [Google Scholar] [CrossRef]
  24. Veerman, J. The Effect of the NaCl Bulk Concentration on the Resistance of Ion Exchange Membranes—Measuring and Modeling. Energies 2020, 13, 1946. [Google Scholar] [CrossRef]
  25. Kamcev, J.; Sujanani, R.; Jang, E.-S.; Yan, N.; Moe, N.; Paul, D.R.; Freeman, B.D. Salt concentration dependence of ionic conductivity in ion exchange membranes. J. Membr. Sci. 2018, 547, 123–133. [Google Scholar] [CrossRef]
  26. Zhu, S.; Kingsbury, R.S.; Call, D.F.; Coronell, O. Impact of solution composition on the resistance of ion exchange membranes. J. Membr. Sci. 2018, 554, 39–47. [Google Scholar] [CrossRef]
  27. Ziv, N.; Dekel, D.R. A practical method for measuring the true hydroxide conductivity of anion exchange membranes. Electrochem. Commun. 2018, 88, 109–113. [Google Scholar] [CrossRef]
  28. Barragán, V.M.; Izquierdo-Gil, M.A.; Maroto, J.C.; Antoranz, P.; Muñoz, S. Estimation of the through-plane thermal conductivity of polymeric ion-exchange membranes using finite element technique. Int. J. Heat Mass Transf. 2021, 176, 121469. [Google Scholar] [CrossRef]
  29. Díaz, J.C.; Kamcev, J. Ionic conductivity of ion-exchange membranes: Measurement techniques and salt concentration dependence. J. Membr. Sci. 2021, 618, 118718. [Google Scholar] [CrossRef]
  30. Pismenskaya, N.D.; Nevakshenova, E.E.; Nikonenko, V.V. Using a Single Set of Structural and Kinetic Parameters of the Microheterogeneous Model to Describe the Sorption and Kinetic Properties of Ion-Exchange Membranes. Pet. Chem. 2018, 58, 465–473. [Google Scholar] [CrossRef]
  31. Agmon, N. The Grotthuss mechanism. Chem. Phys. Lett. 1995, 244, 456–462. [Google Scholar] [CrossRef]
  32. Ponomarev, O.A.; Ionova, I.A. The thermodynamics of ionomers. Polym. Sci. U.S.S.R. 1974, 16, 1181–1190. [Google Scholar] [CrossRef]
  33. Wu, L.; Luo, T.; Yang, X.; Zhao, H.; Wang, X.; Zhang, Z. Impact of the Donnan electrolytes on selectivity of cation exchange membranes evaluated via the ionic membrane conductivity. Sep. Purif. Technol. 2023, 316, 123816. [Google Scholar] [CrossRef]
  34. Logette, S.; Eysseric, C.; Pourcelly, G.; Lindheimer, A.; Gavach, C. Selective permeability of a perfluorosulphonic membrane to different valency cations. Ion-exchange isotherms and kinetic aspects. J. Membr. Sci. 1998, 144, 259–274. [Google Scholar] [CrossRef]
  35. Nikonenko, V.; Nebavsky, A.; Mareev, S.; Kovalenko, A.; Urtenov, M.; Pourcelly, G. Modelling of Ion Transport in Electromembrane Systems: Impacts of Membrane Bulk and Surface Heterogeneity. Appl. Sci. 2019, 9, 25. [Google Scholar] [CrossRef]
  36. Dimov, A.D.; Alexandrova, I. Electrochemical properties of cation-exchange membranes. J. Chem. Soc. Faraday Trans. 1 Phys. Chem. Condens. Phases 1987, 83, 2841–2845. [Google Scholar] [CrossRef]
  37. Zou, Z.; Wu, L.; Luo, T.; Yan, Z.; Wang, X. Assessment of anion exchange membrane selectivity with ionic membrane conductivity, revised with Manning’s theory or the Kohlrausch’s law. J. Membr. Sci. 2021, 635, 119496. [Google Scholar] [CrossRef]
  38. Chaudhury, S.; Agarwal, C.; Goswami, A. Transport Properties of Multivalent Cations in Nafion-117 Membrane with Mixed Ionic Composition. J. Phys. Chem. B 2015, 119, 10566–10572. [Google Scholar] [CrossRef]
Figure 1. The correlation between the total ionic conductivity of cation exchange membranes (CEMs), κ T , and the gel phase ( κ g ) and the inter-gel phase conductivity ( κ i n ) according to the contemporary two-phase model. New attempts are also summarized here. α > 1 results from nonlinear fitting according to Equation (1) and it contradicts the model.
Figure 1. The correlation between the total ionic conductivity of cation exchange membranes (CEMs), κ T , and the gel phase ( κ g ) and the inter-gel phase conductivity ( κ i n ) according to the contemporary two-phase model. New attempts are also summarized here. α > 1 results from nonlinear fitting according to Equation (1) and it contradicts the model.
Membranes 13 00811 g001
Figure 2. Schematical illustration of the experimental procedures for determining the ionic composition of CEMs after equilibrium in mixed-sulfate solutions.
Figure 2. Schematical illustration of the experimental procedures for determining the ionic composition of CEMs after equilibrium in mixed-sulfate solutions.
Membranes 13 00811 g002
Figure 3. Simple illustration of the cell configuration for the measurement of the through-plane membrane resistance (reproduced with permission from Ref. [33], Elsevier 2023).
Figure 3. Simple illustration of the cell configuration for the measurement of the through-plane membrane resistance (reproduced with permission from Ref. [33], Elsevier 2023).
Membranes 13 00811 g003
Figure 4. The linear correlation between the ionic conductivities of three CEMs, (a) 3361BW, (b) CSE and (c) Nafion117, and the conductivity of an external solution with a single sulfate electrolyte in the ln κ T ln κ s coordinates.
Figure 4. The linear correlation between the ionic conductivities of three CEMs, (a) 3361BW, (b) CSE and (c) Nafion117, and the conductivity of an external solution with a single sulfate electrolyte in the ln κ T ln κ s coordinates.
Membranes 13 00811 g004
Figure 5. The variation of the experimental and the predicted total membrane conductivities with the cation composition of the external mixed MgSO4 + Na2SO4 electrolytes of 0.5 M total sulfate, for the 3361BW membrane (a), CSE membrane (b) and Nafion 117 membrane (c). Here, f(Na) (or f(Mg)) represents the volume fraction of the inter-gel phase in the membrane fitted from the single Na2SO4 (or MgSO4) electrolyte solution, and f(Mg/Na) is the arithmetic average value.
Figure 5. The variation of the experimental and the predicted total membrane conductivities with the cation composition of the external mixed MgSO4 + Na2SO4 electrolytes of 0.5 M total sulfate, for the 3361BW membrane (a), CSE membrane (b) and Nafion 117 membrane (c). Here, f(Na) (or f(Mg)) represents the volume fraction of the inter-gel phase in the membrane fitted from the single Na2SO4 (or MgSO4) electrolyte solution, and f(Mg/Na) is the arithmetic average value.
Membranes 13 00811 g005
Figure 6. The variation of the experimental and the predicted total membrane conductivities with the cation composition of the external mixed H2SO4 + Na2SO4 electrolytes of 0.5 M total sulfate, for the 3361BW membrane (a), CSE membrane (b) and Nafion 117 membrane (c). Here, f(Na) (or f(H)) represents the volume fraction of the inter-gel phase in the membrane fitted from the single Na2SO4 (or H2SO4) electrolyte solution, and f(H/Na) is the arithmetic average value.
Figure 6. The variation of the experimental and the predicted total membrane conductivities with the cation composition of the external mixed H2SO4 + Na2SO4 electrolytes of 0.5 M total sulfate, for the 3361BW membrane (a), CSE membrane (b) and Nafion 117 membrane (c). Here, f(Na) (or f(H)) represents the volume fraction of the inter-gel phase in the membrane fitted from the single Na2SO4 (or H2SO4) electrolyte solution, and f(H/Na) is the arithmetic average value.
Membranes 13 00811 g006
Figure 7. Comparison of the total conductivity variations for the CSE membrane in mixed MgSO4 + Na2SO4 electrolytes as a function of the cation composition in the external solution. The inter-gel phase conductivity ( κ i n ) is obtained via two different approaches: calculated with (a) the ionic mobilities in infinite dilution, (b) with the inter-gel phase conductivity data simulated in the OLI Studio.
Figure 7. Comparison of the total conductivity variations for the CSE membrane in mixed MgSO4 + Na2SO4 electrolytes as a function of the cation composition in the external solution. The inter-gel phase conductivity ( κ i n ) is obtained via two different approaches: calculated with (a) the ionic mobilities in infinite dilution, (b) with the inter-gel phase conductivity data simulated in the OLI Studio.
Membranes 13 00811 g007
Figure 8. The variation of the experimental and the predicted total conductivities of (a) 3361BW, (b) CSE and (c) Nafion117 with the cation composition of the external mixed MgSO4 + Na2SO4 electrolytes of 0.5 M total sulfate. The ionic mobilities in the inter-gel phase are the same as those in the external solutions obtained from the OLI Studio data.
Figure 8. The variation of the experimental and the predicted total conductivities of (a) 3361BW, (b) CSE and (c) Nafion117 with the cation composition of the external mixed MgSO4 + Na2SO4 electrolytes of 0.5 M total sulfate. The ionic mobilities in the inter-gel phase are the same as those in the external solutions obtained from the OLI Studio data.
Membranes 13 00811 g008
Figure 9. The variation of the experimental and the predicted total conductivities of (a) 3361BW, (b) CSE and (c) Nafion117 with the cation composition of the external mixed H2SO4 + Na2SO4 electrolytes of 0.5 M total sulfate. The ionic mobilities in the inter-gel phase are the same as those in the external solutions obtained from the OLI Studio data.
Figure 9. The variation of the experimental and the predicted total conductivities of (a) 3361BW, (b) CSE and (c) Nafion117 with the cation composition of the external mixed H2SO4 + Na2SO4 electrolytes of 0.5 M total sulfate. The ionic mobilities in the inter-gel phase are the same as those in the external solutions obtained from the OLI Studio data.
Membranes 13 00811 g009
Figure 10. The correlation between the gel phase conductivity ( κ g ) of CEMs, (a) 3361BW, (b) CSE and (c) Nafion117, and the counter-ions composition in the gel phase under mixed MgSO4 + Na2SO4 electrolytes of different total sulfate concentrations (0.1, 0.5 and 1.0 M).
Figure 10. The correlation between the gel phase conductivity ( κ g ) of CEMs, (a) 3361BW, (b) CSE and (c) Nafion117, and the counter-ions composition in the gel phase under mixed MgSO4 + Na2SO4 electrolytes of different total sulfate concentrations (0.1, 0.5 and 1.0 M).
Membranes 13 00811 g010
Figure 11. The correlation between the gel phase conductivity ( κ g ) of CEMs, (a) 3361BW, (b) CSE and (c) Nafion117, and the counter-ions composition in the gel phase under mixed H2SO4 + Na2SO4 electrolytes of different total sulfate concentrations (0.1, 0.5 and 1.0 M).
Figure 11. The correlation between the gel phase conductivity ( κ g ) of CEMs, (a) 3361BW, (b) CSE and (c) Nafion117, and the counter-ions composition in the gel phase under mixed H2SO4 + Na2SO4 electrolytes of different total sulfate concentrations (0.1, 0.5 and 1.0 M).
Membranes 13 00811 g011
Figure 12. The variation of the experimental and the predicted total membrane conductivities with the cation composition of the mixed external MgSO4 + Na2SO4 electrolytes (total sulfate concentration 0.5 M), for the 3361BW membrane (a), CSE membrane (b) and Nafion 117 membrane (c). The gel phase conductivity ( κ g ) of the CEMs is obtained from the nonlinear correlation method considering the interaction parameter ( α ).
Figure 12. The variation of the experimental and the predicted total membrane conductivities with the cation composition of the mixed external MgSO4 + Na2SO4 electrolytes (total sulfate concentration 0.5 M), for the 3361BW membrane (a), CSE membrane (b) and Nafion 117 membrane (c). The gel phase conductivity ( κ g ) of the CEMs is obtained from the nonlinear correlation method considering the interaction parameter ( α ).
Membranes 13 00811 g012
Table 1. Ionic mobilities in infinite dilution at 298.15 K.
Table 1. Ionic mobilities in infinite dilution at 298.15 K.
IonsNa+Mg2+H+SO42−
u i × 10 8 (m2s−1V−1)5.205.5036.28.27
Table 2. Fundamental properties of commercial cation exchange membranes.
Table 2. Fundamental properties of commercial cation exchange membranes.
Mem.Area
Resistance (Ω cm2)
Density of Dry Polymer
(g cm−3)
Dry Membrane Thickness
(μm)
Water Uptake
(g g−1 Dry Polymer)
IEC
(meq g−1 Dry Polymer)
CSE1.81.16 ± 0.031600.29 ± 0.052.34 ± 0.01
3361BW11.00.93 ± 0.014200.43 ± 0.082.45 ± 0.01
Nafion1171.01.94 ± 0.111800.19 ± 0.010.99 ± 0.01
Table 3. The fitting results of the nonlinear correlation between κ T and κ s according to Equation (1).
Table 3. The fitting results of the nonlinear correlation between κ T and κ s according to Equation (1).
MembranesElectrolytesα [-] f 2 [-]R2 [-]
3361BWNa2SO40.0250.3050.910
MgSO40.2810.3170.868
H2SO4−1.0720.8180.987
CSENa2SO4−0.0060.1680.852
MgSO4−0.7320.9980.924
H2SO4−1.5610.9990.959
Nafion117Na2SO40.4270.1280.971
MgSO41.1800.0580.833
H2SO4−0.7340.1880.750
Table 4. The fitting results of the nonlinear and linear correlations between κ T and κ i n , and the κ i n is calculated from the ionic mobility under infinite dilution (Equation (8)).
Table 4. The fitting results of the nonlinear and linear correlations between κ T and κ i n , and the κ i n is calculated from the ionic mobility under infinite dilution (Equation (8)).
MembranesElectrolytesNonlinearLinear
α [-] f 2 [-]R2 [-] f 2 [-]R2 [-]
3361BWNa2SO4−0.0550.2070.9160.2100.931
MgSO40.2130.1670.7400.1570.801
H2SO4−1.3930.1550.9990.3970.932
CSENa2SO4−0.2820.1410.8680.1530.863
MgSO4−0.3820.0990.8680.0980.879
H2SO4−2.4390.0050.8840.1500.610
Nafion117Na2SO40.2560.1200.9350.0980.931
MgSO40.9400.0860.7850.1250.695
H2SO4−1.2570.0140.7720.0680.581
Table 5. The fitting results of the nonlinear and linear correlations between κ T and κ i n , and the κ i n is obtained from the ionic mobility the same as the external solution (Equation (12)).
Table 5. The fitting results of the nonlinear and linear correlations between κ T and κ i n , and the κ i n is obtained from the ionic mobility the same as the external solution (Equation (12)).
MembranesElectrolytesNonlinearLinear
α [-] f 2 [-]R2 [-] f 2 [-]R2 [-]
3361BWNa2SO40.1620.3100.9190.2690.932
MgSO40.4000.3190.7110.1780.769
H2SO4---0.4400.924
CSENa2SO4−0.3560.1520.8660.2180.858
MgSO4−0.3790.0740.8530.1180.876
H2SO4−1.7030.0060.7950.1590.588
Nafion117Na2SO40.4770.2390.9240.1170.908
MgSO4---0.1630.640
H2SO4−1.8660.0020.7540.0720.569
Table 6. Mobilities of cations in the gel phase obtained with two different sets of fitting results.
Table 6. Mobilities of cations in the gel phase obtained with two different sets of fitting results.
Mem. u N a + g × 105
(cm2 s−1V−1)
u M g 2 + g × 105
(cm2 s−1V−1)
u H + g × 105
(cm2 s−1V−1)
Infinite DilutionOLI SimulationInfinite DilutionOLI SimulationInfinite DilutionOLI Simulation
CSE4.935.149.321.4033.3731.23
3361BW10.855.650.960.4196.0434.60
Nafion11713.2413.521.551.9559.0057.30
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wu, L.; Jiang, H.; Luo, T.; Wang, X. On the Ionic Conductivity of Cation Exchange Membranes in Mixed Sulfates Using the Two-Phase Model. Membranes 2023, 13, 811. https://doi.org/10.3390/membranes13100811

AMA Style

Wu L, Jiang H, Luo T, Wang X. On the Ionic Conductivity of Cation Exchange Membranes in Mixed Sulfates Using the Two-Phase Model. Membranes. 2023; 13(10):811. https://doi.org/10.3390/membranes13100811

Chicago/Turabian Style

Wu, Liansheng, Haodong Jiang, Tao Luo, and Xinlong Wang. 2023. "On the Ionic Conductivity of Cation Exchange Membranes in Mixed Sulfates Using the Two-Phase Model" Membranes 13, no. 10: 811. https://doi.org/10.3390/membranes13100811

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