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.展开更多
This study investigates the use of dynamic a priori error information according to atmospheric moistness and the use of quality controls in temperature and water vapor profile retrievals from hyperspectral infrared ...This study investigates the use of dynamic a priori error information according to atmospheric moistness and the use of quality controls in temperature and water vapor profile retrievals from hyperspectral infrared (IR) sounders. Temperature and water vapor profiles are retrieved from Atmospheric InfraRed Sounder (AIRS) radiance measurements by applying a physical iterative method using regression retrieval as the first guess. Based on the dependency of first-guess errors on the degree of atmospheric moistness, the a priori first-guess errors classified by total precipitable water (TPW) are applied in the AIRS physical retrieval procedure. Compared to the retrieval results from a fixed a priori error, boundary layer moisture retrievals appear to be improved via TPW classification of a priori first-guess errors. Six quality control (QC) tests, which check non-converged or bad retrievals, large residuals, high terrain and desert areas, and large temperature and moisture deviations from the first guess regression retrieval, are also applied in the AIRS physical retrievals. Significantly large errors are found for the retrievals rejected by these six QCs, and the retrieval errors are substantially reduced via QC over land, which suggest the usefulness and high impact of the QCs, especially over land. In conclusion, the use of dynamic a priori error information according to atmospheric moistness, and the use of appropriate QCs dealing with the geographical information and the deviation from the first-guess as well as the conventional inverse performance are suggested to improve temperature and moisture retrievals and their applications.展开更多
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.展开更多
A combined approximation for a kind of compressible miscible displacement problems including molecular diffusion and dispersion in porous media is studied. Mixed finite el- ement method is applied to the ftow equation...A combined approximation for a kind of compressible miscible displacement problems including molecular diffusion and dispersion in porous media is studied. Mixed finite el- ement method is applied to the ftow equation, and the transport one is solved by the symmetric interior penalty discontinuous Galerkin method (SIPG). To avoid the inconve- nience of the cut-off operator in [3, 21], some induction hypotheses different from the ones in [6] are used. Based on interpolation projection properties, a priori hp error estimates are obtained. Comparing with the existing error analysis that only deals with the diffusion case, the current work is more complicated and more significant.展开更多
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 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 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.展开更多
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.展开更多
A new higher-order accurate space-time discontinuous Galerkin(DG)method using the interior penalty flux and discontinuous basis functions,both in space and in time,is pre-sented and fully analyzed for the second-order...A new higher-order accurate space-time discontinuous Galerkin(DG)method using the interior penalty flux and discontinuous basis functions,both in space and in time,is pre-sented and fully analyzed for the second-order scalar wave equation.Special attention is given to the definition of the numerical fluxes since they are crucial for the stability and accuracy of the space-time DG method.The theoretical analysis shows that the DG discre-tization is stable and converges in a DG-norm on general unstructured and locally refined meshes,including local refinement in time.The space-time interior penalty DG discre-tization does not have a CFL-type restriction for stability.Optimal order of accuracy is obtained in the DG-norm if the mesh size h and the time stepΔt satisfy h≅CΔt,with C a positive constant.The optimal order of accuracy of the space-time DG discretization in the DG-norm is confirmed by calculations on several model problems.These calculations also show that for pth-order tensor product basis functions the convergence rate in the L∞and L2-norms is order p+1 for polynomial orders p=1 and p=3 and order p for polynomial order p=2.展开更多
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 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 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 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.展开更多
Under two hypotheses of nonconforming finite elements of fourth order elliptic problems,we present a side-patchwise projection based error analysis method(SPP-BEAM for short).Such a method is able to avoid both the re...Under two hypotheses of nonconforming finite elements of fourth order elliptic problems,we present a side-patchwise projection based error analysis method(SPP-BEAM for short).Such a method is able to avoid both the regularity condition of exact solutions in the classical error analysis method and the complicated bubble function technique in the recent medius error analysis method.In addition,it is universal enough to admit generalizations.Then,we propose a sufficient condition for these hypotheses by imposing a set of in some sense necessary degrees of freedom of the shape function spaces.As an application,we use the theory to design a P3 second order triangular H2 non-conforming element by enriching two P4 bubble functions and,another P4 second order triangular H2 nonconforming finite element,and a P3 second order tetrahedral H2 non-conforming element by enriching eight P4 bubble functions,adding some more degrees of freedom.展开更多
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.展开更多
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 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 GOES-R Algorithm Working Group Program and GOES-R High Impact Weather Project (Grant No NA10NES4400013)supported by the Korea Meteorological Administration Research and Development Program under Grant CATER 2006-2103the BK21 Project of the Korean Government
文摘This study investigates the use of dynamic a priori error information according to atmospheric moistness and the use of quality controls in temperature and water vapor profile retrievals from hyperspectral infrared (IR) sounders. Temperature and water vapor profiles are retrieved from Atmospheric InfraRed Sounder (AIRS) radiance measurements by applying a physical iterative method using regression retrieval as the first guess. Based on the dependency of first-guess errors on the degree of atmospheric moistness, the a priori first-guess errors classified by total precipitable water (TPW) are applied in the AIRS physical retrieval procedure. Compared to the retrieval results from a fixed a priori error, boundary layer moisture retrievals appear to be improved via TPW classification of a priori first-guess errors. Six quality control (QC) tests, which check non-converged or bad retrievals, large residuals, high terrain and desert areas, and large temperature and moisture deviations from the first guess regression retrieval, are also applied in the AIRS physical retrievals. Significantly large errors are found for the retrievals rejected by these six QCs, and the retrieval errors are substantially reduced via QC over land, which suggest the usefulness and high impact of the QCs, especially over land. In conclusion, the use of dynamic a priori error information according to atmospheric moistness, and the use of appropriate QCs dealing with the geographical information and the deviation from the first-guess as well as the conventional inverse performance are suggested to improve temperature and moisture retrievals and their applications.
基金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.
文摘A combined approximation for a kind of compressible miscible displacement problems including molecular diffusion and dispersion in porous media is studied. Mixed finite el- ement method is applied to the ftow equation, and the transport one is solved by the symmetric interior penalty discontinuous Galerkin method (SIPG). To avoid the inconve- nience of the cut-off operator in [3, 21], some induction hypotheses different from the ones in [6] are used. Based on interpolation projection properties, a priori hp error estimates are obtained. Comparing with the existing error analysis that only deals with the diffusion case, the current work is more complicated and more significant.
基金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.
文摘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.
基金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.
基金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.
文摘A new higher-order accurate space-time discontinuous Galerkin(DG)method using the interior penalty flux and discontinuous basis functions,both in space and in time,is pre-sented and fully analyzed for the second-order scalar wave equation.Special attention is given to the definition of the numerical fluxes since they are crucial for the stability and accuracy of the space-time DG method.The theoretical analysis shows that the DG discre-tization is stable and converges in a DG-norm on general unstructured and locally refined meshes,including local refinement in time.The space-time interior penalty DG discre-tization does not have a CFL-type restriction for stability.Optimal order of accuracy is obtained in the DG-norm if the mesh size h and the time stepΔt satisfy h≅CΔt,with C a positive constant.The optimal order of accuracy of the space-time DG discretization in the DG-norm is confirmed by calculations on several model problems.These calculations also show that for pth-order tensor product basis functions the convergence rate in the L∞and L2-norms is order p+1 for polynomial orders p=1 and p=3 and order p for polynomial order p=2.
文摘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.
基金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 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 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.
文摘Under two hypotheses of nonconforming finite elements of fourth order elliptic problems,we present a side-patchwise projection based error analysis method(SPP-BEAM for short).Such a method is able to avoid both the regularity condition of exact solutions in the classical error analysis method and the complicated bubble function technique in the recent medius error analysis method.In addition,it is universal enough to admit generalizations.Then,we propose a sufficient condition for these hypotheses by imposing a set of in some sense necessary degrees of freedom of the shape function spaces.As an application,we use the theory to design a P3 second order triangular H2 non-conforming element by enriching two P4 bubble functions and,another P4 second order triangular H2 nonconforming finite element,and a P3 second order tetrahedral H2 non-conforming element by enriching eight P4 bubble functions,adding some more degrees of freedom.
文摘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 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.