A second order finite difference-spectral method for space fractional diffusion equations

A high order finite difference-spectral method is derived for solving space fractional diffusion equations,by combining the second order finite difference method in time and the spectral Galerkin method in space.The stability and error estimates of the temporal semidiscrete scheme are rigorously discussed,and the convergence order of the proposed method is proved to be O(τ2+Nα-m)in L2-norm,whereτ,N,αand m are the time step size,polynomial degree,fractional derivative index and regularity of the exact solution,respectively.Numerical experiments are carried out to demonstrate the theoretical analysis.     

Numerical simulation for the three-dimension fractional sub-diffusion equation

Fractional sub-diffusion equations have been widely used to model sub-diffusive systems. Most algorithms are designed for one-dimensional problems due to the memory effect in fractional derivative. In this paper, the numerical simulation of the 3D fractional sub-diffusion equation with a time fractional derivative of order $\alpha$ $(0<\alpha <1)$ is considered. A fractional alternating direction implicit scheme (FADIS) is proposed. We prove that FADIS is uniquely solvable, unconditionally stable and convergent in $H^1$ norm by the energy method. A numerical example is given to demonstrate the efficiency of FADIS.     

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.     

An orthogonal spline collocation alternating direction implicit method for second-order hyperbolic problems

On a rectangular region, we consider a linear second-order hyperbolicinitial-boundary value problem involving a mixed derivativeterm, continuous variable coefficients and non-homogeneous Dirichletboundary conditions. In comparison to the alternating directionimplicit Laplace-modified method of Fernandes (1997), we formulateand analyse a new parameter-free alternating direction implicitscheme in which the standard central difference formula is usedfor the time approximation and orthogonal spline collocationis used for the spatial discretization. We establish unconditionalstability of the scheme, and its optimal order in the discretemaximum norm in time and the H¹ norm in space. Numerical experimentsindicate that the new scheme, which has the same order as themethod of Fernandes (1997, Numer. Math., 77, 223–241),is more accurate. We also show that the new scheme is easilygeneralized to the second-order hyperbolic problems on rectangularpolygons. Extensions of the scheme to problems with discontinuouscoefficients, nonlinear problems, and problems with other boundaryconditions are also discussed.     

二维水动力弥散的随机数学模型

首先以二维水动力弥散实验为基础 ,根据实验所提供的信息 ,建立了二维水动力弥散实验的随机水质数学模型 .其次我们对所建模型进行了数值模拟计算 ,计算结果与实测数据比较吻合 ,从而说明所建立的随机模型是合理的 ,有实际应用价值 .     

修正交替方向法求解具有分数阶导数的二维非连续渗流问题

本文考虑在二维均匀介质中带有分数阶导数的非连续渗流问题,此模型修正了众所周知的 Daxcy原理.利用Riemann-Liouvifie和Griinwald-Letnikov分数阶导数之间的关系,提出了求解在二维均匀介质中带有分数阶导数的非连续渗流问题的两种修正的交替方向法:修正的交替方向隐式Euler方法和修正的Peaceman.-Rachford方法.我们讨论了这两种方法的稳定性,相容性和收敛性.最后给出数值例子.     

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.     

A second-order numerical method for two-dimensional two-sided space fractional convection diffusion equation

Space fractional convection diffusion equation describes physical phenomena where particles or energy (or other physical quantities) are transferred inside a physical system due to two processes: convection and superdiffusion. In this paper, we discuss the practical alternating directions implicit method to solve the two-dimensional two-sided space fractional convection diffusion equation on a finite domain. We theoretically prove and numerically verify that the presented finite difference scheme is unconditionally von Neumann stable and second order convergent in both space and time directions.     

IMPLICIT DIFFERENCE APPROXIMATION FOR A TIME FRACTIONAL ADVECTION-DISPERSION EQUATION

In this paper, a time fractional advection-dispersion equation is considered. From the relationship between the Caputo derivative and the Grunwald derivative, the Caputo derivative is approximated by using the Griinwald derivative. An implicit difference approximation for this equation is proposed. We prove that this approximation is unconditionally stable and convergent. Finally, numerical examples are given.