Finite element approximation of a Stefan problem with degenerate Joule heating

Barrett, John W. ; Nürnberg, Robert

ESAIM: Mathematical Modelling and Numerical Analysis , Tome 38 (2004) no. 4, pp. 633-652

Résumé

We consider a fully practical finite element approximation of the following degenerate system

$$\frac{\partial }{\partial t}\rho \left(u\right)-\nabla .\left(\alpha \left(u\right)\nabla u\right)\ni \sigma \left(u\right){\left|\nabla \phi \right|}^2,\nabla .\left(\sigma \left(u\right)\nabla \phi \right)=0$$

subject to an initial condition on the temperature, $u$, and boundary conditions on both $u$ and the electric potential, $\phi$. In the above $\rho \left(u\right)$ is the enthalpy incorporating the latent heat of melting, $\alpha \left(u\right)>0$ is the temperature dependent heat conductivity, and $\sigma \left(u\right)\ge 0$ is the electrical conductivity. The latter is zero in the frozen zone, $u\le 0$, which gives rise to the degeneracy in this Stefan system. In addition to showing stability bounds, we prove (subsequence) convergence of our finite element approximation in two and three space dimensions. The latter is non-trivial due to the degeneracy in $\sigma \left(u\right)$ and the quadratic nature of the Joule heating term forcing the Stefan problem. Finally, some numerical experiments are presented in two space dimensions.

Détail

DOI : 10.1051/m2an:2004030

Classification : 35K55, 35K65, 35R35, 65M12, 65M60, 80A22
Keywords: Stefan problem, Joule heating, degenerate system, finite elements, convergence