Processing math: 78%
Research article

Analysis of local fractional coupled Helmholtz and coupled Burgers' equations in fractal media

  • In this paper, we present a computational algorithm, namely, local fractional natural homotopy analysis method (LFNHAM) to explore the solutions of local fractional coupled Helmholtz and local fractional coupled Burgers' equations (LFCHEs and LFCBEs). This work also investigates the uniqueness and convergence of the solution of a general local fractional partial differential equation (LFPDE) obtained by the suggested method in view of theory of fixed point and Banach spaces. Furthermore, the error analysis of the LFNHAM solution is also discussed. Moreover, the numerical simulations are presented for each of the local fractional coupled equations on the Cantor set. The computational procedure clearly illustrates the validity and reliability of the proposed method for achieving the solutions of local fractional coupled Helmholtz and coupled Burgers' equations. The proposed method also minimizes the computational work unlike other conventional methods while still giving extremely precise results. The implemented combination supplies a more general solution as compared to other methods and assimilates their consequences as a special case. In addition, the acquired solutions are also in excellent match with previously determined solutions.

    Citation: 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[J]. AIMS Mathematics, 2022, 7(5): 8080-8111. doi: 10.3934/math.2022450

    Related Papers:

    [1] Mubashir Qayyum, Efaza Ahmad, Hijaz Ahmad, Bandar Almohsen . New solutions of time-space fractional coupled Schrödinger systems. AIMS Mathematics, 2023, 8(11): 27033-27051. doi: 10.3934/math.20231383
    [2] 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
    [3] 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
    [4] Sunil Kumar, Amit Kumar , Zaid Odibat, Mujahed Aldhaifallah, Kottakkaran Sooppy Nisar . A comparison study of two modified analytical approach for the solution of nonlinear fractional shallow water equations in fluid flow. AIMS Mathematics, 2020, 5(4): 3035-3055. doi: 10.3934/math.2020197
    [5] Mohammed Al-Refai, Dumitru Baleanu . Comparison principles of fractional differential equations with non-local derivative and their applications. AIMS Mathematics, 2021, 6(2): 1443-1451. doi: 10.3934/math.2021088
    [6] Thongchai Botmart, Ravi P. Agarwal, Muhammed Naeem, Adnan Khan, Rasool Shah . On the solution of fractional modified Boussinesq and approximate long wave equations with non-singular kernel operators. AIMS Mathematics, 2022, 7(7): 12483-12513. doi: 10.3934/math.2022693
    [7] Alessandra Jannelli, Maria Paola Speciale . On the numerical solutions of coupled nonlinear time-fractional reaction-diffusion equations. AIMS Mathematics, 2021, 6(8): 9109-9125. doi: 10.3934/math.2021529
    [8] K. Pavani, K. Raghavendar . A novel method to study time fractional coupled systems of shallow water equations arising in ocean engineering. AIMS Mathematics, 2024, 9(1): 542-564. doi: 10.3934/math.2024029
    [9] Changdev P. Jadhav, Tanisha B. Dale, Vaijanath L. Chinchane, Asha B. Nale, Sabri T. M. Thabet, Imed Kedim, Miguel Vivas-Cortez . On solutions of fractional differential equations for the mechanical oscillations by using the Laplace transform. AIMS Mathematics, 2024, 9(11): 32629-32645. doi: 10.3934/math.20241562
    [10] Naveed Iqbal, Muhammad Tajammal Chughtai, Nehad Ali Shah . Numerical simulation of fractional-order two-dimensional Helmholtz equations. AIMS Mathematics, 2023, 8(6): 13205-13218. doi: 10.3934/math.2023667
  • In this paper, we present a computational algorithm, namely, local fractional natural homotopy analysis method (LFNHAM) to explore the solutions of local fractional coupled Helmholtz and local fractional coupled Burgers' equations (LFCHEs and LFCBEs). This work also investigates the uniqueness and convergence of the solution of a general local fractional partial differential equation (LFPDE) obtained by the suggested method in view of theory of fixed point and Banach spaces. Furthermore, the error analysis of the LFNHAM solution is also discussed. Moreover, the numerical simulations are presented for each of the local fractional coupled equations on the Cantor set. The computational procedure clearly illustrates the validity and reliability of the proposed method for achieving the solutions of local fractional coupled Helmholtz and coupled Burgers' equations. The proposed method also minimizes the computational work unlike other conventional methods while still giving extremely precise results. The implemented combination supplies a more general solution as compared to other methods and assimilates their consequences as a special case. In addition, the acquired solutions are also in excellent match with previously determined solutions.



    Helmholtz and Burgers' equations play an important role in various streams of applied physics. The Helmholtz equation frequently occurs in the study of physical phenomena involving elliptic partial differential equations (PDEs) such as wave and diffusion, magnetic fields, seismology, electromagnetic radiation, transmission, vibrating lines, acoustics, and geosciences. This equation is actually derived from the wave equation. The Helmholtz equation is a transformed form of the acoustic wave equation. It is utilized in a stream of seismic wave propagation and imaging. This equation plays a significant role in estimations of acoustic propagation in shallow water at low frequencies and characteristics of geodesic sea floor [1]. Mathematically, the eigenvalue problem for the Laplace operator is called the Helmholtz equation expressed by elliptic type linear PDE 2ϑ=k2ϑ, where 2 denotes the Laplacian differential operator, k2 signifies the eigenvalue, and ϑ is the eigen function. When this equation is used in respect of waves, k is termed as the wave number which measures the spatial frequency of waves. For the first time, Samuel and Thomas [2] suggested the Helmholtz equation with fractional order. Recently, Prakash et al. [3] presented the solution of the space-fractional Helmholtz equation with the q-homotopy analysis transform method (q-HATM). More recently, Shah et al. [4] examined the fractional Helmholtz equation also.

    On the other hand, the Burgers' equations [5,6,7] characterize the nonlinear diffusion phenomenon through the simplest PDEs. Burgers' equations occur mainly in the mathematical model of turbulence, fluid mechanics, and approximation of flow in viscous fluids [5,8,9]. The coupled Burgers' equations in one-dimensional form are described as a sedimentation and/or evolution model of scaled volume concentrations in fluid suspensions. More literature about coupled Burgers' equations can be found in previous works [10,11]. In view of the development of the fractional calculus approach, the Burgers' equation with a fractional derivative was first presented in [12]. After that, many authors investigated the solution for fractional Burgers' equations in past decades using approximate analytical methods (see, for example, [13,14,15,16,17,18,19,20,21,22,23]).

    From the past decade, the concept of local fractional calculus and local fractional derivatives developed in the work of Yang [24,25] has been a centre of attraction among researchers. Further, many authors investigated the equations and models appearing in fractal media through various local fractional methods, for instance, local fractional homotopy perturbation method (LFHPM) for handling local fractional PDEs (LFPDEs) [26,27], local fractional Tricomi equation arising in fractal transonic flow [28], local fractional Klein-Gordon equations [29], local fractional heat conduction equation [30], local fractional wave equation in fractal strings [31], local fractional Laplace equation [32], system of LFPDEs [33], and fractal vehicular traffic flow [34], etc. In this sequence, the 2D local fractional Helmholtz equation (LFHE) was introduced in [35]. Recently, the LFHE was solved by local fractional variational iteration method [36], local fractional series expansion method [37]. The local fractional Helmholtz and coupled Helmholtz equations were handled successfully by Baleanu and Jassim [38,39,40] through various local fractional methods. In recent years, the local fractional coupled Burgers' equations (LFCBEs) were also investigated for solutions through various techniques that can be found in [41,42,43,44,45].

    The 2D local fractional coupled Helmholtz equations (LFCHEs) suggested in [38] are given as follows:

    2σϑ1(γ,τ)τ2σ+2σϑ2(γ,τ)γ2σω12σϑ1(γ,τ)=1(γ,τ),0<σ1, (1.1)
    2σϑ2(γ,τ)τ2σ+2σϑ1(γ,τ)γ2σω22σϑ2(γ,τ)=2(γ,τ),0<σ1, (1.2)

    subject to the initial conditions:

    ϑ1(γ,0)=ϕ1(γ),σϑ1(γ,0)τσ=ψ1(γ),ϑ2(γ,0)=ϕ2(γ),σϑ2(γ,0)τσ=ψ2(γ), (1.3)

    where ϑ1(γ,τ) and ϑ2(γ,τ) are unknown local fractional continuous functions, and 1(γ,τ) and 2(γ,τ) are the nondifferentiable source terms.

    The system of nonlinear coupled Burger's equations with local fractional derivatives can be described as:

    σϑ1τσ+ξ12σϑ1γ2σ+ξ2σϑ1γσϑ1+ρσ[ϑ1ϑ2]γσ=0,ϑ1=ϑ1(γ,τ),0<σ1, (1.4)
    σϑ2τσ+μ12σϑ2γ2σ+μ2σϑ2γσϑ2+ησ[ϑ1ϑ2]γσ=0,ϑ2=ϑ2(γ,τ),0<σ1, (1.5)

    subject to the initial conditions:

    ϑ1(γ,0)=f1(γ),ϑ2(γ,0)=f2(γ), (1.6)

    where ξ1, ξ2, μ1, and μ2 denote real constants, ρ and η specify arbitrary constants that depend on parameters of the system, ϑ1(γ,τ) and ϑ2(γ,τ) are local fractional continuous functions, and γ lies in the computational domain Ω.

    The key purpose of this work is to establish a new coupling of local fractional homotopy analysis method (LFHAM) [43,46] and local fractional natural transform (LFNT) [47], named as local fractional natural homotopy analysis method (LFNHAM) throughout in this paper. The second goal of the paper is to explore the solutions for the LFCHEs and the LFCBEs by utilizing the newly suggested combination LFNHAM. Moreover, the numerical simulations have also been presented for the obtained solutions of LFCHEs and LFCBEs for the fractal order σ=ln2/ln3 of a local fractional derivative by using MATLAB. The originality and novelty of the paper lie in the fact that the LFCHEs and LFCBEs have never been solved by using this newly suggested combination LFNHAM. In addition, convergence and uniqueness of the LFNHAM solution are also examined for the LFNHAM solution of general LFPDE in view of Banach's fixed point theory.

    The notable aspect of the LFNHAM as compared to others is that it offers an extended degree of freedom for analysis and the main ingredient is an auxiliary parameter 0 to ensure the convergence of the acquired series solution. Furthermore, a more appropriate choice of an initial guess & effortless creation of deformation equations are the interesting attributes of this method. The LFNHAM is surely beneficial as it combines two powerful algorithms to attain the solutions for nonlinear LFDEs. The LFNHAM generates a convergent series solution that revolves around a convergence parameter without involving linearization, perturbation, or descretization phenomena. In addition, the LFNHAM also minimizes the numerical work unlike other conventional methods while still giving extremely precise results. The LFNHAM provides a more general solution as compared to LFHPM, local fractional Adomian decomposition method (LFADM) and local fractional natural homotopy perturbation method (LFNHPM) and assimilates their consequences as a special case. In addition, it does not involve the computation of complicated Adomian or He's polynomials. But there is also a point of demerit with this technique. The implementation of LFNHAM can be difficult in the situation of non-evaluation of the LFNT of a function. This work checks the LFNHAM solution regarding uniqueness and convergence for the first time and the error analysis of the LFNHAM solution is also discussed. These points surely illustrate the reliability and validity of the proposed method. The other aspect of the LFNHAM is that the coupling of LFNT with LFHAM performs fast-tracked calculations in comparison to LFHAM and consequently consumes less time and less computer memory.

    Moreover, the LFNT possesses two important attributes, scale property and unit-preserving property, and hence can be utilized to handle LFPDEs without exerting new frequency range. In the light of these facts, the LFNT which possesses the linearity feature, also possesses the feature of linearity of functions, and hence does not involve the changing of units. This transform performs operation similarly as the local fractional Laplace transform (LFLT) and local fractional Sumudu transform (LFST). By virtue of these facts, the LFNT may be used to analyze some complex problems of science and engineering that may be handled hardly with other integral transforms.

    The rest portion of the paper is organized as follows: Section 2 presents definitions and formulae for the local fractional derivative and LFNT. Section 3 illustrates the computational procedure for the suggested scheme LFNHAM. The convergence and uniqueness of the LFNHAM solution is discussed in Section 4. Sections 5 and 6 are devoted to the implementation of the LFNHAM to the LFCHEs and LFCBEs, respectively. In Section 7, numerical simulations have been performed in respect of a fractal value. At the end, Section 8 presents the epilogue.

    The section presents a quick view of some definitions and formulae which have been utilized in this work.

    Definition 2.1.1. [24,25] Let (l1,l2) be the interval and Δt=max{Δt0,Δt1,Δt2,Δt3,...} be a partition of (l1,l2) with (tj,tj+1), j=0,...,N1, t0=l1, tN=l2 with Δtj=tj+1tj. Now, the LFI of ϑ(γ) is formulated as

    l1Iσl2ϑ(γ)=1Γ(1+σ)l2l1ϑ(t)(dt)σ=1Γ(1+σ)limΔt0N1j=0ϑ(tj)(Δtj)σ. (2.1)

    Definition 2.1.2. [24,25] The Mittag-Leffler function is given as

    Eσ(γσ)=q=0γqσΓ(1+qσ),0<σ1. (2.2)

    Definition 2.1.3. [24,25] The fractal sine and cosine functions are given by

    Sinσ(γσ)=q=0(1)qγ(2q+1)σΓ(1+(2q+1)σ),0<σ1, (2.3)
    Cosσ(γσ)=q=0(1)qγ2qσΓ(1+2qσ),0<σ1, (2.4)
    Sinhσ(γσ)=q=0γ(2q+1)σΓ(1+(2q+1)σ),0<σ1, (2.5)
    Coshσ(γσ)=q=0γ2qσΓ(1+2qσ),0<σ1. (2.6)

    Definition 2.1.4. [24,25] The LFD of ϑ(γ)Cσ(l1,l2) of order σ at γ=γ0 is presented as

    Dσγϑ(γ0)=dσϑ(γ0)dγσ=ϑσ(γ0)=Δσ(ϑ(γ)ϑ(γ0))(γγ0)σ,γ(l1,l2), (2.7)

    where Δσ(ϑ(γ)ϑ(γ0))Γ(σ+1)(ϑ(γ)ϑ(γ0)).

    The local fractional partial derivative of ϑ(γ,τ)Cσ(l1,l2) of order σ was provided by Yang [24,25] as follows:

    στσϑ(γ,τ)=Δσ(ϑ(γ,τ)ϑ(γ,τ0))(ττ0)σ,γ(l1,l2), (2.8)

    where Δσ(ϑ(γ,τ)ϑ(γ,τ0))Γ(σ+1)(ϑ(γ,τ)ϑ(γ,τ0)).

    The LFIs and LFDs of special functions used in this study and described in [24,25] are given as follows:

    Dσγaϑ(γ)=aDσγϑ(γ),Dσγ(γqσΓ(1+qσ))=γ(q1)σΓ(1+(q1)σ),qN,
    DσγEσ(γσ)=Eσ(γσ),Dσγ(Sinσ(γσ))=Cosσ(γσ),Dσγ(Cosσ(γσ))=Sinσ(γσ),
    Iσγ(γqσΓ(1+qσ))=γ(q+1)σΓ(1+(q+1)σ),qN,

    where γσ signifies a Cantor function.

    For the first time, Khan & Khan [48] suggested a new integral transform called N-transform. Some years later, Belgacem and Silambarasan [49,50,51] changed its name to Natural transform and also presented a comprehensive study regarding its applications. This transform performs operation similarly as the Laplace and Sumudu transforms.

    Definition 2.2.1. [47] The LFNT of the function ϑ(γ,τ) of order σ is defined as

    LFNσ[ϑ(γ,τ)]=ˉϑσ(γ,s,u)=1Γ(1+σ)0Eσ(sστσuσ)ϑ(γ,τ)uσ(dτ)σ,0<σ1, (2.9)

    and the corresponding inverse LFNT LFN1σ is formulated as

    LFN1σ[ˉϑσ(γ,s,u)]=ϑ(γ,τ)=1(2πi)σρ+iρiEσ(sστσuσ)ˉϑσ(γ,s,u)(ds)σ,0<σ1, (2.10)

    where sσ and uσ signify the LFNT variables and ρ denotes a real constant. The integral in the definition of inverse LFNT is taken along sσ=ρ in the complex plane sσ=xσ+iyσ. It is notable that the LFNT converges to LFLT for u=1 and to LFST for s=1.

    Some properties of the LFNT are being mentioned here:

    Proposition 2.2.2. [47] The LFNT of a LFD is defined by

    LFNσ[ϑ(qσ)(γ,τ)]=sqσuqσˉϑσ(γ,s,u)q1k=0s(qk1)σu(qk)σϑ(kσ)(γ,0). (2.11)

    For q=1,2 and 3, the following expressions are generated

    LFNσ[ϑ(σ)(γ,τ)]=sσuσˉϑσ(γ,s,u)1uσϑ(γ,0),
    LFNσ[ϑ(2σ)(γ,τ)]=s2σu2σˉϑσ(γ,s,u)sσu2σϑ(γ,0)1uσϑ(σ)(γ,0),
    LFNσ[ϑ(3σ)(γ,τ)]=s3σu3σˉϑσ(γ,s,u)s2σu3σϑ(γ,0)sσu2σϑ(σ)(γ,0)1uσϑ(2σ)(γ,0). (2.12)

    Proposition 2.2.3. [47] The linearity property of the LFNT is defined by

    LFNσ[α1ϑ1(γ,τ)+α2ϑ2(γ,τ)]=α1LFNσ[ϑ1(γ,τ)]+α2LFNσ[ϑ2(γ,τ)]=α1ˉϑ1,σ(γ,s,u)+α2ˉϑ2,σ(γ,s,u), (2.13)

    where ˉϑ1,σ(γ,s,u) and ˉϑ2,σ(γ,s,u) denote the LFNT of ϑ1(γ,τ) and ϑ2(γ,τ), respectively.

    Theorem 2.2.4. (Local fractional convolution). If LFNσ{ϑ1(γ,τ)}=ˉϑ1,σ(γ,s,u) and LFNσ{ϑ2(γ,τ)}=ˉϑ2,σ(γ,s,u), we have

    LFNσ{ϑ1(γ,τ)ϑ2(γ,τ)}=uσˉϑ1,σ(γ,s,u)ˉϑ2,σ(γ,s,u), (2.14)

    where

    ϑ1(γ,τ)ϑ2(γ,τ)=1Γ(1+σ)0ϑ1(γ,η)ϑ2(γ,τη)(dη)σ. (2.15)

    Some useful formulae for LFNT are listed in Table 1 [47].

    Table 1.  Formulae of LFNT.
    LFNσ(1)=1sσ LFNσ(Cosσ(τσ))=sσs2σ+u2σ
    LFNσ(τσΓ(1+σ))=uσs2σ LFNσ(Sinσ(τσ))=uσs2σ+u2σ
    LFNσ(τqσΓ(1+qσ))=uqσs(q+1)σ LFNσ(Coshσ(τσ))=sσs2σu2σ

     | Show Table
    DownLoad: CSV

    To explain the basic idea of LFNHAM, the following LFPDE is taken here

    Lσϑ(γ,τ)+Pσϑ(γ,τ)+Qσϑ(γ,τ)=ω(γ,τ),0<γ<1,0<τ<1, (3.1)

    where Lσqστqσ denotes the linear local fractional differential operator (LFDO) of order qσ i.e., a number δ>0 such that Lσϑδϑ, Pσ specifies the linear fractional differential operator of general nature in γ and τ. Here, it is also assumed that Pσ is bounded i.e., |Pσ(ϑϑ)|λ|ϑϑ|. Qσ denotes the nonlinear differential operator which is Lipschitz continuous with ξ>0 fulfilling the criteria |Qσ(ϑϑ)|ξ|ϑϑ|, γ and τ are independent variables, ϑ(γ,τ) and ω(γ,τ) denote local fractional unknown function and nondifferentiable source term, respectively.

    Now, the suggested computational approach recommends the implementation of the LFNT operator PNϵ on Eq (3.1)

    LFNσ[Lσϑ(γ,τ)]+LFNσ[Pσϑ(γ,τ)]+LFNσ[Qσϑ(γ,τ)]=LFNσ[ω(γ,τ)]. (3.2)

    Using the property of LFNT for LFDs, it follows

    ˉϑσ(γ,s,u)=uqσsqσq1k=0s(qk1)σu(qk)σϑ(kσ)(γ,0)+uqσsqσ(LFNσ[ω(γ,τ)])uqσsqσ(LFNσ[Pσϑ(γ,τ)]+LFNσ[Qσϑ(γ,τ)]), (3.3)
    whereˉϑσ(γ,s,u)=LFNσ[ϑ(γ,τ)]. (3.4)

    After simplification, we get

    ˉϑσ(γ,s,u)uqσsqσq1k=0s(qk1)σu(qk)σϑ(kσ)(γ,0)uqσsqσ(LFNσ[ω(γ,τ)])+uqσsqσ(LFNσ[Pσϑ(γ,τ)]+LFNσ[Qσϑ(γ,τ)])=0. (3.5)

    Now on account of Eq (3.5), the nonlinear operator is constituted as

    Φ[φ(γ,τ;κ)]=LFNσ[φ(γ,τ;κ)]uqσsqσq1k=0s(qk1)σu(qk)σφ(kσ)(γ,0;κ)uqσsqσ(LFNσ[ω(γ,τ)])+uqσsqσ(LFNσ[Pσφ(γ,τ;κ)]+LFNσ[Qσφ(γ,τ;κ)])=0, (3.6)

    where κ[0,1] is an embedding parameter, φ(γ,τ;κ) symbolizes the local fractional unknown function of γ,τ and κ, and the symbol LFNσ represents the LFNT operator.

    Now utilizing the traditional approach of LFHAM [43,46] and basic methodology of HAM [52,53,54], the zeroth-order deformation equation is developed in this way:

    (1κ)LFNσ[φ(γ,τ;κ)ϑ0(γ,τ)]=κΦ[φ(γ,τ;κ)], (3.7)

    where 0 is a convergence regulation parameter and ϑ0(γ,τ) symbolizes an initial guess for ϑ(γ,τ).

    It is observed that LFNHAM makes easy the choice of auxiliary parameters, linear operator, and initial guess. The following equations stand firmly for κ=0 and κ=1 in this manner

    φ(γ,τ;0)=ϑ0(γ,τ),φ(γ,τ;1)=ϑ(γ,τ). (3.8)

    Hence, when κ takes values from 0 to 1, φ(γ,τ;κ) deviates from ϑ0(γ,τ) to ϑ(γ,τ). Next, the Taylor's series expansion of φ(γ,τ;κ) about κ generates

    φ(γ,τ;κ)=ϑ0(γ,τ)+μ=1κμϑμ(γ,τ), (3.9)

    where

    ϑμ(γ,τ)=[1Γ(μ+1)μ[ϕ(γ,τ;κ)]κμ]κ=0. (3.10)

    The convergence controller 0 promptly provides the convergence of the series solution (3.9). Thus the series given by Eq (3.9) converges at κ=1 with appropriate pick of ϑ0(γ,τ). Thus, we have

    ϑ(γ,τ)=ϑ0(γ,τ)+μ=1ϑμ(γ,τ). (3.11)

    Equation (3.11) provides a relationship between ϑ0(γ,τ) and the exact solution ϑ(γ,τ) through the terms ϑμ(γ,τ),(μ=1,2,3,...), that will be calculated in upcoming steps. Eq (3.11) provides the solution of Eq (3.1) in the form of a series.

    The vectors are constituted as

    ϑμ={ϑ0(γ,τ),ϑ1(γ,τ),ϑ2(γ,τ),ϑ3(γ,τ),......ϑμ(γ,τ)}. (3.12)

    Now, the μth-order deformation equation is framed as

    LFNσ[ϑμ(γ,τ)χμϑμ1(γ,τ)]=μ[ϑμ1(γ,τ)]. (3.13)

    Operating the inverse of LFNT on Eq (3.13), we get

    ϑμ(γ,τ)=χμϑμ1(γ,τ)+LFN1σ(μ(ϑμ1(γ,τ))). (3.14)

    In Eq (3.14), the value of μ(ϑμ1(γ,τ)) can be written in a new look as

    μ(ϑμ1(γ,τ))=LFNσ[ϑμ1(γ,τ)](1χμ)(uqσsqσq1k=0s(qk1)σu(qk)σϑ(kσ)(γ,0)+uqσsqσ(LFNσ[ω(γ,τ)]))+uqσsqσLFNσ(Pσϑμ1(γ,τ)+Pμ1), (3.15)

    where the value of χμ is presented as

    χμ={0,μ11,μ>1. (3.16)

    In Eq (3.15), Pμ denotes homotopy polynomial suggested in [55] in functioning of LFHAM [43,46], and is formulated as

    Pμ=1Γ(μ)[μκμQσ(φ(γ,τ;κ))]κ=0, (3.17)

    where

    φ=φ0+κφ1+κ2φ2+κ3φ3. (3.18)

    Putting the value of μ(ϑμ1) from Eq (3.15) in Eq (3.14) transforms the Eq (3.14) as follows:

    ϑμ(γ,τ)=(χμ+)ϑμ1(γ,τ)(1χμ)LFN1σ(uqσsqσq1k=0s(qk1)σu(qk)σϑ(kσ)(γ,0)+uqσsqσ(LFNσ[ω(γ,τ)]))+LFN1σ(uqσsqσLFNσ(Pσϑμ1(γ,τ)+Pμ1)). (3.19)

    From Eq (3.19), the components ϑμ(γ,τ) can be evaluated for μ1 and the LFNHAM solution is presented in the following way:

    ϑ(γ,τ)=limNNμ=0ϑμ(γ,τ). (3.20)

    The significant aspect of the LFNHAM is the auxiliary parameter 0 which guarantees the convergence of the series solution of Eq (3.1).

    Theorem 3.1. If a constant 0<<1 can be estimated such that ϑμ+1(γ,τ)ϑμ(γ,τ) for each value of μ. Moreover, if the truncated series Nμ=0ϑμ(γ,τ) is assumed as an approximate solution ϑ then the maximum absolute truncated error is computed as

    ϑ(γ,τ)Nμ=0ϑμ(γ,τ)N+1(1)ϑ0(γ,τ). (3.21)

    Proof. The maximum absolute truncated error is computed in this way:

    ϑ(γ,τ)Nμ=0ϑμ(γ,τ)=μ=N+1ϑμ(γ,τ)μ=N+1ϑμ(γ,τ)μ=N+1μϑ0(γ,τ)()N+1[1++2+]ϑ0(γ,τ)N+1(1)ϑ0(γ,τ).

    This finishes the proof.

    In the upcoming Section 4, we establish the convergence and uniqueness of the LFNHAM solution.

    Theorem 4.1. (Uniqueness theorem). The attainment of solution by implementation of LFNHAM for the LFPDE (3.1) is unique, wherever 0<ρ<1, where

    ρ=(1+)+(λ+ξ)ς. (4.1)

    Proof. The solution of nonlinear LFPDE (3.1) is obtained as

    ϑ(γ,τ)=limNNμ=0ϑμ(γ,τ), (4.2)

    where

    ϑμ(γ,τ)=(χμ+)ϑμ1(γ,τ)(1χμ)LFN1σ(uqσsqσq1k=0s(qk1)σu(qk)σϑ(kσ)(γ,0)+uqσsqσ(LFNσ[ω(γ,τ)]))+LFN1σ(uqσsqσLFNσ(Pσϑμ1(γ,τ)+Pμ1)). (4.3)

    Let ϑ(γ,τ) and ϑ(γ,τ) be two distinct possible solutions for Eq (3.1), then we acquire

    |ϑ(γ,τ)ϑ(γ,τ)|=|(1+)(ϑϑ)+LFN1σ(uqσsqσLFNσ[Pσ(ϑϑ)+Qσ(ϑϑ)])|(1+)|ϑϑ|+LFN1σ(uqσsqσLFNσ|Pσ(ϑϑ)+Qσ(ϑϑ)|). (4.4)

    Utilization of the local fractional convolution theorem for LFNT in Eq (4.4) gives

    |ϑϑ|(1+)|ϑϑ|+Γ(1+σ)τ0|Pσ(ϑϑ)+Qσ(ϑϑ)|(τη)(q1)σΓ(1+(q1)σ)(dη)σ(1+)|ϑϑ|+Γ(1+σ)τ0(λ+ξ)|ϑϑ|(τη)(q1)σΓ(1+(q1)σ)(dη)σ. (4.5)

    Now, with the help of mean value theorem (MVT) of LFI calculus [56,57], inequality (4.5) transforms in the following form

    |ϑ(γ,τ)ϑ(γ,τ)|(1+)|ϑϑ|+(λ+ξ)|ϑϑ|ς=[(1+)+(λ+ξ)ς]|ϑϑ|=ρ|ϑϑ|.(1ρ)|ϑϑ|0, (4.6)

    where ρ=(1+)+(λ+ξ)ς. (4.7)

    Since 0<ρ<1, therefore |ϑϑ|=0, which provides ϑ=ϑ. This ensures the aspect of uniqueness of the solution of Eq (3.1).

    Theorem 4.2. (Convergence theorem). Suppose Ξ is a Banach space and Θ:ΞΞ is a nonlinear mapping. Assume that

    Θ(ω)Θ(ϑ)λωϑ,ϑ,ωΞ. (4.8)

    Then Banach's fixed point theory [58,59] suggests the existence of a fixed point forΘ. Moreover, the sequence constructed by LFNHAM converges to the fixed points of Θ with arbitrary choices of ω0,ϑ0Ξ and

    ϑl1ϑl2λl21λϑ1ϑ0,ϑ,ωΞ. (4.9)

    Proof. It is presumed that (Π[Ω],.), where Π[Ω] signifies the Banach space of continuous functions on real line interval Ω holding the sup norm. Now, it is sufficient to prove that {ϑl2} is a Cauchy sequence in the Banach space Ξ.

    Now, consider

    ϑl1ϑl2=maxτΩ|ϑl1ϑl2|=maxτΩ|(1+)(ϑl11ϑl21)+LFN1σ(uqσsqσLFNσ[Pσ(ϑl11ϑl21)+Qσ(ϑl11ϑl21)])|maxτΩ{(1+)|ϑl11ϑl21|+LFN1σ(uqσsqσLFNσ[|Pσ(ϑl11ϑl21)|+|Qσ(ϑl11ϑl21)|])}. (4.10)

    Employing the local fractional convolution theorem for LFNT in Eq (4.10), we have

    ϑl1ϑl2maxτΩ{(1+)|ϑl11ϑl21|+Γ(1+σ)τ0[|Pσ(ϑl11ϑl21)|+|Qσ(ϑl11ϑl21)|](τη)(q1)σΓ(1+(q1)σ)(dη)σ}maxτΩ{(1+)|ϑl11ϑl21|+Γ(1+σ)τ0(λ+ξ)|ϑl11ϑl21|(τη)(q1)σΓ(1+(q1)σ)(dη)σ}. (4.11)

    Now, application of MVT of LFI calculus [56,57] reduces the inequality (4.11) in the following form

    ϑl1ϑl2maxτΩ{(1+)|ϑl11ϑl21|+(λ+ξ)|ϑl11ϑl21|ς}=maxtΓ{(1+)+(λ+ξ)ς}|ϑl11ϑl21|=λϑl11ϑl21.ϑl1ϑl2λϑl11ϑl21, (4.12)
    whereλ=(1+)+(λ+ξ)ς (4.13)

    Assume l1=l2+1, then it produces

    ϑl2+1ϑl2λϑl2ϑl21λ2ϑl21ϑl22λl2ϑ1ϑ0. (4.14)

    Utilizing the triangular inequality, we have

    ϑl1ϑl2ϑl2+1ϑl2+ϑl2+2ϑl2+1++ϑl1ϑl11(λl2+λl2+1+λl2+2++λl11)ϑ1ϑ0=λl2(1+λ+λ2++λl1l21)ϑ1ϑ0=λl2[1λl1l211λ]ϑ1ϑ0. (4.15)

    Since 0<λ<1, thus 1λl1l21<1, then

    ϑl1ϑl2λl21λϑ1ϑ0. (4.16)

    But ϑ1ϑ0<, thus ϑl1ϑl20 as n, hence {ϑl2} is a Cauchy sequence in Π[Ω] and so {ϑl2} is convergent. This ensures the convergence of the solution ϑ(γ,τ) of LFPDE (3.1). Hence the theorem. □

    In this section, LFNHAM is implemented for deriving the solutions for LFCHEs.

    The following LFCHEs on Cantor set are investigated

    2σϑ1(γ,τ)τ2σ+2σϑ2(γ,τ)γ2σϑ1(γ,τ)=0,0<σ1, (5.1)
    2σϑ2(γ,τ)τ2σ+2σϑ1(γ,τ)γ2σϑ2(γ,τ)=0,0<σ1, (5.2)

    subject to the fractal initial conditions:

    ϑ1(γ,0)=0,σϑ1(γ,0)τσ=Eσ(γσ),
    ϑ2(γ,0)=0,σϑ2(γ,0)τσ=Eσ(γσ), (5.3)

    where ϑ1(γ,τ) and ϑ2(γ,τ) represent the local fractional continuous functions.

    On account of the initial conditions (5.3) and algorithm of LFNHAM, the initial guess are written as

    ϑ1,0(γ,τ)=ϑ1(γ,0)+τσΓ(1+σ)ϑ1,0(σ)(γ,0)=Eσ(γσ)τσΓ(1+σ),
    ϑ2,0(γ,τ)=ϑ2(γ,0)+τσΓ(1+σ)ϑ2,0(σ)(γ,0)=Eσ(γσ)τσΓ(1+σ). (5.4)

    Employing the LFNT operator LFNσ on Eqs (5.1) and (5.2), we get

    LFNσ[2σϑ1(γ,τ)τ2σ]+LFNσ[2σϑ2(γ,τ)γ2σ]LFNσ(ϑ1(γ,τ))=0, (5.5)
    LFNσ[2σϑ2(γ,τ)τ2σ]+LFNσ[2σϑ1(γ,τ)γ2σ]LFNσ(ϑ2(γ,τ))=0. (5.6)

    Now, the implementation of formula of LFNT for local fractional derivatives yields

    s2σu2σˉϑ1(γ,s,u)sσu2σϑ1(γ,0)1uσϑ1(σ)(γ,0)+LFNσ[2σϑ2(γ,τ)γ2σ]LFNσ(ϑ1(γ,τ))=0, (5.7)
    s2σu2σˉϑ2(γ,s,u)sσu2σϑ2(γ,0)1uσϑ2(σ)(γ,0)+LFNσ[2σϑ1(γ,τ)γ2σ]LFNσ(ϑ2(γ,τ))=0. (5.8)

    After rearranging the terms, we get

    ˉϑ1(γ,s,u)=1sσϑ1(γ,0)+uσs2σϑ1(σ)(γ,0)u2σs2σLFNσ[2σϑ2(γ,τ)γ2σ]+u2σs2σLFNσ(ϑ1(γ,τ)), (5.9)
    ˉϑ2(γ,s,u)=1sσϑ2(γ,0)+uσs2σϑ2(σ)(γ,0)u2σs2σLFNσ[2σϑ1(γ,τ)γ2σ]+u2σs2σLFNσ(ϑ2(γ,τ)). (5.10)

    Now, further simplification in view of initial condition (5.3) reduces Eqs (5.9) and (5.10) in the following way

    ˉϑ1(γ,s,u)uσs2σEσ(γσ)+u2σs2σLFNσ[2σϑ2(γ,τ)γ2σ]u2σs2σLFNσ(ϑ1)=0, (5.11)
    ˉϑ2(γ,s,u)+uσs2σEσ(γσ)+u2σs2σLFNσ[2σϑ1(γ,τ)γ2σ]u2σs2σLFNσ(ϑ2)=0. (5.12)

    Now in view of Eqs (5.11) and (5.12), the nonlinear operators are formed as:

    Φ1[φ1(γ,τ;κ)]=LFNσ(φ1(γ,τ;κ))uσs2σEσ(γσ)+u2σs2σLFNσ[2σφ2(γ,τ;κ)γ2σ]u2σs2σLFNσ(φ1(γ,τ;κ)), (5.13)
    Φ2[φ2(γ,τ;κ)]=LFNσ(φ1(γ,τ;κ))+uσs2σEσ(γσ)+u2σs2σLFNσ[2σφ1(γ,τ;κ)γ2σ]u2σs2σLFNσ(φ2(γ,τ;κ)), (5.14)

    where κ is an embedding parameter and φ1(γ,τ;κ) & φ2(γ,τ;κ) are real valued functions of γ,τ and κ.

    Now using the steps of the LFHAM [43,46] and basic methodology of HAM [52,53,54], μth-order deformation equations are constructed as follows:

    ϑ1,μ(γ,τ)=χμϑ1,(μ1)(γ,τ)+LFN1σ(μ(ϑ1,(μ1)(γ,τ))), (5.15)
    ϑ2,μ(γ,τ)=χμϑ2,(μ1)(γ,τ)+LFN1σ(μ(ϑ2,(μ1)(γ,τ))). (5.16)

    In Eqs (5.15) and (5.16), the terms μ(ϑ1,(μ1)(γ,τ)) and μ(ϑ2,(μ1)(γ,τ)) are expressed as

    μ(ϑ1,(μ1)(γ,τ))=LFNσ[ϑ1,(μ1)(γ,τ)](1χμ)uσs2σEσ(γσ)+u2σs2σLFNσ[2σϑ2,(μ1)(γ,τ)γ2σ]u2σs2σLFNσ(ϑ1,(μ1)(γ,τ)), (5.17)
    μ(ϑ2,(μ1)(γ,τ))=LFNσ[ϑ2,(μ1)(γ,τ)]+(1χμ)uσs2σEσ(γσ)+u2σs2σLFNσ[2σϑ1,(μ1)(γ,τ)γ2σ]u2σs2σLFNσ(ϑ2,(μ1)(γ,τ)). (5.18)

    Now implementing the LFNHAM and using Eqs (5.15) to (5.18), we have

    ϑ1,μ(γ,τ)=(χμ+)ϑ1,(μ1)(γ,τ)(1χμ)Eσ(γσ)τσΓ(1+σ)+LFN1σ(u2σs2σLFNσ[2σϑ2,(μ1)(γ,τ)γ2σ]u2σs2σLFNσ(ϑ1,(μ1)(γ,τ))),μ1, (5.19)
    ϑ2,μ(γ,τ)=(χμ+)ϑ2,(μ1)(γ,τ)+(1χμ)Eσ(γσ)τσΓ(1+σ)+LFN1σ(u2σs2σLFNσ[2σϑ1,(μ1)(γ,τ)γ2σ]u2σs2σLFNσ(ϑ2,(μ1)(γ,τ))),μ1. (5.20)

    On account of Eqs (5.19) and (5.20) for μ=1, we have

    ϑ1,1(γ,τ)=ϑ1,0(γ,τ)Eσ(γσ)τσΓ(1+σ)+LFN1σ(u2σs2σLFNσ[2σϑ2,0(γ,τ)γ2σ]u2σs2σLFNσ(ϑ1,0(γ,τ))), (5.21)
    ϑ2,1(γ,τ)=ϑ2,0(γ,τ)+Eσ(γσ)τσΓ(1+σ)+LFN1σ(u2σs2σLFNσ[2σϑ1,0(γ,τ)γ2σ]u2σs2σLFNσ(ϑ2,0(γ,τ))). (5.22)

    Using initial guess values (5.4) and further simplification reduces Eqs (5.21) and (5.22) in this way

    ϑ1,1(γ,τ)=2Eσ(γσ)LFN1σ(u2σs2σuσs2σ)=2Eσ(γσ)τ3σΓ(1+3σ), (5.23)
    ϑ2,1(γ,τ)=2Eσ(γσ)LFN1σ(u2σs2σuσs2σ)=2Eσ(γσ)τ3σΓ(1+3σ). (5.24)

    By means of Eqs (5.19) and (5.20) for μ=2, we have

    ϑ1,2(γ,τ)=(1+)ϑ1,1(γ,τ)+LFN1σ(u2σs2σLFNσ[2σϑ2,1(γ,τ)γ2σ]u2σs2σLFNσ(ϑ1,1(γ,τ))), (5.25)
    ϑ2,2(γ,τ)=(1+)ϑ2,1(γ,τ)+LFN1σ(u2σs2σLFNσ[2σϑ1,1(γ,τ)γ2σ]u2σs2σLFNσ(ϑ2,1(γ,τ))). (5.26)

    Utilizing the values provided by Eqs (5.23) and (5.24) in Eqs (5.25) and (5.26), we obtain

    ϑ1,2(γ,τ)=2(1+)Eσ(γσ)τ3σΓ(1+3σ)+42Eσ(γσ)τ5σΓ(1+5σ), (5.27)
    ϑ2,2(γ,τ)=2(1+)Eσ(γσ)τ3σΓ(1+3σ)42Eσ(γσ)τ5σΓ(1+5σ). (5.28)

    Following the similar procedure, we obtain the rest of the values for ϑ1,μ(γ,τ) and ϑ2,μ(γ,τ) for μ3.

    Setting the convergence-control parameter =1, we attain the following values

    ϑ1,1(γ,τ)=2Eσ(γσ)τ3σΓ(1+3σ),
    ϑ2,1(γ,τ)=2Eσ(γσ)τ3σΓ(1+3σ),
    ϑ1,2(γ,τ)=4Eσ(γσ)τ5σΓ(1+5σ),
    ϑ2,2(γ,τ)=4Eσ(γσ)τ5σΓ(1+5σ), (5.29)

    and so on.

    Hence, the solutions of Eqs (5.1) and (5.2) are obtained as

    ϑ1(γ,τ)=n=0ϑ1,n(γ,τ)=Eσ(γσ)(τσΓ(1+σ)+2τ3σΓ(1+3σ)+4τ5σΓ(1+5σ)+8τ7σΓ(1+7σ)+)=Eσ(γσ)(n=02nτ(2n+1)σΓ(1+(2n+1)σ))=12Eσ(γσ)(n=02n+12τ(2n+1)σΓ(1+(2n+1)σ))=12Eσ(γσ)sinhσ(2τσ). (5.30)

    Similarly, we have

    ϑ2(γ,τ)=n=0ϑ2,n(γ,τ)=Eσ(γσ)(τσΓ(1+σ)2τ3σΓ(1+3σ)4τ5σΓ(1+5σ)8τ7σΓ(1+7σ))=Eσ(γσ)(n=02nτ(2n+1)σΓ(1+(2n+1)σ))=12Eσ(γσ)sinhσ(2τσ). (5.31)

    Finally, the solutions of the coupled Helmholtz Eqs (5.1) & (5.2) are expressed as

    ϑ1(γ,τ)=Eσ(γσ)sinhσ(2τσ)2,
    ϑ2(γ,τ)=Eσ(γσ)sinhσ(2τσ)2. (5.32)

    Thus, the LFNHAM solutions are in complete agreement with the solutions obtained by Yang and Hua [60].

    In this portion, the LFNHAM is executed for deriving the solutions for LFCBEs.

    σϑ1τσ+2σϑ1γ2σ2σϑ1γσϑ1+σ[ϑ1ϑ2]γσ=0,ϑ1=ϑ1(γ,τ),0<σ1, (6.1)
    σϑ2τσ+2σϑ2γ2σ2σϑ2γσϑ2+σ[ϑ1ϑ2]γσ=0,ϑ2=ϑ2(γ,τ),0<σ1, (6.2)

    subject to the fractal initial conditions:

    ϑ1(γ,0)=ϑ2(γ,0)=Eσ(γσ), (6.3)

    where ϑ1(γ,τ) and ϑ2(γ,τ) are local fractional continuous functions.

    In view of initial conditions (6.3) and LFNHAM, the initial guess are expressed as

    ϑ1,0(γ,τ)=ϑ2,0(γ,τ)=Eσ(γσ). (6.4)

    Employing the LFNT operator LFNσ on Eqs (6.1) and (6.2), we get

    LFNσ(σϑ1τσ)+LFNσ(2σϑ1γ2σ)2LFNσ(σϑ1γσϑ1)+LFNσ(σ[ϑ1ϑ2]γσ)=0, (6.5)
    LFNσ(σϑ2τσ)+LFNσ(2σϑ2γ2σ)2LFNσ(σϑ2γσϑ2)+LFNσ(σ[ϑ1ϑ2]γσ)=0. (6.6)

    Now, employing the formula of LFNT for LFDs yields

    sσuσˉϑ1(γ,s,u)1uσϑ1(γ,0)+LFNσ(2σϑ1γ2σ)2LFNσ(σϑ1γσϑ1)+LFNσ(σ[ϑ1ϑ2]γσ)=0, (6.7)
    sσuσˉϑ2(γ,s,u)1uσϑ2(γ,0)+LFNσ(2σϑ2γ2σ)2LFNσ(σϑ2γσϑ2)+LFNσ(σ[ϑ1ϑ2]γσ)=0. (6.8)

    After rearranging the terms, we get

    ˉϑ1(γ,s,u)=1sσϑ1(γ,0)uσsσLFNσ(2σϑ1γ2σ)+2uσsσLFNσ(σϑ1γσϑ1)uσsσLFNσ(σ[ϑ1ϑ2]γσ), (6.9)
    ˉϑ2(γ,s,u)=1sσϑ2(γ,0)uσsσLFNσ(2σϑ2γ2σ)+2uσsσLFNσ(σϑ2γσϑ2)uσsσLFNσ(σ[ϑ1ϑ2]γσ). (6.10)

    Now, further simplification in view of initial conditions (6.3) reduces Eqs (6.9) and (6.10) as follows:

    ˉϑ1(γ,s,u)1sσEσ(γσ)+uσsσLFNσ(2σϑ1γ2σ)2uσsσLFNσ(σϑ1γσϑ1)+uσsσLFNσ(σ[ϑ1ϑ2]γσ)=0 (6.11)
    ˉϑ2(γ,s,u)1sσEσ(γσ)+uσsσLFNσ(2σϑ2γ2σ)2uσsσLFNσ(σϑ2γσϑ2)+uσsσLFNσ(σ[ϑ1ϑ2]γσ)=0. (6.12)

    Now in view of Eqs (6.11) and (6.12), the nonlinear operators are formed as

    Φ1[φ1(γ,τ;κ)]=LFNσ(φ1(γ,τ;κ))1sσEσ(γσ)+uσsσLFNσ[2σφ1(γ,τ;κ)γ2σ]2uσsσLFNσ(σφ1γσφ1)+uσsσLFNσ(σ[φ1(γ,τ;κ)φ2(γ,τ;κ)]γσ), (6.13)
    Φ2[φ2(γ,τ;κ)]=LFNσ(φ2(γ,τ;κ))1sσEσ(γσ)+uσsσLFNσ[2σφ2(γ,τ;κ)γ2σ]2uσsσLFNσ(σφ2γσφ2)+uσsσLFNσ(σ[φ1(γ,τ;κ)φ2(γ,τ;κ)]γσ), (6.14)

    where κ[0,1] is an embedding parameter and φ1(γ,τ;κ) & φ2(γ,τ;κ) are real valued functions of γ,τ, and κ.

    Performing the steps of the LFNHAM, μth-order deformation equations are formed as

    ϑ1,μ(γ,τ)=χμϑ1,(μ1)(γ,τ)+LFN1σ(μ(ϑ1,(μ1)(γ,τ))), (6.15)
    ϑ2,μ(γ,τ)=χμϑ2,(μ1)(γ,τ)+LFN1σ(μ(ϑ2,(μ1)(γ,τ))). (6.16)

    In Eqs (6.15) and (6.16), the terms μ(ϑ1,(μ1)(γ,τ)) and μ(ϑ2,(μ1)(γ,τ)) are computed as

    μ(ϑ1,(μ1)(γ,τ))=LFNσ[ϑ1,(μ1)(γ,τ)](1χμ)1sσEσ(γσ)+uσsσLFNσ[2σϑ1,(μ1)(γ,τ)γ2σ]2uσsσLFNσ(Pμ1)+uσsσLFNσ(σPμ1γσ), (6.17)
    μ(ϑ2,(μ1)(γ,τ))=LFNσ[ϑ2,(μ1)(γ,τ)](1χμ)1sσEσ(γσ)+uσsσLFNσ[2σϑ2,(μ1)(γ,τ)γ2σ]2uσsσLFNσ(Pμ1)+uσsσLFNσ(σPμ1γσ), (6.18)

    where Pμ1, Pμ1, Pμ1 denote the homotopy polynomials [55] and are formulated as

    Pμ1=1Γ(μ)[μκμ(σφ1(γ,τ;κ)γσφ1(γ,τ;κ))], (6.19)
    Pμ1=1Γ(μ)[μκμ(φ1(γ,τ;κ)φ2(γ,τ;κ))], (6.20)
    Pμ1=1Γ(μ)[μκμ(σφ2(γ,τ;κ)γσφ2(γ,τ;κ))], (6.21)

    and

    φ1(γ,t,κ)=φ1,0+κφ1,1+κ2φ1,2+, (6.22)
    φ2(γ,t,κ)=φ2,0+κφ2,1+κ2φ2,2+. (6.23)

    On account of LFNHAM and Eqs (6.15)–(6.18), we have

    ϑ1,μ(γ,τ)=(χμ+)ϑ1,(μ1)(γ,τ)(1χμ)Eσ(γσ)+LFN1σ(uσsσLFNσ[2σϑ1,(μ1)(γ,τ)γ2σ]2uσsσLFNσ(Pμ1)+uσsσLFNσ(σPμ1γσ)),μ1. (6.24)
    {\vartheta _{2, \mu }}\left( {\gamma , \, \tau } \right) = \left( {{\chi _\mu } + \hbar } \right){\vartheta _{2, \, (\mu - 1)}}\left( {\gamma , \, \tau } \right) - \hbar \, \left( {1 - {\chi _\mu }} \right)\, {E_\sigma }\left( {{\gamma ^\sigma }} \right) \\ \;\;\;\;\;\;\;\;\;\; + \hbar {}^{LF}N_\sigma ^{ - 1}\left( {\frac{{{u^\sigma }}}{{{s^\sigma }}}{}^{LF}{N_\sigma }\left[ {\frac{{{\partial ^{2\sigma }}{\vartheta _{2, \, (\mu - 1)}}\left( {\gamma , \tau } \right)}}{{\partial {\gamma ^{2\sigma }}}}} \right] - 2\frac{{{u^\sigma }}}{{{s^\sigma }}}{}^{LF}{N_\sigma }\left( {{{P''}_{\mu - 1}}} \right) + \frac{{{u^\sigma }}}{{{s^\sigma }}}{}^{LF}{N_\sigma }\left( {\frac{{{\partial ^\sigma }{{P'}_{\mu - 1}}}}{{\partial {\gamma ^\sigma }}}} \right)} \right) , \mu \geqslant 1 . (6.25)

    Taking into account the set of iterative schemes (6.24) & (6.25) and initial conditions (6.3), the iterative terms for various values of \mu are computed as follows:

    {\vartheta _{1, 1}}\left( {\gamma , \, \tau } \right) = \hbar \, {E_\sigma }\left( {{\gamma ^\sigma }} \right)\, \frac{{{\tau ^\sigma }}}{{\Gamma \left( {1 + \sigma } \right)}} , \\ {\vartheta _{2, 1}}\left( {\gamma , \, \tau } \right) = \hbar \, {E_\sigma }\left( {{\gamma ^\sigma }} \right)\, \frac{{{\tau ^\sigma }}}{{\Gamma \left( {1 + \sigma } \right)}} , \\ {\vartheta _{1, 2}}\left( {\gamma , \, \tau } \right) = \hbar \left( {1 + \hbar } \right)\, {E_\sigma }\left( {{\gamma ^\sigma }} \right)\, \frac{{{\tau ^\sigma }}}{{\Gamma \left( {1 + \sigma } \right)}} + {\hbar ^2}{E_\sigma }\left( {{\gamma ^\sigma }} \right)\frac{{{\tau ^{2\sigma }}}}{{\Gamma \left( {1 + 2\sigma } \right)}} , \\ {\vartheta _{2, 2}}\left( {\gamma , \, \tau } \right) = \hbar \left( {1 + \hbar } \right)\, {E_\sigma }\left( {{\gamma ^\sigma }} \right)\, \frac{{{\tau ^\sigma }}}{{\Gamma \left( {1 + \sigma } \right)}} + {\hbar ^2}{E_\sigma }\left( {{\gamma ^\sigma }} \right)\frac{{{\tau ^{2\sigma }}}}{{\Gamma \left( {1 + 2\sigma } \right)}} , \\{\vartheta _{1, 3}}\left( {\gamma , \, \tau } \right) = \hbar \left( {1 + \hbar } \right)\, {E_\sigma }\left( {{\gamma ^\sigma }} \right)\, \left( {\left( {1 + \hbar } \right)\frac{{{\tau ^\sigma }}}{{\Gamma \left( {1 + \sigma } \right)}} + \hbar \frac{{{\tau ^{2\sigma }}}}{{\Gamma \left( {1 + 2\sigma } \right)}}} \right)\\ \;\;\;\;\;\;\;\;\;\;+ {\hbar ^2}{E_\sigma }\left( {{\gamma ^\sigma }} \right)\, \left( {\left( {1 + \hbar } \right)\frac{{{\tau ^{2\sigma }}}}{{\Gamma \left( {1 + 2\sigma } \right)}} + \hbar \frac{{{\tau ^{3\sigma }}}}{{\Gamma \left( {1 + 3\sigma } \right)}}} \right) , \\ {\vartheta _{2, 3}}\left( {\gamma , \, \tau } \right) = \hbar \left( {1 + \hbar } \right)\, {E_\sigma }\left( {{\gamma ^\sigma }} \right)\, \left( {\left( {1 + \hbar } \right)\frac{{{\tau ^\sigma }}}{{\Gamma \left( {1 + \sigma } \right)}} + \hbar \frac{{{\tau ^{2\sigma }}}}{{\Gamma \left( {1 + 2\sigma } \right)}}} \right)\\ \;\;\;\;\;\;\;\;\;\;+ {\hbar ^2}{E_\sigma }\left( {{\gamma ^\sigma }} \right)\, \left( {\left( {1 + \hbar } \right)\frac{{{\tau ^{2\sigma }}}}{{\Gamma \left( {1 + 2\sigma } \right)}} + \hbar \frac{{{\tau ^{3\sigma }}}}{{\Gamma \left( {1 + 3\sigma } \right)}}} \right) . (6.26)

    Proceeding in the similar way, we get the remaining values for {\vartheta _{1, \mu }}\left({\gamma, \, \tau } \right) and {\vartheta _{2, \mu }}\left({\gamma, \, \tau } \right) for \mu \geqslant 4 .

    Setting the convergence-control parameter \hbar = - 1 , we attain the following values

    {\vartheta _{1, 1}}\left( {\gamma , \, \tau } \right) = - {E_\sigma }\left( {{\gamma ^\sigma }} \right)\, \frac{{{\tau ^\sigma }}}{{\Gamma \left( {1 + \sigma } \right)}} ,
    {\vartheta _{2, 1}}\left( {\gamma , \, \tau } \right) = - \, {E_\sigma }\left( {{\gamma ^\sigma }} \right)\, \frac{{{\tau ^\sigma }}}{{\Gamma \left( {1 + \sigma } \right)}} ,
    {\vartheta _{1, 2}}\left( {\gamma , \, \tau } \right) = {E_\sigma }\left( {{\gamma ^\sigma }} \right)\frac{{{\tau ^{2\sigma }}}}{{\Gamma \left( {1 + 2\sigma } \right)}} ,
    {\vartheta _{2, 2}}\left( {\gamma , \, \tau } \right) = {E_\sigma }\left( {{\gamma ^\sigma }} \right)\frac{{{\tau ^{2\sigma }}}}{{\Gamma \left( {1 + 2\sigma } \right)}} ,
    {\vartheta _{1, 3}}\left( {\gamma , \, \tau } \right) = - {E_\sigma }\left( {{\gamma ^\sigma }} \right)\, \left( {\frac{{{\tau ^{3\sigma }}}}{{\Gamma \left( {1 + 3\sigma } \right)}}} \right) ,
    {\vartheta _{2, 3}}\left( {\gamma , \, \tau } \right) = - {E_\sigma }\left( {{\gamma ^\sigma }} \right)\, \left( {\frac{{{\tau ^{3\sigma }}}}{{\Gamma \left( {1 + 3\sigma } \right)}}} \right) , (6.27)
    \vdots

    and so on.

    Proceeding in the same way, the rest of the terms of {\vartheta _{1, \mu }}\left({\gamma, \, \tau } \right) and {\vartheta _{2, \mu }}\left({\gamma, \, \tau } \right) for \mu \geqslant 4 are evaluated in a smooth manner, and finally the local fractional series solutions are obtained.

    Hence, the solutions of Eqs (6.1) and (6.2) are obtained as

    {\vartheta _1}\, \left( {\gamma , \tau } \right) = \sum\limits_{n = 0}^\infty {{\vartheta _{1, n}}\, \left( {\gamma , \tau } \right)}\\ \;\;\;\;\;\;\;\;\;\; = {E_\sigma }\left( {{\gamma ^\sigma }} \right)\, \left( {1 - \frac{{{\tau ^\sigma }}}{{\Gamma \left( {1 + \sigma } \right)}} + \frac{{{\tau ^{2\sigma }}}}{{\Gamma \left( {1 + 2\sigma } \right)}} - \frac{{{\tau ^{3\sigma }}}}{{\Gamma \left( {1 + 3\sigma } \right)}} + \frac{{\, {\tau ^{4\sigma }}}}{{\Gamma \left( {1 + 4\sigma } \right)}} - \cdots } \right) \\ \;\;\;\;\;\;\;\;\;\; = {E_\sigma }\left( {{\gamma ^\sigma }} \right)\, {E_\sigma }\left( { - {\tau ^\sigma }} \right) . (6.28)

    Similarly,

    {\vartheta _2}\, \left( {\gamma , \tau } \right) = \sum\limits_{n = 0}^\infty {{\vartheta _{2, n}}\, \left( {\gamma , \tau } \right) = \, } {E_\sigma }\left( {{\gamma ^\sigma }} \right)\, {E_\sigma }\left( { - {\tau ^\sigma }} \right) . (6.29)

    The solutions of LFCHEs and LFCBEs acquired by implementing the LFNHAM are the general form of solutions as compared to the LFADM, LFHPM, and LFNHPM. It is noteworthy that the LFNHAM solution transforms to the LFNHPM solution for \hbar = - 1 . The computational results validate the reliability and accuracy of the proposed method to achieve solutions for LFCHEs and LFCBEs. Moreover, the solutions of LFCHEs are in excellent match with the solutions obtained by Yang and Hua [60]. These facts authenticate the reliability of the solutions obtained by LFNHAM. Conclusively, the suggested hybrid framework can be employed to a variety of local fractional models occurring in a fractal medium.

    In this segment, the numerical simulations are presented for the solutions of the LFCHEs and LFCBEs under fractal initial conditions obtained via LFNHAM. The 3D plots for solutions of LFCHEs and LFCBEs have been generated for the fractal order \sigma = \ln 2/\ln 3 . Here, all the 3D plots on the Cantor sets have been prepared with the aid of MATLAB software. Figures 1 and 2 depict the 3D surface graphics of coupled solutions {\vartheta _1}\left({\gamma, \, \tau } \right) and {\vartheta _2}\left({\gamma, \, \tau } \right) for the LFCHEs. Similarly, Figures 3 and 4 show the 3D surface graphics of coupled solutions {\vartheta _1}\left({\gamma, \, \tau } \right) and {\vartheta _2}\left({\gamma, \, \tau } \right) for LFCBEs. Here, \gamma and \tau have been taken in the closed interval of 0 to 1. The nature of {\vartheta _1}\left({\gamma, \, \tau } \right) and {\vartheta _2}\left({\gamma, \, \tau } \right) have been explored with respect to \gamma and \tau . The fractal solutions of the LFCHEs and LFCBEs show interesting characteristics for \sigma = \ln 2/\ln 3 . The graphical presentation demonstrates that the computed solutions for the LFCHEs and LFCBEs are consistently dependent on the fractal order \sigma . Furthermore, the 3D figures drawn on Cantor sets indicate that the coupled solutions {\vartheta _1}\left({\gamma, \, \tau } \right) and {\vartheta _2}\left({\gamma, \, \tau } \right) are of fractal nature.

    Figure 1.  3D plot of the solution {\vartheta _1}\, \left({\gamma, \tau } \right) with respect to \gamma and \tau in case of LFCHEs for \sigma = \ln 2/\ln 3 .
    Figure 2.  3D plot of the solution {\vartheta _2}\, \left({\gamma, \tau } \right) with respect to \gamma and \tau in case of LFCHEs for \sigma = \ln 2/\ln 3 .
    Figure 3.  3D nature of the solution {\vartheta _1}\, \left({\gamma, \tau } \right) with respect to \gamma and \tau in case of LFCBEs for \sigma = \ln 2/\ln 3 .
    Figure 4.  3D behavior of the solution {\vartheta _2}\, \left({\gamma, \tau } \right) with respect to \gamma and \tau in case of LFCBEs for \sigma = \ln 2/\ln 3 .

    In this paper, the LFNHAM is proposed for computation of solutions for LFCHEs and LFCBEs on Cantor sets. The local fractional series solutions for LFCHEs and LFCBEs have been depicted in terms of Mittag-Leffler function. The 3D plots are presented for solutions of LFCHEs and LFCBEs by using the MATLAB software. It is clearly observed from the surface graphics of the solutions that the figures plotted on the Cantor set for the functions {\vartheta _1}\left({\gamma, \, \tau } \right) and {\vartheta _2}\left({\gamma, \, \tau } \right) are of fractal nature. The computational results authenticate the reliability and accuracy of the implemented method to obtain solutions for LFCHEs and LFCBEs. The combination of LFHAM and LFNT performs faster calculations than LFHAM. The convergence and uniqueness of the LFNHAM solution for a general LFPDE is also discussed in view of Banach's fixed point theory. In a nutshell, the suggested hybrid approach in connection with LFNT can be employed to such types of local fractional models appearing in a fractal media.

    The authors are very grateful to the referees for constructive comments and suggestions towards the improvement of this paper. This research received no external funding.

    The authors declare no conflict of interest.



    [1] J. Y. Liu, S. H. Tsai, C. C. Wang, C. R. Chu, Acoustic wave reflection from a rough seabed with a continuously varying sediment layer overlying an elastic basement, J. Sound Vib., 275 (2004), 739-755. https://doi.org/10.1016/j.jsv.2003.06.012 doi: 10.1016/j.jsv.2003.06.012
    [2] M. S. Samuel, A. Thomas, On fractional Helmholtz equations, Frac. Calc. Appl. Anal., 13 (2010), 295-308.
    [3] A. Prakash, M. Goyal, S. Gupta, Numerical simulation of space-fractional Helmholtz equation arising in seismic wave propagation, imaging and inversion, Pramana-J. Phys., 93 (2019), 28. https://doi.org/10.1007/s12043-019-1773-8 doi: 10.1007/s12043-019-1773-8
    [4] N. A. Shah, E. R. El-Zahar, M. D. Aljoufi, J. D. Chung, An efficient approach for solution of fractional-order Helmholtz equations, Adv. Differ. Equ., 2021 (2021), 14. https://doi.org/10.1186/s13662-020-03167-x doi: 10.1186/s13662-020-03167-x
    [5] J. M. Burgers, A mathematical model illustrating the theory of turbulence, Adv. Appl. Mech., 1 (1948), 171-199. https://doi.org/10.1016/S0065-2156(08)70100-5 doi: 10.1016/S0065-2156(08)70100-5
    [6] J. M. Burgers, The nonlinear diffusion equation, Dordrecht: Springer, 1974. https://doi.org/10.1007/978-94-010-1745-9
    [7] G. B. Whitham, Linear and nonlinear waves, New York: John Wiley & Sons, 2011. https://doi.org/10.1002/9781118032954
    [8] R. A. Oderinua, The reduced differential transform method for the exact solutions of advection, Burgers and coupled Burgers equations, Theory and Applications of Mathematics & Computer Science, 2 (2012), 10-14.
    [9] J. D. Cole, On a quasilinear parabolic equations occurring in aerodynamics, Q. Appl. Math., 9 (1951), 225-236.
    [10] S. E. Esipov, Coupled Burgers equations: a model of polydispersive sedimentation, Phys. Rev. E, 52 (1995), 3711-3718. https://doi.org/10.1103/physreve.52.3711 doi: 10.1103/physreve.52.3711
    [11] J. Nee, J. Duan, Limit set of trajectories of the coupled viscous Burger's equations, Appl. Math. Lett., 11 (1998), 57-61. https://doi.org/10.1016/S0893-9659(97)00133-X doi: 10.1016/S0893-9659(97)00133-X
    [12] N. Sugimoto, Burgers equation with a fractional derivative; hereditary effects on nonlinear acoustic waves, J. Fluid Mech., 225 (1991), 631-653. https://doi.org/10.1017/S0022112091002203 doi: 10.1017/S0022112091002203
    [13] T. Funaki, W. A. Woyczynski, Fractal Burgers equations, J. Differ. Equ., 148 (1998), 9-46. https://doi.org/10.1006/jdeq.1998.3458 doi: 10.1006/jdeq.1998.3458
    [14] S. Momani, Non-perturbative analytical solutions of the space-and time-fractional Burgers equations, Chaos Soliton. Fract., 28 (2006), 930-937. https://doi.org/10.1016/j.chaos.2005.09.002 doi: 10.1016/j.chaos.2005.09.002
    [15] C. H. Chan, M. Czubak, L. Silvestre, Eventual regularization of the slightly supercritical fractional Burgers equation, Discrete Contin. Dyn. Syst., 27 (2010), 847-861. https://doi.org/10.3934/dcds.2010.27.847 doi: 10.3934/dcds.2010.27.847
    [16] G. Karch, C. Miao, X. Xu, On convergence of solutions of fractal Burgers equation toward rarefaction waves, SIAM J. Math. Anal., 39 (2008), 1536-1549. https://doi.org/10.1137/070681776 doi: 10.1137/070681776
    [17] A. Yildirim, S. T. Mohyud-Din, Analytical approach to space-and time-fractional Burgers equations, Chin. Phys. Letts., 27 (2010), 090501. https://doi.org/10.1088/0256-307X/27/9/090501 doi: 10.1088/0256-307X/27/9/090501
    [18] N. Alibaud, C. Imbert, G. Karch, Asymptotic properties of entropy solutions to fractal Burgers equation, SIAM J. Math. Anal., 42 (2010), 354-376. https://doi.org/10.1137/090753449 doi: 10.1137/090753449
    [19] A. Yildirim, A. Kelleci, Homotopy perturbation method for numerical solutions of coupled Burgers equations with time- and space-fractional derivatives, Int. J. Numer. Methods Heat Fluid Flow, 20 (2010), 897-909. https://doi.org/10.1108/09615531011081423 doi: 10.1108/09615531011081423
    [20] J. Liu, G. Hou, Numerical solutions of the space- and time fractional coupled Burgers equations by generalized differential transform method, Appl. Math. Comput., 217 (2011), 7001-7008. https://doi.org/10.1016/j.amc.2011.01.111 doi: 10.1016/j.amc.2011.01.111
    [21] Y. Xu, O. P. Agrawal, Numerical solutions and analysis of diffusion for new generalized fractional Burgers equation, Fract. Calc. Appl. Anal., 16 (2013), 709-736. https://doi.org/10.2478/s13540-013-0045-4 doi: 10.2478/s13540-013-0045-4
    [22] A. Prakash, M. Kumar, K. K. Sharma, Numerical method for solving fractional coupled Burgers equations, Appl. Math. Comput., 260 (2015), 314-320. https://doi.org/10.1016/j.amc.2015.03.037 doi: 10.1016/j.amc.2015.03.037
    [23] J. Singh, D. Kumar, R. Swroop, Numerical solution of time- and space-fractional coupled Burgers' equations via homotopy algorithm, Alex. Eng. J., 55 (2016), 1753-1763. https://doi.org/10.1016/j.aej.2016.03.028 doi: 10.1016/j.aej.2016.03.028
    [24] X. J. Yang, Local fractional functional analysis and its applications, Hong Kong: Asian Academic Publisher, 2011.
    [25] X. J. Yang, Advanced local fractional calculus and its applications, New York: World Science Publisher, 2012.
    [26] X. J. Yang, H. M. Srivastava, C. Cattani, Local fractional homotopy perturbation method for solving fractional partial differential equations arising in mathematical physics, Rom. Rep. Phys., 67 (2015), 752-761.
    [27] Y. Zhang, C. Cattani, X. J. Yang, Local fractional homotopy perturbation method for solving non-homogeneous heat conduction equations in fractal domains, Entropy, 17 (2015), 6753-6764. https://doi.org/10.3390/e17106753 doi: 10.3390/e17106753
    [28] J. Singh, D. Kumar, J. J. Nieto, A reliable algorithm for a local fractional Tricomi equation arising in fractal transonic flow, Entropy, 18 (2016), 206. https://doi.org/10.3390/e18060206 doi: 10.3390/e18060206
    [29] D. Kumar, J. Singh, D. Baleanu, A hybrid computational approach for Klein-Gordon equations on Cantor Sets, Nonlinear Dyn., 87 (2017), 511-517. https://doi.org/10.1007/s11071-016-3057-x doi: 10.1007/s11071-016-3057-x
    [30] D. Zhao, J. Singh, D. Kumar, S. Rathore, X. J. Yang, An efficient computational technique for local fractional heat conduction equations in fractal media, J. Nonlinear Sci. Appl., 10 (2017), 1478-1486. http://dx.doi.org/10.22436/jnsa.010.04.17
    [31] J. Singh, D. Kumar, D. Baleanu, S. Rathore, On the local fractional wave equation in fractal strings, Math. Method. Appl. Sci., 42 (2019), 1588-1595. https://doi.org/10.1002/mma.5458 doi: 10.1002/mma.5458
    [32] V. P. Dubey, J. Singh, A. M. Alshehri, S. Dubey, D. Kumar, A comparative analysis of two computational schemes for solving local fractional Laplace equations, Math. Method. Appl. Sci., 44 (2021), 13540-13559. https://doi.org/10.1002/mma.7642 doi: 10.1002/mma.7642
    [33] H. Jafari, H. K. Jassim, A new approach for solving system of local fractional partial differential equations, Appl. Appl. Math., 11 (2016), 9.
    [34] D. Kumar, F. Tchier, J. Singh, D. Baleanu, An efficient computational technique for fractal vehicular traffic flow, Entropy, 20 (2018), 259. https://doi.org/10.3390/e20040259 doi: 10.3390/e20040259
    [35] Y. J. Hao, H. M. Srivastava, H. Jafari, X. J. Yang, Helmholtz and diffusion equations associated with local fractional derivative operators involving the Cantorian and Cantor-type cylindrical coordinates, Adv. Math. Phys., 2013 (2013), 754248. https://doi.org/10.1155/2013/754248 doi: 10.1155/2013/754248
    [36] X. J. Wang, Y. Zhao, C. Cattani, X. J. Yang, Local fractional variational iteration method for inhomogeneous Helmholtz equation within local fractional derivative operator, Math. Probl. Eng., 2014 (2014), 913202. https://doi.org/10.1155/2014/913202 doi: 10.1155/2014/913202
    [37] A. M. Yang, Z. S. Chen, H. M. Srivastava, X. J. Yang, Application of the local fractional series expansion method and the variational iteration method to the Helmholtz equation involving local fractional derivative operators, Abstr. Appl. Anal., 2013 (2013), 259125. https://doi.org/10.1155/2013/259125 doi: 10.1155/2013/259125
    [38] H. K. Jassim, The approximate solutions of Helmholtz and coupled Helmholtz equations on cantor sets within local fractional operator, J. Zankoy Sulaimani Part A, 17 (2015), 19-25. https://doi.org/10.17656/jzs.10421 doi: 10.17656/jzs.10421
    [39] D. Baleanu, H. K. Jassim, A modification fractional homotopy perturbation method for solving Helmholtz and coupled Helmholtz equations on cantor sets, Fractal Fract., 3 (2019), 30. https://doi.org/10.3390/fractalfract3020030 doi: 10.3390/fractalfract3020030
    [40] D. Baleanu, H. K. Jassim, M. Al Qurashi, Solving Helmholtz equation with local fractional derivative operators, Fractal Fract., 3 (2019), 43. https://doi.org/10.3390/fractalfract3030043 doi: 10.3390/fractalfract3030043
    [41] M. Chen, Y. Xia, Approximate solution for Burger's equation with local fractional derivative by Yang-Laplace decomposition method, Therm. Sci., 21 (2017), 209-216. https://doi.org/10.2298/TSCI17S1209C doi: 10.2298/TSCI17S1209C
    [42] X. J. Yang, J. A. T. Machado, J. Hristov, Nonlinear dynamics for local fractional Burger's equations arising in fractal flow, Nonlinear Dyn., 84 (2016), 3-7. https://doi.org/10.1007/s11071-015-2085-2 doi: 10.1007/s11071-015-2085-2
    [43] D. Ziane, A. Bokhari, R. Belgacem, Local fractional homotopy analysis method for solving coupled nonlinear systems of Burger's equations, Int. J. Open Problems Comput. Math., 12 (2019), 47-57.
    [44] S. O. Edeki, G. O. Akinlabi, Local fractional operator for a one-dimensional coupled Burger equation of non-integer time order parameter, J. Math. Fund. Sci., 50 (2018), 28-39. https://doi.org/10.5614/j.math.fund.sci.2018.50.1.3 doi: 10.5614/j.math.fund.sci.2018.50.1.3
    [45] X. J. Yang, F. Gao, H. M. Srivastava, Exact travelling wave solutions for the local fractional two-dimensional Burgers' type equations, Comput. Math. Appl., 73 (2017), 203-210. https://doi.org/10.1016/j.camwa.2016.11.012 doi: 10.1016/j.camwa.2016.11.012
    [46] S. Maitama, W. Zhao, Local fractional homotopy analysis method for solving nondifferentiable problems on Cantor sets, Adv. Differ. Equ., 2019 (2019), 127. https://doi.org/10.1186/s13662-019-2068-6 doi: 10.1186/s13662-019-2068-6
    [47] S. Maitama, Local fractional natural homotopy perturbation method for solving partial differential equations with local fractional derivative, Prog. Fract. Differ. Appl., 4 (2018), 219-228. http://dx.doi.org/10.18576/pfda/040306
    [48] Z. H. Khan, W. A. Khan, N-transform properties and applications, NUST Journal of Engineering Sciences, 1 (2008), 127-133.
    [49] F. B. M. Belgacem, R. Silambarasan, Theory of natural transform, Mathematics in Engineering, Science and Aerospace, 3 (2012), 105-135.
    [50] F. B. M. Belgacem, R. Silambarasan, Advances in the natural transform, AIP Conference Proceedings, 1493 (2012), 106-110. https://doi.org/10.1063/1.4765477 doi: 10.1063/1.4765477
    [51] R. Silambarasan, F. B. M. Belgacem, Application of the natural transform to Maxwell's equations, In: Progress in Electromagnetics Research Symposium Proceedings (PIERS '11), Suzhou, China, 2011,899-902.
    [52] S. J. Liao, Beyond perturbation: introduction to Homotopy analysis method, Boca Raton: Chapman and Hall/CRC Press, 2003. https://doi.org/10.1201/9780203491164
    [53] S. J. Liao, On the homotopy analysis method for nonlinear problems, Appl. Math. Comput., 147 (2004), 499-513. https://doi.org/10.1016/S0096-3003(02)00790-7 doi: 10.1016/S0096-3003(02)00790-7
    [54] S. J. Liao, Proposed homotopy analysis techniques for the solution of nonlinear problems, Ph.D. Thesis, Shanghai Jiao Tong University, Shanghai, 1992.
    [55] Z. Odibat, S. A. Bataineh, An adaptation of homotopy analysis method for reliable treatment of strongly nonlinear problems: construction of homotopy polynomials, Math. Method. Appl. Sci., 38 (2014), 991-1000. https://doi.org/10.1002/mma.3136 doi: 10.1002/mma.3136
    [56] X. J. Yang, A short note on local fractional calculus of function of one variable, J. Appl. Libr. Inf. Sci., 1 (2012), 1-13.
    [57] G. S. Chen, Mean value theorems for local fractional integrals on fractal space, Adv. Mech. Eng. Appl., 1 (2012), 5-8.
    [58] I. K. Argyros, Convergence and applications of Newton-type iterations, New York: Springer, 2008. https://doi.org/10.1007/978-0-387-72743-1
    [59] A. A. Magrenan, A new tool to study real dynamics: the convergence plane, Appl. Math. Comput., 248 (2014), 215-224. https://doi.org/10.1016/j.amc.2014.09.061 doi: 10.1016/j.amc.2014.09.061
    [60] Y. J. Yang, L. Q. Hua, Variational iteration transform method for fractional differential equations with local fractional derivative, Abstr. Appl. Anal., 2014 (2014), 760957. https://doi.org/10.1155/2014/760957 doi: 10.1155/2014/760957
  • This article has been cited by:

    1. Ved Prakash Dubey, Devendra Kumar, Jagdev Singh, Ahmed M. Alshehri, Sarvesh Dubey, Analysis of local fractional Klein-Gordon equations arising in relativistic fractal quantum mechanics, 2022, 1745-5030, 1, 10.1080/17455030.2022.2112993
    2. Saud Fahad Aldosary, Ram Swroop, Jagdev Singh, Ateq Alsaadi, Kottakkaran Sooppy Nisar, An efficient hybridization scheme for time-fractional Cauchy equations with convergence analysis, 2023, 8, 2473-6988, 1427, 10.3934/math.2023072
    3. Hasib Khan, Jehad Alzabut, Dumitru Baleanu, Ghada Alobaidi, Mutti-Ur Rehman, Existence of solutions and a numerical scheme for a generalized hybrid class of n-coupled modified ABC-fractional differential equations with an application, 2023, 8, 2473-6988, 6609, 10.3934/math.2023334
    4. Ghazala Akram, Saima Arshed, Maasoomah Sadaf, Maham Maqbool, Comparison of fractional effects for Phi-4 equation using beta and M-truncated derivatives, 2023, 55, 0306-8919, 10.1007/s11082-023-04549-7
    5. Ved Prakash Dubey, Devendra Kumar, Hashim M. Alshehri, Sarvesh Dubey, Jagdev Singh, Computational Analysis of Local Fractional LWR Model Occurring in a Fractal Vehicular Traffic Flow, 2022, 6, 2504-3110, 426, 10.3390/fractalfract6080426
    6. Ved Prakash Dubey, Jagdev Singh, Ahmed M. Alshehri, Sarvesh Dubey, Devendra Kumar, Analysis and Fractal Dynamics of Local Fractional Partial Differential Equations Occurring in Physical Sciences, 2023, 18, 1555-1415, 10.1115/1.4056360
    7. Devendra Kumar, Ved Prakash Dubey, Sarvesh Dubey, Jagdev Singh, Ahmed Mohammed Alshehri, Computational analysis of local fractional partial differential equations in realm of fractal calculus, 2023, 167, 09600779, 113009, 10.1016/j.chaos.2022.113009
    8. Hasib Khan, Jehad Alzabut, Haseena Gulzar, Existence of solutions for hybrid modified ABC-fractional differential equations with p-Laplacian operator and an application to a waterborne disease model, 2023, 70, 11100168, 665, 10.1016/j.aej.2023.02.045
    9. Nguyen Tuan, Devendra Kumar, Anh Nguyen, Regularity for a non-local diffusion equation with Riemann-Liouville derivative, 2023, 27, 0354-9836, 373, 10.2298/TSCI23S1373T
    10. Shruti Tomar, Naresh M. Chadha, Santanu Raut, 2023, Chapter 11, 978-3-031-29958-2, 177, 10.1007/978-3-031-29959-9_11
    11. Bhawna Pokhriyal, Pranay Goswami, A generalized local fractional LWR model of vehicular traffic flow and its solution, 2023, 46, 0170-4214, 18899, 10.1002/mma.9598
    12. Jagdev Singh, Ved Prakash Dubey, Devendra Kumar, Sarvesh Dubey, Dumitru Baleanu, Fractal-view analysis of local fractional Fokker–Planck equation occurring in modelling of particle’s Brownian motion, 2024, 56, 1572-817X, 10.1007/s11082-024-06842-5
    13. Usama Ghafoor, Muhammad Abbas, Tayyaba Akram, Emad K. El-Shewy, Mahmoud A. E. Abdelrahman, Noura F. Abdo, An Efficient Cubic B-Spline Technique for Solving the Time Fractional Coupled Viscous Burgers Equation, 2024, 8, 2504-3110, 93, 10.3390/fractalfract8020093
  • 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(2118) PDF downloads(91) Cited by(13)

Figures and Tables

Figures(4)  /  Tables(1)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog