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.展开更多
The purpose of this paper is to investigate the convergence of the mixed finite element method for the initial-boundary value problem for the Sobolev equation Ut-div{aut + b1 u} = f based on the Raviart-Thomas space ...The purpose of this paper is to investigate the convergence of the mixed finite element method for the initial-boundary value problem for the Sobolev equation Ut-div{aut + b1 u} = f based on the Raviart-Thomas space Vh × Wh H(div; × L2(). Optimal order estimates are obtained for the approximation of u, ut, the associated velocity p and divp respectively in L(0,T;L2()), L(0,T;L2()), L(0,T;L2()2), and L2(0, T; L2()). Quasi-optimal order estimates are obtained for the approximations of u, ut in L(0, T; L()) and p in L(0,T; L()2).展开更多
This paper deals with a-posteriori error estimates for piecewise linear finite element approximations of parabolic problems in two space dimensions. The analysis extends previous results for elliptic problems to the p...This paper deals with a-posteriori error estimates for piecewise linear finite element approximations of parabolic problems in two space dimensions. The analysis extends previous results for elliptic problems to the parabolic context.展开更多
An enriched goal-oriented error estimation method with extended degrees of freedom is developed to estimate the error in the continuum-based shell extended finite element method. It leads to high quality local error b...An enriched goal-oriented error estimation method with extended degrees of freedom is developed to estimate the error in the continuum-based shell extended finite element method. It leads to high quality local error bounds in three-dimensional fracture mechanics simulation which involves enrichments to solve the singularity in crack tip. This enriched goal-oriented error estimation gives a chance to evaluate this continuum- based shell extended finite element method simulation. With comparisons of reliability to the stress intensity factor calculation in stretching and bending, the accuracy of the continuum-based shell extended finite element method simulation is evaluated, and the reason of error is discussed.展开更多
Asymptotic expansions for the finite element approximation to the eigenvalues of the mult-igroup diffusion equations on Ω(?)R^n(n=2,3)of reactor theory are given,firstly for a piecewiseuniform triangulation,then for ...Asymptotic expansions for the finite element approximation to the eigenvalues of the mult-igroup diffusion equations on Ω(?)R^n(n=2,3)of reactor theory are given,firstly for a piecewiseuniform triangulation,then for nonuniform quadrilateral meshes,and finally for nonuniformhexahedral meshes.The effect of certain classes of numerical integrations is studied.As applicationsof the expansions,several extrapolation formulas and a posteriori error estimates are obtained.展开更多
A nonlinear Galerkin mixed element (NGME) method and a posteriori error exstimator based on the method are established for the stationary Navier-Stokes equations. The existence and error estimates of the NGME solution...A nonlinear Galerkin mixed element (NGME) method and a posteriori error exstimator based on the method are established for the stationary Navier-Stokes equations. The existence and error estimates of the NGME solution are first discussed, and then a posteriori error estimator based on the NGME method is derived.展开更多
Based on the strain formulation of the quasi-conforming finite element, displacement functions are constructed which have definite physical meaning, and a conclusion can be obtained that the coefficients of the consta...Based on the strain formulation of the quasi-conforming finite element, displacement functions are constructed which have definite physical meaning, and a conclusion can be obtained that the coefficients of the constant and the linear strain are uniquely determined, and the quasi-conforming finite element method is convergent to constant strain. There are different methods for constructing the rigid displacement items, and different methods correspond to different order node errors, and this is different from ordinary displacement method finite element.展开更多
An H^1-Galerkin mixed finite element method is discussed for a class of second order SchrSdinger equation. Optimal error estimates of semidiscrete schemes are derived for problems in one space dimension. At the same t...An H^1-Galerkin mixed finite element method is discussed for a class of second order SchrSdinger equation. Optimal error estimates of semidiscrete schemes are derived for problems in one space dimension. At the same time, optimal error estimates are derived for fully discrete schemes. And it is showed that the H1-Galerkin mixed finite element approximations have the same rate of convergence as in the classical mixed finite element methods without requiring the LBB consistency condition.展开更多
Based on an improved orthogonal expansion in an element, a new error expression of n-degree finite element approximation uh to two-point boundary value problem is derived, and then superconvergence of two order for bo...Based on an improved orthogonal expansion in an element, a new error expression of n-degree finite element approximation uh to two-point boundary value problem is derived, and then superconvergence of two order for both function and derivatives are obtained.展开更多
An error expansion for isoparametric bilinear finite element approximation isestablished with some kind of nonuniform quadrilateral mesh constructed in the followingway:We first decompose the considered polygonal doma...An error expansion for isoparametric bilinear finite element approximation isestablished with some kind of nonuniform quadrilateral mesh constructed in the followingway:We first decompose the considered polygonal domain into several fixed convex macro-quadrilaterals and then link up some equi-proportionate points of the opposite edges in eachmacro-quadrilateral and form a quadrilateral mesh which may not be uniform and so isuseful in the adaptive refinement.展开更多
In this paper, we present the a posteriori error estimate of two-grid mixed finite element methods by averaging techniques for semilinear elliptic equations. We first propose the two-grid algorithms to linearize the m...In this paper, we present the a posteriori error estimate of two-grid mixed finite element methods by averaging techniques for semilinear elliptic equations. We first propose the two-grid algorithms to linearize the mixed method equations. Then, the averaging technique is used to construct the a posteriori error estimates of the two-grid mixed finite element method and theoretical analysis are given for the error estimators. Finally, we give some numerical examples to verify the reliability and efficiency of the a posteriori error estimator.展开更多
This paper proves the saturation assumption for the nonconforming Morley finite ele- ment discretization of the biharmonic equation. This asserts that the error of the Morley approximation under uniform refinement is ...This paper proves the saturation assumption for the nonconforming Morley finite ele- ment discretization of the biharmonic equation. This asserts that the error of the Morley approximation under uniform refinement is strictly reduced by a contraction factor smaller than one up to explicit higher-order data approximation terms. The refinement has at least to bisect any edge such as red refinement or 3-bisections on any triangle. This justifies a hierarchical error estimator for the Morley finite element method, which simply compares the discrete solutions of one mesh and its red-refinement. The related adaptive mesh-refining strategy performs optimally in numerical experiments. A remark for Crouzeix-Raviart nonconforming finite element error control is included.展开更多
Based on the concept of the constitutive relation error along with the residuals of both the origin and the dual problems, a goal-oriented error estimation method with extended degrees of freedom is developed. It lead...Based on the concept of the constitutive relation error along with the residuals of both the origin and the dual problems, a goal-oriented error estimation method with extended degrees of freedom is developed. It leads to the high quality locM error bounds in the problem of the direct-solution steady-state dynamic analysis with a frequency-domain finite element, which involves the enrichments with plural variable basis functions. The solution of the steady-state dynamic procedure calculates the harmonic response directly in terms of the physical degrees of freedom in the model, which uses the mass, damping, and stiffness matrices of the system. A three-dimensional finite element example is carried out to illustrate the computational procedures.展开更多
This paper presents a posteriori residual error estimator for the new mixed el-ement scheme for second order elliptic problem on anisotropic meshes. The reliability and efficiency of our estimator are established with...This paper presents a posteriori residual error estimator for the new mixed el-ement scheme for second order elliptic problem on anisotropic meshes. The reliability and efficiency of our estimator are established without any regularity assumption on the mesh.展开更多
Viscoelastic foundation plays a very important role in civil engineering. It can effectively disperse the structural load into the foundation soil and avoid the damage caused by the concentrated load. The model of Eul...Viscoelastic foundation plays a very important role in civil engineering. It can effectively disperse the structural load into the foundation soil and avoid the damage caused by the concentrated load. The model of Euler-Bernoulli beam on viscoelastic Pasternak foundation can be used to analyze the deformation and response of buildings under complex geological conditions. In this paper, we use Hermite finite element method to get the numerical approximation scheme for the vibration equation of viscoelastic Pasternak foundation beam. Convergence and error estimation are rigourously established. We prove that the fully discrete scheme has convergence order O(τ2+h4), where τis time step size and his space step size. Finally, we give four numerical examples to verify the validity of theoretical analysis.展开更多
In this paper, we propose the nonconforming virtual element method (NCVEM) discretization for the pointwise control constraint optimal control problem governed by elliptic equations. Based on the NCVEM approximation o...In this paper, we propose the nonconforming virtual element method (NCVEM) discretization for the pointwise control constraint optimal control problem governed by elliptic equations. Based on the NCVEM approximation of state equation and the variational discretization of control variables, we construct a virtual element discrete scheme. For the state, adjoint state and control variable, we obtain the corresponding prior estimate in H<sup>1</sup> and L<sup>2</sup> norms. Finally, some numerical experiments are carried out to support the theoretical results.展开更多
For h-lype adaptive finite element method, the local mesh refinement introduces irregular nodes and destroys the standard continuity between elements. The reference nodes of the irregular are used to interpolate eleme...For h-lype adaptive finite element method, the local mesh refinement introduces irregular nodes and destroys the standard continuity between elements. The reference nodes of the irregular are used to interpolate element coordinates and displacements.The improved shape functions, of which the conventional shape functions. are a particular case, are presented to guarantee the continuity, No changes but the shape functions are needed when the mcthod is applied in finite element programs.the computational results the features of the method.higher accuracy,simplicity.fewer degrees of freedom and less computation effort.展开更多
In this paper, it is shown that Quasi-Wilson clement possesses a very special property i.e. the consistency error is of order O(h(2)), one order higher than that of Wilson element.
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 λ.展开更多
An accurate and efficient differential quadrature time element method (DQTEM) is proposed for solving ordi- nary differential equations (ODEs), the numerical dissipation and dispersion of DQTEM is much smaller tha...An accurate and efficient differential quadrature time element method (DQTEM) is proposed for solving ordi- nary differential equations (ODEs), the numerical dissipation and dispersion of DQTEM is much smaller than that of the direct integration method of single/multi steps. Two methods of imposing initial conditions are given, which avoids the tediousness when derivative initial conditions are imposed, and the numerical comparisons indicate that the first method, in which the analog equations of initial displacements and velocities are used to directly replace the differential quadra- ture (DQ) analog equations of ODEs at the first and the last sampling points, respectively, is much more accurate than the second method, in which the DQ analog equations of initial conditions are used to directly replace the DQ analog equations of ODEs at the first two sampling points. On the contrary to the conventional step-by-step direct integration schemes, the solutions at all sampling points can be obtained simultaneously by DQTEM, and generally, one differential quadrature time element may be enough for the whole time domain. Extensive numerical comparisons validate the effi- ciency and accuracy of the proposed method.展开更多
基金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.
文摘The purpose of this paper is to investigate the convergence of the mixed finite element method for the initial-boundary value problem for the Sobolev equation Ut-div{aut + b1 u} = f based on the Raviart-Thomas space Vh × Wh H(div; × L2(). Optimal order estimates are obtained for the approximation of u, ut, the associated velocity p and divp respectively in L(0,T;L2()), L(0,T;L2()), L(0,T;L2()2), and L2(0, T; L2()). Quasi-optimal order estimates are obtained for the approximations of u, ut in L(0, T; L()) and p in L(0,T; L()2).
文摘This paper deals with a-posteriori error estimates for piecewise linear finite element approximations of parabolic problems in two space dimensions. The analysis extends previous results for elliptic problems to the parabolic context.
基金Project supported by the National Natural Science Foundation of China(No.10876100)
文摘An enriched goal-oriented error estimation method with extended degrees of freedom is developed to estimate the error in the continuum-based shell extended finite element method. It leads to high quality local error bounds in three-dimensional fracture mechanics simulation which involves enrichments to solve the singularity in crack tip. This enriched goal-oriented error estimation gives a chance to evaluate this continuum- based shell extended finite element method simulation. With comparisons of reliability to the stress intensity factor calculation in stretching and bending, the accuracy of the continuum-based shell extended finite element method simulation is evaluated, and the reason of error is discussed.
文摘Asymptotic expansions for the finite element approximation to the eigenvalues of the mult-igroup diffusion equations on Ω(?)R^n(n=2,3)of reactor theory are given,firstly for a piecewiseuniform triangulation,then for nonuniform quadrilateral meshes,and finally for nonuniformhexahedral meshes.The effect of certain classes of numerical integrations is studied.As applicationsof the expansions,several extrapolation formulas and a posteriori error estimates are obtained.
文摘A nonlinear Galerkin mixed element (NGME) method and a posteriori error exstimator based on the method are established for the stationary Navier-Stokes equations. The existence and error estimates of the NGME solution are first discussed, and then a posteriori error estimator based on the NGME method is derived.
文摘Based on the strain formulation of the quasi-conforming finite element, displacement functions are constructed which have definite physical meaning, and a conclusion can be obtained that the coefficients of the constant and the linear strain are uniquely determined, and the quasi-conforming finite element method is convergent to constant strain. There are different methods for constructing the rigid displacement items, and different methods correspond to different order node errors, and this is different from ordinary displacement method finite element.
基金Supported by the National Natural Science Foundation of China (10601022)Natural Science Foundation of Inner Mongolia Autonomous Region (200607010106)Youth Science Foundation of Inner Mongolia University(ND0702)
文摘An H^1-Galerkin mixed finite element method is discussed for a class of second order SchrSdinger equation. Optimal error estimates of semidiscrete schemes are derived for problems in one space dimension. At the same time, optimal error estimates are derived for fully discrete schemes. And it is showed that the H1-Galerkin mixed finite element approximations have the same rate of convergence as in the classical mixed finite element methods without requiring the LBB consistency condition.
基金This work was supported by The Special Funds for State Major Basic Research Projects(No. G1999032804)The National Natural Science Foundation of China(Grant No.10471038)
文摘Based on an improved orthogonal expansion in an element, a new error expression of n-degree finite element approximation uh to two-point boundary value problem is derived, and then superconvergence of two order for both function and derivatives are obtained.
文摘An error expansion for isoparametric bilinear finite element approximation isestablished with some kind of nonuniform quadrilateral mesh constructed in the followingway:We first decompose the considered polygonal domain into several fixed convex macro-quadrilaterals and then link up some equi-proportionate points of the opposite edges in eachmacro-quadrilateral and form a quadrilateral mesh which may not be uniform and so isuseful in the adaptive refinement.
文摘In this paper, we present the a posteriori error estimate of two-grid mixed finite element methods by averaging techniques for semilinear elliptic equations. We first propose the two-grid algorithms to linearize the mixed method equations. Then, the averaging technique is used to construct the a posteriori error estimates of the two-grid mixed finite element method and theoretical analysis are given for the error estimators. Finally, we give some numerical examples to verify the reliability and efficiency of the a posteriori error estimator.
文摘This paper proves the saturation assumption for the nonconforming Morley finite ele- ment discretization of the biharmonic equation. This asserts that the error of the Morley approximation under uniform refinement is strictly reduced by a contraction factor smaller than one up to explicit higher-order data approximation terms. The refinement has at least to bisect any edge such as red refinement or 3-bisections on any triangle. This justifies a hierarchical error estimator for the Morley finite element method, which simply compares the discrete solutions of one mesh and its red-refinement. The related adaptive mesh-refining strategy performs optimally in numerical experiments. A remark for Crouzeix-Raviart nonconforming finite element error control is included.
基金Project supported by the National Natural Science Foundation of China (No. 10876100)
文摘Based on the concept of the constitutive relation error along with the residuals of both the origin and the dual problems, a goal-oriented error estimation method with extended degrees of freedom is developed. It leads to the high quality locM error bounds in the problem of the direct-solution steady-state dynamic analysis with a frequency-domain finite element, which involves the enrichments with plural variable basis functions. The solution of the steady-state dynamic procedure calculates the harmonic response directly in terms of the physical degrees of freedom in the model, which uses the mass, damping, and stiffness matrices of the system. A three-dimensional finite element example is carried out to illustrate the computational procedures.
文摘This paper presents a posteriori residual error estimator for the new mixed el-ement scheme for second order elliptic problem on anisotropic meshes. The reliability and efficiency of our estimator are established without any regularity assumption on the mesh.
文摘Viscoelastic foundation plays a very important role in civil engineering. It can effectively disperse the structural load into the foundation soil and avoid the damage caused by the concentrated load. The model of Euler-Bernoulli beam on viscoelastic Pasternak foundation can be used to analyze the deformation and response of buildings under complex geological conditions. In this paper, we use Hermite finite element method to get the numerical approximation scheme for the vibration equation of viscoelastic Pasternak foundation beam. Convergence and error estimation are rigourously established. We prove that the fully discrete scheme has convergence order O(τ2+h4), where τis time step size and his space step size. Finally, we give four numerical examples to verify the validity of theoretical analysis.
文摘In this paper, we propose the nonconforming virtual element method (NCVEM) discretization for the pointwise control constraint optimal control problem governed by elliptic equations. Based on the NCVEM approximation of state equation and the variational discretization of control variables, we construct a virtual element discrete scheme. For the state, adjoint state and control variable, we obtain the corresponding prior estimate in H<sup>1</sup> and L<sup>2</sup> norms. Finally, some numerical experiments are carried out to support the theoretical results.
文摘For h-lype adaptive finite element method, the local mesh refinement introduces irregular nodes and destroys the standard continuity between elements. The reference nodes of the irregular are used to interpolate element coordinates and displacements.The improved shape functions, of which the conventional shape functions. are a particular case, are presented to guarantee the continuity, No changes but the shape functions are needed when the mcthod is applied in finite element programs.the computational results the features of the method.higher accuracy,simplicity.fewer degrees of freedom and less computation effort.
文摘In this paper, it is shown that Quasi-Wilson clement possesses a very special property i.e. the consistency error is of order O(h(2)), one order higher than that of Wilson element.
基金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 (11172028,10772014)
文摘An accurate and efficient differential quadrature time element method (DQTEM) is proposed for solving ordi- nary differential equations (ODEs), the numerical dissipation and dispersion of DQTEM is much smaller than that of the direct integration method of single/multi steps. Two methods of imposing initial conditions are given, which avoids the tediousness when derivative initial conditions are imposed, and the numerical comparisons indicate that the first method, in which the analog equations of initial displacements and velocities are used to directly replace the differential quadra- ture (DQ) analog equations of ODEs at the first and the last sampling points, respectively, is much more accurate than the second method, in which the DQ analog equations of initial conditions are used to directly replace the DQ analog equations of ODEs at the first two sampling points. On the contrary to the conventional step-by-step direct integration schemes, the solutions at all sampling points can be obtained simultaneously by DQTEM, and generally, one differential quadrature time element may be enough for the whole time domain. Extensive numerical comparisons validate the effi- ciency and accuracy of the proposed method.