Superconvergence of a nonconforming brick element for the quad-curl problem

This short note shows the superconvergence of an $H(\mathrm{grad}\,\mathrm{curl})$-nonconforming brick element very recently introduced in [17] for the quad-curl problem. The supercloseness is based on proper modifications for both the interpolation and the discrete formulation, leading to an $O(h^2)$ superclose order in the discrete $H(\mathrm{grad}\,\mathrm{curl})$ norm. Moreover, we propose a suitable postprocessing method to ensure the global superconvergence. Numerical results verify our theory

Conservative flux recovery from the Q1 conforming finite element method on quadrilateral grids

Compared with standard Galerkin finite element methods, mixed methods for second-order elliptic problems give readily available flux approximation, but in general at the expense of having to deal with a more complicated discrete system. This is especially true when conforming elements are involved. Hence it is advantageous to consider a direct method when finding fluxes is just a small part of the overall modeling processes. The purpose of this article is to introduce a direct method combining the standard Galerkin Q1 conforming method with a cheap local flux recovery formula. The approximate flux resides in the lowest order Raviart-Thomas space and retains local conservation property at the cluster level. A cluster is made up of at most four quadrilaterals

A Reissner-Mindlin plate formulation using symmetric Hu-Zhang elements via polytopal transformations

In this work we develop new finite element discretisations of the shear-deformable Reissner--Mindlin plate problem based on the Hellinger-Reissner principle of symmetric stresses. Specifically, we use conforming Hu-Zhang elements to discretise the bending moments in the space of symmetric square integrable fields with a square integrable divergence $\boldsymbol{M} \in \mathcal{HZ} \subset H^{\mathrm{sym}}(\mathrm{Div})$. The latter results in highly accurate approximations of the bending moments $\boldsymbol{M}$ and in the rotation field being in the discontinuous Lebesgue space $\boldsymbol{\phi} \in [L]^2$, such that the Kirchhoff-Love constraint can be satisfied for $t \to 0$. In order to preserve optimal convergence rates across all variables for the case $t \to 0$, we present an extension of the formulation using Raviart-Thomas elements for the shear stress $\mathbf{q} \in \mathcal{RT} \subset H(\mathrm{div})$. We prove existence and uniqueness in the continuous setting and rely on exact complexes for inheritance of well-posedness in the discrete setting. This work introduces an efficient construction of the Hu-Zhang base functions on the reference element via the polytopal template methodology and Legendre polynomials, making it applicable to hp-FEM. The base functions on the reference element are then mapped to the physical element using novel polytopal transformations, which are suitable also for curved geometries. The robustness of the formulations and the construction of the Hu-Zhang element are tested for shear-locking, curved geometries and an L-shaped domain with a singularity in the bending moments $\boldsymbol{M}$. Further, we compare the performance of the novel formulations with the primal-, MITC- and recently introduced TDNNS methods.Comment: Additional implementation material in: https://github.com/Askys/NGSolve_HuZhang_Elemen

Flux Recovery from Primal Hybrid Finite Element Methods

A flux recovery technique is introduced and analyzed for the computed solution of the primal hybrid finite element method for second-order elliptic problems. The recovery is carried out over a single element at a time while ensuring the continuity of the flux across the interelement edges and the validity of the discrete conservation law at the element level. Our construction is general enough to cover all degreesof polynomialsand gridsof triangular or quadrilateral type. We illustrate the principle using the Raviart–Thomas spaces, but other well-known related function spaces such as the Brezzi–Douglas–Marini (BDM) or Brezzi–Douglas–Fortin–Marini (BDFM) space can be used as well. An extension of the technique to the nonlinear case is given. Numerical results are presented to confirm the theoretical results