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.展开更多
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.展开更多
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.展开更多
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.展开更多
A new technique of residual-type a posteriori error analysis is developed for the lowest- order Raviart-Thomas mixed finite element discretizations of convection-diffusion-reaction equations in two- or three-dimension...A new technique of residual-type a posteriori error analysis is developed for the lowest- order Raviart-Thomas mixed finite element discretizations of convection-diffusion-reaction equations in two- or three-dimension. Both centered mixed scheme and upwind-weighted mixed scheme are considered. The a posteriori error estimators, derived for the stress variable error plus scalar displacement error in L_2-norm, can be directly computed with the solutions of the mixed schemes without any additional cost, and are proven to be reliable. Local efficiency dependent on local variations in coefficients is obtained without any saturation assumption, and holds from the cases where convection or reaction is not present to convection- or reaction-dominated problems. The main tools of the analysis are the postprocessed approximation of scalar displacement, abstract error estimates, and the property of modified Oswald interpolation. Numerical experiments are carried out to support our theoretical results and to show the competitive behavior of the proposed posteriori error estimates.展开更多
In this paper, the superconvergence results are derived for a class of boundary control problems governed by Stokes equations. We derive superconvergence results for both the control and the state approximation. Base ...In this paper, the superconvergence results are derived for a class of boundary control problems governed by Stokes equations. We derive superconvergence results for both the control and the state approximation. Base on superconvergence results, we obtain asymptotically exact a posteriori error estimates.展开更多
We construct and analyze conservative local discontinuous Galerkin(LDG)methods for the Generalized Korteweg-de-Vries equation.LDG methods are designed by writing the equation as a system and performing separate approx...We construct and analyze conservative local discontinuous Galerkin(LDG)methods for the Generalized Korteweg-de-Vries equation.LDG methods are designed by writing the equation as a system and performing separate approximations to the spatial derivatives.The main focus is on the development of conservative methods which can preserve discrete versions of the first two invariants of the continuous solution,and a posteriori error estimates for a fully discrete approximation that is based on the idea of dispersive reconstruction.Numerical experiments are provided to verify the theoretical estimates.展开更多
In this paper,an initial boundary value problem of the space-time fractional diffusion equation is studied.Both temporal and spatial directions for this equation are discreted by the Galerkin spectral methods.And then...In this paper,an initial boundary value problem of the space-time fractional diffusion equation is studied.Both temporal and spatial directions for this equation are discreted by the Galerkin spectral methods.And then based on the discretization scheme,reliable a posteriori error estimates for the spectral approximation are derived.Some numerical examples are presented to verify the validity and applicability of the derived a posteriori error estimator.展开更多
In this paper, we study adaptive finite element discretization schemes for an optimal control problem governed by elliptic PDE with an integral constraint for the state. We derive the equivalent a posteriori error est...In this paper, we study adaptive finite element discretization schemes for an optimal control problem governed by elliptic PDE with an integral constraint for the state. We derive the equivalent a posteriori error estimator for the finite element approximation, which particularly suits adaptive multi-meshes to capture different singularities of the control and the state. Numerical examples are presented to demonstrate the efficiency of a posteriori error estimator and to confirm the theoretical results.展开更多
In this paper, we derive optimal order a posteriori error estimates for the local dis- continuous Galerkin (LDC) method for linear convection-diffusion problems in one space dimension. One of the key ingredients in ...In this paper, we derive optimal order a posteriori error estimates for the local dis- continuous Galerkin (LDC) method for linear convection-diffusion problems in one space dimension. One of the key ingredients in our analysis is the recent optimal superconver- gence result in [Y. Yang and C.-W. Shu, J. Comp. Math., 33 (2015), pp. 323-340]. We first prove that the LDG solution and its spatial derivative, respectively, converge in the L2-norm to (p + 1)-degree right and left Radau interpolating polynomials under mesh re- finement. The order of convergence is proved to be p + 2, when piecewise polynomials of degree at most p are used. These results are used to show that the leading error terms on each element for the solution and its derivative are proportional to (p + 1)-degree right and left Radau polynomials. We further prove that, for smooth solutions, the a posteriori LDG error estimates, which were constructed by the author in an earlier paper, converge, at a fixed time, to the true spatial errors in the L2-norm at (.9(hp+2) rate. Finally, we prove that the global effectivity indices in the L2-norm converge to unity at (9(h) rate. These results improve upon our previously published work in which the order of convergence for the a posteriori error estimates and the global effectivity index are proved to be p+3/2 and 1/2, respectively. Our proofs are valid for arbitrary regular meshes using PP polynomials with p ≥ 1. Several numerical experiments are performed to validate the theoretical results.展开更多
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.展开更多
We consider the problem of a posteriori error estimates and adaptivity for three typical force-based atomistic-to-continuum coupling methods.Combining the residual and the stability estimates,we derive computable a po...We consider the problem of a posteriori error estimates and adaptivity for three typical force-based atomistic-to-continuum coupling methods.Combining the residual and the stability estimates,we derive computable a posteriori error estima-tors for the three methods in the energy norm and formulate adaptive algorithms using these estimators.Our numerical experiments show optimal convergence rates of these algorithms.The efficiency of the estimators are also demonstrated numerically.展开更多
This paper is concerned with recovery type a posteriori error estimates of fully discrete finite element approximation for general convex parabolic optimal control problems with pointwise control constraints.The time ...This paper is concerned with recovery type a posteriori error estimates of fully discrete finite element approximation for general convex parabolic optimal control problems with pointwise control constraints.The time discretization is based on the backward Euler method.The state and the adjoint state are approximated by piecewise linear functions and the control is approximated by piecewise constant functions.We derive the superconvergence properties of finite element solutions.By using the superconvergence results,we obtain recovery type a posteriori error estimates.Some numerical examples are presented to verify the theoretical results.展开更多
This article is concerned with the numerical detection of bifurcation points of nonlinear partial differential equations as some parameter of interest is varied.In particular,we study in detail the numerical approxima...This article is concerned with the numerical detection of bifurcation points of nonlinear partial differential equations as some parameter of interest is varied.In particular,we study in detail the numerical approximation of the Bratu problem,based on exploiting the symmetric version of the interior penalty discontinuous Galerkin finite element method.A framework for a posteriori control of the discretization error in the computed critical parameter value is developed based upon the application of the dual weighted residual(DWR)approach.Numerical experiments are presented to highlight the practical performance of the proposed a posteriori error estimator.展开更多
Adaptive methods have been rapidly developed and applied in many fields of scientific and engi- neering computing. Reliable and efficient a posteriori error estimates play key roles for both adaptive finite element ...Adaptive methods have been rapidly developed and applied in many fields of scientific and engi- neering computing. Reliable and efficient a posteriori error estimates play key roles for both adaptive finite element and boundary element methods. The aim of this paper is to develop a posteriori error estimates for boundary element methods. The standard a posteriori error estimates for boundary element methods are obtained from the classical boundary integral equations. This paper presents hyper-singular a posteriori er- ror estimates based on the hyper-singular integral equations. Three kinds of residuals are used as the esti- mates for boundary element errors. The theoretical analysis and numerical examples show that the hyper- singular residuals are good a posteriori error indicators in many adaptive boundary element computations.展开更多
文摘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 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.
基金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 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.
基金The authors are grateful for the anonymous referees for their helpful com- ments. This work was supported in part by The Education Science Foundation of Chongqing (KJ120420), National Natural Science Foundation of China (11171239), The Project-sponsored by Scientific Research Foundation for the Returned Overseas Chinese Scholars and Open Fund of Key Laboratory of Mountain Hazards and Earth Surface Processes, CAS.
文摘A new technique of residual-type a posteriori error analysis is developed for the lowest- order Raviart-Thomas mixed finite element discretizations of convection-diffusion-reaction equations in two- or three-dimension. Both centered mixed scheme and upwind-weighted mixed scheme are considered. The a posteriori error estimators, derived for the stress variable error plus scalar displacement error in L_2-norm, can be directly computed with the solutions of the mixed schemes without any additional cost, and are proven to be reliable. Local efficiency dependent on local variations in coefficients is obtained without any saturation assumption, and holds from the cases where convection or reaction is not present to convection- or reaction-dominated problems. The main tools of the analysis are the postprocessed approximation of scalar displacement, abstract error estimates, and the property of modified Oswald interpolation. Numerical experiments are carried out to support our theoretical results and to show the competitive behavior of the proposed posteriori error estimates.
基金The research was supported by the Special Funds for Major State Basic Research Projects (No. G2000067102), National Natural Science Foundation of China (No. 60474027).
文摘In this paper, the superconvergence results are derived for a class of boundary control problems governed by Stokes equations. We derive superconvergence results for both the control and the state approximation. Base on superconvergence results, we obtain asymptotically exact a posteriori error estimates.
基金The research of O.Karakashian was partially supported by National Science Foundation grant DMS-1216740The research of Y.Xing was partially supported by National Science Foundation grants DMS-1216454 and DMS-1621111.
文摘We construct and analyze conservative local discontinuous Galerkin(LDG)methods for the Generalized Korteweg-de-Vries equation.LDG methods are designed by writing the equation as a system and performing separate approximations to the spatial derivatives.The main focus is on the development of conservative methods which can preserve discrete versions of the first two invariants of the continuous solution,and a posteriori error estimates for a fully discrete approximation that is based on the idea of dispersive reconstruction.Numerical experiments are provided to verify the theoretical estimates.
基金supported by the State Key Program of National Natural Science Foundation of China(No.11931003)National Natural Science Foundation of China(Nos.41974133,11671157 and 11971410)supported by the Innovation Project of Graduate School of South China Normal University(No.2018LKXM008).
文摘In this paper,an initial boundary value problem of the space-time fractional diffusion equation is studied.Both temporal and spatial directions for this equation are discreted by the Galerkin spectral methods.And then based on the discretization scheme,reliable a posteriori error estimates for the spectral approximation are derived.Some numerical examples are presented to verify the validity and applicability of the derived a posteriori error estimator.
基金supported by the National Basic Research Program of P.R.China under the grant 2005CB321703the NSFC under the grants:10441005 and 10571108the Research Fund for Doctoral Program of High Education by China State Education Ministry
文摘In this paper, we study adaptive finite element discretization schemes for an optimal control problem governed by elliptic PDE with an integral constraint for the state. We derive the equivalent a posteriori error estimator for the finite element approximation, which particularly suits adaptive multi-meshes to capture different singularities of the control and the state. Numerical examples are presented to demonstrate the efficiency of a posteriori error estimator and to confirm the theoretical results.
文摘In this paper, we derive optimal order a posteriori error estimates for the local dis- continuous Galerkin (LDC) method for linear convection-diffusion problems in one space dimension. One of the key ingredients in our analysis is the recent optimal superconver- gence result in [Y. Yang and C.-W. Shu, J. Comp. Math., 33 (2015), pp. 323-340]. We first prove that the LDG solution and its spatial derivative, respectively, converge in the L2-norm to (p + 1)-degree right and left Radau interpolating polynomials under mesh re- finement. The order of convergence is proved to be p + 2, when piecewise polynomials of degree at most p are used. These results are used to show that the leading error terms on each element for the solution and its derivative are proportional to (p + 1)-degree right and left Radau polynomials. We further prove that, for smooth solutions, the a posteriori LDG error estimates, which were constructed by the author in an earlier paper, converge, at a fixed time, to the true spatial errors in the L2-norm at (.9(hp+2) rate. Finally, we prove that the global effectivity indices in the L2-norm converge to unity at (9(h) rate. These results improve upon our previously published work in which the order of convergence for the a posteriori error estimates and the global effectivity index are proved to be p+3/2 and 1/2, respectively. Our proofs are valid for arbitrary regular meshes using PP polynomials with p ≥ 1. Several numerical experiments are performed to validate the theoretical results.
基金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.
基金Hao Wang was partially supported by NSFC grant 11501389,11471214 and Sichuan University Starting Up Research Funding No.2082204194117Shaohui Li-u was partially supported by Sichuan University National Students’Platform for Innovation and Entrepreneurship Training Program No.201610610172Stony Brook University PhD Scholarship.
文摘We consider the problem of a posteriori error estimates and adaptivity for three typical force-based atomistic-to-continuum coupling methods.Combining the residual and the stability estimates,we derive computable a posteriori error estima-tors for the three methods in the energy norm and formulate adaptive algorithms using these estimators.Our numerical experiments show optimal convergence rates of these algorithms.The efficiency of the estimators are also demonstrated numerically.
基金supported by Guangdong Province Universities and Colleges Pearl River Scholar Funded Scheme(2008)National Science Foundation of China(10971074)+1 种基金Specialized Research Fund for the Doctoral Program of Higher Education(20114407110009)Hunan Provinical Innovation Foundation for Postgraduate(lx2009 B120)。
文摘This paper is concerned with recovery type a posteriori error estimates of fully discrete finite element approximation for general convex parabolic optimal control problems with pointwise control constraints.The time discretization is based on the backward Euler method.The state and the adjoint state are approximated by piecewise linear functions and the control is approximated by piecewise constant functions.We derive the superconvergence properties of finite element solutions.By using the superconvergence results,we obtain recovery type a posteriori error estimates.Some numerical examples are presented to verify the theoretical results.
基金the financial support of the EPSRC under the grant EP/E013724the support of the EPSRC under the grant EP/F01340X.
文摘This article is concerned with the numerical detection of bifurcation points of nonlinear partial differential equations as some parameter of interest is varied.In particular,we study in detail the numerical approximation of the Bratu problem,based on exploiting the symmetric version of the interior penalty discontinuous Galerkin finite element method.A framework for a posteriori control of the discretization error in the computed critical parameter value is developed based upon the application of the dual weighted residual(DWR)approach.Numerical experiments are presented to highlight the practical performance of the proposed a posteriori error estimator.
基金Supported by the National Key Basic Research and Development(973) Program of China (No. G19990328) and the Knowledge Innovation Program of the Chinese Academy of Sciences
文摘Adaptive methods have been rapidly developed and applied in many fields of scientific and engi- neering computing. Reliable and efficient a posteriori error estimates play key roles for both adaptive finite element and boundary element methods. The aim of this paper is to develop a posteriori error estimates for boundary element methods. The standard a posteriori error estimates for boundary element methods are obtained from the classical boundary integral equations. This paper presents hyper-singular a posteriori er- ror estimates based on the hyper-singular integral equations. Three kinds of residuals are used as the esti- mates for boundary element errors. The theoretical analysis and numerical examples show that the hyper- singular residuals are good a posteriori error indicators in many adaptive boundary element computations.