To read this content please select one of the options below:

Numerical solution of Stefan problem with variable space grid method based on mixed finite element/finite difference approach

Milos Ivanovic (Department of Mathematics and Informatics, University of Kragujevac, Kragujevac, Serbia)
Marina Svicevic (Department of Mathematics and Informatics, University of Kragujevac, Kragujevac, Serbia)
Svetislav Savovic (Department of Mathematics and Informatics, University of Kragujevac, Kragujevac, Serbia and Department of Mathematics, City University of Hong Kong, Kowloon, Hong Kong)

International Journal of Numerical Methods for Heat & Fluid Flow

ISSN: 0961-5539

Article publication date: 4 December 2017

628

Abstract

Purpose

The purpose of this paper is to improve the accuracy and stability of the existing solutions to 1D Stefan problem with time-dependent Dirichlet boundary conditions. The accuracy improvement should come with respect to both temperature distribution and moving boundary location.

Design/methodology/approach

The variable space grid method based on mixed finite element/finite difference approach is applied on 1D Stefan problem with time-dependent Dirichlet boundary conditions describing melting process. The authors obtain the position of the moving boundary between two phases using finite differences, whereas finite element method is used to determine temperature distribution. In each time step, the positions of finite element nodes are updated according to the moving boundary, whereas the authors map the nodal temperatures with respect to the new mesh using interpolation techniques.

Findings

The authors found that computational results obtained by proposed approach exhibit good agreement with the exact solution. Moreover, the results for temperature distribution, moving boundary location and moving boundary speed are more accurate than those obtained by variable space grid method based on pure finite differences.

Originality/value

The authors’ approach clearly differs from the previous solutions in terms of methodology. While pure finite difference variable space grid method produces stable solution, the mixed finite element/finite difference variable space grid scheme is significantly more accurate, especially in case of high alpha. Slightly modified scheme has a potential to be applied to 2D and 3D Stefan problems.

Keywords

Citation

Ivanovic, M., Svicevic, M. and Savovic, S. (2017), "Numerical solution of Stefan problem with variable space grid method based on mixed finite element/finite difference approach", International Journal of Numerical Methods for Heat & Fluid Flow, Vol. 27 No. 12, pp. 2682-2695. https://doi.org/10.1108/HFF-11-2016-0443

Publisher

:

Emerald Publishing Limited

Copyright © 2017, Emerald Publishing Limited

Related articles