A simple criterion is studied for the first time for identifying the discrete energy dissipation of the Crank-Nicolson scheme for Maxwell’s equations in a Cole-Cole dispersive medium.Several numerical formulas that a...A simple criterion is studied for the first time for identifying the discrete energy dissipation of the Crank-Nicolson scheme for Maxwell’s equations in a Cole-Cole dispersive medium.Several numerical formulas that approximate the time fractional derivatives are investigated based on this criterion,including the L1 formula,the fractional BDF-2,and the shifted fractional trapezoidal rule(SFTR).Detailed error analysis is provided within the framework of time domain mixed finite element methods for smooth solutions.The convergence results and discrete energy dissipation law are confirmed by numerical tests.For nonsmooth solutions,the method SFTR can still maintain the optimal convergence order at a positive time on uniform meshes.Authors believe this is the first appearance that a second-order time-stepping method can restore the optimal convergence rate for Maxwell’s equations in a Cole-Cole dispersive medium regardless of the initial singularity of the solution.展开更多
In this article, we report the derivation of high accuracy finite difference method based on arithmetic average discretization for the solution of Un=F(x,u,u′)+∫K(x,s)ds , 0 x s < 1 subject to natural boundary co...In this article, we report the derivation of high accuracy finite difference method based on arithmetic average discretization for the solution of Un=F(x,u,u′)+∫K(x,s)ds , 0 x s < 1 subject to natural boundary conditions on a non-uniform mesh. The proposed variable mesh approximation is directly applicable to the integro-differential equation with singular coefficients. We need not require any special discretization to obtain the solution near the singular point. The convergence analysis of a difference scheme for the diffusion convection equation is briefly discussed. The presented variable mesh strategy is applicable when the internal grid points of the solution space are both even and odd in number as compared to the method discussed by authors in their previous work in which the internal grid points are strictly odd in number. The advantage of using this new variable mesh strategy is highlighted computationally.展开更多
基金supported in part by the Grant No.NSFC 12201322supported in part by Grant No.NSFC 12061053+1 种基金supported in part by the Grant Nos.NSFC 12161063 and the NSF of Inner Mongolia 2021MS01018supported in part by Grant Nos.NSFC 11871092 and NSAF U1930402.
文摘A simple criterion is studied for the first time for identifying the discrete energy dissipation of the Crank-Nicolson scheme for Maxwell’s equations in a Cole-Cole dispersive medium.Several numerical formulas that approximate the time fractional derivatives are investigated based on this criterion,including the L1 formula,the fractional BDF-2,and the shifted fractional trapezoidal rule(SFTR).Detailed error analysis is provided within the framework of time domain mixed finite element methods for smooth solutions.The convergence results and discrete energy dissipation law are confirmed by numerical tests.For nonsmooth solutions,the method SFTR can still maintain the optimal convergence order at a positive time on uniform meshes.Authors believe this is the first appearance that a second-order time-stepping method can restore the optimal convergence rate for Maxwell’s equations in a Cole-Cole dispersive medium regardless of the initial singularity of the solution.
文摘In this article, we report the derivation of high accuracy finite difference method based on arithmetic average discretization for the solution of Un=F(x,u,u′)+∫K(x,s)ds , 0 x s < 1 subject to natural boundary conditions on a non-uniform mesh. The proposed variable mesh approximation is directly applicable to the integro-differential equation with singular coefficients. We need not require any special discretization to obtain the solution near the singular point. The convergence analysis of a difference scheme for the diffusion convection equation is briefly discussed. The presented variable mesh strategy is applicable when the internal grid points of the solution space are both even and odd in number as compared to the method discussed by authors in their previous work in which the internal grid points are strictly odd in number. The advantage of using this new variable mesh strategy is highlighted computationally.
