Loading [MathJax]/jax/output/SVG/jax.js
Research article Special Issues

Travelling wave solutions for fractional Korteweg-de Vries equations via an approximate-analytical method

  • This paper introduces an approximate-analytical method (AAM) for solving nonlinear fractional partial differential equations (NFPDEs) in full general forms. The main advantage of the paper is to apply the proposed AAM to solve the fractional Korteweg-de Vries (KdV) equations. Moreover, the analytical travelling wave solutions for the fractional KdV equation and the modified fractional KdV equation are successfully obtained. The numerical solutions are also obtained in the forms of tables and graphs. The fractional partial derivatives are considered in Caputo sense.

    Citation: Hayman Thabet, Subhash Kendre, James Peters. Travelling wave solutions for fractional Korteweg-de Vries equations via an approximate-analytical method[J]. AIMS Mathematics, 2019, 4(4): 1203-1222. doi: 10.3934/math.2019.4.1203

    Related Papers:

    [1] Musawa Yahya Almusawa, Hassan Almusawa . Numerical analysis of the fractional nonlinear waves of fifth-order KdV and Kawahara equations under Caputo operator. AIMS Mathematics, 2024, 9(11): 31898-31925. doi: 10.3934/math.20241533
    [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] Musong Gu, Chen Peng, Zhao Li . Traveling wave solution of (3+1)-dimensional negative-order KdV-Calogero-Bogoyavlenskii-Schiff equation. AIMS Mathematics, 2024, 9(3): 6699-6708. doi: 10.3934/math.2024326
    [4] 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
    [5] Ihsan Ullah, Aman Ullah, Shabir Ahmad, Hijaz Ahmad, Taher A. Nofal . A survey of KdV-CDG equations via nonsingular fractional operators. AIMS Mathematics, 2023, 8(8): 18964-18981. doi: 10.3934/math.2023966
    [6] Ali Khalouta, Abdelouahab Kadem . A new computational for approximate analytical solutions of nonlinear time-fractional wave-like equations with variable coefficients. AIMS Mathematics, 2020, 5(1): 1-14. doi: 10.3934/math.2020001
    [7] M. Mossa Al-Sawalha, Rasool Shah, Adnan Khan, Osama Y. Ababneh, Thongchai Botmart . Fractional view analysis of Kersten-Krasil'shchik coupled KdV-mKdV systems with non-singular kernel derivatives. AIMS Mathematics, 2022, 7(10): 18334-18359. doi: 10.3934/math.20221010
    [8] Khalid Khan, Amir Ali, Muhammad Irfan, Zareen A. Khan . Solitary wave solutions in time-fractional Korteweg-de Vries equations with power law kernel. AIMS Mathematics, 2023, 8(1): 792-814. doi: 10.3934/math.2023039
    [9] Qasem M. Tawhari . Advanced analytical techniques for fractional Schrödinger and Korteweg-de Vries equations. AIMS Mathematics, 2025, 10(5): 11708-11731. doi: 10.3934/math.2025530
    [10] Yanhua Gu . High-order numerical method for the fractional Korteweg-de Vries equation using the discontinuous Galerkin method. AIMS Mathematics, 2025, 10(1): 1367-1383. doi: 10.3934/math.2025063
  • This paper introduces an approximate-analytical method (AAM) for solving nonlinear fractional partial differential equations (NFPDEs) in full general forms. The main advantage of the paper is to apply the proposed AAM to solve the fractional Korteweg-de Vries (KdV) equations. Moreover, the analytical travelling wave solutions for the fractional KdV equation and the modified fractional KdV equation are successfully obtained. The numerical solutions are also obtained in the forms of tables and graphs. The fractional partial derivatives are considered in Caputo sense.


    Fractional partial differential equations (FPDEs) have emerged from numerous research fields of science, engineering and physics, see [1,3,7,8,15,16,23]. There are many works devoted to the development of efficient methods for solving FPDEs but we should note that finding an analytical or approximate solution is a challenging problem, therefore, accurate methods for finding the solutions of FPDEs are still under investigation (see, e.g., new analytical technique for solving a system of nonlinear fractional partial differential equations [21], analytical solutions for conformable space-time conformable fractional partial differential equations [18], modified least squares homotopy perturbation method for solving FPDEs [19], conformable natural extension of the theory of conformable fractional calculus [22]). Several analytical and numerical methods for solving FPDEs exist in the literature, for example; the Generalized differential transform method[12] and the method was proposed to solve a space- and time-fractional diffusion-wave equation, homotopy perturbation method [6] for fractional Fornberg–Whitham equation, a homotopy perturbation technique [4], and so on. There are, however, a few solution methods for only travelling wave solutions, for example; the multiple exp-function algorithm [10]), fractional complex transformation [25] where the method was employed to solve fractional partial differential equations in the sense of the modified Riemann–Liouville derivative, Lie symmetries, symmetry reductions and conservation laws [2], new modification of Adomian decomposition method for solving a system of nonlinear fractional partial differential equations Research [20] and some references cited therein.

    As a modern mathematical discipline nonlinear wave theory frequently explores asymptotic regimes (such as varying over many scales, high frequency or large amplitude) which are not easily accessible via numerical simulations. Moreover, the nonlinear wave theory plays an important role to the understanding of real water waves, the interaction of light with matter, optical fiber transmission, traffic flow, earthquakes, galaxy formation and the steepening of short gravity waves on the crests of long waves. In this paper, two types of the well-known nonlinear wave equations with a fractional order are studied; one is the nonlinear KdV equation and another one is the modified KdV equation (see [2,13,17,24]).

    Briefly, a nonlinear wave equation of Korteweg de Vries type has the following form:

    ut+auux+bu2ux+uxxx=0 (nonlinear KdV),ut+6auux+6bu2ux+uxxx=0 (modified KdV),

    where the subscripts denote partial derivatives and parameters a,b can be scaled to any real number.

    The goal of this paper is to apply an AAM to solve full general NFPDEs without any assumption and that it gives good results in analytical and numerical experiments. The rest of the paper is organized in as follows: In Section 2, we present basic definitions and preliminaries which are needed in the sequel. In Section 3, we introduce an AAM for solving full general NFPDEs. Approximate analytical and numerical solutions for fractional KdV equations are obtained in Section 4.

    There are various definitions and properties of fractional integrals and derivatives. In this section, we present modifications of some basic definitions and preliminaries of the fractional calculus theory, which are used in this paper and can be found in [5,9,11,14] and some references cited therein.

    Definition 1. Let qRN and q0. Then, the Riemann-Liouville fractional partial integral denoted by Iqt of order q for a function u(x,t)Cμ, μ>1 is defined as:

    {Iqtu(x,t)=1Γ(q)t0(tτ)q1u(x,τ)dτ,q,t>0,I0tu(x,t)=u(x,t),q=0,t>0, (2.1)

    where Γ is the well-known Gamma function.

    Theorem 2.1. Let q1,q2RN,q1,q20 and p>1. Then, for a function u(x,t)Cμ, μ>1, the operator Iqt satisfies the following properties:

    {Iq1tIq2tu(x,t)=Iq1+q2tu(x,t),Iq1tIq2tu(x,t)=Iq2tIq1tu(x,t),Iqttp=Γ(p+1)Γ(p+q+1)tp+q. (2.2)

    Definition 2. Let q,tR,t>0,m1<q<mN and u(x,t)Cmμ. Then, the Caputo fractional partial derivative of order q for a function u(x,t) denoted by Dqtu(x,t) is defined as follows:

    {Dqtu(x,t)=ta(tτ)mq1Γ(mq)mu(x,τ)τmdτ,Dqtu(x,t)=mu(x,t)tm,q=mN. (2.3)

    Theorem 2.2. Let t,qR,t>0 and m1<q<mN. Then

    {IqtDqtu(x,t)=u(x,t)m1k=0tkk!ku(x,0+)tk,DqtIqtu(x,t)=u(x,t). (2.4)

    In this section, we present a AAM to solve NFPDEs with initial values of the following form:

    {Dqtu(ˉx,t)=f(ˉx,t)+Lˉu+Nˉu,m1<q<mN,iuti(ˉx,0)=fi(ˉx),i=0,1,2,,m1, (3.1)

    where L and N are linear and nonlinear operators, respectively, of u(ˉx,t) and its partial derivatives, which might include other fractional partial derivatives of orders less than q; f(ˉx,t) is known analytic function; and Dqt is the Caputo partial derivative of fractional order q, and ˉx=(x1,x2,,xn)Rn. The proposed AAM has much more computational power in obtaining piecewise analytical solutions.

    To introduce the AAM, first we need to present the following results.

    Lemma 3.1. For u(ˉx,t)=k=0pkuk(ˉx,t), the linear operator L(u) satisfies the following property:

    Lu(ˉx,t)=L(k=0pkuk(ˉx,t))=k=0pkL(uk(ˉx,t)). (3.2)

    Theorem 3.2. Let u(x,t)=k=0uk(ˉx,t), for the parameter λ, we define uλ(ˉx,t)=k=0λkuk(ˉx,t), then the nonlinear operator N(uλ) satisfies the following property:

    N(uλ)=N(k=0λkuk)=n=0[1n!nλn[N(nk=0λkuk)]λ=0]λn. (3.3)

    Proof By the help of the Maclaurin expansion of N(k=0λkuk) with respect to λ, we have

    N(uλ)=N(k=0λkuk)=[N(k=0λkuk)]λ=0+[λ[N(k=0λkuk)]λ=0]λ+[12!2λ2[N(k=0λkuk)]λ=0]λ2+=n=0[1n!nλn[N(k=0λkuk)]λ=0]λn=n=0[1n!nλn[N(nk=0λkuk+k=n+1λkuk)]λ=0]λn=n=0[1n!nλn[N(nk=0λkuk)]λ=0]λn.

    Definition 3. The polynomials En(u0,u1,,un) can be defined as:

    En(u0,u1,,un)=1n!nλn[N(nk=0λkuk)]λ=0. (3.4)

    Remark 3.1. Let En=En(u0,u1,,un) be as in Definition 3. Then by using Theorem 3.2, the nonlinear operator N(uλ) can be expressed in terms of En as:

    N(uλ)=n=0λnEn. (3.5)

    The following theorem introduces the approximate analytical solution obtained by the AAM for the nonlinear fractional partial differential equation with initial value given by (3.1).

    Theorem 3.3. Assume that f(ˉx,t),fi(ˉx) are defined as in (3.1) and for m1<q<mN. Then, equation (3.1) admits at least a solution given by

    u(ˉx,t)=f(q)t(ˉx,t)+m1i=0tii!fi(ˉx)+k=1[L(q)t(u(k1))+E(q)(k1)t], (3.6)

    where L(q)t(u(k1)) and E(q)(k1)t denote the fractional partial integral of order q for L(uk1) and E(k1) respectively with respect to t.

    Proof. Let the solution function u(ˉx,t) of equation (3.1) be as in the following analytical expansion:

    u(ˉx,t)=k=0uk(ˉx,t). (3.7)

    To solve the initial value problem (3.1), we consider

    Dqtuλ(ˉx,t)=λ[f(ˉx,t)+L(uλ)+N(uλ)],λ[0,1]. (3.8)

    with initial conditions given by

    iuλ(ˉx,0)ti=gi(ˉx),i=0,1,2,,m1. (3.9)

    Next, we assume that, the initial value problem (3.8) has a solution given by

    uλ(ˉx,t)=k=0λku(ˉx,t). (3.10)

    By applying Riemann-Liouville fractional partial integral of order q with respect to t to both sides of the initial value problem (3.8) and using Theorem 2.2, we obtain

    uλ(ˉx,t)=m1i=0tii!iuλ(ˉx,0)ti+λIqt[f(ˉx,t)+L(uλ)+N(uλ)]. (3.11)

    By using the initial condition given by (3.1), equation (3.11) can be rewritten as:

    uλ(ˉx,t)=m1i=0tii!gi(ˉx)+λ[f(q)t(ˉx,t)+Iqt[L(uλ)]+Iqt[N(uλ)]]. (3.12)

    By substituting (3.10) into (3.12), we obtain

    k=0λkuk(ˉx,t)=m1i=0tii!gi(ˉx)+λ[f(q)t(ˉx,t)+Iqt[L(k=0λkuk)]+Iqt[N(k=0λkuk)]]. (3.13)

    By using Lemma 3.1 and Theorem 3.2, equation (3.13) becomes

    k=0λku(ˉx,t)=m1i=0tii!gi(ˉx)+λf(q)t(ˉx,t)+Iqtλk=0[L(λkuk)]+Iqtλn=0[1n!nλn[N(nk=0λkuk)]λ=0]λn. (3.14)

    Next, we use Definition 3 in equation (3.14), we obtain

    k=0λku(ˉx,t)=m1i=0tii!gi(ˉx)+λf(q)t(ˉx,t)+Iqtλk=0[L(λkuk)]+Iqtλn=0Enλn. (3.15)

    By equating the terms in equation (3.13) with identical powers of λ, we obtain the following components:

    {u0(ˉx,t)=m1i=0tii!gi(ˉx),u1(ˉx,t)=f(q)t(ˉx,t)+L(q)tu0+E(q)0t,uk(ˉx,t)=L(q)tu(k1)+E(q)(k1)t,k=2,3, (3.16)

    Substituting the components from (3.16) into equation (3.10) gives the solution of equation (3.8). Now, from equations (3.7) and (3.10), we obtain

    u(ˉx,t)=limλ1uλ(ˉx,t)=u0(ˉx,t)+u1(ˉx,t)+k=2uk(ˉx,t). (3.17)

    By using the first equations of (3.17), we see that iu(ˉx,0)ti=limλ1iuλ(ˉx,0)ti, which implies that gi(ˉx)=fi(ˉx). Inserting (3.16) into (3.17) completes the proof.

    In the following results, we establish the convergence theorem and maximum absolute error theorem of the analytical solution given by (3.6) for the nonlinear partial differential equation given by (3.1).

    Theorem 3.4. Let B be a Banach space. Then the series solution of equation (3.16) converges to SB for , if there exists γ,0γ<1 such that, unγu(n1) for nN.

    Proof. Define the sequence Sn of partial sums of the series given by equation (3.16) as:

    {S0=u0(ˉx,t),S1=u0(ˉx,t)+u1(ˉx,t),S2=u0(ˉx,t)+u1(ˉx,t)+u2(ˉx,t),Sn=u0(ˉx,t)+u1(ˉx,t)+u2(ˉx,t)++un(ˉx,t), (3.18)

    and we need to show that {Sn} is a Cauchy sequence in Banach space B. For this aim, we introduce

    S(n+1)Sn=u(n+1)(ˉx,t)γun(ˉx,t)γ2u(n1)(ˉx,t)γn+1u0(ˉx,t). (3.19)

    For every n,nN,nn, by using equation (3.19) and triangle inequality successively, we have

    SnSn=S(n+1)Sn+S(n+2)S(n+1)++SnS(n1)S(n+1)Sin+S(n+2)S(n+1)++SnS(n1)γn+1u0(ˉx,t)+γn+2u0(ˉx,t)++γnu0(ˉx,t)=γn+1(1+γ++γnn1)u0(ˉx,t)γn+1(1γnn1γ)u0(ˉx,t). (3.20)

    Since 0<γ<1, so 1γnn1. Then

    SnSnγn+11γu0(ˉx,t). (3.21)

    Since u0(ˉx,t) is bounded, then

    limn,nSnSn=0. (3.22)

    Therefore, {Sn} is a Cauchy sequence in the Banach space B, so the series solution defined in equation (3.7) is convergent and the proof is complete.

    Theorem 3.5. The maximum absolute truncation error of the series solution (3.7) of the nonlinear fractional partial differential equation (3.1) is estimated to be

    sup(ˉx,t)Ω|u(ˉx,t)nk=0u(ˉx,t)|γn+11γsup(ˉx,t)Ω|u0(ˉx,t)|, (3.23)

    where the region ΩRn+1.

    Proof. From Theorem 3.4, we have

    SnSnγn+11γsup(ˉx,t)Ω|u0(ˉx,t)|. (3.24)

    But we assume that Sn=nk=0u(ˉx,t) and since n, we obtain Snu(ˉx,t), so equation (3.24) can be rewritten as:

    u(ˉx,t)Sn=u(ˉx,t)nk=0u(ˉx,t)γn+11γsup(ˉx,t)Ω|u0(ˉx,t)|. (3.25)

    Therefore, the maximum absolute error is

    sup(ˉx,t)Ω|u(ˉx,t)nk=0u(ˉx,t)|γn+11γsup(ˉx,t)Ω|u0(ˉx,t)|, (3.26)

    and this completes the proof.

    In this section, we present examples of some well-known nonlinear fractional partial differential equations, namely, fractional KdV equations. The following examples are chosen due to their closed form solutions which are available in the literature or they have been solved previously by some other well-known methods.

    Example 4.1. Consider the following Korteweg-de Vries (KdV) equation of fractional order with initial value:

    {Dqtu(x,t)+u(x,t)ux(x,t)+uxxx(x,t)=0,u(x,0)=13tanh2(x2+1). (4.1)

    For q=1, the exact solitary wave solution of the mKdV equation (4.1) is given by

    uEx(x,t)=13tanh2(t2+x2+1). (4.2)

    To solve equation (4.1) along with the initial condition, we compare (4.1) with equation (3.1), we obtain

    Dqtu(x,t)=N(u)uxxx(x,t), (4.3)

    where we assume N(u)=u(x,t)ux(x,t).

    Next, we assume equation (4.1) has a solution given by

    u(x,t)=k=0uk(x,t). (4.4)

    To obtain the approximate solution of equation (4.3), we consider the following equation:

    Dqtuλ=λ[N(uλ)uλxxx(x,t)], (4.5)

    subject to the initial conditions given by

    uλ(x,0)=g(x), (4.6)

    and we assume that equation (4.5) has a solution of the following form:

    uλ(x,t)=k=0λkuk(x,t). (4.7)

    By applying Riemann-Liouville fractional partial integral of order q with respect to t to both sides of equation (4.5) and by using Theorem 2.2 and equation (4.6), we obtain

    uλ=g(x)+λIqt[N(uλ)uλxxx(x,t)]. (4.8)

    We use Remark 3.1 and equation (4.7) into equation (4.8), we obtain

    k=0λkuk=g(x)+λIqt[αn=0λnEnk=0λkukxxx]. (4.9)

    We equate the terms in equation (4.9) with identical powers of λ, we obtain the following components:

    {u0=g(x),u1=Iqt[E0u0xxx],uk=Iqt[E(k1)u(k1)xxx], (4.10)

    for k=2,3, where E(k1) can be obtained by using Definition 3.

    By using equations (4.4) and (4.7), we can write

    u(x,t)=limλ1uλ(x,t)=k=0uk(x,t). (4.11)

    By using equation (4.11), we have u(x,0)=limλ1uλ(x,0), it gives that g(x)=u(x,0).

    Consequently, we use the components from (4.10) into (4.11), using Definition 3 and by Mathematica software, we obtain the first few components of the solution for the initial value problem (4.1) as follows:

    u0(x,t)=13tanh2(x2+1),u1(x,t)=24tqsinh4(x2+1)csch3(x+2)Γ(q+1),u2(x,t)=12ex+2(ex+2+1)4Γ(2q+1)(ex+2(ex+24)+1)t2q,u3(x,t)=12ex+2(ex+21)t3q(ex+2+1)7Γ(q+1)2Γ(3q+1)[12ex+2(ex+2×(ex+24)+1)Γ(2q+1)(ex+2(ex+2(ex+2×(ex+2+16)114)+16)+1)Γ(q+1)2]. (4.12)

    Therefore, the approximate solution of third-order term for equation (4.1) is given by

    u(x,t)=13tanh2(t2+x2+1)24tqsinh4(x2+1)csch3(x+2)Γ(q+1)+12ex+2(ex+2(ex+24)+1)t2q(ex+2+1)4Γ(2q+1)+12ex+2(ex+21)t3q(ex+2+1)7Γ(q+1)2Γ(3q+1)×[12ex+2(ex+2(ex+24)+1)Γ(2q+1)(ex+2(ex+2(ex+2(ex+2+16)114)+16)+1)Γ(q+1)2]. (4.13)

    Example 4.2. Consider the following modified Korteweg-de Vries (mKdV) equation of fractional order with initial value:

    {Dqtu(x,t)α(u(x,t))2ux(x,t)+uxxx(x,t)=0,u(x,0)=62αtanh(x2+1), (4.14)

    where α is a constant. When q=1, the exact solitary wave solution of the mKdV equation (4.14) is given by

    uEx(x,t)=62αtanh(t4+x2+1). (4.15)

    To obtain the solution for equation (4.1) along with the initial condition, we compare (4.14) with equation (3.1), we obtain

    Dqtu(x,t)=N(u)uxxx(x,t), (4.16)

    where we assume N(u)=α(u(x,t))2ux(x,t).

    Further, we suppose that equation (4.1) has a solution given by

    u(x,t)=k=0uk(x,t). (4.17)

    To get a solution of equation (4.16), we consider the following equation:

    Dqtuλ=λ[N(uλ)uλxxx(x,t)], (4.18)

    subject to the initial conditions given by

    uλ(x,0)=g(x), (4.19)

    and we assume that equation (4.18) has a solution of the following form:

    uλ(x,t)=k=0λkuk(x,t). (4.20)

    By applying the Riemann-Liouville fractional partial integral of order q with respect to t for both sides of equation (4.18) and by using Theorem 2.2 and equation (4.19), we obtain

    uλ=g(x)+λIqt[N(uλ)uλxxx(x,t)]. (4.21)

    Using Remark 3.1 and equation (4.7), in equation (4.8), we obtain

    k=0λkuk=g(x)+λIqt[n=0λnEnk=0λkukxxx]. (4.22)

    We equate the terms in equation (4.22) with identical powers of λ, we obtain the following components:

    {u0=g(x),u1=Iqt[E0u0xxx],uk=Iqt[E(k1)u(k1)xxx], (4.23)

    for k=2,3,, where the nonlinear polynomial E(k1) can be obtained by using Definition 3.

    From equations (4.17) and (4.20), we can write

    u(x,t)=limλ1uλ(x,t)=k=0uk(x,t). (4.24)

    We use equation (4.11) into (4.24), we have u(x,0)=limλ1uλ(x,0), it gives that g(x)=u(x,0).

    Consequently, we use the components from (4.23) into (4.24) and using Definition 3 with the help of Mathematica software, the first few components of the solution for the initial value problem (4.14) are obtained as follows:

    u0(x,t)=6tanh(x2+1)2α,u1(x,t)=32ex+2tqα(ex+2+1)2Γ(q+1),u2(x,t)=32ex+2(ex+21)t2q2α(ex+2+1)3Γ(2q+1),u3(x,t)=32ex+2t3q4α(ex+2+1)6Γ(q+1)2Γ(3q+1)×[(ex+2(ex+2(ex+2(ex+2+22)78)+22)+1)×Γ(q+1)212ex+2(ex+2(ex+23)+1)Γ(2q+1)],

    and so on. Therefore, the approximate solution of third-order term for equation (4.14) is given by

    u(x,t)=6tanh(x2+1)2α+32ex+2tqα(ex+2+1)2Γ(q+1)32ex+2(ex+21)t2q2α(ex+2+1)3Γ(2q+1)+32ex+2t3q4α(ex+2+1)6Γ(q+1)2Γ(3q+1)×[(ex+2(ex+2(ex+2(ex+2+22)78)+22)+1)Γ(q+1)212ex+2(ex+2(ex+23)+1)Γ(2q+1)].

    The numerical values in Table 1 are of approximate solution uAAM(x,t) obtained by the proposed AAM and the exact solution uEx(x,t) given by equation (4.2) through different values of x,t when q=0.5,0.75 and 1 for Example 4.1. The Figure 1 and Fig. 2 show the graphs of the approximate solution among different values of x,t when q=0.50 and q=0.75 respectively for Example 4.1. In Fig. 3, we plotted the graph of the approximate solution among different values of x,t when q=1 for Example 4.1. The Fig. 4 shows the graph of the exact solution among different values of x,t for Example 4.1. In Fig. 5 and Fig. 6, we plotted the graphs of the approximate and exact solutions respectively among different values of t,q when x is fixed for Example 4.1. From the Fig. 5, we observe that the approximate solution for Example 4.2 decreases when t and q decrease for fixed values of x. In Fig. 6, the exact solution for Example 4.1 decreases when t decreases for fixed values of x.

    Table 1.  Numerical values of the approximate and exact solutions among different values of x,t when q=0.5,0.75 and 1 for Example 4.1.
    x t q=0.5 q=0.75 q=1 Absolute Error
    uAAM uAAM uAAM uEx |uEXuAAM|
    -0.5 0.25 -0.73475 -0.61587 -0.48622 -0.48650 0.000233
    0.50 -0.82575 -0.82914 -0.736403 -0.74008 0.003674
    0.75 -0.81728 -0.94952 -0.94678 -0.96491 0.018129
    0 0.25 -1.17201 -1.06435 -0.96475 -0.96490 0.000159
    0.50 -1.29271 -1.23203 -1.15638 -1.15876 0.002377
    0.75 -1.36508 -1.34165 -1.31108 -1.32229 0.011207
    0.5 0.25 -1.47394 -1.39198 -1.32222 -1.32229 0.000062
    0.50 -1.58180 -1.51238 -1.45700 -1.45788 0.000882
    0.75 -1.66942 -1.60008 -1.56485 -1.56880 0.003946

     | Show Table
    DownLoad: CSV
    Figure 1.  The graph of the approximate solution among different values of x,t when q=0.50 for Example 4.1.
    Figure 2.  The graph of the approximate solution among different values of x,t when q=0.75 for Example 4.1.
    Figure 3.  The graph of the approximate solution among different values of x,t when q=1 for Example 4.1.
    Figure 4.  The graph for the exact solution among different values of x,t for Example 4.1.
    Figure 5.  The graph of the approximate solution among different values of t,q when x is fixed for Example 4.1.
    Figure 6.  The graph for the exact solution among different values of t when x is fixed for Example 4.1.

    In Table 2, we evaluated the numerical values of the approximate solution uAAM(x,t) obtained by the proposed AAM and the exact solution uEx(x,t) given by equation (4.15) through different values of x,t when q=0.5,0.75,1 to make a numerical comparison between uAAM and uEx. The Fig. 7 and Fig. 8 show the graphs of the approximate solution among different values of x,t for α=2 when q=0.50 and q=0.75 respectively for Example 4.2. The Fig. 9 shows the graph of the approximate solution among different values of x,t when q=1 and α=2 for Example 4.1. In Fig. 10, we plotted the graph of the exact solution among different values of x,t when α=2 for Example 4.2. In Fig. 11 and Fig. 12, we plotted the graphs of the approximate and exact solutions respectively among different values of t,q when x is fixed and α=2 for Example 4.2. From the Fig. 11, we observe that the approximate solution for Example 4.2 increases when t increases and q decreases for fixed values of x and α=2. In Fig. 12, the approximate solution for Example 4.2 decreases when t increases for fixed values of x and α=2. From the previous argument and from the numerical values in Table 1 and Table 2, we observe that absolute error is very small which means that the proposed AAM is very effective to obtain the analytical solutions for NFPDEs easily without any assumption.

    Table 2.  Numerical values of the approximate and exact solutions among different values of x,t when q=0.5,0.75,1 and α=2 for Example 4.2.
    x t q=0.5 q=0.75 q=1 Absolute Error
    uAAM uAAM uAAM uEx |uEXuAAM|
    -0.5 0.25 0.61400 0.59605 0.58107 0.58106 1.2653×106
    0.50 0.63635 0.62351 0.60958 0.60960 0.0000194034
    0.75 0.65239 0.64526 0.63563 0.63573 0.0000940212
    0 0.25 0.70393 0.69153 0.68123 0.68123 3.2894×107
    0.50 0.71966 0.71035 0.70087 0.70088 4.6986×106
    0.75 0.73139 0.72529 0.71861 0.71863 0.0000210896
    0.5 0.25 0.76384 0.75578 0.74903 0.74903 1.81425×107
    0.50 0.77414 0.76808 0.76196 0.76195 3.12052×106
    0.75 0.78189 0.77782 0.77355 0.77353 0.0000167879

     | Show Table
    DownLoad: CSV
    Figure 7.  The graph of the approximate solution among different values of x,t when q=0.50 and α=2 for Example 4.2.
    Figure 8.  The graph of the approximate solution among different values of x,t when q=0.75 and α=2 for Example 4.2.
    Figure 9.  The graph of the approximate solution among different values of x,t when q=1 and α=2 for Example 4.2.
    Figure 10.  The graph for the exact solution among different values of x,t when α=2 for Example 4.2.
    Figure 11.  The graph of the approximate solution among different values of t,q when x is fixed and α=2 for Example 4.2.
    Figure 12.  The graph for the exact solution among different values of t when x is fixed and α=2 for Example 4.2.

    An approximate-analytical method for solving full general NFPDEs was introduced in this paper and it was efficaciously applied to solve fractional KdV equations with the aid of the Caputo derivative. The obtained travelling wave solutions to the fractional KdV equations can be applied for further study in analytical fields such as stability analysis in the area of applied sciences for instance mathematics and physics. The real benefit of the AAM is that it provides extra general and a great number of original travelling wave solutions by way of some unrestricted parameters. It is valued to state that putting into practice of the proposed method is very dependable, well-organized and it applicable to solve further fractional nonlinear wave equations. We used Maple and Mathematica software to obtain the numerical values and plotting the graphs.

    Author contributions: Hayman Thabet and Subhash Kendre contributed substantially to this paper. Hayman Thabet wrote this paper; Subhash Kendre Supervised development of the paper and James Peters helped to evaluate and edit the paper.

    The research by James Peters has been supported by the Natural Sciences & Engineering Research Council of Canada (NSERC) discovery grant 185986, Instituto Nazionale di Alta Matematica (INdAM) Francesco Severi, Gruppo Nazionale per le Strutture Algebriche, Geometriche e Loro Applicazioni grant 9 920160 000362, n.prot U 2016/000036 and Scientific and Technological Research Council of Turkey (TÜBİTAK) Scientific Human Resources Development (BIDEB) under grant no: 2221-1059B211301223.

    The authors thank the referees for their valuable comments and suggestions.

    The authors declare no conflict of interest.



    [1] K. A. Touchent, Z. Hammouch, T. Mekkaoui, et al. Implementation and Convergence Analysis of Homotopy Perturbation Coupled With Sumudu Transform to Construct Solutions of Local-Fractional PDEs, Fractal and Fractional, 2 (2018), 22.
    [2] A. Akbulut and F. Taşcan, Lie symmetries, symmetry reductions and conservation laws of time fractional modified korteweg–de vries (mkdv) equation, Chaos Soliton. Fract., 100 (2017), 1-6. doi: 10.1016/j.chaos.2017.04.020
    [3] N. A. Asif, Z. Hammouch, M. B. Riaz, et al. Analytical solution of a Maxwell fluid with slip effects in view of the Caputo-Fabrizio derivative, Eur. Phys. J. Plus, 133 (2018), 272.
    [4] A. M. A. El-Sayed, A. Elsaid, I. L. El-Kalla, et al. A homotopy perturbation technique for solving partial differential equations of fractional order in finite domains, Appl. Math. Comput., 218 (2012), 8329-8340.
    [5] B. Guo, X. Pu and F. Huang, Fractional partial differential equations and their numerical solutions, World Scientific, 2015.
    [6] P. K. Gupta and M. Singh, Homotopy perturbation method for fractional fornberg–whitham equation, Comput. Math. Appl., 61 (2011), 250-254. doi: 10.1016/j.camwa.2010.10.045
    [7] Q. Huang and R. Zhdanov, Symmetries and exact solutions of the time fractional Harry-Dym equation with Riemann–Liouville derivative, Physica A, 409 (2014), 110-118. doi: 10.1016/j.physa.2014.04.043
    [8] A. Keten, M. Yavuz and D. Baleanu, Nonlocal Cauchy Problem via a Fractional Operator Involving Power Kernel in Banach Spaces, Fractal and Fractional, 3 (2019), 27.
    [9] A. A. Kilbas, H. M. Srivastava and J. J. Trujillo, Theory and applications of fractional differential equations, Elsevier Science Limited, 2006.
    [10] W. X. Ma and Z. Zhu, Solving the (3 + 1)-dimensional generalized KP and BKP equations by the multiple exp-function algorithm, Appl. Math. Comput., 218 (2012), 11871-11879.
    [11] K. S. Miller and B. Ross, An introduction to the fractional calculus and fractional differential equations, Wiley-Interscience, 1993.
    [12] S. Momani, Z. Odibat and V. S. Erturk, Generalized differential transform method for solving a space- and time-fractional diffusion-wave equation, Phys. Lett. A, 370 (2007), 379-387. doi: 10.1016/j.physleta.2007.05.083
    [13] G. Paquin and P. Winternitz, Group theoretical analysis of dispersive long wave equations in two space dimensions, Physica D: Nonlinear Phenomena, 46 (1990), 122-138. doi: 10.1016/0167-2789(90)90115-6
    [14] I. Podlubny, Fractional differential equations: an introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications, Academic press, 1998.
    [15] M. B. Riaz, M. A. Imran and K. Shabbir, Analytic solutions of Oldroyd-B fluid with fractional derivatives in a circular duct that applies a constant couple, Alex. Eng. J., 55 (2016), 3267-3275. doi: 10.1016/j.aej.2016.07.032
    [16] M. B. Riaz and A. A. Zafar, Exact solutions for the blood flow through a circular tube under the influence of a magnetic field using fractional Caputo-Fabrizio derivatives, Math. Model. Nat. Phenom., 13 (2018), 8.
    [17] R. Sahadevan and T. Bakkyaraj, Invariant analysis of time fractional generalized Burgers and Korteweg-de Vries equations, J. Math. Anal. Appl., 393 (2012), 341-347. doi: 10.1016/j.jmaa.2012.04.006
    [18] H. Thabet and S. Kendre, Analytical solutions for conformable space-time fractional partial differential equations via fractional differential transform, Chaos Soliton. Fract., 109 (2018), 238-245. doi: 10.1016/j.chaos.2018.03.001
    [19] H. Thabet and S. Kendre, Modified least squares homotopy perturbation method for solving fractional partial differential equations, Malaya Journal of Matematik, 6 (2018), 420-427. doi: 10.26637/MJM0602/0020
    [20] H. Thabet and S. Kendre, New modification of Adomian decomposition method for solving a system of nonlinear fractional partial differential equations, Int. J. Adv. Appl. Math. and Mech, 6 (2019), 1-13.
    [21] H. Thabet, S. Kendre and D. Chalishajar, New Analytical Technique for Solving a System of Nonlinear Fractional Partial Differential Equations, Mathematics, 5 (2017), 47.
    [22] H. Thabet, S. Kendre and J. Peters, Analytical Solutions for Nonlinear Systems of Conformable Space-Time Fractional Partial Differential Equations via Generalized Fractional Differential Transform, Vietnam Journal of Mathematics, (2019), 1-21.
    [23] M. Yavuz and E. Bonyah, New approaches to the fractional dynamics of schistosomiasis disease model, Physica A, 525 (2019), 373-393. doi: 10.1016/j.physa.2019.03.069
    [24] M. Yavuz and N. Özdemir, Comparing the new fractional derivative operators involving exponential and mittag-leffler kernel, Discrete & Continuous Dynamical Systems-S, (2019), 1098-1107.
    [25] E. M. E. Zayed, Y. A. Amer and R. M. A. Shohib, The fractional complex transformation for nonlinear fractional partial differential equations in the mathematical physics, Journal of the Association of Arab Universities for Basic and Applied Sciences, 19 (2016), 59-69.
  • This article has been cited by:

    1. Umar Farooq, Hassan Khan, Fairouz Tchier, Evren Hincal, Dumitru Baleanu, Haifa Bin Jebreen, New approximate analytical technique for the solution of time fractional fluid flow models, 2021, 2021, 1687-1847, 10.1186/s13662-021-03240-z
    2. Hayman Thabet, Subhash Kendre, James Peters, Melike Kaplan, Solitary wave solutions and traveling wave solutions for systems of time-fractional nonlinear wave equations via an analytical approach, 2020, 39, 2238-3603, 10.1007/s40314-020-01163-1
    3. Melike Kaplan, Asma Rashid Butt, Hayman Thabet, Arzu Akbulut, Nauman Raza, Dipankar Kumar, An effective computational approach and sensitivity analysis to pseudo-parabolic-type equations, 2021, 1745-5030, 1, 10.1080/17455030.2021.1989081
    4. Aisha Abdullah Alderremy, Rasool Shah, Nehad Ali Shah, Shaban Aly, Kamsing Nonlaopon, Comparison of two modified analytical approaches for the systems of time fractional partial differential equations, 2023, 8, 2473-6988, 7142, 10.3934/math.2023360
    5. 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, 2022, 7, 2473-6988, 12483, 10.3934/math.2022693
    6. Hayman Thabet, Subhash Kendre, Subhash Unhale, Mawardi Bahri, Numerical Analysis of Iterative Fractional Partial Integro-Differential Equations, 2022, 2022, 2314-4785, 1, 10.1155/2022/8781186
    7. Azam Zahrani, Mashaallah Matinfar, Mostafa Eslami, Arzu Akbulut, Bivariate Chebyshev Polynomials to Solve Time-Fractional Linear and Nonlinear KdV Equations, 2022, 2022, 2314-4785, 1, 10.1155/2022/6554221
    8. Hayman Thabet, Subhash Kendre, James Peters, Advances in solving conformable nonlinear partial differential equations and new exact wave solutions for Oskolkov‐type equations, 2022, 45, 0170-4214, 2658, 10.1002/mma.7945
    9. Jian-Gen Liu, Jian Zhang, A new approximate method to the time fractional damped Burger equation, 2023, 8, 2473-6988, 13317, 10.3934/math.2023674
    10. Hayman Thabet, Subhash Kendre, Conformable mathematical modeling of the COVID‐19 transmission dynamics: A more general study, 2023, 46, 0170-4214, 18126, 10.1002/mma.9549
    11. M. Nagaraja, H. B. Chethan, T. R. Shivamurthy, Mohd Asif Shah, D. G. Prakasha, Semi-analytical approach for the approximate solution of Harry Dym and Rosenau–Hyman equations of fractional order, 2024, 11, 2768-4830, 10.1080/27684830.2024.2401662
    12. Dilek Varol, Solitary and Periodic Wave Solutions of the Space-Time Fractional Extended Kawahara Equation, 2023, 7, 2504-3110, 539, 10.3390/fractalfract7070539
    13. Hegagi Mohamed Ali, Analytical investigation of the fractional nonlinear shallow-water model, 2024, 1598-5865, 10.1007/s12190-024-02172-7
    14. Hayman Thabet, Subhash Kendre, A fractional mathematical model with nonlinear partial differential equations for transmission dynamics of severe acute respiratory syndrome coronavirus 2 infection, 2023, 4, 27724425, 100209, 10.1016/j.health.2023.100209
    15. H. B. Chethan, Rania Saadeh, D. G. Prakasha, Ahmad Qazza, Naveen S. Malagi, M. Nagaraja, Deepak Umrao Sarwe, An efficient approximate analytical technique for the fractional model describing the solid tumor invasion, 2024, 12, 2296-424X, 10.3389/fphy.2024.1294506
    16. Huda Alsaud, Hassan Eltayeb, The Four-Dimensional Natural Transform Adomian Decomposition Method and (3+1)-Dimensional Fractional Coupled Burgers’ Equation, 2024, 8, 2504-3110, 227, 10.3390/fractalfract8040227
  • Reader Comments
  • © 2019 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(5506) PDF downloads(643) Cited by(16)

Figures and Tables

Figures(12)  /  Tables(2)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog