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 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, 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 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 and analyze the interior penalty discontinuous Galerkin method for H(div)-elliptic problem.An optimal a priori error estimate in the energy norm is proved.In addition,a residual-based a poster...In this paper,we propose and analyze the interior penalty discontinuous Galerkin method for H(div)-elliptic problem.An optimal a priori error estimate in the energy norm is proved.In addition,a residual-based a posteriori error estimator is obtained.The estimator is proved to be both reliable and efficient in the energy norm.Some numerical testes are presented to demonstrate the effectiveness of our method.展开更多
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.展开更多
The goal of this paper is to study a mixed finite element approximation of the general convex optimal control problems governed by quasilinear elliptic partial differential equations. The state and co-state are approx...The goal of this paper is to study a mixed finite element approximation of the general convex optimal control problems governed by quasilinear elliptic partial differential equations. The state and co-state are approximated by the lowest order Raviart-Thomas mixed finite element spaces and the control is approximated by piecewise constant functions. We derive a priori error estimates both for the state variables and the control variable. Finally, some numerical examples are given to demonstrate the theoretical results.展开更多
Nearly all inf-sup stable mixed finite elements for the incompressible Stokes equations relax the divergence constraint. The price to pay is that a priori estimates for the ve- locity error become pressure-dependent, ...Nearly all inf-sup stable mixed finite elements for the incompressible Stokes equations relax the divergence constraint. The price to pay is that a priori estimates for the ve- locity error become pressure-dependent, while divergence-free mixed finite elements de- liver pressure-independent estimates. A recently introduced new variational crime using lowest-order Raviart-Thomas velocity reconstructions delivers a much more robust modi- fied Crouzeix-Raviart element, obeying an optimal pressure-independent discrete H1 ve- locity estimate. Refining this approach, a more sophisticated variational crime employing the lowest-order BDM element is proposed, which also allows proving an optimal pressure- independent L2 velocity error. Numerical examples confirm the analysis and demonstrate the improved robustness in the Navier-Stokes case.展开更多
Numerical atomic orbitals have been successfully used in molecular simulations as a basis set,which provides a nature,physical description of the electronic states and is suitable for ■(N)calculations based on the st...Numerical atomic orbitals have been successfully used in molecular simulations as a basis set,which provides a nature,physical description of the electronic states and is suitable for ■(N)calculations based on the strictly localized property.This paper presents a numerical analysis for some simplified atomic orbitals,with polynomial-type and confined Hydrogen-like radial basis functions respectively.We give some a priori error estimates to understand why numerical atomic orbitals are computationally efficient in electronic structure calculations.展开更多
In this paper,we discuss the conforming finite element method for a modified interior transmission eigenvalues problem.We present a complete theoretical analysis for the method,including the a priori and a posteriori ...In this paper,we discuss the conforming finite element method for a modified interior transmission eigenvalues problem.We present a complete theoretical analysis for the method,including the a priori and a posteriori error estimates.The theoretical analysis is conducted under the assumption of low regularity on the solution.We prove the reliability and efficiency of the a posteriori error estimators for eigenfunctions up to higher order terms,and we also analyze the reliability of estimators for eigenvalues.Finally,we report numerical experiments to show that our posteriori error estimator is effective and the approximations can reach the optimal convergence order.The numerical results also indicate that the conforming finite element eigenvalues approximate the exact ones from below,and there exists a monotonic relationship between the conforming finite element eigenvalues and the refractive index through numerical experiments.展开更多
This paper is concerned with an ill-posed problem which results from the area of molecular imaging and is known as BLT problem. Using Tikhonov regularization technique, a quadratic optimization problem can be formulat...This paper is concerned with an ill-posed problem which results from the area of molecular imaging and is known as BLT problem. Using Tikhonov regularization technique, a quadratic optimization problem can be formulated. We provide an improved error estimate for the finite element approximation of the regularized optimization problem. Some numerical examples are presented to demonstrate our theoretical results.展开更多
In this paper,we investigate a priori error estimates for the quadratic optimal control problems governed by semilinear elliptic partial differential equations using higher order triangular mixed finite element method...In this paper,we investigate a priori error estimates for the quadratic optimal control problems governed by semilinear elliptic partial differential equations using higher order triangular mixed finite element methods.The state and the co-state are approximated by the order k Raviart-Thomas mixed finite element spaces and the control is approximated by piecewise polynomials of order k(k≥0).A priori error estimates for the mixed finite element approximation of semilinear control problems are obtained.Finally,we present some numerical examples which confirm our theoretical results.展开更多
In this paper, we present a local discontinuous Galerkin (LDG) method for the AllenCahn equation. We prove the energy stability, analyze the optimal convergence rate of k + 1 in L2 norm and present the (2k+1)-th...In this paper, we present a local discontinuous Galerkin (LDG) method for the AllenCahn equation. We prove the energy stability, analyze the optimal convergence rate of k + 1 in L2 norm and present the (2k+1)-th order negative-norm estimate of the semi- discrete LDG method for the Allen-Cahn equation with smooth solution. To relax the severe time step restriction of explicit time marching methods, we construct a first order semi-implicit scheme based on the convex splitting principle of the discrete Allen-Cahn energy and prove the corresponding unconditional energy stability. To achieve high order temporal accuracy, we employ the semi-implicit spectral deferred correction (SDC) method. Combining with the unconditionally stable convex splitting scheme, the SDC method can be high order accurate and stable in our numerical tests. To enhance the efficiency of the proposed methods, the multigrid solver is adapted to solve the resulting nonlinear algebraic systems. Numerical studies are presented to confirm that we can achieve optimal accuracy of (O(hk+1) in L2 norm and improve the LDG solution from (O(hk+1) to (O(h2k+1) with the accuracy enhancement post-processing technique.展开更多
In this paper, we study variational discretization for the constrained optimal control problem governed by convection dominated diffusion equations, where the state equation is approximated by the edge stabilization G...In this paper, we study variational discretization for the constrained optimal control problem governed by convection dominated diffusion equations, where the state equation is approximated by the edge stabilization Galerkin method. A priori error estimates are derived for the state, the adjoint state and the control. Moreover, residual type a posteriori error estimates in the L^2-norm are obtained. Finally, two numerical experiments are presented to illustrate the theoretical results.展开更多
In this paper,we consider the finite element method for two-dimensional nonlinear modified time-fractional fourth-order diffusion equation.In order to avoid using higher order elements,we introduce an intermediate var...In this paper,we consider the finite element method for two-dimensional nonlinear modified time-fractional fourth-order diffusion equation.In order to avoid using higher order elements,we introduce an intermediate variableσ=∆u and translate the fourth-order derivative of the original problem into a second-order coupled system.We discretize the fractional time derivative terms by using the L1-approximation and discretize the first-order time derivative term by using the second-order backward differentiation formula.In the fully discrete scheme,we implement the finite element method for the spatial approximation.Unconditional stability of the fully discrete scheme is proven and its optimal convergence order is obtained.Numerical experiments are carried out to demonstrate our theoretical analysis.展开更多
In this paper, a hybird approximation scheme for an optimal control problem governed by an elliptic equation with random field in its coefficients is considered. The random coefficients are smooth in the physical spac...In this paper, a hybird approximation scheme for an optimal control problem governed by an elliptic equation with random field in its coefficients is considered. The random coefficients are smooth in the physical space and depend on a large number of random variables in the probability space. The necessary and sufficient optimality conditions for the optimal control problem are obtained. The scheme is established to approximate the optimality system through the discretization by using finite volume element method for the spatial space and a sparse grid stochastic collocation method based on the Smolyak approximation for the probability space, respectively. This scheme naturally leads to the discrete solutions of an uncoupled deterministic problem. The existence and uniqueness of the discrete solutions are proved. A priori error estimates are derived for the state, the co-state and the control variables. Numerical examples are presented to illustrate our theoretical results.展开更多
In this paper,we investigate a priori and a posteriori error estimates of fully discrete H^(1)-Galerkin mixed finite element methods for parabolic optimal control prob-lems.The state variables and co-state variables a...In this paper,we investigate a priori and a posteriori error estimates of fully discrete H^(1)-Galerkin mixed finite element methods for parabolic optimal control prob-lems.The state variables and co-state variables are approximated by the lowest order Raviart-Thomas mixed finite element and linear finite element,and the control vari-able is approximated by piecewise constant functions.The time discretization of the state and co-state are based on finite difference methods.First,we derive a priori error estimates for the control variable,the state variables and the adjoint state variables.Second,by use of energy approach,we derive a posteriori error estimates for optimal control problems,assuming that only the underlying mesh is static.A numerical example is presented to verify the theoretical results on a priori error estimates.展开更多
文摘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 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.
基金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.
文摘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.
文摘In this paper,we propose and analyze the interior penalty discontinuous Galerkin method for H(div)-elliptic problem.An optimal a priori error estimate in the energy norm is proved.In addition,a residual-based a posteriori error estimator is obtained.The estimator is proved to be both reliable and efficient in the energy norm.Some numerical testes are presented to demonstrate the effectiveness of our method.
基金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.
文摘The goal of this paper is to study a mixed finite element approximation of the general convex optimal control problems governed by quasilinear elliptic partial differential equations. The state and co-state are approximated by the lowest order Raviart-Thomas mixed finite element spaces and the control is approximated by piecewise constant functions. We derive a priori error estimates both for the state variables and the control variable. Finally, some numerical examples are given to demonstrate the theoretical results.
文摘Nearly all inf-sup stable mixed finite elements for the incompressible Stokes equations relax the divergence constraint. The price to pay is that a priori estimates for the ve- locity error become pressure-dependent, while divergence-free mixed finite elements de- liver pressure-independent estimates. A recently introduced new variational crime using lowest-order Raviart-Thomas velocity reconstructions delivers a much more robust modi- fied Crouzeix-Raviart element, obeying an optimal pressure-independent discrete H1 ve- locity estimate. Refining this approach, a more sophisticated variational crime employing the lowest-order BDM element is proposed, which also allows proving an optimal pressure- independent L2 velocity error. Numerical examples confirm the analysis and demonstrate the improved robustness in the Navier-Stokes case.
基金The research for this paper has been enabled by the Alexander von Humboldt Foundation,whose support for the long term visit of Huajie Chen at Technische Universit¨at Berlin is gratefully acknowledged.
文摘Numerical atomic orbitals have been successfully used in molecular simulations as a basis set,which provides a nature,physical description of the electronic states and is suitable for ■(N)calculations based on the strictly localized property.This paper presents a numerical analysis for some simplified atomic orbitals,with polynomial-type and confined Hydrogen-like radial basis functions respectively.We give some a priori error estimates to understand why numerical atomic orbitals are computationally efficient in electronic structure calculations.
基金supported by the National Natural Science Foundation of China(Nos.12261024,11561014)Science and Technology Planning Project of Guizhou Province(Guizhou Kehe fundamental research-ZK[2022]No.324).
文摘In this paper,we discuss the conforming finite element method for a modified interior transmission eigenvalues problem.We present a complete theoretical analysis for the method,including the a priori and a posteriori error estimates.The theoretical analysis is conducted under the assumption of low regularity on the solution.We prove the reliability and efficiency of the a posteriori error estimators for eigenfunctions up to higher order terms,and we also analyze the reliability of estimators for eigenvalues.Finally,we report numerical experiments to show that our posteriori error estimator is effective and the approximations can reach the optimal convergence order.The numerical results also indicate that the conforming finite element eigenvalues approximate the exact ones from below,and there exists a monotonic relationship between the conforming finite element eigenvalues and the refractive index through numerical experiments.
基金the National Natural Science Foundation of China (No. 60474027,10301003 and 10771211)the National Basic Research Program under the Grant 2005CB321701
文摘This paper is concerned with an ill-posed problem which results from the area of molecular imaging and is known as BLT problem. Using Tikhonov regularization technique, a quadratic optimization problem can be formulated. We provide an improved error estimate for the finite element approximation of the regularized optimization problem. Some numerical examples are presented to demonstrate our theoretical results.
基金supported by the Foundation for Talent Introduction of Guangdong Provincial Universities and CollegesPearl River Scholar Funded Scheme(2008)National Science Foundation of China(10971074).
文摘In this paper,we investigate a priori error estimates for the quadratic optimal control problems governed by semilinear elliptic partial differential equations using higher order triangular mixed finite element methods.The state and the co-state are approximated by the order k Raviart-Thomas mixed finite element spaces and the control is approximated by piecewise polynomials of order k(k≥0).A priori error estimates for the mixed finite element approximation of semilinear control problems are obtained.Finally,we present some numerical examples which confirm our theoretical results.
文摘In this paper, we present a local discontinuous Galerkin (LDG) method for the AllenCahn equation. We prove the energy stability, analyze the optimal convergence rate of k + 1 in L2 norm and present the (2k+1)-th order negative-norm estimate of the semi- discrete LDG method for the Allen-Cahn equation with smooth solution. To relax the severe time step restriction of explicit time marching methods, we construct a first order semi-implicit scheme based on the convex splitting principle of the discrete Allen-Cahn energy and prove the corresponding unconditional energy stability. To achieve high order temporal accuracy, we employ the semi-implicit spectral deferred correction (SDC) method. Combining with the unconditionally stable convex splitting scheme, the SDC method can be high order accurate and stable in our numerical tests. To enhance the efficiency of the proposed methods, the multigrid solver is adapted to solve the resulting nonlinear algebraic systems. Numerical studies are presented to confirm that we can achieve optimal accuracy of (O(hk+1) in L2 norm and improve the LDG solution from (O(hk+1) to (O(h2k+1) with the accuracy enhancement post-processing technique.
基金support of the Chinese and German Research Foundations through the Sino-German Workshop on Applied Mathematics held in Hangzhou in October 2007support of the German Research Foundation through the grants DFG06-381 and DFG06-382+1 种基金support of the National Basic Research Program under the Grant 2005CB321701the National Natural Science Foundation of China under the Grant 60474027 and 10771211
文摘In this paper, we study variational discretization for the constrained optimal control problem governed by convection dominated diffusion equations, where the state equation is approximated by the edge stabilization Galerkin method. A priori error estimates are derived for the state, the adjoint state and the control. Moreover, residual type a posteriori error estimates in the L^2-norm are obtained. Finally, two numerical experiments are presented to illustrate the theoretical results.
文摘In this paper,we consider the finite element method for two-dimensional nonlinear modified time-fractional fourth-order diffusion equation.In order to avoid using higher order elements,we introduce an intermediate variableσ=∆u and translate the fourth-order derivative of the original problem into a second-order coupled system.We discretize the fractional time derivative terms by using the L1-approximation and discretize the first-order time derivative term by using the second-order backward differentiation formula.In the fully discrete scheme,we implement the finite element method for the spatial approximation.Unconditional stability of the fully discrete scheme is proven and its optimal convergence order is obtained.Numerical experiments are carried out to demonstrate our theoretical analysis.
基金The work is partly supported by the Natural Science Foundation of China (Grant No. 11271231, 11301300, 61572297), by the Shandong Province Outstanding Young Scientists Research Award Fhnd Project (Grant No. BS2013DX010), by the Natural Science Foundation of Shandong Province, China (Grant No. ZR2014FM003), and by the Shandong Academy of Sciences Youth Fund Project (Grant No. 2013QN007).
文摘In this paper, a hybird approximation scheme for an optimal control problem governed by an elliptic equation with random field in its coefficients is considered. The random coefficients are smooth in the physical space and depend on a large number of random variables in the probability space. The necessary and sufficient optimality conditions for the optimal control problem are obtained. The scheme is established to approximate the optimality system through the discretization by using finite volume element method for the spatial space and a sparse grid stochastic collocation method based on the Smolyak approximation for the probability space, respectively. This scheme naturally leads to the discrete solutions of an uncoupled deterministic problem. The existence and uniqueness of the discrete solutions are proved. A priori error estimates are derived for the state, the co-state and the control variables. Numerical examples are presented to illustrate our theoretical results.
基金This work was supported by National Natural Science Foundation of China(11601014,11626037,11526036)China Postdoctoral Science Foundation(2016M 601359)+4 种基金Scientific and Technological Developing Scheme of Jilin Province(20160520108 JH,20170101037JC)Science and Technology Research Project of Jilin Provincial Depart-ment of Education(201646)Special Funding for Promotion of Young Teachers of Beihua University,Natural Science Foundation of Hunan Province(14JJ3135)the Youth Project of Hunan Provincial Education Department(15B096)the construct program of the key discipline in Hunan University of Science and Engineering.
文摘In this paper,we investigate a priori and a posteriori error estimates of fully discrete H^(1)-Galerkin mixed finite element methods for parabolic optimal control prob-lems.The state variables and co-state variables are approximated by the lowest order Raviart-Thomas mixed finite element and linear finite element,and the control vari-able is approximated by piecewise constant functions.The time discretization of the state and co-state are based on finite difference methods.First,we derive a priori error estimates for the control variable,the state variables and the adjoint state variables.Second,by use of energy approach,we derive a posteriori error estimates for optimal control problems,assuming that only the underlying mesh is static.A numerical example is presented to verify the theoretical results on a priori error estimates.