Uniform convergence of Galerkin‐multigrid methods for nonconforming finite elements for second‐order problems with less than full elliptic regularity

In this article we prove uniform convergence estimates for the recently developed Galerkin‐multigrid methods for nonconforming finite elements for second‐order problems with less than full elliptic regularity. These multigrid methods are defined in terms of the “Galerkin approach,” where quadratic forms over coarse grids are constructed using the quadratic form on the finest grid and iterated coarse‐to‐fine intergrid transfer operators. Previously, uniform estimates were obtained for problems with full elliptic regularity, whereas these estimates are derived with less than full elliptic regularity here. Applications to the nonconforming P₁, rotated Q₁, and Wilson finite elements are analyzed. The result applies to the mixed method based on finite elements that are equivalent to these nonconforming elements. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 203–217, 2002; DOI 10.1002/num.10004     

A new rotated nonconforming pyramid element

In this paper, a new nonparametric nonconforming pyramid finite element is introduced. This element takes the five face mean values as the degrees of the freedom and the finite element space is a subspace of P₂. Different from the other nonparametric elements, the basis functions of this new element can be expressed explicitly without solving linear systems locally, which can be achieved by introducing a new reference pyramid. To evaluate the integration, a class of new quadrature formulae with only two/three equally weighted points on pyramid are constructed. We present the error estimation in the presence of quadrature formulae. Numerical results are shown to confirm the optimality of the convergence order for the second order elliptic problems.     

A second order nonconforming rectangular finite element method for approximating Maxwell’s equations

The main objective of this paper is to present a new rectangular nonconforming finite element scheme with the second order convergence behavior for approximation of Maxwell’s equations.Then the corresponding optimal error estimates are derived.The difficulty in construction of this finite element scheme is how to choose a compatible pair of degrees of freedom and shape function space so as to make the consistency error due to the nonconformity of the element being of order O(h 3 ) ,properly one order higher than that of its interpolation error O(h 2 ) in the broken energy norm,where h is the subdivision parameter tending to zero.     

The Morley element for fourth order elliptic equations in any dimensions

In this paper, the well-known nonconforming Morley element for biharmonic equations in two spatial dimensions is extended to any higher dimensions in a canonical fashion. The general n-dimensional Morley element consists of all quadratic polynomials defined on each n-simplex with degrees of freedom given by the integral average of the normal derivative on each (n-1)-subsimplex and the integral average of the function value on each (n-2)-subsimplex. Explicit expressions of nodal basis functions are also obtained for this element on general n-simplicial grids. Convergence analysis is given for this element when it is applied as a nonconforming finite element discretization for the biharmonic equation. The work was supported by the National Natural Science Foundation of China (10571006). This work was supported in part by NSF DMS-0209497 and NSF DMS-0215392 and the Changjiang Professorship through Peking University.     

Ramanujan’s identities and representation of integers by certain binary and quaternary quadratic forms

We revisit old conjectures of Fermat and Euler regarding the representation of integers by binary quadratic form x ²+5y ². Making use of Ramanujan’s ₁ ψ ₁ summation formula, we establish a new Lambert series identity for $\sum_{n,m=-\infty }^{\infty}q^{n^{2}+5m^{2}}We revisit old conjectures of Fermat and Euler regarding the representation of integers by binary quadratic form x ²+5y ². Making use of Ramanujan’s ₁ ψ ₁ summation formula, we establish a new Lambert series identity for ?_n,m=-￥^￥q^n2+5m2\sum_{n,m=-\infty }^{\infty}q^{n^{2}+5m^{2}} . Conjectures of Fermat and Euler are shown to follow easily from this new formula. But we do not stop there. Employing various formulas found in Ramanujan’s notebooks and using a bit of ingenuity, we obtain a collection of new Lambert series for certain infinite products associated with quadratic forms such as x ²+6y ², 2x ²+3y ², x ²+15y ², 3x ²+5y ², x ²+27y ², x ²+5(y ²+z ²+w ²), 5x ²+y ²+z ²+w ². In the process, we find many new multiplicative eta-quotients and determine their coefficients.     

Three‐dimensional quadratic nonconforming brick element

A new nonconforming brick element with quadratic convergence for the energy norm is introduced. The nonconforming element consists of on a cube [?1,1]³, and 14 degree of freedom (DOF). Two types of DOF are introduced. One consists of the value at the eight vertices and six face‐centroids and the other consists of the value at the eight vertices and the integration value of six faces. Error estimates of optimal order are derived in both broken energy and norms for second‐order elliptic problems. If a genuine hexahedron, which is not a parallelepiped, is included in the partition, the proposed element is also convergent, but with a lower order. Copyright © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 158–174, 2014     

Note on the solution of a certain boundary-value problem

A simple algorithm is described for inverting the operatorD _x D _y (D _x andD _y here and subsequently denote partial differentiation with respect tox andy respectively) which occurs in the iterative solution of the equationD _x D _y f (x, y)=g (x, y, f, D _x f, D _x ² f,D _x D _y f, D _y ² f) when boundary values off(x,y) are given along the sides of the rectangle in thexy-plane whose corners are at the points (a,b); (a+(n+1)k,b); (a+(n+1)k,b+(n+1)h); (a,b+(n+1)h).Communication M. R. 43 of the Computation Department of the Mathematical Centre, Amsterdam.     

Guaranteed velocity error control for the pseudostress approximation of the Stokes equations

The pseudostress approximation of the Stokes equations rewrites the stationary Stokes equations with pure (but possibly inhomogeneous) Dirichlet boundary conditions as another (equivalent) mixed scheme based on a stress in H(div) and the velocity in L². Any standard mixed finite element function space can be utilized for this mixed formulation, e.g., the Raviart‐Thomas discretization which is related to the Crouzeix‐Raviart nonconforming finite element scheme in the lowest‐order case. The effective and guaranteed a posteriori error control for this nonconforming velocity‐oriented discretization can be generalized to the error control of some piecewise quadratic velocity approximation that is related to the discrete pseudostress. The analysis allows for local inf‐sup constants which can be chosen in a global partition to improve the estimation. Numerical examples provide strong evidence for an effective and guaranteed error control with very small overestimation factors even for domains with large anisotropy.© 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1411–1432, 2016     

A posteriori error analysis of nonconforming finite volume elements for general second‐order elliptic PDEs

In this article, we study the a posteriori H¹ and L² error estimates for Crouzeix‐Raviart nonconforming finite volume element discretization of general second‐order elliptic problems in ?². The error estimators yield global upper and local lower bounds. Finally, numerical experiments are performed to illustrate the theoretical findings. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

A new highly accurate discretization for three‐dimensional singularly perturbed nonlinear elliptic partial differential equations

We employ a new fourth‐order compact finite difference formula based on arithmetic average discretization to solve the three‐dimensional nonlinear singularly perturbed elliptic partial differential equation ε(u_xx + u_yy + u_zz) = f(x, y, z, u, u_x, u_y, u_z), 0 < x, y, z < 1, subject to appropriate Dirichlet boundary conditions prescribed on the boundary, where ε > 0 is a small parameter. We also describe new fourth‐order methods for the estimates of (?u/?x), (?u/?y), and (?u/?z), which are quite often of interest in many physical problems. In all cases, we require only a single computational cell with 19 grid points. The proposed methods are directly applicable to solve singular problems without any modification. We solve three test problems numerically to validate the proposed derived fourth‐order methods. We compare the advantages and implementation of the proposed methods with the standard central difference approximations in the context of basic iterative methods. Numerical examples are given to verify the fourth‐order convergence rate of the methods. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2006     

A superconvergent nonconforming mixed finite element method for the Navier–Stokes equations

The superconvergence for a nonconforming mixed finite element approximation of the Navier–Stokes equations is analyzed in this article. The velocity field is approximated by the constrained nonconforming rotated Q₁ (CNRQ₁) element, and the pressure is approximated by the piecewise constant functions. Under some regularity assumptions, the superconvergence estimates for both the velocity in broken H¹‐norm and the pressure in L²‐norm are obtained. Some numerical examples are presented to demonstrate our theoretical results. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 646–660, 2016     

Quadratic invariant curves for a planar mapping

A planar mapping was derived from a second order delay differential equation with a piecewise constant argument. Invariant curves for the planar mapping reflects on the dynamics of the differential equation. Results were reported on a planar mapping admitting quadratic invariant curves y=x ²+C, except for the case -3/4≥C≤0. This remaining case is now resolved, and we describe the solutions of the functional equation K(x ²+C)+k(x)=x by iterations of y.     

Convergence and superconvergence analysis of an anisotropic nonconforming finite element methods for semisingularly perturbed reaction–diffusion problems

The numerical approximation by a lower‐order anisotropic nonconforming finite element on appropriately graded meshes are considered for solving semisingular perturbation problems. The quasi‐optimal‐order error estimates are proved in the ε‐weighted H¹‐norm valid uniformly, up to a logarithmic factor, in the singular perturbation parameter. By using the interpolation postprocessing technique, the global superconvergent error estimates in ε‐weighted H¹‐norm are obtained. Numerical experiments are given to demonstrate validity of our theoretical analysis. Copyright © 2007 John Wiley & Sons, Ltd.     

A partially penalized P1/CR immersed finite element method for planar elasticity interface problems

We propose a partially penalized P₁/CR immersed finite element (IFE) method with midpoint values on edges as degrees of freedom for CR elements to solve planar elasticity interface problems. Optimal approximation errors in L² norm and H¹ semi‐norm are obtained for the P₁/CR IFE spaces. Moreover, by adding some stabilization terms on the edges of interface elements, we derive an optimal error estimate for the P₁/CR IFE method. Our method differs from the method with average values on edges as degrees of freedom for P₁/CR elements in Qin et al.'s study, where no approximation theoretical result was presented. Numerical examples confirm our theoretical results.     

A highly accurate numerical solution of a biharmonic equation

The coefficients for a nine–point high–order accuracy discretization scheme for a biharmonic equation ∇ ⁴u = f(x, y) (∇² is the two–dimensional Laplacian operator) are derived. The biharmonic problem is defined on a rectangular domain with two types of boundary conditions: (1) u and ∂²u/∂n² or (2) u and ∂u/part;n (where ∂/part;n is the normal to the boundary derivative) are specified at the boundary. For both considered cases, the truncation error for the suggested scheme is of the sixth-order O(h⁶) on a square mesh (h_x = h_y = h) and of the fourth-order O(h⁴_xh²_xh²_y h⁴_y) on an unequally spaced mesh. The biharmonic equation describes the deflection of loaded plates. The advantage of the suggested scheme is demonstrated for solving problems of the deflection of rectangular plates for cases of different boundary conditions: (1) a simply supported plate and (2) a plate with built-in edges. In order to demonstrate the high–order accuracy of the method, the numerical results are compared with exact solutions. © John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 13: 375–391, 1997     

A new cubic nonconforming finite element on rectangles

A new nonconforming rectangle element with cubic convergence for the energy norm is introduced. The degrees of freedom (DOFs) are defined by the 12 values at the three Gauss points on each of the four edges. Due to the existence of one linear relation among the above DOFs, it turns out the DOFs are 11. The nonconforming element consists of . We count the corresponding dimension for Dirichlet and Neumann boundary value problems of second‐order elliptic problems. We also present the optimal error estimates in both broken energy and norms. Finally, numerical examples match our theoretical results very well. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 691–705, 2015     

An adaptive least squares mixed finite element method for the stress‐displacement formulation of linear elasticity

A least‐squares mixed finite element method for linear elasticity, based on a stress‐displacement formulation, is investigated in terms of computational efficiency. For the stress approximation quadratic Raviart‐Thomas elements are used and these are coupled with the quadratic nonconforming finite element spaces of Fortin and Soulie for approximating the displacement. The local evaluation of the least‐squares functional serves as an a posteriori error estimator to be used in an adaptive refinement algorithm. We present computational results for a benchmark test problem of planar elasticity including nearly incompressible material parameters in order to verify the effectiveness of our adaptive strategy. For comparison, conforming quadratic finite elements are also used for the displacement approximation showing convergence orders similar to the nonconforming case, which are, however, not independent of the Lamé parameters. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Stability of the Quadratic Equation of Pexider Type

We will investigate the stability problem of the quadratic equation (1) and extend the results of Borelli and Forti, Czerwik, and Rassias. By applying this result and an improved theorem of the author, we will also prove the stability of the quadratic functional equation of Pexider type,f ₁ (x +y) + f₂(x -y) =f ₃(x) +f ₄(y), for a large class of functions.     

Remarks on the global regularity and time decay of the 2D MHD equations with partial dissipation

This paper focuses on a system of the two‐dimensional (2D) magnetohydrodynamic (MHD) equations with the partial kinematic dissipation (?_yyu₁,?_xxu₂) and the partial magnetic diffusion (?_yyb₁,?_xxb₂). Based on the basic energy estimates only, we are able to show that this system always possesses a unique global smooth solution when the initial data are sufficiently smooth. Moreover, we obtain optimal large‐time decay rates of both solutions and their higher order derivatives by developing the classic Fourier splitting methods together with the auxiliary decay estimates of the first derivative of solutions and induction technique.