Finite element approximation of a Stefan problem with degenerate Joule heating

John W. Barrett; Robert Nürnberg

ESAIM: Mathematical Modelling and Numerical Analysis (2010)

  • Volume: 38, Issue: 4, page 633-652
  • ISSN: 0764-583X

Abstract

top
We consider a fully practical finite element approximation of the following degenerate system ∂ ∂ t ρ ( u ) - ∇ . ( α ( u ) ∇ u ) ∋ σ ( u ) | ∇ φ | 2 , ∇ . ( σ ( u ) ∇ φ ) = 0 subject to an initial condition on the temperature, u, and boundary conditions on both u and the electric potential, ϕ. In the above p(u) is the enthalpy incorporating the latent heat of melting, α(u) > 0 is the temperature dependent heat conductivity, and σ(u) > 0 is the electrical conductivity. The latter is zero in the frozen zone, u ≤ 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 σ(u) and the quadratic nature of the Joule heating term forcing the Stefan problem. Finally, some numerical experiments are presented in two space dimensions.

How to cite

top

Barrett, John W., and Nürnberg, Robert. "Finite element approximation of a Stefan problem with degenerate Joule heating." ESAIM: Mathematical Modelling and Numerical Analysis 38.4 (2010): 633-652. <http://eudml.org/doc/194231>.

@article{Barrett2010,
abstract = { We consider a fully practical finite element approximation of the following degenerate system $$ \{\frac\{\partial \}\{\partial t\}\} \rho(u) - \nabla . ( \,\alpha(u) \,\nabla u ) \ni \sigma(u)\,|\nabla\phi|^2 , \quad \nabla . (\, \sigma(u) \,\nabla \phi ) = 0 $$ subject to an initial condition on the temperature, u, and boundary conditions on both u and the electric potential, ϕ. In the above p(u) is the enthalpy incorporating the latent heat of melting, α(u) > 0 is the temperature dependent heat conductivity, and σ(u) > 0 is the electrical conductivity. The latter is zero in the frozen zone, u ≤ 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 σ(u) and the quadratic nature of the Joule heating term forcing the Stefan problem. Finally, some numerical experiments are presented in two space dimensions. },
author = {Barrett, John W., Nürnberg, Robert},
journal = {ESAIM: Mathematical Modelling and Numerical Analysis},
keywords = {Stefan problem; Joule heating; degenerate system; finite elements; convergence.; convergence},
language = {eng},
month = {3},
number = {4},
pages = {633-652},
publisher = {EDP Sciences},
title = {Finite element approximation of a Stefan problem with degenerate Joule heating},
url = {http://eudml.org/doc/194231},
volume = {38},
year = {2010},
}

TY - JOUR
AU - Barrett, John W.
AU - Nürnberg, Robert
TI - Finite element approximation of a Stefan problem with degenerate Joule heating
JO - ESAIM: Mathematical Modelling and Numerical Analysis
DA - 2010/3//
PB - EDP Sciences
VL - 38
IS - 4
SP - 633
EP - 652
AB - We consider a fully practical finite element approximation of the following degenerate system $$ {\frac{\partial }{\partial t}} \rho(u) - \nabla . ( \,\alpha(u) \,\nabla u ) \ni \sigma(u)\,|\nabla\phi|^2 , \quad \nabla . (\, \sigma(u) \,\nabla \phi ) = 0 $$ subject to an initial condition on the temperature, u, and boundary conditions on both u and the electric potential, ϕ. In the above p(u) is the enthalpy incorporating the latent heat of melting, α(u) > 0 is the temperature dependent heat conductivity, and σ(u) > 0 is the electrical conductivity. The latter is zero in the frozen zone, u ≤ 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 σ(u) and the quadratic nature of the Joule heating term forcing the Stefan problem. Finally, some numerical experiments are presented in two space dimensions.
LA - eng
KW - Stefan problem; Joule heating; degenerate system; finite elements; convergence.; convergence
UR - http://eudml.org/doc/194231
ER -

References

top
  1. J.W. Barrett and C.M. Elliott, A finite element method on a fixed mesh for the Stefan problem with convection in a saturated porous medium, in Numerical Methods for Fluid Dynamics, K.W. Morton and M.J. Baines Eds., Academic Press (London) (1982) 389–409.  
  2. J.W. Barrett and R. Nürnberg, Convergence of a finite element approximation of surfactant spreading on a thin film in the presence of van der Waals forces. IMA J. Numer. Anal.24 (2004) 323–363.  
  3. C.M. Elliott, On the finite element approximation of an elliptic variational inequality arising from an implicit time discretization of the Stefan problem. IMA J. Numer. Anal.1 (1981) 115–125.  
  4. C.M. Elliott, Error analysis of the enthalpy method for the Stefan problem. IMA J. Numer. Anal.7 (1987) 61–71.  
  5. C.M. Elliott and S. Larsson, A finite element model for the time-dependent Joule heating problem. Math. Comp.64 (1995) 1433–1453.  
  6. R.F. Gariepy, M. Shillor and X. Xu, Existence of generalized weak solutions to a model for in situ vitrification. European J. Appl. Math.9 (1998) 543–559.  
  7. S.S. Koegler and C.H. Kindle, Modeling of the in situ vitrification process. Amer. Ceram. Soc. Bull.70 (1991) 832–835.  
  8. J. Simon, Compact sets in the space Lp(0,T;B). Ann. Math. Pura. Appl.146 (1987) 65–96.  
  9. X. Xu, A compactness theorem and its application to a system of partial differential equations. Differential Integral Equations9 (1996) 119–136.  
  10. X. Xu, Existence for a model arising from the in situ vitrification process. J. Math. Anal. Appl.271 (2002) 333–342.  

NotesEmbed ?

top

You must be logged in to post comments.

To embed these notes on your page include the following JavaScript code on your page where you want the notes to appear.

Only the controls for the widget will be shown in your chosen language. Notes will be shown in their authored language.

Tells the widget how many notes to show per page. You can cycle through additional notes using the next and previous controls.

    
                

                
Note: Best practice suggests putting the JavaScript code just before the closing </body> tag.