Journal of Mathematics
Volume 2018 (2018), Article ID 2391697, 11 pages
https://doi.org/10.1155/2018/2391697

Research Article

Around Chaotic Disturbance and Irregularity for Higher Order Traveling Waves

Emile F. Doungmo Goufo¹ and Ignace Tchangou Toudjeu²

¹Department of Mathematical Sciences, University of South Africa, Florida 0003, South Africa
²Department of Electrical Engineering, Tshwane University of Technology, Pretoria 0183, South Africa

Correspondence should be addressed to Emile F. Doungmo Goufo; franckemile2006@yahoo.ca

Received 27 February 2018; Accepted 29 March 2018; Published 3 June 2018

Academic Editor: S. K. Q. Al-Omari

Copyright © 2018 Emile F. Doungmo Goufo and Ignace Tchangou Toudjeu. 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.

Many unknown features in the theory of wave motion are still captivating the global scientific community. In this paper, we consider a model of seventh order Korteweg–de Vries (KdV) equation with one perturbation level, expressed with the recently introduced derivative with nonsingular kernel, Caputo-Fabrizio derivative (CFFD). Existence and uniqueness of the solution to the model are established and proven to be continuous. The model is solved numerically, to exhibit the shape of related solitary waves and perform some graphical simulations. As expected, the solitary wave solution to the model without higher order perturbation term is shown via its related homoclinic orbit to lie on a curved surface. Unlike models with conventional derivative () where regular behaviors are noticed, the wave motions of models with the nonsingular kernel derivative are characterized by irregular behaviors in the pure factional cases (). Hence, the regularity of a soliton can be perturbed by this nonsingular kernel derivative, which, combined with the perturbation parameter of the seventh order KdV equation, simply causes more accentuated irregularities (close to chaos) due to small irregular deviations.

1. Introduction

In the last decade, a great number of researchers have paid a particular attention to the study of solitary wave equations that undergo the influence of external perturbations. Most of physical dynamics related to the movement of liquids and waves are governed by Korteweg–de Vries (KdV) equation and its variants. Hence, for this particular equation, Cao et al. [1] as well as Grimshaw and Tian [2] have recently shown that a force combined together with dissipation can provoke a chaotic behavior usually detectable by other analysis like phase plane analysis or nonnegative Liapunov exponents. KdV equation and its variants are of infinite dimension and their use to address traveling waves or chaotic dynamics of low dimension is facilitated by numerical approximations, which have proven that correlation dimension established via Grassberger-Procaccia technique and information dimension obtained from formula of Kaplan-Yorke are both between two and three for steady traveling waves [3].

However, many authors (like, e.g., [4–7]) preferred to use numerical approach to analyze the KdV equation or its variants, especially the one with many levels of perturbations. Hence, it was shown in [4] that there is no periodic waves for the autonomous Korteweg–de Vries–Burgers equation of dimension two. We follow, in this paper, the same trend of numerical approach by making use of the recently developed fractional derivative with nonsingular kernel [8–13], to express a seventh order Korteweg–de Vries (KdV) equation with one perturbation level. This is the first instance where such a model is extended to the scope of fractional differentiation and fully investigated. We prove existence and uniqueness of a continuous solution. Before that, we shall give in the following section a brief review of the recent developments done in the theory of fractional differentiation.

2. Around the Nonsingular Kernel Differentiation with Fractional Order

The concept of fractional order derivative is seen by many authors as a great endeavor to ameliorate nonlinear mathematical models, widen their analysis, and expand their interpretation. Today’s literature of the concept has been enriched with many innovative definitions more related to the complexity and diversity of natural phenomena surrounding us. There are fractional order derivatives of local type and also of nonlocal type [8, 14–17]. The Caputo’s definition

remains the most commonly used in the applied science, followed by Riemann-Liouville’s version given by

.

Recent observations by Caputo and Fabrizio [8] stated that the two definitions above better describe physical processes, related to fatigue, damage, and electromagnetic hysteresis, but do not genuinely depict some behavior taking place in multiscale systems and in materials with massive heterogeneities. Hence, the same authors introduced the following new version of fractional order derivative with no singular kernel:

Definition 1 (Caputo-Fabrizio fractional order derivative (CFFD)). Let be a function in ; ; ; then, the Caputo-Fabrizio fractional order derivative (CFFD) is defined as

where is a normalization function such that .

Remark 2. Caputo and Fabrizio [8] substituted the kernel appearing in (1) when by the function and by This immediately removes the singularity at that exists in the previous Caputo’s expression.

For the function that does not belong to , the CFFD is given by

Losada and Nieto [13] upgraded this definition of CFFD by proposing the following:

Unlike the classical version of Caputo fractional order derivative [14, 18], the CFFD with no singular kernel appears to be easier to handle. Furthermore, the CFFD verifies the following equalities:

with any suitable function and the starting point of the integrodifferentiation. The fractional integral related to the CFFD and proposed by Losada and Nieto reads as

This antiderivative represents sort of average between the function and its integral of order one. The Laplace transform of the CFFD reads

where is the Laplace transform of

Definition 3 (New Riemann-Liouville fractional order derivative (NRLFD)). As a response to the CFFD and being aware of the conflicting situations that exist between the classical Riemann-Liouville and Caputo derivatives, the classical Riemann-Liouville definition was modified [9, 10] in order to propose another definition known as the new Riemann-Liouville fractional derivative (NRLFD) without singular kernel and expressed for as

Again, the NRLFD is without any singularity at compared to the classical Riemann-Liouville fractional order derivative and also verifies

Compared to (7), we note here the exact correspondence with at The Laplace transform of the NRLFD reads as [9, 10]

Other versions and innovative definitions of fractional derivatives have since been introduced. This paper however uses the CFFD, so for more details about those recent definitions, please feel free to consult the articles and works mentioned above and also the references mentioned therein.

3. Existence and Uniqueness

In this section we prove the existence and uniqueness results for the seventh order Korteweg–de Vries equation (KdV) with one perturbation level, expressed with the CFFD and given by

assumed to satisfy the initial condition

where is the perturbation parameter and is the Caputo-Fabrizio fractional order derivative (CFFD) given in (5). Existence results for the model (13)-(14) here above are established by making use of the expression of the antiderivative (8). This yields

This can be rewritten as

Set now

The next step is to look for a real constant such that

In fact

Well known properties for the norms give

Keeping in mind that and are bounded functions, then there exists real numbers and such that

Set ; hence,

Therefore, the Lipschitz condition holds for the partial derivatives and and there exists a real constant such that

where the bounded condition (14) has been exploited, whence

with

This proves that satisfies the Lipschitz condition and then, it allows us to state the following proposition.

Proposition 4. If the condition holds, then, there exists a unique and continuous solution to the seventh order Korteweg–de Vries equation with one perturbation level expressed with the CFFD given in (5):

Proof. Let us go back to the model (16) rewritten as

which yields the recurrence formulation given as follows:

Let

and it can be shown that is a continuous solution. Indeed, if we take

then, it is straightforward to see that

More explicitly, we have

Passing this equation to the norm yields

Applying the Lipschitz condition to gives

which can be rewritten as

After integrating, we make use of well known properties of the recursive technique from (35) to have

with , which explicitly shows the existence of the solution and that it is continuous.
The step forward is to prove that the solution of the model (26) is given by the function

For that, let

Making use of (29), we should have In other terms, the gap that exists between and should vanish as Consider

giving

Hence, and from the right hand side, we have

Just take as the solution of (26) that is continuous. Moreover, the Lipschitz condition for yields

This yields

considering the initial condition and taking the limit as gives

Uniqueness. To prove that the solution is unique, we take two different functions and that satisfy the model (26); then,

equivalently

This yields if

where we have used the Lipschitz condition for and this ends the proof.

4. Shape of Solitary Waves via Numerical Approximations

4.1. Shape of Solitary Waves for the Lower-Order Approximation

In this section, we are interested in waves traveling to a specific direction, and then, we consider solutions to the seventh order KdV equation

We start by considering the conventional case where to have, using , the following model:

We investigate the traveling waves taking the form , where is the speed of the wave. We assume that does not depend on independently from but rather depends on the combined variable We also assume that the wave dies at infinity, meaning

Now, it is possible to transform the seventh order KdV equation (49) into an ordinary differential equation (ODE) by making use of the basic properties of differentiation. Then, and which yield the following ODE:

The -integration of this equation once gives

where we have ignored the constant of -integration that is null due to boundary conditions (50). If the higher order perturbation term is ignored, then we have

Then we can solve numerically this equation by transforming it into a system of four ODEs of order one as follows:

Numerical simulations are done in the phase-space () as shown in Figure 1. Let us now come back to the full model (48) with the nonsingular kernel derivative CFFD (given in (5)) and with no higher order perturbation parameter , given as

This fractional equation is solved numerically by making use of the Adams-Bashforth-Moulton type method also known as predictor-corrector (PECE) technique and is fully detailed in the article by Diethelm et al. [19]. Figures 1 and 2 represent the numerical simulations of solutions to the fractional model (55) with different values of the derivative order They clearly point out relative irregularities when the derivative order of the CFFD is

(a)

(b)

Figure 1: Numerical plot in the phase-space for the model KdV model (55) or model (48) with no higher order perturbation term, for (conventional case). As expected in (a), the soliton solution is shown via its related homoclinic orbit to lie on a curved surface. In (b), the projection on the plane .

(a)

(b)

Figure 2: Numerical plot in the phase-space for the model KdV model (55) with no higher order perturbation term, for As expected in (a), even in the pure fractional case, the soliton solution is shown via its related homoclinic orbit to lie on a curved surface, with some irregular movements compared to Figure 1. In (b), the projection on the plane .

4.2. Shape of Solitary Waves for the Higher Order Approximation

Now, the perturbation term is considered to have (48):

subject to localized initial condition

Although we proved existence and uniqueness for this nonlinear problem, it still faces the challenge of providing an explicit expression of exact solution or approximated solution. It is almost impossible to use some analytical methods, like for instance, integral transform methods, the Green function technique, or the technique separation of variables. Hence, a semianalytical method like the Laplace iterative method can be a valuable tool to provide a special solution to the model (56)-(57) as shown in the following lines.

We start by applying the Laplace transform (9) on both sides of (56) to obtain

Let ; then,

Application of inverse Laplace transform on both sides yields

This leads to the following recursive system

then making the solution to be Numerical approximations are performed according to the following steps:

(i) is considered as initial input: Choose as the number of terms to compute. Name to be the approximate solution. Set and (ii)For the other terms use . Compute . to finally get

where

Note we can use , to obtain following results and some numerical approximations with related absolute errors are summarized in Tables 1 and 2 (which present some numerical approximations with related absolute errors as given by the sequences (63) and (65), respectively, representing the pure fractional case () and the standard case ()).

Note that, for , we recover for the conventional model (49)

the following well known result

This result corresponds to the similar one obtained in [20, 21] or in [22] via Adomian decomposition method. The other terms were computed following the same iterative approach and then, the functions admits the closed form found to be

The shape of numerical approximations is depicted in Figures 3–6, plotted for different values of the derivative order and showing a spot of irregular perturbed motion in the pure fractional case () compared to the conventional case (). Before going to conclusions, we can summarize the physical aspects of this work (issued from the mathematical analysis and simulations) as follows: analysis performed on the KdV model with no higher order perturbation term has shown the soliton solution via its related homoclinic orbit lying on a curved surface (Figures 1 and 2), and this remains true in the conventional case as well as the pure fractional case. However, irregular movements are more important in the latter case. Analysis performed on the seventh order KdV equation with one perturbations level has shown the shape of solitary waves characterized by motions with more irregular behaviors tending to become chaotic. The chaos is more likely to happen in pure fractional case compared to the conventional case (Figures 3–6). Hence, the main results of this paper, reflected by the six figures, insinuate that the regularity of a soliton can be perturbed by the nonsingular kernel derivative, which, combined with the perturbation parameter of the KdV model, may lead to chaos.

Table 1: Some values for numerical and exact solutions to the seventh order KdV equation expressed with the CFFD given in (5) at , , and

For

Time ()
Spatial ()
Exact value
Numerical value
Error made

0.5
−15.0

−7.0

0.0
0.100372300
0.100470000

+7.0

+15.0

1.0
−15.0

−7.0

00

0.0
0.100390000
0.100380000

+7.0

+15.0

2.0
−15.0

−7.0

0.0
0.100410000
0.100510000

+7.0

+15.0

4.0
−15.0

−7.0

0.0
0.100430000
0.090430000

+7.0

+15.0

8.0
−15.0

−7.0

0.0
0.090400000
0.100400000

+7.0

+15.0

10.0
−15.0

−7.0

0.0
0.090402000
0.091402000

+7.0

+15.0

Table 2: Some values for numerical and exact solutions to the seventh order KdV equation expressed with the CFFD given in (5) at , , and

For

Time ()
Spatial ()
Exact value
Numerical value
Error made

0.5
−15.0
0.000099315
0.0001 094 20

−7.0
0.593670000
0.593670000

0.0
0.151107000
0.151107000

+7.0
0.06971 8000
0.06971 8000

+15.0
0.000 167 700
0.00009 6770

1.0
−15.0
0.00009 159 2
0.00008 857 2

−7.0
0.05477 9000
0.05477 9000

0.0
0.151907000
0.151907000

+7.0
0.07 554 6000
0.07 554 6000

+15.0
0.0001 266 20
0.00002 6620

2.0
−15.0
0.00007790 1
0.000237900

−7.0
0.04 6635000
0.04663 5000

0.0
0.152007000
0.152140000

+7.0
0.08 868 8000
0.08 868 8000

+15.0
0.0001 488 70
0.0001591 70

4.0
−15.0
0.00005 635 2
0.001316 400

−7.0
0.03 378 5000
0.03 378 5000

0.0
0.152060000
0.162060000

+7.0
0.122 130000
0.120 120000

+15.0
0.0002057 90
0.0002057 90

8.0
−15.0
0.00002 948 8
0.001 289 500

−7.0
0.01 771 1000
0.01 931 1000

0.0
0.152060000
0.162060000

+7.0
0.230 430000
0.231 630000

+15.0
0.0003 932 70
0.0004 033 70

10.0
−15.0
0.00002 133 1
0.00002 007 1

−7.0
0.01 281 9000
0.01283 0000

0.0
0.161760000
0.162760000

+7.0
0.315 320000
0.315 320000

+15.0
0.0005 436 40
0.0005 436 40

Figure 3: The shape of approximated solution in the conventional case (). It is plotted with respect to different fixed values of time and space at , , The motion is quite standard and regular.

Figure 4: The shape of approximated solution in the pure fractional case (), plotted with respect to different fixed values of time and space at , , The motion is shown to behave irregularly compared to the conventional case and the dynamic tending to a chaotic one.

Figure 5: The shape of approximated solution in the conventional case (). It is plotted with respect to different fixed values of time and space at , , The motion is quite standard and regular.

Figure 6: The shape of approximated solution in the pure fractional case (), plotted with respect to different fixed values of time and space at , , The motion is shown to behave irregularly compared to the conventional case and the dynamic tending to a chaotic one.

5. Concluding Remarks

We have made use of the recent version of derivative with nonsingular kernel to prove existence and uniqueness results for a model of seventh order Korteweg–de Vries (KdV) equation with one perturbation level. The unique solution is continuous. We then use numerical approximations to evaluate the behavior of the solution under the influence of external factors. It happened that when applied to equations of wave motion like the seventh order KdV equation, the new derivative acts as one of those external factors, which, combined with the perturbation term of the model, causes the solution to be more irregular and unpredictable. This is the first instance where such a model is fully investigated and such result is exposed. This work differs from the previous ones within introduction of the nonsingular kernel derivative into a powerful model like Korteweg–de Vries’s, which reveals another interesting feature that exists in the domain of wave motion as well as chaos theory.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was partially supported by the National Research Foundation (NRF) of South Africa, Grant no. 105932.

References

Q. Cao, K. Djidjeli, W. G. Price, and E. H. Twizell, “Periodic and chaotic behaviour in a reduced form of the perturbed generalized Korteweg-de Vries and KADomtsev-Petviashvili equations,” Physica D: Nonlinear Phenomena, vol. 125, no. 3-4, pp. 201–221, 1999.

R. Grimshaw and X. Tian, “Periodic and chaotic behaviour in a reduction of the perturbed Korteweg-de Vries equation,” Proceedings of the Royal Society A Mathematical, Physical and Engineering Sciences, vol. 445, no. 1923, pp. 1–21, 1994.

X. Tian and R. H. Grimshaw, “Low-dimensional chaos in a perturbed Korteweg-de Vries equation,” International Journal of Bifurcation and Chaos, vol. 5, no. 4, pp. 1221–1233, 1995.

Z. Feng and R. Knobel, “Traveling waves to a Burgers-Korteweg-de Vries-type equation with higher-order nonlinearities,” Journal of Mathematical Analysis and Applications, vol. 328, no. 2, pp. 1435–1450, 2007.

H. Grad and P. N. Hu, “Unified shock profile in a plasma,” Physics of Fluids, vol. 10, no. 12, pp. 2596–2602, 1967.

P. N. Hu, “Collisional theory of shock and nonlinear waves in a plasma,” Physics of Fluids, vol. 15, no. 5, pp. 854–864, 1972.

R. S. Johnson, “A non-linear equation incorporating damping and dispersion,” Journal of Fluid Mechanics, vol. 42, pp. 49–60, 1970.

M. Caputo and M. Fabrizio, “A new definition of fractional derivative without singular kernel,” Progress in Fractional Differentiation and Applications, vol. 1, no. 2, pp. 1–13, 2015.

E. F. Doungmo Goufo, “Chaotic processes using the two-parameter derivative with non-singular and non-local kernel: basic theory and applications,” Chaos: An Interdisciplinary Journal of Nonlinear Science, vol. 26, no. 8, Article ID 084305, 10 pages, 2016.

E. F. D. Goufo, “Application of the Caputo-Fabrizio fractional derivative without singular kernel to Korteweg-de Vries-Burgers equation,” Mathematical Modelling and Analysis, vol. 21, no. 2, pp. 188–198, 2016.

Y. Khan, S. Panjeh Ali Beik, K. Sayevand, and A. Shayganmanesh, “A numerical scheme for solving differential equations with space and time-fractional coordinate derivatives,” Quaestiones Mathematicae, vol. 38, no. 1, pp. 41–55, 2015.

R. Kandasemy, S. Raut, D. Verma, and G. There, “Design of solar tricycle for handicapped person,” Journal of Mechanical and Civil Engineering, vol. 25, no. 2, pp. 11–24, 2013.

J. Losada and J. J. Nieto, “Properties of the new fractional derivative without singular Kernel,” Progress in Fractional Differentiation and Applications, vol. 1, no. 2, pp. 87–92, 2015.

M. Caputo, “Linear models of dissipation whose Q is almost frequency independent-II,” The Geophysical Journal of the Royal Astronomical Society, vol. 13, no. 5, pp. 529–539, 1967.

E. F. Doungmo Goufo, “A biomathematical view on the fractional dynamics of cellulose degradation,” Fractional Calculus and Applied Analysis, vol. 18, no. 3, pp. 554–564, 2015.

E. F. Doungmo Goufo, “Solvability of chaotic fractional systems with 3D four-scroll attractors,” Chaos, Solitons & Fractals, vol. 104, pp. 443–451, 2017.

Y. Khan, K. Sayevand, M. Fardi, and M. Ghasemi, “A novel computing multi-parametric homotopy approach for system of linear and nonlinear Fredholm integral equations,” Applied Mathematics and Computation, vol. 249, pp. 229–236, 2014.

J. Prüss, Evolutionary Integral Equations and Applications, vol. 87 of Monographs in Mathematics, Birkhäuser, Basel, Switzerland, 1993.

K. Diethelm, N. J. Ford, A. D. Freed, and Y. Luchko, “Algorithms for the fractional calculus: a selection of numerical methods,” Computer Methods Applied Mechanics and Engineering, vol. 194, no. 6-8, pp. 743–773, 2005.

B. R. Duffy and E. J. Parkes, “Travelling solitary wave solutions to a seventh-order generalized KdV equation,” Physics Letters A, vol. 214, no. 5-6, pp. 271-272, 1996.

R.-X. Yao and Z.-B. Li, “On t/x-dependent conservation laws of the generalized nth-order KdV equation,” Chinese Journal of Physics, vol. 42, no. 4, pp. 315–322, 2014.

S. M. El-Sayed and D. Kaya, “An application of the ADM to seven-order Sawada-Kotara equations,” Applied Mathematics and Computation, vol. 157, no. 1, pp. 93–101, 2004.