A detailed study on a new (2+1)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$(2 + 1)$\end{document}-dimensional mKdV equation involving the Caputo–Fabrizio time-fractional derivative

The present article aims to present a comprehensive study on a nonlinear time-fractional model involving the Caputo–Fabrizio (CF) derivative. More explicitly, a new (2+1)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$(2 + 1)$\end{document}-dimensional mKdV (2D-mKdV) equation involving the Caputo–Fabrizio time-fractional derivative is considered and an analytic approximation for it is retrieved through a systematic technique, called the homotopy analysis transform (HAT) method. Furthermore, after proving the Lipschitz condition for the kernel ψ(x,y,t;u)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\psi (x,y, t;u)$\end{document}, the fixed-point theorem is formally utilized to demonstrate the existence and uniqueness of the solution of the new 2D-mKdV equation involving the CF time-fractional derivative. A detailed study finally is carried out to examine the effect of the Caputo–Fabrizio operator on the dynamics of the obtained analytic approximation.


Introduction
The classical KdV equation is a nonlinear partial differential equation to model waves on shallow water surfaces that was established by Korteweg and de Vries in 1895. This exactly solvable model has been the topic of many research works. Nowadays, unique applications of the classical KdV equation have been suggested by many scholars as it can be used to describe long internal waves in a density-stratified ocean, ion-acoustic waves in a plasma, and acoustic waves on a crystal lattice. The mathematical form of the classical KdV equation is given by [1] u t + α 1 u xxx + α 2 uu x = 0.
There are different variations of the classical KdV equation, some of those reported in the literature are [2,3]: u t + α 1 u xxx + α 2 u 2 u x = 0, (modified KdV), u t + α 1 u xxx + α 2 u n u x = 0, (generalized KdV), and the (2 + 1)-dimensional KdV equation Many efforts have been devoted to studying the KdV-type equations, for example, Wazwaz [2] used different reliable methods to obtain solitons and periodic solutions of the KdV, modified KdV, and generalized KdV equations, and Wang [3] derived lump solutions of the (2 + 1)-dimensional KdV equation by means of an ansatz based on the quadratic functions. Recently, the fractional forms of the KdV-type equations have been explored using a series of systematic methods in [4,5].
The classical form of the new 2D-mKdV Eq. (1) was first proposed by Wang and Kara [12] in 2019 using the extended Lax pair. Wang and Kara in [12] extracted a group of solitary wave solutions of the new 2D-mKdV equation (its classical form) through the Lie symmetry method. Recently, nonlinear ODEs/PDEs involving the Caputo-Fabrizio fractional derivative have gained significant attention owing to their potential to describe many complicated physical phenomena. In this respect, Shah et al. [13] analyzed a nonlinear model of dengue fever disease with the CF fractional derivative using the Laplace Adomian decomposition method. Owolabi and Atangana [14] explored a series of nonlinear fractional parabolic differential equations involving the CF derivative thought a numerical scheme. In another work performed by Arshad et al. [15], the CD4 + T-cells model of HIV infection with the CF fractional derivative was studied using an effective numerical scheme. Shaikh et al. [16] employed the iterative Laplace transform method to analyze a group of fractional reaction-diffusion equations involving the CF derivative. More articles can be found in .
The rest of this paper is as follows: In Sect. 2, the Caputo-Fabrizio fractional operators are reviewed in detail. In Sect. 3, the Lipschitz condition for the kernel ψ(x, y, t; u) is proved, then the fixed-point theorem is formally applied to show the existence and uniqueness of the solution of the new 2D-mKdV equation involving the CF time-fractional derivative. In Sect. 4, an analytic approximation for the new 2D-mKdV equation involving the Caputo-Fabrizio time-fractional derivative is acquired using the HAT method. The results of this paper are summarized in the last section.

Basic definitions and features
This section presents the basic definitions and features of the Caputo-Fabrizio fractional operators. In this respect, first the Caputo-Fabrizio fractional derivative and integral are defined, then the Laplace transform of the Caputo-Fabrizio fractional derivative is given. H 1 (a, b), b > a and α ∈ (0, 1]. Then, the Caputo-Fabrizio fractional derivative of u(t) of order α is given by [17] CF where M(α) is a normalized function satisfying M(0) = M(1) = 1.

Definition 2
The Caputo-Fabrizio fractional integral of u(t) of order α (α ∈ (0, 1]) is given by [18] CF and in the general case, we have

Theorem 1 The following Lipschitz condition holds for the Caputo-Fabrizio fractional derivative given in Definition
Proof In a similar manner as in [23], it is easy to show that

The model and the existence and uniqueness of its solution
To start, let us consider This suggests that Eq. (1) can be rewritten as Applying the CF fractional integral to both sides of Eq. (2) results in In order to show that the kernel ψ(x, y, t; u) satisfies the Lipschitz condition, we first consider bounded functions u(x, y, t) and v(x, y, t). Using the triangle property of norms, one can find This confirms that the Lipschitz condition is satisfied for the kernel ψ(x, y, t; u). Now, based on the Eq. (3) and the fixed-point theorem, an iterative scheme is established as follows:

Theorem 2 If the function u(x, y, t) is bounded, then
Proof by induction Suppose that n = 1. Then, one can write e 1 (x, y, t) = u 1 (x, y, t)u 0 (x, y, t) u(x, y, 0) . Now, if the relation holds for n -1, namely then, we will prove that To show this, we proceed as follows:

Theorem 3 If at t = t 0 we have
For t = t 0 , one obtains λt 0 < 1, the above series is convergent, and therefore, u n (x, y, t) exists and is bounded for any n.
Besides, by assuming one can easily prove that Proof by contradiction To start, let us consider two solutions, say u(x, y, t) and v(x, y, t), for the model (1). One can write Consequently,

The new CF time-fractional 2D-mKdV equation and its analytical solutions
In the present section, first soliton solutions of the classical form of the model are extracted using an ansatz method, then the HAT method is used to acquire an analytic approximation for the new CF time-fractional 2D-mKdV equation.

Soliton solutions of the classical form of the model
To find soliton solutions, a test function is considered as follows: u(x, y, t) = A 0 + A 1 sech(α 1 x + α 2 y + α 3 t).

The model and its analytic approximation
To obtain an analytic approximation, we apply the Laplace transform to both sides of Eq. (1). Such an operation results in -6u 2 (x, y, t)u y (x, y, t) + u xxx (x, y, t)u yyy (x, y, t) -3u xxy (x, y, t) + 3u xyy (x, y, t) = 0.

Therefore, the series solution of the new CF time-fractional 2D-mKdV equation is derived as
u(x, y, t) = -4 e x-y 1 + e 2(x-y) -32 It is noteworthy that for α = 1 and h = -1, the above series converges to the following exact solution: u(x, y, t) = -4 e x-y+8t 1 + e 2(x-y+8t) . Figure 1 presents the 3-order approximation of the new CF time-fractional 2D-mKdV equation for α = 1, 0.99, and 0.98 against the exact solution. From this figure, a full agreement between the 3-order approximation (when α = 1) and the exact solution is obviously observed. The absolute error of the 3-order approximation (when α = 1) and the exact solution has been presented in Table 1. The results confirm the efficiency of the HAT method   in deriving an analytic approximation with high accuracy. Finally, Fig. 2 shows the 3D plots of the exact solution and the 3-order approximation for α = 1, 0.99, and 0.98. It is believed that the analytic approximation given by the HAT method can precisely predict the dynamics of the dark soliton solution of the new 2D-mKdV equation involving the CF time-fractional derivative.

Conclusion
A thorough study on a nonlinear model involving the Caputo-Fabrizio time-fractional derivative was carried out successfully in the current paper. In this respect, a new 2D-mKdV equation designed with the CF time-fractional derivative was considered, and an analytic approximation for it was formally derived using a systematic approach, named the homotopy analysis transform method. The existence and uniqueness of the solution of the new 2D-mKdV equation involving the Caputo-Fabrizio time-fractional derivative were studied by proving the Lipschitz condition for the kernel ψ(x, y, t; u) and applying the fixed-point theorem. A detailed study was finally accomplished to investigate the effect of the Caputo-Fabrizio operator on the dynamics of the obtained analytic approximation. The results presented herein confirm the efficiency of the HAT method in deriving an analytic approximation with high accuracy for nonlinear models involving the Caputo-Fabrizio time-fractional derivative. Our future work is to obtain an analytic approximation of the new 2D-mKdV equation with the Atangana-Baleanu time-fractional derivative and study the effect of the Atangana-Baleanu operator on the dynamics of the approximate solution.