Research Papers

Isoflux Nusselt Number and Slip Length Formulae for Superhydrophobic Microchannels

[+] Author and Article Information
Ryan Enright

Stokes Institute,
University of Limerick,
Limerick, Ireland
e-mail: ryan.enright@alcatel-lucent.com

Marc Hodes

Department of Mechanical Engineering,
Tufts University,
Medford, MA 02155

Todd Salamon

Bell Labs,
600 Mountain Avenue,
Murray Hill, NJ 07974

Yuri Muzychka

Faculty of Engineering and Applied Science,
Memorial University of Newfoundland,
St. John's, NL A1B 3X5

1Present address: Thermal Management Research Group, Efficient Energy Transfer (ηET) Dept., Bell Labs Ireland, Alcatel-Lucent Ireland, Blanchardstown Business & Technology Park, Blanchardstown, Snugborough Road, Dublin 15, Ireland.

Contributed by the Heat Transfer Division of ASME for publication in the Journal of Heat Transfer. Manuscript received September 10, 2012; final manuscript received May 8, 2013; published online October 17, 2013. Assoc. Editor: Jose L. Lage.

J. Heat Transfer 136(1), 012402 (Oct 17, 2013) (9 pages) Paper No: HT-12-1493; doi: 10.1115/1.4024837 History: Received September 10, 2012; Revised May 08, 2013

We analytically and numerically consider the hydrodynamic and thermal transport behavior of fully developed laminar flow through a superhydrophobic (SH) parallel-plate channel. Hydrodynamic slip length, thermal slip length and heat flux are prescribed at each surface. We first develop a general expression for the Nusselt number valid for asymmetric velocity profiles. Next, we demonstrate that, in the limit of Stokes flow near the surface and an adiabatic and shear-free liquid–gas interface, both thermal and hydrodynamic slip lengths can be found by redefining existing solutions for conduction spreading resistances. Expressions for the thermal slip length for pillar and ridge surface topographies are determined. Comparison of fundamental half-space solutions for the Laplace and Stokes equations facilitate the development of expressions for hydrodynamic slip length over pillar-structured surfaces based on existing solutions for the conduction spreading resistance from an isothermal source. Numerical validation is performed and an analysis of the idealized thermal transport behavior suggests conditions under which superhydrophobic microchannels may enhance heat transfer.

Copyright © 2014 by ASME
Your Session has timed out. Please sign back in to continue.


Quéré, D., 2008, “Wetting and Roughness,” Annu. Rev. Mater. Res., 38, pp. 71–99. [CrossRef]
Ou, J., Perot, B., and Rothstein, J. P., 2004, “Laminar Drag Reduction in Microchannels Using Ultrahydrophobic Surfaces,” Phys. Fluids, 16(12), p. 4635. [CrossRef]
Ou, J., and Rothstein, J. P., 2005, “Direct Velocity Measurements of the Flow Past Drag-Reducing Ultrahydrophobic Surfaces,” Phys. Fluids, 17(10), p. 3606. [CrossRef]
Woolford, B., Jeffs, K., Maynes, D., and Webb, B. W., 2005, “Laminar Fully-Developed Flow in a Microchannel With Patterned Ultrahydrophobic Walls,” Proceedings of HT2005, ASME.
Maynes, D., Jeffs, K., Woolford, B., and Webb, B. W., 2007, “Laminar Flow in a Microchannel With Hydrophobic Surface Patterned Microribs Oriented Parallel to the Flow Direction,” Phys. Fluids, 19, p. 093603. [CrossRef]
Joseph, P., 2005, “Étude expérimentale du glissement sur surfaces lisses et texturées,” Ph.D. thesis, Université Paris 6, Paris, France.
Joseph, P., Cottin-Bizonne, C., Benoit, J. M., Ybert, C., Journet, C., Tabeling, P., and Bocquet, L., 2006, “Slippage of Water Past Superhydrophobic Carbon Nanotube Forests in Microchannels,” Phys. Rev. Lett., 97(15), p. 6104. [CrossRef]
Enright, R., Krupenkin, T. N., Kolodner, P., Hodes, M., Salamon, T., and Dalton, T., 2008, “Effects of Interfacial Position on Drag Reduction in a Superhydrophobic Microchannel,” Proceedings of ICNMM2008, ASME.
Lee, C., Choi, C. H., and Kim, C. J., 2008, “Structured Surfaces for a Giant Liquid Slip,” Phys. Rev. Lett., 101(6), p. 4501. [CrossRef]
Maynes, D., Webb, B. W., and Davies, J., 2008, “Thermal Transport in a Microchannel Exhibiting Ultrahydrophobic Microribs Maintained at Constant Temperature,” ASME J. Heat Transfer, 130(2), p. 2402. [CrossRef]
Maynes, D., Webb, B. W., Crockett, J., and Solovjov, V., 2013, “Analysis of Laminar Slip-Flow Thermal Transport in Microchannels With Transverse Ribs and Cavity Structured Superhydrophobic Walls at Constant Heat Flux,” ASME J. Heat Transfer, 135(2), p. 021701. [CrossRef]
Baier, T., Steffes, C., and Hardt, S., 2010, “Thermocapillary Flow on Superhydrophobic Surfaces,” Phys. Rev. E, 82, p. 037301. [CrossRef]
Rosengarten, G., Stanley, C., and Kwok, F., 2008, “Superinsulating Heat Transfer Surfaces for Microfluidic Channels,” Int. J. Transp. Phenom., 10, pp. 293–306.
Ybert, C., Barentin, C., Cottin-Bizonne, C., Joseph, P., and Bocquet, L., 2007, “Achieving Large Slip With Superhydrophobic Surfaces: Scaling Laws for Generic Geometries,” Phys. Fluids, 19(12), p. 3601. [CrossRef]
Lauga, E., and Stone, H. A., 2003, “Effective Slip in Pressure-Driven Stokes Flow,” J. Fluid Mech., 489, pp. 55–77. [CrossRef]
Philip, J. R., 1972, “Flows Satisfying Mixed No-Slip and No-Shear Conditions,” Z. Agnew. Math. Phys., 23, pp. 353–370. [CrossRef]
Philip, J. R., 1972, “Integral Properties of Flows Satisfying Mixed No-Slip and No-Shear Conditions,” Z. Agnew. Math. Phys., 23, pp. 960–968. [CrossRef]
Navier, C. L. M. H., 1823, “Mémoire sur les lois du movement des fluids,” Mémoires de l'Académie Royale des Sciences de l'Institut de France, 6, pp. 389–440.
Cottin-Bizonne, C., Barentin, C., Charlaix, E., Bocquet, L., and Barrat, J. L., 2004, “Dynamics of Simple Liquids at Heterogenous Surfaces: Molecular-Dynamics Simulations and Hydrodynamic Description,” Eur. Phys. J. E, 15, pp. 427–438. [CrossRef]
Choi, C. H., Ulmanella, U., Kim, J., Ho, C. M., and Kim, C. J., 2006, “Effective Slip and Friction Reduction in Nanograted Superhydrophobic Microchannels,” Phys. Fluids, 18(8), p. 087105. [CrossRef]
Davis, A. M. J., and Lauga, E., 2010, “Hydrodynamic Friction of Fakir-Like Superhydrophobic Surfaces,” J. Fluid Mech., 661, pp. 402–411. [CrossRef]
Davis, A. M. J., and Lauga, E., 2009, “Geometric Transition in Friction for Flow Over a Bubble Mattress,” Phys. Fluids, 21(1), p. 011701. [CrossRef]
Ng, C. O., and Wang, C. Y., 2010, “Apparent Slip Arising From Stokes Shear Flow Over a Bidimensional Patterned Surface,” Microfluid. Nanofluid., 8, pp. 361–371. [CrossRef]
Ng, C. O., and Wang, C. Y., 2009, “Stokes Shear Flow Over a Grating: Implications for Superhydrophobic Slip,” Phys. Fluids, 21(1), p. 3602. [CrossRef]
Woolford, B., Maynes, D., and Webb, B. W., 2008, “Liquid Flow Through Microchannels With Grooved Walls Under Wetting and Superhydrophobic Conditions,” Microfluid. Nanofluid., 7(1), pp. 121–135. [CrossRef]
Davies, J., Maynes, D., Webb, B. W., and Woolford, B., 2006, “Laminar Flow in a Microchannel With Superhydrophobic Walls Exhibiting Transverse Ribs,” Phys. Fluids, 18(8), p. 7110. [CrossRef]
Sparrow, E. M., Lin, S. H., and Lundgren, T. S., 1964, “Flow Development in the Hydrodynamic Entrance Region of Tubes and Ducts,” Phys. Fluids, 7(3), pp. 338–347. [CrossRef]
Sparrow, E. M., and Lin, S. H., 1962, “Laminar Heat Transfer in Tubes Under Slip-Flow Conditions,” ASME J. Heat Transfer, 84(4), pp. 363–369. [CrossRef]
Bocquet, L., and Barrat, J.-L., 2007, “Flow Boundary Conditions From Nano- to Micro-Scales,” Soft Matter, 3(6), pp. 685–693. [CrossRef]
Nield, D. A., 2004, “Forced Convection in a Parallel Plate Channel With Asymmetric Heating,” Int. J. Heat Mass Transfer, 47, pp. 5609–5612. [CrossRef]
Nield, D. A., 2008, “Erratum to “Forced Convection in a Parallel Plate Channel With Asymmetric Heating,”” Int. J. Heat Mass Transfer, 51, p. 2108. [CrossRef]
Kays, W. M., and Crawford, M. E., 1993, Convective Heat and Mass Transfer, 3rd ed., McGraw-Hill, New York.
Inman, R. M., 1964, “Laminar Slip Flow Heat Transfer in a Parallel-Plate Channel or a Round Tube With Uniform Wall Heating,” NASA Technical Note D-2393, Lewis Research Center, Cleveland, OH.
Rohsenow, W. M., Hartnett, J. P., and Ganić, E. N., 1985, Handbook of Heat Transfer Fundamentals, 2nd ed., McGraw-Hill, New York.
Enright, R., Hodes, M., Salamon, T. R., and Muzychka, Y., 2010, “Analysis and Simulation of Heat Transfer in a Superhydrophobic Microchannel,” Proceedings of 14th International Heat Transfer Conference, Vol. 6, August, ASME, pp. 157–168.
Cottin-Bizonne, C., Barentin, C., and Bocquet, L., 2012, “Scaling Laws for Slippage on Superhydrophobic Fractal Surfaces,” Physics of Fluids, 24(1), p. 012001. [CrossRef]
Bejan, A., and Kraus, A. D., eds., 2003, Heat Transfer Handbook, John Wiley & Sons, Hoboken, NJ.
Mikic, B. B., and Rohsenow, W. M., 1966, “Thermal Contact Resistance,” MIT Mechanical Engineering Report No. DSR 74542-41.
Negus, K. J., Yovanovich, M. M., and Beck, J. V., 1989, “On the Nondimensionalization of Constriction Resistance for Semi-Infinate Heat Flux Tubes,” ASME J. Heat Transfer, 111, pp. 804–807. [CrossRef]
Veziroglu, T. N., and Chandra, S., 1969, “Thermal Conductance of Two-Dimensional Constrictions,” Prog. Astronaut. Aeronaut., 21, pp. 591–615.
Carslaw, H. S., and Jaeger, J. C., eds., 1959, Conduction of Heat in Solids, 2nd ed., Oxford University Press, London.
Lamb, H., ed., 1932, Hydrodynamics, 6th ed., Cambridge University Press, Cambridge, UK.
Yovanovich, M. M., 1975, “General Expressions for Constriction Resistances Due To Arbitrary Flux Distributions at Non-Symmetric, Coaxial Contacts,” 13th Aerospace Sciences Meeting, AIAA.
Negus, K. J., and Yovanovich, M. M., 1984, “Constriction Resistance of Circular Flux Tubes With Mixed Boundary Conditions by Linear Superposition of Neumann Solutions,” ASME National Heat Transfer Conference, ASME.
Lobaton, E. J., and Salamon, T., 2007, “Computation of Constant Mean Curvature Surfaces: Application to the Gas-Liquid Interface of a Pressurized Fluid on a Superhydrophobic Surface,” J. Colloid Interface Sci., 314, pp. 184–198. [CrossRef] [PubMed]
Sbragaglia, M., and Prosperetti, A., 2007, “A Note on the Effective Slip Properties for Microchannel Flows With Ultrahydrophobic Surfaces,” Phys. Fluids, 19(04), p. 3603. [CrossRef]
Teo, C. J., and Khoo, B. C., 2009, “Flow Past Superhydrophobic Surfaces Containing Longitudinal Grooves: Effects of Interface Curvature,” Microfluid. Nanofluid., 9(2), pp. 499–511. [CrossRef]
Salamon, T., Lee, W., Krupenkin, T., Hodes, M., Kolodner, P., Enright, R., and Salinger, A., 2005, “Numerical Simulation of Fluid Flow in microchannels With Superhydrophobic Walls,” Proceedings of IMECE 2005, ASME.
Tuckerman, D. B., and Pease, R. F. W., 1981, “High Performance Heat Sinking for VLSI,” IEEE Electron. Device Lett., 2(5), pp. 126–129. [CrossRef]
Hodes, M., Zhang, R., Wilcoxon, R., and Lower, N., 2012, “On the Cooling Potential of Galinstan-Based Minichannel Heat Sinks,” ITherm 2012, IEEE.


Grahic Jump Location
Fig. 1

Geometry of a parallel-plate channel where H is the channel height. Both the channel width and length are considered much larger than H. Each of the surfaces exhibits arbitrary hydrodynamic and thermal slip.

Grahic Jump Location
Fig. 2

Isometric (a) and plan (b) views of the solution domain. The origin of the coordinate system is defined at the center of the pillar tip. The composite interface at y = 0 is composed of a liquid–solid interface, Bs, subjected to no-slip and constant heat flux. The liquid–gas interface, Bi, was specified as shear-free and adiabatic. The domain extended to a symmetric boundary, Bsymy=H/2, at the channel center-line. Flow imposed by a pressure gradient entered the domain from the left (x-direction) through Bperx=-l/2 with a specified mean temperature and exited to the right through Bperx=l/2 in a periodic manner. Bsymz=0 and Bsymz=l/2 represent symmetry boundaries spaced l/2 apart in the z-direction.

Grahic Jump Location
Fig. 3

Dimensionless temperature profile for φs = 0.01 (a) in the channel and (b) a magnified view near the surface

Grahic Jump Location
Fig. 4

The numerical thermal slip length data (•) are compared to the analytical conduction solution, Eq. (29) (solid curve). The predicted behavior of ridges, Eq. (31), is shown by the dashed curve.

Grahic Jump Location
Fig. 5

Comparison of our numerical hydrodynamic slip length data (▲) to the asymptotic expression of Davis and Lauga [21], b/l=(3π/16)φs-1/2-0.421 (dotted–dashed curve), Eq. (39) (dashed curve), and Eq. (41) (solid curve) as a function of (a) φs-1/2 and (b) φs. The inset of (a) shows the ratio of the slip length predicted by Davis and Lauga [21] (dashed curve) and Eq. (37) (solid curve) to Eq. (41). The numerical data of Ng and Wang [23] (open circles/squares corresponding to circular/square pillars) and Ybert et al. [14] (filled squares corresponding to square pillars) are also shown in (b).

Grahic Jump Location
Fig. 6

(a) Relationship between the thermal and hydrodynamic slip lengths. The numerical data (▲) are compared to the solid line which corresponds to the analytical prediction for pillars. The dashed and dashed–dotted curves correspond to the analytical prediction for ridges aligned parallel and transverse to the channel flow, respectively. (b) Nusselt number Nu and (c) the thermohydraulic measure Φ/Φc plotted as a function of b˜ for pillars (b˜t≈1.5b˜) (solid curve), parallel ridges (b˜t≈1.05b˜) (dashed curve) and transverse ridges (b˜t≈2.1b˜) (dotted–dashed curve) in the Stokes flow limit. Also shown is the limit of b˜t=0 (blue dashed curve) and a range of 0≤b˜t≤2.5b* in increments of 0.25b˜t (light dotted curves).




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In