Research Papers

Two-Phase Flow With Phase Change in Porous Channels

[+] Author and Article Information
Gaurav Tomar

Assistant Professor
Department of Mechanical Engineering,
Indian Institute of Science,
Bangalore 560012, India
e-mail: gtom@mecheng.iisc.ernet.in

Contributed by the Heat Transfer Division of ASME for publication in the JOURNAL OF THERMAL SCIENCE AND ENGINEERING APPLICATIONS. Manuscript received July 27, 2014; final manuscript received December 15, 2014; published online January 21, 2015. Assoc. Editor: Ranganathan Kumar.

J. Thermal Sci. Eng. Appl 7(2), 021006 (Jun 01, 2015) (8 pages) Paper No: TSEA-14-1172; doi: 10.1115/1.4029458 History: Received July 27, 2014; Revised December 15, 2014; Online January 21, 2015

Phase change heat transfer in porous media finds applications in various geological flows and modern heat pipes. We present a study to show the effect of phase change on heat transfer in a porous channel. We show that the ratio of Jakob numbers based on wall superheat and inlet fluid subcooling governs the liquid–vapor interface location in the porous channel and below a critical value of the ratio, the liquid penetrates all the way to the extent of the channel in the flow direction. In such cases, the Nusselt number is higher due to the proximity of the liquid–vapor interface to the heat loads. For higher heat loads or lower subcooling of the liquid, the liquid–vapor interface is pushed toward the inlet, and heat transfer occurs through a wider vapor region thus resulting in a lower Nusselt number. This study is relevant in the designing of efficient two-phase heat exchangers such as capillary suction based heat pipes where a prior estimation of the interface location for the maximum heat load is required to ensure that the liquid–vapor interface is always inside the porous block for its operation.

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


Lasance, C. M., 1997, “Advances in High-Performance Cooling for Electronics,” Electron. Cool., 11(4), pp. 22–39.
Incropera, F. P., 1988, “Convection Heat Transfer in Electronic Equipment Cooling,” ASME J. Heat Transfer, 110, pp. 1097–1111. [CrossRef]
Pamula, V., and Chakrabarty, K., 2003, “Cooling of Integrated Circuits Using Droplet-Based Microfluidics,” ACM Great Lakes Symposium on VLSI, pp. 84–87. [CrossRef]
Lee, D., and Vafai, K., 1999, “Comparative Analysis of Jet Impingement and Microchannel Cooling for High Heat Flux Applications,” Int. J. Heat Mass Transfer, 42(9), pp. 1555–1568. [CrossRef]
Mohapatra, S., and Loikitis, D., 2005, “Advances in Liquid Coolant Technologies for Electronics Cooling,” 21st SemiTherm Symposium, Mar. 15–17, pp. 354–360. [CrossRef]
Alkam, M. K., Al-Nimr, M. A., and Hamdan, M. O., 2001, “Enhancing Heat Transfer in Parallel-Plate Channels by Using Porous Inserts,” Int. J. Heat Mass Transfer, 44(5), pp. 931–938. [CrossRef]
North, M. T., and Cho, W.-L., 2003, “High Heat Flux Liquid-Cooled Porous Metal Heat Sink,” ASME Paper No. IPACK2003-35320. [CrossRef]
Ko, K.-H., 2004, “Heat Transfer Enhancement in a Channel With Porous Baffles,” Ph.D. thesis, Texas A&M University, College Station, TX.
Vasiliev, L. L., 2005, “Heat Pipes in Modern Heat Exchangers,” Appl. Therm. Eng., 25(1), pp. 1–19. [CrossRef]
Wang, G., Mishkinis, D., and Nikanpour, D., 2008, “Capillary Heat Loop Technology: Space Applications and Recent Canadian Activities,” Appl. Therm. Eng., 28(4), pp. 284–303. [CrossRef]
Figus, C., Bray, Y. L., Bories, S., and Prat, M., 1999, “Heat and Mass Transfer With Phase Change in a Porous Structure Partially Heated: Continuum Model and Pore Network Simulations,” Int. J. Heat Mass Transfer, 42(14), pp. 2557–2569. [CrossRef]
Demidov, A. S., and Yatsenko, E. S., 1994, “Investigation of Heat and Mass Transfer in the Evaporation Zone of a Heat Pipe Operating by the ‘Inverted Meniscus’ Principle,” Int. J. Heat Mass Transfer, 37(14), pp. 2155–2163. [CrossRef]
Meakin, P., and Xu, Z., 2008, “Disspative Particle Dynamics and Other Particle Methods for Multiphase Fluid Flow in Fractured and Porous Media,” 6th International Conference on CFD in Oil and Gas, Metallurgical and Process Industries, Paper No. CFD08-52.
Liu, M., Meakin, P., and Huang, H., 2006, “Dissipative Particle Dynamics With Attractive and Repulsive Particle–Particle Interactions,” Phys. Fluids, 18(1), p. 017101. [CrossRef]
Carciofi, B. A. M., Prat, M., and Laurindo, J. B., 2011, “Homogeneous Volume of Fluid Method for Simulating Imbibition in Porous Media Saturated by Gas,” Energy Fuels, 25(5), pp. 2267–2273. [CrossRef]
Riaz, A., and Tchelepi, H. A., 2006, “Numerical Simulations of Immiscible Two-Phase Flow in Porous Media,” Phys. Fluids, 18(1), p. 014104. [CrossRef]
Aziz, K., and Settari, A., 1979, Petroleum Reservoir Simulation, Applied Science Publishers, London.
Yang, Z., Peng, X. F., and Ye, P., 2008, “Numerical Experimental Investigation of Two Phase Flow During Boiling in Coiled Tube,” Int. J. Heat Mass Transfer, 51(5), pp. 1003–1016. [CrossRef]
Buckley, S. E., and Leverett, M. C., 1942, “Mechanisms of Fluid Displacements in Sands,” Trans. AIME, 146(1), pp. 107–116. [CrossRef]
Leveque, R. J., 2004, Finite-Volume Method for Hyperbolic Problems, Cambridge University Press, Cambridge, UK. [CrossRef]
Bastian, P., 1999, “Numerical Computation of Multiphase Flows in Porous Media,” Ph.D. thesis, Technischen Fakultät der Christian-Albrechts-Universität Kiel, Heidelberg, Germany.
Bear, J., 1972, Dynamics of Fluids in Porous Media, Dover Publications, Mineola, NY.
Chavent, G., and Jaffre, J., 1978, Mathematical Models and Finite Elements for Reservoir Simulations (Studies in Mathematics and Its Applications, Vol. 17), Elsevier, Philadelphia, PA.
Fang, C., David, M., Rogacs, A., and Goodson, K., 2010, “Volume of Fluid Simulations of Boiling Two-Phase Flow in a Vapor-Venting Microchannel,” Front. Heat Mass Transfer, 1, p. 013002. [CrossRef]
Chavent, G., 1975, “A New Formulation of Diphysic Incompressible Flows in Porous Media,” Applications of Methods of Functional Analysis to Problem in Mechanics, Springer, Berlin, Germany, pp. 258–270. [CrossRef]
Brinkman, H. C., 1947, “Calculation of the Viscous Force Extended by a Flowing Fluid on a Dense Swarm of Particles,” Appl. Sci. Res., A1, pp. 27–34.
Vafai, K., and Tien, C. L., 1981, “Boundary and Inertia Effects on Flow and Heat Transfer in Porous Media,” Int. J. Heat Mass Transfer, 24(2), pp. 195–203. [CrossRef]
Falgout, R., Jones, J., and Yang, U., 2006, “The Design and Implementation of Hypre: A Library of Parallel High Performance Preconditioners,” Numerical Solution of Partial Differential Equations on Parallel Computers (Lecture Notes in Computational Science and Engineering, Vol. 51), A. Bruaset and A. Tveito, eds., Springer, Berlin, pp. 267–294. [CrossRef]
Sazhenkov, S. A., 2008, “Studying the Darcy–Stefan Problem on Phase Transition in a Saturated Porous Soil,” ASME J. Appl. Mech. Tech. Phys., 49(4), pp. 587–597. [CrossRef]


Grahic Jump Location
Fig. 1

Schematic showing a porous channel with liquid being imbibed from the left and vapor exiting from the right end. Heating elements, shown by two rectangular blocks at the top and bottom walls of the channel, cause evaporation of the liquid.

Grahic Jump Location
Fig. 2

Typical curves for the fractional flow fw as a function of saturation of the wetting phase. The dashed line is the tangent on the curve for M = 1/4 which passes through origin. The slope, dxs/dt, is the shock speed for an imbibition flow with M = 1/4. The remaining portion of the curve shows the behavior of the expansion wave.

Grahic Jump Location
Fig. 3

Comparison of the numerical results with the solution obtained from Buckley–Leverett analysis for M = 1/4 for two grid-sizes

Grahic Jump Location
Fig. 4

Saturation profiles for different viscosity ratios obtained using the IMPES method

Grahic Jump Location
Fig. 5

Comparison of the analytical profiles (Eq. (39)) with the numerical results at different times

Grahic Jump Location
Fig. 6

Saturation profiles corresponding to the temperature field at different times shown in Fig. 5

Grahic Jump Location
Fig. 7

Steady-state saturation field for Jaw/Jaf = 0.025. Region on the right from the interface represents vapor whereas the region on the left of the interface marks the saturated liquid. The gray rectangular blocks represent the location of the heating elements at temperature Tw. The length of the heating elements is one-third of the channel length, which are placed at the center of the top and bottom walls.

Grahic Jump Location
Fig. 8

Steady-state temperature field (θ = (T-Tsat)/(Tw-Tsat)) for Jaw/Jaf = 0.025

Grahic Jump Location
Fig. 9

Steady-state saturation field for Jaw/Jaf = 0.00525. Region on the right of the interface represents vapor whereas on the left of the interface marks the imbibed liquid.

Grahic Jump Location
Fig. 10

Steady-state temperature field (θ = (T-Tsat)/(Tw-Tsat)) for Jaw/Jaf = 0.00525

Grahic Jump Location
Fig. 11

Variation in the maximum penetration of the liquid in the porous channel with the ratio Jaw/Jaf

Grahic Jump Location
Fig. 12

Variation in the space averaged Nusselt number at the heating elements with the ratio Jaw/Jaf



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