Convergence results for 3D sparse grid approaches

textabstractThe convergence behaviour is investigated of solution algorithms for the anisotropic Poisson problem on partially ordered, sparse families of regular grids in 3D. In order to study multilevel techniques on sparse families of grids, first we consider the convergence of a two-level algorithm that applies semi-coarsening successively in each of the coordinate directions. This algorithm shows good convergence, but recursive application of the successive semi-coarsening is not sufficiently efficient. Therefore we introduce another algorithm, which uses collective 3D semi-coarsened coarse grid corrections. The convergence behaviour of this collective version is worse, due to the lack of correspondence between the solutions on the different grids. By solving for the trivial solution we demonstrate that a good convergence behaviour of the collective version of the algorithm can be retained when the different solutions are sufficiently coherent. In order to solve also non-trivial problems, we develop a defect correction process. This algorithm makes use of hierarchical smoothing in order to deal with the problems related to the lack of coherence between the solutions on the different grids. Now good convergence rates are obtained also for non-trivial solutions. All convergence results are obtained for two-level processes. The results show convergence rates which are bounded, independent of the discretisation level and of the anisotropy in the problem.

[1]  W. Mulder.  A new multigrid approach to convection problems , 1989 .

[2]  Michael Griebel,et al.  Parallel multigrid methods on sparse grids , 1990, Forschungsberichte, TU Munich.

[3]  Michael Griebel,et al.  A combination technique for the solution of sparse grid problems , 1990, Forschungsberichte, TU Munich.

[4]  Hans-joachim Bungartz,et al.  An adaptive poisson solver using hierarchical bases and sparse grids , 1991, Forschungsberichte, TU Munich.

[5]  W. Mulder.  A high-resolution Euler solver based on multigrid, semi-coarsening, and defective correction , 1992 .

[6]  Hans-Joachim Bungartz,et al.  Dünne Gitter und deren Anwendung bei der adaptiven Lösung der dreidimensionalen Poisson-Gleichung , 1992 .

[7]  P. Oswald,et al.  On additive Schwarz preconditioners for sparse grid discretizations , 1993 .

[8]  Michael Griebel,et al.  Tensor product type subspace splittings and multilevel iterative methods for anisotropic problems , 1995, Adv. Comput. Math..

[9]  C. Liem,et al.  The Splitting Extrapolation Method: A New Technique in Numerical Solution of Multidimensional Problems , 1995 .

[10]  P. Oswald,et al.  Remarks on the Abstract Theory of Additive and Multiplicative Schwarz Algorithms , 1995 .

[11]  Gero Nieeen.  On the stability of sparse grid splittingsGero , 1995 .

[12]  Pieter W. Hemker.  Sparse-grid finite-volume multigrid for 3D-problems , 1995, Adv. Comput. Math..

[13]  Christoph Zenger,et al.  The Solution of Multidimensional Real Helmholtz Equations on Sparse Grids , 1996, SIAM J. Sci. Comput..

[14]  P. Hemker,et al.  Approximation on partially ordered sets of regular grids , 1997 .

[15]  Panayot S. Vassilevski,et al.  On Two Ways of Stabilizing the Hierarchical Basis Multilevel Methods , 1997, SIAM Rev..

[16]  P. Hemker,et al.  Semi-coarsening in three directions for Euler-flow computations in three dimensions , 1997 .

[17]  C. Pflaum.  A multilevel algorithm for the solution of second order elliptic differential equations on sparse grids , 1998 .