Abstract

After the recent introduction of the Caputo-Fabrizio derivative by authors of the same names, the question was raised about an eventual comparison with the old version, namely, the Caputo derivative. Unlike Caputo derivative, the newly introduced Caputo-Fabrizio derivative has no singular kernel and the concern was about the real impact of this nonsingularity on real life nonlinear phenomena like those found in shallow water waves. In this paper, a nonlinear Sawada-Kotera equation, suitable in describing the behavior of shallow water waves, is comprehensively analyzed with both types of derivative. In the investigations, various fixed-point theories are exploited together with the concept of Piccard -stability. We are then able to obtain the existence and uniqueness results for the models with both versions of derivatives. We conclude the analysis by performing some numerical approximations with both derivatives and graphical simulations being presented for some values of the derivative order Similar behaviors are pointed out and they concur with the expected multisoliton solutions well known for the Sawada-Kotera equation. This great observation means either of both derivatives is suitable to describe the motion of shallow water waves.

1. Introduction

The concept of derivative with generalized order has been intensively used to model many real life phenomena with the general formulation given in many works [14]. The version of Caputo fractional derivative (CD) [5] remains the most popular thanks to its useful properties suitable to describe real life phenomena [510], especially in those phenomena related to the dynamics of diseases [1114]. However, to address the issues caused by the singular kernel noticed in Caputo fractional derivative, Caputo and Fabrizio recently proposed another fractional derivative with no singularity called the Caputo-Fabrizio fractional derivative (CFD). Since then, many authors used the CFD to solve some real world problems as shown in some works [7, 15, 16]. The CFD was modified in [7] to generate the related Riemann-Liouville version. Recall that the Caputo-Fabrizio fractional derivative reads as the convolution of the local derivative and the function exponential. The concept was generalized by introducing the two-parameter fractional derivative with nonlocal and nonsingular kernel [17, 18], like, for instance, the Goufo-Caputo derivative with two parameters and the Goufo-Riemann-Liouville derivative with two parameters proven to have a huge impact in Chaos theory and image processing. Investigations on fractional derivatives like the CFD are still ongoing around the world and this paper follows the same momentum by providing an analysis of the Sawada-Kotera equation within the scope of both Caputo and Caputo-Fabrizio derivatives.

It is important to notice that the Sawada-Kotera equation is a nonlinear differential equation that falls into the class of Kortewegde-de Vries- (KdV-) type equations. It has applications in waves theory, like, for instance, in shallow water waves, ion-acoustic waves, or stratified internal waves. It is also classified in the hierarchy of equations suitable to be solved via Lax operator [1921]. Hence, the nonlinear Sawada-Kotera type model is comprehensively investigated in this work using the CD and CFD defined, respectively, by, and, where is a normalization function such that

The well-known antiderivative associated with the CD is given aswhile the antiderivative associated with the CFD is defined as [15, 16] The latter antiderivative when applied to a function gives a sort of average between the function and its first-order integral.

As we mentioned earlier, the CFD is different form the CD in the sense that it presents a kernel with no singularity at Furthermore, the following relations hold for the CFD:

2. Analysis Using the Caputo Fractional Derivative (CD)

Existence of the exact solution of the nonlinear Sawada-Kotera equation modeled with Caputo time fractional derivative is proven in this section. Recall that similar analysis was performed in other works [22, 23]. We focus our attention on the following Sawada-Kotera model expressed assubject to the initial conditionwhere represents the CD as defined in (1). The first step to follow in this analysis is to transform (6) into an integral equation. Hence, we apply on both sides the antiderivative (3) to getFor reasons of simplicity, the following operator is considered:The next stage is to prove that the Lipschitz condition with respect to holds for the operator Thus,Exploiting the norm’s properties, the above yieldsExploiting the Lipschitz condition for the standard derivative function, there are positive constants , , , , , , , ; and and such that Substituting (12) into (11) leads toTaking thus Hence, the desired Lipschitz conditions are proven. Making use of the results above, it is possible to evaluate the quantity

Considering the fact that we are analyzing a real world phenomenon, we must assume existence of initial conditions that are bounded, whencewhich gives following theorem.

Theorem 1. Take and assume that given in (7) is a bounded function. Let (with ) be a continuous with respect to for any fixed in ,
If , then satisfies model (6) if and only if satisfies the corresponding Volterra integral model (8).

Proof. We start with the proof of the necessity condition by assuming that verifies model (6). Since for all , then (6) means that it is possible to find the Caputo fractional derivative of in . HoweverExploiting existing results in [9] for , then Again applying the results in [9] for to then we obtainwith Integration by parts performed in (21), followed by differentiating the resulting expression, leads toChanging of variable with yieldsObviously, and . Then, models (23) and (21) take the formBecause of and exploiting the Lipschitz condition of , we get However, applying on both sides of (6) and considering the initial condition, we obtain the desired Volterra version, which proves the necessity condition.
Conversely, we prove the sufficient condition by assuming that verifies the Volterra version of (6). We have to show that the initial condition of (6) holds for If we differentiate both sides of the Volterra version, then The similar change of variable applied in the Volterra version for yields Considering the limit as and exploiting the fact that is continuous prove that the required the initial conditions hold, which shows the sufficiency condition.

Theorem 2. Since the condition is satisfied and the Lipschitz condition holds for the operator J, there is a unique solution for model (8) in the space

Proof. The proof begins by proving that it is possible to find a unique solution of model (6) in . By Theorem 2 it will be sufficient to point out the existence of a unique solution of the nonlinear Volterra integral equation (8). Recall that (8) makes sense when it is applied in any interval Let us then choose an appropriate such that and then prove the existence of a unique solution To perform this, we can exploit the method of successive approximation by puttingObviously, , and the differentiation of (30) with respect to leads to Due to the fact that is differentiable with respect to , then
At this stage, we can estimate for ; thus, Also we have that Proceeding by repeating the above steps -times, we get the inequality Hence, it happens that the sequence approaches a limit function . Also thanks to the assumption However, choosing , we estimate Nevertheless, considering the Lipschitz condition of , we have The equations obtained above prove that is the solution of (8) in the space .
Uniqueness. Here we assume that it is possible to find two different solutions and verifying model (8) on : The above leads to and is a contradiction compared to the assumption. Thus, there is a unique solution
Now, let us consider the closed interval so that with and Hence, for , we get Considering the above proof of uniqueness on , we certainly have where is a known function.
The whole analysis and discussion presented above lead to the conclusion that there exists a unique solution of (8). Considering the interval so that with and , it is possible to repeat the same steps as above and finally prove the existence of a unique solution of (8) and, therefore, existence for (6) in the space , which proves the theorem.

Theorem 3. Since satisfies Theorem 1, then if holds, then any sequence , as defined above is an approximation of the exact solution with

Proof. The proof is done by induction on .
For recursively, we obtain butand the proof is complete.

The above theorem also shows a rapid convergence of to

3. Analysis Using the Caputo-Fabrizio Fractional Derivative (CFD)

3.1. Definition and Preliminary Concepts

Let be a Banach space and consider as a self-map of . Let be a recursive method where we assume that is a set representing all the fixed points of the map and is assumed to have at least one element. Assume that tends to an element Let and define .

implies that Then, with the above settings, we say that the recursive relation is -stable. Moreover, with the assumption that has an upper boundary and if all the conditions above hold for , then it is called Picard’s iteration and is -stable.

Lemma 4 (see [24]). Let be a Banach space and a self-map of satisfying for all , where , With the assumption that has a fixed point , then is Piccard -stable.

3.2. Analysis

Theorem 5. Consider the self-map and redefine model (6) as the recurrence relation: Then, it is -stable in .

Proof. We start by showing that has a fixed-point. Let ; then, where is the fractional integral (antiderivative) associated with the Caputo-Fabrizio fractional derivative defined in (4). Then, where the Lipschitz condition for the differential operator has been used together with the following positive constants ; ; and and Then, where and this proves that the Lipschitz condition holds for the nonlinear operator and hence, it has a fixed point.
Lastly, if we take and , then the conditions of Lemma 4 hold for which is therefore Piccard -stable and the proof is complete.

Remark 6. The previous theorem shows the existence of a unique solution for model (6) obtained via the fixed point of in the recurrence relation (47).

4. Numerical Schemes

Numerical methods are performed in this section for model with both derivatives. Before that, we show how the solutions can be derived for both types of models with the use of some recursive techniques together with Laplace transform (for the model with CD) and Sumudu transform (for the model with CFD).

4.1. Model with Caputo Derivative

To proceed with the solvability, consider model (6) and apply the Laplace transform on both sides. We obtain equivalently, The inverse Laplace transform gives At this stage we can introduce the following iterative system: This system is exploited to determine an approximation for the solution of the model with Caputo derivative. The stability of the method is well known since it uses the fixed-point theorem. Whence, this approximation reads as

4.2. Model with Caputo-Fabrizio Derivative

Similarly, to proceed with the solvability, we will consider model (6) and apply the Sumudu transform on both sides. But first of all, we have to recall the well-known relation that exists between the Sumudu transform and the Caputo-Fabrizio derivative given as where is the Sumudu transform of Applying the Sumudu transform on both sides of (6) gives After rearrangement The inverse Sumudu transform leads to Similarly as above, we can introduce the following iterative system:

Again, this system is exploited to determine an approximation for the solution of the Caputo-Fabrizio model. The stability of the method is well known and is established using the fixed-point theorem and the approximation is given by

Exploiting the recurrence schemes above, numerical simulations for models with both derivatives, namely, the Caputo derivative (CD) and Caputo-Fabrizio derivative (CFD), are depicted in Figure 1 for different values of the derivative order equal to and There are clear similarities in the observed behaviors in Figure 1(a) (model with the CD) and Figure 1(b) (model with the CFD) both represented for We note existence of the standard multisoliton solutions well known to be related to Sawada-Kotera equation. The same analysis holds for Figure 1(c) (model with the CD) and Figure 1(d) (model with the CFD) all drawn for

5. Concluding Remarks

The prime goal of this work has been to exploit the concept of Piccard -stability as well as various fixed-point theories, Sumudu and Laplace transforms to address the solvability of the nonlinear Sawada-Kotera model expressed with both Caputo and Caputo-Fabrizio derivatives with fractional order. The analysis performed in this paper has enabled us to show existence and uniqueness results for both models and compare the behaviors of their solutions. As expected, numerical simulations done for some values of the derivative order show similar behaviors of the solutions and point out existence of multisoliton solutions for both equations. This comparative analysis is more than useful in the sense that it validates the use of any of both derivatives in the description of natural phenomena like those related to shallow water waves or ion-acoustic waves and will surely yield more similar and complex investigations within the scientific community.

Competing Interests

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