Sine‐cosine wavelet method for fractional oscillator equations

Purpose In this article, a novel computational method is introduced for solving the fractional nonlinear oscillator differential equations on the semi‐infinite domain. The purpose of the proposed method is to get better and more accurate results. Design/methodology/approach The proposed method is the combination of the sine‐cosine wavelets and Picard technique. The operational matrices of fractional‐order integration for sine‐cosine wavelets are derived and constructed. Picard technique is used to convert the fractional nonlinear oscillator equations into a sequence of discrete fractional linear differential equations. Operational matrices of sine‐cosine wavelets are utilized to transformed the obtained sequence of discrete equations into the systems of algebraic equations and the solutions of algebraic systems lead to the solution of fractional nonlinear oscillator equations. Findings The convergence and supporting analysis of the method are investigated. The operational matrices contains many zero entries, which lead to the high efficiency of the method, and reasonable accuracy is achieved even with less number of collocation points. Our results are in good agreement with exact solutions and more accurate as compared with homotopy perturbation method, variational iteration method, and Adomian decomposition method. Originality/value Many engineers can utilize the presented method for solving their nonlinear fractional models.     

Application of the Caputo‐Fabrizio and Atangana‐Baleanu fractional derivatives to mathematical model of cancer chemotherapy effect

In this paper, we obtain approximate‐analytical solutions of a cancer chemotherapy effect model involving fractional derivatives with exponential kernel and with general Mittag‐Leffler function. Laplace homotopy perturbation method and the modified homotopy analysis transform method were applied. The first method is based on a combination of the Laplace transform and homotopy methods, while the second method is an analytical technique based on homotopy polynomial. The cancer chemotherapy effect equations are solved numerically and analytically using the aforesaid methods. Illustrative examples are included to demonstrate the validity and applicability of the presented technique with new fractional‐order derivatives with exponential decay law and with general Mittag‐Leffler law.     

Numerical analysis nonlinear multi‐term time fractional differential equation with collocation method via fractional B‐spline

In this work, we present numerical analysis for nonlinear multi‐term time fractional differential equation which involve Caputo‐type fractional derivatives for . The proposed method is based on utilization of fractional B‐spline basics in collocation method. The scheme can be readily obtained efficient and quite accurate with less computational work numerical result. The proposal approach transform nonlinear multi‐term time fractional differential equation into a suitable linear system of algebraic equations which can be solved by a suitable numerical method. The numerical experiments will be verify to demonstrate the effectiveness of our method for solving one‐ and two‐dimensional multi‐term time fractional differential equation.     

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.     

A new analysis for fractional model of regularized long‐wave equation arising in ion acoustic plasma waves

The key purpose of the present work is to constitute a numerical scheme based on q‐homotopy analysis transform method to examine the fractional model of regularized long‐wave equation. The regularized long‐wave equation explains the shallow water waves and ion acoustic waves in plasma. The proposed technique is a mixture of q‐homotopy analysis method, Laplace transform, and homotopy polynomials. The convergence analysis of the suggested scheme is verified. The scheme provides and n‐curves, which show that the range convergence of series solution is not a local point effects and elucidate that it is superior to homotopy analysis method and other analytical approaches. Copyright © 2017 John Wiley & Sons, Ltd.     

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.     

Solving 2D time‐fractional diffusion equations by a pseudospectral method and Mittag‐Leffler function evaluation

Two‐dimensional time‐fractional diffusion equations with given initial condition and homogeneous Dirichlet boundary conditions in a bounded domain are considered. A semidiscrete approximation scheme based on the pseudospectral method to the time‐fractional diffusion equation leads to a system of ordinary fractional differential equations. To preserve the high accuracy of the spectral approximation, an approach based on the evaluation of the Mittag‐Leffler function on matrix arguments is used for the integration along the time variable. Some examples along with numerical experiments illustrate the effectiveness of the proposed approach. Copyright © 2016 John Wiley & Sons, Ltd.     

Synchronization of nonlinear master–slave systems under input delay and slope‐restricted input nonlinearity

This article addresses the synchronization of nonlinear master–slave systems under input time‐delay and slope‐restricted input nonlinearity. The input nonlinearity is transformed into linear time‐varying parameters belonging to a known range. Using the linear parameter varying (LPV) approach, applying the information of delay range, using the triple‐integral‐based Lyapunov–Krasovskii functional and utilizing the bounds on nonlinear dynamics of the nonlinear systems, nonlinear matrix inequalities for designing a simple delay‐range‐dependent state feedback control for synchronization of the drive and response systems is derived. The proposed controller synthesis condition is transformed into an equivalent but relatively simple criterion that can be solved through a recursive linear matrix inequality based approach by application of cone complementary linearization algorithm. In contrast to the conventional adaptive approaches, the proposed approach is simple in design and implementation and is capable to synchronize nonlinear oscillators under input delays in addition to the slope‐restricted nonlinearity. Further, time‐delays are treated using an advanced delay‐range‐dependent approach, which is adequate to synchronize nonlinear systems with either higher or lower delays. Furthermore, the resultant approach is applicable to the input nonlinearity, without using any adaptation law, owing to the utilization of LPV approach. A numerical example is worked out, demonstrating effectiveness of the proposed methodology in synchronization of two chaotic gyro systems. © 2015 Wiley Periodicals, Inc. Complexity 21: 220–233, 2016     

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     

Analytical approximate solutions of Riesz fractional diffusion equation and Riesz fractional advection–dispersion equation involving nonlocal space fractional derivatives

In this paper, we consider the analytical solutions of fractional partial differential equations (PDEs) with Riesz space fractional derivatives on a finite domain. Here we considered two types of fractional PDEs with Riesz space fractional derivatives such as Riesz fractional diffusion equation (RFDE) and Riesz fractional advection–dispersion equation (RFADE). The RFDE is obtained from the standard diffusion equation by replacing the second‐order space derivative with the Riesz fractional derivative of order α∈(1,2]. The RFADE is obtained from the standard advection–dispersion equation by replacing the first‐order and second‐order space derivatives with the Riesz fractional derivatives of order β∈(0,1] and of order α∈(1,2] respectively. Here the analytic solutions of both the RFDE and RFADE are derived by using modified homotopy analysis method with Fourier transform. Then, we analyze the results by numerical simulations, which demonstrate the simplicity and effectiveness of the present method. Here the space fractional derivatives are defined as Riesz fractional derivatives. Copyright © 2015 John Wiley & Sons, Ltd.     

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.     

Numerical solutions of time fractional degenerate parabolic equations by variational iteration method with Jumarie‐modified Riemann–Liouville derivative

In this article, the fractional variational iteration method is employed for computing the approximate analytical solutions of degenerate parabolic equations with fractional time derivative. The time‐fractional derivatives are described by the use of a new approach, the so‐called Jumarie modified Riemann–Liouville derivative, instead in the sense of Caputo. The approximate solutions of our model problem are calculated in the form of convergent series with easily computable components. Moreover, the numerical solution is compared with the exact solution and the quantitative estimate of accuracy is obtained. The results of the study reveal that the proposed method with modified fractional Riemann–Liouville derivatives is efficient, accurate, and convenient for solving the fractional partial differential equations in multi‐dimensional spaces without using any linearization, perturbation or restrictive assumptions. Copyright © 2011 John Wiley & Sons, Ltd.     

Multi‐switching combination–combination synchronization of non‐identical fractional‐order chaotic systems

In this paper, multi‐switching combination–combination synchronization scheme has been investigated between a class of four non‐identical fractional‐order chaotic systems. The fractional‐order Lorenz and Chen's systems are taken as drive systems. The combination–combination of multi drive systems is then synchronized with the combination of fractional‐order Lü and Rössler chaotic systems. In multi‐switching combination–combination synchronization, the state variables of two drive systems synchronize with different state variables of two response systems simultaneously. Based on the stability of fractional‐order chaotic systems, the multi‐switching combination–combination synchronization of four fractional‐order non‐identical systems has been investigated. For the synchronization of four non‐identical fractional‐order chaotic systems, suitable controllers have been designed. Theoretical analysis and numerical results are presented to demonstrate the validity and feasibility of the applied method. Copyright © 2017 John Wiley & Sons, Ltd.     

Variational iteration method for solving the space‐ and time‐fractional KdV equation

This paper presents numerical solutions for the space‐ and time‐fractional Korteweg–de Vries equation (KdV for short) using the variational iteration method. The space‐ and time‐fractional derivatives are described in the Caputo sense. In this method, general Lagrange multipliers are introduced to construct correction functionals for the problems. The multipliers in the functionals can be identified optimally via variational theory. The iteration method, which produces the solutions in terms of convergent series with easily computable components, requiring no linearization or small perturbation. The numerical results show that the approach is easy to implement and accurate when applied to space‐ and time‐fractional KdV equations. The method introduces a promising tool for solving many space–time fractional partial differential equations. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2007     

Synchronization of memristor‐based delayed BAM neural networks with fractional‐order derivatives

This article deals with the problem of synchronization of fractional‐order memristor‐based BAM neural networks (FMBNNs) with time‐delay. We investigate the sufficient conditions for adaptive synchronization of FMBNNs with fractional‐order 0 < α < 1. The analysis is based on suitable Lyapunov functional, differential inclusions theory, and master‐slave synchronization setup. We extend the analysis to provide some useful criteria to ensure the finite‐time synchronization of FMBNNs with fractional‐order 1 < α < 2, using Mittag‐Leffler functions, Laplace transform, and linear feedback control techniques. Numerical simulations with two numerical examples are given to validate our theoretical results. Presence of time‐delay and fractional‐order in the model shows interesting dynamics. © 2016 Wiley Periodicals, Inc. Complexity 21: 412–426, 2016     

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.     

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.     

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     

A three‐dimensional finite element model for the control of certain non‐linear bioreactors

This paper deals with the dynamics of non‐linear distributed parameter fixed‐bed bioreactors. The model consists of a pair of non‐linear partial differential (evolution) equations. The true spatially three‐dimensional situation is considered instead of the usual one‐dimensional approximation. This enables one to take into account the effects of flow profiles and the true location of the measurement transducer. The (output) evolution of the corresponding open‐loop control system is simulated. Furthermore, the associated closed‐loop system with respect to the relevant output function is considered. Especially, the asymptotic output tracking is found to be successful by applying the usual process based on the state feedback linearization. Copyright © 2000 John Wiley & Sons, Ltd.