An enthalpy-based finite element method for nonlinear heat problems involving phase change

Abstract Within the framework of the finite element method, we present in this paper an efficient algorithm for solving nonlinear heat problems involving phase change. In particular, two problems are considered in this work: the stationary convection–diffusion problem, and the classical transient heat problem. The mathematical models used to solve these problems are based upon enthalpy formulations. The algorithmic design is based on a Newton-type iterative procedure for the stationary problem, and for the transient one, a combination with classical finite difference schemes in time is performed. The proposed phase change algorithmic treatment is applicable for both the situations in which the latent heat takes place over a temperature range or at fixed temperature. Hence, for this latter situation, a regularization over a narrow temperature range is not necessary. The numerical implementation is discussed in detail where we include also the other possible nonlinearities, namely, the temperature dependence of the thermal conductivity and the radiation-type boundary conditions. Finally, a set of representative numerical simulations to illustrate the effectiveness of the proposed method is given.

[1]  J. R Ockendon,et al.  Moving boundary problems in heat flow and diffusion : being the proceedings of the conference held at the University of Oxford, 25-27 March 1974 , 1976 .

[2]  Boris Rubinsky,et al.  Numerical Computation Using Finite Elements for the Moving Interface in Heat Transfer Problems with Phase Transformation , 1983 .

[3]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[4]  P. Humbert.  CESAR-LCPC UN CODE GENERAL DE CALCUL PAR ELEMENTS FINIS , 1989 .

[5]  Maurizio Paolini,et al.  Finite element approximations of singular parabolic problems , 1988 .

[6]  Roland W. Lewis,et al.  An improved algrorithm for heat conduction problems with phase change , 1978 .

[7]  M. Cross,et al.  An enthalpy method for convection/diffusion phase change , 1987 .

[8]  Bernhard A. Schrefler,et al.  Numerical methods in heat transfer , 1983 .

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

[10]  Pekka Neittaanmäki.  Numerical Methods for Free Boundary Problems , 1991 .

[11]  E. Isaacson,et al.  Stefan-like problems , 1950 .

[12]  R. E. White.  A modified finite difference scheme for the Stefan problem , 1983 .

[13]  Antonio Fasano,et al.  Numerical solution of phase-change problems , 1973 .

[14]  K. Bathe.  Finite Element Procedures , 1995 .

[15]  Kumar K. Tamma,et al.  Recent advances, trends and new perspectives via enthalpy-based finite element formulations for applications to solidification problems , 1990 .

[16]  D. R. Atthey.  A Finite Difference Scheme for Melting Problems , 1974 .

[17]  A. Segal,et al.  Comparison of finite element techniques for solidification problems , 1986 .

[18]  Paul T. Boggs,et al.  Moving boundary problems , 1978 .

[19]  J. Z. Zhu,et al.  The finite element method , 1977 .

[20]  Roland W. Lewis,et al.  Finite element solution of non‐linear heat conduction problems with special reference to phase change , 1974 .

[21]  Boris Rubinsky,et al.  A finite element method for the solution of one-dimensional phase change problems , 1981 .

[22]  Antonio Fasano,et al.  Free boundary problems : theory and applications , 1983 .

[23]  Ricardo H. Nochetto,et al.  Energy error estimates for a linear scheme to approximate nonlinear parabolic problems , 1987 .

[24]  Klaus-Jürgen Bathe,et al.  AN EFFICIENT ALGORITHM FOR ANALYSIS OF NONLINEAR HEAT TRANSFER WITH PHASE CHANGES , 1982 .

[25]  Sergio Idelsohn,et al.  A temperature‐based finite element solution for phase‐change problems , 1986 .