Finite Element Analysis of Solidification Process Using the Temperature-Enthalpy Relationship

온도-엔탈피 관계를 이용한 응고과정의 유한요소 해석

A finite element method is developed for calculating the temperature and enthalpy distribution and accordingly the solid, liquid and mushy zone in a three-dimensional body subjected to any heat boundary conditions. The method concurrently consider both temperature and enthalpy for consideration of the latent heat effect, differently from other methods of using a special energy balance equation for solving a mushy zone. The developed brick element has eight nodes with one degree of freedom at each node. The numerical method and procedure are verified using the results of one and two dimensional analytic solutions and by other researchers. It is shown that the present method presents a consistent and stable results in either abrupt or ranged phase change problems. Moreover, the numerical results by the present method are hardly effected by the calculation time steps which otherwise are difficult to determine in most phase change problems. Finally, as a three-dimensional application, a T-shaped body of a phase change is presented and the temperature and enthalpy variation along the time are solved.