Technical Brief

Assessment of Basic Approaches to Numerical Modeling of Phase Change Problems—Accuracy, Efficiency, and Parallel Decomposition

[+] Author and Article Information
Tomas Mauder

Department of Thermodynamics,
Brno University of Technology,
Brno 616 00, Czech Republic
e-mail: mauder@fme.vutbr.cz

Pavel Charvat

Department of Thermodynamics,
Brno University of Technology,
Brno 616 00, Czech Republic
e-mail: charvat@fme.vutbr.cz

Josef Stetina

Department of Thermodynamics,
Brno University of Technology,
Brno 616 00, Czech Republic
e-mail: stetina@fme.vutbr.cz

Lubomir Klimes

Department of Thermodynamics,
Brno University of Technology,
Brno 616 00, Czech Republic
e-mail: klimes@fme.vutbr.cz

1Corresponding author.

Contributed by the Heat Transfer Division of ASME for publication in the JOURNAL OF HEAT TRANSFER. Manuscript received June 27, 2016; final manuscript received January 27, 2017; published online April 11, 2017. Editor: Portonovo S. Ayyaswamy.

J. Heat Transfer 139(8), 084502 (Apr 11, 2017) (5 pages) Paper No: HT-16-1425; doi: 10.1115/1.4036081 History: Received June 27, 2016; Revised January 27, 2017

The fast and accurate modeling of phase change is of a significant importance in many processes from steel casting to latent heat thermal energy storage. The paper presents a numerical case study on the transient 3D heat diffusion problem with phase change. Three different approaches to modeling of the solid–liquid phase change in combination with four commonly used numerical schemes are compared for their efficiency, accuracy, applicability, simplicity of implementation, and robustness. The possibility of parallel decomposition of the approaches is also discussed. The results indicate that the best accuracy was achieved with the second-order implicit methods, and the best efficiency was reached with the simple explicit methods.

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


Ockenden, J. R. , and Hodgkin, W. R. , 1975, Moving Boundary Value Problems in Heat Flow and Diffusion, Oxford University Press, Oxford, UK.
Tavakoli, R. , and Davami, P. , 2007, “ Unconditionally Stable Fully Explicit Finite Difference Solution of Solidification Problems,” Metall. Mater. Trans. B, 38(1), pp. 121–142. [CrossRef]
Douglas, J. , and Gunn, J. E. , 1964, “ A General Formulation of Alternating Direction Methods—Part I: Parabolic and Hyperbolic Problems,” Numer. Math., 6(1), pp. 428–453. [CrossRef]
Muhieddine, M. , Canot, É. , and March, R. , 2009, “ Various Approaches for Solving Problems in Heat Conduction With Phase Change,” Int. J. Finite, 3, pp. 1–20.
Liu, Y.-Ch. , and Chao, L.-S. , 2006, “ Modified Effective Specific Heat Method of Solidification Problems,” Mater. Trans. J., 47(11), pp. 2737–2744. [CrossRef]
Zhang, Y. , Du, K. , He, J. P. , Yang, L. , and Li, Y. J. , 2014, “ Impact Factors Analysis of the Enthalpy Method and the Effective Heat Capacity Method on the Transient Nonlinear Heat Transfer in Phase Change Materials (PCMs),” Numer. Heat Transfer A, 65(1), pp. 66–83. [CrossRef]
Li, Ch.-Y. , Garimella, S. V. , and Simpson, J. E. , 2003, “ Fixed-Grid Front-Tracking Algorithm for Solidification Problems—Part I: Method and Validation,” Numer. Heat Transfer, Part B, 43(2), pp. 117–141. [CrossRef]
Voller, V. , 1983, “ An Explicit Numerical Method to Track a Moving Phase Change Front,” Int. J. Heat Mass Transfer, 26(1), pp. 147–150. [CrossRef]
Hashemi, T. H. , and Sliepcevich, C. M. , 1968, “ A Numerical Method for Solving Two-Dimensional Problems of Heat Conduction With Change of Phase,” Chem. Eng. Prog. Symp. Ser., 63, pp. 34–41.
Lewisa, R. W. , and Ravindran, K. , 2000, “ Finite Element Simulation of Metal Casting,” Int. J. Numer. Methods Eng., 47(1–3), pp. 29–59. [CrossRef]
Dusinberre, G. M. , 1945, “ Numerical Methods for Transient Heat Flows,” Trans. ASME, 67, pp. 703–712.
Mochnacki, B. , and Lara, S. , 2003, “ Application of Generalized Finite Difference Method in Numerical Modelling of Moving Boundary Problems,” Sci. Res. Inst. Math. Comput. Sci., 2(1), pp. 129–143.
Minkowycz, W. J. , Sparrow, E. M. , Schneider, G. E. , and Pletcher, R. H. , 1988, Handbook of Numerical Heat Transfer, Wiley, New York.
Voller, V. R. , 1996, “ An Overview of Numerical Methods for Solving Phase Change Problems,” Advances in Numerical Heat Transfer, Vol. 1, W. J. Minkowycz and E. M. Sparrow , eds., Taylor & Francis, London, pp. 341–375.
Hahn, D. W. , and Ozisik, M. N. , 2012, Heat Conduction, 3rd ed., Wiley, Hoboken, NJ. [PubMed] [PubMed]
Yang, W. Y. , Cao, W. , Chung, T.-S. , and Morris, J. , 2005, Applied Numerical Methods Using Matlab, Wiley, Hoboken, NJ.
Otto, S. R. , and Denier, J. P. , 2005, An Introduction to Programming and Numerical Methods in MATLAB, Springer-Verlag, London.
Patankar, S. V. , 1980, Numerical Heat Transfer and Fluid Flow, McGraw-Hill, New York.
Larkin, B. K. , 1964, “ Some Stable Explicit Difference Approximations to the Diffusion Equation,” Math. Comput., 18(86), pp. 196–202. [CrossRef]
Barakat, H. Z. , and Clark, J. A. , 1966, “ On the Solution of Diffusion Equation by Numerical Methods,” ASME J. Heat Transfer, 88(4), pp. 421–427. [CrossRef]
Warming, R. F. , and Beam, R. M. , 1979, “ An Extension of A-Stability to Alternating Direction Implicit Methods,” BIT Numer. Math., 19(3), pp. 395–417. [CrossRef]
Wang, T.-Y. , and Chen, Ch. C.-P. , 2002, “ 3-D Thermal-ADI: A Linear-Time Chip Level Transient Thermal Simulator,” IEEE Trans. Comput.-Aided Des. Integr. Circuits Syst., 21(12), pp. 1434–1445. [CrossRef]
Rubitherm Technologies GmbH, 2016, Data sheet RT28HC, Berlin, accessed Dec. 8, 2016, http://www.rubitherm.eu/
Mauder, T. , Sandera, C. , and Stetina, J. , 2015, “ Optimal Control Algorithm for Continuous Casting Process by Using Fuzzy Logic,” Steel Res. Int., 86(7), pp. 785–798. [CrossRef]
Fu, W. S. , Wang, W. H. , and Huang, S. H. , 2012, “ An Investigation of Natural Convection of Three Dimensional Horizontal Parallel Plates From a Steady to an Unsteady Situation by a CUDA Computation Platform,” Int. J. Heat Mass Transfer, 55(17–18), pp. 4638–4650. [CrossRef]


Grahic Jump Location
Fig. 1

Results of 1D simulations



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