Heterogeneous nucleation capability of conical microstructures for water droplets

Wei Xu, Zhong Lan, Benli Peng, Rongfu Wen and Xuehu Ma*
Liaoning Provincial Key Laboratory of Clean Utilization of Chemical Resources, Institute of Chemical Engineering, Dalian University of Technology, Dalian, China. E-mail: xuehuma@dlut.edu.cn; Fax: +86 411 84707892; Tel: +86 411 84707892

Received 14th October 2014 , Accepted 14th November 2014

First published on 17th November 2014


Abstract

The presence of microstructures on a substrate has a great effect on the heterogeneous nucleation of water droplets. A circular conical apex and a cavity are adopted as the physical model to represent the typical defects which exist widely on substrates, and classic nucleation theory is used to quantitatively analyze the nucleation capability of different microstructures at different condensation conditions. The results indicate that conical cavities with narrower cone angles can reduce the nucleation free energy barrier as compared with apexes and a planar substrate, yielding a relatively higher nucleation capability. With the vapor pressure and supersaturation increasing, the nucleation rate increases rapidly, and some of the cavities that are originally not preferred for nucleation gradually translate into active nucleation sites. Consequently, the activated nucleation sites are finite for practical substrates under certain nucleation conditions, and the nucleation sites number density can be affected by the condensation conditions and the distribution of micro cavities on the substrate. The analysis also indicated that it is possible to realize spatial control of nucleation sites by the construction of micro cavities, and the nucleation sites number density can be intensified by increasing the amount of micro cavities on the substrate.


1. Introduction

It can be observed that, when dew drops form, although they may be positioned randomly on flat leaves, they often tend to accumulate in the direction of the leaf tip as time continues. The mechanism of this behavior has been reported by Shanahan1 in terms of surface free energy minimization. Another question that naturally arises is: why exactly do the initial dew drop doses not preferably form at the leaf tip in the first place?

The formation of the initial dew drop is a typical heterogeneous nucleation process of liquid droplets from a bulk vapor phase, which is a typical process that exists in nature and industrial applications, such as the formation of rain drops and hailstones, crystal formation,2–5 chemical vapor deposition,6–8 and nucleation of initial droplets in dropwise condensation. As far as the nucleation process is concerned, the thermodynamic model was usually adopted and referred to as the Classic Nucleation Theory (CNT).9,10 Subsequently, the CNT model has been modified to obtain more accurate results,11–14 and the effect of a planar substrate and the heterogeneous nucleation model were considered.15–24 Experimental results revealed that the nucleation process can be affected by artificially distributed hydrophilic areas to realize the so called “spatial control” of nucleation sites.25 On the other hand, the relationship between the substrate structure and the nucleation process,26–33 and the effect of micro particles on the nucleation rate, were also widely investigated.34,35 The results suggest that the initial droplets tend to appear on the substrate defects and deposited heterogeneous particles,26,27,33 and the nucleation sites number density increases with the surface roughness.28,36 These results all indicated that the substrate properties, such as the wettability of the substrate material, and the distribution of microstructures, have a great effect on the nucleation of the initial nucleus. However, most of the works concentrated on the physico-chemical properties that are obtained from statistics, such as the apparent contact angle and surface roughness, and the underlying mechanism of how the substrate structures affect the nucleation process is still not very clear. Hence, a better understanding of the behavior of initial dew drop formation on different structures is very significant to explain natural phenomenon and predict the effect of substrate structures on nucleation processes.

In this paper, the behavior of heterogeneous nucleation of water droplets on conical microstructures will be explored. A circular conical apex and a cavity are adopted as the physical model to represent the typical defects widely existing on the substrates, and the CNT model is adopted to quantitatively analyze the nucleation capability of the microstructures under different condensation conditions.

2. Physical model

A circular conical apex and a cavity are shown in Fig. 1 to represent the defects which exist in actual substrates, and a planar substrate is also considered for comparison. The structures are characterized by the intrinsic wetting angle α of the substrate material and the cone angle β. As the initial nucleus is relatively small, gravitational force can be neglected as compared with the surface tension effect. According to the classic capillary approximation16 and molecular dynamics simulation results,37 the nucleus will be deposited on three different substrate structures in the form of a spherical cap as shown in Fig. 1, with the local contact angle equal to the intrinsic wetting angle of the substrate material. It has to be pointed out that, for structures that fulfill the condition of α + β/2 < 90°, the meniscus within micro cavities will be concave instead of that described in Fig. 1c. The additional pressure provided by the concave meniscus is favorable for the stability of embryos of any sizes and hence is favorable for the formation of the nucleus. Here, the term of embryo is adopted to denote a molecular cluster that has not reach up to the critical size, in order to differentiate from nucleus. Based on this fact, the structures fulfilling the α + β/2 < 90° condition will not be considered in this paper. The sizes of the microstructures are comparable to that of the initial nucleus, which is on the nanoscale under typical nucleation conditions of water droplets.
image file: c4ra12352b-f1.tif
Fig. 1 Schematic representation of an axisymmetric nucleus on three microstructures. (a) Circular conical apex, (b) planar substrate, (c) circular conical cavity.

In the following thermodynamic model section, it can be found that the nucleation process is greatly affected by the volume and liquid–vapor and liquid–solid interfacial areas of the initial nucleus. For a nucleus deposited on an apex as described by Fig. 1a, the related parameters can be calculated from the geometry:

 
R = r* × cos(αβ/2) (1)
 
hcap = r* × [1 + sin(αβ/2)] (2)
 
hcone = R × cot(β/2) (3)
where r* is the critical curvature radius of the initial nucleus, R is the bottom radius of the cone, and hcap and hcone are the heights of the spherical cap and the cone, respectively.

The volume of the nucleus Vdrop, the liquid–vapor and liquid–solid interfacial areas of Slv and Sls can be obtained as follows:

 
image file: c4ra12352b-t1.tif(4)
 
Slv = 2πr*2[1 + sin(αβ/2)] (5)
 
Sls = πr*2[thin space (1/6-em)]cos2(αβ/2)/sin(β/2) (6)

In the case of the planar substrate, β can be simply set as 180°. For a cavity structure, the cone angle β* can be converted into β by the formula β = 360 − β*, and the formulas above are also applicable. Consequently, eqn (1)–(6) can be used to calculate the parameters related to all three physical structures.

3. Thermodynamic model

Generally, the nucleation capability at certain conditions can be evaluated either by the nucleation free energy barrier ΔG(r*) or the nucleation rate J. Here, ΔG(r*) is defined as the Gibbs free energy barrier that has to be overcome to form a nucleus of critical size r*; J is defined as the number of initial nuclei formed within a unit time period and a unit volume or area for homogeneous or heterogeneous nucleation processes. Based on classic nucleation theory, the nucleation rate J can be expressed as follows:19,25
 
image file: c4ra12352b-t2.tif(7)
where J0 is a kinetic pre-factor that is directly connected with vapor conditions and nucleus configurations (spherical shape for homogeneous nucleation, and spherical cap for heterogeneous nucleation on a planar substrate).19,38 kB is the Boltzmann constant and T is the temperature.

Using the CNT model, the capillary approximation16,19,21 and Young’s equation, the nucleation free energy barrier of critical nucleus can be expressed as:

 
image file: c4ra12352b-t3.tif(8)
where ρ, σlv and M are the density, surface tension and molecular mass of the condensate liquid, respectively. According to the capillary assumption,16,19,21 ρ and σlv are obtained from the physico-chemical properties of the bulk liquid. NA is Avogadro’s number, and S is the supersaturation defined as the ratio between vapor pressure Pv and the equilibrium pressure at nucleation temperature.

The critical radius of the nucleus can be obtained as:25

 
image file: c4ra12352b-t4.tif(9)
where nl is the number of molecules per unit volume of condensate liquid.

By substituting eqn (9) into eqn (4)–(6) to calculate Vdrop, Slv and Sls, and then substituting them into eqn (8), a general formula of ΔG(r*) can be obtained as follows with respect to the microstructure configurations:

 
image file: c4ra12352b-t5.tif(10)
with F(α,β), the form factor that is connected with substrate structure parameters, given by the following expression:
 
image file: c4ra12352b-t6.tif(11)

For nucleation processes of liquid droplets from the bulk vapor phase, the nucleation rate J can be deduced from the following simplified general formula:21

 
image file: c4ra12352b-t7.tif(12)
where N(n) is the number density of embryos containing n molecules per unit area for heterogeneous nucleation, je(n) is the attachment rate of vapor molecules onto an embryo of size n, and A(n) is the vapor–liquid interfacial area of an embryo.

The integral formula in eqn (12) means that the nuclei are developed from embryos with various sizes. The number distribution of these embryos N(n) can be obtained from the classic cluster size distribution model:11,21,24

 
image file: c4ra12352b-t8.tif(13)
where ρN,v is the number of vapor molecules per unit volume of the bulk vapor. The pre-factor of the above formula is either ρN,v or (ρN,v)2/3 for homogeneous or heterogeneous nucleation.21 In the present study, the nucleation occurs on the substrate, and the pre-factor of (ρN,v)2/3 is used. The analysis is conducted under steady state conditions, and there are no depletion effects on the finite number of microstructures. ΔG(r) is the free energy change for the formation of embryo of radius r. The difference between ΔG(r) and ΔG(r*) is whether the critical radius is achieved.

The radius of the embryo r is related to the number of molecules n in the embryo by:

 
image file: c4ra12352b-t9.tif(14)
where m is the mass of one molecule, and vl is the specific volume of the condensate, which is also adopted from the bulk liquid.16,19,21

ΔG(r) can be obtained by expanding the right hand side term of eqn (10) in a Taylor power series in terms of rr* about the equilibrium radius r*:21

 
image file: c4ra12352b-t10.tif(15)

je(n) can be deduced from the kinetic theory of gases:

 
image file: c4ra12352b-t11.tif(16)

With the consideration of microstructure configurations, the interfacial area of A(n) can be expressed as:

 
A(n) = 2πr2[1 + sin(αβ/2)] (17)

Substituting eqn (13), (16) and (17) into eqn (12), the nucleation rate J is written:

 
image file: c4ra12352b-t12.tif(18)

Meanwhile, the relationship between cluster size n and r can be described by eqn (14). Substituting eqn (14) into the above formula, J can be organized as:

 
image file: c4ra12352b-t13.tif(19)

Substituting eqn (15) into the integration term of the above formula, we can get:

 
image file: c4ra12352b-t14.tif(20)

Setting intermediate variable B as:

 
image file: c4ra12352b-t15.tif(21)

eqn (20) then can be simplified as follows:

 
image file: c4ra12352b-t16.tif(22)

Finally, substituting eqn (22) into eqn (19), a general form of J can be obtained as:

 
image file: c4ra12352b-t17.tif(23)

If β is taken to be 180°, eqn (23) becomes identical to the expression for heterogeneous nucleation on a planar substrate.

Meanwhile, the kinetic pre-factor J0 in eqn (7) also can be obtained as:

 
image file: c4ra12352b-t18.tif(24)

As indicated by eqn (10) and (23), the nucleation rate J is a function of nucleation conditions (Pv and S) and substrate structure parameters (α,β).

4. Results and discussion

4.1 Effect of β on the nucleation rate

According to eqn (23), the nucleation rate of the heterogeneous nucleation process is determined by J0 and ΔG(r*). The calculated results of J0 and ΔG(r*) under various structure parameters are shown in Fig. 2. As β increases, the substrate structures translate from apexes to cavities, with J0 and ΔG(r*) decreasing rapidly in the same manner. As indicated by the schematic diagram of the physical model, a relatively larger liquid–vapor interfacial area is expected for an initial nucleus deposited on the top of an apex, providing a relatively higher probability for vapor molecular attachment. As a result, J0 for an apex is larger than that of planar substrate or cavities. For instance, the calculated results of J0 for an apex with β = 60° is 1.37 times that of a planar substrate (β = 180°) and 3.73 times that of a cavity with β = 300° (β* = 60°). On the other hand, as β increases, the volume and interfacial area of the initial nucleus decrease accordingly due to the space-confining effect of cavities. Considering that the main part of the nucleation free energy barrier is caused by the formation of new interfaces, the decrease of the interfacial area is thus favorable for the decrease of ΔG(r*). For instance, ΔG(r*) for an apex with β = 60° is 1.87 times that of a planar substrate and 13.93 times that of a cavity with β = 300°.
image file: c4ra12352b-f2.tif
Fig. 2 ΔG(r*), J0 and J under various structure parameters (Pv = 100 kPa, S = 1.5, α = 90°).

It is necessary to point out that, as the nucleation rate is an exponential function of ΔG(r*), the nucleation free energy barrier is thus decisive to the heterogeneous nucleation process. According to eqn (23), the nucleation rate for a cavity with β* = 60° is 9.0 × 1032 times that of a planar substrate and 1.6 × 1066 times that of an apex with β = 60° under the considered nucleation condition. As the microstructures translate from apexes to cavities, comparably lower nucleation free energy barriers are required to form the initial nucleus due to the space-confining effect of cavities, and the nucleation rates are thus increased. Based on the discussions above, the order of the nucleation capability of the three structure configurations is cavities, planar substrate and apexes. The presence of cavities with narrower cone angles is favorable for nucleation processes.

The nucleation rate J calculated from eqn (23) is shown in Fig. 3, under the nucleation condition of Pv = 100 kPa and S = 1.5. Apexes are not considered as they are not preferred for nucleation compared with a planar substrate and cavities. The structure parameter β is translated into β* using the relation of β* = 360 − β, and a smaller β* denotes a narrower cavity. The intrinsic wetting angle α is restricted within 60–110° with the consideration of most practical substrate materials.21 As expected, the nucleation rate decreases with α sharply, indicating that the nucleation rate on a hydrophilic surface was higher than a hydrophobic one for any substrate structures. Meanwhile, the nucleation rates for narrower cavities are obviously greater than the planar substrate for the same α, indicating relatively higher nucleation capabilities for micro cavities. It has been reported by Varanasi25 that the nucleation rate on a hydrophilic surface with α ∼ 25° is about 10129 times higher than that on a hydrophobic surface with α ∼ 110°, and the nucleation sites can be artificially controlled by a hydrophilic–hydrophobic hybrid surface. According to the model analysis above, so called spatial control of nucleation sites also can be realized by the appropriate substrate structure constructions, except for the regulation of surface wettability.


image file: c4ra12352b-f3.tif
Fig. 3 Nucleation rate as a function of wetting angle α and structure parameter β* (Pv = 100 kPa, S = 1.5).

It also can be found that the effect of structure parameter is so great that the nucleation rates under certain conditions are extremely low and almost no nuclei can be formed under these circumstances. This behavior actually provides a threshold of α and β for heterogeneous nucleation processes. In the contour map of Fig. 3, a threshold of 1 m−2 s−1 is chosen following Carey’s analysis.21 As a result, only those micro cavities that fulfill the threshold can be activated to form initial nuclei (see the lower-right part of the contour map with color fill).

For water vapor condensation, the kinetic pre-factor J0 is usually in the range of 1023–1026 m−2 s−1, which means that the exponential part of eqn (23) has to be greater than 10−26–10−23 to fulfill the threshold. The threshold proposed here actually provides an upper limit for the nucleation free energy barrier. To ensure effective nucleation, ΔG(r*) has to be low enough, and the upper limit of ΔG(r*) can be calculated from eqn (23) at different condensation conditions. In principle, a relatively lower wetting angle and narrower cavity are favorable for nucleation processes. This also explains why randomly distributed cavities, grooves, scratches (smaller β*) and heterogeneous particles (smaller α) can act as nucleation sites.26,27,39

One of the important conclusions obtained from the above analysis is that the number of active nucleation sites will be finite on a practical condensation substrate with randomly distributed micro cavities. The number of nucleation sites per unit surface area is usually defined as the nucleation sites number density (Ns), an important parameter in dropwise condensation heat transfer theory.40 Considering that the actual condensation substrate is composed of randomly distributed micro cavities, grooves and apexes of different structure parameters, the nucleation capabilities of different areas are inherently different depending on whether the above thresholds are well fulfilled. During the initial condensation stage, the micro cavities with higher nucleation capabilities will be rapidly occupied by initial nuclei. As condensation continues, more cavities that fulfill the threshold will be gradually activated, and the number of initial nuclei increases accordingly until all of the possible nucleation sites are occupied, yielding a maximum value of nuclei numbers. After that, the subsequent condensation process will be realized by the growth of pre-existing nuclei from critical size to micro droplets, while no nuclei could form on the blank surface between adjacent droplets. The blank area is inherently not preferred for nucleation basically due to the failure to fulfill the threshold.

4.2 Nucleation capability of micro cavities at various condensation conditions

As indicated by eqn (9) and eqn (23), the nucleation rate is also a function of vapor pressure and supersaturation. The calculated results of nucleation rate at various condensation conditions are presented in Fig. 4 and 5, with S varying from 1.2 to 1.6 under the vapor pressure of 100 kPa, and Pv varying from 10 kPa to 100 kPa under the supersaturation of 1.5, respectively. As described earlier, an arbitrary threshold of J = 1 m−2 s−1 is chosen to determine whether the micro cavity can be activated as an nucleation site, and the combination of α and β that fulfils the threshold is shown in Fig. 4 and 5 by color fill.
image file: c4ra12352b-f4.tif
Fig. 4 Nucleation capability of micro cavities under various supersaturations.

image file: c4ra12352b-f5.tif
Fig. 5 Nucleation capability of micro cavities under various vapor pressures.

As S increases from 1.2 to 1.6, the nucleation rate for the same structure increases accordingly. According to eqn (9), the critical size of the nucleus decreases sharply as S increases, yielding a comparably lower nucleation free energy barrier as indicated by eqn (10). Consequently, the nucleation rate is greatly increased for the same micro cavity. Meanwhile, some of the micro cavities that are originally not preferred for nucleation under low supersaturations can translate into effective nucleation sites when S is increased to a certain degree.

On the other hand, the nucleation rate also increases with Pv under a constant supersaturation, as shown in Fig. 5. According to eqn (9) and (10), as vapor pressure increases, the critical radius of the nucleus decreases slightly, yielding a relatively lower nucleation free energy barrier that is preferred for nucleation. Meanwhile, as indicated by eqn (24), the kinetic pre-factor J0 also increases with Pv. As a result, J increases rapidly with Pv for the same S and αβ. Similar to the behavior observed in Fig. 4, some of the micro cavities that are originally not preferred for nucleation under low vapor pressures can translate into effective nucleation sites when Pv is increased to a certain degree, as shown in Fig. 5. As Pv and S increase, more micro cavities with wider cone angles can translate into active nucleation sites, suggesting that Ns may increase with Pv and S accordingly. The relationship between Ns and the condensation conditions has been noticed by different researchers,41,42 and the results also suggest that Ns increases with Pv and S.

5. Conclusions

A circular conical apex and a cavity are proposed as physical models to represent the typical defects which widely exist on substrates, and classic nucleation theory is adopted to quantitatively analyze the nucleation capability of different microstructures at different water vapor nucleation conditions.

The results indicate that the kinetic pre-factor and nucleation free energy barrier all decrease when the substrate structures translate from apexes to cavities, and the nucleation rate of narrower cavities is higher than for a planar substrate and apexes. The cavities that are distributed on a substrate can act as nucleation sites, and the activated nucleation sites are finite for practical substrates due to the different nucleation capabilities of different surface areas. The spatial control of nucleation sites and the intensification of the nucleation sites number density can be realized by substrates with relatively lower wetting angle and the presence of narrower cavities.

The nucleation capability can also be affected by condensation conditions. As vapor pressure and supersaturation increase, the nucleation rate increases rapidly for the same structures, and some of the cavities that are originally not preferred for nucleation gradually translate into active nucleation sites, suggesting that the nucleation sites number density may increase with vapor pressure and supersaturation.

Acknowledgements

The authors greatly appreciate the financial support from the National Natural Science Foundation of China (no. 51236002, no. 51476018).

Notes and references

  1. M. E. R. Shanahan, Langmuir, 2011, 27, 14919–14922 CrossRef CAS PubMed.
  2. V. G. Karpov, Y. A. Kryukov, M. Mitra and I. V. Karpov, J. Appl. Phys., 2008, 104, 054507 CrossRef PubMed.
  3. A. T. Appapillai, C. Sachs and E. M. Sachs, J. Appl. Phys., 2011, 109, 084916 CrossRef PubMed.
  4. Z. Kožíšek and P. Demo, J. Cryst. Growth, 1998, 194, 239–246 CrossRef.
  5. S. Klein and D. M. Herlach, J. Appl. Phys., 2013, 114, 183510 CrossRef PubMed.
  6. S. A. Moshkalev and C. Verssimo, J. Appl. Phys., 2007, 102, 044303 CrossRef PubMed.
  7. Y. K. Chae and H. Komiyama, J. Appl. Phys., 2001, 90, 3610–3613 CrossRef CAS PubMed.
  8. J. Robertson, J. Appl. Phys., 2003, 93, 731–735 CrossRef CAS PubMed.
  9. M. Volmer and A. Weber, Z. Phys. Chem., 1926, 119, 277–301 CAS.
  10. J. C. Fisher, J. H. Hollomon and D. Turnbull, J. Appl. Phys., 1948, 19, 775–784 CrossRef CAS PubMed.
  11. A. Dillmann and G. E. A. Meier, J. Chem. Phys., 1991, 94, 3872–3884 CrossRef CAS PubMed.
  12. C. F. Delale and G. E. A. Meier, J. Chem. Phys., 1993, 98, 9850–9858 CrossRef CAS PubMed.
  13. I. J. Ford, A. Laaksonen and M. Kulmala, J. Chem. Phys., 1993, 99, 764–765 CrossRef CAS PubMed.
  14. V. I. Kalikmanov and M. E. H. Van Dongen, J. Chem. Phys., 1995, 103, 4250–4255 CrossRef CAS PubMed.
  15. D. Turnbull, J. Chem. Phys., 1950, 18, 198–203 CrossRef CAS PubMed.
  16. N. H. Fletcher, J. Chem. Phys., 1958, 29, 572–576 CrossRef CAS PubMed.
  17. A. Laaksonen, I. J. Ford and M. Kulmala, Phys. Rev. E: Stat. Phys., Plasmas, Fluids, Relat. Interdiscip. Top., 1994, 49, 5517–5524 CrossRef CAS.
  18. V. Talanquer and D. W. Oxtoby, Phys. A, 1995, 220, 74–84 CrossRef.
  19. D. Kashchiev, Nucleation: Basic Theory with Applications, Butterworth-Heinemann, Burlington, MA, 2000, ch. 3, pp. 32–33, ch. 9, pp. 124–125 Search PubMed.
  20. M. Qian, Acta Mater., 2007, 55, 943–953 CrossRef CAS PubMed.
  21. V. P. Carey, Liquid-vapor Phase-Change Phenomena, Taylor and Francis, New York, USA, 2008 Search PubMed.
  22. C. M. L. V. Kroon and I. J. Ford, Atmos. Res., 2011, 101, 553–561 CrossRef PubMed.
  23. A. S. Abyzov and J. W. P. Schmelzer, J. Chem. Phys., 2013, 138, 164504 CrossRef PubMed.
  24. W. Xu, Z. Lan, B. L. Peng, R. F. Wen and X. H. Ma, RSC Adv., 2014, 4, 31692–31699 RSC.
  25. K. K. Varanasi, M. Hsu, N. Bhate, W. Yang and T. Deng, Appl. Phys. Lett., 2009, 95, 094101 CrossRef PubMed.
  26. J. L. McCormick and J. W. Westwater, Chem. Eng. Sci., 1965, 20, 1021–1036 CrossRef CAS.
  27. T. Y. Song, Z. Lan, X. H. Ma and T. Bai, Int. J. Therm. Sci., 2009, 48, 2228–2236 CrossRef CAS PubMed.
  28. C. F. Mu, J. J. Pang, Q. Y. Lu and T. Q. Liu, Chem. Eng. Sci., 2008, 63, 874–880 CrossRef CAS PubMed.
  29. X. M. Chen, J. Wu, R. Y. Ma, M. Hua, N. Koratkar, S. H. Yao and Z. K. Wang, Adv. Funct. Mater., 2011, 21, 4617–4623 CrossRef CAS.
  30. E. Ruckenstein and G. O. Berim, J. Colloid Interface Sci., 2010, 351, 277–282 CrossRef CAS PubMed.
  31. M. Qian and J. Ma, J. Cryst. Growth, 2012, 355, 73–77 CrossRef CAS PubMed.
  32. K. Y. Li, S. Xu, W. X. Shi, M. He, H. L. Li, S. Z. Li, X. Zhou, J. J. Wang and Y. L. Song, Langmuir, 2012, 28, 10749–10754 CrossRef CAS PubMed.
  33. C. W. Lo, C. C. Wang and M. C. Lu, Adv. Funct. Mater., 2014, 24, 1211–1217 CrossRef CAS.
  34. M. Lazaridis and Ø. Hov, Atmos. Res., 2000, 55, 103–113 CrossRef CAS.
  35. V. Y. Smorodin and P. K. Hopke, Atmos. Res., 2006, 82, 591–604 CrossRef CAS PubMed.
  36. I. O. Ucar and H. Y. Erbil, Appl. Surf. Sci., 2012, 259, 515–523 CrossRef CAS PubMed.
  37. A. Shahraz, A. Borhan and K. A. Fichthorn, Langmuir, 2013, 29, 11632–11639 CrossRef CAS PubMed.
  38. J. W. P. Schmelzer, J. Non-Cryst. Solids, 2010, 356, 2901–2907 CrossRef CAS PubMed.
  39. K. Rykaczewski, Langmuir, 2012, 28, 7720–7729 CrossRef CAS PubMed.
  40. S. Kim and K. J. Kim, J. Heat Transfer, 2011, 133, 081502 CrossRef.
  41. C. Graham and P. Griffith, Int. J. Heat Mass Transfer, 1973, 16, 337–346 CrossRef CAS.
  42. J. W. Rose, Int. J. Heat Mass Transfer, 1976, 19, 1363–1370 CrossRef CAS.

This journal is © The Royal Society of Chemistry 2015
Click here to see how this site uses Cookies. View our privacy policy here.