Localized modes in time-fractional modified coupled Korteweg-de Vries equation with singular and non-singular kernels

Abstract: In this article, the modified coupled Korteweg-de Vries equation with Caputo and CaputoFabrizio time-fractional derivatives are considered. The system is studied by applying the modified double Laplace transform decomposition method which is a very effective tool for solving nonlinear coupled systems. The proposed method is a composition of the double Laplace and decomposition method. The results of the problems are obtained in the form of a series solution for 0 < α ≤ 1, which is approaching to the exact solutions when α = 1. The precision and effectiveness of the considered method on the proposed model are confirmed by illustrated with examples. It is observed that the proposed model describes the nonlinear evolution of the waves suffered by the weak dispersion effects. It is also observed that the coupled system forms the wave solution which reveals the evolution of the shock waves because of the steeping effect to temporal evolutions. The error analysis is performed, which is comparatively very small between the exact and approximate solutions, which signifies the importance of the proposed method.


Introduction
The time-fractional singular and non-singular operator are widely used in modern sciences and technology to study the behavior and applications of nonlinear ordinary and partial differential equations [1][2][3][4]. Fractional-order models are useful to study numerous real-world problems for longtime memory, and chaotic behaviour [5]. Due to these characteristics, the nonlinear models having fractional-order derivatives have been widely attracted in many areas of fractional calculus including image processing and signal, mechanics, biophysics and bioengineering, electrical engineering, biology, viscoelasticity, rheology, and control theory [6][7][8].
The Korteweg-de Vries (KdV) equation has extensively studied for nonlinear models describing the evaluation in time of long, unidirectional shallow water waves [9]. It was mainly presented by Boussinesq in 1877 and then retrieved by Diederik Korteweg and Gustav-de Vries in 1895 [10]. The existence solution to the a KdV equation can be seen in [11,12]. The KdV system having time-fractional derivatives attained remarkable attention in plasma physics especially in electrons acoustic waves (EAWs) propagation, due to its important role in studying diverse forms of mutual developments experimentally [13,14]. It has observed that the time-fractional operators in the system dramatically change the soliton amplitudes of the electron-acoustic solitary waves. This effects have been particularly equated with the structures of the broadband electrostatic noise experienced in the dayside auroral zone [15]. The system of coupled KdV system plays a leading part in various areas of sciences and engineering, particularly in water waves, quantum field theory, hydrodynamics and plasma physics [16,17]. It also represents the relations in extended waves with altered dispersion relations and describes iterations of water waves [18,19].
Here, we consider the mCKdV in the form [20] ∂ α φ ∂t α + 3φ 2 ∂φ ∂x − 3 2 where λ is a real number to be chosen accordingly. The proposed mCKdV Eq (1.1) has more interesting features than classical, because the operators will be defined by an integral which play a vital role in modern technology, engineering, plasma physics, hydrodynamics and quantum theory [21]. The nonlinear differential equations contain numerous fractional differential singular and non-singular operators such as Hilfer, Caputo, Caputo-Fabrizio, Riemann-Liouville, Antagana-Balenau in Caputo's sense [22][23][24]. These operators can be reduced in Caputo's form after some parametric addition. One can assume that the fractional operator could provide a power-law estimate of the local conduct of non-differentiable functions [25]. The Caputo operator possesses a power-law kernel and has restrictions to apply in modeling physical phenomena.
A modified laplace decomposition method (MLDM) is applied to Schrodinger-KdV equation in the sense of Atangana-Baleanu derivative in [26]. To deal successfully in such a situation, an alternative fractional operator possess a kernel with exponential decay has been introduced [27]. This novel operator is called the Caputo-Fabrizio (CF) operator which has a non-singular kernel. This operator is broadly applicable for modeling particular type of physical problems which follows the exponential decay law. Currently, mathematical and physical models having the CF operator have a remarkable development. The characteristics and applications of the above derivatives has been extensively studied (see [28][29][30][31][32][33][34][35] and the reference therein).
Some resent contribution of the time fractional order KdV equations have been studied by using different techniques [53][54][55][56][57]. The proposed method is an essential and effective approach to finding approximate solutions of the nonlinear models having time-fractional derivatives. In the non-linear model (1.1), the time fractional derivative has been taken in the form of Caputo's and Caputo-Fabrizio operator form which have a remarkable advantages in such physical models. The model is solved by using an affective method called modified double Laplace decomposition method (MDLDM) to obtain an approximate solution. The numerical solution of the model has been discuss in the form tables follow by its plots.
The rest of the paper is organized as follows: In Section 2, some main definitions, remarks, important results and a brief discussion of the proposed method is included. The convergence and uniqueness of the proposed method with the help of Banach contraction principle theorems are also studied in Section 3. In Section 4, the problem is discussed in Caputo's and Fabrizio forms and its solutions are obtained in a series form. We consider two numerical examples in Caputo's and Fabrizio form with the application of the proposed method, approximate solutions are obtained in the same section. This section also includes numerical discussion and figures for both the examples. Section 5 concludes the article followed by the Section 6 contain the future work in the manuscript, and Appendix contains the numerical values of the examples in the form of tables and also includes some parameters values.

Some basic definitions
In this section, we provide some basic definitions, lemmas and remarks regarding to the proposed method. Here also given some basics rules and definition related to double Laplace transform and decomposition method. Definition 1. [5,35,58] Caputo's fractional derivative of positive real order α > 0 of a function φ(x, t) is given by function φ : (0, ∞) → IR is defined by where n = [α] + 1, [α] denotes the integer part of a real number α, provided the right-hand side is point-wise defined on (0, ∞).
Definition 2. [35,59] Let φ(x, t) ∈ H 1 (a, b), with b > a and α ∈ (0, 1], then the fractional order in sense of Caputo-Fabrizio is define as where p and s are complex numbers.
Definition 4. [61,62] Applying the definition of double Laplace transform on fractional order derivative with respect t and x in Caputo sense is given by Definition 5. [61,62] Applying the definition of double Laplace transform on fractional order derivative with respect t and x in Caputo-Fabrizio sense is given by From the above definitions, we can conclude that whereφ (p, s) is an analytic function ∀p and s defined in the region by the inequalities Re (p) ≥ c and Re (s) ≥ d, where c, d ∈ R to be considered accordingly.

Series solutions of mCKdV using the proposed method
The proposed method MDLDM is hybrid method used widely for non-linear ordinary or partial differential equations. It is the combination of double Laplace and Adomian decomposition method to obtain an approximate solution for problems in hand. To discuss the analysis of the proposed method, we consider the following coupled non-linear problems: where L = c D α t (.) is the fractional-time derivative in Caputo's form, R 1 , R 2 are the operators contains the linear terms of Eq (1.1), N 1 , N 2 are the non-linear operators and f 1 (x, t), f 2 (x, t) are some external functions. Applying double Laplace on both sides to Eq (2.21), we obtain (2.22) Using the scheme defined above of double Laplace on the n th -derivative in Caputo's form, we obtain (2.24) Consider the series solution of the form 25) let N 1 and N 2 be the non-linear operators in the model are defined by where A n and B n are well known Adomian polynomials [50] of the functions φ 0 , φ 1 , φ 2 , . . . and ψ 0 , ψ 1 , ψ 2 , . . . respectively. For convince consider A n can be described by the following formula: Applying an inverse double Laplace on both sides to Eq (2.23): equating terms on both sides, we obtain In general, the following recursive formulas can be obtained, the final solutions can be obtained as for numerical purposes we take α = 1, the exact solution can be obtained as Similar procedures can be used to a problem having Caputo-Fabrizio fractional-time derivative to obtain a recursive relation.

Convergence analysis of the proposed method
In this section, we are discussing the existence and uniqueness of the solutions (2.29) and (2.30) of the general coupled time fractional non-linear partial differential equation (2.21). The following theorems are taken from [63].
where κ 1 , κ 2 and κ 3 are Lipschitz constants, φ andφ are the function's distinct values. Now we proceed as follow: The mapping is contraction under the condition 0 < σ < 1. As a result of the Banach fixed point theorem for contraction, Eq (2.21) has a unique solution. This marks the end of the proof.
Next we discuss convergence analysis of the problem.
Theorem 2. (Convergence theorem) The solution of Eq (2.21) in general forum will be convergence.
Proof. The Banach space of all continuous functions on the interval J with the norm φ =max t∈J φ is denoted as (C[J], . ). Define the {S n} sequence of partial sums i.e S n = n j=0 . With n ≥ m, let S n and S m be arbitrary partial sums. In this Banach space, we will show that S n is a Cauchy sequence. We obtain by employing a new formulation of Adomian polynomials.
Choosing n = m + 1, then , by using the following triangular inequality Now by definition 0 < σ < 1, we have 1 − σ n−m < 1, thus we have and also as |φ| < ∞ (φ is bounded), therefore, S n − S m → 0, hence S n is a Cauchy sequence in the Banach space B, so the series n j=0 φ j is convergent. Remark 1. Theorems 1 and 2 can also be applicable to the solution (2.30) for its uniqueness and convergent.

Applications of the proposed method
In this section, we illustrate some examples on Eq (1.1) with time-fractional derivatives in the form of Caputo and Caputo-Fabrizio operators in each case, and apply the proposed method discussed in this section with ( f 1 , f 2 ) = 0. Example 1. Consider the following (mCKdV) [20] in Caputo's form The exact solution of Eq (4.1) for α = 1 is [20] φ(x, t) = r 1 2k where Apply the modified double Laplace transform decomposition (MDLDM) scheme discussed in this section to Eq (4.1) and single Laplace transform to Eq (4.2), we obtain the following recurrence relations: for n ≥ 0. The Adomian polynomials A n , · · · F n , n = 1, 2, 3, for nonlinearity appears in Eq (4.4) are calculated by using the general rules shown in Eq (2.27) For more details, one can see [51,52,64]. Plugging Eq (4.5) in Eq (4.4), we obtain φ 0 = r 1 2k + k tanh(kx), ψ 0 = λ(r 1 + k) 2r 1 + r 1 tanh(kx), a 1 tanh(kx)(a 2 cosh(2kx)) + a 3 tanh(kx)) sech 4 (kx) , where the coefficients given above can be seen in Appendix. It should be noted that, other terms can be calculated in the similar way. The final solutions can be obtained For numerical purposes we take α = 1, the exact solution can be obtained as lim n→∞ (φ n , ψ n ) = (φ, ψ).
Example 2. Consider the following example having time fractional derivative in Caputo-Fabrizio (CF) form With the help of Eq (2.17) and the procedure used in Example 1, we can obtain the following series solutions in Caputo-Fabrizio form a 1 tanh(kx)(a 2 cosh(2kx)) + a 3 tanh(kx)) sech 4 (kx) , Other terms can be calculated in the similar way. The general solutions can be obtained in the form For numerical purposes we take α = 1, the exact solution can be obtained as lim n→∞ (φ n , ψ n ) = (φ, ψ).
The shock wave solutions φ(x, t) and ψ(x, t) in Eq (4.8) for the modified coupled Korteweg-de Vries (mCKdV) Eq 4.1 with Caputo's Fabrizio operators (CF) are depicted in Figure 4 Table 2.  [20] with the same parametric data for Example 2.

Conclusions
We have studied analytically the mCKdV with Caputo's and Caputo-Fabrizio (CF) derivatives using MDLDM. The proposed method is applied to the mCKdV system and approximate results are obtained in the form of series solutions by considering examples. The error analysis of the proposed method is also discussed. It is observed that the suggested scheme is one of the vigorous tools to investigate nonlinear problems. The main advantage of the suggested method is to analyze analytical solutions of the considered problem without using linearization and discretization. For validation, two examples in Caputo's and CF form are studied numerically and the results are compared with the numerical results with a physical interpretation. When non-integer order was employed in models including noninteger order derivatives, many scientists were less interested. For non-integer order fractional time derivative with singular and non-singular kernels of the mKdV equations, the MDLDM can be used. The MDLDM approach for handling nonlinear evolutionary equations was shown to have a broader applicability in this work. Physical meaning was given to the parameters in the resulting travelling wave solution.