Crank-Nicolson Quasi-Compact Scheme for the Nonlinear Two-Sided Spatial Fractional Advection-Diffusion Equations

The higher-order numerical scheme of nonlinear advection-diffusion equations is studied in this article, where the space fractional derivatives are evaluated by using weighted and shifted Grünwald difference operators and combining the compact technique, in the time direction is discretized by the Crank-Nicolson method. Through the energy method, the stability and convergence of the numerical scheme in the sense of L2-norm are proved, and the convergence order is . Some examples are given to show that our numerical scheme is effective.

ALTERNATING DIRECTION IMPLICIT OSC SCHEME FOR THE TWO-DIMENSIONAL FRACTIONAL EVOLUTION EQUATION WITH A WEAKLY SINGULAR KERNEL

In this paper, a new kind of alternating direction implicit （ADI） Crank-Nicolson-type orthogonal spline collocation （OSC） method is formulated for the two-dimensional frac-tional evolution equation with a weakly singular kernel arising in the theory of linear viscoelas-ticity. The novel OSC method is used for the spatial discretization, and ADI Crank-Nicolson-type method combined with the second order fractional quadrature rule are considered for thetemporal component. The stability of proposed scheme is rigourously established, and nearlyoptimal order error estimate is also derived. Numerical experiments are conducted to supportthe predicted convergence rates and also exhibit expected super-convergence phenomena.

GENERALIZED UPWIND SCHEME WITH FRACTIONAL STEPS FOR 3-D PROBLEM OF CONVECTION DOMINATING GROUNDWATER TRANSPORT

A generalized upwind scheme with fractional steps for 3-D mathematical models of convection dominating groundwater quality is suggested. The mass transport equation is split into a convection equation and a dispersive equation. The generalized upwind scheme is used to solve the convection equation and the finite element method is used to compute the dispersive equation.These procedures which not only overcome the phenomenon of the negative concentration and numerical dispersion appear frequently with normal FEM or FDM to solve models of convection dominating groundwater transport but also avoid the step for computing each node velocity give a more suitable method to calculate the concentrations of the well points.

Control Scheme Discussion for Bottom Oil Slurry System of FCCU Fractionator

流化催化裂化装置（FCCU）分馏塔底油浆系统的控制效果对整个装置的平稳安全运行具有重要意义。近年来国内新建了许多大型FCCU，既包含国内各设计院自己的技术，也有部分引进国外的工艺包，其分馏塔底油浆系统的控制方案也有多种形式。总结了近年来国内新建FCCU常用的三种较为典型的油浆系统控制方案，阐述了各方案的特点，对其进行了简要的对比和分析。在此基础上尝试性地提出了这三种方案的适用场合，方便设计人员根据每个项目的特点设计合理的控制方案，达到最佳控制效果。

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.

Compact finite difference schemes for the backward fractional Feynman–Kac equation with fractional substantial derivative

The fractional Feynman-Kac equations describe the distributions of functionals of non-Brownian motion, or anomalous diffusion, including two types called the forward and backward fractional Feynman-Kac equations, where the nonlocal time-space coupled fractional substantial derivative is involved. This paper focuses on the more widely used backward version. Based on the newly proposed approximation operators for fractional substantial derivative, we establish compact finite difference schemes for the backward fractional Feynman-Kac equation. The proposed difference schemes have the q-th(q = 1, 2, 3, 4) order accuracy in temporal direction and fourth order accuracy in spatial direction, respectively. The numerical stability and convergence in the maximum norm are proved for the first order time discretization scheme by the discrete energy method, where an inner product in complex space is introduced. Finally, extensive numerical experiments are carried out to verify the availability and superiority of the algorithms. Also, simulations of the backward fractional Feynman-Kac equation with Dirac delta function as the initial condition are performed to further confirm the effectiveness of the proposed methods.

A Compact Finite Volume Scheme for the Multi-Term Time Fractional Sub-Diffusion Equation

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∞-norm. The convergence order is O(τ2-α + h4). Finally, two numerical examples are given to confirm the theoretical results. Some tables listed also can explain the stability and convergence of the scheme.

基于L2-1_(σ)格式逼近时间分数阶扩散方程的差分方法及其收敛性分析

针对时间分数阶扩散方程,在时间方向上结合L2-1_(σ)格式,空间上采用二阶中心差分方法进行离散,并对离散格式进行了收敛性和稳定性分析,离散格式和分析方法可以很容易推广到空间高维情形。最后,通过数值算例对L2-1_(σ)格式和L1格式进行了误差和收敛阶的对比,显示出L2-1_(σ)格式在时间分数阶导数逼近上的优势。

空间分数阶KGS方程组的辛差分格式

该文研究分数阶KGS方程组的辛差分格式.首先,作者给出了无穷维分数阶Hamilton系统,并将KGS方程组转化为Hamilton系统.然后,基于分数阶中心差分格式对分数阶KGS方程组进行空间离散,得到的半离散系统是一个有限维Hamilton系统.接着,利用辛中点格式对时间进行离散得到全离散格式,并且对该格式进行了守恒性分析.最后,通过数值实验验证了该数值格式的有效性.

多项时间分数阶扩散方程类Carey非协调元的误差分析

基于L^(1)全离散格式,针对具有Caputo导数的二维多项时间分数阶扩散方程给出了类Carey非协调有限元方法.利用该单元的特殊性质和分数阶导数巧妙的处理技巧导出了L^(2)模和H^(1)模意义下的最优误差估计.

Discrete Chaos in Fractional Henon Map

In this study, a discrete fractional Henon map is proposed in the Caputo discrete delta’s sense. The results show that the discrete fractional calculus is an efficient tool and the maps derived in this way have simpler forms but hold rich dynamical behaviors.

初值奇异性非线性分数阶常微分方程的高阶数值方法

考虑非线性分数阶常微方程高阶格式的精确解具有初值奇异性,从而引入初值变量和逐块方法,再利用拉格朗日插值公式,提出一种新的高阶数值格式。该高阶数值格式为非光滑解条件下的5+α阶。

Investigation of Uncertainties of Establishment Schemes in Dynamic Global Vegetation Models

In Dynamic Global Vegetation Models （DGVMs）, the establishment of woody vegetation refers to flowering, fertiliza- tion, seed production, germination, and the growth of tree seedlings. It determines not only the population densities but also other important ecosystem structural variables. In current DGVMs, establishments of woody plant functional types （PFTs） are assumed to be either the same in the same grid cell, or largely stochastic. We investigated the uncertainties in the competition of establishment among coexisting woody PFTs from three aspects： the dependence of PFT establishments on vegetation states; background establishment; and relative establishment potentials of different PFTs. Sensitivity experi- ments showed that the dependence of establishment rate on the fractional coverage of a PFT favored the dominant PFT by increasing its share in establishment. While a small background establishment rate had little impact on equilibrium states of the ecosystem, it did change the timescale required for the establishment of alien species in pre-existing forest due to their disadvantage in seed competition during the early stage of invasion. Meanwhile, establishment purely fiom background （the scheme commonly used in current DGVMs） led to inconsistent behavior in response to the change in PFT specification （e.g., number of PFTs and their specification）. Furthermore, the results also indicated that trade-off between irtdividual growth and reproduction/colonization has significant influences on the competition of establishment. Hence, further development of es- tablishment parameterization in DGVMs is essential in reducing the uncertainties in simulations of both ecosystem structures and successions.

Development of an Establishment Scheme for a DGVM

Environmental changes are expected to shift the distribution and abundance of vegetation by determining seedling estab- lishment and success. However, most current ecosystem models only focus on the impacts of abiotic factors on biogeophysics （e.g., global distribution, etc.）, ignoring their roles in the population dynamics （e.g., seedling establishment rate, mortality rate, etc.） of ecological communities. Such neglect may lead to biases in ecosystem population dynamics （such as changes in population density for woody species in forest ecosystems） and characteristics. In the present study, a new establishment scheme for introducing soil water as a function rather than a threshold was developed and validated, using version 1.0 of the IAP-DGVM as a test bed. The results showed that soil water in the establishment scheme had a remarkable influence on forest transition zones. Compared with the original scheme, the new scheme significantly improved simulations of tree population density, especially in the peripheral areas of forests and transition zones. Consequently, biases in forest fractional coverage were reduced in approximately 78.8% of the global grid cells. The global simulated areas of tree, shrub, grass and bare soil performed better, where the relative biases were reduced from 34.3% to 4.8%, from 27.6% to 13.1%, from 55.2% to 9.2%, and from 37.6% to 3.6%, respectively. Furthermore, the new scheme had more reasonable dependencies of plant functional types （PFTs） on mean annual precipitation, and described the correct dominant PFTs in the tropical rainforest peripheral areas of the Amazon and central Africa.

A Class of Semi-Implicit Parallel Difference Method for Time Fractional Diffusion Equations

In this paper, we construct a class of semi-implicit difference method for time fractional diffusion equations—the group explicit (GE) difference scheme, which is a difference scheme with good parallelism constructed using Saul’yev asymmetric scheme. The stability and convergence of the GE scheme of time fractional diffusion equation are analyzed by mathematical induction. Then, the theoretical analysis is verified by numerical experiments, which shows that the GE scheme is effective for solving the time fractional diffusion equation.

High-order numerical methods of fractional-order Stokes' first problem for heated generalized second grade fluid

The high-order implicit finite difference schemes for solving the fractional- order Stokes＇ first problem for a heated generalized second grade fluid with the Dirichlet boundary condition and the initial condition are given. The stability, solvability, and convergence of the numerical scheme are discussed via the Fourier analysis and the matrix analysis methods. An improved implicit scheme is also obtained. Finally, two numerical examples are given to demonstrate the effectiveness of the mentioned schemes

Square-Root Dynamics of a SIR-Model in Fractional Order

In this paper, we consider an SIR-model for which the interaction term is the square root of the susceptible and infected individuals in the form of fractional order differential equations. First the non-negative solution of the model in fractional order is presented. Then the local stability analysis of the model in fractional order is discussed. Finally, the general solutions are presented and a discrete-time finite difference scheme is constructed using the nonstandard finite difference (NSFD) method. A comparative study of the classical Runge-Kutta method and ODE45 is presented in the case of integer order derivatives. The solutions obtained are presented graphically.

时间分数阶扩散方程有限体积法的隐式差分格式

对于时间Caputo型导数的扩散方程,根据Caputo型导数和Grunwald-Letnikov型导数的关系,利用Grunwald-Letnikov型导数的离散格式离散分数阶导数,构造有限体积法的隐式差分格式,并证明差分格式的无条件稳定性和无条件收敛性,并给出数值例子。

Explicit Approximation Solutions and Proof of Convergence of the Space-Time Fractional Advection Dispersion Equations

The space-time fractional advection dispersion equations are linear partial pseudo-differential equations with spatial fractional derivatives in time and in space and are used to model transport at the earth surface. The time fractional order is denoted by β∈ and ?is devoted to the space fractional order. The time fractional advection dispersion equations describe particle motion with memory in time. Space-fractional advection dispersion equations arise when velocity variations are heavy-tailed and describe particle motion that accounts for variation in the flow field over entire system. In this paper, I focus on finding the precise explicit discrete approximate solutions to these models for some values of ?with ?, ?while the Cauchy case as ?and the classical case as ?with ?are studied separately. I compare the numerical results of these models for different values of ?and ?and for some other related changes. The approximate solutions of these models are also discussed as a random walk with or without a memory depending on the value of . Then I prove that the discrete solution in the Fourierlaplace space of theses models converges in distribution to the Fourier-Laplace transform of the corresponding fractional differential equations for all the fractional values of ?and .

分数阶多延迟抛物方程不同差分格式的分析

为获得分数阶多延迟抛物方程精确的解析表达式及其数值解,结合分数阶微积分的定义和数值方法,对分数阶多延迟抛物方程构造一类线性化的Crank-Nicolson差分格式和紧致差分格式.通过数值算例对差分格式的可解性、稳定性和收敛性进行验证.结果表明,该分数阶多延迟抛物方程的Crank-Nicolson差分格式和紧致差分格式具有良好的精确性和有效性.