A grid based ADI method for the problem of two phase solidification

Abstract This paper proposes a transformation based, unconditionally stable, Alternating Direction Implicit (ADI) scheme for solving two-phase Stefan problems of solidification in arbitrary bounded domains. The governing equations of each phase are transformed, from a complex physical domain to a fixed rectangular domain, using body-fitted coordinates. ADI method is used to solve the transformed equations of each phase separately. The unconditional stability of the proposed ADI scheme is discussed numerically using von-Neumann method. Several numerical experiments are carried out for the case of stable solidification to verify the applicability of the proposed method. An excellent agreement has been found between the numerically generated values and the exact/existing solutions. Further, the developed scheme has also been tested on the problems of unstable solidification with mild surface tension and kinetic mobility. Once again the interface location with time has been computed very accurately.

[1]  J. Caldwell,et al.  Numerical methods for one‐dimensional Stefan problems , 2004 .

[2]  W. Shyy,et al.  Computation of Solid-Liquid Phase Fronts in the Sharp Interface Limit on Fixed Grids , 1999 .

[3]  Abhishek Kumar Singh,et al.  A Stefan problem with variable thermal coefficients and moving phase change material , 2019, Journal of King Saud University - Science.

[4]  D. Juric,et al.  A Front-Tracking Method for Dendritic Solidification , 1996 .

[5]  J. Strain.  Linear stability of planar solidification fronts , 1988 .

[6]  Svetislav Savović,et al.  Finite difference solution of one-dimensional Stefan problem with periodic boundary conditions , 2003 .

[7]  J. Otero,et al.  Constant mass model for the liquid-solid phase transition on a one-dimensional Stefan problem: Transient and steady state regimes , 2017 .

[8]  M. D. Simone,et al.  Analytical model for solidification and melting in a finite PCM in steady periodic regime , 2015 .

[9]  Roland Herzog,et al.  Optimal Control of the Classical Two-Phase Stefan Problem in Level Set Formulation , 2011, SIAM J. Sci. Comput..

[10]  R. Sekerka,et al.  Stability of a Planar Interface During Solidification of a Dilute Binary Alloy , 1964 .

[11]  Suresh V. Garimella,et al.  FIXED-GRID FRONT-TRACKING ALGORITHM FOR SOLIDIFICATION PROBLEMS, PART I: METHOD AND VALIDATION , 2003 .

[12]  K. A. Rathjen,et al.  Heat Conduction With Melting or Freezing in a Corner , 1971 .

[13]  Baskar Ganapathysubramanian,et al.  Modelling dendritic solidification with melt convection using the extended finite element method , 2006, J. Comput. Phys..

[14]  M. Hennessy,et al.  The one-dimensional Stefan problem with non-Fourier heat conduction , 2019, International Journal of Thermal Sciences.

[15]  R. Almgren.  Variational algorithms and pattern formation in dendritic solidification , 1993 .

[16]  James Caldwell,et al.  Spherical solidification by the enthalpy method and the heat balance integral method , 2000 .

[17]  Exact solution for Stefan problem with general power-type latent heat using Kummer function , 2015 .

[18]  S. Osher,et al.  A Simple Level Set Method for Solving Stefan Problems , 1997, Journal of Computational Physics.

[19]  J. Caldwell,et al.  Starting solutions for the boundary immobilization method , 2004 .

[20]  D. Mazzeo,et al.  Thermal field and heat storage in a cyclic phase change process caused by several moving melting and solidification interfaces in the layer , 2018, International Journal of Thermal Sciences.

[21]  Sarah L. Mitchell,et al.  Application of Standard and Refined Heat Balance Integral Methods to One-Dimensional Stefan Problems , 2010, SIAM Rev..

[22]  M. Oxley,et al.  The boundary element method applied to moving boundary problems , 1990 .

[23]  A third order accurate discontinuous finite element method for the one-dimensional Stefan problem , 1979 .

[24]  V. Voller,et al.  An analytical solution for a Stefan problem with variable latent heat , 2004 .

[25]  Luisa F. Cabeza,et al.  Experimental validation of the exact analytical solution to the steady periodic heat transfer problem in a PCM layer , 2017 .

[26]  Bruno Welfert,et al.  Stability of ADI schemes applied to convection-diffusion equations with mixed derivative terms , 2007 .

[27]  X. Duan,et al.  An improved layered thermal resistance model for solid-liquid phase change time estimation , 2020 .

[28]  Willem Hundsdorfer,et al.  Accuracy and stability of splitting with stabilizing corrections , 2002 .

[29]  Joe F. Thompson,et al.  Automatic numerical generation of body-fitted curvilinear coordinate system for field containing any number of arbitrary two-dimensional bodies , 1974 .

[30]  George Beckett,et al.  A moving mesh finite element method for the two-dimensional Stefan problems , 2001 .

[31]  Rizwan-uddin.  A Nodal Method for Phase Change Moving Boundary Problems , 1999 .

[32]  Sarah L. Mitchell,et al.  On the numerical solution of two-phase Stefan problems with heat-flux boundary conditions , 2014, J. Comput. Appl. Math..

[33]  Mustafa Turkyilmazoglu,et al.  Stefan problems for moving phase change materials and multiple solutions , 2018 .

[34]  Guillermo Marshall.  A front tracking method for one dimensional moving boundary problems , 1986 .

[35]  Sarah L. Mitchell,et al.  Finite-difference methods with increased accuracy and correct initialization for one-dimensional Stefan problems , 2009, Appl. Math. Comput..