Novel Computations of the Time-Fractional Coupled Korteweg–de Vries Equations via Non-Singular Kernel Operators in Terms of the Natural Transform

: In the present research, we establish an effective method for determining the time-fractional coupled Korteweg–de Vries (KdV) equation’s approximate solution employing the fractional derivatives of Caputo–Fabrizio and Atangana–Baleanu. KdV models are crucial because they can accurately represent a variety of physical problems, including thin-ﬁlm ﬂows and waves on shallow water surfaces. Some theoretical physical features of quantum mechanics are also explained by the KdV model. Many investigations have been conducted on this precisely solvable model. Numerous academics have proposed new applications for the generation of acoustic waves in plasma from ions and crystal lattices. Adomian decomposition and natural transform decomposition techniques are combined in the natural decomposition method (NDM). We ﬁrst apply the natural transform to examine the fractional order and obtain a recurrence relation. Second, we use the Adomian decomposition approach to the recurrence relation, and then, using successive iterations and the initial conditions, we can establish the series solution. We note that the proposed fractional model is highly accurate and valid when using this technique. The numerical outcomes demonstrate that only a small number of terms are required to arrive at an approximation that is exact, efﬁcient, and trustworthy. Two examples are given to illustrate how the technique performs. Tables and 3D graphs display the best current numerical and analytical results. The suggested method provides a series form solution, which makes it quite easy to understand the behavior of the fractional models.


Introduction
Fractional calculus (FC) is a 17th-century invention that generalizes integer-order calculus to arbitrary-order calculus.The primary advantage of FC is that it describes a beneficial technique for researching memory and genetic characteristics in a wide range of phenomena.Furthermore, ordinary calculus is a subset of FC.The fundamental research of fractional derivatives has advanced rapidly in recent decades.Fractional calculus has received much attention in the last thirty years or more.Several academics have noticed that developing unique fractional derivatives (FDs) with distinct singular or nonsingular kernels is critical to address the demand for modeling real-world problems in various areas.Because most FDs do not have perfect solutions, approximations and numerical techniques must be used.More information on the definitions and properties of fractional derivatives can be checked at [1].Over the last five decades, research on the theory and application of differential equations (DEs) in terms of Caputo FD has been achieved [2,3].However, the Caputo FD has a unique kernel.Caputo and Fabrizio offered a solution to the solitary kernel problem using an exponential function in the last decade [4].Although this operator is local, it has several concerns.The corresponding integral is not fractional in the fractional order derivative.Atangana and Baleanu work hard to overcome local issues [5].They designed the Liouville-Caputo and Riemann-Liouville derivatives of the modified Mittag-Leffler function.Actually, this derivative is not only a differential operator; it may also be thought of as a filter regulator.This interesting derivative also has the benefit of explaining some of the macroscopic behavior of certain materials.Numerous researchers have been closely observing these derivatives' stimulating behaviors in recent years [6][7][8].
Non-singular FDs have been employed in several models.Several models have recently been studied for non-singular FDs [9,10].The authors investigated the dynamics of several physical issues using the Atangana-Baleanu (AB) and Caputo-Fabrizio (CF) FDs.Wang et al. used non-singular FD to develop a model for bank data with real field data from 2004 to 2014 [11].They illustrate the fractional ABC operator's better accuracy and adaptability and how it may be used with comfort to simulate such real-world events.Saifullah et al. used the non-singular FD in the CF sense to show the complicated progression of HIV infection [12].Khan [15].The works listed in the citations [16][17][18] are all helpful.The fractional calculus has many uses; see [19][20][21][22][23][24] for a few examples.
FC has been used to model physical and engineering processes best represented by differential fractional equations (FDEs).In recent decades, FDEs have been extensively used in numerous engineering and applied science fields.Nonlinear differential equations define most of the phenomena in nature.Thus, researchers pay more attention to the different branches of science and engineering to try to solve them.However, finding an exact solution is hard due to the involvement of nonlinear parts in these equations.The solutions of nonlinear FPDEs are of great concern in mathematics and useful applications [25][26][27].Consequently, knowing how to build a reliable technique to obtain the approximate or the exact solution of FPDEs is of great interest in the field of research of fractional models.Numerous analytical techniques have been used to find the solution of these problems.Such as iterative Laplace transform method [28], Laplace variational iteration method [29], approximate-analytical method [30], Laplace Adomian decomposition method [31], optimal homotopy asymptotic method [32], reduced differential transform method [33], homotopy analysis method [34], Natural transform decomposition method [35], Adomian decomposition method [36] and many more [37][38][39][40][41].
The Korteweg-de Vries (KdV) equation is a sort of partial differential equation that has recently been employed in defining several physical occurrences as an example of the formation and association of nonlinear waves.It was constructed as a modified equation controlling the movement of one-dimensional, large, small-amplitude surface gravity waves in a shallow water channel, as demonstrated in [42].The KdV equation is currently being investigated in various physical science fields, such as stratified internal waves, plasma physics, ion-acoustic waves, lattice dynamics, and collision-free hydromagnetic waves [43].A KdV model has been employed in quantum physics to explain several hypothetical physical phenomena.It is applied in fluid dynamics, aerodynamics, and continuum mechanics as a model for shock wave production, turbulence, solitons, mass transport, and boundary layer behavior [44].In this study, we derive an analytical solution to the nonlinear coupled time-fractional KdV equations that are provided by: subjected to initial sources ( where a, d < 0; b, c and e stand for constant parameters.Numerous academics have explored the KdV equation, which has applications in analyzing shallow-water waves and many other physical phenomena.The Adomian decomposition method, which yields accurate solutions in the form of a convergent series, is a well-known technique for solving homogeneous and nonhomogeneous, linear and nonlinear, homogeneous and nonlinear differential and partial differential equations, as well as integro differential and fractional differential equations.Without the need for linearization or disturbance, the Adomian decomposition method has been successfully and efficaciously applied to examine issues that have arisen in science and technology.On the other hand, the Adomian decomposition method requires a significant quantity of computer memory and more time for computational work.Therefore, it is inevitable that this method will be used with already-existing transform methods.Differential equations were solved by combining a number of transforms with additional methods.In [45,46], the natural decomposition method (NDM), a linked natural transform and Adomian decomposition approach, was developed for solving differential equations.It provides an approximation solution in series form.The central theme of this work is to solve time-fractional coupled KdV equations NTDM.Numerous studies have employed the NTDM to obtain approximate analytical solutions; it generated accurate and closely convergent outcomes.
This work is structured as follows.Section 2 defines and describes the natural transform's properties.Section 3 describes the overall implementation of the proposed technique.Section 4 covers the new technique and compares it to two different ways using two examples and presents tables and graphs to validate the NTDM.Section 5 contains the manuscript's conclusion.
(3) Definition 2. The non-integer derivative in Caputo sense is defined as [47] D Definition 3. The non-integer derivative in CF sense is stated as [47] D Definition 4. The non-integer derivative in the ABC sense is stated as [47] D and for ζ ∈ (0, ∞), the NT of F(ζ) is as where H(ζ) is the Heaviside function.
Definition 6.The inverse NT of F(ς, κ) is as having c 1 and c 2 constants.N[D Definition 8.The NT of D with M[β] denoting a normalization function.

The Proposed Scheme
This section focuses on an analytical approach for obtaining the solution of the differential equation having fractional-order as given below [49]: having the initial guess where L, N demonstrates the linear, non-linear functions, respectively, and h(ω, ζ) is an indicated source function.

Case I (NTDM CF )
In terms of NT and fractional CF derivative, Equation (1) transformed into with ). and By operating the inverse NT, we have The solution of F(ω, ζ) is expanded in series form as and N(F(ω, ζ)) is illustrated as with the Adomian polynomials A i as By switching Equations ( 23)-( 24) into ( 22), we have From (25), we obtain, By utilizing ( 26) into (23), we obtain the solution to Equation (1) in the NTDM CF sense as

Case II (NTDM ABC )
In terms of NT and fractional ABC derivative, Equation ( 1) is transformed into with By operating the inverse NT, we have After, we have From (25), we obtain Thus, we acquire the outcomes of (1) in terms of NTDM ABC as

Numerical Results
Example 1.Let us assume the fractional coupled KdV Equation ( 1) with a = −ς, b = −6ς, c = 2ν, d = −ψ, and e = 3ψ, having the initial guess In terms of NT, we obtain After, we obtain which simplifies to By operating the inverse NT, we have

NDM CF solution
The solution of F(ω, ζ) and G(ω, ζ) are expanded in series form as The nonlinear terms according to Adomian polynomials are as follows now by putting these terms in Equation (37), we obtain By equating both sides of Equation ( 39), we acquire Finally, we obtain the analytical solution of F(ω, ζ) and G(ω, ζ) as

NDM ABC solution
The solutions of F(ω, ζ) and G(ω, ζ) are expanded in series form as The nonlinear terms according to Adomian polynomials are now by putting these terms in Equation (37), we obtain By equating both sides of Equation ( 43), we acquire , Finally, we obtain the analytical solution of F(ω, ζ) and G(ω, ζ) as At β = 1, we obtain the exact solution as Example 2. Let us assume fractional coupled KdV Equation ( 1) with a = −1, b = −6, c = 3, d = −1, and e = 3, having initial guess In terms of NT, we obtain After, we obtain By operating the inverse NT, we have

NDM CF solution
The solution of F(ω, ζ) and G(ω, ζ) are expanded in series form as The nonlinear terms according to Adomian polynomials are now by putting these terms in Equation ( 50), we obtain By equating both sides of Equation ( 52), we acquire Finally, we obtain the analytical solution of F(ω, ζ) and G(ω, ζ) as

NDM ABC solution
The solutions of F(ω, ζ) and G(ω, ζ) are expanded in series form as The nonlinear terms according to Adomian polynomials are now by putting these terms in Equation ( 50), we obtain By equating both sides of Equation ( 56), we acquire , Finally, we obtain the analytical solution of F(ω, ζ) and G(ω, ζ) as At β = 1, we obtain the exact solution as

Results Discussion
The graphical and numerical analysis presented in this section offers valuable insights into the behavior and accuracy of our proposed solution method for the coupled Kortewegde Vries (KdV) equations using non-singular kernel operators in conjunction with the natural transform across varying values of the fractional parameter β.
In Figure 1, we depict the behavior of the exact solution of F(ω, ζ) alongside our approach's solutions at different values of β, including β = 1, β = 0.80, and β = 0.60 of F(ω, ζ) for Example 1.These graphs allow us to visually compare how well our method approximates the exact solution as β varies.Our approach evidently provides a reasonably accurate representation of the exact solution, with deviations becoming more noticeable as β decreases.
Figure 2 follows a similar pattern as Figure 1 but for G(ω, ζ).We observe the behavior of the exact solution and our approach's solutions at β = 1, β = 0.80, and β = 0.60 of G(ω, ζ).Again, these visualizations enable us to assess the performance of our method in approximating the exact solution.As β decreases, some deviation from the exact solution is observed, but our approach remains a promising approximation method.
We continue our analysis in Figure 3, but now for Example 2. We explore the behavior of the exact solution and our approach's solutions at different β values, including β = 1, β = 0.80, and β = 0.60 of F(ω, ζ) and G(ω, ζ).These graphs highlight the ability of our method to adapt to varying fractional parameters and provide reasonable approximations of the exact solution.Tables 1 and 2 present a quantitative analysis of the accuracy of our method by comparing the results of our approach with the exact solutions for F(ω, ζ) and G(ω, ζ) at various orders of β.The comparison between Tables 3 and 4 shows that the solutions obtained in this paper are more accurate than those obtained in [50].Table 5 presents a quantitative analysis of the accuracy of our method by comparing the results of our approach with the exact solutions for F(ω, ζ) and G(ω, ζ) at various orders of β.These tables offer numerical evidence of our method's performance and consistency across different orders of β.In summary, our graphical and numerical analysis demonstrates the effectiveness of our proposed method in approximating the solutions of coupled KdV equations via non-singular kernel operators within the framework of the natural transform.While some deviation from the exact solution is observed as β decreases, our method consistently provides reasonably accurate results, making it a valuable tool for solving these equations across various fractional parameter values.

Conclusions
Our study presents a practical methodology for approximating solutions to the timefractional coupled Korteweg-de Vries (KdV) equation, leveraging the power of fractional derivatives through Caputo-Fabrizio and Atangana-Baleanu formulations.We employed a natural decomposition method (NDM), which amalgamates the natural transform and Adomian decomposition methods.Initially, we utilized the natural transform to scrutinize the fractional order and derive a recurrence relation.Subsequently, the Adomian decomposition method was employed to process this recurrence relation.We derived a series solution by iteratively applying this approach and incorporating initial conditions.Our findings demonstrate that the proposed fractional model is highly accurate and robust when utilizing this method, and it remains valid even under extensive computational loads or limitations.To exemplify the technique's performance, we provided two illustrative examples, complemented by tables and 3D graphs, showcasing the excellence of our numerical and analytical results.This method furnishes a series-based solution that significantly enhances our comprehension of the behavior of fractional models, making it a valuable tool for analyzing and interpreting their dynamics.

having c 1
and c 2 constants.Definition 7. The NT of D β ζ F(ζ) in a Caputo manner is given by[47]
et al. used non-singular FDs to investigate the nonlinear Schrodinger equation's wave propagation [13].Rahman et al. examined the Φ 4 −equation with nonsingular FDs [14].Khan et al. examined the KdV-Burger equation's wave dynamics assuming non-singular FDs

Table 1 .
Comparison between the proposed method and exact solutions for F(ω, ζ) at numerous orders of β of Example 1. Solution at β = 0.8 (NT DM ABC ) at β = 1 (NT DM CF ) at β = 1 Exact Solution

Table 2 .
Comparison between the exact solution and our solution for G(ω, ζ) at numerous orders of β of Example 1.

Table 3 .
[50]r comparison between our solution for F(ω, ζ) and the results obtained in[50]of Example 1.

Table 4 .
[50]r comparison between our solution for G(ω, ζ) and the results obtained in[50]of Example 1.

Table 5 .
Comparison between the exact solution and our solution for F(ω, ζ) and G(ω, ζ) at numerous orders of β of Example 2.