Numerical simulation of moving boundary problem with moving phase change material and size-dependent thermal conductivity

Abstract

The Keller box method is applied to the one-phase moving boundary problem with moving phase change material, size-dependent thermal conductivity, and periodic boundary conditions. The phase transition process is allowed to occur when the material is forced to move in one direction or the other at a constant speed. The boundary immobilization method is applied to immobilize the moving boundary, and for the numerical approximation of moving boundary problem, we proposed Keller box method. Keller box method also accommodates the non-linearity in thermal conductivity and Stefan condition. Using the convergence analysis, the proposed scheme obtains stability, as well as second-order accuracy for both spatial and temporal directions, under reasonable conditions. To validate the proposed numerical scheme, we have considered a particular case of this problem having a similarity solution. It is found that the numerical results obtained by the Keller box method have good agreement with the similarity solution and also verified that the computational rate of convergence of our scheme is two. The effects of various parameters and size-dependent thermal conductivity on the position of moving boundary are also investigated.

Keywords

Moving boundary problem moving phase change material size-dependent thermal conductivity boundary immobilization method Keller box method

Get full access to this article

View all access options for this article.

References

Stefan

. Über die theorie der eisbildung, insbesondere lquouber die eisbildung im polarmeere. Ann Phys 1891; 278: 269–286.

Crank

. Free and moving boundary problems. New York: Oxford University Press, 1987.

Turkyilmazoglu

. Stefan problems for moving phase change materials and multiple solutions. Int J Thermal Sci 2018; 126: 67–73.

Font

. A one-phase Stefan problem with size-dependent thermal conductivity. Appl Math Model 2018; 63: 172–178.

Chaurasiya

Kumar

Rai

, et al. A computational solution of a phase-change material in the presence of convection under the most generalized boundary condition. Thermal Sci Eng Progr 2020; 20: 100664.

Jitendra, Rai KN and Singh J. Wavelet based numerical approach of non-classical moving boundary problem with convection effect and variable latent heat under the most generalized boundary conditions. Eur J Mech-B/Fluids 2021; 87: 1–11.

Bollati

Briozzo

. Stefan problems for the diffusion–convection equation with temperature-dependent thermal coefficients. Int J Non Linear Mech 2021; 134: 103732.

Parhizi

Jain

. Eigenfunction-based solution for one-dimensional solid–liquid phase change heat transfer problems with advection. Int J Thermal Sci 2022; 172: 107262.

Nandi

Sanyasiraju

YVS

. A second-order accurate fixed-grid method for multi-dimensional Stefan problem with moving phase change materials. Appl Math Comput 2022; 416: 126719.

10.

Chaurasiya

Rai

Singh

. Heat transfer analysis for the solidification of a binary eutectic system under imposed movement of the material. J Therm Anal Calorim 2021; 147: 1–18.

11.

Kumar A and Rajeev. A Stefan problem with moving phase change material, variable thermal conductivity and periodic boundary condition. Appl Math Comput 2020; 386: 125490.

12.

Landau

. Heat conduction in a melting solid. Q Appl Math 1950; 8: 81–94.

13.

Friedman

. Partial differential equations of parabolic type. New York: Courier Dover Publications, 2008.

14.

Ladyvzenskaja

Solonnikov

Ural’ceva

. Linear and quasi-linear equations of parabolic type, volume 23. American Mathematical Soc., 1988.

15.

Keller

. A new difference scheme for parabolic problems. In Numerical Solution of Partial Differential Equations–II, pages 327–350. Elsevier, 1971.

16.

Strikwerda

. Finite difference schemes and partial differential equations. Philadelphia: SIAM, 2004.