Space‐time spectral collocation method for the one‐dimensional sine‐Gordon equation

In this article, we introduce a new space‐time spectral collocation method for solving the one‐dimensional sine‐Gordon equation. We apply a spectral collocation method for discretizing spatial derivatives, and then use the spectral collocation method for the time integration of the resulting nonlinear second‐order system of ordinary differential equations (ODE). Our formulation has high‐order accurate in both space and time. Optimal a priori error bounds are derived in the L²‐norm for the semidiscrete formulation. Numerical experiments show that our formulation have exponential rates of convergence in both space and time. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 670–690, 2015     

Fractional stochastic Burgers‐type equation in Hölder space—Wellposedness and approximations

In this work, we use the spectral Galerkin method to prove the existence of a pathwise unique mild solution of a fractional stochastic partial differential equation of Burgers type in a Hölder space. We get the temporal regularity, and using a combination of Galerkin and exponential‐Euler methods, we obtain a full discretization scheme of the solution. Moreover, we calculate the rates of convergence for both approximations (Galerkin and full discretization) with respect to time and to space.     

Space‐time spectral method for two‐dimensional semilinear parabolic equations

In this paper, a high‐order accurate numerical method for two‐dimensional semilinear parabolic equations is presented. We apply a Galerkin–Legendre spectral method for discretizing spatial derivatives and a spectral collocation method for the time integration of the resulting nonlinear system of ordinary differential equations. Our formulation can be made arbitrarily high‐order accurate in both space and time. Optimal a priori error bound is derived in the L²‐norm for the semidiscrete formulation. Extensive numerical results are presented to demonstrate the convergence property of the method, show our formulation have spectrally accurate in both space and time. John Wiley & Sons, Ltd.     

Finite difference spectral approximation for the time‐space fractional telegraph equation and its parameter estimation

The object of this paper is to present the numerical solution of the time‐space fractional telegraph equation. The proposed method is based on the finite difference scheme in temporal direction and Fourier spectral method in spatial direction. The fast Fourier transform (FFT) technique is applied to practical computation. The stability and convergence analysis are strictly proven, which shows that this method is stable and convergent with (2?α) order accuracy in time and spectral accuracy in space. Moreover, the Levenberg‐Marquardt (L‐M) iterative method is employed for the parameter estimation. Finally, some numerical examples are given to confirm the theoretical analysis.     

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.     

Analytical and numerical solutions of a two‐dimensional multi‐term time‐fractional Oldroyd‐B model

In this paper, we consider a two‐dimensional multi‐term time‐fractional Oldroyd‐B equation on a rectangular domain. Its analytical solution is obtained by the method of separation of variables. We employ the finite difference method with a discretization of the Caputo time‐fractional derivative to obtain an implicit difference approximation for the equation. Stability and convergence of the approximation scheme are established in the L_∞ ‐norm. Two examples are given to illustrate the theoretical analysis and analytical solution. The results indicate that the present numerical method is effective for this general two‐dimensional multi‐term time‐fractional Oldroyd‐B model.     

Convergence analysis for parallel‐in‐time solution of hyperbolic systems

Parallel‐in‐time algorithms have been successfully employed for reducing time‐to‐solution of a variety of partial differential equations, especially for diffusive (parabolic‐type) equations. A major failing of parallel‐in‐time approaches to date, however, is that most methods show instabilities or poor convergence for hyperbolic problems. This paper focuses on the analysis of the convergence behavior of multigrid methods for the parallel‐in‐time solution of hyperbolic problems. Three analysis tools are considered that differ, in particular, in the treatment of the time dimension: (a) space–time local Fourier analysis, using a Fourier ansatz in space and time; (b) semi‐algebraic mode analysis, coupling standard local Fourier analysis approaches in space with algebraic computation in time; and (c) a two‐level reduction analysis, considering error propagation only on the coarse time grid. In this paper, we show how insights from reduction analysis can be used to improve feasibility of the semi‐algebraic mode analysis, resulting in a tool that offers the best features of both analysis techniques. Following validating numerical results, we investigate what insights the combined analysis framework can offer for two model hyperbolic problems, the linear advection equation in one space dimension and linear elasticity in two space dimensions.     

Gauss‐Lobatto‐Legendre‐Birkhoff pseudospectral approximations for the multi‐term time fractional diffusion‐wave equation with Neumann boundaryconditions

A second‐order finite difference/pseudospectral scheme is proposed for numerical approximation of multi‐term time fractional diffusion‐wave equation with Neumann boundary conditions. The scheme is based upon the weighted and shifted Grünwald difference operators approximation of the time fractional calculus and Gauss‐Lobatto‐Legendre‐Birkhoff (GLLB) pseudospectral method for spatial discretization. The unconditionally stability and convergence of the scheme are rigorously proved. Numerical examples are carried out to verify theoretical results.     

Fractional quasi AKNS‐technique for nonlinear space–time fractional evolution equations

This paper aims to formulate the fractional quasi‐inverse scattering method. Also, we give a positive answer to the following question: can the Ablowitz‐Kaup‐Newell‐Segur (AKNS) method be applied to the space–time fractional nonlinear differential equations? Besides, we derive the Bäcklund transformations for the fractional systems under study. Also, we construct the fractional quasi‐conservation laws for the considered fractional equations from the defined fractional quasi AKNS‐like system. The nonlinear fractional differential equations to be studied are the space–time fractional versions of the Kortweg‐de Vries equation, modified Kortweg‐de Vries equation, the sine‐Gordon equation, the sinh‐Gordon equation, the Liouville equation, the cosh‐Gordon equation, the short pulse equation, and the nonlinear Schrödinger equation.     

A new numerical approach for single rational soliton solution of Chen‐Lee‐Liu equation with Riesz fractional derivative in optical fibers

In this paper, time‐splitting spectral approximation technique has been proposed for Chen‐Lee‐Liu (CLL) equation involving Riesz fractional derivative. The proposed numerical technique is efficient, unconditionally stable, and of second‐order accuracy in time and of spectral accuracy in space. Moreover, it conserves the total density in the discretized level. In order to examine the results, with the aid of weighted shifted Grünwald‐Letnikov formula for approximating Riesz fractional derivative, Crank‐Nicolson weighted and shifted Grünwald difference (CN‐WSGD) method has been applied for Riesz fractional CLL equation. The comparison of results reveals that the proposed time‐splitting spectral method is very effective and simple for obtaining single soliton numerical solution of Riesz fractional CLL equation.     

Numerical solutions of the space‐time fractional advection‐dispersion equation

Fractional advection‐dispersion equations are used in groundwater hydrologhy to model the transport of passive tracers carried by fluid flow in a porous medium. In this paper we present two reliable algorithms, the Adomian decomposition method and variational iteration method, to construct numerical solutions of the space‐time fractional advection‐dispersion equation in the form of a rabidly convergent series with easily computable components. The fractional derivatives are described in the Caputo sense. Some examples are given. Numerical results show that the two approaches are easy to implement and accurate when applied to space‐time fractional advection‐dispersion equations. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008     

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.     

High‐order difference scheme for the solution of linear time fractional klein–gordon equations

In this article, we apply a high‐order difference scheme for the solution of some time fractional partial differential equations (PDEs). The time fractional Cattaneo equation and the linear time fractional Klein–Gordon and dissipative Klein–Gordon equations will be investigated. The time fractional derivative which has been described in the Caputo's sense is approximated by a scheme of order , and the space derivative is discretized with a fourth‐order compact procedure. We will prove the solvability of the proposed method by coefficient matrix property and the unconditional stability and ‐convergence with the energy method. Numerical examples demonstrate the theoretical results and the high accuracy of the proposed scheme. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1234–1253, 2014     

Local radial basis function interpolation method to simulate 2D fractional‐time convection‐diffusion‐reaction equations with error analysis

In this article, a kind of meshless local radial point interpolation (MLRPI) method is proposed to two‐dimensional fractional‐time convection‐diffusion‐reaction equations and satisfactory agreements are archived. This method is based on meshless methods and benefits from collocation ideas but it does not belong to the traditional global meshless collocation methods. In MLRPI method, it does not need any kind of integration locally or globally over small quadrature domains which is essential in the finite element method and those meshless methods based on Galerkin weak form. Also, it is not needed to determine shape parameter which plays important role in collocation method based on the radial basis functions (Kansa's method). Therefore, computational costs of this kind of MLRPI method is less expensive. The stability and convergence of this meshless approach are discussed and theoretically proven. It is proved that the present meshless formulation is very effective for modeling and simulation of fractional differential equations. Furthermore, the numerical studies on sensitivity analysis and convergence analysis show the stability and reliable rates of convergence. © 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 974–994, 2017     

Finite‐time flocking problem of a Cucker–smale‐type self‐propelled particle model

In this article, we propose a Cucker–Smale‐type self‐propelled particle model with continuous non‐Lipschitz protocol. We show that the flocking can occur in finite time if the communication rate function satisfies a lower bound condition. Both our theoretical and numerical results uncover a power‐law relationship between the convergence time and the number of individuals. Our result implies that the individuals in groups with high density can transit rapidly to ordered collective motion. We also investigate the influence of control parameter on the convergence speed. © 2015 Wiley Periodicals, Inc. Complexity 21: 354–361, 2016     

Stability analysis for fractional‐order partial differential equations by means of space spectral time Adams‐Bashforth Moulton method

In this article, a new numerical scheme space Spectral time Fractional Adam Bashforth Moulton method for the solution of fractional partial differential equations is offered. The proposed method is obtained by modifying, in a suitable way; the spectral technique and the method of lines. The attention is focused on the stability properties and hence an elegant stability analysis for the current approach is also provided. Finally, two examples are presented to illustrate the effectiveness of the reported method. Obtained results confirm the convergence and spectral accuracy of the proposed method in both space and time. In addition, a comparison with the existing studies is also made as a limiting case of the considered problem at the end and found in good agreement.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 34: 19–29, 2018     

A high‐order B‐spline collocation scheme for solving a nonhomogeneous time‐fractional diffusion equation

A high‐accuracy numerical approach for a nonhomogeneous time‐fractional diffusion equation with Neumann and Dirichlet boundary conditions is described in this paper. The time‐fractional derivative is described in the sense of Riemann‐Liouville and discretized by the backward Euler scheme. A fourth‐order optimal cubic B‐spline collocation (OCBSC) method is used to discretize the space variable. The stability analysis with respect to time discretization is carried out, and it is shown that the method is unconditionally stable. Convergence analysis of the method is performed. Two numerical examples are considered to demonstrate the performance of the method and validate the theoretical results. It is shown that the proposed method is of order O(Δx⁴ + Δt^{2 ? α}) convergence, where α ∈ (0,1) . Moreover, the impact of fractional‐order derivative on the solution profile is investigated. Numerical results obtained by the present method are compared with those obtained by the method based on standard cubic B‐spline collocation method. The CPU time for present numerical method and the method based on cubic B‐spline collocation method are provided.     

Asymptotic behavior of solutions to space‐time fractional diffusion‐reaction equations

This article discusses the analyticity and the long‐time asymptotic behavior of solutions to space‐time fractional diffusion‐reaction equations in . By a Laplace transform argument, we prove that the decay rate of the solution as t→∞ is dominated by the order of the time‐fractional derivative. We consider the decay rate also in a bounded domain. Copyright © 2016 John Wiley & Sons, Ltd.     

Finite element method for two‐dimensional time‐fractional tricomi‐type equations

In this article, we consider the finite element method (FEM) for two‐dimensional linear time‐fractional Tricomi‐type equations, which is obtained from the standard two‐dimensional linear Tricomi‐type equation by replacing the first‐order time derivative with a fractional derivative (of order α, with 1 <α< 2 ). The method is based on finite element method for space and finite difference method for time. We prove that the method is unconditionally stable, and the error estimate is presented. The comparison of the FEM results with the exact solutions is made, and numerical experiments reveal that the FEM is very effective. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013