Natural Convection Numeric Simulation on Metal Freezing Using Differential Method

Heri Suprianto, Eko Prasetyo Budiana, Purwadi Joko Widodo


The research of modeling of natural convection in metal solidification process with finite different method was conducted to determine temperature distribution and fluid flow profil with variations value Rayleigh number. The research conducted by solving governing equation of natural convection with finite difference approximation. Governing equation of natural convection consist of continuity equation, momentum equations, and energy equation. The ADI (Alternating Directional Implicit) method was used to discriteze for governing equation of natural convection. Finite difference method was written in Fortran language whereas the temperature distribution and fluid flow profile were visualized with Matlab software. The results of this research was validated by comparing the results obtained with Rajiv Sampath research. Comparison of the results of research showed good agreement. The result showed that solidification process occurs faster at Ra 10^4 compared with 10^5 and 10^6

Full Text:



Belhamadia, Y., Kane, A.S., Fortin, A., 2012, An Enhanced Mathematical Model for Phase Change Problem With Natural Convection, International Journal of Numerical Analysis and Modeling, Vol.3, 192-206.

Brandes, E.A., Brook, G.B., 1992, Smithells Metals Reference Book, 7th edition, Butterworth-Heinemann, Oxford.

Cengel, A.Y., 2008, Heat Transfer A Practical Approach Second Edition, McGraw-Hill, New York.

Chen, Y., Im, Y.T., Yoo, J., 1995, Finite Element Analaysis of Solidification of Aluminum with Natrural Convection, Journal of Material Processing Technology, Vol.52, 592-609, Elsevier.

Hoffman, K.A., 2000, Computation Fluid Dynamics For Engineering Volume I, Engineering System TM , Kansas USA.

Holman, J.P., 2010, Heat Transfer Tenth Edition, McGraw-Hill, New York.

Incropera, F.M., 1996, Introduction to Heat Transfer, John Wiley & Sons, USA.

McDaniel, D. J., Zabaras, N., 1994, A least-squares front-tracking finite element method analysis of phase change with natural convection, International Journal Numerical Method Engineering, Vol. 37, 2755-2777.


  • There are currently no refbacks.

Mekanika: Majalah Ilmiah Mekanika

ISSN: 1412-7962 || eISSN: 2579-3144

Address :  Jl. Ir Sutami no 36 A, Building I, Faculty of Engineering, Universitas Sebelas Maret, Surakarta.

Phone    :  +62271632163 

email     :


Creative Commons License
This work is licensed under a Creative Commons Attribution 4.0 International License