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,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.展开更多
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.展开更多
The subject of this work is to propose adaptive finite element methods based on an optimal maximum norm error control estimate.Using estimators of the local regularity of the unknown exact solution derived from comput...The subject of this work is to propose adaptive finite element methods based on an optimal maximum norm error control estimate.Using estimators of the local regularity of the unknown exact solution derived from computed approximate solutions,the proposed procedures are analyzed in detail for a non-trivial class of corner problems and shown to be efficient in the sense that they generate the correct type of refinement and lead to the desired control under consideration.展开更多
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.展开更多
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.展开更多
This paper considers the finite element approximation to parabolic optimal control problems with measure data in a nonconvex polygonal domain.Such problems usually possess low regularity in the state variable due to t...This paper considers the finite element approximation to parabolic optimal control problems with measure data in a nonconvex polygonal domain.Such problems usually possess low regularity in the state variable due to the presence of measure data and the nonconvex nature of the domain.The low regularity of the solution allows the finite element approximations to converge at lower orders.We prove the existence,uniqueness and regularity results for the solution to the control problem satisfying the first order optimality condition.For our error analysis we have used piecewise linear elements for the approximation of the state and co-state variables,whereas piecewise constant functions are employed to approximate the control variable.The temporal discretization is based on the implicit Euler scheme.We derive both a priori and a posteriori error bounds for the state,control and co-state variables.Numerical experiments are performed to validate the theoretical rates of convergence.展开更多
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.展开更多
A mixed finite element method is presented for the Biot consolidation problem in poroe-lasticity.More precisely,the displacement is approximated by using the Crouzeix-Raviart nonconforming finite elements,while the fu...A mixed finite element method is presented for the Biot consolidation problem in poroe-lasticity.More precisely,the displacement is approximated by using the Crouzeix-Raviart nonconforming finite elements,while the fuid pressure is approximated by using the node conforming finite elements.The well-posedness of the fully discrete scheme is established,and a corresponding priori error estimate with optimal order in the energy norm is also derived.Numerical experiments are provided to validate 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.展开更多
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 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.
基金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.
基金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.
文摘The subject of this work is to propose adaptive finite element methods based on an optimal maximum norm error control estimate.Using estimators of the local regularity of the unknown exact solution derived from computed approximate solutions,the proposed procedures are analyzed in detail for a non-trivial class of corner problems and shown to be efficient in the sense that they generate the correct type of refinement and lead to the desired control under consideration.
基金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.
文摘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.
文摘This paper considers the finite element approximation to parabolic optimal control problems with measure data in a nonconvex polygonal domain.Such problems usually possess low regularity in the state variable due to the presence of measure data and the nonconvex nature of the domain.The low regularity of the solution allows the finite element approximations to converge at lower orders.We prove the existence,uniqueness and regularity results for the solution to the control problem satisfying the first order optimality condition.For our error analysis we have used piecewise linear elements for the approximation of the state and co-state variables,whereas piecewise constant functions are employed to approximate the control variable.The temporal discretization is based on the implicit Euler scheme.We derive both a priori and a posteriori error bounds for the state,control and co-state variables.Numerical experiments are performed to validate the theoretical rates of convergence.
文摘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.
基金supported by the Guangdong Basic and Applied Basic Research Foundation(Grant No.2020A1515011032)The work of M.Cai is supported in part by the NIH-BUILD(Grant No.UL1GM118973)+2 种基金by the NIH-RCMI(Grant No.U54MD013376)the National Science Foundation awards(Grant Nos.1700328,1831950)The work of L.Zhong is supported by the National Natural Science Foundation of China(Grant No.12071160)。
文摘A mixed finite element method is presented for the Biot consolidation problem in poroe-lasticity.More precisely,the displacement is approximated by using the Crouzeix-Raviart nonconforming finite elements,while the fuid pressure is approximated by using the node conforming finite elements.The well-posedness of the fully discrete scheme is established,and a corresponding priori error estimate with optimal order in the energy norm is also derived.Numerical experiments are provided to validate 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.
基金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.