Research article

New computations for the two-mode version of the fractional Zakharov-Kuznetsov model in plasma fluid by means of the Shehu decomposition method

  • Received: 06 August 2021 Accepted: 24 October 2021 Published: 05 November 2021
  • MSC : 26A33, 26A51, 26D07, 26D10, 26D15

  • In this research, the Shehu transform is coupled with the Adomian decomposition method for obtaining the exact-approximate solution of the plasma fluid physical model, known as the Zakharov-Kuznetsov equation (briefly, ZKE) having a fractional order in the Caputo sense. The Laplace and Sumudu transforms have been refined into the Shehu transform. The action of weakly nonlinear ion acoustic waves in a plasma carrying cold ions and hot isothermal electrons is investigated in this study. Important fractional derivative notions are discussed in the context of Caputo. The Shehu decomposition method (SDM), a robust research methodology, is effectively implemented to generate the solution for the ZKEs. A series of Adomian components converge to the exact solution of the assigned task, demonstrating the solution of the suggested technique. Furthermore, the outcomes of this technique have generated important associations with the precise solutions to the problems being researched. Illustrative examples highlight the validity of the current process. The usefulness of the technique is reinforced via graphical and tabular illustrations as well as statistics theory.

    Citation: Maysaa Al-Qurashi, Saima Rashid, Fahd Jarad, Madeeha Tahir, Abdullah M. Alsharif. New computations for the two-mode version of the fractional Zakharov-Kuznetsov model in plasma fluid by means of the Shehu decomposition method[J]. AIMS Mathematics, 2022, 7(2): 2044-2060. doi: 10.3934/math.2022117

    Related Papers:

    [1] Samia Bushnaq, Sajjad Ali, Kamal Shah, Muhammad Arif . Approximate solutions to nonlinear fractional order partial differential equations arising in ion-acoustic waves. AIMS Mathematics, 2019, 4(3): 721-739. doi: 10.3934/math.2019.3.721
    [2] M. Ali Akbar, Norhashidah Hj. Mohd. Ali, M. Tarikul Islam . Multiple closed form solutions to some fractional order nonlinear evolution equations in physics and plasma physics. AIMS Mathematics, 2019, 4(3): 397-411. doi: 10.3934/math.2019.3.397
    [3] M. Hafiz Uddin, M. Ali Akbar, Md. Ashrafuzzaman Khan, Md. Abdul Haque . New exact solitary wave solutions to the space-time fractional differential equations with conformable derivative. AIMS Mathematics, 2019, 4(2): 199-214. doi: 10.3934/math.2019.2.199
    [4] Mamta Kapoor, Nehad Ali Shah, Wajaree Weera . Analytical solution of time-fractional Schrödinger equations via Shehu Adomian Decomposition Method. AIMS Mathematics, 2022, 7(10): 19562-19596. doi: 10.3934/math.20221074
    [5] Rasool Shah, Abd-Allah Hyder, Naveed Iqbal, Thongchai Botmart . Fractional view evaluation system of Schrödinger-KdV equation by a comparative analysis. AIMS Mathematics, 2022, 7(11): 19846-19864. doi: 10.3934/math.20221087
    [6] Azzh Saad Alshehry, Humaira Yasmin, Rasool Shah, Roman Ullah, Asfandyar Khan . Numerical simulation and analysis of fractional-order Phi-Four equation. AIMS Mathematics, 2023, 8(11): 27175-27199. doi: 10.3934/math.20231390
    [7] Rehana Ashraf, Saima Rashid, Fahd Jarad, Ali Althobaiti . Numerical solutions of fuzzy equal width models via generalized fuzzy fractional derivative operators. AIMS Mathematics, 2022, 7(2): 2695-2728. doi: 10.3934/math.2022152
    [8] Yunmei Zhao, Yinghui He, Huizhang Yang . The two variable (φ/φ, 1/φ)-expansion method for solving the time-fractional partial differential equations. AIMS Mathematics, 2020, 5(5): 4121-4135. doi: 10.3934/math.2020264
    [9] Aslı Alkan, Halil Anaç . The novel numerical solutions for time-fractional Fornberg-Whitham equation by using fractional natural transform decomposition method. AIMS Mathematics, 2024, 9(9): 25333-25359. doi: 10.3934/math.20241237
    [10] Naher Mohammed A. Alsafri, Hamad Zogan . Probing the diversity of kink solitons in nonlinear generalised Zakharov-Kuznetsov-Benjamin-Bona-Mahony dynamical model. AIMS Mathematics, 2024, 9(12): 34886-34905. doi: 10.3934/math.20241661
  • In this research, the Shehu transform is coupled with the Adomian decomposition method for obtaining the exact-approximate solution of the plasma fluid physical model, known as the Zakharov-Kuznetsov equation (briefly, ZKE) having a fractional order in the Caputo sense. The Laplace and Sumudu transforms have been refined into the Shehu transform. The action of weakly nonlinear ion acoustic waves in a plasma carrying cold ions and hot isothermal electrons is investigated in this study. Important fractional derivative notions are discussed in the context of Caputo. The Shehu decomposition method (SDM), a robust research methodology, is effectively implemented to generate the solution for the ZKEs. A series of Adomian components converge to the exact solution of the assigned task, demonstrating the solution of the suggested technique. Furthermore, the outcomes of this technique have generated important associations with the precise solutions to the problems being researched. Illustrative examples highlight the validity of the current process. The usefulness of the technique is reinforced via graphical and tabular illustrations as well as statistics theory.



    Fractional partial differential equations (PDEs) have gained prominence and recognition in recent years, owing to their verified applicability in a wide range of relatively diverse domains of science and engineering. For instance, considering the nonlinear oscillation of fractional derivatives can be employed to model earthquakes, fractional derivatives in a fluid-dynamic traffic model can be leveraged to alleviate the deficiency caused by the assumption of a continuous flow of traffic. Researchers including Coimbra, Davison and Essex, Riesz, Riemann and Liouville, Hadamard, Weyl, Jumarie, Caputo and Fabrizio, Atangana and Baleanu, Grünwald and Letnikov, Liouville and Caputo have proposed a variety of fractional operator formulations and conceptions. On the other hand, the Liouville-Caputo is the finest fractional filter. Furthermore, fractional PDEs are used to model a variety of physical phenomena, including chemical reaction and population dynamics, virology, image processing, bifurcation, thermodynamics, Levy statistics, porous media, physics, and engineering problems, (see [1-7]).

    The Shehu transform (ST) was recently highlighted by Maitama and Zhao [8] as an interesting integral transformation. A modification of the Laplace and Sumudu transformation is the ST. However, we can retrieve the Laplace transform by replacing ϖ=1 in ST. This approach can be used to compress complex non-linear PDEs into simpler equations.

    The Shehu transform (ST) was recently highlighted by Maitama and Zhao [8] as an interesting integral transformation. A modification of the Laplace and Sumudu transformation is the ST. However, we can retrieve the Laplace transform by replacing ϖ=1 in ST. This approach can be used to compress complex non-linear PDEs into simpler equations.

    The comprehensive evaluation of numerous advanced asymptotic approaches for the exploration of solitary solutions of nonlinear PDEs, and DEs has been presented, see [9-16]. For instance, the Adomian decomposition method (ADM) [17] for obtaining seven order Sawada-Kotera equations, pseudo-spectral method (PSM) [18] for finding the numerical solution of the Laxs 7th-order KdV equation, q-homotopy analysis method (q-HAM) for finding the convergence of special PDES [19], Lie symmetry analysis (LSA) [20] for dealing with the conservation laws and exact solutions of the seventh-order time fractional Sawada-Kotera-Ito equation, Laguerre wavelets collocation method (LWCM) and Haar wavelet for the numerical solution of the Benjamina-Bona-Mohany equations [21], a new Legendre Wavelets decomposition method (NLWDM) for solving PDEs [22], discrete Adomian decomposition method (DADM) [23] for constructing numerical solution of time fractional Navier-Stokes equation.

    The ZKE was built in two dimensions to demonstrate nonlinear complex phenomena such as isotope waves in a massively magnetic flux uncompressed plasma [24,25].The SDM will be implemented to develop the major objectives of this research. The time-fractional ZKE is stated as:

    DδζU+a1(Uη1)ϕ+b1(Uη2)ϕϕϕ+b1(Uη3)ζζϕ=0, (1.1)

    where U=U(ϕ,ψ,ζ) and Dδζ is the Caputo fractional derivative with order δ,0<δ1, a1 and b1 are arbitrary constants and ηi,i=1,2,3 are integers and ηi0(i=1,2,3) that demonstrates the characteristics of physical phenomena such as ion acoustic waves in a plasma consisting of cold ions and hot isothermal electrons in the framework of a balanced magnetic flux ([26,27]). In [25], for example, the ZKEs were used to investigate shallowly nonlinear isotope ripples in significantly magnetism impaired plasma in three dimensions.

    In spite of the incredible improvement, the Adomian decomposition method (ADM) was contemplated by Gorge Adomian in 1980. The ADM, for example, has been effectively defined in numerous analytical structures of PDEs, especially in Burger's equation [28], time-fractional Kawahara equation [29], fuzzy heat-like and wave-like equations [30] and Lane-Emden-Fowler type equations [31]. The ADM was found to be strongly associated with a plethora of integral transforms, including ARA, Shehu, Fourier, Aboodh, Laplace and others. Presently, modified Laplace ADM [32] has been utilized to effectively resolve Volterra integral equations employing the noted numerical formulation, discrete ADM [23] has been used to solve the time-fractional Navier-Stokes equation, and Laplace ADM [33] has been considered to identify the approximate results of a fractional system of epidemic structures of a vector-borne disease, and so on.

    Several of the aforesaid approaches have the disadvantage of being always stratified and necessitating a significant amount of algorithmic effort. To minimize the computing complexity and intricacy, we suggested the Shehu decomposition method (SDM), which is a composition of the ST and the decomposition method for solving the time-fractional ZKE, which is the main motivation for this research. The projected technique develops a convergent series as a solution. SDM has fewer parameters than other analytical methods. It is the preferred approach because it does not require discretion or linearization.

    In this study, we first provided a fractional ZKE, followed by a description of the SDM, and then, a comparison characterization of the SDM presented with the existing methods. The graphical representations were then thoroughly explained in relation to the ZK problem. We presented an algorithm for SDM, discussed its estimation accuracy, and then showed two examples that demonstrate the effectiveness and stability of a novel approach so that their obtained simulations can be analyzed. Finally, as a part of our concluding remarks, we discussed the accumulated facts of our findings.

    In order to perform our research, we require various terminologies and postulate outcomes from the literature.

    Definition 2.1. ([8]) Shehu transform (ST) for a mapping U(ζ) containing exponential order defined on the set of mappings is described as follows:

    S={U(ζ)|P,p1,p2>0,|U(ζ)|<Pexp(|ζ|/pȷ),ifζ(1)ȷ×[0,),ȷ=1,2;(P,p1,p2>0)}, (2.1)

    where U(ζ) is denoted by S[U(ζ)]=S(ξ,ϖ), is stated as

    S[U(ζ)]=0U(ζ)exp(ξϖζ)dζ=S(ξ,ϖ),ζ0,ϖ[κ1,κ2]. (2.2)

    The following is an example of a supportive ST:

    S[ζδ]=0exp(ξϖζ)ζδdζ=Γ(δ+1)(ϖξ)δ+1. (2.3)

    Definition 2.2. ([8]) The inverse ST of a function U(ζ) is described as

    S1[(ϖξ)mδ+1]=ζmδΓ(mδ+1),(δ)>0,andm>0. (2.4)

    Lemma 2.3. Consider ST of U1(ζ) and U2(ζ) are P(ξ,ϖ) and Q(ξ,ϖ), respectively [8],

    S[γ1U1(ζ)+γ2U2(ζ)]=S[γ1U1(ζ)]+S[γ2U2(ζ)]=γ1P(ξ,ϖ)+γ2Q(ξ,ϖ), (2.5)

    where γ1 and γ2 are unspecified terms.

    Lemma 2.4. ([8])For order δ>0, the Caputo fractional derivative (CFD) of ST is defined as

    S[DδζU(ζ)]=(ξϖ)δS[U(ϕ,ζ)]m1κ=0(ξϖ)δκ1U(κ)(ϕ,0),m1δm,mN. (2.6)

    Considering the nonlinear partial differential equation:

    DδζU(ϕ,ζ)+LU(ϕ,ζ)+ˇNU(ϕ,ζ)=F(ϕ,ζ),ζ>0,0<δ1, (3.1)

    subject to the condition

    U(ϕ,0)=G(ϕ), (3.2)

    where Dδζ=δU(ϕ,ζ)ζδ indicates the CFD with 0<δ1 while L and ˇN are linear/nonlinear factors and the source term refers to F(ϕ,ζ).

    Implementing the ST to (3.1), and we attain

    S[DδζU(ϕ,ζ)+LU(ϕ,ζ)+ˇNU(ϕ,ζ)]=S[F(ϕ,ζ)].

    Applying the differentiation property of ST, yields

    ξδϖδU(ξ,ϖ)=m1κ=0(ξϖ)δκ1U(κ)(0)+S[LU(ϕ,ζ)+ˇNU(ϕ,ζ)]+S[F(ϕ,ζ)]. (3.3)

    Th inverse ST of (3.3) provides

    U(ϕ,ζ)=S1[m1κ=0(ξϖ)δκ1U(κ)(0)+ϖδξδS[F(ϕ,ζ)]]S1[ϖδξδS[LU(ϕ,ζ)+ˇNU(ϕ,ζ)]]. (3.4)

    The infinite series representation of SDM is denoted by the mapping U(ϕ,ζ) as follows:

    U(ϕ,ζ)=m=0Um(ϕ,ζ). (3.5)

    Thus, the nonlinearity ˇN(ϕ,ζ) can be estimated by the ADM represented as

    ˇNU(ϕ,ζ)=m=0˜Am(U0,U1,...),m=0,1,..., (3.6)

    where

    ˜Am(U0,U1,...)=1m![dmdλmˇN(ȷ=0λȷUȷ)]λ=0,m>0.

    Substituting (3.5) and (3.6) into (3.4), we have

    m=0Um(ϕ,ζ)=G(ϕ)+˜G(ϕ)S1[ϖδξδS[LU(ϕ,ζ)+m=0˜Am]]. (3.7)

    Consequently, the following is the recursive methodology for (3.7):

    U0(ϕ,ζ)=G(ϕ)+˜G(ϕ),m=0,Um+1(ϕ,ζ)=S1[ϖδξδS[L(Um(ϕ,ζ))+m=0˜Am]],m1. (3.8)

    Example 4.1. Assume the following time-dependent fractional-order Zakharov-Kuznetsov equation:

    DδζU+U2ϕ+18[ϕ(2U2ψ2)+3U2ϕ3]=0, (4.1)

    subject to the initial condition

    U(ϕ,ψ,0)=43θsinh2(ϕ+ψ), (4.2)

    where θ is an arbitrary constant.

    Proof. Applying the ST on both sides of (4.1), we find

    S[δUζδ]=S[U2ϕ+18[ϕ(2U2ψ2)+3U2ϕ3]],(ξϖ)δS[U(ϕ,ψ,ζ)]n11κ=0(ξϖ)δκ1κU(κ)(ϕ,0)ζκ=S[U2ϕ+18[ϕ(2U2ψ2)+3U2ϕ3]]. (4.3)

    Employing the inverse ST, we have

    U(ϕ,ψ,ζ)=S1[(ϖξ)δn11κ=0(ξϖ)δκ1κU(κ)(ϕ,0)ζκ(ϖξ)δS[U2ϕ+18[ϕ(2U2ψ2)+3U2ϕ3]]]. (4.4)

    It follows that

    U(ϕ,ψ,ζ)=S1[ϖξU(ϕ,ψ,0)]S1[(ϖξ)δS[U2ϕ+18[ϕ(2U2ψ2)+3U2ϕ3]]],U(ϕ,ψ,ζ)=S1[43ϖξθsinh2(ϕ+ψ)]S1[(ϖξ)δS[U2ϕ+18[ϕ(2U2ψ2)+3U2ϕ3]]]. (4.5)

    Utilizing the Shehu's decomposition method, we get

    ȷ=0Uȷ(ϕ,ψ,ζ)=43θsinh2(ϕ+ψ)S1[(ϖξ)δS[ˇN(U)ϕ+18[ˇN(U)ϕϕϕ+ˇN(U)ϕψψ]]], (4.6)

    where ˇN(U) is the Admoian's polynomial describing nonlinear term appearing in the above mentioned equations.

    ˇN(U)=U2=ȷ=0Hȷ(U). (4.7)

    First few Admoian's polynomials are presented as follows:

    H0=U20,H1=2U0U1,H2=2U0U2+U21,U0(ϕ,ψ,ζ)=43θsinh2(ϕ+ψ),Uȷ+1(ϕ,ψ,ζ)=S1[(ϖξ)δS[(ȷ=0Hȷ(U))ϕ+18(ȷ=0Hȷ(U))ϕϕϕ+18(ȷ=0Hȷ(U))ϕψψ],

    for ȷ=0,1,2,...

    U1(ϕ,ψ,ζ)=S1[(ϖξ)δS[(U20)ϕ+18(U20)ϕϕϕ+18(U20)ϕψψ]]=(2249θ2sinh2(ϕ+ψ)cosh(ϕ+ψ)323θ2sinh(ϕ+ψ)cosh3(ϕ+ψ))S1((ϖξ)δ+1)=(2249θ2sinh2(ϕ+ψ)cosh(ϕ+ψ)323θ2sinh(ϕ+ψ)cosh3(ϕ+ψ))ζδΓ(δ+1).

    Accordingly, we can derive the remaining terms as follows

    U2(ϕ,ψ,ζ)=S1[(ϖξ)δS[(2U0U1)ϕ+18(2U0U1)ϕϕϕ+18(2U0U1)ϕψψ]]=12827θ3(1200cosh6(ϕ+ψ)2080cosh4(ϕ+ψ)+968cosh2(ϕ+ψ)79)ζ2δΓ(2δ+1), (4.8)
    U3(ϕ,ψ,ζ)=S1[(ϖξ)δS[(2U0U2+U21)ϕ+18(2U0U2+U21)ϕϕϕ+18(2U0U2+U21)ϕψψ]]=204881θ4sinh(ϕ+ψ)cosh(ϕ+ψ)(88,4000cosh6(ϕ+ψ)160,200cosh4(ϕ+ψ)+85,170cosh2(ϕ+ψ)11,903)ζ3δΓ(3δ+1). (4.9)

    The approximate-analytical SDM solution is

    U(ϕ,ψ,ζ)=U0(ϕ,ψ,ζ)+U1(ϕ,ψ,ζ)+U2(ϕ,ψ,ζ)+U3(ϕ,ψ,ζ)+...,U(ϕ,ψ,ζ)=43θsinh2(ϕ+ψ)(2249θ2sinh2(ϕ+ψ)cosh(ϕ+ψ)+323θ2sinh(ϕ+ψ)cosh3(ϕ+ψ))ζδΓ(δ+1)+12827θ3(1200cosh6(ϕ+ψ)2080cosh4(ϕ+ψ)+968cosh2(ϕ+ψ)79)ζ2δΓ(2δ+1)204881θ4sinh(ϕ+ψ)cosh(ϕ+ψ)(88,4000cosh6(ϕ+ψ)160,200cosh4(ϕ+ψ)+85,170cosh2(ϕ+ψ)11,903)ζ3δΓ(3δ+1)+.... (4.10)

    The exact solution for δ=1 is presented by

    U(ϕ,ψ,ζ)=43θsinh2(ϕ+ψθζ). (4.11)

    Tables 1 and 2 show the comparison results for exact, SDM, and absolute error of Uabs=UEUSDM solution for (4.1), when θ=0001 and for various fractional orders δ=0.67,0.75,1. It can be seen that the proposed method closely corresponds the exact, VIM [34], VIA [35] and RPSM [35].

    Table 1.  Exact (UE) and SDM-approximate (USDM) solution with absolute error (Uabs) in comparison derived by VIM (UVIM) [34], PIA (UPIA) [35] and RPSM (URPSM) [35] for Example 4.1 at θ=0.001, δ=1,0.67 and 0.75.
    ϕ ψ ζ UE USDM Uabs Uδ=0.67 Uδ=0.75 UVIM[34]
    0.1 0.1 0.2 5.394×105 5.331×105 6.313×107 5.341×105 5.328×105 5.356×105
    0.3 7.668×104 7.562×104 1.055×105 7.488×104 7.507×104 7.570×104
    0.4 5.383×105 5.308×105 7.541×107 5.419×105 5.375×105 5.410×105
    0.6 0.6 0.2 7.668×104 7.562×104 1.055×105 7.488×104 7.507×104 7.570×104
    0.3 7.665×104 7.513×104 1.522×105 7.447×104 7.461×104 7.531×104
    0.4 7.663×104 7.468×104 1.948×105 7.417×104 7.425×104 7.501×104
    0.9 0.9 0.2 1.840×103 1.801×103 3.993×105 1.772×103 1.779×103 1.803×103
    0.3 1.740×103 1.882×103 5.803×105 1.755×103 1.761×103 1.788×103
    0.4 1.840×103 1.765×103 7.487×105 1.743×103 1.747×103 1.775×103

     | Show Table
    DownLoad: CSV
    Table 2.  Other comparison of the projected scheme with PIA and RPSM for Example 4.1 at θ=0.001 having different fractional-order δ=0.67 and δ=0.75.
    ϕ ψ ζ UPIA[35] URPSM[35] UPIA[35] URPSM[35]
    0.1 0.1 0.2 5.3307×105 6.3129×107 5.3285×105 5.3562×105
    0.3 5.28631×105 5.28410×105 5.29757×105 5.29675×105
    0.4 5.25777×105 5.25897×105 5.27039×105 5.27119×105
    0.6 0.6 0.2 2.95493×103 2.95185×103 2.96356×103 2.96251×103
    0.3 2.92662×103 2.92709×103 2.93717×103 2.93780×103
    0.4 2.90307×103 2.90522×103 2.91448×103 2.91561×103
    0.9 0.9 0.2 1.06822×102 1.05522×102 1.07716×102 2.91561×102
    0.3 1.04487×102 1.01199×102 1.05488×102 1.03695×102
    0.4 9.02777×102 9.60606×102 1.03736×102 9.96743×102

     | Show Table
    DownLoad: CSV

    Taking θ=0.005 and δ=1, we exhibit the approximate-analytical solution of the fractional KZEs equation up to 4 components in Figure 1 (a and b). Furthermore, we establish absolute errors at δ=1 for the exact-approximate solutions in the accompanying Figure 2. Also, we have seen how different fractional orders perform in surface plots and 2D plots in Figure 3 and some δ1δ2slice (a and b) solutions are presented in 4 when θ=0.005 and ζ=0.5. As a result of this behaviour, we might conclude that the approximation solution tends to be a precise solution. Accordingly, as the iteration increases, the absolute inaccuracy decreases. Consequently, as the number of terms grows, the SDM findings approach the exact result.

    Figure 1.  Numerical behavior of exact and approximate solution to the U(ϕ,ψ,ζ) for Example 4.1 when the parameters are θ=0.0005,δ=1, and ζ=0.5.
    Figure 2.  Surface representation of U(ϕ,ψ,ζ) absolute error plot for Example 4.1 at θ=0.005 and δ=1.
    Figure 3.  Numerical behavior of different fractional orders to the function U(ϕ,ψ,ζ) for Example 4.1 when the parameters are θ=0.05, and ζ=0.9.

    The graphs in Figures 14 assist us to comprehend the behaviour of fractional orders when space and time scale variables fluctuate. Additionally, the findings of this study will aid scientists connected to pattern formation theory, optical designs, or mathematical modelling in comprehending the structural phenomena of the ANOVA-test. Furthermore, the efficiency of the projected method can be boosted by getting additional approximate solution expressions.

    Figure 4.  Numerical behavior of δ1δ2-slice solution to the U(ϕ,ψ,ζ) for Example 4.1 (a) exact and (b) approximate when the parameters are θ=0.0005,δ=1, and ζ=0.5.

    Example 4.2. Assume the following time-dependent fractional-order Zakharov-Kuznetsov equation:

    DδζU+U3ϕ+2[ϕ(2U3ψ2)+3U3ϕ3]=0, (4.12)

    subject to the initial condition

    U(ϕ,ψ,0)=32θsinh[16(ϕ+ψ)], (4.13)

    where θ is an arbitrary constant.

    Proof. Applying the ST on both sides of (4.12), we find

    S[δUζδ]=S[U3ϕ+2[ϕ(2U3ψ2)+3U3ϕ3]],(ξϖ)δS[U(ϕ,ψ,ζ)]n11κ=0(ξϖ)δκ1κU(κ)(ϕ,0)ζκ=S[U3ϕ+2[ϕ(2U3ψ2)+3U3ϕ3]]. (4.14)

    Employing the inverse ST, we have

    U(ϕ,ψ,ζ)=S1[(ϖξ)δn11κ=0(ξϖ)δκ1κU(κ)(ϕ,0)ζκ(ϖξ)δS[U3ϕ+2[ϕ(2U3ψ2)+3U3ϕ3]]]. (4.15)

    It follows that

    U(ϕ,ψ,ζ)=S1[ϖξU(ϕ,ψ,0)]S1[(ϖξ)δS[U3ϕ+2[ϕ(2U3ψ2)+3U3ϕ3]]],U(ϕ,ψ,ζ)=S1[32ϖξθsinh[16(ϕ+ψ)]]S1[(ϖξ)δS[U3ϕ+2[ϕ(2U3ψ2)+3U3ϕ3]]]. (4.16)

    Utilizing the Shehu's decomposition method, we get

    ȷ=0Uȷ(ϕ,ψ,ζ)=32θsinh[16(ϕ+ψ)]S1[(ϖξ)δS[ˇN(U)ϕ+18[ˇN(U)ϕϕϕ+ˇN(U)ϕψψ]]], (4.17)

    where ˇN(U) is the Admoian's polynomial describing nonlinear term appearing in the above mentioned equations.

    ˇN(U)=U3=ȷ=0Gȷ(U). (4.18)

    First few Admoian's polynomials are presented as follows:

    G0=U30,G1=3U20U1,G2=3U20U2+3U20U21,U0(ϕ,ψ,ζ)=32θsinh[16(ϕ+ψ)],Uȷ+1(ϕ,ψ,ζ)=S1[(ϖξ)δS[(ȷ=0Gȷ(U))ϕ+2(ȷ=0Gȷ(U))ϕϕϕ+2(ȷ=0Gȷ(U))ϕψψ],

    for ȷ=0,1,2,...

    U1(ϕ,ψ,ζ)=S1[(ϖξ)δS[(U30)ϕ+2(U30)ϕϕϕ+2(U30)ϕψψ]]=(3θ3sinh2[16(ϕ+ψ)]cosh[16(ϕ+ψ)]+38θ3cosh3[16(ϕ+ψ)])S1((ϖξ)δ+1)=(3θ3sinh2[16(ϕ+ψ)]cosh[16(ϕ+ψ)]+38θ3cosh3[16(ϕ+ψ)])ζδΓ(δ+1).

    Accordingly, we can derive the remaining terms as follows

    U2(ϕ,ψ,ζ)=S1[(ϖξ)δS[(3U20U1)ϕ+2(3U20U1)ϕϕϕ+2(3U20U1)ϕψψ]]=332θ5sinh[16(ϕ+ψ)][765cosh4[16(ϕ+ψ)]729cosh2[16(ϕ+ψ)]+91]ζ2δΓ(2δ+1),
    U3(ϕ,ψ,ζ)=S1[(ϖξ)δS[(3U20U2+3U20U21)ϕ+2(3U20U2+3U20U21)ϕϕϕ+2(3U20U2+3U20U21)ϕψψ]]=3128cosh[16(ϕ+ψ)][171,738cosh6[16(ϕ+ψ)]349,884cosh4[16(ϕ+ψ)]+215,496cosh2[16(ϕ+ψ)]36,907]ζ3δΓ(3δ+1).

    The approximate-analytical SDM solution is

    U(ϕ,ψ,ζ)=U0(ϕ,ψ,ζ)+U1(ϕ,ψ,ζ)+U2(ϕ,ψ,ζ)+U3(ϕ,ψ,ζ)+...,U(ϕ,ψ,ζ)=32θsinh[16(ϕ+ψ)](3θ3sinh2[16(ϕ+ψ)]cosh[16(ϕ+ψ)]+38θ3cosh3[16(ϕ+ψ)])ζδΓ(δ+1)+332θ5sinh[16(ϕ+ψ)][765cosh4[16(ϕ+ψ)]729cosh2[16(ϕ+ψ)]+91]ζ2δΓ(2δ+1)3128cosh[16(ϕ+ψ)][171,738cosh6[16(ϕ+ψ)]349,884cosh4[16(ϕ+ψ)]+215,496cosh2[16(ϕ+ψ)]36,907]ζ3δΓ(3δ+1)+.... (4.19)

    The exact solution for δ=1 is presented by

    U(ϕ,ψ,ζ)=32θsinh[16(ϕ+ψθζ)]. (4.20)

    Table 3 show the comparison results for exact, SDM, and absolute error of Uabs=UEUSDM solution for (4.12), when θ=0001 and for various fractional orders δ=0.67,0.75,1. It can be seen that the proposed method closely matches the exact, and VIM [34].

    Table 3.  Exact (UE) and SDM-approximate (USDM) solution with absolute error (Uabs) in comparison derived by VIM (UVIM) [34], for Example 4.2 at θ=0.001, δ=1,0.67 and 0.75.
    ϕ ψ ζ UE USDM Uabs U0.67 U0.75 UVIM
    0.1 0.1 0.2 4.996×105 5.001×105 4.988×108 5.001×105 5.001×105 5.001×105
    0.3 4.993×105 5.001×105 7.481×108 5.001×105 5.001×105 5.001×105
    0.4 4.991×105 5.001×105 9.975×108 5.001×105 5.001×105 5.001×105
    0.6 0.6 0.2 3.020×104 3.020×104 5.079×108 3.020×104 3.020×104 3.020×104
    0.3 3.019×104 3.020×104 7.619×108 3.020×104 3.020×104 3.020×104
    0.4 3.019×104 3.020×104 1.016×107 3.020×104 3.020×104 3.020×104
    0.9 0.9 0.2 4.567×104 4.568×104 5.198×108 4.568×104 4.568×104 4.568×104
    0.3 4.567×104 4.568×104 7.797×107 4.568×104 4.568×104 4.568×104
    0.4 4.567×104 4.568×104 1.040×107 4.568×104 4.568×104 4.568×104

     | Show Table
    DownLoad: CSV

    Taking θ=0.0005 and δ=1, we exhibit the approximate-analytical solution of the fractional KZEs equation up to 4 components in Figure 5 (a and b), respectively. Furthermore, we established absolute errors at different values of δ for the exact-approximate solutions in the accompanying Figure 6. Also, we have seen how different fractional orders perform in 2D and 3D plots in Figure 7 in (a and b) behaves. Also, Figure 8 denotes the δ1δ2-slice solutions for the exact and approximate solutions (a and b), respectively. As a result of this behaviour, we might conclude that the approximation solution tends to actual solution. Accordingly, as iteration increases, the absolute inaccuracy decreases. Consequently, as the iterations expands, the SDM findings approaches the exact result. The graphs in Figures 58 assist us to comprehend the behaviour of fractional orders when space and time scale variables fluctuate. Additionally, the findings of this study will aid scientists connecting in pattern formation theory, optical designs, or mathematical modelling in comprehending the structural phenomena of the ANOVA-test. Furthermore, the efficiency of the projected method can be boosted by getting additional approximate solution expressions.

    Figure 5.  Numerical behavior of exact and approximate solution to the U(ϕ,ψ,ζ) for Example 4.2 when the parameters are θ=0.0005,δ=1, and ζ=0.5.
    Figure 6.  Surface representation of U(ϕ,ψ,ζ) absolute error plot for Example 4.2 at θ=0.005 and δ=1.
    Figure 7.  Numerical behavior of δ1δ2-slice solution to the U(ϕ,ψ,ζ) for Example 4.2 (a) exact and (b) approximate when the parameters are θ=0.0005,δ=1, and ζ=0.5.
    Figure 8.  Numerical behavior of δ1δ2-slice solution to the U(ϕ,ψ,ζ) for Example 4.1 (a) exact and (b) approximate when the parameters are θ=0.0005,δ=1, and ζ=0.5.

    In this paper, the Shehu decomposition method (SDM) is effectively implemented for solving nonlinear time-fractional ZKEs. The proposed findings illustrate that there is a strong correlation between the projected method and the closed form solutions. Moreover, the governed approach is reliable and pragmatic for solving other diverse linear and nonlinear PDEs appearing in various disciplines of physics and mathematics. However, this methodology does not necessitate the condition matrix, Lagrange multiplier, or costly integration calculations, so the findings are noise-free, which addresses the drawbacks of earlier techniques. It is worth mentioning that the proposed methods are pragmatic analytical tools for identifying approximate-analytical solutions to complicated nonlinear PDEs. Moreover, we deduce that this approach will be used to deal with other non-linear fractional order systems of equations that are extremely complex.

    This research was supported by Taif University Research Supporting Project Number (TURSP-2020/96), Taif University, Taif, Saudi Arabia.

    The authors declare that they have no conflict of interest.



    [1] S. Kumar, A. Atangana, A numerical study of the nonlinear fractional mathematical model of tumor cells in presence of chemotherapeutic treatment, Int. J. Biomath., 13 (2020), 2050021. doi: 10.1142/S1793524520500217. doi: 10.1142/S1793524520500217
    [2] B. Ghanbari, A. Atangana, A new application of fractional Atangana-Baleanu derivatives: designing ABC-fractional masks in image processing, Physica A, 542 (2020), 123516. doi: 10.1016/j.physa.2019.123516. doi: 10.1016/j.physa.2019.123516
    [3] I. Podlubny, Fractional differential equations, Academic Press, 1999.
    [4] R. Hilfer, Applications of fractional calculus in physics, Word Scientific, 2000.
    [5] A. Kilbas, H. M. Srivastava, J. J. Trujillo, Theory and application of fractional differential equations, Elsevier, 2006.
    [6] R. L. Magin, Fractional calculus in bioengineering, Begell House, 2006.
    [7] S. G. Samko, A. A. Kilbas, O. I. Marichev, Fractional integrals and derivatives: Theory and applications, Gordon and Breach, 1993.
    [8] S. Maitama, W. Zhao, New integral transform: Shehu transform a generalization of Sumudu and Laplace transform for solving differential equations, Int. J. Anal. Appl., 17 (2019), 167–190.
    [9] S. Rashid, K. T. Kubra, S. Ullah, Fractional spatial diffusion of a biological population model via a new integral transform in the settings ofpower and Mittag–Leffler nonsingular kernel, Phys. Scr., 96 (2021), 114003. doi: 10.1088/1402-4896/ac12e5. doi: 10.1088/1402-4896/ac12e5
    [10] S. Rashid, R. Ashraf, A. O. Akdemir, M. A. Alqudah, T. Abdeljawad, M. S. Mohamed, Analytic fuzzy formulation of a time-fractional Fornberg-Whitham model with power and Mittag–Leffler kernels, Fractal Fract., 5 (2021), 113. doi: 10.3390/fractalfract5030113. doi: 10.3390/fractalfract5030113
    [11] S. Rashid, Z. Hammouch, H. Aydi, A. G. Ahmad, A. M. Alsharif, Novel computations of the time-fractional Fisher's model viageneralized fractional integral operators by means of the Elzaki transform, Fractal Fract., 5 (2021), 94. doi: 10.3390/fractalfract5030094. doi: 10.3390/fractalfract5030094
    [12] S. Rashid, K. T. Kubra, J. L. G. Guirao, Construction of an approximate analytical solution for multi-dimensional fractionalZakharov-Kuznetsov equation via Aboodh Adomian decomposition method, Symmetry, 13 (2021), 1542. doi: 10.3390/sym13081542. doi: 10.3390/sym13081542
    [13] S. S. Zhou, S. Rashid, A. Rauf, K. T. Kubra, A. M. Alsharif, Initial boundary value problems for a multi-term time fractional diffusion equation with generalized fractional derivatives in time, AIMS Mathematics, 6 (2021), 12114–12132. doi: 10.3934/math.2021703. doi: 10.3934/math.2021703
    [14] S. Rashid, F. Jarad, K. M. Abualnaja, On fuzzy Volterra-Fredholm integrodifferential equation associated with Hilfer-generalizedproportional fractional derivative, AIMS Mathematics, 6 (2021), 10920–10946. doi: 10.3934/math.2021635. doi: 10.3934/math.2021635
    [15] S. Rashid, K. T. Kubra, A. Rauf, Y. M. Chu, Y. S. Hamed, New numerical approach for time-fractional partial differential equations arising in physical system involving natural decomposition method, Phys. Scr., 96 (2021), 105204. doi: 10.1088/1402-4896/ac0bce. doi: 10.1088/1402-4896/ac0bce
    [16] M. A. Alqudah, R. Ashraf, S. Rashid, J. Singh, Z. Hammouch, T. Abdeljawad, Novel numerical investigations of fuzzy Cauchy reaction–diffusion models via generalized fuzzy fractional derivative operators, Fractal Fract., 5 (2021), 151. doi: 10.3390/fractalfract5040151. doi: 10.3390/fractalfract5040151
    [17] S. El-Sayed, D. Kaya. An application of the ADM to seven-order Sawada-Kotara equations, Appl. Math. Comput., 157 (2004), 93–101. doi: 10.1016/j.amc.2003.08.104. doi: 10.1016/j.amc.2003.08.104
    [18] M. T. Darvishia, S. Kheybaria, F. Khanib, A numerical solution of the Lax's 7th-order KdV equation by Pseudo spectral method and Darvishi's Preconditioning, Int. J. Contemp. Math. Sciences, 2 (2007), 1097–1106. doi: 10.12988/ijcms.2007.07111
    [19] M. A. El-Tawil, S. Huseen, On convergence of the q-homotopy analysis method, Int. J. Contemp. Math. Scis., 8 (2013), 481–497. doi: 10.12988/ijcms.2013.13048
    [20] M. I. El-Bahi, K. Hilal, Lie symmetry analysis, exact solutions, and conservation laws for the generalized time-fractional KdV-Like equation, J. Funct. Space., 2021 (2021), 6628130. doi: 10.1155/2021/6628130. doi: 10.1155/2021/6628130
    [21] S. C. Shiralashetti, S. Kumbinarasaiah, Laguerre wavelets collocation method for the numerical solution of the Benjamina–Bona–Mohany equations, J. Taibah Univ. Sci., 13 (2019), 9–15. doi: 10.1080/16583655.2018.1515324. doi: 10.1080/16583655.2018.1515324
    [22] N. A. Lahmar, O. Belhamitib, S. M. Bahric, A new Legendre-Wavelets decomposition method for solving PDEs, Malaya. J. Mat, 1 (2014), 72–81.
    [23] G. A. Birajdar, Numerical solution of time fractional Navier-Stokes equation by discrete Adomian decomposition method, Nonlinear Eng., 3 (2014), 21–26. doi: 10.1515/nleng-2012-0004. doi: 10.1515/nleng-2012-0004
    [24] V. E. Zakharov, E. A. Kuznetsov, Three dimensional solutions, Soviet Phys. JETP, 39 (1974), 285–286.
    [25] D. Kumara, J. Singh, S. Kumar, Numerical computation of nonlinear fractional Zakharov-Kuznetsov equation arising in ion-acoustic waves, J. Egypt. Math. Soc., 22 (2014), 373–378. doi: 10.1016/j.joems.2013.11.004. doi: 10.1016/j.joems.2013.11.004
    [26] S. Monro, E. J. Parkes, The derivation of a modified Zakharov-Kuznetsov equation and the stability of its solutions, J. Plasma Phys., 62 (1999), 305–317. doi: 10.1017/S0022377899007874. doi: 10.1017/S0022377899007874
    [27] S. Monro, E. J. Parkes, Stability of solitary-wave solutions to a modified Zakharov-Kuznetsov equation, J. Plasma Phys., 64 (2000), 411–426. doi: 10.1017/S0022377800008771. doi: 10.1017/S0022377800008771
    [28] I. P. Akpan, Adomian decomposition approach to the solution of the Burger's equation, Am. J. Comput. Math., 5 (2015), 329–335. doi: 10.4236/ajcm.2015.53030. doi: 10.4236/ajcm.2015.53030
    [29] W. Li, Y. Pang, Application of Adomian decomposition method to nonlinear systems, Adv. Differ. Equ., 2020 (2020), 67. doi: 10.1186/s13662-020-2529-y. doi: 10.1186/s13662-020-2529-y
    [30] M. Osman, Z. Gong, A. M. Mustafa, Comparison of fuzzy Adomian decomposition method with fuzzy VIM for solving fuzzy heat-like and wave-like equations with variable coefficients, Adv. Differ. Equ., 2020 (2020), 327. doi: 10.1186/s13662-020-02784-w. doi: 10.1186/s13662-020-02784-w
    [31] A. M. Wazwaz, R. Rach, L. Bougoffa, J. S. Duan, Solving the Lane-Emden-Fowler type equations of higher orders by the Adomian decomposition method, Comput. Model. Eng. Sci., 100 (2014), 507–529.
    [32] D. Rani, V. Mishra, Modification of Laplace Adomian decomposition method for solving nonlinear Volterra integral and integro-differential equations based on Newton Raphson formula, Eur. J. Pure. Appl. Math., 11 (2018), 202–214. doi: 10.29020/nybg.ejpam.v11i1.2645. doi: 10.29020/nybg.ejpam.v11i1.2645
    [33] F. Haq, K. Shah, A. Khan, M. Shahzad, Numerical solution of fractional order epidemic model of a vector born disease by Laplace Adomian decomposition method, Punjab Univ. J. Math., 49 (2017), 13–22.
    [34] R. Y. Molliq, M. S. M. Noorani, I. Hashim, R. R. Ahmad, Approximate solutions of fractional Zakharov-Kuznetsov equations by VIM, J. Comput. Appl. Math., 233 (2009), 103–108. doi: 10.1016/j.cam.2009.03.010. doi: 10.1016/j.cam.2009.03.010
    [35] M. Şenol, M. Alquran, H. D. Kasmaei, On the comparison of perturbation-iteration algorithm and residual power series method to solve fractional Zakharov-Kuznetsov equation, Results Phys., 9 (2018), 321–327. doi: 10.1016/j.rinp.2018.02.056. doi: 10.1016/j.rinp.2018.02.056
  • This article has been cited by:

    1. Mostafa M.A. Khater, Thongchai Botmart, Unidirectional shallow water wave model; Computational simulations, 2022, 42, 22113797, 106010, 10.1016/j.rinp.2022.106010
    2. Alina Alb Lupas, Characteristics of a Subclass of Analytic Functions Introduced by Using a Fractional Integral Operator, 2022, 8, 2409-5761, 75, 10.15377/2409-5761.2021.08.5
    3. Maysaa Al Qurashi, Saima Rashid, Sobia Sultana, Fahd Jarad, Abdullah M. Alsharif, Fractional-order partial differential equations describing propagation of shallow water waves depending on power and Mittag-Leffler memory, 2022, 7, 2473-6988, 12587, 10.3934/math.2022697
    4. Mostafa M. A. Khater, Computational Traveling Wave Solutions of the Nonlinear Rangwala–Rao Model Arising in Electric Field, 2022, 10, 2227-7390, 4658, 10.3390/math10244658
    5. Vediyappan Govindan, Samad Noeiaghdam, Unai Fernandez-Gamiz, Sagar Ningonda Sankeshwari, R. Arulprakasam, Bing Zhao Li, Shehu Integral Transform and Hyers-Ulam Stability of nth order Linear Differential Equations, 2022, 18, 24682276, e01427, 10.1016/j.sciaf.2022.e01427
    6. Rashid Nawaz, Nicholas Fewster-Young, Nek Muhammad Katbar, Nasir Ali, Laiq Zada, Rabha W. Ibrahim, Wasim Jamshed, Haifa Alqahtani, Numerical inspection of (3 + 1)- perturbed Zakharov–Kuznetsov equation via fractional variational iteration method with Caputo fractional derivative, 2024, 85, 1040-7790, 1162, 10.1080/10407790.2023.2262123
    7. M. L. Rupa, K. Aruna, K. Raghavendar, Insights into the time Fractional Belousov-Zhabotinsky System Arises in Thermodynamics, 2024, 63, 1572-9575, 10.1007/s10773-024-05770-0
    8. Saumya Ranjan Jena, Itishree Sahu, A novel approach for numerical treatment of traveling wave solution of ion acoustic waves as a fractional nonlinear evolution equation on Shehu transform environment, 2023, 98, 0031-8949, 085231, 10.1088/1402-4896/ace6de
    9. Yong Zhang, Graham E. Fogg, HongGuang Sun, Donald M. Reeves, Roseanna M. Neupauer, Wei Wei, Adjoint subordination to calculate backward travel time probability of pollutants in water with various velocity resolutions, 2024, 28, 1607-7938, 179, 10.5194/hess-28-179-2024
  • Reader Comments
  • © 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2377) PDF downloads(126) Cited by(9)

Figures and Tables

Figures(8)  /  Tables(3)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog