Open AccessArticle

Analysis of Membrane Transport Equations for Reverse Electrodialysis (RED) Using Irreversible Thermodynamics

Wojciech Kujawski

^1,2,*

Andriy Yaroshchuk

Emiliy Zholkovskiy

⁴,

Izabela Koter

¹ and

Stanislaw Koter

^1,*

Faculty of Chemistry, Nicolaus Copernicus University in Toruń, 87-100 Toruń, Poland

MEPhI, National Research Nuclear University, 115409 Moscow, Russia

ICREA & Polytechnic University of Catalonia—Barcelona Tech, 08034 Barcelona, Spain

⁴

Institute of Bio-Colloid Chemistry, National Academy of Sciences of Ukraine, 03680 Kyiv-142, Ukraine

Authors to whom correspondence should be addressed.

Int. J. Mol. Sci. 2020, 21(17), 6325; https://doi.org/10.3390/ijms21176325

Submission received: 31 July 2020 / Revised: 26 August 2020 / Accepted: 27 August 2020 / Published: 31 August 2020

(This article belongs to the Special Issue Ion and Molecule Transport in Membrane Systems 2.0)

Download

Browse Figures

Figure 1
The system: ion-exchange membrane with adjacent boundary diffusion layers. lp—thickness of concentration polarization layer; lm—thickness of the membrane; c′, c″—concentrations of ions in the bulk solution; c′(m), c″(m)—concentrations of ions at the membrane surface. "> Figure 2
Electric power vs. −I calculated using Equation (5) (full ITE) for the real thickness of Nafion 120 equal to lm = 0.3 mm, and the assumed value lm = 0.1 mm, for 0.1–1.0 M and 0.05–3.7 M NaCl, T = 298 K. lp—thickness of the concentration polarization layer, lm—thickness of the membrane. "> Figure 3
Pressure profile across the membrane calculated from Equation (5) (full ITE) for lm = 0.3 mm and lm = 0.1 mm, for 0.1–1.0 M and 0.05–3.7 M NaCl, T = 298 K. Thickness of the concentration polarization layer lp = 0.05 mm. Thickness of the membrane—lm. "> Figure 4
Profiles of electroosmotic (el), osmotic (osm), and pressure-driven (p), components of Jv calculated from Equation (5) (full ITE) for the current density –I = 100 and 200 A/m2; one can notice a dominant role of electroosmotic contribution to Jv at higher current density. "> Figure 5
Dependence of the maximum of power density, Pmax, and the current density, −I, at Pmax on the membrane thickness, lm, for 0.1–1.0 M and 0.05–3.7 M NaCl, T = 298 K, Equation (5), full ITE. "> Figure 6
Dependence of the maximum of power density, Pmax, and the current density, −I, at Pmax on the diffusion boundary layer thickness, lp, for 0.1–1.0 M and 0.05–3.7 M NaCl, T = 298 K, Equation (5), full ITE. "> Figure 7
Power density vs. current density for full model and approximate approaches: TCP—3-coefficients approach (P1, P2, and L12 calculated from <math display="inline"><semantics> <mrow> <msub> <mover accent="true"> <mi>t</mi> <mo>¯</mo> </mover> <mrow> <mn>1</mn> <mo>,</mo> <mi>a</mi> <mi>p</mi> <mi>p</mi> </mrow> </msub> </mrow> </semantics></math>, km, and Ps), (t1app, Ps); TP, CP, TC—2-coefficients approach (P1 and P2 calculated from <math display="inline"><semantics> <mrow> <msub> <mover accent="true"> <mi>t</mi> <mo>¯</mo> </mover> <mrow> <mn>1</mn> <mo>,</mo> <mi>a</mi> <mi>p</mi> <mi>p</mi> </mrow> </msub> </mrow> </semantics></math> and Ps, km and Ps, <math display="inline"><semantics> <mrow> <msub> <mover accent="true"> <mi>t</mi> <mo>¯</mo> </mover> <mrow> <mn>1</mn> <mo>,</mo> <mi>a</mi> <mi>p</mi> <mi>p</mi> </mrow> </msub> </mrow> </semantics></math> and km, respectively); the Nafion 120 membrane separates (a) 0.1 and 1 M NaCl, (b) 0.05 and 3.7 M NaCl; lm = 0.1 mm, lp = 0.05 mm, T = 298 K. "> Figure 7 Cont.
Power density vs. current density for full model and approximate approaches: TCP—3-coefficients approach (P1, P2, and L12 calculated from <math display="inline"><semantics> <mrow> <msub> <mover accent="true"> <mi>t</mi> <mo>¯</mo> </mover> <mrow> <mn>1</mn> <mo>,</mo> <mi>a</mi> <mi>p</mi> <mi>p</mi> </mrow> </msub> </mrow> </semantics></math>, km, and Ps), (t1app, Ps); TP, CP, TC—2-coefficients approach (P1 and P2 calculated from <math display="inline"><semantics> <mrow> <msub> <mover accent="true"> <mi>t</mi> <mo>¯</mo> </mover> <mrow> <mn>1</mn> <mo>,</mo> <mi>a</mi> <mi>p</mi> <mi>p</mi> </mrow> </msub> </mrow> </semantics></math> and Ps, km and Ps, <math display="inline"><semantics> <mrow> <msub> <mover accent="true"> <mi>t</mi> <mo>¯</mo> </mover> <mrow> <mn>1</mn> <mo>,</mo> <mi>a</mi> <mi>p</mi> <mi>p</mi> </mrow> </msub> </mrow> </semantics></math> and km, respectively); the Nafion 120 membrane separates (a) 0.1 and 1 M NaCl, (b) 0.05 and 3.7 M NaCl; lm = 0.1 mm, lp = 0.05 mm, T = 298 K. "> Figure 8
Power density vs. current density for full model—comparison of the TC const approach with logarithmic mean with the TC approach and the full ITE for the concentrations 0.1–1 M and 0.05–3.7 M NaCl; lm = 0.1 mm, lp = 0.05 mm, T = 298 K. ">

Versions Notes

Abstract

Reverse electrodialysis (RED) is an electro-membrane process for the conversion of mixing energy into electricity. One important problem researchers’ face when modeling the RED process is the choice of the proper membrane transport equations. In this study, using experimental data that describe the membrane Nafion 120 in contact with NaCl aqueous solutions, the linear transport equation of irreversible thermodynamics was applied to calculate the power density of the RED system. Various simplifying assumptions about transport equation (i.e., four-, three-, and two-coefficients approaches) are proposed and discussed. We found that the two-coefficients approach, using the membrane conductivity and the apparent transport number of ions, describes the power density with good accuracy. In addition, the influence of the membrane thickness and the concentration polarization on the power density is also demonstrated.

Keywords:

reverse electrodialysis; irreversible thermodynamics; power density; Nafion 120 membrane

1. Introduction

Global population growth contributes to the escalation of energy demand for households, transportation, and industry. The use of fossil fuels results in various environmental threats. Therefore, the search is increasingly widespread for alternative methods of energy conversion, e.g., osmotic into electrical [1] and mechanical energy (so called “osmotic pump” [2,3], or the conversion of electrokinetic energy [4]. Reverse electrodialysis (RED) is an electro-membrane process for the conversion of mixing energy into electricity. The concept of RED was invented in 1950 [5,6]; however, it was not heavily investigated until the 21st century, when this technique regained the attention of researchers looking for novel, unconventional, and renewable energy sources [7,8,9].

One of the important RED issues is the cost-effectiveness of the process. Turek and Bandura [10] suggested that energy conversion using RED is not profitable; however, Giacalone et al. [11] presented a much more optimistic perspective. There are a number of parameters that influence the efficiency of RED, like ionic shortcut currents [12], type of electrodes [13], kind of membrane [14], membrane thickness [15], fouling issues [16,17], efficiency of solutions mixing in the stack [18,19], and geometry of spacers [20].

In this section, we present a description of membrane transport in RED. In the literature discussing the modeling of RED, apart from very simplified membrane transport (i.e., considering only the transport numbers of ions [21]), usually the Nernst–Planck equation or its transformations are used [22,23,24,25,26,27,28]. In some research, the volume flow caused by the osmotic and electroosmotic transport of water through the membrane was considered [25,29]. Whereas in other research, only osmotic transport was considered [22,24]. In other studies, no volume transport was anticipated [27]. The problem of membrane transport is probably neglected because of the many other important aspects of RED. A general (i.e., not exactly related to RED) review of transport through ion-exchange membranes can be found elsewhere [30].

The Nernst–Planck equation is a simplified version of the linear transport equation of irreversible thermodynamics. The transport equation is comprised of six transport coefficients for the system membrane in contact with a single electrolyte aqueous solution, under the condition that the Onsager symmetry of the coefficients is fulfilled [31]. The determination of all these coefficients requires extensive work especially because their concentration dependence should be considered [32,33,34,35]. Consequently, the practical use of the full transport equation is limited.

Therefore, the aim of this work was to discuss which of the simplifications of the irreversible thermodynamics transport equation (ITE) is useful for calculating the power density in RED. The following three approaches based on four, three, and two coefficients (to be determined experimentally) were analyzed and discussed in detail:

(1): the four-coefficients approach—the pressure difference in the full set of ITE is neglected,
(2): the three-coefficients approach—the coupling of ion flows with the gradient of chemical potential of water is neglected,
(3): the two-coefficients approach (i.e., Nernst–Planck equation)—the ion flux is a function of the gradient of its electrochemical potential only.

The results of these simplifications (i.e., values of power density vs. electric current density) are compared with the values obtained from the full set of ITE. The transport coefficients used in these latter equations were determined for the cation-exchange membrane Nafion 120 in contact with NaCl solutions [34,35] in a wide range of concentrations (0.05–3.7 M). To consider the concentration dependence, the differential form of the transport equation was used. This method was applied for the first time by McCallum and Meares to describe transport through an ion-exchange membrane [36]. The Nafion 120 membrane is relatively thick (0.3 mm); therefore, the transport coefficients are rescaled to smaller thickness values to demonstrate the influence of that parameter on the power density.

2. Theory

In the following, one basic part of the RED module is considered, i.e., one ion-exchange membrane with adjacent concentration polarization layers (Figure 1). The steady-state conditions in the system are assumed, i.e., the fluxes of ions in the membrane and both polarization layers are the same.

2.1. Transport Equations of Irreversible Thermodynamics (ITE)

The starting point for the ITE approach to isothermal conditions is the dissipation function, Φ, the product of entropy production, and temperature, which can be expressed as the sum of driving forces, X_i, and coupled with fluxes, J_i:

Φ = T σ = \sum_{i} J_{i} X_{i} > 0 .

For discontinuous systems, the forces are represented by the differences in chemical potentials of transported species. In continuous systems, the negative values of gradients of these potentials are used.

According to irreversible thermodynamics, it is assumed that fluxes and forces are linearly related:

J_{i} = L_{i 1} X_{1} + L_{i 2} X_{2} + \dots

(1)

and the Onsager reciprocity relation L_ik = L_ki holds if the forces are not too high.

One can choose many sets of fluxes and forces; all of them have to fulfill the condition of the entropy production invariance:

Φ = \sum_{i} J_{i}^{'} X_{i}^{'} = \sum_{i} J_{i}^{″} X_{i}^{″} .

(2)

From this condition, we obtain various forms of the fluxes and forces, among which is the set of practical equations useful for experimental determination of transport coefficients [37]:

\begin{matrix} J_{1} / ν_{1} = & L_{π π} Δ π / {\tilde{c}}_{s} + L_{π p} (Δ p - Δ π) + L_{π E} Δ E \\ J_{v} = & L_{p π} Δ π / {\tilde{c}}_{s} + L_{p p} (Δ p - Δ π) + L_{p E} Δ E \\ I = & L_{E π} Δ π / {\tilde{c}}_{s} + L_{E p} (Δ p - Δ π) + L_{E E} Δ E \end{matrix}

(3)

where J₁ is the flux of ions 1 (cations), ν₁ is the number of ions 1 in the electrolyte molecule, s is volume flux, I is electric current density, Δp is the pressure difference, Δπ is the osmotic pressure difference,

{\tilde{c}}_{s}

is the mean concentration of electrolyte solutions contacting the membrane (c_s′, c_s″), and

{\tilde{L}}_{α β}

is the mean transport coefficient for that concentration range. ΔE is the electric potential difference measured with electrodes reversible to anions (2):

Δ E = Δ {\tilde{μ}}_{2} / z_{2} F = Δ φ + \frac{R T}{z_{2} F} Δ \ln a_{2} .

(4)

The power density is calculated from the electric potential difference Δφ, determined from Equation (4) assuming that

Δ \ln a_{2} \approx Δ \ln a_{\pm}

The determination of all L_αβ coefficients is a laborious task. In Table 1, we list the experiments we performed to determine L_αβ in Equation (3).

Having these coefficients determined for the vanishingly small concentration differences (Δπ, Δμ_s → 0) allowed us to assume the Onsager reciprocity relation (L_αβ = L_βα) and determine the concentration dependence of six L_αβ (α, β = π, p, E) [34,35]. The method of such determination was described in detail elsewhere [32]. The finite concentration differences (L_αβ) in Equation (3) are not symmetric because of their strong concentration dependence [32,38].

Calculating the ion flux, volume flux, and electric potential difference, considering the concentration dependence of L_αβ is the most exact approach. To calculate these quantities for the assumed conditions (concentrations on both sides of membrane, electric current density, membrane thickness), we applied the method proposed by McCallum and Meares [36]. According to this method, the membrane is divided into n_m slices and for each slice the transport is given by the set of Equation (3) with local variables (electrolyte 1:1):

\begin{matrix} J_{1} = & (L_{π π} ({\tilde{c}}_{s}) Δ π / {\tilde{c}}_{s} + L_{π p} ({\tilde{c}}_{s}) (Δ p - Δ π) + L_{π E} ({\tilde{c}}_{s}) Δ E) l_{m} / w_{s l} \\ J_{v} = & (L_{p π} ({\tilde{c}}_{s}) Δ π / {\tilde{c}}_{s} + L_{p p} ({\tilde{c}}_{s}) (Δ p - Δ π) + L_{p E} ({\tilde{c}}_{s}) Δ E) l_{m} / w_{s l} \\ I = & (L_{E π} ({\tilde{c}}_{s}) Δ π / {\tilde{c}}_{s} + L_{E p} ({\tilde{c}}_{s}) (Δ p - Δ π) + L_{E E} ({\tilde{c}}_{s}) Δ E) l_{m} / w_{s l} \end{matrix}

(5)

where:

Δ Y = Y_{k - 1} - Y_{k}, Y = π, p, E

{\tilde{c}}_{s} = (c_{s, k - 1} + c_{s, k}) / 2

where l_m is the thickness of the real membrane for which the coefficients L_αβ = f(c_s) were determined and w_sl is the width of slice equal to the assumed membrane thickness (l_m_,assumed) divided by the number of slices (w_sl = l_m_,assumed/n_m). The concentration difference is expressed as Δπ and the formula

Δ π = \partial π / \partial c_{s} (\tilde{c}) (c_{s, k - 1} - c_{s, k})

The fluxes across each slice are the same, assuming the stationary state in the considered system. Therefore, for n_m slices, we have 3n equations and 3n_m unknowns: c_s,k, p_k, E_k, k = 1, …, n − 1, J₁, and J_v, which we solve for given values of parameters (electric current, concentration, pressure on both sides of membrane, and the assumed membrane thickness). All of these refer to the hypothetical (virtual) solution being in equilibrium with the membrane at coordinate x = kw_sl, k = 1, …, n − 1; such a solution is defined by equality of chemical potentials of the transported species in that solution and in the membrane at a given x.

The concentration polarization layers adjacent to membrane surfaces were treated in a similar way as the membrane. The experimental data for NaCl solutions were taken from [39]. The used equations are provided in Appendix A.

A comment should be made on the extended Nernst–Planck equation, frequently used in the literature for describing the membrane transport of ions. Written in terms of hypothetical solution parameters, it takes the following form:

J_{i} = \frac{P_{i} c_{i}}{R T} (- \frac{d {\tilde{μ}}_{i}}{d x}) + c_{i} K_{p, i} α_{i} J_{v}

where K_p,i is the partition coefficient of ion I and α_i is the convective coupling coefficient [40,41]. There are 2 coefficients per ion (K_p,iα_i is treated as one coefficient); thus, for two ions and J_v described by 3 coefficients (L_pπ, L_pp, L_pE), there are 7 coefficients—more than in the set of Equation (3) assuming the symmetry of coefficients. For this reason, the extended Nernst–Planck equation will not be discussed here.

2.2. The Four-Coefficients Approach

In this case, the pressure difference across each slice is assumed to be zero. Although on both sides of membrane the pressure is the same, it may change across the membrane to maintain the volume flux across each slice the same. As the number of unknowns is reduced, the equation for J_v in Equation (5) can be neglected (J₁ is needed in the concentration polarization layers). Thus Equation (5) is reduced to:

\begin{matrix} J_{1} = & ({L^{'}}_{π π} ({\tilde{c}}_{s}) Δ π + L_{π E} ({\tilde{c}}_{s}) Δ E) l_{m} / w_{s l} \\ I = & ({L^{'}}_{E π} ({\tilde{c}}_{s}) Δ π + L_{E E} ({\tilde{c}}_{s}) Δ E) l_{m} / w_{s l} \end{matrix}

(6)

where

{L^{'}}_{π π} ({\tilde{c}}_{s}) = L_{π π} ({\tilde{c}}_{s}) / {\tilde{c}}_{s} - L_{π p} ({\tilde{c}}_{s})

{L^{'}}_{E π} ({\tilde{c}}_{s}) = L_{E π} ({\tilde{c}}_{s}) / {\tilde{c}}_{s} - L_{E p} ({\tilde{c}}_{s})

The four coefficients (L′_ππ, L_πE, L′_Eπ, L_EE) can be determined from P_s,dif,

{\bar{t}}_{1}

, k_m, and

{\bar{t}}_{1, a p p}

—the Electromotive Force (EMF) measurements. However, practically, it is better to determine the electroosmotic coefficient W (Table 1) and to calculate

{\bar{t}}_{1}

from W and

{\bar{t}}_{1, a p p}

2.3. The Three-Coefficients Approach

In the equation of ion transport, the coupling with water, i.e., the term L_i₀X₀ is neglected. Thus, Equation (1) with L_ii replaced by P_ic_i/RT reduces to:

\begin{array}{l} J_{1} = \frac{P_{1} c_{1}}{R T} X_{1} + L_{12} X_{2} \\ J_{2} = L_{21} X_{1} + \frac{P_{2} c_{2}}{R T} X_{2} \end{array}

(7)

where X_i is:

X_{i} \equiv - \frac{d {\tilde{μ}}_{i}}{d x} = - (R T \frac{d \ln a_{i}}{d x} + z_{i} F \frac{d ϕ}{d x}) .

(8)

We assume that L₂₁ = L₁₂; therefore, there are 3 parameters that can be obtained from (1) specific conductivity, k_m, (2) electrolyte permeability P_s, and (3) the apparent transport number,

{\bar{t}}_{1, a p p}

. Instead of

{\bar{t}}_{1, a p p}

, we could choose the real transport number of cation,

{\bar{t}}_{1}

; however, because

{\bar{t}}_{1, a p p}

includes the electroosmotic number of water and is much easier to determine than

{\bar{t}}_{1}

, we choose

{\bar{t}}_{1, a p p}

. Notably,

{\bar{t}}_{1, a p p}

and

{\bar{t}}_{1}

obtained from Equation (7) are the same because the flux of water is neglected here (see formulae for

{\bar{t}}_{1}

and

{\bar{t}}_{1, a p p}

in Table 1). The EMF measurements allow for the determination of the apparent transport number of ion 1. Using Equation (7), we obtain the following set of practical equations (for the 1:1 electrolyte):

\begin{matrix} J_{1} = & (P_{1} ({\tilde{c}}_{s}) (1 - {\tilde{c}}_{s} {\bar{v}}_{s}) Δ π / R T + ({\bar{t}}_{1, a p p} κ_{m} / F) Δ E) / w_{s l} \\ I = & (({\bar{t}}_{1, a p p} κ_{m} / F) (1 - {\tilde{c}}_{s} {\bar{v}}_{s}) Δ π / {\tilde{c}}_{s} + κ_{m} Δ E) / w_{s l} \end{matrix}

(9)

where P₁ is given by:

P_{1} = P_{s} + R T κ_{m} {({\bar{t}}_{1, a p p} / F)}^{2} / c_{s},

and P_s by:

P_{s} \equiv \frac{R T}{(1 - c_{s} {\bar{v}}_{s})} {(\frac{J_{1}}{- d π / d x})}_{I = 0} .

2.4. The Two-Coefficients Approach

In this approach, we assume that the flux of ion I depends only on the gradient of electrochemical potential of that ion:

J_{i} = L_{i i} X_{i} = - \frac{P_{i} c_{i}}{R T} \frac{d {\tilde{μ}}_{i}}{d x} i = 1, 2 .

(10)

The following practical equations are derived from Equation (10):

\begin{matrix} J_{1} = & P_{1} ({\tilde{c}}_{s}) ((1 - {\tilde{c}}_{s} {\bar{v}}_{s}) Δ π + {\tilde{c}}_{s} F Δ E) / (R T w_{s l}) \\ I = & P_{1} ({\tilde{c}}_{s}) (F (1 - {\tilde{c}}_{s} {\bar{v}}_{s}) Δ π + F^{2} {\tilde{c}}_{s} (1 + P_{2} ({\tilde{c}}_{s}) / P_{1} ({\tilde{c}}_{s})) Δ E) / (R T w_{s l}) \end{matrix}

(11)

There are 2 transport coefficients only, P₁ and P₂, with four experimentally available parameters: (1) specific conductivity k_m, (2) transport number of cation

{\bar{t}}_{1}

, (3) electrolyte permeability P_s, and (4)

{\bar{t}}_{1, a p p}

. The relation between these quantities and P₁, P₂ is as follows (assuming 1:1 electrolyte, and denoting ion 1 as a cation):

κ_{m} \equiv {(\frac{I}{- d E / d x})}_{c} = \frac{F^{2}}{R T} c_{s} (P_{1} + P_{2}),

(12)

{\bar{t}}_{1} \equiv {(\frac{F J_{1}}{I})}_{c} = \frac{P_{1}}{P_{1} + P_{2}},

(13)

P_{s} \equiv \frac{R T}{(1 - c_{s} {\bar{v}}_{s})} {(\frac{J_{1}}{- d π / d x})}_{I = 0} = \frac{P_{1} P_{2}}{P_{1} + P_{2}},

(14)

and

{\bar{t}}_{1, a p p} = {(\frac{d {\tilde{μ}}_{2}}{d μ_{s}})}_{I = 0} = - F {(\frac{d E}{d μ_{s}})}_{I = 0} = \frac{P_{1}}{P_{1} + P_{2}} .

(15)

Similar to the 3-coefficients approach, the expressions for

{\bar{t}}_{1}

and

{\bar{t}}_{1, a p p}

are the same because the water transport is also neglected.

The determination of

{\bar{t}}_{1, a p p}

by the EMF method is much easier and therefore much more popular than the determination of the real transport number,

{\bar{t}}_{1}

. We chose the following pairs of experiments to determine P₁ and P₂:

(1) k_m,

{\bar{t}}_{1, a p p}

, (2) k_m, P_s, (3)

{\bar{t}}_{1, a p p}

, P_s.

Additionally, contrary to

{\bar{t}}_{1}

, the value of

{\bar{t}}_{1, a p p}

includes the electroosmotic transport of water; thus, to some extent, the transport of water molecules is considered.

3. Results and Discussion

As mentioned above, in this work, only one membrane (cation-exchange membrane) with adjacent concentration polarization layers is discussed.

In Figure 2, a typical dependence of the power density P on the current density I, calculated as P = −IΔφ, based on the full transport equation (Equation (5)), is shown for two concentration differences (0.1–1.0 M and 0.05–3.7 M) and two values of membrane thickness (0.1 and 0.3 mm). Two common trends were observed: the higher the concentrations difference, the higher the value of P; the lower the membrane thickness, the higher the value of P.

The four-coefficients approach, neglecting the pressure changes (Equation (6)), yields similar results. For 0.1–1.0 M NaCl, the difference in P_max is below 1% irrespective of the membrane thickness. For higher concentration differences (0.05–3.7 M), this difference increases to ca. 5–6%.

Changes of pressure in the hypothetical/virtual solution affect the ion fluxes. Some examples of pressure variations are shown in Figure 3. The negative values of pressure indicated that there must have been an underpressure inside the membrane needed to keep the volume flow constant. Similar results were obtained by McCallum and Meares [36]. They suggested that the membrane thickness could decrease in that case. At higher I values, the pressure becomes positive. This can be explained by different changes in the profiles of electroosmotic and osmotic parts of J_v across the membrane caused by the change of concentration profiles (Figure 4). To keep J_v constant, the pressure part of J_v should change accordingly, which demands various profiles of pressure. Assuming the constant pressure value (the four-coefficients approach), the volume flux across the membrane is not constant—it can deviate from the actual steady state value even several times, which is unrealistic.

The influence of the membrane thickness on the maximum of P and the current density at P_max is shown in Figure 5. On both curves (P_max = f(l_m) and −I(P_max) = f(l_m)), maxima are observed at l_m < 0.1 mm, and the maximum on P_max is at a slightly higher l_m value than that of −I(P_max). The higher concentration difference demands a thicker membrane (maxima are shifted toward thicker membranes). The thinner the membrane, the higher I, but also the higher the permeability of electrolytes; the impact of diffusion layers (of the assumed constant thickness) also becomes relatively stronger.

The dependence of the maximum of power density P_max and the current density at P_max on the diffusion boundary layer thickness l_p is a monotonically decreasing function (Figure 6). Due to a substantial decrease in P_max with l_p, the reduction of the thickness of that layer by ensuring adequate stirring of the solutions at the membrane surface is an important practical issue. Even more important is the decrease in l_p on the lower concentration side of membrane. To reduce the concentration polarization effect, an extra energy for pumping is needed.

The Two-Coefficients Approach

The results for the two-coefficients approach, for which P₁ and P₂ are calculated either from

{\bar{t}}_{1, a p p}

and k_m (TC), or from

{\bar{t}}_{1, a p p}

and P_s (TP), or from k_m and P_s (CP), are shown in Figure 7. Among these three combinations of coefficients, the smallest deviation from the power density calculated from the full ITE was observed for the TC combination (

{\bar{t}}_{1, a p p}

, k_m). It yielded even better results than the three-coefficients approach based on all three parameters (TCP:

{\bar{t}}_{1, a p p}

, k_m and P_s).

Because of the quite good fit of the TC approach, another question arose: is it necessary to determine the concentration dependence of

{\bar{t}}_{1, a p p}

and k_m in the whole range c′–c″, or is it enough to determine

{\tilde{\bar{t}}}_{1, a p p}

in that range and k_m for the mean concentration of c′ and c″?

In that case (P₁, P₂ = const), we can determine the integrated form for J₁ from Equation (10):

J_{1} = \frac{P_{1}}{P_{1} + P_{2}} j / F - \frac{P_{1} P_{2}}{(P_{1} + P_{2}) R T l_{m}} \int_{{c^{'}}_{s}}^{{c^{″}}_{s}} (1 - c_{s} {\bar{v}}_{s}) \frac{\partial π}{\partial c_{s}} d c_{s} .

(16)

ΔE can be expressed as a function of a certain mean concentration of c′ and c″:

Δ E \equiv E^{'} - E^{″} = \frac{R T}{F} \frac{(j / F - J_{1}) l_{m}}{{\tilde{c}}_{s} P_{2}} .

(17)

P₁ and P₂ are determined from k_m (Equation (12)), with c_s replaced by

{\tilde{c}}_{s}

, and

{\tilde{\bar{t}}}_{1, a p p}

calculated from (15) with dE/dμ_s replaced by (E(c″) − E(c′))/(μ_s(c″) − μ_s(c′)). To check this simplification, the calculations were performed for

{\tilde{c}}_{s}

being the arithmetic and logarithmic means of c′, c″. The boundary diffusion layers were treated as in the previous cases to keep the same conditions for all discussed approaches. To simplify calculations, in the calculations of P₁, P₂, and

{\tilde{c}}_{s}

, the external bulk concentrations were used, not those at the membrane surfaces. A comparison of this approach (TC const.) with the full ITE equation and with the TC approach is shown in Figure 8.

This simplification of the logarithmic mean of concentration (marked in Figure 8 as “TC const ln”) yields comparable results to the other approaches for the lower concentration difference (0.1–1.0 M). For 0.05–3.7 M the predictions of “TC const ln” are too low at the current density range 0–j(P_max).

In this case, the concentration dependence of

{\bar{t}}_{1, a p p}

and k_m is needed to improve the predictions of power density. The arithmetic mean of concentration yields much worse results especially for high concentration difference and is not shown in Figure 8.

4. Conclusions

Among the tested simplifications of the full set of transport equations of irreversible thermodynamics (six independent coefficients), the four-coefficients approach (k_m,

{\bar{t}}_{1}

{\bar{t}}_{1, a p p}

, P_s_,dif) with only the pressure effect neglected yields almost the same results for 0.1–1 M as the full ITE approach.

The three-coefficients approach (k_m,

{\bar{t}}_{1, a p p}

, P_s_,dif), where only ion fluxes are considered with their couplings, slightly overestimates the power density at lower concentration difference. However, at higher concentration differences, it strongly overestimates and becomes much less accurate than the two-coefficients approach (equation of the Nernst–Planck type, assuming no coupling between ion fluxes) based on k_m and

{\bar{t}}_{1, a p p}

. Other combinations of the two-coefficients approach are substantially worse. The model based on k_m and

{\bar{t}}_{1, a p p}

can be further simplified by using constants P₁ and P₂ calculated from the mean values of k_m and

{\bar{t}}_{1, a p p}

in a given concentration range c′–c″ and using the logarithmic mean of these concentrations. These coefficients are easily accessible experimentally.

None of the discussed approximations considered the volume flow effect on the concentrations of the solutions contacting the membrane. To do so, we need to estimate the osmotic volume permeability and the electroosmotic coefficient W.

To obtain the optimal power density for given concentrations of feed solutions, the appropriate membrane thickness should be chosen. The higher the concentration difference, the thicker the membrane should be. As the position of P_max depends on the membrane properties, the optimal thickness should be independently determined for cation- and anion-exchange membranes.

The thickness of the diffusion boundary layer should be as low as possible especially on the low concentration side of the membrane—it strongly decreases the power density.

Author Contributions

Conceptualization, S.K.; Formal analysis, A.Y., E.Z., and S.K.; Funding acquisition, W.K.; Methodology, S.K.; Project administration, W.K.; Software, S.K.; Visualization, I.K. and S.K.; Writing—original draft, S.K.; Writing—review & editing, W.K., A.Y., and I.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by funds for fundamental research from Nicolaus Copernicus University in Toruń, Poland (Faculty of Chemistry—“Membranes and membrane separation processes—fundamental and applied research”). The APC was funded by International Journal of Molecular Science.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Transport Equation in the Concentration Polarization Layers

In [39], the transport coefficients, l_ik, are determined for the set of equations:

J_{i}^{(0)} = l_{i 1} (- \frac{d {\tilde{μ}}_{1}}{d x}) + l_{i 2} (- \frac{d {\tilde{μ}}_{2}}{d x}) i = 1, 2

(A1)

in the solvent-fixed reference frame. Equation (A1) can be transformed with forces similar to those in Equation (3) (electrolyte 1:1):

\begin{array}{l} J_{1}^{(0)} = l_{11} a (- \frac{d π}{d x}) + l_{1 E} (- \frac{d E}{d x}) \\ I = l_{1 E} a (- \frac{d π}{d x}) + κ_{e} (- \frac{d E}{d x}) \end{array}

(A2)

where

l_{1 E} = F (l_{11} - l_{12})

κ_{e} = F^{2} (l_{11} - 2 l_{12} + l_{22})

, and

a = (1 - c_{s} {\bar{v}}_{s}) / c_{s}

. Due to volume flow through the membrane, we assume that the flux of ion 1 in the boundary layer is given by:

J_{1} = J_{1}^{(0)} + \frac{c_{s}}{c_{0}} J_{0} .

(A3)

In the practical equations (Equations (3) and (5)), J_v is defined as

J_{v} = {\bar{v}}_{0} J_{0} + {\bar{v}}_{s} J_{1} / ν_{1}

; thus, for our case (electrolyte 1:1), Equation (A3) can be rewritten as:

J_{1} = J_{1}^{(0)} (1 - c_{s} {\bar{v}}_{s}) + c_{s} J_{v} .

To simultaneously calculate the c and E profiles in the boundary layers with the membrane, we divide both layers into n_e slices for which the transport equations are as follows:

\begin{matrix} J_{1} = & (l_{11} ({\tilde{c}}_{s}) a ({\tilde{c}}_{s}) Δ π + L_{1 E} ({\tilde{c}}_{s}) Δ E) (1 - {\tilde{c}}_{s} {\bar{v}}_{s}) / w_{s l, e} + {\tilde{c}}_{s} J_{v} \\ I = & (L_{1 E} ({\tilde{c}}_{s}) a ({\tilde{c}}_{s}) Δ π + κ_{e} ({\tilde{c}}_{s}) Δ E) / w_{s l, e} \end{matrix}

(A4)

where w_sl,e is the width of the slices in those layers. Dividing each layer into n_p slices we have (4n_p + 3n_m) equations (including equations for the membrane) and the same number of unknowns: (2n_p + n_m − 1) of c_s,k, (2n_p + n_m) of E_k, (n_m − 1) of p_k, and J₁, J_v.

References

Lee, K.H.; Baker, R.W.; Lonsdale, H.K. Membranes for power generation by pressure retarded osmosis. J. Membr. Sci. 1981, 8, 141–171. [Google Scholar] [CrossRef]
Levenspiel, O.; De Nevers, N. The osmotic pump. Science 1974, 183, 157–160. [Google Scholar] [CrossRef] [PubMed]
Narebska, A.; Koter, S.; Kujawski, W. Conversion of Osmotic into Mechanical Energy in Systems with Charged Membranes. J. Non Equilib. Thermodyn. 1990, 15, 1–10. [Google Scholar] [CrossRef]
Kilsgaard, B.S.; Haldrup, S.; Catalano, J.; Bentien, A. High figure of merit for electrokinetic energy conversion in Nafion membranes. J. Power Source 2014, 247, 235–242. [Google Scholar] [CrossRef]
Manecke, G. Membranakkumulator. Z. Physik. Chem. 1952, 201, 1–15. [Google Scholar] [CrossRef]
Pattle, R.F. Production of Electric Power by mixing Fresh and Salt Water in the Hydroelectric Pile. Nature 1954, 174, 660. [Google Scholar] [CrossRef]
Dlugolecki, P.; Gambier, A.; Nijmeijer, K.; Wessling, M. Practical Potential of Reverse Electrodialysis As Process for Sustainable Energy Generation. Environ. Sci. Technol. 2009, 43, 6888–6894. [Google Scholar] [CrossRef]
Veerman, J.; de Jong, R.M.; Saakes, M.; Metz, S.J.; Harmsen, G.J. Reverse electrodialysis: Comparison of six commercial membrane pairs on the thermodynamic efficiency and power density. J. Membr. Sci. 2009, 343, 7–15. [Google Scholar] [CrossRef]
Veerman, J.; Saakes, M.; Metz, S.J.; Harmsen, G.J. Electrical Power from Sea and River Water by Reverse Electrodialysis: A First Step from the Laboratory to a Real Power Plant. Environ. Sci. Technol. 2010, 44, 9207–9212. [Google Scholar] [CrossRef]
Turek, M.; Bandura, B. Renewable energy by reverse electrodialysis. Desalination 2007, 205, 67–74. [Google Scholar] [CrossRef]
Giacalone, F.; Papapetrou, M.; Kosmadakis, G.; Tamburini, A.; Micale, G.; Cipollina, A. Application of reverse electrodialysis to site-specific types of saline solutions: A techno-economic assessment. Energy 2019, 181, 532–547. [Google Scholar] [CrossRef]
Veerman, J.; Post, J.W.; Saakes, M.; Metz, S.J.; Harmsen, G.J. Reducing power losses caused by ionic shortcut currents in reverse electrodialysis stacks by a validated model. J. Membr. Sci. 2008, 310, 418–430. [Google Scholar] [CrossRef] [Green Version]
Veerman, J.; Saakes, M.; Metz, S.J.; Harmsen, G.J. Reverse electrodialysis: Evaluation of suitable electrode systems. J. Appl. Electrochem. 2010, 40, 1461–1474. [Google Scholar] [CrossRef] [Green Version]
Guler, E.; Elizen, R.; Vermaas, D.A.; Saakes, M.; Nijmeijer, K. Performance-determining membrane properties in reverse electrodialysis. J. Membr. Sci. 2013, 446, 266–276. [Google Scholar] [CrossRef]
Tedesco, M.; Hamelers, H.V.M.; Biesheuvel, P.M. Nernst-Planck transport theory for (reverse) electrodialysis: III. Optimal membrane thickness for enhanced process performance. J. Membr. Sci. 2018, 565, 480–487. [Google Scholar] [CrossRef] [Green Version]
Vermaas, D.A.; Kunteng, D.; Saakes, M.; Nijmeijer, K. Fouling in reverse electrodialysis under natural conditions. Water Res. 2013, 47, 1289–1298. [Google Scholar] [CrossRef]
Pintossi, D.; Saakes, M.; Borneman, Z.; Nijmeijer, K. Electrochemical impedance spectroscopy of a reverse electrodialysis stack: A new approach to monitoring fouling and cleaning. J. Power Source 2019, 444, 227302. [Google Scholar] [CrossRef]
Vermaas, D.A.; Saakes, M.; Nijmeijer, K. Enhanced mixing in the diffusive boundary layer for energy generation in reverse electrodialysis. J. Membr. Sci. 2014, 453, 312–319. [Google Scholar] [CrossRef]
Pawlowski, S.; Rijnaarts, T.; Saakes, M.; Nijmeijer, K.; Crespo, J.G.; Velizarov, S. Improved fluid mixing and power density in reverse electrodialysis stacks with chevron-profiled membranes. J. Membr. Sci. 2017, 531, 111–121. [Google Scholar] [CrossRef]
Mehdizadeh, S.; Yasukawa, M.; Abo, T.; Kakihana, Y.; Higa, M. Effect of spacer geometry on membrane and solution compartment resistances in reverse electrodialysis. J. Membr. Sci. 2019, 572, 271–280. [Google Scholar] [CrossRef]
Jeong, H.I.; Kim, H.J.; Kim, D.K. Numerical analysis of transport phenomena in reverse electrodialysis for system design and optimization. Energy 2014, 68, 229–237. [Google Scholar] [CrossRef]
Veerman, J.; Saakes, M.; Metz, S.J.; Harmsen, G.J. Reverse electrodialysis: A validated process model for design and optimization. Chem. Eng. J. 2011, 166, 256–268. [Google Scholar] [CrossRef]
Kim, D.K. Numerical study of power generation by reverse electrodialysis in ion-selective nanochannels. J. Mech. Sci. Technol. 2011, 25, 5–10. [Google Scholar] [CrossRef]
Tedesco, M.; Cipollina, A.; Tamburini, A.; van Baak, W.; Micale, G. Modelling the Reverse ElectroDialysis process with seawater and concentrated brines. Desalin. Water Treat. 2012, 49, 404–424. [Google Scholar] [CrossRef] [Green Version]
Tedesco, M.; Cipollina, A.; Tamburini, A.; Bogle, I.D.L.; Micale, G. A simulation tool for analysis and design of reverse electrodialysis using concentrated brines. Chem. Eng. Res. Des. 2015, 93, 441–456. [Google Scholar] [CrossRef] [Green Version]
Tedesco, M.; Hamelers, H.V.M.; Biesheuvel, P.M. Nernst-Planck transport theory for (reverse) electrodialysis: I. Effect of co-ion transport through the membranes. J. Membr. Sci. 2016, 510, 370–381. [Google Scholar] [CrossRef]
Kim, H.; Jeong, N.; Yang, S.; Choi, J.; Lee, M.S.; Nam, J.Y.; Jwa, E.; Kim, B.; Ryu, K.S.; Choi, Y.W. Nernst-Planck analysis of reverse-electrodialysis with the thin-composite pore-filling membranes and its upscaling potential. Water Res. 2019, 165, 114970. [Google Scholar] [CrossRef]
Moya, A.A. Uphill transport in improved reverse electrodialysis by removal of divalent cations in the dilute solution: A Nernst-Planck based study. J. Membr. Sci. 2020, 598, 117784. [Google Scholar] [CrossRef]
Tedesco, M.; Hamelers, H.V.M.; Biesheuvel, P.M. Nernst-Planck transport theory for (reverse) electrodialysis: II. Effect of water transport through ion-exchange membranes. J. Membr. Sci. 2017, 531, 172–182. [Google Scholar] [CrossRef] [Green Version]
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] [Green Version]
Kedem, O.; Katchalsky, A. Permeability of composite membranes. Part 1. Electric current, volume flow and flow of solute through membranes. Trans. Faraday Soc. 1963, 59, 1918–1930. [Google Scholar] [CrossRef]
Foley, T.; Klinowski, J.; Meares, P. Differential Conductance Coefficients in a Cation-Exchange Membrane. Proc. R. Soc. Lond. Ser. A 1974, 336, 327–354. [Google Scholar]
Kumamoto, E.; Kimizuka, H. Transport Properties of the Barium Form of a Poly(styrenesu1fonic acid) Cation-Exchange Membrane. J. Phys. Chem. 1981, 85, 635–642. [Google Scholar] [CrossRef]
Narebska, A.; Koter, S.; Kujawski, W. Irreversible Thermodynamics of Transport across Charged Membranes. Part, I. Macroscopic Resistance Coefficients for a System with Nafion 120 Membrane. J. Membr. Sci. 1985, 25, 153–170. [Google Scholar] [CrossRef]
Narebska, A.; Kujawski, W.; Koter, S. Irreversible Thermodynamics of Transport across Charged Membranes. Part II. Ion-Water Interactions at Permeation of Alkalis. J. Membr. Sci. 1987, 30, 125–140. [Google Scholar] [CrossRef]
McCallum, C.; Meares, P. Computer prediction of stationary states of membranes from differential permeabilities. J. Membr. Sci. 1976, 1, 65–98. [Google Scholar] [CrossRef]
Krämer, H.; Meares, P. Correlation of electrical and permeability properties of ion-selective membranes. Biophys. J. 1969, 9, 1006–1028. [Google Scholar] [CrossRef] [Green Version]
Koter, S. Transport of electrolytes across cation-exchange membranes. Test of Onsager reciprocity in zero-current processes. J. Membr. Sci. 1993, 78, 155–162. [Google Scholar] [CrossRef]
Miller, D.G. Application of Irreversible Thermodynamics to Electrolyte Solutions. I. J. Phys. Chem. 1966, 70, 2639–2659. [Google Scholar] [CrossRef]
Meares, P. Coupling of ion and water fluxes in synthetic membranes. J. Membr. Sci. 1981, 8, 295–307. [Google Scholar] [CrossRef]
Dresner, L. Stability of the Extended Nernst-Planck Equations in the Description of Hyperfiltration through Ion-Exchange Membranes. J. Phys. Chem. 1972, 76, 2256–2267. [Google Scholar] [CrossRef]

Figure 1. The system: ion-exchange membrane with adjacent boundary diffusion layers. l_p—thickness of concentration polarization layer; l_m—thickness of the membrane; c′, c″—concentrations of ions in the bulk solution; c′^(m), c″^(m)—concentrations of ions at the membrane surface.

Figure 2. Electric power vs. −I calculated using Equation (5) (full ITE) for the real thickness of Nafion 120 equal to l_m = 0.3 mm, and the assumed value l_m = 0.1 mm, for 0.1–1.0 M and 0.05–3.7 M NaCl, T = 298 K. l_p—thickness of the concentration polarization layer, l_m—thickness of the membrane.

Figure 3. Pressure profile across the membrane calculated from Equation (5) (full ITE) for l_m = 0.3 mm and l_m = 0.1 mm, for 0.1–1.0 M and 0.05–3.7 M NaCl, T = 298 K. Thickness of the concentration polarization layer l_p = 0.05 mm. Thickness of the membrane—l_m.

Figure 4. Profiles of electroosmotic (el), osmotic (osm), and pressure-driven (p), components of J_v calculated from Equation (5) (full ITE) for the current density –I = 100 and 200 A/m²; one can notice a dominant role of electroosmotic contribution to J_v at higher current density.

Figure 5. Dependence of the maximum of power density, P_max, and the current density, −I, at P_max on the membrane thickness, l_m, for 0.1–1.0 M and 0.05–3.7 M NaCl, T = 298 K, Equation (5), full ITE.

Figure 6. Dependence of the maximum of power density, P_max, and the current density, −I, at P_max on the diffusion boundary layer thickness, l_p, for 0.1–1.0 M and 0.05–3.7 M NaCl, T = 298 K, Equation (5), full ITE.

Figure 7. Power density vs. current density for full model and approximate approaches: TCP—3-coefficients approach (P₁, P₂, and L12 calculated from

{\bar{t}}_{1, a p p}

, k_m, and P_s), (t_1app, P_s); TP, CP, TC—2-coefficients approach (P₁ and P₂ calculated from

{\bar{t}}_{1, a p p}

and P_s, k_m and P_s,

{\bar{t}}_{1, a p p}

and k_m, respectively); the Nafion 120 membrane separates (a) 0.1 and 1 M NaCl, (b) 0.05 and 3.7 M NaCl; l_m = 0.1 mm, l_p = 0.05 mm, T = 298 K.

Figure 7. Power density vs. current density for full model and approximate approaches: TCP—3-coefficients approach (P₁, P₂, and L12 calculated from

{\bar{t}}_{1, a p p}

, k_m, and P_s), (t_1app, P_s); TP, CP, TC—2-coefficients approach (P₁ and P₂ calculated from

{\bar{t}}_{1, a p p}

and P_s, k_m and P_s,

{\bar{t}}_{1, a p p}

and k_m, respectively); the Nafion 120 membrane separates (a) 0.1 and 1 M NaCl, (b) 0.05 and 3.7 M NaCl; l_m = 0.1 mm, l_p = 0.05 mm, T = 298 K.

Figure 8. Power density vs. current density for full model—comparison of the TC const approach with logarithmic mean with the TC approach and the full ITE for the concentrations 0.1–1 M and 0.05–3.7 M NaCl; l_m = 0.1 mm, l_p = 0.05 mm, T = 298 K.

Table 1. List of experiments, equation for 1:1 electrolyte.

Experiment	Equation
conductivity	$κ \equiv l_{m} {(\frac{I}{Δ E})}_{c, p} = l_{m} L_{E E}$
the real transport number of ion 1	${\bar{t}}_{1} \equiv {(\frac{F J_{1}}{I})}_{Δ c, Δ p = 0} = \frac{L_{π E}}{L_{E E}}$
electroosmotic coefficient	$W \equiv {(\frac{F J_{v}}{I})}_{Δ c, Δ p = 0} = {\bar{v}}_{0} {\bar{t}}_{0} + {\bar{v}}_{s} {\bar{t}}_{1} = \frac{L_{p E}}{L_{E E}}$
the apparent transport number of ion 1 (electromotive force)	${\bar{t}}_{1, a p p} = {(\frac{F Δ E}{Δ μ_{s}})}_{I, Δ p = 0} = \frac{F}{1 - {\tilde{c}}_{s} {\bar{v}}_{s}} (\frac{L_{E π}}{L_{E E}} - {\tilde{c}}_{s} \frac{L_{E p}}{L_{E E}})$ ${({\bar{t}}_{1, a p p})}_{Δ c \to 0} = {\bar{t}}_{1} - \frac{c_{s}}{c_{0}} {\bar{t}}_{0}$
electrolyte permeability	$L_{s} \equiv {(\frac{J_{s}}{Δ π})}_{I, Δ p = 0} = L_{π π} / {\tilde{c}}_{s} - L_{π p} + \frac{L_{E p} L_{π E} - L_{π E} L_{E π} / {\tilde{c}}_{s}}{L_{E E}}$
osmotic permeability	$L_{v, o s m} \equiv - {(\frac{J_{v}}{Δ π})}_{I, Δ p = 0} = L_{p p} - L_{π p} / {\tilde{c}}_{s} + \frac{L_{p E} L_{E π} / {\tilde{c}}_{s} - L_{p E} L_{E p}}{L_{E E}}$
hydrodynamic permeability	$L_{p} \equiv {(\frac{J_{v}}{Δ p})}_{I, Δ c = 0} = L_{p p} - \frac{L_{p E} L_{E p}}{L_{E E}}$

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Kujawski, W.; Yaroshchuk, A.; Zholkovskiy, E.; Koter, I.; Koter, S. Analysis of Membrane Transport Equations for Reverse Electrodialysis (RED) Using Irreversible Thermodynamics. Int. J. Mol. Sci. 2020, 21, 6325. https://doi.org/10.3390/ijms21176325

AMA Style

Kujawski W, Yaroshchuk A, Zholkovskiy E, Koter I, Koter S. Analysis of Membrane Transport Equations for Reverse Electrodialysis (RED) Using Irreversible Thermodynamics. International Journal of Molecular Sciences. 2020; 21(17):6325. https://doi.org/10.3390/ijms21176325

Chicago/Turabian Style

Kujawski, Wojciech, Andriy Yaroshchuk, Emiliy Zholkovskiy, Izabela Koter, and Stanislaw Koter. 2020. "Analysis of Membrane Transport Equations for Reverse Electrodialysis (RED) Using Irreversible Thermodynamics" International Journal of Molecular Sciences 21, no. 17: 6325. https://doi.org/10.3390/ijms21176325

APA Style

Kujawski, W., Yaroshchuk, A., Zholkovskiy, E., Koter, I., & Koter, S. (2020). Analysis of Membrane Transport Equations for Reverse Electrodialysis (RED) Using Irreversible Thermodynamics. International Journal of Molecular Sciences, 21(17), 6325. https://doi.org/10.3390/ijms21176325

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

Article Menu