Jun Hu

Jun Hu, Shangyou Zhang

This paper presents a family of mixed finite elements on triangular grids for solving the classical Hellinger-Reissner mixed problem of the elasticity equations. In these elements, the matrix-valued stress field is approximated by the full $C^0$-$P_k$ space enriched by $(k-1)$ $H(\d)$ edge bubble functions on each internal edge, while the displacement field by the full discontinuous $P_{k-1}$ vector-valued space, for the polynomial degree $k\ge 3$. The main challenge is to find the correct stress finite element space matching the full $C^{-1}$-$P_{k-1}$ displacement space. The discrete stability analysis for the inf-sup condition does not rely on the usual Fortin operator, which is difficult to construct. It is done by characterizing the divergence of local stress space which covers the $P_{k-1}$ space of displacement orthogonal to the local rigid-motion. The well-posedness condition and the optimal a priori error estimate are proved for this family of finite elements. Numerical tests are presented to confirm the theoretical results.

Long Chen, Jun Hu, Xuehai Huang

A block diagonal preconditioner with the minimal residual method and a block triangular preconditioner with the generalized minimal residual method are developed for Hu-Zhang mixed finite element methods of linear elasticity. They are based on a new stability result of the saddle point system in mesh-dependent norms. The mesh-dependent norm for the stress corresponds to the mass matrix which is easy to invert while the displacement it is spectral equivalent to Schur complement. A fast auxiliary space preconditioner based on the $H^1$ conforming linear element of the linear elasticity problem is then designed for solving the Schur complement. For both diagonal and triangular preconditioners, it is proved that the conditioning numbers of the preconditioned systems are bounded above by a constant independent of both the crucial Lamé constant and the mesh-size. Numerical examples are presented to support theoretical results. As byproducts, a new stabilized low order mixed finite element method is proposed and analyzed and superconvergence results of Hu-Zhang element are obtained.

Long Chen, Jun Hu, Xuehai Huang

In this paper, we design two classes of stabilized mixed finite element methods for linear elasticity on simplicial grids. In the first class of elements, we use $\boldsymbol{H}(\mathbf{div}, Ω; \mathbb{S})$-$P_k$ and $\boldsymbol{L}^2(Ω; \mathbb{R}^n)$-$P_{k-1}$ to approximate the stress and displacement spaces, respectively, for $1\leq k\leq n$, and employ a stabilization technique in terms of the jump of the discrete displacement over the faces of the triangulation under consideration; in the second class of elements, we use $\boldsymbol{H}_0^1(Ω; \mathbb{R}^n)$-$P_{k}$ to approximate the displacement space for $1\leq k\leq n$, and adopt the stabilization technique suggested by Brezzi, Fortin, and Marini. We establish the discrete inf-sup conditions, and consequently present the a priori error analysis for them. The main ingredient for the analysis is two special interpolation operators, which can be constructed using a crucial $\boldsymbol{H}(\mathbf{div})$ bubble function space of polynomials on each element. The feature of these methods is the low number of global degrees of freedom in the lowest order case. We present some numerical results to demonstrate the theoretical estimates.

Snorre H. Christiansen, Jun Hu, Kaibo Hu

We construct 2D and 3D finite element de Rham sequences of arbitrary polynomial degrees with extra smoothness. Some of these elements have nodal degrees of freedom (DoFs) and can be considered as generalisations of scalar Hermite and Lagrange elements. Using the nodal values, the number of global degrees of freedom is reduced compared with the classical Nédélec and Brezzi-Douglas-Marini (BDM) finite elements, and the basis functions are more canonical and easier to construct. Our finite elements for ${H}(\mathrm{div})$ with regularity $r=2$ coincide with the nonstandard elements given by Stenberg (Numer Math 115(1): 131-139, 2010). We show how regularity decreases in the finite element complexes, so that they branch into known complexes. The standard de Rham complexes of Whitney forms and their higher order version can be regarded as the family with the lowest regularity. The construction of the new families is motivated by the finite element systems.%, and we also establish local exact sequences (geometric decomposition) for the new elements.

Long Chen, Jun Hu, Xuehai Huang

A V-cycle multigrid method for the Hellan-Herrmann-Johnson (HHJ) discretization of the Kirchhoff plate bending problems is developed in this paper. It is shown that the contraction number of the V-cycle multigrid HHJ mixed method is bounded away from one uniformly with respect to the mesh size. The uniform convergence is achieved for the V-cycle multigrid method with only one smoothing step and without full elliptic regularity. The key is a stable decomposition of the kernel space which is derived from an exact sequence of the HHJ mixed method, and the strengthened Cauchy Schwarz inequality. Some numerical experiments are provided to confirm the proposed V-cycle multigrid method. The exact sequences of the HHJ mixed method and the corresponding commutative diagram is of some interest independent of the current context.

Long Chen, Jun Hu, Xuehai Huang et al.

A posteriori error estimators for the symmetric mixed finite element methods for linear elasticity problems of Dirichlet and mixed boundary conditions are proposed. Stability and efficiency of the estimators are proved. Finally, we provide numerical examples to verify the theoretical results.

Jun Hu, Rui Ma

In this paper, a new method is proposed to prove the superconvergence of both the Crouzeix-Raviart and Morley elements. The main idea is to fully employ equivalences with the first order Raviart-Thomas element and the first order Hellan-Herrmann-Johnson element, respectively. In this way, some special conformity of discrete stresses is explored and superconvergence of mixed elements can be used to analyze superconvergence of nonconforming elements. Finally, a half order superconvergence by postprocessing is proved for both nonconforming elements.

Wei Cai, Jun Hu, Shangyou Zhang

In this paper, we will use the interior functions of an hierarchical basis for high order $BDM_p$ elements to enforce the divergence-free condition of a magnetic field $B$ approximated by the H(div) $BDM_p$ basis. The resulting constrained finite element method can be used to solve magnetic induction equation in MHD equations. The proposed procedure is based on the fact that the scalar $(p-1)$-th order polynomial space on each element can be decomposed as an orthogonal sum of the subspace defined by the divergence of the interior functions of the $p$-th order $BDM_p$ basis and the constant function. Therefore, the interior functions can be used to remove element-wise all higher order terms except the constant in the divergence error of the finite element solution of $B$-field. The constant terms from each element can be then easily corrected using a first order H(div) basis globally. Numerical results for a 3-D magnetic induction equation show the effectiveness of the proposed method in enforcing divergence-free condition of the magnetic field.

Jun Hu, Guozhu Yu

In this paper, we present a unified analysis of both convergence and optimality of adaptive mixed finite element methods for a class of problems when the finite element spaces and corresponding a posteriori error estimates under consideration satisfy five hypotheses. We prove that these five conditions are sufficient for convergence and optimality of the adaptive algorithms under consideration. The main ingredient for the analysis is a new method to analyze both discrete reliability and quasi-orthogonality. This new method arises from an appropriate and natural choice of the norms for both the discrete displacement and stress spaces, namely, a mesh-dependent discrete $H^1$ norm for the former and a $L^2$ norm for the latter, and a newly defined projection operator from the discrete stress space on the coarser mesh onto the discrete divergence free space on the finer mesh. As applications, we prove these five hypotheses for the Raviart--Thomas and Brezzi--Douglas--Marini elements of the Poisson and Stokes problems in both 2D and 3D.

Jun Hu, Yunqing Huang, Qun Lin

The aim of the paper is to introduce a new systematic method that can produce lower bounds for eigenvalues. The main idea is to use nonconforming finite element methods. The general conclusion herein is that if local approximation properties of nonconforming finite element spaces $V_h$ are better than global continuity properties of $V_h$, corresponding methods will produce lower bounds for eigenvalues. More precisely, under three conditions on continuity and approximation properties of nonconforming finite element spaces we first show abstract error estimates of approximate eigenvalues and eigenfunctions. Subsequently, we propose one more condition and prove that it is sufficient to guarantee nonconforming finite element methods to produce lower bounds for eigenvalues of symmetric elliptic operators. As one application, we show that this condition hold for most nonconforming elements in literature. As another important application, this condition provides a guidance to modify known nonconforming elements in literature and to propose new nonconforming elements. In fact, we enrich locally the Crouzeix-Raviart element such that the new element satisfies the condition; we propose a new nonconforming element for second order elliptic operators and prove that it will yield lower bounds for eigenvalues. Finally, we prove the saturation condition for most nonconforming elements.

Jun Hu, Zhongci Shi

This paper is devoted to the $L^2$ norm error estimate of the Adini element for the biharmonic equation. Surprisingly, a lower bound is established which proves that the $ L^2$ norm convergence rate can not be higher than that in the energy norm. This proves the conjecture of [Lascaux and Lesaint, Some nonconforming finite elements for the plate bending problem, RAIRO Anal. Numer. 9 (1975), pp. 9--53.] that the convergence rates in both $L^2$ and $H^1$ norms can not be higher than that in the energy norm for this element.

Jun Hu, Zhongci Shi, Xueqin Yang

In the present paper, superconvergence of second order, after an appropriate postprocessing, is achieved for both the two and three dimensional first order rectangular Morley elements of biharmonic equations. The analysis is dependent on superconvergence of second order for the consistency error and a corrected canonical interpolation operator, which help to establish supercloseness of second order for the corrected canonical interpolation. Then the final superconvergence follows a standard postprocessing. For first order nonconforming finite element methods of both two and three dimensional fourth order elliptic problems, it is the first time that full superconvergence of second order is obtained without an extra boundary condition imposed on exact solutions. It is also the first time that superconvergence is established for nonconforming finite element methods of three dimensional fourth order elliptic problems. Numerical results are presented to demonstrate the theoretical results.

Jun Hu, Mira Schedensack

In this paper, we propose two low order nonconforming finite element methods (FEMs) for the three-dimensional Stokes flow that generalize the non-conforming FEM of Kouhia and Stenberg (1995, Comput. Methods Appl. Mech. Engrg.). The finite element spaces proposed in this paper consist of two globally continuous components (one piecewise affine and one enriched component) and one component that is continuous at the midpoints of interior faces. We prove that the discrete Korn inequality and a discrete inf-sup condition hold uniformly in the meshsize and also for a non-empty Neumann boundary. Based on these two results, we show the well-posedness of the discrete problem. Two counterexamples prove that there is no direct generalization of the Kouhia-Stenberg FEM to three space dimensions: The finite element space with one non-conforming and two conforming piecewise affine components does not satisfy a discrete inf-sup condition with piecewise constant pressure approximations, while finite element functions with two non-conforming and one conforming component do not satisfy a discrete Korn inequality.

Jun Hu, Limin Ma

In this paper, three high-accuracy methods for eigenvalues of second order elliptic operators are proposed by using the nonconforming Crouzeix-Raviart(CR for short) element and the nonconforming enriched Crouzeix-Raviart(ECR for short) element. They are based on a crucial full one order superconvergence of the first order mixed Raviart-Thomas(RT for short) element. The main ingredient of such a superconvergence analysis is to employ a discrete Helmholtz decomposition of the difference between the canonical interpolation and the finite element solution of the RT element. In particular, it allows for some vital cancellation between terms in one key sum of boundary terms. Consequently, a full one order superconvergence follows from a special relation between the CR element and the RT element, and the equivalence between the ECR element and the RT element for these two nonconforming elements. These superconvergence results improve those in literature from a half order to a full one order for the RT element, the CR element and the ECR element. Based on the aforementioned superconvergence of the RT element, asymptotic expansions of eigenvalues are established and employed to achieve high accuracy extrapolation methods for these two nonconforming elements. In contrast to a classic analysis in literature, the novelty herein is to use not only the canonical interpolations of these nonconforming elements but also that of the RT element to analyze such asymptotic expansions. Based on the superconvergence of these nonconforming elements, asymptotically exact a posteriori error estimators of eigenvalues are constructed and analyzed for them. Finally, two post-processing methods are proposed to improve accuracy of approximate eigenvalues by employing these a posteriori error estimators.Numerical tests are provided to justify and compare the performance of the aforementioned methods.

Jun Hu, Rui Ma

We propose two families of mixed finite elements for solving the classical Hellinger-Reissner mixed problem of the linear elasticity equations in three dimensions. First, a family of conforming mixed triangular prism elements is constructed by product of elements on triangular meshes and elements in one dimension. The well-posedness is established for all elements with $k\geq1$, which are of $k+1$ order convergence for both the stress and displacement. Besides, a family of reduced stress spaces is proposed by dropping the degrees of polynomial functions associated with faces. As a result, the lowest order conforming mixed triangular prism element has 93 plus 33 degrees of freedom on each element. Second, we construct a new family of nonconforming mixed tetrahedral elements. The shape function spaces of our stress spaces are different from those of the elements in literature.

Jun Hu, Pengzhan Jin

We propose a hybrid iterative method based on MIONet for PDEs, which combines the traditional numerical iterative solver and the recent powerful machine learning method of neural operator, and further systematically analyze its theoretical properties, including the convergence condition, the spectral behavior, as well as the convergence rate, in terms of the errors of the discretization and the model inference. We show the theoretical results for the frequently-used smoothers, i.e. Richardson (damped Jacobi) and Gauss-Seidel. We give an upper bound of the convergence rate of the hybrid method w.r.t. the model correction period, which indicates a minimum point to make the hybrid iteration converge fastest. Several numerical examples including the hybrid Richardson (Gauss-Seidel) iteration for the 1-d (2-d) Poisson equation are presented to verify our theoretical results, and also reflect an excellent acceleration effect. As a meshless acceleration method, it is provided with enormous potentials for practice applications.

Jun Hu, Limin Ma

In this paper we propose a penalized Crouzeix-Raviart element method for eigenvalue problems of second order elliptic operators. The key idea is to add a penalty term to tune the local approximation property and the global continuity property of the discrete eigenfunctions. The feature of this method is that by adjusting the penalty parameter, the resulted discrete eigenvalues can be in a state of "chaos", and consequently a large portion of them can be reliable and approximate the exact ones with high accuracy. Furthermore, we design an algorithm to select such a quasi-optimal penalty parameter. Finally, we provide numerical tests to demonstrate the performance of the proposed method.

Jun Hu, Xueqin Yang, Shuo Zhang

This paper is devoted to the convergence analysis of the Adini element scheme for the fourth order problem in arbitrary dimension. We prove that, the Adini element scheme is $\mathcal{O}(h^2)$ order convergent in energy norm provided the exact solution is in $H^4$, and the convergence rate in $L^2$ norm can not be nontrivially higher than $\mathcal{O}(h^2)$ order. Numerical verifications are presented.

Jun Hu

The design of mixed finite element methods in linear elasticity with symmetric stress approximations has been a longstanding open problem until Arnold and Winther designed the first family of mixed finite elements where the discrete stress space is the space of $H(\d,Ω; \mathbb {S})$---$P_{k+1}$ tensors whose divergence is a $P_{k-1}$ polynomial on each triangle for $k\geq 2$.Such a two dimensional family was extended, by Arnold, Awanou and Winther, to a three dimensional family of mixed elements where the discrete stress space is the space of $H(\d,Ω; \mathbb {S})$---$P_{k+2}$ tensors, whose divergence is a $P_{k-1}$ polynomial on each tetrahedron for $k\geq 2$. In this paper, we are able to construct, in a unified fashion, mixed finite element methods with symmetric stress approximations on an arbitrary simplex in $\mathbb{R}^n$ for any space dimension. On the contrary, the discrete stress space here is the space of $H(\d,Ω;\mathbb {S})$---$P_k$ tensors, and the discrete displacement space here is the space of $L^2(Ω; \mathbb{R}^n)$---$P_{k-1}$ vectors for $k\geq n+1$. These finite element spaces are defined with respect to an arbitrary simplicial triangulation of the domain, and can be regarded as extensions to any dimension of those in two and three dimensions by Hu and Zhang.

Jun Hu, Shangyou Zhang

In this paper, we construct, in a unified fashion, lower order finite element subspaces of spaces of symmetric tensors with square-integrable divergence on a domain in any dimension. These subspaces are essentially the symmetric H(div)-Pk (1=<k<=n) tensor spaces, enriched, for each n-1 dimensional simplex, by (n+1)n/2 H(div)-Pn+1 bubble functions when 1=< k<= n-1, and by (n-1)n/2 H(div)-P n+1 bubble functions when k= n. These spaces can be used to approximate the symmetric matrix field in a mixed formulation problem where the other variable is approximated by discontinuous piecewise Pk-1 polynomials. This in particular leads to first order mixed elements on simplicial grids with total degrees of freedom per element $18$ plus $3$ in 2D, 48 plus 6 in 3D. The previous record of the degrees of freedom of first order mixed elements is, 21 plus 3 in 2D, and 156 plus 6 in 3D, on simplicial grids. We also derive, in a unified way and without using any tools like differential forms, a family of auxiliary mixed finite elements in any dimension. One example in this family is the Raviart-Thomas elements in one dimension, the second example is the mixed finite elements for linear elasticity in two dimensions due to Arnold and Winther, the third example is the mixed finite elements for linear elasticity in three dimensions due to Arnold, Awanou and Winther.

Jun Hu, Xueqin Yang

In this paper, we analyze the lower bound property of the discrete eigenvalues by the rectangular Morley elements of the biharmonic operators in both two and three dimensions. The analysis relies on an identity for the errors of eigenvalues. We explore a refined property of the canonical interpolation operators and use it to analyze the key term in this identity. In particular, we show that such a term is of higher order for two dimensions, and is negative and of second order for three dimensions, which causes a main difficulty. To overcome it, we propose a novel decomposition of the first term in the aforementioned identity. Finally, we establish a saturation condition to show that the discrete eigenvalues are smaller than the exact ones. We present some numerical results to demonstrate the theoretical results.