Next Article in Journal
Runge–Kutta Numerical Method Followed by Richardson’s Extrapolation for Efficient Ion Rejection Reassessment of a Novel Defect-Free Synthesized Nanofiltration Membrane
Next Article in Special Issue
Theoretical Analysis of a Mathematical Relation between Driving Pressures in Membrane-Based Desalting Processes
Previous Article in Journal
Debye-Hückel Free Energy of an Electric Double Layer with Discrete Charges Located at a Dielectric Interface
Previous Article in Special Issue
Removal of Sulfadiazine by Polyamide Nanofiltration Membranes: Measurement, Modeling, and Mechanisms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Need for Accurate Osmotic Pressure and Mass Transfer Resistances in Modeling Osmotically Driven Membrane Processes

1
Chemical and Biochemical Procedures Laboratory, Institute of Biomolecular and Chemical Engineering, University of Pannonia, H-8200 Veszprem, Hungary
2
Department of Biophysics and Radiation Biology, Semmelweis University, Tüzoltó u. 37-47, H-1094 Budapest, Hungary
3
Rohsenow Kendall Heat Transfer Laboratory, Department of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
*
Authors to whom correspondence should be addressed.
Membranes 2021, 11(2), 128; https://doi.org/10.3390/membranes11020128
Submission received: 28 December 2020 / Revised: 27 January 2021 / Accepted: 1 February 2021 / Published: 14 February 2021
(This article belongs to the Special Issue Numerical Modelling in Membrane Processes)

Abstract

:
The widely used van ’t Hoff linear relation for predicting the osmotic pressure of NaCl solutions may result in errors in the evaluation of key system parameters, which depend on osmotic pressure, in pressure-retarded osmosis and forward osmosis. In this paper, the linear van ’t Hoff approach is compared to the solutions using OLI Stream Analyzer, which gives the real osmotic pressure values. Various dilutions of NaCl solutions, including the lower solute concentrations typical of river water, are considered. Our results indicate that the disparity in the predicted osmotic pressure of the two considered methods can reach 30%, depending on the solute concentration, while that in the predicted power density can exceed over 50%. New experimental results are obtained for NanoH2O and Porifera membranes, and theoretical equations are also developed. Results show that discrepancies arise when using the van ’t Hoff equation, compared to the OLI method. At higher NaCl concentrations (C > 1.5 M), the deviation between the linear approach and the real values increases gradually, likely indicative of a larger error in van ’t Hoff predictions. The difference in structural parameter values predicted by the two evaluation methods is also significant; it can exceed the typical 50–70% range, depending on the operating conditions. We find that the external mass transfer coefficients should be considered in the evaluation of the structural parameter in order to avoid overestimating its value. Consequently, measured water flux and predicted structural parameter values from our own and literature measurements are recalculated with the OLI software to account for external mass transfer coefficients.

1. Introduction

The osmotically driven membrane separation processes, such as pressure-retarded osmosis (PRO) and forward osmosis (FO), play a critical role in energy generation [1,2,3], water purification, and dewatering/desalination [4,5,6,7]. Consequently, these processes have been investigated experimentally [1,8,9], theoretically [10,11,12,13,14], and numerically [15,16,17] in the literature for aqueous sodium chloride solutions. Intensive research has also been carried out by hybrid PRO processes that combine reverse osmosis desalination and pressure-retarded osmosis [17,18,19,20] in order to reduce the energy consumption of water desalination. The osmotic pressure with brine solution as draw [19,21] and wastewater retentate as feed [22] has also been studied. With the increasing application of hollow-fiber modules, several studies have numerically evaluated and optimized scaled-up modules [15,16,23,24,25]. In the last decade, many experiments have been performed on the nanocomposite (used graphene oxide, nanotube, carbon quantum or other nanomaterial) membranes [26,27,28,29,30,31], which may be the harbinger of the next generation of membranes [16].
The accurate determination of the osmotic pressure is central to analyzing these processes, for key system parameters, such as water flux, power density, and membrane structural parameter, depend on osmotic pressure. The non-ideality of sodium chloride solutions complicates the evaluation of osmotic pressure. As discussed by Mistry and Lienhard [32,33], the rational activity coefficient, osmotic coefficient and osmotic pressure of NaCl and mixed electrolyte solutions can deviate significantly from those of an ideal solution (the rational activity coefficient varies from 0.67–1.2 for a molality of 0–6 molNaCl/kgwater). Nonetheless, the literature often uses the ideal van ’t Hoff linear equation (π = iCRT), valid only for dilute solutions [34], to predict osmotic pressure in aqueous solutions [2,8,12,13,16,24]. We compare the van ’t Hoff prediction for osmotic pressure to that of the OLI software (OLI Stream Analyzer 2.0 software [35]), which gives the real osmotic pressure values as a function of the solute concentration [36,37], for aqueous NaCl solutions. The error in the van ’t Hoff prediction, in some cases over 10–15% depending on the draw solution concentration, then propagates to other system parameters (water flux, power density, membrane structural parameter [38,39,40]) derived from the predicted osmotic pressure. These parameters are necessary to evaluate the separation efficiency of asymmetric [41] or thin film composite membranes [2,42,43] in PRO and FO. In particular, the membrane structural parameter S crucially determines membrane performance [8,10,12,38,40]. In addition, many papers have studied the mass transport in hollow fiber membrane modules [44,45,46,47], but transport models are not the topic of this paper.
Several reviews have also been conducted on models and methods for determining the structural parameter, S [48,49,50]. Cath et al. [2] worked on a method to standardize membrane characterization in forward osmosis. Lee et al. [38,51] showed that membrane surface porosity can significantly affect water and solute flux values, which serve as the basis for structural parameter calculations. Manickam and McCutcheon [40] discuss how the values of the predicted structural parameter can deviate by one to two orders of magnitude from those obtained through direct measurements [52]. The soft properties of the membrane layer prevent accurate direct measurements of this parameter [40]. Thus, mass transport models using fitted parameters based on measured water flux data [2,8,10,11,24,38,40,46,53,54] are used to predict S and are well accepted in the literature [12,24,46,54,55,56,57].
These models require accurate values of transport parameters, including external mass transfer coefficients, the bulk diffusion coefficient, and/or water (A) and salt (B) permeabilities. Loeb et al. [10] developed basic expressions for determining the structural parameter in PRO and FO. These expressions do not account for external concentration polarization (ECP) layers. Achilli et al. [1] extended this model to account for ECP. Yip et al. [42] later combined ECP with internal concentration polarization (ICP) to calculate solute transport. Nagy [54] and Bui et al. [12] accounted for all four transport layers (active layer, ICP, ECP on both sides of the membrane) in their transport models. Although mass transport models accounting for all layers exist, authors typically use Loeb’s model [10] for predicting the structural parameter [2,8,37,42,58,59].
Experimentally, two sets of commercially available membranes were investigated, namely, a thin-film nanocomposite membrane (NanoH2O), with inorganic nanoparticles incorporated in the selective layer [60], and the asymmetric Porifera membrane [61,62]. The main objective of these experiments was twofold. Firstly, the experiments were conducted to determine the main characteristic properties (A, B, R) of the membranes. And secondly, these characteristic properties were used to calculate values of the structural parameters to show the effect of the osmotic pressure’s prediction on their values, for comparison with data from the literature [37,40].
This paper aims to determine whether the discrepancy between van ’t Hoff’s approach and OLI software predictions in osmotic pressure leads to inaccurate prediction of the water flux, power density, membrane structural parameter, and separation efficiency for NaCl solutions or whether these errors are negligible in system analyses of PRO and FO. For the conditions under which the van ’t Hoff approach yields misleading values of system parameters, such as water flux and power density, the value of the structural parameter will be determined. The impact of external mass transfer coefficients on the predicted value of the structural parameter will also be investigated, using mass transport models that account for all transport layers [12,24,46,54,55,56,57].

2. Materials and Methods

This section is an overview of the applied membranes and the experimental, theoretical methods.

2.1. Materials

2.1.1. Membranes

Two different flat-sheet membrane layers were chosen for investigation: the Porifera (PFO) membrane, which is commonly used for FO and PRO applications (manufactured by Porifera Inc., Hayward, CA, USA); and the NanoH2O membrane, which is commonly used for Reverse Osmosis (RO) (NanoH2O Inc. Torrance, CA, USA, recently LG NanoH2O Inc., Torrance, CA, USA). The NanoH2O membrane has a thin film nanocomposite polyamide active layer containing inorganic nanoparticles. The active layer has a thickness of 500–900 nm with a polysulfone support layer. This two-layer membrane is supported by a non-woven fabric layer to strengthen the mechanical stability. The thickness of the NanoH2O flat-sheet membrane is around 200 μm.
Porifera (PFO) is an asymmetric, two-layer flat-sheet membrane, consisting of a Nomex® polymer, which is durable and heat resistant. The polymer contains aromatic groups in its molecular structure [61]; the overall thickness of the membrane was measured to be between 70–90 μm.

2.1.2. A, B and R Values of Investigated Membranes

Permeabilities can be determined experimentally for reverse osmosis (RO) or, e.g., through simultaneously fitting A, B, and S [8]. In this study, the values of A [63] and B [64], measured by conventional methods, are considered as intrinsic membrane parameters. Cross-flow experiments on RO using deionized water DI water water at velocities ranging from 0.20–0.25 m/s and pressures ranging from 3–9 bar were conducted to determine the water permeability A, solute permeability B, and percentage of solute rejection R (Table 1). Porifera and NanoH2O membranes were investigated using the PRO operating mode (draw solution faces the selective layer) and FO operating mode (draw solution faces the support layer). However, the results for Porifera membranes in FO mode are not shown in this study due to the low obtained values of R. The real cause of this behavior needs further investigations; it is not the aim of this paper. Three trials tests were carried out for both membranes. The solute rejection is determined from the concentration difference between the bulk feed and permeate salt concentration: R = 1 − Cp/Cb, while the B value is obtained by B = Jw(1 − R)/R.

2.2. Methods

2.2.1. Membrane Characterization

Microstructure and pore size distributions of PFO and NanoH2O membranes were tested by an FEI/ThermoFisher Apreo S scanning electron microscope (ThermoFisher Scientific, Waltham, MA, USA). Observations made by SEM were carried out in a low vacuum with an accelerating voltage of 2.0 and 5.0 kV, depending on the electron-beam sensitivity of the samples. In order to acquire the best analytical results for the microstructure, the samples were mounted in a methyl methacrylate casting resin (Dentacryl Technicky, SpofaDental Inc., Jicin, Czeh Republic) and cut with a PowerTome-PCZ (RMC) type ultramicrotom. Since the resin is non-conductive by nature, the resin mounted samples conductive for SEM analysis, a JEOL IB-29510VET-type carbon evaporator (JEOL, USA Inc., Ltd., Peabody, MA, USA) was used, preparing a thin electron transparent layer (approximately 15–25 nm) by carbon (Figure 1a–c,c1,e,f,f1,f2) and gold (Figure 1d) onto the sample surface.

2.2.2. Determination of the Pure Water and the Solute (NaCl) Permeability Coefficients

The pure water permeability, A, was determined by measuring the water permeate rate in the RO mode, applying hydraulic transmembrane pressure differences, ranging. The applied hydraulic pressures ranged from 2 to 10 bar, at a cross-flow velocity of 0.1–0.25 m/s, at 22 °C. In addition, the pure water permeance was also measured by a dead-end filtration method by pressurizing the liquid with an inert gas. The intrinsic salt rejection and the salt permeability coefficient were measured using 2 g/L NaCl solution, in the hydraulic pressure ranges between 2 and 9 bar. Values of A, B and R were calculated using expressions given in the literature [41,46].

2.2.3. Determination of the Osmotic Pressure

The osmotic pressure of the NaCl solution was predicted using two different methods. Firstly, the van ’t Hoff approach (π = iRTC, where R is the universal gas constant, T is the absolute temperature, and i is the van ’t Hoff factor) was used to calculate the osmotic pressure. In addition, the OLI software (OLI Stream Analyzer 2.0 software) [35] was also employed for a relative comparison. The approach uses the following curve-fit expression for the osmotic pressure, π = 5.94028C2 + 37.4521C, based on the values produced by the OLI software. Here, C is the actual solute concentration in mol/L, and π is the osmotic pressure in bar, as given by the OLI software when C > 0.6 mol/L. For values of concentration below 0.6 M, values were linearly interpolated at intervals of 0.01 g/g using tabulated OLI data shown in Ref. [35].

2.2.4. Evaluation of the Osmotic Water Flux and Solute Flux

The water flux and the reverse salt fluxes on both the NanoH2O and Porifera membranes were determined in parallel using devices from Beroplan GmbH (St. Inberg, Germany) and a custom-made, cross-flow FO set-up. The membrane areas used across the Beroplan GmbH and the in-house setup were 50 and 90 cm2 membrane surfaces, respectively. This membrane module is comprised of two geometrically similar channels on both sides of the membrane. Both liquid phases were recirculated separately through 2L-2L reservoirs, in closed loops. In certain cases (measuring the pure water permeability and the salt permeability by the RO method), a high-pressure peristaltic pump was used to circulate the draw solution. The change in the solutions’ weight was monitored over time. The concentration change in the liquid phases was measured using a conductivity meter. Experiments were performed by orienting the membrane in both PRO (active membrane layer facing the draw solution) and FO (active layer facing the feed solution) modes. The cross-flow velocity was kept between 0.2 and 0.25 m/s at 22 °C, without the presence of a hydraulic pressure difference. The draw solution concentration was chosen to be between 0.25 and 1.5 M NaCl, while DI water was used as the feed solution. The conductivity of both the solutions was monitored to determine the reverse salt flux across both membranes. The external mass transfer coefficients were predicted using data from [12].

2.2.5. Prediction of the Structural Parameter

Known theoretical expressions, i.e., Equations (13) and (15), for cases of kd→∞ and kf→∞, and the presented ones, Equations (16) and (17), were applied for prediction of the values of the structural parameter, S. These were conducted using the measured water flux data under the previously specified operating conditions. There has been no direct, acceptable measurement’s method for the determination of the structural parameter until now. The main aim of this paper is to compare the values of the structural parameter obtained by the two different methods of the osmotic pressure’s calculations and their alignment with experiments. The in-house computer program written by us for these calculations uses QuickBASIC software. The high accuracy of this computer program (with accuracy of 16 decimals) was applied for the prediction. The stepwise variations of the Jw values were compared to those measured until the relative difference between them was less than 0.01%.

3. Theoretical Development

The power density in PRO is determined by calculating the differences between the osmotic pressure and hydraulic pressure on both sides of the asymmetric membrane. The concentrations and osmotic pressures on both sides of the membrane active layer can be expressed explicitly. The individual concentration values of the solute on either side of the membrane selective layer were defined by Nagy [46,54] and Nagy et al. [36,55] for both PRO and FO operation modes. This enables the user to determine the osmotic pressure separately on both sides of the active layer, and thus the osmotic pressure difference. The concentration difference obtained using these individual interface values should be the same as that given by the conventional expressions (Equation (3) or (6) below). Knowledge of the individual osmotic pressures can be advantageous when the linear relationship between the ratio of concentration and the osmotic pressure (e.g., as assumed by Lee et al. [59]) is not valid.
The notation established for the interface concentrations and other nomenclature used during the course of this study is illustrated in Appendix A (Figure A1): Cs denotes the concentration between the selective and support layers in both operation modes (PRO and FO), while Cm denotes the interface concentration between the selective layer and draw side boundary layer (in the case of PRO) or feed side boundary layer (in the case of FO).
The interface concentrations are given by the following expressions ([46], pp. 519–521).
Case A: Pressure-retarded osmosis (PRO).
C m = C d 1 + B J w e J w S / D + 1 / k f 1 e J w / k d + C f B J w 1 e J w / k d e J w S / D + 1 / k f 1 + B J w e J w S / D + 1 / k f e J w / k d
C s = C d B J w e J w S / D + 1 / k f 1 e J w / k d + C f 1 + B J w 1 e J w / k d e J w S / D + 1 / k f 1 + B J w e J w S / D + 1 / k f e J w / k d
The concentration difference between the two sides of the selective layer is determined by subtracting Equation (2) from Equation (1). This provides the following mathematical expression for the pressure-retarded osmosis case:
Δ C m C m C s = C d e J w / k d C f e J w 1 / k f + S / D 1 + B J w e J w 1 / k f + S / D e J w / k d
Case B: forward osmosis (FO).
The interfacial concentrations are provided below ([46], pp. 447–456, [55]):
C s = C d 1 + B J w e J w / k f 1 e J w S / D + 1 / k d + C f B J w 1 e J w S / D + 1 / k d e J w / k f 1 + B J w e J w / k f e J w S / D + 1 / k d
And
C m = C d B J w e J w / k f 1 e J w S / D + 1 / k d + C f 1 + B J w 1 e J w S / D + 1 / k d e J w / k f 1 + B J w e J w / k f e J w S / D + 1 / k d
The concentration difference between the two sides of the selective layer is similar, calculated by taking the difference between Equation (4) and Equation (5):
Δ C m C s C m = C d e J w 1 / k d + S / D C f e J w 1 / k f 1 + B J w e J w 1 / k f e J w 1 / k d + S / D

3.1. Prediction of the Osmotic Pressure Difference, Δπm

The osmotic pressure difference can be determined using the interface concentrations at each side of the membrane selective layer for either PRO and FO. Accordingly, the PRO osmotic pressure difference is:
Δ π = π m C m π s C s
The OLI software determines the osmotic pressure of the NaCl aqueous solution as a function of the water activity in conjunction with an activity coefficient model:
π C = R T v ¯ l n a w
where v ¯ is the partial molal volume of water, and aw is the activity of water. The activity can be rewritten as the product of the fugacity coefficient of the solvent and its mole fraction. Knowing the concentrations on both sides of the active layer makes it possible to determine the osmotic pressure difference using Equation (7), using the concentrations given by Equations (1) and (2).
Using the van ’t Hoff approach, namely π = iRT C, where R is the universal gas constant, T is the absolute temperature, and i is the van ’t Hoff factor (number of ions). Then, by multiplying Equation (1) by iRT, the osmotic pressure at the draw side of the skin layer can be predicted via the van ’t Hoff approach.
π m C m = π d 1 + B J w e J w S / D + 1 / k f 1 e J w / k d + π f B J w 1 e J w / k d e J w S / D + 1 / k f 1 + B J w e J w S / D + 1 / k f e J w / k d
Similarly, for the other side of the active layer, the linearized osmotic pressure takes the following form:
π s C s = π d B J w e J w S / D + 1 / k f 1 e J w / k d + π f 1 + B J w 1 e J w / k d e J w S / D + 1 / k f 1 + B J w e J w S / D + 1 / k f e J w / k d
The osmotic pressure of NaCl and MgCl2 solutions as a function of concentration [4,35,46] is illustrated in Figure A2. Here, osmotic pressure has been plotted using both the van ’t Hoff approach and by the OLI Stream Analyzer. The individual data points represent the measured data. The continuous line is a fitted equation of the following form: π = 5.94028C2 + 37.4521C, where C is the actual solute concentration in mol/L, and π is the osmotic pressure in bar, as given by the OLI software. At lower solute concentrations (C ≤ 0.6 M), the tabular data were linearly interpolated at intervals of 0.01 g/g. This method provides additional resolution and accuracy in determining the osmotic pressure. Knowing the individual osmotic pressures, π(Cm) and π(Cs), the PRO osmotic pressure difference can be calculated by solving for the difference between both terms, as seen by Equation (7). For FO mode, the terms on the right-hand side of Equation (7) are reversed.
For the van ’t Hoff approach, the PRO-mode osmotic pressure difference can be calculated by taking the product of Equation (3), the universal gas constant, the van ’t Hoff factor, and the temperature (seen in Equation (11).
Δ π π m π s = π d e J w / k d π f e J w 1 / k f + S / D 1 + B J w e J w 1 / k f + S / D e J w / k d
The osmotic pressure difference can be similarly calculated for the FO operating mode, by applying the van ’t Hoff linear expression to Equation (6). The result is:
Δ π π s π m = π d e J w 1 / k d + S / D π f e J w 1 / k f 1 + B J w e J w 1 / k f e J w 1 / k d + S / D
The water flux can be obtained by calculating the product of the osmotic pressure difference and the water permeability, A [ J w = A Δ π Δ P ] The flux takes the following equation form for FO mode, when there is no hydraulic pressure difference: J w = A Δ π . Similarly, the solute flux can be determined using the concentration differences at the membrane surface and the salt permeability, B. The solute flux equation follows the same format, J s = B Δ C m , and is consistent for both operating modes (PRO and FO).
The solution procedure requires an initial guess for the water flux to compute the concentrations on both sides of the membrane. The guess is updated by recalculating the flux through the product of the osmotic pressure and the water permeability. This process is iterative and stopped when the difference between the previously iterated flux value and the updated one is less than 0.01%.

3.2. Prediction of the Membrane Structural Parameter

Determining intrinsic membrane properties such as water (A) and salt permeability (B) is generally conducted through reverse osmosis (RO) experiments that independently vary hydraulic pressure, as illustrated in the literature [8,56,57,63,64]. The value of A is obtained from the gradient of water flux vs. the hydraulic pressure difference. B is obtained by means of the rejection value of the solute component and the water flux ([46], p. 525). The structural parameter is often predicted by measured water flux data in a system run in both the pressure-retarded osmosis (PRO) and/or FO operating mode, neglecting the effect of the external mass transfer resistances [8,12,37,40]. Loeb et al. developed an expression for the prediction of the membrane structural parameter [10]. Accordingly, its value can be obtained for both the PRO and FO modes.
For PRO mode (the draw solution faces the membrane skin layer), in case of kd→∞ and kf→∞ [2,40,45,52,58,65,66,67,68,69] (Jw,∞ represents the water flux while S means the value of structural parameter, both obtained without external solute transfer resistances):
S = D J w l n B J w + A π d B + A π f
Or the water flux can be expressed from Equation (13) as:
J w = D S l n B J w + A π d B + A π f .
At a given S value, the water flux can be predicted by this relatively simple expression. The water flux predicted by Equation (14) should give the same results as those obtained using Equation (11) in the limiting case, i.e., kd→∞ and kf→∞.
For FO mode (the draw solution faces the membrane support layer), for the case of kd→∞ and kf→∞ [4,10,58,70,71]:
S = D J w l n B + A π d B + J w + A π f
Equations (13) and (15) define the value of the structural parameter as a function of more easily determinable intrinsic parameters. These include A and B as well as the water flux rate, Jw without external transfer resistances, and the diffusion coefficient in the bulk fluid phase, D. This procedure does not involve any external mass transport resistances. Since Loeb’s paper [10,58] was published, transport models have regularly taken into account one of the external boundary layer resistances, i.e., the one facing the active membrane layer [1,8,42,72]. More recently developed models also take into account both of the external fluid phase resistances [12,46,54,55,56]. These most recent transport expressions involve the effect of all four transport layers (see Figure A1). These resistances can be introduced into Loeb et al. equations resulting in the expressions in Equations (16) and (17) for the case of kd ≠ ∞ and kf ≠ ∞.
For the PRO mode, applying Equation (7) or Equation (8) with J w = A ( Δ π Δ P ) :
S = D J w l n B J w e J w / k d + A π d B + A π f e J w 1 / k f + 1 / k d
For the FO mode, applying, e.g., Equation (12) and J W = A Δ π :
S = D J w l n B + A π d B + J w e J w / k f + A π f e J w 1 / k f + 1 / k d
Equations (16) and (17) clearly show that the external resistance adjacent to the support layer, kf, can also have a significant impact in determining the structural parameter. The degree to which this effect can influence the predicted value of the structural parameter is similar in the extent to the external resistance of the active layer. Thus, neglecting this external resistance might also cause significant errors in prediction of the actual values of S, as discussed in Section 3.2.
It can be noted that relatively accurate predictions of the external mass transfer coefficients, kd, kf, as well as A, B, and D, are necessary to obtain a precise value for the structural parameter. Therefore, if these relevant parameters are calculated to a high degree of accuracy, the correctness of the applied membrane mass transfer model used for both the PRO and FO processes should provide accurate and real values for the S structural parameter. Considering the measured values of A and B as intrinsic parameters, their determined values can be regarded as relatively accurate [35]. As mentioned previously, the external mass transfer resistances are also important for determination of the salt permeability [2,12,68,71]. It can, however, be stated that the predicted kf and kd values are usually practically well-determined ones, and, therefore, their predicted values should not introduce too much error into the values of the membrane structural parameter. The value of the diffusion coefficient in the pores of the support layer remains uncertain. It is generally acceptable to assume that the diffusion coefficient value within the support layer is equal to its bulk fluid diffusion coefficient value.
Conversely, the measured values of A and B should also be independent of the RO process used for their determinations [35]. Thus, the previously experimentally determined A and B values can really be regarded as accurate. Moreover, a reliable prediction of the external mass transfer coefficients is still required for correct determination of the S value. This can be done using the measured water flux data. Taking into account the effect of the external mass transfer resistance, during determination of the salt permeability [2,12,68,71] is also critical for accuracy. The predicted kf and kd values are already well determined. Consequently, their measured/predicted values do not lead to large errors in the values of the membrane structural parameter.

4. Results and Discussion

4.1. Characterization of the NanoH2O and the Porifera (PFO) Membranes

Figure 1 presents selected micrographs of the NanoH2O (Figure 1a–c,c1,c2) and Porifera (Figure 1d–f,f1,f2) membrane layers. The NanoH2O membrane consists of three layers with the following thicknesses: a polyamide active layer (500–900 nm), a microporous polysulfone support layer (about 50 μm), and a structural support, or non-woven fabric, layer comprised of cylindrical fibers (15–25 μm) [73]. Figure 1a,b show the surface of the top (active) and bottom (non-woven fabric layer) of the NanoH2O membrane, respectively. Given the limitations of the electron microscopy we employed, the nano-scale pore size distribution on the active layer surface is not visible. Figure 1c represents the complete cross-section of this membrane, with magnified images shown in Figure 1c1,c2. The polyamide active layer contains approximately 50–100 nm thick cavities [74], reflecting its greater porosity relative to the polysulfone support layer. However, these cavities are covered and connected by a nonporous polyamide layer of 20 nm thickness, which is permeable only to water molecules [75]. The non-woven fabric layer provides the strong mechanical stability of the membrane. The three membranes layers, particularly the non-woven fabric layer, may significantly hinder solute transport relative to water transport, i.e., increase solute transport resistance against osmotic pressure difference [60]. This phenomenon is further investigated in Section 4.2.
Similarly, the PFO membrane consists of a thin and compact active layer (500–900 nm) and a thick porous support layer (70–80 μm). The active layer surface does not contain detectable pores (<20–50 nm) at micrometric or sub-micrometric scales (Figure 1d), unlike the support layer surface that contains microscopic pores (30–100 nm) (Figure 1e). Figure 1f shows the cross-section structure of the PFO membrane. Cylindrical fibers of Nomex polymer (about 50 μm thick) are detectable in the support layer. The support layer does not have the usual capillary microstructure. The wall of the capillary-like structure contains numerous holes, giving it a sponge-like character (Figure 1f2).

4.2. Osmotic Performance of Membranes

Figure 2a shows the water flux of the NanoH2O membrane measured in PRO and FO modes, with five NaCl dilutions acting as the draw solutions and DI water acting as the feed solution. The water flux in FO mode is higher at larger draw concentrations than PRO mode. This finding is unexpected, because draw solute concentration, in PRO mode, typically generates larger water flux in asymmetric membranes. However, the NanoH2O membrane is an atypical asymmetric membrane due to the nanosized particles on its selective layer. The water flux data are approximately one-order of magnitude lower than those of the cellulose triacetate membrane from [40] (see Section 4.4.2, it discussed later) and the Porifera membrane in PRO mode (Figure 3). This low water flux, despite the NanoH2O membrane’s permeability, may be caused by the fabric layer as well as the rather dense support layer hindering water transfer, thereby reducing the effect of osmotic pressure relative to the hydraulic pressure difference ([62] reached a similar conclusion). The salt fluxes of the NanoH2O membrane in both operating modes are shown in Figure 2b. The salt flux and the uncertainty of the measured flux is somewhat higher in FO mode compared to PRO mode.
Figure 3 demonstrates the osmotic water flux as a function of the draw solute concentration, with DI water on the feed side, for the Porifera membrane. Our values match those obtained by Motsa and Bamba [62]. Variations in measured water flux even at similar draw concentrations support the high error of approximately ±20%, similar to those given by [40]. The substantially higher Porifera water fluxes relative to those of the NanoH2O membrane suggest that the Porifera membrane is more suitable for use in the PRO process.

4.3. Theoretical Analyses of the Membrane Performance with Typical Parameter Values

When evaluating experimental data, to what extent does the choice of osmotic pressure prediction method affect the predicted values of the structural parameter, S? In this section, we compare results based on two prediction methods, the van ’t Hoff approach and the OLI stream analyzer, using NaCl as the solute compound.
The water flux and power density are plotted and briefly discussed as a function of draw solute concentration, hydraulic pressure difference, membrane structural parameter, water permeability, and solute permeability. Other parameter values, listed in Table 2, are held constant. The membrane parameters chosen here reflect the performance of today’s more typical CTA membranes, which differ substantially from the Porifera and NanoH2O measurements. The effect of the external mass transfer coefficients on the predicted values of the structural parameter is also shown in some figures, using data of [36], and discussed in Section 4.4.2 of the main text, using measured water flux data from [37,40], demonstrating the significance of their effect on the evaluation of the measured data.
The solution pair applied here is seawater as the draw solution (Cd = 0.6 M, i.e., ~35 g/L) and river water as the feed solution (Cf = 0.015 M, i.e., ~0.9 g/L). The effect of values Cd, kd and kf are briefly discussed for PRO and FO systems. Accurate characterization of the solute bulk diffusion coefficient is important because it proportionally affects the predicted value of the membrane structural parameter (Equations (13)–(17)). The value of the fluid diffusion coefficient varies by less than 5% across the range of relevant solute concentrations, namely, between 1.323 and 1.383 × 10−9 m2/s at 20 °C and C = 0–1.5 M, as given by [12] or 1.510–1.556 × 10−9 m2/s at C = 0–60 g/L according to [1,8,76]. The solute diffusivity was chosen to be 1.367 × 10−9 m2/s (Table 2), according to the prediction by Tow and Lienhard [77].
In the next Section 4.3.1, Section 4.3.2 and Section 4.3.3, the data of Table 2 are used for the prediction of the water flux and the effect of the membrane or operating parameters as a function of one of the varied parameters.

4.3.1. Effect of Osmotic Pressure and Transport Parameters on Water Flux and Power Density

The osmotic pressure is predicted by using both the van ’t Hoff linear approach and the OLI System Analyzer [35], which gives real values of the osmotic pressure as a function of the solute concentration. The osmotic pressure of NaCl solution as a function of concentration according to both prediction methods is shown in Figure A2 (see Appendix B). The water flux and power density are compared, using both osmotic pressure prediction methods, as a function of the draw solute concentration, hydraulic pressure difference, membrane structural parameter, solute permeability, and water permeability.

Effect of Draw Concentration on Power Density

In cases of most solute components, the osmotic pressure increases super-linearly with solute concentration [4]. This function can strongly depend on the solute component and the number of ions. For example, the osmotic pressure of the MgCl2 is significantly higher than that of sodium chloride [4], because the van ’t Hoff factor I, which is related to the number of solute ions, is equal to three for MgCl2 and two for sodium chloride. Figure A2 shows the predicted osmotic pressure using both the van ’t Hoff equation and the OLI software in cases of both NaCl and MgCl2 solute compounds. The left-hand ordinate shows the osmotic pressure of the NaCl solution, while the right-hand ordinate the MgCl2 solution. Values obtained by the van ’t Hoff approach fall on the same curve in the two cases because the ratio of their values is equal to 3/2. This figure clearly illustrates the high difference between the osmotic pressure values (more detailed discussion is given in the Appendix B). Accordingly, significantly different results can be expected of the two prediction methods when evaluating the water flux or the power density. This deviation will also affect the calculated value of the structural parameter when based upon water flux measurements.
Figure 4 illustrates the change in power density (note that according to E = Jw ΔP, 1 W/m2 power density is obtained by Jw = 10 × 10−6 m/s water flux in case of ΔP = 10 bar), as a function of the draw concentration (values of other parameters are listed in Table 2; values of A, B, the bulk diffusion coefficient were chosen according to [36]; chosen values of S and the external mass transfer coefficients, kd, kf, are typical). Let us first look at the curves obtained for NaCl solutions, in both PRO mode (draw solution facing the skin layer) and FO mode (draw solution facing the support layer). At a solute concentration greater than about 1.9 M, the predicted power density is higher when using the OLI Stream Analyzer. The difference in power densities obtained by the two prediction modes gradually increases and reaches 14% at Cd = 6 M. On the other hand, the van ’t Hoff approach predicts higher power densities at low concentrations with the seawater-river water pair (Figure A2 illustrates the osmotic pressure in the concentration range less than 0.6 M, in the upper part of this figure, on the left-hand side). In that concentration range (up to about 1.9 M draw concentration) the linear approach overestimates the osmotic pressure, which is shown in the inset of Figure A2.
Additionally, the power density is also plotted in this figure in the case of MgCl2 draw solution in PRO mode. The relative difference between the power densities of these two solutes decreases from 350% to 35%, in the concentration range of 0.6 and 3 M (see curves obtained by OLI software in Figure 1). Looking at the osmotic pressures, e.g., at Cd = 3 M, the differences between the van ’t Hoff and OLI models are about 15 and 275 bar for NaCl and MgCl2, respectively (see Figure A2). Against the huge difference in osmotic pressure, e.g., at Cd = 3 M, the difference in power densities of the two components (~35%) is rather moderate. This is the consequence of the solute transport resistances in the membrane support layer and in the external boundary layers.
The relative values of deviation are illustrated by Figure 5, as a function of the solute concentration. (The difference is related to the higher osmotic or hydraulic pressure values.) This figure well illustrates the error in osmotic pressure caused by applying the van ’t Hoff linear equation, comparing its data to those of OLI software ones. If one assumes that the acceptable error should be falling below 5%, then this is valid in a concentration range of 1 to 2.4 M. The intersection point, where the error is equal to unity, is at about 1.9 M. The deviation can reach even the 30% above 2.4 M solute concentration, where its value increases with the concentration. Figure 5 clearly illustrates the concentration range, in that the linear approach can be acceptable, perhaps between 1 and 2.4 M NaCl concentration range as draw solution. However, it should be noted that the real effect of this deviation on the value of the mass transport parameter can vary, so we consider each parameter in the following sections.

Effect of Hydraulic Pressure Difference on Water Flux and Power Density

The effect of the hydraulic pressure difference on water flux and power density is illustrated by the commonly used river water and seawater pair, applying again the data listed in Table 2 (Figure 6). The shape of these curves is well known from the literature [1]. It is well known that there is a linear relationship between the water flux, Jw, and the hydraulic pressure difference, ΔP: Jw = A(Δπ − ΔP). There exists an optimal hydraulic pressure difference that provides the maximum power density. In our case, it is important to see the differences in power density between the two prediction modes! It varies, increasing strongly with the increase in power density. Let us show the difference, e.g., at ΔP = 12 bar, where the power density has approximately its maximum value under the van ’t Hoff approach. The harvestable energy values are 7.1 and 4.9 W/m2, as obtained by the linear approach and OLI software, respectively. Thus, the deviation between the two prediction methods reaches 30%, which is significant. At higher osmotic pressures, the difference further increases. This clearly illustrates that the application of the linear or van ’t Hoff approach to describe the solute transport and membrane performance can lead to significant errors (more than 100% at lower ΔP values, not shown here) in the evaluation of transport in the river water-seawater pair.

Effect of the Membrane Structural Parameter, S on Water Flux

The membrane structural parameter is crucially important to membrane performance. It is generally the greatest mass transport resistance of the four different mass transfer layers, illustrated in Figure A1. Lowering its value is one of the most important tasks for the industrial producers of membranes. The effect of the S value on water flux is plotted in Figure 7 using the parameter values listed in Table 2. Two different ΔP values are chosen, i.e., ΔP = 0 and 10 bar. Generally, the deviation obtained by the two prediction methods increases with increasing hydraulic pressure difference, which is in harmony with data given in Figure 6. On the other hand, the data obtained by OLI software are remarkably lower than those obtained by the linear approach. Let us emphasize water flux data, e.g., at S = 500 μm and ΔP = 10, as examples. The difference in predicted water flux is about 30%. At increasing values of the membrane structural parameter, the relative error in predicted water flux increases. These results show that the usage of the van ’t Hoff linear approach can give larger errors in water flux/energy density as the structural parameter increases, even for low values of the draw concentrations, such as those paired with seawater.

Effect of the Water Permeability, A, on Power Density

In this subsection and the next subsection, we analyze the effect of water permeability and solute permeability on power density. Together, these two key parameters (so called ‘intrinsic membrane parameters’) decisively determine membrane performance. The value of the water permeability determines the power density of a PRO membrane process. Figure 8 illustrates the strong effect of A on the power density (note that according to E = Jw ΔP, 1 W/m2 power density is obtained by Jw = 1 × 10−6 m/s water flux, at ΔP = 10 bar; accordingly, the value of E in Figure 5 corresponds to that water flux). Other transport parameters are kept constant as listed in Table 2. It can be clearly seen in Figure 8 that the difference in power density, as obtained by the two prediction methods, gradually increases with water permeability. On the other hand, the relative error does not increase that much, e.g., at A = 1.9 × 10−7 m/s-bar and 10 × 10−7 m/s-bar, the relative difference of the value of E were obtained to be 29% and 23%, respectively, due to the stronger increase in the absolute value in water fluxes, which values are 1.36 × 10-6 and 4.4 × 10−6 m/s, at these A values.
We also show how the active layers’ solute concentrations, Cm and Cs, vary as a function of the water permeability. The role of the internal and external mass transfer resistances gradually increases with the increase in water flux (these are related to values of Cm/Cd and partly to value of Cs/Cd), which then can moderate the increase in the water flux or power density. This is a consequence of the water and solute transport occurring in opposite directions. Due to this, the increased water flux increases the overall resistance of the boundary layers [46,55]. Though the value of ΔCm (=CmCs) gradually decreases, the power density increases due to the more rapidly increasing value of the water permeability or water flux. On the other hand, the two prediction methods give remarkably different concentration values. It is also worth noting that the value of E is higher for the van ’t Hoff approach, though the relative value of ΔCm will be lower with increasing value of A (dotted lines of Cm/Cd and Cs/Cd in FO mode), but its absolute value gradually increases, due to the larger increase in the difference in their absolute values (not shown here).

Effect of the Solute Permeability, B, on Power Density

The solute permeability, B, also has a huge effect on the power density (Figure 9). A relatively high value was chosen for the water permeability here: A = 10 × 10−7 m/s·bar. The power density decreases strongly, with an accelerated tendency, as a function of the value of B. The difference between the two prediction methods for the calculation of the osmotic pressure here is also rather essential. The relative difference in values of E, obtained by the two prediction methods gradually increases with solute permeability, though the absolute difference between them somewhat lowers due to their lowering absolute values. At solute permeability values of B = 5 × 10−7 and 10 × 10−7 m/s, the differences in power density obtained by the prediction methods are 24% and 33%, respectively. Increasing B value lowers the concentration difference across the active membrane layer and thus decreases the water flux and power density, as illustrated in Figure 9.
Summarizing the results discussed in the above subsections, it can be stated that the van ’t Hoff equation and the OLI Stream Analyzer software give remarkably different results using the seawater and river water pair. Figure A2 excellently illustrates the difference in the predicted osmotic pressure in the lower concentration range. Accordingly, for the evaluation of the measured results in PRO systems, the OLI software is recommended if one wants to get more realistic membrane performance data on the energy-producing processes.

4.3.2. Membrane Structural Parameter, S

It seems there is no direct measuring method for determination of the value of the membrane structural parameter. Its real value depends on the thickness of the membrane, on the membrane porosity, ε, and on the membrane tortuosity, τ:S = δτ/ε.
In this subsection, measured water flux, Jw,∞, [37,40] is applied for the prediction of the value of the structural parameter, i.e., the values of S/S.
The effect of kd was calculated by measured values of water flux, for both of PRO (Figure 10) and FO (Figure 11) systems, taking into account the error of the measured data, given by Manickam and McCutcheon [40] (it is about ±20%), as well as data of Tang et al. [37]. The effect of the low concentration side mass transfer resistance (1/kf→0) was neglected here for the sake of simplification. Additionally, the effect of the external mass transfer resistances on the S/S was plotted in Figure 12. The effect of changing A and B are also illustrated (Figure 13 and Figure 14), for constant values of S. Then, the measured water flux data were used for prediction of the values of S and discussed both without and with the draw side mass transfer coefficient. Their results are listed in Table 3, Table 4, Table 5 and Table 6.

Effect of the Draw Side Mass Transfer Coefficient in PRO Mode

Let us develop a suitable expression for the prediction of the S/S values, first for a PRO system. Equations (13) and (16) clearly show that the fluid phase mass transfer coefficients can strongly affect the value of the membrane structural parameter. The ratio of these equations can define the real effect of these parameters. This ratio can then easily be obtained, e.g., for a PRO system, as (it is worth noting that the measured water flux, Jw,exp, is replaced in the following expression; its value is assumed to be equal to Jw,∞ without external mass transfer resistance, i.e., Jw,expJw,∞):
S S = l n B J w , e x p e J w , e x p / k d + π d A / B + π f A J w , e x p 1 / k d + 1 / k f l n B J w , e x p + A π d / B + A π f
where S in the numerator of Equation (18) means the value of the structural parameter of the support membrane layer in the presence of the external mass transfer coefficients, while S is that without taking into account the external mass transfer resistances, i.e., kd→∞,kf→∞, i.e., values of these mass transfer coefficients are enough high, that the value of, e.g., exp(1/kf) ≅ 1 and that of kd, as well. Equation (18) is a simple, explicit expression for prediction of the effect on the membrane structural parameter in presence of external mass transfer resistance. The value of Jw,exp simply represents the measured water flux, independent of whether or not external mass transfer resistances are present. In order to get water flux without external mass transfer resistance, such operating conditions should be provided where the assumption that kd→∞, kf→∞ is fulfilled. During our recalculation of data, the value of the structural parameter was predicted by the trial-error method (see Section 4.4). Value of S/S expressed by Equation (18) is equal to unity, when kd→∞, kf→∞, and its value gradually decreases with the increase of any value of the external mass transfer coefficients, though to different degrees. Equation (18) clearly shows that both external mass transfer resistances can affect the predicted values of S, although not in the same manner. The water flux also varies with the change of the kd and/or kf values, when other parameters are constant.
Figure 10 illustrates the effect of the draw side mass transfer coefficient on the value of S/S as a function of the water flux. The water flux was changed, taking into account its measured value (Jw = 4.9 × 10−6 m/s, dotted line in Figure 10; see also Table 4, at Cd = 1.0 M) and its error, in the range of ± 20%, i.e., between 3.8 × 10−6 and 5.5 × 10−6 m/s. As can be seen, the value of the draw side external mass transfer coefficient depends strongly on the water flux. In principle, any of the water flux data, used in this calculation, could have been the measured one. It is obvious that the value of the structural parameter varies with the change of the water flux. Values of S are listed on the upper part of this figure. Let us look at values of S/S at Jw,exp (≡ Jw,∞; dotted, vertical line). It varies between 0.95 and 0.42 when kd values are decreased from 10 × 10−5 m/s down to 1.5 × 10−5 m/s. The measured values of the kd can change between about (2 and 4) × 10−5 m/s [12,77] or even higher [1], at cross-flow velocities of 20–30 cm/s. These S/S data underscore the requirement that the measured water flux data be very carefully evaluated regarding the draw side and feed side mass transfer resistances. More precise determination of the external mass transfer coefficients is needed and one should avoid neglecting the external resistances.

Effect of the Draw Side Mass Transfer Coefficient in FO Mode

The effect of the external mass transfer coefficients under PRO mode is essentially different from that obtained in the FO system. This behavior is caused by the fact that the solute concentrations on the two sides of the active layer are much lower in FO mode due to the lower interface concentrations caused by the resistance of the membrane support layer. Accordingly, the concentration difference in the active layer can also differ from those obtained in the case of PRO systems ([46], pp. 450–453). Against that, analysis of the FO operation mode can also be useful. Let us express the value of S/S for this separation mode, using Equations (15) and (17) for FO mode:
S S = l n B + π d A / B + J w , e x p e J w , e x p / k f + π f A J w , e x p 1 / k d + 1 / k f l n B + A π d / B + J w , e x p + A π f
It is worth noting that the measured water flux is assumed to be equal to that obtained without external mass transfer resistance, i.e., Jw,expJw,∞. It is worth noting that kf represents the external mass transfer coefficient on the feed side, which faces the skin membrane layer in FO operation mode (Figure A1). S/S as a function of the draw side water flux, taking into account the about 20% error of the measured Jw,∞ value, is shown in Figure 8, considering the measured value published in [40], at Cd = 1.5 M. Other data can be found in Table 2. Values of the structural parameter are also given in this figure. Different values of the draw side fluid mass transfer coefficient are plotted (value of kd is varied between 1.5 × 10−5 m/s up to 10 × 10−5 m/s), without transfer resistance on the feed side, kf→∞. The effect of the kd external mass transfer coefficient on the S/S value is significantly less than in PRO mode. Its effect may be neglected, at lower values of kd, only. On the other hand, the water flux affects the value of the structural parameter strongly, in the water flux range investigated.

Value of S/S as a Function of Draw Side Mass Transfer Coefficient by Both of PRO and FO Operation

How S/S varies as a function of the draw side mass transfer coefficient is shown in Figure 12, applying parameters’ data of Manickam and McCutcheon again [40], at Jw,exp = 4.9 × 10−6 m/s, Cd = 1.0 M, S = 803 μm for PRO and Jw, = 3.0 × 106 m/s; Cd = 1.0 M; S = 434 μm, for FO operation (see Table 5). Other constant parameter values are given in the caption of Figure 10. Generally, it can be stated that the effect of the external mass transfer coefficients can practically be neglected when values of these coefficients are larger than 10 × 10−5 m/s. Thus, e.g., the value of S/S is equal to 0.92 for PRO and 0.97 for FO mode, at kd = 10 × 10−4 m/s, while there is no resistance on the permeate side. On the other hand, S/S quickly decreases with the kd mass transfer rate between 1 × 10−5 m/s and 10 × 10−5 m/s (note here that the draw side mass transfer resistance faces the membrane support layer). It is worth noting that the value of the structural parameter without external resistance, S∞ and Jw,exp were kept constant during the calculation; thus, in order to reach the measure water flux, the actual value of S should gradually be lowered. Zero value of S/S, at the given external mass transfer resistances, means that, to reach the value of Jw,exp, the diffusion path of the membrane support layer, S, should be equal to zero. The sum of the external fluid layers and the active layer mass transfer resistances reaches the maximum value of the possible transport resistance. This figure also confirms that the effect of the draw side mass transfer resistance is much stronger; therefore, its effect should be taken into account in order to avoid essential error in the predicted value of the structural parameter. It can also be seen that the relative deviation of results between the two modes, with and without permeate side mass transfer resistance, continuously increases in PRO mode. Thus, the relative deviations between the two modes of osmotic pressure calculation, e.g., between kd = 2 × 10−5 m/s and kf = kd or kf→∞, as well as kd = 1.5 × 10−5 m/s and kf = kd or kf→∞, are 14% and 28%, respectively. The continuous lines represent data with external transfer resistance on both sides of the membrane (kd = kf), while the broken lines show data without external transfer resistance on the feed side.
Working under the usual operating conditions, i.e., with cross-flow velocity at about 0.2 and 0.5 m/s, the external mass transfer coefficients are mostly falling in the mass transfer coefficient range of 2 × 10−5 and about 8 × 10−5 m/s [1,12,77]. This should mean that the predicted value of the structural parameter will be significantly higher than its real values, which are obtained taking into account the external mass transfer resistances, at a given value of the water flux. The overestimation of the structural parameter then can lead to a faulty evaluation of membrane properties. Accordingly, the accurate prediction of values of the external mass transfer coefficients is crucially important for the correct evaluation of the membrane properties.

4.3.3. The Effect of the External Mass Transfer Coefficient as a Function of A and B

Considering that both the water and solute permeability can significantly affect the water flux and, consequently, the effect of the external mass transfer resistances, as it is plotted in Figure 10, Figure 11 and Figure 12, the change of the water and the solute permeability were also studied and briefly discussed here. Figure 13 shows the variation of the S/S value as a function of the water permeability, using the parameters in the caption of Figure 10 and at Cd = 1.0 M; Sme s= S = 803 µm; kf→∞. The value of S/S has a minimum as a function of A; the minimum occurs at larger A values as the external mass transfer coefficient decreases. The reason for this behavior in the value of S/S, in the A range of (0.1–1) × 10−7 m/s·bar, is that the value of Jaw/kd increases by about an order of magnitude, and accordingly the increase in the numerator’s value is much less than that of the denominator in Equation (18), due to the increasing value of the second term of the right-hand side in this expression. On the other hand, with further increase in the value of A beyond 1 × 10−7 m/s·bar, the value of πdA in the numerator also gradually increases, pushing the value of S/S up to unity. Values of S/S vary significantly across the practical value range of the water permeability, i.e., between 1 × 107 m/s·bar and 8 × 10−7 m/s·bar. This figure clearly shows the importance of the A value in the real predicted values of the S/S in the presence of an external mass transfer resistance. The water flux also strongly increases in the function of the A value (dotted line), as is expected.
How the solute permeability affects the predicted values of S/S is illustrated by Figure 14 at Cd = 1.0 M; S = 803μm, kf∞ (other constant parameter values as in the caption of Figure 10). A great variation in values of S/S∞ was obtained in Figure 14, especially at low values of the solute permeability. The external draw side mass transfer coefficient significantly affects the predicted value of S/S, in this value range of B. The water flux also changed as a function of the B value: its value lowers, close to linearly, with the increasing value of the solute permeability parameter. Comparing the change of A values in range between 1 × 10−7 m/s·bar and 100 × 10−7 m/s·bar (Figure 13), the change of the water flux is much higher than those in case of B, in the range of (1–100) × 10−7 m/s.
The data discussed in this and previous sections show that the change in value of all parameters can affect the measure of the effect of the external mass transfer resistance and on the predicted S/S. values. Accordingly, neglecting their effects without any critical considerations can cause an essential error in predicted S values. This should suggest that a prediction of the structural parameter should always carefully consider the possible effect of the external mass transfer rates, which depend mainly on the cross-flow velocity.

4.4. Evaluation of Measured Data Focusing on Values of the Membrane Structural Parameter

The value of the structural parameter cannot be directly measured. Thus, solute transport models are generally used to predict the S values. These models depend on measured water flux; membrane parameters A, B, and D; external mass transfer coefficients, kd, kf; and solute concentrations, or osmotic pressure π©. This section uses experimental water flux data measured by the authors and those published by Manickam and McCutcheon [40] and Tang et al. [37] in PRO and FO systems to investigate the effect of the commonly neglected external mass transfer coefficients on the structural parameter. We conclude that accounting for external mass transfer coefficients is necessary to accurately determine the structural parameter.

4.4.1. Prediction of S Values Based on Measured Data Using NanoH2O and Porifera Membranes

The correct prediction of the structural parameter of the support layer is extremely important for the correct prediction of membrane performance, which gives essential information on the membrane’s quality. That is why not only experimental data are analyzed but detailed discussion of typical theoretical operating conditions will be shown. Accordingly, this paper offers a surveying evaluation of the membrane performance for both the PRO and FO systems. When evaluating experimental data, to what extent does the choice of osmotic pressure prediction method affect the predicted values of the structural parameter, S? In this section, we compare results based on two prediction methods, namely, the van ’t Hoff approach and the OLI stream analyzer, using NaCl as the solute.
Table 3 lists the predicted data of the structural parameter values of NanoH2O membrane using Equation (13) or Equation (16) for PRO and Equation (15) or Equation (17) for FO, in cases of kd = kf∞ and kd = kf = 2 × 10−5 m/s, respectively, at the different draw concentrations. Values of kd = kf = 2 × 10−5 m/s are regarded as actual ones for cross-flow velocity between 0.2 and 0.25 m/s, taking into account the measured data of [12]. It is remarkable the rather low values of the water flux compared to the typical commercial membranes, e.g., cellulose triacetate or thin-film composite membranes [12,42,75]. This low water flux data can be caused by the structure of the membrane. This may be caused partly by the compact character of the support layer and partly by the non-woven fabric layer, which can hinder the transport of the water molecules. This conclusion is in harmony with [62]. On the one hand, the FO system serves somewhat higher values of water flux at larger draw solute concentration. On the other hand, the rejection of the solute component (90.2%R; see Table 1) is remarkable in case of the FO operating mode.
As can be seen in Table 3, values of the structural parameter are unusually high, especially in the PRO system. On the other hand, values of S are significantly lower in the FO system than those in case of the PRO one. In the case of the PRO operating mode, these S values are surprisingly close two times higher than those obtained for PRO membranes, which are mostly below 1000 μm [37,40]. This large discrepancy is obviously caused by the compact structure of the support layer and also the non-woven fabric layer. Further detailed investigation is needed to clarify exactly the cause of this behavior.
The predicted values of the structural parameter for the Porifera membrane are listed in Table 4. All measured points, plotted in Figure 3, are taken into account. The average values do not involve the highest and the lowest measured points. Remarkably, the average data obtained by the van ’t Hoff approach and the OLI software give minor deviations. This difference is less than 4–5%. The average values of the structural parameter are also rather high, more than 1000 μm. Basically, these values are similar to those predicted in the literature [8,37,40].

4.4.2. Prediction of S Values Based on PRO and FO Data in Manickam and McCutcheon

Manickam and McCutcheon [40] investigated the effect of the draw side concentration on the structural parameter for an FO system (ΔP = 0) operating in PRO (draw solution faces the active membrane layer) and FO (draw solution faces the membrane support layer) modes for two types of membranes. In our analysis, we consider cellulose triacetate membrane from Hydration Technology, USA, (HTI-CTA). Values of A and B were determined experimentally by the standard methodologies using RO. Several studies predict the S value of this membrane without accounting for external mass transfer resistances: Tiraferri et al. [8] predicted S = 481–541 µm, Kim et al. [9] S = 503–560 µm, Cath et al. [4] S = 450–560 µm, Achilli et al. [1] S = 675 µm, and Kim et al. [9] predicted S = 560–590 µm. The value of S for the HTI-CTA membrane, when accounting for external mass transfer resistances, is lower than that without resistance at a given measured value of Jw. For example, Bui et al. [12] predicts S = 200–500 µm depending on the cross-flow velocity. The error in S increases as the external mass transfer resistances increase (see Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11).
An iterative method using stepwise interpolation (0.01 g/g intervals) to determine C vs. π(C) at low concentrations (C < 0.6 M) and the expression π(C) = 5.94028 C2 + 37.4521 C for high concentrations (C > 0.6 M) are employed, for OLI prediction, to recalculate water flux without external mass transfer resistances, Jw,∞, and with external mass transfer resistances, Jw. Our objective is to ensure that our calculated water flux values are within 0.01% accuracy of the published water flux data. The S value is then fitted to these values of water flux using Equations (7) and (10) or (13). The S/S values are also calculated using Equation (18) for PRO and Equation (19) for FO (another approach is to calculate the S value using Equation (13) for PRO or Equation (15) for FO in order to predict water flux from Equations (7) and (11) or (12)).
The experimental water flux data and predicted structural parameter data using the OLI software can be found in Table 5. For the sake of completeness, we also have predicted the values of structural parameter, applying the van ’t Hoff approach, as well and its error related to the OLI data. The deviation between results obtained by the two prediction methods is relatively high. Its value gradually decreases with an increase in the draw concentration, in harmony with the tendency of the osmotic pressure function, plotted in Figure 5. The predicted values of S (kd = kf∞) can deviate from those published by Manickam and McCutcheon, because S is sensitive to changes in water flux. The average value of S from Manickam and McCutcheon [40], 793 μm, matches the average value determined in this study, 794 μm, for PRO systems. The external resistance values, for kd = 5.0 × 10−5 and kd=2.0 × 10−5 m/s, are chosen based on those measured in the literature. (Bui et al. [12] finds kd = 1.74–1.84 × 105 and 2.0–2.1 × 105 m/s for 0–1.5 M NaCl concentrations and kf = 1.8–2.1 × 105 m/s for DI water with cross-flow velocities between 0.21 and 0.31 m/s. Manickam and McCutcheon [40] similarly measured kd = 2.0 × 105 m/s at 0.26 m/s cross-flow velocity. In contrast, other studies predicted or measured higher kd values, e.g., 4.8 × 105 m/s [77] or 8.5 × 105 m/s [1]).
For various water fluxes and draw concentrations, the ratio of our predicted S for these resistances relative to S ranges from 0.45 to 0.92 in PRO (Columns 5–6 in Table 5) and 0.77 to 0.93 in FO (Columns 5–6 in Table 5). These results illustrate that the external mass transfer coefficient can significantly impact the predicted values of the structural parameter, depending on the draw concentration. For example, at a water flux of 2.8 × 10−6 m/s and a draw concentration of 0.5 mol/L in PRO, the structural parameter, with kd = 2.0 × 10−5 and kd = 5.0 × 10−5 m/s, is 649 and 356 μm, respectively. These values are 82% and 45% of the structural parameter without external mass transfer resistance, which is 794 μm. Figure 10, Figure 11 and Figure 12 match the results in Table 5. The increase in the internal polarization layer’s resistance, i.e., an increase in S, may significantly reduce the effect of the external mass transfer coefficients. S/S varies between 0.2 and 0.8 for S values of 636 to 1000 µm, respectively, at kd = kf = 2.0 × 10−5 m/s and for a decreasing Jw,∞ (Figure 10). As Jw,∞ increases, S/S increases at approximately constant S (Table 5). Accordingly, the effect of the external polarization layer may depend on mass transport (kd, kf, S) and membrane (A,B) parameters (see Equations (3) and (7)). Consequently, neglecting this effect requires careful consideration.

4.4.3. Prediction of S Values Based on PRO and FO Data in Tang et al.

Tang et al. [37] uses an FO system (ΔP = 0) operating in PRO and FO modes. They apply the van ’t Hoff linear approach for predicting osmotic pressure and mass transfer models without external mass transfer resistances for predicting structural parameter. However, Equations (14) and value of Jaw expressed from Equation (15) should have been applied rather than Equations (7), (11) and (12) used in the literature, which are used by multiplying A for the calculation of the water fluxes, for the prediction of the water flux without external mass transfer resistances, Jw,, at a given value of S. Knowing the value of Jw,∞, the value of S/S was predicted by Equations (18) or (19). We have re-evaluated the measured water flux based on osmotic pressures on both sides of the active membrane layer that were calculated using the OLI Stream Analyzer software [35]. Equations (16) and (17), accounting for the external mass transfer resistances, are used to predict the membrane structural parameters. Table 6 includes S values from Tang based on van ’t Hoff and our results based on OLI for various draw concentrations and measured water fluxes in PRO and FO. Results in Table 6 match with data shown in Figure 10 for PRO mode, and Figure 11; Figure 12 in FO mode. Our calculated S values without external resistances differ from those obtained by Tang et al., depending on draw solute concentration. This difference decreases with increasing draw concentration, in line with osmotic pressure deviations. For example, the S values largely vary from one another at the lower Cd = 0.5 M, due to differences in osmotic pressure calculations using van ’t Hoff (24.3 bar) and OLI (20.2 bar). Beginning at approximately Cd = 1.5 M, van ’t Hoff underestimates osmotic pressure, i.e., π(OLI) > π(van ’t Hoff). The effect of external mass transfer resistance on structural parameter similarly depends on Cd, as S/S increases with increasing draw concentration. Consistent with findings in Sec. 4.5.1, neglecting external resistance for external mass transfer coefficients larger than 4–5 × 10−5 m/s causes less error in the prediction of the membrane structural parameter.
We clearly show that the van ’t Hoff approach and the OLI software can yield different results, especially in the case of seawater–river water pairing, which has a low draw side concentration. This deviation strongly depends on the values of all transport parameters affecting solute and water transport rates.

5. Conclusions

The membrane structural parameter is a key indicator of membrane performance. This study investigates how the structural parameter calculation is affected by the osmotic pressure calculation and by the external mass transfer coefficients. Our analysis shows the importance of the careful application of these transport parameters in order to accurately predict membrane performance. We show that the van ’t Hoff approach does not agree with results based on the true osmotic pressure over a range of solute concentration. This difference can strongly affect the measured water flux and the predicted membrane performance, especially in the case of seawater-river water pairing, for which the draw side concentration is low. Accordingly, we recommend that the van ’t Hoff linear approach is not used to determine the osmotic pressure. In addition, the external mass transfer coefficients should not be neglected in predicting structural parameter, as is typical in the literature, unless their values are greater than approximately 8–10 × 10−5 m/s. Failing to consider these resistances can lead to the overestimation of the membrane structural parameter at typical cross-flow velocities of 20–35 cm/s in PRO and FO. Consequently, the accurate prediction of external mass transfer coefficients is crucially important for the accurate evaluation of the membrane properties.

Author Contributions

E.N. performed the model development, formal analysis, and writing. I.H. performed visualization and formal analysis. D.R., Q.J.W., Y.D.A., and J.H.L. provided conceptual guidance and editing of the manuscript. D.R. also performed formal analysis. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Hungaryan National Development Agency grants OTKA 116727 and NKFIH-1158/2019. This research did not receive any specific grant from funding agencies in the public, commercial or not-for-profit sector.

Acknowledgments

Thanks to the Hungarian National Development Office.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Awater permeability, m/(s bar)
Bsalt permeability, m/s
Csalt concentration, mol/L
Dfluid diffusion coefficient, m2/s
Epower density, W/m2
Jssolute transport rate, kg/m2s
Jwwater flux, m/s
Jw,∞water flux at kd = kf→∞, m/s
kdiffusive mass transfer coefficient, m/s
Phydraulic pressure, bar
Rrejection coefficient
Sstructural parameter, m
Sstructural parameter without external mass transfer resistances, m
vconvective fluid velocity perpendicular to the membrane surface, m/s
Greek
βdiffusive plus convective transport coefficient, m/s
πosmotic pressure, bar
δthickness of the fluid boundary layer, m
εporosity,
τtortuosity,
ΔCmconcentration difference on the active layer, mol/L
ΔPhydraulic pressure difference, bar
Δπeffective osmotic pressure difference on the active layer, bar
Subscripts
bbulk solute concentration
ddraw solution
ffeed solution
mmembrane active layer
iinterface between selective and the support layer (with P, π, β, or C)
ppermeate
ssupport layer (with S, τ, ε, δ, or D)
spexternal interface of the support layer
wwater

Appendix A. Concentration Distribution and Nomination of Transport Parameters in PRO and FO

Figure A1. Illustration of the concentration distribution and nomenclature of parameters for PRO (left-hand side figure) and FO (figure on the right-hand side) processes, including external mass transfer resistances. Cm and Cs correspond to the interface concentration of the active layer facing the boundary layer and facing the porous support layer, respectively. Δπeff is the effective osmotic pressure difference on the active layer.
Figure A1. Illustration of the concentration distribution and nomenclature of parameters for PRO (left-hand side figure) and FO (figure on the right-hand side) processes, including external mass transfer resistances. Cm and Cs correspond to the interface concentration of the active layer facing the boundary layer and facing the porous support layer, respectively. Δπeff is the effective osmotic pressure difference on the active layer.
Membranes 11 00128 g0a1
The interface concentrations and the external mass transfer coefficients in the PRO and FO membrane processes are illustrated in Figure A1. The subscripts d and f correspond to the draw and the feed side, respectively. The interface concentrations of the asymmetric membrane are Cm, Cs and Csp; those of the fluid phase are Cd and Cf; and the diffusive mass transfer coefficients of the membrane transfer layers are kd, B, ksp and kf, while the corresponding convective ones are βd, βsp, βf.

Appendix B. Osmotic Pressure as a Function of Draw Concentration

The osmotic pressure as a function of NaCl (left-hand vertical axis) and MgCl2 (right-hand vertical axis) solute concentrations is illustrated in Figure A2. Differences in the osmotic pressure values based on the van ’t Hoff linear approach and OLI Stream Analyzer [4,35,46] are shown. The deviation between the two osmotic pressure curves increases as a function of draw concentration, especially in the case of MgCl2. For example, the osmotic pressure difference between the van ’t Hoff and OLI data is about 284 bar at C = 3 M [4]. Though the difference is lower for NaCl, particularly when C < 2 M, it is still not negligible, as shown in this study.
Figure A2. Osmotic pressure, predicted by the van ’t Hoff linear approach and the OLI Streamline Analyzer software, as a function of the solute concentration for NaCl and MgCl2. The vertical axis corresponds to the osmotic pressure.
Figure A2. Osmotic pressure, predicted by the van ’t Hoff linear approach and the OLI Streamline Analyzer software, as a function of the solute concentration for NaCl and MgCl2. The vertical axis corresponds to the osmotic pressure.
Membranes 11 00128 g0a2

Appendix C. Estimation of the Uncertainty of Predicted Data

This study used an in-house computer program written in Quick-Basic computer language. Measured data from the literature were used to calculate the water flux and power density as a function of different membrane parameters. The error in measured water flux from the literature can reach up to 10%. We estimate a measurement error of 20% to ensure that the reading error does not increase the uncertainty of our calculations. Accuracy of our calculated water flux data is 16 decimals using Equations (7) and (11) for PRO and Equation (12) for FO; it is fitted to the measured data within 0.01% accuracy for all calculations of structural parameter. The structural parameter was iteratively calculated using Equations (16) or (17), and/or (18) and (19) until the error was less than 0.1%.
Uncertainty of the PRO and FO experimental data: every data point for both the water flux and the solute concentration was measured in triplicate; also, some of these triplicate measurements were repeated. The average values of data were then plotted in Figure 2; Figure 3. Generally, it can be stated that the scatter of these data reached the value of ±20%. The error of the results of RO experiments for determination of the solute permeability, applying triplicate measurements again, was obtained to be ±40%.

References

  1. Achilli, A.; Cath, T.Y.; Childress, A.E. Power generation with pressure retarded osmosis: An experimental and theoretical investigation. J. Membr. Sci. 2009, 343, 42–52. [Google Scholar] [CrossRef]
  2. Cath, T.Y.; Elimelech, M.; McCutcheon, J.R.; McGinnis, R.L.; Achilli, A.; Anastasio, D.; Brady, A.R.; Childress, A.E.; Farr, I.V.; Hancock, N.T.; et al. Standard Methodology for Evaluating Membrane Performance in Osmotically Driven Membrane Processes. Desalination 2013, 312, 31–38. [Google Scholar] [CrossRef] [Green Version]
  3. Altaee, A.; Sharif, A. Pressure retarded osmosis: Advancement in the process applications for power generation and desalination. Desalination 2015, 356, 31–46. [Google Scholar] [CrossRef]
  4. Cath, T.; Childress, A.; Elimelech, M. Forward osmosis: Principles, applications, and recent developments. J. Membr. Sci. 2006, 281, 70–87. [Google Scholar] [CrossRef]
  5. McGovern, R.K.; Lienhard, J.H. On the potential of forward osmosis to energetically outperform reverse osmosis desalination. J. Membr. Sci. 2014, 469, 245–250. [Google Scholar] [CrossRef] [Green Version]
  6. Shaffer, D.L.; Werber, J.R.; Jaramillo, H.; Lin, S.; Elimelech, M. Forward osmosis: Where are we now? Desalination 2015, 356, 271–284. [Google Scholar] [CrossRef]
  7. Mohammadifakhr, M.; De Grooth, J.; Roesink, H.D.W.; Kemperman, A.J.B. Forward Osmosis: A Critical Review. Processes 2020, 8, 404. [Google Scholar] [CrossRef] [Green Version]
  8. Tiraferri, A.; Yip, N.Y.; Straub, A.P.; Castrillon, S.R.-V.; Elimelech, M. A method for the simultaneous determination of transport and structural parameters of forward osmosis membranes. J. Membr. Sci. 2013, 444, 523–538. [Google Scholar] [CrossRef]
  9. Kim, D.I.; Kim, J.; Hong, S. Changing membrane orientation in pressure retarded osmosis for sustainable power generation with low fouling. Desalination 2016, 389, 197–206. [Google Scholar] [CrossRef]
  10. Loeb, S. Effect of porous support fabric on osmosis through a Loeb-Sourirajan type asymmetric membrane. J. Membr. Sci. 1997, 129, 243–249. [Google Scholar] [CrossRef]
  11. Wang, J.; Dlamini, D.S.; Mishra, A.K.; Pendergast, M.T.M.; Wong, M.C.; Mamba, B.B.; Freger, V.; Verliefde, A.R.; Hoek, E.M. A critical review of transport through osmotic membranes. J. Membr. Sci. 2014, 454, 516–537. [Google Scholar] [CrossRef]
  12. Bui, N.-N.; Arena, J.T.; McCutcheon, J.R. Proper accounting of mass transfer resistances in forward osmosis: Improving the accuracy of model predictions of structural parameter. J. Membr. Sci. 2015, 492, 289–302. [Google Scholar] [CrossRef] [Green Version]
  13. Cheng, Z.L.; Chung, T.-S. Mass transport of various membrane configurations in pressure retarded osmosis (PRO). J. Membr. Sci. 2017, 537, 160–176. [Google Scholar] [CrossRef]
  14. Anissimov, Y. Aspects of Mathematical Modelling of Pressure Retarded Osmosis. Membranes 2016, 6, 13. [Google Scholar] [CrossRef]
  15. Kishimoto, M.; Tanaka, Y.; Yasukawa, M.; Goda, S.; Higa, M.; Matsuyama, H. Optimization of Pressure-Retarded Osmosis with Hollow-Fiber Membrane Modules by Numerical Simulation. Ind. Eng. Chem. Res. 2019, 58, 6687–6695. [Google Scholar] [CrossRef]
  16. Chae, S.H.; Kim, Y.M.; Park, H.; Seo, J.; Lim, S.J.; Kim, J.H. Modeling and Simulation Studies Analyzing the Pressure-Retarded Osmosis (PRO) and PRO-Hybridized Processes. Energies 2019, 12, 243. [Google Scholar] [CrossRef] [Green Version]
  17. Wang, Q.; Zhou, Z.; Li, J.; Tang, Q.; Hu, Y. Investigation of the reduced specific energy consumption of the RO-PRO hybrid system based on temperature-enhanced pressure retarded osmosis. J. Membr. Sci. 2019, 581, 439–452. [Google Scholar] [CrossRef]
  18. Kurihara, M.; Takeuchi, H. SWRO-PRO System in “Mega-ton Water System” for Energy Reduction and Low Environmental Impact. Water 2018, 10, 48. [Google Scholar] [CrossRef] [Green Version]
  19. Lee, S.; Choi, J.; Park, Y.-G.; Shon, H.; Ahn, C.H.; Kim, S.-H. Hybrid desalination processes for beneficial use of reverse osmosis brine: Current status and future prospects. Desalination 2019, 454, 104–111. [Google Scholar] [CrossRef]
  20. Bargiacchi, E.; Orciuolo, F.; Ferrari, L.; Desideri, U. Use of Pressure-Retarded-Osmosis to reduce Reverse Osmosis energy consumption by exploiting hypersaline flows. Energy 2020, 211, 118969. [Google Scholar] [CrossRef]
  21. Altaee, A.; AlZainati, N. Novel Thermal Desalination Brine Reject-Sewage Effluent Salinity Gradient for Power Generation and Dilution of Brine Reject. Energies 2020, 13, 1756. [Google Scholar] [CrossRef] [Green Version]
  22. Yang, T.; Wan, C.F.; Xiong, J.Y.; Chung, T.-S. Pre-treatment of wastewater retentate to mitigate fouling on the pressure retarded osmosis (PRO) process. Sep. Purif. Technol. 2019, 215, 390–397. [Google Scholar] [CrossRef]
  23. Sanahuja-Embuena, V.; Khensir, G.; Yusuf, M.; Andersen, M.F.; Nguyen, X.T.; Trzaskus, K.; Pinelo, M.; Nielsen, C. Role of Operating Conditions in a Pilot Scale Investigation of Hollow Fiber Forward Osmosis Membrane Modules. Membranes 2019, 9, 66. [Google Scholar] [CrossRef] [Green Version]
  24. Chen, Y.; Loh, C.H.; Zhang, L.; Setiawan, L.; She, Q.; Fang, W.; Hu, X.; Wang, R. Module scale-up and performance evaluation of thin film composite hollow fiber membranes for pressure retarded osmosis. J. Membr. Sci. 2018, 548, 398–407. [Google Scholar] [CrossRef]
  25. Chen, Y.; AlAnezi, A.A.; Zhou, J.; Altaee, A.; Shaheed, M.H. Optimization of module pressure retarded osmosis membrane for maximum energy extraction. J. Water Process. Eng. 2019, 32, 100935. [Google Scholar] [CrossRef]
  26. Yadav, S.; Saleem, H.; Ibrar, I.; Naji, O.; Hawari, A.A.; AlAnezi, A.A.; Zaidi, S.J.; Altaee, A.; Zhou, J. Recent developments in forward osmosis membranes using carbon-based nanomaterials. Desalination 2020, 482, 114375. [Google Scholar] [CrossRef]
  27. Lim, S.; Park, M.J.; Phuntsho, S.; Mai-Prochnow, A.; Murphy, A.B.; Seo, D.; Shon, H.K. Dual-layered nanocomposite membrane incorporating graphene oxide and halloysite nanotube for high osmotic power density and fouling resistance. J. Membr. Sci. 2018, 564, 382–393. [Google Scholar] [CrossRef]
  28. Akther, N.; Phuntsho, S.; Chen, Y.; Ghaffour, N.; Shon, H.K. Recent advances in nanomaterial-modified polyamide thin-film composite membranes for forward osmosis processes. J. Membr. Sci. 2019, 584, 20–45. [Google Scholar] [CrossRef]
  29. Gai, W.; Zhao, D.L.; Chung, T.-S. Novel thin film composite hollow fiber membranes incorporated with carbon quantum dots for osmotic power generation. J. Membr. Sci. 2018, 551, 94–102. [Google Scholar] [CrossRef]
  30. Park, M.J.; Lim, S.; Gonzales, R.R.; Phuntsho, S.; Han, D.S.; Abdel-Wahab, A.; Adham, S.; Shon, H.K. Thin-film composite hollow fiber membranes incorporated with graphene oxide in polyethersulfone support layers for enhanced osmotic power density. Desalination 2019, 464, 63–75. [Google Scholar] [CrossRef]
  31. Ren, J.; Chowdhury, M.R.; Qi, J.; Xia, L.; Huey, B.D.; McCutcheon, J. Relating osmotic performance of thin film composite hollow fiber membranes to support layer surface pore size. J. Membr. Sci. 2017, 540, 344–353. [Google Scholar] [CrossRef]
  32. Mistry, K.H.; Lienhard, J.H. Effect of Nonideal Solution Behavior on Desalination of a Sodium Chloride Solution and Comparison to Seawater. J. Energy Resour. Technol. 2013, 135, 042003. [Google Scholar] [CrossRef]
  33. Mistry, K.H.; Hunter, H.A.; Lienhard, J.H. Effect of composition and nonideal solution behavior on desalination calculations for mixed electrolyte solutions with comparison to seawater. Desalination 2013, 318, 34–47. [Google Scholar] [CrossRef] [Green Version]
  34. Ives, D.J.G. Chemical Thermodynamics: With Special Reference to Inorganic Chemistry; (University chemistry series); Macdonald Technical and Scientific: London, UK, 1971. [Google Scholar]
  35. OLI Stream Analyzer 2.0, OLI Systems Inc. Available online: https://silo.tips/download/analyzer-20-installation-guide-contents (accessed on 13 February 2021).
  36. Nagy, E.; Hegedus, I.; Tow, E.W.; Lienhard, J.H. Effect of fouling on performance of pressure retarded osmosis (PRO) and forward osmosis (FO). J. Membr. Sci. 2018, 565, 450–462. [Google Scholar] [CrossRef]
  37. Tang, C.Y.; She, Q.; Lay, W.C.; Wang, R.; Fane, A.G. Coupled effects of internal concentration polarization and fouling on flux behavior of forward osmosis membranes during humic acid filtration. J. Membr. Sci. 2010, 354, 123–133. [Google Scholar] [CrossRef]
  38. Park, M.; Lee, J.J.; Lee, S.; Kim, J.H. Determination of a constant membrane structure parameter in forward osmosis processes. J. Membr. Sci. 2011, 375, 241–248. [Google Scholar] [CrossRef]
  39. Ghanbari, M.; Emadzadeh, D.; Lau, W.; Riazi, H.; Almasi, D.; Ismail, A. Minimizing structural parameter of thin film composite forward osmosis membranes using polysulfone/halloysite nanotubes as membrane substrates. Desalination 2016, 377, 152–162. [Google Scholar] [CrossRef]
  40. Manickam, S.S.; McCutcheon, J.R. Model thin film composite membranes for forward osmosis: Demonstrating the inaccuracy of existing structural parameter models. J. Membr. Sci. 2015, 483, 70–74. [Google Scholar] [CrossRef]
  41. Park, H.B.; Kamcev, J.; Robeson, L.M.; Elimelech, M.; Freeman, B.D. Maximizing the right stuff: The trade-off between membrane permeability and selectivity. Science 2017, 356, eaab0530. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Yip, N.Y.; Tiraferri, A.; Phillip, W.A.; Schiffman, J.D.; Hoover, L.A.; Kim, Y.C.; Elimelech, M. Thin-Film Composite Pressure Retarded Osmosis Membranes for Sustainable Power Generation from Salinity Gradients. Environ. Sci. Technol. 2011, 45, 4360–4369. [Google Scholar] [CrossRef]
  43. Wan, C.F.; Li, B.; Yang, T.; Chung, N.T.-S. Design and fabrication of inner-selective thin-film composite (TFC) hollow fiber modules for pressure retarded osmosis (PRO). Sep. Purif. Technol. 2017, 172, 32–42. [Google Scholar] [CrossRef]
  44. Sivertsen, E.; Holt, T.; Thelin, W.R.; Brekke, G. Modelling mass transport in hollow fibre membranes used for pressure retarded osmosis. J. Membr. Sci. 2012, 69–79. [Google Scholar] [CrossRef]
  45. Cheng, Z.L.; Li, X.; Chung, T.-S. The forward osmosis-pressure retarded osmosis (FO-PRO) hybrid system: A new process to mitigate membrane fouling for sustainable osmotic power generation. J. Membr. Sci. 2018, 559, 63–74. [Google Scholar] [CrossRef]
  46. Nagy, E. Basic Equation of Mass Transport through a Membrane Layer, 2nd ed.; Elsevier: Amsterdam, The Netherlands, 2019. [Google Scholar]
  47. Nagy, E.; Hegedüs, I. Diffusive Plus Convective Mass Transport, Accompanied by Biochemical Reaction, Across Capillary Membrane. Catalysts 2020, 10, 1115. [Google Scholar] [CrossRef]
  48. Alvarez-Fernandez, A.; Reid, B.; Fornerod, M.J.J.J.; Taylor, A.; Divitini, G.; Guldin, S. Structural Characterization of Mesoporous Thin Film Architectures: A Tutorial Overview. ACS Appl. Mater. Interfaces 2020, 12, 5195–5208. [Google Scholar] [CrossRef] [Green Version]
  49. Kim, B.; Gwak, G.; Hong, S. Review on methodology for determining forward osmosis (FO) membrane characteristics: Water permeability (A), solute permeability (B), and structural parameter (S). Desalination 2017, 422, 5–16. [Google Scholar] [CrossRef]
  50. Holt, T.; Sivertsen, E.; Thelin, W.R.; Brekke, G. Pressure dependency of the membrane structure parameter and implication in pressure retarded osmosis (PRO). In Osmotically Driven Membrane Processes. Approach, Development and Current State, 1st ed.; Du, H., Thompson, A., Wang, X., Eds.; InTech: Zagreb, Chroatia, 2018; pp. 111–130. [Google Scholar]
  51. Lee, W.; Kang, P.K.; Kim, A.S.; Lee, S. Impact of surface porosity on water flux and structural parameter in forward osmosis. Desalination 2018, 439, 46–57. [Google Scholar] [CrossRef]
  52. Manickam, S.S.; Gelb, J.; McCutcheon, J.R. Pore structure characterization of asymmetric membranes: Non-destructive characterization of porosity and tortuosity. J. Membr. Sci. 2014, 454, 549–554. [Google Scholar] [CrossRef]
  53. Maisonneuve, J.; Pillay, P.; Laflamme, C.B. Pressure-retarded osmotic power system model considering non-ideal effects. Renew. Energy 2015, 75, 416–424. [Google Scholar] [CrossRef]
  54. Nagy, E. A general, resistance-in-series, salt- and water flux models for forward osmosis and pressure-retarded osmosis for energy generation. J. Membr. Sci. 2014, 460, 71–81. [Google Scholar] [CrossRef]
  55. Nagy, E.; Dudás, J.; Hegedus, I. Improvement of the energy generation by pressure retarded osmosis. Energy 2016, 116, 1323–1333. [Google Scholar] [CrossRef]
  56. Tiraferri, A.; Yip, N.Y.; Phillip, W.A.; Schiffman, J.D.; Elimelech, M. Relating performance of thin-film composite forward osmosis membranes to support layer formation and structure. J. Membr. Sci. 2011, 367, 340–352. [Google Scholar] [CrossRef] [Green Version]
  57. Yip, N.Y.; Elimelech, M. Performance Limiting Effects in Power Generation from Salinity Gradients by Pressure Retarded Osmosis. Environ. Sci. Technol. 2011, 45, 10273–10282. [Google Scholar] [CrossRef] [PubMed]
  58. Loeb, S. Production of energy from concentrated brines by pressure-retarded osmosis. J. Membr. Sci. 1976, 1, 49–63. [Google Scholar] [CrossRef]
  59. Lee, K.; Baker, R.; Lonsdale, H. Membranes for power generation by pressure-retarded osmosis. J. Membr. Sci. 1981, 8, 141–171. [Google Scholar] [CrossRef]
  60. Hofs, B.; Schurer, R.; Harmsen, D.; Ceccarelli, C.; Beerendonk, E.; Cornelissen, E. Characterization and performance of a commercial thin film nanocomposite seawater reverse osmosis membrane and comparison with a thin film composite. J. Membr. Sci. 2013, 446, 68–78. [Google Scholar] [CrossRef]
  61. Abdullah, W.N.A.S.; Tiandee, S.; Lau, W.; Aziz, F.; Ismail, A.F. Potential use of nanofiltration like-forward osmosis membranes for copper ion removal. Chin. J. Chem. Eng. 2020, 28, 420–428. [Google Scholar] [CrossRef]
  62. Motsa, M.M.; Mamba, B.B. Forward Osmosis as a Pre-Treatment Step for Seawater Dilution and Wastewater Reclamation. In Osmotically Driven Membrane Processes. Approach, Developmet and Current States; Du, H., Thompson, A., Wang, X., Eds.; InTech: Zagreb, Croatia, 2018; Volume i, pp. 207–234. [Google Scholar]
  63. She, Q.; Jin, X.; Tang, C.Y. Osmotic power production from salinity gradient resource by pressure retarded osmosis: Effects of operating conditions and reverse solute diffusion. J. Membr. Sci. 2012, 401–402, 262–273. [Google Scholar] [CrossRef]
  64. Roy, Y.; Lienhard, J.H. On the presence of solute-solvent transport coupling in reverse osmosis. J. Membr. Sci. 2020, 611, 118272. [Google Scholar] [CrossRef]
  65. Altaee, A.; Zhou, J.; AlAnezi, A.A.; Zaragoza, G. Pressure retarded osmosis process for power generation: Feasibility, energy balance and controlling parameters. Appl. Energy 2017, 206, 303–311. [Google Scholar] [CrossRef]
  66. Bilal, M.; Asgher, M.; Cheng, H.; Yan, Y.; Iqbal, H.M. Multi-point enzyme immobilization, surface chemistry, and novel platforms: A paradigm shift in biocatalyst design. Crit. Rev. Biotechnol. 2019, 39, 202–219. [Google Scholar] [CrossRef] [PubMed]
  67. Kim, Y.C.; Lee, J.H.; Park, S.-J. Novel crossflow membrane cell with asymmetric channels: Design and pressure-retarded osmosis performance test. J. Membr. Sci. 2015, 476, 76–86. [Google Scholar] [CrossRef]
  68. Wan, C.F.; Chung, N.T.-S. Osmotic power generation by pressure retarded osmosis using seawater brine as the draw solution and wastewater retentate as the feed. J. Membr. Sci. 2015, 479, 148–158. [Google Scholar] [CrossRef]
  69. Gao, Y.; Wang, Y.-N.; Li, W.; Tang, C.Y. Characterization of internal and external concentration polarizations during forward osmosis processes. Desalination 2014, 338, 65–73. [Google Scholar] [CrossRef]
  70. Yip, N.Y.; Tiraferri, A.; Phillip, W.A.; Schiffman, J.D.; Elimelech, M. High Performance Thin-Film Composite Forward Osmosis Membrane. Environ. Sci. Technol. 2010, 44, 3812–3818. [Google Scholar] [CrossRef] [PubMed]
  71. Wang, R.; Shi, L.; Tang, C.Y.; Chou, S.; Qiu, C.; Fane, A.G. Characterization of novel forward osmosis hollow fiber membranes. J. Membr. Sci. 2010, 355, 158–167. [Google Scholar] [CrossRef]
  72. Manickam, S.S.; McCutcheon, J.R. Understanding mass transfer through asymmetric membranes during forward osmosis: A historical perspective and critical review on measuring structural parameter with semi-empirical models and characterization approaches. Desalination 2017, 421, 110–126. [Google Scholar] [CrossRef]
  73. Lee, K.P.; Arnot, T.C.; Mattia, D. A review of reverse osmosis membrane materials for desalination—Development to date and future potential. J. Membr. Sci. 2011, 370, 1–22. [Google Scholar] [CrossRef] [Green Version]
  74. Fujioka, T.; Oshima, N.; Suzuki, R.; Price, W.E.; Nghiem, L.D. Probing the internal structure of reverse osmosis membranes by positron annihilation spectroscopy: Gaining more insight into the transport of water and small solutes. J. Membr. Sci. 2015, 486, 106–118. [Google Scholar] [CrossRef] [Green Version]
  75. Yan, H.; Miao, X.; Xu, J.; Pan, G.; Zhang, Y.; Shi, Y.; Guo, M.; Liu, Y. The porous structure of the fully-aromatic polyamide film in reverse osmosis membranes. J. Membr. Sci. 2015, 475, 504–510. [Google Scholar] [CrossRef]
  76. Vitagliano, V.; Lyons, P.A. Diffusion Coefficients for Aqueous Solutions of Sodium Chloride and Barium Chloride. J. Am. Chem. Soc. 1956, 78, 1549–1552. [Google Scholar] [CrossRef]
  77. Tow, E.W.; Lienhard, J.H. Quantifying osmotic membrane fouling to enable comparisons across diverse processes. J. Membr. Sci. 2016, 511, 92–107. [Google Scholar] [CrossRef] [Green Version]
Figure 1. SEM images of the surface and cross-section structure of two membranes used in experiments: QuantumFlux (NanoH2O) membrane (ac) and PFO (Porifera) membrane (df). (a) Surface structure of the active layer of NanoH2O membrane; (b) surface structure of the support, or non-woven fabric, layer of NanoH2O membrane with observable circular fibers (15–25 μm thick); (c) the cross-section of NanoH2O membrane contains a thin (500–900 nm) active layer and a thick (70–80 μm) porous support layer, which can result in unique separation features; (d) the active surface of PFO membrane; (e) the surface of the support layer of PFO membrane; (f) the cross-section of PFO membrane contains a very thin active layer (f1) and porous sponge-like support layer (f2). PFO membrane fibers of Nomex polymer (50 μm thick) are also well detectable.
Figure 1. SEM images of the surface and cross-section structure of two membranes used in experiments: QuantumFlux (NanoH2O) membrane (ac) and PFO (Porifera) membrane (df). (a) Surface structure of the active layer of NanoH2O membrane; (b) surface structure of the support, or non-woven fabric, layer of NanoH2O membrane with observable circular fibers (15–25 μm thick); (c) the cross-section of NanoH2O membrane contains a thin (500–900 nm) active layer and a thick (70–80 μm) porous support layer, which can result in unique separation features; (d) the active surface of PFO membrane; (e) the surface of the support layer of PFO membrane; (f) the cross-section of PFO membrane contains a very thin active layer (f1) and porous sponge-like support layer (f2). PFO membrane fibers of Nomex polymer (50 μm thick) are also well detectable.
Membranes 11 00128 g001
Figure 2. (a) Osmotic water fluxes of the NanoH2O membrane in pressure-retarded osmosis (PRO) and forward osmosis (FO) modes. Experiments were performed at 22 °C and 0.20–0.25 m/s cross-flow velocity with DI water at the feed side and without a hydraulic pressure difference (1 × 10−6 m/s = 3.6 Lmh). (b) Salt fluxes of the NanoH2O membrane in PRO and FO modes. Experiments were performed at 22 °C and 0.20–0.25 m/s cross-flow velocity with DI water at the feed side and without a hydraulic pressure difference.
Figure 2. (a) Osmotic water fluxes of the NanoH2O membrane in pressure-retarded osmosis (PRO) and forward osmosis (FO) modes. Experiments were performed at 22 °C and 0.20–0.25 m/s cross-flow velocity with DI water at the feed side and without a hydraulic pressure difference (1 × 10−6 m/s = 3.6 Lmh). (b) Salt fluxes of the NanoH2O membrane in PRO and FO modes. Experiments were performed at 22 °C and 0.20–0.25 m/s cross-flow velocity with DI water at the feed side and without a hydraulic pressure difference.
Membranes 11 00128 g002
Figure 3. Osmotic water flux of the Porifera membrane as a function of the draw NaCl concentration in PRO mode. Experiments were performed at 22 °C and at 0.20–0.25 m/s cross-flow velocity with DI water at the feed side and without a hydraulic pressure difference (the continuous line is a fitted curve to the measured points).
Figure 3. Osmotic water flux of the Porifera membrane as a function of the draw NaCl concentration in PRO mode. Experiments were performed at 22 °C and at 0.20–0.25 m/s cross-flow velocity with DI water at the feed side and without a hydraulic pressure difference (the continuous line is a fitted curve to the measured points).
Membranes 11 00128 g003
Figure 4. Predicted power density/water flux as a function of the solute concentration for PRO and FO systems comparing the linear approach and the OLI Stream Analyzer’s data [35] for the osmotic pressure. (Other parameters are listed in Table 2).
Figure 4. Predicted power density/water flux as a function of the solute concentration for PRO and FO systems comparing the linear approach and the OLI Stream Analyzer’s data [35] for the osmotic pressure. (Other parameters are listed in Table 2).
Membranes 11 00128 g004
Figure 5. Change of the relative error in percent ( π l i n π O L I / π x , where πx indicates the higher pressure in each case) as a function of the solute concentration, applying for prediction of the π value by the van ’t Hoff linear approach and the OLI software.
Figure 5. Change of the relative error in percent ( π l i n π O L I / π x , where πx indicates the higher pressure in each case) as a function of the solute concentration, applying for prediction of the π value by the van ’t Hoff linear approach and the OLI software.
Membranes 11 00128 g005
Figure 6. Power density/water flux as a function of the hydraulic pressure difference applying the van ’t Hoff approach and the OLI software for prediction of the osmotic pressure (other parameters are listed in Table 2).
Figure 6. Power density/water flux as a function of the hydraulic pressure difference applying the van ’t Hoff approach and the OLI software for prediction of the osmotic pressure (other parameters are listed in Table 2).
Membranes 11 00128 g006
Figure 7. Water flux as a function of the membrane structural parameter at two different values of the hydraulic pressure (Other parameters are listed in Table 2) (1 Lmh = 0.2778 × 10−6 m/s).
Figure 7. Water flux as a function of the membrane structural parameter at two different values of the hydraulic pressure (Other parameters are listed in Table 2) (1 Lmh = 0.2778 × 10−6 m/s).
Membranes 11 00128 g007
Figure 8. Power density and concentration of the membrane active layer as a function of the water permeability (constant parameters are given in Table 2).
Figure 8. Power density and concentration of the membrane active layer as a function of the water permeability (constant parameters are given in Table 2).
Membranes 11 00128 g008
Figure 9. Power density and concentration of the membrane active layer as a function of the solute permeability (Remained parameter are given in Table 2).
Figure 9. Power density and concentration of the membrane active layer as a function of the solute permeability (Remained parameter are given in Table 2).
Membranes 11 00128 g009
Figure 10. Change of the S/S values as a function of the water flux in a PRO system, taking into account their predicted error of the measured data (±20%), given by Manickam and McCutcheon [40]. Parameter is the external mass transfer coefficient in the boundary layer facing the active membrane layer, at Cd = 1.0 M, and kf∞ (A = 1.71 × 10−7 m/s·bar; B = 1.94 × 10−7 m/s; Cf = 0; S = 803 × 10−6 m; ΔP = 0; D = 1.5 × 10−9 m2/s). On the upper horizontal axis, the calculated value of the corresponding structural parameter is given for kd→∞, kf→∞.
Figure 10. Change of the S/S values as a function of the water flux in a PRO system, taking into account their predicted error of the measured data (±20%), given by Manickam and McCutcheon [40]. Parameter is the external mass transfer coefficient in the boundary layer facing the active membrane layer, at Cd = 1.0 M, and kf∞ (A = 1.71 × 10−7 m/s·bar; B = 1.94 × 10−7 m/s; Cf = 0; S = 803 × 10−6 m; ΔP = 0; D = 1.5 × 10−9 m2/s). On the upper horizontal axis, the calculated value of the corresponding structural parameter is given for kd→∞, kf→∞.
Membranes 11 00128 g010
Figure 11. Change of the S/S values as a function of the water flux in the FO system, taking into account ±20% error of the measured data, at Cd = 1.5 M, kf→∞ applied data of Manickam and McCutcheon [40] (other parameters are given in caption of Figure 10). On the upper horizontal axis, the calculated value of the corresponding structural parameter, S, is given for kd∞, kf→∞.
Figure 11. Change of the S/S values as a function of the water flux in the FO system, taking into account ±20% error of the measured data, at Cd = 1.5 M, kf→∞ applied data of Manickam and McCutcheon [40] (other parameters are given in caption of Figure 10). On the upper horizontal axis, the calculated value of the corresponding structural parameter, S, is given for kd∞, kf→∞.
Membranes 11 00128 g011
Figure 12. Value of S/S as a function of the draw side mass transfer coefficient with both operating modes, i.e., PRO and FO; Cd = 1.0 M; kf∞; for PRO: Jw, = 4.92 × 10−6 m/s; S = 803 μm; for FO: J = 3.0 × 10−6 m/s; S = 395 μm; other parameters as they are given in caption of Figure 10.
Figure 12. Value of S/S as a function of the draw side mass transfer coefficient with both operating modes, i.e., PRO and FO; Cd = 1.0 M; kf∞; for PRO: Jw, = 4.92 × 10−6 m/s; S = 803 μm; for FO: J = 3.0 × 10−6 m/s; S = 395 μm; other parameters as they are given in caption of Figure 10.
Membranes 11 00128 g012
Figure 13. Effect of the water permeability on the value of S/S and on the water flux, Jw,∞, as a function of draw side mass transfer coefficient, in a PRO system; Cd = 1.0 M; S = 803 μm; πd = 43.3 bar; kf∞; other parameters as they are given in caption of Table 5.
Figure 13. Effect of the water permeability on the value of S/S and on the water flux, Jw,∞, as a function of draw side mass transfer coefficient, in a PRO system; Cd = 1.0 M; S = 803 μm; πd = 43.3 bar; kf∞; other parameters as they are given in caption of Table 5.
Membranes 11 00128 g013
Figure 14. Effect of the solute permeability on the value of S/S and on the water flux, Jw,∞, as a function of draw side mass transfer coefficient, in a PRO system; Cd = 1.0 M; S = 803 μm; kf→∞; other parameters as in Table 5.
Figure 14. Effect of the solute permeability on the value of S/S and on the water flux, Jw,∞, as a function of draw side mass transfer coefficient, in a PRO system; Cd = 1.0 M; S = 803 μm; kf→∞; other parameters as in Table 5.
Membranes 11 00128 g014
Table 1. Selective layer characteristics of the Porifera and NanoH2O membranes determined by cross-flow RO experiments at 22 °C and 0.20–0.25 m/s cross-flow velocity. deionized water (DI water) water was fed at 3, 5, 7, 9 bar to measure the A values, and 2 g/L NaCl solution was used to measure the B and R values at 3, 5, 7 bar.
Table 1. Selective layer characteristics of the Porifera and NanoH2O membranes determined by cross-flow RO experiments at 22 °C and 0.20–0.25 m/s cross-flow velocity. deionized water (DI water) water was fed at 3, 5, 7, 9 bar to measure the A values, and 2 g/L NaCl solution was used to measure the B and R values at 3, 5, 7 bar.
MembraneWater Permeability, A,
10−7 m/s·bar
Salt Permeability, B
10−7 m/s
NaCl Rejection Coefficient, %R
Porifera5.0 ± 0.151.5 ± 0.690 ± 4.6
NanoH2O1.0 ± 0.030.5± 0.1290.2 ± 2.5
Table 2. General operating and transport parameters.
Table 2. General operating and transport parameters.
ParametersValues
Feed solute concentration0.015 M
Solute concentration in the draw solution0.6 M, or varies
External mass transfer coefficients, draw side3.85 × 10−5 m/s
External mass transfer coefficient, feed side3.85 × 10−5 m/s
Diffusion coefficient at high salinity1.367 × 10−9 m2/s
Diffusion coefficient at low salinity1.294 × 10−9 m2/s
Hydraulic pressure difference10 bar (or 0 and 15 bar)
Membrane transport parametersA = 1.9 × 10−7 m/s·bar or varies
B = 5.02 × 10−7 m/s or varies
S = 5 × 10−4 m or varies
Table 3. Predicted values of the structural parameter for NanoH2O membrane as a function of draw salt concentration. DI is used for the feed phase; the osmotic pressure was predicted by the van ’t Hoff approach and by the OLI software for calculation of the structural parameter, S, with (kd = kf = 2 × 10−5 m/s) and without external mass transfer resistances (kd = kf∞).
Table 3. Predicted values of the structural parameter for NanoH2O membrane as a function of draw salt concentration. DI is used for the feed phase; the osmotic pressure was predicted by the van ’t Hoff approach and by the OLI software for calculation of the structural parameter, S, with (kd = kf = 2 × 10−5 m/s) and without external mass transfer resistances (kd = kf∞).
PROFO
Cd, mol/L,
PRO, FO
Jw,
10−6 m/s
S, 10−6 m
NanoH2O
Jw,
10−6 m/s
S, 10−6 m
NanoH2O
van ’t HoffOLI * van ’t HoffOLI *
1 *2 *1 *2 * 1 *2 *1 *2 *
0.270.2815,95015,78814,54014,3640.268438835072677178
0.620.3715,39615,23514,61014,4300.329898981290868999
1.030.4016,68716,53116,27016,0490.477262717968876804
1.530.4715,51015,35515,30615,2220.527531744874187335
Average 15,88615,72718,18215,016 8282819776657579
1 *: kd = kf∞; 2 *: kd = kf = 2 × 10−5 m/s; OLI *: π = 5.94028C2 + 37.4521C.
Table 4. Predicted values of the structural parameter for the Porifera membrane as a function of draw salt concentration. DI is used for feed phase; the osmotic pressure was predicted by the van ’t Hoff approach and by the OLI software, using the curve-fit expression (π = 5.94028C2 + 37.4521C for C > 0.6) and also the linear interpolation of the OLI data given in Ref. [35], for calculation of the structural parameter, S, with (kd = kf = 2 × 10−5 m/s) and without external mass transfer resistances (kd = kf→∞).
Table 4. Predicted values of the structural parameter for the Porifera membrane as a function of draw salt concentration. DI is used for feed phase; the osmotic pressure was predicted by the van ’t Hoff approach and by the OLI software, using the curve-fit expression (π = 5.94028C2 + 37.4521C for C > 0.6) and also the linear interpolation of the OLI data given in Ref. [35], for calculation of the structural parameter, S, with (kd = kf = 2 × 10−5 m/s) and without external mass transfer resistances (kd = kf→∞).
Draw Concentration, MS (van ’t Hoff), 10−6 mS (OLI, Interpolated from Table), 10−6 mS (OLI, Curve Fit), 10−6 m;
1 *2 *1 *2 *1 *2 *
0.26205718581976177018151591
0.4813708931118132712571034
0.501114116910788471019766
0.60120099911699621108885
0.82125010581228104111911003
0.88104985610298321002799
1.28112695011639891102923
1.33105587710568791150973
Average:1169991113510241138936
1 * kd = kf→∞; 2 * kd = kf = 2 × 10−5 m/s.
Table 5. Application of measured data of Manickam and McCutcheon [40], for illustration of the effect of the external transfer coefficient on the predicted values of the structural parameter. The first two columns contain data published in [40]. The measured water flux data is shown to one decimal point in order to reflect that the reading error can reach 10% [40]. The third column (kd = kf→∞) is data published in [40]. The fourth and fifth columns are calculations conducted in this study. (A = 0.606 L/m2h = 1.71 × 10−7 m/s·bar; B = 0.699 L/m2h = 1.94 × 10−7 m/s; D = 1.5 × 10−9 m2/s; Cf = 0; ΔP = 0).
Table 5. Application of measured data of Manickam and McCutcheon [40], for illustration of the effect of the external transfer coefficient on the predicted values of the structural parameter. The first two columns contain data published in [40]. The measured water flux data is shown to one decimal point in order to reflect that the reading error can reach 10% [40]. The third column (kd = kf→∞) is data published in [40]. The fourth and fifth columns are calculations conducted in this study. (A = 0.606 L/m2h = 1.71 × 10−7 m/s·bar; B = 0.699 L/m2h = 1.94 × 10−7 m/s; D = 1.5 × 10−9 m2/s; Cf = 0; ΔP = 0).
Cd
mol/L
Jw,
Measured
10−6 m/s
S, (OLI Software)
10−6 m
S, (van ’t Hoff), 10−6 mError, %
kd = kfkd = 5 × 10−5 m/s, kfkd = 2 × 10−5 m/s, kfkd = kf
P0.52.8791649(S/S = (0.82)356 (0.45)111729.1
R1.04.9803708 (0.88)514 (0.64)89710.5
O1.56.4789726 (0.92)584 (0.74)8163.3
F0.52.1327298 (0.91)252 (0.77)45928.8
O1.03.0434404 (0.93)360 (0.83)49111.6
1.54.0395360 (0.92)320 (0.81)4154.8
Table 6. Application of measured data of Tang et al. [37] to recalculate the effect of the measured water flux on the prediction of structural parameter, using OLI software instead of the van ’t Hoff linear approach, and to illustrate the effect of the external transfer coefficient on the predicted values of S. The first two columns (water flux, predicted structural parameter) are data published in [37], and the other data are calculated by authors of this paper. AS indicates that it is less than zero, i.e., the measured water flux is higher than expected for the tested conditions. (A = 2.22 × 10−7 m/s·bar; B = 1.7 × 10−7 m/s; D = 1.5 × 10−9 m2/s; Cf = 0.01 M; ΔP = 0).
Table 6. Application of measured data of Tang et al. [37] to recalculate the effect of the measured water flux on the prediction of structural parameter, using OLI software instead of the van ’t Hoff linear approach, and to illustrate the effect of the external transfer coefficient on the predicted values of S. The first two columns (water flux, predicted structural parameter) are data published in [37], and the other data are calculated by authors of this paper. AS indicates that it is less than zero, i.e., the measured water flux is higher than expected for the tested conditions. (A = 2.22 × 10−7 m/s·bar; B = 1.7 × 10−7 m/s; D = 1.5 × 10−9 m2/s; Cf = 0.01 M; ΔP = 0).
Cd
mol/L
Jw,
Measured
10−6 m/s
S, (van ‘t Hoff)
10−6 m
kd = kf
S, (OLI Software), 10−6 m,
kd = kfkd = 5 × 10−5 m/s
kf→∞
kd = 2 × 10−5 m/s
kf→∞
0.54.2581210--
P1.07.3546474320 (0.675)-
R2.011.6471481407 (0.846)142 (0.295)
O3.013.6431499447 (0.896)332 (0.665)
4.015.3460494449 (0.91)361 (0.730)
0.251.6467227197 (0.87)152 (0.67)
F0.753.0460381352 (0.92)309 (0.81)
O1.54.2460441412 (0.93)369 (0.84)
2.04.8456461431 (0.93)385 (0.84)
3.05.6456490458 (0.93)410 (0.84)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nagy, E.; Hegedüs, I.; Rehman, D.; Wei, Q.J.; Ahdab, Y.D.; Lienhard, J.H. The Need for Accurate Osmotic Pressure and Mass Transfer Resistances in Modeling Osmotically Driven Membrane Processes. Membranes 2021, 11, 128. https://doi.org/10.3390/membranes11020128

AMA Style

Nagy E, Hegedüs I, Rehman D, Wei QJ, Ahdab YD, Lienhard JH. The Need for Accurate Osmotic Pressure and Mass Transfer Resistances in Modeling Osmotically Driven Membrane Processes. Membranes. 2021; 11(2):128. https://doi.org/10.3390/membranes11020128

Chicago/Turabian Style

Nagy, Endre, Imre Hegedüs, Danyal Rehman, Quantum J. Wei, Yvana D. Ahdab, and John H. Lienhard. 2021. "The Need for Accurate Osmotic Pressure and Mass Transfer Resistances in Modeling Osmotically Driven Membrane Processes" Membranes 11, no. 2: 128. https://doi.org/10.3390/membranes11020128

APA Style

Nagy, E., Hegedüs, I., Rehman, D., Wei, Q. J., Ahdab, Y. D., & Lienhard, J. H. (2021). The Need for Accurate Osmotic Pressure and Mass Transfer Resistances in Modeling Osmotically Driven Membrane Processes. Membranes, 11(2), 128. https://doi.org/10.3390/membranes11020128

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

Article Metrics

Back to TopTop