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.展开更多
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.展开更多
Under the hypotheses that the second-order and third-order derivatives of a function are bounded, an estimate of the radius of the convergence ball of Ostrowski-Traub's method is obtained. An error analysis is given ...Under the hypotheses that the second-order and third-order derivatives of a function are bounded, an estimate of the radius of the convergence ball of Ostrowski-Traub's method is obtained. An error analysis is given which matches the convergence order of the method. Finally, two examples are provided to show applications of our theorem.展开更多
Wavelet has rapid development in the current mathematics new areas. It also has a double meaning of theory and application. In signal and image compression, signal analysis, engineering technology has a wide range of ...Wavelet has rapid development in the current mathematics new areas. It also has a double meaning of theory and application. In signal and image compression, signal analysis, engineering technology has a wide range of applications. In this paper, we use wavelet method, for estimating the density function for censoring data. We evaluate the mean integrated squared error, convergence ratio of given estimator. Also, we obtain empirical distribution of given estimator and verify the conclusion by two simulation examples.展开更多
This paper presents derivation of a priori error estimates and convergence rates of finite element processes for boundary value problems (BVPs) described by self adjoint, non-self adjoint, and nonlinear differential o...This paper presents derivation of a priori error estimates and convergence rates of finite element processes for boundary value problems (BVPs) described by self adjoint, non-self adjoint, and nonlinear differential operators. A posteriori error estimates are discussed in context with local approximations in higher order scalar product spaces. A posteriori error computational framework (without the knowledge of theoretical solution) is presented for all BVPs regardless of the method of approximation employed in constructing the integral form. This enables computations of local errors as well as the global errors in the computed finite element solutions. The two most significant and essential aspects of the research presented in this paper that enable all of the features described above are: 1) ensuring variational consistency of the integral form(s) resulting from the methods of approximation for self adjoint, non-self adjoint, and nonlinear differential operators and 2) choosing local approximations for the elements of a discretization in a subspace of a higher order scalar product space that is minimally conforming, hence ensuring desired global differentiability of the approximations over the discretizations. It is shown that when the theoretical solution of a BVP is analytic, the a priori error estimate (in the asymptotic range, discussed in a later section of the paper) is independent of the method of approximation or the nature of the differential operator provided the resulting integral form is variationally consistent. Thus, the finite element processes utilizing integral forms based on different methods of approximation but resulting in VC integral forms result in the same a priori error estimate and convergence rate. It is shown that a variationally consistent (VC) integral form has best approximation property in some norm, conversely an integral form with best approximation property in some norm is variationally consistent. That is best approximation property of the integral form and the VC of the integral form is equivalent, one cannot exist without the other, hence can be used interchangeably. Dimensional model problems consisting of diffusion equation, convection-diffusion equation, and Burgers equation described by self adjoint, non-self adjoint, and nonlinear differential operators are considered to present extensive numerical studies using Galerkin method with weak form (GM/WF) and least squares process (LSP) to determine computed convergence rates of various error norms and present comparisons with the theoretical convergence rates.展开更多
In this paper we give an almost sharp error estimate of Halley’s iteration for the majorizing sequence. Compared with the corresponding results in [6,14], it is far better. Meanwhile,the convergence theorem is establ...In this paper we give an almost sharp error estimate of Halley’s iteration for the majorizing sequence. Compared with the corresponding results in [6,14], it is far better. Meanwhile,the convergence theorem is established .for Halley’s iteration in Banach spaces.展开更多
In this paper we discuss the convergence of a modified Newton’s method presented by A. Ostrowski [1] and J.F. Traub [2], which has quadratic convergence order but reduces one evaluation of the derivative at every two...In this paper we discuss the convergence of a modified Newton’s method presented by A. Ostrowski [1] and J.F. Traub [2], which has quadratic convergence order but reduces one evaluation of the derivative at every two steps compared with Newton’s method. A convergence theorem is established by using a weak condition a≤3-2(2<sup>1/2</sup>) and a sharp error estimate is given about the iterative sequence.展开更多
This work deals with the numerical solution of singular perturbation system of ordinary differential equations with boundary layer. For the numerical solution of this problem fitted finite difference scheme on a unifo...This work deals with the numerical solution of singular perturbation system of ordinary differential equations with boundary layer. For the numerical solution of this problem fitted finite difference scheme on a uniform mesh is constructed and analyzed. The uniform error estimates for the approximate solution are obtained.展开更多
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.展开更多
This article is concerned with the pointwise error estimates for vanishing vis- cosity approximations to scalar convex conservation laws with boundary.By the weighted error function and a bootstrap extrapolation techn...This article is concerned with the pointwise error estimates for vanishing vis- cosity approximations to scalar convex conservation laws with boundary.By the weighted error function and a bootstrap extrapolation technique introduced by Tadmor-Tang,an optimal pointwise convergence rate is derived for the vanishing viscosity approximations to the initial-boundary value problem for scalar convex conservation laws,whose weak entropy solution is piecewise C 2 -smooth with interaction of elementary waves and the ...展开更多
In this paper two kernel density estimators are introduced and investigated. In order to reduce bias, we intuitively subtract an estimated bias term from ordinary kernel density estimator. The second proposed density ...In this paper two kernel density estimators are introduced and investigated. In order to reduce bias, we intuitively subtract an estimated bias term from ordinary kernel density estimator. The second proposed density estimator is a geometric extrapolation of the first bias reduced estimator. Theoretical properties such as bias, variance and mean squared error are investigated for both estimators. To observe their finite sample performance, a Monte Carlo simulation study based on small to moderately large samples is presented.展开更多
In this paper, we give a convergence theorem and error estimates for an iteration method under new Kantorovitch-Ostrowski type condition using the information of higher derivatives at initial points. Compared with the...In this paper, we give a convergence theorem and error estimates for an iteration method under new Kantorovitch-Ostrowski type condition using the information of higher derivatives at initial points. Compared with the corresponding study in [3], the convergence determination is established under one global condition, instead of two, on the function.展开更多
In this paper we computed the data of pore structure with Lagrange interpolation method, analyzed the convergence of the method and deduced the corresponding error estimation formula.
In this paper, we prove the convergence of the nodal expansion method, a new numerical method for partial differential equations and provide the error estimates of approximation solution.
In this paper we deal with the convergence analysis of the finite element method for an elliptic penalized unilateral obstacle optimal control problem where the control and the obstacle coincide.Error estimates are es...In this paper we deal with the convergence analysis of the finite element method for an elliptic penalized unilateral obstacle optimal control problem where the control and the obstacle coincide.Error estimates are established for both state and control variables.We apply a fixed point type iteration method to solve the discretized problem.To corroborate our error estimations and the eficiency of our algorithms,the convergence results and numerical experiments are illustrated by concrete examples.展开更多
The fractional optimal control problem leads to significantly increased computational complexity compared to the corresponding classical integer-order optimal control problem,due to the global properties of fractional...The fractional optimal control problem leads to significantly increased computational complexity compared to the corresponding classical integer-order optimal control problem,due to the global properties of fractional differential operators.In this paper,we focus on an optimal control problem governed by fractional differential equations with an integral constraint on the state variable.By the proposed first-order optimality condition consisting of a Lagrange multiplier,we design a spectral Galerkin discrete scheme with weighted orthogonal Jacobi polynomials to approximate the resulting state and adjoint state equations.Furthermore,a priori error estimates for state,adjoint state and control variables are discussed in details.Illustrative numerical tests are given to demonstrate the validity and applicability of our proposed approximations and theoretical results.展开更多
In this paper, we discuss virtual element method (VEM) approximation of optimal control problem governed by Brinkman equations with control constraints. Based on the polynomial projections and variational discretizati...In this paper, we discuss virtual element method (VEM) approximation of optimal control problem governed by Brinkman equations with control constraints. Based on the polynomial projections and variational discretization of the control variable, we build up the virtual element discrete scheme of the optimal control problem and derive the discrete first order optimality system. A priori error estimates for the state, adjoint state and control variables in L<sup>2</sup> and H<sup>1</sup> norm are derived. The theoretical findings are illustrated by the numerical experiments.展开更多
In this paper, we study a weakly over-penalized interior penalty method for non-self- adjoint and indefinite problems. An optimal a priori error estimate in the energy norm is derived. In addition, we introduce a resi...In this paper, we study a weakly over-penalized interior penalty method for non-self- adjoint and indefinite problems. An optimal a priori error estimate in the energy norm is derived. In addition, we introduce a residual-based a posteriori error estimator, which is proved to be both reliable and efficient in the energy norm. Some numerical testes are presented to validate our 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.
基金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 (10871178)the Natural Science Foundation(Y606154)the Foundation of the Eduction Department of Zhejiang Province of China (Y200804008)
文摘Under the hypotheses that the second-order and third-order derivatives of a function are bounded, an estimate of the radius of the convergence ball of Ostrowski-Traub's method is obtained. An error analysis is given which matches the convergence order of the method. Finally, two examples are provided to show applications of our theorem.
文摘Wavelet has rapid development in the current mathematics new areas. It also has a double meaning of theory and application. In signal and image compression, signal analysis, engineering technology has a wide range of applications. In this paper, we use wavelet method, for estimating the density function for censoring data. We evaluate the mean integrated squared error, convergence ratio of given estimator. Also, we obtain empirical distribution of given estimator and verify the conclusion by two simulation examples.
文摘This paper presents derivation of a priori error estimates and convergence rates of finite element processes for boundary value problems (BVPs) described by self adjoint, non-self adjoint, and nonlinear differential operators. A posteriori error estimates are discussed in context with local approximations in higher order scalar product spaces. A posteriori error computational framework (without the knowledge of theoretical solution) is presented for all BVPs regardless of the method of approximation employed in constructing the integral form. This enables computations of local errors as well as the global errors in the computed finite element solutions. The two most significant and essential aspects of the research presented in this paper that enable all of the features described above are: 1) ensuring variational consistency of the integral form(s) resulting from the methods of approximation for self adjoint, non-self adjoint, and nonlinear differential operators and 2) choosing local approximations for the elements of a discretization in a subspace of a higher order scalar product space that is minimally conforming, hence ensuring desired global differentiability of the approximations over the discretizations. It is shown that when the theoretical solution of a BVP is analytic, the a priori error estimate (in the asymptotic range, discussed in a later section of the paper) is independent of the method of approximation or the nature of the differential operator provided the resulting integral form is variationally consistent. Thus, the finite element processes utilizing integral forms based on different methods of approximation but resulting in VC integral forms result in the same a priori error estimate and convergence rate. It is shown that a variationally consistent (VC) integral form has best approximation property in some norm, conversely an integral form with best approximation property in some norm is variationally consistent. That is best approximation property of the integral form and the VC of the integral form is equivalent, one cannot exist without the other, hence can be used interchangeably. Dimensional model problems consisting of diffusion equation, convection-diffusion equation, and Burgers equation described by self adjoint, non-self adjoint, and nonlinear differential operators are considered to present extensive numerical studies using Galerkin method with weak form (GM/WF) and least squares process (LSP) to determine computed convergence rates of various error norms and present comparisons with the theoretical convergence rates.
基金Jointly supported by China Major Key Project for Basic Researcher and Provincial Natrual Science Foundation.
文摘In this paper we give an almost sharp error estimate of Halley’s iteration for the majorizing sequence. Compared with the corresponding results in [6,14], it is far better. Meanwhile,the convergence theorem is established .for Halley’s iteration in Banach spaces.
基金Jointly supported by China Major Key Project for Basic Researcher and Provincial Natrural Science Foundation.
文摘In this paper we discuss the convergence of a modified Newton’s method presented by A. Ostrowski [1] and J.F. Traub [2], which has quadratic convergence order but reduces one evaluation of the derivative at every two steps compared with Newton’s method. A convergence theorem is established by using a weak condition a≤3-2(2<sup>1/2</sup>) and a sharp error estimate is given about the iterative sequence.
文摘This work deals with the numerical solution of singular perturbation system of ordinary differential equations with boundary layer. For the numerical solution of this problem fitted finite difference scheme on a uniform mesh is constructed and analyzed. The uniform error estimates for the approximate solution are obtained.
基金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 NSF China#10571075NSF-Guangdong China#04010473+1 种基金The research of the second author was supported by Jinan University Foundation#51204033the Scientific Research Foundation for the Returned Overseas Chinese Scholars,State education Ministry#2005-383
文摘This article is concerned with the pointwise error estimates for vanishing vis- cosity approximations to scalar convex conservation laws with boundary.By the weighted error function and a bootstrap extrapolation technique introduced by Tadmor-Tang,an optimal pointwise convergence rate is derived for the vanishing viscosity approximations to the initial-boundary value problem for scalar convex conservation laws,whose weak entropy solution is piecewise C 2 -smooth with interaction of elementary waves and the ...
文摘In this paper two kernel density estimators are introduced and investigated. In order to reduce bias, we intuitively subtract an estimated bias term from ordinary kernel density estimator. The second proposed density estimator is a geometric extrapolation of the first bias reduced estimator. Theoretical properties such as bias, variance and mean squared error are investigated for both estimators. To observe their finite sample performance, a Monte Carlo simulation study based on small to moderately large samples is presented.
文摘In this paper, we give a convergence theorem and error estimates for an iteration method under new Kantorovitch-Ostrowski type condition using the information of higher derivatives at initial points. Compared with the corresponding study in [3], the convergence determination is established under one global condition, instead of two, on the function.
文摘In this paper we computed the data of pore structure with Lagrange interpolation method, analyzed the convergence of the method and deduced the corresponding error estimation formula.
基金This project is supported by the National Science Foundation of China
文摘In this paper, we prove the convergence of the nodal expansion method, a new numerical method for partial differential equations and provide the error estimates of approximation solution.
文摘In this paper we deal with the convergence analysis of the finite element method for an elliptic penalized unilateral obstacle optimal control problem where the control and the obstacle coincide.Error estimates are established for both state and control variables.We apply a fixed point type iteration method to solve the discretized problem.To corroborate our error estimations and the eficiency of our algorithms,the convergence results and numerical experiments are illustrated by concrete examples.
基金This work was partly supported by National Natural Science Foundation of China(Grant Nos.:12101283,12271233 and 12171287)Natural Science Foundation of Shandong Province(Grant Nos.:ZR2019YQ05,2019KJI003,and ZR2016JL004).
文摘The fractional optimal control problem leads to significantly increased computational complexity compared to the corresponding classical integer-order optimal control problem,due to the global properties of fractional differential operators.In this paper,we focus on an optimal control problem governed by fractional differential equations with an integral constraint on the state variable.By the proposed first-order optimality condition consisting of a Lagrange multiplier,we design a spectral Galerkin discrete scheme with weighted orthogonal Jacobi polynomials to approximate the resulting state and adjoint state equations.Furthermore,a priori error estimates for state,adjoint state and control variables are discussed in details.Illustrative numerical tests are given to demonstrate the validity and applicability of our proposed approximations and theoretical results.
文摘In this paper, we discuss virtual element method (VEM) approximation of optimal control problem governed by Brinkman equations with control constraints. Based on the polynomial projections and variational discretization of the control variable, we build up the virtual element discrete scheme of the optimal control problem and derive the discrete first order optimality system. A priori error estimates for the state, adjoint state and control variables in L<sup>2</sup> and H<sup>1</sup> norm are derived. The theoretical findings are illustrated by the numerical experiments.
基金We thank the anonymous referees for their valuable comments and suggestions which lead to an improved presentation of this paper. This work was supported by NSFC under the grant 11371199, 11226334 and 11301275, the Jiangsu Provincial 2011 Program (Collaborative Innovation Center of Climate Change), the Program of Natural Science Research of Jiangsu Higher Education Institutions of China (Grant No. 12KJB110013), Natural Science Foundation of Guangdong Province of China (Grant No. S2012040007993) and Educational Commission of Guangdong Province of China (Grant No. 2012LYM0122).
文摘In this paper, we study a weakly over-penalized interior penalty method for non-self- adjoint and indefinite problems. An optimal a priori error estimate in the energy norm is derived. In addition, we introduce a residual-based a posteriori error estimator, which is proved to be both reliable and efficient in the energy norm. Some numerical testes are presented to validate our theoretical analysis.