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 maximum a posteriori( MAP) algorithm is proposed to improve the accuracy of super resolution( SR) reconstruction in traditional methods. The algorithm applies both joints image registration and SR reconstruction...A maximum a posteriori( MAP) algorithm is proposed to improve the accuracy of super resolution( SR) reconstruction in traditional methods. The algorithm applies both joints image registration and SR reconstruction in the framework,but separates them in the process of iteratiion. Firstly,we estimate the shifting parameters through two lowresolution( LR) images and use the parameters to reconstruct initial HR images. Then,we update the shifting parameters using HR images. The aforementioned steps are repeated until the ideal HR images are obtained. The metrics such as PSNR and SSIM are used to fully evaluate the quality of the reconstructed image. Experimental results indicate that the proposed method can enhance image resolution efficiently.展开更多
In this paper, a posteriori error estimates for the generalized Schwartz method with Dirichlet boundary conditions on the interfaces for advection-diffusion equation with second order boundary value problems are prove...In this paper, a posteriori error estimates for the generalized Schwartz method with Dirichlet boundary conditions on the interfaces for advection-diffusion equation with second order boundary value problems are proved by using the Euler time scheme combined with Galerkin spatial method. Furthermore, an asymptotic behavior in Sobolev norm is de- duced using Benssoussau-Lions' algorithm. Finally, the results of some numerical experiments are presented to support the theory.展开更多
Residual based on a posteriori error estimates for conforming finite element solutions of incompressible Navier-Stokes equations with stream function form which were computed with seven recently proposed two-level met...Residual based on a posteriori error estimates for conforming finite element solutions of incompressible Navier-Stokes equations with stream function form which were computed with seven recently proposed two-level method were derived. The posteriori error estimates contained additional terms in comparison to the error estimates for the solution obtained by the standard finite element method. The importance of these additional terms in the error estimates was investigated by studying their asymptotic behavior. For optimal scaled meshes, these bounds are not of higher order than of convergence of discrete solution.展开更多
In this paper, a posteriori error estimates were derived for piecewise linear finite element approximations to parabolic obstacle problems. The instrumental ingredient was introduced as a new interpolation operator wh...In this paper, a posteriori error estimates were derived for piecewise linear finite element approximations to parabolic obstacle problems. The instrumental ingredient was introduced as a new interpolation operator which has optimal approximation properties and preserves positivity. With the help of the interpolation operator the upper and lower bounds were obtained.展开更多
The present study regards the numerical approximation of solutions of systems of Korteweg-de Vries type,coupled through their nonlinear terms.In our previous work[9],we constructed conservative and dissipative finite ...The present study regards the numerical approximation of solutions of systems of Korteweg-de Vries type,coupled through their nonlinear terms.In our previous work[9],we constructed conservative and dissipative finite element methods for these systems and presented a priori error estimates for the semidiscrete schemes.In this sequel,we present a posteriori error estimates for the semidiscrete and fully discrete approximations introduced in[9].The key tool employed to effect our analysis is the dispersive reconstruction devel-oped by Karakashian and Makridakis[20]for related discontinuous Galerkin methods.We conclude by providing a set of numerical experiments designed to validate the a posteriori theory and explore the effectivity of the resulting error indicators.展开更多
The paper presents an algorithm of automatic target detection in Synthetic Aperture Radar(SAR) images based on Maximum A Posteriori(MAP). The algorithm is divided into three steps. First, it employs Gaussian mixture d...The paper presents an algorithm of automatic target detection in Synthetic Aperture Radar(SAR) images based on Maximum A Posteriori(MAP). The algorithm is divided into three steps. First, it employs Gaussian mixture distribution to approximate and estimate multi-modal histogram of SAR image. Then, based on the principle of MAP, when a priori probability is both unknown and learned respectively, the sample pixels are classified into different classes c = {target,shadow, background}. Last, it compares the results of two different target detections. Simulation results preferably indicate that the presented algorithm is fast and robust, with the learned a priori probability, an approach to target detection is reliable and promising.展开更多
Interferogram noise reduction is a very important processing step in Interferometric Synthetic Aperture Radar(InSAR) technique. The most difficulty for this step is to remove the noises and preserve the fringes simult...Interferogram noise reduction is a very important processing step in Interferometric Synthetic Aperture Radar(InSAR) technique. The most difficulty for this step is to remove the noises and preserve the fringes simultaneously. To solve the dilemma, a new interferogram noise reduction algorithm based on the Maximum A Posteriori(MAP) estimate is introduced in this paper. The algorithm is solved under the Total Generalized Variation(TGV) minimization assumption, which exploits the phase characteristics up to the second order differentiation. The ideal noise-free phase consisting of piecewise smooth areas is involved in this assumption, which is coincident with the natural terrain. In order to overcome the phase wraparound effect, complex plane filter is utilized in this algorithm. The simulation and real data experiments show the algorithm can reduce the noises effectively and meanwhile preserve the interferogram fringes very well.展开更多
Based on the auxiliary subspace techniques,a posteriori error estimator of nonconforming weak Galerkin finite element method(WGFEM)for Stokes problem in two and three dimensions is presented.Without saturation assumpt...Based on the auxiliary subspace techniques,a posteriori error estimator of nonconforming weak Galerkin finite element method(WGFEM)for Stokes problem in two and three dimensions is presented.Without saturation assumption,we prove that the WGFEM approximation error is bounded by the error estimator up to an oscillation term.The computational cost of the approximation and the error problems is considered in terms of size and sparsity of the system matrix.To reduce the computational cost of the error problem,an equivalent error problem is constructed by using diagonalization techniques,which needs to solve only two diagonal linear algebraic systems corresponding to the degree of freedom(d.o.f)to get the error estimator.Numerical experiments are provided to demonstrate the effectiveness and robustness of the a posteriori error estimator.展开更多
In this paper,we study the a posteriori error estimator of SDG method for variable coefficients time-harmonic Maxwell's equations.We propose two a posteriori error estimators,one is the recovery-type estimator,and...In this paper,we study the a posteriori error estimator of SDG method for variable coefficients time-harmonic Maxwell's equations.We propose two a posteriori error estimators,one is the recovery-type estimator,and the other is the residual-type estimator.We first propose the curl-recovery method for the staggered discontinuous Galerkin method(SDGM),and based on the super-convergence result of the postprocessed solution,an asymptotically exact error estimator is constructed.The residual-type a posteriori error estimator is also proposed,and it's reliability and effectiveness are proved for variable coefficients time-harmonic Maxwell's equations.The efficiency and robustness of the proposed estimators is demonstrated by the numerical experiments.展开更多
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 addresses tensile shock physics in thermoviscoelastic (TVE) solids without memory. The mathematical model is derived using conservation and balance laws (CBL) of classical continuum mechanics (CCM), incorpo...This paper addresses tensile shock physics in thermoviscoelastic (TVE) solids without memory. The mathematical model is derived using conservation and balance laws (CBL) of classical continuum mechanics (CCM), incorporating the contravariant second Piola-Kirchhoff stress tensor, the covariant Green’s strain tensor, and its rates up to order n. This mathematical model permits the study of finite deformation and finite strain compressible deformation physics with an ordered rate dissipation mechanism. Constitutive theories are derived using conjugate pairs in entropy inequality and the representation theorem. The resulting mathematical model is both thermodynamically and mathematically consistent and has closure. The solution of the initial value problems (IVPs) describing evolutions is obtained using a variationally consistent space-time coupled finite element method, derived using space-time residual functional in which the local approximations are in hpk higher-order scalar product spaces. This permits accurate description problem physics over the discretization and also permits precise a posteriori computation of the space-time residual functional, an accurate measure of the accuracy of the computed solution. Model problem studies are presented to demonstrate tensile shock formation, propagation, reflection, and interaction. A unique feature of this research is that tensile shocks can only exist in solid matter, as their existence requires a medium to be elastic (presence of strain), which is only possible in a solid medium. In tensile shock physics, a decrease in the density of the medium caused by tensile waves leads to shock formation ahead of the wave. In contrast, in compressive shocks, an increase in density and the corresponding compressive waves result in the formation of compression shocks behind of the wave. Although these are two similar phenomena, they are inherently different in nature. To our knowledge, this work has not been reported in the published literature.展开更多
This paper aims at a general guideline to obtain a posteriori error estimates for the finite element error control in computational partial differential equations.In the abstract setting of mixed formulations,a genera...This paper aims at a general guideline to obtain a posteriori error estimates for the finite element error control in computational partial differential equations.In the abstract setting of mixed formulations,a generalised formulation of the corresponding residuals is proposed which then allows for the unified estimation of the respective dual norms.Notably,this can be done with an approach which is applicable in the same way to conforming,nonconforming and mixed discretisations.Subsequently,the unified approach is applied to various model problems.In particular,we consider the Laplace,Stokes,Navier-Lamé,and the semi-discrete eddy current equations.展开更多
Four primal discontinuous Galerkin methods are applied to solve reactive transport problems, namely, Oden-BabuSka-Baumann DG (OBB-DG), non-symmetric interior penalty Galerkin (NIPG), symmetric interior penalty Gal...Four primal discontinuous Galerkin methods are applied to solve reactive transport problems, namely, Oden-BabuSka-Baumann DG (OBB-DG), non-symmetric interior penalty Galerkin (NIPG), symmetric interior penalty Galerkin (SIPG), and incomplete interior penalty Galerkin (IIPG). A unified a posteriori residual-type error estimation is derived explicitly for these methods. From the computed solution and given data, explicit estimators can be computed efficiently and directly, which can be used as error indicators for adaptation. Unlike in the reference [10], we obtain the error estimators in L^2 (L^2) norm by using duality techniques instead of in L^2(H^1) norm.展开更多
In this paper,we investigate the Legendre Galerkin spectral approximation of quadratic optimal control problems governed by parabolic equations.A spectral approximation scheme for the parabolic optimal control problem...In this paper,we investigate the Legendre Galerkin spectral approximation of quadratic optimal control problems governed by parabolic equations.A spectral approximation scheme for the parabolic optimal control problem is presented.We obtain a posteriori error estimates of the approximated solutions for both the state and the control.展开更多
A unified a posteriori error analysis has been developed in [18, 21-23] to analyze the finite element error a posteriori under a universal roof. This paper contributes to the finite element meshes with hanging nodes w...A unified a posteriori error analysis has been developed in [18, 21-23] to analyze the finite element error a posteriori under a universal roof. This paper contributes to the finite element meshes with hanging nodes which are required for local mesh-refining. The twodimensional 1-irregular triangulations into triangles and parallelograms and their combinations are considered with conforming and nonconforming finite element methods named after or by Courant, Q1, Crouzeix-Raviart, Poisson, Stokes and Navier-Lamé equations Han, Rannacher-Turek, and others for the The paper provides a unified a priori and a posteriori error analysis for triangulations with hanging nodes of degree ≤ 1 which are fundamental for local mesh refinement in self-adaptive finite element discretisations.展开更多
Two residual-based a posteriori error estimators of the nonconforming Crouzeix-Raviart element are derived for elliptic problems with Dirac delta source terms.One estimator is shown to be reliable and efficient,which ...Two residual-based a posteriori error estimators of the nonconforming Crouzeix-Raviart element are derived for elliptic problems with Dirac delta source terms.One estimator is shown to be reliable and efficient,which yields global upper and lower bounds for the error in piecewise W1,p seminorm.The other one is proved to give a global upper bound of the error in Lp-norm.By taking the two estimators as refinement indicators,adaptive algorithms are suggested,which are experimentally shown to attain optimal convergence orders.展开更多
This paper discusses conforming mixed finite element approximations for the Stokes eigenvalue problem. Firstly, several mixed finite element identities are proved. Based on these identities, the following new results ...This paper discusses conforming mixed finite element approximations for the Stokes eigenvalue problem. Firstly, several mixed finite element identities are proved. Based on these identities, the following new results are given: (1) It is proved that the numerical eigenvalues obtained by mini-element, P1-P1 element and Q1-Q1 element approximate the exact eigenvalues from above. (2) As for the P1-P1, Q1-Q1 and Q1-Po element eigenvalues, the asymptotically exact a posteriori error indicators are presented. (3) The reliable and efficient a posteriori error estimator proposed by Verfiirth is applied to mini-element eigenfunctions. Finally, numerical experiments are carried out to verify the theoretical analysis.展开更多
This paper develops a posteriori error estimates of residual type for conforming and mixed finite element approximations of the fourth order Cahn-Hilliard equation ut + △(ε△Au-ε^-1f(u)) = 0. It is shown that ...This paper develops a posteriori error estimates of residual type for conforming and mixed finite element approximations of the fourth order Cahn-Hilliard equation ut + △(ε△Au-ε^-1f(u)) = 0. It is shown that the a posteriori error bounds depends on ε^-1 only in some low polynomial order, instead of exponential order. Using these a posteriori error estimates, we construct at2 adaptive algorithm for computing the solution of the Cahn- Hilliard equation and its sharp interface limit, the Hele-Shaw flow. Numerical experiments are presented to show the robustness and effectiveness of the new error estimators and the proposed adaptive algorithm.展开更多
In this paper, we discuss the a posteriori error estimate of the finite element approximation for the boundary control problems governed by the parabolic partial differential equations. Three different a posteriori er...In this paper, we discuss the a posteriori error estimate of the finite element approximation for the boundary control problems governed by the parabolic partial differential equations. Three different a posteriori error estimators are provided for the parabolic boundary control problems with the observations of the distributed state, the boundary state and the final state. It is proven that these estimators are reliable bounds of the finite element approximation errors, which can be used as the indicators of the mesh refinement in adaptive finite element methods.展开更多
基金supported by the National Basic Research Program under the Grant 2005CB321701the National Natural Science Foundation of China under the Grants 60474027 and 10771211.
文摘In this paper,we investigate a streamline diffusion finite element approxi- mation scheme for the constrained optimal control problem governed by linear con- vection dominated diffusion equations.We prove the existence and uniqueness of the discretized scheme.Then a priori and a posteriori error estimates are derived for the state,the co-state and the control.Three numerical examples are presented to illustrate our theoretical results.
基金Supported by the National Natural Science Foundation of China(61405191)
文摘A maximum a posteriori( MAP) algorithm is proposed to improve the accuracy of super resolution( SR) reconstruction in traditional methods. The algorithm applies both joints image registration and SR reconstruction in the framework,but separates them in the process of iteratiion. Firstly,we estimate the shifting parameters through two lowresolution( LR) images and use the parameters to reconstruct initial HR images. Then,we update the shifting parameters using HR images. The aforementioned steps are repeated until the ideal HR images are obtained. The metrics such as PSNR and SSIM are used to fully evaluate the quality of the reconstructed image. Experimental results indicate that the proposed method can enhance image resolution efficiently.
文摘In this paper, a posteriori error estimates for the generalized Schwartz method with Dirichlet boundary conditions on the interfaces for advection-diffusion equation with second order boundary value problems are proved by using the Euler time scheme combined with Galerkin spatial method. Furthermore, an asymptotic behavior in Sobolev norm is de- duced using Benssoussau-Lions' algorithm. Finally, the results of some numerical experiments are presented to support the theory.
文摘Residual based on a posteriori error estimates for conforming finite element solutions of incompressible Navier-Stokes equations with stream function form which were computed with seven recently proposed two-level method were derived. The posteriori error estimates contained additional terms in comparison to the error estimates for the solution obtained by the standard finite element method. The importance of these additional terms in the error estimates was investigated by studying their asymptotic behavior. For optimal scaled meshes, these bounds are not of higher order than of convergence of discrete solution.
基金Project supported by National Natural Science Foundation ofChina (Grant No .10471089)
文摘In this paper, a posteriori error estimates were derived for piecewise linear finite element approximations to parabolic obstacle problems. The instrumental ingredient was introduced as a new interpolation operator which has optimal approximation properties and preserves positivity. With the help of the interpolation operator the upper and lower bounds were obtained.
基金This work was supported in part by the National Science Foundation under grant DMS-1620288。
文摘The present study regards the numerical approximation of solutions of systems of Korteweg-de Vries type,coupled through their nonlinear terms.In our previous work[9],we constructed conservative and dissipative finite element methods for these systems and presented a priori error estimates for the semidiscrete schemes.In this sequel,we present a posteriori error estimates for the semidiscrete and fully discrete approximations introduced in[9].The key tool employed to effect our analysis is the dispersive reconstruction devel-oped by Karakashian and Makridakis[20]for related discontinuous Galerkin methods.We conclude by providing a set of numerical experiments designed to validate the a posteriori theory and explore the effectivity of the resulting error indicators.
文摘The paper presents an algorithm of automatic target detection in Synthetic Aperture Radar(SAR) images based on Maximum A Posteriori(MAP). The algorithm is divided into three steps. First, it employs Gaussian mixture distribution to approximate and estimate multi-modal histogram of SAR image. Then, based on the principle of MAP, when a priori probability is both unknown and learned respectively, the sample pixels are classified into different classes c = {target,shadow, background}. Last, it compares the results of two different target detections. Simulation results preferably indicate that the presented algorithm is fast and robust, with the learned a priori probability, an approach to target detection is reliable and promising.
文摘Interferogram noise reduction is a very important processing step in Interferometric Synthetic Aperture Radar(InSAR) technique. The most difficulty for this step is to remove the noises and preserve the fringes simultaneously. To solve the dilemma, a new interferogram noise reduction algorithm based on the Maximum A Posteriori(MAP) estimate is introduced in this paper. The algorithm is solved under the Total Generalized Variation(TGV) minimization assumption, which exploits the phase characteristics up to the second order differentiation. The ideal noise-free phase consisting of piecewise smooth areas is involved in this assumption, which is coincident with the natural terrain. In order to overcome the phase wraparound effect, complex plane filter is utilized in this algorithm. The simulation and real data experiments show the algorithm can reduce the noises effectively and meanwhile preserve the interferogram fringes very well.
基金the Natural Science Foundation of Jiangsu Province(No.BK20210540)the Natural Science Foundation of The Jiangsu Higher Education Institutions of China(No.21KJB110015)the National Key Research and Development Program of China(grant no.2020YFA0713601).
文摘Based on the auxiliary subspace techniques,a posteriori error estimator of nonconforming weak Galerkin finite element method(WGFEM)for Stokes problem in two and three dimensions is presented.Without saturation assumption,we prove that the WGFEM approximation error is bounded by the error estimator up to an oscillation term.The computational cost of the approximation and the error problems is considered in terms of size and sparsity of the system matrix.To reduce the computational cost of the error problem,an equivalent error problem is constructed by using diagonalization techniques,which needs to solve only two diagonal linear algebraic systems corresponding to the degree of freedom(d.o.f)to get the error estimator.Numerical experiments are provided to demonstrate the effectiveness and robustness of the a posteriori error estimator.
基金supported by NSFC Projects(Nos.11771371,12171411,11971410)Project of Scientific Research Fund of Hunan Provincial Science and Technology Department(No.2018WK4006)+1 种基金Project of Scientific Research Fund of Hunan Provincial Science and Technology Department,China(No.2020ZYT003)National defense basic scientific research program JCKY2019403D001.
文摘In this paper,we study the a posteriori error estimator of SDG method for variable coefficients time-harmonic Maxwell's equations.We propose two a posteriori error estimators,one is the recovery-type estimator,and the other is the residual-type estimator.We first propose the curl-recovery method for the staggered discontinuous Galerkin method(SDGM),and based on the super-convergence result of the postprocessed solution,an asymptotically exact error estimator is constructed.The residual-type a posteriori error estimator is also proposed,and it's reliability and effectiveness are proved for variable coefficients time-harmonic Maxwell's equations.The efficiency and robustness of the proposed estimators is demonstrated by the numerical experiments.
基金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.
文摘This paper addresses tensile shock physics in thermoviscoelastic (TVE) solids without memory. The mathematical model is derived using conservation and balance laws (CBL) of classical continuum mechanics (CCM), incorporating the contravariant second Piola-Kirchhoff stress tensor, the covariant Green’s strain tensor, and its rates up to order n. This mathematical model permits the study of finite deformation and finite strain compressible deformation physics with an ordered rate dissipation mechanism. Constitutive theories are derived using conjugate pairs in entropy inequality and the representation theorem. The resulting mathematical model is both thermodynamically and mathematically consistent and has closure. The solution of the initial value problems (IVPs) describing evolutions is obtained using a variationally consistent space-time coupled finite element method, derived using space-time residual functional in which the local approximations are in hpk higher-order scalar product spaces. This permits accurate description problem physics over the discretization and also permits precise a posteriori computation of the space-time residual functional, an accurate measure of the accuracy of the computed solution. Model problem studies are presented to demonstrate tensile shock formation, propagation, reflection, and interaction. A unique feature of this research is that tensile shocks can only exist in solid matter, as their existence requires a medium to be elastic (presence of strain), which is only possible in a solid medium. In tensile shock physics, a decrease in the density of the medium caused by tensile waves leads to shock formation ahead of the wave. In contrast, in compressive shocks, an increase in density and the corresponding compressive waves result in the formation of compression shocks behind of the wave. Although these are two similar phenomena, they are inherently different in nature. To our knowledge, this work has not been reported in the published literature.
基金supported by the German National Science Foundation DFG within the Research Center MATHEON and,with the first author,by the WCU program through KOSEF(R31-2008-000-10049-0)support by the NSF grants DMS-0707602,DMS-0810176,DMS-0811153,DMS-0914788+2 种基金the German National Science Foundation within the Priority Programs SPP 1253,SPP 1506the German Federal Ministry for Education and Research(BMBF)within the projects BMBF-FROPT and BMBF-MeFreSimthe European Science Foundation(ESF)within the ESF Program OPTPDE.
文摘This paper aims at a general guideline to obtain a posteriori error estimates for the finite element error control in computational partial differential equations.In the abstract setting of mixed formulations,a generalised formulation of the corresponding residuals is proposed which then allows for the unified estimation of the respective dual norms.Notably,this can be done with an approach which is applicable in the same way to conforming,nonconforming and mixed discretisations.Subsequently,the unified approach is applied to various model problems.In particular,we consider the Laplace,Stokes,Navier-Lamé,and the semi-discrete eddy current equations.
基金This work is supported by Program for New Century Excellent Talents in University of China State Education Ministry NCET-04-0776, National Science Foundation of China, the National Basic Research Program under the Grant 2005CB321703, and the key project of China State Education Ministry and Hunan Education Commission.
文摘Four primal discontinuous Galerkin methods are applied to solve reactive transport problems, namely, Oden-BabuSka-Baumann DG (OBB-DG), non-symmetric interior penalty Galerkin (NIPG), symmetric interior penalty Galerkin (SIPG), and incomplete interior penalty Galerkin (IIPG). A unified a posteriori residual-type error estimation is derived explicitly for these methods. From the computed solution and given data, explicit estimators can be computed efficiently and directly, which can be used as error indicators for adaptation. Unlike in the reference [10], we obtain the error estimators in L^2 (L^2) norm by using duality techniques instead of in L^2(H^1) norm.
基金the National Basic Research Programthe National Natural Science Foundation of China(Grant No.2005CB321703)+2 种基金Scientific Research Fund of Hunan Provincial Education Departmentthe Outstanding Youth Scientist of the National Natural Science Foundation of China(Grant No.10625106)the National Basic Research Program of China(Grant No.2005CB321701)
文摘In this paper,we investigate the Legendre Galerkin spectral approximation of quadratic optimal control problems governed by parabolic equations.A spectral approximation scheme for the parabolic optimal control problem is presented.We obtain a posteriori error estimates of the approximated solutions for both the state and the control.
基金supported by DFG Research Center MATHEON"Mathematics for key technologies" in Berlinsupported by the NSFC under Grant 10601003 and A Foundation for the Author of National Excellent Doctoral Dissertation of PR China 200718support of two Sino-German workshops on Applied and Computational Mathematics held in 2005 and 2007 through the Sino-German office in Beijing.
文摘A unified a posteriori error analysis has been developed in [18, 21-23] to analyze the finite element error a posteriori under a universal roof. This paper contributes to the finite element meshes with hanging nodes which are required for local mesh-refining. The twodimensional 1-irregular triangulations into triangles and parallelograms and their combinations are considered with conforming and nonconforming finite element methods named after or by Courant, Q1, Crouzeix-Raviart, Poisson, Stokes and Navier-Lamé equations Han, Rannacher-Turek, and others for the The paper provides a unified a priori and a posteriori error analysis for triangulations with hanging nodes of degree ≤ 1 which are fundamental for local mesh refinement in self-adaptive finite element discretisations.
基金the National Natural Science Foundation of China(Grant No.10771150)the National Basic Research Program of China(Grant No.2005CB321701)the Program for New Century Excellent Talents in University(Grant No.NCET-07-0584)
文摘Two residual-based a posteriori error estimators of the nonconforming Crouzeix-Raviart element are derived for elliptic problems with Dirac delta source terms.One estimator is shown to be reliable and efficient,which yields global upper and lower bounds for the error in piecewise W1,p seminorm.The other one is proved to give a global upper bound of the error in Lp-norm.By taking the two estimators as refinement indicators,adaptive algorithms are suggested,which are experimentally shown to attain optimal convergence orders.
基金supported by National Natural Science Foundation of China (Grant No.10761003)Science and Technology Foundation of Guizhou Province of China (Grant No. [2011] 2111)
文摘This paper discusses conforming mixed finite element approximations for the Stokes eigenvalue problem. Firstly, several mixed finite element identities are proved. Based on these identities, the following new results are given: (1) It is proved that the numerical eigenvalues obtained by mini-element, P1-P1 element and Q1-Q1 element approximate the exact eigenvalues from above. (2) As for the P1-P1, Q1-Q1 and Q1-Po element eigenvalues, the asymptotically exact a posteriori error indicators are presented. (3) The reliable and efficient a posteriori error estimator proposed by Verfiirth is applied to mini-element eigenfunctions. Finally, numerical experiments are carried out to verify the theoretical analysis.
基金the NSF grants DMS-0410266 and DMS-0710831the China National Basic Research Program under the grant 2005CB321701+1 种基金the Program for the New Century Outstanding Talents in Universities of Chinathe Natural Science Foundation of Jiangsu Province under the grant BK2006511
文摘This paper develops a posteriori error estimates of residual type for conforming and mixed finite element approximations of the fourth order Cahn-Hilliard equation ut + △(ε△Au-ε^-1f(u)) = 0. It is shown that the a posteriori error bounds depends on ε^-1 only in some low polynomial order, instead of exponential order. Using these a posteriori error estimates, we construct at2 adaptive algorithm for computing the solution of the Cahn- Hilliard equation and its sharp interface limit, the Hele-Shaw flow. Numerical experiments are presented to show the robustness and effectiveness of the new error estimators and the proposed adaptive algorithm.
基金National Nature Science Foundation under Grants 60474027 and 10771211the National Basic Research Program under the Grant 2005CB321701
文摘In this paper, we discuss the a posteriori error estimate of the finite element approximation for the boundary control problems governed by the parabolic partial differential equations. Three different a posteriori error estimators are provided for the parabolic boundary control problems with the observations of the distributed state, the boundary state and the final state. It is proven that these estimators are reliable bounds of the finite element approximation errors, which can be used as the indicators of the mesh refinement in adaptive finite element methods.