A New Unconditionally Stable Method for Telegraph Equation Based on Associated Hermite Orthogonal Functions

The present paper proposes a new unconditionally stable method to solve telegraph equation by using associated Hermite (AH) orthogonal functions. Unlike other numerical approaches, the time variables in the given equation can be handled analytically by AH basis functions. By using the Galerkin’s method, one can eliminate the time variables from calculations, which results in a series of implicit equations. And the coefficients of results for all orders can then be obtained by the expanded equations and the numerical results can be reconstructed during the computing process. The precision and stability of the proposed method are proved by some examples, which show the numerical solution acquired is acceptable when compared with some existing methods

High-order Compact Difference Schemes for the Modified Anomalous Subdiffusion Equation

In this paper, two kinds of high-order compact finite difference schemes for second-order derivative are developed. Then a second-order numerical scheme for Riemann-Liouvile derivative is established based on fractional center difference operator. We apply these methods to fractional anomalous subdiffusion equation to construct two kinds of novel numerical schemes. The solvability, stability and convergence analysis of these difference schemes are studied by Fourier method in details. The convergence orders of these numerical schemes are $\mathcal {O}(\tau^2+h^6)$ and $\mathcal {O}(\tau^2+h^8)$, respectively. Finally, numerical experiments are displayed which are in line with the theoretical analysis.Comment:

Approximate Analytical Solutions of Space-Fractional Telegraph Equations by Sumudu Adomian Decomposition Method

The main goal in this work is to establish a new and efficient analytical scheme for space fractional telegraph equation (FTE) by means of fractional Sumudu decomposition method (SDM). The fractional SDM gives us an approximate convergent series solution. The stability of the analytical scheme is also studied. The approximate solutions obtained by SDM show that the approach is easy to implement and computationally very much attractive. Further, some numerical examples are presented to illustrate the accuracy and stability for linear and nonlinear cases

High-Order Algorithms for Riesz Derivative and Their Applications (

We firstly develop the high-order numerical algorithms for the left and right Riemann-Liouville derivatives. Using these derived schemes, we can get high-order algorithms for the Riesz fractional derivative. Based on the approximate algorithm, we construct the numerical scheme for the space Riesz fractional diffusion equation, where a fourth-order scheme is proposed for the spacial Riesz derivative, and where a compact difference scheme is applied to approximating the first-order time derivative. It is shown that the difference scheme is unconditionally stable and convergent. Finally, numerical examples are provided which are in line with the theoretical analysis

B-Splines Based Finite Difference Schemes For Fractional Partial Differential Equations

Fractional partial differential equations (FPDEs) are considered to be the extended formulation of classical partial differential equations (PDEs). Several physical models in certain fields of sciences and engineering are more appropriately formulated in the form of FPDEs. FPDEs in general, do not have exact analytical solutions. Thus, the need to develop new numerical methods for the solutions of space and time FPDEs. This research focuses on the development of new numerical methods. Two methods based on B-splines are developed to solve linear and non-linear FPDEs. The methods are extended cubic B-spline approximation (ExCuBS) and new extended cubic B-spline approximation (NExCuBS). Both methods have the same basis functions but for the NExCuBS, a new approximation is used for the second order space derivative

A finite difference method for fractional diffusion equations with Neumann boundary conditions

A finite difference numerical method is investigated for fractional order diffusion problems in one space dimension. The basis of the mathematical model and the numerical approximation is an appropriate extension of the initial values, which incorporates homogeneous Dirichlet or Neumann type boundary conditions. The well-posedness of the obtained initial value problem is proved and it is pointed out that each extensions is compatible with the original boundary conditions. Accordingly, a finite difference scheme is constructed for the Neumann problem using the shifted Grünwald--Letnikov approximation of the fractional order derivatives, which is based on infinite many basis points. The corresponding matrix is expressed in a closed form and the convergence of an appropriate implicit Euler scheme is proved