We propose a family of preconditioners for linear systems of equations arising from a piecewise polynomial symmetric interior penalty discontinuous Galerkin discretization of H(curl,?)-elliptic boundary value problems on conforming meshes. The design and analysis of the proposed preconditioners rely on the auxiliary space method (ASM) employing an auxiliary space of H(curl,?)-conforming finite element functions together with a relaxation technique (local smoothing). On simplicial meshes, the proposed preconditioner enjoys asymptotic optimality with respect to mesh refinement. It is also robust with respect to jumps in the coefficients ? and b in the second-and zeroth-order parts of the operator, respectively, except when the problem changes from curl-dominated to reaction-dominated and vice versa. On quadrilateral/hexahedral meshes some of the proposed ASM solvers may fail, since the related H(curl,?)-conforming finite element space does not provide a spectrally accurate discretization. Extensive numerical experiments are included to verify the theory and assess the performance of the preconditioners.

Auxiliary space preconditioners for SIP-DG discretizations of H(curl)-elliptic problems with discontinuous coefficients

2017

Abstract

We propose a family of preconditioners for linear systems of equations arising from a piecewise polynomial symmetric interior penalty discontinuous Galerkin discretization of H(curl,?)-elliptic boundary value problems on conforming meshes. The design and analysis of the proposed preconditioners rely on the auxiliary space method (ASM) employing an auxiliary space of H(curl,?)-conforming finite element functions together with a relaxation technique (local smoothing). On simplicial meshes, the proposed preconditioner enjoys asymptotic optimality with respect to mesh refinement. It is also robust with respect to jumps in the coefficients ? and b in the second-and zeroth-order parts of the operator, respectively, except when the problem changes from curl-dominated to reaction-dominated and vice versa. On quadrilateral/hexahedral meshes some of the proposed ASM solvers may fail, since the related H(curl,?)-conforming finite element space does not provide a spectrally accurate discretization. Extensive numerical experiments are included to verify the theory and assess the performance of the preconditioners.
2017
Istituto di Matematica Applicata e Tecnologie Informatiche - IMATI -
auxiliary space preconditioning
discontinuous coefficients
discontinuous Galerkin methods
H(curl?)-elliptic problems
File in questo prodotto:
File Dimensione Formato  
prod_371708-doc_124157.pdf

solo utenti autorizzati

Descrizione: Auxiliary space preconditioners for SIP-DG discretizations of H(curl)-elliptic problems
Tipologia: Versione Editoriale (PDF)
Dimensione 1.64 MB
Formato Adobe PDF
1.64 MB Adobe PDF   Visualizza/Apri   Richiedi una copia

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/20.500.14243/330281
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 7
  • ???jsp.display-item.citation.isi??? 6
social impact