1. Introduction
Salinity gradient energy (SGE) can be generated from the reversible mixing of river water and seawater. It is potentially a clean and sustainable form of energy, producing brackish water and electrical energy only. The theoretical energy content of mixing 1 m
3 river water with a large surplus of is 2.5 MJ or 1.7 MJ when mixed with 1 m
3 [
1]. The global potential of SGE is estimated to be 2.6 TW [
2] which is comparable with the current world electricity consumption of about 3 TW [
3] The two main membrane-based technologies that can be used to harvest this energy are reverse electrodialysis (RED) and pressure retarded osmosis (PRO) [
1]. This paper is focused on RED.
A RED stack consists of a large number of cells, each containing an anion exchange membrane (AEM), a seawater compartment, a cation exchange membrane (CEM) and a river water compartment. The stack is terminated on each side with an electrode compartment and an electrode. On one side, there is an extra ion exchange membrane (IEM), the same type as on the other side (
Figure 1). The driving force in the cell is the chemical potential difference between the sea- and river water. Positive ions diffuse through the CEM from the seawater compartment to the river water compartment in one direction; the negative ions through the AEM in the other direction. The internal ionic current is converted to an external electron current on the electrodes by proper redox reactions.
The feed consists of two solutions with different salt content: i.e., seawater with river water, Red Sea water with Dead Sea water, an industrial brine waste with seawater and so on. In this paper, we focus on the combination of seawater with river water.
Table 1 shows that the main ions in seawater are Na
+ and Cl
− and we restrict our model to these ions.
The feasibility of RED is described by two response parameters: power density P
d (power per m
2 total active membrane) and energy efficiency Y (the generated amount of energy as part of the exergy from the inlets). In earlier experiments, we studied the effect of different parameters on these response parameters [
5]. Experimental power densities are typical 1.2 W/m
2. If operated under conditions of maximal power density, the energy efficiency is about 20%, increasing to more than 85% with low power density and using circulating feed waters [
6]. Membrane properties are decisive for optimal operating and influence power density as well as fuel efficiency by osmosis, electro-osmosis, co-ion transport, monovalent/multivalent competition and electrical as well hydrodynamic losses.
There are two main contributors to the performance of a RED stack: the feed water supply through the spacer filled compartments (parallel to the membranes) and the ion and water transport through the membranes (perpendicular to the membranes). The first movement is described by hydrodynamics; the driving force is the pressure difference as applied by the pumps. The driving forces for ion and water transport are the chemical and electrical potentials; resistance against the particle movements is caused by friction between the different particles and between particles and immobilized charges within the membranes. This complicated process is described by the Maxwell–Stefan theory and is the subject of this paper [
7,
8].
A membrane has primary and secondary properties. Primary properties can be seen as the causes and include type (homogeneous or heterogeneous), ion exchange capacity, swelling degree, composition, degree of cross-linking, etc. Secondary properties are the consequences and examples are electrical resistance, counter ion permselectivity, water permeability, permselectivity for monovalent ions, density and tensile strength.
- -
Electrical resistance: With pure NaCl solutions, the RED process is clear. A low electrical resistance is important for high power density; high permselectivity together with low osmosis is necessary for good energy efficiency. The areal resistance of common membranes as used in RED is about 2–3 Ω∙cm
2 for CEMs and 1–3 Ω∙cm
2 for AEMs [
9].
- -
Permselectivity for counter-ions: A second membrane parameter is the permselectivity for counter-ions. A CEM is also slightly permeable for chloride ions and an AEM for sodium ions. This transport of co-ions affects not only the generated voltage (and therewith power density), but also fuel efficiency [
10]. Most common IEMs as used in RED have permselectivities of 95–100% for CEMs and 90–95% for AEMs [
9].
- -
Water permeability: Besides ions, also water can be transported through the membrane. If the RED stack is not connected to an electrical load, the ion flow through the membranes is low and the water will flow from the river water compartment to the seawater side. In fact, the water flows from high to low water concentration— a phenomenon called osmosis. However, if the stack produces power, there is an ionic flow inside the stack from the seawater to the river water compartments. Ions in aqueous solution are hydrated with a hydration number of 5.4–5.5 for Na
+ and 7.1–7.3 for Cl
−, depending slightly on the concentration [
11]. These water molecules—in direct contact with the ion—form the inner hydration shell. The moving hydrated ions can drag even more water molecules that are situated in the so-called outer hydration shell. All attached water causes a water flow opposite the osmotic flow, a phenomenon known as electro-osmosis. Osmosis and electro-osmosis are opposite in direction and can even completely cancel each other out at a certain electrical current density [
10]. In normal ED, as used for desalination, osmosis and electro-osmosis have the same direction and here too this water transport is detrimental to the energy demand of the process.
Permselectivity for monovalent ions: Most commercial membranes do not discriminate between monovalent and divalent ions. However, it was found that divalent ions in the feed waters have a negative influence on the power density of a RED stack equipped with these membranes [
12,
13,
14,
15]. The main reason that divalent ions adversely affect the operation of a RED system is their low membrane potential E. This is shown by the Donnan equation for a CEM, as follows:
where
aCEM stands for the permselectivity of the CEM,
R for the gas constant (8.314 J mol
−1K
−1),
T for the temperature (K),
z for the valence of the ion and
F for the Faraday constant (96,485 C mol
−1); activities are denoted by
a, activity coefficients by
γ and concentrations (mol/L) by
C. Subscripts
S and
R are used for sea- and river water. A similar equation holds for the membrane potential over an AEM. The membrane potential is determined by the monovalent ions (
z = 1) because those are in large excess and amounts to 78 mV for pure NaCl solutions of 30 g/L and 1 g/L [
14]. Supposing that the membrane is part of a RED stack operating at maximum power, the voltage will drop to half, namely 39 mV. Divalent ions (
z = 2) generate half of the potential of monovalent ions—assuming the same ratio of activities—and this potential is equal to that 39 mV. In that case, those ions experience no net electrical force and therefore remain at rest and do not contribute to the power output. If there are relatively more divalent ions in the river water, then there will be an uphill transport of these ions. It is clear that monovalent and divalent ions have different velocities within the membrane. A measure to mitigate the effect of divalent ions is the application of mono-selective membranes. However, these membranes have a higher resistance and ultimately yield little or no gain in power density efficiency or efficiency.
To get grip on these complex systems, Higa et al. [
16,
17] applied the Nernst–Planck equation to such systems. However, this method takes no account of the mutual interactions between the different particles. In the next decade, many articles were published on the multicomponent problem in ion exchange membranes using the Maxwell–Stefan (MS) theory [
7,
8,
18,
19,
20]. This theory takes all interactions between ions, water molecules and fixed membrane ions into account. Movement of particles is derived from driving forces and from friction with other particles. In a RED stack—operating on pure NaCl solutions—the driving forces are the chemical and electrical potential while the process is slowed down by friction between H
2O–Na
+, H
2O–Cl
−, Na
+–Cl
−, Na
+–~SO
3−, Cl
−–~NR
3+, H
2O–~SO
3−, H
2O–~NR
3+, Na
+–~NR
3+, and Cl
−–~SO
3− (~NR
3+ and ~SO
3− are fixed charges in the AEM and CEM).The secondary properties resistance and permselectivity are dependent on these friction coefficients.
To describe the drag between particles, the MS theory uses diffusivities, inverse friction coefficients. In a simplified model of seawater as shown in
Table 1, there are 8 ions; the MS approach of such a system in RED uses a total of 11 particles: the 8 mentioned for the solved ions, 1 for water and 2 for the two types of fixed charges on the membranes (AEM and CEM). All diffusivities can be written as a single 11 × 11 matrix. However, the diffusivities on the main diagonal are not interesting and the matrix is symmetrical, and therefore only 55 elements remain. Because there is no interaction between the fixed ions on the AEM and the CEM, the concerning element can be omitted, remaining 54 elements. In general—following this reasoning—with
n different ions, the number of MS interactions (
NMS) is as follows:
Each interaction has an associated diffusivity and with multiple ions NMS increases rapidly. Thus for pure NaCl solutions, NMS = 9 and for mixtures of NaCl and MgSO4, NMS = 40. MS-diffusivities are in most cases barely obtained. Only the system consisting of Na+, Cl− and water in membranes is well described, but it appears that the diffusivities are dependent on concentration and type of membrane. Furthermore, published data of other ions are exceptionally scare. Therefore our model and experiments are restricted to NaCl solutions only and transport phenomena through compartments and diffusion layers will not be taken in account. In this case only 2 ions remain, resulting in 9 diffusivities. The model in this paper is validated by data achieved experimentally; the transport of ions and water and electrical properties are measured and compared with the values of the model.
3. Application of the Maxwell–Stefan Theory
The transport through a membrane is considered as a one-dimensional event in the z-direction. The MS description of species i in a mixture of n different particles in an ion exchange membrane in a RED stack, is as follows:
The left side of the equation, the gradient of the chemical potential, represents the driving force (N/mol) of species i. The right side shows the friction between species i and all other species (j = 1, …, n; j ≠ i) and, depending on the difference in average velocity (vi − vj), and the mole fraction of each other species (xj). For each combination of i with j, there is a special Maxwell–Stefan diffusivity Ðij.
The most important driving forces in the RED process are the activity gradient and the electric gradient. Used pressures in RED are low (<1 bar) with respect to the osmotic pressure of seawater (about 25 bar) and can be neglected. If the these driving forces are substituted in Equation (3) the result is:
The first term on the left side indicates the activity contribution where the mole fraction x is corrected by the activity coefficient γ. The second term is the electric part of the driving force. This equation in our model is approximated by:
The superscripts R and S indicate the boundaries in the membrane on the river water side and on the seawater side.
Maxwell–Stefan diffusivities Đ of particles in ion exchange membranes are collected from different sources as shown in
Table 2. The amount of found data is rather poor; moreover, data are related to different membranes, achieved with different methods at different temperatures. The last column in the table lists the median of the different diffusivities; these are used as a starting value in the fit of the MS model.
4. Modelling
Figure 2 shows the design of the MS model. Our model is limited to transport in the membranes; mass transfer in the intermediate feed water compartments was not considered for reasons of simplicity. The variables are the velocities of the particles (v1, v2 and v3) through the CEM and (w1, w2 and w3) through the AEM. Indices 1, 2 and 3 stand for H
2O, Na
+ and Cl
−; U
C and U
A are the voltages over the CEM and AEM. With the solver the 6 velocities and 2 voltages are obtained as functions of the electrical current I. The Maxwell–Stefan model consists of 5 steps which are summarized below. The program code is written in Mathcad Express Prime 9.0.0.0.0 and is available from the author.
4.1. Feed Water Concentrations
Because the experiment was performed with 7 different current densities, 7 inlet concentrations were imported for both inlet streams (S and R) and 7 for both outlet concentrations. The model used the mean value of each stream. For each of the 4 mean concentrations the activity of all ions was calculated using the extended Debye–Hückel equation. The mole fraction of each particle (ions and water) was also calculated.
4.2. Membrane Properties
The relevant membrane properties are the concentration of the fixed charges (
Cn) and the total concentration of all particles (
Ct); they were obtained from the reported values of the swelling degree (
S), density (
ρ) and the ion exchange capacity (
IEC) as shown:
Here, MW is the molar mass of water (0.018 kg/mol).
4.3. Boundary Concentrations within the Membranes
The following compounds are involved in the model: H
2O (1), Na
+ (2), Cl
− (3), for the dissolved particles and ~SO
3− (4) and ~NR
3+ (5) for the fixed charges in the membranes. The number between parentheses is used as an index in the notation of mole fractions, etc. Concentrations in the river and seawater compartments were taken independently of the distance to the membranes. Within the membranes there are four interfaces, indicated by the superscripts CS (CEM–Sea), CR (CEM–River), AS (AEM–Sea) and AR (AEM–River). On each interface the mole fraction of the ions in the membrane phase was calculated by applying the theory of the Donnan equilibrium [
17]. For example, on the sea-CEM interface the equilibria for all charged particles are as follows:
In these expressions
KCS represents the Donnan equilibrium constant for the relevant interface (
CS); the power applied on
K (1 and −1) is the valence of the relevant ion. To solve the system, an extra equation should be added. For reasons of electroneutrality, this is as follows:
Using Equations (8) and (9) the mole fractions x2CS and x3CS and the Donnan equilibrium constant KCS can be solved. Such systems are also solved for the other interfaces, giving the mole fractions xiCS, xiCR, xiAS, xiAR (i = 2, 3) and the equilibrium constants KCS, KCR, KAS and KAR. In this particular case, with only three unknowns, the solver function does not need to be used and the equilibrium constant K and the mole fractions of Na+ and Cl− can be solved with a quadratic equation.
4.4. The Solver Box
The solver box contains 6 equations similar to Equation (5), namely 3 for the transport of H2O, Na+ and Cl− through the CEM and 3 for the transport through the AEM. Additional equations are added for the mass balance through the membranes and for the relation between mass transport and electrical current. With these 8 equations, 8 variables are solved (v1, v2, v3, w1, w2, w3, and the voltages over the two membranes UC and UA).
4.5. Comparison of MS Model and Experiment
Power density was determined from membrane potentials and current density. The internal fluxes of ions and water were derived from the velocities and concentrations of the particles. These data were compared with the experimental values. This is discussed in the next section.
6. Conclusions and Outlook
The Maxwell–Stefan methodology provided credible values for three of the nine diffusivities. The other diffusivities have such high values that they play no role in the model. Nevertheless, after the fitting procedure, the model provided good agreements with the experimental values for (electro)osmosis, salt transport and power density which is promising for further developments in the application of the MS methodology to RED.
To achieve this, it is important to refine both experiment and model. The experiment could be improved by carrying out very accurate flow measurements, not only at the outlets, as in the experiments described in this paper and explained in
Appendix A, but also at the inlet. By working with metering pumps (e.g., piston pumps) instead of peristaltic pumps, internal water transport can be measured directly. In this way, more reliable values can be achieved for the water and salt transport through membranes, parameters that are of decisive importance for the adapted procedure. Another threat is internal fluid leakage from the seawater to the river water compartment or vice versa. This may perhaps be checked with high molecular weight dyes that cannot cross the membranes themselves but can pass through the leaks in solution.
There is still a lot to be done regarding the model. At higher current intensities, internal salt transport increases and therefore the outlet concentrations change. The results described in this paper were produced using average values for inlet and outlet concentrations. As a result, the calculation of the mole fractions at the membrane interfaces is not entirely correct. A remedy could be to perform the adapted procedure for each individual amperage. That would then lead to a set of current dependent diffusivities that could be averaged afterwards.
With better coordination between experiment and model, it would be possible to include divalent ions in experiments and model. This is particularly interesting because divalent ions in a RED setup can move oppositely to monovalent ions with the same charge sign [
14]. Both seawater and river water contain significant amounts of divalent ions, especially Mg
2+, Ca
2+ and SO
42−. A better understanding of all interactions between water and free and fixed ions would be a boost for the development of RED technology as a source of renewable energy. Furthermore, this knowledge would also contribute to the development of ED, especially if this technique is used for the desalination of seawater.