Research Papers: Micro/Nanoscale Heat Transfer

Simulation of Single Bubble Evaporation in a Microchannel in Zero Gravity With Thermocapillary Effect

[+] Author and Article Information
Wei Li

Fellow ASME
Department of Energy Engineering,
Zhejiang University,
Hangzhou 310027, China
e-mail: weili96@zju.edu.cn

Yang Luo

Department of Energy Engineering,
Zhejiang University,
Hangzhou 310027, China;
Department of Energy Engineering,
Collaborative Innovation Center of Advanced
Zhejiang University,
Hangzhou 310027, China

Jingzhi Zhang

School of Energy and Power Engineering,
Shandong University,
Jinan 250061, China

W. J. Minkowycz

Department of Mechanical and
Industrial Engineering,
University of Illinois at Chicago,
Chicago, IL 60607

Contributed by the Heat Transfer Division of ASME for publication in the JOURNAL OF HEAT TRANSFER. Manuscript received November 3, 2017; final manuscript received April 27, 2018; published online July 23, 2018. Assoc. Editor: Guihua Tang.

J. Heat Transfer 140(11), 112403 (Jul 23, 2018) (9 pages) Paper No: HT-17-1654; doi: 10.1115/1.4040147 History: Received November 03, 2017; Revised April 27, 2018

This paper presents fundamental research on the hydrodynamics and heat transfer surrounding a single elongated bubble during flow boiling in a circular microchannel. A continuum surface force (CSF) model based on the volume of fluid (VOF) method is combined with the thermocapillary force to explore the effects of thermocapillarity for flow boiling in microchannels. To validate the self-defined codes, a two-phase thermocapillary-driven flow and a Taylor bubble growing in a capillary tube are studied. Results of both test cases show good convergence and agreement with data from the earlier literature. The bubble motion and the local heat transfer coefficient (HTC) on the heated wall with respect to time are discussed. It is found that for large Marangoni number (case 3), variation of surface tension has affected the bubble shape and temperature profile. The thermocapillary effect induces convection in a thin liquid film region, which augments the HTCs at specified positions. The numerical investigation also shows that the average HTC increased by 6.7% in case 3 when compared with case 1. Thus, it is very important to study further the effects of themocapillarity and the Marangoni effect on bubble growth in microchannels.

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


Kandlikar, S. G. , and Grande, W. J. , 2003, “ Evolution of Microchannel Flow Passages-Thermohydraulic Performance and Fabrication Technology,” Heat Transfer Eng., 24(1), pp. 3–17. [CrossRef]
Li, W. , and Wu, Z. , 2010, “ A General Criterion for Evaporative Heat Transfer in Micro/Mini-Channels,” Int. J. Heat Mass Transfer, 53(9–10), pp. 1967–1976. [CrossRef]
Kandlikar, S. G. , 2002, “ Fundamental Issues Related to Flow Boiling in Minichannels and Microchannels,” Exp. Therm. Fluid Sci., 26(2–4), pp. 389–407. [CrossRef]
Thome, J. R. , 2004, “ Boiling in Microchannels: A Review of Experiment and Theory,” Int. J. Heat Fluid Flow, 25(2), pp. 128–139. [CrossRef]
Garimella, S. V. , and Sobhan, C. B. , 2003, “ Transport in Microchannels—A Critical Review,” Annu. Rev. Heat Transfer, 13(13), pp. 1–50. [CrossRef]
Leong, K. C. , Ho, J. Y. , and Wong, K. K. , 2017, “ A Critical Review of Pool and Flow Boiling Heat Transfer of Dielectric Fluids on Enhanced Surfaces,” Appl. Therm. Eng., 112, pp. 999–1019. [CrossRef]
Guo, Z. , Fletcher, D. F. , and Haynes, B. S. , 2014, “ A Review of Computational Modelling of Flow Boiling in Microchannels,” J. Comput. Multiphase Flow, 6(2), pp. 79–110. [CrossRef]
Jacobi, A. M. , and Thome, J. R. , 2002, “ Heat Transfer Model for Evaporation of Elongated Bubble Flows in Microchannels,” ASME J. Heat Transfer, 124(6), pp. 1131–1131. [CrossRef]
Mukherjee, A. , and Kandlikar, S. G. , 2005, “ Numerical Simulation of Growth of a Vapor Bubble During Flow Boiling of Water in a Microchannel,” Microfluid. Nanofluid, 1(2), pp. 137–145. [CrossRef]
Magnini, M. , Pulvirenti, B. , and Thome, J. R. , 2013, “ Numerical Investigation of the Influence of Leading and Sequential Bubbles on Slug Flow Boiling Within a Microchannel,” Int. J. Therm. Sci., 71(9), pp. 36–52. [CrossRef]
Ling, K. , Son, G. , Sun, D.-L. , and Tao, W.-Q. , 2015, “ Three Dimensional Numerical Simulation on Bubble Growth and Merger in Microchannel Boiling Flow,” Int. J. Therm. Sci., 98, pp. 135–147. [CrossRef]
Magnini, M. , and Thome, J. R. , 2017, “ An Updated Three-Zone Heat Transfer Model for Slug Flow Boiling in Microchannels,” Int. J. Multiphase Flow, 91, pp. 296–314. [CrossRef]
Thome, J. R. , Dupont, V. , and Jacobi, A. M. , 2004, “ Heat Transfer Model for Evaporation in Microchannels—Part I: Presentation of the Model,” Int. J. Heat Mass Transfer, 47(14–16), pp. 3375–3385. [CrossRef]
Karbalaei, A. , Kumar, R. , and Cho, H. J. , 2016, “ Thermocapillarity in Microfluidics—A Review,” Micromachines, 7(1), p. 13.
Zhuan, R. , and Wang, W. , 2010, “ Simulation on Nucleate Boiling in Micro-Channel,” Int. J. Heat Mass Transfer, 53(1–3), pp. 502–512. [CrossRef]
Wang, Z. , Li, S. , Chen, R. , Zhu, X. , and Liao, Q. , 2017, “ Simulation on the Dynamic Flow and Heat and Mass Transfer of a Liquid Column Induced by the IR Laser Photothermal Effect Actuated Evaporation in a Microchannel,” Int. J. Heat Mass Transfer, 113, pp. 975–983. [CrossRef]
Zhang, J. , Li, W. , and Minkowycz, W. J. , 2017, “ Numerical Simulation of R410A Condensation in Horizontal Microfin Tubes,” Numer. Heat Transfer, Part B, 71(4), pp. 361–376. [CrossRef]
Zhang, J. , Li, W. , and Sherif, S. A. , 2016, “ A Numerical Study of Condensation Heat Transfer and Pressure Drop in Horizontal Round and Flattened Minichannels,” Int. J. Therm. Sci., 106, pp. 80–93. [CrossRef]
Zhang, J. , and Li, W. , 2016, “ Numerical Study on Heat Transfer and Pressure Drop Characteristics of R410A Condensation in Horizontal Circular Mini/Micro-Tubes,” Can. J. Chem. Eng., 94(9), pp. 1809–1819. [CrossRef]
Zhang, J. , Fletcher, D. F. , and Li, W. , 2016, “ Heat Transfer and Pressure Drop Characteristics of Gas–Liquid Taylor Flow in Mini Ducts of Square and Rectangular Cross-Sections,” Int. J. Heat Mass Transfer, 103, pp. 45–56. [CrossRef]
Zhang, J. , Li, W. , Sherif, S. A. , Fletcher, D. F. , Li, W. , and Minkowycz, W. J. , 2016, “ Investigation of Hydrodynamic and Heat Transfer Characteristics of Gas-Liquid Taylor Flow in Vertical Capillaries,” Int. Commun. Heat Mass, 74, pp. 1–10. [CrossRef]
Sasmal, G. P. , and Hochstein, J. I. , 1994, “ Marangoni Convection With a Curved and Deforming Free Surface in a Cavity,” ASME J. Fluids Eng., 116(3), pp. 577–582. [CrossRef]
Rattner, A. S. , and Garimella, S. , 2014, “ Simple Mechanistically Consistent Formulation for Volume-Of-Fluid Based Computations of Condensing Flows,” ASME J. Heat Transfer, 136(7), p. 071501.
Hirt, C. W. , and Nichols, B. D. , 1981, “ Volume of Fluid (VOF) Method for the Dynamics of Free Boundaries,” J. Comput. Phys., 39(1), pp. 201–225. [CrossRef]
Brackbill, J. U. , Kothe, D. B. , and Zemach, C. , 1992, “ A Continuum Method for Modeling Surface Tension,” J. Comput. Phys., 100(2), pp. 335–354. [CrossRef]
Nabil, M. , and Rattner, A. S. , 2016, “ InterThermalPhaseChangeFoam—A Framework for Two-Phase Flow Simulations With Thermally Driven Phase Change,” SoftwareX, 5, pp. 216–226. [CrossRef]
Issa, R. I. , Ahmadi-Befrui, B. , Beshay, K. R. , and Gosman, A. D. , 1991, “ Solution of the Implicitly Discretised Reacting Flow Equations by Operator-Splitting,” J. Comput. Phys., 93(2), pp. 388–410. [CrossRef]
Weller, H. G. , Tabor, G. , Jasak, H. , and Fureby, C. , 1998, “ A Tensorial Approach to Computational Continuum Mechanics Using Object-Oriented Techniques,” Comput. Phys., 12(6), pp. 620–631. [CrossRef]
Yang, Z. , Peng, X. F. , and Ye, P. , 2008, “ Numerical and Experimental Investigation of Two Phase Flow During Boiling in a Coiled Tube,” Int. J. Heat Mass Transfer, 51(5–6), pp. 1003–1016. [CrossRef]
Sen, A. K. , and Davis, S. H. , 1982, “ Steady Thermocapillary Flows in Two-Dimensional Slots,” J. Fluid Mech., 121(1), pp. 163–186. [CrossRef]
Saldi, Z. S. , 2011, “ Marangoni Driven Free Surface Flows in Liquid Weld Pools,” Ph.D. dissertation, Delft University of Technology, Delft, The Netherlands. https://repository.tudelft.nl/islandora/object/uuid:8401374b-9e9c-4d25-86b7-fc445ec73d27?collection=research
Francois, M. M. , Sicilian, J. M. , and Kothe, D. B. , 2006, “ Modeling of Thermocapillary Forces Within a Volume Tracking Algorithm,” Modeling of Casting, Welding and Advanced Solidification Processes–XI, Opio, France, May 28–June 2, pp. 935–942. https://www.researchgate.net/publication/276919120_Modeling_of_Thermocapillary_Forces_within_a_Volume_Tracking_Algorithm
Son, G. , 2014, “ A Sharp-Interface Level-Set Method for Analysis of Marangoni Effect on Microdroplet Evaporation,” Int. Commun. Heat Mass, 58, pp. 156–165. [CrossRef]
Yamamoto, T. , Okano, Y. , and Dost, S. , 2017, “ Validation of the S-CLSVOF Method With the Density-Scaled Balanced Continuum Surface Force Model in Multiphase Systems Coupled With Thermocapillary Flows,” Int. J. Numer. Methods Fluids, 83(3), pp. 223–244. [CrossRef]
Magnini, M. , Pulvirenti, B. , and Thome, J. R. , 2013, “ Numerical Investigation of Hydrodynamics and Heat Transfer of Elongated Bubbles During Flow Boiling in a Microchannel,” Int. J. Heat Mass Transfer, 59(1), pp. 451–471. [CrossRef]
Pan, Z. , Weibel, J. A. , and Garimella, S. V. , 2016, “ A Saturated-Interface-Volume Phase Change Model for Simulating Flow Boiling,” Int. J. Heat Mass Transfer, 93, pp. 945–956. [CrossRef]
Li, S. , Chen, R. , Zhu, X. , and Liao, Q. , 2016, “ Numerical Investigation of the Marangoni Convection During the Liquid Column Evaporation in Microchannels Caused by IR Laser Heating,” Int. J. Heat Mass Transfer, 101, pp. 970–980. [CrossRef]
Han, Y. , and Shikazono, N. , 2009, “ Measurement of the Liquid Film Thickness in Micro Tube Slug Flow,” Int. J. Heat Fluid Flow, 30(5), pp. 842–853. [CrossRef]
Luo, Y. , Zhang, J. , Li, W. , Sokolova, E. , Li, Y. , and Minkowycz, W. J. , 2017, “ Numerical Investigation of the Bubble Growth in Horizontal Rectangular Microchannels,” Numer. Heat Transfer, Part A, 71(12), pp. 1175–1188. [CrossRef]
Han, Y. , Shikazono, N. , and Kasagi, N. , 2012, “ The Effect of Liquid Film Evaporation on Flow Boiling Heat Transfer in a Micro Tube,” Int. J. Heat Mass Transfer, 55(4), pp. 547–555. [CrossRef]


Grahic Jump Location
Fig. 1

Surface tension coefficient versus temperature for R113; the properties of R113 are obtained from the software refprop 9.11

Grahic Jump Location
Fig. 2

Algorithm to solve phase change problem with thermocapillary effect

Grahic Jump Location
Fig. 3

Illustration of tangential and normal components of surface tension force

Grahic Jump Location
Fig. 4

(a) Problem sketch of free surface thermocapillary-driven flows in a cavity and (b) temperature field (solid line) and streamline (dashed line) obtained at the steady-state (d = 0.01). The dimensionless interval between temperature contours is 0.1. The dimensionless interval between streamline contours is 1.6 × 10−3 in both phases.

Grahic Jump Location
Fig. 5

Schematic illustration of: (a) an elongated vapor bubble flowing through the microchannel in saturated liquid with an upstream adiabatic section and downstream heating section and (b) a spherical vapor bubble growing in a microchannel of diameter D = 0.2 mm

Grahic Jump Location
Fig. 6

(a) Bubble nose instantaneous location versus time after entering the heating section and (b) bubble shape and temperature profile of the case with thermocapillary effect at tm; Comparison of HTCs on the heated wall at tm. The results of our simulation are compared with those in the literature [35,36].

Grahic Jump Location
Fig. 7

Bubble interface evolution and temperature distribution of case 1 (top) and case 3 (bottom)

Grahic Jump Location
Fig. 8

(a) Comparison of bubble nose and tail location against time for different Marangoni numbers and (b) dimensionless liquid film thickness (δ/D) and superheat temperature of the heated wall along the axial direction at 6 ms

Grahic Jump Location
Fig. 9

Time evolution of HTC. tin denotes the time instance when the bubble nose enters the position of x/R = 5, 10, and 15.

Grahic Jump Location
Fig. 10

Time evolution of average HTC on the heated wall



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