Next Article in Journal
Removal of Typical PPCPs by Reverse Osmosis Membranes: Optimization of Treatment Process by Factorial Design
Next Article in Special Issue
Tuning the Properties of Polyvinylidene Fluoride/Alkali Lignin Membranes to Develop a Biocatalytic Membrane Reactor for an Organophosphorus Pesticide Degradation
Previous Article in Journal
Interactions of Surfactants with Biomimetic Membranes—2. Generation of Electric Potential with Non-Ionic Surfactants
Previous Article in Special Issue
Development and Implementation of MBR Monitoring: Use of 2D Fluorescence Spectroscopy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison of Energy Efficiency between Atmospheric Batch Pressure-Retarded Osmosis and Single-Stage Pressure-Retarded Osmosis

1
School of Civil and Environmental Engineering, Nanyang Technological University, 50 Nanyang Avenue, Singapore 639798, Singapore
2
Singapore Membrane Technology Centre, Nanyang Environment and Water Research Institute, Nanyang Technological University, 1 Cleantech Loop, Clean Tech One, #06-08, Singapore 637141, Singapore
3
Interdisciplinary Graduate Programme, Nanyang Technological University, 50 Nanyang Avenue, Singapore 639798, Singapore
*
Author to whom correspondence should be addressed.
Membranes 2023, 13(3), 354; https://doi.org/10.3390/membranes13030354
Submission received: 14 February 2023 / Revised: 2 March 2023 / Accepted: 16 March 2023 / Published: 19 March 2023
(This article belongs to the Special Issue Honorary Issue for Professor Anthony Fane)

Abstract

:
Batch pressure-retarded osmosis (PRO) with varied-pressure and multiple-cycle operation using a pressurized variable-volume tank has been proposed as a high-efficiency osmotic energy harvesting technology, but it suffers scalability constraints. In this study, a more scalable batch PRO, namely, atmospheric batch PRO (AB-PRO), was proposed, utilizing an atmospheric tank to receive and store the intermediate diluted draw solution (DS) and a pressure exchanger to recover the pressure energy from the diluted DS before being recycled into the tank. Its performance was further compared with single-stage PRO (SS-PRO) at different flow schemes via analytic models. The results show that the AB-PRO with an infinitesimal per-cycle water recovery (r) approaches the thermodynamic maximum energy production under ideal conditions, outperforming the SS-PRO with lower efficiencies caused by under-pressurization (UP). However, when considering inefficiencies, a ~40% efficiency reduction was observed in AB-PRO owing to UP and entropy generation as the optimal r is no-longer infinitesimal. Nonetheless, AB-PRO is still significantly superior to SS-PRO at low water recoveries (R) and maintains a stable energy efficiency at various R, which is conducive to meeting the fluctuating demand in practice by flexibly adjusting R. Further mitigating pressure losses and deficiencies of energy recovery devices can significantly improve AB-PRO performance.

1. Introduction

Renewable energy, as an alternative to fossil fuels in power production to reduce carbon emissions, has gained considerable attention in recent decades [1,2]. Among various types of renewable energy, osmotic energy, also called salinity-gradient energy, which originates from the mixing of two solutions with different salinities, has raised interest owing to its huge global capacity and potential accessibility from both natural and industrial sources [3,4,5,6]. The global osmotic power generated from the mixing of river water with seawater was estimated to be 1.6–2.6 TW, of which around 980 GW could be extracted as electricity, which is comparable to that generated from hydropower [7]. Despite the enormous potential, osmotic energy has not been tapped at large scales.
Among various technologies for the harvesting of osmotic energy, membrane-based pressure-retarded osmosis (PRO) [8,9] and reverse electrodialysis (RED) [10,11] have the highest potential for large-scale applications. Simulated and experimental results from previous studies indicate that PRO outperforms RED in terms of higher membrane power density [12,13,14]. In PRO, a hydraulic pressure, lower than the osmotic pressure difference between a higher concentration draw solution (DS, e.g., seawater) and a lower concentration feed solution (FS, e.g., river water), is applied on the DS, which allows DS to “draw” water molecules from FS and results in a pressurized diluted DS with greater volume. In theory, net electric energy can be produced by depressurizing the increased volume of DS through a hydro-turbine. The first PRO plant with a single-stage and co-current flow configuration, as shown in Figure 1a, was designed and installed in Norway by Statkraft [15,16], in which a constant hydraulic pressure is applied throughout the operation. However, it was difficult to obtain a positive net energy production in such a PRO system in practice owing to fouling, extra energy inputs in pre-treatment stages, and inherent inefficiencies of its system design [17,18]. Figure 1a shows the sectional view of the membrane module of this co-current single-stage PRO (SS-PRO), and the dotted curves indicate the osmotic pressure profiles of DS and FS. It reveals the significant reduction of osmotic pressure difference between DS and FS at the outlet because of the simultaneous dilution of DS and concentration of FS caused by water permeation from FS to DS along the same flow direction [19]. For example, using the same volume of a 0.6 M NaCl solution (equivalent to the salinity of seawater) and a 0.01 M NaCl solution (equivalent to the salinity of river water) as the DS and FS, respectively, and assuming 90% of the water in the FS permeates to the DS during the co-current SS-PRO process, the osmotic pressure gradient shrinks from 29.2 bar at the inlet to 10.68 bar at the outlet even if salt leakage is not considered. The highest applied hydraulic pressure is constrained by the resulting low osmotic pressure difference across the membrane at the end of membrane modules. Therefore, the insufficient hydraulic pressure, also called under-pressurization (UP), especially in the first few elements, incurs irreversible energy losses.
To reduce the energy loss caused by UP, the counter-current flow mode with the opposite flow direction of the DS and FS streams is applied to PRO. As illustrated by the dotted curves in Figure 1b, the different directions of FS concentration and DS dilution make the osmotic pressure difference distributed more evenly along the membrane module, which is conducive to the application of a higher pressure, thus improving energy efficiency [20,21]. Under the same conditions discussed in the previous example, the highest hydraulic pressure that can be applied theoretically improves to 15.14 bar in the counter-current SS-PRO. More recently, a multi-stage PRO was developed to reduce energy loss by tailoring the hydraulic pressure close to the osmotic pressure difference [22,23]. A plurality of membrane modules are successively linked by the connection of their FS and/or DS inlets and the outlets of the front stage where the hydraulic pressure decreases gradually from the first to the last stage in accordance with the reducing osmotic pressure difference. However, increasing the number of stages implies additional inter-stage devices and accessories which significantly magnifies the capital cost [24,25]. To reduce energy loss and avoid increasing stages, a batch PRO (BPRO) process was proposed in which a pressurized variable-volume tank is installed on the DS side to allow the recirculation of the diluted DS to residual DS in the tank and realize the multi-cycle operation with the single-stage configuration. The osmotic pressure difference in the membrane module is more evenly distributed due to the smaller water recovery in each cycle of BPRO. The applied hydraulic pressure in BPRO varies with the osmotic pressure difference. The ideal performance of BPRO was evaluated by Li Mingheng [26], and the results demonstrated that BPRO outperforms single-stage and two- or three-stage PRO in terms of both energy production and power density, showing its potential for high-efficiency osmotic energy harvesting without additional stages and expenditures. However, in practical scenarios, scaling up the pressurized variable-volume tank remains challenging, restricting the application of BPRO [27]. Moreover, inefficiency factors, such as the pump, energy recovery device (ERD) deficiencies, and pressure drops, have yet to be considered when evaluating BPRO performance.
In this study, an alternative design of the BPRO by incorporating atmospheric tanks for storing the intermediate DS and FS and a pressure exchanger for recovering pressure energy from the DS effluent as shown in Figure 1c, namely, atmospheric batch PRO (AB-PRO), was proposed, which avoids the use of intricate pressurized variable-volume draw tank in the BPRO reported previously. Specific energy production (SEP) and energy production efficiency (EPE) of AB-PRO were calculated and compared with other PRO technologies operated in a single-stage system (i.e., co-current SS-PRO and counter-current SS-PRO) in both ideal and practical scenarios via analytical modeling.

2. Methods

2.1. Configuration and Operation

As shown in Figure 1a,b, SS-PRO is operated in an open-loop, continuous mode, where both the DS and the FS are discharged out of the system after passing through the membrane module once. The DS influent is pressurized through exchanging hydraulic pressure with the DS effluent inside a pressure exchanger (PX). Owing to the imperfect efficiency of the PX, a circulation (or booster) pump is installed to further pressurize the DS influent to the target value before it enters the membrane module. Under a constant hydraulic pressure below the osmotic pressure difference (∆π), water molecules permeate through the membrane from the FS to the DS, which results in a diluted DS effluent with an increased flow rate. The diluted DS effluent is split into two streams. The first stream, at the same flow rate as the DS influent, is directed to the PX to exchange pressure with the DS influent flowing to the module. The other stream at the same flow rate as the water permeating through the membrane is depressurized by an energy recover device (ERD), such as a hydro-turbine (HT) or PX, to produce electricity or pressurize a fluid in anther system [28,29,30]. The DS and the FS flowing to the membrane module can follow two different flow schemes, the co-current mode (Figure 1a) and the counter-current mode (Figure 1b).
In contrast, as shown in Figure 1c, AB-PRO is operated in a closed-loop, variable-pressure mode with two atmospheric tanks to store DS and FS. In this case, the intermediate diluted DS and the concentrated FS are recycled and mixed with the residual solutions in the DS and FS tanks, respectively, instead of being discharged out of the system. It allows AB-PRO to be operated in multiple cycles with a smaller water recovery in each cycle (i.e., per-cycle water recovery) to achieve the total water recovery requirement. Therefore, the AB-PRO process can start at a higher hydraulic pressure because only a small portion of water passing through the membrane from FS to DS in the first cycle brings about a milder decline in the osmotic pressure difference along the membrane module than that of the SS-PRO, which is also demonstrated by the osmotic pressure profiles (dotted curves) in Figure 1. As DS is diluted and FS is concentrated from cycle to cycle, the applied hydraulic pressure in AB-PRO is gradually reduced to the same level as the constant pressure applied in SS-PRO. Owing to a higher average applied pressure, AB-PRO features less energy loss due to under-pressurization compared to SS-PRO. Moreover, the AB-PRO proposed in this study adopts a PX for the energy exchange between the influents and effluents of DS, which allows the practical atmospheric tank to be used instead of the less scalable pressurized variable-volume tank proposed in the previous study. Both co-current and counter-current flow orientations can be employed in AB-PRO, but the impact of flow schemes is insignificant since the change of osmotic pressure difference in each cycle of AB-PRO is not obvious due to the small per-cycle water recovery (r) as illustrated by the osmotic pressure profile in Figure 1c.

2.2. Derivation of Energy Production

To have an insight into the practical performance of each PRO process, the impacts of inefficiencies of devices (e.g., pump, PX, and ERD) and pressure losses on both the DS and FS sides were systematically assessed. Analytic expressions of the practical maximum specific energy production (SEP) and energy production efficiency (EPE) were derived to evaluate the energy production performance of different PRO modes under ideal and practical conditions. SEP is defined as the energy generated per total volume of DS and FS, and EPE refers to the ratio of the energy extracted by PRO to the Gibbs free energy generated from the mixing of DS and FS. To simplify these analyses, it was assumed that the salt rejection of the membranes is 100%. The osmotic pressure (π) was assumed to be linearly proportional to the salt concentration (c), which can be expressed by the van’t Hoff equation [27,31]:
π = v R g T c
where v is the number of ionic species each solute molecule dissociates into, R g is the ideal gas constant, and T is the absolute temperature, which is assumed to be 298 K in the following calculations.

2.2.1. Specific Gibbs Free Energy

The Gibbs free energy of mixing ( G m ) normalized by the total volume of DS and FS defines the thermodynamic maximum specific energy production ( S E P t h e o m o ). As illustrated in Equations (2)–(4) [19,32], the S E P t h e o r e is related to the concentration of FS and DS, the water recovery, and the volumetric fraction of DS and FS.
S E P t h e o m o = G m V F , 0 + V D , 0 = π D , 0 ( 1 ) ln ( 1 + 1 R ) + π F , 0 ln ( 1 R )
= V F , 0 V F , 0 + V D , 0
R = V P , T V F , 0
where V D , 0 and V F , 0 are the initial volumes of DS and FS, respectively, π D , 0 and π F , 0 are the initial osmotic pressures of DS and FS, respectively, R is the water recovery defined as the ratio of total water permeation volume ( V P , T ) to V F , 0 , which ranges from 0 to 1 (when FS is pure water), and is the volumetric fraction of FS. It should be noted that was optimized according to Equation (A1) in Appendix A in the following calculations.

2.2.2. Energy Production Performance of SS-PRO

In SS-PRO, the highest constant hydraulic pressure ( P C ) equals to the lowest osmotic pressure difference ( π ) along the membrane module. For a given R and in SS-PRO, the practical maximum SEP ( S E P S S P R O ) is obtained at the highest applied constant pressure. Therefore, S E P S S P R O can be computed via Equation (5), and the SEP under ideal conditions ( S E P S S P R O , i d e a l ) can be obtained by Equation (6).
S E P S S P R O = η E R D R P C ( 1 ) η P P C + ( 1 ) η P X η P P C [ η P P F , l o s s + ( 1 ) η P P D , l o s s ]
S E P S S P R O , i d e a l = R P C
where η P , η P X , and η E R D are the efficiencies of pump PX and ERD, respectively, and P D , l o s s and P F , l o s s are the pressure losses on the DS side and the FS side, respectively. Equation (5) also reveals the four energy components contributing to SEP, the specific energy production through the ERD ( η E R D R P C ), the specific energy consumption for pumping ( ( 1 ) η P P C ), the specific energy recovered by the PX ( ( 1 ) η P X η P P C ), and the specific energy loss due to pressure drop ( η P P F , l o s s + ( 1 ) η P P D , l o s s ).
P C is different with various flow schemes in SS-PRO. As shown in Figure 1a, the highest pressure in the co-current SS-PRO ( P C , c o ) depends on the π between the final diluted DS and the final concentrated FS at the end point of the membrane module (Equation (7)):
P C , co = π D , 0 1 + R ( 1 ) π F , 0 ( 1 R )
However, the highest pressure in the counter-current SS-PRO ( P C , c o u n t e r ) is determined by the smaller one between the π at the two ends of the membrane module (i.e., at the outlet of FS or DS).
P C , c o u n t e r = { P 1 = π D , 0 π F , 0 ( 1 R ) ,   a t   t h e   o u t l e t   o f   F S   a n d   w h e n   P 1 < P 2 P 2 = π D , 0 1 + R ( 1 ) π F , 0 ,   a t   t h e   o u t l e t   o f   D S   a n d   w h e n   P 1 > P 2

2.2.3. Energy Production Performance of AB-PRO

Unlike the constant pressure applied in SS-PRO, AB-PRO features variable-pressure according to the changing π during the recirculation of DS and FS. Therefore, the hydraulic pressure applied in AB-PRO is primarily related to the per-cycle water recovery ( r ), rather than the total target water recovery ( R ). In this study, r was assumed to be constant during a batch of AB-PRO. Moreover, the spatial effects in AB-PRO were ignored by assuming a spatially invariant concentration throughout the pressure vessels and tanks at a given moment. Such an assumption facilitates the derivation of the analytic expression of SEP of the AB-PRO without compromised accuracy since the length of the pressure vessel in AB-PRO can be very short and the change of π can be finished in a short time [27]. The impact of flow schemes on SEP was also ignored based on this assumption. The maximum variable pressure applied in AB-PRO is a function of the time ( t ) as follows:
P t = π D , 0 ( 1 + r θ ) [ 1 + ( 1 ) t τ ] π F , 0 ( 1 r ) ( 1 t τ )
r = Q P Q F
θ = Q D Q F
τ = V F , 0 Q P
where Q F and Q D are the FS and DS circulation rates, respectively, Q P is the water permeation rate, τ is defined as the maximum retention time of the FS under the permeate flow rate of Q P , θ is the ratio of Q F to Q D in AB-PRO, which was assumed to be the same as that in SS-PRO.
The approximation of practical maximum SEP of AB-PRO can be obtained via Equation (13). It should be noted that the pressure loss on both sides of DS and FS in AB-PRO is assumed to be 1/7 times that of SS-PRO due to the shorter module length of AB-PRO.
S E P A B P R O = [ ( η P X 1 ) θ η P r + η E R D ] [ π D , 0 ( 1 ) ( 1 + r θ ) ln ( 1 + 1 R ) + π F , 0 ( 1 r ) ln ( 1 R ) ] R η P r ( θ P D , l o s s + P F , l o s s )
By differentiating the S E P A B P R O with respect to r , the optimal r can be computed by the iterative method. Under ideal conditions, the maximum SEP (Equation (A14) in Appendix B) can be obtained with r 0 , which equals the specific Gibbs free energy of mixing (Equation (2)).
The SEP of different PRO processes was compared under different scenarios. The values of parameters, including η P , η P X , η E R D , P D , l o s s , and P F , l o s s , are listed in Table 1.

3. Results and Discussions

3.1. Ideal Energy Production Performance

Figure 2 shows the SEP and EPE of different PRO processes with various pairs of DS and FS. The results indicate that the SEP increases with a higher salinity gradient between the DS and the FS in all three PRO processes.
As illustrated in Figure 2a,b, both co-current flow and counter-current flow SS-PRO processes show imperfect energy production efficiency (EPE) even under ideal conditions. The constant-pressure operating mode of SS-PRO results in unavoidable energy loss caused by under-pressurization (UP), which is indicated by the yellow areas in Figure 2d,e. The energy loss increases with the increase of R , leading to the decrease of EPE. However, there are still differences between the two SS-PRO processes, which have different flow schemes. For co-current SS-PRO, as R increases, the SEP first increases to a peak then decreases. According to Equation (6), without considering the inefficiencies, the SEP of SS-PRO only depends on (1) the applied pressure ( P C ), (2) the water recovery ( R ), and (3) the volumetric fraction of FS ( ). For a specific pair of DS and FS at the optimal (Equation (A1) in Appendix A), when R increases, there is a tradeoff between the increased loss of energy owing to the drop of P C and the increased release of Gibbs free energy of mixing ( G m ) due to more water permeation. When the adverse impact of UP dominates, both S E P and E P E drop significantly with the increasing R as a result of the declining P C owing to the rapid decrease of osmotic pressure difference ( π ) along the module in the co-current SS-PRO. When R is maximized, the osmotic pressure of the final concentrated FS is the same as that of the final diluted DS, which means that no hydraulic pressure can be applied to extract energy although G m reaches maximum. Unlike the trend of SEP, the EPE keeps decreasing as R rises. Theoretically, to obtain a high SEP and EPE simultaneously, it is necessary to operate the co-current SS-PRO at low R using a pair of DS and FS with a high salinity gradient. For example, the SEP of 0.10 kWh·m−3 and the EPE of 65% at R = 0.50 using a 0.6 M NaCl solution as the DS and a 0.05 M NaCl solution as the FS can be improved to 0.14 kWh·m−3 and 86%, respectively, by changing the DS into a 1.2 M NaCl solution and reducing the R to 0.20. For the counter-current SS-PRO, the flow scheme induces a more even distribution of π along the membrane module. As illustrated by the red dotted curve in Figure 2e describing the change of π in the counter-current mode, it starts from a lower value and ends at a higher value compared to the black dash curve showing the change of π in the co-current mode. It can be explained by the osmotic pressure difference profile depicted by the dotted curves in Figure 1c that the opposite directions of DS and FS flows in the counter-current SS-PRO result in a lower π of ( π D , i n π D , o u t ) at the DS inlet than that ( π = π D , i n π F , i n ) in the co-current SS-PRO, and a higher π of ( π D , o u t π D , i n ) at the DS outlet than that ( π = π D , o u t π F , o u t ) in the co-current SS-PRO. Since the hydraulic pressure applied in SS-PRO is constrained by the lowest π , the counter-current flow mode contributes to a higher appliable hydraulic pressure and more extractable osmotic energy at the same R . Equations (7) and (8) also demonstrate a higher P C , c o u n t e r than P C , c o from a mathematical perspective. Therefore, the counter-current SS-PRO outperforms the co-current SS-PRO on both SEP and EPE especially at high R . In addition, as a considerable concentration change is required to bring π to 0 (e.g., π D , i n = π F , o u t , or π D , o u t = π F , i n ) in the counter-current SS-PRO, the situation where no pressure can be applied to harvest energy will not occur in most cases. As shown in Figure 1b, the SEP of counter-current SS-PRO is monotonically increasing with R ,   and the EPE remains above 70% at any R .
In contrast, as shown in Figure 2c, AB-PRO can approach the thermodynamic maximum SEP and 100% of EPE under ideal conditions, which are significantly higher than those of the two SS-PRO processes. The mechanism illustrated in Figure 2f reveals that the energy loss caused by UP can be eliminated in the AB-PRO by operating with an infinitesimal per-cycle water recovery ( r ), where the variable hydraulic pressure applied on the DS side is as close as possible to the changing osmotic pressure difference. Using the 1.2 M NaCl solution and 0.05 M NaCl solution as DS and FS, respectively, 0.34 kWh·m−3 of power can be produced in AB-PRO at R = 0.5, which is equivalent to that of a 120 m-high hydropower system.
Theoretical results suggest that AB-PRO is potentially superior to SS-PRO benefiting from the reversible mixing process and variable-pressure operation. The following sections will further assess its energy production when considering practical inefficiencies.

3.2. Impact of Inefficiencies

Inefficiencies as listed in Table 1, including pump efficiency ( η P ), pressure exchanger efficiency ( η P X ), energy recovery device efficiency ( η E R D ), and pressure loss ( P l o s s ) on the DS and FS sides, were considered when evaluating the practical performance of different PRO processes. The pressure drops on DS and FS sides were assumed to be the same (i.e., P l o s s = P F , l o s s = P D , l o s s ), and the P l o s s in AB-PRO was assumed to be 1/7 of that in SS-PRO as a shorter membrane module is allowed to be utilized in AB-PRO because of the smaller water recovery in each cycle. The parameters of η P of 0.85, η P X of 0.98, η E R D of 0.90, P l o s s of 0.2 bar in AB-PRO and P l o s s of 1.4 bar in SS-PRO were set as the baseline conditions. The impact of each parameter among η P , η E R D , and P l o s s was investigated by the control variate method.

3.2.1. Overall Impact of Inefficiencies

As revealed by Figure 3, inefficiencies induce significant reductions in the SEP and EPE of all three PRO processes. For SS-PRO with both flow modes, the energy production cannot compensate the unavoidable energy loss caused by device inefficiencies and pressure drops in practice at low R . An R of at least 0.1 is required to gain the net energy production as indicated by Figure 3a,b. When R ranges from 0.1 to 0.25, the SEP and EPE have an insignificant difference between the co-current SS-PRO and counter-current SS-PRO, which rises with the increase of R in both SS-PRO processes. At a higher   R , the EPE of the co-current SS-PRO decreases due to the rapid decrease of ∆ π at the outlet of pressure vessels while that of the counter-current SS-PRO remains stable (above 50%) benefitting from the more uniform distribution of ∆ π along the membrane module.
A significant energy loss of over 40% also occurs in AB-PRO when considering practical inefficiencies (baseline). As the r is no longer approaching 0 in practical scenarios (e.g., the optimal r is ~0.15 in the baseline case), the maximum hydraulic pressure that can be applied will be lower according to Equation (9), leading to more energy loss caused by under-pressurization (UP). Moreover, there is a larger salinity gradient between the recirculated solution and residual solution in both DS and FS tanks when r increases, hence increasing the amount of entropy generation by mixing and energy loss. Furthermore, there is a tradeoff in AB-PRO: where a higher r leads to an increase of energy loss caused by UP and mixing, but a reduction of energy loss by pressure drops and device inefficiencies as the total volume of solution passing through the devices and membrane module is diminished (Equations (A8) and (A9) in Appendix B). Therefore, the SEP of AB-PRO was optimized by optimizing r at each R in this study, including the data in Figure 3c and other figures showing SEP of AB-PRO in the following sections. Although AB-PRO features a similar EPE (within a 5% difference) as the counter-current SS-PRO, its SEP and EPE are significantly higher than those of both SS-PRO processes at R < 0.4. In addition, the performance of AB-PRO is relatively stable. The change of EPE with R does not exceed 10%, which facilitates the flexible adjustment of R to meet the fluctuating electricity or pressure demand in application. When R = 0.5 in the baseline case, SEP and EPE are 0.19 kWh·m−3 and 56%, respectively. It should be noted that the value of SEP is related to the normalization method. The total energy production is normalized by the total volume of DS and FS in this study, but the SEP can also be defined as the energy production per volume of DS or FS to show the osmotic energy harvesting capacity if FS or DS is not limited. The results normalized by the initial DS volume (Figure A1 in Appendix C) exhibit a higher SEP of 0.48 kWh·m−3 at R = 0.5 under the same conditions.

3.2.2. Impact of Pump Efficiency

By comparing the PRO performances at the same mode but different η P , the impact of pump on SEP and EPE is not significant as demonstrated in Figure 4. In all three PRO processes, an improvement of pump efficiency from 0.80 to 0.95 only results in an increase in EPE by less than 5%. The results are mainly contributed by the configuration of a high-efficiency pressure exchanger ( η P X = 0.98) which recovers most of the pressure generated by pumps and alleviates the negative effect of pump inefficiency, while the pump used in PRO is only for boosting the pressure after the PX, typically for several bars (<5 bar). The insignificant impact of η P implies more freedom in the choice of pump quality in practice. Moreover, as η P generally varies at various flow rates, the results also indicate a strong stability of the three PRO systems under different operating conditions.

3.2.3. Impact of Energy Recovery Device Efficiency

As shown in Figure 1, an energy recovery device (ERD), such as a hydro-turbine (HT) or a pressure exchanger (PX), is installed in PRO processes to produce electricity or transfer the pressure to another system. The impact of η E R D on the SEP and EPE is significant, as demonstrated by Figure 5. With an increase in η E R D , the energy production performance is substantially improved in all three PRO processes, while the trends of performance curves for all the processes remain consistent with those in the baseline case. The maximum EPE increment in both SS-PRO and AB-PRO can approach 9% when η E R D increases from 0.85 to 0.95. By using a HT with a high efficiency of 0.95, totals of 0.19 kWh·m−3, 0.22 kWh·m−3, and 0.23 kWh·m−3 of electricity can be generated by the co-current SS-RO, counter-current SS-PRO, and AB-PRO, respectively, at DS = 1.2 M NaCl solution, FS = 0.05 M NaCl solution, and R = 0.60 . Moreover, if there is another system requiring pressure, η E R D can be further enhanced to 0.98 by replacing the HT by a PX. In this case, the EPE of ~59% can be maintained in both the AB-PRO and counter-current SS-PRO processes when R > 0.8. A higher EPE ranging from 60 to 66% can be realized in AB-PRO at R < 0.8, which is superior to that in the two SS-PRO processes, especially at low R .

3.2.4. Impact of Pressure Loss

Pressure loss ( P l o s s ) caused by the friction of fluid in the flow channel (e.g., in membrane modules) of both DS and FS induces energy loss and lessens the net energy extractable by PRO. When assuming the same optimal hydrodynamic conditions (e.g., hydraulic diameter, flow velocity, and regimes) for SS-PRO and AB-PRO, the P l o s s is proportional to the module length according to the Darcy-Weisbach equation [18,27]. In the following investigation, the P l o s s of AB-PRO was assumed to be 1/7 of that of SS-PRO (Table 1) since batch PRO can use much shorter membrane modules to achieve a similar recovery in multiple-cycle operation. Moreover, the membrane module was assumed to have high mechanical strength, hence there was no membrane deformation in the PRO processes [38].
The results in Figure 6 indicate a significant decay of the energy harvesting performance when P l o s s increases. In SS-PRO, a higher P l o s s requires a higher R threshold to gain net energy production and leads to lower SEP and EPE at the same R . For example, if P l o s s rises from 1.4 bar to 3.5 bar, the minimal R for net energy generation increases from ~0.1 to ~0.2 in both SS-PRO processes. Meanwhile, the reduction of SEP ranges from 15% to 42% at different recoveries, and the highest EPE is only ~30% at R = 0.5 in the co-current mode and ~40% at R = 0.9 in the counter-current mode. On the contrary, the EPE can increase by 5–25% if P l o s s can be reduced to 0.7 bar in SS-PRO through optimizing hydrodynamic conditions, feed spacer geometry, and so on [38,39]. The counter-current SS-PRO always outperforms the co-current SS-PRO at various P l o s s and is capable of maintaining a relatively stable EPE of ~60% at R > 0.30 and P l o s s = 0.7 bar. At the P l o s s of 0.5 bar, AB-PRO presents a better performance than SS-PRO (at the P l o s s of 3.5 bar) and realizes an EPE of 40–50% benefitting from the shorter module length. However, alleviating P l o s s from 0.2 bar (baseline) to 0.1 bar only results in a 2–5% increase in EPE in AB-PRO. Although the decrease of per-cycle water recovery ( r ) from ~0.22 to ~0.15 because of the reduction in P l o s s brings about less energy loss caused by UP and entropy generation, it also induces more energy loss owing to inefficiencies in AB-PRO, resulting in a slight improvement of EPE. When P l o s s = 0.1 bar in AB-PRO and 0.7 bar in SS-PRO, the counter-current SS-PRO shows a small advantage compared to AB-PRO at R > 0.3, but the EPE difference between them does not exceed 3%. AB-PRO has an absolute predominance at R < 0.30.

4. Conclusions

The performance of atmospheric batch PRO (AB-PRO) was explored and compared with conventional single-stage PRO (SS-PRO) with different flow schemes. Variable-pressure AB-PRO with an infinitesimal per-cycle water recovery ( r ) can approach the thermodynamic maximum SEP and 100% of EPE under ideal conditions, while the efficiencies of two SS-PRO processes decrease with the increasing R owing to the irreversible energy loss of under-pressurization caused by the constant-pressure operation. The impact of inefficiencies, including device deficiencies and pressure losses, was also investigated for all three PRO processes. In the practical case, a significant decay of performance was observed in all three PRO processes. Although counter-current SS-PRO shows a comparable performance with AB-PRO at high R , the SEP and EPE of AB-PRO are significantly higher than that of SS-PRO at low R . AB-PRO is capable of maintaining a relatively stable and high efficiency in the entire range of R , which facilitates meeting the fluctuating energy or pressure demand in application by the flexible adjustment of R . When utilizing a 1.2 M NaCl solution as the DS and a 0.05 M NaCl solution as the FS, the SEP and EPE of AB-PRO reach 0.19 kWh·m−3 and 56%, respectively, in the baseline case at R = 0.5. For either AB-PRO or SS-PRO process, the advancement in pump and ERDs together with process design are critical to the further enhancement of energy production capacity. With a fixed high pressure exchanger efficiency of 0.98, pressure loss plays the most important role in the overall efficiency of PRO, followed by the ERD efficiency and pump efficiency.

Author Contributions

Conceptualization, D.L. and Q.S.; methodology, D.L. and Q.S.; software, D.L.; validation, D.L. and Z.M.; formal analysis, D.L. and Q.S.; investigation, D.L. and Q.S.; data curation, D.L.; writing—original draft preparation, D.L.; writing—review and editing, D.L., Z.M. and Q.S.; visualization, D.L.; supervision, Q.S.; project administration, Q.S.; funding acquisition, Q.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Education, Singapore, under the Academic Research Fund Tier 1 [RG123/21] and the Singapore Energy Centre [SgEC-Core2021-44].

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the Ministry of Education, Singapore and the Singapore Energy Centre for supporting and funding this research work.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

AbbreviationMeaning
PROPressure-retarded osmosis
REDReverse electrodialysis
DSDraw solution
FSFeed solution
SS-PROSingle-stage pressure-retarded osmosis
BPROBatch pressure-retarded osmosis
AB-PROAtmospheric batch pressure-retarded osmosis
SEPSpecific energy production
EPEEnergy production efficiency
ERDEnergy recovery device
PXPressure exchanger
HTHydro-turbine
UPUnder-pressurization

Nomenclature

SymbolMeaningUnit
S E P t h e o m o Thermodynamic maximum specific energy productionkWh·m−3
S E P S S P R O Specific energy production of single-stage PROkWh·m−3
S E P A B P R O Specific energy production of atmospheric batch PROkWh·m−3
S E P S S P R O , i d e a l Ideal specific energy production of single-stage PROkWh·m−3
S E P A B P R O , i d e a l Ideal specific energy production of atmospheric batch PROkWh·m−3
G m Gibbs free energy of mixingkWh·m−3
R g Ideal gas constantJ·K−1·mol−1
v Number of ionic species-
T Absolute temperatureK
Volumetric fraction of feed solution-
o p t i m a l Optimal volumetric fraction of feed solution-
R Water recovery-
r Per-cycle water recovery-
V D , 0 Initial volume of draw solutionm3
V F , 0 Initial volume of feed solutionm3
V P , T Total water permeation volumem3
Q D Circulation rate of draw solutionm3·s−1
Q F Circulation rate of feed solutionm3·s−1
Q P Flow rate of permeatem3/s
π Osmotic pressurebar
π D , 0 Initial osmotic pressure of draw solutionbar
π F , 0 Initial osmotic pressure of feed solutionbar
π Osmotic pressure differencebar
η P Pump efficiency-
η P X Pressure exchanger efficiency-
η E R D Energy recovery device efficiency-
P C Constant hydraulic pressurebar
P C , c o Constant hydraulic pressure in co-current flow modebar
P C , c o u n t e r Constant hydraulic pressure in counter-current flow modebar
P t Hydraulic pressure at time tbar
P l o s s Pressure lossbar
P D , l o s s Pressure loss in draw solution sidebar
P F , l o s s Pressure loss in feed solution sidebar
θ Ratio of to Q D to Q F -
τ Ratio of V F , 0 to Q P
E p u m p Energy consumption by pumpingkWh
E P X Energy recovered by pressure exchangerkWh
E l o s s Energy consumption by compensating pressure losskWh
E E R D Energy production/recovery by energy recovery devicekWh

Appendix A. The Optimal Volumetric Fraction

The volumetric fraction of FS was optimized to achieve the highest energy production from the mixing of FS and DS in this study. Based on Equation (1) in the main text, the optimal can be obtained by the following equation:
o p t i m a l = e [ π F , 0 ln ( π F , 0 ) π D , 0 ln ( π D , 0 ) ( π F , 0 π D , 0 ) 1 ] π D , 0 ( π F , 0 π D , 0 )

Appendix B. Derivations of Specific Energy Production

Appendix B.1. SS-PRO

SS-PRO is a continuous, constant-pressure process; hence the derivation of energy production is straightforward. The energy components in SS-PRO include:
(1)
Energy consumption by pumping ( E p u m p )
E p u m p = 1 η P P C V D , 0
(2)
Energy recovered by the pressure exchanger ( E P X )
E P X = η P X η P P C V D , 0
(3)
Energy consumption by compensating pressure loss ( E l o s s )
E l o s s = 1 η P P F , l o s s V F , 0 + 1 η P P D , l o s s V D , 0
(4)
Energy production by an energy recovery device ( E E R D )
E E R D = η E R D P C V P
The SEP is the net energy production normalized by the total volume of DS and FS.
S E P = E p u m p + E P X E l o s s + E E R D V F , 0 + V D , 0
Equation (A6) can be simplified to Equation (5) in the main text.
In co-current flow SS-PRO, P c is the osmotic pressure difference between the final diluted DS ( π D , f ) and the final concentration FS ( π F , f ) at the end point of the membrane module with P C , c o = ( π D , f π F , f ) , π D , f = π D , 0 1 + R ( 1 ) and π F , f = π F , 0 ( 1 R ) . In counter-current flow SS-PRO, P C is lower between the osmotic pressure difference at the FS outlet ( P 1 = π D , 0 π F , f = π D , 0 π F , 0 ( 1 R ) ) and at the DS outlet ( P 2 = π D , f π F , 0 = π D , 0 1 + R ( 1 ) π F , 0 ). In the ideal case, η P = η P X = η E R D = 1, P F , l o s s = P D , l o s s = 0 bar.

Appendix B.2. AB-PRO

AB-PRO is a continuous but variable-pressure ( P t ) process, hence the derivation of its energy production is based on the integration with respect to time ( t ). The duration time ( t B ) of a batch of AB-PRO can be obtained by Equation (A7).
t B = R V F , 0 r Q F
It was assumed that the total water permeation rate ( Q P ) is constant during the AB-PRO process. Thus, regarding a constant recirculation flow rate of FS ( Q F ), r is also constant.
The volume of solution passing through the devices and membrane module on the FS side ( V p a s s , F ) and DS side ( V p a s s , D ) in t B can be obtained by Equations (A8) and (A9), respectively.
V p a s s , F = Q F t B = R V F , 0 r
V p a s s , D = Q D t B = θ R V F , 0 r
The energy components in AB-PRO include:
(1)
Energy consumption by pumping ( E p u m p )
E p u m p = 1 η P 0 t B Q D P t d t
(2)
Energy recovered by pressure exchanger ( E P X )
E P X = η P X η P 0 t B Q D P t d t
(3)
Energy consumption by compensating pressure loss ( E l o s s )
E l o s s = 1 η P P D , l o s s Q D t B + 1 η P P F , l o s s Q F t B
(4)
Energy production by energy recovery device ( E E R D )
E E R D = η E R D 0 t B Q P P t d t
The SEP of AB-PRO can also be obtained by Equation (A6) combining with Equations (9)–(12) and further simplified to Equation (13) in the main text.
In the ideal case, not only η P = η P X = η E R D = 1 and P F , l o s s = P D , l o s s = 0 bar, but also r 0 . The highest SEP of AB-PRO can be approached by the following equation, which is the same as the thermodynamic maximum SEP (Equation (2) in the main text).
S E P A B P R O , i d e a l = π D , 0 ( 1 ) ln ( 1 + 1 R ) + π F , 0 ln ( 1 R )

Appendix C. Specific Energy Production Normalized by the Initial Volume of Draw Solution

Figure A1. Specific energy production (SEP) normalized by initial volume of draw solution (DS) and energy production efficiency (EPE) of (a) co-current single-stage PRO (SS-PRO), (b) counter-current SS-PRO, and (c) atmospheric batch PRO (AB-PRO) under ideal and practical baseline conditions. The feed solution (FS) is a 0.05 M NaCl solution, and the DS is a 1.2 M NaCl solution. Practical baseline conditions: η P = 0.85 , η E R D = 0.90 , η P X = 0.98 , and P l o s s = 0.2   bar in AB-PRO and 1.4 bar in SS-PRO.
Figure A1. Specific energy production (SEP) normalized by initial volume of draw solution (DS) and energy production efficiency (EPE) of (a) co-current single-stage PRO (SS-PRO), (b) counter-current SS-PRO, and (c) atmospheric batch PRO (AB-PRO) under ideal and practical baseline conditions. The feed solution (FS) is a 0.05 M NaCl solution, and the DS is a 1.2 M NaCl solution. Practical baseline conditions: η P = 0.85 , η E R D = 0.90 , η P X = 0.98 , and P l o s s = 0.2   bar in AB-PRO and 1.4 bar in SS-PRO.
Membranes 13 00354 g0a1

References

  1. Dincer, I. Renewable energy and sustainable development: A crucial review. Renew. Sustain. Energy Rev. 2000, 4, 157–175. [Google Scholar] [CrossRef]
  2. Solomon, S.; Daniel, J.S.; Sanford, T.J.; Murphy, D.M.; Plattner, G.K.; Knutti, R.; Friedlingstein, P. Persistence of climate changes due to a range of greenhouse gases. Proc. Natl. Acad. Sci. USA 2010, 107, 18354–18359. [Google Scholar] [CrossRef] [Green Version]
  3. Wick, G.L. Power from salinity gradients. Energy 1978, 3, 95–100. [Google Scholar] [CrossRef]
  4. Yip, N.Y.; Brogioli, D.; Hamelers, H.V.; Nijmeijer, K. Salinity Gradients for Sustainable Energy: Primer, Progress, and Prospects. Environ. Sci. Technol. 2016, 50, 12072–12094. [Google Scholar] [CrossRef]
  5. Alvarez-Silva, O.A.; Osorio, A.F.; Winter, C. Practical global salinity gradient energy potential. Renew. Sustain. Energy Rev. 2016, 60, 1387–1395. [Google Scholar] [CrossRef]
  6. Liu, S.; Song, W.; Meng, M.; Xie, M.; She, Q.; Zhao, P.; Wang, X. Engineering pressure retarded osmosis membrane bioreactor (PRO-MBR) for simultaneous water and energy recovery from municipal wastewater. Sci. Total Environ. 2022, 826, 154048. [Google Scholar] [CrossRef] [PubMed]
  7. Ramon, G.Z.; Feinberg, B.J.; Hoek, E.M.V. Membrane-based production of salinity-gradient power. Energy Environ. Sci. 2011, 4, 4423–4434. [Google Scholar] [CrossRef]
  8. Gonzales, R.R.; Abdel-Wahab, A.; Adham, S.; Han, D.S.; Phuntsho, S.; Suwaileh, W.; Hilal, N.; Shon, H.K. Salinity gradient energy generation by pressure retarded osmosis: A review. Desalination 2021, 500, 114841. [Google Scholar] [CrossRef]
  9. 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]
  10. Post, J.W.; Hamelers, H.V.M.; Buisman, C.J.N. Energy recovery from controlled mixing salt and fresh water with a reverse electrodialysis system. Environ. Sci. Technol. 2008, 42, 5785–5790. [Google Scholar] [CrossRef]
  11. Mei, Y.; Tang, C.Y. Recent developments and future perspectives of reverse electrodialysis technology: A review. Desalination 2018, 425, 156–174. [Google Scholar] [CrossRef]
  12. Yip, N.Y.; Elimelech, M. Comparison of Energy Efficiency and Power Density in Pressure Retarded Osmosis and Reverse Electrodialysis. Environ. Sci. Technol. 2014, 48, 11002–11012. [Google Scholar] [CrossRef] [PubMed]
  13. Post, J.W.; Veerman, J.; Hamelers, H.V.M.; Euverink, G.J.W.; Metz, S.J.; Nymeijer, K.; Buisman, C.J.N. Salinity-gradient power: Evaluation of pressure-retarded osmosis and reverse electrodialysis. J. Membr. Sci. 2007, 288, 218–230. [Google Scholar] [CrossRef] [Green Version]
  14. Ju, J.; Choi, Y.; Lee, S.; Jeong, N. Comparison of fouling characteristics between reverse electrodialysis (RED) and pressure retarded osmosis (PRO). Desalination 2021, 497, 114648. [Google Scholar] [CrossRef]
  15. Crown Princess of Norway to Open the World’s First Osmotic Power Plant. Available online: http://www.statkraft.com/presscentre/press-releases/crown-princess-mette-marit-to-open-the-worlds-first-osmotic-power-plant.aspx (accessed on 1 March 2021).
  16. Achilli, A.; Childress, A.E. Pressure retarded osmosis: From the vision of Sidney Loeb to the first prototype installation—Review. Desalination 2010, 261, 205–211. [Google Scholar] [CrossRef]
  17. Straub, A.P.; Lin, S.; Elimelech, M. Module-scale analysis of pressure retarded osmosis: Performance limitations and implications for full-scale operation. Environ. Sci. Technol. 2014, 48, 12435–12444. [Google Scholar] [CrossRef]
  18. Kim, M.; Kim, S. Practical limit of energy production from seawater by full-scale pressure retarded osmosis. Energy 2018, 158, 373–382. [Google Scholar] [CrossRef]
  19. Straub, A.P.; Deshmukh, A.; Elimelech, M. Pressure-retarded osmosis for power generation from salinity gradients: Is it viable? Energy Environ. Sci. 2016, 9, 31–48. [Google Scholar] [CrossRef]
  20. He, W.; Wang, Y.; Shaheed, M.H. Modelling of osmotic energy from natural salt gradients due to pressure retarded osmosis: Effects of detrimental factors and flow schemes. J. Membr. Sci. 2014, 471, 247–257. [Google Scholar] [CrossRef]
  21. Yang, W.; Song, L.; Zhao, J.; Chen, Y.; Hu, B. Numerical analysis of performance of ideal counter-current flow pressure retarded osmosis. Desalination 2018, 433, 41–47. [Google Scholar] [CrossRef]
  22. Touati, K.; Tadeo, F.; Elfil, H. Osmotic energy recovery from Reverse Osmosis using two-stage Pressure Retarded Osmosis. Energy 2017, 132, 213–224. [Google Scholar] [CrossRef] [Green Version]
  23. He, W.; Wang, Y.; Shaheed, M.H. Enhanced energy generation and membrane performance by two-stage pressure retarded osmosis (PRO). Desalination 2015, 359, 186–199. [Google Scholar] [CrossRef]
  24. 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. [Google Scholar] [CrossRef]
  25. Lin, S.; Elimelech, M. Staged reverse osmosis operation: Configurations, energy efficiency, and application potential. Desalination 2015, 366, 9–14. [Google Scholar] [CrossRef]
  26. Li, M. Dynamic Operation of Batch Reverse Osmosis and Batch Pressure Retarded Osmosis. Ind. Eng. Chem. Res. 2020, 59, 3097–3108. [Google Scholar] [CrossRef]
  27. Werber, J.R.; Deshmukh, A.; Elimelech, M. Can batch or semi-batch processes save energy in reverse-osmosis desalination? Desalination 2017, 402, 109–122. [Google Scholar] [CrossRef] [Green Version]
  28. Sim, V.S.T.; She, Q.; Chong, T.H.; Tang, C.Y.; Fane, A.G.; Krantz, W.B. Strategic Co-Location in a Hybrid Process Involving Desalination and Pressure Retarded Osmosis (PRO). Membranes 2013, 3, 98–125. [Google Scholar] [CrossRef] [Green Version]
  29. 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]
  30. Lee, S.; Park, T.-S.; Park, Y.-G.; Lee, W.-I.; Kim, S.-H. Toward scale-up of seawater reverse osmosis (SWRO)–pressure retarded osmosis (PRO) hybrid system: A case study of a 240 m3/day pilot plant. Desalination 2020, 491, 114429. [Google Scholar] [CrossRef]
  31. Mo, Z.; Li, D.; She, Q. Semi-closed reverse osmosis (SCRO): A concise, flexible, and energy-efficient desalination process. Desalination 2022, 544, 116147. [Google Scholar] [CrossRef]
  32. Lin, S.; Straub, A.P.; Elimelech, M. Thermodynamic limits of extractable energy by pressure retarded osmosis. Energy Environ. Sci. 2014, 7, 2706–2714. [Google Scholar] [CrossRef]
  33. Madsen, H.T.; Bruun Hansen, T.; Nakao, T.; Goda, S.; Søgaard, E.G. Combined geothermal heat and pressure retarded osmosis as a new green power system. Energy Convers. Manag. 2020, 226, 113504. [Google Scholar] [CrossRef]
  34. Manzoor, H.; Selam, M.A.; Abdur Rahman, F.B.; Adham, S.; Castier, M.; Abdel-Wahab, A. A tool for assessing the scalability of pressure-retarded osmosis (PRO) membranes. Renew. Energy 2020, 149, 987–999. [Google Scholar] [CrossRef]
  35. Huang, B.; Pu, K.; Wu, P.; Wu, D.; Leng, J. Design, Selection and Application of Energy Recovery Device in Seawater Desalination: A Review. Energies 2020, 13, 4150. [Google Scholar] [CrossRef]
  36. Bharadwaj, D.; Struchtrup, H. Large scale energy storage using multistage osmotic processes: Approaching high efficiency and energy density. Sustain. Energy Fuels 2017, 1, 599–614. [Google Scholar] [CrossRef]
  37. 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]
  38. She, Q.; Hou, D.; Liu, J.; Tan, K.H.; Tang, C.Y. Effect of feed spacer induced membrane deformation on the performance of pressure retarded osmosis (PRO): Implications for PRO process operation. J. Membr. Sci. 2013, 445, 170–182. [Google Scholar] [CrossRef]
  39. Benjamin, J.; Al Mashrafi, S.; Tejada-Martinez, A.; Diaz-Elsayed, N.; Arias, M.E.; Zhang, Q. Optimizing pressure retarded osmosis spacer geometries: An experimental and CFD modeling study. J. Membr. Sci. 2022, 647, 120284. [Google Scholar] [CrossRef]
Figure 1. Schematic diagrams and sectional views of (a) co-current single-stage PRO (SS-PRO), (b) counter-current SS-PRO, (c) atmospheric batch PRO (AB-PRO). The dotted curves in the sectional views show the osmotic pressure profile along the membrane module. DS: draw solution, FS: feed solution, CP: circulation pump, PX: pressure exchanger, ERD: energy recovery device (i.e., PX or hydro-turbine), π D , i n and π D , o u t : osmotic pressure of DS influents and effluents, respectively, π F , i n and π F , o u t : osmotic pressure of FS influents and effluents, respectively.
Figure 1. Schematic diagrams and sectional views of (a) co-current single-stage PRO (SS-PRO), (b) counter-current SS-PRO, (c) atmospheric batch PRO (AB-PRO). The dotted curves in the sectional views show the osmotic pressure profile along the membrane module. DS: draw solution, FS: feed solution, CP: circulation pump, PX: pressure exchanger, ERD: energy recovery device (i.e., PX or hydro-turbine), π D , i n and π D , o u t : osmotic pressure of DS influents and effluents, respectively, π F , i n and π F , o u t : osmotic pressure of FS influents and effluents, respectively.
Membranes 13 00354 g001
Figure 2. Specific energy production (SEP) and energy production efficiency (EPE) of (a) co-current single-stage PRO (SS-PRO), (b) counter-current SS-PRO, and (c) atmospheric batch PRO (AB-PRO), and energy components of (d) co-current SS-PRO, (e) counter-current SS-PRO, and (f) AB- PRO. In figure (ac), various salt solutions with different concentrations were used for osmotic energy harvesting including a 0.6 M NaCl solution (e.g., equivalent to the salinity of seawater) and a 1.2 M NaCl solution (e.g., equivalent to the salinity of SWRO brine) as draw solution (DS) alternatives, a 0.01 M NaCl solution (e.g., equivalent to the salinity of river water) and a 0.05 M NaCl solution (e.g., equivalent to the salinity of wastewater concentrate) as feed solution (FS) alternatives. It should be noted that the x axis, water recovery, in (ac) is different from that in (df). (ac) show the SEP and EPE performance at various target total water recoveries, while (df) indicate the change of osmotic pressure difference and hydraulic pressure at the real-time water recovery from 0 to the target value (R’).
Figure 2. Specific energy production (SEP) and energy production efficiency (EPE) of (a) co-current single-stage PRO (SS-PRO), (b) counter-current SS-PRO, and (c) atmospheric batch PRO (AB-PRO), and energy components of (d) co-current SS-PRO, (e) counter-current SS-PRO, and (f) AB- PRO. In figure (ac), various salt solutions with different concentrations were used for osmotic energy harvesting including a 0.6 M NaCl solution (e.g., equivalent to the salinity of seawater) and a 1.2 M NaCl solution (e.g., equivalent to the salinity of SWRO brine) as draw solution (DS) alternatives, a 0.01 M NaCl solution (e.g., equivalent to the salinity of river water) and a 0.05 M NaCl solution (e.g., equivalent to the salinity of wastewater concentrate) as feed solution (FS) alternatives. It should be noted that the x axis, water recovery, in (ac) is different from that in (df). (ac) show the SEP and EPE performance at various target total water recoveries, while (df) indicate the change of osmotic pressure difference and hydraulic pressure at the real-time water recovery from 0 to the target value (R’).
Membranes 13 00354 g002
Figure 3. Specific energy production (SEP) and energy production efficiency (EPE) of (a) co-current single-stage PRO (SS-PRO), (b) counter-current SS-PRO, and (c) atmospheric batch PRO (AB-PRO) under ideal and practical baseline conditions. The feed solution (DS) is a 0.05 M NaCl solution, and the draw solution (DS) is a 1.2 M NaCl solution. Practical baseline conditions: η P = 0.85 , η E R D = 0.90 , η P X = 0.98 , and P l o s s = 0.2   bar in AB-PRO and 1.4 bar in SS-PRO.
Figure 3. Specific energy production (SEP) and energy production efficiency (EPE) of (a) co-current single-stage PRO (SS-PRO), (b) counter-current SS-PRO, and (c) atmospheric batch PRO (AB-PRO) under ideal and practical baseline conditions. The feed solution (DS) is a 0.05 M NaCl solution, and the draw solution (DS) is a 1.2 M NaCl solution. Practical baseline conditions: η P = 0.85 , η E R D = 0.90 , η P X = 0.98 , and P l o s s = 0.2   bar in AB-PRO and 1.4 bar in SS-PRO.
Membranes 13 00354 g003
Figure 4. The impact of pump efficiency ( η P ) on specific energy production (SEP) and energy production efficiency (EPE) of (a) co-current single-stage PRO (SS-PRO), (b) counter-current SS-PRO, and (c) atmospheric batch PRO (AB-PRO). The feed solution (DS) is a 0.05 M NaCl solution, and the draw solution (DS) is a 1.2 M NaCl solution. In addition to the different η P as shown in the figure, other conditions are fixed at the baseline level: η P X = 0.98 , η E R D = 0.90 , and P l o s s = 0.2   bar in AB-PRO and 1.4 bar in SS-PRO.
Figure 4. The impact of pump efficiency ( η P ) on specific energy production (SEP) and energy production efficiency (EPE) of (a) co-current single-stage PRO (SS-PRO), (b) counter-current SS-PRO, and (c) atmospheric batch PRO (AB-PRO). The feed solution (DS) is a 0.05 M NaCl solution, and the draw solution (DS) is a 1.2 M NaCl solution. In addition to the different η P as shown in the figure, other conditions are fixed at the baseline level: η P X = 0.98 , η E R D = 0.90 , and P l o s s = 0.2   bar in AB-PRO and 1.4 bar in SS-PRO.
Membranes 13 00354 g004
Figure 5. The impact of energy recovery device efficiency ( η E R D ) on specific energy production (SEP) and energy production efficiency (EPE) of (a) co-current single-stage PRO (SS-PRO), (b) counter-current SS-PRO, and (c) atmospheric batch PRO (AB-PRO). The feed solution (DS) is a 0.05 M NaCl solution, and the draw solution (DS) is 1.2 M NaCl solution. In addition to the different η E R D as shown in the figure, other conditions are fixed at the baseline level: η P = 0.85 , η P X = 0.98 , and P l o s s = 0.2   bar in AB-PRO and 1.4 bar in SS-PRO.
Figure 5. The impact of energy recovery device efficiency ( η E R D ) on specific energy production (SEP) and energy production efficiency (EPE) of (a) co-current single-stage PRO (SS-PRO), (b) counter-current SS-PRO, and (c) atmospheric batch PRO (AB-PRO). The feed solution (DS) is a 0.05 M NaCl solution, and the draw solution (DS) is 1.2 M NaCl solution. In addition to the different η E R D as shown in the figure, other conditions are fixed at the baseline level: η P = 0.85 , η P X = 0.98 , and P l o s s = 0.2   bar in AB-PRO and 1.4 bar in SS-PRO.
Membranes 13 00354 g005
Figure 6. The impact of pressure loss ( P l o s s )   on specific energy production (SEP) and energy production efficiency (EPE) of (a) co-current single-stage PRO (SS-PRO), (b) counter-current SS-PRO, and (c) atmospheric batch PRO (AB-PRO). The feed solution (DS) is a 0.05 M NaCl solution, and the draw solution (DS) is a 1.2 M NaCl solution. In addition to the different P l o s s as shown in the figure, other conditions are fixed at the baseline level: η P = 0.85 , η P X = 0.98 , and η E R D = 0.90 .
Figure 6. The impact of pressure loss ( P l o s s )   on specific energy production (SEP) and energy production efficiency (EPE) of (a) co-current single-stage PRO (SS-PRO), (b) counter-current SS-PRO, and (c) atmospheric batch PRO (AB-PRO). The feed solution (DS) is a 0.05 M NaCl solution, and the draw solution (DS) is a 1.2 M NaCl solution. In addition to the different P l o s s as shown in the figure, other conditions are fixed at the baseline level: η P = 0.85 , η P X = 0.98 , and η E R D = 0.90 .
Membranes 13 00354 g006
Table 1. Values of parameters for the energy production calculation.
Table 1. Values of parameters for the energy production calculation.
Term aIdeal CasePractical CaseReference
BaselineOther Scenarios
η P 10.850.80–0.95[31,33]
η P X 10.980.98[27,34]
η E R D  b10.90HT: 0.85–0.95[35,36,37]
PX: 0.98
P l o s s
( P D , l o s s = P F , l o s s ) c
SS-PRO0 bar1.4 bar0.7–3.5 bar[18,38]
AB-PRO0 bar0.7 bar0.1–0.5 bar
a  η P is the pump efficiency, η P X is the pressure exchanger (PX) efficiency, η E R D is the energy recovery device (ERD) efficiency, P l o s s is the pressure loss. b the configured ERD can be either a hydro-turbine (HT) to generate electricity or a PX to exchange pressure with another system. c the pressure loss on both DS and FS sides is assumed to be the same (i.e., P D , l o s s = P F , l o s s ).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, D.; Mo, Z.; She, Q. Comparison of Energy Efficiency between Atmospheric Batch Pressure-Retarded Osmosis and Single-Stage Pressure-Retarded Osmosis. Membranes 2023, 13, 354. https://doi.org/10.3390/membranes13030354

AMA Style

Li D, Mo Z, She Q. Comparison of Energy Efficiency between Atmospheric Batch Pressure-Retarded Osmosis and Single-Stage Pressure-Retarded Osmosis. Membranes. 2023; 13(3):354. https://doi.org/10.3390/membranes13030354

Chicago/Turabian Style

Li, Dan, Zijing Mo, and Qianhong She. 2023. "Comparison of Energy Efficiency between Atmospheric Batch Pressure-Retarded Osmosis and Single-Stage Pressure-Retarded Osmosis" Membranes 13, no. 3: 354. https://doi.org/10.3390/membranes13030354

APA Style

Li, D., Mo, Z., & She, Q. (2023). Comparison of Energy Efficiency between Atmospheric Batch Pressure-Retarded Osmosis and Single-Stage Pressure-Retarded Osmosis. Membranes, 13(3), 354. https://doi.org/10.3390/membranes13030354

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