A Unified Model-Based Analysis and Optimization ... - ACS Publications

Nov 18, 2013 - Process intensification of seawater reverse osmosis through enhanced train capacity and module size – Simulation on Lanzarote IV SWRO...
0 downloads 0 Views 921KB Size
Article pubs.acs.org/IECR

A Unified Model-Based Analysis and Optimization of Specific Energy Consumption in BWRO and SWRO Mingheng Li* Department of Chemical and Materials Engineering, California State Polytechnic University, Pomona, California 91768, United States ABSTRACT: A unified mathematical model is developed to describe both brackish water (BW) and seawater (SW) reverse osmosis (RO) desalination processes. The model contains two dimensionless parameters: κ, quantifying the effect of retentate pressure drop, and γ, reflecting the membrane capacity demand ratio. A γ vs. κ map derived from data obtained from RO desalination plants worldwide is used to show that BWRO and SWRO are in completely different regions. The minimization of specific energy consumption (SEC) is formulated and solved as a nonlinear optimization problem using representative values of γ and κ in both BWRO and SWRO. Several computational case studies are carried out to demonstrate the differences in operating BWRO and SWRO. Discussions are made on the effect of fractional recovery, RO configurations (single-stage, two-stage, and two-stage with booster pump) as well as brine recirculation on SEC.





INTRODUCTION Specific energy consumption (SEC), or pump energy consumption divided by water production rate, is a very important topic in RO desalination.1,2 There have been significant research and development efforts in model-based analysis and optimization of RO design and operation.3−9 Model-based control has also been developed recently to reduce SEC in a pilot-scale RO system10,11 In two previous papers, 12,13 the author provided a comprehensive analysis of single- and multistage SWRO with/without energy recovery device (ERD) from firstprinciples, based on the assumption of negligible retentate pressure drop. Several dimensionless parameters (γ, α, and Y) were introduced to show the coupled behavior between membrane property (area and permeability), feed conditions (flow rate and salinity), and operating conditions (hydraulic pressure difference and fractional recovery). Later on, the retentate pressure drop effect was explicitly accounted for in the optimization model for both single-train and multitrain BWRO processes.14,15 A validation of model-based optimization results was carried out in a BWRO plant in Chino, California, and a 10% decrease in SEC and a substantial reduction in brine disposal cost were demonstrated.16 As a continuation of the author’s previous work, this paper aims to utilize a unified dimensionless model to explain the similarities and differences between SWRO and BWRO. The model is based on several dimensionless parameters, including κ, quantifying the effect of retentate pressure drop, and γ, a previously derived parameter reflecting the membrane capacity demand ratio.12 Data from both BWRO and SWRO plants worldwide are used to show the typical ranges of γ and κ. Using representative values of γ and κ in both BWRO and SWRO, a mathematical model is formulated to reduce SEC by optimizing a dimensionless pressure parameter α. The optimal SECs in both BWRO and SWRO using three different RO configurations (single-stage, two-stage, and two-stage with booster pump) are compared. The effect of brine recirculation on SEC is also discussed. © 2013 American Chemical Society

DIMENSIONLESS RO MODEL On the basis of assumptions of Darcy’s law for mass transfer and quadratic dependence of retentate pressure drop on velocity, the following mathematical model was derived to describe a single-stage RO process14 −

dQ (x) dx

⎛ ⎞ Q = ALp⎜ΔP − 0 Δπ0⎟ Q ⎝ ⎠

d(ΔP(x)) =−kQ 2 dx Q (x )

=Q 0@x = 0

ΔP(x)

= ΔP0@x = 0

(1)

where Q is the retentate flow rate, −dQ is the flow rate of water across membrane of area dA (dA = Adx; x is the distance from entrance normalized by the pressure vessel length, and A is the total area of all RO elements in this stage), Lp is the membrane hydraulic permeability, ΔP and Δπ are the differences in hydraulic pressures and osmotic pressures across the membrane, respectively. Δπ0, Q0, and ΔP0 are Δπ, Q, and ΔP at the stage entrance. k is a coefficient describing the pressure drop along the pressure vessel in the retentate stream. Note that the pressure differential in the permeate stream is typically small. More detailed assumptions and justifications of eq 1 are discussed in the author’s previous work.14 This coupled differential equation captures the profile of driving force and water flux, and is more advanced than standard membrane projection model in literature, which uses averaged flux and driving force. Defining α = Δπ0/ΔP0, γ = ALpΔπ0/Q0, κ = kQ20/Δπ0, q = Q/Q0, p = ΔP/ΔP0, the following equation is derived: Received: Revised: Accepted: Published: 17241

September 18, 2013 November 8, 2013 November 15, 2013 November 18, 2013 dx.doi.org/10.1021/ie4031083 | Ind. Eng. Chem. Res. 2013, 52, 17241−17248

Industrial & Engineering Chemistry Research

Article

relationships between p and q at the stage level and those at the system level are shown as follows

dp(x) = −καq2(x) dx ⎛ p(x ) dq(x) 1 ⎞ = −γ ⎜ − ⎟ dx q(x) ⎠ ⎝ α p(x )

= 1, @x = 0

q(x)

= 1, @x = 0

i−1

p(x) = pi (x) ∏ pk (k), i − 1 < x ≤ i , i = 1, 2, ..., N k=0 i−1

q(x) = qi(x) ∏ qk (k), i − 1 < x ≤ i , i = 1, 2, ..., N

(2)

k=0

(5)

An integration of eq 2 would provide the profiles of p(x) and q(x) along the RO stage. The fractional recovery of water is then calculated as Y = 1 − q(1), where q(1) is q evaluated at x = 1, or the end of this stage. It is seen from eq 2 that the RO system is characterized by three different dimensionless parameters α, γ, and κ, all described in the form of ratios. As discussed in previous publications,12,13 α is the ratio of osmotic pressure difference to hydraulic pressure difference across the membrane at the entrance of the stage, and γ is the ratio of flow rate across the membrane assuming a constant drive force of Δπ0 to the retentate flow rate at the entrance of the stage. γ may be considered as a measure of the membrane capacity demand ratio. κ introduced in this work is the ratio of maximum retentate pressure drop (i.e., the theoretical pressure drop when there is no water recovery, or the retentate flow rate is always the same as the intake flow rate) to the osmotic pressure difference at the entrance of the stage. It is a measure of driving force loss that adversely affects the water production. According to fluid mechanics of channel flow, κ can be described as follows κ = kQ 02/Δπ0 = 2f (Lpv /DH)ρ(Q 0/Ac)2 /Δπ0

where p0(0) ≡ q0(0) ≡ 1. The system recovery at the end of all stages Y = 1 − Πi N= 1qi(i). The parameters αi, γi, and κi should be updated as eq 4 is solved sequentially from stage to stage. For example, in a typical two-stage RO process with a 2:1 array to achieve a high system system recovery and to maintain a relatively uniform feed flow in each stage,14,17 A2 = A1/2, Ac2 = Ac1/2, k2 = 22k1, Δπ1 = Δπ0/ q1(1), ΔP1 = ΔP0/p1(1), Q1 = Q0q1(1). As a result, γ2 = γ1/ 2q21(1), α2 = α1/q1(1)/p1(1), κ2 = 4q31(1)κ1. Once the first stage is solved, α2, γ2, and κ2 are determined so that the second stage can be solved in a similar way. The system recovery at the end of the second stage is calculated as Y = 1−q1(1)q2(2). The SEC normalized by the feed osmotic pressure, or NSEC, is a dimensionless number12

dx dqi(x) dx

(3)

NSEC=

min J α1

qi(x)

= 1, @x = i − 1, i = 1, 2, ..., N

1/α − (Pbp − Pp)/Δπ0 N

ηpump(1 − ∏i = 1 qi(i))

(7)

=

1/α1 − (Pbp − Pp)/Δπ0 N

ηpump(1 − ∏i = 1 qi(i))

s. t . dpi (x) dx dqi(x) dx

⎛ p (x ) 1 ⎞⎟ , i = 1, 2, ..., N = −γi⎜⎜ i − qi(x) ⎟⎠ ⎝ αi = 1, @x = i − 1, i = 1, 2, ..., N

(6)

The minimization of NSEC in a single- or multistage RO is then formulated as an optimization problem as follows

= −κiαiqi2(x), i = 1, 2, ..., N

pi (x)

ηpump(Q 0Y )Δπ0

where ΔPpump is the pressure increase across the pump, ηpump is the pump efficiency. Note that ΔP0 = Pbp + ΔPpump − Pp, where Pbp is the feed pressure before pump, and Pp is the permeate pressure. Therefore

where f is the Fanning friction coefficient, DH is the hydraulic diameter between membrane sheets, Ac is the total crosssectional area in the RO stage for retentate flow, and Lpv is the total length of the pressure vessel that houses the membrane elements in series. Apparently, κ will reduce if more pressure vessels are put in parallel in a RO stage, or fewer membrane elements are enclosed in a pressure vessel. The introduction of dimensionless parameters simplifies the analysis of RO process as the results are scalable. Because BWRO and SWRO differ in Δπ0 (around 10−20 psi in BWRO and 390 psi in SWRO), they may exhibit different process characteristics and require different design and operation strategies. Equation 2 can be easily extended to describe RO networks with complicated configurations. For example, a RO system with N stages in series may be modeled as follows dpi (x)

Q 0ΔPpump

NSEC =

= −κiαiqi2(x), i = 1, 2, ..., N ⎛ p (x ) 1 ⎞⎟ , i = 1, 2, ..., N = −γi⎜⎜ i − qi(x) ⎟⎠ ⎝ αi

pi (x)

= 1, @x = i − 1, i = 1, 2, ..., N

qi(x)

= 1, @x = i − 1, i = 1, 2, ..., N

(8)

The optimization of a RO network with more complicated configurations can be formulated using methods developed in the author’s previous work.13 For a two-stage RO with interstage booster pump, the optimization problem (neglecting pump efficiency and assuming pbp = pp for simplicity) may be formulated as follows

(4)

where pi and qi are defined based on the ratio of ΔP(x) and Q(x) in RO stage i to the ones at the entrance of stage i. The 17242

dx.doi.org/10.1021/ie4031083 | Ind. Eng. Chem. Res. 2013, 52, 17241−17248

Industrial & Engineering Chemistry Research min J

=

Article

in parentheses. In all calculations, osmotic pressure is estimated based on TDS data using an online calculator (http://www. lenntech.com/calculators/osmotic/osmotic-pressure.htm), if it is not provided in the literature. A significant difference between BWRO and SWRO is observed. For BWRO, γ1 is typically about 0.1 (even γtotal is less than 0.2 if multiple stages are used), whereas κ1 is between 1 and 6. For SWRO, γ1 is between 0.5 and 1.5, whereas κ1 is less than 0.2. The differences are presented in the γ vs. κ map in Figure 1. These differences are mainly caused by the fact that

1/α1 + 1/α2 − p1 (1)q1(1) 1 − q1(1)q2(2)

α1, α2

s. t . dpi (x) dx dqi(x) dx

= −κiαiqi2(x), i = 1, 2 ⎛ p (x ) 1 ⎞⎟ , i = 1, 2 = −γi⎜⎜ i − qi(x) ⎟⎠ ⎝ αi

pi (x)

= 1, @x = i − 1, i = 1, 2

qi(x)

= 1, @x = i − 1, i = 1, 2

γ1

= (2/3)γtotal

γ2

= γ1/2/q12(1)

κ2

= 4κ1q13(1)

α2

≤ α1/p1 (1)/q1(1)

−αi

≤ 0, i = 1, 2

αi

≤ 1, i = 1, 2

(9)

In eq 9, γtotal is defined based on the total area of all RO elements in both stages (γtotal = AtotalLpΔπ0/Q0). Note that while Atotal = A1 + A2, γtotal ≠ γ1 + γ2. Inequality constraint α2 ≤ α1/p1(1)/q1(1) is included to ensure non-negative pressure increase across the interstage pump. The optimization models described by eqs 8 and 9 can be solved using common constrained multivariable nonlinear optimization packages, e.g., fmincon in Matlab.

Figure 1. γ and κ in BWRO and SWRO.

the osmotic pressure of seawater is at least 1 order of magnitude higher than the one of brackish water. Current BWRO plants are not built with a large γ because of constraints in capital cost. For example, the WMZM plant in Jordan has three BWRO stages with a total of 791 elements, but γtotal is estimated to be 0.12 base on the operational data in one of the skids.27 The differences in γ and κ imply different process characteristics as well as design and operation strategies in SWRO and BRWO. It is worth noting that in a two-pass seawater treatment plant, the first pass is SWRO, whereas the second-pass is BWRO.30−32 Optimization of NSEC in Single-Stage BWRO and SWRO. The optimization problem described by eq 8 is solved for single-stage RO with representative values of γ and κ in BWRO and SWRO. For demonstration purposes, it is assumed that Pp = Pbp. Moreover, the effect of pump efficiency is neglected. Interested readers are referred to the author’s previous work14 for optimization accounting for pump characteristic maps and other process details. The optimal NSEC and the corresponding α and Y are shown in Figures 2, 3, and 4. It is seen that the optimal NSEC decreases as κ decreases and γ increases. When both κ and 1/γ go to zero, NSEC reaches 4 and Y = 50% by applying ΔP0 = 2Δπ0. This is well-known in the RO literature.12 However, in real-world RO plants, γ is always finite (and may be very small in BWRO) because of considerations in capital cost. Even though κ is small in SWRO, it is substantial in BWRO. A larger γ allows the system to be operated closer to the thermodynamic limit, leading to a lower NSEC.12,13 A larger κ, however, reduces the driving force for permeate flow, leading to a higher NSEC. When γ decreases (less membrane capacity) or κ increases (more pressure drop), a higher recovery and a greater applied pressure are suggested in order to achieve the NSEC at



RESULTS AND DISCUSSION Analysis of BWRO and SWRO Plant Data. Data of BWRO and SWRO plants worldwide14,18−29 are analyzed and the parameters α, γ, and κ are either calculated directly based on their definitions or correlated to match plant measurements including pressure differentials and recoveries. The results are summarized in Table 1. γ, κ, and α reported here are based on the first stage for plants employing multistage operation as they change from stage to stage. For these plants, γtotal are provided Table 1. Correlation of Plant Data with Mathematical Model plant

RO type 18

Carlsbad Pilot, USA Tuas, Singapore19 Tuas, Singapore20 Jeddah, Saudi Arabia21,22 North Obhor, Saudi Arabia23 Perth, Australia24 ADC Pilot, USA25 WMZM, Jordan (threestage)26,27 Chino, USA (twostage)14 James City (two-stage), USA28 Ocean Reef Club, USA29

γ1 (γtotal)

κ1

α1

Y

SWRO SWRO SWRO SWRO

1.43 0.66 0.86 0.86

0.13 0.13 0.12 0.08

0.49 0.46 0.50 0.53

0.50 0.45 0.45 0.41

SWRO

0.88

0.06

0.57

0.36

SWRO SWRO BWRO

0.61 0.56 0.07 (0.12)

0.05 0.04 4.15

0.46 0.45 0.09

0.45 0.45 0.88

BWRO

0.04 (0.05)

5.44

0.05

0.81

BWRO

0.10 (0.15)

1.24

0.13

0.78

BWRO

0.08

1.06

0.10

0.65

17243

dx.doi.org/10.1021/ie4031083 | Ind. Eng. Chem. Res. 2013, 52, 17241−17248

Industrial & Engineering Chemistry Research

Article

its lowest attainable level. It is worth noting that the maximum recovery is limited by the ability to produce potable quality permeate in SWRO33 and the occurrence of silica precipitation in BWRO.34 From Figures 2, 3, and 4, one can read optimal operating conditions for typical BWRO and SWRO processes using representative values of γ and κ. The results are shown in Table 2. Apparently, NSEC is only slightly above 4 in SWRO but much larger in BWRO. To achieve the optimal NSECs, ΔP0 is 2.3 and 17 times Δπ0 in SWRO and in BWRO, leading to a significant difference in water recoveries. The optimal recovery in BWRO is much larger. The dimensionless driving forces at the entrance and exit of the stage (p/α − 1/q) are provided in the same table. It appears that SWRO is operated near the thermodynamic limit, but BWRO is far away from it. Comparison of SEC between Single-Stage and TwoStage Designs. Both single-stage and two-stage designs have been used in RO operation, and there are still debates on which one is more efficient from a viewpoint of energy consumption.35 In this subsection, several case studies will be used to clarify the effect of RO configurations on SEC in both BWRO and SWRO. For BWRO, the effect of retentate pressure drop on water recovery is first quantified before the detailed comparisons. Equation 2 was applied to an industrial two-stage BWRO train in Chino, California.14 This BWRO train has 28 pressure vessels in parallel in the first stage and 14 in the second stage. Each pressure vessel houses 7 Dow FILMTEC BW30−400 RO elements in series. The BWRO train does not have an interstage booster pump or ERD. It was derived from plant data that k1 = 2.12 × 10−5 psi/gpm2 and Lp is 0.11 gfd/psi.14 The model predictions using these two parameters were able to match plant data collected in a wide range of operating conditions.16 On the basis of these parameters, it is calculated that κ1 = 5.5 and γ1 = 0.035. A comparison of simulation results with/without retentate pressure drop is shown in Figure 5. When the pressure drop effect is completely ignored (or assuming κ1 = κ2 = 0), the calculated recovery is 91%, leading to an over prediction of about 12% compared to its measured value of 81%. This case study indicates that the retentate pressure drop should not be ignored in BWRO where κ is large. With the above information, one would expect that if all 42 pressure vessels are laid out in parallel in just one stage, this BWRO plant could be more energy-efficient because of reduced retentate pressure drop. The detailed comparison is shown in Figures 6 and 7. The overall retentate pressure drop ratio is only 4% in the single-stage design as compared to 24% in the two-stage design. There are two factors that contribute to such a significant difference. First, the channel area for retentate flow in the single-stage configuration is 1.5 and 3 times of the one in the first and second stage in the two-stage design, respectively. The increased channel area leads to a reduced retentate velocity along the pressure vessel. Second, the total pressure vessel length in the single-stage design is only a half of

Figure 2. Optimal NSEC as a function of κ and γ in single-stage RO without ERD.

Figure 3. Optimal applied pressure as a function of κ and γ in singlestage RO without ERD.

Figure 4. Optimal recovery as a function of κ and γ in single-stage RO without ERD.

Table 2. Comparison between Single-Stage BWRO and SWRO under typical Operation Conditions RO type

γ

κ

αopt

NESCopt

Yopt

p(0)/α−1/q(0)

p(1)/α−1/q(1)

BWRO SWRO Thermodynamic Limit12

0.05 1 ∞

1 0.1 0

0.059 0.438 0.5

22.9 4.2 4

0.75 0.54 0.5

16.1 1.28 1

12.7 0.05 0

17244

dx.doi.org/10.1021/ie4031083 | Ind. Eng. Chem. Res. 2013, 52, 17241−17248

Industrial & Engineering Chemistry Research

Article

calculation shows that a ΔP0 of 163 psi (10% lower than 180 psi in the current two-stage design) could keep the 81% recovery with the same Q0. Therefore, the single-stage configuration is more effective than its two-stage counterpart in BWRO from a viewpoint of SEC (about 10% lower in this case). However, it should be noted that the two-stage design has a less drastic change in retentate flow per stage and less tendency for concentration polarization36 because of a higher velocity. The same conclusion may be drawn for SWRO, even though the difference between single-stage and two-stage configurations is small, given the fact that κ is close to zero in SWRO. However, it is worth pointing out that when an interstage booster pump is used in conjunction with the two-stage design in SWRO, it would be advantageous over the single-stage design, according to the author’s published theoretical analysis.12,13 The computational analysis is based on the assumption κ = 0. The works also show that when γtotal < 0.2, the optimal solution is corresponding to a zero pressure increase in the interstage booster pump, or the two-stage configuration is essentially one-stage. There are industrial cases to substantialize reduced SEC in SWRO brine conversion using interstage booster pumps.37 For example, in the SWRO plant in Las Palmas III, Spain, an increase of 8% in the recovery of and an improvement of 12.65% in the specific energy consumption were observed after interstage booster pumps were installed.38 To theoretically validate whether a booster pump would always reduce SEC for a nonzero κ in SWRO and BWRO, the optimization problems described by eqs 8 and 9 are solved for three different configurations (single-stage, two-stage, and twostage with booster pump, see Figure 8) using representative

Figure 5. Effect of retentate pressure drop on water recovery in an industrial two-stage BWRO train without interstage booster pump.

Figure 6. Comparison of dimensionless applied pressure and retentate flow between single-stage and two-stage designs in BWRO.

Figure 8. Schematic of three different RO configurations used in industry.

values of γ and κ. In order to have a fair comparison, all three configurations have exactly the same membrane properties and flow characteristics. These mean that γtotal = AtotalLpΔπ/Q0 is the same in all configurations. However, it should be noted that if the pressure loss parameter is κ in the first stage of a twostage RO, it becomes (2/3)2κ in a single-stage RO because of increased cross-sectional channel area for retentate flow (based on eq 3). The optimized NSECs and corresponding recoveries for BWRO and SWRO are shown in Figures 9 and 10. For BWRO, the optimal NSECs with/without a booster pump are almost the same. Both two-stage designs have a higher NSEC than the one in single-stage design at their optimal conditions. Therefore, installing a booster pump would not help reduce SEC in BWRO. In SWRO, the two-stage with booster pump

Figure 7. Comparison of dimensionless driving force between singlestage and two-stage designs in BWRO.

the one in the two-stage design. A less retentate pressure drop results in a higher driving force for water production (see Figure 7). It is shown that the predicted water recovery in the single-stage design is 89% with the same Q0 and ΔPpump. To maintain the same production, Q0 and/or ΔPpump may be lowered in the single-stage configuration. For example, the 17245

dx.doi.org/10.1021/ie4031083 | Ind. Eng. Chem. Res. 2013, 52, 17241−17248

Industrial & Engineering Chemistry Research

Article

Figure 10. Optimal (a) NSECs and (b) recoveries in SWRO for single-stage, two-stage, and two-stage with booster pump using κ = 0.1.

Figure 9. Optimal (a) NSECs and (b) recoveries in BWRO for singlestage, two-stage, and two-stage with booster pump using κ = 1.

Figure 11 shows an example in which an ERD (PX by Energy Recovery Inc.) is intentionally unbalanced in a BWRO system

design outperforms the other two designs at a large γtotal (e.g., 0.8 or greater in this case). However, the single-stage may be the most energy efficient at a γtotal smaller than 0.8. Calculation also shows that this critical γtotal reduces from 0.8 to 0.7 when κ decreases from 0.1 to 0.05. In both BWRO and SWRO studied here, the two-stage RO without booster pump is always no better than the other two configurations from a viewpoint of SEC. The above comparisons among three RO configurations are made at their best operating conditions. It is possible that one design at its optimal condition is more energy-effective than another at a nonoptimal condition. Moreover, because the twostage designs (with/without booster pumps) have higher recoveries at its optimal condition than the one in the singlestage, adding a constraint of maximum recovery may also change the comparison results. Effect of Brine Recirculation on SEC in BWRO. The use of ERD is a proven technology to reduce SEC in SWRO.20,21,39 It has also been recently demonstrated in BWRO.29 However, in some cases there is leakage of brine into the feed in the ERD, which may increase the feed salinity and reduce the driving force for separation. The same issue exists in brine recirculation systems where energy contained in the brine is recovered. In this subsection, an example is used to quantify the effect of brine feed mixing on SEC in BWRO.

Figure 11. Diagram of RO with ERD.

which yields an overall system recovery (F divided by A) of 85% and 2000 TDS feedwater, whereas the membrane recovery (F divided by E) is at 65% and 4886 TDS feedwater.40 The pressures and flows in all locations are shown in Table 3. Based on these data, it can be derived that α = 0.232, γ = 0.268, and κ = 0.426. If RO properties are the same when ERD is used, α = 0.092, γ = 0.138, and κ = 0.638 because γ ∝ Δπ0/Q0 17246

dx.doi.org/10.1021/ie4031083 | Ind. Eng. Chem. Res. 2013, 52, 17241−17248

Industrial & Engineering Chemistry Research

Article

Table 3. Comparison of Single-Stage BWRO with/without ERD A flow (gpm) pressure (psi) TDS (mg/L)

1774 25 2000

flow (gpm) pressure (psi) TDS (mg/L)

1774 25 2000

B 250 25 2000

C 1524 241 2000

D with PX40 775 221 10563 without PX

1774 241 2000

and κ∝Q20/Δπ0. The model-predicted pressures and flow rates based on the new values of α, γ, and κ are shown in Table 3 for a comparison. Without ERD, the recovery predicted by the model is 90% (see Figure 12), implying that the pump power may be reduced



E

F

G

H

2299 241 4886

1500 5 92

799 231 14000

274 10 14000

1774 241 2000

1604 5 92

170 237 20000

170 237 20000

retentate pressure drop. However, with an interstage booster pump, two-stage may be more energy efficient than single-stage in SWRO. In both BWRO and SWRO studied in this paper, the two-stage RO without booster pump is always no better than the other two configurations from a viewpoint of SEC. • In brine recirculation and/or ERD, mixing of brine and feed may increase salinity, adversely affecting SEC. A trade-off between increased feed salinity and increased energy recovery should be carefully evaluated.

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Phone: +1-909-8693668. Fax: +1-909-869-6920. Notes

The authors declare no competing financial interest.



NOMENCLATURE α Δπ0/ΔP, dimensionless Δπ0 osmotic pressure difference across the membrane at entrance, dimensionless γ ALpΔπ0/Q0, dimensionless κ κQ20/Δπ0, dimensionless π osmotic pressure, bar A membrane area, m2 Ac cross-sectional area for retentate flow, m2 DH hydraulic diameter between membrane sheets, m f Fanning friction coefficient, dimensionless k pressure drop coefficient, bar/(m3/s)2 Lp hydraulic permeability of membrane, m/s/bar Lpv length of pressure vessel, m P pressure, bar P ΔP/ΔP0, dimensionless Pb permeate pressure, bar Pbp pressure before pump, bar Q flow rate, m3 s−1 q Q/Q0, dimensionless Y recovery, dimensionless

Figure 12. Effect of brine leakage in ERD on retentate pressure drop and recovery.

to maintain the same production. Without blending the fresh feed with the highly concentrated brine stream, a lower osmotic pressure is maintained in the retentate. Moreover, a smaller feed flow rate (1774 gpm instead of 2299 gpm) leads to a smaller pressure drop in the retentate stream. Both factors contribute to a larger driving force for permeate production. This case indicates that a trade-off between reduced driving force and increased energy recovery should be carefully evaluated from a viewpoint of SEC when brine recirculation is used. This is especially important for BWRO because the pressure drop effect is significant.



CONCLUSIONS The following conclusions are drawn based on the computational studies of RO desalination in this work: • Even though BWRO and SWRO may be described by the same mathematical model, their process characteristics are different because they are in completely different regions in the γ vs. κ map. The pressure drop effect on RO performance cannot be ignored for BWRO. • In current RO plants, SWRO may be operated near the thermodynamic limit to reduce SEC. However, BWRO is operated far away from the thermodynamic limit because of a small γ. • Without interstage booster pumps, single-stage is more energy efficient than two-stage in BWRO because of less



REFERENCES

(1) Semiat, R. Energy issues in desalination processes. Environ. Sci. Technol. 2008, 42, 8193−8201. (2) Elimelech, M.; Phillip, W. A. The future of seawater desalination: Energy, technology, and the environment. Science 2011, 333, 712−717. (3) Zhu, A.; Christofides, P. D.; Cohen, Y. Effect of thermodynamic restriction on energy cost optimization of RO membrane water desalination. Ind. Eng. Chem. Res. 2009, 48, 6010−6021. (4) Oh, H.; Hwang, T.; Lee, S. A simplified simulation model of RO systems for seawater desalination. Desalination 2009, 238, 128−139.

17247

dx.doi.org/10.1021/ie4031083 | Ind. Eng. Chem. Res. 2013, 52, 17241−17248

Industrial & Engineering Chemistry Research

Article

(28) 5MGD Ground Water Reverse Osmosis Treatment; James City Service Authority Five Forks Water Treatment Facility: James City County, VA; http://www.jccegov.com/pdf/jcsapdfs/Reverse Osmosis Plant Booklet.pdf (29) MacHarg, J. P.; McClellan, S. A. Pressure Exchanger Helps Reduce Energy Costs in Brackish Water RO System. J. AWWA 2004, 96, 44−47. (30) Redondo, J.; Busch, M.; Witte, J. D. Boron removal from seawater using FILMTEC high rejection SWRO membranes. Desalination 2003, 156, 229−238. (31) Hassan, A. M.; Al-Jarrah, S.; Al-Lohabi, T.; Al-Hamdan, A.; Bakheet, L. M.; Amri, M. M. I. Performance Evaluation of SWCC SWRO Plants. In Proceedings of IDA World Conference; Washington, D.C. ; International Desalination Association: Topsfield, MA, 1991; pp 208−224 (32) Hassan, A. M.; Abanmy, A. M.; Al-Thobiety, M.; Mani, T.; AlLuhibi, T.; Al-Masudi, I.; Al-Gherier, A. A.; Bakheet, L.; Amri, M. M. I.; Atiya, K. Performance Evaluation of SWCC SWRO Plants Part II. In Proceedings of IDA World Conference; Washington, D.C. ; International Desalination Association: Topsfield, MA, 1991; pp 412−439 (33) Wilf, M.; Klinko, K. Effective new pretreatment for seawater reverse osmosis systems. Desalination 1998, 117, 323−331. (34) Tarquin, A. J. Desalination and Water Purification Research and Development Program Report No. 108: Volume Reduction of High-Silica RO Concentrate Using Membranes and and Lime Treatment; ,U.S. Department of the Interior Bureau of Reclamation: Denver, CO, 2005. (35) Zhu, A. Energy and Cost Optimization of Reverse Osmosis Desalination. Ph.D. thesis, UCLA, Los Angeles, 2012. (36) Kim, S.; Hoek, E. M. V. Modeling concentration polarization in reverse osmosis processes. Desalination 2005, 186, 111−128. (37) Kurihara, M.; Yamamura, H.; Nakanishi, T. High recovery/high pressure membranes for brine conversion SWRO process development and its performance data. Desalination 1999, 125, 9−15. (38) Lemes, R.; Talavera, J. L. P.; Falcon, R.; Arocha, R.; Curbelo, J.; Platas, V.; Lorenzo, L. D.; Zarzo, D. Evolution of production and energy savings in SWRO plant of Las Palmas III. In Proceedings of IDA World Conference; Taipei, Taiwan ; International Desalination Association: Topsfield, MA, 2011. (39) Penate, B.; de la Fuente, J. A.; Barreto, M. Operation of the RO kinetic energy recovery system: Description and real experiences. Desalination 2010, 252, 179−185. (40) MacHarg, J. P. Energy Optimization of Brackish Groundwater Reverse Osmosis; Texas Water Development Board Contract Report Number 0804830845 ; Texas Water Development Board: Austin, TX, 2011.

(5) Skiborowski, M.; Mhamdi, A.; Kraemer, K.; Marquardt, W. Model-based structural optimization of seawater desalination plants. Desalination 2012, 292, 30−44. (6) Esfahani, I. J.; Ataei, A.; Kim, M.; Kang, O.; Yoo, C. Parametric analysis and optimization of combined gas turbine and reverse osmosis system using refrigeration cycle. Desalin. Water Treat. 2012, 43, 149− 158. (7) Altaee, A. Computational model for estimating reverse osmosis system design and performance: Part-one binary feed solution. Desalination 2012, 291, 101−105. (8) Yechiel, A.; Shevah, Y. Optimization of energy costs for SWRO desalination plants. Desalin. Water Treat. 2012, 46, 304−311. (9) Feinberg, B. J.; Ramon, G. Z.; Hoek, E. M. V. Thermodynamic Analysis of Osmotic Energy Recovery at a Reverse Osmosis Desalination Plant. Environ. Sci. Technol. 2013, 47, 2982−2989. (10) Bartman, A.; Christofides, P. D.; Cohen, Y. Nonlinear Modelbased control of an experimental reverse osmosis water desalination system. Ind. Eng. Chem. Res. 2009, 48, 6126−6136. (11) Bartman, A.; Zhu, A.; Christofides, P. D.; Cohen, Y. Minimizing energy consumption in reverse osmosis membrane desalination using optimization-based control. J. Process Control 2010, 20, 1261−1269. (12) Li, M. Minimization of energy in reverse osmosis water desalination using constrained nonlinear optimization. Ind. Eng. Chem. Res. 2010, 49, 1822−1831. (13) Li, M. Reducing specific energy consumption in reverse osmosis (RO) water desalination: An analysis from first principles. Desalination 2011, 276, 128−135. (14) Li, M. Optimal Plant Operation of Brackish Water Reverse Osmosis (BWRO) Desalination. Desalination 2012, 293, 61−68. (15) Li, M. Optimization of multitrain brackish water reverse osmosis (BWRO) desalination. Ind. Eng. Chem. Res. 2012, 51, 3732−3739. (16) Li, M.; Noh, B. Validation of model-based optimization of brackish water reverse osmosis (BWRO) plant operation. Desalination 2012, 304, 20−24. (17) FILMTEC Membranes System Design: Multi-Stage System; The Dow Chemical Company: Midland, MI; Form No. 609−02050−604. (18) Franks, R.; Wilf, M.; Voutchkov, N.; Murkute, P.; Kizer, J. A Pilot Study Using Seawater Reverse Osmosis Membranes in Combination with Various Pretreatments to Meet the Challenges of Pacific Seawater Desalination; Technical Report, Hydranautics: Oceanside, CA. (19) Schneider, B. Selection, operation and control of a work exchanger energy recovery system based on the Singapore project. Desalination 2005, 184, 197−210. (20) Stover, R. L. Seawater reverse osmosis with isobaric energy recovery devices. Desalination 2007, 203, 168−175. (21) Guirguis, M. J. Energy Recovery Devices in Seawater Reverse Osmosis Desalination Plants with Emphasis on Efficiency and Economical Analysis of Isobaric versus Centrifugal Devices. M.Sc. thesis, University of South Florida, Tampa, FL, 2011. (22) Ayyash, Y.; Imai, H.; Yamada, T.; Fukuda, T.; Yanaga, Y.; Taniyama, T. Performance of reverse osmosis membrane in Jeddah Phase I plant. Desalination 1994, 96, 215−224. (23) Kammourie, N.; Dajanib, T. F. F.; Cioffib, S.; Rybarb, S. SAWACO − North Obhor SWRO Plant Operational Experience. Desalination 2008, 221, 101−106. (24) Sanz, M. A.; Stover, R. L. Low Energy Consumption in the Perth Seawater Desalination Plant. In Proceedings of IDA World Congress; Maspalomas, Gran Canaria, Spain ; International Desalination Association: Topsfield, MA, 2007 (25) Dundorf, S.; MacHarg, J.; Sessions, B.; Seacord, T. F. Optimizing Lower Energy Seawater Desalination. In Proceedings of IDA World Conference; Dubai, United Arab Emirates ; International Desalination Association: Topsfield, MA, 2009. (26) Mohsen, M. S.; Gammoh, S. Performance evaluation of reverse osmosis desalination plant: A case study of Wadi Ma’in, Zara and Mujib Plant. Desalin. Water Treat. 2010, 14, 265−272. (27) Yassin, A. Performance Evaluation of Five Years Operation Experience of WMZM RO Desalination Plant; Taylor & Francis: Abingdon, U.K., 2012. 17248

dx.doi.org/10.1021/ie4031083 | Ind. Eng. Chem. Res. 2013, 52, 17241−17248