Solving a system of nonlinear fractional partial differential equations using homotopy analysis method

In this article, the homotopy analysis method (HAM) has been employed to obtain solutions of a System of nonlinear fractional partial differential equations. This indicates the validity and great potential of the homotopy analysis method for solving system of fractional partial differential equations. The fractional derivative is described in the Caputo sense.     

Comparison between homotopy analysis method and optimal homotopy asymptotic method for nth‐order integro‐differential equation

This paper presents general framework for solving the nth‐order integro‐differential equation using homotopy analysis method (HAM) and optimal homotopy asymptotic method (OHAM). OHAM is parameter free and can provide better accuracy over the HAM at the same order of approximation. Furthermore, in OHAM the convergence region can be easily adjusted and controlled. Comparison, via two examples, between our solution using HAM and OHAM and the exact solution shows that the HAM and the OHAM are effective and accurate in solving the nth‐order integro‐differential equation. Copyright © 2011 John Wiley & Sons, Ltd.     

The homotopy perturbation method applied to the nonlinear fractional Kolmogorov–Petrovskii–Piskunov equations

The fractional derivatives in the sense of Caputo, and the homotopy perturbation method are used to construct approximate solutions for nonlinear Kolmogorov–Petrovskii–Piskunov (KPP) equations with respect to time and space fractional derivatives. Also, we apply complex transformation to convert a time and space fractional nonlinear KPP equation to an ordinary differential equation and use the homotopy perturbation method to calculate the approximate solution. This method is efficient and powerful in solving wide classes of nonlinear evolution fractional order equations.     

Analysis of nonlinear fractional partial differential equations with the homotopy analysis method

In this paper, the time fractional partial differential equations are investigated by means of the homotopy analysis method. This technique is extended to study the partial differential equations of fractal order for the first time. The accurate series solutions are obtained. This indicates the validity and great potential of the homotopy analysis method for solving nonlinear fractional partial differential equations.     

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.     

Approximations of the nonlinear Volterra's population model by an efficient numerical method

In this paper, we apply the new homotopy perturbation method to solve the Volterra's model for population growth of a species in a closed system. This technique is extended to give solution for nonlinear integro‐differential equation in which the integral term represents the total metabolism accumulated fromtime zero. The approximate analytical procedure only depends on two components. The newhomotopy perturbationmethodwas applied to nonlinear integro‐differential equations directly and by converting the problem into nonlinear ordinary differential equation. We also compare this method with some other numerical results and show that the present approach is less computational and is applicable for solving nonlinear integro‐differential equations and ordinary differential equations as well. Copyright © 2011 John Wiley & Sons, Ltd.     

On the homotopy analysis method for backward/forward-backward stochastic differential equations

In this paper, an analytic approximation method for highly nonlinear equations, namely the homotopy analysis method (HAM), is employed to solve some backward stochastic differential equations (BSDEs) and forward-backward stochastic differential equations (FBSDEs), including one with high dimensionality (up to 12 dimensions). By means of the HAM, convergent series solutions can be quickly obtained with high accuracy for a FBSDE in a 6-dimensional case, within less than 1 % CPU time used by a currently reported numerical method for the same case [34]. Especially, as dimensionality enlarges, the increase of computational complexity for the HAM is not as dramatic as this numerical method. All of these demonstrate the validity and high efficiency of the HAM for the backward/forward-backward stochastic differential equations in science, engineering, and finance.     

A new method for homoclinic solutions of ordinary differential equations

Consideration is given to the homoclinic solutions of ordinary differential equations. We first review the Melnikov analysis to obtain Melnikov function, when the perturbation parameter is zero and when the differential equation has a hyperbolic equilibrium. Since Melnikov analysis fails, using Homotopy Analysis Method (HAM, see [Liao SJ. Beyond perturbation: introduction to the homotopy analysis method. Boca Raton: Chapman & Hall/CRC Press; 2003; Liao SJ. An explicit, totally analytic approximation of Blasius’ viscous flow problems. Int J Non-Linear Mech 1999;34(4):759–78; Liao SJ. On the homotopy analysis method for nonlinear problems. Appl Math Comput 2004;147(2):499–513] and others [Abbasbandy S. The application of the homotopy analysis method to nonlinear equations arising in heat transfer. Phys Lett A 2006;360:109–13; Hayat T, Sajid M. On analytic solution for thin film flow of a forth grade fluid down a vertical cylinder. Phys Lett A, in press; Sajid M, Hayat T, Asghar S. Comparison between the HAM and HPM solutions of thin film flows of non-Newtonian fluids on a moving belt. Nonlinear Dyn, in press]), we obtain homoclinic solution for a differential equation with zero perturbation parameter and with hyperbolic equilibrium. Then we show that the Melnikov type function can be obtained as a special case of this homotopy analysis method. Finally, homoclinic solutions are obtained (for nontrivial examples) analytically by HAM, and are presented through graphs.     

Approximate analytical solutions of Schnakenberg systems by homotopy analysis method

In this paper, the homotopy analysis method (HAM) has been employed to obtain analytical solution of a two reaction–diffusion systems of fractional order (fractional Schnakenberg systems) which has been modeling morphogen systems in developmental biology. Different from all other analytic methods, HAM provides us with a simple way to adjust and control the convergence region of solution series by choosing proper values for auxiliary parameter h. The fractional derivative is described in the Caputo sense. The reason of using fractional order differential equations (FOD) is that FOD are naturally related to systems with memory which exists in most biological systems. Also they are closely related to fractals which are abundant in biological systems. The results derived of the fractional system are of a more general nature. Respectively, solutions of FOD spread at a faster rate than the classical differential equations, and may exhibit asymmetry. However, the fundamental solutions of these equations still exhibit useful scaling properties that make them attractive for applications.     

Explicit series solutions of some linear and nonlinear Schrodinger equations via the homotopy analysis method

In this paper, by means of the homotopy analysis method (HAM), the solutions of some Schrodinger equations are exactly obtained in the form of convergent Taylor series. The HAM contains the auxiliary parameter ?, that provides a convenient way of controlling the convergent region of series solutions. This analytical method is employed to solve linear and nonlinear examples to obtain the exact solutions. HAM is a powerful and easy-to-use analytic tool for nonlinear problems.     

Semi-exact solution for thermo-mechanical analysis of functionally graded elastic-strain hardening rotating disks

In this paper, distributions of stress and strain components of rotating disks with non-uniform thickness and material properties subjected to thermo-elasto-plastic loading are obtained by semi-exact method of Liao’s homotopy analysis method (HAM) and finite element method (FEM). The materials are assumed to be elastic-linear strain hardening and isotropic. The analysis of rotating disk is based on Von Mises’ yield criterion. A two dimensional plane stress analysis is used. The distribution of temperature is assumed to have power forms with the hotter point located at the outer surface of the disk. A mathematical technique of transformation has been proposed to solve the homotopy equations which are originally hard to be handled. The domain of the solution has been substituted by a new domain through which the unknown variable has been taken out from the argument of the function. This makes the solution much easier. A numerical solution of the governing differential equations is also presented based on the Runge–Kutta’s method. The results of three methods are presented and compared which shows good agreements. This verifies the implementation of the HAM and demonstrates its applicability to provide accurate solution for a very complicated case of strongly high nonlinear differential equations with no exact solution. It is important to notice that compared with other methods, HAM needs significant more computation time and computer hardware requirements which limit its application for those problems that other methods can easily handle them.     

Approximate solutions for the Burger and regularized long wave equations by means of the homotopy analysis method

In this work, the homotopy analysis method (HAM), one of the most effective method, is implemented for finding approximate solutions of the Burger and regularized long wave (RLW) equations. Comparisons are made between the results of the proposed method and homotopy perturbation method (HPM). It illustrates the validity and the great potential of the homotopy analysis method in solving nonlinear partial differential equations.     

Solution of the laminar viscous flow in a semi-porous channel in the presence of a uniform magnetic field by using the homotopy analysis method

In this paper, the problem of laminar viscous flow in a semi-porous channel in the presence of a transverse magnetic field is presented and the homotopy analysis method (HAM) is employed to compute an approximation to the solution of the system of nonlinear differential equations governing the problem. It has been attempted to show the capabilities and wide-range applications of the homotopy analysis method in comparison with the numerical method in solving this problem. The obtained solutions, in comparison with the numeric solutions admit a remarkable accuracy. A clear conclusion can be drawn from the numerical method’s (NM) results that the HAM provides highly accurate solutions for nonlinear differential equations.     

Solving nonlinear fractional partial differential equations using the homotopy analysis method

In this article, the homotopy analysis method is applied to solve nonlinear fractional partial differential equations. On the basis of the homotopy analysis method, a scheme is developed to obtain the approximate solution of the fractional KdV, K(2,2), Burgers, BBM‐Burgers, cubic Boussinesq, coupled KdV, and Boussinesq‐like B(m,n) equations with initial conditions, which are introduced by replacing some integer‐order time derivatives by fractional derivatives. The homotopy analysis method for partial differential equations of integer‐order is directly extended to derive explicit and numerical solutions of the fractional partial differential equations. The solutions of the studied models are calculated in the form of convergent series with easily computable components. The results of applying this procedure to the studied cases show the high accuracy and efficiency of the new technique. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Numerical solution of nonlinear partial quadratic integro‐differential equations of fractional order via hybrid of block‐pulse and parabolic functions

In this paper, an effective numerical approach based on a new two‐dimensional hybrid of parabolic and block‐pulse functions (2D‐PBPFs) is presented for solving nonlinear partial quadratic integro‐differential equations of fractional order. Our approach is based on 2D‐PBPFs operational matrix method 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 nonlinear partial quadratic integro‐differential equations is achieved. Convergence analysis and an error estimate associated with the proposed method is obtained, and it is proved that the numerical convergence order of the suggested numerical method is O(h³) . The validity and applicability of the method are demonstrated by solving three numerical examples. Numerical examples are presented in the form of tables and graphs to make comparisons with the exact solutions much easier.     

Two‐dimensional shifted Legendre polynomial collocation method for electromagnetic waves in dielectric media via almost operational matrices

In this paper, a numerical solution of fractional partial differential equations (FPDEs) for electromagnetic waves in dielectric media will be discussed. For the solution of FPDEs, we developed a numerical collocation method using an algorithm based on two‐dimensional shifted Legendre polynomials approximation, which is proposed for electromagnetic waves in dielectric media. By implementing the partial Riemann–Liouville fractional derivative operators, two‐dimensional shifted Legendre polynomials approximation and its operational matrix along with collocation method are used to convert FPDEs first into weakly singular fractional partial integro‐differential equations and then converted weakly singular fractional partial integro‐differential equations into system of algebraic equation. Some results concerning the convergence analysis and error analysis are obtained. Illustrative examples are included to demonstrate the validity and applicability of the technique. Copyright © 2017 John Wiley & Sons, Ltd.     

Three-dimensional flow of a Jeffery fluid over a linearly stretching sheet

This investigation reports the three-dimensional flow of Jeffrey fluid over a linearly stretching surface. Transformation method has been utilized for the reduction of partial differential equations into the ordinary differential equations. The solutions of the nonlinear systems are presented by a homotopy analysis method (HAM). The reported graphical results are analyzed. A comparative study with the previous results of viscous fluid in the literature is made.     

The scaled boundary FEM for nonlinear problems

The traditional scaled boundary finite-element method (SBFEM) is a rather efficient semi-analytical technique widely applied in engineering, which is however valid mostly for linear differential equations. In this paper, the traditional SBFEM is combined with the homotopy analysis method (HAM), an analytic technique for strongly nonlinear problems: a nonlinear equation is first transformed into a series of linear equations by means of the HAM, and then solved by the traditional SBFEM. In this way, the traditional SBFEM is extended to nonlinear differential equations. A nonlinear heat transfer problem is used as an example to show the validity and computational efficiency of this new SBFEM.     

Approximate limit cycles of coupled nonlinear oscillators with fractional derivatives

In this paper, the homotopy analysis method (HAM) is presented to establish the accurate approximate analytical solutions for multi-degree-of-freedom (MDOF) coupled nonlinear oscillators with fractional derivatives. Approximate limit cycles (LCs) of two systems of the coupled fractional van der Pol (VDP) oscillators and the fractional damped Duffing resonator driven by a fractional VDP oscillator are exampled for illustrating the validity and great potential of the HAM. The presented approach can provide approximate LCs very accurately and efficiently compared with some direct simulation results. This method can keep high accuracy and efficiency for both weakly and strongly nonlinear problems with any given fractional order. Furthermore, it is capable of tracking unstable LCs which cannot be generated by some time-marching numerical algorithm. Based on the obtained results, we analyze effect of different fractional orders, coupling coefficient, and nonlinear coefficient of the coupled equations on amplitudes and frequencies of the LCs.     

Solutions of nonlinear thickness-shear vibrations of an infinite isotropic plate with the homotopy analysis method

As a preliminary attempt for the study on nonlinear vibrations of a finite crystal plate, the thickness-shear mode of an infinite and isotropic plate is investigated. By including nonlinear constitutive relations and strain components, we have established nonlinear equations of thickness-shear vibrations. Through further assuming the mode shape of linear vibrations, we utilized the standard Galerkin approximation to obtain a nonlinear ordinary differential equation depending only on time. We solved this nonlinear equation and obtained its amplitude–frequency relation by the homotopy analysis method (HAM). The accuracy of the present results is shown by comparison between our results and the perturbation method. Numerical results show that the homotopy analysis solutions can be adjusted to improve the accuracy. These equations and results are useful in verifying the available methods and improving our further solution strategy for the coupled nonlinear vibrations of finite piezoelectric plates.