Two new locking-free nonconforming finite elements for the pure displacement planar elasticity problem are presented. Convergence rates of the elements are uniformly optimal with respect to A. The energy norm and L2 n...Two new locking-free nonconforming finite elements for the pure displacement planar elasticity problem are presented. Convergence rates of the elements are uniformly optimal with respect to A. The energy norm and L2 norm errors are proved to be O(h2) and O(h3), respectively. Numerical tests confirm the theoretical analysis.展开更多
We prove that the error estimates of a large class of nonconforming finite elements are dominated by their approximation errors, which means that the well-known Cea's lemma is still valid for these nonconforming f...We prove that the error estimates of a large class of nonconforming finite elements are dominated by their approximation errors, which means that the well-known Cea's lemma is still valid for these nonconforming finite element methods. Furthermore, we derive the error estimates in both energy and L2 norms under the regularity assumption u ∈ H1+s(Ω) with any s > 0. The extensions to other related problems are possible.展开更多
Two nonconforming penalty methods for the two-dimensional stationary Navier-Stokes equations are studied in this paper.These methods are based on the weakly continuous P1 vector fields and the locally divergence-free(...Two nonconforming penalty methods for the two-dimensional stationary Navier-Stokes equations are studied in this paper.These methods are based on the weakly continuous P1 vector fields and the locally divergence-free(LDF)finite elements,which respectively penalize local divergence and are discontinuous across edges.These methods have no penalty factors and avoid solving the saddle-point problems.The existence and uniqueness of the velocity solution are proved,and the optimal error estimates of the energy norms and L^(2)-norms are obtained.Moreover,we propose unified pressure recovery algorithms and prove the optimal error estimates of L^(2)-norm for pressure.We design a unified iterative method for numerical experiments to verify the correctness of the theoretical analysis.展开更多
We consider, in this paper, the trace averaging domain decomposition method for the second order self-adjoint elliptic problems discretized by a class of nonconforming finite elements, which is only continuous at the ...We consider, in this paper, the trace averaging domain decomposition method for the second order self-adjoint elliptic problems discretized by a class of nonconforming finite elements, which is only continuous at the nodes of the quasi-uniform mesh. We show its geometric convergence and present the dependence of the convergence factor on the relaxation factor, the subdomain diameter H and the mesh parameter h. In essence;, this method is equivalent to the simple iterative method for the preconditioned capacitance equation. The preconditioner implied in this iteration is easily invertible and can be applied to preconditioning the capacitance matrix with the condition number no more than O((1 + In H/h)max(1 + H-2, 1 + In H/h)).展开更多
An optimal order of the multigrid method is given in energy-norm for the nonconforming finite element for solving the biharmonic equation, by using the nodal interpolation operator as the transfer operator between grids.
In this paper, a class of rectangular finite elements for 2m-th-oder elliptic boundary value problems in n-dimension (m, n ≥1) is proposed in a canonical fashion, which includes the (2m - 1)-th Hermite interpolat...In this paper, a class of rectangular finite elements for 2m-th-oder elliptic boundary value problems in n-dimension (m, n ≥1) is proposed in a canonical fashion, which includes the (2m - 1)-th Hermite interpolation element (n = 1), the n-linear finite element (m = 1) and the Adini element (m = 2). A nonconforming triangular finite element for the plate bending problem, with convergent order (O(h2), is also proposed.展开更多
In this paper, we study the explicit expressions of the constants in the error estimates of the lowest order mixed and nonconforming finite element methods. We start with an explicit relation between the error constan...In this paper, we study the explicit expressions of the constants in the error estimates of the lowest order mixed and nonconforming finite element methods. We start with an explicit relation between the error constant of the lowest order Raviart-Thomas interpolation error and the geometric characters of the triangle. This gives an explicit error constant of the lowest order mixed finite element method. Furthermore, similar results can be ex- tended to the nonconforming P1 scheme based on its close connection with the lowest order Raviart-Thomas method. Meanwhile, such explicit a priori error estimates can be used as computable error bounds, which are also consistent with the maximal angle condition for the optimal error estimates of mixed and nonconforming finite element methods.展开更多
Under two hypotheses of nonconforming finite elements of fourth order elliptic problems,we present a side-patchwise projection based error analysis method(SPP-BEAM for short).Such a method is able to avoid both the re...Under two hypotheses of nonconforming finite elements of fourth order elliptic problems,we present a side-patchwise projection based error analysis method(SPP-BEAM for short).Such a method is able to avoid both the regularity condition of exact solutions in the classical error analysis method and the complicated bubble function technique in the recent medius error analysis method.In addition,it is universal enough to admit generalizations.Then,we propose a sufficient condition for these hypotheses by imposing a set of in some sense necessary degrees of freedom of the shape function spaces.As an application,we use the theory to design a P3 second order triangular H2 non-conforming element by enriching two P4 bubble functions and,another P4 second order triangular H2 nonconforming finite element,and a P3 second order tetrahedral H2 non-conforming element by enriching eight P4 bubble functions,adding some more degrees of freedom.展开更多
This paper is a generalization of some recent results concerned with the lower bound property of eigenvalues produced by both the enriched rotated Q_(1) and Crouzeix-Raviart elements of the Stokes eigenvalue problem.T...This paper is a generalization of some recent results concerned with the lower bound property of eigenvalues produced by both the enriched rotated Q_(1) and Crouzeix-Raviart elements of the Stokes eigenvalue problem.The main ingredient are a novel and sharp L^(2) error estimate of discrete eigenfunctions,and a new error analysis of nonconforming finite element methods.展开更多
The main aim of this article is to study the approximation of a locking-free anisotropic nonconforming finite element for the pure displacement boundary value problem of planar linear elasticity. The optimal error est...The main aim of this article is to study the approximation of a locking-free anisotropic nonconforming finite element for the pure displacement boundary value problem of planar linear elasticity. The optimal error estimates are obtained by using some novel approaches and techniques. The method proposed in this article is robust in the sense that the convergence estimates in the energy and L^2-norms are independent-of the Lame parameter λ.展开更多
A lumped mass approximation scheme of a low order Crouzeix-Raviart type noncon- forming triangular finite element is proposed to a kind of nonlinear parabolic integro-differential equations. The L2 error estimate is d...A lumped mass approximation scheme of a low order Crouzeix-Raviart type noncon- forming triangular finite element is proposed to a kind of nonlinear parabolic integro-differential equations. The L2 error estimate is derived on anisotropic meshes without referring to the traditional nonclassical elliptic projection.展开更多
A modified penalty scheme is discussed for solving the Stokes problem with the Crouzeix-Raviart type nonconforming linear triangular finite element. By the L^2 projection method, the superconvergence results for the v...A modified penalty scheme is discussed for solving the Stokes problem with the Crouzeix-Raviart type nonconforming linear triangular finite element. By the L^2 projection method, the superconvergence results for the velocity and pressure are obtained with a penalty parameter larger than that of the classical penalty scheme. The numerical experiments are carried out to confirm the theoretical results.展开更多
This paper deals with a new nonconforming anisotropic rectangular finite element approximation for the planar elasticity problem with pure displacement boundary condition. By use of the special properties of this elem...This paper deals with a new nonconforming anisotropic rectangular finite element approximation for the planar elasticity problem with pure displacement boundary condition. By use of the special properties of this element, and by introducing the complementary space and a series of novel techniques, the optimal error estimates of the energy norm and the L^2-norm are obtained. The restrictions of regularity assumption and quasi-uniform assumption or the inverse assumption on the meshes required in the conventional finite element methods analysis are to be got rid of and the applicable scope of the nonconforming finite elements is extended.展开更多
This paper studies a low order mixed finite element method (FEM) for nonstationary incompressible Navier-Stokes equations. The velocity and pressure are approximated by the nonconforming constrained Q1^4ot element a...This paper studies a low order mixed finite element method (FEM) for nonstationary incompressible Navier-Stokes equations. The velocity and pressure are approximated by the nonconforming constrained Q1^4ot element and the piecewise constant, respectively. The superconvergent error estimates of the velocity in the broken H^1-norm and the pressure in the L^2-norm are obtained respectively when the exact solutions are reasonably smooth. A numerical experiment is carried out to confirm the theoretical results.展开更多
For transient Naiver-Stokes problems, a stabilized nonconforming finite element method is presented, focusing on two pairs inf-sup unstable finite element spaces, i.e., pNC/pNC triangular and pNQ/pNQ quadrilateral fin...For transient Naiver-Stokes problems, a stabilized nonconforming finite element method is presented, focusing on two pairs inf-sup unstable finite element spaces, i.e., pNC/pNC triangular and pNQ/pNQ quadrilateral finite element spaces. The semi- and full-discrete schemes of the stabilized method are studied based on the pressure projection and a variational multi-scale method. It has some attractive features: avoiding higher-order derivatives and edge-based data structures, adding a discrete velocity term only on the fine scale, being effective for high Reynolds number fluid flows, and avoiding increased computation cost. For the full-discrete scheme, it has second-order estimations of time and is unconditionally stable. The presented numerical results agree well with the theoretical results.展开更多
A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretizatio...A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretization. Crouzeix-Raviart nonconforming finite element approximation, namely, nonconforming (P1)2 - P0 element, is used for the velocity and pressure fields with the streamline diffusion technique to cope with usual instabilities caused by the convection and time terms. Stability and error estimates are derived with suitable norms.展开更多
Based on combination of two variational principles, a nonconforming stabilized finite element method is presented for the Reissner-Mindlin plates. The method is convergent when the finite element space is energy-compa...Based on combination of two variational principles, a nonconforming stabilized finite element method is presented for the Reissner-Mindlin plates. The method is convergent when the finite element space is energy-compatible. Error estimates are derived. In particular, three finite element spaces are applied in the computation. Numerical results show that the method is insensitive to the mesh distortion and has better performence than the MITC4 and DKQ methods. With properly chosen parameters, high accuracy can be obtained at coarse meshes.展开更多
In this paper, the quadratic nonconforming brick element (MSLK element) intro- duced in [10] is used for the 3D Stokes equations. The instability for the mixed element pair MSLK-P1 is analyzed, where the vector-valu...In this paper, the quadratic nonconforming brick element (MSLK element) intro- duced in [10] is used for the 3D Stokes equations. The instability for the mixed element pair MSLK-P1 is analyzed, where the vector-valued MSLK element approximates the velocity and the piecewise P1 element approximates the pressure. As a cure, we adopt the piecewise P1 macroelement to discretize the pressure instead of the standard piecewise P1 element on cuboid meshes. This new pair is stable and the optimal error estimate is achieved. Numerical examples verify our theoretical analysis.展开更多
The superconvergence in the finite element method is a phenomenon in which the fi-nite element approximation converges to the exact solution at a rate higher than the optimal order error estimate. Wang proposed and an...The superconvergence in the finite element method is a phenomenon in which the fi-nite element approximation converges to the exact solution at a rate higher than the optimal order error estimate. Wang proposed and analyzed superconvergence of the conforming finite element method by L2-projections. However, since the conforming finite element method (CFEM) requires a strong continuity, it is not easy to construct such finite elements for the complex partial differential equations. Thus, the nonconforming finite element method (NCFEM) is more appealing computationally due to better stability and flexibility properties compared to CFEM. The objective of this paper is to establish a general superconvergence result for the nonconforming finite element approximations for second-order elliptic problems by L2-projection methods by applying the idea presented in Wang. MATLAB codes are published at https://github.com/annaleeharris/Superconvergence-NCFEM for anyone to use and to study. The results of numerical experiments show great promise for the robustness, reliability, flexibility and accuracy of superconvergence in NCFEM by L2- projections.展开更多
In this article, we study the explicit expressions of the constants in the error estimate of the nonconforming finite element method. We explicitly obtain the approximation error estimate and the consistency error est...In this article, we study the explicit expressions of the constants in the error estimate of the nonconforming finite element method. We explicitly obtain the approximation error estimate and the consistency error estimate for the Wilson's element without the regular assumption, respectively, which implies the final finite element error estimate. Such explicit a priori error estimates can be used as computable error bounds.展开更多
基金Project supported by the National Natural Science Foundation of China (Nos. 10771198 and 11071226)the Foundation of International Science and Technology Cooperation of Henan Province
文摘Two new locking-free nonconforming finite elements for the pure displacement planar elasticity problem are presented. Convergence rates of the elements are uniformly optimal with respect to A. The energy norm and L2 norm errors are proved to be O(h2) and O(h3), respectively. Numerical tests confirm the theoretical analysis.
基金supported by the Special Funds for Major State Basic Research Project(Grant No.2005CB321701)
文摘We prove that the error estimates of a large class of nonconforming finite elements are dominated by their approximation errors, which means that the well-known Cea's lemma is still valid for these nonconforming finite element methods. Furthermore, we derive the error estimates in both energy and L2 norms under the regularity assumption u ∈ H1+s(Ω) with any s > 0. The extensions to other related problems are possible.
基金National Nature Science Foundation of China(No.11971337,No.11801387)。
文摘Two nonconforming penalty methods for the two-dimensional stationary Navier-Stokes equations are studied in this paper.These methods are based on the weakly continuous P1 vector fields and the locally divergence-free(LDF)finite elements,which respectively penalize local divergence and are discontinuous across edges.These methods have no penalty factors and avoid solving the saddle-point problems.The existence and uniqueness of the velocity solution are proved,and the optimal error estimates of the energy norms and L^(2)-norms are obtained.Moreover,we propose unified pressure recovery algorithms and prove the optimal error estimates of L^(2)-norm for pressure.We design a unified iterative method for numerical experiments to verify the correctness of the theoretical analysis.
文摘We consider, in this paper, the trace averaging domain decomposition method for the second order self-adjoint elliptic problems discretized by a class of nonconforming finite elements, which is only continuous at the nodes of the quasi-uniform mesh. We show its geometric convergence and present the dependence of the convergence factor on the relaxation factor, the subdomain diameter H and the mesh parameter h. In essence;, this method is equivalent to the simple iterative method for the preconditioned capacitance equation. The preconditioner implied in this iteration is easily invertible and can be applied to preconditioning the capacitance matrix with the condition number no more than O((1 + In H/h)max(1 + H-2, 1 + In H/h)).
文摘An optimal order of the multigrid method is given in energy-norm for the nonconforming finite element for solving the biharmonic equation, by using the nodal interpolation operator as the transfer operator between grids.
基金The work was supported by the National Natural Science Foundation of China (10871011).
文摘In this paper, a class of rectangular finite elements for 2m-th-oder elliptic boundary value problems in n-dimension (m, n ≥1) is proposed in a canonical fashion, which includes the (2m - 1)-th Hermite interpolation element (n = 1), the n-linear finite element (m = 1) and the Adini element (m = 2). A nonconforming triangular finite element for the plate bending problem, with convergent order (O(h2), is also proposed.
基金supported by the Special Funds for Major State Basic Research Project(No.2005CB321701)
文摘In this paper, we study the explicit expressions of the constants in the error estimates of the lowest order mixed and nonconforming finite element methods. We start with an explicit relation between the error constant of the lowest order Raviart-Thomas interpolation error and the geometric characters of the triangle. This gives an explicit error constant of the lowest order mixed finite element method. Furthermore, similar results can be ex- tended to the nonconforming P1 scheme based on its close connection with the lowest order Raviart-Thomas method. Meanwhile, such explicit a priori error estimates can be used as computable error bounds, which are also consistent with the maximal angle condition for the optimal error estimates of mixed and nonconforming finite element methods.
文摘Under two hypotheses of nonconforming finite elements of fourth order elliptic problems,we present a side-patchwise projection based error analysis method(SPP-BEAM for short).Such a method is able to avoid both the regularity condition of exact solutions in the classical error analysis method and the complicated bubble function technique in the recent medius error analysis method.In addition,it is universal enough to admit generalizations.Then,we propose a sufficient condition for these hypotheses by imposing a set of in some sense necessary degrees of freedom of the shape function spaces.As an application,we use the theory to design a P3 second order triangular H2 non-conforming element by enriching two P4 bubble functions and,another P4 second order triangular H2 nonconforming finite element,and a P3 second order tetrahedral H2 non-conforming element by enriching eight P4 bubble functions,adding some more degrees of freedom.
基金The author would like to thank Prof.Shangyou Zhang for helping the numerical experiments.The author was supported by the NSFC under Grants Nos.11571023 and 11401015.
文摘This paper is a generalization of some recent results concerned with the lower bound property of eigenvalues produced by both the enriched rotated Q_(1) and Crouzeix-Raviart elements of the Stokes eigenvalue problem.The main ingredient are a novel and sharp L^(2) error estimate of discrete eigenfunctions,and a new error analysis of nonconforming finite element methods.
基金The research is supported by NSF of China (10371113 10471133)
文摘The main aim of this article is to study the approximation of a locking-free anisotropic nonconforming finite element for the pure displacement boundary value problem of planar linear elasticity. The optimal error estimates are obtained by using some novel approaches and techniques. The method proposed in this article is robust in the sense that the convergence estimates in the energy and L^2-norms are independent-of the Lame parameter λ.
基金Supported by the National Natural Science Foundation of China (10671184)
文摘A lumped mass approximation scheme of a low order Crouzeix-Raviart type noncon- forming triangular finite element is proposed to a kind of nonlinear parabolic integro-differential equations. The L2 error estimate is derived on anisotropic meshes without referring to the traditional nonclassical elliptic projection.
基金supported by the National Natural Science Foundation of China (Nos. 10971203 and 11271340)the Research Fund for the Doctoral Program of Higher Education of China (No. 20094101110006)
文摘A modified penalty scheme is discussed for solving the Stokes problem with the Crouzeix-Raviart type nonconforming linear triangular finite element. By the L^2 projection method, the superconvergence results for the velocity and pressure are obtained with a penalty parameter larger than that of the classical penalty scheme. The numerical experiments are carried out to confirm the theoretical results.
基金Supported by the National Natural Science Foundation of China(10371113,10671184)
文摘This paper deals with a new nonconforming anisotropic rectangular finite element approximation for the planar elasticity problem with pure displacement boundary condition. By use of the special properties of this element, and by introducing the complementary space and a series of novel techniques, the optimal error estimates of the energy norm and the L^2-norm are obtained. The restrictions of regularity assumption and quasi-uniform assumption or the inverse assumption on the meshes required in the conventional finite element methods analysis are to be got rid of and the applicable scope of the nonconforming finite elements is extended.
基金Project supported by the National Natural Science Foundation of China(No.11271340)
文摘This paper studies a low order mixed finite element method (FEM) for nonstationary incompressible Navier-Stokes equations. The velocity and pressure are approximated by the nonconforming constrained Q1^4ot element and the piecewise constant, respectively. The superconvergent error estimates of the velocity in the broken H^1-norm and the pressure in the L^2-norm are obtained respectively when the exact solutions are reasonably smooth. A numerical experiment is carried out to confirm the theoretical results.
基金supported by the National Natural Science Foundation of China(No.11271273)
文摘For transient Naiver-Stokes problems, a stabilized nonconforming finite element method is presented, focusing on two pairs inf-sup unstable finite element spaces, i.e., pNC/pNC triangular and pNQ/pNQ quadrilateral finite element spaces. The semi- and full-discrete schemes of the stabilized method are studied based on the pressure projection and a variational multi-scale method. It has some attractive features: avoiding higher-order derivatives and edge-based data structures, adding a discrete velocity term only on the fine scale, being effective for high Reynolds number fluid flows, and avoiding increased computation cost. For the full-discrete scheme, it has second-order estimations of time and is unconditionally stable. The presented numerical results agree well with the theoretical results.
基金supported by the National Natural Science Foundation of China(No.10771150)the National Basic Research Program of China(No.2005CB321701)+1 种基金the Program for New Century Excellent Talents in University(No.NCET-07-0584)the Natural Science Foundation of Sichuan Province(No.07ZB087)
文摘A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretization. Crouzeix-Raviart nonconforming finite element approximation, namely, nonconforming (P1)2 - P0 element, is used for the velocity and pressure fields with the streamline diffusion technique to cope with usual instabilities caused by the convection and time terms. Stability and error estimates are derived with suitable norms.
基金supported by the Key Technologies R&D Program of Sichuan Province of China(No. 05GG006-006-2)
文摘Based on combination of two variational principles, a nonconforming stabilized finite element method is presented for the Reissner-Mindlin plates. The method is convergent when the finite element space is energy-compatible. Error estimates are derived. In particular, three finite element spaces are applied in the computation. Numerical results show that the method is insensitive to the mesh distortion and has better performence than the MITC4 and DKQ methods. With properly chosen parameters, high accuracy can be obtained at coarse meshes.
基金Supported by the National Natural Science Foundation of China(11171052,11301053,61328206 and 61272371)the Fundamental Research Funds for the Central Universities
文摘In this paper, the quadratic nonconforming brick element (MSLK element) intro- duced in [10] is used for the 3D Stokes equations. The instability for the mixed element pair MSLK-P1 is analyzed, where the vector-valued MSLK element approximates the velocity and the piecewise P1 element approximates the pressure. As a cure, we adopt the piecewise P1 macroelement to discretize the pressure instead of the standard piecewise P1 element on cuboid meshes. This new pair is stable and the optimal error estimate is achieved. Numerical examples verify our theoretical analysis.
文摘The superconvergence in the finite element method is a phenomenon in which the fi-nite element approximation converges to the exact solution at a rate higher than the optimal order error estimate. Wang proposed and analyzed superconvergence of the conforming finite element method by L2-projections. However, since the conforming finite element method (CFEM) requires a strong continuity, it is not easy to construct such finite elements for the complex partial differential equations. Thus, the nonconforming finite element method (NCFEM) is more appealing computationally due to better stability and flexibility properties compared to CFEM. The objective of this paper is to establish a general superconvergence result for the nonconforming finite element approximations for second-order elliptic problems by L2-projection methods by applying the idea presented in Wang. MATLAB codes are published at https://github.com/annaleeharris/Superconvergence-NCFEM for anyone to use and to study. The results of numerical experiments show great promise for the robustness, reliability, flexibility and accuracy of superconvergence in NCFEM by L2- projections.
基金supported by National Natural Science Foundation of China (11071226 11201122)
文摘In this article, we study the explicit expressions of the constants in the error estimate of the nonconforming finite element method. We explicitly obtain the approximation error estimate and the consistency error estimate for the Wilson's element without the regular assumption, respectively, which implies the final finite element error estimate. Such explicit a priori error estimates can be used as computable error bounds.