Next Article in Journal
Surfactant-Enhanced Extraction of Lutein from Marigold Petals using an Aqueous Two-Phase System
Next Article in Special Issue
Modeling, Simulation, and Optimization of Membrane Processes
Previous Article in Journal
Potential of Advanced Oxidation as Pretreatment for Microplastics Biodegradation
Previous Article in Special Issue
A Study of Copper (II) Ions Removal by Reverse Osmosis under Various Operating Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multiscale Analysis of Permeable and Impermeable Wall Models for Seawater Reverse Osmosis Desalination

1
School of Materials Science and Engineering, Sun Yat-sen University, Guangzhou 510006, China
2
School of Computer Science and Engineering, Sun Yat-sen University, Guangzhou 510006, China
3
National Supercomputing Center in Guangzhou (NSCC-GZ), Guangzhou 510006, China
4
Guangdong Province Key Laboratory of Computational Science, Guangzhou 510006, China
*
Author to whom correspondence should be addressed.
Separations 2023, 10(2), 134; https://doi.org/10.3390/separations10020134
Submission received: 23 January 2023 / Revised: 7 February 2023 / Accepted: 10 February 2023 / Published: 15 February 2023
(This article belongs to the Special Issue Modeling, Simulation, and Optimization of Membrane Processes)

Abstract

:
In recent years, high permeability membranes (HPMs) have attracted wide attention in seawater reverse osmosis (SWRO) desalination. However, the limitation of hydrodynamics and mass transfer characteristics for conventional spiral wound modules defeats the advantage of HPMs. Feed spacer design is one of the effective ways to improve module performance by enhancing permeation flux and mitigating membrane fouling. Herein, we propose a multiscale modeling framework that integrates a three-dimensional multi-physics model with a permeable wall and an impermeable wall, respectively, at a sub-millimeter scale and a system-level model at a meter scale. Using the proposed solution framework, a thorough quantitative analysis at different scales is conducted and it indicates that the average errors of the friction coefficient and the Sherwood number using the impermeable wall model are less than 2% and 9%, respectively, for commercial SWRO membrane (water permeability 1 L m−2 h−1 bar−1) and HPMs (3 L m−2 h−1 bar−1, 5 L m−2 h−1 bar−1 and 10 L m−2 h−1 bar−1) systems, compared to the predictions using the permeable wall model. Using both the permeable and impermeable wall models, the system-level simulations, e.g., specific energy consumption, average permeation flux, and the maximum concentration polarization factor at the system inlet are basically the same (error < 2%), while the impermeable wall model has a significant advantage in computational efficiency. The multiscale framework coupling the impermeable wall model can be used to guide the efficient and accurate optimal spacer design and system design for HPMs using, e.g., a machine learning approach.

1. Introduction

Membrane-based separation technologies, such as micro- [1,2], ultra- [3,4], or nanofiltration (NF) [5,6], reverse osmosis (RO) [7,8], electrodialysis [9,10] and membrane distillation [11,12], have increasingly attracted extensive attention in industry and academia. Due to its energy-saving and cost-competitiveness, RO is one of the most advanced membrane-based technologies and has widespread applications in desalination [13,14], wastewater treatment [15,16], pharmaceutical production [17,18], electronics industries [19,20] etc. Scaling and fouling are common challenging problems in RO, especially in the condition of high flux (e.g., more than 50 L m−2 h−1 [21]) with the use of a high permeability membrane (HPM). Spacer design in spiral wound membrane (SWM) modules has wide application prospects to mitigate membrane scaling and fouling and improve the performance of a SWM module [22].
Lin et al. [23] summarized the research development on feed spacer in the last twenty years and emphasized that it is critical for feed spacer design to balance the feed channel pressure drop and mass transfer on the membrane surfaces. Among these works, the computational fluid dynamics (CFD) simulation technology was widely applied to capture the local fluid flow and transport phenomena in membrane-based processes. The CFD models in a spacer-filled channel are divided into a permeable wall model and an impermeable (or impermeable-dissolving) wall model [24], according to the type of boundary conditions enforced on the membrane walls. For the permeable wall model, Anqi et al. [25] adopted the three-dimensional (3D) CFD technology to estimate the effect of the angle (α = 60°, 90°, 120°) between the spacer filaments and the Reynolds numbers (100, 400, 800). The results indicate that increasing crossflow contributes to enhance mass transfer, and reducing the angle α will decrease the friction loss. The experimental and simulation results by Lin et al. [26] showed that the feed spacer geometries with non-uniform filament have a significant effect on both anti-fouling performance and hydraulics. The study by Foo et al. [27] showed that the mass transfer in the spacer-filled channel can be markedly enhanced by an unsteady forced-slip combined with a spacer design. In this case, the optimization of the geometrical parameters can greatly affect the resonant frequency and thereby reduce the concentration polarization (CP). The study by Singh et al. [28] showed that the feed spacer configurations with the higher water permeation are accompanied with a higher flow resistance. Li et al. [29] proposed a hybrid model that couples a CFD model at a small-scale and a system level model for brackish water RO (BWRO) desalination. The simulated results fit well with the measurements from industry. Lin et al. [30] designed a novel membrane module with a diagonal-flow feed channel. The CFD simulations and experimental results indicate that the novel membrane module has better performance (e.g., higher permeation flux and higher salt rejection) compared to the conventional membrane module. This is because the former one has a higher average cross-flow in the feed channel, which can reduce CP and mitigate membrane scaling. Apart from the RO process, spacer design also has widespread applications in other membrane-based processes, such membrane distillation [31].
However, the typical used permeable wall model is fully coupled and is computationally expensive. Because the cross velocity is several orders of magnitude higher than the permeation flux, the boundary condition on the membrane surface barely affects the fluid flow and solute transport. Some researchers used impermeable wall boundary conditions to obtain reasonable simulation results in their CFD simulations [32,33,34,35,36,37]. The impermeable wall model is a one-way coupling that greatly reduces the nonlinear effect compared with the commonly used permeable wall model. Geraldes and Afonso [38] proposed a generalized correction factor to correct the mass transfer coefficients at impermeable walls with an average error of 3.6% compared with the CFD predictions with a permeable wall, which is independent of the flow regime and feed spacer geometry. The predicted results using the impermeable wall model are in good agreement with the experimental data [32,33] with an error of approximately 10%. Toh et al. [37] presented a multi-scale model that couples an impermeable wall model at a small-scale and a large-scale model for RO desalination. The simulated results for BWRO and SWRO systems with the use of an advanced feed spacer and HPMs showed that the significance of improving the spacer performance is even greater than increasing the membrane permeability in terms of cost-effectiveness. Luo et al. [39] proposed a heuristic optimization framework based on a field synergy principle and a high throughput computing strategy for BWRO systems. Using the optimized feed spacer, the axial pressure drop for the entire system is negligibly small compared to a commercial feed spacer. Chen et al. [40] further developed several novel feed spacers based on triply periodic minimal surfaces, and the average water permeation flux can be enhanced by about 8.6% in comparison to a conventional feed spacer. Gu et al. [41] used a multilayer artificial neural network to establish a surrogate pressure drop model for spacer-filled channels in BWRO desalination that enables a quantitative description between the axial pressure drop and design parameters, including feed spacer geometry and inlet velocity magnitude. The surrogate model can be guided to the optimal design of feed spacer. Recently, Luo et al. [42] further proposed a supercomputing-based machine learning-driven multiscale optimization design framework for SWRO desalination systems with HPMs. Using the optimization framework, the specific energy consumption (SEC) and the required membrane area for a designed two-stage SWRO system reduced by 27.5% (to 1.66 kWh m−3) and 37.2%, respectively, under typical conditions (feed salinity 35,000 ppm, recovery rate 50%). By spacer optimization and system design, the CP at the system inlet could be suppressed to no more than 1.20.
Although both the permeable and impermeable wall models have been widely applied in membrane-based processes, the quantitative analyses at different scales for both models were rarely reported, especially in the conditions of high flux with the use of HPMs. In this paper, we establish a multiscale modeling framework that integrates the permeable and impermeable wall models (Section 2.1), respectively, with a system-level model (Section 2.3) by using relations of the friction factor and the Sherwood number with respect to various Reynolds numbers (Section 2.2). A thorough quantitative analysis of the framework is conducted for CFD simulations at a sub-millimeter scale (Section 3.1) and system-level simulations at a meter scale (Section 3.2).

2. Mathematical Modeling

In spacer-filled RO channels, the “periodic fully-developed” transport phenomena were confirmed in previous work [29,43] using CFD techniques, which makes it possible to estimate the cell average mass transfer coefficient at different locations along the axis of the SWM module by changing the average inlet velocity in the CFD model. Here, we present a multiscale modeling framework that integrates the permeable wall model and the impermeable wall model, respectively, as shown in Figure 1.

2.1. The Permeable and Impermeable Wall Models at a Sub-Millimeter Scale

2.1.1. Problem Description

Both the permeable and impermeable wall models in the SWRO desalination consider fluid flow and solute transport, which can be mathematically formulated below [42].
{ ρ ( u · ) u = · [ P I + μ ( u + ( u ) T ) ] , in Ω , ρ · u = 0 , in Ω , D s 2 c = u · c , in Ω , u in = u out , ( u ) in = ( u ) out , u ¯ = U , on Γ I , ( P ) in = ( P ) out , c = c 0 , on Γ I , P = 0 , n · D s c = 0 , on Γ O , J W = w = G 1 ( L P , P , c ) , J S = G 2 ( B , J W , c ) , on Γ U Γ B , u L = u R , P L = P R , c L = c R , on Γ L Γ R ,
where the governing equations, the Navier–Stokes (N–S) equation, the continuity equation under a laminar flow condition and the convection–diffusion equation are defined in a 3D computational domain ( Ω ), or a narrow spacer-filled channel (Figure 2). The variables consist of the velocity vector u ( u , v , w ) , hydraulic pressure (P) and molar concentration (c). The density ( ρ ) and viscosity ( μ ) of fluid and the diffusivity ( D s ) of salt are assumed as constants [42]. L P and B denote water permeability and salt permeability, respectively.

2.1.2. Boundary and Initial Conditions

For the permeable wall model, the permeation flux of water, or water flux ( J W ), is proportional to the difference ( Δ P Δ π ) between the transmembrane pressure ( Δ P = Δ P 0 Δ P c ) and transmembrane osmotic pressure ( Δ π ). The specific value is L P . Δ π is the function of the local concentration on the membrane surfaces. For the impermeable wall model, the permeation flux of water is assumed as zero, and a constant concentration boundary condition is applied on the membrane surfaces. Thus, J W can be expressed as
J W = G 1 ( L P , P , c ) = { L P [ Δ P 0 Δ P c Δ π ( c ) ] , f o r   permeable   wall   modeo 0 , f o r   impermeable   wall   modeo
where Δ P 0 is the inlet transmembrane pressure. The pressure drop ( Δ P c ) along the feed direction (x-direction, see Figure 2b) is the difference of the average pressures at the inlet ( Γ I ) and outlet ( Γ O ), respectively. In the convection–diffusion equation, the boundary conditions for a permeable flux of salt ( J S ) and a constant concentration ( c w , 0 ) are enforced on the membrane walls ( Γ U Γ B ) for the permeable wall model and impermeable wall model, respectively, as below.
{ G 2 ( B , J W , c ) = n ( D s c + c u ) = J W c w ( 1 R salt ) , f o r   permeable   wall   model c w = c w , 0 . f o r   impermeable   wall   model
c w is the concentration on the membrane walls. The salt rejection ( R salt ) is defined by
R salt = J W J W + B .
A fully developed boundary condition with an average velocity magnitude and a constant concentration boundary condition ( c = c 0 ) are applied at the inlet boundary ( Γ I ) for both the permeable and impermeable wall models. The periodic boundary conditions are applied at the lateral boundaries ( Γ L Γ R ), which are commonly used in numerical simulation involving structural repeating units (see Figure 2a). For the impermeable wall model, the N–S equation and continuity equation are solved first. The initial velocity and pressure conditions are assumed as constants. The obtained velocity and pressure profiles are used to further calculate the concentration profile by solving the convection– diffusion equation. For the permeable wall model, the calculated velocity, pressure and concentration profiles are solved simultaneously. The simulated results for Case 1 (Commercial SWRO membrane: L P = 1 L m−2 h−1 bar−1) are considered as the appropriate initial results to accelerate the solutions for Cases 2–4 (HPMs: L P = 3 L m−2 h−1 bar−1; 5 L m−2 h−1 bar−1; 10 L m−2 h−1 bar−1).

2.2. The Relations Coupling the CFD Model and System-Level Model

After solving Equations (1)–(4), the 3D velocity, pressure and concentration profiles in the spacer-filled channel ( Ω ) can be estimated for both the permeable and impermeable wall models. The relations of the friction factor ( f ), Sherwood number ( S h ) and Reynolds numbers (Re = D H u ¯ ρ μ ) are obtained by the fitting curves (f-Re, Sh-Re), namely,
{ f per = k 1 Re t 1 , f o r   permeable   wall   model S h per = k 2 Re t 2 , f o r   permeable   wall   model
and
{ f imp = k 1 Re t 1 , f o r   impermeable   wall   model S h imp = k 2 Re t 2 . f o r   impermeable   wall   model
The friction factor ( f i ) and Sherwood number ( S h i ) are calculated by
f i = Δ P c ,   i L D H / ( 1 2 u ¯ 2 ) ,               i = per ,   imp
and
S h i = k ¯ m ,   i D H D S ,                                             i = per ,   imp
respectively. For a given spacer geometry, the hydraulic diameter ( D H ) can be identified [29]. Then the pressure drop per unit length along the feed direction (x-direction, see Figure 2b can be calculated by
Δ P c ,   i L = P ¯ ( x = 0 ) P ¯ ( x = L ) L ,                                       i = per ,   imp
P ¯ ( x = L ) and P ¯ ( x = 0 ) denote the average pressure at the inlet ( x = 0 ) and outlet ( x = L ). The cell average mass transfer coefficient ( k ¯ m ,   i ) is calculated by
k ¯ m ,   i = 0 L d x 0 W ( D s c r c w c z | z = H / 2 ) d y 0 L d x 0 W d y ,                         i = per ,   imp
L , W and H denote the length (x-direction), width (y-direction) and height (z-direction) of the computational domain (Figure 2b). It should be pointed out that the estimated Sherwood number using the impermeable wall model ( k ¯ m ,   imp ) needs to be corrected using the reported relations [38]
k ¯ m ,   per = k ¯ m ,   imp [ ψ + ( 1 + 0.26 ψ 1.4 ) 1.7 ] ,           ( ψ < 20 )
ψ = J W k ¯ m ,   imp ,
in which J W is calculated using the system-level model, Equations (13)–(15).

2.3. The System-Level Model for SWRO Desalination at a Meter Scale

Furthermore, we established a system-level model [36,37] that enables us to quantify the effect of spacer geometry on the entire performance of RO desalination.
{ d Q d X = J W A tot X = 0 ,   Q = Q 0 , d ( Δ P ) d X = ρ u ¯ 2 2 D H f ( n mem l x ) X = 0 ,   Δ P = Δ P 0 , d w b d X = J W A tot Q ( w b w p ) X = 0 ,   w b = w b ,   0 , J W = L P ( Δ P σ φ R salt w w ) Q = N PV n sp H l y ε u ¯ ,
where Q , Δ P and J W are the flow rate, transmembrane pressure and water flux, respectively. N PV , n mem and n sp denote the number of pressure vessels, the number of membrane module per pressure vessel and the number of feed spacers per membrane module, respectively. l x and l y are the module length parallel to flow and the module length perpendicular to flow, respectively. The porosity, ε is calculated by the formula in previous work [29]. In this paper, the average inlet velocity ( U 0 ) is considered as a constant. w b , w p and w w denote bulk salinity, permeate salinity and salinity on the membrane walls, respectively. These variables change with respect to various dimensionless length, X = x n mem l x . The total membrane area, A tot , is calculated by A tot =   N PV n mem A 0 ( A 0 : membrane area per membrane module). σ and φ denote the reflection coefficient and osmotic pressure coefficient, respectively. The salinity on the membrane walls ( w w ) is calculated by
w w = w p 1 R salt ,
in which w p can be expressed by
w p = w b / [ exp ( ln J W B J W k ¯ m ,   per ) + 1 ] .
A more detailed derivation for Equation (15) can be found in our latest work [42].

3. Results and Discussion

3.1. CFD Simulations

In this paper, the detailed geometry parameters of the feed spacer are listed in Table 1 that are referred to previous literature [44]. The spacer thickness is about 34 mil (863 μm). The number of spatial discrete elements is 1,477,393 (Figure 3), which has comparable precision compared with our previous work [39]. The detailed mesh independence analysis was reported there. Based on the finite element method, Comsol Multiphysics was used to establish and solve the CFD models in this paper. More detailed model parameters can be found in our latest work [42]. The calculated velocity magnitude distributions using the permeable wall model (Figure 4a) and impermeable wall model (Figure 4b) demonstrate that the effect of the permeation flux of water (or water flux) on the cross velocity in the spacer-filled channel is negligible. This is because the water flux (≈10−5 m/s, see Figure 5a) is much smaller than the cross velocity in the feed channel (Figure 4a), which is consistent with the reported literature [29]. Even for the HPMs systems (e.g., L P = 5 L m−2 h−1 bar−1, Figure 5b), this conclusion still holds true. Similarly, the calculated pressure magnitude distribution is basically the same using the permeable wall model (Figure 6a) and impermeable wall model (Figure 6b).
Using the permeable and impermeable wall models, the calculated friction coefficient ( f ) and Sherwood number ( S h ) are shown in Figure 7a–d and Figure 7e–h, respectively, with the use of a commercial SWRO membrane (Case 1: L p = 1 L m−2 h−1 bar−1) and HPMs (Case 2: L p = 3 L m−2 h−1 bar−1; Case 3: L p = 5 L m−2 h−1 bar−1); Case 4: L p = 10 L m−2 h−1 bar−1). It should be noted that the calculated Sherwood number using the impermeable wall model in Figure 7e–h is corrected using Equations (11)–(15). The average errors of the calculated friction coefficient using the impermeable wall model are 0.12%, 0.44%, 0.69% and 1.17%, respectively, for Cases 1–4 compared with that of using the permeable wall model. Correspondingly, the average errors of the estimated Sherwood number are 6.47%, 7.96%, 7.39% and 8.46%, which basically coincide with the reported results in previous work [38]. In an RO system, it is more likely that fouling will occur in the lead elements [45], especially in the condition of high flux. Therefore, we further evaluated the Sherwood number using the impermeable wall model at the system inlet with the errors of 4.25%, 0.27%, 3.49% and 9.51%, respectively. Overall, the accuracy of the impermeable wall model is comparable and acceptable on the prediction of the friction coefficient and Sherwood number versus the permeable wall model. In computational time, the latter one (77.8 h) for each case is about three times than that of the former one (25.3 h). It should be emphasized that the impermeable wall model is only associated with the spacer geometry and inlet velocity, and is independent of other parameters, such as water permeability and operating pressure. Therefore, the Sherwood number with respect to various water permeabilities (Cases 1–4) can be directly calculated using Equations (11)–(15) based on the obtained relations ( f imp Re , S h imp Re ). However, the permeable wall model must be recalculated to obtain the Sherwood number when the water permeability or the other conditions change. Thus, the impermeable wall model has a significant advantage in computational efficiency in comparison to the permeable wall model, which can greatly reduce the computational cost, and can be used for the optimal design of the feed spacer.

3.2. Performance Evaluations at a System-Level

Using the system-level model, Equations (13)–(15) coupled with the permeable wall model, Equations (1)–(5), (7)–(10) and the impermeable wall model, Equations (1)–(4), (6)–(12) respectively, the effect of both CFD models on the RO system performance was evaluated. The input parameters for the SWRO system model (Cases 1–4) are listed in Table 2. The simulated transmembrane pressure ( Δ P ), flow rate ( Q ), brine salinity ( w b ), permeate salinity ( w p ), permeation flux of water ( J W ) and C P F along the dimensionless length ( X ) (Inlet: X = 0; Outlet: X = 1) for Cases 1–4 are shown in Figure 8a–f, Figure 9a–f, Figure 10a–f and Figure 11a–f, respectively. Obviously, the simulated results at a system-level based on the impermeable wall model are highly consistent with that of the permeable wall model, even for the HPMs systems. Furthermore, the errors of the calculated average permeate salinity ( w ¯ p ), average permeation flux ( J ¯ w ), maximum CP factor, (max (CPF)), recovery rate ( R r ), permeate rate ( Q p ) and S E C are all less than 2% using both the permeable and impermeable wall models. More detailed results are listed in Table 3. As the water permeability increases from 1 L m−2 h−1 bar−1 to 10 L m−2 h−1 bar−1, the maximum CP factor at the system inlet can reach up to about 1.50, which will intensify the membrane fouling. Therefore, it is crucial to enhance the mass transfer with a pressure loss penalty when the SWRO system is operated at higher flux. To achieve this, the spacer design and improvement of the cross velocity or unsteady-state shear technologies [46] are alternative methods.

4. Conclusions and Outlook

In this paper, we propose a multiscale model framework that couples the CFD models with a permeable wall and an impermeable wall, respectively, and the SWRO model at the system level. The simulated results indicate that the CFD simulations (e.g., velocity, pressure) and the system simulations (e.g., SEC, average permeation flux) using the permeable and impermeable wall models are basic unanimous (error < 2%). The error of the estimated Sherwood number at the system inlet using both models is smaller than 10%, even under extremely high-flux conditions (≈200 L m−2 h−1 at system inlet). Furthermore, the impermeable wall model is one-way coupled, which can greatly reduce the computational time compared to the fully coupled permeable wall model. The proposed multiscale framework coupling the impermeable wall model enables the feed spacer design and system design simultaneously for high permeability RO/NF membrane systems. Nevertheless, the impermeable wall model still has some limitations. On the one hand, it cannot accurately predict the local CP due to the assumed constant concentration boundary condition on membrane walls. On the other hand, the correct relations from k ¯ m ,   imp to k ¯ m ,   per must satisfy the condition of J W / k ¯ m ,   imp < 20 .

Author Contributions

Q.Y.: Methodology, Software, Validation, Formal analysis, Investigation, Visualization, Writing-original draft. Y.H.: Methodology, Writing-review and editing, Supervision, Resources. Y.J.: Methodology, Writing-review and editing. J.L.: Conceptualization, Methodology, Visualization, Writing-original draft, Supervision, Date curation. All authors have read and agreed to the published version of the manuscript.

Funding

The author Prof. Dr. Yi Heng acknowledges support provided by Key-Area Research and Development Program of Guangdong Province (No.2021B0101190003), Natural Science Foundation of Guangdong Province (No. 2022A1515011514) and Zhujiang Talent Program of Guangdong Province (2017GC010576). The corresponding author Dr. Jiu Luo thanks support by the project funded by China Postdoctoral Science Foundation (2022M723674).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author (Jiu Luo) upon reasonable request.

Conflicts of Interest

All authors in this paper declare no competing interests.

References

  1. Anis, S.; Hashaikeh, R.; Hilal, N. Microfiltration membrane processes: A review of research trends over the past decade. J. Water Process. Eng. 2019, 32, 100941. [Google Scholar] [CrossRef]
  2. Mouiya, M.; Abourriche, A.; Bouazizi, A.; Benhammou, A.; El Hafiane, Y.; Abouliatim, Y.; Nibou, L.; Oumam, M.; Ouammou, M.; Smith, A.; et al. Flat ceramic microfiltration membrane based on natural clay and Moroccan phosphate for desalination and industrial wastewater treatment. Desalination 2018, 427, 42–50. [Google Scholar] [CrossRef]
  3. Al Aani, S.; Mustafa, T.; Hilal, N. Ultrafiltration membranes for wastewater and water process engineering: A comprehensive statistical review over the past decade. J. Water Process. Eng. 2020, 35, 101241. [Google Scholar] [CrossRef]
  4. Brover, S.; Lester, Y.; Brenner, A.; Sahar-Hadar, E. Optimization of ultrafiltration as pre-treatment for seawater RO desalination. Desalination 2022, 524, 115478. [Google Scholar] [CrossRef]
  5. Zhou, D.; Zhu, L.; Fu, Y.; Zhu, M.; Xue, L. Development of lower cost seawater desalination processes using nanofiltration technologies—A review. Desalination 2015, 376, 109–116. [Google Scholar] [CrossRef]
  6. Wafi, M.; Hussain, N.; Abdalla, O.E.-S.; Al-Far, M.; Al-Hajaj, N.; Alzonnikah, K. Nanofiltration as a cost-saving desalination process. SN Appl. Sci. 2019, 1, 751. [Google Scholar] [CrossRef] [Green Version]
  7. Qasim, M.; Badrelzaman, M.; Darwish, N.; Darwish, N.; Hilal, N. Reverse osmosis desalination: A state-of-the-art review. Desalination 2019, 459, 59–104. [Google Scholar] [CrossRef] [Green Version]
  8. Lim, Y.; Goh, K.; Kurihara, M.; Wang, R. Seawater desalination by reverse osmosis: Current development and future challenges in membrane fabrication—A review. J. Membr. Sci. 2021, 629, 119292. [Google Scholar] [CrossRef]
  9. Galama, A.; Saakes, M.; Bruning, H.; Rijnaarts, H.; Post, J. Seawater predesalination with electrodialysis. Desalination 2014, 342, 61–69. [Google Scholar] [CrossRef]
  10. Al-Amshawee, S.; Yunus, M.; Azoddein, A.; Hassell, D.; Dakhil, I.; Hasan, H. Electrodialysis desalination for water and wastewater: A review. Chem. Eng. J. 2020, 380, 122231. [Google Scholar] [CrossRef]
  11. Yadav, A.; Labhasetwar, P.; Shahi, V. Membrane distillation using low-grade energy for desalination: A review. J. Environ. Chem. Eng. 2021, 9, 105818. [Google Scholar] [CrossRef]
  12. Yadav, A.; Patel, R.; Singh, C.; Labhasetwar, P.; Shahi, V. Experimental study and numerical optimization for removal of methyl orange using polytetrafluoroethylene membranes in vacuum membrane distillation process. Colloids Surf. A Physicochem. Eng. Asp. 2022, 635, 128070. [Google Scholar] [CrossRef]
  13. Feinberg, B.; Ramon, G.; Hoek, E. Thermodynamic analysis of osmotic energy recovery at a reverse osmosis desalination plant. Environ. Sci. Technol. 2013, 47, 2982–2989. [Google Scholar] [CrossRef] [PubMed]
  14. Lim, Y.; Ma, Y.; Chew, J.; Wang, R. Assessing the potential of highly permeable reverse osmosis membranes for desalination: Specific energy and footprint analysis. Desalination 2022, 533, 115771. [Google Scholar] [CrossRef]
  15. Zhang, Z.; Wu, Y.; Luo, L.; Li, G.; Li, Y.; Hu, H. Application of disk tube reverse osmosis in wastewater treatment: A review. Sci. Total Environ. 2021, 792, 148291. [Google Scholar] [CrossRef] [PubMed]
  16. Colla, V.; Branca, T.; Rosito, F.; Lucca, C.; Vivas, B.; Delmiro, V. Sustainable Reverse Osmosis application for wastewater treatment in the steel industry. J. Clean. Prod. 2016, 130, 103–115. [Google Scholar] [CrossRef]
  17. Licona, K.; Geaquinto, L.; Nicolini, J.; Figueiredo, N.; Chiapetta, S.; Habert, A.; Yokoyama, L. Assessing potential of nanofiltration and reverse osmosis for removal of toxic pharmaceuticals from water. J. Water Process. Eng. 2018, 25, 195–204. [Google Scholar] [CrossRef]
  18. Lin, Y.; Chiou, J.; Lee, C. Effect of silica fouling on the removal of pharmaceuticals and personal care products by nanofiltration and reverse osmosis membranes. J. Hazard. Mater. 2014, 277, 102–109. [Google Scholar] [CrossRef]
  19. Bejaoui, I.; Mnif, A.; Hamrouni, B. Performance of Reverse Osmosis and Nanofiltration in the Removal of Fluoride from Model Water and Metal Packaging Industrial Effluent. Sep. Sci. Technol. 2014, 49, 1135–1145. [Google Scholar] [CrossRef]
  20. Ndiaye, P.; Moulin, P.; Dominguez, L.; Millet, J.; Charbit, F. Removal of fluoride from electronic industrial effluent by RO membrane separation. Desalination 2005, 173, 25–32. [Google Scholar] [CrossRef]
  21. Zhao, D.; Japip, S.; Zhang, Y.; Weber, M.; Maletzko, C.; Chung, T. Emerging thin-film nanocomposite (TFN) membranes for reverse osmosis: A review. Water Res. 2020, 173, 115557. [Google Scholar] [CrossRef]
  22. Rahmawati, R.; Bilad, M.; Nawi, N.; Wibisono, Y.; Suhaimi, H.; Shamsuddin, N.; Arahman, N. Engineered spacers for fouling mitigation in pressure driven membrane processes: Progress and projection. J. Environ. Chem. Eng. 2021, 9, 106285. [Google Scholar] [CrossRef]
  23. Lin, W.; Zhang, Y.; Li, D.; Wang, X.; Huang, X. Roles and performance enhancement of feed spacer in spiral wound membrane modules for water treatment: A 20-year review on research evolvement. Water Res. 2021, 198, 117146. [Google Scholar] [CrossRef]
  24. Fimbres-Weihs, G.; Wiley, D. Review of 3D CFD modeling of flow and mass transfer in narrow spacer-filled channels in membrane modules. Chem. Eng. Process. 2010, 49, 759–781. [Google Scholar] [CrossRef]
  25. Anqi, A.; Alkhamis, N.; Oztekin, A. Computational study of desalination by reverse osmosis—Three-dimensional analyses. Desalination 2016, 388, 38–49. [Google Scholar] [CrossRef]
  26. Lin, W.; Shao, R.; Wang, X.; Huang, X. Impacts of non-uniform filament feed spacers characteristics on the hydraulic and anti-fouling performances in the spacer-filled membrane channels: Experiment and numerical simulation. Water Res. 2020, 185, 116251. [Google Scholar] [CrossRef] [PubMed]
  27. Foo, K.; Liang, Y.; Weihs, G. CFD study of the effect of SWM feed spacer geometry on mass transfer enhancement driven by forced transient slip velocity. J. Membr. Sci. 2020, 597, 117643. [Google Scholar] [CrossRef]
  28. Singh, C.; Yadav, A.; Kumar, A. Numerical simulations of the effect of spacer filament geometry and orientation on the performance of the reverse osmosis process. Colloid. Surface A 2022, 650, 129664. [Google Scholar] [CrossRef]
  29. Li, M.; Bui, T.; Chao, S. Three-dimensional CFD analysis of hydrodynamics and concentration polarization in an industrial RO feed channel. Desalination 2016, 397, 194–204. [Google Scholar] [CrossRef]
  30. Lin, W.; Lei, J.; Wang, Q.; Wang, X.-m.; Huang, X. Performance enhancement of spiral-wound reverse osmosis membrane elements with novel diagonal-flow feed channels. Desalination 2022, 523, 115447. [Google Scholar] [CrossRef]
  31. Yadav, A.; Singh, C.; Patel, R.; Kumar, A.; Labhasetwar, P. Investigations on the effect of spacer in direct contact and air gap membrane distillation using computational fluid dynamics. Colloids Surf. A Physicochem. Eng. Asp. 2022, 654, 130111. [Google Scholar] [CrossRef]
  32. Koutsou, C.; Yiantsios, S.; Karabelas, A. A numerical and experimental study of mass transfer in spacer-filled channels: Effects of spacer geometrical characteristics and Schmidt number. J. Membr. Sci. 2009, 326, 234–251. [Google Scholar] [CrossRef]
  33. Completo, C.; Semiao, V.; Geraldes, V. Efficient CFD-based method for designing cross-flow nanofiltration small devices. J. Membr. Sci. 2016, 500, 190–202. [Google Scholar] [CrossRef]
  34. Fimbres-Weihs, G.; Wiley, D. Numerical study of mass transfer in three-dimensional spacer-filled narrow channels with steady flow. J. Membr. Sci. 2007, 306, 228–243. [Google Scholar] [CrossRef]
  35. Shakaib, M.; Hasani, S.; Mahmood, M. CFD modeling for flow and mass transfer in spacer-obstructed membrane feed channels. J. Membr. Sci. 2009, 326, 270–284. [Google Scholar] [CrossRef]
  36. Liang, Y.; Toh, K.; Weihs, G.F. 3D CFD study of the effect of multi-layer spacers on membrane performance under steady flow. J. Membr. Sci. 2019, 580, 256–267. [Google Scholar] [CrossRef]
  37. Toh, K.; Liang, Y.; Lau, W.; Weihs, G. The techno-economic case for coupling advanced spacers to high-permeance RO membranes for desalination. Desalination 2020, 491, 114534. [Google Scholar] [CrossRef]
  38. Geraldes, V.; Afonso, M. Generalized mass-transfer correction factor for nanofiltration and reverse osmosis. AIChE J. 2006, 52, 3353–3362. [Google Scholar] [CrossRef]
  39. Luo, J.; Li, M.; Heng, Y. A hybrid modeling approach for optimal design of non-woven membrane channels in brackish water reverse osmosis process with high-throughput computation. Desalination 2020, 489, 114463. [Google Scholar] [CrossRef]
  40. Chen, Q.; Luo, J.; Heng, Y. High-Throughput Optimal Design of Spacers Using Triply Periodic Minimal Surfaces in BWRO. Separations 2022, 9, 62. [Google Scholar] [CrossRef]
  41. Gu, J.; Luo, J.; Li, M.; Huang, C.; Heng, Y. Modeling of pressure drop in reverse osmosis feed channels using multilayer artificial neural networks. Chem. Eng. Res. Des. 2020, 159, 146–156. [Google Scholar] [CrossRef]
  42. Luo, J.; Li, M.; Hoek, E.; Heng, Y. Supercomputing and machine learning-aided optimal design of high permeability seawater reverse osmosis membrane systems. Sci. Bull. 2023, in press. [Google Scholar] [CrossRef] [PubMed]
  43. Li, F.; Meindersma, W.; Haan, A.; Reith, T. Optimization of commercial net spacers in spiral wound membrane modules. J. Membr. Sci. 2002, 208, 289–302. [Google Scholar] [CrossRef] [Green Version]
  44. Bucs, S.; Radu, A.; Lavric, V.; Vrouwenvelder, J.; Picioreanu, C. Effect of different commercial feed spacers on biofouling of reverse osmosis membrane systems: A numerical study. Desalination 2014, 343, 26–37. [Google Scholar] [CrossRef]
  45. Al-Amoudi, A.; Lovitt, R. Fouling strategies and the cleaning system of NF membranes and factors affecting cleaning efficiency. J. Membr. Sci. 2007, 303, 4–28. [Google Scholar] [CrossRef]
  46. Zamani, F.; Chew, J.; Akhondi, E.; Krantz, W.; Fane, A. Unsteady-state shear strategies to enhance mass-transfer for the implementation of ultrapermeable membranes in reverse osmosis: A review. Desalination 2015, 356, 328–348. [Google Scholar] [CrossRef]
Figure 1. The multiscale modeling framework for SWRO desalination that integrates the permeable wall model and the impermeable wall model.
Figure 1. The multiscale modeling framework for SWRO desalination that integrates the permeable wall model and the impermeable wall model.
Separations 10 00134 g001
Figure 2. (a) Spacer sheet; (b) The computational domain (spacer-filled channel); (c) A spacer unit; (d) Filament A; and (e) Filament B.
Figure 2. (a) Spacer sheet; (b) The computational domain (spacer-filled channel); (c) A spacer unit; (d) Filament A; and (e) Filament B.
Separations 10 00134 g002
Figure 3. The geometric model discretization for (a) the computational domain and (b) the zoomed spacer section.
Figure 3. The geometric model discretization for (a) the computational domain and (b) the zoomed spacer section.
Separations 10 00134 g003
Figure 4. The velocity magnitude distributions with an inlet average velocity magnitude of 0.15 m s−1 using (a) the permeable wall model (Case 1: L P = 1 L m−2 h−1 bar−1) and (b) the impermeable wall model.
Figure 4. The velocity magnitude distributions with an inlet average velocity magnitude of 0.15 m s−1 using (a) the permeable wall model (Case 1: L P = 1 L m−2 h−1 bar−1) and (b) the impermeable wall model.
Separations 10 00134 g004
Figure 5. The water flux with inlet average velocity magnitude of 0.15 m s−1 using the permeable wall model for (a) the commercial membrane ( L P = 1 L m−2 h−1 bar−1) and (b) the HPM ( L P = 5 L m−2 h−1 bar−1) systems.
Figure 5. The water flux with inlet average velocity magnitude of 0.15 m s−1 using the permeable wall model for (a) the commercial membrane ( L P = 1 L m−2 h−1 bar−1) and (b) the HPM ( L P = 5 L m−2 h−1 bar−1) systems.
Separations 10 00134 g005
Figure 6. The pressure magnitude distributions with an inlet average velocity magnitude of 0.15 m s−1 using (a) the permeable wall model and (b) the impermeable wall model.
Figure 6. The pressure magnitude distributions with an inlet average velocity magnitude of 0.15 m s−1 using (a) the permeable wall model and (b) the impermeable wall model.
Separations 10 00134 g006
Figure 7. The friction coefficient ( f ) and and Sherwood number ( Sh ) with respect to various Reynolds numbers (Re) using the impermeable wall model and the permeable wall model for Case 1 (a,e) ( L P = 1 L m−2 h−1 bar−1), (b,f) Case 2 ( L P = 3 L m−2 h−1 bar−1), (c,g) Case 3 ( L P = 5 L m−2 h−1 bar−1) and (d,h) Case 4 ( L P = 10 L m−2 h−1 bar−1).
Figure 7. The friction coefficient ( f ) and and Sherwood number ( Sh ) with respect to various Reynolds numbers (Re) using the impermeable wall model and the permeable wall model for Case 1 (a,e) ( L P = 1 L m−2 h−1 bar−1), (b,f) Case 2 ( L P = 3 L m−2 h−1 bar−1), (c,g) Case 3 ( L P = 5 L m−2 h−1 bar−1) and (d,h) Case 4 ( L P = 10 L m−2 h−1 bar−1).
Separations 10 00134 g007
Figure 8. The simulated results at the system level using the impermeable and permeable wall models with the feed salinity ( w ¯ b ,   0 = 35,000 ppm), water recovery ( R r ,   0 = 50%), pump efficiency (0.85), energy recovery device efficiency (0.95) and water permeability ( L p = 1 L m−2 h−1 bar−1). (a) Transmembrane pressure, Δ P (bar), (b) Flow rate, Q (m3 h−1), (c) Brine salinity, w b (ppm), (d) Permeate salinity, w p (ppm), (e) Permeation flux of water, J W (L m−2 h−1), (f) CP factor, C P F (−) along the dimensionless length ( X ).
Figure 8. The simulated results at the system level using the impermeable and permeable wall models with the feed salinity ( w ¯ b ,   0 = 35,000 ppm), water recovery ( R r ,   0 = 50%), pump efficiency (0.85), energy recovery device efficiency (0.95) and water permeability ( L p = 1 L m−2 h−1 bar−1). (a) Transmembrane pressure, Δ P (bar), (b) Flow rate, Q (m3 h−1), (c) Brine salinity, w b (ppm), (d) Permeate salinity, w p (ppm), (e) Permeation flux of water, J W (L m−2 h−1), (f) CP factor, C P F (−) along the dimensionless length ( X ).
Separations 10 00134 g008
Figure 9. The simulated results at the system level using the impermeable and permeable wall models with the feed salinity ( w ¯ b ,   0 = 35,000 ppm), water recovery ( R r ,   0 = 50%), pump efficiency (0.85), energy recovery device efficiency (0.95) and water permeability ( L p = 3 L m−2 h−1 bar−1). (a) Transmembrane pressure, Δ P (bar), (b) Flow rate, Q (m3 h−1), (c) Brine salinity, w b (ppm), (d) Permeate salinity, w p (ppm), (e) Permeation flux of water, J W (L m−2 h−1), (f) CP factor, C P F (−) along the dimensionless length ( X ).
Figure 9. The simulated results at the system level using the impermeable and permeable wall models with the feed salinity ( w ¯ b ,   0 = 35,000 ppm), water recovery ( R r ,   0 = 50%), pump efficiency (0.85), energy recovery device efficiency (0.95) and water permeability ( L p = 3 L m−2 h−1 bar−1). (a) Transmembrane pressure, Δ P (bar), (b) Flow rate, Q (m3 h−1), (c) Brine salinity, w b (ppm), (d) Permeate salinity, w p (ppm), (e) Permeation flux of water, J W (L m−2 h−1), (f) CP factor, C P F (−) along the dimensionless length ( X ).
Separations 10 00134 g009
Figure 10. The simulated results at the system level using the impermeable and permeable wall models with the feed salinity ( w ¯ b ,   0 = 35,000 ppm), water recovery ( R r ,   0 = 50%), pump efficiency (0.85), energy recovery device efficiency (0.95) and water permeability ( L p = 5 L m−2 h−1 bar−1). (a) Transmembrane pressure, Δ P (bar), (b) Flow rate, Q (m3 h−1), (c) Brine salinity, w b (ppm), (d) Permeate salinity, w p (ppm), (e) Permeation flux of water, J W (L m−2 h−1), (f) CP factor, C P F (−) along the dimensionless length ( X ).
Figure 10. The simulated results at the system level using the impermeable and permeable wall models with the feed salinity ( w ¯ b ,   0 = 35,000 ppm), water recovery ( R r ,   0 = 50%), pump efficiency (0.85), energy recovery device efficiency (0.95) and water permeability ( L p = 5 L m−2 h−1 bar−1). (a) Transmembrane pressure, Δ P (bar), (b) Flow rate, Q (m3 h−1), (c) Brine salinity, w b (ppm), (d) Permeate salinity, w p (ppm), (e) Permeation flux of water, J W (L m−2 h−1), (f) CP factor, C P F (−) along the dimensionless length ( X ).
Separations 10 00134 g010
Figure 11. The simulated results at the system level using the impermeable and permeable wall models with the feed salinity ( w ¯ b ,   0 = 35,000 ppm), water recovery ( R r ,   0 = 50%), pump efficiency (0.85), energy recovery device efficiency (0.95) and water permeability ( L p = 10 L m−2 h−1 bar−1). (a) Transmembrane pressure, Δ P (bar), (b) Flow rate, Q ( m3 h−1), (c) Brine salinity, w b (ppm), (d) Permeate salinity, w p (ppm), (e) Permeation flux of water, J W (L m−2 h−1), (f) CP factor, C P F (−) along the dimensionless length ( X ).
Figure 11. The simulated results at the system level using the impermeable and permeable wall models with the feed salinity ( w ¯ b ,   0 = 35,000 ppm), water recovery ( R r ,   0 = 50%), pump efficiency (0.85), energy recovery device efficiency (0.95) and water permeability ( L p = 10 L m−2 h−1 bar−1). (a) Transmembrane pressure, Δ P (bar), (b) Flow rate, Q ( m3 h−1), (c) Brine salinity, w b (ppm), (d) Permeate salinity, w p (ppm), (e) Permeation flux of water, J W (L m−2 h−1), (f) CP factor, C P F (−) along the dimensionless length ( X ).
Separations 10 00134 g011
Table 1. The detailed geometrical parameters of the CFD models [44].
Table 1. The detailed geometrical parameters of the CFD models [44].
Geometrical ParametersValue
Spacer unit L 1 ,   A (μm)234
L 2 ,   A (μm)389
L 3 ,   A (μm)500
L 4 ,   A (μm)696
L 5 ,   A (μm)566
L 1 ,   B (μm)281
L 2 ,   B (μm)535
L 3 ,   B (μm)92
L 4 ,   B (μm)878
L 5 ,   B (μm)599
D 1 ,   A (μm)422
D 2 ,   A (μm)223
D 1 ,   B (μm)445
D 2 ,   B (μm)223
D tot (μm)863
α (°)90
Computational domainLength, L (mm)16.864
Width, W (mm)3.3729
Height, H (mm)0.853
Porosity, ε (-)0.90
Hydraulic diameter, D H (mm)1.058
Table 2. The input parameters for SWRO systems in Cases 1–4.
Table 2. The input parameters for SWRO systems in Cases 1–4.
Input ParametersValue
Operating conditionsInlet transmembrane pressure, Δ P 0 (bar)60
Inlet average velocity magnitude, U 0 (m s−1)0.2
Properties of feed soluteFeed salinity, w b ,   0 (ppm)35,000
Density, ρ (kg m−3)1021
Viscosity, μ (Pa s)9.41 × 10−4
Diffusion coefficient, D s (m2 s−1)1.45 × 10−9
Reflection coefficient, σ (−)1
Osmotic pressure coefficient, φ (bar)805.1
Membrane propertiesWater permeability, L P (L m−2 h−1 bar−1)1 (Case1); 3 (Case 2); 5 (Case 3); 10 (Case 4)
Salt permeability, B (L m−2 h−1)0.05
Membrane module parametersModule length parallel to flow, l x (m)1
Module length perpendicular to flow, l y (m)1
Number of feed spacers per element, n sp 23
Membrane area per membrane module, A 0 (m2)37.2
Module configurationsNumber of membrane elements per pressure vessel, n mem 8 (Case1); 5 (Case 2); 3 (Case 3); 2 (Case 4)
Number of pressure vessels, N pv 30
EfficienciesEfficiency of high-pressure pump, η pump (−)0.85
Efficiency of energy recovery device, η R (−)0.95
Table 3. Output results and errors analysis at the system-level coupling with the impermeable wall model (M2) compared with that of using the permeable wall model (M1) with respect to Case 1 ( L P = 1 L m−2 h−1 bar−1), Case 2 ( L P = 3 L m−2 h−1 bar−1), Case 3 ( L P = 5 L m−2 h−1 bar−1) and Case 4 ( L P = 10 L m−2 h−1 bar−1).
Table 3. Output results and errors analysis at the system-level coupling with the impermeable wall model (M2) compared with that of using the permeable wall model (M1) with respect to Case 1 ( L P = 1 L m−2 h−1 bar−1), Case 2 ( L P = 3 L m−2 h−1 bar−1), Case 3 ( L P = 5 L m−2 h−1 bar−1) and Case 4 ( L P = 10 L m−2 h−1 bar−1).
Output ResultsValue
Case 1Case 2Case 3Case 4
M1M2ErrorM1M2ErrorM1M2ErrorM1M2Error
w ¯ p (ppm)1381390.55%88890.79%56561.22%40401.47%
J ¯ W (L m−2 h−1)18.418.30.38%34.133.90.64%54.554.00.99%81.079.91.26%
max (CPF) (−)1.081.080.32%1.211.210.04%1.321.310.65%1.501.481.86%
R r (−)0.430.430.38%0.500.490.64%0.480.470.99%0.470.471.26%
Q p (m3 h−1)1641640.38%1901890.64%1821810.99%1811781.26%
S E C (kWh m−3)2.172.170.06%2.112.110.08%2.122.120.12%2.112.120.14%
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

Yang, Q.; Heng, Y.; Jiang, Y.; Luo, J. Multiscale Analysis of Permeable and Impermeable Wall Models for Seawater Reverse Osmosis Desalination. Separations 2023, 10, 134. https://doi.org/10.3390/separations10020134

AMA Style

Yang Q, Heng Y, Jiang Y, Luo J. Multiscale Analysis of Permeable and Impermeable Wall Models for Seawater Reverse Osmosis Desalination. Separations. 2023; 10(2):134. https://doi.org/10.3390/separations10020134

Chicago/Turabian Style

Yang, Qingqing, Yi Heng, Ying Jiang, and Jiu Luo. 2023. "Multiscale Analysis of Permeable and Impermeable Wall Models for Seawater Reverse Osmosis Desalination" Separations 10, no. 2: 134. https://doi.org/10.3390/separations10020134

APA Style

Yang, Q., Heng, Y., Jiang, Y., & Luo, J. (2023). Multiscale Analysis of Permeable and Impermeable Wall Models for Seawater Reverse Osmosis Desalination. Separations, 10(2), 134. https://doi.org/10.3390/separations10020134

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