Error estimates are established for the finite element methods to solve a class of second order nonlinear parabolic equations. Optimal rates of convergence in L 2 and H 1 norms are derived.Meanwhile,the schemes are se...Error estimates are established for the finite element methods to solve a class of second order nonlinear parabolic equations. Optimal rates of convergence in L 2 and H 1 norms are derived.Meanwhile,the schemes are second order correct in time.展开更多
The effect of numerical integration in finite element methods applied to a class of nonlinear parabolic equations is considered and some sufficient conditions on the quadrature scheme to ensure that the order of conve...The effect of numerical integration in finite element methods applied to a class of nonlinear parabolic equations is considered and some sufficient conditions on the quadrature scheme to ensure that the order of convergence is unaltered in the presence of numerical integration are given. Optimal L 2 and H 1 estimates for the error and its time derivative are established.展开更多
This paper presents alternating direction finite volume element methods for three-dimensional parabolic partial differential equations and gives four computational schemes, one is analogous to Douglas finite differenc...This paper presents alternating direction finite volume element methods for three-dimensional parabolic partial differential equations and gives four computational schemes, one is analogous to Douglas finite difference scheme with second-order splitting error, the other two schemes have third-order splitting error, and the last one is an extended LOD scheme. The L2 norm and H1 semi-norm error estimates are obtained for the first scheme and second one, respectively. Finally, two numerical examples are provided to illustrate the efficiency and accuracy of the methods.展开更多
In this paper we introduce two kinds of parallel Schwarz domain decomposition me thods for general, selfadjoint, second order parabolic equations and study the dependence of their convergence rates on parameters of ti...In this paper we introduce two kinds of parallel Schwarz domain decomposition me thods for general, selfadjoint, second order parabolic equations and study the dependence of their convergence rates on parameters of time-step and space-mesh. We prove that the, approximate solution has convergence independent of iteration times at each time-level. And the L^2 error estimates are given.展开更多
In this work, system of parabolic equations with discontinuous coefficients is studied. The domain decomposition method modified by a characteristic finite element procedure is applied. A function is defined to approx...In this work, system of parabolic equations with discontinuous coefficients is studied. The domain decomposition method modified by a characteristic finite element procedure is applied. A function is defined to approximate the fluxes on inner boundaries by using the solution at the previous level. Thus the parallelism is achieved. Convergence analysis and error estimate are also presented.展开更多
In this paper we investigate the existence, uniqueness and regularity of the solution of semilinear parabolic equations with coefficients that are discontinuous across the interface, some prior estimates are obtained....In this paper we investigate the existence, uniqueness and regularity of the solution of semilinear parabolic equations with coefficients that are discontinuous across the interface, some prior estimates are obtained. A net shape of the finite elements around the singular points was designed in [7] to solve the linear elliptic problems, by means of that net, we prove that the approximate solution has the same convergence rate as that without singularity.展开更多
In this paper,we study the initial-boundary value problem for the semilinear parabolic equations ut-△Xu=|u|p-1u,where X=(X1,X2,…,Xm) is a system of real smooth vector fields which satisfy the H?rmander’s conditio...In this paper,we study the initial-boundary value problem for the semilinear parabolic equations ut-△Xu=|u|p-1u,where X=(X1,X2,…,Xm) is a system of real smooth vector fields which satisfy the H?rmander’s condition,and △X=∑j=1m Xj2 is a finitely degenerate elliptic operator.Using potential well method,we first prove the invariance of some sets and vacuum isolating of solutions.Finally,by the Galerkin method and the concavity method we show the global existence and blow-up in finite time of solutions with low initial energy or critical initial energy,and also we discuss the asymptotic behavior of the global solutions.展开更多
In this paper, we present the a posteriori error estimate of two-grid mixed finite element methods by averaging techniques for semilinear elliptic equations. We first propose the two-grid algorithms to linearize the m...In this paper, we present the a posteriori error estimate of two-grid mixed finite element methods by averaging techniques for semilinear elliptic equations. We first propose the two-grid algorithms to linearize the mixed method equations. Then, the averaging technique is used to construct the a posteriori error estimates of the two-grid mixed finite element method and theoretical analysis are given for the error estimators. Finally, we give some numerical examples to verify the reliability and efficiency of the a posteriori error estimator.展开更多
A nonlinear parabolic system is derived to describe incompressible nuclear waste-disposal contamination in porous media. A sequential implicit tirne-stepping is defined, in which the pressure and Darcy velocity of the...A nonlinear parabolic system is derived to describe incompressible nuclear waste-disposal contamination in porous media. A sequential implicit tirne-stepping is defined, in which the pressure and Darcy velocity of the mixture are approximated simultaneously by a mixed finite element method and the brine, radionuclid and heat are treated by a combination of a Galerkin finite element method and the method of characteristics. Optimal-order convergence in L2 is proved. Time-truncation errors of standard procedures are reduced by time stepping along the characteristics of the hyperbolic part of the brine, radionuclide and heal equalios, temporal and spatial error are lossened by direct compulation of the velocity in the mixed method, as opposed to differentiation of the pressure.展开更多
In this paper, we extend the applications of proper orthogonal decomposition (POD) method, i.e., apply POD method to a mixed finite element (MFE) formulation naturally satisfied Brezz-Babu^ka for parabolic equatio...In this paper, we extend the applications of proper orthogonal decomposition (POD) method, i.e., apply POD method to a mixed finite element (MFE) formulation naturally satisfied Brezz-Babu^ka for parabolic equations, establish a reduced-order MFE formulation with lower dimensions and sufficiently high accuracy, and provide the error estimates between the reduced-order POD MFE solutions and the classical MFE solutions and the implementation of algorithm for solving reduced-order MFE formulation. Some numerical examples illustrate the fact that the results of numerical computation are consis- tent with theoretical conclusions. Moreover, it is shown that the new reduced-order MFE formulation based on POD method is feasible and efficient for solving MFE formulation for parabolic equations.展开更多
We study spatially semidiscrete and fully discrete two-scale composite nite element method for approximations of the nonlinear parabolic equations with homogeneous Dirichlet boundary conditions in a convex polygonal d...We study spatially semidiscrete and fully discrete two-scale composite nite element method for approximations of the nonlinear parabolic equations with homogeneous Dirichlet boundary conditions in a convex polygonal domain in the plane.This new class of nite elements,which is called composite nite elements,was rst introduced by Hackbusch and Sauter[Numer.Math.,75(1997),pp.447-472]for the approximation of partial di erential equations on domains with complicated geometry.The aim of this paper is to introduce an effcient numerical method which gives a lower dimensional approach for solving partial di erential equations by domain discretization method.The composite nite element method introduces two-scale grid for discretization of the domain,the coarse-scale and the ne-scale grid with the degrees of freedom lies on the coarse-scale grid only.While the ne-scale grid is used to resolve the Dirichlet boundary condition,the dimension of the nite element space depends only on the coarse-scale grid.As a consequence,the resulting linear system will have a fewer number of unknowns.A continuous,piecewise linear composite nite element space is employed for the space discretization whereas the time discretization is based on both the backward Euler and the Crank-Nicolson methods.We have derived the error estimates in the L^(∞)(L^(2))-norm for both semidiscrete and fully discrete schemes.Moreover,numerical simulations show that the proposed method is an efficient method to provide a good approximate solution.展开更多
The numerical approximations of the dynamical systems governed by semilinear parabolic equations are considered. An abstract framework for long time error estimates is established. When applied to reaction diffusion...The numerical approximations of the dynamical systems governed by semilinear parabolic equations are considered. An abstract framework for long time error estimates is established. When applied to reaction diffusion equation, Navier Stokes equations and Chan Hilliard equation, approximated by Galerkin and nonlinear Galerkin methods in space and by Runge Kutta method in time, our framework yields error estimates uniform in time.展开更多
This paper deals with a-posteriori error estimates for piecewise linear finite element approximations of parabolic problems in two space dimensions. The analysis extends previous results for elliptic problems to the p...This paper deals with a-posteriori error estimates for piecewise linear finite element approximations of parabolic problems in two space dimensions. The analysis extends previous results for elliptic problems to the parabolic context.展开更多
Poisson-Nernst-Planck equations are a coupled system of nonlinear partial differential equations consisting of the Nernst-Planck equation and the electrostatic Poisson equation with delta distribution sources,which de...Poisson-Nernst-Planck equations are a coupled system of nonlinear partial differential equations consisting of the Nernst-Planck equation and the electrostatic Poisson equation with delta distribution sources,which describe the electrodiffusion of ions in a solvated biomolecular system.In this paper,some error bounds for a piecewise finite element approximation to this problem are derived.Several numerical examples including biomolecular problems are shown to support our analysis.展开更多
This paper proposes a reliable and efficient a posteriori error estimator for the finite element methods for the beam problem. It is proved that the error can be bounded by the computable error estimator from above an...This paper proposes a reliable and efficient a posteriori error estimator for the finite element methods for the beam problem. It is proved that the error can be bounded by the computable error estimator from above and below up to multiplicative constants that do neither depend on the meshsize nor on the thickness of the beam.展开更多
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.展开更多
A proper orthogonal decomposition (POD) method is applied to a usual finite element (FE) formulation for parabolic equations so that it is reduced into a POD FE formulation with lower dimensions and enough high accura...A proper orthogonal decomposition (POD) method is applied to a usual finite element (FE) formulation for parabolic equations so that it is reduced into a POD FE formulation with lower dimensions and enough high accuracy. The errors between the reduced POD FE solution and the usual FE solution are analyzed. It is shown by numerical examples that the results of numerical computations are consistent with theoretical conclusions. Moreover, it is also shown that this validates the feasibility and efficiency of POD method.展开更多
We study an indirect finite element approximation for two-sided space-fractional diffusion equations in one space dimension.By the representation formula of the solutions u(x)to the proposed variable coefficient model...We study an indirect finite element approximation for two-sided space-fractional diffusion equations in one space dimension.By the representation formula of the solutions u(x)to the proposed variable coefficient models in terms of v(x),the solutions to the constant coefficient analogues,we apply finite element methods for the constant coefficient fractional diffusion equations to solve for the approximations vh(x)to v(x)and then obtain the approximations uh(x)of u(x)by plugging vh(x)into the representation of u(x).Optimal-order convergence estimates of u(x)−uh(x)are proved in both L2 and Hα∕2 norms.Several numerical experiments are presented to demonstrate the sharpness of the derived error estimates.展开更多
Adaptive space-time finite element method, continuous in space but discontinuous in time for semi-linear parabolic problems is discussed. The approach is based on a combination of finite element and finite difference ...Adaptive space-time finite element method, continuous in space but discontinuous in time for semi-linear parabolic problems is discussed. The approach is based on a combination of finite element and finite difference techniques. The existence and uniqueness of the weak solution are proved without any assumptions on choice of the spacetime meshes. Basic error estimates in L-infinity (L-2) norm, that is maximum-norm in time, L-2-norm in space are obtained. The numerical results are given in the last part and the analysis between theoretic and experimental results are obtained.展开更多
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.展开更多
基金National Natural Science Foundation and Doctoral Foundation of Education Ministry of Stat
文摘Error estimates are established for the finite element methods to solve a class of second order nonlinear parabolic equations. Optimal rates of convergence in L 2 and H 1 norms are derived.Meanwhile,the schemes are second order correct in time.
文摘The effect of numerical integration in finite element methods applied to a class of nonlinear parabolic equations is considered and some sufficient conditions on the quadrature scheme to ensure that the order of convergence is unaltered in the presence of numerical integration are given. Optimal L 2 and H 1 estimates for the error and its time derivative are established.
文摘This paper presents alternating direction finite volume element methods for three-dimensional parabolic partial differential equations and gives four computational schemes, one is analogous to Douglas finite difference scheme with second-order splitting error, the other two schemes have third-order splitting error, and the last one is an extended LOD scheme. The L2 norm and H1 semi-norm error estimates are obtained for the first scheme and second one, respectively. Finally, two numerical examples are provided to illustrate the efficiency and accuracy of the methods.
基金This work was supported by Natural Science Foundation of China and Shandong Province.
文摘In this paper we introduce two kinds of parallel Schwarz domain decomposition me thods for general, selfadjoint, second order parabolic equations and study the dependence of their convergence rates on parameters of time-step and space-mesh. We prove that the, approximate solution has convergence independent of iteration times at each time-level. And the L^2 error estimates are given.
基金This work was supported by the Youth Development Foundation of Shandong University at Weihai (Grant No. Z200607).
文摘In this work, system of parabolic equations with discontinuous coefficients is studied. The domain decomposition method modified by a characteristic finite element procedure is applied. A function is defined to approximate the fluxes on inner boundaries by using the solution at the previous level. Thus the parallelism is achieved. Convergence analysis and error estimate are also presented.
文摘In this paper we investigate the existence, uniqueness and regularity of the solution of semilinear parabolic equations with coefficients that are discontinuous across the interface, some prior estimates are obtained. A net shape of the finite elements around the singular points was designed in [7] to solve the linear elliptic problems, by means of that net, we prove that the approximate solution has the same convergence rate as that without singularity.
基金supported by National Natural Science Foundation of China(11631011 and 11626251)
文摘In this paper,we study the initial-boundary value problem for the semilinear parabolic equations ut-△Xu=|u|p-1u,where X=(X1,X2,…,Xm) is a system of real smooth vector fields which satisfy the H?rmander’s condition,and △X=∑j=1m Xj2 is a finitely degenerate elliptic operator.Using potential well method,we first prove the invariance of some sets and vacuum isolating of solutions.Finally,by the Galerkin method and the concavity method we show the global existence and blow-up in finite time of solutions with low initial energy or critical initial energy,and also we discuss the asymptotic behavior of the global solutions.
文摘In this paper, we present the a posteriori error estimate of two-grid mixed finite element methods by averaging techniques for semilinear elliptic equations. We first propose the two-grid algorithms to linearize the mixed method equations. Then, the averaging technique is used to construct the a posteriori error estimates of the two-grid mixed finite element method and theoretical analysis are given for the error estimators. Finally, we give some numerical examples to verify the reliability and efficiency of the a posteriori error estimator.
基金The research was supported by the Natural Science Foundation of China
文摘A nonlinear parabolic system is derived to describe incompressible nuclear waste-disposal contamination in porous media. A sequential implicit tirne-stepping is defined, in which the pressure and Darcy velocity of the mixture are approximated simultaneously by a mixed finite element method and the brine, radionuclid and heat are treated by a combination of a Galerkin finite element method and the method of characteristics. Optimal-order convergence in L2 is proved. Time-truncation errors of standard procedures are reduced by time stepping along the characteristics of the hyperbolic part of the brine, radionuclide and heal equalios, temporal and spatial error are lossened by direct compulation of the velocity in the mixed method, as opposed to differentiation of the pressure.
基金supported by the National Science Foundation of China(11271127 and 11061009)Science Research Program of Guizhou(GJ[2011]2367)the Co-Construction Project of Beijing Municipal Commission of Education
文摘In this paper, we extend the applications of proper orthogonal decomposition (POD) method, i.e., apply POD method to a mixed finite element (MFE) formulation naturally satisfied Brezz-Babu^ka for parabolic equations, establish a reduced-order MFE formulation with lower dimensions and sufficiently high accuracy, and provide the error estimates between the reduced-order POD MFE solutions and the classical MFE solutions and the implementation of algorithm for solving reduced-order MFE formulation. Some numerical examples illustrate the fact that the results of numerical computation are consis- tent with theoretical conclusions. Moreover, it is shown that the new reduced-order MFE formulation based on POD method is feasible and efficient for solving MFE formulation for parabolic equations.
基金The author gratefully acknowledges valuable support provided by the Department of Mathematics,NIT Calicut and the DST,Government of India,for providing support to carry out this work under the scheme TIST(No.SR/FST/MS-I/2019/40).
文摘We study spatially semidiscrete and fully discrete two-scale composite nite element method for approximations of the nonlinear parabolic equations with homogeneous Dirichlet boundary conditions in a convex polygonal domain in the plane.This new class of nite elements,which is called composite nite elements,was rst introduced by Hackbusch and Sauter[Numer.Math.,75(1997),pp.447-472]for the approximation of partial di erential equations on domains with complicated geometry.The aim of this paper is to introduce an effcient numerical method which gives a lower dimensional approach for solving partial di erential equations by domain discretization method.The composite nite element method introduces two-scale grid for discretization of the domain,the coarse-scale and the ne-scale grid with the degrees of freedom lies on the coarse-scale grid only.While the ne-scale grid is used to resolve the Dirichlet boundary condition,the dimension of the nite element space depends only on the coarse-scale grid.As a consequence,the resulting linear system will have a fewer number of unknowns.A continuous,piecewise linear composite nite element space is employed for the space discretization whereas the time discretization is based on both the backward Euler and the Crank-Nicolson methods.We have derived the error estimates in the L^(∞)(L^(2))-norm for both semidiscrete and fully discrete schemes.Moreover,numerical simulations show that the proposed method is an efficient method to provide a good approximate solution.
文摘The numerical approximations of the dynamical systems governed by semilinear parabolic equations are considered. An abstract framework for long time error estimates is established. When applied to reaction diffusion equation, Navier Stokes equations and Chan Hilliard equation, approximated by Galerkin and nonlinear Galerkin methods in space and by Runge Kutta method in time, our framework yields error estimates uniform in time.
文摘This paper deals with a-posteriori error estimates for piecewise linear finite element approximations of parabolic problems in two space dimensions. The analysis extends previous results for elliptic problems to the parabolic context.
基金supported by the China NSF(NSFC 11001062,NSFC 11161014)the fund from Education Department of Guangxi Province under grant 201012MS094B.Z.Lu was supported by the National Center for Mathematics and Interdisciplinary Sciences,Chinese Academy of Sciences and the China NSF(NSFC10971218).
文摘Poisson-Nernst-Planck equations are a coupled system of nonlinear partial differential equations consisting of the Nernst-Planck equation and the electrostatic Poisson equation with delta distribution sources,which describe the electrodiffusion of ions in a solvated biomolecular system.In this paper,some error bounds for a piecewise finite element approximation to this problem are derived.Several numerical examples including biomolecular problems are shown to support our analysis.
基金the fortieth postdoctoral foundation of Chinathe NSFC under Grant 10601003+1 种基金the Scientific Research Foundation for the Returned Overseas Chinese Scholars,State Education MinistryA Foundation for the Author of National Excellent Doctoral Dissertation of PRC 200718
文摘This paper proposes a reliable and efficient a posteriori error estimator for the finite element methods for the beam problem. It is proved that the error can be bounded by the computable error estimator from above and below up to multiplicative constants that do neither depend on the meshsize nor on the thickness of the beam.
基金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.
基金supported by National Natural Science Foundation of China (Grant Nos. 10871022, 10771065,and 60573158)Natural Science Foundation of Hebei Province (Grant No. A2007001027)
文摘A proper orthogonal decomposition (POD) method is applied to a usual finite element (FE) formulation for parabolic equations so that it is reduced into a POD FE formulation with lower dimensions and enough high accuracy. The errors between the reduced POD FE solution and the usual FE solution are analyzed. It is shown by numerical examples that the results of numerical computations are consistent with theoretical conclusions. Moreover, it is also shown that this validates the feasibility and efficiency of POD method.
基金the OSD/ARO MURI Grant W911NF-15-1-0562the National Science Foundation under Grant DMS-1620194.
文摘We study an indirect finite element approximation for two-sided space-fractional diffusion equations in one space dimension.By the representation formula of the solutions u(x)to the proposed variable coefficient models in terms of v(x),the solutions to the constant coefficient analogues,we apply finite element methods for the constant coefficient fractional diffusion equations to solve for the approximations vh(x)to v(x)and then obtain the approximations uh(x)of u(x)by plugging vh(x)into the representation of u(x).Optimal-order convergence estimates of u(x)−uh(x)are proved in both L2 and Hα∕2 norms.Several numerical experiments are presented to demonstrate the sharpness of the derived error estimates.
文摘Adaptive space-time finite element method, continuous in space but discontinuous in time for semi-linear parabolic problems is discussed. The approach is based on a combination of finite element and finite difference techniques. The existence and uniqueness of the weak solution are proved without any assumptions on choice of the spacetime meshes. Basic error estimates in L-infinity (L-2) norm, that is maximum-norm in time, L-2-norm in space are obtained. The numerical results are given in the last part and the analysis between theoretic and experimental results are obtained.
基金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.