%0 Journal Article %1 bochev2009rehabilitation %A Bochev, Pavel B. %A Ridzal, Denis %D 2009 %I Society for Industrial & Applied Mathematics (SIAM) %J SIAM Journal on Numerical Analysis %K 35j25-bvps-2nd-order-elliptic-equations 65f10-iterative-methods-for-linear-systems 65f30-other-matrix-algorithms 65n12-pdes-bvps-stability-and-convergence-of-numerical-methods 65n30-pdes-bvps-finite-elements 65n50-pdes-bvps-mesh-generation-and-refinement 78a30-electro-and-magnetostatics %N 1 %P 487--507 %R 10.1137/070704265 %T Rehabilitation of the Lowest-Order Raviart--Thomas Element on Quadrilateral Grids %U https://epubs.siam.org/doi/10.1137/070704265 %V 47 %X A recent study D. N. Arnold, D. Boffi, and R. S. Falk, SIAM J. Numer. Anal., 42 (2005), pp. 2429–2451 reveals that convergence of finite element methods using $H(div\,,Ømega)$-compatible finite element spaces deteriorates on nonaffine quadrilateral grids. This phenomena is particularly troublesome for the lowest-order Raviart–Thomas elements, because it implies loss of convergence in some norms for finite element solutions of mixed and least-squares methods. In this paper we propose reformulation of finite element methods, based on the natural mimetic divergence operator M. Shashkov, Conservative Finite Difference Methods on General Grids, CRC Press, Boca Raton, FL, 1996, which restores the order of convergence. Reformulations of mixed Galerkin and least-squares methods for the Darcy equation illustrate our approach. We prove that reformulated methods converge optimally with respect to a norm involving the mimetic divergence operator. Furthermore, we prove that standard and reformulated versions of the mixed Galerkin method lead to identical linear systems, but the two versions of the least-squares method are veritably different. The surprising conclusion is that the degradation of convergence in the mixed method on nonaffine quadrilateral grids is superficial, and that the lowest-order Raviart–Thomas elements are safe to use in this method. However, the breakdown in the least-squares method is real, and there one should use our proposed reformulation.

@article{bochev2009rehabilitation, abstract = { A recent study [D. N. Arnold, D. Boffi, and R. S. Falk, SIAM J. Numer. Anal., 42 (2005), pp. 2429–2451] reveals that convergence of finite element methods using $H(\mathrm{div}\,,\Omega)$-compatible finite element spaces deteriorates on nonaffine quadrilateral grids. This phenomena is particularly troublesome for the lowest-order Raviart–Thomas elements, because it implies loss of convergence in some norms for finite element solutions of mixed and least-squares methods. In this paper we propose reformulation of finite element methods, based on the natural mimetic divergence operator [M. Shashkov, Conservative Finite Difference Methods on General Grids, CRC Press, Boca Raton, FL, 1996], which restores the order of convergence. Reformulations of mixed Galerkin and least-squares methods for the Darcy equation illustrate our approach. We prove that reformulated methods converge optimally with respect to a norm involving the mimetic divergence operator. Furthermore, we prove that standard and reformulated versions of the mixed Galerkin method lead to identical linear systems, but the two versions of the least-squares method are veritably different. The surprising conclusion is that the degradation of convergence in the mixed method on nonaffine quadrilateral grids is superficial, and that the lowest-order Raviart–Thomas elements are safe to use in this method. However, the breakdown in the least-squares method is real, and there one should use our proposed reformulation.}, added-at = {2021-12-19T04:45:39.000+0100}, author = {Bochev, Pavel B. and Ridzal, Denis}, biburl = {https://www.bibsonomy.org/bibtex/275cc1a92016606833cf2fa5e6720a126/gdmcbain}, doi = {10.1137/070704265}, interhash = {10039a6a3bc29e0854a3eba1588f2a45}, intrahash = {75cc1a92016606833cf2fa5e6720a126}, journal = {{SIAM} Journal on Numerical Analysis}, keywords = {35j25-bvps-2nd-order-elliptic-equations 65f10-iterative-methods-for-linear-systems 65f30-other-matrix-algorithms 65n12-pdes-bvps-stability-and-convergence-of-numerical-methods 65n30-pdes-bvps-finite-elements 65n50-pdes-bvps-mesh-generation-and-refinement 78a30-electro-and-magnetostatics}, month = jan, number = 1, pages = {487--507}, publisher = {Society for Industrial {\&} Applied Mathematics ({SIAM})}, timestamp = {2021-12-19T04:45:39.000+0100}, title = {Rehabilitation of the Lowest-Order Raviart--Thomas Element on Quadrilateral Grids}, url = {https://epubs.siam.org/doi/10.1137/070704265}, volume = 47, year = 2009 }