Microsoft Word - ETASR_V11_N2_pp6902-6906 Engineering, Technology & Applied Science Research Vol. 11, No. 2, 2021, 6902-6906 6902 www.etasr.com Maddah: Predicting Flux Rates against Pressure via Solution-Diffusion in Reverse Osmosis Membranes Predicting Flux Rates against Pressure via Solution- Diffusion in Reverse Osmosis Membranes Hisham A. Maddah Department of Chemical Engineering King Abdulaziz University Rabigh, Saudi Arabia hmaddah@kau.edu.sa Abstract-This paper suggests a new method of predicting flux values at Reverse Osmosis (RO) desalination plants. The solution- diffusion model is utilized to determine the osmotic pressure drops for seawater sources. The same technique was applied to the groundwater source at the Abqaiq plant (500 RO plant) to calculate the osmotic pressure. The calculated osmotic pressures were utilized to determine the appropriate flux rates and membrane resistances of different BWRO Toray membranes and a performance comparison between various membranes has been established. The model results confirm an inverse relationship between membrane thickness and water flux rate. Also, a proportional linear relation between the overall water flux and the applied pressure is identified. Higher flux rates and lower salinity indicate lower membrane resistance yielding higher production. The modeled data predict that BWRO Toray TM720D-440 with an 8" membrane is the optimal choice for treating waters from the three water sources at the Abqaiq plant. Keywords-reverse osmosis; treatment; desalination; modeling I. INTRODUCTION The solution-diffusion model is a popular expression used to explain the transport in dialysis, reverse osmosis, gas permeation, and pervaporation. Previous experimental data and modeling results verified that the flux rate is proportional to the gradient in the chemical potential [1]. There are two different models that describe and control the permeation in membranes for better separation. The first model is the solution-diffusion model where permeants dissolve (sorption) in the membrane material at the upstream interface in the presence of a concentration gradient that allows permeants to diffuse through the membrane and desorbed on the downstream interface side. The separation between different permeants occurs because each material has a different diffusion rate in the membrane. The solution-diffusion model has been used since 1940 to explain the transport of gases across polymeric membranes. A second model, called the pore-flow model, depends on the presence of a pressure gradient for a convection flow of permeants through the membrane's tiny pores, and is more limited compared to the first model. Exclusion or filtration of larger permeant's pores is the separation technique explained via the pore-flow model [1, 2]. There is a major difference between the solution-diffusion model and the pore-flow model in expressing the chemical potential. In the solution-diffusion model, the pressure within a membrane is uniform and the chemical potential gradient is expressed only as a concentration gradient. Solution-diffusion membranes transmit pressure in the same way as liquids, which is the reason for expressing the pressure difference across the membrane as a concentration gradient only. On the other hand, the chemical potential gradient in the pore-flow model is expressed only as a pressure gradient since the concentrations of both solvent and solute within a membrane are uniform. Comparisons between the two models for a one-component solution in a pressure-driven permeation system were conducted in [1, 2]. The objective of this work is to estimate the osmotic pressure drop value of the high rejection brackish water RO membrane (Toray TM720D-400 with 8") by using the solution- diffusion model that is applied to the Abqaiq plant (500 RO plant) for Shedgum/Abqaiq groundwater at Saudi Aramco, Dhahran, Saudi Arabia. Osmotic pressure drops have been calculated for the groundwater, the Arabian Gulf water, and the Red Sea water at the same plant configuration and operating conditions of the Abqaiq plant in Aramco. The calculated osmotic pressures are utilized to determine the applied pressure drop across the membrane and the applicability of using different BWRO Toray membrane types for the treatment of seawater. The maximum achievable water flux values are determined for the various suggested BWRO membranes for the three water sources. Also, the membrane resistance values have been investigated for comparison purposes. The ideal membrane for the treatment of various water sources at a RO plant with the same configuration of the Abqaiq plant has been selected. The feasibility of using BWRO membranes in the desalination of Red Sea water in Jeddah, Saudi Arabia is studied at the same flux rate of the Arabian Gulf water source and the same plant conditions. Osmotic pressure drop, applied pressure drop, flux rates, and membrane resistance values for the Red Sea water source were compared with those of Shedgum/Abqaiq groundwater and Arabian Gulf water. II. REVERSE OSMOSIS In reverse osmosis, water flows from the salt solution to the pure waterside by applying pressure (∆�) that is greater than the osmotic pressure (∆�) [1]. Generally, in reverse osmosis, the condition ∆� > ∆� must be satisfied all the time to allow water to pass through the membrane and reach the permeate Corresponding author: Hisham A. Maddah Engineering, Technology & Applied Science Research Vol. 11, No. 2, 2021, 6902-6906 6903 www.etasr.com Maddah: Predicting Flux Rates against Pressure via Solution-Diffusion in Reverse Osmosis Membranes side [1, 2]. Reverse osmosis membranes are preferred over ultrafiltration and nanofiltration since they are capable of removing 90 to 99% of TDS in water [3]. The osmotic pressure ∆π is defined as the pressure difference �� − �ℓ across the membrane. If a pressure higher than the osmotic pressure is applied to the feed side of the membrane [1]. This process is called the reverse osmosis. The driving forces in a reverse osmosis membrane according to the solution-diffusion and the pore-flow models are visualized in [1]. �� and � are the chemical potential and activity coefficient, respectively, of component i [1]. III. METHOD AND EQUATIONS Collected Abqaiq 500 RO plant data (Table I) have been used to determine osmotic pressure drop values for the RO membrane (Toray TM720D-400 with 8") from (1) and (2). However, in order to calculate the osmotic pressure for seawater sources, the same information of Shedgum/Abqaiq groundwater at Abqaiq 500 RO plant was applied, except for the flux and salinity values, for the treatment of either the Arabian Gulf or the Red Sea waters as listed in Table I [1, 4]. Water permeability is approximately determined to be 9.5×10 -7 cm 2 /s [8]. For water-salt solution, reverse osmosis permeation expression can be simplified as [1, 5]: �� = �(�� − ��) (1) � = ���������ℓ (2) where �� is the membrane flux of component �, water (gfd), Δ� is the applied pressure drop across the membrane (psi), Δ� is the osmotic pressure drop across the membrane (psi), A is the water permeability constant (cm/atm×s), �� is the permeability of component � , water (cm2/s), ��� is the initial mole concentration of water (ppm), �� is the water molar volume (cm 3 /mol), T is the w ater temperature (K), R is the gas constant, and ℓ is the membrane thickness which is assumed to be similar to spacer thickness (mil). Membrane resistance [8] constants for each BWRO Toray membrane have been calculated by using (3). �� = ∆%& ℛ( (3) where �� is the membrane flux of component �, water (gfd), Δ� is the applied pressure across the membrane (psi), ) is the dynamic viscosity of water (lb s/ft 2 ), and ℛ* is the membrane resistance (1/ft). π = ℳ-. (4) In (4), ℳ is the molar concentration of dissolved species (mol/L), R is the ideal gas constant (0.08206L.atm/mol.K), and T is the water temperature (K). Equation (5) defines the ability of a membrane to separate salt from the feed solution which is known as membrane removal percentage (/) and it increases with the applied pressure. The feed TDS concentration is taken from the three studied sources, as shown in Table I, while the outlet TDS concentration is determined by using (5) at a similar removal percentage of Toray TM720D-400 with 8’’ membrane that is 99.8% (Table IV). The water molecular weight (18g/mol) should be used to convert the ppm values to molar concentrations of TDS. / = 0�1�2�1ℓ�1� 3 × 100 (5) where / is the membrane removal percentage, �7� is the initial concentration of component 8, salt (ppm), and �7ℓ is the final concentration of component 8 (ppm). Table II shows the applied pressure drop must be at 20psi or below per element (RO module) and 60psi or below per vessel [4, 6]. The assumption of having equal pressure on membranes per vessel would simplify our calculations. Altaee’s study showed that permeate flow, pressure and recovery rate are distributed almost equally to membranes per RO vessel [10]. A field study confirmed an improved performance by rearranging the elements in pressure vessels in order to reduce the pressure drop and permeate conductivity across the vessel [11]. Typical flux rates and maximum recovery values for the groundwater and the two studied water source scenarios (the Arabian Gulf and the Red Sea waters) at Abqaiq 500 RO plant are given in Table III. TABLE I. DATA OF RO MEMBRANE PROCESS AT ABQAIQ 500 RO PLANT AND THE TWO SEAWATER STUDIED SCENARIOS [1, 4, 6, 7] Parameter Shedgum/Abqaiq groundwater Arabian Gulf water Red Sea water Membrane type Toray TM720D-400 with 8" RO module 72 parallel membranes × 8 units Membrane thickness (ℓ) [3] Assumed to be similar to spacer thickness of 34 ;�< Membrane area (�=>?) [3] 400ft2 Max pressure drop per vessel (∆�) ~60psi Max pressure drop per membrane (∆�) ~20psi Water salinity (���)* ~2800 [4] ~41070 [6] ~42070 [7] Membrane water flux (�� )* ~18gfd ~12 gfd ~12 gfd Water temperature (.) 25C Water permeability constant (��)** 9.5×10-7cm2/s Water molar volume (��) 18cm3/mol Gas constant (-) 8.2057×10-5m3.atm/mol.K *Averaged values ** Taken from [2], regardless of the temperature effect on permeability. Can be calculated at different temperatures [16] Engineering, Technology & Applied Science Research Vol. 11, No. 2, 2021, 6902-6906 6904 www.etasr.com Maddah: Predicting Flux Rates against Pressure via Solution-Diffusion in Reverse Osmosis Membranes TABLE II. OPERATING DESIGN LIMITS OF THE OVERALL RO MODULE AT ABQAIQ 500 RO PLANT [4, 12, 13] Operating limits Maximum operating pressure 600psi (4.1MPa) Maximum feed water temperature 113°F (45°C) Maximum feed water SDI15 5 Feed water chlorine concentration Not detectable Feed water pH range, continuous operation 2-11 Feed water pH range, chemical cleaning 1-12 Maximum pressure drop per element 20psi (0.14MPa) Maximum pressure drop per vessel 60psi (0.4MPa) TABLE III. CHARACTERISTICS OF GROUNDWATER SOURCE AND STUDIED WATER SOURCES AT THE ABQAIQ 500 RO PLANT [4] Water source Shedgum/Abqaiq groundwater Arabian Gulf Red Sea Feed silt density index @AB < 3 @AB < 4 @AB < 4 Typical target flux (gfd) 18 12 12 Max. element recovery (%) 19 14 14 The determined osmotic pressure values for the RO membrane (Toray TM720D-400 with 8") of the groundwater and the two studied water sources are used again in (1) to calculate the applied pressure drop and suggested flux values. The same osmotic pressure drop for each case is utilized to determine the results of different Toray BWRO membrane types at high, low, and standard operating pressure as shown in Table IV. TM-720-370 and TM720-440 are standard BWRO membranes and TM720C-440, TM720L-400, and TM720L- 440 are low-pressure BWRO membranes whereas TM720DA400, TM720D-400, and TM720D-440 are high- pressure BWRO membranes. It is worth mentioning that our applied pressure drop must be higher than the calculated osmotic pressure in order to have a positive flux. TS-diagrams [7] are used to determine the exact value of water densities at different feed sources from the average water temperature and water salinity (Table V). The exact water densities allow us to convert gas constant values from m 3 .atm/mol.K to kg.atm/mol.K to progress calculations. TABLE IV. TORAY BRACKISH WATER RO 8’’ DIAMETER MEMBRANES [13, 14] Category Type Rejection (%) Thickness (mil)* Standard BWRO TM720-370 99.7 31 TM720-440 99.7 28 High- pressure BWRO TM720DA400 99.8 31 TM720D-400 99.8 34 TM720D-440 99.8 28 Low-pressure BWRO TM720C-440 99.2 28 TM720L-400 99.5 31 TM720L-440 99.5 28 * The membrane thickness is assumed to be the same as spacer thickness TABLE V. WATER DENSITIES FROM TS-DIAGRAMS [4, 6, 7, 15] Water Source Temperature (°C) Salinity (ppm) Density (kg/m 3 ) Shedgum/Abqaiq groundwater 25 2800 [2] 999.19 Arabian Gulf 25 41070 [4] 1027.97 Red Sea 25 42070 [6] 1028.67 IV. RESULTS AND DISCUSSION Equations (1) and (2) allow us to calculate the osmotic pressure drop for each water source (Table VI). The osmotic pressure of the groundwater source is less than the Arabian Gulf and the Red Sea water sources which is related to the flux rates and water salinity. Flux rates for the Arabian Gulf and the Red Sea waters are approximately the half of the groundwater source. However, the water salinity of the groundwater source is much lower than the other sources. Therefore, the required applied pressure drop must be larger in the case of seawater sources due to their higher determined osmotic pressure values. Since the plant configuration has 8 elements per vessel, we should have a maximum osmotic pressure of 60psi or less per vessel which is equivalent to a max pressure of 7.5psi per membrane, assuming that the pressure is distributed equally on membranes per vessel. The selected applied pressure range for our study is 6.5 to 7.5psi. Maximum pressure values are assigned to the different membranes based on their category as illustrated in Table VII. TABLE VI. CALCULATED OSMOTIC PRESSURE DROP (∆�) FOR EACH WATER SOURCE Water source A (cm/atm.s) Ji (cm/s) Ji/A (atm) ∆π (atm) ∆π (psi) ∆π per vessel < 60 (psi) Shedgum/ Abqaiq groundwater 0.00808 0.00083 0.1028 0.44 6.48 51.84 Arabian Gulf 0.00755 0.00056 0.0742 0.47 6.90 55.21 Red Sea 0.00754 0.00056 0.0743 0.47 6.90 55.20 TABLE VII. ASSIGNED PRESSURE VALUES FOR TORAY BWRO MEMBRANES Category Type ∆D Eange (DFG)* Standard BWRO TM720-370 6.50 - 7.25 TM720-440 6.50 - 7.25 High-pressure BWRO TM720DA400 6.50 - 7.50 TM720D-400 6.50 - 7.50 TM720D-440 6.50 - 7.50 Low-pressure BWRO TM720C-440 6.50 - 7.00 TM720L-400 6.50 - 7.00 TM720L-440 6.50 - 7.00 * High and low-pressure values are taken relative to the standard pressure The relationship between the applied pressure drops and the overall water flux rates for the groundwater source are obtained in Figure 1(a)-(c) for standard, high-pressure, and low-pressure Toray BWRO membranes. Figure 1 shows that the maximum possible flux for the groundwater in the standard membranes is around 11gfd for TM720-440 membrane. In Figures 3(b) and 3(c) the highest observed groundwater flux in the high-pressure and low-pressure membranes are 14.7gfd for TM720D-440 and 7.5gfd for TM720C-440 and TM720L-440 respectively (blue and green lines overlap). This observation is associated with the membrane thickness in which the least membrane thickness (28mils) has been capable to achieve the highest flux. This confirms an inverse relationship between the membrane thickness and the water flux rate. Further, there is a linear relationship between the applied pressure drop and the overall water flux. Engineering, Technology & Applied Science Research Vol. 11, No. 2, 2021, 6902-6906 6905 www.etasr.com Maddah: Predicting Flux Rates against Pressure via Solution-Diffusion in Reverse Osmosis Membranes Fig. 1. Effect of different applied pressures on the groundwater flux for (a) Toray standard BWRO membranes, (b) Toray high-pressure BWRO membranes, and (c) Toray low-pressure BWRO membranes. Figure 2 identifies a proportional relationship between the water flux and the applied pressure across the membrane. The highest recorded flux is accounted for TM720D-440 for Shedgum/Abqaiq groundwater because water TDS is low for groundwater and TM720D-440 has the lowest thickness and the highest pressure range. The Arabian Gulf and the Red Sea water sources almost have similar flux rates at the same applied pressures due to the similarities in their water salinity levels. TM720C-440, TM720L-400, and TM720L-440 membranes reserved the lowest flux values since they are categorized as low-pressure BWRO membranes. Fig. 2. Observed water flux for various water sources at different applied pressures using Toray BWRO membranes. Figure 3 demonstrates the membrane resistance for the three studied water sources. Seawater sources have higher membrane resistances than the groundwater source because of their lower flux and higher TDS. TM720L-400 has the highest membrane resistance since it is in the low-pressure category and has the highest membrane thickness of 31mils. Equation (4) calculations are shown in Table VIII. The study predictions estimated that the overall osmotic pressure drops required for seawater and groundwater treatment plants are approximately 55psi and 830psi respectively. The higher the salinity difference between the fed and the produced water, the more the osmotic pressure drop we need to overcome in order to produce treated water (positive flux). Fig. 3. Observed membrane resistance of various water sources in Toray BWRO membranes. TABLE VIII. VAN'T HOFF CALCULATIONS FOR THE REQUIRED OSMOTIC PRESSURES Water source Concentration (mol/L) Membrane removal (%) Osmotic pressure (atm) |∆I| (psi) TDSin TDSout πin πout |∆I| Shedgum/ Abqaiq groundwater 0.156 0.00031 99.8 3.81 0.01 3.80 55.80 Arabian Gulf 2.282 0.005 99.8 55.82 0.11 55.71 818.41 Red Sea 2.337 0.005 99.8 57.18 0.11 57.07 838.34 V. CONCLUSION The application of the solution-diffusion model to the Abqaiq plant (500 RO plant) is initiated by using various parameters to calculate the osmotic pressure of Toray TM720D-400 with 8" membrane for Shedgum/Abqaiq groundwater treatment. For the same membrane, the osmotic pressure values are determined for the Arabian Gulf and the Red Sea waters to predict flux rates in other membranes for seawater situations. Low, standard, and high pressure BWRO Toray membranes performances have been compared to identify the optimal membrane for treating saline water from the three studied water sources at the Abqaiq 500 RO plant. Engineering, Technology & Applied Science Research Vol. 11, No. 2, 2021, 6902-6906 6906 www.etasr.com Maddah: Predicting Flux Rates against Pressure via Solution-Diffusion in Reverse Osmosis Membranes The assumption of having a membrane thickness that is similar to its spacer thickness may not seem very accurate. However, it is true that we should have a proportional relation between both thicknesses which suggests that our results are still valid. A linear relationship has been observed between the water flux and the applied pressure drops. It is proved that membrane flux decreases with the increase in membrane thickness at constant pressure drop. Modeling results endorse that BWRO Toray TM720D-440 with 8"membrane is the optimum membrane choice for the water treatment from the three water sources at Abqaiq 500 RO plant since it has the lowest membrane resistance and the highest overall water flux. REFERENCES [1] J. G. Wijmans and R. W. Baker, "The solution-diffusion model: a review," Journal of Membrane Science, vol. 107, no. 1, pp. 1–21, Nov. 1995, https://doi.org/10.1016/0376-7388(95)00102-I. [2] D. R. Paul, "Reformulation of the solution-diffusion theory of reverse osmosis," Journal of Membrane Science, vol. 241, no. 2, pp. 371–386, Oct. 2004, https://doi.org/10.1016/j.memsci.2004.05.026. [3] H. A. Maddah and A. M. Chogle, "Applicability of low pressure membranes for wastewater treatment with cost study analyses," Membrane Water Treatment, vol. 6, no. 6, pp. 477–488, Nov. 2015, http://dx.doi.org/10.12989/mwt.2015.6.6.477. [4] Abqaiq plants: 500 RO plant, Operator Handbook. Dhahran, Saudi Arabia: Saudi Aramco Company: Training & Career Development Center, 2008. [5] A. Sagle and B. Freeman, "Fundamentals of membranes for water treatment," The future of desalination in Texas, vol. 2, no. 363, 2004, Art. no. 137. [6] R. Smith, A. Purnama, and H. H. Al-Barwani, "Sensitivity of hypersaline Arabian Gulf to seawater desalination plants," Applied Mathematical Modelling, vol. 31, no. 10, pp. 2347–2354, Oct. 2007, https://doi.org/10.1016/j.apm.2006.09.010. [7] E. A. Abdel-Aal, M. E. Farid, F. S. M. Hassan, and A. E. Mohamed, "Desalination of Red Sea water using both electrodialysis and reverse osmosis as complementary methods," Egyptian Journal of Petroleum, vol. 24, no. 1, pp. 71–75, Mar. 2015, https://doi.org/10.1016/j.ejpe. 2015.02.007. [8] J. C. Crittenden, R. R. Trussell, D. W. Hand, K. J. Howe, and G. Tchobanoglous, MWH’s Water Treatment: Principles and Design, 3rd ed. Hoboken, NJ, USA: John Wiley & Sons, 2012. [9] J. Prausnitz, R. Lichtenthaler, and E. G. de Azevedo, Molecular Thermodynamics of Fluid-Phase Equilibria, 3rd edition. Upper Saddle River, NJ, USA: Pearson, 1998. [10] A. Altaee, "A Computational Model to Estimate the Performance of 8 inches RO Membranes in Pressure Vessel," Journal of Membrane and Separation Technology, vol. 1, no. 1, pp. 60–71, Oct. 2012. [11] S. S. Madaeni, M. Afshar, N. Jaafarzadeh, F. Tarkian, and K. Ghasemipanah, "Rearrangement of membrane elements in the pressure vessels for optimum utilization of reverse osmosis process," Chemical Engineering Research and Design, vol. 89, no. 1, pp. 48–54, Jan. 2011, https://doi.org/10.1016/j.cherd.2010.04.021. [12] "TorayTM720D-400 with 8"," in ROMEMBRA® Toray RO Membrane, Tokyo, Japan: Toray Corporation, 2011. [13] "Brackish Water & Sea Water Membranes, TM (inner connecter type) Datasheets," in ROMEMBRA® Toray RO Membrane, Tokyo, Japan: Toray Corporation, 2016. [14] Maintenance and handling manual for membrane elements (rev. 105, 106). Poway, CA, USA: Toray Membrane USA, 2016. [15] T. Beer, Environmental Oceanography, 2nd ed. Boca Raton, FL, USA: CRC Press, 1996. [16] H. A. Maddah, "Application of Finite Fourier Transform and Similarity Approach in a Binary System of the Diffusion of Water in a Polymer," Journal of Materials Science and Chemical Engineering, vol. 4, no. 4, pp. 20–30, Apr. 2016, https://doi.org/10.4236/msce.2016.44003. AUTHORS PROFILE Hisham Maddah is an independent researcher and a faculty member working with the Chemical Engineering Department at King Abdulaziz University in Rabigh (KAU-Rabigh). He has been awarded a full scholarship and received his PhD in "Naturally-Sensitized Photoanodes for Molecular Photovoltaics" from the University of Illinois at Chicago (UIC) by Summer 2020. He has completed his MS and BS in Chemical Engineering in the University of Southern California (USC) in 2017 and KAU in 2012 respectively. He is an expert in organic-inorganic dye-sensitized solar cells for solar energy harvesting. His research interests include mini-passive solar stills, machine learning for solar-desalination systems, membrane separation, capacitive deionization, and activated carbons for water treatment, separation technologies, and higher education administration. He has published more than 50 articles and 3 book chapters in reputed international peer-reviewed journals and conferences and co-arranged 5 outreach activities on third- generation solar cells.