High-Order Accurate Runge-Kutta (Local) Discontinuous Galerkin Methods for One- and Two-Dimensional Fractional Diffusion Equations

As the generalization of the integer order partial differential equations (PDE), the fractional order PDEs are drawing more and more attention for their applications in fluid flow, finance and other areas. This paper presents high-order accurate Runge-Kutta local discontinuous Galerkin (DG) methods for one- and two-dimensional fractional diffusion equations containing derivatives of fractional order in space. The Caputo derivative is chosen as the representation of spatial derivative, because it may represent the fractional derivative by an integral operator. Some numerical examples show that the convergence orders of the proposed local $P^k$-DG methods are $O(h^{k+1})$ both in one and two dimensions, where $P^k$ denotes the space of the real-valued polynomials with degree at most $k$.     

Fast second-order accurate difference schemes for time distributed-order and Riesz space fractional diffusion equations

The aim of this paper is to develop fast second-order accurate difference schemes for solving one- and two-dimensional time distributed-order and Riesz space fractional diffusion equations. We adopt the same measures for one- and two-dimensional problems as follows: we first transform the time distributed-order fractional diffusion problem into the multi-term time-space fractional diffusion problem with the composite trapezoid formula. Then, we propose a second-order accurate difference scheme based on the interpolation approximation on a special point to solve the resultant problem. Meanwhile, the unconditional stability and convergence of the new difference scheme in $L_2$-norm are proved. Furthermore, we find that the discretizations lead to a series of Toeplitz systems which can be efficiently solved by Krylov subspace methods with suitable circulant preconditioners. Finally, numerical results are presented to show the effectiveness of the proposed difference methods and demonstrate the fast convergence of our preconditioned Krylov subspace methods.     

A Short Note on the Quasilinearization Method for Fractional Differential Equations

Recent literature shows that for certain classes of fractional differential equations the monotone iterative technique fails to guarantee the quadratic convergence of the quasilinearization method. The present work proves the quadratic convergence of the quasilinearization method and the existence and uniqueness of the solution of such a class of fractional differential equations. Our analysis depends upon the classical Kantorovich theorem on Newton's method. Various examples are discussed in order to illustrate our approach.     

分数(k,q)阶差分方程的解

本文首次提出了一种分数阶差分,分数阶和分以及分数阶差分方程的定义,并利用Z变换理论,给出(k,q)阶常系数分数阶差分方程的具体解法.     

(2,q)阶分数差分方程的解

首次提出了一种分数阶差分,分数阶和分以及分数阶差分方程的定义,并给出(2,q)阶常系数分数阶差分方程的具体解法.     

时间分数阶扩散方程的一个新的高阶有限差分/谱逼近

研究时间分数阶扩散方程,结合时间方向的有限差分格式和空间方向的Legendre Collocation谱方法,构造了一个高阶稳定数值格式.数值算例表明该格式是无条件稳定和长时间稳定的,其收敛阶为O(△t3-α+N-m),其中△t,N和m分别是时间步长,空间多项式阶数以及精确解的正则度.     

Stability and Convergence of an Implicit Difference Approximation for the Space Riesz Fractional Reaction-Dispersion Equation

In this paper,we consider a Riesz space-fractional reaction-dispersion equation (RSFRDE).The RSFRDE is obtained from the classical reaction-dispersion equation by replacing the second-order space derivative with a Riesz derivative of orderβ∈(1,2]. We propose an implicit finite difference approximation for RSFRDE.The stability and convergence of the finite difference approximations are analyzed.Numerical results are found in good agreement with the theoretical analysis.     

A Note on Krylov Methods for Fractional Evolution Problems

The solution of linear fractional-order differential problems is addressed. For this purpose, rational approximations obtained by projections on resolvent Krylov subspaces are considered. Their convergence properties in Hilbert spaces are investigated.     

Analysis of an Implicit Finite Difference Scheme for Time Fractional Diffusion Equation

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.     

Diagonal and Toeplitz splitting iteration methods for diagonal‐plus‐Toeplitz linear systems from spatial fractional diffusion equations

The finite difference discretization of the spatial fractional diffusion equations gives discretized linear systems whose coefficient matrices have a diagonal‐plus‐Toeplitz structure. For solving these diagonal‐plus‐Toeplitz linear systems, we construct a class of diagonal and Toeplitz splitting iteration methods and establish its unconditional convergence theory. In particular, we derive a sharp upper bound about its asymptotic convergence rate and deduct the optimal value of its iteration parameter. The diagonal and Toeplitz splitting iteration method naturally leads to a diagonal and circulant splitting preconditioner. Analysis shows that the eigenvalues of the corresponding preconditioned matrix are clustered around 1, especially when the discretization step‐size h is small. Numerical results exhibit that the diagonal and circulant splitting preconditioner can significantly improve the convergence properties of GMRES and BiCGSTAB, and these preconditioned Krylov subspace iteration methods outperform the conjugate gradient method preconditioned by the approximate inverse circulant‐plus‐diagonal preconditioner proposed recently by Ng and Pan (M.K. Ng and J.‐Y. Pan, SIAM J. Sci. Comput. 2010;32:1442‐1464). Moreover, unlike this preconditioned conjugate gradient method, the preconditioned GMRES and BiCGSTAB methods show h‐independent convergence behavior even for the spatial fractional diffusion equations of discontinuous or big‐jump coefficients.     

Spectral Optimization Methods for the Time Fractional Diffusion Inverse Problem

An inverse problem of reconstructing the initial condition for a time fractional diffusion equation is investigated. On the basis of the optimal control framework, the uniqueness and first order necessary optimality condition of the minimizer for the objective functional are established, and a time-space spectral method is proposed to numerically solve the resulting minimization problem. The contribution of the paper is threefold: 1) a priori error estimate for the spectral approximation is derived; 2) a conjugate gradient optimization algorithm is designed to efficiently solve the inverse problem; 3) some numerical experiments are carried out to show that the proposed method is capable to find out the optimal initial condition, and that the convergence rate of the method is exponential if the optimal initial condition is smooth.     

变时间分数阶扩散方程的非一致时间网格有限差分方法

本文在非一致时间网格上,使用有限差分方法求解变时间分数阶扩散方程?α(x,t)u(x,t)/tα(x,t)-2u(x,t)/x2=f(x,t),0α(x,t)q≤1,证明了该方法在最大范数下的稳定性与收敛性,收敛阶为C(Δt2-q+h2).数值实例验证了理论分析的结果.     

Poisson和公式在分数阶热方程中的应用

本文利用Poisson和公式,证明了如下分数阶热方程(D_t~αlu=D_x~2u u(x1 0)=f(x))当f分别为周期函数和f∈S(■)时(速降函数空间),它们的热核满足关系H_t~α(x)=∑n=-∞H_t~α(x+n)进一步,我们把结论推广到更一般的分数阶微分方程和高维情形     

反应扩散方程解的渐近性态

研究了一类反应扩散方程初始边值问题. 利用微分不等式理论得到了问题解的渐近性态.     

High Accuracy Spectral Method for the Space-Fractional Diffusion Equations

In this paper, a high order accurate spectral method is presented for the space-fractional diffusion equations. Based on Fourier spectral method in space and Chebyshev collocation method in time, three high order accuracy schemes are proposed. The main advantages of this method are that it yields a fully diagonal representation of the fractional operator, with increased accuracy and efficiency compared with low-order counterparts, and a completely straightforward extension to high spatial dimensions. Some numerical examples, including Allen-Cahn equation, are conducted to verify the effectiveness of this method.     

The Numerical Simulation of Space-Time Variable Fractional Order Diffusion Equation

Many physical processes appear to exhibit fractional order behavior that may vary with time or space. The continuum of order in the fractional calculus allows the order of the fractional operator to be considered as a variable. Numerical methods and analysis of stability and convergence of numerical scheme for the variable fractional order partial differential equations are quite limited and difficult to derive. This motivates us to develop efficient numerical methods as well as stability and convergence of the implicit numerical methods for the space-time variable fractional order diffusion equation on a finite domain. It is worth mentioning that here we use the Coimbra-definition variable time fractional derivative which is more efficient from the numerical standpoint and is preferable for modeling dynamical systems. An implicit Euler approximation is proposed and then the stability and convergence of the numerical scheme are investigated. Finally, numerical examples are provided to show that the implicit Euler approximation is computationally efficient.     

On the convergence of two‐level Krylov methods for singular symmetric systems

We discuss the convergence of a two‐level version of the multilevel Krylov method for solving linear systems of equations with symmetric positive semidefinite matrix of coefficients. The analysis is based on the convergence result of Brown and Walker for the Generalized Minimal Residual method (GMRES), with the left‐ and right‐preconditioning implementation of the method. Numerical results based on diffusion problems are presented to show the convergence.     

Numerical Solution of Fractional Diffusion-Wave Equation

In this article, a novel compact finite difference scheme is mboxconstructed to solve the fractional diffusion-wave equation based on its equivalent integro-differential equation. In the temporal direction, the product trapezoidal scheme is employed to treat the fractional integral term. The convergence and stability of the scheme are proved. Numerical examples are also provided to verify the theoretical analysis.     

一类分数阶非线性波方程的精确解及应用

基于分离变量的思想构造了分数阶非线性波方程含常系数的解的形式.在用待定系数法求解时,根据原方程确定假设解中的待定参数,得到具体解的表达式.利用该方法求解了3个非线性波方程,即分数阶CH(Camassa-Holm)方程、时间分数阶空间五阶Kdv-like方程、分数阶广义Ostrovsky方程.比较简便地得到了这些方程的精确解.文献中关于整数阶非线性波方程的结果成为本文结果的特例.通过数值模拟给出了部分解的图像.对能够通过待定系数法求出精确解的分数阶微分方程所应满足的条件进行了阐述.     

Eikonal型方程粘性解的表达式

对于圆锥型和棱锥型Hamiltonian的Eikonal型方程,本文给出了一种几何方法,得出其初值问题解的表达式并且说明由此式给出的解为原初值问题的粘性解．首先用一个凸函数序列逼近Eikonal型方程中的Hamiltonian,再由Hopf-Lax公式给出方程序列的粘性解,最后证明了该粘性解序列会收敛到Eikonal方程的粘性解．