An H^1-Galerkin mixed finite element method is discussed for a class of second order SchrSdinger equation. Optimal error estimates of semidiscrete schemes are derived for problems in one space dimension. At the same t...An H^1-Galerkin mixed finite element method is discussed for a class of second order SchrSdinger equation. Optimal error estimates of semidiscrete schemes are derived for problems in one space dimension. At the same time, optimal error estimates are derived for fully discrete schemes. And it is showed that the H1-Galerkin mixed finite element approximations have the same rate of convergence as in the classical mixed finite element methods without requiring the LBB consistency condition.展开更多
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.展开更多
<div style="text-align:justify;"> In this paper, we study the error estimates for direct discontinuous Galerkin methods based on the upwind-biased fluxes. We use a newly global projection to obtain the...<div style="text-align:justify;"> In this paper, we study the error estimates for direct discontinuous Galerkin methods based on the upwind-biased fluxes. We use a newly global projection to obtain the optimal error estimates. The numerical experiments imply that <em>L</em><sup>2 </sup>norms error estimates can reach to order <em>k</em> + 1 by using time discretization methods. </div>展开更多
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.展开更多
Allen and Liu (1995) introduced a new method for a time-dependent convection dominated diffusion problem, which combines the modified method of characteristics and method of streamline diffusion. But they ignored the ...Allen and Liu (1995) introduced a new method for a time-dependent convection dominated diffusion problem, which combines the modified method of characteristics and method of streamline diffusion. But they ignored the fact that the accuracy of time discretization decays at half an order when the characteristic line goes out of the domain. In present paper, the author shows that, as a remedy, a simple lumped scheme yields a full accuracy approximation. Forthermore, some local error bounds independent of the small viscosity axe derived for this scheme outside the boundary layers.展开更多
In the present paper, a new numerical method for solving initial-boundary value problems of evolutionary equations is proposed and studied, combining difference method with high accuracy with boundary integral equatio...In the present paper, a new numerical method for solving initial-boundary value problems of evolutionary equations is proposed and studied, combining difference method with high accuracy with boundary integral equation method. The numerical approximate schemes for both problems on a bounded or unbounded domain in R3 are proposed and their prior error estimates are obtained.展开更多
In this paper,the uniform error estimates with respect to t∈[0, ∞ ) of the nonlinear Galerkin method are given for the long time integration of the Kuramoto-Sivashinsky equation. The nonlinear Galerkin method is use...In this paper,the uniform error estimates with respect to t∈[0, ∞ ) of the nonlinear Galerkin method are given for the long time integration of the Kuramoto-Sivashinsky equation. The nonlinear Galerkin method is used to study the asymptotic behaviour of Kuramoto-Sivashinsky equation and to construct the bifurcation diagrams.展开更多
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 purpose of this paper is to investigate the convergence of the mixed finite element method for the initial-boundary value problem for the Sobolev equation Ut-div{aut + b1 u} = f based on the Raviart-Thomas space ...The purpose of this paper is to investigate the convergence of the mixed finite element method for the initial-boundary value problem for the Sobolev equation Ut-div{aut + b1 u} = f based on the Raviart-Thomas space Vh × Wh H(div; × L2(). Optimal order estimates are obtained for the approximation of u, ut, the associated velocity p and divp respectively in L(0,T;L2()), L(0,T;L2()), L(0,T;L2()2), and L2(0, T; L2()). Quasi-optimal order estimates are obtained for the approximations of u, ut in L(0, T; L()) and p in L(0,T; L()2).展开更多
Although the two-grid finite element decoupled scheme for mixed Navier-Stokes/ Darcy model in literatures has given the numerical results of optimal convergence order, the theoretical analysis only obtain the optimal ...Although the two-grid finite element decoupled scheme for mixed Navier-Stokes/ Darcy model in literatures has given the numerical results of optimal convergence order, the theoretical analysis only obtain the optimal error order for the porous media flow and a non-optimal error order for the fluid flow. In this article, we give a more rigorous of the error analysis for the fluid flow and obtain the optimal error estimates of the velocity and the pressure.展开更多
A Fourier spectral method for the generalized Korteweg-de Vries equation with periodic boundary conditions is analyzed, and a corresponding optimal error estimate in L^2-norm is obtained. It improves the result presen...A Fourier spectral method for the generalized Korteweg-de Vries equation with periodic boundary conditions is analyzed, and a corresponding optimal error estimate in L^2-norm is obtained. It improves the result presented by Maday and Quarteroni. A modified Fourier pseudospectral method is also presented, with the same convergence properties as the Fourier spectral method.展开更多
In this paper, a spectral method to analyze the generalized Benjamin Bona Mahony equations is used. The existence and uniqueness of global smooth solution of these equations are proved. The large time error estimati...In this paper, a spectral method to analyze the generalized Benjamin Bona Mahony equations is used. The existence and uniqueness of global smooth solution of these equations are proved. The large time error estimation between the spectral approximate solution and the exact solution is obtained.展开更多
The proper orthogonal decomposition (POD) is a model reduction technique for the simulation Of physical processes governed by partial differential equations (e.g., fluid flows). It has been successfully used in th...The proper orthogonal decomposition (POD) is a model reduction technique for the simulation Of physical processes governed by partial differential equations (e.g., fluid flows). It has been successfully used in the reduced-order modeling of complex systems. In this paper, the applications of the POD method are extended, i.e., the POD method is applied to a classical finite difference (FD) scheme for the non-stationary Stokes equation with a real practical applied background. A reduced FD scheme is established with lower dimensions and sufficiently high accuracy, and the error estimates are provided between the reduced and the classical FD solutions. Some numerical examples illustrate that the numerical results are consistent with theoretical conclusions. Moreover, it is shown that the reduced FD scheme based on the POD method is feasible and efficient in solving the FD scheme for the non-stationary Stokes equation.展开更多
In this paper,we prove the optimal error estimates in L2 norm of the semidiscrete local discontinuous Galerkin(LDG)method for the thin film epitaxy problem without slope selection.To relax the severe time step restric...In this paper,we prove the optimal error estimates in L2 norm of the semidiscrete local discontinuous Galerkin(LDG)method for the thin film epitaxy problem without slope selection.To relax the severe time step restriction of explicit time marching methods,we employ a class of exponential time differencing(ETD)schemes for time integration,which is based on a linear convex splitting principle.Numerical experiments of the accuracy and long time simulations are given to show the efficiency and capability of the proposed numerical schemes.展开更多
In this paper we deal with the convergence analysis of the finite element method for an elliptic penalized unilateral obstacle optimal control problem where the control and the obstacle coincide.Error estimates are es...In this paper we deal with the convergence analysis of the finite element method for an elliptic penalized unilateral obstacle optimal control problem where the control and the obstacle coincide.Error estimates are established for both state and control variables.We apply a fixed point type iteration method to solve the discretized problem.To corroborate our error estimations and the eficiency of our algorithms,the convergence results and numerical experiments are illustrated by concrete examples.展开更多
The fractional optimal control problem leads to significantly increased computational complexity compared to the corresponding classical integer-order optimal control problem,due to the global properties of fractional...The fractional optimal control problem leads to significantly increased computational complexity compared to the corresponding classical integer-order optimal control problem,due to the global properties of fractional differential operators.In this paper,we focus on an optimal control problem governed by fractional differential equations with an integral constraint on the state variable.By the proposed first-order optimality condition consisting of a Lagrange multiplier,we design a spectral Galerkin discrete scheme with weighted orthogonal Jacobi polynomials to approximate the resulting state and adjoint state equations.Furthermore,a priori error estimates for state,adjoint state and control variables are discussed in details.Illustrative numerical tests are given to demonstrate the validity and applicability of our proposed approximations and theoretical results.展开更多
In this paper, we propose the nonconforming virtual element method (NCVEM) discretization for the pointwise control constraint optimal control problem governed by elliptic equations. Based on the NCVEM approximation o...In this paper, we propose the nonconforming virtual element method (NCVEM) discretization for the pointwise control constraint optimal control problem governed by elliptic equations. Based on the NCVEM approximation of state equation and the variational discretization of control variables, we construct a virtual element discrete scheme. For the state, adjoint state and control variable, we obtain the corresponding prior estimate in H<sup>1</sup> and L<sup>2</sup> norms. Finally, some numerical experiments are carried out to support the theoretical results.展开更多
This paper is devoted to the error estimate for the iterative discontinuous Galerkin(IDG)method introduced in[P.Yin,Y.Huang and H.Liu.Commun.Comput.Phys.16:491-515,2014]to the nonlinear Poisson-Boltzmann equation.The ...This paper is devoted to the error estimate for the iterative discontinuous Galerkin(IDG)method introduced in[P.Yin,Y.Huang and H.Liu.Commun.Comput.Phys.16:491-515,2014]to the nonlinear Poisson-Boltzmann equation.The total error includes both the iteration error and the discretization error of the direct DG method to linear elliptic equations.For the DDG method,the energy error is obtained by a constructive approach through an explicit global projection satisfying interface conditions dictated by the choice of numerical fluxes.The L^(2) error of order O(h^(m+1))for polynomials of degree m is further recovered.The bounding constant is also shown to be independent of the iteration times.Numerical tests are given to validate the established convergence theory.展开更多
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 will investigate the error estimates and the superconvergence property of mixed finite element methods for a semilinear elliptic control problem with an integral constraint on control.The state and co...In this paper,we will investigate the error estimates and the superconvergence property of mixed finite element methods for a semilinear elliptic control problem with an integral constraint on control.The state and co-state are approximated by the lowest order Raviart-Thomas mixed finite element and the control variable is approximated by piecewise constant functions.We derive some superconvergence properties for the control variable and the state variables.Moreover,we derive L∞-and H−1-error estimates both for the control variable and the state variables.Finally,a numerical example is given to demonstrate the theoretical results.展开更多
基金Supported by the National Natural Science Foundation of China (10601022)Natural Science Foundation of Inner Mongolia Autonomous Region (200607010106)Youth Science Foundation of Inner Mongolia University(ND0702)
文摘An H^1-Galerkin mixed finite element method is discussed for a class of second order SchrSdinger equation. Optimal error estimates of semidiscrete schemes are derived for problems in one space dimension. At the same time, optimal error estimates are derived for fully discrete schemes. And it is showed that the H1-Galerkin mixed finite element approximations have the same rate of convergence as in the classical mixed finite element methods without requiring the LBB consistency condition.
文摘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.
文摘<div style="text-align:justify;"> In this paper, we study the error estimates for direct discontinuous Galerkin methods based on the upwind-biased fluxes. We use a newly global projection to obtain the optimal error estimates. The numerical experiments imply that <em>L</em><sup>2 </sup>norms error estimates can reach to order <em>k</em> + 1 by using time discretization methods. </div>
文摘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.
文摘Allen and Liu (1995) introduced a new method for a time-dependent convection dominated diffusion problem, which combines the modified method of characteristics and method of streamline diffusion. But they ignored the fact that the accuracy of time discretization decays at half an order when the characteristic line goes out of the domain. In present paper, the author shows that, as a remedy, a simple lumped scheme yields a full accuracy approximation. Forthermore, some local error bounds independent of the small viscosity axe derived for this scheme outside the boundary layers.
基金This research was supported by the National Natural Science Foundation of China
文摘In the present paper, a new numerical method for solving initial-boundary value problems of evolutionary equations is proposed and studied, combining difference method with high accuracy with boundary integral equation method. The numerical approximate schemes for both problems on a bounded or unbounded domain in R3 are proposed and their prior error estimates are obtained.
文摘In this paper,the uniform error estimates with respect to t∈[0, ∞ ) of the nonlinear Galerkin method are given for the long time integration of the Kuramoto-Sivashinsky equation. The nonlinear Galerkin method is used to study the asymptotic behaviour of Kuramoto-Sivashinsky equation and to construct the bifurcation diagrams.
基金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 purpose of this paper is to investigate the convergence of the mixed finite element method for the initial-boundary value problem for the Sobolev equation Ut-div{aut + b1 u} = f based on the Raviart-Thomas space Vh × Wh H(div; × L2(). Optimal order estimates are obtained for the approximation of u, ut, the associated velocity p and divp respectively in L(0,T;L2()), L(0,T;L2()), L(0,T;L2()2), and L2(0, T; L2()). Quasi-optimal order estimates are obtained for the approximations of u, ut in L(0, T; L()) and p in L(0,T; L()2).
基金Subsidized by NSFC(11571274 and 11171269)the Ph.D.Programs Foundation of Ministry of Education of China(20110201110027)
文摘Although the two-grid finite element decoupled scheme for mixed Navier-Stokes/ Darcy model in literatures has given the numerical results of optimal convergence order, the theoretical analysis only obtain the optimal error order for the porous media flow and a non-optimal error order for the fluid flow. In this article, we give a more rigorous of the error analysis for the fluid flow and obtain the optimal error estimates of the velocity and the pressure.
基金Project supported by the National Natural Science Foundation of China (No. 60874039)Shanghai Leading Academic Discipline Project (No. J50101)
文摘A Fourier spectral method for the generalized Korteweg-de Vries equation with periodic boundary conditions is analyzed, and a corresponding optimal error estimate in L^2-norm is obtained. It improves the result presented by Maday and Quarteroni. A modified Fourier pseudospectral method is also presented, with the same convergence properties as the Fourier spectral method.
文摘In this paper, a spectral method to analyze the generalized Benjamin Bona Mahony equations is used. The existence and uniqueness of global smooth solution of these equations are proved. The large time error estimation between the spectral approximate solution and the exact solution is obtained.
基金Project supported by the National Natural Science Foundation of China (Nos. 10871022, 11061009, and 40821092)the National Basic Research Program of China (973 Program) (Nos. 2010CB428403, 2009CB421407, and 2010CB951001)the Natural Science Foundation of Hebei Province of China (No. A2010001663)
文摘The proper orthogonal decomposition (POD) is a model reduction technique for the simulation Of physical processes governed by partial differential equations (e.g., fluid flows). It has been successfully used in the reduced-order modeling of complex systems. In this paper, the applications of the POD method are extended, i.e., the POD method is applied to a classical finite difference (FD) scheme for the non-stationary Stokes equation with a real practical applied background. A reduced FD scheme is established with lower dimensions and sufficiently high accuracy, and the error estimates are provided between the reduced and the classical FD solutions. Some numerical examples illustrate that the numerical results are consistent with theoretical conclusions. Moreover, it is shown that the reduced FD scheme based on the POD method is feasible and efficient in solving the FD scheme for the non-stationary Stokes equation.
基金This work is supported by NSFC grants No.11601490.
文摘In this paper,we prove the optimal error estimates in L2 norm of the semidiscrete local discontinuous Galerkin(LDG)method for the thin film epitaxy problem without slope selection.To relax the severe time step restriction of explicit time marching methods,we employ a class of exponential time differencing(ETD)schemes for time integration,which is based on a linear convex splitting principle.Numerical experiments of the accuracy and long time simulations are given to show the efficiency and capability of the proposed numerical schemes.
文摘In this paper we deal with the convergence analysis of the finite element method for an elliptic penalized unilateral obstacle optimal control problem where the control and the obstacle coincide.Error estimates are established for both state and control variables.We apply a fixed point type iteration method to solve the discretized problem.To corroborate our error estimations and the eficiency of our algorithms,the convergence results and numerical experiments are illustrated by concrete examples.
基金This work was partly supported by National Natural Science Foundation of China(Grant Nos.:12101283,12271233 and 12171287)Natural Science Foundation of Shandong Province(Grant Nos.:ZR2019YQ05,2019KJI003,and ZR2016JL004).
文摘The fractional optimal control problem leads to significantly increased computational complexity compared to the corresponding classical integer-order optimal control problem,due to the global properties of fractional differential operators.In this paper,we focus on an optimal control problem governed by fractional differential equations with an integral constraint on the state variable.By the proposed first-order optimality condition consisting of a Lagrange multiplier,we design a spectral Galerkin discrete scheme with weighted orthogonal Jacobi polynomials to approximate the resulting state and adjoint state equations.Furthermore,a priori error estimates for state,adjoint state and control variables are discussed in details.Illustrative numerical tests are given to demonstrate the validity and applicability of our proposed approximations and theoretical results.
文摘In this paper, we propose the nonconforming virtual element method (NCVEM) discretization for the pointwise control constraint optimal control problem governed by elliptic equations. Based on the NCVEM approximation of state equation and the variational discretization of control variables, we construct a virtual element discrete scheme. For the state, adjoint state and control variable, we obtain the corresponding prior estimate in H<sup>1</sup> and L<sup>2</sup> norms. Finally, some numerical experiments are carried out to support the theoretical results.
基金The authors thank the referees for valuable suggestionswhich led to significant improvements in this revised version.This work was supported by the National Science Foundation of USA under Grant DMS1312636by NSF Grant RNMS(Ki-Net)1107291.Huang’s work was supported by National Science Foundation of China under Grant 91430213.
文摘This paper is devoted to the error estimate for the iterative discontinuous Galerkin(IDG)method introduced in[P.Yin,Y.Huang and H.Liu.Commun.Comput.Phys.16:491-515,2014]to the nonlinear Poisson-Boltzmann equation.The total error includes both the iteration error and the discretization error of the direct DG method to linear elliptic equations.For the DDG method,the energy error is obtained by a constructive approach through an explicit global projection satisfying interface conditions dictated by the choice of numerical fluxes.The L^(2) error of order O(h^(m+1))for polynomials of degree m is further recovered.The bounding constant is also shown to be independent of the iteration times.Numerical tests are given to validate the established convergence theory.
基金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.
基金supported by National Science Foundation of China(11271145)Foundation for Talent Introduction of Guangdong Provincial University,Specialized Research Fund for the Doctoral Program of Higher Education(20114407110009)+1 种基金the Project of Department of Education of Guangdong Province(2012KJCX0036)supported by China Postdoctoral Science Foundation funded project(2013M542188).
文摘In this paper,we will investigate the error estimates and the superconvergence property of mixed finite element methods for a semilinear elliptic control problem with an integral constraint on control.The state and co-state are approximated by the lowest order Raviart-Thomas mixed finite element and the control variable is approximated by piecewise constant functions.We derive some superconvergence properties for the control variable and the state variables.Moreover,we derive L∞-and H−1-error estimates both for the control variable and the state variables.Finally,a numerical example is given to demonstrate the theoretical results.