Published in Int. J. Numer. Meth. Engng. Vol. 23 (1), pp. 99-119, 1986
doi: 10.1002/nme.1620230109


A finite element procedure for solving multidimensional phase change problems is described. The algorithm combines a temperature formulation with a finite element treatment of the differential equation and discontinuous integration within the two‐phase elements to avoid the necessity of regularization. A new criterion for the computation of the iteration matrix is proposed. It is based on a quasi‐Newton correction of the Jacobian matrix for conduction problems without change of phase. A set of test problems with exact solution is analysed and demonstrates that the procedure can accurately evaluate the front position and temperature history with a reasonable computational effort.

Back to Top

Document information

Published on 23/02/19
Submitted on 15/02/19

DOI: 10.1002/nme.1620230109
Licence: CC BY-NC-SA license

Document Score


Times cited: 41
Views 4
Recommendations 0

Share this document

claim authorship

Are you one of the authors of this document?