YTUP
Journals
About
Services
Guides
Sign InSubmit Article
HomeJournalsJournal of Thermal Engineering10.18186/thermal.1374673
JoJournal of Thermal Engineering
Get Alerted Download PDF
AbstractKeywordsIntroductionModeling ApproachResults And Discussion16. The maximum interface temperature Tsl for the contactConclusionNomenclatureReferencesShare and CiteRelated Articles
Article Open Access1 January 2023

Heat transfer model for dropwise condensation on hydrophobic and superhydrophobic interfaces

Order Reprints Cite Share

R. YUVARAJ1, and D. SENTHILKUMAR1

1Department of Mechanical Engineering, Sona College of Technology, Salem, Tamilnadu, 636005, India

Journal of Thermal Engineering 2023, Vol. 9, Issue 5, pp. 1245-1259; doi.org/10.18186/thermal.1374673

Download PDF View DOI record

Abstract

Heat transfer models for condensation on hydrophobic and superhydrophobic interfaces are broadly available based on thermal resistance correlations. In the previous studies, very few models are presented based on the scaling factor or Nusselt number, and no model is available that directly correlates Biot number. This study develops a heat transfer model for dropwise condensation underneath a horizontal surface. The present model correlates with the Biot number to predict the heat transfer, temperature variation at the interfaces, solid-liquid, and liquid-vapor, and the growth rate of droplet condensate on the hydrophobic and superhydro-phobic interfaces by using Archimedes’ hat-box theorem. The present model is validated with analytical and experimental results against hydrophobic and superhydrophobic contact angles of similar working parameters made excellent agreements. The analytical model for dropwise condensation produces inaccurate results due to discrepancies and discontinuities due to mul-tiple correlations in the modeling. The present model is modified to obtain a continuous result using experimental data. The modified model is used for analyzing heat transfer by varying Biot numbers from 0.0001 to 1000 using Python 3.6.1 with an accuracy of 10-4. Simulation of the present model results in constant heat transfer at Bi = 4, irrespective of the contact angle. A negligible amount of coating resistance and interface resistance when Bi > 0.1, curvature effect when Bi > 0.04, droplet resistance when Bi < 0.02, the maximum liquid-vapor interface tem-perature at Bi ≈ 10, and maximum solid-liquid interface temperature at Bi ≈ 5, are presented.

Keywords: Dropwise Condensation; Biot Number; Heat Transfer; Hydrophobic Surface; Superhydrophobic Surface; Interface Temperature

Introduction

Condensation is a phase change mechanism widely used in many applications such as power generation [1-3], desalination [4], electronic thermal management [5], heat exchangers [6], building energy management [7-8], refrigeration, and air conditioning system [9]. Condensation on

non-wetting surfaces offers dropwise condensation due to low surface energy. The focus on dropwise condensation has increased widely since its discovery in 1930 [10]. This mechanism yields heat transfer rates up to an order of magnitude higher than filmwise condensation [11-12]. The low surface energy of a substrate offers increases in contact

*Corresponding author. *E-mail address: yuvarajr@sonatech.ac.in This paper was recommended for publication in revised form by Regional Editor Omid Mahian Published by Yıldız Technical University Press, İstanbul, Turkey Copyright 2021, Yıldız Technical University. This is an open access article under the CC BY-NC license (http://creativecommons.org/licenses/by-nc/4.0/).

angle (θ) from hydrophilic (θ < 90°) to hydrophobic (θ > 90°) surfaces. Several works have since fabricated superhydrophobic (θ ≥ 150°) θ ≥ 150° surfaces to achieve a variety of applications, including self-cleaning, defrosting, and condensation heat transfer enhancement. The initial droplet size of a diameter < 1 μm condensate on the dropwise condensation process at discrete nucleation sites distributed over a cooled surface [13]. The minimum thermodynamically viable radius (rmin) [14] of the condensed droplet grows due to condensation and coalescence with neighboring droplets until they become large enough (rmax) to be removed by gravity. The Rolling and merging of large droplets grow up to the maximum diameter and depart from the surface. This allows new primary droplets of radius rmin to form and grow from nucleation sites to the maximum radius rmax. In dropwise condensation, heat from the vapor transfers through the liquid-vapor interface and then transfers through the droplet by conduction to the cooled wall. Though many experimental studies [15-31] and numerical studies [32-35] of dropwise condensation on hydrophobic and superhydrophobic surfaces have been conducted in the past decades, accurate modeling of the heat transfer process has lagged. Kim and Kim proposed an analytical model to predict the heat transfer rate (Qd) through a single droplet by considering heat drop due to interfacial resistance (ΔTi ), curvature resistance (ΔTc), droplet resistance (ΔTdrop), and the coating resistance (ΔTcoat) [36]. Chavan et al. proposed an individual heat transfer model with the Nusselt number Nu defined as the ratio of convective resistance to conductive resistance and the Biot number Bi defined as the ratio of internal resistance to surface resistance. Adhikari et al. presented the average heat flux through the base of the droplet as [35] with the scaling factor f (Bi, θ)is a function of Biot number Bi and contact angle θ, for the conduction contribution to the overall thermal resistance. Most of the previous studies’ assumed f = 0 for an infinitely conductive droplet [14,39]. Several studies have been performed to determine the formulae for f (Bi, θ). Sadhal and Martin proposed a model to find the scaling factor for contact angle (θ < 90°) [41]. Yuvaraj et al. [42] reported dropwise condensation on the superhydrophobic surface with regular complete sphere-like droplets, increasing the contact surface area between the vapor and the plate surface. There is no room for the formation of films on the superhydrophobic surface. Instead, the water particles are formed as droplets, resulting decrease in resistance to the flow of heat and the fall-off diameter of the droplet from the superhydrophobic surface. This causes an increase in the heat transfer coefficient of dropwise condensation and the rate of condensation. Chavan et al. have reported numerical and experimental results for hydrophobic surfaces (θ > 90°) along with three different correlations to predict f (Bi, θ) [21]. Generally, this model would provide discontinuity in the simulation for dropwise condensation analyses. Accurate modeling of the continuous heat transfer process

from nano-size to the macro size of the condensate droplet on hydrophobic and superhydrophobic interfaces is currently unavailable. In the present study, a simple heat transfer model for dropwise condensation in terms of Biot number is developed to provide continuous results for heat transfer, temperature variation at the interfaces, and the growth rate of droplet condensate on the hydrophobic and superhydrophobic interfaces by using Archimedes’ hat-box theorem. The present model is validated with analytical results obtained using the Kim and Kim model [36] and experimental results given by Chavan et al. [21] of similar working parameters made excellent agreements. The continuous results from the nano to the macro level are obtained from the present model and compared with the experimental data presented by Chavan et al. [21]. The modified model is used for analyzing heat transfer, heat flux, the growth rate of the droplet, and heat drops due to various thermal resistances and interface temperatures by varying Biot numbers Bi = 0.001 → 1000 using Python 3.6.1 with an accuracy of 10-4. Python is an open-source high-level computer language and very simple compared to other high-level languages. The present dropwise condensation heat transfer model is developed and executed by creating code using Python version 3.6.1. Python contains many libraries, including NumPy, MatPlotLib, Xlsxwriter, etc., for solving linear algebraic equations, plotting the results, and storing the data. Python 3.6.1 version also supports the plotting of contours and heat maps.

Modeling Approach

The analytical heat transfer model presented by Kim and Kim is given as [36] (1)

where hi is the interfacial heat transfer coefficient, r is the radius of the droplet, θ is the contact angle, k is the thermal conductivity of the water droplet, δ is promoter coating layer thickness, kcoat is the thermal conductivity of the coating, and rmin is the minimum thermodynamic viable radius of the droplet. The heat transfer through a droplet also be calculated from the enthalpy dissipated by the condensing vapor [34] (6) where ρ is the density of water, hfg is the latent heat and is the growth rate of droplet condensate. The interfacial heat transfer coefficient is determined by using the kinetic theory [37-38]. (7) where α is the condensation coefficient, Rg is the gas constant, vg is the specific volume of vapor and Tsat is the saturation temperature of the vapor. Condensation coefficient α is defined as the ratio of vapor molecules captured by the liquid phase to the total number of vapor molecules reaching the liquid surface (α = 0 → 1). In the present work, a clean environment is assumed, and α is taken as 1. Chavan et al. proposed an individual heat transfer model as [21]

(8) where Nu is the Nusselt number defined as the ratio of convective resistance to conductive resistance, and Bi is the

Biot number defined as the ratio of internal resistance to surface resistance. Adhikari et al. presented the average heat flux through the base of the droplet as [35] (9) where, f (Bi, θ) is the scaling factor, a function of Biot number Bi and contact angle θ, for the conduction contribution to the overall thermal resistance. Most of the previous studies’ assumed f = 0 for an infinitely conductive droplet [14,39]. For f = 0.25 as assumed by Graham and Griffith [40], the conduction-front propagation velocity is greater than the growth velocity . Several studies have been performed to determine the formulae for f (Bi, θ). Figure 1 shows the water droplet underneath a flat surface of radius r, contact radius rc, droplet height h, contact angle θ with negligible hysteresis, equivalent height H, the volume of the droplet V, and its equivalent volume Ve is obtained from Archimedes’ hat-box theorem [43]. By equating V = Ve, the relation between h and H can be written as (10) Rose and Glicksman reported that most heat is transferred through tiny droplets of constant condensation rate with time [46]. Droplet diameter d < 150 μm transfer 90% of total heat for low-pressure condensation, which occupies 35% of surface area, and d < 40 μm transfer 90% of heat at atmospheric pressure with 23% of the surface area [40]. Droplets in this size range are sufficiently small such that Marangoni circulation and other convection effects can be neglected [45-46]. Hence, conduction is the primary

mode of heat transfer in dropwise condensation. As droplets grow, the internal temperature distribution varies, and conduction-front propagation velocity is much greater than the droplet growth rate. Therefore, the heat transfer process can be modeled as quasi-steady, and the in-droplet condensation heat transfer process reduces to conduction in a solid medium [35]. With this assumption, the conductive resistance of the droplet can be rewritten as (11) The Biot number Bi and the Nusselt number Nu are defined in terms of contact radius rc as (12)

(13) After substituting and rearranging the eq (10) and eq (11) in eq (5), the heat transfer model in terms of Biot number can be expressed as

(14) where c = 6πk2/hi is the condensation coefficient, a is the empirical value and b = Bicoat/Bi is the Biot number ratio between a coating to the droplet. The following parameters are assumed as similar to the previous works [34-36], interface heat transfer coefficient hi = 400 MW/ m2 K , the thermal conductivity of the water droplet k =

0.6. W/mK, the thermal conductivity of the coating layer

kc = 0.2 W/mK, coating thickness δ = 1 μm, the minimum radius rmin ≈ 10 nm, Subcooling temperature ΔT = 10 K, the density of the liquid ρ = 998 kg/m3, latent heat of condensation hfg = 2.3 MJ/kg, and the saturation temperature Tsat = 373.15 K. The Biot number for coating is obtained as Bicoat = 2 from the parameters mentioned above. The value of a is obtained from the experimental values given by Chavan et al. After substituting the values of a and b in eq (14) further, the equation reduced to

through the droplet, predict interface temperatures, heat drop due to various resistances, and growth rate of condensate on hydrophobic and superhydrophobic interfaces for varying Biot numbers (Bi = 0.0001 → 1000). (16) Eq (16) is used to find the growth rate of condensate on hydrophobic and superhydrophobic interfaces for varying Biot numbers (Bi = 0.0001 → 1000 [34].

Results And Discussion

Validation of the Present Model with Analytical Model [36] and Numerical Model [21] The present model is validated by an analytical model proposed by Kim and Kim [36] and a numerical model presented by Chavan et al. [21] for the growth rate of the droplet as a function of the radius (r) of the droplet given in eq (16). The analytical model produces continuous results for different contact angles (θ = 90°, 120°, 150° & 170°) as shown in Figure 2. Whereas, in the numerical model, the discontinuities are noted in Figure 2. due to three different correlations, i) Bi < 0.5 ii) 0.5 < Bi < 2, and iii) Bi > 2, for predicting the Nu are used by Chavan et al. [21] Chavan et al. [21] compared the experimental growth of the droplet with numerical results and Kim’s analytical results and reported analytical model underestimates the growth rate due to the lack of taking into account local heat transfer effects at the three-phase contact line. Also, it was observed in the experiments that the large droplet diameters

(15) Eq (15) is the modified heat transfer model used in the present work without noncondensable gases during condensation. The present model can simulate heat transfer

Figure 2. Validation of the present model with Kim’s analytical model [36] and Chavan’s numerical model [21] regarding growth rate.

growing on the superhydrophobic surface lower the growth rates due to added conduction thermal resistance and low droplet diameters growing faster due to interfacial thermal resistance dominated growth [21]. When the contact angle θ = 90°, the present model deviates from the numerical model at a radius of the droplet r ≈ 50 μm and it produces an identical result at r = 1 mm. The present model made good agreements with the numerical results when the contact angle increases from θ = 90° to θ = 170° as shown in Figure 2. Hence, the present model is best suitable for dropwise condensation on hydrophobic and superhydrophobic interfaces. Validation of the Present Model with the Experimental works [21,24] The present model is also validated by experimental results obtained by Chavan et al. [21] and Enright et al. [24] for the radius of the droplet r as a function of time t. The experimental results are shown in Figures 3 (a) and

(b) obtained by Chavan et al. for the subcooling temperature ΔTexperimental = 5°C and corresponding contact angles θ = 140° and θ = 170° respectively. In the present work, the ΔTnumerical = 0.0545°C for θ = 140° and ΔTnumerical = 0.092°C for θ = 170° are taken by mapping the experimental results and made an excellent agreement for the case θ = 140° and good agreement made for the θ = 170°. The difference between the present model and the experiments can be recognized as the presence of noncondensable gases (NCG). Rose reported that the condensation of water vapor in the presence of NGCs (air) acts as a diffusion barrier for water vapor [15]. The present model is also validated by another experimental result obtained by Enright et al. [24] best fits with ΔTnumerical = 0.064°C as shown in Figure 3 (c). To obtain the complete simulation of the present model by varying contact angles from θ = 90° → 180°, the ΔTnumerical = 0.0545°C is used due to its best fit with experimental results.

Figure 3. Validation of present work with Kim’s analytical [36] and Chavan’s experimental work [21] for (a) contact angle θ = 140°, (b) contact angle θ = 170° and Enright’s experimental work [24] for (c) contact angle θ = 120°.

Simulation of Heat Transfer Through a Single Droplet In the present work, the heat transfer, heat drop, interface temperature variation, and growth rates are analyzed by simulating the eq (14) for different values of Biot number and contact angles. Since the Biot number is a function of the contact angle and radius of the droplet, Bi = f (θ,r) the variation of droplet radius is shown in Figure 4. The droplet radius measures the macro scale for the Biot numbers in the range of 10 > Bi > 1000, microscale for the Biot numbers in the range of 0.01 > Bi > 1, and nanoscale for the Biot numbers 0.00001 > Bi > 0.001. Since the minimum thermodynamically viable radius rmin ≈ 10 nm size, the range of Biot numbers in the present work is chosen from Bi = 0.0001 → 1000. Because of the spherical nature of the droplet, the radius of the droplet varies in a similar pattern against the contact angles θ = 90° to θ ≈ 180° for different Biot numbers, as shown in Figure 4.

Figure 4. Variation of droplet radius with contact angle for different Biot numbers (Bi = 0.001 → 1000).

The heat transfer through a single droplet Qd is shown in Figure 5 for different Biot numbers Bi = 0.01 → 1 in Figure (a) and Bi = 4 → 1000 in Figure 5 (b). The thermodynamically viable radius of a droplet can be expressed in the form as (17) Where σ is the surface tension of the water, in the literature, the subcooling temperature has been taken in the range ΔT = 1K - 10K. The subcooling temperature ΔT = 10K and the corresponding minimum viable radius of the droplet rmin ≈ 10 nm are used to solve the present model [21,34,36]. In the present work, the numerical subcooling temperature is taken as ΔTnumerical = 0.0545 K for the best fit of experimental data and the respective viable radius rmin ≈ 400 nm is used for simulation. The curvature effect is used in eq (14) becomes negative when the radius of the droplet r is smaller than the rmin value i.e. 400 nm. From this point forward, whatever the values are in the negative sign means that r < rmin. The heat transfer through a single droplet Qd is plotted as a function of contact angles varying from θ = 90° → 180° as shown in Figure 5. When the Bi < 0.3, the heat transfer through the droplet is very small due to high interface resistance and curvature resistance. The heat transfer rate increases with an increase in contact angle for Bi > 0.3 as shown in Figure 5 (a). Further, with an increase in Biot number towards Bi = 1 the variation of heat transfer follows linearly and approaches a constant line when Bi = 4 as shown in Figure 5 (b). The heat transfer remains constant for all the values of contact angles in both hydrophobic and superhydrophobic sides at Bi = 4. The heat transfer rate decreases with an increase in contact angle when the Biot number is increased from Bi = 4 → 1000 as shown in Figure 5 (b). It is essential to study the variations in the region of Bi < 0.3, to 1000.

(a) (b) Figure 5. Heat transfer rate for different contact angles (a) Bi = 0.01 → 1 and (b) Bi = 4 → 1000.

Figure 6. Heat transfer contours for different contact angles (a) Bi = 0.001 → 1, (b) Bi = 1 → 1000 , (c) Bi = 0.01 → 1 , (d) Bi = 0.1 → 0.3, (e) Bi = 1 → 4, and (f) Bi = 4 → 10. The heat transfer variations from Bi = 0.001 → 1000 are predicted by using the modified model eq (14), and contours are plotted as shown in Figure 6 (a-b). The maximum heat transfer rate in the superhydrophobic region for lower Biot numbers shown in Figure 6 (a) is shifted to the hydrophobic region for higher Biot numbers as shown in Figure 6 (b). Figure 6 (c-f) shows the higher

heat transfer shifting behavior for different ranges of Biot numbers Bi = 0.01 → 0.1, Bi = 0.1 → 0.3, Bi = 1 → 4, and Bi = 4 → 10 respectively. The cold spot (lower heat transfer rate) in the hydrophobic region and the hot spot (higher heat transfer rate) in the superhydrophobic region for Bi = 0.01 → 0.1 are shown in Figure 6 (c). The cold spot propagates hydrophobic to the superhydrophobic region,

(a) (b) Figure 7. Heat flux with different contact angles for (a) Bi = 0.01 → 1 and (b) Bi = 1 → 1000. and Bi = 0.1 → 0.3 both regions become cold spots, as shown in Figure 6 (d). An important finding is observed in Figure 6 (e) that the hot spot is almost the same for all contact angles at Bi = 4. In Figure 6 (f), the hot spot is shifted from the superhydrophobic to the hydrophobic region at Bi = 10 and the heat transfer rate decreases as the contact angle increases. The heat flux qd is the ratio between the heat transfer rate of a single droplet Qd and the contact area of the droplet between the solid and liquid interface. Figure 7 (a) shows the heat flux varies nearly similar to the heat transfer rate up to Bi = 4. In Figure 7 (b), the heat flux is constant throughout the hydrophobic and superhydrophobic region at Bi = 4, and afterward, increase in Biot number decreases the heat flux due to the high conductive resistance offered by the solid-liquid interface area. When the Bi > 4, the heat flux decreases for increases in contact angle.

The heat drop due to various resistances, interface resistance ΔTInt, curvature resistance ΔTCurv, droplet resistance ΔTDrop, and coating resistance ΔTCoat, are analyzed and shown in Figures 8-11 for different Biot numbers. In the condensation process, the first and foremost resistance is interfacial resistance. Eq (1) is used to predict the heat drop due to interfacial resistance for Bi = 0.01 → 1 as shown in Figure 8 (a), the maximum heat drop is attained at Bi = 1, and for Bi = 1 → 1000 as shown in Figure 8 (b), the heat drop decreases when the Biot number increases. For all the Biot numbers, when Bi > 0.3, minimum heat drop in the superhydrophobic region is occurred due to the large interface area offered by the superhydrophobic surface. When the contact angle exceeds 170°, Figures 8 (a) and (b) show a sudden heat drop. The increase in the contact angle decreases the heat drop due to the greater interfacial resistance offered by the liquid-vapor

(a) (b) Figure 8. Heat drop due to interfacial resistance between droplet and vapor with contact angles (θ = 90° → 179°) for different Biot numbers (a) Bi = 0.01 → 1 and (b) Bi = 1 → 1000.

Figure 9. Heat drop due to curvature effect with contact angles (θ = 90° → 179°) for different Biot numbers (Bi = 0.0001 → 1000).

Figure 10. Heat drop due to conductive resistance inside the droplet with contact angles (θ = 90° → 179°) for different Biot numbers (Bi = 0.01 → 1000).

interface. For a lower Biot number, the liquid-vapor surface area approaches zero, resulting in heat drop due to interfacial resistance, as shown in Figure 8 (a). The maximum heat drop of ΔTInt = 0.008 K is obtained for Bi = 1 and θ = 90°, which is 15% of the total heat drop ΔTInt = 0.0545 K. The heat drop due to curvature resistance ΔTCurv is predicted using the eq (2) for different Biot numbers shown in Figure 9. The curvature resistance only depends on the droplet’s radius, and the variation shows an identical pattern for all the Biot numbers. The maximum heat drop in Figure 9 is noted for the lowest Biot number Bi = 0.0001 and decreases with an increase in Biot numbers. For larger values of Biot numbers, the curvature effect can be neglected. The curvature resistance decreases for increases in contact angles, and the maximum heat drop for the Bi = 1 is obtained as ΔTCurv = 0.02, which is 37% of the total heat drop.

Eq (11) is the modified heat drop ΔTDrop due to droplet resistance and predicts the heat drops for different Biot numbers, as shown in Figure 10. The conductive resistance offered by the droplet is the primary source of heat drop for Bi > 0.1. The heat drop due to conductive resistance increases with an increase in contact angles for the Biot numbers Bi = 0.1 → 10 and closely follows linear variation. When the Biot number is greater than 100, the heat drop decreases in the hydrophobic region from contact angle θ = 90° → 140° and increases in the superhydrophobic region. The maximum heat drop occurs at a contact angle θ ≈ 180° for Bi = 1 is ΔTDrop = 0.02 K and at θ = 90° ΔTDrop = 0.005 K which is 9% of the total heat drop. The coating resistance is predicted by using Eq (4) for different Biot numbers shown in Figures 11 (a) and (b). The

(a) (b) Figure 11. Heat drop due to coating resistance with contact angles (θ = 90° → 179°) for different Biot numbers (a) Bi =

0.01. → 1 and (b) Bi = 1 → 1000.

Figure 12. Total heat drop through the droplet with contact angles (θ = 90° → 179°) for different Biot numbers (Bi = 0.001 → 1000).

Figure 13. Variation of Nusselt number as a function of Biot number (Bi = 0.01 → 1000) for different contact angles (θ = 90° → 170°).

coating thickness δ = 1 μm is assumed to predict heat drop due to the coating resistance ΔTCoat. Figure 11 (a) shows the heat drop variation for Bi = 0.01 → 1 with an increase in contact angle increases in heat drop obtained. At Bi = 3, the interpretation of heat drop closely follows constantly for all the contact angles in the hydrophobic and superhydrophobic regions, as shown in Figure 11 (b). The contact area of liquid and solid increases when the radius of the droplets increases from micro size to macro size. Since the droplet radius is under micron size for Biot number Bi = 1, the heat drop due to coating resistance decreases when the Biot number increases greater than one. Further increase in Biot number offers a decrease in heat drop due to coating resistance. When Bi > 3, the heat drop decrease with an increase in contact angles. The total heat drop is shown in Figure 12 for the Biot numbers varies from Bi = 0.001 → 1000. The total heat drop ΔTTotal decreases with an increase in contact angle when the Bi < 1. There is no variation in the total heat drop for Bi = 1, and it remains constant for all the contact angles in both the hydrophobic and superhydrophobic regions. The inner resistance to the surface resistances of the droplet is the same at Biot number Bi = 1. The heat drop increases when the contact angle increases for the Bi > 1. For more significant value Biot numbers, superhydrophobic regions offer higher total heat drops, as shown in Figure 12. The Nusselt number, the ratio of conductive resistance to convective resistance, and variation for different Biot numbers are shown in Figure 13. The superhydrophobic interface gives a more prominent value Nusselt number compared with the hydrophobic interface up to Bi ≈ 5 with a minimal magnitude due to low conductive resistances than convection. Whereas Biot number Bi > 5, the shift has occurred, and the hydrophobic interface offers a more

Figure 14. Growth rate variation with Biot number (Bi = 0.0001 → 1000) for different contact angles θ = 90°, 120°, 150° & 170°.

considerable value of Nusslet number than the superhydrophobic interface due to higher conductive resistance than the convection. The growth rate dr/dt as a function of Biot number Bi = 0.0001 → 1000 for different contact angles is shown in Figure 14. For all contact angles, the growth rate remains constant for a specific range of Biot numbers. The content of growth rate increase for superhydrophobic interfaces with lower magnitude and higher growth rate with a comparatively short degree for hydrophobic interfaces are shown in Figure 14 with an accuracy of 10-4.

Comparison of Hydrophobic and Superhydrophobic Heat Transfer In the previous section, the simulation of heat transfer is predicted by using the present modified heat transfer model given by eq (14) which correlates Biot number and contact angle. For complete simulation for the range of contact angle θ = 90° → 180° and Biot numbers Bi = 0.00O1 → 1000, the numerical subcooling temperature

was taken as ΔTnumerical = 0.0545. The present modified model validated with the experimental work, as shown in Figure 3 (a) and (b) are directly used to compare the heat transfer for the contact angles θ = 140° and θ = 170°. The numerical subcooling temperature for the contact angle θ = 170° is taken as ΔTnumerical = 0.092 for the best fit, as shown in Figure 3 (b). The crucial findings are noted for different resistances, and Nusselt number variations are

Figure 15. Comparison of hydrophobic (θ = 140°) and superhydrophobic θ = 170°) heat transfer for different Biot numbers (Bi = 0.0001 → 1000) (a) interfacial heat drop (b) curvature heat drop, (c) droplet heat drop, (d) coating heat drop, (e) total heat drop and (f) Nusselt number.

Figure 16. Interface temperature as a function of the Biot number for different contact angles (θ = 140° & 170°).

shown in Figure 15 (a-f). The interfacial resistance affects the heat drop intensively when the Biot number Bi < 0.1 due to low conductive resistance and remains constant for both contact angles, as shown in Figure 15 (a). The curvature effect is significant for the Biot number Bi < 0.04 and remains constant, as shown in Figure 15 (b). The vital finding observed in Figure 15 (c) is that a massive heat drop occurred in the range of Biot numbers Bi = 1 → 100 due to higher heat transfer through the droplet by conduction. The heat drop due to conduction for the superhydrophobic interface is higher than for the hydrophobic interface. When the Bi < 0.01, the heat transfer through the droplet is negligible, and most of the heat drop is carried by the curvature effect and interface resistance. The coating resistance is also significant for the Boit number

Bi < 0.05, as shown in Figure 15 (d). For Biot number in the range Bi = 0.01 → 1, the total heat drop produces a marginal difference, significant variation for Biot number Bi < 0.01, and a moderate difference is noted for higher Biot numbers as shown in Figure 15 (e). Figure 15 (f) shows the variation of the Nusselt number as a function of the Biot number in which the variation is shifted at Bi ≈ 5, and the Nusselt number decreases for higher contact angles. The interface temperatures, liquid-vapor interface, and solid-liquid are also predicted by using the Eqns (1, 2, 4, 11) for the contact angles θ = 140° and θ = 170°. The liquid-vapor interface temperatures Tly decreases immensely for the Biot number Bi < 2 due to a large amount of heat transfer through tiny droplets and remains unchanged for the higher value of Biot numbers due to higher conductive resistance offered by the droplets. The maximum interface temperature Tly is obtained at Biot number Bi > 2 for both contact angles θ = 140°and θ = 170°. The solid-liquid interface temperature Tsl exactly follows its variation up to Bi ≈

16. The maximum interface temperature Tsl for the contact

angle, θ = 140° and θ = 170°are obtained at Bi ≈ 0.5 and Bi ≈

0.4. respectively. The interface temperature Tsl decreases for

the higher values of Biot number. Heat Transfer Rate Variation with Model Sensitivity The present model sensitivity is observed by changing the subcooling temperature ΔT and minimum viable radius rmin. Figure 17 (a) and (b) show that the heat transfer rate variation for different Biot numbers varies from Bi = 0.001 → 1000 by changing ΔT and rmin by 5% and 10% respectively. The present model is low sensitive for lower size droplets of radius r < 10 μm and slightly higher sensitive for larger droplets of radius r > 10 μm as shown in Figure 17 (a). At Bi

(a) (b) Figure 17. Heat transfer rate variation with model sensitivity by changing ΔT & rmin by (a) 5% decrement and (b) 10% decrement.

= 10, the hydrophobic interface produces a 4% reduced heat transfer rate, and the superhydrophobic interface makes a 5% reduced heat transfer rate for a 5% decrement in the values of ΔT and rmin. At Bi = 1000, the maximum reduction of 6% heat transfer obtained in the superhydrophobic interface and hydrophobic interface provides a 5% reduction in heat transfer. Whereas for a 10% decrement, the maximum reduction in heat transfer of 11% is obtained from the superhydrophobic interface at Biot number Bi = 1000 as shown in Figure 17 (b).

Conclusion

An unique heat transfer model for dropwise condensation on the hydrophobic and superhydrophobic interface is developed by correlating the Biot number and contact angle. The developed model is validated with analytical, numerical, and experimental work, offering excellent agreement. The present heat transfer model is modified for the best fit of observed values. The modified heat transfer model is used to predict the heat transfer through the single droplet, heat flux, heat drops due to various resistances such as interfacial resistance, curvature resistance, droplet resistance, coating resistance, the growth rate of the droplet, and interfacial temperatures for different Biot numbers varies from Bi = 0.0001 → 1000 and different contact angles in the range θ = 90° → 180°. Besides, heat transfer contours are plotted for other Biot numbers as a function of contact angles. From the simulation of the present modified model, the following are the findings for dropwise condensation on hydrophobic and superhydrophobic interfaces, • A continuous single correlation model predicts the heat transfer for a vast range of Biot numbers Bi = 0.0001 → 1000 • For the best fit, the numerical value for subcooling temperature is found as ΔTnumerical = 0.0545 K against the experimental value ΔTexperimental = 5 K. • At Bi = 4, the heat transfer rate remains unchanged for all the contact angles in the hydrophobic and superhydrophobic regions. • A wide variation of the growth rate for all the Biot numbers in the range Bi = 0.0001 → 1000 is presented. • The heat drop due to interfacial resistance reaches a maximum value at Bi = 1. • For Bi = 1, the total heat drops are almost constant for all the contact angles in the hydrophobic and superhydrophobic regions. • The Nusselt number shift occurs at Bi ≈ 5. • The interfacial resistance and curvature effect, and coating resistance are significant for Bi < 0.01, and droplet resistance is substantial for Bi > 0.05. • The liquid-vapor interface temperature reaches its maximum value at Bi ≈ 5 and remains constant afterward. • The solid-liquid interface temperature reaches a maximum at Bi ≈ 0.4 → 0.5 and decreases afterward.

Nomenclature

Bi Biot number Nu Nusselt number θ Contact angle, degree r Droplet radius, m rmin Minimum viable radius, m rmax Maximum radius, m Q Heat transfer rate, W ΔT Temperature difference, K hi Interface heat transfer coefficient, W/m2K δ Coating thickness, m k Thermal conductivity, W/mK ρ Density, kg/m2 hfg Latent heat, J/kg dr/dt Growth rate, m/s rc Contact radius, m h Droplet height, m V Volume, m3 H Effective height, m a & b Empirical value c Condensation coefficient

References

  1. REFERENCES
  2. [1] McNeil DA, Burnside BM, Cuthbertson G. Dropwise condensation of steam on a small tube bundle at turbine condenser conditions. Exp Heat Transf 2000;13:89105. [CrossRef]
  3. [2] Kurshakov AV, Ryzhenkov AV, Bodrov AA, Ryzhenkov OV, Patakin AA, Chernov EF. Heat transfer enhancement in steam-turbine condensers with the use of surface-active substances. Therm Eng 2014;61:785789. [CrossRef]
  4. [3] Beér JM. High efficiency electric power generation: The environmental role. Prog Energy Combust Sci 2007;33:107134. [CrossRef]
  5. [4] Lara JR, Noyes G, Holtzapple MT. An investigation of high operating temperatures in mechanical vapor-compression desalination. Desalination 2008;227:217232. [CrossRef]
  6. [5] Peters TB, McCarthy M, Allison J, Dominguez-Espinosa FA, Jenicek D, Kariya HA, et al. Design of an integrated loop heat pipe air-cooled heat exchanger for high performance electronics. IEEE Trans Components Packag Manuf Technol 2012;2:16371648. [CrossRef]
  7. [6] Kim MH, Bullard CW. Air-side performance of brazed aluminum heat exchangers under dehumidifying conditions. Int J Refrig 2002;25:924934. [CrossRef]
  8. [7] Li B, Yao R. Urbanisation and its impact on building energy consumption and efficiency in China. Renew Energy 2009;34:19941998. [CrossRef]
  9. [8] Pérez-Lombard L, Ortiz J, Pout C. A review on buildings energy consumption information. Energy Build 2008;40:394398. [CrossRef]
  10. [9] Sett S, Sokalski P, Boyina K, Li L, Rabbi KF, Auby H, et al. Stable Dropwise Condensation of Ethanol and Hexane on Rationally Designed Ultrascalable Nanostructured Lubricant-Infused Surfaces. Nano Lett 2019;19:52875296. [CrossRef]
  11. [10] Schmidt E, Schurig W, Sellschopp W. Versuche über die Kondensation von Wasserdampf in Film- und Tropfenform. Techn Mech Thermodyn 1930;1:5363. [Deutsch] [CrossRef]
  12. [11] Le Fevre EJ, Rose JW. An experimental study of heat transfer by dropwise condensation. Int J Heat Mass Transf 1965;8:11171133. [CrossRef]
  13. [12] Rose JW. On the mechanism of dropwise condensation. Int J Heat Mass Transf 1967;10:755762. [CrossRef]
  14. [13] Stylianou SA, Rose JW. Drop-to-filmwise condensation transition: Heat transfer measurements for ethanediol. Int J Heat Mass Transf 1983;26:747760. [CrossRef]
  15. [14] Burnside BM, Hadi HA. Digital computer simulation of dropwise condensation from equilibrium droplet to detectable size. Int J Heat Mass Transf 1999;42:31373146. [CrossRef]
  16. [15] Rose JW. Dropwise condensation theory and experiment: A review. Proc Inst Mech Eng Part A J Power Energy. 2002;216:115128. [CrossRef]
  17. [16] Vemuri S, Kim KJ. An experimental and theoretical study on the concept of dropwise condensation. Int J Heat Mass Transf 2006;49:649657. [CrossRef]
  18. [17] Leach RN, Stevens F, Langford SC, Dickinson JT. Dropwise condensation: Experiments and simulations of nucleation and growth of water drops in a cooling system. Langmuir 2006;22:88648872. [CrossRef]
  19. [18] Zhong L, Xuehu M, Sifang W, Mingzhe W, Xiaonan L. Effects of surface free energy and nanostructures on dropwise condensation. Chem Eng J 2010;156:546552. [CrossRef]
  20. [19] Miljkovic N, Enright R, Wang EN. Effect of droplet morphology on growth dynamics and heat transfer during condensation on superhydrophobic nanostructured surfaces. ACS Nano 2012;6:17761785. [CrossRef]
  21. [20] Peng B, Ma X, Lan Z, Xu W, Wen R. Experimental investigation on steam condensation heat transfer enhancement with vertically patterned hydrophobic–hydrophilic hybrid surfaces. Int J Heat Mass Transf 2015;83:2738. [CrossRef]
  22. [21] Chavan S, Cha H, Orejon D, Nawaz K, Singla N, Yeung YF, et al. Heat transfer through a condensate droplet on hydrophobic and nanostructured superhydrophobic surfaces. Langmuir 2016;32:77747787. [CrossRef]
  23. [22] Macner AM, Daniel S, Steen PH. Simulating heat transfer during transient dropwise condensation on a low-thermal-conductivity substrate. Langmuir 2019;35:1156611578. [CrossRef]
  24. [23] Enright R, Miljkovic N, Al-Obeidi A, Thompson CV, Wang EN. Condensation on superhydrophobic surfaces: the role of local energy barriers and structure length scale. Langmuir 2012;28:1442414432. [CrossRef]
  25. [24] Enright R, Miljkovic N, Dou N, Nam Y, Wang EN. Condensation on superhydrophobic copper oxide nanostructures. Proceedings of the ASME 2012 Third International Conference on Micro/Nanoscale Heat and Mass Transfer 2013:135. [CrossRef]
  26. [25] Ölçeroğlu E, Hsieh CY, Rahman MM, Lau KKS, McCarthy M. Full-field dynamic characterization of superhydrophobic condensation on biotemplated nanostructured surfaces. Langmuir 2014;30:75567566. [CrossRef]
  27. [26] Li G, Alhosani MH, Yuan S, Liu H, Ghaferi AA, Zhang T. Microscopic droplet formation and energy transport analysis of condensation on scalable superhydrophobic nanostructured copper oxide surfaces. Langmuir 2014;30:1449814511. [CrossRef]
  28. [27] Baba S, Sawada K, Tanaka K, Okamoto A. Dropwise condensation on a hierarchical nanopillar structured surface. Langmuir 2020;36:1003310042. [CrossRef]
  29. [28] Cha H, Vahabi H, Wu A, Chavan S, Kim MK, Sett S, et al. Dropwise condensation on solid hydrophilic surfaces. Sci Adv 2020;6:eaax0746. [CrossRef]
  30. [29] Song Z, Lu M, Chen X. Investigation of dropwise condensation heat transfer on laser-ablated superhydrophobic/hydrophilic hybrid copper surfaces. ACS Omega 2020;5:2358823595. [CrossRef]
  31. [30] Wang R, Wu F, Xing D, Yu F, Gao X. Density maximization of one-step electrodeposited copper nanocones and dropwise condensation heat-transfer performance evaluation. ACS Appl Mater Interfaces 2020;12:2451224520. [CrossRef]
  32. [31] Ludwicki JM, Robinson FL, Steen PH. Switchable wettability for condensation heat transfer. ACS Appl Mater Interfaces 2020;12:2211522119. [CrossRef]
  33. [32] Orazzo A, Tanguy S. Direct numerical simulations of droplet condensation. Int J Heat Mass Trans 2019;129:432448. [CrossRef]
  34. [33] Stevens KA, Crockett J, Maynes D, Iverson BD. Simulation of drop-size distribution during dropwise and jumping drop condensation on a vertical surface: ımplications for heat transfer modeling. Langmuir 2019;35:1285812875. [CrossRef]
  35. [34] Birbarah P, Chavan S, Miljkovic N. Numerical simulation of jumping droplet condensation. Langmuir 2019;35:1030910321. [CrossRef]
  36. [35] Adhikari S, Nabil M, Rattner AS. Condensation heat transfer in a sessile droplet at varying Biot number and contact angle. Int J Heat Mass Transf 2017;115:926931. [CrossRef]
  37. [36] Kim S, Kim KJ. Dropwise condensation modeling suitable for superhydrophobic surfaces. J Heat Transf 2011;133:081502. [CrossRef]
  38. [37] Miljkovic N, Enright R, Wang EN. Modeling and Optimization of Superhydrophobic Condensation. J Heat Transf 2013;135:111004. [CrossRef]
  39. [38] Tanasawa I. Advances in Condensation Heat Transfer. In: Hartnett JP, Irvine TF, Cho YI (editors). Advances in Heat Transfer. Amsterdam: Elsevier.;1991. Vol. 21, p. 55139. [CrossRef]
  40. [39] Glicksman LR, Hunt AW. Numerical simulation of dropwise condensation. Int J Heat Mass Transf 1972;15:22512269. [CrossRef]
  41. [40] Graham C, Griffith P. Drop size distributions and heat transfer in dropwise condensation. Int J Heat Mass Transf 1973;16:337346. [CrossRef]
  42. [41] Sadhal SS, Martin WW. Heat transfer through drop condensate using differential inequalities. Int J Heat Mass Transf 1977;20:14011407. [CrossRef]
  43. [42] Yuvaraj R, Senthil Kumar D. Study of droplet dynamics and condensation heat transfer on superhydrophobic copper surface. Therm Sci 2021;25:653664. [CrossRef]
  44. [43] Wolfram MathWorld. Archimedes' Hat-Box Theorem. Available at: https://mathworld.wolfram.com/ArchimedesHat-BoxTheorem.html Last Accessed Date: 26.09.2023.
  45. [44] Rose JW, Glicksman LR. Dropwise condensation—The distribution of drop sizes. Int J Heat Mass Transf 1973;16:411425. [CrossRef]
  46. [45] Tanaka H, Tsuruta T. A microscopic study of dropwise condensation. Int J Heat Mass Transf 1984;27:327335. [CrossRef]
  47. [46] Xu Z, Zhang L, Wilke K, Wang EN. Multiscale dynamic growth and energy transport of droplets during condensation. Langmuir 2018;34:90859095. [CrossRef]

Share and Cite

YUVARAJ, R.; SENTHILKUMAR, D. Heat transfer model for dropwise condensation on hydrophobic and superhydrophobic interfaces. Journal of Thermal Engineering 2023, Vol. 9, pp. 1245-1259. https://doi.org/10.18186/thermal.1374673

Export:

Related Articles

Conceptualization thermal analysis and manufacturing of nano-textured micro-structured surfaces forMete Budakli, 1 January 2021Numerical analysis on heat transfer flow structure and exergy loss of combined truncated and circulaS.Mohamed ILLYAS, Kumaresan VELLAISAMY et al., 1 January 2023Numerical simulation to study mixing vane spacer effects on heat transfer performance of supercriticSatish Kumar DHURANDHAR, Shobha Lata SINHA et al., 1 January 2023Heat transfer enhancement and applications of thermal energy storage techniques on solar air collectKafel AZEEZ, Riyadh Ibraheem AHMED et al., 1 January 2023
Publication History
Published1 January 2023
Versionv1
AccessOpen Access
10.18186/thermal.1374673
Article Figures (9)
Figure 1Figure 2Figure 3Figure 4Figure 5Figure 6Figure 7Figure 8Figure 9
Related Articles
Conceptualization thermal analysis and manufacturing of nano-textured micro-structured surfaces forMete BudakliJournal of Thermal Engineering, 1 January 2021Numerical analysis on heat transfer flow structure and exergy loss of combined truncated and circulaS.Mohamed ILLYAS, Kumaresan VELLAISAMY et al.Journal of Thermal Engineering, 1 January 2023Numerical simulation to study mixing vane spacer effects on heat transfer performance of supercriticSatish Kumar DHURANDHAR, Shobha Lata SINHA et al.Journal of Thermal Engineering, 1 January 2023
Journal of Thermal Engineering coverJournal of Thermal Engineering Download PDF

Subscribe to YTUP

Stay connected and receive the latest research updates directly in your inbox.

YTUP — Yıldız Technical University Publishing

Advancing knowledge and fostering innovation through high-quality, peer-reviewed academic publications.

About YTU

Discover

  • ›Articles
  • ›Journals
  • ›Research Topics
  • ›Open Access Policy

Guidelines

  • ›Author guidelines
  • ›Services for authors
  • ›Policies and publication ethics
  • ›Editor guidelines
  • ›Fee policy

Explore

  • ›Articles
  • ›Research Topics
  • ›Journals
  • ›How we publish

Support

  • ›Help center
  • ›Emails and alerts
  • ›Contact us
  • ›Submit
  • ›Career opportunities
YTU Logo

© 2026 Yıldız Technical University (Istanbul, Turkey)

Terms and ConditionsTerms of UsePrivacy PolicyPrivacy SettingsDisclaimer
Like this platform? Join our teamHave feedback or questions?
Supervisor