We aim to approximate contrast problems by means of a numerical scheme which does not require that the computational mesh conforms with the discontinuity between coefficients. We focus on the approximation of diffusion–reaction equations in the framework of finite elements. In order to improve the unsatisfactory behavior of Lagrangian elements for this particular problem, we resort to an enriched approximation space, which involves elements cut by the interface. Firstly, we analyze the H1-stability of the finite element space with respect to the position of the interface. This analysis, applied to the conditioning of the discrete system of equations, shows that the scheme may be ill posed for some configurations of the interface. Secondly, we propose a stabilization strategy, based on a scaling technique, which restores the standard properties of a Lagrangian finite element space and results to be very easily implemented. We also address the behavior of the scheme with respect to large contrast problems ending up with a choice of Nitscheʼs penalty terms such that the extended finite element scheme with penalty is robust for the worst case among small sub-elements and large contrast problems. The theoretical results are finally illustrated by means of numerical experiments.

An unfitted interface penalty method for the numerical approximation of contrast problems

ZUNINO, PAOLO;CATTANEO, LAURA;
2011

Abstract

We aim to approximate contrast problems by means of a numerical scheme which does not require that the computational mesh conforms with the discontinuity between coefficients. We focus on the approximation of diffusion–reaction equations in the framework of finite elements. In order to improve the unsatisfactory behavior of Lagrangian elements for this particular problem, we resort to an enriched approximation space, which involves elements cut by the interface. Firstly, we analyze the H1-stability of the finite element space with respect to the position of the interface. This analysis, applied to the conditioning of the discrete system of equations, shows that the scheme may be ill posed for some configurations of the interface. Secondly, we propose a stabilization strategy, based on a scaling technique, which restores the standard properties of a Lagrangian finite element space and results to be very easily implemented. We also address the behavior of the scheme with respect to large contrast problems ending up with a choice of Nitscheʼs penalty terms such that the extended finite element scheme with penalty is robust for the worst case among small sub-elements and large contrast problems. The theoretical results are finally illustrated by means of numerical experiments.
File in questo prodotto:
File Dimensione Formato  
APNUM-ZCC.pdf

Accesso riservato

: Post-Print (DRAFT o Author’s Accepted Manuscript-AAM)
Dimensione 801.64 kB
Formato Adobe PDF
801.64 kB Adobe PDF   Visualizza/Apri

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: http://hdl.handle.net/11311/608728
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 26
  • ???jsp.display-item.citation.isi??? 24
social impact