On Heat Conduction with Phase Change: Accurate Explicit Numerical Method


1 Service Transferts, Interfaces et Procédés-TIPs, Faculté des Sciences Appliquées/Ecole Polytechnique, Université Libre de Bruxelles-ULB, CP165/67, avenue F.D. Roosevelt 50, 1050 Bruxelles, Belgique

2 Laboratoire de Mécanique des Fluides Théorique et Appliquée, Faculté de Physique, Université des Sciences et de la Technologie Houari Boumediene-USTHB, B.P. 32, El-Alia, 16111, Alger, Algérie

3 LUNAM Université, Université de Nantes, CNRS, GEPEA UMR- 6144, CRTT, BP 406 44602 Saint-Nazaire Cedex – France


The one phase Stefan problem is discussed using the Goodman HBI method and an explicit numerical method including modified boundary immobilization scheme. The main advantage of the HBI method lie in the remarkable association of simplicity, flexibility and acceptable accuracy which an error less than 2.5% in predicting the moving front location for Stefan number less than unity which covers most usual isothermal phase change material. An accurate explicit numerical model to track the moving front in Stefan-like problems is provided. The scheme is obtained using the variable space step method based on variable domain. The method is developed using central difference approximations to replace spatial and temporal derivatives. Furthermore, iterative procedure, in numerical calculation, is avoided by introducing simple assumptions. The numerical results show that the accuracy of the method has been considerably improved without additional computational cost.