Maxwell–Stefan model of multicomponent ion transport inside a monolayer Nafion membrane for intensified chlor-alkali electrolysis

A mathematical model based on a generalized Maxwell–Stefan equation has been developed to describe multicomponent ion and water transport inside a cation-exchange membrane. This model has been validated using experimental data and has been used to predict concentration profiles, membrane potential drop, and transport numbers of ions and water for the chlor-alkali process at increased current densities. Several improvements have been made to previously developed Maxwell– Stefan models. In our model, the generalized Maxwell–Stefan equation is written in terms of concentration instead of mole fraction and the fixed group (membrane) concentration is assumed to be constant. We have adapted the Augmented matrix method using the built-in partial differential equation parabolic elliptic (pdepe) solver in Matlab®, and both the concentration and the electrical potential gradients have been solved simultaneously. The boundary conditions are determined with the Donnan equilibrium at the membrane–solution interface. We have also employed semi-empirical correlations to define the Maxwell–Stefan diffusivities inside the membrane. For the bulk diffusivities, we applied the correlations for the concentrated solution instead of the values at infinite dilution. With the diffusivities presented in this work, the model shows a better fit to the experimental data than with previously reported fitted diffusivities. Prediction of the sodium transport number and water transport number is generally good, whereas the deviations with regard to membrane potential might also be related to issues with the experimental data. The model predicts an increase in both sodium and water transport numbers at increased current density operation of chlor-alkali production.


Introduction
Renewable energy sources are increasingly available as an alternative electricity source for large-scale electrolysis processes, such as chlorine production and water electrolysis. Due to the intermittent nature of the renewables, it is crucial to develop new intensified electrolysis cells that can be operated with high flexibility. This can be achieved by increasing the maximum current density. For the chlor-alkali process, an ambitious target would be 30 kA m −2 (currently the process is limited to 7 kA m −2 ) [1].
A key aspect of operation at higher current densities is the membrane performance. Unfortunately, the information on membrane performance at high current densities remains scarce. Ion transport across a membrane under current load is not completely understood despite the fact that there have been several attempts to model the behavior by either using Maxwell-Stefan (MS) or Nernst-Planck (NP) models [2][3][4][5][6][7][8]. The Nernst-Planck approach assumes an ideal solution and neglects ion-ion interactions. This model is known to be valid for dilute ionic systems [9,10], but chlor-alkali electrolysis involves highly concentrated solutions, typically around 5 M sodium chloride (NaCl) and 10 M sodium hydroxide (NaOH) [11]. In this case, the Maxwell-Stefan is considered more reliable since the interactions of different components and the non-ideal solutions are taken into account [9,[12][13][14]. Moreover, the Maxwell-Stefan approach includes the water transport via the solvent-ion interactions, whereas the Nernst-Planck model has to introduce a separate equation (i.e., the Schlögl equation) to account for the water transport [10,15].
The process intensification of the chlor-alkali process would benefit from a mathematical model that can predict the membrane performance. Up to this date, only Van der Stegen et al. [2] have developed a Maxwell-Stefan model for the chlor-alkali system. However, this model has some limitations. It was derived in mole fraction which assumed that the total concentration is known. The ionic fixed groups of the membrane were regarded as one of the mobile components in the aqueous mixture to obtain convergence, whereas the fixed charge group concentration is generally considered to be constant inside the membrane. This model has also simplified the calculation of the membrane potential gradient by neglecting the concentration gradient and by using Ohm's law to derive the potential gradient explicitly. The neutrality condition is broken by this simplification, which has been numerically proven during the investigation of the extended Nernst-Planck model [8].
The main drawback in applying the Maxwell-Stefan approach is the lack of reliable data on diffusivities at high concentrations. The Maxwell-Stefan diffusivities ( i, j ) are required for the interaction between the components in a mixture. In a mixture with n components, the number of Maxwell-Stefan diffusivities should be 0.5 × n × (n − 1) based on the Onsager relations [16]. The chlor-alkali system contains at least five components (Na + , Cl − , OH − , H 2 O, and -SO 3 − ) and hence requires ten binary diffusivities. Wesselingh et al. [17] proposed that the diffusivities inside the membrane can be related to the diffusivities in the bulk using the tortuosity factor (τ), (Eq. 1). 0 i, w is the diffusion coefficient in infinitely diluted aqueous solution (Eq. 2). The values for 0 i, w at 25 °C and 90 °C are given in Table 1 for sodium, chloride, and hydroxide in water. Kraaijeveld et al. [4] proposed a correlation for the diffusivity of positive ions and sulfonate groups in the membrane ( m in Eq. (3). As can be seen from the equation, it is suggested that this diffusivity is related to the diffusivity of the same ion with water in the membrane. The diffusivity of negative ions and sulfonate groups in the membrane were fitted in the model based on the experimental data of dialysis.
According to Van der Stegen et al. [2], the low values of diffusivities from the correlations of Wesselingh et al. resulted in an unreasonably high membrane potential for the chloralkali system. Instead they opted for estimating diffusivities using a sensitivity analysis based on four output parameters in the chlor-alkali process: current efficiency, cell potential, relative water transport number, and the chloride concentration in caustic (catholyte solution). They assumed that the membrane potential was about 15% of the total cell potential. The input parameters for the model were anolyte (NaCl) and catholyte (NaOH) concentration, membrane thickness, membrane equivalent weight (EW), current density, temperature and initial values of diffusivities. They excluded four out of ten Maxwell-Stefan diffusivities by choosing high values (1 × 10 −8 m 2 s −1 ) to neglect the interaction between these components. They indicated that the values of the diffusivities are a function of current density (I) as given in Eq.  Table 1.
Other groups have looked at diffusivities in ion-exchange membranes for other systems than chlor-alkali. Visser et al. [3] encountered a similar problem with the low values of diffusivities using the semi-empirical equations of Wesselingh et al. Unlike Van der Stegen et al., they investigated the interactions of different electrolytes: HCl, H 2 SO 4 , NaCl, NaOH, and Na 2 SO 4 . The diffusivities were fitted from several partial diffusion experiments: diffusion dialysis (salt diffusion flux and osmotic water flux), electro-osmotic, membrane resistance, pressure driven volume flow, and electrodialysis. They performed in total 26 experiments at 25 °C for a Nafion 450 membrane to define 21 binary diffusivities. The experiments used current densities from 0 to 1 kA m −2 . The diffusivity values for NaCl and NaOH in the membrane are listed in Table 1. Despite their concern about low diffusivity issues from the semi-empirical equation, they applied Eq. (1) to estimate the chloride and hydroxide interaction with water inside the membrane ( m OH − , w , m Cl − , w ) . It is observed in Table 1 that the value of the binary diffusivities of sodium and hydroxide ( m Na + , OH − = 100 × 10 −10 m 2 s −1 ) is remarkably different from the diffusivities of sodium and chloride ( m Na + , Cl − = 0.580 × 10 −10 m 2 s −1 ) , for which no explanation was given. It is important to note that they have excluded the Donnan potential as the boundary condition because of the convergence issue in the model. The co-ion concentration inside the membrane was estimated based on the algebraic relations as a function of the external composition instead of the Donnan potential. Chapman et al. [18] investigated Maxwell-Stefan diffusivities for concentrated electrolyte systems and made correlations to relate the diffusivities for high concentrations to diffusivities at infinite dilution as given in Eqs. (5)(6)(7)(8). The interaction of negative and positive ions increases substantially with increasing concentration as shown in Eq. (9). The value of q in Eq. (9) is defined in Eq. (10). C is the concentration of the negative ion and C o represents the bulk concentration. G (in K 3/2 M −1/2 ) is considered to be a correction for the idealized theory. It is strongly related to the concentration and the type of electrolyte. There was no systematic trend found for the effect of the temperature. Sodium chloride of 5 M at 50 °C has a G value of 979 K 3/2 M −1/2 . The concentration of sodium hydroxide investigated is limited to 1.5 M with the temperature at 25 °C, which has the value of 2606 K 3/2 M −1/2 . For the high concentration above 3 M, the average value for uni-univalent electrolyte is reported to be 3000 ± 1000 K 3/2 M −1/2 .
This work aims to improve the Maxwell-Stefan modeling for the chlor-alkali process and to increase its applicability to higher current densities. In contrast with the work of Van der Stegen, the ionic fixed group (membrane) concentration is defined based on the known membrane properties. As proposed by Krishna [19], both the concentration (chemical potential) gradient and the electrical potential gradient can be calculated simultaneously using an augmented matrix method. By adopting this method, no further assumption about the potential gradient is needed. This paper also investigates the influence of the Maxwell-Stefan diffusivities on the membrane performance in OH − , w (NaOH) = 0 Na + , w + 3.747 × 10 −6 × C 1∕2 NaOH + 7.090 × 10 −5 × C NaOH + 1.708 × 10 −4 × C 3∕2 NaOH + 9.499 × 10 −5 × C 2 Not available terms of membrane permselectivity (sodium transport number), membrane potential, and relative water transport number. It compares the fitted values of diffusivities with the more general correlations which are applicable for different operating conditions such as concentration, temperature, and current density. The model is then validated with the experimental data reported in the literature for the chlor-alkali system at 2 and 3 kA m −2 [20][21][22][23]. Lastly, the model predicts the membrane performance in terms of membrane permselectivity and membrane potential for high current densities up to 30 kA m −2 .

Maxwell-Stefan equation
Maxwell-Stefan theory is a steady-state force balance between driving forces and friction forces acting on a certain component in the mixture. Equation (11) presents the relation between the driving force on a component i in the mixture and the sum of the friction forces between i and the other component j in terms of mole fraction [9].
The derivation of the driving forces of the generalized Maxwell-Stefan is based on irreversible thermodynamics [9,24]. If the system contains ions, the partial molar Gibbs free energy depends not only on the chemical potential (µ i ), but also on the electrical potential (ϕ). The combination of these potentials is called the electrochemical potential (η i ) as given in Eq. (12) [25].
where  is the Faraday constant, V i is the partial molar volume of the solvent, P is the pressure, R is the universal gas constant, T is the temperature, and z i is ionic charge of component i. The activity of component i (a i ) is defined by Eq. (13) using the activity coefficient i to account for the non-ideal solution.
The contribution of the pressure gradient is negligible in electrochemical cells when compared to the concentration and the electrical potential gradients [2,9,26]. For an ideal system, the generalized Maxwell-Stefan equation from Eq. (11) can be written in terms of the concentration (Eq. 14). The relation between the ionic fluxes and the current density is shown in Eq. (15). Hence, the electroneutrality condition needs to be met according to Eq. (16). (11)

Input parameters
The input parameters used in the model are presented in Table 2. The generalized Maxwell-Stefan equation in Eq. (14) is derived for an ideal solution with constant pressure and temperature. The concentration of the stationary fixed charges per void fraction in the membrane is calculated by Eq. (17) [8,11,27]. The EW is defined as the dry weight of polymer in gram per mole of sulfonic acid groups. The ratio of weight fraction of fixed ionic groups and electrolyte in the ion cluster (f m ∕f e ) is around 0.4-1.32 depending on the membrane properties and electrolyte concentration [11]. The density of the sodium hydroxide solution in equilibrium with the membrane is used to represent the density of the electrolyte adsorbed in the membrane (ρ e ).
The water uptake (W w ) in the membrane depends on the EW and it is a function of the sodium hydroxide concentration. The water content decreases with increasing EW. The correlations for the water uptake in weight percentage of dry polymer as a function of sodium hydroxide up to 10 M are given in Eqs. (18,19), for sulfonate EW1100 and sulfonate EW1200 [11].

Boundary conditions
In this paper, we focus on the investigation of the mass transfer behavior inside the membrane. Therefore, the mass transfer resistance is only considered inside the membrane by assuming that a high mass transfer takes place in the bulk solution. The Donnan equilibrium theory is applied to define the concentration of ions at the interface ( Table 3). The Donnan equilibrium for all ions at the membrane surfaces is expressed in Eq. (31) using the same distribution ratio (K) shown in Eq. (32) [8,27,29]. Figure 1 depicts the concentration jump of the ionic species at the solution and the membrane interface for both anolyte and catholyte sides. The ion concentration at the solution interface is assumed to be the same as that of the bulk concentration.
The concentration of water is calculated using the density correlations given in Eqs. (20)(21)(22)(23)(24)(25)(26)(27)(28)(29)(30). The density correlation is available for the mixed electrolyte of sodium hydroxide and sodium chloride but not for the sulfonate group of the membrane. Considering the simplification of the model, the total sodium ion, calculated from the Donnan equilibrium at the left boundary, is used to define the weight percentage 1 3 of sodium chloride as presented in Eq. (33). Similarly, the total sodium ion at the right boundary determines the weight percentage of sodium hydroxide shown in Eq. (34).

Augmented matrix
The Maxwell-Stefan equation given in Eq. (14) forms a system of non-linear differential algebraic equations (DAEs). The built-in pdepe solver in Matlab® can solve DAEs using the ordinary differential equations (ODE15s) solver, only applicable for DAE index 1 or index 0. The index of DAEs is defined from the number of differentiations needed to reduce DAEs to ODEs. Equation (14) consists of both concentration and potential gradients (DAE index 2), which cannot be solved by the built-in pdepe solver in Matlab. The DAE index 2 is reduced to index 1 by applying the augmented matrix method as proposed by Krishna et al. [19] shown in Eq. (35).
The non-linear flux equations need to be arranged in a matrix format (Eqs. [36][37][38][39][40]. The derivation of augmented formulation is explained in Appendix. The concentration gradient (dC i ∕dz) as driving force in (b i ) is a linear combination of both flux (N i ) and the potential gradient (d ∕dz).   We apply Eq.  correlation is available for the diffusion coefficient of negative ion-negative ion ( m OH − , Cl − ) and a base case value of 1.0 × 10 −10 m 2 s −1 is used in this model. The calculated values for three different temperatures are listed in Table 5. Table 5 contains four sets of Maxwell-Stefan diffusivities used in the model simulation. It should be noted that the fitted values of diffusivities obtained by Visser et al. were used at 25 °C and a sodium concentration up to 4 M. Therefore, this set of diffusivities is less applicable for the chloralkali system. However, due to very limited availability of data for Maxwell-Stefan diffusivities, these values are used in the simulation for comparison. Figure 2 depicts the concentration profiles of sodium, hydroxide, chloride, and water inside the membrane at 2 kA m −2 using different values of Maxwell-Stefan diffusivities as listed in Table 5. The conditions are 80 °C, EW = 1150, membrane thickness = 0.25 mm, 25 wt% NaCl as anolyte, and 32 wt% NaOH as catholyte similar to the experimental work by Yeager [20,21]. All sets of Maxwell-Stefan diffusivities generate non-linear concentration profiles for each component. The concentration profiles of the ions using diffusivities from Visser, Van der Stegen, and this work show similar trends, but clearly differ from the base case. The water concentration profile is similar for the diffusivities from Visser and this work, but is clearly different for the base case and van der Stegen. The profiles of the potential gradient in Fig. 3a also show a non-linear behavior and are depending on the values of diffusivities. This confirms the  Despite the clearly different concentration profiles for the different sets of diffusivities, all models predict similar pH profiles inside the membrane due to the fact that pH is logarithmic (Fig. 3b). This demonstrates that the hydroxide ion concentration is relatively high, which results in a pH above 12 throughout the membrane. Similar results have been reported in the previously developed models    [7,31]. It should be mentioned that in our model we use pKw = 12.60 at 80 °C [28] instead of the pKw = 14 at 25 °C.

Model validation for low current densities
The membrane permselectivity (sodium transport number) and the membrane potential are key performance criteria in the chlor-alkali process. To validate the models with the different diffusivity sets, the model was run with the same conditions as experiments by Yeager et al. and Berzins [20,23] in which the catholyte concentration was varied from 10 to 35 wt% (≈ 2.5-12 M). Figure 4a shows the sodium transport number as a function of the catholyte concentration at 2 kA m −2 using a monolayer sulfonate membrane with EW = 1150, membrane thickness = 0.25 mm, and temperature = 80 °C. Similar profiles are also observed in Fig. 4b for a monolayer membrane with EW = 1100 and a membrane thickness of 0.1 mm at 3 kA m −2 . It is observed that the influence of the diffusivities on the calculated sodium transport number is even more pronounced than for the concentration profiles. The base case diffusivities overpredict the sodium transport number, while the opposite is observed for both the Visser and Van der Stegen cases. The diffusivities based on the empirical correlation used in our model show the best fit to the experimental work. Sensitivity analysis was performed for diffusion coefficient of negative ion-negative ion ( m OH − , Cl − ) for the range 0.1 × 10 −10 − 50 × 10 −10 m 2 s −1 , and no significant impact was found on the ion transport number and the membrane potential for this system. It could be because the chloride ion flux through the membrane is negligible. Therefore, a base case value of 1.0 × 10 −10 m 2 s −1 is used throughout the simulation.
Water is also transported together with sodium to the catholyte site. The Maxwell-Stefan model includes the water transport via water-ion interactions. Therefore, a separate semi-empirical equation (i.e., the Schlögl equation) is no longer required. Figure 5 shows water transport calculated with the different models and compares these to experimental data. It can be seen that both Van der Stegen and Visser predict too low water transport, especially at high caustic concentrations, whereas the base case and our work show a good match. It is important to note that water concentration at the membrane interface is calculated based on the density correlations of sodium chloride and sodium hydroxide, which does not completely reflect reality since it does not take into account the density difference compared to sulfonate groups in the membrane.
The performance of the models can be explained by considering the diffusivities. Visser uses a very high binary diffusivity of sodium ion-hydroxide ion inside the membrane ( m Na + , OH − = 100 × 10 −10 m 2 s −1 ) , which results in a high hydroxide ion transport inside the membrane and a low sodium transport number. The diffusivity exceeds the value calculated using the correlation suggested by Chapman for 10 M (≈ 32 wt%) NaOH ( bulk Na + , OH − = 7.5 × 10 −10 m 2 s −1 at 25 °C). The diffusivity values suggested by Van der Stegen result in even lower sodium transport numbers and negative water transport numbers. In this case, it does not seem logical that the diffusivities of chloride inside the membrane are higher than the value for infinite dilution as listed in Table 1. The diffusivities presented in this work seem more realistic and this is reflected in both sodium and relative water transport numbers.
The values of diffusivities also affect the membrane potential drop as shown in Fig. 6. The suggested values  Table 5 (see legend). EW = 1150, membrane thickness = 0.25 mm, temperature = 80 °C, 25 wt% NaCl, and 32 wt% NaOH 1 3 of diffusivities by Van der Stegen as a function of current density as given in Eq. (4) show an unrealistic profile. The model was run with the operating conditions reported by Bergner et al. [32] with a typical thickness of 0.29 mm for Nafion N954 [33] and EW = 1100. Our model shows a reasonable match with the experimental values of the membrane potential at 80 °C around 0.51 V at 6 kA m −2 [34] and 0.291 V at 3 kA m −2 [11]. However, the linear correlations proposed by Bergner et al. for both their experimental data and the experimental data of Nidola [35] show higher values than the predicted model at lower current densities. The fact that they still show a significant membrane potential at zero current density is not easy to explain. Osmotic pressure differences between the anolyte and catholyte could play a role, but it might also be related to experimental issues. They mentioned that during the membrane potential measurement, two Luggin capillaries could not be placed with sufficient accuracy (the zero gap configuration had a gap of 0.5-1 mm). Chandran et al. [36] also reported similar experimental difficulties, resulting in large measured membrane potentials.  [20,21] and Berzins [23]) of the sodium transport number as a function of catholyte concentration. For the models, the different values of Maxwell-Stefan diffusivities as listed in Table 5 [20,21]) of the relative water transport number as a function of catholyte concentration using different values of Maxwell-Stefan diffusivities as listed in Table 5 (see legend) at 2 kA m −2 . EW = 1150, membrane thickness = 0.25 mm, temperature = 80 °C and 25 wt% NaCl  Table 5. The experimental data of Bergner et al. [32]: NaCl = 18 wt%, NaOH = 33 wt%, temperature = 90 °C, EW = 1100, membrane thickness = 0.29 mm (properties of Nafion N954 [33]). The typical values of membrane potential at 80 °C for current densities of 3.5 kA m −2 and 6 kA m −2 are 0.291 V [11] and 0.51 V [34], respectively All in all one can conclude that the model with the diffusivities presented in this work shows the best fit compared to experimental data. Prediction of the sodium transport number and water transport number is good, whereas the deviations with regard to membrane potential might also be related to issues with the experimental data. Given this good performance, we can now consider using the model to predict membrane performance at conditions for which experimental data are not yet available such as high current densities.

Model prediction for high current densities
The Maxwell-Stefan diffusivities suggested in this work are applied for further simulation. The model uses the operating conditions of the chlor-alkali process: 25 wt% NaCl as anolyte, 32 wt% NaOH as catholyte, and a Nafion single-layer membrane, N-1110 (EW = 1100 and thickness = 0.27 mm) [8,33]. Figure 7 presents the transport number of sodium, hydroxide, water, and membrane potential drop as a function of current density from 2 to 30 kA m −2 for low and high temperature (25 °C and 90 °C). The sodium transport number shows an increasing trend with increasing current density, which is also experimentally observed in chlor-alkali [11]. This can be explained by the fact that at low current densities, diffusion is still important compared to migration. This leads to a higher back diffusion of sodium from catholyte (around 10 M NaOH) to anolyte (around 5 M NaCl). At higher current densities, the contribution of diffusion diminishes compared to the increasing migration term. Increased temperature results in decreased membrane permselectivity, which can be explained by the fact that higher temperature leads to higher mobility of ions and enhances ion diffusion. The model also predicts that the chloride ion transport is negligible for both high current density and temperature.
As shown in Fig. 7, the membrane potential is significantly lower for higher temperature. This is related to the increased conductivity of the membrane. In the model, the conductivity is not explicitly present as a variable but it is taken into account in the diffusivities as shown in Eq. (2). The infinite dilution of diffusivities is related to the limiting ionic conductivity. The ionic conductivity increases with increasing temperature which results in a lower resistance, thus lower membrane potential drop. Figure 8 illustrates the effects of the fixed ionic group concentration and the current densities on the membrane permselectivity and the membrane potential. The model predicts a higher sodium transport number with higher fixed ionic group concentration and current densities. The latter is in line with the results shown in Fig. 7. Higher sulfonate groups in the membrane prevents the transport of the hydroxide ion as co-ion based on the Donnan exclusion and enhances the sodium ion transport as counter-ion. The water concentration decreases with increasing sulfonate groups concentration, and the model predicts a lower water Fig. 7 The transport number of sodium, hydroxide, water, and membrane potential as a function of current density using the Maxwell-Stefan diffusivities based on this work for two different temperatures (see legend). The model simulation used the Nafion 1110 with properties: EW = 1100, membrane thickness = 0.27 mm, 25 wt% NaCl as anolyte, and 32 wt% NaOH as catholyte transport for higher fixed ionic group concentration. A slight increase is observed in the membrane potential and this can be related to the decreased water concentration for higher membrane concentration.

Conclusion
A Maxwell-Stefan model has been developed to investigate the non-linear behavior of multicomponent ion and water transport inside a cation-exchange membrane. The non-linear concentration profiles, the membrane potential drop, and the transport number of ions and water strongly depend on the values of the Maxwell-Stefan diffusivities. To conclude, the combined correlations proposed by Wesselingh et al., Kraaijveld et al., and Chapman et al. show a good agreement with the available experimental data of chlor-alkali for sodium transport number and relative water transport number. Thus, it is our considered opinion that the semi-empirical correlations are suitable for defining the Maxwell-Stefan diffusivities and therefore are suitable for further simulations.
The improvement of our model compared to the previously developed models is the ability to calculate both fluxes and the membrane potential drop simultaneously by adopting the augmented matrix method. A typical expected membrane potential for chlor-alkali process around 0.5 V at 6 kA m −2 and 80 °C is well predicted by our model. OpenAccess This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creat iveco mmons .org/licen ses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Appendix: Augmented matrix formulation
The fixed charged groups of the cation-selective membrane are treated as a component in the aqueous electrolyte solution. The chlor-alkali process contains in total 5 components including the fixed charged groups (n = 5). The flux of fixed charged groups of the membrane is zero. This results in a matrix of (n − 1) fluxes. In addition, the concentration of the fixed charged groups is constant. The fixed charged group is therefore set as the nth (5th) component.
Every component is defined a number: The concentration of fixed charged group is known from the membrane property. The total concentration is the sum of the concentration of every component: (47) Fig. 8 The transport number of sodium, hydroxide, water, and membrane potential as a function of ionic fixed group concentration using the Maxwell-Stefan diffusivities suggested in this work for three dif-ferent current densities (see legend). Membrane thickness = 0.27 mm, temperature = 90 °C, NaCl as anolyte = 25 wt%, and NaOH as catholyte = 32 wt% Based on the literature written by Krishna (1987), the DAE index 2 can be reduced to index 1 by applying the augmented matrix method [19].
Except in the region close the electrode surface, where there occurs charge separation due to double-layer phenomena, the electroneutrality condition needs to be met. This means that there is no net electrical body force acting on the mixture as a whole: Since the water does not carry either positive or negative charge (z 3 = 0) , one concentration variable can be eliminated: The concentration of the fourth component can therefore be calculated: Both the flux and the driving force of the fixed charge groups are zero. This leads to (n − 1) equations: (51) There are in total four equations with five unknowns (N 1 , N 2 , N 3 , N 4 and d ∕dx).
The current density can be defined in terms of flux: The fixed charged groups of the cation-exchange membrane contain zero flux (N 5 = N SO − 3 = 0) because these ions are not transported (kept in place).
Using the Eq. (53), there are now in total five equations with five unknowns: (52)  The augmented vector of the driving forces: The augmented vector of the unknown variables: The fluxes and the electro potential gradient can be therefore calculated: The negative sign is added for the physical meaning of the fluxes to be a positive value (driving forces (b) are positive, while the diagonal matrix [B] contains negative values). (59)