Next Article in Journal
Effect of Short-Term Contact with C1–C4 Monohydric Alcohols on the Water Permeance of MPD-TMC Thin-Film Composite Reverse Osmosis Membranes
Next Article in Special Issue
Characterization and Assessment of a Novel Plate and Frame MD Module for Single Pass Wastewater Concentration–FEED Gap Air Gap Membrane Distillation
Previous Article in Journal
Hybrid Hollow Fiber Nanofiltration–Calcite Contactor: A Novel Point-of-Entry Treatment for Removal of Dissolved Mn, Fe, NOM and Hardness from Domestic Groundwater Supplies
Previous Article in Special Issue
Evaluation of Permeate Quality in Pilot Scale Membrane Distillation Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

CFD Investigation of Spacer-Filled Channels for Membrane Distillation

1
Dipartimento di Ingegneria, Università degli Studi di Palermo, viale delle Scienze ed.6, 90128 Palermo, Italy
2
Center of Excellence in Desalination Technology, King Abdulaziz University, Jeddah 21589, Saudi Arabia
3
Mechanical Engineering Department, King Abdulaziz University, Jeddah 21589, Saudi Arabia
*
Author to whom correspondence should be addressed.
Membranes 2019, 9(8), 91; https://doi.org/10.3390/membranes9080091
Submission received: 11 June 2019 / Revised: 12 July 2019 / Accepted: 15 July 2019 / Published: 25 July 2019
(This article belongs to the Special Issue Membrane Distillation Process)

Abstract

:
The membrane distillation (MD) process for water desalination is affected by temperature polarization, which reduces the driving force and the efficiency of the process. To counteract this phenomenon, spacer-filled channels are used, which enhance mixing and heat transfer but also cause higher pressure drops. Therefore, in the design of MD modules, the choice of the spacer is crucial for process efficiency. In the present work, different overlapped spacers are investigated by computational fluid dynamics (CFD) and results are compared with experiments carried out with thermochromic liquid crystals (TLC). Results are reported for different flow attack angles and for Reynolds numbers (Re) ranging from ~200 to ~800. A good qualitative agreement between simulations and experiments can be observed for the areal distribution of the normalized heat transfer coefficient. Trends of the average heat transfer coefficient are reported as functions of Re for the geometries investigated, thus providing the basis for CFD-based correlations to be used in higher-scale process models.

1. Introduction

In regions where the existing freshwater sources do not meet the water demand, production of potable water from saline and brackish waters is needed [1]. Among the desalination processes, membrane distillation (MD) has recently attracted significant interest due to its features, which make MD a process easily scalable to small–medium size [2,3,4]. MD is a thermally driven membrane-based process where the transport of water vapor is realized through a micro-porous hydrophobic membrane, which leads to a theoretical 100% rejection of salts. The separation process is driven by the vapor pressure difference at the membrane/liquid interfaces due to the temperature gradient existing between the two sides of the membrane. MD is characterized by its compactness and low operational pressure requirement as well as low working temperature. These features make it possible to power MD by low grade heat (waste heat) or renewable energy and boost its applicability in remote areas and small scale production [5,6,7]. Another important distinct feature of MD is its ability to desalinate highly saline brines and wastewater. MD is the most appropriate option to increase the recovery ratio of multi-stage flash (MSF) and reverse osmosis (RO) plants due to it its ability to treat feeds with high salinity and impurities. One main advantage of MD technology is that evaporation and condensation surfaces are tightly packed, thus leading to a compact system with low capital cost per unit product. The classification of MD systems is related to the adopted condensation methods. In general, MD systems are classified into four basic simple configurations designed relative to the permeate vapor condensation method used. These include direct contact membrane distillation (DCMD), air gap membrane distillation (AGMD), sweeping gas membrane distillation (SGMD), and vacuum membrane distillation (VMD) [5]. The succession of these different configurations was dictated by a need to improve the performance of the process. DCMD is the simplest configuration that has liquid phases (hot feed and cold permeate) in direct contact with both sides of the porous hydrophobic membrane. The vapor diffusion path is limited to the thickness of the membrane, therefore reducing the mass and heat transfer resistances and thus enhancing the process [8,9].
Despite these advantages and the many experimental and theoretical studies conducted in recent years, MD is not yet fully developed from the commercial point of view because it is still highly energy demanding (typical thermal energy consumption ~100 kWh/m3) compared, for example, to multiple effect distillation, or MED (thermal energy consumption ~40 kWh/m3 for large plants) and reverse osmosis, or RO (mechanical energy consumption 2–4 kWh/m3) [10,11,12]. Other reasons why MD is still on a pilot scale are related to membrane, module design, and low water flux. One of the main causes responsible for low water flux is the temperature polarization phenomenon, while the concentration polarization has a negligible effect on the permeate flow rate reduction [13]. Temperature polarization can lead to a dramatic reduction of the actual transmembrane temperature driving force [10,13]. A reduction of these effects and a performance improvement could be achieved by the optimization of the geometry of the channels and of the adopted spacers acting as mixing promoters.
Computational fluid dynamics (CFD) is suitable to investigate flow field and transport phenomena in spacer-filled channels. Since the 1990s, CFD has been used to study channels with spacers, by 2D [14,15,16] or 3D simulations [17,18,19,20], in the steady regime [16,21,22,23,24] or considering unsteady flows [25,26,27]. In the unsteady regime, very accurate direct numerical simulations (DNS) become computationally expensive, and turbulence models may be more suitable [28,29]. In regard to CFD results validation, in most of the literature only global quantities from CFD and experiments are compared, such as the friction coefficient and the average Nusselt or Sherwood numbers [25,30,31]. However, local data (i.e., flow fields and temperature distributions) are more suitable than global quantities because they better allow one to check if the assumptions made in CFD simulations (e.g., assuming simple boundary conditions or neglecting geometrical irregularities) reproduce sufficiently well the actual physics of the system. Experiments based on thermochromic liquid crystals (TLC) and digital image processing (DIP) can provide local results, such as temperature and heat transfer coefficient distributions, which can also be processed to obtain global quantities [27,32,33,34,35].
The aim of the present work is to characterize and design novel prototypes of spacers for membrane separation modules. We performed 3D simulations from the laminar to the unsteady regime by using, in this latter case, the shear stress transport (SST) turbulence model. For selected test cases, CFD results were compared with experimental data obtained in a purpose built test facility where the TLC technique was used. The validated CFD model was then used to investigate several spacers in different conditions.

2. Materials and Methods

2.1. Geometry, Operating Conditions, and Performance Parameters

An overlapped spacer of the kind examined in the present study is characterized by a few geometric parameters (Figure 1), namely: the filament diameter (d); the orthogonal distance between filaments (pitch, P); the intrinsic angle between filaments (θ); and the flow attack angle (α) between the directions of the main flow and of the filament layer adjacent to a specified (e.g., top) wall. Note that d and P are assumed to be the same for the two overlapped layers of filaments. In the present study, for the above parameters the values reported in Table 1 were investigated by both experiments and CFD simulations. In order to exactly reproduce the geometry of the experimental module, the geometrical model for CFD was created from the same CAD files used for the 3D printing of the spacer. The channel height H, chosen to be within the range of typical spacer thickness values in MD [18,36,37], was not exactly equal to 2d (twice the filament diameter), i.e., 4 mm, but slightly less (3.8 mm) because a certain amount of interpenetration existed. Thus, the pitch to channel height ratio was P/H = 2.63.
Both in the experiments and in the computations, the bottom wall was assumed adiabatic and heat transfer occurring only through the top one (active wall), in which distributions of temperature, wall heat flux, and heat transfer coefficient were assessed. This last quantity was locally defined as:
h = q T b u l k T w a l l
in which Tbulk is the local fluid bulk temperature, while q″ and Twall are the local wall heat flux and wall temperature, respectively. For the reasons discussed in [23], an average heat transfer coefficient is better defined not as the surface average 〈h〉 of the local h, but rather as
h a v g = q T b u l k T w a l l
in which the symbol 〈 〉 denotes averaging operation on the active wall.
The Reynolds number (Re) is defined as:
Re = U D h ν
where Dh is the hydraulic diameter, conventionally assumed to be equal to twice the channel height (7.6 mm) as in a void plane channel of infinite width. U is the approach, or superficial, velocity, defined as the velocity that the fluid would have if the channel were void of any spacer; it is equal to the ratio between the flow rate and the passage area U = Q/Ap (in the experiments Ap = 9.12 × 10−4 m2). ν is the kinematic viscosity, which was assumed equal to 6.78 × 10−7 m2/s in the cases investigated (water at ~43 °C).
For each configuration, the values of the flow rate investigated and the corresponding values of the Reynolds number are those given in Table 2. Flow rates were selected based on the experiments’ limitations (the range of the flowmeters, measurement uncertainty at low flow rates, and pressure build up at the Plexiglas channel), and to enable comparison of results with several studies conducted in the literature by other researchers within the same range of values.

2.2. Experimental Set-Up and Procedures

The experimental results were obtained in a previous study [35] using the setup shown in Figure 2. It is composed of two narrow (~4 mm height) hot and cold water channels separated by a 1 mm thick polycarbonate sheet. These channels are enclosed by an upper and a lower 2 cm thick transparent Plexiglass sheet. Test spacers were custom designed and fabricated using 3D printing with overall dimensions of 24 cm × 24 cm × 0.38 cm. The spacers were designed to fit and fill the hot channel while touching a thermochromic liquid crystals (TLC) sheet, which was made to adhere to the lower surface of the polycarbonate sheet. TLC sheets adopted had a color play from 35 to 40 °C. The hot and cold water temperatures were measured using RTD sensors located at the inlet and outlet of the two channels, while the flow rates were measured using flow sensors and rotameters. Julabo temperature-controlled heating and cooling circulators were used to control the hot and cold water temperatures at the channel inlets. Further details on experimental procedure and data processing can be found in [35]. The color distribution of the TLC at the upper surface of the hot channel was monitored by a high resolution camera, and images (1920 × 1080 pixels) were captured for different feed temperatures and flow rates (from 1 to 4 L/min), while the flow rate in the cold channel was kept at 8 L/min.
These images were processed by first converting them to the HSV (hue, saturation, value) format. It is known that the hue component of this format correlates with temperature but requires prior calibration. At the beginning of each spacer experiment, a TLC calibration test was preliminarily performed and a best fit polynomial curve was obtained for the temperature as a function of hue. During the calibration test both the upper and lower channels were fed by the same hot water source. In the real experiments, the cold water was set to 30 °C and the hot water was set to 43 °C, allowing a 13 °C difference between the two channels for heat transfer. Hot and cold water inlet temperatures were selected after several pre-tests to ensure that the temperature distributions on the TCL surface were clearly visible and fell within the color play range for all the flow rates.
It should be observed that in order to characterize the temperature polarization phenomenon, it is sufficient to reproduce the thermal field and convective heat transfer occurring in the spacer-filled channel only. This can be achieved by imposing at the active wall a realistic thermal boundary condition, representative of those occurring in MD modules, independent of whether there are or not an actual membrane and a vapor flux.

2.3. CFD Investigations

The problem was described by the continuity, momentum (Navier–Stokes), and energy equations for water at T = 43 °C and p = 1 bar, which were solved by the ANSYS-CFX® code (ANSYS, USA). In the simulations, the unit periodic cell approach was used [22,34,38], which simulates flow and heat exchange phenomena in periodic lattices under fully developed conditions (i.e., at sufficient distance from inlets and outlets). A source term, accounting for the large-scale temperature gradient, was implemented in the energy equation, and a body force per unit volume, accounting for the large-scale pressure gradient, in the momentum equations. In each run, this latter term was dynamically adjusted to obtain the required Re value (corresponding to one of those achieved in the experiments).
As reported in the literature, for geometries similar to those examined here [19,26], the fluid flow becomes unsteady for Re > 350. For this reason, steady-state laminar simulations were carried out for Re = 205 and 305, while, for 410 ≤ Re ≤ 820, the shear stress transport (SST) turbulence model was used, which is a blend between the k-ω model near the walls and the k-ε model in the outer region. As demonstrated in previous works [29], for the present transitional flows, ω-based models, which fully resolve the near-wall layer, are preferable to ε-based models, which make use of wall functions; among them, the SST model gives the most accurate predictions in terms of both distributions and average values of the heat transfer coefficient.
The computational domains (different according to the intrinsic angle) are shown in Figure 3, while details of the computational finite volume grid for case (c) are shown in Figure 4. Only four corner blocks were meshed with tetrahedra; the remaining volume was meshed with mapped hexahedra or with sweep mode in the regions adjacent to the cylindrical surfaces. In regard to the boundary conditions, periodicity was imposed at the walls perpendicular to the filaments; at the other walls a no slip condition was set. The fluid–filament interfaces and the bottom wall were assumed adiabatic, while at the top wall a mixed condition (Robin) was imposed by prescribing an external temperature of 30 °C and an external heat transfer coefficient equal to λ/δ, λ being the fluid’s thermal conductivity and δ the channel thickness.

3. Results and Discussion

3.1. CFD Prediction of Temperature and Heat Transfer Coefficient Distributions and Comparison with Experiments

For the sake of brevity, results are illustrated here in detail only for the two geometrical configurations in Figure 5. Configuration 1 is characterized by an intrinsic angle θ between filaments of 60° and a flow attack angle α of 30°, while configuration 2 is characterized by θ = 90° and α = 60°. The solid profiles indicate the spacer filaments, which touch the upper, thermally active, wall. The results of the whole computational campaign (aiming at the determination of the overall dependence of 〈h〉 on Re for each geometrical configuration, to be used in process modelling applications) are reported in the final section.
For each configuration, three Reynolds numbers were considered: the lowest (205), an intermediate one (410), and the highest (820) of those reported in Table 2. Distributions of h, normalized by the mean value 〈h〉, are shown. Experimental distributions were obtained by post-processing the TLC images, as described in Section 2.2, and averaging the results over 3 × 3 unit cells to reduce dispersion.

3.1.1. Configuration 1

Results for configuration 1 are reported in Figure 6. For clarity purposes, the insets show the direction of the flow and the arrangement of the filaments.
The comparison shows a fair overall agreement between simulations and experiments in the heat transfer coefficient distribution. At the lowest Re, the h distribution was roughly symmetric between the upstream (left) and downstream (right) halves of the cell, showing that inertial effects were negligible. As Re increased, a stripe of low h values grew shortly after the upstream upper filament, associated with a region of separated flow (as confirmed by an analysis of the CFD-predicted flow fields). At the same time, a band of high h values appeared shortly before the downstream upper filament, in correspondence with the reattachment of the separated shear layer. Some discrepancy existed between experiments and CFD results, which predicted stripes of high h, narrower, longer, and inclined with respect to the filaments. Note also that experimental results showed residual fluctuations, not eliminated by the spatial averaging performed (which only involved nine unit cells); CFD simulations, by their nature themselves, yielded time-averaged values of h (and of any other variable) and thus appeared smooth and regular.

3.1.2. Configuration 2

Results for configuration 2 are reported in Figure 7. As in the previous Figure 6, the inset shows the direction of the flow and the arrangement of the filaments.
In this case only a partial agreement between simulation and experiments could be registered. The numerical simulations predicted, both under laminar flow conditions (Figure 7b) and using the turbulence model (Figure 7d,f), a band where the most intense heat transfer occurs, located shortly before the downstream (upper) spacer filament, and a secondary band located shortly upstream of it (i.e., closer to the center of the unit cell). These features were not present in the experimental results, which showed a broad maximum of h only slightly upstream of the upper filament and centered on the opposite filament (the one that did not touch the active upper wall).
General levels and the overall distribution of h/havg, however, were satisfactorily reproduced. As in the previous configuration, experimental results exhibited a considerable amount of residual fluctuations, not completely damped by the 9-cell spatial averaging performed, while CFD results were smooth and regular.

3.2. CFD Prediction of Average Heat Transfer Coefficient and Comparison with Experiments

Figure 8 and Figure 9 show the average heat transfer coefficient havg (Equation (2)) as a function of the Reynolds number for configurations 1 (θ = 60°, α = 30°) and 2 (θ = 90°, α = 60°), respectively.
A fair agreement between the experiments and the CFD simulations could be noticed. CFD results for configuration 1 (Figure 8) overpredicted havg by ~13% when Re was lower than ~600; at higher Re, a slight underprediction was registered (~14%). In Figure 9, results for configuration 2 showed that CFD overestimates havg by ~14%, with the maximum difference between experiments and simulations (~25% overprediction) observed at Re ≈ 500–600.

3.3. Summary Results for All Geometries Investigated

In Figure 10 and Figure 11, the CFD results obtained for all the configurations reported in Table 1 are shown. In Figure 10, havg is plotted as a function of Reynolds number for the three configurations in which the flow direction bisects the intrinsic angle θ, i.e., θ = 30°, α = 15°; θ = 60°, α = 30°; and θ = 90°, α = 45°. Curves of havg against Re, for the configuration with θ = 90°, are reported in Figure 11 considering flow attack angles values lower (graph a) or higher (graph b) than 45°.
From Figure 10 and Figure 11, it is possible to identify the configuration with θ = 90°, α = 60° as that providing the highest heat transfer coefficient.

4. Conclusions

Flow and heat transfer were predicted by computational fluid dynamics for channels provided with spacers consisting of two overlapped layers with a fixed pitch to height ratio (~2.63) and different values of the intrinsic angle between the filaments and of the flow attack angle. The Reynolds number ranged between 200 and 820.
In a previous work, experiments were carried out for the same configurations and Reynolds numbers by using thermochromic liquid crystals (TLC), which provide the distribution of the heat transfer coefficient h on the active wall.
For two selected cases, distributions of h obtained by CFD were compared with the experimental distributions, obtained by post-processing the real images of the TLC. From the comparison, a fair agreement concerning the shape of the h distributions was registered; the match was better at low Reynolds number, when no turbulence model was used, and worse at high Re, when simulations relied on the SST turbulence model.
For the same test cases, the average heat transfer coefficient was also fairly well predicted by CFD simulations, with discrepancies of the order of 10–20%, which can be regarded as only minor in view of the complex geometry and flow structure.
Finally, the results of the whole computational campaign were reported in terms of average heat transfer coefficient as a function of the Reynolds number and can be used in process modelling applications involving spacers with the same geometrical features as those investigated here.
In MD systems, geometrical configurations yielding the highest possible heat transfer coefficient h should be adopted in order to minimize the membrane surface (the largest component of plant cost). With reference to the geometries investigated in this work, for Reynolds numbers above 400–500, this optimal configuration seems to be characterized by an intrinsic angle θ of 90° and a flow attack angle α of 60°.

Author Contributions

Conceptualization, A.C., M.C., G.M., and M.A.; methodology, M.C. and M.L.C.; software, M.C. and M.L.C.; validation, M.A., N.T., and S.B.; formal analysis, M.C. and M.A.; investigation, M.L.C., N.T., and S.B.; resources, M.A. and A.C.; data curation, M.L.C. and M.C.; writing—original draft preparation, M.L.C.; writing—review and editing, A.C., M.C., and M.A.; visualization, M.L.C., N.T., and S.B.; supervision, A.C., M.C., and M.A.; project administration, M.A.; funding acquisition, M.A.

Funding

This project was supported by the NSTIP strategic technologies program in the Kingdom of Saudi Arabia, Project No. 11WAT1576-03.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gude, V.G.; Nirmalakhandan, N. Sustainable desalination using solar energy. Energy Convers. Manag. 2010, 51, 2245–2251. [Google Scholar] [CrossRef]
  2. Summers, E.K.; Arafat, H.A.; Lienhard, V.J.H. Energy efficiency comparison of single-stage membrane distillation (MD) desalination cycles in different configurations. Desalination 2012, 290, 54–66. [Google Scholar] [CrossRef]
  3. Zuo, G.; Guan, G.; Wang, R. Numerical modeling and optimization of vacuum membrane distillation module for low-cost water production. Desalination 2014, 339, 1–9. [Google Scholar] [CrossRef]
  4. Dong, G.; Kim, J.F.; Kim, J.H.; Drioli, E.; Lee, Y.M. Open-source predictive simulators for scale-up of direct contact membrane distillation modules for seawater desalination. Desalination 2017, 402, 72–87. [Google Scholar] [CrossRef]
  5. Drioli, E.; Ali, A.; Macedonio, F. Membrane distillation: Recent developments and perspectives. Desalination 2015, 356, 56–84. [Google Scholar] [CrossRef]
  6. González-Bravo, R.; Elsayed, N.A.; Ponce-Ortega, J.M.; Nápoles-Rivera, F.; El-Halwagi, M.M. Optimal design of thermal membrane distillation systems with heat integration with process plants. Appl. Therm. Eng. 2015, 75, 154–166. [Google Scholar] [CrossRef]
  7. Subramani, A.; Jacangelo, J.G. Emerging desalination technologies for water treatment: A critical review. Water Res. 2015, 75, 164–187. [Google Scholar] [CrossRef]
  8. Pangarkar, B.L.; Deshmukh, S.K.; Sapkal, V.S.; Sapkal, R.S. Review of membrane distillation process for water purification. Desalin. Water Treat. 2014, 3994, 1–23. [Google Scholar] [CrossRef]
  9. Long, R.; Lai, X.; Liu, Z.; Liu, W. Direct contact membrane distillation system for waste heat recovery: Modelling and multi-objective optimization. Energy 2018, 148, 1060–1068. [Google Scholar] [CrossRef]
  10. Wang, P.; Chung, T.S. Recent advances in membrane distillation processes: Membrane development, configuration design and application exploring. J. Membr. Sci. 2015, 474, 39–56. [Google Scholar] [CrossRef]
  11. Fritzmann, C.; Lowenberg, J.; Wintgens, T.; Melin, T. State-of-the-art of reverse osmosis desalination. Desalination 2007, 216, 1–76. [Google Scholar] [CrossRef]
  12. Ghaffour, N.; Lattemann, S.; Missimer, T.; Ng, K.C.; Sinha, S.; Amy, G. Renewable energy-driven innovative energy-efficient desalination technologies. Appl. Energy 2014, 136, 1155–1165. [Google Scholar] [CrossRef] [Green Version]
  13. Martínez-Díez, L.; Vázquez-González, M. Temperature and concentration polarization in membrane distillation of aqueous salt solutions. J. Membr. Sci. 1999, 156, 265–273. [Google Scholar] [CrossRef]
  14. Cao, Z. CFD simulations of net-type turbulence promoters in a narrow channel. J. Membr. Sci. 2001, 185, 157–176. [Google Scholar] [CrossRef]
  15. Sousa, P.; Soares, A.; Monteiro, E.; Rouboa, A. A CFD study of the hydrodynamics in a desalination membrane filled with spacers. Desalination 2014, 349, 22–30. [Google Scholar] [CrossRef]
  16. Amokrane, M.; Sadaoui, D.; Koutsou, C.; Karabelas, A.; Dudeck, M.; Karabelas, A. A study of flow field and concentration polarization evolution in membrane channels with two-dimensional spacers during water desalination. J. Membr. Sci. 2015, 477, 139–150. [Google Scholar] [CrossRef]
  17. Koutsou, C.; Yiantsios, S.; Karabelas, A.; 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]
  18. Taamneh, Y.; Bataineh, K. Improving the performance of direct contact membrane distillation utilizing spacer-filled channel. Desalination 2017, 408, 25–35. [Google Scholar] [CrossRef]
  19. Gurreri, L.; Tamburini, A.; Cipollina, A.; Micale, G.; Ciofalo, M. Flow and mass transfer in spacer-filled channels for reverse electrodialysis: A CFD parametrical study. J. Membr. Sci. 2016, 497, 300–317. [Google Scholar] [CrossRef]
  20. La Cerva, M.F.; Di Liberto, M.; Gurreri, L.; Tamburini, A.; Cipollina, A.; Micale, G.; Ciofalo, M. Coupling CFD with a one-dimensional model to predict the performance of reverse electrodialysis stacks. J. Membr. Sci. 2017, 541, 595–610. [Google Scholar] [CrossRef] [Green Version]
  21. Saeed, A.; Vuthaluru, R.; Vuthaluru, H.B. Investigations into the effects of mass transport and flow dynamics of spacer filled membrane modules using CFD. Chem. Eng. Res. Des. 2015, 93, 79–99. [Google Scholar] [CrossRef]
  22. Gurreri, L.; Tamburini, A.; Cipollina, A.; Micale, G.; Ciofalo, M. CFD simulation of mass transfer phenomena in spacer filled channels for reverse electrodialysis applications. Chem. Eng. Trans. 2013, 32, 1879–1884. [Google Scholar]
  23. La Cerva, M.; Ciofalo, M.; Gurreri, L.; Tamburini, A.; Cipollina, A.; Micale, G. On some issues in the computational modelling of spacer-filled channels for membrane distillation. Desalination 2017, 411, 101–111. [Google Scholar] [CrossRef] [Green Version]
  24. Al-sharif, S.; Albeirutty, M.; Cipollina, A.; Micale, G. Modelling flow and heat transfer in spacer-filled membrane distillation channels using open source CFD code. Desalination 2013, 311, 103–112. [Google Scholar] [CrossRef]
  25. Koutsou, C.P.; Yiantsios, S.G.; Karabelas, A.J. Direct numerical simulation of flow in spacer-filled channels: Effect of spacer geometrical characteristics. J. Membr. Sci. 2007, 291, 53–69. [Google Scholar] [CrossRef]
  26. Mojab, S.M.; Pollard, A.; Pharoah, J.G.; Beale, S.B.; Hanff, E.S. Unsteady laminar to turbulent flow in a spacer-filled channel. Flow Turbul. Combust. 2014, 92, 563–577. [Google Scholar] [CrossRef]
  27. Ponzio, F.N.; Tamburini, A.; Cipollina, A.; Micale, G.; Ciofalo, M. Experimental and computational investigation of heat transfer in channels filled by woven spacers. Int. J. Heat Mass Transf. 2017, 104, 163–177. [Google Scholar] [CrossRef] [Green Version]
  28. Shakaib, M.; Hasani, S.M.F.; Ahmed, I.; Yunus, R.M. A CFD study on the effect of spacer orientation on temperature polarization in membrane distillation modules. Desalination 2012, 284, 332–340. [Google Scholar] [CrossRef]
  29. Ciofalo, M.; La Cerva, M.; Di Liberto, M. Selection of Turbulence Models for Low-Reynolds Number Turbulent Flow in Spacer-Filled Channels. In Proceedings of the 36th UIT Heat Transfer Conference, Catania, Italy, 25–27 June 2018. [Google Scholar]
  30. Chang, H.; Hsu, J.A.; Chang, C.L.; Ho, C.D.; Cheng, T.W. Simulation study of transfer characteristics for spacer-filled membrane distillation desalination modules. Appl. Energy 2017, 185, 2045–2057. [Google Scholar] [CrossRef]
  31. Katsandri, A. A theoretical analysis of a spacer filled flat plate membrane distillation modules using CFD: Part I: Velocity and shear stress analysis. Desalination 2017, 408, 145–165. [Google Scholar] [CrossRef]
  32. Tamburini, A.; Pitò, P.; Cipollina, A.; Micale, G.; Ciofalo, M. A Thermochromic Liquid Crystals Image Analysis technique to investigate temperature polarization in spacer-filled channels for Membrane Distillation. J. Membr. Sci. 2013, 447, 260–273. [Google Scholar] [CrossRef] [Green Version]
  33. Tamburini, A.; Cipollina, A.; Al-Sharif, S.; Albeirutty, M.; Gurreri, L.; Micale, G.; Ciofalo, M. Assessment of temperature polarization in membrane distillation channels by liquid crystal thermography. Desalin. Water Treat. 2014, 55, 37–41. [Google Scholar] [CrossRef]
  34. Tamburini, A.; Renda, M.; Cipollina, A.; Micale, G.; Ciofalo, M. Investigation of heat transfer in spacer-filled channels by experiments and direct numerical simulations. Int. J. Heat Mass Transf. 2016, 93, 1190–1205. [Google Scholar] [CrossRef]
  35. Albeirutty, M.; Turkmen, N.; Al-Sharif, S.; Bouguecha, S.; Malik, A.; Faruki, O.; Cipollina, A.; Ciofalo, M.; Micale, G. An experimental study for the characterization of fluid dynamics and heat transport within the spacer-filled channels of membrane distillation modules. Desalination 2018, 430, 136–146. [Google Scholar] [CrossRef]
  36. Phattaranawik, J.; Jiraratananon, R.; Fane, A. Heat transport and membrane distillation coefficients in direct contact membrane distillation. J. Membr. Sci. 2003, 212, 177–193. [Google Scholar] [CrossRef]
  37. Phattaranawik, J.; Jiraratananon, R.; Fane, A.G.; Halim, C. Mass flux enhancement using spacer filled channels in direct contact membrane distillation. J. Membr. Sci. 2001, 187, 193–201. [Google Scholar] [CrossRef]
  38. Gurreri, L.; Tamburini, A.; Cipollina, A.; Micale, G.; Ciofalo, M. CFD prediction of concentration polarization phenomena in spacer-filled channels for reverse electrodialysis. J. Membr. Sci. 2014, 468, 133–148. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Sketch of a generic overlapped spacer with the main geometric parameters. The red arrow indicates the flow direction.
Figure 1. Sketch of a generic overlapped spacer with the main geometric parameters. The red arrow indicates the flow direction.
Membranes 09 00091 g001
Figure 2. (a) Schematic representation and (b) picture of the experimental setup.
Figure 2. (a) Schematic representation and (b) picture of the experimental setup.
Membranes 09 00091 g002
Figure 3. Computational domain (unit cell): (a) θ = 30°; (b) θ = 60°; (c) θ = 90°.
Figure 3. Computational domain (unit cell): (a) θ = 30°; (b) θ = 60°; (c) θ = 90°.
Membranes 09 00091 g003
Figure 4. Mesh of a unit cell for case (c); one of the regions meshed with tetrahedra is highlighted.
Figure 4. Mesh of a unit cell for case (c); one of the regions meshed with tetrahedra is highlighted.
Membranes 09 00091 g004
Figure 5. Sketch of the configurations for which experimental and computational distributions of the normalized heat transfer coefficient are compared. The intrinsic angle θ between the filaments and the flow attack angle α are reported.
Figure 5. Sketch of the configurations for which experimental and computational distributions of the normalized heat transfer coefficient are compared. The intrinsic angle θ between the filaments and the flow attack angle α are reported.
Membranes 09 00091 g005
Figure 6. Configuration 1: distributions of the normalized heat transfer coefficient for increasing flow rates. (a,b) Q = 1 L/min. Re ≈ 205; (c,d) Q = 2 L/min. Re ≈ 410; (e,f) Q = 4 L/min. Re ≈ 820. Left column (a,c,e) experimental results; right column (b,d,f) computational fluid dynamics (CFD) predictions. Coordinates are in mm.
Figure 6. Configuration 1: distributions of the normalized heat transfer coefficient for increasing flow rates. (a,b) Q = 1 L/min. Re ≈ 205; (c,d) Q = 2 L/min. Re ≈ 410; (e,f) Q = 4 L/min. Re ≈ 820. Left column (a,c,e) experimental results; right column (b,d,f) computational fluid dynamics (CFD) predictions. Coordinates are in mm.
Membranes 09 00091 g006
Figure 7. Configuration 2: distributions of the normalized heat transfer coefficient for increasing flow rates. (a,b) Q = 1 L/min. Re ≈ 205; (c,d) Q = 2 L/min. Re ≈ 410; (e,f) Q = 4 L/min. Re ≈ 820. Left column (a,c,e) experimental results; right column (b,d,f) CFD predictions. Coordinates are in mm.
Figure 7. Configuration 2: distributions of the normalized heat transfer coefficient for increasing flow rates. (a,b) Q = 1 L/min. Re ≈ 205; (c,d) Q = 2 L/min. Re ≈ 410; (e,f) Q = 4 L/min. Re ≈ 820. Left column (a,c,e) experimental results; right column (b,d,f) CFD predictions. Coordinates are in mm.
Membranes 09 00091 g007
Figure 8. havg vs. Re for θ = 60°, α = 30° (configuration 1). Symbols denote the experimental results, while CFD results are represented by the thick solid line.
Figure 8. havg vs. Re for θ = 60°, α = 30° (configuration 1). Symbols denote the experimental results, while CFD results are represented by the thick solid line.
Membranes 09 00091 g008
Figure 9. havg vs. Re for θ = 90°, α = 60° (configuration 2). Symbols denote the experimental results, while CFD results are represented by the thick solid line.
Figure 9. havg vs. Re for θ = 90°, α = 60° (configuration 2). Symbols denote the experimental results, while CFD results are represented by the thick solid line.
Membranes 09 00091 g009
Figure 10. CFD results: havg vs. Re for θ = 30°,α = 15°; θ = 60°, α = 30°; and θ = 90°, α = 45°.
Figure 10. CFD results: havg vs. Re for θ = 30°,α = 15°; θ = 60°, α = 30°; and θ = 90°, α = 45°.
Membranes 09 00091 g010
Figure 11. CFD results: havg vs. Re for θ = 90°: (a) α ≤ 45° and (b) α ≥ 45°.
Figure 11. CFD results: havg vs. Re for θ = 90°: (a) α ≤ 45° and (b) α ≥ 45°.
Membranes 09 00091 g011
Table 1. Investigated values of geometric parameters of spacers.
Table 1. Investigated values of geometric parameters of spacers.
d (mm)P (mm)θ (°)α (°)
2103015
6030
900
15
30
45
60
75
90
Table 2. Investigated values of the flow rate and corresponding Reynolds numbers (Re).
Table 2. Investigated values of the flow rate and corresponding Reynolds numbers (Re).
Q (L/min)Re
1205
1.5307
2410
2.5512
3615
3.5717
4820

Share and Cite

MDPI and ACS Style

La Cerva, M.; Cipollina, A.; Ciofalo, M.; Albeirutty, M.; Turkmen, N.; Bouguecha, S.; Micale, G. CFD Investigation of Spacer-Filled Channels for Membrane Distillation. Membranes 2019, 9, 91. https://doi.org/10.3390/membranes9080091

AMA Style

La Cerva M, Cipollina A, Ciofalo M, Albeirutty M, Turkmen N, Bouguecha S, Micale G. CFD Investigation of Spacer-Filled Channels for Membrane Distillation. Membranes. 2019; 9(8):91. https://doi.org/10.3390/membranes9080091

Chicago/Turabian Style

La Cerva, Mariagiorgia, Andrea Cipollina, Michele Ciofalo, Mohammed Albeirutty, Nedim Turkmen, Salah Bouguecha, and Giorgio Micale. 2019. "CFD Investigation of Spacer-Filled Channels for Membrane Distillation" Membranes 9, no. 8: 91. https://doi.org/10.3390/membranes9080091

APA Style

La Cerva, M., Cipollina, A., Ciofalo, M., Albeirutty, M., Turkmen, N., Bouguecha, S., & Micale, G. (2019). CFD Investigation of Spacer-Filled Channels for Membrane Distillation. Membranes, 9(8), 91. https://doi.org/10.3390/membranes9080091

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