The Performance of an AGMD Unit Coupled with a Water Solar Collector the Experimental Validation for Desalination of Seawater
Received Date: Jan 24, 2018 / Accepted Date: Apr 18, 2018 / Published Date: Apr 25, 2018
Our work consists in presenting the results of an invention for a membrane distillation system coupled to an efficient and robust water solar collector which produces potable water with high quality and a small percentage of brackish discharge independent of salinity of the water source. The air gap membrane distillation (AGMD) model for the modules has been developed. This model is based on mathematical equations that describe the heat and mass transfer mechanisms of a single-stage AGMD process. It can simulate AGMD modules in regimes. The theoretical model was validated using in AGMD under different operating conditions and parameters. The predicted water vapor flux was compared to the flux measured at five different feed water temperatures, two different feed water salinities. The model was then used to study and analyze the parameters that have significant effect on scaling-up the AGMD process such as the effect of increasing the membrane length, and feed and coolant flow rates. The model was also used to analyze the maximum thermal efficiency of the AGMD process
Keywords: Desalination; Membrane distillation; Air gap membrane distillation; Solar energy; Coupling
A: area (m2); B: effective absorptivity; Cp: heat capacity (J/kg K); D: diffusion coefficient (m2/s) ; Fj : objective functionj; H: height (m); HVL: heat transfer rate of phase change (J/ m2s); h: heat transfer coefficient (W/m2K); I: intensity of solar radiation (W/m2); K: thermal conductivity (W/m K); Kc: proportional gain; Ku: ultimate gain; k: mass transfer coefficient (m/s); L: length (m); M: mass (kg); Mw: molecular weight of water kg/kmol); mf: fluid flowrate (kg/s); mf,c: fluid collector flowrate (kg/s); N: Mass flux (kmol/m2s); OP: controller output; P: pressure (Pa); Q: heat transfer rate (J/s); QN: sensible heat transfer rate (J/s); R: gas constant (J/kmol K); S: collector absorber surface area (m2); Si: stream number I; T: temperature (K); U: overall heat transfer coefficient of the heat exchanger (W/m2K); U’: overall heat loss coefficient between the collector absorber and the surroundings (W/m2K); W: width (m); DHvap: heat of vaporization (J/kmol); D: thickness (m); e: membrane porosity; m: viscosity (kg/m s); r: density (kg/m3); s: membrane tortuosity; sI: integral time (s); su: ultimate period (s); L: liquid; G: gas; Sat: saturated; AG: air gap ; A: Ambient ; Air: air ; Avg: average; c: solar collector; CL: cold liquid; CONL: condensate liquid; F: condensing film surface in the air gap; F: circulation fluid in the solar collector; fh: circulation fluid between the coil and the heat exchanger; fs: circulation fluid between the internal coil and the solar; Tf: fluid temperature (K); Tc: collector temperature (K); GM: gas in the membrane; HL: hot liquid; MD: membrane module; MEM: membrane; MET: metal; Nu: Nusselt number; Pr: Prandtl number; Re: Reynolds number; T: total; w: water; wa: watereair
Membrane distillation is a separation process that involves transport of water vapour through porous hydrophobic membranes. A variety of methods may be employed to impose this vapour pressure difference. In the present work, the air gap membrane distillation method (AGMD) is considered. The principal advantage of AGMD against other configurations of membrane distillation arises from the possibility to condensing the permeate vapours on a cold surface rather than directly in a cold liquid. In this configuration, the mass transfer steps involve movement within the liquid feed toward the membrane surface, evaporation at the membrane interface and transport of the vapour through the membrane pores and air gap prior to condensation. So, the separation mechanism of membrane distillation and its performance is based on vapour-liquid equilibrium. The driving force of the process is supplied by the vapour pressure difference caused by the existing temperature difference between the liquid-vapour interfaces. The benefits of membrane distillation over conventional distillation processes are found in the lower operating temperatures and pressures that reduce the equipment surface and costs, the compact modules and the possibility of overcoming corrosion problems by using plastic equipment. The lower operating temperatures can use available energy sources such as solar and geothermal energies or waste energy in industrial processes. That is why transmembrane evaporation applied to the desalination of seawater is a very promising technology to reduce the production costs of water from fresh water sources  (lakes, solar pond, etc.). However, at this moment little knowledge is available about the optimal design of such a transmembrane evaporation process. In fact, the coupling and non-linearity of the equations, the interaction between the flow and the thermal plume in the cavity and the complicated geometries involved have forced experimental solutions. So, the evident lack of data for thermal and dynamic interactions in natural convection within the cavity prompted the present study. This work presents a model for preliminary design calculations that is used to evaluate the effect of the main design parameters in the transmembrane evaporation module performance for water desalination. Developing mathematical models for prediction of membrane separation processes is an important tool in the field of membrane science and technology. The models play a meaningful role in simulation and optimization of membrane systems leading to efficient and economical designs of separation processes [1-3].This paper presents the modeling of phenomena (mass and heat transfers) and the experimental and simulated results of the parametric study of an installation of type PTFE on the pilot scale (500 L) is shown in Figure 1.
Materials And Methods
The different desalination methods
Figure 2 illustrates the desalination techniques classified in tow broad categories: membrane processes and thermal processes. The processes acting on the chemical bonds and processes are being performed by phase change. A method for separating salt water desalination in two parts: fresh water containing a low concentration of dissolved salts and concentrate brine. This process is energy-consuming various desalination techniques have been implemented over the years on the basis of the available energy [4,5].
Principles of MD: MD is a thermal process in which water vapor is transported through a hydrophobic porous membrane. The liquid phase to be treated must be kept in contact with one face of the membrane without penetrating the pores unless the trans membrane pressure is greater than the inlet pressure. The hydrophobicity of the membrane prevents the liquid entering the pores due to the surface tension. Thus, liquid/vapor interfaces are created in the vicinity of the pores.
Different configurations of MD
Direct contact membrane distillation (DCMD): A colder aqueous solution than the feed solution is maintained in direct contact with the permeate membrane (Figure 3). In this case, the transmembrane temperature difference induces a vapor pressure difference. Consequently, the volatile molecules evaporate at the liquid / vapor interface of the feed and condensate at the liquid / vapor interface of the colder permeate.
Air gap membrane distillation (AGMD): In Figure 4, there was an air knife is interposed between the membrane and the condensation surface. As a result, evaporated volatile molecules will pass through the pores of the membrane and the air knife to finally condense on a colder surface inside the membrane module.
Sweep gas membrane distillation (SGMD): In Figure 5 is presenting cold procedure side of membrane swept air in order to carry away the transferred steam. An inert cold gas sweeps the membrane on the side of the permeate carrying with it the vaporized molecules. The condensation then takes place outside the membrane module.
Vacuum membrane distillation (VMD): The Figure 6 shows a vacuum is applied to the permeate side by means of a vacuum pump. The applied vacuum pressure is less than the saturating vapor pressure of the volatile molecules to be separated. In this case, the condensation of the molecules takes place outside the membrane module.
Modeling of the membrane distillation unit
The 1-D model for individual AGMD is formulated and resolved numerically on the platform, where the membrane module can be developed via a graphical interface and the model components are linked and solved in an installation on the equations (Figure 7). The specifications of the components analyzed in this study are defined in Table 1 and by the following equations (4,5,10,11).
|Membrane area (m2)||10|
|Pore diameter (mm)||0.2|
|Thermal conductivity (W/m K)||0.173|
|Hot liquid channel thickness (mm)||0.77|
|Cold liquid channel thickness (mm)||0.77|
|Air gap channel thickness (mm)||0.43|
|Metal foil thickness (mm)||0.98|
|Thermal conductivity (W/m K)||398|
Table 1: Lists of membrane parameters and solar collector.
The 1-D AGMD mathematical model developed was for flat sheet membranes typically used in commercial MD pilot units. The model is based on dividing the AGMD module longitudinally into small elements. Within each element, different zones exist where significant mass and energy exchange occurs along the boundaries of these zones. As depicted in Figure 8, the zones, in order from left to right, are as follows:
• The hot fluid channel
• The polymeric membrane layer
• The air gap space
• The condensate film on the cooling plate
• The cooling plate sheet
• The cold fluid channel
In Figure 9, it is assumed that flow direction, x, is the same as that of the hot feed flow, and each small element is assumed to have a length of dx and a constant width W. Moreover, the mathematical model calculations were simplified according to the following assumptions:
• The system is at steady state condition.
• The hot and cold fluids are assumed to flow in the x direction only.
• The pressure inside the air gap is constant (no pressure drop along the air gap zone).
• The condensation on the cooling plate is film-wise and the thickness of the falling film inside the air gap is small in comparison with the width of the air gap.
• Within the air gap, there is no bulk velocity of the air-vapor mixture. Heat is transferred by conduction while mass is transferred through diffusion.
• Pure water vapor is only transported through membrane pores.
• There is no heat being exchanged with the surrounding.
The main components are the solar collector and the MD module, but you can add a heat exchanger with a slack storage if we will work day and night. The MD module is featured with a design of energy [4-6] recovery. The System’s performance is determined by the profiles of solar radiation, the design of individual system components, component integration, the business model and the control strategy. Modeling and control methods of solar desalination, Ben Bacha et al.  and Roca et al.  presented studies for a solar cycle system condensation and evaporation multiple hybrid fossil fuel powered solar distillation system, respectively. Both groups have developed reduced process for integrating their control algorithms proposed on the basis of linear technology control and feedback linearization technique, models respectively. Due to the intermittent and unpredictable nature of solar radiation, the steady state operation of the solar desalination process is not easy to achieve and the application of modern control algorithms is difficult. The purpose of this study is to develop a dynamic model including models for key components, and the latter to investigate the overall system optimization. The model is built on platform , which allows the analysis and control system design.
The hot fluid and hot fluid is against the current. For simplicity, seawater is used, along with hot and cold fluids in the model [10,11]. The purpose of this study is to investigate the optimization and control of the entire system, this simplification will not cause significant differences in the test results. The mass transfers resistance of the hot fluid side which is insignificant in a previous analysis  to the hollow fiber module. Transferring the mass flow is determined by taking into account the mass transfer resistance in the gap and membrane. However, the heat transfer resistance of the whole layers is taken into account. Mass and energy flows for all layers, including the hot fluid, the membrane gap, a cold fluid and metal sheet are illustrated in Figure 10, the model equations are summarized below.
The mass balances:
Balancing the energy:
The heat transfer coefficients for hot and cold fluid sides are estimated using the correlations reported by Schock and Miquel for the module flat plate wound membrane.
For the condensing heat transfer to the film, the following relationship is employed :
Results And Discussion
Simulation for membrane distillation unit
Model validation at different operating parameters: The first set of experiments was conducted to test the reproducibility and to determine the experimental errors. The measured water vapor flux at different feed water temperatures was repeatable and the variation in flux was a maximum of ±0.12 kg/m2•hr (2%). The mathematical model results were then validated against different experimental data. Figures 11a and 11b.
Using AGMD model, effects of heat and mass of the membrane and the air gap are studied by varying the heat transfer coefficients and mass simultaneously by a factor of 0.1 or 10. The results are shown in Figure 12.
Model validation at different operating parameters: The first set of experiments was conducted to test the reproducibility and to determine the experimental errors. The measured water vapor flux at different feed water temperatures was repeatable and the variation in flux was a maximum of ±0.12 kg/m2•hr (2%). The mathematical model results were then validated against different experimental data. Figure 13 shows a comparison between the predicted mass fluxes and the measured water vapor fluxes for a range of deionized feed water temperatures (40°C-80°C). The model predicted an exponential behavior of the AGMD flux as a function of feed water temperature. Such behavior is not only supported by our experimental data but also reported in published AGMD literature [14-17]. However, the validity of the mathematical model should not be judged based on predicting the trend of the process but also on how closely it predicts the absolute experimental data. Our current purpose of developing this model is to utilize it as a tool for further analyzing the AGMD process and for scale-up. Such a goal may require relaxed a criterion toward which we may judge the validity of our module. Nonetheless, the prediction of the model was within the range of the experimental error.
To validate the model further we replaced the deionized water (feed) with Red Sea water to see how the model predicts the water vapor flux for a seawater salinity of 4.2 wt%. The distillate conductivity was continuously measured to check for any pore wetting that may took place and the distillate conductivity was always below 20 μS. As shown in Figure 14 the predicted water vapor flux was also within the range of experimental error.
The effect of air gap width was also investigated. As shown in Figure 15 the model predicted decay in flux as the air gap increased. However, the model predictions for water vapor flux at different air gap widths was not as good as were the predictions for variations in feed temperature. Analysis of the results showed that the water vapor flux was very sensitive to the change of air gap width, especially when it is very small. A reduction in air gap width results in higher production capacity and higher errors. These errors are more significant when the air gap width is very small. Therefore, any small error in measuring the gap width (i.e., by 0.1 mm) will affect the water vapor flux significantly. The error of our measurements to the gap width was about ± 0.5 mm. Our investigation showed that this was due to the deformation of the parafilm tape used in sealing the module. Further experimental tests with a modified module are required in the future to better evaluate the model prediction at small air gap width.
Finally, the model was validated against experimental data using different membrane pore sizes. The model prediction was good enough (± 10%), although it didn’t predict well the data (15%) at feed temperature of 70°C for the 0.45 μm membrane (Figure 16). In this region the flux is increasing significantly as feed temperature is increased, so variations in the inlet temperature will have a larger effect on the measured flux compared to measurements at lower feed temperatures, and the error of 15% appears reasonable.
SEM images and membrane properties
The SEM observation indicated that the PDFE membrane (Table 1) studied is characterized by a foam structure and is symmetrical. Some differences in pore size occurred only on the outer surface (Figure 17). The maximum pore size observed on the inner membrane surface did not exceed a few microns.
Mass transfer in the MD process is diffuse. Therefore, the permeate flux is strongly affected by the wall thickness of the membrane and the pore diameter. The results of the study confirmed that the type of membrane used has a significant influence on the efficiency of the MD process (Figure 15).
Given the membranes with similar wall thickness, a higher flux was obtained, having a larger pore size in a membrane. Molecular diffusion and Knudsen influence mass transfer in the MD process, therefore, the permeate flux increases with an increase in pore diameter.
The mathematical model of a membrane module was developed from mass equations and heat transfer. The model calculations were based on the division of the AGMD module in different longitudinal areas. Normal for these areas of the module was cut into small cells. The mass and energy equation were applied to these slices and resolved by iterative procedures. The model has been validated under different conditions such as water temperatures, nutrition salinity, membrane pore sizes, and gap widths.
The model used in the analysis of complex and interdependent AGMD processes, which are essential elements for the extension of processes. The current analysis of the residence time in the indoor module AGMD is very important for the extension of processes because it has a direct effect on the process flow and its thermal efficiency. The flux decreases as the length of the membrane increases and increases with flow rate.
In addition, the efficiency of the process increases so that the surface of the membrane that increases the AGMD processes to operate at low temperature difference across the membrane. The effect of the type of flow is more visible with the variations of the flow rates than with the variation of the temperature of supply of the hot fluid. This effect disappears when the thickness of the air gap is high. The maximum permeate flux obtained was 7.4 kg / m² h with a fluid temperature of 80°C., air gap of 1.04 mm and hot and cold flow rates of 5 l / min. This work reveals that even at the low hot fluid supply temperature of 25°C, the AGMD configuration is capable of producing desalinated water. This aspect of the process may be useful in coupling with low temperature heat sources.
The membrane was then characterized by the SEM and FTIR technique to locate the presence of salt on the membrane.
Other tests will be carried out soon on the long-term use and the clogging effect of the membrane which is the main lock to the development of this kind of desalination technique.
- Lloyd DR, Lawson KW (1997) Membrane distillation. J Memb Sci 124:1-25.
- Bourawi MS, Ding Z, Ma R, Khayet M (2006) A framework for better understanding membrane distillation separation process. JMemb Sci 285:4-29.
- Kalogirou SA (2005) Seawater desalination using renewable energy sources. ProgrEnergy Combust Sci 31:242-281.
- Banat F, Jwaied N, Rommel M, Koschikowski J, Wieghaus M (2007) Performance evaluation of the “large SMADES” autonomous desalination solar-driven membrane distillation plant in Aqaba, Jordan. Desalination 217:17-28
- Banat F, Jwaied N, Rommel M, Koschikowski J, Wieghaus M(2006) Desalination by a“compact SMADES”autonomous solar-powered membrane distillation unit.Desalination 217:29-37.
- van Medevoort J, Jansen A, Hanemaaijer JH, Dotremont C, Nelemnas B, et al. (2008)Memstill: Seawater desalination a solution to water scarcity.
- Ben Bacha H, Bouzguenda M, Damak T, Abid MS, Maalej AY (2000) A study of a water desalination station using the SMCEC technique: production optimization. Renewable Energy 21:523-536.
- Roca L, Berenguel M, Yebra L, Alarcón DC (2008) Preliminary modeling and control studies in AQUASOL project. Desalination 222:466-473.
- Ben Bacha H, Damak T, Bouzguenda M, Maalej AY, Ben Dhia H (2003) A methodology to design and predict operation of a solar collector for a solar-powered desalination unit using the SMCEC principle. Desalination 156:305-313.
- Avlonitis S, Hanbury WT, Ben Boudinar M (1993) Spiral wound modules performance an analytical solution: part II. Desalination 89:227-246.
- Chang H, Liau JS, Ho CD, Wang WH(2009) Simulation of membrane distillation modules for desalination by developing user’s model on Aspen Plus platform. Desalination 249: 380-387.
- Holman JP (1989) Heat transfer. McGraw-Hill, New York.
- Koshchikowski J, Wieghaus M, Rommel M (2003) Solar thermal-driven desalination plants based on membrane distillation. Desalination 156:295-304.
- Schock G, Miquel A (1987) Mass transfer and pressure loss in spiral wound modules. Desalination 64:339-352.
- Chen CX, Yu LX, Dai TY (1996) Review of new membrane and membrane process. Technology of Water Treatment 22:307-313.
- Liu LY, Ding CW, Chang LJ (2008) Progress in membrane separation enhanced by ultrasound. Chemical Industry and Engineering Progess 27:32-37.
Citation: Boukhriss M, Bacha HB (2018) The Performance of an AGMD Unit Coupled with a Water Solar Collector the Experimental Validation for Desalination of Seawater. J Fundam Renewable Energy Appl 8: 256. DOI: 10.4172/2090-4541.1000257
Copyright: © 2018 Boukhriss M, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Select your language of interest to view the total content in your interested language
Share This Article
8th International Conference on Environment and Climate Change
November 22-23, 2018 Bucharest, Romania
April 01-02, 2019 Paris, France
- Total views: 382
- [From(publication date): 0-2018 - Oct 20, 2018]
- Breakdown by view type
- HTML page views: 347
- PDF downloads: 35