The ranges of the physical parameters.
The main objective of this chapter is to study the liquid film condensation in a thermal desalination process, which is based on the phase change phenomenon. The external tube wall is subjected to a constant temperature. The set of the non-linear and coupled equations expressing the conservation of mass, momentum and energy in the liquid and gas mixtures is solved numerically. An implicit finite difference method is employed to solve the coupled governing equations for liquid film and gas flow together with the interfacial matching conditions. Results include radial direction profiles of axial velocity, temperature and vapour mass fraction, as well as axial variation of the liquid film thickness. Additionally, the effects of varying the inlet conditions on the phase change phenomena are examined. It was found that increasing the inlet-to-wall temperature difference improves the condensate film thickness. Decreasing the radius of the tube increased the condensation process. Additionally, non-condensable gas is a decisive factor in reducing the efficiency of the heat and mass exchanges. Overall, these parameters are relevant factors to improve the effectiveness of the thermal desalination units.
- thermal process
- vapour-gas mixtures
- heat and mass transfer
- phase change
The demand of fresh water supply is increasing due to the economic development and the fast population growth. With limited resources of fresh water, desalination of seawater and brackish water offers the potential to encounter the increasing water demands around the world. Generally, the reverse osmosis has about the great part of the market share in the world compared to thermal desalination technologies. Consequently, the necessity to improve the thermal processes, which are based on the phase change phenomenon of evaporation and condensation, continues to receive a high interest. Condensation on the cooling surfaces is a phenomenon of major significance in the chemical industries, refrigeration, heat exchangers and desalination units, including thermal desalination.
The mechanism of condensation can be classified by various ways: geometric configurations like tube, channel, internal, external, horizontal or vertical; species of fluid such as steam, refrigerant or mixture with the presence of non-condensable gas; condensing phenomena as filmwise, dropwise or fog; and flow regime like laminar and turbulent. Since the first analysis of Nusselt  for film condensation on a vertical plate, a numerous number of studies have been done on improving film condensation modelling and to contribute to the comprehension of this complex phenomenon. Lebedev et al.  performed experimentally a combined study of heat and mass transfer from water vapour on a flat plate. They observed an enhancement of the condensation heat transfer with the increase of the inlet relative humidity. Dobran and Thorsen  studied the laminar filmwise condensation of a saturated vapour inside a vertical tube. They found that the mechanism of condensation is governed by ratio of vapour to liquid viscosity, Froude number to Reynolds number ratio, subcooling number and Prandtl number of liquid. Siow et al. [4, 5] presented a numerical study of the laminar film condensation with the presence of non-condensable gas in horizontal and then in vertical channels. They analysed the effect of the inlet Reynolds number, the inlet pressure and the inlet-to-wall temperature difference on the condensation mechanism. They studied also the liquid film condensation from steam-air mixtures inside a vertical channel. Results indicate that a higher concentration of non-condensable gas caused substantial reduces in the local Nusselt number, the pressure gradient and the film thickness. Belhadj et al.  conducted a numerical analysis to improve the condensation process of water vapour inside a vertical channel. Their results show that the phenomenon of phase change is sensitive to the inlet temperature of liquid film. For different values of the system parameters at the inlet of the tube, Dharma et al.  estimated from a numerical study the local and average values of Nusselt number, the pressure drop, the condensate Reynolds number and the gas-liquid interface temperature. Lee and Kim  carried out experimental and analytical studies to analyse the effect of the non-condensable gas (nitrogen) on the condensation of water vapour along a vertical tube with a small diameter. The experimental results demonstrate that the heat transfer coefficients become important with a high inlet vapour flow and the reduction of mass fraction of nitrogen. In addition, the authors developed a new correlation to evaluate the heat transfer coefficient regardless the diameter of the condenser tube. Nebuloni and Thome  developed a numerical and theoretical model to predict the laminar film condensation inside various channel shapes. They showed that the channel shape strongly affects the overall thermal performance. Chantana and Kumar  investigated experimentally and theoretically the heat transfer characteristics of steam-air during condensation inside a vertical tube. They observed that a higher Reynolds number and mass fraction of vapour improve the process of condensation. Dahikar et al.  conducted an experimental and CFD studies in the case of the film condensation with downward steam inside a vertical pipe. They found that a larger interfacial shear affects the momentum transfer because of the great velocity gradient especially at the gas-liquid interface. Merouani et al.  presented a numerical analysis during the condensation of steam-gas mixtures between two coaxial cylinders. They observed that a higher vapour concentration at the inlet and molar mass of the non-condensable gas increases the heat flux at the inner wall. Qiujie et al.  presented a numerical study in the case of steam-air condensation on isothermal vertical plate by using volume of fluid (VOF) method. Their results indicate that the mass fraction variation of the non-condensable gas directly affects the liquid film condensation and then influences the heat transfer.
In the thermal desalination unit, the condenser is used for producing freshwater from the saline water sources. In fact, in order to enhance the condensation process with the presence of non-condensable gas in thermal seawater desalination processes, many studies have been conducted. Semiat and Galperin  found from steam condensation that even a small air mass fraction decreases the heat transfer coefficient in seawater desalination plant. Al-Shammari et al.  have shown from an experimental study that non-condensable gas has a negative effect on the heat transfer. An experimental study on the role of non-condensable gases in the condensation of steam inside slightly inclined tubes was presented by Caruso et al. . The experiments were carried out under the following conditions: inner diameter of the tube 12.6 and 26.8 mm, inclination of the tube 7°, mass fraction of the non-condensable gas between 5 and 42%, and vapour saturated at atmospheric pressure. Due to the gravity, the condensate is collected mainly in the lower part of the tube. They also developed a correlation to calculate the local condensation heat transfer coefficient. Hassaninejadfarahani et al.  investigated numerically a liquid film condensation with high amount of non-condensable gas inside vertical tube. They examined the effects of varying the inlet air mass fraction, the inlet relative humidity, the inlet Reynolds number and the radius of the tube on the simultaneous heat and mass transfer during condensation. Recently, Charef et al.  investigated the condensation process of water vapour-air into liquid film inside a vertical tube under two different boundary conditions: imposed temperature and imposed heat flux. The results indicated a better condensation process under imposed heat flux. It was found that the presence of non-condensable gas affects negatively the system efficiency.
The purpose of this study is to numerically develop and investigate the problem of water vapour condensation in the presence of non-condensable gas in a vertical tube. In order to improve the effectiveness of the steam condensation in desalination process, special attention is addressed to examine the effects of the tube geometry and the inlet conditions on the condensation process. In the following, we present the studied problem, the numerical method and the main results.
2. Mathematical model
2.1. Physical model and assumptions
The geometry under consideration is a vertical tube with length
For the mathematical formulation of the problem, it has been assumed that the gas flow is laminar, incompressible and two-dimensional. The vapour and liquid phases are in thermodynamic equilibrium at the interface. In addition, viscous dissipation and other secondary effects are negligible, and the humid air is assumed to be a perfect gas.
2.2. Mathematical formulation
With respect to the mentioned assumptions, the governing equations for the conservation of mass, momentum and energy, respectively, in the liquid region are written as
Conservation of mass:
Conservation of momentum:
Conservation of energy:
Similarly, the mass conservation, momentum, energy and diffusion equations for the gas phase are written as follows:
Conservation of mass:
Conservation of momentum:
Conservation of energy:
Species diffusion equation:
2.3. Boundary and interfacial conditions
The governing equations are subjected to the following boundary conditions:
At the tube inlet
At the centre line of the tube
At the wall of the tube
At the interface vapour-liquid
Continuities of velocity and temperature:
Continuity of shear stress:
Heat balance at the interface:
where hfg is the latent heat of condensation and is the mass flux at the interface .
The radial velocity of water vapour-air mixture is calculated by considering that the interface is semipermeable  and that the solubility of air in the liquid film is negligibly small, which implies that the air velocity in the radial direction is zero at the interface. The velocity of the steam-air mixture at the interface can be written as
The governing Eqs. (1)–(7) with interfacial conditions (8)–(13) are used to determine the field of variables . To complete the mathematical model, two equations are used. At every axial location, the overall mass balance in the liquid phase and the gas flow should be satisfied:
A dimensionless accumulated condensation is introduced to estimate the mass transfer along the tube:
A transformation of coordinates was performed to ensure that the computational grid would clearly define the gas-liquid interface at each station along the tube. The coordinates are transformed into as follows:
The pure component data (in previous formulations) is approached by polynomials in terms of mass fraction and temperature. For more information, the thermo-physical properties are available in [20, 21].
3. Numerical solution method
The set of non-linear governing equations are discretized using a finite difference numerical scheme. The radial diffusion and the axial convection terms are approximated by the central and the backward differences, respectively. Hence, we arrange the system of discretized algebraic equations coupled with the boundary conditions into a matrix. Finally, the matrix resolution is carried out using the tri-diagonal matrix algorithm (TDMA) . Besides that, a special care was made to ensure accuracy of the numerical computation, by generating a non-uniform grid in both directions. Accordingly, the grid is refined at the interface. In fact, it is important to note that as the liquid goes to the outlet, the film thickness varies along the tube. For that reason, during the downstream marching at each iteration, our finite difference computational grid deals with the variation of the liquid and gas computational domain.
3.1. Marching procedure
A set of non-linear algebraic equations is realized for and the two scalars and . The computational solution is advanced as follows:
For any axial position x, guess an arbitrary values of and .
Calculate the error of the liquid film mass balance using Eq. (15).
The best approximation to the thickness of the liquid film is then obtained using the secant method . ThusE21
The convergence criteria used is . Usually, six to seven iterations suffice to get converged solution.
Calculate the error in the gas flow balance using Eq. (16).
Check the satisfaction of the convergence of velocity, temperature and species concentrations. If the relative error between two consecutive iterations is small enough, that is
The solution for the actual axial position is complete. If not, repeat procedures (1) to (7), where represents the variables .
3.2. Velocity and pressure coupling
Owing to satisfy the global mass flow constraint, the pressure correction gradient and axial velocity profile are performed applying a method proposed by Raithby and Schneider , described by Anderson et al. . To fully explain, we let . Due to an initial guesses for , we calculate provisional velocities and a mass flow rate of gas . Because of the linearity of the equation of momentum with frozen coefficients, the correct velocity at each point from an application of Newton’s method is as follow:
is the change in the gradient of the pressure required to satisfy the global mass flow constraint. In addition, we specify . The difference equations are indeed differentiated with respect to the pressure gradient (H) to have difference equations for , which have a tridiagonal form. The coefficients for the unknowns in these equations will be the same as for the original implicit difference equations. The system of algebraic equations for is resolved by Thomas algorithm. Furthermore, the boundary conditions on must be coherent with the boundary conditions of the velocity. At boundaries, where the velocity is specified, . The solution of is used to calculate , noting that to satisfy the constraint of the global mass flow, is the correction in velocity at each point. So, we can write
where the integral is estimated using numerical means. The in Eq. (24) is the known value specified in the initial conditions. The required value of is given by Eq. (24), whereas the correct values of the velocity may be determined from Eq. (23). Besides, the continuity equation then allows to calculate .
3.3. Mesh stability and validation of the numerical model
To validate the grid independency of results and to avoid convergence problems due to the use of thin grids, it is helpful to choose an optimum solution between computational time and result precision. Several grid sizes have been examined to ensure that the results are grid independent. Figure 2 shows that in all grid arrangements, the difference in local Nusselt number of sensible heat is always less than 3%. The grid with
In order to check the accuracy and validity of the numerical method, the obtained results were first compared to those reported by Hassaninejadfarahani et al.  in the case of laminar condensation of a steam and non-condensable gas in a vertical tube, in which the tube wall is maintained at a constant temperature. A good agreement was found between the current computational study and the results provided by Hassaninejadfarahani et al.  as shown in Figure 3a, b, which illustrates the vapour mass fraction evolution and mixture temperature, respectively.
The computations have been also compared with experimental results of Lebedev et al. . It is important to indicate that Lebedev et al.  examined the simultaneous heat and mass transfer during humid air condensation in a vertical duct. So, to obtain the case of Lebedev et al. , equivalent hydraulic diameter
4. Distribution of axial velocity, temperature and mass fraction profiles along the vertical tube
This chapter investigates the process of the liquid film condensation from the water vapour and non-condensable gas mixtures inside a vertical tube. The results of this study have been obtained for the case of inlet gas temperature
|Tube length (||3.0, 4.5, 6.0|
|Tube radius ||0.008, 0.01, 0.012|
|Inlet vapour mass fraction (||0.05, 0.125, 0.2|
|wall temperature (||5, 20, 35|
|Non-condensable gas||Oxygen, air, nitrogen|
Figures 5–7 Illustrate the profiles of velocity, temperature and the mass fraction of water vapour at different axial locations of the tube. From the distribution of velocity in Figure 5, it is observed that the variation of the velocity in the gas mixture is higher than that in the liquid region. Moreover, as the gas flow progresses along the tube, the velocity in the mixture decreases, while the velocity in the liquid film slightly rises. This behaviour is due to the mass transfer from the mixture to the liquid film. In fact, when the gas mixture loses the mass, it loses velocity too, however, the liquid film gaining mass as well as acceleration. Figure 6 Presents the evolution of the temperature profiles in both mixture and liquid phases at different tube sections. It can be seen that in the liquid phase, the temperature profiles are close to the temperature of the wall and nearly linear. This indicates that the interface temperature decreases from the inlet to the tube exit, resulting in the reduction of the heat transfer across the condensate film. Also, the slopes of the temperature of the mixture are decreasing along the tube owing to the absorbed energy transferred from the gas flow to the liquid film. The distributions of the mass fraction of water vapour in the gas region are illustrated in Figure 7. It is interesting to observe that the vapour mass fraction w0 decreases from the entrance to the tube exit, which implies that the condensation rate is decreasing along the tube. Consequently, w0 is reduced from the Centre line (η = 1) to the liquid-vapour interface (η = 0).
5. Effect of the tube geometry (length
Land radius R)
In the desalination units, the tube geometry of which the water vapour condenses (whether it is the length or the radius) contributes positively to the improvement of the condensation process if they are well dimensioned.
To reveal their impacts, we first examined the impact of the length of the tube on the liquid film thickness and the condensing mass flux at the interface along the tube. Figure 8 shows the influence of the tube length on the thickness of the film and condensing mass flux at the interface. It is noted that increases with increasing the tube length
The effect of changing the tube radius
6. Effect of water vapour mass fraction
In the majority of thermal desalination units, the water vapour that does not condense at the first effect, with all the non-condensable gases content, is transferred to the second effect, and this produces gas accumulation up to inadmissible concentrations. These gases cause a reduction in the performance of the system.
This result is confirmed in Figure 10, which shows that the thickness of the film increases considerably from the inlet to the exit of the tube. It is also observed that the increase in the mass fraction of water vapour
7. Effect of the wall temperature
The effect of the wall temperature
8. Effect of the non-condensable gas type
In thermal desalination units, when the water vapour condenses, the presence of a non-condensable gas hinders this phenomenon. The accumulation of non-condensable components at the vapour-liquid interface plays the role of an obstacle for heat and mass transfer. This causes a reduction in the efficiency of the system and therefore an increase in costs in most desalination units using phase change.
In this section, we analyse the influence of the non-condensable gas type during water vapour condensation. We considered mixed mixtures of water-oxygen, water-air and water- nitrogen. The molar mass of oxygen, air and nitrogen are equal to 31.99, 28.95 and 28.01 g/mol, respectively. Since the liquid film thickness
A numerical analysis has been carried out to investigate the liquid film condensation of the water vapour with the presence of non-condensable gas inside a vertical tube. The main conclusions drawn from this study are as follows:
The efficiency of the system is enhanced by increasing the tube length and decreasing the radius, which allows condensing the maximum of the water vapour.
A small amount of non-condensable gas improves the heat and mass exchanges.
Decreasing the wall temperature enhances the liquid film thickness and the accumulated condensation.
The non-condensable gas type has a great effect on the condensation process.
|Cp||specific heat (J. kg−1. K−1)|
|D||diffusion coefficient (m2. s−1)|
|d||diameter of the tube (m)|
|R||radius of the tube (d/2) (m)|
|L||tube length (m)|
|g||gravitational acceleration (m. s−2)|
|Nus||sensible Nusselt number|
|P||atmospheric pressure (Pa)|
|u||axial velocity (m. s−1)|
|v||radial velocity (m. s−1)|
|w||mass fraction of vapour|
|J”||mass flux at the interface (kg.m−2.s−1)|
|hfg||latent heat of condensation (J. kg−1)|
|r||radial coordinate (m)|
|δ||liquid film thickness (m)|
|λ||thermal conductivity (w.k−1. m−1)|
|μ||dynamic viscosity (kg. m−1. s−1)|
|ρ||density (kg. m−3)|
|ɸ||relative humidity (%)|
|L||referring to the liquid|
|G||referring to the gas mixture|
|0||condition at inlet of the tube|
|W||condition at wall of the tube|
|a||referring to the air|