By the standard theory,the stable Qk+1,k−Qk,k+1/Qdck divergence-free element converges with the optimal order of approximation for the Stokes equations,but only order k for the velocity in H1-norm and the pressure in...By the standard theory,the stable Qk+1,k−Qk,k+1/Qdck divergence-free element converges with the optimal order of approximation for the Stokes equations,but only order k for the velocity in H1-norm and the pressure in L2-norm.This is due to one polynomial degree less in y direction for the first component of velocity,which is a Qk+1,k polynomial of x and y.In this manuscript,we will show by supercloseness of the divergence free element that the order of convergence is truly k+1,for both velocity and pressure.For special solutions(if the interpolation is also divergence-free),a two-order supercloseness is shown to exist.Numerical tests are provided confirming the accuracy of the theory.展开更多
The Wilson finite element method is considered to solve a class of two- dimensional second order elliptic boundary value problems. By using of the particular structure of the element and some new techniques, we obtain...The Wilson finite element method is considered to solve a class of two- dimensional second order elliptic boundary value problems. By using of the particular structure of the element and some new techniques, we obtain the superclose and global superconvergence on anisotropic meshes. Numerical example is also given to confirm our theoretical analysis.展开更多
This paper is concerned with convergence and superconvergence properties of the local discontinuous Galerkin(LDG)method for two-dimensional semilinear second-order elliptic problems of the form−Δu=f(x,y,u)on Cartesia...This paper is concerned with convergence and superconvergence properties of the local discontinuous Galerkin(LDG)method for two-dimensional semilinear second-order elliptic problems of the form−Δu=f(x,y,u)on Cartesian grids.By introducing special GaussRadau projections and using duality arguments,we obtain,under some suitable choice of numerical fuxes,the optimal convergence order in L2-norm of O(h^(p+1))for the LDG solution and its gradient,when tensor product polynomials of degree at most p and grid size h are used.Moreover,we prove that the LDG solutions are superconvergent with an order p+2 toward particular Gauss-Radau projections of the exact solutions.Finally,we show that the error between the gradient of the LDG solution and the gradient of a special Gauss-Radau projection of the exact solution achieves(p+1)-th order superconvergence.Some numerical experiments are performed to illustrate the theoretical results.展开更多
In this paper,a two-grid mixed finite element method(MFEM)of implicit Backward Euler(BE)formula is presented for the fourth order time-dependent singularly perturbed Bi-wave problem for d-wave superconductors by the n...In this paper,a two-grid mixed finite element method(MFEM)of implicit Backward Euler(BE)formula is presented for the fourth order time-dependent singularly perturbed Bi-wave problem for d-wave superconductors by the nonconforming EQ_(1)^(rot) element.In this approach,the original nonlinear system is solved on the coarse mesh through the Newton iteration method,and then the linear system is computed on the fine mesh with Taylor’s expansion.Based on the high accuracy results of the chosen element,the uniform superclose and superconvergent estimates in the broken H^(1)-norm are derived,which are independent of the negative powers of the perturbation parameter appeared in the considered problem.Numerical results illustrate that the computing cost of the proposed two-grid method is much less than that of the conventional Galerkin MFEM without loss of accuracy.展开更多
An anisotropic nonconforming finite element method is presented for a class of nonlinear Sobolev equations. The optimal error estimates and supercloseness are obtained for both semi-discrete and fully-discrete approxi...An anisotropic nonconforming finite element method is presented for a class of nonlinear Sobolev equations. The optimal error estimates and supercloseness are obtained for both semi-discrete and fully-discrete approximate schemes, which are the same as the traditional finite element methods. In addition, the global superconvergence is derived through the postprocessing technique. Numerical experiments are included to illustrate the feasibility of the proposed method.展开更多
Quasi-Wilson nonconforming finite element approximation for a class of nonlinear Sobolev equa- tions is discussed on rectangular meshes. We first prove that this element has two special characters by novel approaches....Quasi-Wilson nonconforming finite element approximation for a class of nonlinear Sobolev equa- tions is discussed on rectangular meshes. We first prove that this element has two special characters by novel approaches. One is that (Vh(U -- Ihu),VhVh)h may be estimated as order O(h2) when u E H3(Ω), where Iuu denotes the bilinear interpolation of u, vh is a polynomial belongs to quasi-Wilson finite element space and △h denotes the piecewise defined gradient operator, h is the mesh size tending to zero. The other is that the consistency error of this element is of order O(h2)/O(h3) in broken Hi-norm, which is one/two order higher than its interpolation error when u ε Ha(Ω)/H4 ((1). Then we derive the optimal order error estimate and su- perclose property via mean-value method and the known high accuracy result of bilinear element. Furthermore, we deduce the global superconvergence through interpolation post processing technique. At last, an extrapola- tion result of order O(h3), two order higher than traditional error estimate, is obtained by constructing a new suitable extrapolation scheme.展开更多
In this paper, nonconforming quasi-Wilson finite element approximation to a class of nonlinear sine-Gordan equations is discussed. Based on the known higher accuracy results of bilinear element and different technique...In this paper, nonconforming quasi-Wilson finite element approximation to a class of nonlinear sine-Gordan equations is discussed. Based on the known higher accuracy results of bilinear element and different techniques from the existing literature, it is proved that the inner product △↓(u - Ih^1u), △↓vh) and the consistency error can be estimated as order O(h^2) in broken H^1 - norm/L^2 - norm when u ∈ H^3(Ω)/H^4(Ω), where Ih^1u is the bilinear interpolation of u, Vh belongs to the quasi-Wilson finite element space. At the same time, the superclose result with order O(h^2) for semi-discrete scheme under generalized rectangular meshes is derived. Furthermore, a fully-discrete scheme is proposed and the corresponding error estimate of order O(h^2 + τ^2) is obtained for the rectangular partition when u ∈ H^4(Ω), which is as same as that of the bilinear element with ADI scheme and one order higher than that of the usual analysis on nonconforming finite elements.展开更多
In this paper, we consider the linear elasticity problem based on the Hellinger-Reissner variational principle. An O(h2) order superclose property for the stress and displacement and a global superconvergence result...In this paper, we consider the linear elasticity problem based on the Hellinger-Reissner variational principle. An O(h2) order superclose property for the stress and displacement and a global superconvergence result of the displacement are established by employing a Clement interpolation, an integral identity and appropriate postprocessing techniques.展开更多
EQrot nonconforming finite element approximation to a class of nonlinear dual phase lagging heat conduction equations is discussed for semi-discrete and fully-discrete schemes. By use of a special property, that is, t...EQrot nonconforming finite element approximation to a class of nonlinear dual phase lagging heat conduction equations is discussed for semi-discrete and fully-discrete schemes. By use of a special property, that is, the consistency error of this element is of order O(h2) one order higher than its interpolation error O(h), the superclose results of order O(h2) in broken Hi-norm are obtained. At the same time, the global superconvergence in broken Hi-norm is deduced by interpolation postprocessing technique. Moreover, the extrapolation result with order O(h4) is derived by constructing a new interpolation postprocessing operator and extrapolation scheme based on the known asymptotic expansion formulas of EQrot element. Finally, optimal error estimate is gained for a proposed fully-discrete scheme by different approaches from the previous literature.展开更多
In this paper, a new splitting positive definite nonconforming mixed finite element method is proposed for pseudo-hyperbolic equations, in which a quasi-Wilson quadrilateral element is used for the flux p, and the bil...In this paper, a new splitting positive definite nonconforming mixed finite element method is proposed for pseudo-hyperbolic equations, in which a quasi-Wilson quadrilateral element is used for the flux p, and the bilinear element is used for u. Superconvergence results in ||·||div,h norm for p and optimal error estimates in L2 norm for u are derived for both semi-discrete and fully discrete schemes under almost uniform meshes.展开更多
This paper is concerned with the superconvergent points of the continuous Galerkin solutions for delay differential equations of pantograph type. We prove the local nodal superconvergence of continuous Galerkin soluti...This paper is concerned with the superconvergent points of the continuous Galerkin solutions for delay differential equations of pantograph type. We prove the local nodal superconvergence of continuous Galerkin solutions under uniform meshes and locate all the superconvergent points based on the supercloseness between the continuous Galerkin solution U and the interpolation Hhu of the exact solution u. The theoretical results are illustrated by numerical examples.展开更多
The lowest order H1-Galerkin mixed finite element method (for short MFEM) is proposed for a class of nonlinear sine-Gordon equations with the simplest bilinear rectangular element and zero order Raviart- Thomas elem...The lowest order H1-Galerkin mixed finite element method (for short MFEM) is proposed for a class of nonlinear sine-Gordon equations with the simplest bilinear rectangular element and zero order Raviart- Thomas element. Base on the interpolation operator instead of the traditional Ritz projection operator which is an indispensable tool in the traditional FEM analysis, together with mean-value technique and high accuracy analysis, the superclose properties of order O(h2)/O(h2 + r2) in Hi-norm and H(div; Ω)-norm axe deduced for the semi-discrete and the fully-discrete schemes, where h, r- denote the mesh size and the time step, respectively, which improve the results in the previous literature.展开更多
In this paper, we first introduce interpolation operator of projection type in three dimen- sions, from which we derive weak estimates for tricubic block finite elements. Then using the estimate for the W 2, 1-seminor...In this paper, we first introduce interpolation operator of projection type in three dimen- sions, from which we derive weak estimates for tricubic block finite elements. Then using the estimate for the W 2, 1-seminorm of the discrete derivative Green's function and the weak estimates, we show that the tricubic block finite element solution uh and the tricubic interpolant of projection type Πh3u have superclose gradient in the pointwise sense of the L∞-norm. Finally, this supercloseness is applied to superconvergence analysis, and the global superconvergence of the finite element approximation is derived.展开更多
In this paper,nonconforming finite element methods(FEMs)are proposed for the constrained optimal control problems(OCPs)governed by the nonsmooth elliptic equations,in which the popular EQr1 ot element is employed to a...In this paper,nonconforming finite element methods(FEMs)are proposed for the constrained optimal control problems(OCPs)governed by the nonsmooth elliptic equations,in which the popular EQr1 ot element is employed to approximate the state and adjoint state,and the piecewise constant element is used to approximate the control.Firstly,the convergence and superconvergence properties for the nonsmooth elliptic equation are obtained by introducing an auxiliary problem.Secondly,the goal-oriented error estimates are obtained for the objective function through establishing the negative norm error estimate.Lastly,the methods are extended to some other well-known nonconforming elements.展开更多
In this paper,we present a two-grid discretization scheme for semilinear parabolic integro-differential equations by H1-Galerkin mixed finite element methods.We use the lowest order Raviart-Thomas mixed finite element...In this paper,we present a two-grid discretization scheme for semilinear parabolic integro-differential equations by H1-Galerkin mixed finite element methods.We use the lowest order Raviart-Thomas mixed finite elements and continuous linear finite element for spatial discretization,and backward Euler scheme for temporal discretization.Firstly,a priori error estimates and some superclose properties are derived.Secondly,a two-grid scheme is presented and its convergence is discussed.In the proposed two-grid scheme,the solution of the nonlinear system on a fine grid is reduced to the solution of the nonlinear system on a much coarser grid and the solution of two symmetric and positive definite linear algebraic equations on the fine grid and the resulting solution still maintains optimal accuracy.Finally,a numerical experiment is implemented to verify theoretical results of the proposed scheme.The theoretical and numerical results show that the two-grid method achieves the same convergence property as the one-grid method with the choice h=H^(2).展开更多
In this paper,the superconvergence properties of the time-dependent Navier-Stokes equations are investigated by a low order nonconforming mixed finite element method(MFEM).In terms of the integral identity technique,t...In this paper,the superconvergence properties of the time-dependent Navier-Stokes equations are investigated by a low order nonconforming mixed finite element method(MFEM).In terms of the integral identity technique,the superclose error estimates for both the velocity in broken H-norm and the pressure in L2-norm are first obtained,which play a key role to bound the numerical solution in Lx-norm.Then the corresponding global superconvergence results are derived through a suitable interpolation postprocessing approach.Finally,some numerical results are provided to demonstrated the theoretical analysis.展开更多
文摘By the standard theory,the stable Qk+1,k−Qk,k+1/Qdck divergence-free element converges with the optimal order of approximation for the Stokes equations,but only order k for the velocity in H1-norm and the pressure in L2-norm.This is due to one polynomial degree less in y direction for the first component of velocity,which is a Qk+1,k polynomial of x and y.In this manuscript,we will show by supercloseness of the divergence free element that the order of convergence is truly k+1,for both velocity and pressure.For special solutions(if the interpolation is also divergence-free),a two-order supercloseness is shown to exist.Numerical tests are provided confirming the accuracy of the theory.
基金Project supported by the National Natural Science Foundation of China (No. 10371113)
文摘The Wilson finite element method is considered to solve a class of two- dimensional second order elliptic boundary value problems. By using of the particular structure of the element and some new techniques, we obtain the superclose and global superconvergence on anisotropic meshes. Numerical example is also given to confirm our theoretical analysis.
基金This research was supported by the NASA Nebraska Space Grant(Federal Grant/Award Number 80NSSC20M0112).
文摘This paper is concerned with convergence and superconvergence properties of the local discontinuous Galerkin(LDG)method for two-dimensional semilinear second-order elliptic problems of the form−Δu=f(x,y,u)on Cartesian grids.By introducing special GaussRadau projections and using duality arguments,we obtain,under some suitable choice of numerical fuxes,the optimal convergence order in L2-norm of O(h^(p+1))for the LDG solution and its gradient,when tensor product polynomials of degree at most p and grid size h are used.Moreover,we prove that the LDG solutions are superconvergent with an order p+2 toward particular Gauss-Radau projections of the exact solutions.Finally,we show that the error between the gradient of the LDG solution and the gradient of a special Gauss-Radau projection of the exact solution achieves(p+1)-th order superconvergence.Some numerical experiments are performed to illustrate the theoretical results.
基金supported by the National Natural Science Foundation of China(Grant Nos.12201640,12071443).
文摘In this paper,a two-grid mixed finite element method(MFEM)of implicit Backward Euler(BE)formula is presented for the fourth order time-dependent singularly perturbed Bi-wave problem for d-wave superconductors by the nonconforming EQ_(1)^(rot) element.In this approach,the original nonlinear system is solved on the coarse mesh through the Newton iteration method,and then the linear system is computed on the fine mesh with Taylor’s expansion.Based on the high accuracy results of the chosen element,the uniform superclose and superconvergent estimates in the broken H^(1)-norm are derived,which are independent of the negative powers of the perturbation parameter appeared in the considered problem.Numerical results illustrate that the computing cost of the proposed two-grid method is much less than that of the conventional Galerkin MFEM without loss of accuracy.
基金supported by the National Natural Science Foundation of China No.10671184
文摘An anisotropic nonconforming finite element method is presented for a class of nonlinear Sobolev equations. The optimal error estimates and supercloseness are obtained for both semi-discrete and fully-discrete approximate schemes, which are the same as the traditional finite element methods. In addition, the global superconvergence is derived through the postprocessing technique. Numerical experiments are included to illustrate the feasibility of the proposed method.
基金Supported by National Natural Science Foundation of China (Grant Nos. 10971203 11101381+4 种基金 11271340)Specialized Research Fund for the Doctoral Program of Higher Education (Grant No. 20094101110006)Tianyuan Mathematics Foundation of the National Natural Science Foundation of China(Grant No. 11026154)the Natural Science Foundation of Henan Province (Grant Nos. 112300410026 122300410266)
文摘Quasi-Wilson nonconforming finite element approximation for a class of nonlinear Sobolev equa- tions is discussed on rectangular meshes. We first prove that this element has two special characters by novel approaches. One is that (Vh(U -- Ihu),VhVh)h may be estimated as order O(h2) when u E H3(Ω), where Iuu denotes the bilinear interpolation of u, vh is a polynomial belongs to quasi-Wilson finite element space and △h denotes the piecewise defined gradient operator, h is the mesh size tending to zero. The other is that the consistency error of this element is of order O(h2)/O(h3) in broken Hi-norm, which is one/two order higher than its interpolation error when u ε Ha(Ω)/H4 ((1). Then we derive the optimal order error estimate and su- perclose property via mean-value method and the known high accuracy result of bilinear element. Furthermore, we deduce the global superconvergence through interpolation post processing technique. At last, an extrapola- tion result of order O(h3), two order higher than traditional error estimate, is obtained by constructing a new suitable extrapolation scheme.
文摘In this paper, nonconforming quasi-Wilson finite element approximation to a class of nonlinear sine-Gordan equations is discussed. Based on the known higher accuracy results of bilinear element and different techniques from the existing literature, it is proved that the inner product △↓(u - Ih^1u), △↓vh) and the consistency error can be estimated as order O(h^2) in broken H^1 - norm/L^2 - norm when u ∈ H^3(Ω)/H^4(Ω), where Ih^1u is the bilinear interpolation of u, Vh belongs to the quasi-Wilson finite element space. At the same time, the superclose result with order O(h^2) for semi-discrete scheme under generalized rectangular meshes is derived. Furthermore, a fully-discrete scheme is proposed and the corresponding error estimate of order O(h^2 + τ^2) is obtained for the rectangular partition when u ∈ H^4(Ω), which is as same as that of the bilinear element with ADI scheme and one order higher than that of the usual analysis on nonconforming finite elements.
文摘In this paper, we consider the linear elasticity problem based on the Hellinger-Reissner variational principle. An O(h2) order superclose property for the stress and displacement and a global superconvergence result of the displacement are established by employing a Clement interpolation, an integral identity and appropriate postprocessing techniques.
基金Supported by the National Natural Science Foundation of China (Nos. 10971203 11101381)+3 种基金Tianyuan Mathe-matics Foundation of National Natural Science Foundation of China (No. 11026154)Natural Science Foundation of Henan Province (No. 112300410026)Natural Science Foundation of the Education Department of Henan Province (Nos. 2011A110020 12A110021)
文摘EQrot nonconforming finite element approximation to a class of nonlinear dual phase lagging heat conduction equations is discussed for semi-discrete and fully-discrete schemes. By use of a special property, that is, the consistency error of this element is of order O(h2) one order higher than its interpolation error O(h), the superclose results of order O(h2) in broken Hi-norm are obtained. At the same time, the global superconvergence in broken Hi-norm is deduced by interpolation postprocessing technique. Moreover, the extrapolation result with order O(h4) is derived by constructing a new interpolation postprocessing operator and extrapolation scheme based on the known asymptotic expansion formulas of EQrot element. Finally, optimal error estimate is gained for a proposed fully-discrete scheme by different approaches from the previous literature.
基金Supported by the National Natural Science Foundation of China(No.10971203,11271340,11101384)the Specialized Research Fund for the Doctoral Program of Higher Education(No.20094101110006)
文摘In this paper, a new splitting positive definite nonconforming mixed finite element method is proposed for pseudo-hyperbolic equations, in which a quasi-Wilson quadrilateral element is used for the flux p, and the bilinear element is used for u. Superconvergence results in ||·||div,h norm for p and optimal error estimates in L2 norm for u are derived for both semi-discrete and fully discrete schemes under almost uniform meshes.
基金Acknowledgments. The second author is supported by NSFC (Nos. 11571027, 91430215), by Beijing Nova Program (No. 2151100003150140) and by the Importation and Development of High-Caliber Talents Project of Beijing Municipal Institutions (No. CIT&TCD201504012). The third author is supported by the Natural Science Foundation of Fujian Province of China (No.2013J05015), by NSFC (No.11301437), and by the Fundamental Research ~nds for the Central Universities (No. 20720150004).
文摘This paper is concerned with the superconvergent points of the continuous Galerkin solutions for delay differential equations of pantograph type. We prove the local nodal superconvergence of continuous Galerkin solutions under uniform meshes and locate all the superconvergent points based on the supercloseness between the continuous Galerkin solution U and the interpolation Hhu of the exact solution u. The theoretical results are illustrated by numerical examples.
基金Supported in part by the National Natural Science Foundation of China under Grant Nos.11671369,11271340the Natural Science Foundation of the Education Department of Henan Province under Grant Nos.14A110009,16A110022
文摘The lowest order H1-Galerkin mixed finite element method (for short MFEM) is proposed for a class of nonlinear sine-Gordon equations with the simplest bilinear rectangular element and zero order Raviart- Thomas element. Base on the interpolation operator instead of the traditional Ritz projection operator which is an indispensable tool in the traditional FEM analysis, together with mean-value technique and high accuracy analysis, the superclose properties of order O(h2)/O(h2 + r2) in Hi-norm and H(div; Ω)-norm axe deduced for the semi-discrete and the fully-discrete schemes, where h, r- denote the mesh size and the time step, respectively, which improve the results in the previous literature.
基金supported by Natural Science Foundation of Ningbo City (Grant No. 2008A610020)National Natural Science Foundation of China (Grant No. 10671065)the Scientific Research Fund of Hunan Provincial Education Department (Grant Nos. 07C576, 03C212)
文摘In this paper, we first introduce interpolation operator of projection type in three dimen- sions, from which we derive weak estimates for tricubic block finite elements. Then using the estimate for the W 2, 1-seminorm of the discrete derivative Green's function and the weak estimates, we show that the tricubic block finite element solution uh and the tricubic interpolant of projection type Πh3u have superclose gradient in the pointwise sense of the L∞-norm. Finally, this supercloseness is applied to superconvergence analysis, and the global superconvergence of the finite element approximation is derived.
基金supported by the National Natural Science Foundation of China(Nos.11501527,11671369).
文摘In this paper,nonconforming finite element methods(FEMs)are proposed for the constrained optimal control problems(OCPs)governed by the nonsmooth elliptic equations,in which the popular EQr1 ot element is employed to approximate the state and adjoint state,and the piecewise constant element is used to approximate the control.Firstly,the convergence and superconvergence properties for the nonsmooth elliptic equation are obtained by introducing an auxiliary problem.Secondly,the goal-oriented error estimates are obtained for the objective function through establishing the negative norm error estimate.Lastly,the methods are extended to some other well-known nonconforming elements.
基金Science and Technology Research Project of Jilin Provincial Department of Education(JJKH20190634KJ)The work of C.M.Liu was supported by the National Natural Science Foundation of China(11901189)+5 种基金the Key Project of Hunan Provincial Education Department(19A191)L.P.Chen was supported by Natural Science Foundation of China(11501473)the Fundamental Research Funds of the Central Universities of China(2682016CX108)The work of Y.Yang was supported by National Natural Science Foundation of China Project(11671342,11771369,11931003)the Project of Scientific Research Fund of Hunan Provincial Science and Technology Department(2018JJ2374,2018WK4006,2019YZ3003)the Key Project of Hunan Provincial Department of Education(17A210).
文摘In this paper,we present a two-grid discretization scheme for semilinear parabolic integro-differential equations by H1-Galerkin mixed finite element methods.We use the lowest order Raviart-Thomas mixed finite elements and continuous linear finite element for spatial discretization,and backward Euler scheme for temporal discretization.Firstly,a priori error estimates and some superclose properties are derived.Secondly,a two-grid scheme is presented and its convergence is discussed.In the proposed two-grid scheme,the solution of the nonlinear system on a fine grid is reduced to the solution of the nonlinear system on a much coarser grid and the solution of two symmetric and positive definite linear algebraic equations on the fine grid and the resulting solution still maintains optimal accuracy.Finally,a numerical experiment is implemented to verify theoretical results of the proposed scheme.The theoretical and numerical results show that the two-grid method achieves the same convergence property as the one-grid method with the choice h=H^(2).
基金This work is supported by National Natural Science Foundation of China(Nos.11671369,11271340).
文摘In this paper,the superconvergence properties of the time-dependent Navier-Stokes equations are investigated by a low order nonconforming mixed finite element method(MFEM).In terms of the integral identity technique,the superclose error estimates for both the velocity in broken H-norm and the pressure in L2-norm are first obtained,which play a key role to bound the numerical solution in Lx-norm.Then the corresponding global superconvergence results are derived through a suitable interpolation postprocessing approach.Finally,some numerical results are provided to demonstrated the theoretical analysis.