A significant obstacle impeding the advancement of the time fractional Schrodinger equation lies in the challenge of determining its precise mathematical formulation.In order to address this,we undertake an exploratio...A significant obstacle impeding the advancement of the time fractional Schrodinger equation lies in the challenge of determining its precise mathematical formulation.In order to address this,we undertake an exploration of the time fractional Schrodinger equation within the context of a non-Markovian environment.By leveraging a two-level atom as an illustrative case,we find that the choice to raise i to the order of the time derivative is inappropriate.In contrast to the conventional approach used to depict the dynamic evolution of quantum states in a non-Markovian environment,the time fractional Schrodinger equation,when devoid of fractional-order operations on the imaginary unit i,emerges as a more intuitively comprehensible framework in physics and offers greater simplicity in computational aspects.Meanwhile,we also prove that it is meaningless to study the memory of time fractional Schrodinger equation with time derivative 1<α≤2.It should be noted that we have not yet constructed an open system that can be fully described by the time fractional Schrodinger equation.This will be the focus of future research.Our study might provide a new perspective on the role of time fractional Schrodinger equation.展开更多
The Time Fractional Burger equation was solved in this study using the Mabel software and the Variational Iteration approach. where a number of instances of the Time Fractional Burger Equation were handled using this ...The Time Fractional Burger equation was solved in this study using the Mabel software and the Variational Iteration approach. where a number of instances of the Time Fractional Burger Equation were handled using this technique. Tables and images were used to present the collected numerical results. The difference between the exact and numerical solutions demonstrates the effectiveness of the Mabel program’s solution, as well as the accuracy and closeness of the results this method produced. It also demonstrates the Mabel program’s ability to quickly and effectively produce the numerical solution.展开更多
In this paper,a local discontinuous Galerkin(LDG)scheme for the time-fractional diffusion equation is proposed and analyzed.The Caputo time-fractional derivative(of orderα,with 0<α<1)is approximated by a finit...In this paper,a local discontinuous Galerkin(LDG)scheme for the time-fractional diffusion equation is proposed and analyzed.The Caputo time-fractional derivative(of orderα,with 0<α<1)is approximated by a finite difference method with an accuracy of order3-α,and the space discretization is based on the LDG method.For the finite difference method,we summarize and supplement some previous work by others,and apply it to the analysis of the convergence and stability of the proposed scheme.The optimal error estimate is obtained in the L2norm,indicating that the scheme has temporal(3-α)th-order accuracy and spatial(k+1)th-order accuracy,where k denotes the highest degree of a piecewise polynomial in discontinuous finite element space.The numerical results are also provided to verify the accuracy and efficiency of the considered scheme.展开更多
A finite difference/spectral scheme is proposed for the time fractional Ito equation.The mass conservation and stability of the numerical solution are deduced by the energy method in the L^(2)norm form.To reduce the c...A finite difference/spectral scheme is proposed for the time fractional Ito equation.The mass conservation and stability of the numerical solution are deduced by the energy method in the L^(2)norm form.To reduce the computation costs,the fast Fourier transform technic is applied to a pair of equivalent coupled differential equations.The effectiveness of the proposed algorithm is verified by the first numerical example.The mass conservation property and stability statement are confirmed by two other numerical examples.展开更多
In this paper, the Adomian decomposition method was used to solve the Time Fractional Burger equation using Mabel program. This method was applied to a number of examples of the Time Fractional Burger Equation. The ob...In this paper, the Adomian decomposition method was used to solve the Time Fractional Burger equation using Mabel program. This method was applied to a number of examples of the Time Fractional Burger Equation. The obtained numerical results were presented in the form of tables and graphics. The difference between the exact solutions and the numerical solutions shows us the effectiveness of the solution using the Mabel program and that this method gave accurate results and was close to the exact solution, in addition to its ability to obtain the numerical solution quickly and efficiently using the Mabel program.展开更多
It is an important issue to numerically solve the time fractional Schrödinger equation on unbounded domains, which models the dynamics of optical solitons propagating via optical fibers. The perfectly matched lay...It is an important issue to numerically solve the time fractional Schrödinger equation on unbounded domains, which models the dynamics of optical solitons propagating via optical fibers. The perfectly matched layer approach is applied to truncate the unbounded physical domain, and obtain an initial boundary value problem on a bounded computational domain, which can be efficiently solved by the finite difference method. The stability of the reduced initial boundary value problem is rigorously analyzed. Some numerical results are presented to illustrate the accuracy and feasibility of the perfectly matched layer approach. According to these examples, the absorption parameters and the width of the absorption layer will affect the absorption effect. The larger the absorption width, the better the absorption effect. There is an optimal absorption parameter, the absorption effect is the best.展开更多
In this work, we propose a new approach, namely ansatz method, for solving fractional differential equations based on a fractional complex transform and apply it to the nonlinear partial space-time fractional modified...In this work, we propose a new approach, namely ansatz method, for solving fractional differential equations based on a fractional complex transform and apply it to the nonlinear partial space-time fractional modified Benjamin-Bona- Mahoney (mBBM) equation, the time fractional mKdV equation and the nonlinear fractional Zoomeron equation which gives rise to some new exact solutions. The physical parameters in the soliton solutions: amplitude, inverse width, free parameters and velocity are obtained as functions of the dependent model coefficients. This method is suitable and more powerful for solving other kinds of nonlinear fractional PDEs arising in mathematical physics. Since the fractional deriva- tives are described in the modified Riemann-Liouville sense.展开更多
Time fractional diffusion equation is usually used to describe the problems involving non-Markovian random walks. This kind of equation is obtained from the standard diffusion equation by replacing the first-order tim...Time fractional diffusion equation is usually used to describe the problems involving non-Markovian random walks. This kind of equation is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α∈(0, 1). In this paper, an implicit finite difference scheme for solving the time fractional diffusion equation with source term is presented and analyzed, where the fractional derivative is described in the Caputo sense. Stability and convergence of this scheme are rigorously established by a Fourier analysis. And using numerical experiments illustrates the accuracy and effectiveness of the scheme mentioned in this paper.展开更多
This work is concerned with the application of a redefined set of extended uniform cubic B-spline(RECBS)functions for the numerical treatment of time-fractional Telegraph equation.The presented technique engages finit...This work is concerned with the application of a redefined set of extended uniform cubic B-spline(RECBS)functions for the numerical treatment of time-fractional Telegraph equation.The presented technique engages finite difference formulation for discretizing the Caputo time-fractional derivatives and RECBS functions to interpolate the solution curve along the spatial grid.Stability analysis of the scheme is provided to ensure that the errors do not amplify during the execution of the numerical procedure.The derivation of uniform convergence has also been presented.Some computational experiments are executed to verify the theoretical considerations.Numerical results are compared with the existing schemes and it is concluded that the present scheme returns superior outcomes on the topic.展开更多
In this article,time fractional Fornberg-Whitham equation of He’s fractional derivative is studied.To transform the fractional model into its equivalent differential equation,the fractional complex transform is used ...In this article,time fractional Fornberg-Whitham equation of He’s fractional derivative is studied.To transform the fractional model into its equivalent differential equation,the fractional complex transform is used and He’s homotopy perturbation method is implemented to get the approximate analytical solutions of the fractional-order problems.The graphs are plotted to analysis the fractional-order mathematical modeling.展开更多
The present paper deals with the numerical solution of time-fractional partial differential equations using the element-free Galerkin (EFG) method, which is based on the moving least-square approximation. Compared w...The present paper deals with the numerical solution of time-fractional partial differential equations using the element-free Galerkin (EFG) method, which is based on the moving least-square approximation. Compared with numerical methods based on meshes, the EFG method for time-fractional partial differential equations needs only scattered nodes instead of meshing the domain of the problem. It neither requires element connectivity nor suffers much degradation in accuracy when nodal arrangements are very irregular. In this method, the first-order time derivative is replaced by the Caputo fractional derivative of order α(0 〈 α≤ 1). The Galerkin weak form is used to obtain the discrete equations, and the essential boundary conditions are enforced by the penalty method. Several numerical examples are presented and the results we obtained are in good agreement with the exact solutions.展开更多
A fast Time Domain Integral Equation(TDIE) solver is presented for analysis of transient scattering from electrically large conducting complex objects.The numerical process of Marching-On-in-Time(MOT) method based TDI...A fast Time Domain Integral Equation(TDIE) solver is presented for analysis of transient scattering from electrically large conducting complex objects.The numerical process of Marching-On-in-Time(MOT) method based TDIE encounters high computational cost and exorbitant memory requirements.A group-style accelerated method-Plane Wave Time Domain(PWTD) algorithm,which permits rapid evaluation of transient wave field generated by temporally bandlimited sources,is employed to reduce the computational cost of MOT-based TDIE solvers.An efficient compressed storage technique for sparse matrix is adopted to decrease the enormous memory requirements of MOT.The scheme of the Multi-Level PWTD(MLPWTD)-enhanced MOT with compressed storage for sparse matrix is presented for analysis of transient scattering from electrically large complex objects in this paper.The numerical simulation results demonstrate the validity and efficiency of the presented scheme.展开更多
For two-dimensional(2D)time fractional diffusion equations,we construct a numerical method based on a local discontinuous Galerkin(LDG)method in space and a finite differ-ence scheme in time.We investigate the numeric...For two-dimensional(2D)time fractional diffusion equations,we construct a numerical method based on a local discontinuous Galerkin(LDG)method in space and a finite differ-ence scheme in time.We investigate the numerical stability and convergence of the method for both rectangular and triangular meshes and show that the method is unconditionally stable.Numerical results indicate the effectiveness and accuracy of the method and con-firm the analysis.展开更多
In this paper, We show for isentropic equations of gas dynamics with adiabatic exponent gamma=3 that approximations of weak solutions generated by large time step Godunov's scheme or Glimm's scheme give entrop...In this paper, We show for isentropic equations of gas dynamics with adiabatic exponent gamma=3 that approximations of weak solutions generated by large time step Godunov's scheme or Glimm's scheme give entropy solution in the limit if Courant number is less than or equal to 1.展开更多
The generalized master equation for the space-time coupled continuous time random walk is derived analytically, in which the space-time coupling is considered through the correlated function 9(t) ~ t^γ, 0 ≤ γ 〈...The generalized master equation for the space-time coupled continuous time random walk is derived analytically, in which the space-time coupling is considered through the correlated function 9(t) ~ t^γ, 0 ≤ γ 〈 2, and the probability density function ω(t) of a particle's waiting time t follows a power law form for large t: ω(t) ~t^-(1+α), 0 〈 α 〈 1. The results indicate that the expressions of the generalized master equation are determined by the correlation exponent 7 and the long-tailed index α of the waiting time. Moreover, the diffusion results obtained from the generalized master equation are in accordance with the previous known results and the numerical simulation results.展开更多
In this paper, we use the Mittag-Leffler addition formula to solve the Green function of generalized time fractional diffusion equation in the whole plane and prove the convergence of the Green function.
In this paper, we introduce high-order finite volume methods for the multi-term time fractional sub-diffusion equation. The time fractional derivatives are described in Caputo’s sense. By using some operators, we obt...In this paper, we introduce high-order finite volume methods for the multi-term time fractional sub-diffusion equation. The time fractional derivatives are described in Caputo’s sense. By using some operators, we obtain the compact finite volume scheme have high order accuracy. We use a compact operator to deal with spatial direction;then we can get the compact finite volume scheme. It is proved that the finite volume scheme is unconditionally stable and convergent in L<sub>∞</sub>-norm. The convergence order is O(τ<sup>2-α</sup> + h<sup>4</sup>). Finally, two numerical examples are given to confirm the theoretical results. Some tables listed also can explain the stability and convergence of the scheme.展开更多
We propose a solution method of Time Dependent Schr?dinger Equation (TDSE) and the advection equation by quantum walk/quantum cellular automaton with spatially or temporally variable parameters. Using numerical method...We propose a solution method of Time Dependent Schr?dinger Equation (TDSE) and the advection equation by quantum walk/quantum cellular automaton with spatially or temporally variable parameters. Using numerical method, we establish the quantitative relation between the quantum walk with the space dependent parameters and the “Time Dependent Schr?dinger Equation with a space dependent imaginary diffusion coefficient” or “the advection equation with space dependent velocity fields”. Using the 4-point-averaging manipulation in the solution of advection equation by quantum walk, we find that only one component can be extracted out of two components of left-moving and right-moving solutions. In general it is not so easy to solve an advection equation without numerical diffusion, but this method provides perfectly diffusion free solution by virtue of its unitarity. Moreover our findings provide a clue to find more general space dependent formalisms such as solution method of TDSE with space dependent resolution by quantum walk.展开更多
The Cauchy problem for some parabolic fractional partial differential equation of higher orders and with time delays is considered. The existence and unique solution of this problem is studied. Some smoothness propert...The Cauchy problem for some parabolic fractional partial differential equation of higher orders and with time delays is considered. The existence and unique solution of this problem is studied. Some smoothness properties with respect to the parameters of these delay fractional differential equations are considered.展开更多
In this work,a novel time-stepping L1 formula is developed for a hidden-memory variable-order Caputo’s fractional derivative with an initial singularity.This formula can obtain second-order accuracy and an error esti...In this work,a novel time-stepping L1 formula is developed for a hidden-memory variable-order Caputo’s fractional derivative with an initial singularity.This formula can obtain second-order accuracy and an error estimate is analyzed strictly.As an application,a fully discrete difference scheme is established for the initial-boundary value problem of a hidden-memory variable-order time fractional diffusion model.Numerical experiments are provided to support our theoretical results.展开更多
基金Project supported by the National Natural Science Foun dation of China(Grant No.11274398).
文摘A significant obstacle impeding the advancement of the time fractional Schrodinger equation lies in the challenge of determining its precise mathematical formulation.In order to address this,we undertake an exploration of the time fractional Schrodinger equation within the context of a non-Markovian environment.By leveraging a two-level atom as an illustrative case,we find that the choice to raise i to the order of the time derivative is inappropriate.In contrast to the conventional approach used to depict the dynamic evolution of quantum states in a non-Markovian environment,the time fractional Schrodinger equation,when devoid of fractional-order operations on the imaginary unit i,emerges as a more intuitively comprehensible framework in physics and offers greater simplicity in computational aspects.Meanwhile,we also prove that it is meaningless to study the memory of time fractional Schrodinger equation with time derivative 1<α≤2.It should be noted that we have not yet constructed an open system that can be fully described by the time fractional Schrodinger equation.This will be the focus of future research.Our study might provide a new perspective on the role of time fractional Schrodinger equation.
文摘The Time Fractional Burger equation was solved in this study using the Mabel software and the Variational Iteration approach. where a number of instances of the Time Fractional Burger Equation were handled using this technique. Tables and images were used to present the collected numerical results. The difference between the exact and numerical solutions demonstrates the effectiveness of the Mabel program’s solution, as well as the accuracy and closeness of the results this method produced. It also demonstrates the Mabel program’s ability to quickly and effectively produce the numerical solution.
基金supported by the State Key Program of National Natural Science Foundation of China(11931003)the National Natural Science Foundation of China(41974133)。
文摘In this paper,a local discontinuous Galerkin(LDG)scheme for the time-fractional diffusion equation is proposed and analyzed.The Caputo time-fractional derivative(of orderα,with 0<α<1)is approximated by a finite difference method with an accuracy of order3-α,and the space discretization is based on the LDG method.For the finite difference method,we summarize and supplement some previous work by others,and apply it to the analysis of the convergence and stability of the proposed scheme.The optimal error estimate is obtained in the L2norm,indicating that the scheme has temporal(3-α)th-order accuracy and spatial(k+1)th-order accuracy,where k denotes the highest degree of a piecewise polynomial in discontinuous finite element space.The numerical results are also provided to verify the accuracy and efficiency of the considered scheme.
基金the National Natural Science Foundation of China(No.11701103)the Young Top-notch Talent Program of Guangdong Province of China(No.2017GC010379)+4 种基金the Natural Science Foundation of Guangdong Province of China(No.2022A1515012147)the Project of Science and Technology of Guangzhou of China(No.202102020704)the Opening Project of Guangdong Province Key Laboratory of Computational Science at the Sun Yat-sen University of China(2021023)the Science and Technology Development Fund,Macao SAR(File No.0005/2019/A)the University of Macao of China(File Nos.MYRG2020-00035-FST,MYRG2018-00047-FST).
文摘A finite difference/spectral scheme is proposed for the time fractional Ito equation.The mass conservation and stability of the numerical solution are deduced by the energy method in the L^(2)norm form.To reduce the computation costs,the fast Fourier transform technic is applied to a pair of equivalent coupled differential equations.The effectiveness of the proposed algorithm is verified by the first numerical example.The mass conservation property and stability statement are confirmed by two other numerical examples.
文摘In this paper, the Adomian decomposition method was used to solve the Time Fractional Burger equation using Mabel program. This method was applied to a number of examples of the Time Fractional Burger Equation. The obtained numerical results were presented in the form of tables and graphics. The difference between the exact solutions and the numerical solutions shows us the effectiveness of the solution using the Mabel program and that this method gave accurate results and was close to the exact solution, in addition to its ability to obtain the numerical solution quickly and efficiently using the Mabel program.
文摘It is an important issue to numerically solve the time fractional Schrödinger equation on unbounded domains, which models the dynamics of optical solitons propagating via optical fibers. The perfectly matched layer approach is applied to truncate the unbounded physical domain, and obtain an initial boundary value problem on a bounded computational domain, which can be efficiently solved by the finite difference method. The stability of the reduced initial boundary value problem is rigorously analyzed. Some numerical results are presented to illustrate the accuracy and feasibility of the perfectly matched layer approach. According to these examples, the absorption parameters and the width of the absorption layer will affect the absorption effect. The larger the absorption width, the better the absorption effect. There is an optimal absorption parameter, the absorption effect is the best.
文摘In this work, we propose a new approach, namely ansatz method, for solving fractional differential equations based on a fractional complex transform and apply it to the nonlinear partial space-time fractional modified Benjamin-Bona- Mahoney (mBBM) equation, the time fractional mKdV equation and the nonlinear fractional Zoomeron equation which gives rise to some new exact solutions. The physical parameters in the soliton solutions: amplitude, inverse width, free parameters and velocity are obtained as functions of the dependent model coefficients. This method is suitable and more powerful for solving other kinds of nonlinear fractional PDEs arising in mathematical physics. Since the fractional deriva- tives are described in the modified Riemann-Liouville sense.
基金Supported by the Discipline Construction and Teaching Research Fund of LUTcte(20140089)
文摘Time fractional diffusion equation is usually used to describe the problems involving non-Markovian random walks. This kind of equation is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α∈(0, 1). In this paper, an implicit finite difference scheme for solving the time fractional diffusion equation with source term is presented and analyzed, where the fractional derivative is described in the Caputo sense. Stability and convergence of this scheme are rigorously established by a Fourier analysis. And using numerical experiments illustrates the accuracy and effectiveness of the scheme mentioned in this paper.
文摘This work is concerned with the application of a redefined set of extended uniform cubic B-spline(RECBS)functions for the numerical treatment of time-fractional Telegraph equation.The presented technique engages finite difference formulation for discretizing the Caputo time-fractional derivatives and RECBS functions to interpolate the solution curve along the spatial grid.Stability analysis of the scheme is provided to ensure that the errors do not amplify during the execution of the numerical procedure.The derivation of uniform convergence has also been presented.Some computational experiments are executed to verify the theoretical considerations.Numerical results are compared with the existing schemes and it is concluded that the present scheme returns superior outcomes on the topic.
基金supported by the National Natural Science Foundation of China under Grant No.11561051。
文摘In this article,time fractional Fornberg-Whitham equation of He’s fractional derivative is studied.To transform the fractional model into its equivalent differential equation,the fractional complex transform is used and He’s homotopy perturbation method is implemented to get the approximate analytical solutions of the fractional-order problems.The graphs are plotted to analysis the fractional-order mathematical modeling.
基金Project supported by the National Natural Science Foundation of China(Grant No.11072117)the Natural Science Foundationof Zhejiang Province,China(Grant Nos.Y6110007and Y6110502)the K.C.Wong Magna Fund in Ningbo University,China
文摘The present paper deals with the numerical solution of time-fractional partial differential equations using the element-free Galerkin (EFG) method, which is based on the moving least-square approximation. Compared with numerical methods based on meshes, the EFG method for time-fractional partial differential equations needs only scattered nodes instead of meshing the domain of the problem. It neither requires element connectivity nor suffers much degradation in accuracy when nodal arrangements are very irregular. In this method, the first-order time derivative is replaced by the Caputo fractional derivative of order α(0 〈 α≤ 1). The Galerkin weak form is used to obtain the discrete equations, and the essential boundary conditions are enforced by the penalty method. Several numerical examples are presented and the results we obtained are in good agreement with the exact solutions.
文摘A fast Time Domain Integral Equation(TDIE) solver is presented for analysis of transient scattering from electrically large conducting complex objects.The numerical process of Marching-On-in-Time(MOT) method based TDIE encounters high computational cost and exorbitant memory requirements.A group-style accelerated method-Plane Wave Time Domain(PWTD) algorithm,which permits rapid evaluation of transient wave field generated by temporally bandlimited sources,is employed to reduce the computational cost of MOT-based TDIE solvers.An efficient compressed storage technique for sparse matrix is adopted to decrease the enormous memory requirements of MOT.The scheme of the Multi-Level PWTD(MLPWTD)-enhanced MOT with compressed storage for sparse matrix is presented for analysis of transient scattering from electrically large complex objects in this paper.The numerical simulation results demonstrate the validity and efficiency of the presented scheme.
文摘For two-dimensional(2D)time fractional diffusion equations,we construct a numerical method based on a local discontinuous Galerkin(LDG)method in space and a finite differ-ence scheme in time.We investigate the numerical stability and convergence of the method for both rectangular and triangular meshes and show that the method is unconditionally stable.Numerical results indicate the effectiveness and accuracy of the method and con-firm the analysis.
基金Supported in part by the National Natural Science of China, NSF Grant No. DMS-8657319.
文摘In this paper, We show for isentropic equations of gas dynamics with adiabatic exponent gamma=3 that approximations of weak solutions generated by large time step Godunov's scheme or Glimm's scheme give entropy solution in the limit if Courant number is less than or equal to 1.
基金Supported by the National Natural Science Foundation of China under Grant Nos 11605003 and 11547231
文摘The generalized master equation for the space-time coupled continuous time random walk is derived analytically, in which the space-time coupling is considered through the correlated function 9(t) ~ t^γ, 0 ≤ γ 〈 2, and the probability density function ω(t) of a particle's waiting time t follows a power law form for large t: ω(t) ~t^-(1+α), 0 〈 α 〈 1. The results indicate that the expressions of the generalized master equation are determined by the correlation exponent 7 and the long-tailed index α of the waiting time. Moreover, the diffusion results obtained from the generalized master equation are in accordance with the previous known results and the numerical simulation results.
文摘In this paper, we use the Mittag-Leffler addition formula to solve the Green function of generalized time fractional diffusion equation in the whole plane and prove the convergence of the Green function.
文摘In this paper, we introduce high-order finite volume methods for the multi-term time fractional sub-diffusion equation. The time fractional derivatives are described in Caputo’s sense. By using some operators, we obtain the compact finite volume scheme have high order accuracy. We use a compact operator to deal with spatial direction;then we can get the compact finite volume scheme. It is proved that the finite volume scheme is unconditionally stable and convergent in L<sub>∞</sub>-norm. The convergence order is O(τ<sup>2-α</sup> + h<sup>4</sup>). Finally, two numerical examples are given to confirm the theoretical results. Some tables listed also can explain the stability and convergence of the scheme.
基金supported in part by TUT Programs on Advanced Simulation Engineering,Toyohashi University of Technology.
文摘We propose a solution method of Time Dependent Schr?dinger Equation (TDSE) and the advection equation by quantum walk/quantum cellular automaton with spatially or temporally variable parameters. Using numerical method, we establish the quantitative relation between the quantum walk with the space dependent parameters and the “Time Dependent Schr?dinger Equation with a space dependent imaginary diffusion coefficient” or “the advection equation with space dependent velocity fields”. Using the 4-point-averaging manipulation in the solution of advection equation by quantum walk, we find that only one component can be extracted out of two components of left-moving and right-moving solutions. In general it is not so easy to solve an advection equation without numerical diffusion, but this method provides perfectly diffusion free solution by virtue of its unitarity. Moreover our findings provide a clue to find more general space dependent formalisms such as solution method of TDSE with space dependent resolution by quantum walk.
文摘The Cauchy problem for some parabolic fractional partial differential equation of higher orders and with time delays is considered. The existence and unique solution of this problem is studied. Some smoothness properties with respect to the parameters of these delay fractional differential equations are considered.
基金supported by the National Natural Science Foundation of China(No.12201076)the China Postdoctoral Science Foundation(No.2023M732180)。
文摘In this work,a novel time-stepping L1 formula is developed for a hidden-memory variable-order Caputo’s fractional derivative with an initial singularity.This formula can obtain second-order accuracy and an error estimate is analyzed strictly.As an application,a fully discrete difference scheme is established for the initial-boundary value problem of a hidden-memory variable-order time fractional diffusion model.Numerical experiments are provided to support our theoretical results.