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.展开更多
The present study regards the numerical approximation of solutions of systems of Korteweg-de Vries type,coupled through their nonlinear terms.In our previous work[9],we constructed conservative and dissipative finite ...The present study regards the numerical approximation of solutions of systems of Korteweg-de Vries type,coupled through their nonlinear terms.In our previous work[9],we constructed conservative and dissipative finite element methods for these systems and presented a priori error estimates for the semidiscrete schemes.In this sequel,we present a posteriori error estimates for the semidiscrete and fully discrete approximations introduced in[9].The key tool employed to effect our analysis is the dispersive reconstruction devel-oped by Karakashian and Makridakis[20]for related discontinuous Galerkin methods.We conclude by providing a set of numerical experiments designed to validate the a posteriori theory and explore the effectivity of the resulting error indicators.展开更多
In this paper,we investigate a streamline diffusion finite element approxi- mation scheme for the constrained optimal control problem governed by linear con- vection dominated diffusion equations.We prove the existenc...In this paper,we investigate a streamline diffusion finite element approxi- mation scheme for the constrained optimal control problem governed by linear con- vection dominated diffusion equations.We prove the existence and uniqueness of the discretized scheme.Then a priori and a posteriori error estimates are derived for the state,the co-state and the control.Three numerical examples are presented to illustrate our theoretical results.展开更多
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.展开更多
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.展开更多
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).展开更多
On the basis of composition duality principles, augmented three-field macrohybrid mixed variational problems and finite element schemes are analyzed. The compatibility condition adopted here, for compositional dualiza...On the basis of composition duality principles, augmented three-field macrohybrid mixed variational problems and finite element schemes are analyzed. The compatibility condition adopted here, for compositional dualization, is the coupling operator surjectivity, property that expresses in a general operator sense the Ladysenskaja-Babulka-Brezzi inf-sup condition. Variational macro-hybridization is performed under the assumption of decomposable primal and dual spaces relative to nonoverlapping domain decompositions. Then, through compositional dualization macro-hybrid mixed problems are obtained, with internal boundary dual traces as Lagrange multipliers. Also, "mass" preconditioned aug- mentation of three-field formulations are derived, stabilizing macro-hybrid mixed finite element schemes and rendering possible speed up of rates of convergence. Dual mixed incompressible Darcy flow problems illustrate the theory throughout the paper.展开更多
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.展开更多
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.展开更多
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.展开更多
The combined hybrid finite element method is of an intrinsic mechanism of enhancing coarse-mesh-accuracy of lower order displacement schemes. It was confirmed that the combined hybrid scheme without energy error leads...The combined hybrid finite element method is of an intrinsic mechanism of enhancing coarse-mesh-accuracy of lower order displacement schemes. It was confirmed that the combined hybrid scheme without energy error leads to enhancement of accuracy at coarse meshes, and that the combination parameter plays an important role in the enhancement. As an improvement of conforming bilinear Q(4)-plane element, the combined hybrid method adopted the most convenient quadrilateral displacements-stress mode, i.e.,the mode of compatible isoparametric bilinear displacements and pure constant stresses. By adjusting the combined parameter, the optimized version of the combined hybrid element was obtained and numerical tests indicated that this parameter-adjusted version behaves much better than Q(4)-element and is of high accuracy at coarse meshes. Due to elimination of stress parameters at the elemental level, this combined hybrid version is of the same computational cost as that of Q(4)-element.展开更多
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 develops a posteriori error estimates of residual type for conforming and mixed finite element approximations of the fourth order Cahn-Hilliard equation ut + △(ε△Au-ε^-1f(u)) = 0. It is shown that ...This paper develops a posteriori error estimates of residual type for conforming and mixed finite element approximations of the fourth order Cahn-Hilliard equation ut + △(ε△Au-ε^-1f(u)) = 0. It is shown that the a posteriori error bounds depends on ε^-1 only in some low polynomial order, instead of exponential order. Using these a posteriori error estimates, we construct at2 adaptive algorithm for computing the solution of the Cahn- Hilliard equation and its sharp interface limit, the Hele-Shaw flow. Numerical experiments are presented to show the robustness and effectiveness of the new error estimators and the proposed adaptive algorithm.展开更多
This paper is concerned with recovery type a posteriori error estimates of fully discrete finite element approximation for general convex parabolic optimal control problems with pointwise control constraints.The time ...This paper is concerned with recovery type a posteriori error estimates of fully discrete finite element approximation for general convex parabolic optimal control problems with pointwise control constraints.The time discretization is based on the backward Euler method.The state and the adjoint state are approximated by piecewise linear functions and the control is approximated by piecewise constant functions.We derive the superconvergence properties of finite element solutions.By using the superconvergence results,we obtain recovery type a posteriori error estimates.Some numerical examples are presented to verify the theoretical results.展开更多
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 λ.展开更多
In this paper, a constrained distributed optimal control problem governed by a first- order elliptic system is considered. Least-squares mixed finite element methods, which are not subject to the Ladyzhenkaya-Babuska-...In this paper, a constrained distributed optimal control problem governed by a first- order elliptic system is considered. Least-squares mixed finite element methods, which are not subject to the Ladyzhenkaya-Babuska-Brezzi consistency condition, are used for solving the elliptic system with two unknown state variables. By adopting the Lagrange multiplier approach, continuous and discrete optimality systems including a primal state equation, an adjoint state equation, and a variational inequality for the optimal control are derived, respectively. Both the discrete state equation and discrete adjoint state equation yield a symmetric and positive definite linear algebraic system. Thus, the popular solvers such as preconditioned conjugate gradient (PCG) and algebraic multi-grid (AMG) can be used for rapid solution. Optimal a priori error estimates are obtained, respectively, for the control function in L2 (Ω)-norm, for the original state and adjoint state in H1 (Ω)-norm, and for the flux state and adjoint flux state in H(div; Ω)-norm. Finally, we use one numerical example to validate the theoretical findings.展开更多
In this paper,we study the mathematical formulation for an optimal control problem governed by a linear parabolic integro-differential equation and present the optimality conditions.We then set up its weak formulation...In this paper,we study the mathematical formulation for an optimal control problem governed by a linear parabolic integro-differential equation and present the optimality conditions.We then set up its weak formulation and the finite element approximation scheme.Based on these we derive the a priori error estimates for its finite element approximation both in H1 and L^(2)norms.Furthermore some numerical tests are presented to verify the theoretical results.展开更多
In this paper,we present an a posteriori error estimates of semilinear quadratic constrained optimal control problems using triangular mixed finite element methods.The state and co-state are approximated by the orde...In this paper,we present an a posteriori error estimates of semilinear quadratic constrained optimal control problems using triangular mixed finite element methods.The state and co-state are approximated by the order k≤1 RaviartThomas mixed finite element spaces and the control is approximated by piecewise constant element.We derive a posteriori error estimates for the coupled state and control approximations.A numerical example is presented in confirmation of the theory.展开更多
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.展开更多
文摘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 work was supported in part by the National Science Foundation under grant DMS-1620288。
文摘The present study regards the numerical approximation of solutions of systems of Korteweg-de Vries type,coupled through their nonlinear terms.In our previous work[9],we constructed conservative and dissipative finite element methods for these systems and presented a priori error estimates for the semidiscrete schemes.In this sequel,we present a posteriori error estimates for the semidiscrete and fully discrete approximations introduced in[9].The key tool employed to effect our analysis is the dispersive reconstruction devel-oped by Karakashian and Makridakis[20]for related discontinuous Galerkin methods.We conclude by providing a set of numerical experiments designed to validate the a posteriori theory and explore the effectivity of the resulting error indicators.
基金supported by the National Basic Research Program under the Grant 2005CB321701the National Natural Science Foundation of China under the Grants 60474027 and 10771211.
文摘In this paper,we investigate a streamline diffusion finite element approxi- mation scheme for the constrained optimal control problem governed by linear con- vection dominated diffusion equations.We prove the existence and uniqueness of the discretized scheme.Then a priori and a posteriori error estimates are derived for the state,the co-state and the control.Three numerical examples are presented to illustrate our theoretical results.
基金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 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.
文摘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).
文摘On the basis of composition duality principles, augmented three-field macrohybrid mixed variational problems and finite element schemes are analyzed. The compatibility condition adopted here, for compositional dualization, is the coupling operator surjectivity, property that expresses in a general operator sense the Ladysenskaja-Babulka-Brezzi inf-sup condition. Variational macro-hybridization is performed under the assumption of decomposable primal and dual spaces relative to nonoverlapping domain decompositions. Then, through compositional dualization macro-hybrid mixed problems are obtained, with internal boundary dual traces as Lagrange multipliers. Also, "mass" preconditioned aug- mentation of three-field formulations are derived, stabilizing macro-hybrid mixed finite element schemes and rendering possible speed up of rates of convergence. Dual mixed incompressible Darcy flow problems illustrate the theory throughout the paper.
基金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.
基金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.
基金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.
文摘The combined hybrid finite element method is of an intrinsic mechanism of enhancing coarse-mesh-accuracy of lower order displacement schemes. It was confirmed that the combined hybrid scheme without energy error leads to enhancement of accuracy at coarse meshes, and that the combination parameter plays an important role in the enhancement. As an improvement of conforming bilinear Q(4)-plane element, the combined hybrid method adopted the most convenient quadrilateral displacements-stress mode, i.e.,the mode of compatible isoparametric bilinear displacements and pure constant stresses. By adjusting the combined parameter, the optimized version of the combined hybrid element was obtained and numerical tests indicated that this parameter-adjusted version behaves much better than Q(4)-element and is of high accuracy at coarse meshes. Due to elimination of stress parameters at the elemental level, this combined hybrid version is of the same computational cost as that of Q(4)-element.
基金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.
基金the NSF grants DMS-0410266 and DMS-0710831the China National Basic Research Program under the grant 2005CB321701+1 种基金the Program for the New Century Outstanding Talents in Universities of Chinathe Natural Science Foundation of Jiangsu Province under the grant BK2006511
文摘This paper develops a posteriori error estimates of residual type for conforming and mixed finite element approximations of the fourth order Cahn-Hilliard equation ut + △(ε△Au-ε^-1f(u)) = 0. It is shown that the a posteriori error bounds depends on ε^-1 only in some low polynomial order, instead of exponential order. Using these a posteriori error estimates, we construct at2 adaptive algorithm for computing the solution of the Cahn- Hilliard equation and its sharp interface limit, the Hele-Shaw flow. Numerical experiments are presented to show the robustness and effectiveness of the new error estimators and the proposed adaptive algorithm.
基金supported by Guangdong Province Universities and Colleges Pearl River Scholar Funded Scheme(2008)National Science Foundation of China(10971074)+1 种基金Specialized Research Fund for the Doctoral Program of Higher Education(20114407110009)Hunan Provinical Innovation Foundation for Postgraduate(lx2009 B120)。
文摘This paper is concerned with recovery type a posteriori error estimates of fully discrete finite element approximation for general convex parabolic optimal control problems with pointwise control constraints.The time discretization is based on the backward Euler method.The state and the adjoint state are approximated by piecewise linear functions and the control is approximated by piecewise constant functions.We derive the superconvergence properties of finite element solutions.By using the superconvergence results,we obtain recovery type a posteriori error estimates.Some numerical examples are presented to verify the theoretical results.
基金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 λ.
基金Acknowledgements. The authors would like to thank the editor and the anonymous referee for their valuable comments and suggestions on an earlier version of this paper. The work of Hongxing Rui (corresponding author) was supported by the NationM Natural Science Founda- tion of China (No. 11171190). The work of Hongfei Fu was supported by the National Natural Science Foundation of China (No. 11201485), the Promotive Research Fund for Excellent Young and Middle-aged Scientists of Shandong Province (No. BS2013NJ001), and the Fundamental Research Funds for the Central Universities (No. 14CX02217A).
文摘In this paper, a constrained distributed optimal control problem governed by a first- order elliptic system is considered. Least-squares mixed finite element methods, which are not subject to the Ladyzhenkaya-Babuska-Brezzi consistency condition, are used for solving the elliptic system with two unknown state variables. By adopting the Lagrange multiplier approach, continuous and discrete optimality systems including a primal state equation, an adjoint state equation, and a variational inequality for the optimal control are derived, respectively. Both the discrete state equation and discrete adjoint state equation yield a symmetric and positive definite linear algebraic system. Thus, the popular solvers such as preconditioned conjugate gradient (PCG) and algebraic multi-grid (AMG) can be used for rapid solution. Optimal a priori error estimates are obtained, respectively, for the control function in L2 (Ω)-norm, for the original state and adjoint state in H1 (Ω)-norm, and for the flux state and adjoint flux state in H(div; Ω)-norm. Finally, we use one numerical example to validate the theoretical findings.
基金W.F.Shen was supported by National Natural Science Foundation of China(Grant:11326226)Nature Science Foundation of Shandong Province(No.ZR2012GM018)D.P.Yang partially was supported by National Natural Science Foundation of China,Grant:11071080.
文摘In this paper,we study the mathematical formulation for an optimal control problem governed by a linear parabolic integro-differential equation and present the optimality conditions.We then set up its weak formulation and the finite element approximation scheme.Based on these we derive the a priori error estimates for its finite element approximation both in H1 and L^(2)norms.Furthermore some numerical tests are presented to verify the theoretical results.
基金supported by Guangdong Provincial‘Zhujiang Scholar Award Project’National Science Foundation of China 10671163+2 种基金the National Basic Research Program under the Grant 2005CB321703Scientific Research Fund of Hunan Provincial Education Department 06A069Hunan Provincial Innovation Foundation for Postgraduate S2008yjscx04。
文摘In this paper,we present an a posteriori error estimates of semilinear quadratic constrained optimal control problems using triangular mixed finite element methods.The state and co-state are approximated by the order k≤1 RaviartThomas mixed finite element spaces and the control is approximated by piecewise constant element.We derive a posteriori error estimates for the coupled state and control approximations.A numerical example is presented in confirmation of the theory.
基金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.