The present invention relates to heterogeneous surfaces.
Superhydrophobic surfaces have received significant interest for dropwise condensation to increase the efficiency of energy applications such as heat exchangers, power plants, and solar thermal energy conversion systems. However, nucleation densities on regular superhydrophobic surfaces are difficult to achieve due to the high energy barrier for nuclei formation and hence enhancement of heat transfer can be limited.
In general, heterogeneous surface structures can be made by infusing microstructured surfaces with low-surface tension oil, which lead to nucleation densities that were increased by over an order of magnitude while maintaining low droplet adhesion. The approach offers a simple and scalable approach to create surfaces that can be tailored for enhanced heat transfer.
In one aspect, a superhydrophobic surface can include a patterned substrate having a surface including a plurality of first regions distributed in a second hydrophobic region, the first regions including a surface modifying layer and the second hydrophobic region including a material infused into regions of the substrate. The first regions can have hydrophobic features and hydrophilic features.
In another aspect, a method of increasing nucleation density on a surface can include infusing a material into regions of a patterned substrate to form a surface including a plurality of first regions distributed in a second hydrophobic region, the first regions including a surface modifying layer and the second hydrophobic region including a material infused into regions of the substrate.
In certain embodiments, the patterned substrate can include a periodic structure on the surface of the substrate that form the regions of the substrate into which the material is infused. The first regions can be associated with the pattern of the patterned substrate. The periodic structure can be a micropillar or microcolumn.
In certain embodiments, the surface modifying layer can include a functionalized silane. In certain circumstances, the surface modifying layer can include a plurality of scattered hydrophilic sites while exhibiting overall hydrophobicity.
In other embodiments, the material can be an oil or wax. The oil can be a fluorinated oil.
Other features, objects, and advantages will be apparent from the description and drawings, and from the claims.
Condensation heat transfer has wide applications in various systems such as heat exchangers, heat pipes and power plants. The heat transfer coefficient of condensation is of great significance to the efficiency of such systems. Dropwise condensation, where the condensate forms discrete droplets rather than continuous films covering the substrate, is considered as one of the most promising approaches to enhance the heat transfer coefficient. Previous work has demonstrated the application of nanostructured superhydrophobic surfaces where condensate can be spontaneously removed via a surface-tension-driven mechanism. See, for example, J. B. Boreyko and C.-H. Chen, “Self-Propelled Dropwise Condensate on Superhydrophobic Surfaces”, PRL, 2009. 103(18): p. 184501, which is incorporated by reference in its entirety. However, the nucleation density on these surfaces is relatively low since the phase change process relies on high energy active sites to initiate nucleation at low supersaturations (low ΔT), limiting the overall heat transfer performance. Furthermore, air pockets trapped beneath the droplets during growth reduce the contact area between the condensing droplet and substrate, which increases the thermal resistance and reduces the heat transfer coefficient. See, for example, N. Miljkovic, R. Enright, and E. N. Wang, “Effect of Droplet Morphology on Growth Dynamics and Heat Transfer during Condensation on Superhydrophobic Nanostructured Surfaces”, ACS Nano, 2012 6(2): p. 1776-1785, which is incorporated by reference in its entirety.
Recently, Wong et al. demonstrated a liquid-solid composite surface created by infusing a porous fluoropolymer with water-immiscible, low-surface-tension Krytox oil. See, for example, T.-S. Wong, S. H. Kang, S. K. Y. Tang, E. J. Smythe, B. D. Hatton, A. Grinthal and J. Aizenberg, “Bioinspired self-repairing slippery surfaces with pressure-stable omniphobicity”, Nature, 2011. 477: p. 443-447, which is incorporated by reference in its entirety. On such a composite surface, the contact area between droplet and substrate can be large while contact line pinning remains very low allowing easy removal of droplets. These properties make the surfaces potentially suitable for enhanced condensation heat transfer. While such behavior is possible with the proper choice of silane, e.g., dichlorodimethylsilane on SiO2, the nucleation density is limited by the presence of high surface energy defects and contaminants at low supersaturations (low ΔT). The nucleation density on oil-infused, silane-coated structured surfaces can be significantly increased by the use of disordered long-chain silane coatings that result in nucleation sites limited only by the density of pillar structures comprising the surface. The increase in the nucleation density can be explained by heterogeneity in the surface energy of the silane coating and the reduced water-oil interfacial energy. This effect could potentially be used to significantly improve the heat transfer coefficient in condensation by controlling the nucleation density.
The surfaces described herein can increase the nucleation density during condensation process while maintaining the easy removal of condensate to enhance condensation rates. In order to achieve this, the surface can include three components. The first is a surface coating which is overall hydrophobic with local nanometer or micrometer scale hydrophilic sites, which create a heterogeneous surface structure. The purpose of the heterogeneity is to provide nucleation sites for condensation to happen while the overall hydrophobicity allows the easy removal of the condensate. The second is a filling fluid which is immiscible with the condensate and has low interfacial tension with the condensate. The filling serves for two purposes: providing a reduced interfacial tension between the condensate and the oil to reduce the energy cost of condensation; and help to remove the condensate from the substrate. The third component is micrometer or nanometer scale roughness to enhance the wetting of the filling fluid.
In particular, a superhydrophobic surface can be formed on a substrate from a pattern on the surface of a substrate. The pattern can be formed from a plurality of structures on the substrate. The structures can have nanometer sized or micrometer sized features. The features can be bumps, columns, pillars, channel, or trough. The features can be periodically spaced on the surface. For example, the features can be less than 10 micrometers, less than 5 micrometers, less than 2 micrometers, less than 1 micrometer, less than 0.5 micrometers, or less than 0.1 micrometers in width. The features can be spaced in intervals of about 0.5 micrometer, 1 micrometer, 2 micrometers, 5 micrometer, 10 micrometers or 20 micrometers, or more, from each other on the surface. The features can have a height of about 0.5 micrometer, 1 micrometer, 2 micrometers, 5 micrometer, 10 micrometers or 20 micrometers. In certain circumstances, the features are etched or machined from the substrate. In other circumstances, one or more feature can be grown or deposited on the surface of the substrate.
Once the patterned substrate has been formed, the surface can be coated with one or more coating layers. The coating layer can be selected to impart desired properties on the surface, such as, for example, mechanical robustness or increased hydrophobicity, or both. For example, the superhydrophobic surface can include a surface modifying layer on at least a portion of the nanostructures. The surface modifying layer can be a single layer or a multilayer. For example, an initial coating layer, e.g., a metallic layer can be deposited by (for example) electroless plating, chemical vapor deposition or atomic layer deposition. The initial coating layer can be a polymer or a metal. The surface modifying layer can be a hydrophobic material, such as a polymer or self-assembled monolayer, directly on the nanostructure or on the initial coating layer. For example, a silane or a thiol can be assembled on a surface. The hydrophobic material; e.g., a hydrophobic polymer, hydrophobic thiol, hydrophobic carboxylic acid or hydrophobic silane, can include hydrocarbon (e.g., a saturated hydrocarbon) groups, halohydrocarbon groups (e.g., a saturated fluorohydrocarbon), or halocarbon groups (e.g., a perfluorinated alkyl group). In certain examples, the hydrophobic material can be trichloro(1H,1H,2H,2H-perfluorooctyl) silane, (tridecafluoro-1,1,2,2-tetrahydrooctyl)-1-trichlorosilane, (1H,1H,2H,2H-perfluorodecyl acrylate), a Teflon amorphous fluoropolymer resin, or an alkyl or fluoroalkyl thiol deposited by appropriate techniques. The hydrophobic material can have C2-C18 groups that can be fluorinated to varying degrees. The trifluoromethyl or difluoromethyl groups on the surface can allow the surface properties to be tailored, for example, to have advancing wetting angles of 105, 110, 115 or 120 degrees, depending on the choice of fluorinated alkyl group and base structure. The coatings can have a plurality of hydrophilic sites scattered in the overall hydrophobic background. The size of the hydrophilic sites can be 10 nanometer, 100 nanometer or 500 nanometer. The fraction of hydrophilic sites can be 1%, 5%, 10% or 40%. The local contact angles of the hydrophilic sites can be 50-70 degrees.
The surface modified patterned substrate can then be infused with a material. The infusion of the material can include drop coating, dip coating or roll coating the surface with the material. The material can be an oil, for example a fluorinated oil, or low melting point solid, such as a wax. The material can be a low surface tension material, which can allow it to infuse the pattern readily and completely.
For example,
Oil-impregnated surfaces have been recently reported as a promising approach to enhance condensation heat transfer surfaces due to the ultra-low droplet adhesion. However, easy droplet removal is not the only desired property for high heat transfer performance. Low contact angle and high nucleation densities are also essential to further enhance condensation heat transfer. By combining surface heterogeneity and oil-infusion, the nucleation density in condensation can be increased by over an order of magnitude via immersion condensation while maintaining low droplet adhesion. The increase in nucleation densities via the a combined effect of heterogeneity and the reduced oil-water interfacial tension was explained by this disclosure based on classical nucleation theory, which were also corroborated with control experiments using silane-coated silicon micropillar arrays. With improved understanding of the physics, oil-infused superhydrophobic copper oxide surfaces as a platform for condensation enhancement in practical systems were investigated. The condensation heat transfer coefficient on such oil-infused heterogeneous surfaces can be enhanced by approximately 100% compared to state-of-the-art dropwise surfaces in the presence of non-condensables gases. An order of magnitude increase in nucleation density could contribute to approximately 80% increase in the overall heat transfer coefficient. Meanwhile, the low departure radii and low contact angle also assisted in the total improvement. Achieving the three key aspects of condensation simultaneously can be important to realize heat transfer enhancement by as high as 100%. Further work is needed to tailor oil and coating properties, as well as surface geometry to minimize oil loss during operation and maximize condensing surface area. With continued development, immersion condensation promises to be an important condensation mode for a variety of heat transfer and resource conserving applications.
Fabrication
Well-defined silicon micro/nanopillar arrays with diameters, d, ranging from 0.4 μm to 5 μm, periods, l, ranging from 4 μm to 25 μm, and heights, h ranging from 10 μm to 25 μm were used in these experiments. The silicon surfaces were functionalized with three different chemicals: 1) (Tridecafluoro-1,1,2,2-tetrahydrooctyl)-1-trichlorosilane (TFTS) (UCT Specialties), which forms a self-assembled coating (SAC) by chemical vapor deposition (CVD) with a relatively long carbon chain (MW=481.54 g/mol), 2) Dimethyldichlorosilane (DMCS) (Sigma-Aldrich), which forms a self-assembled monolayer (SAM) by CVD with a short carbon chain (MW=129.06 g/mol), and 3) Poly(1H,1H,2H,2H-perfluorodecyl acrylate) (PFDA) polymer, which was deposited using initiated chemical vapor deposition (iCVD) with a typical film thickness of 35 nm. Goniometric measurements on smooth functionalized silicon surfaces showed advancing and receding contact angles of: θa/θr=122°±1.3°/78°±1.3° (equilibrium contact angle θe≈102.1°±0.9°); θa/θr=103.8°±0.5°/102.7°±0.4° (θe≈103.2°±0.3°); and θa/θr=121.1°±2.2°/106.3°±2.4° (θe≈113.5°±1.6°) for deposited films of TFTS, DMCS, and PFDA respectively. A small droplet of Krytox GPL 100 oil (DuPont) was applied to the functionalized silicon pillar arrays. The surface tension of Krytox oil is ˜17-19 mN/m, allowing the oil to spread on the surface. A dry nitrogen stream was used to assist spreading and remove excess oil. Typical scanning electron micrographs (SEM) of the silicon pillar arrays without and with the oil are as shown in
The nucleation behavior on the surfaces with and without the oil were investigated under white light optical microscopy (OM). The samples were horizontally mounted on a thermal stage (Instec Inc.) inside an enclosure and cooled to Tw=283.1±0.1K in a dry nitrogen atmosphere. Following thermal equilibration (˜5 min), nucleation was initiated by flowing water-saturated nitrogen into the enclosure. The supersaturation, defined as the ratio of the vapor pressure to the saturation pressure at the stage temperature (pv/pw), was controlled by the temperature of the water reservoir through which the nitrogen carrier gas was sparged and measured using a humidity probe (Hygroclip, Rotronic) located ˜1 cm above the sample. Typical values of the supersaturation were around S≈1.6. The nucleation density and subsequent growth behavior was recorded at a frame rate of 30 fps using a CMOS camera (Phantom V7.1, Vision Research) attached to the optical microscope.
Referring to
As shown in
Meanwhile, on pillar arrays coated with DMCS and PFDA, similar increases in nucleation density were not observed, as shown in
AFM Imaging and Contact Angle Analysis
To investigate the mechanism for this drastic change in nucleation density, atomic force microscopy (AFM) was performed in tapping mode on a smooth TFTS-coated silicon surface and observed the presence of micelle structures, as shown in
The local contact angles on the hydrophobic substrate and the micelle structures can be determined based on a modified Cassie-Baxter model. Assuming the local contact angles on the hydrophobic substrate and the micelle structures to be θ1 and θ2, respectively, the macroscopic advancing and receding contact angles are determined as
θa=θ1, (1)
cos θr=√{square root over (f)} cos θ2+(1−√{square root over (f)})cos θ1, (2)
where f is the area fraction of the micelles.
Based on the macroscopically measured advancing and receding angles, θa=122°±1.3° and θr=78°±1.3°, and the fraction of the micelles determined as f≈0.4 from AFM, the local contact angles on the hydrophobic substrate and the hydrophilic micelles were found to be θ1=122°±1.3° and θ2=60°±1.5°.
Nucleation Theory
The nucleation rate, J, can be determined by classical nucleation theory (CNT) as (D. Kashchiev, Nucleation: Basic Theory with Applications. 1 ed. 2000, Oxford: Butterworth-Heinemann, which is incorporated by reference in its entirety)
J=zf*exp(−G*) (3)
In Eqn. (3), z is the Zeldovich factor and G* is the dimensionless energy barrier, given by
z=(kT ln S)2/8πνo√{square root over (kTψ(θ)γ2)} (4)
G*=16πψ(θ)νo2γ3/3(kT)2(ln S)2 (5)
where S is the supersaturation and ψ(θ) is the activity that accounts for the effect of contact angle. f* is the frequency of monomer attachment to the critical droplet nucleus dependent on the nature of the nucleus growth. The main modes of growth during heterogeneous nucleation are limited via surface diffusion or direct impingement of monomers to the nucleus. See, for example, G. M. Pound, M. T. Simnad, and L. Yang, “Heterogeneous nucleation of crystals from vapor” J. Chem. Phys., 1954. 22(1215) and R. A. Sigsbee, “Atom capture and growth rates of nuclei”, JAP, 1971. 42(10): p. 3904-3915, each of which is incorporated by reference in its entirety. Volumetric diffusion is a third growth limiting step, which is only considered important for nucleation taking place in liquid or solid solutions. See, for example, D. Kashchiev, Nucleation: Basic Theory with Applications. 1 ed. 2000, Oxford: Butterworth-Heinemann, which is incorporated by reference in its entirety. However, all three mechanisms were included when calculating the nucleation rates.
The frequency of monomer attachment due to direct vapor impingement is given by
f*I=γn[(1−cos(θw))/2φ2/3(θ)](36πνo2)1/3 ln2/3 (6)
where γn is the sticking coefficient (0<γn<1), I is the classical Hertz-Knudsen impingement rate (I=P/√{square root over (2πmokT)})), n is the number of molecules in the nucleated cluster, and νo is the volume of an individual water molecule (νo=3×10−29 m3). To determine an upper bound on the nucleation rate, a sticking coefficient of one was assumed (γn=1).
The frequency of monomer attachment due to surface diffusion is given by f*sd=γnc*λs2l, where c* is the capture number due to surface diffusion (1<c*<5), and λs is the mean surface diffusion distance of an adsorbed monomer on the substrate. The capture number c* is size independent and approximately equal to 1.9 for heterogeneous condensation of water vapor. See, for example, D. J. Pocker, and S. J. Hruska, “Detailed calculations of the number of distinct sites visited in random walk on several two-dimensional substrate lattices.” J. Vac. Sci. Tech., 1971. 8(6): p. 700-707, which is incorporated by reference in its entirety. The mean surface diffusion distance is dependent on the wettability of the substrate and is given by λs=√{square root over (Dsdτd)} where Dsd is the surface diffusion coefficient (Dsd=ds2νsexp[−Esd/kT]), τd is the desorption time (τd=(1/νs)exp[−Edes/kT]), νs is the adsorbed molecule vibration frequency determined using the Debye approximation (νs=VDα/2), ds is the length of a molecular jump along the substrate surface approximated by the lattice constant of the substrate (ds=5.4 Å) (J. P. Hirth, and G. M. Pound, Condensation and evaporation—nucleation and growth kinetics. 1963, England: Pergamon Press, which is incorporated by reference in its entirety) and VD is the speed of sound in the substrate (VD=8433 m/s). The desorption and surface diffusion energies are given by Edes=EI+σsvao and Esd=0.5Edes (P. A. Thiel, and T. E. Madey, “The interaction of water with solid surfaces: Fundamental aspects.” Surface Science Reports, 1987. 7(6-8): p. 211-385, which is incorporated by reference in its entirety), respectively, where EI is the binding energy of an n=1 sized cluster, σsv is the solid vapour interfacial energy and ao is the water molecule surface area (ao=4.67×10−19 m2). The calculated energies of desorption show excellent agreement with that of experiment and molecular dynamics simulations (Edes,SiO2=0.9 eV). See, for example, J. N. Israelachvili, Intermolecular and surface forces. 2nd ed. 1991, Amsterdam: Academic Press and Y. Ma, A. S. Foster, and R. M. Nieminen, “Reactions and clustering of water with silica surface”. J. Chem. Phys., 2005. 122(144709), which is incorporated by reference in its entirety.
The frequency of monomer attachment due to volumetric diffusion is given by
where D is the self diffusion coefficient of water vapor (D=(3/8πnodo2)√{square root over (kT/πmo)}), C is the equilibrium concentration of monomers (C=(1/ao)exp(−W1/kT)), do, mo, and no are the water molecule diameter (do=3.0 Å) (J. N. Israelachvili, Intermolecular and surface forces. 2nd ed. 1991, Amsterdam: Academic Press, which is incorporated by reference in its entirety), mass (mo=3×10−26 kg) (D. Kashchiev, Nucleation: Basic Theory with Applications. 1 ed. 2000, Oxford: Butterworth-Heinemann, which is incorporated by reference in its entirety) and number density (no=NA/νM), respectively.
By adding the nucleation rate from the three mechanisms together, the nucleation rate, J, can be determined as a function of the contact angle and surface tension of the condensate at given supersaturations, as shown in
The surface tension of water in air is 72 mN/m and the interfacial tension between water and the Krytox oil was found to be 58 mN/m by measuring the contact angle of water droplet on oil film. From
In order to validate this assumption, condensation experiments were carried out using silicon pillar arrays coated with 3-(trimethoxysilyl)propyl methacrylate (3-TMPM). The advancing and receding contact angle of water on a smooth silicon surface coated with 3-TMPM are 65°±1.5° and 53°±1.1°, respectively. The contact angle is in the range where the nucleation rate will be almost zero with a surface tension of 72 mN/m and nucleation should occur on every tip of pillars with an interfacial tension of 58 mN/m. Optical images of the condensation experiments are as shown in
Heat Transfer Coefficient
In practice, filmwise condensation, where a thin liquid film covers the surface, is the most prevalent condensation mode due to the high wettability of common heat transfer materials. In this condensation mode, the heat transfer coefficient is limited by the thermal resistance associated with the condensate film which insulates the surface. See Mills, A. F. Heat and Mass Transfer. 2 edn, (Prentice-Hall, 1999), which is incorporated by reference in its entirety. Accordingly, efforts spanning eight decades have been devoted to the realization of non-wetting surfaces for dropwise condensation where shedding droplets clear the surface for droplet re-nucleation/re-growth, leading to enhanced heat transfer rates. See Schmidt, E., Schurig, W. & Sellschopp, W. Versuche über die Kondensation von Wasserdampf in Film- and Tropfenform. Forschung im Ingenieurwesen 1, 53-63, (1930), Tanner, D. W., Potter, C. J., Pope, D. & West, D. Heat transfer in dropwise condensation—Part I The effects of heat flux, steam velocity and non-condensable gas concentration. International Journal of Heat and Mass Transfer 8, 419-426, (1965), O'Neill, G. A. & Westwater, J. W. Dropwise condensation of steam on electroplated silver surfaces. International Journal of Heat and Mass Transfer 27, 1539-1549, (1984), Boreyko, J. B. & Chen, C.-H. Self-Propelled Dropwise Condensate on Superhydrophobic Surfaces. Phys Rev Lett 103, 184501 (2009), Chen, C.-H. et al. Dropwise condensation on superhydrophobic surfaces with two-tier roughness. Appl Phys Lett 90, 173108-173103 (2007), and Le Fevre, E. J. & Rose, J. W. An experimental study of heat transfer by dropwise condensation. International Journal of Heat and Mass Transfer 8, 1117-1133, (1965), each of which is incorporated by reference in its entirety. One order of magnitude higher heat transfer coefficients compared to filmwise condensation have been reported using dropwise condensation in pure vapor environments. See Daniel, S., Chaudhury, M. K. & Chen, J. C. Fast Drop Movements Resulting from the Phase Change on a Gradient Surface. Science 291, 633-636, (2001), which is incorporated by reference in its entirety. In order to maximize the heat transfer coefficient, a high performance dropwise condensation surface should simultaneously achieve three properties: low contact angle hysteresis to minimize droplet departure radii, low contact angle to reduce the conduction resistance of the droplet, and high nucleation density (see Miljkovic, N., Enright, R. & Wang, E. N. Effect of Droplet Morphology on Growth Dynamics and Heat Transfer during Condensation on Superhydrophobic Nanostructured Surfaces. Acs Nano 6, 1776-1785, (2012), which is incorporated by reference in its entirety), as shown in
Immersion condensation, a new approach to enhance condensation heat transfer by introducing heterogeneous surface chemistry composed of discrete hydrophilic domains on a hydrophobic background in oil-infused micro and nanostructured surfaces is disclosed. This approach allows water droplets to nucleate immersed within the oil to achieve high nucleation densities while maintaining easy droplet removal and low contact angles (
Surface Heterogeneity by Self-Assembled Coatings
A self-assembled coating (SAC) of (tridecafluoro-1,1,2,2-tetrahydrooctyl)-1-trichlorosilane (TFTS) was deposited from the vapor phase (See Methods for the deposition process). The SAC coating method is capable of forming heterogeneity by agglomeration. See Bunker, B. C. et al. The Impact of Solution Agglomeration on the Deposition of Self-Assembled Monolayers. Langmuir 16, 7742-7751, (2000), which is incorporated by reference in its entirety. The SAC method was chosen due to its simplicity and scalability, but alternative methods are also available to generate heterogeneity at the appropriate length scale, e.g., block copolymer or nano-imprinting. See Park, M., Harrison, C., Chaikin, P. M., Register, R. A. & Adamson, D. H. Block Copolymer Lithography: Periodic Arrays of ˜1011 Holes in 1 Square Centimeter. Science 276, 1401-1404, (1997) and Guo, L. J., Cheng, X. & Chou, C.-F. Fabrication of Size-Controllable Nanofluidic Channels by Nanoimprinting and Its Application for DNA Stretching. Nano Lett 4, 69-73, (2003), each is which is incorporated by reference in its entirety. Height and phase atomic force microscope (AFM) images of the TFTS coating on a smooth silicon surface were obtained and are shown in
Immersion Condensation on Silicon Micropillars
The SAC was deposited on silicon micropillar arrays to fundamentally investigate nucleation behavior on oil-infused surfaces. Silicon micropillar arrays were fabricated with diameters, d, ranging from 0.4-5 μm, periods, l, ranging from 4-25 μm, and heights, h, ranging from 10-25 μm using contact lithography and deep reactive ion etching (DRIE) processes. The geometries were chosen to satisfy the imbibition condition to enable oil spreading and to stabilize the oil film. See Bico, J., Thiele, U. & Quéré, D. Wetting of textured surfaces. Colloids and Surfaces A: Physicochemical and Engineering Aspects 206, 41-46, (2002), which is incorporated by reference in its entirety. The pillar surfaces were subsequently functionalized with the TFTS SAC, and infused with a fluorinated oil, Krytox GPL 100. The low surface tension of Krytox oil (≈17-19 mN/m) allowed it to spread on the surface and form a stable film via capillarity. A dry N2 stream was used to assist spreading and remove excess oil. Typical scanning electron microscope (SEM) images of the coated pillar arrays without and with oil-infusion are shown in
The increase in nucleation density on the oil-infused TFTS surfaces was achieved via the combination of the high-surface-energy sites and reduced water-oil interfacial energy. Based on classical nucleation theory, the nucleation rate can be determined as a function of the contact angle and the surface energy of the condensate at a given supersaturation, as shown in
Immersion Condensation on Scalable Copper Oxide Nanostructures
The overall heat transfer performance of an immersion condensation surface is disclosed. While studies on well-defined silicon micropillar arrays can provide physical insight into immersion condensation behavior, they are not practical due to cost and challenges in interfacing the silicon substrate and the heat transfer measurement apparatus with minimum uncertainties. Therefore, immersion condensation heat transfer measurements on oil-infused copper oxide (CuO) nanostructures functionalized with TFTS was performed, which promises a scalable, low cost platform for condensation surfaces. See Nam, Y. & Ju, Y. S. Comparative Study of Copper Oxidation Schemes and Their Effects on Surface Wettability. ASME Conference Proceedings 2008, 1833-1838 (2008), which is incorporated by reference in its entirety. SEM images of representative copper oxide nanostructures without and with Krytox oil-infusion are shown in
Overall heat transfer coefficients were measured to evaluate the performance on three different CuO-based surfaces: a hydrophobic surface for typical dropwise condensation, a superhydrophobic TFTS-coated copper oxide surface, and a Krytox oil-infused, TFTS-coated CuO surface (
Methods
Surface Fabrication
The silicon micropillar arrays were fabricated using contact lithography followed by deep reactive ion etching. For copper oxide surfaces, commercially available oxygen-free Cu tubes (99.9% purity) with outer diameters, DOD=6.35 mm, inner diameters, DID=3.56 mm, and lengths L=131 mm as the test samples were used for the experiments. Each Cu tube was cleaned in an ultrasonic bath with acetone for 10 minutes and rinsed with ethanol, isopropyl alcohol and de-ionized (DI) water. The tubes were then dipped into a 2.0 M hydrochloric acid solution for 10 minutes to remove the native oxide film on the surface, then triple-rinsed with DI water, and dried with clean nitrogen gas.
Nanostructured CuO films were formed by immersing the cleaned tubes into a hot (96±3° C.) alkaline solution composed of NaClO2, NaOH, Na3PO4.12H2O, and DI water (3.75:5:10:100 wt. %). See Enright, R., Dou, N., Miljkovic, N., Nam, Y. & Wang, E. N. Condensation on Superhydrophobic Copper Oxide Nanostructures. 3rd Micro/Nanoscale Heat & Mass Transfer International Conference (2012), which is incorporated by reference in its entirety. During the oxidation process, a thin (<200 nm) Cu2O layer was formed that then re-oxidized to form sharp, knife-like CuO structures with heights of h≈1 μm, solid fraction φ≈0.023 and roughness factor r≈10. To verify the independence of oxide thickness on chemical oxidation time (see Nam, Y. & Ju, Y. S. Comparative Study of Copper Oxidation Schemes and Their Effects on Surface Wettability. Imece 2008: Heat Transfer, Fluid Flows, and Thermal Systems, Vol 10, Pts a-C, 1833-1838 (2009), which is incorporated by reference in its entirety), four separate samples were made using oxidation times, τ=5, 10, 20, and 45 minutes. The sharp CuO structures were then coated with silane SAC to create SHP surfaces.
In addition to SHP surfaces, cleaned copper tubes were also immersed into hydrogen peroxide solutions at room temperature to form a thin smooth layer of Cu2O. The smooth surfaces were also coated with TFTS to achieve typical hydrophobic surfaces for dropwise condensation (DHP).
Surface Coating Deposition
The self-assembled coatings (SAC) were formed using a vapor deposition process. First, the silicon surfaces were cleaned using a Piranha solution (H2O2:H2SO4=1:3) to remove possible organic contamination and to create a large number of OH bonds on the surface, which enables the bonding between silane molecules and the silicon surface. For the copper oxide surfaces, the surfaces were cleaned by intensive plasma (≈1 hr). The samples were then placed in a desiccator (Cole-Palmer) together with a small petri dish containing ≈1 mL of the silane liquid. The desiccator was pumped down to ≈10 kPa. The pump was then shut off and the valve was closed so that the silane liquid could evaporate in the low-pressure environment of the desiccator and attach to the surfaces to form the SAC via the following reaction,
Si—OH+R—Si—Cl→Si—O—Si—R+HCl.
During the self-assembly process, the silane molecule form nanoscale agglomerates with diameters of ≈200-500 nm shown in
Surface Characterization
Advancing and receding contact angles for all samples were measured and analyzed using a micro-goniometer (MCA-3, Kyowa Interface Science Co., Japan). Field emission electron microscopy was performed on a Zeiss Ultra Plus FESEM (Carl Zeiss GMBH) at an imaging voltage of 3 kV.
OM Imaging Procedure
The samples were horizontally mounted on a thermal stage inside an enclosure and cooled to Tw=283.1±0.1 K in a dry nitrogen atmosphere. Following thermal equilibration (≈5 minutes), nucleation was initiated by flowing water-saturated nitrogen into the enclosure. The humidity of the gas flow was measured using a humidity probe located 1 cm above the sample to determine the supersaturation, S, defined as the ratio of the vapor pressure to the saturation pressure at the stage temperature (S=pv/pw). Typical values of supersaturation were S≈1.6. The nucleation density and subsequent growth behavior was recorded at a frame rate of 10 frames per second using a high speed camera (Phantom V7.1, Vision Research) attached to the optical microscope. The observable nucleation density during each experiment was determined by counting the number of nuclei in the captured images and dividing the number of nuclei by the imaging area. Multiple experiments were performed to determine the average nucleation densities on the different surfaces.
ESEM Imaging Procedure
Condensation nucleation and growth were studied on these fabricated surfaces using an environmental scanning electron microscope (EVO 55 ESEM, Carl Zeiss GMBH). Backscatter detection mode was used with a high gain. The water vapor pressure in the ESEM chamber was 800±80 Pa. Typical image capture was obtained with a beam potential of 20 kV and variable probe current depending on the stage inclination angle. To limit droplet heating effects, probe currents were maintained below 2.0 nA and the view area was kept above 400 μm×300 μm. See Rykaczewski, K., Scott, J. H. J. & Fedorov, A. G. Electron beam heating effects during environmental scanning electron microscopy imaging of water condensation on superhydrophobic surfaces. Appl Phys Lett 98 (2011), which is incorporated by reference in its entirety. A 500 μm lower aperture was used in series with a 100 μm variable pressure upper aperture to obtain greater detail. The sample temperature was initially set to 4±1.5° C. and was allowed to equilibrate for 5 minutes. The surface temperature was subsequently decreased to 3±1.5° C., resulting in nucleation of water droplets on the sample surface. Accordingly, the supersaturation, S, during the imaging process was in the range of 1<S<1.29. Images and recordings were obtained at an inclination angle of 45° from the horizontal to observe droplet growth. Cu tape was used for mounting the sample to the cold stage to ensure good thermal contact.
Heat Transfer Measurements
The test samples, 6.35 mm diameter tubes with different surface treatments, were placed in an environmental chamber (Kurt J. Lesker) for the heat transfer measurements. A water reservoir, which was connected to the chamber via a vapor valve, was heated to >95° C. to produce steam. The vapor valve was opened to allow steam to flow into the chamber after the chamber was pumped down to the targeted non-condensable pressure (≈30 Pa). Chilled water flowed along the inside of the tube where the inlet temperature and outlet temperature were both measured by thermocouples so that the heat flux could be determined by the temperature rise. The temperature difference, ΔT was determined as the log-mean temperature difference (LMTD) between the vapor and the chilled water. Each data point in
Parameters Affecting the Condensation Heat Transfer Coefficient
Based on the model developed by Miljkovic et al. (Miljkovic, N.; Enright, R.; Wang, E. N., Effect of Droplet Morphology on Growth Dynamics and Heat Transfer during Condensation on Superhydrophobic Nanostructured Surfaces. Acs Nano 2012, 6 (2), 1776-1785, which is incorporated by reference in its entirety), on a dropwise condensation surface, the heat transfer rate through a single growing droplet can be determined as
where Rtot is the total thermal resistance through the droplet, R is the droplet radius, ρw is the liquid water density, hfg is the latent heat of vaporization, Tsat is the vapor saturation temperature, σ is the water surface tension, ΔT is the temperature difference between the saturated vapor and substrate (Tsat−Ts), δHC and h are the hydrophobic coating thickness (˜1 nm) and pillar height, respectively, kHC, kw, and kP are the hydrophobic coating, water, and pillar thermal conductivities, respectively, and h, is the interfacial condensation heat transfer coefficient. See Umur, A.; Griffith, P., Mechanism of Dropwise Condensation. J Heat Transf 1965, 87 (2), 275-&, which is incorporated by reference in its entirety. φ is the solid fraction of the micro/nanostructures. In the special case of a flat surface, φ=1 and h=0.
Droplet size distribution theory was considered to determine the fraction of droplets with a given radius, R, in the droplet heat transfer model. For small droplets, the droplet distribution is determined by
{circumflex over (R)} is the average departure radius, R* is the critical droplet size for nucleation, τ is the droplet sweeping period, and Re is the radius when droplets begin to merge and grow by droplet coalescence afterwards, Re=lC/2 with lc being the coalescence length determined by nucleation density, N
lC=(4N)−2 (16)
See Kim, S.; Kim, K. J., Dropwise Condensation Modeling Suitable for Superhydrophobic Surfaces. J Heat Trans-T Asme 2011, 133 (8), and Rose, J. W., On the mechanism of dropwise condensation. International Journal of Heat and Mass Transfer 1967, 10 (6), 755-762, each of which is incorporated by reference in its entirety.
For large droplets growing mainly due to coalescence, the droplet distribution can be determined as
The total surface condensation heat flux, q″, can be obtained by incorporating the individual droplet heat transfer rate (Eqn. 8) with the droplet size distributions (Eqns. 9 and 17)
q″=∫R″R
The total condensation heat transfer coefficient is determined as
hC=q″/ΔT (19)
Therefore, the sensitivity of hc on the departure radius, advancing contact angle and nucleation density can be obtained as shown in
Estimation of Local Contact Angle on High-Surface-Energy Sites with the Existence of Oil
As shown in
where σvs is the interfacial energy between the surface and vapor, σws is the interfacial energy between water and the surface, and σwv is the interfacial energy between water and vapor, which is 72 mJ/m2.
Similarly, with the introduction of oil which surrounds the water droplet on a surface, the local contact angle, θws(o), can be determined as
where σos is the interfacial energy between the surface and oil, σws is the interfacial energy between water and the surface, and σwo is the interfacial energy between water and oil, which is 49 mJ/m2. See Anand, S.; Paxson, A. T.; Dhiman, R.; Smith, J. D.; Varanasi, K. K., Enhanced Condensation on Lubricant Impregnated Nanotextured Surfaces. Acs Nano 2012, which is incorporated by reference in its entirety.
Since σos is experimentally difficult to obtain for the system, bounds for the local contact angle for the water-oil-substrate system, θws(o), are provided as follows. The contact angle of oil on the high-surface-energy domain is considered using
where σov is the interfacial energy between vapor and oil, which is 17 mJ/m2. Since the oil wets the TFTS-coated surface, which means θos(v)<90°. Therefore, it can be determined that 0<σvs−σos<17 mJ/m2.
As a result, the local contact angle of the water droplet on a surface surrounded by oil can be bounded as
Derivation of Nucleation Rate as a Function of Contact Angle and Interfacial Energy
The nucleation rate, J, can be determined by classical nucleation theory (CNT) as
J=zf*exp(−G*) (24)
See, Kashchiev., D., Nucleation: Basic Theory with Applications. 1 ed.; Oxford: Butterworth-Heinemann.: 2000, which is incorporated by reference in its entirety. In Eqn. (20), z is the Zeldovich factor and G* is the dimensionless energy barrier, given by
z=(kT ln S)2/8πνo√{square root over (kTψ(θ)γ3)} (25)
G*=16πω(θ)νo2γ3/3(kT)2(ln S)2 (26)
where S is the supersaturation and ψ(θ) is the activity that accounts for the effect of contact angle. f* is the frequency of monomer attachment to the critical droplet nucleus dependent on the nature of the nucleus growth. The main modes of growth during heterogeneous nucleation are limited via direct impingement of monomers to the nucleus or surface diffusion. See Pound, G. M.; Simnad, M. T.; Yang, L., Heterogeneous Nucleation of Crystals from Vapor. The Journal of Chemical Physics 1954, 22 (7), 1215-1219 and Sigsbee, R. A., Adatom Capture and Growth Rates of Nuclei. J Appl Phys 1971, 42 (10), 3904-3915, each of which is incorporated by reference in its entirety.
The frequency of monomer attachment due to direct vapor impingement is given by
f*I=γn[(1−cos(θw))/2ψ2/3(θ)](36πνo2)1/3 ln2/3 (27)
where γn is the sticking coefficient (0<γn<1), I is the classical Hertz-Knudsen impingement rate (I=P/√{square root over (2πmokT)})), n is the number of molecules in the nucleated cluster, and νo is the volume of an individual water molecule (νo=3×10−29 m3). To determine an upper bound on the nucleation rate, a sticking coefficient of one was assumed (γn=1).
The frequency of monomer attachment due to surface diffusion is given by
f*sd=γnc*λs2I (28)
where c* is the capture number due to surface diffusion (1<c*<5), and λs is the mean surface diffusion distance of an adsorbed monomer on the substrate. The capture number c* is size independent and approximately equal to 1.9 for heterogeneous condensation of water vapor. See Pocker, D. J.; Hruska, S. J., Detailed Calculations of the Number of Distinct Sites Visited in Random Walk on Several Two-Dimensional Substrate Lattices. Journal of Vacuum Science and Technology 1971, 8 (6), 700-707, which is incorporated by reference in its entirety. The mean surface diffusion distance is dependent on the wettability of the substrate and is given by λs=√{square root over (Dsdτd)} where Dsd is the surface diffusion coefficient (Dsd=ds2νs exp[−Esd/kT]), τd is the desorption time (τd=(1/νs)exp[−Edes/kT]), νs is the adsorbed molecule vibration frequency determined using the Debye approximation (νs=VDα/2), ds is the length of a molecular jump along the substrate surface approximated by the lattice constant of the substrate (ds=5.4 Å) (see J. P. Hirth; G. M. Pound, Condensation and evaporation—nucleation and growth kinetics England: Pergamon Press.: 1963, which is incorporated by reference in its entirety) and VD is the speed of sound in the substrate (VD=8433 m/s). The desorption and surface diffusion energies are given by Edes=E1+σsvao and Esd=0.5Edes (See Thiel, P. A.; Madey, T. E., The interaction of water with solid surfaces: Fundamental aspects. Surface Science Reports 1987, 7 (6-8), 211-385, which is incorporated by reference in its entirety), respectively, where EI is the binding energy of an n=1 sized cluster, σsv is the solid-vapor interfacial energy and ao is the water molecule surface area (ao=4.67×10−19 m2). The calculated energies of desorption show excellent agreement with that of the experiments and molecular dynamics simulations (Edes,SiO2=0.9 eV). See Israelachvili, J. N., Intermolecular and surface forces. 2nd ed.; Academic Press: Amsterdam, 1991, and Ma, Y.; Foster, A. S.; Nieminen, R. M., Reactions and clustering of water with silica surface. The Journal of Chemical Physics 2005, 122 (14), 144709-9, each of which is incorporated by reference in its entirety.
By adding the nucleation rate from the two mechanisms together, the nucleation rate, J, can be determined as a function of the contact angle and interfacial energy of the condensate at given supersaturations, as shown in
Control Experiments on Homogeneous Hydrophobic Surfaces
Dimethyldichlorosilane (DMCS), which is a homogeneous hydrophobic coating, was used in the studies for control experiments. DMCS can be deposited on silicon surfaces using the vapor deposition process as described in the Methods section.
An atomic force microscope (AFM) image of a flat silicon surface coated by DMCS is shown in
Condensation experiments were performed on micropillar arrays coated by DMCS with and without oil-infusion using the same experimental setup for the condensation experiment on TFTS-coated micropillar arrays, as described in the Methods section. The results are summarized and compared to the TFTS-coated surfaces in
Heat Transfer Measurement Apparatus and Experimental Procedure
A custom environmental chamber was built to test the heat transfer performance of each sample for the study. The vacuum chamber was made of stainless steel with two viewing windows. Resistive heater lines were wrapped around the exterior of the chamber walls to prevent condensation at the inside walls, and the chamber was wrapped with insulation on the exterior walls. Two insulated stainless steel water flow lines (Swagelok) were fed into the chamber via a KF flange port to supply cooling water to the chamber from a large capacity chiller (System III, Neslab). A flow meter (7 LPM MAX, Hedland) having an accuracy off ±2% was integrated along the water inflow line.
A secondary stainless steel tube line was fed into the chamber via a KF adapter port that served as the flow line for the incoming water vapor supplied from a heated steel water reservoir. The vapor line was wrapped with a rope heater (60 W, Omega) and controlled by a power supply (Agilent). The vapor reservoir was wrapped with another independently-controlled rope heater (120 W, Omega) and insulated to limit heat losses to the environment. The access tubes were welded to the vapor reservoir, each with independently-controlled valves. The first valve (Diaphragm Type, Swagelok), connecting the bottom of the reservoir to the ambient, was used to fill the reservoir with water. The second valve (BK-60, Swagelok), connecting the top of the reservoir to the inside of the chamber, was used to provide a path for vapor inflow. K-type thermocouples were located along the length of the water vapor reservoir to monitor temperature. To obtain the temperatures within the chamber, K-type thermocouple bundles were connected through the chamber apertures via a thermocouple feed through (Kurt J. Lesker). A pressure transducer (925 Micro Pirani, MKS) was attached to monitor pressure within the chamber. The thermocouple bundles and the pressure transducer were both connected to an analog input source (RAQ DAQ, National Instruments), which was interfaced to a computer to record and store data. A second bellows valve (Kurt J. Lesker) on the chamber was connected to a vacuum pump to bring the chamber down to vacuum conditions prior to vapor filling. A liquid nitrogen cold trap was placed between the chamber and vacuum pump which served to remove any moisture from the pump-down process.
To run the test samples inside the chamber, the stainless steel bellows tube lines (¼″, Swagelok) were connected to the external water flow lines. T-connection adapters (Swagelok) with bore through Ultra-Torr fittings (Swagelok) were used to adapt K-type thermocouple probes (Omega) at the water inlet and outlet. Prior to experimentation, the thermocouple probes were calibrated using a high precision temperature controlled bath (Lauda Brinkman) to an accuracy of ±0.2 K. The test samples, 6.35 mm diameter tubes with different surface treatments, were connected via a Swagelok compression fitting onto the T-connection. Chilled water flows through the inlet bellows tube, along the inside of the tube sample and through the outlet. Two supports were used to hold the sample and the entire configuration in place. Two separate pieces of insulation were embedded with K-type thermocouple leads and used for wet bulb temperature measurements during experimental runs. A third thermocouple was placed beside the sample to measure the reference temperature inside the chamber. As the experiment progressed, the wet-bulb insulating wick collected water from the bottom of the chamber to the embedded thermocouple. The temperature measured by this thermocouple was compared to the reference temperature calculated from the saturation pressure. This allowed for a high accuracy secondary measurement of saturation conditions inside the chamber. F Figure F3 shows the schematic of the test setup for the heat transfer performance measurement.
For each experimental trial, a set of strict procedures were used to ensure consistency throughout the experiments. The water vapor reservoir was filled with approximately 3.5 liters of DI water (99% full) using a syringe through the vapor release valve. After opening the vapor inflow valve and closing the vapor release valve, the rope heater around the water vapor reservoir was turned on and the heater controller set to maximum output. Then the rope heater connected to the vapor inflow valve was turned on. The temperature of the water reservoir was monitored with the installed thermocouples. Once boiling was achieved and all thermocouples on the reservoir reached >95° C. for at least 10 minutes, the vapor inflow valve was closed.
The next step was to begin the vacuum pump-down procedure. Valves connecting the chamber with the ambient, and valves connecting the chamber and the vacuum pump were both closed while the valve connected to the liquid nitrogen cold trap was opened. The vacuum pump was then turned on, initiating the pump-down process where the pressure inside the chamber was carefully monitored. This process took ≈30 minutes in order to achieve the target non-condensable gases pressure (≈30 Pa).
After pumping down, the vapor inflow valve was opened to allow steam flow into the chamber and condensation occurred on the surface of the tube. The heat flux was determined by the rise in the temperature of the chilled water from the inlet to the outlet. The temperature difference, ΔT was determined as the log-mean temperature difference (LMTD) between the vapor and the chilled water. See Mills, A. F., Heat and Mass Transfer. 2 ed.; Prentice-Hall: 1999, which is incorporated by reference in its entirety. Each data point in
Droplet Shedding Radius
The experimentally-determined average droplet shedding radii ({circumflex over (R)}) for a typical dropwise hydrophobic surface (
Examples of Condensation Behavior
In case of TFTS-coated silicon micropillar array where the pillar diameters are 5 μm, periods are 15 μm, and the supersaturation in the experiment is S=1.6, almost no nucleation was observed except on sparse defects in the TFTS coating where hydrophilic silicon oxide substrate was exposed. In case of immersion condensation behavior on oil-infused TFTS-coated silicon micropillar array, where the pillar diameters are 5 μm, periods are 15 μm, and the supersaturation in the experiment is S=1.6, Nucleation occurred on every tip of the pillars, which yields over an order of magnitude higher nucleation density compared to TFTS-coated silicon micropillar array.
When a regular hydrophobic copper tube is horizontally placed with chilled water flowing inside with flow rate of 5 L/min and the vapor pressure in the experiment is ≈2.4 kPa, droplets grow and coalesce before removed by gravity at diameters around 2 mm. When an oil-infused TFTS-coated copper oxide tube is horizontally placed with chilled water flowing inside with flow rate of 5 L/min, and the vapor pressure in the experiment is ≈2.4 kPa, higher droplet density was observed compared to the regular hydrophobic copper tube while the departure diameter is reduced to approximately 0.98±0.13 mm.
In summary, over an order of magnitude increase in the nucleation density on hydrophobic silicon pillar arrays coated with a long-chain silane molecule was observed when hydrophobic oil was introduced on the surface. AFM imaging revealed the existence of locally hydrophilic micelles despite the overall hydrophobicity of the silane self-assembled coating (SAC). The increased nucleation density is explained in the context of classical nucleation theory as the combined effect of the hydrophilic micelles and the reduction in interfacial energy between water and oil. Control experiments on silicon pillar arrays with hydrophobic coatings without micelles and hydrophilic coatings were performed to support these findings. Such phenomena could potentially be used to create surfaces for enhanced condensation heat transfer for a variety of thermal and energy systems.
Other embodiments are within the scope of the following claims.
This application is a continuation of claims U.S. application Ser. No. 13/908,439, filed Jun. 3, 2013, which claims priority to U.S. Provisional Application Ser. No. 61/654,945, filed Jun. 3, 2012, each of which is incorporated by reference in its entirety.
This invention was made with government support under Contract No. N00014-09-1-1000 awarded by the U.S. Navy. The government has certain rights in this invention.
Number | Name | Date | Kind |
---|---|---|---|
20120052241 | King et al. | Mar 2012 | A1 |
20130244001 | Wang et al. | Sep 2013 | A1 |
20140011013 | Jin et al. | Jan 2014 | A1 |
Number | Date | Country | |
---|---|---|---|
20160097606 A1 | Apr 2016 | US |
Number | Date | Country | |
---|---|---|---|
61654945 | Jun 2012 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 13908439 | Jun 2013 | US |
Child | 14519931 | US |