AN ODYSSEY INTO LOCAL REFINEMENT AND MULTILEVEL PRECONDITIONING II: STABILIZING HIERARCHICAL BASIS METHODS

In this article, we examine the wavelet modified (or stabilized) hierarchical basis (WHB) methods of Vassilevski and Wang, and extend their original quasiuniformity-based framework and results to local 2D and 3D red-green refinement procedures. The concept of a stable Riesz basis plays a critical role in the original work on WHB, and in the design of efficient multilevel preconditioners in general. We carefully examine the impact of local mesh refinement on Riesz bases and matrix conditioning. In the analysis of WHB methods, a critical first step is to establish that the BPX preconditioner is optimal for the refinement procedures under consideration. Therefore, the first article in this series was devoted to extending the results of Dahmen and Kunoth on the optimality of BPX for 2D local red-green refinement to 3D local red-green refinement procedures introduced by Bornemann-Erdmann-Kornhuber (BEK). These results from the first article, together with the local refinement extension of the WHB analysis framework presented here, allow us to establish optimality of the WHB preconditioner on locally refined meshes in both 2D and 3D. In particular, with the minimal smoothness assumption that the PDE coefficients are in L∞, we establish optimality for the additive WHB preconditioner on locally refined 2D and 3D meshes. An interesting implication of the optimality of WHB preconditioner is the a priori H1-stability of the L2-projection. The existing a posteriori approaches in the literature dictate a reconstruction of the mesh if such conditions cannot be satisfied. The proof techniques employed throughout the paper allow extension of the optimality results, the H1-stability of L2-projection results, and the various supporting results to arbitrary spatial dimension d ≥ 1.

[1]  Panayot S. Vassilevski,et al.  Stabilizing the Hierarchical Basis by Approximate Wavelets, I: Theory , 1997, Numer. Linear Algebra Appl..

[2]  Panayot S. Vassilevski,et al.  Wavelet-Like Methods in the Design of Efficient Multilevel Preconditioners for Elliptic PDEs , 1997 .

[3]  R. Skeel,et al.  Unified Multilevel Adaptive Finite Element Methods for Elliptic Problems , 1988 .

[4]  Joseph E. Pasciak,et al.  On the stability of the L2 projection in H1(Omega) , 2002, Math. Comput..

[5]  Randolph E. Bank,et al.  Hierarchical bases and the finite element method , 1996, Acta Numerica.

[6]  Harry Yserentant,et al.  A basic norm equivalence for the theory of multilevel methods , 1993 .

[7]  Barry Joe,et al.  Quality Local Refinement of Tetrahedral Meshes Based on Bisection , 1995, SIAM J. Sci. Comput..

[8]  Michael J. Holst,et al.  An Odyssey into Local Refinement and Multilevel Preconditioning III: Implementation and Numerical Experiments , 2003, SIAM J. Sci. Comput..

[9]  Randolph E. Bank,et al.  Analysis Of A Two-Level Scheme For Solving Finite Element Equations , 1980 .

[10]  E. G. Sewell,et al.  Automatic generation of triangulations for piecewise polynomial approximation , 1972 .

[11]  Maria Elizabeth G. Ong,et al.  Hierarchical Basis Preconditioners in Three Dimensions , 1997, SIAM J. Sci. Comput..

[12]  F. Bornemann,et al.  Adaptive multivlevel methods in three space dimensions , 1993 .

[13]  Rob Stevenson,et al.  A robust hierarchical basis preconditioner on general meshes , 1997 .

[14]  M. E. Go Ong,et al.  Hierachical basis preconditioners for second order elliptic problems in three dimensions , 1989 .

[15]  Panayot S. Vassilevski,et al.  Stabilizing the Hierarchical Basis by Approximate Wavelets II: Implementation and Numerical Results , 1998, SIAM J. Sci. Comput..

[16]  Panayot S. Vassilevski,et al.  Stabilizing the Hierarchical Basis by Approximate Wavelets, I: Theory , 1997 .

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

[18]  S. Jaffard Wavelet methods for fast resolution of elliptic problems , 1992 .

[19]  Rob P. Stevenson Robustness of the additive and multiplicative frequency decomposition multi-level method , 2005, Computing.

[20]  Carsten Carstensen,et al.  Merging the Bramble-Pasciak-Steinbach and the Crouzeix-Thomée criterion for H1-stability of the L2-projection onto finite element spaces , 2002, Math. Comput..

[21]  Burak Aksoylu,et al.  AN ODYSSEY INTO LOCAL REFINEMENT AND MULTILEVEL PRECONDITIONING I: OPTIMALITY OF THE BPX PRECONDITIONER , 2005 .

[22]  H. Yserentant Old and new convergence proofs for multigrid methods , 1993, Acta Numerica.

[23]  R. Bank,et al.  The hierarchical basis multigrid method , 1988 .

[24]  H. Yserentant On the multi-level splitting of finite element spaces , 1986 .

[25]  GermanyNumerische Mathematik,et al.  Multilevel Preconditioning , 1992 .

[26]  Jinchao Xu,et al.  Some Estimates for a Weighted L 2 Projection , 1991 .

[27]  Jinchao Xu,et al.  Iterative Methods by Space Decomposition and Subspace Correction , 1992, SIAM Rev..