An adaptive piecewise linear finite element approximation to a linear method, the so-called nonlinear Chernoff formula, for the simplest two-phase Stefan problem in two dimensions are discussed. The full discretization amounts to solving linear positive definite symmetric systems, followed by an explicit nodewise algebraic correction. Consecutive meshes are highly graded and noncompatible. The pointwise information needed to generate a new mesh is extracted from both the discrete enthalpy U(n) and temperature THETA(n). For well-behaved discrete transition layers T(n) = { 0 < U(n) - THETA(n) < 1}, triangles are expected to be O(tau) within a thin strip of width O(tau1/2), the so-called refined region, which in turn must contain the transition layer. The local meshsize increases up to O(tau1/2) elsewhere; here T stands for the uniform time-step. The expected number of spatial degrees of freedom becomes O(tau-3/2), which compares quite favorably with that required without adaptivity, namely O(tau-2). Stability is examined in a number of Sobolev norms and then is used to derive a quasi-optimal rate of convergence of order essentially O(tau1/2) in the natural energy spaces. Numerical experiments illustrate the performance of the proposed method.

A fully discrete adaptive nonlinear Chernoff formula / R.H. Nochetto, M. Paolini, C. Verdi. - In: SIAM JOURNAL ON NUMERICAL ANALYSIS. - ISSN 0036-1429. - 30:4(1993), pp. 991-1014.

A fully discrete adaptive nonlinear Chernoff formula

C. Verdi
Ultimo
1993

Abstract

An adaptive piecewise linear finite element approximation to a linear method, the so-called nonlinear Chernoff formula, for the simplest two-phase Stefan problem in two dimensions are discussed. The full discretization amounts to solving linear positive definite symmetric systems, followed by an explicit nodewise algebraic correction. Consecutive meshes are highly graded and noncompatible. The pointwise information needed to generate a new mesh is extracted from both the discrete enthalpy U(n) and temperature THETA(n). For well-behaved discrete transition layers T(n) = { 0 < U(n) - THETA(n) < 1}, triangles are expected to be O(tau) within a thin strip of width O(tau1/2), the so-called refined region, which in turn must contain the transition layer. The local meshsize increases up to O(tau1/2) elsewhere; here T stands for the uniform time-step. The expected number of spatial degrees of freedom becomes O(tau-3/2), which compares quite favorably with that required without adaptivity, namely O(tau-2). Stability is examined in a number of Sobolev norms and then is used to derive a quasi-optimal rate of convergence of order essentially O(tau1/2) in the natural energy spaces. Numerical experiments illustrate the performance of the proposed method.
Settore MAT/08 - Analisi Numerica
Article (author)
File in questo prodotto:
Non ci sono file associati a questo prodotto.
Pubblicazioni consigliate

Caricamento pubblicazioni consigliate

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/2434/178595
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 11
  • ???jsp.display-item.citation.isi??? 7
social impact