• Convergence analysis of a nested iterative scheme proposed by Bank,Welfert and Yserentant (BWY) ([Numer. Math., 666: 645-666, 1990]) for solving saddle point system is presented. It is shown that this scheme converges under weaker conditions: the contraction rate for solving the $(1,1)$ block matrix is bound by $(\sqrt{5}-1)/2$. Similar convergence result is also obtained for a class of inexact Uzawa method with even weaker contraction bound $\sqrt{2}/2$. Preconditioned generalized minimal residual method using BWY method as a preconditioner is shown to converge with realistic assumptions.
  • Finite element exterior calculus (FEEC) has been developed as a systematical framework for constructing and analyzing stable and accurate numerical method for partial differential equations by employing differential complexes. This paper is devoted to analyze convergence of adaptive mixed finite element methods for Hodge Laplacian equations based on FEEC without considering harmonic forms. More precisely, a residual type posteriori error estimates is obtained by using the Hodge decomposition, the regular decomposition and bounded commuting quasi-interpolants. An additional marking strategy is added to ensure the quasi-orthogonality. Using this quasi-orthogonality, the uniform convergence of adaptive mixed finite element methods is obtained without assuming the initial mesh size is small enough.
  • Due to the indefiniteness and poor spectral properties, the discretized linear algebraic system of the vector Laplacian by mixed finite element methods is hard to solve. A block diagonal preconditioner has been developed and shown to be an effective preconditioner by Arnold, Falk, and Winther [Acta Numerica, 15:1--155, 2006]. The purpose of this paper is to propose alternative and effective block diagonal and block triangular preconditioners for solving this saddle point system. A variable V-cycle multigrid method with the standard point-wise Gauss-Seidel smoother is proved to be a good preconditioner for a discrete vector Laplacian operator. This multigrid solver will be further used to build preconditioners for the saddle point systems of the vector Laplacian and the Maxwell equations with divergent free constraint. The major benefit of our approach is that the point-wise Gauss-Seidel smoother is more algebraic and can be easily implemented as a black-box smoother.
  • Superconvergence and a posteriori error estimators of recovery type are analyzed for the 4-node hybrid stress quadrilateral finite element method proposed by Pian and Sumihara (Int. J. Numer. Meth. Engrg., 1984, 20: 1685-1695) for linear elasticity problems. Uniform superconvergence of order $O(h^{1+\min\{\alpha,1\}})$ with respect to the Lam\'{e} constant $\lambda$ is established for both the recovered gradients of the displacement vector and the stress tensor under a mesh assumption, where $\alpha>0$ is a parameter characterizing the distortion of meshes from parallelograms to quadrilaterals. A posteriori error estimators based on the recovered quantities are shown to be asymptotically exact. Numerical experiments confirm the theoretical results.