An approximate solution for a fractional model of generalized Harry Dym equation

The nonlinear partial differential equation of Harry Dym is generalized by replacing the first-order time derivative by a fractional derivative of order $$\alpha ,\; 0 \le \alpha \le 2$$α,0≤α≤2. The aim of the present paper is to obtain an approximate solution of time fractional generalized Harry Dym equation using Adomian Decomposition Method (ADM). The fractional derivative is described in the Caputo sense. Numerical examples are given to show the application of the present technique. The results show that the solution of ADM is in good agreement with the exact solution when $$\alpha =1$$α=1, also reveal that the method is very simple and effective.


Introduction
Nonlinear partial differential equations appear in many branches of physics, engineering and applied mathematics. In recent years, there has been a growing interest in the field of fractional calculus. Oldham and Spanier [1], Miller and Ross [2] and Podlubny [3] provide the history and a comprehensive treatment of this subject. Fractional calculus is the field of mathematical analysis, which deals with the investigation and applications of integrals and derivatives of arbitrary order, which can be real or complex. The idea appeared in a letter by Leibniz to L'Hospital in (1695). The subject of fractional calculus have gained importance during the past three decades and popularity, mainly due to its demonstrated applications in different areas of physics and engineering. Several fields of applications of fractional differentiation and fractional integration are already well established, some others have just started. Many applications of fractional calculus can be found in turbulence and fluid dynamics, stochastic dynamical systems, plasma physics and controlled thermonuclear fusion, nonlinear control theory, image processing, nonlinear biological systems; for more details see [4] and the references therein. Indeed, it provides several potentially useful tools for solving differential equations. The most important advantage of using fractional differential equations is their nonlocal property. It is known [5] that the integer order differential operator is nonlocal. This means that the next state of a system depends not only upon its current state but also upon all of its historical states. Hence, the importance of investigating fractional equations arises from the necessity to sharpen the concepts of equilibrium, stability states, and time evolution in the long time limit. There have been few attempts to solve linear problems with multiple fractional derivatives. In [6], an approximate solution based on the decomposition method is given for the generalized fractional diffusionwave equation. Not much work has been done for nonlinear problems, and only a few numerical schemes have been proposed to solve nonlinear fractional differential equations. More recently, applications have included classes of nonlinear equations with multi-order fractional derivative, and this motivates us to develop a numerical scheme for their solutions.
The main objective of the present paper is to use the Adomian decomposition method [7,8] to calculate the approximate solutions of the Dym equation that is named for Harry Dym and occurs in the study of Solitons. The Dym equation first appeared in Kruskal [9] and is attributed to an unpublished paper by Harry Dym. Harry Dym is a completely integrable nonlinear evolution equation that may be solved by means of the inverse scattering transform. It is interesting because it obeys an infinite number of conservation laws and it does not possess the Painleve property. The Dym equation has strong links to the Korteweg-de Vries equation. As mentioned in [10], the Harry Dym has been recently generalized in different ways. The authors in [11] introduced and investigated the coupled Harry Dym equations, and showed that these equations are concerned with the new iso-spectral flows. Different generalizations of the Harry Dym equation have been constructed in [12]. In [13] group analysis of time fractional Harry Dym equation was performed and some group invariant solutions are obtained. An efficient approach based on Homotopy perturbation method and Sumudu transform were proposed in [14] to solve fractional Harry Dym equation for 0 a 1. Beside these generalizations, we would like to present a new generalization of the nonlinear partial differential equation of Harry Dym by replacing the first-order time derivative by a fractional derivative of order a; 0 a 2, and takes the form where a is a parameter describing the order of the fractional derivative, x and t are the space and time variables, and uðx; tÞ is the field defined in the space domain ðÀ1; 1Þ. Theoretically, a can be any positive number. Note that for a ¼ 1, Eq. (1.1) represents the standard Harry Dym equation, which has exact solution (see, [5]) , where a, b and c are suitable constants. This type of equation has been investigated by several authors. In [5], an approximate analytical solution of time fractional Harry Dym equation is obtained using Homotopy perturbation method only for 0 a 1.

Explicit solutions for Harry Dym equation was obtained by
Fuchssteinert et al. [12]. Soliton solutions of the 2 þ 1 dimensional Harry Dym equation was found by Halim in [17]. The plan of this work is the following. A brief review of the fractional calculus theory is given in Sect. 2. In Sect. 3, we use the decomposition method to construct our numerical solutions for Eq. (1.1). The general response expression contains a parameter describing the order of the fractional derivative that can be varied to obtain various responses. In Sect. 4, we present some numerical results (Graphs and Tables) to show the nature of the solution as the fractional derivative parameter is changed.

Basic definition of fractional calculus
This section is devoted to a description of the operational properties of the purpose of acquainting with sufficient fractional calculus theory, to enable us to follow the solution of the generalized Harry Dym equation. Many definitions and studies of fractional calculus have been proposed in the last two centuries. These definitions include, Riemman-Liouville, Weyl, Reize, Campos, Caputa, and Nishimoto fractional operator. Mainly, in this paper, we will re-introduce section 2 of [6]. The Riemann-Liouville definition of fractional derivative operator J a a is defined as follows: Properties of the operator J a can be found in [1], we mention the following: for f 2 L 1 ½a; b; a; b ! 0 and Cðaþcþ1Þ ðx À aÞ aþc .  The Riemann-Liouville derivative has certain disadvantages when trying to model real-world phenomena with fractional differential equations. Therefore, we shall introduce now a modified fractional differentiation operator D a proposed by Caputo in his work on the theory of visco-elasticity [19].
The fractional derivative of f ðxÞ in the Caputo sense is defined as The fractional derivative is considered in the Caputo sense. To solve differential equations, we need to specify additional conditions to produce a unique solution. For the case of Caputo fractional differential equations, these additional conditions are just the traditional conditions, which are taken to those of classical differential equations, and are, therefore, familiar to us. In contrast, for Riemann-Liouville fractional differential equations, these additional conditions constitute certain fractional derivatives of the unknown solution at the initial point x ¼ 0, which are functions of x. The unknown function u ¼ uðx; tÞ is assumed to be a causal function of time, i.e., vanishing for t\0. Also, the initial conditions are not physical; furthermore, it is not clear how much quantities are to be measured from experiment, say, so that they can be appropriately assigned in an analysis. For more details on the geometric and physical interpretation for fractional derivatives of both Riemann-Liouville and Caputo types see [18]. For mathematical properties of fractional derivatives and integrals, one can consult the above-mentioned references.

Analysis of the method
To illustrate the basic idea of the ADM for fractional Harry Dym equation, in an operator form, Eq. (1.1) can be written as where the fractional differential operator D a t is defined as D a t ¼ o a ot a , so that D a t is the operator defined in (2.1). In this study, we shall consider Eq. (1.1) subject to the initial conditions uðx; 0Þ ¼ g 1 ðxÞ; u t ðx; 0Þ ¼ g 2 ðxÞ; x 2 IR ð3:2Þ To solve the nonlinear fractional Eq. (1.1), we apply the operator J a , the inverse of the operator D a t , on both sides of Eq. (3.1) and using the initial condition (3.2) yields where the components u n ðx; tÞ will be determined recursively. The nonlinear function /ðuÞ is then written in the decomposed form /ðuÞ ¼ X 1 n¼0 A n ðu 0 ; u 1 ; . . .; u n Þ; ð3:5Þ where A n are called the Adomian polynomials, these polynomials can be calculated for all forms of nonlinearity according to specific algorithms constructed by Adomian [7,8]. In this specific nonlinearity, we use the general form of formula for A n polynomials as This formula is easy to set computer code to get as many polynomials as we need in the calculation of the numerical solution as well as explicit solutions. To be easy to follow by the reader, we can give the first few Adomian polynomials for /ðuÞ ¼ u 3 u xxx of the nonlinearity as and so on, the rest of the polynomials can be constructed in a similar manner. Wazwaz, in [21] developed an alternative approach for the construction of these polynomials. Substituting (3.4) and (3.5) into both sides of Eq. (3.3) gives Following the decomposition method, we introduce the recursive relations as give an approximate solution to the nonlinear fractional Harry Dym equation. Moreover, the decomposition series solutions (3.11) generally converge very rapidly in real physical problems. The convergence of the decomposition series have been investigated by several authors [15,16].

Numerical applications and results
In this section, we present some numerical results to demonstrate the behavior of the solution as the order of the time fractional derivative is changed. We shall illustrate the numerical scheme by two cases. The first case is somewhat artificial in the sense that the exact answer, for the special case a ¼ 1, is known in advance, and the initial condition is directly taken from this answer. Nonetheless, such an approach is needed to evaluate the accuracy of the numerical scheme, and to examine the effect of varying the order of the time fractional derivative on the behavior of the solution. To calculate the terms of the decomposition series (3.11) for uðx; tÞ, we shall mention that a second initial condition u t ðx; 0Þ ¼ gðxÞ, for 1\a 2 is assumed to ensure continuous dependence of the solution on the parameter a in the transition from a ¼ 1 À to a ¼ 1 þ . Therefore, in (3.9), we need to distinguish two cases: 1. Case I: For 0\a 1, in this case, we choose m ¼ 1 in (3.9), and so upon using Mathematica, the solution can be obtained as It can be seen from Fig. 1 that the solution obtained by the present method is nearly identical to the exact. Also, Fig. 1 shows the exact solution uðx; tÞ of the regular Harry Dym equation (a ¼ 1) given by Mokhtari [20] for constants value of a ¼ 4; b ¼ 1; c ¼ 1.
It is known that using ADM, the series solution converges very rapidly. The rapid convergence means only few terms are required to get analytical function. From Fig. 1, it is to be noted that only four terms of the decomposition series were used for our approximations. It is evident that the solution can be improved by adding more terms of the decomposition series. 2. Case II: For 1\a 2, in this case, we choose m ¼ 2 in (3.9) and the components are determined by   Figure 4, shows the obtained approximate solution to this case when a ¼ 2. Figure 5, shows profile solutions when t ¼ 1 for different values of a, i.e. a ¼ 1:1; 1:2; 1:4; 1:9; 2.

Conclusions
In this paper, the application of Adomian decomposition method was extended to explicit and numerical solution of time fractional nonlinear PDEs in mathematical physics with initial conditions. It may be concluded that the decomposition methodology is a very powerful and efficient technique that provides more realistic solutions. It also provides series solutions which generally converge very rapidly in real physical problems. The decomposition method, which is proved to be efficient in solving differential equations, can be applied to various types of fractional differential equations [6]. Finally, the recent appearance of nonlinear fractional partial differential equations as adequate models in science and engineering makes it necessary to investigate the method of solutions for such equations.
Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.