Continuously deforming finite elements for the solution of parabolic problems, with and without phase change

Abstract

AbstractA number of transport problems are complicated by the presence of physically important transition zones where quantities exhibit steep gradients and special numerical care is required. When the location of such a transition zone changes as the solution evolves through time, use of a deforming numerical mesh is appropriate in order to preserve the proper numerical features both within the transition zone and at its boundaries.A general finite element solution method is described wherein the elements are allowed to deform continuously, and the effects of this deformation are accounted for exactly. The method is based on the Galerkin approximation in space, and uses finite difference approximations for the time derivatives. In the absence of element deformation, the method reduces to the conventional Galerkin formulation.The method is applied to the two‐phase Stefan problem associated with the melting and solidification of A substance. The interface between the solid and liquid phase form an internal moving boundary, and latent heat effects are accounted for in the associated boundary condition. By allowing continuous mesh deformation, as dictated by this boundary condition, the moving boundary always lies on element boundaries. This circumvents the difficulties inherent in interpolation of parameters and dependent variables across regions where those quantities change abruptly.Basis functions based on Hermite polynomials are used, to allow exact specification of the flux‐latent heat balance condition at the phase boundary. Analytic solutions for special cases provide tests of the method.