A wavelet operational method for solving fractional partial differential equations numerically

Fractional calculus is an extension of derivatives and integrals to non-integer orders, and a partial differential equation involving the fractional calculus operators is called the fractional PDE. They have many applications in science and engineering. However not only the analytical solution existed for a limited number of cases, but also the numerical methods are very complicated and difficult. In this paper, we newly establish the simulation method based on the operational matrices of the orthogonal functions. We formulate the operational matrix of integration in a unified framework. By using the operational matrix of integration, we propose a new numerical method for linear fractional partial differential equation solving. In the method, we (1) use the Haar wavelet; (2) establish a Lyapunov-type matrix equation; and (3) obtain the algebraic equations suitable for computer programming. Two examples are given to demonstrate the simplicity, clarity and powerfulness of the new method.     

Exact solutions for nonlinear partial differential equations via Exp‐function method

In this article, the Exp‐function method is applied to nonlinear Burgers equation and special fifth‐order partial differential equation. Using this method, we obtain exact solutions for these equations. The method is straightforward and concise, and its applications are promising. This method can be used as an alternative to obtain analytical and approximate solutions of different types of nonlinear differential equations. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

A new approach for solving highly nonlinear partial differential equations by successive differentiation method

In this work successive differentiation method is applied to solve highly nonlinear partial differential equations (PDEs) such as Benjamin–Bona–Mahony equation, Burger's equation, Fornberg–Whitham equation, and Gardner equation. To show the efficacy of this new technique, figures have been incorporated to compare exact solution and results of this method. Wave variable is used to convert the highly nonlinear PDE into ordinary differential equation with order reduction. Then successive differentiation method is utilized to obtain the numerical solution of considered PDEs in this paper. Copyright © 2017 John Wiley & Sons, Ltd.     

The alternative Legendre Tau method for solving nonlinear multi-order fractional differential equations

In this paper, the alternative Legendre polynomials (ALPs) are used to approximate the solution of a class of nonlinear multi-order fractional differential equations (FDEs). First, the operational matrix of fractional integration of an arbitrary order and the product operational matrix are derived for ALPs. These matrices together with the spectral Tau method are then utilized to reduce the solution of the mentioned equations into the one of solving a system of nonlinear algebraic equations with unknown ALP coefficients of the exact solution. The fractional derivatives are considered in the Caputo sense and the fractional integration is described in the Riemann-Liouville sense. Numerical examples illustrate that the present method is very effective for linear and nonlinear multi-order FDEs and high accuracy solutions can be obtained only using a small number of ALPs.     

The convergence and stability analysis of the Jacobi collocation method for solving nonlinear fractional differential equations with integral boundary conditions

In this paper, we apply the Jacobi collocation method for solving nonlinear fractional differential equations with integral boundary conditions. Due to existence of integral boundary conditions, after reformulation of this equation in the integral form, the method is proposed for solving the obtained integral equation. Also, the convergence and stability analysis of the proposed method are studied in two main theorems. Furthermore, the optimum degree of convergence in the L₂ norm is obtained for this method. Furthermore, some numerical examples are presented in order to illustrate the performance of the presented method. Finally, an application of the model in control theory is introduced. Copyright © 2015 John Wiley & Sons, Ltd.     

Generalized fractional‐order Jacobi functions for solving a nonlinear systems of fractional partial differential equations numerically

In this paper, a collocation spectral numerical algorithm is presented for solving nonlinear systems of fractional partial differential equations subject to different types of conditions. A proposed error analysis investigates the convergence of the mentioned algorithm. Some numerical examples confirm the efficiency and accuracy of the method.     

A Haar wavelet collocation approach for solving one and two-dimensional second-order linear and nonlinear hyperbolic telegraph equations

We have developed a new numerical method based on Haar wavelet (HW) in this article for the numerical solution (NS) of one- and two-dimensional hyperbolic Telegraph equations (HTEs). The proposed technique is utilized for one- and two-dimensional linear and nonlinear problems, which shows its advantage over other existing numerical methods. In this technique, we approximated both space and temporal derivatives by the truncated Haar series. The algorithm of the method is simple and we can implement easily in any other programming language. The technique is tested on some linear and nonlinear examples from literature. The maximum absolute errors (MAEs), root mean square errors (RMSEs), and computational convergence rate are calculated for different number of collocation points (CPs) and also some 3D graphs are also drawn. The results show that the proposed technique is simply applicable and accurate.     

An exponential regula falsi method for solving nonlinear equations

A class of modified regula falsi iterative formulae for solving nonlinear equations is presented in this paper. This method is shown to be quadratically convergent for the sequence of diameters and the sequence of iterative points. The numerical experiments show that new method is effective and comparable to well-known methods.     

Efficient analytic method for solving nonlinear fractional differential equations

In this paper, we propose a new approach of the generalized differential transform method (GDTM) for solving nonlinear fractional differential equations. In GDTM, it is a key to derive a recurrence relation of generalized differential transform (GDT) associated with the solution in the given fractional equation. However, the recurrence relations of complex nonlinear functions such as exponential, logarithmic and trigonometry functions have not been derived before in GDTM. We propose new algorithms to construct the recurrence relations of complex nonlinear functions and apply the GDTM with the proposed algorithms to solve nonlinear fractional differential equations. Several illustrative examples are demonstrated to show the effectiveness of the proposed method. It is shown that the proposed technique is robust and accurate for solving fractional differential equations.     

Implicit Euler method for numerical solution of nonlinear stochastic partial differential equations with multiplicative trace class noise

In this paper, we consider the numerical approximation of stochastic partial differential equations with nonlinear multiplicative trace class noise. Discretization is obtained by spectral collocation method in space, and semi‐implicit Euler method is used for the temporal approximation. Our purpose is to investigate the convergence of the proposed method. The rate of convergence is obtained, and some numerical examples are included to illustrate the estimated convergence rate.     

A BFGS algorithm for solving symmetric nonlinear equations

In this article, we propose a BFGS method for solving symmetric nonlinear equations. The presented method possesses some favourable properties: (a) the generated sequence of iterates is norm descent; (b) the generated sequence of the quasi-Newton matrix is positive definite and (c) this method possesses the global convergence and superlinear convergence. Numerical results show that the presented method is interesting.     

Oscillation for systems of nonlinear neutral type parabolic partial functional differential equations

This paper discusses the oscillation of solutions for systems of nonlinear neutral type parabolic partial fuctional differential equations of the form     

Multiquadric quasi‐interpolation methods for solving partial differential algebraic equations

In this article, we propose two meshless collocation approaches for solving time dependent partial differential algebraic equations (PDAEs) in terms of the multiquadric quasi‐interpolation schemes. In presenting the process of the solution, the error is estimated. Furthermore, the comparisons on condition numbers of the collocation matrices using different methods and the sensitivity of the shape parameter c are given. With the use of the appropriate collocation points, the method for PDAEs with index‐2 is improved. The results show that the methods have some advantages over some known methods, such as the smaller condition numbers or more accurate solutions for PDAEs which has an modal index‐2 or an impulse solution with index‐2. Copyright © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 95–119, 2014     

Comparison between variational iteration method and homotopy perturbation method for linear and nonlinear partial differential equations with the nonhomogeneous initial conditions

In this study, linear and nonlinear partial differential equations with the nonhomogeneous initial conditions are considered. We used Variational iteration method (VIM) and Homotopy perturbation method (HPM) for solving these equations. Both methods are used to obtain analytic solutions for different types of differential equations. Four examples are presented to show the application of the present techniques. In these schemes, the solution takes the form of a convergent series with easily computable components. The present methods perform extremely well in terms of efficiency and simplicity. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

A new technique of using homotopy analysis method for solving high‐order nonlinear differential equations

In this paper, a new technique of homotopy analysis method (HAM) is proposed for solving high‐order nonlinear initial value problems. This method improves the convergence of the series solution, eliminates the unneeded terms and reduces time consuming in the standard homotopy analysis method (HAM) by transform the nth‐order nonlinear differential equation to a system of n first‐order equations. Second‐ and third‐ order problems are solved as illustration examples of the proposed method. Copyright © 2010 John Wiley & Sons, Ltd.     

A new smoothing Newton method for solving constrained nonlinear equations

In this paper, a new smoothing Newton method is proposed for solving constrained nonlinear equations. We first transform the constrained nonlinear equations to a system of semismooth equations by using the so-called absolute value function of the slack variables, and then present a new smoothing Newton method for solving the semismooth equations by constructing a new smoothing approximation function. This new method is globally and quadratically convergent. It needs to solve only one system of unconstrained equations and to perform one line search at each iteration. Numerical results show that the new algorithm works quite well.     

An integrated numerical method with error analysis for solving fractional differential equations of quintic nonlinear type arising in applied sciences

In this study, fractional differential equations having quintic nonlinearity are considered by proposing an accurate numerical method based on the matching polynomial and matrix‐collocation system. This method provides an integration between matrix and fractional derivative, which makes it fast and efficient. A hybrid computer program is designed by making use of the fast algorithmic structure of the method. An error analysis technique consisting of the fractional‐based residual function is constructed to scrutinize the precision of the method. Some error tests are also performed. Figures and tables present the consistency of the approximate solutions of highly stiff model problems. All results point out that the method is effective, simple, and eligible.     

A numerical scheme based on Bernoulli wavelets and collocation method for solving fractional partial differential equations with Dirichlet boundary conditions

In this paper, an efficient and accurate numerical method is presented for solving two types of fractional partial differential equations. The fractional derivative is described in the Caputo sense. Our approach is based on Bernoulli wavelets collocation techniques together with the fractional integral operator, described in the Riemann‐Liouville sense. The main characteristic behind this approach is to reduce such problems to those of solving systems of algebraic equations, which greatly simplifies the problem. By using Newton's iterative method, this system is solved and the solution of fractional partial differential equations is achieved. Some results concerning the error analysis are obtained. The validity and applicability of the method are demonstrated by solving four numerical examples. Numerical examples are presented in the form of tables and graphs to make comparisons with the results obtained by other methods and with the exact solutions much easier.     

Generalized wavelet quasilinearization method for solving population growth model of fractional order

The primary aim of this study is to introduce and develop a generalized wavelet method together with the quasilinearization technique to solve the Volterra's population growth model of fractional order. Unlike the existing operational matrix methods based on orthogonal functions, we formulate the wavelet operational matrices of general order integration without using the block pulse functions. Consequently, the governing problem is transformed into an equivalent system of algebraic equations, which can be tackled with any classical method. The applicability of the proposed method is demonstrated via an illustrative comparison of the numerical outcomes with those found by other known methods. The experimental outcomes demonstrate that the proposed method is fast, accurate, simple, and computationally reliable.