Research article Special Issues

Numerical simulation of fractional-order two-dimensional Helmholtz equations

  • In this paper, we investigate the exact solutions of several fractional-order Helmholtz equations using the homotopy perturbation transform method. We specify sufficient requirements for its convergence and provide error estimations. The homotopy perturbation transform method yields a quickly converging succession of solutions. Solutions for various fractional space derivatives are compared to present approaches and explained using figures. Appropriate parameter selection produces approximations identical to the exact answer. Test examples are provided to demonstrate the proposed approach's precision and competence. The results demonstrate that our system is appealing, user-friendly, dependable, and highly effective.

    Citation: Naveed Iqbal, Muhammad Tajammal Chughtai, Nehad Ali Shah. Numerical simulation of fractional-order two-dimensional Helmholtz equations[J]. AIMS Mathematics, 2023, 8(6): 13205-13218. doi: 10.3934/math.2023667

    Related Papers:

    [1] Naveed Iqbal, Saleh Alshammari, Thongchai Botmart . Evaluation of regularized long-wave equation via Caputo and Caputo-Fabrizio fractional derivatives. AIMS Mathematics, 2022, 7(11): 20401-20419. doi: 10.3934/math.20221118
    [2] Manal Alqhtani, Khaled M. Saad, Rasool Shah, Thongchai Botmart, Waleed M. Hamanah . Evaluation of fractional-order equal width equations with the exponential-decay kernel. AIMS Mathematics, 2022, 7(9): 17236-17251. doi: 10.3934/math.2022949
    [3] Aslı Alkan, Halil Anaç . A new study on the Newell-Whitehead-Segel equation with Caputo-Fabrizio fractional derivative. AIMS Mathematics, 2024, 9(10): 27979-27997. doi: 10.3934/math.20241358
    [4] Muath Awadalla, Abdul Hamid Ganie, Dowlath Fathima, Adnan Khan, Jihan Alahmadi . A mathematical fractional model of waves on Shallow water surfaces: The Korteweg-de Vries equation. AIMS Mathematics, 2024, 9(5): 10561-10579. doi: 10.3934/math.2024516
    [5] Zeliha Korpinar, Mustafa Inc, Dumitru Baleanu . On the fractional model of Fokker-Planck equations with two different operator. AIMS Mathematics, 2020, 5(1): 236-248. doi: 10.3934/math.2020015
    [6] Ved Prakash Dubey, Jagdev Singh, Ahmed M. Alshehri, Sarvesh Dubey, Devendra Kumar . Analysis of local fractional coupled Helmholtz and coupled Burgers' equations in fractal media. AIMS Mathematics, 2022, 7(5): 8080-8111. doi: 10.3934/math.2022450
    [7] Azzh Saad Alshehry, Naila Amir, Naveed Iqbal, Rasool Shah, Kamsing Nonlaopon . On the solution of nonlinear fractional-order shock wave equation via analytical method. AIMS Mathematics, 2022, 7(10): 19325-19343. doi: 10.3934/math.20221061
    [8] 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
    [9] Shabir Ahmad, Aman Ullah, Ali Akgül, Fahd Jarad . A hybrid analytical technique for solving nonlinear fractional order PDEs of power law kernel: Application to KdV and Fornberg-Witham equations. AIMS Mathematics, 2022, 7(5): 9389-9404. doi: 10.3934/math.2022521
    [10] Manoj Singh, Ahmed Hussein, Msmali, Mohammad Tamsir, Abdullah Ali H. Ahmadini . An analytical approach of multi-dimensional Navier-Stokes equation in the framework of natural transform. AIMS Mathematics, 2024, 9(4): 8776-8802. doi: 10.3934/math.2024426
  • In this paper, we investigate the exact solutions of several fractional-order Helmholtz equations using the homotopy perturbation transform method. We specify sufficient requirements for its convergence and provide error estimations. The homotopy perturbation transform method yields a quickly converging succession of solutions. Solutions for various fractional space derivatives are compared to present approaches and explained using figures. Appropriate parameter selection produces approximations identical to the exact answer. Test examples are provided to demonstrate the proposed approach's precision and competence. The results demonstrate that our system is appealing, user-friendly, dependable, and highly effective.



    In recent years, fractional differential equations have gained prominence due to their proven usefulness in several unrelated scientific and engineering fields. For example, the nonlinear oscillations of an earthquake can be characterized by a fractional derivative, and the fractional derivative of the traffic fluid dynamics model can solve the insufficiency resulting from the assumption of continuous traffic flows [1,3]. Numerous chemical processes, mathematical biology, engineering, and scientific problems [4,5,6,7] are also modeled with fractional differential equations. Nonlinear partial differential equations (NPDEs) characterize various physical, biological, and chemical phenomena. Current research is focused on developing precise traveling wave solutions for such equations. Exact and explicit solutions help scientists understand the complicated physical phenomena and dynamic processes portrayed by NPDEs [8,9,10]. In the past four decades, numerous essential methodologies for attaining accurate solutions to NPDEs have been proposed [11,12].

    The Helmholtz equation (HE) derives from the elliptic and wave equations. In a multi-dimensional nonhomogeneous isotropic standard with velocity c, the wave result is υ(ξ,ψ), which corresponds to a source of harmonic (ξ,ψ) vibrating at a given frequency and satisfying the Helmholtz equation in the area R. The classical order HE is

    D2ξυ(ξ,ψ)+D2ψυ(ξ,ψ)+ευ(ξ,ψ)=υ(ξ,ψ). (1.1)

    Here, υ is a suitable boundary differentiable term of R, is a known function, and the wave number with wavelength 2/ξ=0 renders Eq (1.1) homogeneous. If (1.1) is expressed as

    D2ξυ(ξ,ψ)+D2ψυ(ξ,ψ)ευ(ξ,ψ)=υ(ξ,ψ).

    Then it explains mass transfer with density biochemical processes of the 1st order. Equation (1.1) is investigated using the decomposition method [13], the finite element approach [14], the differential transform method [15], the Trefftz method [16], and the spectral collocation method [17], among others [3,4,5].

    The Helmholtz equation is a partial differential equation that describes wave phenomena in various fields of physics, such as electromagnetism, acoustics, and fluid mechanics. Traditionally, the Helmholtz equation has been formulated using integer-order derivatives. However, in recent years, there has been a growing interest in the use of fractional-order derivatives to describe complex phenomena more accurately. In particular, fractional-order space Helmholtz equations are derived directly from mathematical formulas that involve fractional derivatives, rather than being generalized from integer-order space derivative Helmholtz equations. These equations can provide a more accurate description of wave propagation in complex media, such as porous materials, biological tissues, and fractal structures. Fractional-order space Helmholtz equations have attracted significant attention due to their potential applications in a wide range of fields, including medical imaging, geophysics, and telecommunications. They offer a promising avenue for understanding the behavior of waves in complex media and developing new technologies for wave-based sensing and imaging [6,7,8].

    It is advantageous to utilize fractional differential equations in physical problems due to their nonlocal features. Non-locality characterizes fractional-order derivatives, whereas locality characterizes integer-order derivatives [24,25,26,27]. It demonstrates that the future state of the physical system depends on all of its previous states in addition to its current state. Consequently, fractional models are more accurate. In fractional differential equations, the response expression has a parameter that specifies the fractional derivative of the variable order, which may vary to achieve many responses [9,10,11].

    Standard HEs can be generalized to fractional-order Helmholtz equations by extending the Caputo fractional-order space derivative to the integer-order space derivative. The fractional Helmholtz equation in space is

    Dϱξυ(ξ,ψ)+D2ψυ(ξ,ψ)+ευ(ξ,ψ)=ψ(ξ,ψ),

    with υ(0,ψ)=g(ψ) as the initial condition (IC). Gupta et al. [31] solved the multi-dimensional fractional Helmholtz equation using the homotopy perturbation approach. In contrast, Abuasad et al. [14] recently solved a fractional model of the Helmholtz problem using the reduced differential transform method.

    This section describes the properties of the fractional derivatives and a few essential details concerning the Yang transform.

    Definition 2.1. The fractional derivative in terms of Caputo is as follows

     Dϱψυ(ξ,ψ)=1Γ(kϱ)ψ0(ψϱ)kϱ1υ(k)(ξ,ϱ)dϱ,k1<ϱk,kN. (2.1)

    Definition 2.2. The YT is represented as follows

    Y{υ(ψ)}=M(u)=0eψuυ(ψ)dψ,  ψ>0,  u(ψ1,ψ2), (2.2)

    having inverse YT as follows

    Y1{M(u)}=υ(ψ). (2.3)

    Definition 2.3. The nth derivative YT is stated as follows

    Y{υn(ψ)}=M(u)unn1k=0υk(0)unk1,    n=1,2,3, (2.4)

    Definition 2.4. The YT of derivative having fractional-order is stated as follows

    Y{υϱ(ψ)}=M(u)uϱn1k=0υk(0)uϱ(k+1),  0<ϱn. (2.5)

    Consider the general fractional partial differential equations,

    Dϱψυ(ξ,ψ)+Mυ(ξ,ψ)+Nυ(ξ,ψ)=h(ξ,ψ),ψ>0,0<ϱ1,υ(ξ,0)=g(ξ),ν. (3.1)

    Using Yang transform of Eq (3.1), we get

     Y[Dϱyυ(ξ,ψ)+Mυ(ξ,ψ)+Nυ(ξ,ψ)]=Y[h(ξ,ψ)],ψ>0,0<ϱ1,υ(ξ,ψ)=sg(ξ)+sϱY[h(ξ,ψ)]sϱY[Mυ(ξ,ψ)+Nυ(ξ,ψ)]. (3.2)

    Now, applying inverse Yang transform, we have

    υ(ξ,ψ)=F(ξ,ψ)Y1[sϱY{Mυ(ξ,ψ)+Nυ(ξ,ψ)}], (3.3)

    where

     F(ξ,ψ)=Y1[sg(ξ)+sϱY[h(ξ,ψ)]]=g(ν)+Y1[sϱY[h(ξ,ψ)]]. (3.4)

    The parameter p is perturbation technique and p[0,1] defined as

    υ(ξ,ψ)=k=0pkυk(ξ,ψ), (3.5)

    The nonlinear function is expressed as

    Nυ(ξ,ψ)=k=0pkHk(υk), (3.6)

    where Hn are He's polynomials in term of υ0,υ1,υ2,,υn, and can be calculated as

    Hn(υ0,υ1,,υn)=1ϱ(n+1)Dıp[N(ı=0pıυı)]p=0, (3.7)

    where Dıp=ıpı.

    Putting Eqs (3.6) and (3.7) in Eq (3.3), we achieved as

     ı=0pıυı(ξ,ψ)=F(ξ,ψ)p×[Y1{sϱY{Mı=0pıυı(ξ,ψ)+ı=0pıHı(υı)}}]. (3.8)

    Comparison both sides of coefficient p, we get

    p0:υ0(ξ,ψ)=F(ξ,ψ),p1:υ1(ξ,ψ)=Y1[sϱY(Mυ0(ξ,ψ)+H0(υ))],p2:υ2(ξ,ψ)=Y1[sϱY(Mυ1(ξ,ψ)+H1(υ))],pı:υı(ξ,ψ)=Y1[sϱY(Mυı1(ξ,ψ)+Hı1(υ))],ı>0,ıN. (3.9)

    Finally, present the obtained solution and check it with any available analytical or numerical solutions for the given PDE. The υı(ξ,ψ) components can be calculated easily which quickily converges to series form. We can get p1,

     υ(ξ,ψ)=limMMı=1υı(ξ,ψ). (3.10)

    Problem 4.1. Consider the space fractional Helmholtz equation

    ϱυ(ξ,ψ)ξϱ+2υ(ξ,ψ)ψ2υ(ξ,ψ)=0,1<ϱ2, (4.1)

    with the ICs

     υ(0,ψ)=ψ  andυξ(0,ψ)=0. (4.2)

    Using the Yang transform of Eq (4.1), we obtained as

     1sϱY[υ(ξ,ψ)]=υ(0,ψ)s1ϱY{2υ(ξ,ψ)ψ2υ(ξ,ψ)}, (4.3)
     Y[υ(ξ,ψ)]=sυ(0,ψ)sϱY{2υ(ξ,ψ)ψ2υ(ξ,ψ)}, (4.4)

    Taking inverse Yang Transformation, we have

    Y[υ(ξ,ψ)]=ψY1[sϱY{2υ(ξ,ψ)ψ2υ(ξ,ψ)}], (4.5)

    Implemented HPM in Eq (4.5), we can achieve as

     ı=0pıυı(ξ,ψ)=ψp[Y1{sϱY{(ı=0pıυı(ξ,ψ))ψψı=0pıυı(ξ,ψ)}}]. (4.6)

    On both sides comparing coefficients of p, we get

     p0:υ0(ξ,ψ)=ψ,p1:υ1(ξ,ψ)=Y1[sϱY{2υ0(ξ,ψ)ψ2υ0(ξ,ψ)}]=ξϱΓ(ϱ+1)ψ,p2:υ2(ξ,ψ)=Y1[sϱY{2υ1(ξ,ψ)ψ2υ1(ξ,ψ)}]=ξ2ϱΓ(2ϱ+1)ψ,p3:υ3(ξ,ψ)=Y1[sϱY{2υ2(ξ,ψ)ψ2υ2(ξ,ψ)}]=ξ3ϱΓ(3ϱ+1)ψ,p4:υ4(ξ,ψ)=Y1[sϱY{2υ3(ξ,ψ)ψ2υ3(ξ,ψ)}]=ξ4ϱΓ(4ϱ+1)ψ, (4.7)

    The series type result of the first problem example is

     υ(ξ,ψ)=υ0(ξ,ψ)+υ1(ξ,ψ)+υ2(ξ,ψ)+υ3(ξ,ψ)+υ4(ξ,ψ)+υ(ξ,ψ)=ψ[1+ξϱΓ(ϱ+1)+ξ2ϱΓ(2ϱ+1)+x3ϱΓ(3ϱ+1)+ξ4ϱΓ(4ϱ+1)+]. (4.8)

    The exact solution is

    υ(ξ,ψ)=ψcoshξ.

    Similarly y-space can be calculated as:

    ϱυ(ξ,ψ)ψϱ+2υ(ξ,ψ)ξ2υ(ξ,ψ)=0, (4.9)

    with the IC

    υ(ξ,0)=ξ. (4.10)

    Thus, the solution of the above Eq (4.9) is obtain as

    υ(ξ,ψ)=ξ(1+ψϱΓ(ϱ+1)+ψ2ϱΓ(2ϱ+1)+ψ3ϱΓ(3ϱ+1)+ψ4ϱΓ(4ϱ+1)+),

    in the case when ϱ=2, then the solution through HPTM is

    υ(ξ,ψ)=ξcoshψ. (4.11)

    Figure 1 illustrates the exact and HPTM solutions in two-dimensional plots for various values of ϱ ranging from 2 to 1.5, with ξ values ranging from 0 to 1 and ψ set to 1. The solutions are presented in Figure 1a and b for the exact and HPTM methods, respectively. Figure 2 displays the 3-dimensional plots of the exact and HPTM solutions for ϱ=2 and analyzes the point of intersection between the two solutions. Figure 2c and d depict the HPTM solutions at ϱ=1.8 and 1.6 respectively for Problem 4.1. The fractional results were also evaluated for their convergence towards an integer-order result for each problem. Similarly, the figures for the ψ-space can also be generated using the same approach.

    Figure 1.  The first graph show that the exact and approximate solution and second HPTM solution at the different fractional-order graph of Problem 4.1.
    Figure 2.  Exact and proposed method solution at various fractional orders of Problem 4.1.
    Table 1.  Exact and proposed method solution of Problem 4.1 at various fractional orders.
    (ξ,ψ) υ(ξ,ψ) at ϱ=1.5 υ(ξ,ψ) at ϱ=1.75 (HPTM) at ϱ= 2 Exact result
    (0.2, 0.1) 0.2102371 0.2100072 0.2100000 0.2100000
    (0.4, 0.1) 0.4104630 0.4100141 0.4100000 0.4100000
    (0.6, 0.1) 0.6106889 0.6100209 0.6100000 0.6100000
    (0.2, 0.2) 0.2103355 0.2100121 0.2100000 0.2100000
    (0.4, 0.2) 0.4106550 0.4100237 0.4100000 0.4100000
    (0.6, 0.2) 0.6109746 0.6100353 0.6100000 0.6100000
    (0.2, 0.3) 0.2104110 0.2100164 0.2100000 0.2100000
    (0.4, 0.3) 0.4108025 0.4100321 0.4100000 0.4100000
    (0.6, 0.3) 0.6111940 0.6100478 0.6100000 0.6100000
    (0.2, 0.4) 0.2104747 0.2100204 0.2100000 0.2100000
    (0.4, 0.4) 0.4109269 0.4100399 0.4100000 0.4100000
    (0.6, 0.4) 0.6113790 0.6100593 0.6100000 0.6100000
    (0.2, 0.5) 0.2105309 0.2100241 0.2100000 0.2100000
    (0.4, 0.5) 0.4110365 0.4100471 0.4100000 0.4100000
    (0.6, 0.5) 0.6115421 0.6100701 0.6100000 0.6100000

     | Show Table
    DownLoad: CSV

    Problem 4.2. Consider the space-fractional HE

    ϱυ(ξ,ψ)ξϱ+2υ(ξ,ψ)ψ2+5υ(ξ,ψ)=0,1<ϱ2, (4.12)

    with the ICs

     υ(0,ψ)=ψ  andυξ(0,ψ)=0. (4.13)

    Using the Yang transform of Eq (4.12), we obtain as

     1sϱY[υ(ξ,ψ)]=υ(0,ψ)s1ϱY{2υ(ξ,ψ)ψ2+5υ(ξ,ψ)}, (4.14)
     Y[υ(ξ,ψ)]=sυ(0,ψ)sϱY{2υ(ξ,ψ)ψ2+5υ(ξ,ψ)}. (4.15)

    Applying the inverse Yang Transform, we get

    Y[υ(ξ,ψ)]=ψY1[sϱY{2υ(ξ,ψ)ψ2+5υ(ξ,ψ)}], (4.16)

    Using the HPM in Eq (4.16), we obtained as

     ı=0pıυı(ξ,ψ)=ψp[Y1{sϱY{(ı=0pıυı(ξ,ψ))ψψ+5ı=0pıυı(ξ,ψ)}}]. (4.17)

    On both sides comparing coefficients of p, we get

     p0:υ0(ξ,ψ)=ψ,p1:υ1(ξ,ψ)=Y1[sϱY{2υ0(ξ,ψ)ψ2+5υ0(ξ,ψ)}]=5ψξϱΓ(ϱ+1),p2:υ2(ξ,ψ)=Y1[sϱY{2υ1(ξ,ψ)ψ2+5υ1(ξ,ψ)}]=25ψξ2ϱΓ(2ϱ+1),p3:υ3(ξ,ψ)=Y1[sϱY{2υ2(ξ,ψ)ψ2+5υ2(ξ,ψ)}]=125ξ3ϱΓ(3ϱ+1),p4:υ4(ξ,ψ)=Y1[sϱY{2υ3(ξ,ψ)ψ2+5υ3(ξ,ψ)}]=625ψξ4ϱΓ(4ϱ+1), (4.18)

    The series type result of second problem as

     υ(ξ,ψ)=υ0(ξ,ψ)+υ1(ξ,ψ)+υ2(ξ,ψ)+υ3(ξ,ψ)+υ4(ξ,ψ)+υ(ξ,ψ)=ψ[15ξϱΓ(ϱ+1)+25ξ2ϱΓ(2ϱ+1)125x3ϱΓ(3ϱ+1)+625ξ4ϱΓ(4ϱ+1)+]. (4.19)

    The exact solution is

    υ(ξ,ψ)=ψcos5ξ.

    Now similarly, the result of y-space can be calculated with the help of homotopy perturbation

    ϱυ(ξ,ψ)ψϱ+2υ(ξ,ψ)ξ2+5υ(ξ,ψ)=0, (4.20)

    with the IC

    υ(ξ,0)=ξ. (4.21)

    The solution of the Eq (4.20) is expressed as

    υ(ξ,ψ)=ξ(15ψϱΓ(ϱ+1)+25ψ2ϱΓ(2ϱ+1)125ψ3ϱΓ(3ϱ+1)+625ψ4ϱΓ(4ϱ+1)+).

    The exact solution is

    υ(ξ,ψ)=ξcos5ψ. (4.22)

    Figure 3 illustrates the solutions of exact and HPTM in two-dimensional plots, as shown in Figure 3a and b for different values of ϱ, ranging from 2 to 1.5, respectively. The interval considered for ξ is [0, 1], while ψ is constant at 1. The results obtained from the fractional-order model converge to the integer-order solution of the problem. In Figure 4, the 3-dimensional plots of exact and HPTM solutions are presented in Figures (a) and (b), respectively, for ϱ=2. The closed contact of the two solutions is analyzed. Additionally, Figure 4c and d depict the HPTM solutions at ϱ=1.8 and 1.6, respectively, for Problem 4.2. Similarly, graphs for ψ-space can also be generated.

    Figure 3.  Exact and proposed method solution at various fractional orders of Problem 4.2.
    Figure 4.  Exact and proposed method solution at various fractional orders of Problem 4.2.

    Problem 4.3. Consider the space-fractional HE

    ϱυ(ξ,ψ)ξϱ+2υ(ξ,ψ)ψ22υ(ξ,ψ)=(12ξ23ξ4)sinψ,1<ϱ2,0ψ2π, (4.23)

    with the ICs

     υ(0,ψ)=0  andυξ(0,ψ)=0. (4.24)

    Applying the Yang transform of Eq (4.23), we achieve

     1sϱY[υ(ξ,ψ)]=υ(0,ψ)s1ϱY{2υ(ξ,ψ)ψ22υ(ξ,ψ)}, (4.25)
     Y[υ(ξ,ψ)]=sυ(0,ψ)sϱY{2υ(ξ,ψ)ψ22υ(ξ,ψ)}. (4.26)

    Implementing inverse Yang transform, we get

    Y[υ(ξ,ψ)]=(ξ4ξ610)sinψY1[sϱY{2υ(ξ,ψ)ψ22υ(ξ,ψ)}]. (4.27)

    Applying Homotopy perturbation method in Eq (4.27), we achieved as

     ı=0pıυı(ξ,ψ)=ψp[Y1{sϱY{(ı=0pıυı(ξ,ψ))ψψ2ı=0pıυı(ξ,ψ)}}]. (4.28)

    Both sides on comparison coefficients of p, we obtain

     p0:υ0(ξ,ψ)=(ξ4ξ610)sinψ,p1:υ1(ξ,ψ)=Y1[sϱY{2υ0(ξ,ψ)ψ22υ0(ξ,ψ)}]=3(ξϱ+4Γ(ϱ+5)72ξϱ+6Γ(ϱ+7))sinψ,p2:υ2(ξ,ψ)=Y1[sϱY{2υ1(ξ,ψ)ψ22υ1(ξ,ψ)}]=3(ξ2ϱ+4Γ(2ϱ+5)216ξ2ϱ+6Γ(2ϱ+7))sinψ,p3:υ3(ξ,ψ)=Y1[sϱY{2υ2(ξ,ψ)ψ22υ2(ξ,ψ)}]=3(ξ3ϱ+4Γ(3ϱ+5)648ξ3ϱ+6Γ(3ϱ+7))sinψ,p4:υ4(ξ,ψ)=Y1[sϱY{2υ3(ξ,ψ)ψ22υ3(ξ,ψ)}]=3(ξ4ϱ+4Γ(4ϱ+5)1944ξ2ϱ+6Γ(2ϱ+7))sinψ, (4.29)

    The series type result of the third problem is

     υ(ξ,ψ)=υ0(ξ,ψ)+υ1(ξ,ψ)+υ2(ξ,ψ)+υ3(ξ,ψ)+υ4(ξ,ψ)+υ(ξ,ψ)=(ξ4ξ610)sinψ+3(ξϱ+4Γ(ϱ+5)72ξϱ+6Γ(ϱ+7))sinψ+3(ξ2ϱ+4Γ(2ϱ+5)216ξ2ϱ+6Γ(2ϱ+7))sinψ+3(ξ3γ+4Γ(3ϱ+5)648ξ3ϱ+6Γ(3ϱ+7))sinψ+3(ξ4ϱ+4Γ(4ϱ+5)1944ξ2ϱ+6Γ(2ϱ+7))sinψ+. (4.30)

    The exact solution is

    υ(ξ,ψ)=ξ4sinψ.

    Figure 5a and b display the exact and HPTM solutions, respectively, in a 3-dimensional plot at ϱ=2. The closed contact between the exact and HPTM solutions is examined. Figure 6 depicts the exact and HPTM solutions in two-dimensional plot for various values of ϱ=2,1.9,1.8,1.7,1.6,1.5 for ξ[0,1] and ψ=1. The fractional results are observed to approach an integer-order solution of the problem. Similarly, the graphs for ψ-space fractional-order derivative can also be plotted.

    Figure 5.  Exact and proposed method solution at various fractional orders of Problem 4.3.
    Figure 6.  Exact and proposed method solution at various fractional orders of Problem 4.3.

    In this study, fractional-order Helmholtz equations were solved using the Homotopy Perturbation Yang transform method. Due to the great agreement between the generated approximative solution and the precise solution, the homotopy perturbation Yang transform method was demonstrated to be a successful method for solving partial differential equations with Caputo operators. The computation size of the approach was compared to those required by other numerical methods to demonstrate how tiny it is. Additionally, the procedure's quick convergence demonstrates its dependability and marks a notable advancement in the way linear and non-linear fractional-order partial differential equations are solved.

    This research has been funded by Deputy for Research & Innovation, Ministry of Education through Initiative of Institutional Funding at University of Ha'il–Saudi Arabia through project number IFP-22 064.

    The authors declare that they have no competing interests.



    [1] J. Li, X. N. Su, K. Y. Zhao, Barycentric interpolation collocation algorithm to solve fractional differential equations, Math. Comput. Simulat., 205 (2023), 340–367. https://doi.org/10.1016/j.matcom.2022.10.005 doi: 10.1016/j.matcom.2022.10.005
    [2] M. M. Al-Sawalha, R. P. Agarwal, R. Shah, O. Y. Ababneh, W. Weera, A reliable way to deal with fractional-order equations that describe the unsteady flow of a polytropic gas, Mathematics, 10 (2022), 2293. https://doi.org/10.3390/math10132293 doi: 10.3390/math10132293
    [3] L. A. Said, A. H. Madian, A. G. Radwan, A. M. Soliman, Fractional order oscillator with independent control of phase and frequency, In 2014 2nd International Conference on Electronic Design (ICED), 2014,224–229. https://doi.org/10.1109/ICED.2014.7015803
    [4] J. T. Machado, V. Kiryakova, F. Mainardi, Recent history of fractional calculus, Commun. Nonlinear Sci. Numer. Simul., 16 (2011), 1140–1153. https://doi.org/10.1016/j.cnsns.2010.05.027 doi: 10.1016/j.cnsns.2010.05.027
    [5] J. Sabatier, O. P. Agrawal, J. A. Tenreiro Machado, Advances in Fractional Calculus, Dordrecht: Springer, 2007.
    [6] D. Baleanu, K. Diethelm, E. Scalas, J. J. Trujillo, Fractional Calculus: Models and Numerical Methods, World Scientific, 2012. https://doi.org/10.1142/10044
    [7] L. Debnath, Recent applications of fractional calculus to science and engineering, Int. J. Math. Math. Sci., 2003 (2003), 1–30. https://doi.org/10.1155/S0161171203301486 doi: 10.1155/S0161171203301486
    [8] M. M. Al-Sawalha, A. S. Alshehry, K. Nonlaopon, R. Shah, O. Y. Ababneh, Fractional view analysis of delay differential equations via numerical method, AIMS Mathematics, 7 (2022), 20510–20523. https://doi.org/10.3934/math.20221123 doi: 10.3934/math.20221123
    [9] S. Mukhtar, R. Shah, S. Noor, The numerical investigation of a fractional-order multi-dimensional model of Navier-Stokes equation via novel techniques, Symmetry, 14 (2022), 1102. https://doi.org/10.3390/sym14061102 doi: 10.3390/sym14061102
    [10] M. M. Al-Sawalha, A. S. Alshehry, K. Nonlaopon, R. Shah, O. Y. Ababneh, Approximate analytical solution of time-fractional vibration equation via reliable numerical algorithm, AIMS Mathematics, 7 (2022), 19739–19757. https://doi.org/10.3934/math.20221082 doi: 10.3934/math.20221082
    [11] M. M. Al-Sawalha, R. Shah, A. Khan, O. Y. Ababneh, T. Botmart, Fractional view analysis of Kersten-Krasil'shchik coupled KdV-mKdV systems with non-singular kernel derivatives, AIMS Mathematics, 7 (2022), 18334–18359. https://doi.org/10.3934/math.20221010 doi: 10.3934/math.20221010
    [12] Y. Kai, S. Q. Chen, K. Zhang, Z. X. Yin, Exact solutions and dynamic properties of a nonlinear fourth-order time-fractional partial differential equation, Wave. Random Complex, 2022. https://doi.org/10.1080/17455030.2022.2044541 doi: 10.1080/17455030.2022.2044541
    [13] F. Ihlenburg, I. Babuska, Finite element solution of the Helmholtz equation with high wave number part Ⅱ: The h-p version of the FEM, SIAM J. Numer. Anal., 34 (1997), 315–358. https://doi.org/10.1137/S0036142994272337 doi: 10.1137/S0036142994272337
    [14] S. M. El-Sayed, D. Kaya, Comparing numerical methods for Helmholtz equation model problem, Appl. Math. Comput., 150 (2004), 763–773. https://doi.org/10.1016/S0096-3003(03)00305-9 doi: 10.1016/S0096-3003(03)00305-9
    [15] Y. K. Cheung, W. G. Jin, O. C. Zienkiewicz, Solution of Helmholtz equation by Trefftz method, Int. J. Numer. Meth. Eng., 32 (1991), 63–78. https://doi.org/10.1002/nme.1620320105 doi: 10.1002/nme.1620320105
    [16] A. Prakash, M. Goyal, S. Gupta, Numerical simulation of space-fractional Helmholtz equation arising in seismic wave propagation, imaging and inversion, Pramana, 93 (2019), 28. http://doi.org/10.1007/s12043-019-1773-8 doi: 10.1007/s12043-019-1773-8
    [17] S. Nguyen, C. Delcarte, A spectral collocation method to solve Helmholtz problems with boundary conditions involving mixed tangential and normal derivatives, J. Comput. Phys., 200 (2004), 34–49. https://doi.org/10.1016/j.jcp.2004.03.004 doi: 10.1016/j.jcp.2004.03.004
    [18] X. Li, Z. Q. Dong, L. P. Wang, X. D. Niu, H. Yamaguchi, D. C. Li, et al., A magnetic field coupling fractional step lattice Boltzmann model for the complex interfacial behavior in magnetic multiphase flows, Appl. Math. Model., 117 (2023), 219–250. https://doi.org/10.1016/j.apm.2022.12.025 doi: 10.1016/j.apm.2022.12.025
    [19] X. L. Xie, T. F. Wang, W. Zhang, Existence of solutions for the (p, q)-Laplacian equation with nonlocal Choquard reaction, Appl. Math. Lett., 135 (2023), 108418. https://doi.org/10.1016/j.aml.2022.108418 doi: 10.1016/j.aml.2022.108418
    [20] L. Wang, H. Zhao, X. Liu, Z. L. Zhang, X. H. Xia, S. Evans, Optimal remanufacturing service resource allocation for generalized growth of retired mechanical products: Maximizing matching efficiency, IEEE Access, 9 (2021), 89655–89674. https://doi.org/10.1109/ACCESS.2021.3089896 doi: 10.1109/ACCESS.2021.3089896
    [21] S. Abuasad, K. Moaddy, I. Hashim, Analytical treatment of two-dimensional fractional Helmholtz equations, J. King Saud Univ. Sci., 31 (2019), 659–666. https://doi.org/10.1016/j.jksus.2018.02.002 doi: 10.1016/j.jksus.2018.02.002
    [22] X. Wang, X. J. Lyu, Experimental study on vertical water entry of twin spheres side-by-side, Ocean Eng., 221 (2021), 108508. https://doi.org/10.1016/j.oceaneng.2020.108508 doi: 10.1016/j.oceaneng.2020.108508
    [23] Y. Hu, J. X. Qing, Z. H. Liu, Z. J. Conrad, J. N. Cao, X. P. Zhang, Hovering efficiency optimization of the ducted propeller with weight penalty taken into account, Aerosp. Sci. Technol., 117 (2021), 106937. https://doi.org/10.1016/j.ast.2021.106937 doi: 10.1016/j.ast.2021.106937
    [24] H. Y. Jin, Z. A. Wang, Global stabilization of the full attraction-repulsion Keller-Segel system, Discrete Contin. Dyn. Syst., 40 (2020), 3509–3527. https://doi.org/10.3934/dcds.2020027 doi: 10.3934/dcds.2020027
    [25] H. Y. Jin, Z. A. Wang, Boundedness, blowup and critical mass phenomenon in competing chemotaxis, J. Differ. Equ., 260 (2016), 162–196. https://doi.org/10.1016/j.jde.2015.08.040 doi: 10.1016/j.jde.2015.08.040
    [26] H. Y. Jin, Z. A. Wang, Asymptotic dynamics of the one-dimensional attraction-repulsion Keller-Segel model, Math. Meth. Appl. Sci., 38 (2015), 444–457. https://doi.org/10.1002/mma.3080 doi: 10.1002/mma.3080
    [27] L. Liu, S. Zhang, L. Ch. Zhang, G. Pan, J. Z. Yu, Multi-UUV maneuvering counter-game for dynamic target scenario based on fractional-order recurrent neural network, IEEE Trans. Cybernetics, 2022, 1–14. https://doi.org/10.1109/TCYB.2022.3225106 doi: 10.1109/TCYB.2022.3225106
    [28] N. Iqbal, M. T. Chughtai, R. Ullah, Fractional study of the non-linear Burgers' equations via a semi-analytical technique, Fractal Fract., 7 (2023), 103. https://doi.org/10.3390/fractalfract7020103 doi: 10.3390/fractalfract7020103
    [29] P. Liu, J. P. Shi, Z. A. Wang, Pattern formation of the attraction-repulsion Keller-Segel system, Discrete Contin. Dyn. Syst. Ser. B, 18 (2013), 2597–2625. https://doi.org/10.3934/dcdsb.2013.18.2597 doi: 10.3934/dcdsb.2013.18.2597
    [30] M. Alesemi, N. Iqbal, N. Wyal, Novel evaluation of fuzzy fractional Helmholtz equations, J. Funct. Spaces, 2022 (2022), 8165019. https://doi.org/10.1155/2022/8165019 doi: 10.1155/2022/8165019
    [31] P. K. Gupta, A. Yildirim, K. N. Rai, Application of He's homotopy perturbation method for multi-dimensional fractional Helmholtz equation, Internat. J. Numer. Methods Heat Fluid Flow, 22 (2012), 424–435. https://doi.org/10.1108/09615531211215738 doi: 10.1108/09615531211215738
  • This article has been cited by:

    1. Muhammad Nadeem, Zitian Li, Devendra Kumar, Yahya Alsayaad, A robust approach for computing solutions of fractional-order two-dimensional Helmholtz equation, 2024, 14, 2045-2322, 10.1038/s41598-024-54870-8
    2. Jinxing Liu, Muhammad Nadeem, Asad Islam, Sorin Mureşan, Loredana Florentina Iambor, A Modified Residual Power Series Method for the Approximate Solution of Two-Dimensional Fractional Helmholtz Equations, 2023, 15, 2073-8994, 2152, 10.3390/sym15122152
    3. Jinxing Liu, Muhammad Nadeem, Yahya Alsayyad, António M. Lopes, Approximate solution of Newell-Whitehead-Segel model with time-fractional derivative, 2023, 18, 1932-6203, e0288740, 10.1371/journal.pone.0288740
    4. Muhammad Nadeem, Mohamed Sharaf, Saipunidzam Mahamad, Numerical investigation of two-dimensional fractional Helmholtz equation using Aboodh transform scheme, 2024, 0961-5539, 10.1108/HFF-07-2024-0543
  • Reader Comments
  • © 2023 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(1346) PDF downloads(93) Cited by(4)

Figures and Tables

Figures(6)  /  Tables(1)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog