Numerical investigation of fractional-order Kersten–Krasil’shchik coupled KdV–mKdV system with Atangana–Baleanu derivative

In this article, we present a fractional Kersten–Krasil’shchik coupled KdV-mKdV nonlinear model associated with newly introduced Atangana–Baleanu derivative of fractional order which uses Mittag-Leffler function as a nonsingular and nonlocal kernel. We investigate the nonlinear behavior of multi-component plasma. For this effective approach, named homotopy perturbation, transformation approach is suggested. This scheme of nonlinear model generally occurs as a characterization of waves in traffic flow, multi-component plasmas, electrodynamics, electromagnetism, shallow water waves, elastic media, etc. The main objective of this study is to provide a new class of methods, which requires not using small variables for finding estimated solution of fractional coupled frameworks and unrealistic factors and eliminate linearization. Analytical simulation represents that the suggested method is effective, accurate, and straightforward to use to a wide range of physical frameworks. This analysis indicates that analytical simulation obtained by the homotopy perturbation transform method is very efficient and precise for evaluation of the nonlinear behavior of the scheme. This result also suggests that the homotopy perturbation transform method is much simpler and easier, more convenient and effective than other available mathematical techniques.


Introduction
Many researchers have been working on various aspects of fractional derivatives in recent years. Caputo and Fabrizio modified the existing Caputo derivative to develop the Caputo-Fabrizio fractional derivative [1][2][3][4][5] based on a nonsingular kernel. Because of its advantages, numerous researchers utilized this operator to investigate various types of fractional-order partial differential equations [6][7][8][9]. To address this issue, Atangana and Baleanu proposed a new fractional operator called Atangana-Baleanu derivative, which combines Caputo and Riemann-Liouville derivatives. Because of the existence of the Mittag-Leffler kernel, which is a generalization of the exponential kernel, this new Atangana-Baleanu derivative has a long memory. Moreover, the Atangana-Baleanu operator outperforms other operators, and different scientific models have been successfully solved. Many advances have been made in fractional calculus over the last few years by borrowing ideas from classical calculus, but it does not remain easy. Scholars have the main concern to obtain a numerical solution; for this, numerous efficient methodologies have been constructed for fractional differential equations, such as the Adomian decomposition transform method [10], variational iteration transform method [11,12], optimal homotopy asymptotic method [13], the homotopy perturbation method [14,15], etc.
Fractional coupled systems are widely applied to study complex behavior of plasma contains multi components such as ions, free electrons, atoms, etc. Many researchers made efforts to study this behavior numerically. In this direction, recently Paul Kersten and Joseph Krasil'shchik studied KdV equation and modified KdV equation and proposed absolute complexity between coupled KdV-mKdV nonlinear systems for studying the behavior of nonlinear systems [16][17][18][19]. Numerous variations of this Kersten-Krasil'shchik coupled KdV-mKdV nonlinear system have been introduced by many researchers [20][21][22][23][24][25][26][27]. Among these variations, the mathematical model for describing the behavior of multi-component plasma for waves propagating in positive χ axis, known as nonlinear fractional Kersten-Krasil'shchik coupled KdV-mKdV system, is given by where η is temporal coordinate and χ is spatial coordinate. The factor σ represents order of the fractional operator. This operator is studied in the Caputo form. When σ = 1, a fractional coupled system converts to the classical system as follows: If we put φ = 0, then the Kersten-Krasil'shchik coupled KdV-mKdV system converts into the well-known KdV system If we put ψ = 0, then the Kersten-Krasil'shchik coupled KdV-mKdV system converts into the well-known modified KdV system In view of that, the Kersten-Krasil'shchik coupled KdV-mKdV system can be assumed to be a combination of a KdV system and a mKdV system represented by (2) to (4).
In this study, we also consider a fractional nonlinear two component homogeneous time fractional coupled third-order KdV system as follows: where η is temporal coordinate and χ is spatial coordinate, σ is a factor representing the order of the fractional operator. This operator is studied in the Caputo form. When σ = 1, a fractional coupled system converts to the classical system He [28][29][30] proposed the homotopy perturbation method for solving numerous linear and nonlinear initial and boundary value problems. Several researchers have looked into using the homotopy perturbation method to solve nonlinear equations in engineering and science [31][32][33]. Because of the difficulties caused by nonlinear terms, the Laplace transform is completely incapable of handling nonlinear equations. Recently, various methods for dealing with such nonlinearities, such as the Laplace decomposition algorithm [34,35] and the homotopy perturbation transform method (HPTM) [36,37], have been proposed to produce highly effective techniques for solving many nonlinear problems. The primary goal of this paper is to use an effective homotopy perturbation method modification to overcome the deficiency. For solving the system of KdV equations, we use the homotopy perturbation transform method. All conditions can be met using this method. One or two iteration steps also yield very accurate results over a wide range. The proposed homotopy perturbation transform method provides the solution in a rapid convergent series, which may lead to a closed solution [38][39][40][41].

Basic definitions
Definition 2.1 The fractional-order Caputo derivative is defined by

Definition 2.2
The Laplace transformation connected with fractional Caputo derivative LC D σ {f ( )} is expressed by

Definition 2.3 In the Caputo sense the Atangana-Baleanu derivative is defined as
and E σ represents the Mittag-Leffler function.

Definition 2.4
The Atangana-Baleanu derivative in the Riemann-Liouville sense is defined as

Definition 2.5 The Laplace transform connected with the Atangana-Baleanu operator is defined as
. Definition 2.6 Consider 0 < σ < 1, and f is a function of σ , then the fractional-order integral operator of σ is given as

The general methodology of HPTM
In this section, the HPTM for the general form of FPDEs with the initial condition where is D σ η ψ(χ, η) = ∂ σ ∂η σ the Caputo fractional derivative of order σ , M, and N , are linear and nonlinear functions, respectively, and h is a source operator. Using the Laplace transform of Eq. (7), we have Now, by taking the inverse Laplace transform, we get where Now, the perturbation method parameter p is defined as where perturbation term p ∈ [0, 1].
The nonlinear functions can be defined as where H i are He's polynomials of ψ 0 , ψ 1 , ψ 2 , . . . , ψ i and can be determined as Putting equations (13) and (14) in equation (11), we have Both sides comparison coefficient of p, we have . . .

Hence a series solution is given by
and Putting σ = 1 in (29) and (30), we get the solution of the problem The solution represented by Eqs. (31) and (32) is similar to the exact solution in a closed form as follows: In Figs. 1 and 3, the actual and HPTM solutions of φ(χ, η) and ψ(χ, η) are calculated at σ = 1. In Figs. 2 and 4, the 3D graphs for φ(χ, η) and ψ(χ, η) for different fractional order show that the HPTM solutions derived are in a strong agreement with the actual    and the approximate solution. This comparison shows that the HPTM and the actual solutions are very close. As a result, the HPTM is a dependable new study that requires less computation of computations, is adaptable, and simple to use. In Tables 1 and 2  KdV system as follows: with the initial condition φ(χ, 0) = 3 -6 tanh 2 χ 2 , Using the Laplace transform on Eq. (34) by the application of initial condition given by Eq. (35), we get Applying the inverse Laplace transform, we get Using HPM on Eq. (37), we get Nonlinear steps given by He's polynomials H i (φ) and H i (ψ) are given as follows: Values of factors of He's polynomials are given as follows: Hence a series solution is given by Putting σ = 1 in (45), we get the solution of the problem φ(χ, η) = 3 -6 tanh 2 χ 2 + 6η sech 2 χ 2 tanh χ 2 (46) The solution given by Eq. (46) is similar to a closed form solution  Tables 3 and   4, the exact result and HPTM solution of different fractional order of φ(χ, η) and ψ(χ, η) at η = 0.01 are given.

Conclusions
In this paper, we calculated the fractional-order Kersten-Krasil'shchik coupled KdV-mKdV nonlinear system using the Laplace transform and the Atangana-Baleanu derivative. The suggested method is applied to obtain the solution of the given two problems.
The HPTM solution is in close contact with the exact result of the given problems. The present scenario also calculated the results of the given problems with fractional-order derivatives. The figures of the fractional-order results achieved have shown the convergence towards the results of integer order. Furthermore, the present method is simple, straightforward and required less computational cost; the current technique can be modified to solve other fractional-order partial differential equations.