Fractional Numerical Simulation of Coupled Approximate Long Wave and Modiﬁed Boussinesq Equations Involving Mittag-Lefﬂer Kernel

: This study examines approximate long wave and the modiﬁed Boussinesq equations, as well as their complexities with the Atangana–Baleanu fractional derivative operator in the Caputo sense. The analytical solution of the aforementioned model is discussed using the Elzaki transform and the Adomian decomposition method. These problems are indispensable for deﬁning the characteristics of surface water waves by applying a particular relationship of dispersion. We used Elzaki transformation on time-fractional approximate long wave and modiﬁed Boussinesq equations, followed by inverse Elzaki transformation, to achieve the results of the equations. To validate the methodology, we concentrated on two systems and compared them to the actual solutions. The numerical and graphical results demonstrate that the proposed method is computationally precise and straightforward for investigating and resolving fractionally coupled nonlinear phenomena that occur in scientiﬁc and technological.


Introduction
In recent years, fractional partial differential equations (PDEs) have achieved popularity and validation due to their demonstrated applicability in a broad range of applied science fields of study. Evaluating that the nonlinear oscillation of fractional derivatives can be used to model earthquakes, the fractional fluid-dynamic traffic model can be used to compensate for the shortcoming caused by the presumption of continued traffic flow. Numerous researchers, such as Coimbra, Riemann-Liouville, Riesz, Weyl, Hadamard, Caputo-Fabrizio, Atangana-Baleanu, Jumarie, Grunwald-Letnikov, Liouville, and Caputo, have proposed various fractional operator formulations and concepts. In addition, fractional partial differential equations are used to model numerous physical phenomena, such as population dynamics, chemical reaction, image processing, virology, thermodynamics, bifurcation, porous media, Levy statistics, physics, and engineering problems [1][2][3][4][5][6]. Symmetry analysis is beautiful, especially when it comes to the study of partial differential equations and, more specifically, those equations that come from the Mathematics of Finance. The key to comprehending nature is symmetry; however, the majority of natural phenomena lack symmetry. The technique of spontaneous symmetry-breaking is a sophisticated method for concealing symmetry. There are two types of symmetries: infinitesimal and finite. On the other hand, space is in a constant state of flux. Mathematicians have always been intrigued by patterns [7,8]. Classifications of both planar and spatial patterns became more serious in the nineteenth century. In physics and applied mathematics, the necessity of determining novel approximate or exact solutions for nonlinear differential equations remain a fundamental problem requiring the development of new methodologies [9][10][11][12][13].
In the 20th century, Whitham [14], Broer [15], and Kaup [16] investigated the Whitham-Broer-Kaup (WBK) equations, which explain the propagation of shallow water with a different dispersion comparison. Consider the system of fractional-order Whitham-Broer-Kaup equations [17]: where ∂ β ∂τ β is the fractional derivative, u = u(ψ, τ) is the horizontal velocity and v = v(ψ, τ) is the height that is different from the liquid's equilibrium position. In addition, since a and b are constants representing various diffusion strengths, Equation (1) becomes a modified Boussinesq equation when a = 1 and b = 0. Similarly, for a = 0 and b = 1, the scheme represents the standard long wave model. These relations are derived from fluid mechanics to demonstrate the propagation of waves in dissipative and nonlinear media. They are suggested for models via the water leakage through porous groundwater strata and are broadly used in coastal and marine ecosystems technology [18][19][20]. In addition, Equation (1) is the basis for different designs used to represent the unconfined groundwater, such as drains and surface water fluid models.
The Adomian decomposition method was developed by the American scientist G. Adomian. It emphasizes the search for a set of solutions and the decomposition of the nonlinear operator into a sequence in which Adomian polynomials are repeatedly computed using the terms. This method is enhanced by the Elzaki transformation, and the enhanced method is called the Elzaki decomposition method [21][22][23][24]. Elzaki transform is a contemporary integral transform that Tarig Elzaki introduced in 2010. Elzaki transformation is an altered combination of the Sumudu and Laplace transforms [25][26][27].

Theorem 1.
Elzaki transformation convolution theorem, The following equality holds: where the Elzaki transform is indicated by E {.}.
Proof. From Definition 2, we have: Then, taking into account the definition and convolution of the Elzaki transform, we get

Methodology
Consider the general fractional partial differential equation is given as with initial condition here L and M are linear and non-linear terms and source function is h(ψ, τ).
Applying the Elzaki transform of Equation (11), we get Using inverse Elzaki transform (13), we have The nonlinear term M(ν(ψ, τ)) is defined as where A i represents the Adomian polynomials and they are expressed as By putting Equations (15) and (16) into (14), we get From (8), we get we achieved the EDM result of (11) by put (18) into (16),

Numerical Examples
Example 1. Consider the fractional modified Boussinesq equations with the initial conditions Using the Elzaki transformation to Equation (20), we get The nonlinear term can be expressed as By making it easier, the above equation can be written as We obtain by using inverse Elzaki transform to Equation (24), we get Now we are implementing EDM Suppose that the unknown terms u(ψ, τ) and v(ψ, τ) obtain the infinite series solution is given as The Adomian polynomials represented the nonlinear functions and are expressed by On both sides, comparing Equation (27) The remaining terms of the EDM solution u ı and v ı for (ı ≥ 3) can be achieved simply. The series type solution is determined as For Equation (20), the exact solutions are achieved at β = 1 In Figure 1, the actual and analytical results for u(ψ, τ) at β = 1 of Example 1 with a show that they closed contact with each other. Figure 2, analytic result figure for u(ψ, τ) at β = 0.8 and 0.6 and Figure 3, analytic result figure for u(ψ, τ) at different value of β of Example 1. Figure 4, the actual and analytic result of v(ψ, τ) at β = 1 of Example 1 with a show that they closed contact with each other and Figure 5 with the initial conditions Using the Elzaki transformation to Equation (33), we get The nonlinear term can be written as By making it easier, the above equation can be written as (37) Applying inverse Elzaki transform to Equation (37), we get Now we apply EDM Suppose that the unknown terms u(ψ, τ) and v(ψ, τ) obtain the infinite series solution as The nonlinear functions expressed with the help of Adomian polynomials uu ψ = ∑ ∞ ı=0 A ı , uv ψ = ∑ ∞ ı=0 B ı , and vu ψ = ∑ ∞ ı=0 C ı . We can be written the Equation (38) as a solution for using these functions, On both sides, comparing the Equation (40), we get .
The remaining terms of the EDM solution u ı and v ı (ı ≥ 3) can be obtained simply. The series type solution is determined as For Equation (33), the exact solutions is achieved at β = 1 In Figure 6, the actual and analytical results for u(ψ, τ) at β = 1 of Example 2 with a show that they closed contact with each other. Figure 7, analytic result figure for u(ψ, τ) at β = 0.8 and 0.6 and Figure 8, analytic result figure for u(ψ, τ) at different value of β of Example 2. Figure 9, the actual and analytic result of v(ψ, τ) at β = 1 of Example 2 with show that the closed contact with each other and Figure 10, the analytic result graph for v(ψ, τ) at different value of β of Example 2. In Tables 3 and 4, the exact and EDM solutions of u(ψ, τ) and v(ψ, τ) of Example 2 at different fractional-order of β having various values of ψ and τ.

Conclusions
In this study, the Elzaki decomposition method is applied to analyze the system of approximate long wave and modified Boussinesq equations. The two problems are investigated to calculate and validate the accuracy of the proposed method. Compared to the current analytic method for finding approximate results in the system of fractional nonlinear partial differential equations, the proposed methodology is quick and easy. The solutions have been demonstrated in the shape of tables and graphs. The suggested technique makes a list of results in the form of a given equation, which is more accurate and requires less work. For both fractionally coupled systems, calculations were done to figure out the absolute error. At β = 1, a number of mathematical results are evaluated by comparing what is known as the analytic method and the actual solutions. The current methods are better because they require less math and are more accurate. Moreover, it is shown that the presented method is very simple to use, it works, and can be used to solve a different fractional system of differential equations. In the future, we plan to extend the application of the Elzaki decomposition method to analyze higher dimensional physical applications. In particular, we will consider the recent implementations of fractionalorder problems.

Data Availability Statement:
The numerical data used to support the findings of this study are included within the article.