Filmwise condensation in a stagnant pure vapor reservoir
From ThermalFluidsPedia
Contents 
Laminar Flow Regime
The classical analysis of laminar film condensation on a vertical or inclined wall was performed by Nusselt (1916). The physical conditions of laminar film condensation have been shown in Fig. 1. As is the case in any heat transfer analysis, the final goal is to obtain the heat transfer coefficient and the corresponding Nusselt number for the heat transfer device under consideration. Therefore, the objective of this section is to find the heat transfer coefficient and the Nusselt number for the laminar flow regime in film condensation on a vertical surface. The classical Nusselt analysis requires many assumptions in order to achieve a closedform solution of the boundary layer type momentum equations. These assumptions include the following:
1.The flow is laminar.
2.Fluid properties are constant.
3.Subcooling of the liquid is negligible in the energy balance, i.e., all condensation occurs at the saturation temperature corresponding to the pressure in the liquid film near the wall.
4.Inertia and convection effects are negligible in the boundary layer momentum and energy equations, respectively.
5.The vapor is assumed stagnant, and therefore, shear stress is considered to be negligible at the liquidvapor interface.
6.The liquidvapor interface is smooth, i.e., condensate film is laminar and not in the wavy or turbulent stages.
Since the vapor phase is stationary, the governing equations for the vapor phase are no longer needed. The continuity and momentum equations for the liquid phase ard:
The pressure gradient in the momentum equation may be approximated in terms of conditions outside the liquid film, as required by boundary layer approximation. It follows that pressure in the liquid film satisfies

Substituting eq. (1) into the momentum equation, one obtains

The lefthand side term represents the inertia effects in the slender film region, while the righthand terms are the effects due to friction and the sinking effect. As stated by Assumption #4 of the above list, the inertia term is negligible in this analysis. Therefore, the above equation simplifies to

Integrating this equation twice with respect to y and using the boundary conditions of nonslip at the wall (u=0 at y=0) and zero shear at the liquidvapor interface (δu/δy=0 at y = δ, Assumption #5) yields the following:

It can be seen from this equation that the downward velocity component is directly dependent on both the x and y coordinates due to the varying (increasingly thick in the xdirection) film thickness. The mass flow rate per unit width of surface, Γ, of this liquid film at any point can be found by integrating the above velocity profile across the liquid film thickness and multiplying by the liquid density, i.e.,

where it can now be seen that the mass flow rate is a function of the x coordinate.
Recalling that the heat transfer across the liquid film is by conduction only, and assuming that no subcooling exists at the liquidvapor interface (Assumption #3), Fourier’s Law can be used to obtain the heat flux across the film thickness:

where T_{sat} and T_{w} are the saturation temperature and wall temperature, respectively.
The heat transfer rate per unit width for the control volume shown in Fig. 1 is given as

where ΔT = T_{sat} − T_{w}. Since no subcooling in the liquid film exists, the latent heat effects of condensation dominate the process. Therefore, it can be said that

where dΓ is found by differentiating the expression for mass flow rate per unit surface, eq. (5). It is found to be

This expression and the expression for heat flow across the control volume are then substituted into the energy equation, eq. (8), and the following is found:

Finally, δ can be found by integrating the above and using the boundary condition that δ = 0 at x = 0:

The local heat transfer coefficient is found to be

The local Nusselt number, Nux, of the laminar film condensation process on a vertical plate is found to be

It is also desirable to obtain the average heat transfer coefficient and Nusselt number for a plate of length L. The mean heat transfer coefficient can be defined as

Substituting the expression for , eq. (12), into eq. (14) and integrating, the following expressions are obtained:

The local and average heat transfer coefficient can also be nondimensionalized in term of Reynolds number Re_{δ}, defined as

Substituting eq. (5) into eq. (16), the Reynolds number for laminar film condensation becomes

Substituting eqs. (16) and (11) into eqs. (13) and (15), the following nondimensional correlations are obtained:

where the bracket term on the lefthand side, together with its exponent (1/3), is the characteristic length.
The above Nusselt analysis assumes no subcooling in the liquid condensate. This assumption can be removed by making an energy balance at the interface that takes subcooling of the liquid into account, as follows:

The final term of the righthand side takes into account a temperature gradient across the liquid condensate film. Substituting the velocity profile in eq. (4) and using a linear temperature profile,

to evaluate eq. (19), the energy balance can be written as

where

Equation (22) is identical to the energy balance in eqs. (7) and (8) except that has been replaced by . Rohsenow (1956) improved this analysis even more by including the effects of convection in the liquid along with liquid subcooling, and thereby developing the following:

Finally, if it is desirable to include the effects of vapor superheat in the above analysis, the latent heat obtained by eqs. (22) or (23) can be further modified by adding c_{pv}(T_{v} − T_{sat}).
The above analysis can also be applied for condensation outside a vertical tube if δ/D = 1, where D is the diameter of the tube. For an inclined wall with an inclination angle of θ (the angle between the wall and the vertical direction), the component of the gravitational acceleration along the inclined wall is gcosθ. It follows from eqs. (12) and (15) that and .
Wavy Condensate Regime
Wavy flows of thin liquid films have higher heat transfer coefficients than smooth thin films, due to mixing action and an increase in the interfacial surface area. It has been shown by experiments that even during well defined laminar film flow, the film surface (liquidvapor interface) can be wavy. These waves can occur on film flow over a rough wall or a polished wall. This wave formation can even occur when the vapor reservoir is stagnant; however, the waves can be more pronounced when the vapor reservoir has an average velocity, which leads to a higher shear stress at the liquidvapor interface. These small disturbances can amplify under specific circumstances, i.e., when the film reaches a critical thickness and produces full waves that are regular with respect to time and therefore are not considered turbulent. These waves can lead to an improvement in the heat transfer coefficient by as much as 50%, compared to that of the above laminar Nusselt analysis. The mechanisms of heat transfer enhancement by wave include enlargement of liquidvapor interfacial area and decrease in mean film thickness.
The flow regime of the liquid condensate is determined by the Reynolds number defined by eq. (17). While experimental observations of the condensate indicate that the laminar flow regime usually becomes wavy in the general vicinity of Re_{δ}=30, Brauer (1956) suggested that the Reynolds number for onset of waves is related to the Archimedes number by

where Archimedes number is defined as

The wavy laminar film remains so until approximately Re_{δ}=1800, where it becomes turbulent.
In heat transfer analysis of the wavy film regime, there is no one theory reliable for calculating the heat transfer across the wavy film. Many numerical analyses have been performed and different coefficients have been found for the heat transfer coefficient and corresponding Nusselt number. Kutateladze (1982) gave the following correlation for the mean Nusselt number of film condensation on a vertical plate where wave effects are present:

To use eq. (26) to determine the heat transfer coefficient, it is first necessary to know the Reynolds number, which depends on the mass flow rate of condensate per unit width, Γ, as indicated by eq. (17). The mass flow rate of the condensate can be determined by

Substituting eq. (27) into eq. (16), one obtains

Rearranging eq. (28) yields

Combining eq. (26) and (29) yields the Reynolds number for film condensation with waves as follows:

After the Reynolds number is determined by eq. (30), the heat transfer coefficient can be determined by eq. (26).
Turbulent Film Regime
If one follows the film condensate farther down the vertical wall, the film reaches a critical thickness where the waves become irregular with respect to both time and space. This is where the turbulent film regime exists and can exist even with a stagnant vapor reservoir. It follows that this regime occurs when the Reynolds number, Re_{δ}, becomes a critical value that is much higher than that of the wavy regime.
As stated in the previous section, the condensate film becomes turbulent at approximately Re_{δ}=1800.The heat transfer rate in this turbulent regime is much larger than in laminar and wavy flow. The turbulent region in film condensation, as in any turbulent flow, is extremely difficult to model, and all accurate results come from empirical correlations or from detailed numerical modeling. For turbulent flow of condensation on a vertical plate, Labuntsov (1957) recommended the following empirical correlation:

Chun and Seban (1971) obtained the following experimental correlation for the local turbulent heat transfer coefficient for evaporation of water from a vertical wall, which is applicable for local heat transfer coefficient for condensation:

Butterworth (1983) obtained the average heat transfer coefficient for film condensation that covers laminar, wavy laminar, and turbulent flow by combining eqs. (19), (26), (31) and as follows

The Reynolds number, Re_{δ}, is needed in order to use eq. (33) to determine the heat transfer coefficient for turbulent film condensation. Equation (29) was obtained by energy balance and it is valid for all film condensation regimes. Combining eqs. (33) and (29), the Reynolds number for turbulent flow is obtained:

which can be used together with eq. (33) to determine the heat transfer coefficient for turbulent film condensation.
References
Brauer, H., 1956, “Stromung und Warmeubergang bei Reiselfilmen,” VDI Forschung, Vol. 22, pp. 140.
Butterworth, D., 1983, “Film Condensation of Pure Vapor,” Heat Exchanger Handbook, Chapter 2.6.2, Hemisphere, Washington, DC.
Faghri, A., and Zhang, Y., 2006, Transport Phenomena in Multiphase Systems, Elsevier, Burlington, MA
Faghri, A., Zhang, Y., and Howell, J. R., 2010, Advanced Heat and Mass Transfer, Global Digital Press, Columbia, MO.
Kutateladze, S.S., 1982, “Semiempirical Theory of Film Condensation of Pure Vapors,” International Journal of Heat and Mass Transfer, Vol. 25, pp. 653660.
Labuntsov, D.A., 1957, “Heat Transfer in Film Condensation of Pure Steam on Vertical Surfaces and Horizontal Tubes,” Teploenergetika, Vol. 4, pp. 7280.
Nusselt, W., 1916, “Die Oberflächenkondensation des Wasserdampfes,” Z. Vereins deutscher Ininuere, Vol. 60, pp. 541575.