A simple and efficient adaptive local mesh refinement algorithm is devised and analyzed for two-phase Stefan problems in 2D. A typical triangulation is coarse away from the discrete interface, where discretization parameters satisfy a parabolic relation, whereas it is locally refined in the vicinity of the discrete interface so that the relation becomes hyperbolic. Several numerical tests are performed on the computed temperature to extract information about its first and second derivatives as well as to predict discrete free boundary locations. Mesh selection is based upon equidistributing pointwise interpolation errors between consecutive meshes and imposing that discrete interfaces belong to the so-called refined region. Consecutive meshes are not compatible in that they are not produced by enrichment or coarsening procedures but rather regenerated. A general theory for interpolation between noncompatible meshes is set up in L(p)-based norms. The resulting scheme is stable in various Sobolev norms and necessitates fewer spatial degrees of freedom than previous practical methods on quasi-uniform meshes, namely O(tau-3/2) as opposed to O(tau-2), to achieve the same global asymptotic accuracy; here tau > 0 is the (uniform) time step. A rate of convergence of essentially O(tau-1/2) is derived in the natural energy spaces provided the total number of mesh changes is restricted to O(tau-1/2), which in turn is compatible with the mesh selection procedure. An auxiliary quasi-optimal pointwise error estimate for the Laplace operator is proved as well. Numerical results illustrate the scheme's efficiency in approximating both solutions and interfaces.

An adaptive finite element method for two-phase stefan problems in two space dimensions. part i: stability and error estimates / R. Nochetto, M. Paolini, C. Verdi. - In: MATHEMATICS OF COMPUTATION. - ISSN 0025-5718. - 57:195(1991), pp. 73-108.

An adaptive finite element method for two-phase stefan problems in two space dimensions. part i: stability and error estimates

C. Verdi
Ultimo
1991

Abstract

A simple and efficient adaptive local mesh refinement algorithm is devised and analyzed for two-phase Stefan problems in 2D. A typical triangulation is coarse away from the discrete interface, where discretization parameters satisfy a parabolic relation, whereas it is locally refined in the vicinity of the discrete interface so that the relation becomes hyperbolic. Several numerical tests are performed on the computed temperature to extract information about its first and second derivatives as well as to predict discrete free boundary locations. Mesh selection is based upon equidistributing pointwise interpolation errors between consecutive meshes and imposing that discrete interfaces belong to the so-called refined region. Consecutive meshes are not compatible in that they are not produced by enrichment or coarsening procedures but rather regenerated. A general theory for interpolation between noncompatible meshes is set up in L(p)-based norms. The resulting scheme is stable in various Sobolev norms and necessitates fewer spatial degrees of freedom than previous practical methods on quasi-uniform meshes, namely O(tau-3/2) as opposed to O(tau-2), to achieve the same global asymptotic accuracy; here tau > 0 is the (uniform) time step. A rate of convergence of essentially O(tau-1/2) is derived in the natural energy spaces provided the total number of mesh changes is restricted to O(tau-1/2), which in turn is compatible with the mesh selection procedure. An auxiliary quasi-optimal pointwise error estimate for the Laplace operator is proved as well. Numerical results illustrate the scheme's efficiency in approximating both solutions and interfaces.
degenerate parabolic problems; plane polygonal domains; approximation
Settore MAT/08 - Analisi Numerica
Article (author)
File in questo prodotto:
File Dimensione Formato  
S0025-5718-1991-1079028-X.pdf

accesso aperto

Tipologia: Publisher's version/PDF
Dimensione 4.05 MB
Formato Adobe PDF
4.05 MB Adobe PDF Visualizza/Apri
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/180721
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 43
  • ???jsp.display-item.citation.isi??? 42
social impact