Construction of an Approximate Analytical Solution for Multi-Dimensional Fractional Zakharov–Kuznetsov Equation via Aboodh Adomian Decomposition Method

: In this paper, the Aboodh transform is utilized to construct an approximate analytical solution for the time-fractional Zakharov–Kuznetsov equation (ZKE) via the Adomian decomposition method. In the context of a uniform magnetic ﬂux, this framework illustrates the action of weakly nonlinear ion acoustic waves in plasma carrying cold ions and hot isothermal electrons. Two compressive and rarefactive potentials (density fraction and obliqueness) are illustrated. With the aid of the Caputo derivative, the essential concepts of fractional derivatives are mentioned. A powerful research method, known as the Aboodh Adomian decomposition method, is employed to construct the solution of ZKEs with success. The Aboodh transform is a reﬁnement of the Laplace transform. This scheme also includes uniqueness and convergence analysis. The solution of the projected method is demonstrated in a series of Adomian components that converge to the actual solution of the assigned task. In addition, the ﬁndings of this procedure have established strong ties to the exact solutions to the problems under investigation. The reliability of the present procedure is demonstrated by illustrative examples. The present method is appealing, and the simplistic methodology indicates that it could be straightforwardly protracted to solve various nonlinear fractional-order partial differential equations.


Introduction
In recent years, fractional calculus has sparked a wave of interest, and it has been successfully tested and applied in a variety of real-world problems in science and technology [1][2][3][4][5][6][7][8]. Furthermore, it has been the subject of numerous investigations in many domains: for instance, signal processing, random walks, Levy statistics, chaos, porous media, electromagnetic flux, thermodynamics, circuits theory, optical fibre, and solid state physics. Moreover, a systematic attempt has been conducted to derive explicit solutions of partial differential equations (PDEs) [9][10][11][12][13].
The development of an integral transform to locate solutions in science can be connected back to P. S. Laplace's (1749-1827) work on statistical mechanics in the 1780s, in addition to J. B. Fourier's (1768-1830) treatise "La Théorie Analytique de La chaleur" (1822) reported in [14]. In 2013, K. S. Aboodh [15] introduced a new integral transform which is a modification of the Laplace transform. Aboodh transform (AT) is a valuable tool for solving certain DEs that the Sumudu transform cannot solve. Ever since, researchers have been particularly interested in the formation and acquisition of new integral transforms for numerous enhancements [16][17][18][19][20][21][22][23].
where F = F (x 1 , x 2 , t), D ρ t is the Caputo fractional derivative with order ρ, a 1 and b 1 are arbitrary constants and σ i , i = 1, 2, 3 are integers, and σ i = 0 (i = 1, 2, 3) shows the nature of nonlinear phenomena such as ion acoustic waves in the context of a symmetrical magnetic field in a plasma containing cold ions and hot isothermal electrons [39,40]. For example, in [38], the ZKEs were proposed to analyse a shallowly nonlinear isotope ripple in substantially magnetization impairment plasma in three dimensions. The approximate analytical solutions of fractional ZKEs are examined by the variation iteration method [41] and HPM [42], respectively. The detriment of many of the above mentioned strategies is that they are always hierarchical and require a lot of computational complexity. To mitigate computational cost and difficulty, we proposed a new approach called the Aboodh Adomian decomposition method (AADM), which is an amalgamation of the AT and the ADM for solving the time-fractional ZKE, which is the innovation of this research. The suggested technique generates a convergent series as a solution. AADM has fewer parameters than other analytical methods, and it is the preferred approach because it does not require discretion or linearization In this study, we first provide a fractional ZKE, followed by a description of the AADM, and then a uniqueness characterization of the AADM is presented. The convergence analysis is then explained in order to be applied to the ZK problem. We present an algorithm for AADM, discuss its estimation accuracy, and then show two examples that demonstrate the effectiveness and stability of a novel approach so that their obtained simulations can be analysed. Rarefaction curves are drawn for a graphical representation of variations in density fraction and obliqueness, which are associated with the derived results of electron superthermality. Finally, as a part of our concluding remarks, we discuss the accumulated facts of our findings.

Prelude
Several definitions and axiom outcomes from the literature are prerequisites in our analysis.
(2) Definition 2 ([15]). Aboodh transform (AT) for a function F (t) having exponential order over the set of functions is stated as where F (t) is represented by A F (t) = A(ω) and is described as Definition 3 ([43]). The inverse AT of a mapping F (t) is stated as Lemma 1. (Linearity property of AT) Let AT of F 1 (t) and F 2 (t) be P (ω) and Q(ω), respectively [44]: where γ 1 and γ 2 are arbitrary constants.

Lemma 2 ([45]
). The AT of Caputo fractional derivative of order ρ is stated as

Configuration for Aboodh Adomian Decomposition Method
In this note, we state the fundamental concept of AADM. The transform being utilized here is the refinement of the Laplace transform, and it is assumed for the time domain t ≥ 0. The AADM is addressed to the solution of the time-fractional KZE with the fractional time-derivative of the order ρ presented as follows: with the initial condition where D ρ = ∂ ρ ∂t ρ is the Caputo operator, while L and N are linear and nonlinear terms, respectively, andh(x 1 , x 2 , t) is the source term.
Employing the AT on (8) and utilizing the initial condition, we have The following infinite series demonstrates the AADM solution of F (x 1 , x 2 , t) as The following are Adomian polynomial forms for the nonlinear term in the given problem: where H  is represented as Substituting (12) and (13) in (11), we obtain In view of the linearity property of AT, we have Transforming the inverse AT into (15) yields

Qualitative Aspects of Aboodh-Adomian Decomposition Method
In what follows, we will demonstrate that the sufficient conditions assure the existence of a unique solution. Our desired existence of solutions in the case of AADM follows [46].
Proof. Assume that K = C[I], . represents all continuous mappings on the Banach space, defined on I = [0, T ] having the norm . . For this, we introduce a mapping Q : K → K, and we have and M F (x 1 , t) are also Lipschitzian with RF − R F <Ľ 1 F − F and LF − L F < L 2 F − F , whereĽ 1 andĽ 2 are Lipschitz constants, respectively, and F , F are various values of the mapping.
Under the assumption 0 < < 1, the mapping is contraction. Thus, by Banach contraction fixed point theorem, there exists a unique solution to (8). Hence, this completes the proof.

Theorem 2. (Convergence Analysis)
The general form solution of (8) will be convergent.
Firstly, we show that {S n } is a Cauchy sequence in Banach space in K. Taking into consideration a new representation of Adomian polynomials, we obtain Now, Consider n = q + 1; then, Analogously, from the triangular inequality, we have Hence, {S 1 } is a Cauchy sequence in K. As a result, the series ∞ ∑ n=0 F n is convergent, and this completes the proof.

Numerical Illustrations
Problem 1. Assume the following time-dependent fractional-order Zakharov-Kuznetsov equation [41,42]: subject to the initial condition where λ is an arbitrary constant.
Proof. Applying the AT on both sides of (21), we find Employing the inverse AT, we have It follows that Utilizing the Adomian decomposition method, we obtain where N (F ) is the He's polynomial describing a nonlinear term appearing in the abovementioned equations.
Accordingly, we can derive the remaining terms as follows: The approximate analytical AADM solution is The exact solution for ρ = 1 is presented by Table 1 and Table 2 demonstrates the exact AADM solution and the absolute error E abs = E exact − E approx for Problem 1. Figure 1 represents the comparison between the exact (left) and the approximate (right) solution, while Figure 2 describes the surface plot of the absolute error of the solution when ρ = 1, and λ = 0.001. Figure 3 represents a surface plot of approximate solutions for various fractional orders, ρ = 0.55, 0.67, 0.75, 0.85, 0.95, and 1. In addition, Figure 4 addresses approximate solutions for various fractional orders: ρ = 0.55, 0.67, 0.75, 0.77, 0.95, and 1 converge very rapidly to exact solutions, implying that approximate solutions are almost similar to exact solutions. As a result, the VIM [41] and HPM [42] demanded the evaluation of the Lagrangian multiplier, but the AADM demanded the evaluation of the Adomian polynomials, which entails less computation algebraic work. By obtaining further expressions of approximate solutions, the reliability of the analysis can be strengthened.
subject to the initial condition where λ is an arbitrary constant.
Proof. Applying the AT on both sides of (32), we find Employing the inverse AT, we have It follows that Utilizing the Adomian decomposition method, we obtain where N (F ) is the He's polynomial describing a nonlinear term appearing in the abovementioned equations.
Accordingly, we can derive the remaining terms as follows: The approximate analytical AADM solution is The exact solution for ρ = 1 is presented by Table 3 and Table 4 demonstrates the exact AADM solution and the absolute error E abs = E exact − E approx for Problem 2. Figure 5 represents the comparison between the exact (left) and the approximate (right) solution, while Figure 6 describes the surface plot of the absolute error of the solution when ρ = 1, and λ = 0.001. Figure 7 represents a surface plot of approximate solutions for various fractional orders ρ = 0.55, 0.67, 0.75, 0.85, 0.95, and 1. In addition, Figure 8 addresses approximate solutions for various fractional orders: ρ = 0.55, 0.67, 0.75, 0.77, 0.95, and 1 converge very rapidly to exact solutions, implying that approximate solutions are almost similar to exact solutions. As a comparison, the VIM [41] and HPM [42] necessitated the evaluation of the Lagrangian multiplier, but the AADM required the evaluation of the Adomian polynomials, which involved less algebraic computation. By obtaining further expressions of approximate solutions, the reliability of the analysis can be strengthened.

Other Aspects of ZKEs
Firstly, considering the fractional order to be 1 and rotating the coordinate axes (t, ζ) through an angle ϑ, maintaining the -axis stationary, in order to evaluate the temperature dependence of solitary waves in a direction making an angle ϑ with the t-axis, i.e., with the magnetization and lying in the (ζ − t) plane, the independent variables are adjusted in the following manner: Utilizing the aforesaid Scheme (41) in the ZKE (1), yields where Now, the steady state solution of the ZKE (42) in the form is investigated as follows: where Λ = x 1 − U t, whereas U is a constant velocity normalized to C. Employing (44) in (42), then, the steady state formulation is represented as Utilizing the suitable boundary assumptions, viz., (F 0 , F 0 and F 0 ) tends to 0 when Λ → ±∞, then, the solution of (45) is derived as where F m = 3U /η 1 denotes the peak amplitude, and L = 4η 2 /U is the width of solitons, respectively. Since the amplitude and width of ion acoustic waves in plasma are influenced by a variety of factors and physical parameters, it is fascinating to quantitatively determine their consequences on plasma carrying superthermality of cold and hot electrons. Figure 9a,b exhibited symmetric behaviour for positive and negative pressure structures with varied values of density fraction depending on the unperturbed cold electron to fluid ion concentration ratio, in order to see the influence of cold electron superthermality. It is remarkable that with fluctuations in the value of the superthermality of electrons, the wave profile is revealed to be dramatically altered by the superthermality of electrons.
The impact of obliqueness ϑ on both positive and negative potential is represented in Figure 10a,b. As a result, the increment in obliqueness ϑ strengthened the amplitude and width, respectively.

Conclusions
In this study, the AADM was proposed to investigate the time-fractional Zakharov-Kuznetsov equation regulating the nonlinear evolution of ion acoustic waves in a magnetised plasma having cold and hot temperature electrons. For the various physical characteristics, both positive (compressive) and negative (rarefactive) potential structures are generated that are symmetric with respect to origin. The methodology of the suggested technique has been considered to be more effective than other analytical schemes due to its confined number of estimations. The technique is clearly understood by the researchers because it involves implementing the AT explicitly to the projected problem and then adapt-ing the ADM. The inverse Aboodh transform is then employed to derive the approximate solution for the projected problem. To demonstrate the conformity of the developed model and precise solutions to the problems, we have shown 2D and 3D graphs, respectively. The findings acquired by the current report are in excellent accordance with the actual solution of Example 1 and 2 in the paper. Furthermore, the manuscript includes a graph of absolute errors and tabular results which have already been presented and addressed. This demonstrates that the proposed model provided adequate accuracy to the problem solution even though two terms of the series solution were considered. The simulation process reveals that the AADM has achieved an excellent agreement. It may be assumed that the AADM is extremely efficient and easy to implement in determining approximate analytical solutions of several fractional physical and biological models.