A high‐order compact scheme for the nonlinear fractional Klein–Gordon equation

In this article, a high‐order finite difference scheme for a kind of nonlinear fractional Klein–Gordon equation is derived. The time fractional derivative is described in the Caputo sense. The solvability of the difference system is discussed by the Leray–Schauder fixed point theorem, while the stability and L_∞ convergence of the finite difference scheme are proved by the energy method. Numerical examples are provided to demonstrate the theoretical results. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 706–722, 2015     

An efficient numerical approach to solve a class of variable‐order fractional integro‐partial differential equations

The main purpose of this work is to investigate an initial boundary value problem related to a suitable class of variable order fractional integro‐partial differential equations with a weakly singular kernel. To discretize the problem in the time direction, a finite difference method will be used. Then, the Sinc‐collocation approach combined with the double exponential transformation is employed to solve the problem in each time level. The proposed numerical algorithm is completely described and the convergence analysis of the numerical solution is presented. Finally, some illustrative examples are given to demonstrate the pertinent features of the proposed algorithm.     

A new approach by two‐dimensional wavelets operational matrix method for solving variable‐order fractional partial integro‐differential equations

In this paper, a new computational scheme based on operational matrices (OMs) of two‐dimensional wavelets is proposed for the solution of variable‐order (VO) fractional partial integro‐differential equations (PIDEs). To accomplish this method, first OMs of integration and VO fractional derivative (FD) have been derived using two‐dimensional Legendre wavelets. By implementing two‐dimensional wavelets approximations and the OMs of integration and variable‐order fractional derivative (VO‐FD) along with collocation points, the VO fractional partial PIDEs are reduced into the system of algebraic equations. In addition to this, some useful theorems are discussed to establish the convergence analysis and error estimate of the proposed numerical technique. Furthermore, computational efficiency and applicability are examined through some illustrative examples.     

Long time error analysis of the fourth‐order compact finite difference methods for the nonlinear Klein–Gordon equation with weak nonlinearity

We present the fourth‐order compact finite difference (4cFD) discretizations for the long time dynamics of the nonlinear Klein–Gordon equation (NKGE), while the nonlinearity strength is characterized by ?^p with a constant p ∈ ?⁺ and a dimensionless parameter ? ∈ (0, 1] . Based on analytical results of the life‐span of the solution, rigorous error bounds of the 4cFD methods are carried out up to the time at O(?^?p) . We pay particular attention to how error bounds depend explicitly on the mesh size h and time step τ as well as the small parameter ? ∈ (0, 1] , which indicate that, in order to obtain ‘correct’ numerical solutions up to the time at O(?^?p) , the ? ‐scalability (or meshing strategy requirement) of the 4cFD methods should be taken as: h = O(?^p/4) and τ = O(?^p/2) . It has better spatial resolution capacity than the classical second order central difference methods. By a rescaling in time, it is equivalent to an oscillatory NKGE whose solution propagates waves with wavelength at O(1) in space and O(?^p) in time. It is straightforward to get the error bounds of the oscillatory NKGE in the fixed time. Finally, numerical results are provided to confirm our theoretical analysis.     

An accurate numerical method for solving the linear fractional Klein–Gordon equation

In this article, an implementation of an efficient numerical method for solving the linear fractional Klein–Gordon equation (LFKGE) is introduced. The fractional derivative is described in the Caputo sense. The method is based upon a combination between the properties of the Chebyshev approximations and finite difference method (FDM). The proposed method reduces LFKGE to a system of ODEs, which is solved using FDM. Special attention is given to study the convergence analysis and deduce an error upper bound of the proposed method. Numerical example is given to show the validity and the accuracy of the proposed algorithm. Copyright © 2013 John Wiley & Sons, Ltd.     

Numerical simulation of the nonlinear generalized time‐fractional Klein–Gordon equation using cubic trigonometric B‐spline functions

In this paper, an efficient numerical procedure for the generalized nonlinear time‐fractional Klein–Gordon equation is presented. We make use of the typical finite difference schemes to approximate the Caputo time‐fractional derivative, while the spatial derivatives are discretized by means of the cubic trigonometric B‐splines. Stability and convergence analysis for the numerical scheme are discussed. We apply our scheme to some typical examples and compare the obtained results with the ones found by other numerical methods. The comparison shows that our scheme is quite accurate and can be applied successfully to a variety of problems of applied nature.     

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.     

A novel attempt for finding comparatively accurate solution for sine‐Gordon equation comprising Riesz space fractional derivative

In this paper, a numerical procedure involving Chebyshev wavelet method has been implemented for computing the approximate solution of Riesz space fractional sine‐Gordon equation (SGE). Two‐dimensional Chebyshev wavelet method is implemented to calculate the numerical solution of space fractional SGE. The fractional SGE is considered as an interpolation between the classical SGE (corresponding to α = 2) and nonlocal SGE (corresponding to α = 1). As a consequence, the approximate solutions of fractional SGE obtained by using Chebyshev wavelet approach were compared with those derived by using modified homotopy analysis method with Fourier transform. Copyright © 2015 John Wiley & Sons, Ltd.     

Fractional‐order orthogonal Bernstein polynomials for numerical solution of nonlinear fractional partial Volterra integro‐differential equations

In this paper, a new two‐dimensional fractional polynomials based on the orthonormal Bernstein polynomials has been introduced to provide an approximate solution of nonlinear fractional partial Volterra integro‐differential equations. For this aim, the fractional‐order orthogonal Bernstein polynomials (FOBPs) are constructed, and its operational matrices of integration, fractional‐order integration, and derivative in the Caputo sense and product operational matrix are derived. These operational matrices are utilized to reduce the under study problem to a nonlinear system of algebraic equations. Using the approximation of FOBPs, the convergence analysis and error estimate associated to the proposed problem have been investigated. Finally, several examples are included to clarify the validity, efficiency, and applicability of the proposed technique via FOBPs approximation.     

A numerical method for computing radially symmetric solutions of a dissipative nonlinear modified Klein‐Gordon equation

In this article we develop a finite‐difference scheme to approximate radially symmetric solutions of a dissipative nonlinear modified Klein‐Gordon equation subject to smooth initial conditions ? and ψ in an open sphere D around the origin, with constant internal and external damping coefficients—β and γ, respectively—, and nonlinear term of the form G′(w) = w^p, with p > 1 an odd number. The functions ? and ψ are radially symmetric in D, and ?, ψ, r?, and rψ are assumed to be small at infinity. We prove that our scheme is consistent order ??(Δt²) + ??(Δr²) for G′ identically equal to zero and provide a necessary condition for it to be stable order n. Part of our study will be devoted to compare the physical effects of β and γ. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Numerical solutions of nonlinear fractional Schrödinger equations using nonstandard discretizations

In this article, numerical study for both nonlinear space‐fractional Schrödinger equation and the coupled nonlinear space‐fractional Schrödinger system is presented. We offer here the weighted average nonstandard finite difference method (WANSFDM) as a novel numerical technique to study such kinds of partial differential equations. The space fractional derivative is described in the sense of the quantum Riesz‐Feller definition. Stability analysis of the proposed method is studied. To show that this method is reliable and computationally efficient different numerical examples are provided. We expect that the proposed schemes can be applicable to different systems of fractional partial differential equations. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1399–1419, 2017     

A numerical method for the one‐dimensional sine‐Gordon equation

A numerical method based on a predictor–corrector (P‐C) scheme arising from the use of rational approximants of order 3 to the matrix‐exponential term in a three‐time level recurrence relation is applied successfully to the one‐dimensional sine‐Gordon equation, already known from the bibliography. In this P‐C scheme a modification in the corrector (MPC) has been proposed according to which the already evaluated corrected values are considered. The method, which uses as predictor an explicit finite‐difference scheme arising from the second order rational approximant and as corrector an implicit one, has been tested numerically on the single and the soliton doublets. Both the predictor and the corrector schemes are analyzed for local truncation error and stability. From the investigation of the numerical results and the comparison of them with other ones known from the bibliography it has been derived that the proposed P‐C/MPC schemes at least coincide in terms of accuracy with them. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Stability analysis and a numerical scheme for fractional Klein‐Gordon equations

Fractional order nonlinear Klein‐Gordon equations (KGEs) have been widely studied in the fields like; nonlinear optics, solid state physics, and quantum field theory. In this article, with help of the Sumudu decomposition method (SDM), a numerical scheme is developed for the solution of fractional order nonlinear KGEs involving the Caputo's fractional derivative. The coupled method provides us very efficient numerical scheme in terms of convergent series. The iterative scheme is applied to illustrative examples for the demonstration and applications.     

Conservative compact finite difference scheme for the N‐coupled nonlinear Klein–Gordon equations

In this article, a compact finite difference method is developed for the periodic initial value problem of the N‐coupled nonlinear Klein–Gordon equations. The present scheme is proved to preserve the total energy in the discrete sense. Due to the difficulty in obtaining the priori estimate from the discrete energy conservation law, the cut‐off function technique is employed to prove the convergence, which shows the new scheme possesses second order accuracy in time and fourth order accuracy in space, respectively. Additionally, several numerical results are reported to confirm our theoretical analysis. Lastly, we apply the reliable method to simulate and study the collisions of solitary waves numerically.     

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     

On the numerical solution of the Klein‐Gordon equation

A predictor–corrector (P–C) scheme based on the use of rational approximants of second‐order to the matrix‐exponential term in a three‐time level reccurence relation is applied to the nonlinear Klein‐Gordon equation. This scheme is accelerated by using a modification (MPC) in which the already evaluated values are used for the corrector. Both the predictor and the corrector scheme are analyzed for local truncation error and stability. The proposed method is applied to problems possessing periodic, kinks and single, double‐soliton waves. The accuracy as well as the long time behavior of the proposed scheme is discussed. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

A multigrid compact finite difference method for solving the one‐dimensional nonlinear sine‐Gordon equation

The aim of this paper is to propose a multigrid method to obtain the numerical solution of the one‐dimensional nonlinear sine‐Gordon equation. The finite difference equations at all interior grid points form a large sparse linear system, which needs to be solved efficiently. The solution cost of this sparse linear system usually dominates the total cost of solving the discretized partial differential equation. The proposed method is based on applying a compact finite difference scheme of fourth‐order for discretizing the spatial derivative and the standard second‐order central finite difference method for the time derivative. The proposed method uses the Richardson extrapolation method in time variable. The obtained system has been solved by V‐cycle multigrid (VMG) method, where the VMG method is used for solving the large sparse linear systems. The numerical examples show the efficiency of this algorithm for solving the one‐dimensional sine‐Gordon equation. Copyright © 2014 John Wiley & Sons, Ltd.     

A local radial basis function collocation method to solve the variable‐order time fractional diffusion equation in a two‐dimensional irregular domain

The local radial basis function (RBF) method is a promising solver for variable‐order time fractional diffusion equation (TFDE), as it overcomes the computational burden of the traditional global method. Application of the local RBF method is limited to Fickian diffusion, while real‐world diffusion is usually non‐Fickian in multiple dimensions. This article is the first to extend the application of the local RBF method to two‐dimensional, variable‐order, time fractional diffusion equation in complex shaped domains. One of the main advantages of the local RBF method is that only the nodes located in the subdomain, surrounding the local point, need to be considered when calculating the numerical solution at this point. This approach can perform well with large scale problems and can also mitigate otherwise ill‐conditioned problems. The proposed numerical approach is checked against two examples with curved boundaries and known analytical solutions. Shape parameter and subdomain node number are investigated for their influence on the accuracy of the local RBF solution. Furthermore, quantitative analysis, based on root‐mean‐square error, maximum absolute error, and maximum error of the partial derivative indicates that the local RBF method is accurate and effective in approximating the variable‐order TFDE in two‐dimensional irregular domains.     

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.