Efficient preconditioning of variational time discretization methods for parabolic Partial Differential Equations

Weller, Stephan; Basting, Steffen

ESAIM: Mathematical Modelling and Numerical Analysis , Tome 49 (2015) no. 2, pp. 331-347.

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

A. Aziz and P. Monk, Continuous finite elements in space and time for the heat equation. Math. Comput. 52 (1989) 255–274. | DOI | Zbl

E. Bänsch, P. Morin and R.H. Nochetto, Preconditioning a class of fourth order problems by operator splitting. Numer. Math. 118 (2011) 197–228. | DOI | Zbl

R.D. Falgout, J.E. Jones and U.M. Yang, The design and implementation of hypre, a library of parallel high performance preconditioners. In Numerical solution of partial differential equations on parallel computers. Springer (2006) 267–294. | Zbl

S. Hussain, F. Schieweck and S. Turek, Higher order galerkin time discretizations and fast multigrid solvers for the heat equation. J. Numer. Math. 19 (2011) 41–61. | DOI | Zbl

J.L. Lions and E. Magenes, Non-homogeneous Boundary Value Problems and Applications. Vol. I. Springer, New York (1972). | Zbl

A. Logg, K.-A. Mardal and G. Wells, Automated solution of differential equations by the finite element method: The fenics book. Vol. 84. Springer (2012). | Zbl

P. Lesaint and P.A. Raviart, On a Finite Element Method for Solving the Neutron Transport Equation. Analyse Numérique. University Paris VI, Labo (1974). | Zbl

T. Richter, A. Springer and B. Vexler, Efficient numerical realization of discontinuous galerkin methods for temporal discretization of parabolic problems. Numer. Math. (2012) 1–32. | Zbl

F. Schieweck, A-stable discontinuous Galerkin–Petrov time discretization of higher order. J. Numer. Math. 18 (2010) 25–57. | DOI | Zbl

F. Schieweck and G. Matthies, Higher order variational time discretizations for nonlinear systems of ordinary differential equations. Preprint 23/2011, Otto-von-Guericke-Universität Magdeburg (2011).

D. Schötzau and C. Schwab, Time discretization of parabolic problems by the hp-version of the discontinuous galerkin finite element method. SIAM J. Numer. Anal. 38 (2000) 837–875. | DOI | Zbl

D. Schötzau and C. Schwab, hp-discontinuous galerkin time-stepping for parabolic problems. C. R. Acad. Sci. Ser. I Math. 333 (2001) 1121–1126. | Zbl

V. Thomée, Galerkin Finite Element Methods for Parabolic Problems. Number 1054 in Springer Lect. Notes Math. 2nd edition. Springer (1984). | Zbl

T. Werder, K. Gerdes, D. Schötzau and C. Schwab, hp-discontinuous galerkin time stepping for parabolic problems. Comput. Methods Appl. Mech. Eng. 190 (2001) 6685–6708. | DOI | Zbl