Numerical Details of Convective Heat Transfer by Micro-Encapsulated PCM Case Study: Annular Slurry Flow

Document Type : Research Paper


Department of Mechanical Engineering, Faculty of Engineering, Razi University, P. Box: 6714414971, Kermanshah, Iran


Due to the high latent heat value, using microencapsulated PCMs increases the heat transfer coefficient in the heat sinks in mini electronic devices, chilled celling, … . In this paper, convective heat transfer by mixed PCM particles in a fluid as slurry, has been studied by the Eulerian-Lagrangian two-phase method. In this method, the fluid phase is studied by the Eulerian and the particle phase is studied using the Lagrangian view. In this paper, the base fluid is water and the particles made of encapsulated micro-size paraffin wax which has covered by a thin layer of Fe3O4. The fluid phase is solved by a control volume method (SIMPLE) and the velocities of the particle phase are solved by the 4th order of the Runge-Kutta method. Due to high Biot number for particles, the lumped temperature assumption for particles is not valid and the transient one dimensional conduction equation has been solved. In this paper details of solving the energy equation inside the particles has been presented. The results include the local and mean Nusselt numbers for different Reynolds numbers including 200, 350 and 500, wide range of the volume fraction from 0-5% for PCM particle with 10 micro-meter diameter, inside the mini annular tube with inner diameter of 1 mm and outer diameter of 3 mm. The results show for  and Re=200, 500, the Nusselt number increases by 10 and 12.5%, while the pressure loss increases by 2 and 5.5% respectively. The maximum performance coefficient is 1.078 and occurs for Re=200 at .


[1]        X. Wang, J. Niu, Performance of cooled-ceiling operating with MPCM slurry, Energy Conversion and Management, 50 (2009), 583-591.
[2]        P.W. Griffiths, P.C. Eames, Performance of chilled ceiling panels using phase change material slurries as the heat transport medium, Applied Thermal Engineering 27 (2007) 1756–1760.
[3]        M. Safdari, R. Ahmadi, S. Sadeghzadeh, Numerical investigation on PCM encapsulation shape used in the passive-active battery thermal management, Energy (2020), doi:
[4]        W. M. Yan, C. J. Ho, Y. T. Tseng, C. Qin, S. Rashidi, Experimental and numerical investigation of the latent heat thermal storage unit with PCM packing at the inner side of a tube, International Journal of Heat and Mass Transfer, 152 (2020), 119480.
[5]        M. Ghalambaz, A. J. Chamkha, D. Wen, Natural convection flow and heat transfer of nano-Encapsulated phase change Materials (NEPCMs) in a cavity, International Journal of Heat and Mass Transfer, 138 (2019), 738-749.
[6]        J. Rostami, A. Abbassi, J. Harting, Heat Transfer by Nanofluids in Wavy Microchannels, Advanced Powder Technology, 29 (4), (2018), 925-933.
[7]        M. Kalteh, A., Abbassi, M., Saffar-avval, J., Harting, Eulerian- Eulerian Two-Phase Numerical Simulation of Nanofluid Laminar Forced Convection in a Microchannel, International Journal of Heat and Fluid Flow, 32 (2011), 107-116.
[8]        M. Mirzaei, M. Saffar-Avval, H. Naderan, Heat Transfer Investigation of Laminar Developing Flow of Nanofluids in a Microchannel Based on Eulerian-Lagrangian Approach, The Canadian Journal of Chemical Engineering, 92 (2014), 1139-1149.
[9]        W. J. Minkowycz, E. M. Sparrow, J. Y. Murthy, 2006, Handbook of Numerical Heat Transfer, second ed., John Wiley & Sons, Hoboken, NJ.
[10]     M. Zeneli, I. Malgarinos, A. Nikolopoulos, N. Nikolopoulos, P. Grammelis, S. Karellas, E. Kakaras, Numerical simulation of a silicon-based latent heat thermal energy storage system operating at ultra-high temperatures, Applied Energy, 242 (2019), 837-853.
[11]     D. Wen, L. Zhang, Y. He, Flow and Migration of Nanoparticle in a Single Channel, Heat and Mass Transfer, 45 (2009), 1061-1067.
[12]     J. Rostami, Convective Heat Transfer by Micro-Encapsulated PCM in a mini-Duct, International Journal of Thermal Sciences, 161 (2021), 106737
[13]     S. V. Patankar and D.B. Spalding, A Calculation Procedure for Heat, Mass and Momentum Transfer in Three-Dimensional Parabolic Flows, International Journal of Heat and Mass Transfer, 15 (1972), 1787-1806.
[14]     C. M. Rhie, W.L. Chow, Numerical Study of the Turbulent Flow Past an Airfoil with Trading Edge Separation, AIAA J., 21 (11) (1983) 1525-1535.
[15]     D. B. Spalding, A Novel Finite Difference Formulation for Differential Expressions Involving Both First and Second Derivatives, International-Journal for Numerical Methods in Engineering, 4 (1972), 551-559.
[16]     W. E. Ranz, W. R. Marshall, JR., Evaporation from Drops, Chemical Engineering Progress, 48 (3), (1952), 141-146
[17]     P. Bhattacharya, S. K. Saha, A. Yadav, P. E. Phelan, R. S. Prasher, Brownian dynamics simulation to determine the effective thermal conductivity of nanofluids, Journal of Applied Physics 95(11) (2004), 6492–6494.
[18]     W. M. Kays, M. E. Crawford, 1980, Convection Heat and Mass Transfer, 2nd ed., Mc Graw-Hill, New York
[19]      A. Bejan, Convection Heat Transfer, John Wiley& Sons, Hoboken, New Jersey, 4th ed. (2013)
[20]      J. P. Holman, in: Heat Transfer, eighth ed. McGraw-Hill Inc., New York (1997), pp. 218-282.
[21]     J. Y. Jung, H. S. Oh, H. Y. Kwak, Forced Convection Heat Transfer of Nanofluids in Microchannels, International Journal of Heat and Mass Transfer, 52 (2009), 466-472.
Volume 53, Issue 1
March 2022
Pages 66-82
  • Receive Date: 26 December 2021
  • Revise Date: 03 April 2022
  • Accept Date: 04 April 2022
  • First Publish Date: 04 April 2022