Advances in Mathematical Physics
Volume 2015 (2015), Article ID 746327, 8 pages
http://dx.doi.org/10.1155/2015/746327

Research Article

On the Nonlinear Perturbation Rosenau-Hyman Equation: A Model of Nonlinear Scattering Wave

Abdon Atangana,¹ Dumitru Baleanu,2,3 Maysaa’ Mohamed Al Qurashi,⁴ and Xiao-Jun Yang⁵

¹Institute for Groundwater Studies, Faculty of Natural and Agricultural Sciences, University of the Free State, Bloemfontein 9300, South Africa
²Institute of Space Sciences, P.O. Box MG-23, Măgurele, 077125 Bucharest, Romania
³Department of Mathematics and Computer Sciences, Faculty of Arts and Sciences, Cankaya University, 06530 Ankara, Turkey
⁴Mathematics Department, King Saud University, P.O. Box 22452, Riyadh 11495, Saudi Arabia
⁵Department of Mathematics and Mechanics, China University of Mining and Technology, Xuzhou, Jiangsu 221008, China

Received 18 May 2015; Accepted 4 November 2015

Academic Editor: Fawang Liu

Copyright © 2015 Abdon Atangana et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We investigate a nonlinear wave phenomenon described by the perturbation Rosenau-Hyman equation within the concept of derivative with fractional order. We used the Caputo fractional derivative and we proposed an iteration method in order to find a particular solution of the extended perturbation equation. We proved the stability and the convergence of the suggested method for solving the extended equation without any restriction on and also on the perturbations terms. Using the inner product we proved the uniqueness of the special solution. By choosing randomly the fractional orders and , we presented the numerical solutions.

1. Introduction

Fractional calculus is referring to calculus of integrals and derivatives with any real or complex order and has increased significantly attractiveness throughout the past three decades, owing for the most part to its established relevance in many fields of science and engineering. It does without a doubt make quite a lot of potentially advantageously helpful equipment available for finding the solutions of differential and integral equations and further problems connecting exceptional functions of mathematics physics in addition to their additional rooms and generalizations in one and more variables. No wonder researchers nowadays are interested in modelling several physical problems within the scope of the fractional calculus. This leads us to have a look at the nonlinear scattering waves.

Generalized Korteweg-de Vries equations with nonlinear dispersion can propagate compactly supported solitary waves, referred to as compactons [1–7]. Numerical simulations give you an idea about an original pulse wider than a compacton with small amount of radiation; in addition compactons collide elastically suffering only a phase shift after the collision and generating a small amplitude, zero-mass, compact (see [8–13]). Primarily discovered in the Rosenau-Hyman equation for the modeling of pattern developments in liquid drops, compactons have more than a few functions in physics and science [1], for instance, pattern configuration on liquid plane [14]. The equation is in addition the permanent boundary of the disconnected equations of a nonlinear lattice [1–15] and has been generalized to advanced dimensions [16]. Let us also note that the equation has other kinds of solutions, for instance, the elliptic compactons [1, 17, 18]. At long last, several generalizations of the equations have also been well thought-out in the literature, such as the insertion of time-dependent damping and dispersion [19] or the addition of fifth-order dispersion [20]. In this work we are very much interested in investigating the special solution of the perturbation Rosenau-Hyman equation, which we will generalize as follows:

where are small parameters. The above equation will be called the perturbation fractional Rosenau-Hyman equation. Because of the plentiful profits offered by these fractional derivatives, many researchers have been stressing on proposing new definitions of fractional derivatives. We shall present the brief summary of these derivatives in the following section.

2. Some Fractional Calculus Definitions

It is perhaps important to mention that the concept of noninteger order derivative is an older concept, since it is considered to have stanched after an interrogation rose in the year 1695 by Marquis de L’Hopital. We shall present some of these definitions here.

Definition 1 (see [21–26]). The Riemann-Liouville fractional derivative is as follows: according to Riemann-Liouville the fractional derivative of a function says is given as

Definition 2. The Riemann-Liouville fractional integral is as follows: according to Riemann-Liouville, the fractional integral that is considered as antifractional derivative of a function is given as

Definition 3. Caputo fractional derivative is as follows: according to Caputo, the fractional derivative of a continuous and -time differentiable function is given as

Definition 4. The modified Riemann-Liouville fractional derivative of a function is given as

There are other definitions that are not mentioned here. We shall present the derivation of the special solution in the next section. However, in this paper we shall use the Caputo type.

3. Derivation of the Special Solution

One of the great challenges in the field of partial of differential equation is to derive the solution of especially nonlinear equations, not to mention nonlinear equations described with the fractional order derivative, for example, the perturbation fractional Rosenau-Hyman equation that has stronger nonlinearity. No wonder scholars of this area are always trying to find suitable and easier methods to derive at least approximate solution of these difficult equations. In this paper we will make use of a simple iteration method to propose a special solution to the perturbation fractional Rosenau-Hyman equation; the methods is called the new variational iteration method (NVIM) [27], which uses the idea of Lagrange multiplier. We shall first show the methodology to accommodate readers that are not used to this method. Now consider a general partial differential equation with high order () with respect to time; then, in this new VIM, the first step is to apply the Laplace transform on both sides of equation to obtain

The recursive formula of (6) can now be used to put forward the main recursive method connecting the Lagrange multiplier as

Now considering as the restricted term, the Lagrange multiplier can be obtained as [27]

Now applying the inverse Laplace transform on both sides of (6), we obtain the following iteration:

with the first term to be

Now following the above methodology, we can obtain the Lagrange multiplier of (1) to be

and the iteration formula associate is given by

and then

Theorem 5 (convergence analysis). Let us consider

and consider the initial and boundary condition for (1); then the new variation iteration method leads to a special solution of (1).

Proof. To achieve this we shall think about the following fractional sub-Hilbert space of the Hilbert space [28] that can be defined as the set of those functions:

We harmoniously assume that the differential operators are restricted under the norms. Using the definition of the operator, , we have the following:

Our next task is to evaluate with being the inner product. Then, we have the following:

To evaluate the above expression, we shall consider case by case: we shall start with the following:

With the benefit of the Cauchy-Schwartz inequality, we have the following relation

The most important part in the above inequality is ; let us handle this part; first making use of the continuity properties of the derivative, it is possible for us to find a positive constant such that

Therefore

Now if one makes use of the triangular inequality, we transform the above to

However to proceed with our demonstration, we must assume that are bounded, implying that we can find a positive constant, say , such that . Therefore,

Thus, we have the following:

We shall now consider the following:

Again, using the Cauchy-Schwartz inequality, we have the following relationship:

Using the properties of the inner product the above can further be converted to

However,

using the same argument as before, we obtain

Also

Now replacing these relations into (17), we obtain the following relation:

Take

such that

The next step in this proof will be to evaluate

Nevertheless, following the discussion presented, we shall have

where is defined as

This completes the proof.

Our next concern will consist of investigating on the uniqueness of the special solution. We shall assume that there exist two special solutions, say and , satisfying (2); then, using (33), we have

Since both solutions satisfy (2), thus ; then, . Nonetheless,

Thus if is the exact solution of (2) we have the following relation,

with being a very small positive parameter closer to zero. Then the inequality (37) can be converted to

And this implied, with the benefits of the norm, that

3.1. Application of the Scheme

We shall present the special solutions for some examples using the scheme presented earlier.

Example 1. Here we assume that the perturbation fractional Rosenau-Hyman equation is such that ; then we have the following equation:

Using the presented scheme and also using the fractional integral proposed as antiderivative of the conformable fractional derivative, we obtain the following iteration formula:

Now using the algorithm, we reach the following:

Using the iteration formula the remaining terms can be calculated; however, if we go to infinity, the special solution for this version can be given as

which of course satisfies that, for alpha equal to 1, we obtain the exact solution of the Rosenau-Hyman equation

We shall present the numerical results in Figures 1–4 for different value of alpha and . We have depicted the numerical results for different alpha and in Figures 1, 2, 3, and 4; for example, if = 50 and alpha = 0.5, we have Figure 1.
For example if alpha is 0.55 and is 50 we will observe nonlinear wave as represented in Figure 2.
Figures 1–4 as a result of introduction of fractional order derivative help us to understand better the role of nonlinear dispersion in prototype configuration. The solitary wave solutions of these equations have extraordinary assets according to the variation of the order of fractional derivative together with ; in addition the solitary wave collides elastically and they have compact support. The figures also reveal that, when two compactons collide, the interface location is indicated by the beginning of low-amplitude compacton-anticompacton pairs. These equations seem to have only a finite number of local conservation laws [1]. Nonetheless, the behaviour and the stability of these compactons are very similar to those observed in completely integrable systems [1].

Figure 1: The special solution for alpha = 0.5 and = 50.

Figure 2: The special solution for alpha = 0.55 and = 50.

Figure 3: The special solution for alpha = 0.90 and = 20.

Figure 4: The special solution for alpha = 0.95 and = 20.

4. Conclusion

One of the most difficult tasks in the area of applied mathematics is perhaps to make use of mathematical equations in order to explain adequately the physical phenomenon. Partial differential equations have been intensively used for this purpose in the past decades. However researchers have encountered some limitation by modelling physical phenomenon using the integer derivative order. In the way of trying to extend these limitations, the concept of noninteger order derivative has been introduced. In the contemporary day numerous objective occurrences were clarified with pronounced accomplishment in the light of the perception of noninteger order derivatives. Exclusively, the compensations of fractional calculus and fractional order replicas and their diligences in the field of nonlinear wave motion have beforehand been intensively reexamined throughout the last few epochs with exceptional conclusion. We have therefore investigated within that scope of fractional calculus the nonlinear wave phenomenon that is usually described via the Perturbation Rosenau-Hyman equation. We derived the special solution of the extended equation using the so-called Laplace transform and the Lagrange multiplier. We proved with great success the convergence of this method for solving the extended equation. On the other hand we showed the uniqueness of the special solution using some properties of the inner product within a well-constructed Hilbert space. Some numerical results were depicted from Figures 1–4. The special solution seems to describe the real world problem better than the standard solution.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This research project was supported by a grant from the “Research Center of the Center for Female Scientific and Medical Colleges,” Deanship of Scientific Research, King Saud University. The authors are thankful to Visiting Professor Program at King Saud University for support.

References

P. Rosenau and J. M. Hyman, “Compactons: solitons with finite wavelength,” Physical Review Letters, vol. 70, no. 5, pp. 564–567, 1993.

F. Cooper, H. Shepard, and P. Sodano, “Solitary waves in a class of generalized Korteweg-de Vries equations,” Physical Review E, vol. 48, no. 5, pp. 4027–4032, 1993.

P. Rosenau, “On nonanalytic solitary waves formed by a nonlinear dispersion,” Physics Letters A, vol. 230, no. 5-6, pp. 305–318, 1997.

J. de Frutos, M. A. López-Marcos, and J. M. Sanz-Serna, “A finite difference scheme for the K(2, 2) compacton equation,” Journal of Computational Physics, vol. 120, no. 2, pp. 248–252, 1995.

F. Rus and F. R. Villatoro, “Padé numerical method for the Rosenau-Hyman compacton equation,” Mathematics and Computers in Simulation, vol. 76, no. 1–3, pp. 188–192, 2007.

B. Mihaila, A. Cardenas, F. Cooper, and A. Saxena, “Stability and dynamical properties of Rosenau-Hyman compactons using Padé approximants,” Physical Review E, vol. 81, no. 5, Article ID 056708, 2010.

A. Cardenas, B. Mihaila, F. Cooper, and A. Saxena, “Properties of compacton-anticompacton collisions,” Physical Review E, vol. 83, no. 6, Article ID 066705, 2011.

A. Ludu and J. P. Draayer, “Patterns on liquid surfaces: cnoidal waves, compactons and scaling,” Physica D, vol. 123, no. 1–4, pp. 82–91, 1998.

A. S. Kovalev and M. V. Gvozdikova, “Bose gas with nontrivial particle interaction and semiclassical interpretation of exotic solitons,” Low Temperature Physics, vol. 24, no. 7, article 484, 1998.

A. L. Bertozzi and M. Pugh, “The lubrication approximation for thin viscous films: regularity and long-time behavior of weak solutions,” Communications on Pure and Applied Mathematics, vol. 49, no. 2, pp. 85–123, 1996.

V. Kardashov, S. Einav, Y. Okrent, and T. Kardashov, “Nonlinear reaction-diffusion models of self-organization and deterministic chaos: theory and possible applications to description of electrical cardiac activity and cardiovascular circulation,” Discrete Dynamics in Nature and Society, vol. 2006, Article ID 98959, 16 pages, 2006.

Y. S. Kivshar, “Intrinsic localized modes as solitons with a compact support,” Physical Review E, vol. 48, no. 1, pp. R43–R45, 1993.

F. Cooper, J. M. Hyman, and A. Khare, “Compacton solutions in a class of generalized fifth-order korteweg-de vries equations,” Physical Review E, vol. 64, no. 2, Article ID 026608, 2001.

J. Kevorkian and J. D. Cole, Multiple Scale and Singular Perturbation Methods, Springer, New York, NY, USA, 1996.

Y. S. Kivshar and B. A. Malomed, “Dynamics of solitons in nearly integrable systems,” Reviews of Modern Physics, vol. 61, no. 4, pp. 763–915, 1989.

D. W. McLaughlin and A. C. Scott, “Perturbation analysis of fluxon dynamics,” Physical Review A, vol. 18, no. 4, pp. 1652–1680, 1978.

M. S. Ismail and T. R. Taha, “A numerical study of compactons,” Mathematics and Computers in Simulation, vol. 47, no. 6, pp. 519–530, 1998.

A. Chertock and D. Levy, “Particle methods for dispersive equations,” Journal of Computational Physics, vol. 171, no. 2, pp. 708–730, 2001.

F. Rus and F. R. Villatoro, “Adiabatic perturbations for compactons under dissipation and numerically-induced dissipation,” Journal of Computational Physics, vol. 228, no. 11, pp. 4291–4302, 2009.

F. Rus and F. R. Villatoro, “A repository of equations with cosine/sine compactons,” Applied Mathematics and Computation, vol. 215, no. 5, pp. 1838–1851, 2009.

A. Atangana and A. Secer, “A note on fractional order derivatives and table of fractional derivatives of some special functions,” Abstract and Applied Analysis, vol. 2013, Article ID 279681, 8 pages, 2013.

A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo, Theory and Applications of Fractional Differential Equations, Elsevier, Amsterdam, The Netherlands, 2006.

K. B. Oldham and J. Spanier, The Fractional Calculus, Academic Press, New York, NY, USA, 1974.

G. Jumarie, “On the solution of the stochastic differential equation of exponential growth driven by fractional Brownian motion,” Applied Mathematics Letters, vol. 18, no. 7, pp. 817–826, 2005.

A. Atangana and A. Secer, “The time-fractional coupled-Korteweg-de-Vries equations,” Abstract and Applied Analysis, vol. 2013, Article ID 947986, 8 pages, 2013.

R. Khalil, M. Al Horani, A. Yousef, and M. Sababheh, “A new definition of fractional derivative,” Journal of Computational and Applied Mathematics, vol. 264, pp. 65–70, 2014.

G. C. Wu, “New trends in the variational iteration method,” Communications in Fractional Calculus, vol. 2, no. 3, pp. 59–75, 2011.

C. M. Bender and S. A. Orszag, Advanced Mathematical Methods for Scientists and Engineers, Springer, New York, NY, USA, 1999.