MIXED FINITE ELEMENT METHODS FOR FRACTIONAL NAVIER-STOKES EQUATIONS

This paper gives the detailed numerical analysis of mixed finite element method for fractional Navier-Stokes equations.The proposed method is based on the mixed finite element method in space and a finite difference scheme in time.The stability analyses of semi-discretization scheme and fully discrete scheme are discussed in detail.Furthermore,We give the convergence analysis for both semidiscrete and flly discrete schemes and then prove that the numerical solution converges the exact one with order O(h2+k),where h and k:respectively denote the space step size and the time step size.Finally,numerical examples are presented to demonstrate the effectiveness of our numerical methods.     

Gegenbauer spectral method for time‐fractional convection–diffusion equations with variable coefficients

In this paper, we study the numerical solution to time‐fractional partial differential equations with variable coefficients that involve temporal Caputo derivative. A spectral method based on Gegenbauer polynomials is taken for approximating the solution of the given time‐fractional partial differential equation in time and a collocation method in space. The suggested method reduces this type of equation to the solution of a linear algebraic system. Finally, some numerical examples are presented to illustrate the efficiency and accuracy of the proposed method. Copyright © 2014 John Wiley & Sons, Ltd.     

HIGH ORDER FINITE DIFFERENCE/SPECTRAL METHODS TO A WATER WAVE MODEL WITH NONLOCAL VISCOSITY

In this paper, efficient numerical schemes are proposed for solving the water wave model with nonlocal viscous term that describe the propagation of surface water wave. By using the Caputo fractional derivative definition to approximate the nonlocal fractional operator, finite difference method in time and spectral method in space are constructed for the considered model. The proposed method employs known 5/2 order scheme for fractional derivative and a mixed linearization for the nonlinear term. The analysis shows that the proposed numerical scheme is unconditionally stable and error estimates are provided to predict that the second order backward differentiation plus 5/2 order scheme converges with order 2 in time, and spectral accuracy in space. Several numerical results are provided to verify the efficiency and accuracy of our theoretical claims. Finally, the decay rate of solutions is investigated.     

Numerical Solution of the Two-Sided Space–Time Fractional Telegraph Equation Via Chebyshev Tau Approximation

The operational matrices of left Caputo fractional derivative, right Caputo fractional derivative, and Riemann–Liouville fractional integral, for shifted Chebyshev polynomials, are presented and derived. We propose an accurate and efficient spectral algorithm for the numerical solution of the two-sided space–time Caputo fractional-order telegraph equation with three types of non-homogeneous boundary conditions, namely, Dirichlet, Robin, and non-local conditions. The proposed algorithm is based on shifted Chebyshev tau technique combined with the derived shifted Chebyshev operational matrices. We focus primarily on implementing the novel algorithm both in temporal and spatial discretizations. This algorithm reduces the problem to a system of algebraic equations greatly simplifying the problem. This system can be solved by any standard iteration method. For confirming the efficiency and accuracy of the proposed scheme, we introduce some numerical examples with their approximate solutions and compare our results with those achieved using other methods.     

Fast evaluation and high accuracy finite element approximation for the time fractional subdiffusion equation

In this article, an efficient algorithm for the evaluation of the Caputo fractional derivative and the superconvergence property of fully discrete finite element approximation for the time fractional subdiffusion equation are considered. First, the space semidiscrete finite element approximation scheme for the constant coefficient problem is derived and supercloseness result is proved. The time discretization is based on the L1‐type formula, whereas the space discretization is done using, the fully discrete scheme is developed. Under some regularity assumptions, the superconvergence estimate is proposed and analyzed. Then, extension to the case of variable coefficients is also discussed. To reduce the computational cost, the fast evaluation scheme of the Caputo fractional derivative to solve the fractional diffusion equations is designed. Finally, numerical experiments are presented to support the theoretical results.     

Galerkin finite element method for nonlinear fractional Schrödinger equations

In this paper, a class of nonlinear Riesz space-fractional Schrödinger equations are considered. Based on the standard Galerkin finite element method in space and Crank-Nicolson difference method in time, the semi-discrete and fully discrete systems are constructed. By Brouwer fixed point theorem and fractional Gagliardo-Nirenberg inequality, we prove the fully discrete system is uniquely solvable. Moreover, we focus on a rigorous analysis and consideration of the conservation and convergence properties for the semi-discrete and fully discrete systems. Finally, a linearized iterative finite element algorithm is introduced and some numerical examples are given to confirm the theoretical results.     

分数阶Cable方程的有限点法分析

通过采用中心差分格式离散Riemann-Liouville时间分数阶导数和用有限点法建立离散代数系统,提出了数值求解分数阶Cable方程的无网格有限点法,详细推导了该方法的理论误差估计。数值算例证实了该方法的有效性和收敛性,并验证了理论分析结果。     

A fully discrete H 1-Galerkin method with quadrature for nonlinear advection–diffusion–reaction equations

We propose and analyze a fully discrete H ¹-Galerkin method with quadrature for nonlinear parabolic advection–diffusion–reaction equations that requires only linear algebraic solvers. Our scheme applied to the special case heat equation is a fully discrete quadrature version of the least-squares method. We prove second order convergence in time and optimal H ¹ convergence in space for the computer implementable method. The results of numerical computations demonstrate optimal order convergence of scheme in H ^k for k = 0, 1, 2. Support of the Australian Research Council is gratefully acknowledged.     

An Explicit Pseudo-Spectral Scheme with Almost Unconditional Stability for the Cahn-Hilliard Equation

In this paper, an explicit fully discrete three-level pseudo-spectral scheme with almost unconditional stability for the Cahn-Hilliard equation is proposed. Stability and convergence of the scheme are proved by Sobolev's inequalities and the bounded extensive method of the nonlinear function (B.N. Lu (1995)). The scheme possesses the almost same stable condition and convergent accuracy as the Creak-Nicloson scheme but it is an explicit scheme. Thus the iterative method to solve nonlinear algebraic system is avoided. Moreover, the linear stability of the critical point $u_0$ is investigated and the linear dispersive relation is obtained. Finally, the numerical results are supplied, which check the theoretical results.     

非局部Allen-Cahn模型全离散数值格式的误差估计

在这项工作中,我们研究了求解非局部体积守恒Allen-Cahn (AC)方程的全离散傅里叶伪谱数值格式的误差估计.该数值格式的时间行军方法基于著名的不变能量平方法(IEQ). 我们证明了所提出的全离散数值方法是唯一可解,无条件能量稳定的,并获得了该方案在空间和时间上的最优误差估计.此外,我们还进行了一些数值检验来验证理论结果.     

Galerkin-Legendre spectral schemes for nonlinear space fractional Schrödinger equation

In the paper, we first propose a Crank-Nicolson Galerkin-Legendre (CN-GL) spectral scheme for the one-dimensional nonlinear space fractional Schrödinger equation. Convergence with spectral accuracy is proved for the spectral approximation. Further, a Crank-Nicolson ADI Galerkin-Legendre spectral method for the two-dimensional nonlinear space fractional Schrödinger equation is developed. The proposed schemes are shown to be efficient with second-order accuracy in time and spectral accuracy in space which are higher than some recently studied methods. Moreover, some numerical results are demonstrated to justify the theoretical analysis.     

分数阶Cahn-Hilliard方程的高效数值算法

给出了时空分数阶Cahn-Hilliard方程的一个高效数值算法.首先,利用Laplace变换将时空分数阶Cahn-Hilliard方程转化为空间分数阶Cahn-Hilliard方程;然后,结合Fourier谱方法和有限差分法得到一个时间二阶、空间谱精度的高效数值格式;最后,通过数值实验验证本文数值算法的有效性,并验证其满足能量耗散性质和质量守恒定律.     

Convergence analysis of a fast second-order time-stepping numerical method for two-dimensional nonlinear time–space fractional Schrödinger equation

In this article, we consider the two-dimensional nonlinear time–space fractional Schrödinger equation with space described by the fractional Laplacian. A second-order fractional backward difference formula in the temporal direction while Fourier spectral method in the spatial direction is proposed to solve the model numerically. In the numerical implementation, a fast method is applied based on a globally uniform approximation of the trapezoidal rule for the integral on the real line to decrease the memory requirement and computational cost. By using the generalized discrete Gronwall inequality developed by Dixon and McKee and the temporal–spatial error splitting argument, the convergence of the fast time-stepping numerical method is also proved in a simple manner without imposing the Courant-Friedrichs-Lewy (CFL) condition. Finally, some numerical results are provided to support the theoretical analysis.     

非线性分数阶反应扩散方程组的间断时空有限元方法

利用时间间断空间连续的时空有限元方法构造了空间分数阶反应扩散方程组的可以逐时间层求解的全离散格式.在时间离散区间上,采用Radau积分公式,将插值理论与有限元理论相结合,给出了全离散格式解的存在唯一性结果,并证明了所给格式是无条件稳定的,进而详细给出最优阶L~∞(L~2)模误差估计过程.最后用数值算例验证了理论分析的正确性.     

An efficient linear second order unconditionally stable direct discretization method for the phase-field crystal equation on surfaces

We develop an unconditionally stable direct discretization scheme for solving the phase-field crystal equation on surfaces. The surface is discretized by using an unstructured triangular mesh. Gradient, divergence, and Laplacian operators are defined on triangular meshes. The proposed numerical method is second-order accurate in space and time. At each time step, the proposed computational scheme results in linear elliptic equations to be solved, thus it is easy to implement the algorithm. It is proved that the proposed scheme satisfies a discrete energy-dissipation law. Therefore, it is unconditionally stable. A fast and efficient biconjugate gradients stabilized solver is used to solve the resulting discrete system. Numerical experiments are conducted to demonstrate the performance of the proposed algorithm.     

A new numerical algorithm to solve fractional differential equations based on operational matrix of generalized hat functions

In this paper, we propose a new numerical algorithm for solving linear and non linear fractional differential equations based on our newly constructed integer order and fractional order generalized hat functions operational matrices of integration. The linear and nonlinear fractional order differential equations are transformed into a system of algebraic equations by these matrices and these algebraic equations are solved through known computational methods. Further some numerical examples are given to illustrate and establish the accuracy and reliability of the proposed algorithm. The results obtained, using the scheme presented here, are in full agreement with the analytical solutions and numerical results presented elsewhere.     

Application of the modified operational matrices in multiterm variable‐order time‐fractional partial differential equations

In this paper, we present a novel discrete scheme based on Genocchi polynomials and fractional Laguerre functions to solve multiterm variable‐order time‐fractional partial differential equations (M‐V‐TFPDEs) in the large interval. In this purpose, the accurate modified operational matrices are constructed to reduce the problems into a system of algebraic equations. Also, the computational algorithm based on the method and modified operational matrices in the large interval is easily implemented. Furthermore, we discuss the error estimation of the proposed method. Ultimately, to confirm our theoretical analysis and accuracy of numerical approach, several examples are presented.     

Locally linearized fractional step methods for nonlinear parabolic problems

This work deals with the efficient numerical solution of a class of nonlinear time-dependent reaction-diffusion equations. Via the method of lines approach, we first perform the spatial discretization of the original problem by applying a mimetic finite difference scheme. The system of ordinary differential equations arising from that process is then integrated in time with a linearly implicit fractional step method. For that purpose, we locally decompose the discrete nonlinear diffusion operator using suitable Taylor expansions and a domain decomposition splitting technique. The totally discrete scheme considers implicit time integrations for the linear terms while explicitly handling the nonlinear ones. As a result, the original problem is reduced to the solution of several linear systems per time step which can be trivially decomposed into a set of uncoupled parallelizable linear subsystems. The convergence of the proposed methods is illustrated by numerical experiments.     

A fast second‐order difference scheme for the space–time fractional equation

In this paper, a fast second‐order accurate difference scheme is proposed for solving the space–time fractional equation. The temporal Caputo derivative is approximated by ?L2 ‐1_σ formula which employs the sum‐of‐exponential approximation to the kernel function appeared in Caputo derivative. The second‐order linear spline approximation is applied to the spatial Riemann–Liouville derivative. At each time step, a fast algorithm, the preconditioned conjugate gradient normal residual method with a circulant preconditioner (PCGNR), is used to solve the resulting system that reduces the storage and computational cost significantly. The unique solvability and unconditional convergence of the difference scheme are shown by the discrete energy method. Numerical examples are given to verify numerical accuracy and efficiency of the difference schemes.     

An efficient numerical approach to solve Schrödinger equations with space fractional derivative

We design and analyze an efficient numerical approach to solve the coupled Schrödinger equations with space‐fractional derivative. The numerical scheme is based on leap‐frog in time direction and Fourier method in spatial direction. The advantage of the numerical scheme is that only a linear equation needs to be solved for each time step size, and we proved that the energy and mass of space‐fractional coupled Schrödinger equations (SFCSEs) are conserved in the case of full‐discrete scheme. Moreover, we also analyze the error estimate of the numerical scheme, and numerical solutions converge with the order in L² norm. Numerical examples are illustrated to verify the theoretical results.