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.| 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.


