Stability and convergence of two discrete schemes for a degenerate solutal non-isothermal phase-field model
Dpto. E.D.A.N., University of Sevilla,
Aptdo. 1160, 41080 Sevilla, Spain. email@example.com; firstname.lastname@example.org
Revised: 21 May 2008
Revised: 1 October 2008
We analyze two numerical schemes of Euler type in time and C0 finite-element type with -approximation in space for solving a phase-field model of a binary alloy with thermal properties. This model is written as a highly non-linear parabolic system with three unknowns: phase-field, solute concentration and temperature, where the diffusion for the temperature and solute concentration may degenerate. The first scheme is nonlinear, unconditionally stable and convergent. The other scheme is linear but conditionally stable and convergent. A maximum principle is avoided in both schemes, using a truncation operator on the L2 projection onto the finite element for the discrete concentration. In addition, for the model when the heat conductivity and solute diffusion coefficients are constants, optimal error estimates for both schemes are shown based on stability estimates.
Mathematics Subject Classification: 35Q72 / 35K65 / 65M12 / 65M60
Key words: Phase-field models / diffuse interface model / solidification process / degenerate parabolic systems / backward Euler schemes / finite elements / stability / convergence / error estimates.
© EDP Sciences, SMAI, 2009