Issue |
ESAIM: M2AN
Volume 49, Number 2, March-April 2015
|
|
---|---|---|
Page(s) | 331 - 347 | |
DOI | https://doi.org/10.1051/m2an/2014036 | |
Published online | 05 February 2015 |
Efficient preconditioning of variational time discretization methods for parabolic Partial Differential Equations∗,∗∗
Applied Mathematics III, Dept. of Mathematics,
Cauerstr. 11, 91058
Erlangen,
Germany
weller@math.fau.de
Received: 19 March 2014
Revised: 21 August 2014
This paper is concerned with the design of efficient preconditioners for systems arising from variational time discretization methods for parabolic partial differential equations. We consider the first order discontinuous Galerkin method (dG(1)) and the second order continuous Galerkin Petrov method (cGP(2)). The time-discrete formulation of these methods leads to a coupled 2 × 2 block system whose efficient solution strongly depends on efficient preconditioning strategies. The preconditioner proposed in this paper is based on a Schur complement formulation for the so called essential unknown. By introducing an inexact factorization of this ill-conditioned fourth order operator, we are able to circumvent complex arithmetic and prove uniform bounds for the condition number of the preconditioned system. In addition, the resulting preconditioned operator is symmetric and positive definite, therefore allowing for the usage of efficient Krylov subspace solvers such as the conjugate gradient method. For both the dG(1) and cGP(2) method, we provide optimal choices for the sole parameter of the preconditioner and deduce corresponding upper bounds for the condition number of the resulting preconditioned system. Several numerical experiments including the heat equation and a convection-diffusion example confirm the theoretical findings.
Mathematics Subject Classification: 65M12 / 65M60
Key words: Finite element method / time discretization / discontinuous Galerkin / preconditioning
Parts of this work have been funded by the German Research Foundation (DFG) as part of the Priority Program (SPP) 1506 “Transport processes at fluidic interfaces”.
The authors would like to thank Friedhelm Schieweck for introducing them to the topic of variational time discretization and Eberhard Bänsch for pointing out known efficient preconditioners for fourth order problems. The authors would also like to thank the anonymous reviewers of the manuscript for their helpful remarks and suggestions.
© EDP Sciences, SMAI, 2015
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.