 Research article
 Open Access
 Published:
On the application of the SpieglerKedem model to forward osmosis
BMC Chemical Engineering volume 1, Article number: 15 (2019)
Abstract
In Forward Osmosis the diffusion of the solute is counter to that of the solvent i.e. there is socalled “reverse salt diffusion”. Furthermore, the ratio of the two fluxes is generally taken to be a constant because of the assumption of ideal semipermeability. However with the SpieglerKedem (SK) model there is an allowance for a minor deviation from ideal semipermeability and the ratio of the solute flux and solvent flux is no longer constant. The theoretical variation of the solute flux with increasing draw solution concentration is illustrated for various degrees of deviation from ideal semipermeability. A novel variant of the SK model is also introduced and predictions compared with those obtained using the standard form. With the acceptance that the form of “breakthrough” involving cocurrent flow is impossible, a limitation is imposed upon the SK model but even with this limitation the theoretically predicted variation of solvent flux with increasing draw concentration is for certain sets of parameters of an unexpected form for minor deviation from ideal semipermeability. That intriguing counterintuitive outcomes can result from application of the SK model indicates a need to rethink its formulation of the equations and the expressions for the coefficients. This will have implications for forward osmosis and possibly reverse osmosis modelling.
Introduction
Energy production from the controlled mixing of seawater and fresh water has the advantage that it can be operated continuously unlike processes based on solar and wind energy. The main process in this area and the process that motivated this work is PressureRetarded Osmosis (PRO). The pioneers of reverse osmosis recognised the conceptual simplicity of PRO [1] but concern about costs has always been a problem [2,3,4]. The first pilot plant seeking to develop osmotic energy was opened by the Norwegian energy company Statkraft in 2009 but the pilot plant was closed in 2014 due to the low power density (power per unit of membrane area) exhibited by the prototypes and the questionable economic feasibility of the process [5]. Today it is generally accepted that if PRO is to be commercially viable then it will be necessary to use resources with a higher salinity than seawater, for example brine from a reverse osmosis desalination plant [6, 7], but the power output will still be modest. However in a recent purely theoretical paper a thought provoking theoretical analysis by Yaroshchuk [8] suggested that under some conditions and with a membrane displaying ‘leakiness’, a ‘breakthrough’ mode might occur even to the extent that with the appropriate membranes reverse solute diffusion would be eliminated. Indeed there was the prediction that with these appropriate membranes there would be cocurrent flow of both solvent and solute against a concentration gradient.
Previous exploration of the implications of the effect of minor deviations from ideal semipermeability using the SpieglerKedem (SK) model have confined themselves to a single value of the reflection coefficient [8]. In this article a range of values are used. As noted elsewhere [9], whilst irreversible thermodynamic arguments were used to derive the solute and solvent transport equations of the SK model, the membrane itself was treated as a black box". The novel variants of the SK model introduced here in the context of forward osmosis are not fundamentally based but one form does avoid the incongruity (indeed thermodynamically inappropriateness) arising from the other formulations of the SK model.
It has been suggested that the total solute flux can change sign, and that the “dramatic change in the behaviour is ultimately caused by the change in the direction of solute flow through the membrane” [8]. This specific possibility has been refuted elsewhere where it was also shown (contrary to statements in [8]) that the use of thick support layers is undesirable.
Theory
Following the SK equations given in [8], these can be simplified for FO systems by to omitting the pressure term because there should be no overall pressure drop across the barrier layer. The resulting equations are:
where c is the reference (virtual) solute concentration (as given in [10] and used in [8]), ω is the solute permeability, σ is the solute reflection coefficient, J_{s} is the solute flux, J_{v} is the solvent flux and χ is the hydraulic permeability. β is the van’t Hoff factor, R the universal gas constant and T the absolute temperature of the system. The implicit assumption in adopting this approach will be discussed later.
Now using Eq. 1 to define \( \frac{dc}{dx}/{J}_v \) and substituting the result into Eq.2 it is found that the concentration in the solute flux equation cannot be taken to be a variable because the other terms are all invariant. Consequentially all terms in (3) are fixed.
In [8] the concentration c in Eq.2 was a variable but the final term in Eq. 3 would be better written as: (1 − σ)c_{i} where c_{i} is the concentration of solute at the boundary between the support layer and the barrier layer. Consequently a sounder alternative for the barrier layer would have been:
where A and B are the standard ‘A’ parameter and ‘B’ parameter for water flux and salt flux respectively.
The consequences that follow from the adoption of Eq. 4 are explored below. The layout of the system is illustrated in Fig. 1. For the study of a potential breakthrough mode, this orientation is the one of interest [8]. The arrows show the direction of the actual fluxes but mathematically they are taken as positive in the positive xdirection which coinsides with the water flux, J_{v}. Hence the set of equations for the support, the barrier layer and the draw side boundary layer are related to each other as follows:
where Pe_{s} is the Peclet number for the support layer, Pe_{bl} is the Peclet number for the drawside layer and c_{m} is the concentration at the interface between the membrane and the draw solution. c_{d} is the bulk concentration on the draw side. The external concentration on the feedside (which is at the low concentration c_{f}) is insignificant compared with the internal concentration polarisation within the support and a separate term for the mass transfer coefficient on the feedside has not been included.
As components of Eq. 5 are referred to individually below they are listed out as:
The use of Eq. (4) is consistent with the ‘leakiness’ giving rise to a convective term that does not vary across the membrane. Initially this paper was going to have an exclusive focus upon exploring (i) the set of equations given by Eq. (5) as an alternative to those in [8]; and (ii) the effect of varying the reflection coefficient. However a new threeparameter model was briefly explored due to increasing concern about the validity of applying the SK model to FO with concentration invariant parameters. Now FO is the simplest of setting for the application of the SpieglerKedem model so when this model yields irrational results this indicates that some detail in the SK model is wrong. This might be the use of concentration invariant parameters. Given an increasing concern about whether the SK model captures the basic physics, a purely empirical equation was developed as an alternative.
Combining Eqs. (1) and (2) one obtains:
Now the concentration in the barrier layer at its interface with the draw solution is higher than that at the interface with the support layer and so for constant J_{s} Eq. (6) suggests that \( \frac{dc}{dx} \) will vary with position within the membrane (unless σ = 1) but the invariant value of the volumetric flux J_{v} with position will via Eq. (1) suggest that \( \frac{dc}{dx} \) is constant. There is thus (unless the product of concentration and hydraulic permeability, cχ is constant) an apparent contradiction at the heart of the SK model when it is applied to FO. This point illustrates the cautioning comments given elsewhere [9] and noted above. Furthermore the permeability of salt in polymers is not a constant and for sulfonated polymers, the salt diffusion coefficient increases markedly as the external salt concentration increases [11]. So notwithstanding the established nature of the SpielgerKedem (SK) model [12] it should not be seen as so superior as to be above challenge.
An additional reason for asserting that the SK model, as particularised in [8], is invalid is that it is of the form that one obtains when two resistances are in parallel. If one considers an electrical circuit where there is a current, I, flowing through two parallel resistances (R_{1} and R_{2}) driven by a potential difference V, then
Comparing Eqs. 6 and 7, the parallels are obvious. However if the second channel associated with (1 − σ) is considered to consist of pores, the pores will be in contact with the bulk fluids. In reverse osmosis the flow through such channels is sustained by pressure but in FO there is no pressure difference between the two sides and osmosis can only be sustained by permselectivity. Whilst there might be a mitigating effect due to convective flow of solvent it is poorly modelled by the above equation. Caution must be taken in the use of the SK model and it should not be overinterpreted.
The SK model is a threeparameter model well established for reverse osmosis but rarely used for FO. Previously it has been used to describe the transport across the active layer of a FO/PRO membrane [13]. A comparison with the standard twoparameter solutiondiffusion (SD) model was made and little difference was found [12]. Having noted concerns regarding the SK model, at least as it is applied to FO the scope of the paper has been extended to include a purely empirical equation. This equation (i) tends of the SD model at low flux and (ii) \( \left(\frac{J_s}{J_v}\right) \) declines with increasing J_{v} but the ratio is never positive i.e. it tends to zero and does not transition into a regime of supposed cocurrent flow of solute and solvent.
where ‘a’ would be determined from experiments. Herein the concern is to establish the influence of ‘a’ upon the fluxdraw concentration relationship.
Methodology
The results displayed in Figs. 2, 3, 4, and 5, are obtained by solving Eqs. 5 and 9. The difference c_{m} − c_{i} is eliminated by recognising (as in [8]) that the transmembrane solvent flux is also given by:
It will be readily appreciated that the four Eqs. 5a, b, c and 9 contain four unknowns: \( \left(\frac{{J}_s}{J_v}\right) \), c_{i}, c_{m} and either the flux J_{v} (if c_{d} is known) or c_{d} if flux is specified. The procedure adopted was to specify values of the solvent flux, J_{v}, for fixed values of feed salinity, c_{f}, and the relevant membrane and system parameters. This enabled one to calculate the Peclet numbers. Then equating Eq. 5a and 5b, one can find c_{i}. As the difference (c_{m} − c_{i}) can be calculated from Eq. 9 one then obtains the other boundary concentration, c_{m}. Also with a known value of c_{i}, the flux ratio, \( \left(\frac{{J}_s}{J_v}\right) \), can be found from Eq. 5b. Finally from Eq. 5c, the draw solution concentration, c_{d}, is found.
Results and discussion
Unless indicated otherwise the parameters used are the standard conditions given in Table 1. Fig. 2 illustrates the influence of the draw concentration, c_{d}, upon solvent flux, J_{v}, for the solutiondiffusion model and two versions of the SpieglerKedem model, the one introduced herein and the one given in [8]. This implementation of the SK model to Forward Osmosis has the following equation for the barrier layer:
Figure 2 has used dashed lines for those part of the predicted J_{v} c_{d}, curves that correspond with predicted cocurrent flow of solute and solvent. This is counter to the experimentally observed reverse salt flux that is characteristic of forward osmosis (FO). Herein the possibility of cocurrent fluxes in FO is refuted. For the version of the SK model introduced herein and the set of conditions indicated, the model is not valid for draw concentrations higher than 1 M. Overall this figure indicates that for physically realistic conditions there would essentially be no benefit from developing membranes that displayed a degree of ‘leakiness’ even if such a development were possible.
It has been suggested that if allowance is made for minor deviations from ideal semipermeability then operation in an overlooked mode of “breakthrough” osmosis would be possible and importantly it would yield relatively large rates of osmosis [8]. Whilst this has been refuted for realistic conditions, it is of interest to check whether there are any conditions under which minor deviations from ideal semipermeability would give higher fluxes. As shown by Eq. 5a a solute free feed (c_{f} = 0) has to have a positive value for the term (−J_{s}/J_{v}) i.e. there is countercurrent flow of solute and solvent. Thus it is of interest to compare the SpieglerKedem model with the solutiondiffusion model and with the solutiondiffusion model with B = 0 (labelled ‘ideal’) in order to understand their respective behaviours. As shown in Fig. 4, the use of the three parameter SK model with σ = 0.99 predicts fluxes that are intermediate between those of a standard solutiondiffusion model and those of the solutiondiffusion model with B = 0. Clearly it is beneficial to reduce reverse salt diffusion. In qualitative terms this is not a surprising result, but the relative importance in quantitative terms can be startling. From Eq. 5a and Eq. 5c the concentration difference across the active layer can be written as:
This equation clearly shows the important of solute diffusion, a matter which has been discuss in [14].
To conclude the exploration of the importance of the \( \left(\frac{J_s}{J_v}\right) \) term, the predictions given by Eq. 8 were analysed using standard conditions as given in Table 1 except for the ‘B’ parameter. As Eq. 8 hypothesises that the magnitude of \( \left(\frac{J_s}{J_v}\right) \) declines with increasing J_{v} the initial value of ‘B’ was 0.4 μm s^{− 1}. If it were found that the ratio \( \left(\frac{J_s}{J_v}\right) \) declines with increasing solvent flux then ‘a’ would be determined from experiments. Herein the concern is to establish the influence of ‘a’ upon the fluxdraw concentration relationship and this is illustrated in Fig. 5. It shows that whilst the term \( \left(\frac{J_s}{J_v}\right) \) can rapidly reduce to zero (but does not go negative) there is very little effect upon the fluxdraw solution relationship for the realistic set of values used.
Finally, a remark is made regarding the omission of the pressure term from the volumetric flux equation (Eq. 1). The implicit assumption is that one can apply the maximum principle for pressure [16]. When this principle is valid, the maximum and minimum pressure values are located at the domain’s boundary. Given that in FO the pressure can be taken to be the same at each boundary the application of this principle would imply that the pressure throughout the barrier layer is the same as that at the boundary. Strictly speaking, this has been assumed and has not been proven.
Conclusions
From the above theoretical analysis there is no reason to believe that it would be of any benefit to explore practically various degrees of deviation from ideal semipermeability. Whether 3 parameter models are worth pursuing is an open question but with regard to the 3 parameter SK model, it should be applied to FO and PRO with great prudence. That intriguing counterintuitive outcomes can result from its application clearly indicates a need to rethink the formulation of the equations and the expressions for the coefficients especially when it is applied outside of the sphere of RO.
The possibility of different transport behaviours in RO versus FO has been debated from as early as 1976, when Loeb [17] mentioned the need for verifying if membrane water permeance (A) was the same in RO and PRO, and a more recent review reported that the ratio of (A/B)_{FO} to (A/B)_{RO}, has been reported to be very different [18].
Nomenclature
A ‘A’ parameter for water flux (μm s^{− 1} MPa^{− 1})
a Empirically determined parameter in Eq. 8 (s μm^{− 1})
B ‘B’ parameter for salt flux (μm s^{− 1})
c solute concentration (kg m^{− 3} or M)
D diffusion coefficient (m^{2} s^{− 1})
k mass transfer coefficient (μm s^{− 1})
J_{s} solute flux (kg m^{− 2} s^{− 1})
J_{v} volumetric flux (μm s^{− 1} or m s^{− 1})
R universal gas constant (J kg^{− 1} K^{− 1})
S structural parameter (μm)
T absolute temperature (K)
Pe boundary layer (modified) Peclet number
β van ‘t Hoff factor
σ solute reflection coefficient
χ hydraulic permeability
ω solute permeability
Subscripts
b bulk
bl boundary layer on draw side
d draw side (bulk)
f feed (bulk)
i boundary between support layer and barrier layer
m surface of barrier layer on draw side
s support layer
v volumetric
Abbreviations
 SD:

Solution Diffusion model
 SK:

SpieglerKedem model
 SKW:

SpieglerKedem model
 SKY:

SpieglerKedem model
References
 1.
Loeb S. Production of energy from concentrated brines by pressure retarded osmosis. I Preliminary technical and economic correlations. J Membr Sci. 1976;1:49–63.
 2.
Loeb S, Norman RS. Osmotic power plants. Science. 1975;189:654–5.
 3.
Lee KL, Baker R, Lonsdale H. Membranes for power generation by pressureretarded osmosis. J Membr Sci. 1981;8:141–71.
 4.
Loeb S. Largescale power production by pressureretarded osmosis using river water and sea water passing through spiral wound modules. Desalination. 2002;143:115–22.
 5.
http://www.powermag.com/statkraftshelvesosmoticpowerproject/ Accessed 27 Oct 2017.
 6.
Straub AP, Deshmukh A, Elimelech M. Pressureretarded osmosis for power generation from salinity gradients: is it viable? Energy Environ Sci. 2016;9:31–48.
 7.
Lin S, Straub AP, Elimelech M. Thermodynamic limits of extractable energy by pressure retarded osmosis. Energy Environ Sci. 2014;7:276–2714.
 8.
Yaroshchuk A. “Breakthrough” osmosis and unusually high power densities in pressureretarded osmosis in nonideally semipermeable supported membranes. Sci Rep. 2017;7:45168. https://doi.org/10.1038/srep45168.
 9.
Wang J, et al. A critical review of transport through osmotic membranes. J Membr Sci. 2014;454:516–37.
 10.
Yaroshchuk AE. Osmosis and reverse osmosis in fineporous charged diaphragms and membranes. Adv Colloid Interf Sci. 1995;60:1–93.
 11.
Geise GM, Freeman BD, Paul DR. Sodium chloride diffusion in sulfonated polymers for membrane applications. J Membr Sci. 2013;427:186–96.
 12.
Spiegler KS, Kedem O. Thermodynamics of hyperfiltration (reverse osmosis): criteria for efficient membranes. Desalination. 1966;1:311–26.
 13.
Attarde D, Jain M, Gupta SK. Modeling of a forward osmosis and a pressureretarded osmosis spiral wound module using the SpieglerKedem model and experimental validation. Sep Purif Technol. 2016;164:182–97.
 14.
Field RW, Wu JJ. On boundary layers and the attenuation of driving forces in forward osmosis and other membrane processes. Desalination. 2018;429:167–74.
 15.
Field RW, Siddiqui FA, Ang P, Wu JJ. Analysis of the influence of module construction upon forward osmosis performance. Desalination. 2018;431:151–6.
 16.
Golubkin VN, Kovalev VP, Sizykh GB. Maximum principle for pressure in ideal incompressible fluid flows. TsAGI Sci J. 2016;47:599–609.
 17.
Loeb S. Production of energy from concentrated brines by pressureretarded osmosis: I. preliminary technical and economic correlations. J Membr Sci. 1976;1:49–63.
 18.
Manickam SS, McCutcheon JR. Understanding mass transfer through asymmetric membranes during forward osmosis: a historical perspective and critical review on measuring structural parameter with semiempirical models and characterization approaches. Desalination. 2017;421:110–26.
Acknowledgements
For this published article, none.
Funding
The author is grateful to EPSRC and Durham University for GCRF Seedcorn funding.
Availability of data and materials
Not applicable; all data generated or analysed during this study are included in this published article.
Author information
Affiliations
Contributions
Singleauthor paper. The author read and approved the final manuscript.
Corresponding author
Correspondence to Jun Jie Wu.
Ethics declarations
Competing interests
The author declares that he has no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/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. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Forward osmosis
 SpieglerKedem model
 Reverse osmosis
 Pressureretarded osmosis
 Renewable energy