An enriched finite element algorithm for the implicit simulation of the Stefan problem

An implicit enriched finite element algorithm is proposed to simulate heat transfer involving isothermal phase changes. This technique is based on a mixed variational formulation discretized by means of an enriched finite element approximation of the enthalpy in space. The interface is implicitly described without coupling with an interface-capturing technique. The time integration is carried out with an implicit (backward) Euler algorithm in time. Two examples in 1D and 2D clearly evidence the efficiency of the method developed.

RésuméUn algorithme implicite de type éléments finis enrichis est proposé pour simuler des transferts de chaleur avec des changements de phase isothermes. Cette technique repose sur une formulation variationnelle mixte, discrétisée au moyen dʼune approximation de lʼenthalpie par éléments finis enrichis. Lʼinterface est décrite de manière implicite sans couplage avec une technique spécifique de capture dʼinterface. Lʼintégration temporelle est réalisée par un algorithme dʼEuler purement implicite en temps. Deux exemples 1D et 2D mettent clairement en évidence lʼefficacité de la méthode développée.