An analytical method for the solution of two phase Stefan problem in cylindrical geometry

Two phase Stefan problem was solved using analytical method in cylindrical domain. To solve governing equations Eigen conditions were formulated by using separation of variable technique. Eigenvalues of the eigencondition were obtained by applying corresponding boundary conditions for liquid and solid phase. Eigenvalues are graphically validated by using window size method in Mathematica. It is noted radial eigenvalues are free from imaginary values. Interface equation obtained from this method were solved and analyzed by varying the Stefan number and introducing the forced and natural convection. Conduction and convection heat transfer mechanism was studied and results obtained by varying thermal diffusivity, thermal conductivity and Stefan number were discussed. Natural convection effects were studied by introducing Rayleigh number and results showed Stefan number has significant effect than Rayleigh number during Phase transition process. Furthermore, eigen function expansion Method was compared with exact solution of Exponential Integral function method and results showed good agreement for Q = 1.

[1]  M. Rahimi,et al.  Experimental and computational evolution of a shell and tube heat exchanger as a PCM thermal storage system , 2014 .

[2]  Yogesh Jaluria,et al.  A COMPARISON OF DIFFERENT SOLUTION METHODOLOGIES FOR MELTING AND SOLIDIFICATION PROBLEMS IN ENCLOSURES , 1993 .

[3]  Prashant K. Jain,et al.  An exact analytical solution for two-dimensional, unsteady, multilayer heat conduction in spherical coordinates , 2010 .

[4]  Kacem Mansouri,et al.  Analytical solution of heat transfer in a shell-and-tube latent thermal energy storage system , 2015 .

[5]  T. Saitoh,et al.  Numerical Method for Multi-Dimensional Freezing Problems in Arbitrary Domains , 1978 .

[6]  Ephraim M Sparrow,et al.  Analysis of melting in the presence of natural convection in the melt region , 1977 .

[7]  J. Prusa,et al.  Melting and Freezing , 1989 .

[8]  Blackmore,et al.  Discrete-ordinate method of solution of Fokker-Planck equations with nonlinear coefficients. , 1985, Physical review. A, General physics.

[9]  Stewart Paterson,et al.  Propagation of a Boundary of Fusion , 1952, Proceedings of the Glasgow Mathematical Association.

[10]  H. G. Landau,et al.  Heat conduction in a melting solid , 1950 .

[11]  A. M. Meirmanov,et al.  The Stefan Problem , 1992 .

[12]  Majid Ali,et al.  A novel method for analytical solution of transient heat conduction and Stefan problem in cylindrical coordinate , 2016, 2016 13th International Bhurban Conference on Applied Sciences and Technology (IBCAST).

[13]  R. M. Furzeland,et al.  A Comparative Study of Numerical Methods for Moving Boundary Problems , 1980 .

[14]  Michael F. Malone,et al.  Analysis of Two-Dimensional Diffusion-Controlled Moving Boundary Problems , 1975 .

[15]  L. S. Yao,et al.  Effects of Natural Convection in the Melted Region Around a Heated Horizontal Cylinder , 1980 .

[16]  Sergei Levendorskii,et al.  THE EIGENFUNCTION EXPANSION METHOD IN MULTI‐FACTOR QUADRATIC TERM STRUCTURE MODELS , 2007 .

[17]  L. Burmeister.  Convective heat transfer , 1983 .

[18]  Ephraim M Sparrow,et al.  Effect of Subcooling on Cylindrical Melting , 1978 .

[19]  Suneet Singh,et al.  Analytical solution to transient heat conduction in polar coordinates with multiple layers in radial direction , 2008 .

[20]  Bernie D. Shizgal,et al.  A Gaussian quadrature procedure for use in the solution of the Boltzmann equation and related problems , 1981 .

[21]  Kamal Abdel Radi Ismail,et al.  Numerical solution of the phase change problem around a horizontal cylinder in the presence of natural convection in the melt region , 2003 .