Open access peer-reviewed chapter

Computational Study of Liquid Film Condensation with the Presence of Non-Condensable Gas in a Vertical Tube

By Adil Charef, M’barek Feddaoui, Abderrahman Nait Alla and Monssif Najim

Submitted: November 20th 2017Reviewed: March 23rd 2018Published: November 5th 2018

DOI: 10.5772/intechopen.76753

Downloaded: 600

Abstract

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.

Keywords

  • thermal process
  • vapour-gas mixtures
  • condensation
  • heat and mass transfer
  • phase change

1. Introduction

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 [1] 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. [2] 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 [3] 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. [6] 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. [7] 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 [8] 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 [9] 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 [10] 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. [11] 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. [12] 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. [13] 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 [14] found from steam condensation that even a small air mass fraction decreases the heat transfer coefficient in seawater desalination plant. Al-Shammari et al. [15] 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. [16]. 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. [17] 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. [18] 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 L and radius R (Figure 1). The tube wall is subjected to a constant temperature. A mixture of water vapour and non-condensable gas enters the tube with a uniform velocity u0, vapour mass fraction w0, temperature T0 and pressure P0. The vapour condenses and forms a liquid film thickness as the mixture flowing downwards.

Figure 1.

Geometry of the problem.

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:

xρLuL+1rrrρLvL=0E1
  • Conservation of momentum:

xρLuLuL+1rrrρLvLuL=dpdx+1rrrμLuLr+ρLgE2
  • Conservation of energy:

xρLuLCp,LTL+1rrrρLvLCp,LTL=1rrrλLTLrE3

Similarly, the mass conservation, momentum, energy and diffusion equations for the gas phase are written as follows:

  • Conservation of mass:

xρGuG+1rrrρGvG=0E4
  • Conservation of momentum:

xρGuGuG+1rrrρGvGuG=dpdx+1rrrμGuGr+ρGgE5
  • Conservation of energy:

xρGuGCp,GTG+1rrrρGvGCp,GTG=1rrrλGTGrE6
  • Species diffusion equation:

xρGuGw+1rrrρGvGw=1rrrρGDwrE7

2.3. Boundary and interfacial conditions

The governing equations are subjected to the following boundary conditions:

  • At the tube inlet x=0

uG=u0;TG=T0;PG=P0;wG=w0E8
  • At the centre line of the tube r=0

vG=0;uGr=TGr=wr=0E9
  • At the wall of the tube r=R

uL=vL=0;TL=TWE10
  • At the interface vapour-liquid r=Rδx

Continuities of velocity and temperature:

uIx=uG,I=uL,I;TIx=TG,I=TL,IE11

Continuity of shear stress:

τI=μurL,I=μurG,IE12

Heat balance at the interface:

λLTLr=λGTGrJhfgE13

where hfg is the latent heat of condensation and Jis the mass flux at the interface J=ρGvI.

The radial velocity of water vapour-air mixture is calculated by considering that the interface is semipermeable [19] 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

vI=i=12DG,imwGir1i=12wGiE14

The governing Eqs. (1)(7) with interfacial conditions (8)–(13) are used to determine the field of variables uL,vL,TL,uG,vG,TG,w. 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:

m0L2π=RδxRrρudrL0xρGvIRδxdxE15
Rδ022ρ0u0=0RδxrρudrG+0xρGvIRδxdxE16

A dimensionless accumulated condensation is introduced to estimate the mass transfer along the tube:

mcd=2π0xρGvIRδxdxE17

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 r,xcoordinates are transformed into η,Xas follows:

η=RδxrRδx0rRδxE18
η=RδxrδxRδxrRE19
X=xLE20

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) [22]. 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 uL,vL,TL,uG,vG,TG,wand the two scalars dpdxand δx. The computational solution is advanced as follows:

  1. For any axial position x, guess an arbitrary values of dpdxand δx.

  2. Solve the finite difference forms of Eqs. (2)(3) and (5)(7) simultaneously for uL,TL,uG,TG,w.

  3. Numerically, integrate the continuities of Eqs. (1) and (5) to find vLand vG.

  4. The interfacial conditions of velocity, temperature, shear stress and heat balance are obtained from Eqs. (11)(13).

  5. Calculate the error of the liquid film mass balance ELittusing Eq. (15).

  6. The best approximation to the thickness of the liquid film is then obtained using the secant method [23]. Thus

    δxitt+1=δxittδxitt+δxitt1ELittELitt1ELittE21

The convergence criteria used is ELitt=105. Usually, six to seven iterations suffice to get converged solution.

  • Calculate the error in the gas flow balance EGittusing 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

    Err=maxϒi,jnϒi,jn1maxϒi,jn<105E22

    The solution for the actual axial position is complete. If not, repeat procedures (1) to (7), where ϒrepresents the variables uL,TL,uG,TG,w.

  • 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 [24], described by Anderson et al. [25]. To fully explain, we let H=dpdx. Due to an initial guesses for dpdx=dpdx, we calculate provisional velocities ujn+1and a mass flow rate of gas M˙jn+1. 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:

    ujn+1=ujn+1+ujn+1HΔHE23

    ΔHis the change in the gradient of the pressure required to satisfy the global mass flow constraint. In addition, we specify up,jn+1=ujn+1H. The difference equations are indeed differentiated with respect to the pressure gradient (H) to have difference equations for up,jn+1, 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 up,jn+1is resolved by Thomas algorithm. Furthermore, the boundary conditions on up,jn+1must be coherent with the boundary conditions of the velocity. At boundaries, where the velocity is specified, up,jn+1=0. The solution of up,jn+1is used to calculate ΔH, noting that to satisfy the constraint of the global mass flow, up,jn+1ΔHis the correction in velocity at each point. So, we can write

    M.in+1M.in=2πΔH0Rδxup,jn+1drE24

    where the integral is estimated using numerical means. The M.in+1in Eq. (24) is the known value specified in the initial conditions. The required value of ΔHis given by Eq. (24), whereas the correct values of the velocity ujn+1may be determined from Eq. (23). Besides, the continuity equation then allows to calculate vjn+1.

    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 NI × (NJ + NL) = 131 × (81 + 31) is chosen because it gives results close enough to those of the thin grid and sufficiently accurate to describe the heat and mass transfer. Note that NI is the total grid points in the axial direction, NJ is the total grid points in the radial direction at the gas region, and NL is the total grid points in the radial direction at the liquid region.

    Figure 2.

    Comparison of sensible heat Nusselt number.

    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. [17] 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. [17] as shown in Figure 3a, b, which illustrates the vapour mass fraction evolution and mixture temperature, respectively.

    Figure 3.

    Comparison with numerical study of Hassaninejadfarahani et al. [17] for (a) dimensionless mass fraction at the tube exit, and (b) dimensionless mixture temperature.

    The computations have been also compared with experimental results of Lebedev et al. [2]. It is important to indicate that Lebedev et al. [2] examined the simultaneous heat and mass transfer during humid air condensation in a vertical duct. So, to obtain the case of Lebedev et al. [2], equivalent hydraulic diameter de is chosen [7, 26]. Figure 4 is a plot of local condensate heat transfer coefficient compared with the study of Lebedev et al. [2] for T0 = 60°C, P0 = 1 atm, L = 0.6 m, TW = 5°C and de = 0.02 m. A good agreement between our computations and the experiment curves is found with a maximum relative error of 4.7% for both curves (u0 = 1.4 m/s and u0 = 0.7 m/s).

    Figure 4.

    Comparison with experimental data of Lebedev et al. [2] for local condensate heat transfer coefficient.

    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 T0=70°C, inlet pressure P0=1 atm and inlet Reynolds number is fixed at Re0 = 2000. The range of each parameter for this study is listed in Table 1. At first, the air is used as non-condensable gas.

    Tube length (L (m))3.0, 4.5, 6.0
    Tube radius (R(m))0.008, 0.01, 0.012
    Inlet vapour mass fraction (w0)0.05, 0.125, 0.2
    wall temperature (TW(°C))5, 20, 35
    Non-condensable gasOxygen, air, nitrogen

    Table 1.

    The ranges of the physical parameters.

    Figures 57 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).

    Figure 5.

    Distributions of axial velocity profile in both the liquid and vapour phases at different tube sections.

    Figure 6.

    Distributions of axial temperature profile in both the liquid and vapour phases at different tube sections.

    Figure 7.

    Distributions of axial vapour mass fraction profile in gas phases at different tube sections.

    5. Effect of the tube geometry (length L and 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 δxincreases with increasing the tube length L. The results also indicate that, for a high tube length, the mechanism of the condensation is important when the distance is less than (X = 0.8), especially for L = 4.5 m and 6.0 m, because for a fixed Reynolds gas number (the inlet velocity is fixed too), the condensed vapour decreases with the increase of tube length. This means that near to the tube exit, the condensation process becomes almost unimportant. It is also observed that the condensing mass flux reduces from the inlet to the outlet, particularly for the two large values of L approaching zero, because the total amount of the water vapour is transferred to liquid film and remains only the air at the vapour-liquid interface.

    Figure 8.

    Effect of the tube length on the variation of the liquid film thickness and the condensing mass flux at the interface along the tube.

    The effect of changing the tube radius R is shown in Figure 9. From the curves in this figure, it is interesting to observe that for a fixed value of T0, w0 and Re0, a smaller tube radius corresponds to a thicker liquid film thickness. This trend is true for every position X from the inlet to the outlet of the tube. These results are directly related to the velocity of the gas mixture at the inlet. Obviously, for a fixed Reynolds number, δxincreases according to the velocity increases with a weak tube radius. This implies that a higher inlet velocity tends to move the air away from the interface and thus maintains its lower fraction, leading to the increase of the heat transfer coefficient, which improves the condensation process. It is also found that the condensing mass flux increases with the tube radius only near the inlet when X < 4 due to a high interfacial shear stress. This tendency is reversed as the gas mixture progresses along the tube.

    Figure 9.

    Effect of the tube radius on the variation of the liquid film thickness and the condensing mass flux at the interface.

    6. Effect of water vapour mass fraction w0

    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 w0 considerably improves the condensation mechanism along the tube. Indeed, for a constant T0, an increase of w0 affects the thermo-physical properties of the gas mixture at the inlet, which leads to an augmentation of the vapour partial pressure and the temperature at the vapour-liquid interface. Consequently, the condensing mass flux at the interface J“ increases significantly with w0 leading to an increase in the rate of condensation, which improves the thickness of the liquid film. On the other hand, a small amount of w0 (inversely proportional to the mass fraction of the non-condensable gas) causes a remarkable reduction of the condensed mass flux rate and the axial variation of the thickness of the film along the tube. This is due to the presence of air, which plays the role of thermal and mass transfer resistance at the vapour-liquid interface.

    Figure 10.

    Effect of the inlet vapour mass fraction on the variation of the liquid film thickness and the condensing mass flux at the interface.

    7. Effect of the wall temperature TW

    The effect of the wall temperature TW on the liquid film thickness and the condensing mass flux is presented in Figure 11. It is noted that the thickness of the liquid film varies inversely with the imposed temperature of the wall. Clearly, there is a significant growth of the liquid film thickness when TW reduces from 35°C to 5°C because the amount of the condensed vapour is enhanced along the tube by the increase of the heat transfer and hence the thickness of the liquid film becomes thicker. From this figure, it is also observed that the condensing mass flux is important for a large value of the inlet-to-wall temperature difference and then decreases along the tube. The results indicate that J“ decreases as the gas flow progresses along the tube because the condensation is accompanied by a diminution in the temperature of the vapour phase and the heat transferred by the latent mode during the condensation to the liquid film. The effect of wall temperature on the accumulated condensation mcd is illustrated in Figure 12. It can be seen that mcd becomes important for low temperatures. When the temperature difference (T0 - TW) increases, the heat transfer increases, and consequently the density of the condensed flux increases. This explains why the condensation process is favoured for a higher-temperature difference.

    Figure 11.

    Effect of TW on the evolution of the liquid film thickness and the condensing mass flux at the interface.

    Figure 12.

    Effect of TW on the accumulated condensation along the tube.

    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 δ and the condensing mass flux at the interface J“ determine the condensation efficiency, Figure 13 illustrates the variation of δ and J“ along the tube. For a fixed w0, P0 and Re0, an increase in the molar mass of the gas leads to a growth in the density of the gas mixture, the temperature of the gas mixture, as well as a decrease in the saturation concentration especially for water-oxygen mixture. This causes a strong vapour concentration gradient and condensing mass flux at the vapour-liquid interface. The results obtained show that the evolution of the film thickness and condensing mass flux in the water-oxygen mixture are significantly greater than those of the other mixtures. In addition, the condensing mass flux decreases along the tube and tends to a lower value especially near to the tube exit, which means the end of the condensation process.

    Figure 13.

    Effect of the type of non-condensable gas on the evolution of the liquid film thickness and the condensing mass flux at the interface along the tube.

    9. Conclusion

    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:

    1. 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.

    2. A small amount of non-condensable gas improves the heat and mass exchanges.

    3. Decreasing the wall temperature enhances the liquid film thickness and the accumulated condensation.

    4. The non-condensable gas type has a great effect on the condensation process.

    Nomenclature

    Cpspecific heat (J. kg−1. K−1)
    Ddiffusion coefficient (m2. s−1)
    ddiameter of the tube (m)
    Rradius of the tube (d/2) (m)
    Ltube length (m)
    ggravitational acceleration (m. s−2)
    Nussensible Nusselt number
    mcdaccumulated condensation
    Patmospheric pressure (Pa)
    Ttemperature (°C)
    uaxial velocity (m. s−1)
    vradial velocity (m. s−1)
    wmass fraction of vapour
    J”mass flux at the interface (kg.m−2.s−1)
    hfglatent heat of condensation (J. kg−1)
    rradial coordinate (m)
    δliquid film thickness (m)
    λthermal conductivity (w.k−1. m−1)
    μdynamic viscosity (kg. m−1. s−1)
    ρdensity (kg. m−3)
    τshear stress
    ɸrelative humidity (%)
    Lreferring to the liquid
    Greferring to the gas mixture
    Iinterface
    0condition at inlet of the tube
    Wcondition at wall of the tube
    areferring to the air
    xaxial

    © 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

    How to cite and reference

    Link to this chapter Copy to clipboard

    Cite this chapter Copy to clipboard

    Adil Charef, M’barek Feddaoui, Abderrahman Nait Alla and Monssif Najim (November 5th 2018). Computational Study of Liquid Film Condensation with the Presence of Non-Condensable Gas in a Vertical Tube, Desalination and Water Treatment, Murat Eyvaz and Ebubekir Yüksel, IntechOpen, DOI: 10.5772/intechopen.76753. Available from:

    chapter statistics

    600total chapter downloads

    1Crossref citations

    More statistics for editors and authors

    Login to your personal dashboard for more detailed statistics on your publications.

    Access personal reporting

    Related Content

    This Book

    Next chapter

    Distilled Water Production by Vacuum Heat Pump

    By Liu Bin, Cai Ling, Li Tianyin and Sajid Muhammad

    Related Book

    First chapter

    Hydrogen Generation by Water Electrolysis

    By Youssef Naimi and Amal Antar

    We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

    More About Us