This paper aims to numerically study the generalized time-fractional Burgers equation in two spatial dimensions using the L1/LDG method. Here the L1 scheme is used to approximate the time-fractional derivative, i.e., ...This paper aims to numerically study the generalized time-fractional Burgers equation in two spatial dimensions using the L1/LDG method. Here the L1 scheme is used to approximate the time-fractional derivative, i.e., Caputo derivative, while the local discontinuous Galerkin (LDG) method is used to discretize the spatial derivative. If the solution has strong temporal regularity, i.e., its second derivative with respect to time being right continuous, then the L1 scheme on uniform meshes (uniform L1 scheme) is utilized. If the solution has weak temporal regularity, i.e., its first and/or second derivatives with respect to time blowing up at the starting time albeit the function itself being right continuous at the beginning time, then the L1 scheme on non-uniform meshes (non-uniform L1 scheme) is applied. Then both uniform L1/LDG and non-uniform L1/LDG schemes are constructed. They are both numerically stable and the \(L^2\) optimal error estimate for the velocity is obtained. Numerical examples support the theoretical analysis.展开更多
Fractional initial-value problems(IVPs) and time-fractional initial-boundary value problems(IBVPs), each with a Caputo temporal derivative of order α ∈(0, 1), are considered. An averaged variant of the well-known L1...Fractional initial-value problems(IVPs) and time-fractional initial-boundary value problems(IBVPs), each with a Caputo temporal derivative of order α ∈(0, 1), are considered. An averaged variant of the well-known L1 scheme is proved to be O(N^(-2)) convergent for IVPs on suitably graded meshes with N points, thereby improving the O(N^(-(2-α))) convergence rate of the standard L1 scheme. The analysis relies on a delicate decomposition of the temporal truncation error that yields a sharp dependence of the order of convergence on the degree of mesh grading used. This averaged L1 scheme can be combined with a finite difference or piecewise linear finite element discretization in space for IBVPs, and under a restriction on the temporal mesh width, one gets again O(N^(-2)) convergence in time, together with O(h^(2)) convergence in space,where h is the spatial mesh width. Numerical experiments support our results.展开更多
We study the l^1-stability of a Haxniltonian-preserving scheme, developed in [Jin and Wen, Comm. Math. Sci., 3 （2005）, 285-315], for the Liouville equation with a discontinuous potential in one space dimension. We p...We study the l^1-stability of a Haxniltonian-preserving scheme, developed in [Jin and Wen, Comm. Math. Sci., 3 （2005）, 285-315], for the Liouville equation with a discontinuous potential in one space dimension. We prove that, for suitable initial data, the scheme is stable in the l^1-norm under a hyperbolic CFL condition which is in consistent with the l^1-convergence results established in [Wen and Jin, SIAM J. Numer. Anal., 46 （2008）, 2688-2714] for the same scheme. The stability constant is shown to be independent of the computational time. We also provide a counter example to show that for other initial data, in particular, the measure-valued initial data, the numerical solution may become l^1-unstable.展开更多
Numerical approximation for a linearized time fractional KdV equation with initial singularity using L1 scheme on graded mesh is considered.It is proved that the L1 scheme can attain order 2−αconvergence rate with ap...Numerical approximation for a linearized time fractional KdV equation with initial singularity using L1 scheme on graded mesh is considered.It is proved that the L1 scheme can attain order 2−αconvergence rate with appropriate choice of the grading parameter,whereα(0<α<1)is the order of temporal Caputo fractional derivative.A fully discrete spectral scheme is constructed combing a Petrov-Galerkin spectral method for the spatial discretization,and its stability and convergence are theoretically proved.Some numerical results are provided to verify the theoretical analysis and demonstrated the sharpness of the error analysis.展开更多
This paper is concerned with numerical solutions of time-fractional nonlinear parabolic problems by a class of L1-Galerkin finite element methods.The analysis of L1 methods for time-fractional nonlinear problems is li...This paper is concerned with numerical solutions of time-fractional nonlinear parabolic problems by a class of L1-Galerkin finite element methods.The analysis of L1 methods for time-fractional nonlinear problems is limited mainly due to the lack of a fundamental Gronwall type inequality.In this paper,we establish such a fundamental inequality for the L1 approximation to the Caputo fractional derivative.In terms of the Gronwall type inequality,we provide optimal error estimates of several fully discrete linearized Galerkin finite element methods for nonlinear problems.The theoretical results are illustrated by applying our proposed methods to the time fractional nonlinear Huxley equation and time fractional Fisher equation.展开更多
By employing EQ^(ROT)_(1) nonconforming finite element,the numerical approximation is presented for multi-term time-fractional mixed sub-diffusion and diffusion-wave equation on anisotropic meshes.Comparing with the m...By employing EQ^(ROT)_(1) nonconforming finite element,the numerical approximation is presented for multi-term time-fractional mixed sub-diffusion and diffusion-wave equation on anisotropic meshes.Comparing with the multi-term time-fractional sub-diffusion equation or diffusion-wave equation,the mixed case contains a special time-space coupled derivative,which leads to many difficulties in numerical analysis.Firstly,a fully discrete scheme is established by using nonconforming finite element method(FEM)in spatial direction and L1 approximation coupled with Crank-Nicolson(L1-CN)scheme in temporal direction.Furthermore,the fully discrete scheme is proved to be unconditional stable.Besides,convergence and superclose results are derived by using the properties of EQ^(ROT)_(1) nonconforming finite element.What's more,the global superconvergence is obtained via the interpolation postprocessing technique.Finally,several numerical results are provided to demonstrate the theoretical analysis on anisotropic meshes.展开更多
In this paper,a nonlinear time-fractional coupled diffusion system is solved by using a mixed finite element(MFE)method in space combined with L1-approximation and implicit second-order backward difference scheme in t...In this paper,a nonlinear time-fractional coupled diffusion system is solved by using a mixed finite element(MFE)method in space combined with L1-approximation and implicit second-order backward difference scheme in time.The stability for nonlinear fully discrete finite element scheme is analyzed and a priori error estimates are derived.Finally,some numerical tests are shown to verify our theoretical analysis.展开更多
A two-grid finite element method with L1 scheme is presented for solving two-dimen-sional time-fractional nonlinear Schrodinger equation.The finite element solution in the L-norm are proved bounded without any time-st...A two-grid finite element method with L1 scheme is presented for solving two-dimen-sional time-fractional nonlinear Schrodinger equation.The finite element solution in the L-norm are proved bounded without any time-step size conditions(dependent on spatial-step size).The classical L1 scheme is considered in the time direction,and the two-grid finite element method is applied in spatial direction.The optimal order error estimations of the two-grid solution in the LP-norm is proved without any time-step size conditions.It is shown,both theoretically and numerically,that the coarse space can be extremely coarse,with no loss in the order of accuracy.展开更多
In this paper,we develop a two-grid method(TGM)based on the FEM for 2D nonlinear time fractional two-term mixed sub-diffusion and diffusion wave equations.A two-grid algorithm is proposed for solving the nonlinear sys...In this paper,we develop a two-grid method(TGM)based on the FEM for 2D nonlinear time fractional two-term mixed sub-diffusion and diffusion wave equations.A two-grid algorithm is proposed for solving the nonlinear system,which consists of two steps:a nonlinear FE system is solved on a coarse grid,then the linearized FE system is solved on the fine grid by Newton iteration based on the coarse solution.The fully discrete numerical approximation is analyzed,where the Galerkin finite element method for the space derivatives and the finite difference scheme for the time Caputo derivative with orderα∈(1,2)andα1∈(0,1).Numerical stability and optimal error estimate O(h^(r+1)+H^(2r+2)+τ^(min{3−α,2−α1}))in L^(2)-norm are presented for two-grid scheme,where t,H and h are the time step size,coarse grid mesh size and fine grid mesh size,respectively.Finally,numerical experiments are provided to confirm our theoretical results and effectiveness of the proposed algorithm.展开更多
The main contents of this paper are to establish a finite element fully-discrete approximate scheme for multi-term time-fractional mixed sub-diffusion and diffusionwave equation with spatial variable coefficient,which...The main contents of this paper are to establish a finite element fully-discrete approximate scheme for multi-term time-fractional mixed sub-diffusion and diffusionwave equation with spatial variable coefficient,which contains a time-space coupled derivative.The nonconforming EQ^(rot)_(1)element and Raviart-Thomas element are employed for spatial discretization,and L1 time-stepping method combined with the Crank-Nicolson scheme are applied for temporal discretization.Firstly,based on some significant lemmas,the unconditional stability analysis of the fully-discrete scheme is acquired.With the assistance of the interpolation operator I_(h)and projection operator Rh,superclose and convergence results of the variable u in H^(1)-norm and the flux~p=k_(5)(x)ru(x,t)in L^(2)-norm are obtained,respectively.Furthermore,the global superconvergence results are derived by applying the interpolation postprocessing technique.Finally,the availability and accuracy of the theoretical analysis are corroborated by experimental results of numerical examples on anisotropic meshes.展开更多
基金the National Natural Science Foundation of China(Nos.11671251 and 12101266).
文摘This paper aims to numerically study the generalized time-fractional Burgers equation in two spatial dimensions using the L1/LDG method. Here the L1 scheme is used to approximate the time-fractional derivative, i.e., Caputo derivative, while the local discontinuous Galerkin (LDG) method is used to discretize the spatial derivative. If the solution has strong temporal regularity, i.e., its second derivative with respect to time being right continuous, then the L1 scheme on uniform meshes (uniform L1 scheme) is utilized. If the solution has weak temporal regularity, i.e., its first and/or second derivatives with respect to time blowing up at the starting time albeit the function itself being right continuous at the beginning time, then the L1 scheme on non-uniform meshes (non-uniform L1 scheme) is applied. Then both uniform L1/LDG and non-uniform L1/LDG schemes are constructed. They are both numerically stable and the \(L^2\) optimal error estimate for the velocity is obtained. Numerical examples support the theoretical analysis.
基金supported by National Natural Science Foundation of China (Grant Nos. 12101509, 12171283, 12171025 and NSAF-U1930402)the Science Foundation Program for Distinguished Young Scholars of Shandong (Overseas) (Grant No. 2022HWYQ-045)。
文摘Fractional initial-value problems(IVPs) and time-fractional initial-boundary value problems(IBVPs), each with a Caputo temporal derivative of order α ∈(0, 1), are considered. An averaged variant of the well-known L1 scheme is proved to be O(N^(-2)) convergent for IVPs on suitably graded meshes with N points, thereby improving the O(N^(-(2-α))) convergence rate of the standard L1 scheme. The analysis relies on a delicate decomposition of the temporal truncation error that yields a sharp dependence of the order of convergence on the degree of mesh grading used. This averaged L1 scheme can be combined with a finite difference or piecewise linear finite element discretization in space for IBVPs, and under a restriction on the temporal mesh width, one gets again O(N^(-2)) convergence in time, together with O(h^(2)) convergence in space,where h is the spatial mesh width. Numerical experiments support our results.
基金Innovation Project of the Chinese Academy of Sciences grants K5501312S1,K5502212F1,K7290312G7 and K7502712F7NSFC grant 10601062+1 种基金NSF grant DMS-0608720NSAF grant 10676017
文摘We study the l^1-stability of a Haxniltonian-preserving scheme, developed in [Jin and Wen, Comm. Math. Sci., 3 （2005）, 285-315], for the Liouville equation with a discontinuous potential in one space dimension. We prove that, for suitable initial data, the scheme is stable in the l^1-norm under a hyperbolic CFL condition which is in consistent with the l^1-convergence results established in [Wen and Jin, SIAM J. Numer. Anal., 46 （2008）, 2688-2714] for the same scheme. The stability constant is shown to be independent of the computational time. We also provide a counter example to show that for other initial data, in particular, the measure-valued initial data, the numerical solution may become l^1-unstable.
基金The work of Hu Chen is supported in part by NSF of China(No.11801026)and China Postdoctoral Science Foundation Under No.2018M631316the work of Xiaohan Hu and Yifa Tang is supported in part by NSF of China(No.11771438)+3 种基金the work of Jincheng Ren is supported in part by NSF of China(No.11601119)sponsored by Program for HASTIT(No.18HASTIT027)Young talents Fund of HUELthe work of Tao Sun is supported in part by NSF of China(No.11401380).
文摘Numerical approximation for a linearized time fractional KdV equation with initial singularity using L1 scheme on graded mesh is considered.It is proved that the L1 scheme can attain order 2−αconvergence rate with appropriate choice of the grading parameter,whereα(0<α<1)is the order of temporal Caputo fractional derivative.A fully discrete spectral scheme is constructed combing a Petrov-Galerkin spectral method for the spatial discretization,and its stability and convergence are theoretically proved.Some numerical results are provided to verify the theoretical analysis and demonstrated the sharpness of the error analysis.
基金This work is supported by NSFC(Grant Nos.11771035,11771162,11571128,61473126,91430216,91530204,11372354 and U1530401),a grant from the RGC of HK 11300517,China(Project No.CityU 11302915),China Postdoctoral Science Foundation under grant No.2016M602273,a grant DRA2015518 from 333 High-level Personal Training Project of Jiangsu Province,and the USA National Science Foundation grant DMS-1315259the USA Air Force Office of Scientific Research grant FA9550-15-1-0001.Jiwei Zhang also thanks the hospitality of Hong Kong City University during the period of his visiting.
文摘This paper is concerned with numerical solutions of time-fractional nonlinear parabolic problems by a class of L1-Galerkin finite element methods.The analysis of L1 methods for time-fractional nonlinear problems is limited mainly due to the lack of a fundamental Gronwall type inequality.In this paper,we establish such a fundamental inequality for the L1 approximation to the Caputo fractional derivative.In terms of the Gronwall type inequality,we provide optimal error estimates of several fully discrete linearized Galerkin finite element methods for nonlinear problems.The theoretical results are illustrated by applying our proposed methods to the time fractional nonlinear Huxley equation and time fractional Fisher equation.
基金National Natural Science Foundation of China(No.11971416)Scientific Research Innovation Team of Xuchang University(No.2022CXTD002)+3 种基金Foundation for University Key Young Teacher of Henan Province(No.2019GGJS214)Key Scientific Research Projects in Universities of Henan Province(Nos.21B110007,22A110022)National Natural Science Foundation of China(International cooperation key project:No.12120101001)Australian Research Council via the Discovery Project(DP190101889).
文摘By employing EQ^(ROT)_(1) nonconforming finite element,the numerical approximation is presented for multi-term time-fractional mixed sub-diffusion and diffusion-wave equation on anisotropic meshes.Comparing with the multi-term time-fractional sub-diffusion equation or diffusion-wave equation,the mixed case contains a special time-space coupled derivative,which leads to many difficulties in numerical analysis.Firstly,a fully discrete scheme is established by using nonconforming finite element method(FEM)in spatial direction and L1 approximation coupled with Crank-Nicolson(L1-CN)scheme in temporal direction.Furthermore,the fully discrete scheme is proved to be unconditional stable.Besides,convergence and superclose results are derived by using the properties of EQ^(ROT)_(1) nonconforming finite element.What's more,the global superconvergence is obtained via the interpolation postprocessing technique.Finally,several numerical results are provided to demonstrate the theoretical analysis on anisotropic meshes.
基金the National Natural Science Fund(11661058,11301258,11361035)the Natural Science Fund of Inner Mongolia Autonomous Region(2016MS0102,2015MS0101)+1 种基金the Scientific Research Projection of Higher Schools of Inner Mongolia(NJZZ12011)the National Undergraduate Innovative Training Project(201510126026).
文摘In this paper,a nonlinear time-fractional coupled diffusion system is solved by using a mixed finite element(MFE)method in space combined with L1-approximation and implicit second-order backward difference scheme in time.The stability for nonlinear fully discrete finite element scheme is analyzed and a priori error estimates are derived.Finally,some numerical tests are shown to verify our theoretical analysis.
基金supported by the State Key Program of National Natural Science Foundation of China(Grant No.11931003)by the National Natural Science Foundation of China(Grant Nos.41974133,12271233)+1 种基金The work of first author was supported by the Guangdong Basic and Applied Basic Research Foundation(Grant Nos.2024A1515012430,2020A1515011032)by the Educational Commission of Guangdong Province,China(Grant No.2019KTSCX174).
文摘A two-grid finite element method with L1 scheme is presented for solving two-dimen-sional time-fractional nonlinear Schrodinger equation.The finite element solution in the L-norm are proved bounded without any time-step size conditions(dependent on spatial-step size).The classical L1 scheme is considered in the time direction,and the two-grid finite element method is applied in spatial direction.The optimal order error estimations of the two-grid solution in the LP-norm is proved without any time-step size conditions.It is shown,both theoretically and numerically,that the coarse space can be extremely coarse,with no loss in the order of accuracy.
基金This work is supported by the State Key Program of National Natural Science Foundation of China(11931003)National Natural Science Foundation of China(41974133,11971410)+2 种基金Project for Hunan National Applied Mathematics Center of Hunan Provincial Science and Technology Department(2020ZYT003)Hunan Provincial Innovation Foundation for Postgraduate,China(XDCX2020B082,XDCX2021B098)Postgraduate Scientific Research Innovation Project of Hunan Province(CX20210597).
文摘In this paper,we develop a two-grid method(TGM)based on the FEM for 2D nonlinear time fractional two-term mixed sub-diffusion and diffusion wave equations.A two-grid algorithm is proposed for solving the nonlinear system,which consists of two steps:a nonlinear FE system is solved on a coarse grid,then the linearized FE system is solved on the fine grid by Newton iteration based on the coarse solution.The fully discrete numerical approximation is analyzed,where the Galerkin finite element method for the space derivatives and the finite difference scheme for the time Caputo derivative with orderα∈(1,2)andα1∈(0,1).Numerical stability and optimal error estimate O(h^(r+1)+H^(2r+2)+τ^(min{3−α,2−α1}))in L^(2)-norm are presented for two-grid scheme,where t,H and h are the time step size,coarse grid mesh size and fine grid mesh size,respectively.Finally,numerical experiments are provided to confirm our theoretical results and effectiveness of the proposed algorithm.
基金The work is supported by the National Natural Science Foundation of China(Nos.11971416 and 11871441)the Scientific Research Innovation Team of Xuchang University(No.2022CXTD002)the Foundation for University Key Young Teacher of Henan Province(No.2019GGJS214).
文摘The main contents of this paper are to establish a finite element fully-discrete approximate scheme for multi-term time-fractional mixed sub-diffusion and diffusionwave equation with spatial variable coefficient,which contains a time-space coupled derivative.The nonconforming EQ^(rot)_(1)element and Raviart-Thomas element are employed for spatial discretization,and L1 time-stepping method combined with the Crank-Nicolson scheme are applied for temporal discretization.Firstly,based on some significant lemmas,the unconditional stability analysis of the fully-discrete scheme is acquired.With the assistance of the interpolation operator I_(h)and projection operator Rh,superclose and convergence results of the variable u in H^(1)-norm and the flux~p=k_(5)(x)ru(x,t)in L^(2)-norm are obtained,respectively.Furthermore,the global superconvergence results are derived by applying the interpolation postprocessing technique.Finally,the availability and accuracy of the theoretical analysis are corroborated by experimental results of numerical examples on anisotropic meshes.
