MPEMathematical Problems in Engineering1563-51471024-123XHindawi10.1155/2021/66980286698028Research ArticleNew Explicit Solutions to the Fractional-Order Burgers’ Equationhttps://orcid.org/0000-0003-3725-5472UddinM. Hafiz1https://orcid.org/0000-0002-2892-1683ArefinMohammad Asif1https://orcid.org/0000-0001-5688-6259AkbarM. Ali2https://orcid.org/0000-0003-4996-8373IncMustafa345PeraMaria Patrizia1Department of MathematicsJashore University of Science and TechnologyJashore 7408Bangladesh2Department of Applied MathematicsUniversity of RajshahiRajshahi 6205Bangladeshru.ac.bd3Department of Computer EngineeringBiruni UniversityIstanbulTurkeybiruni.edu.tr4Department of MathematicsFirat UniversityElazigTurkeyfirat.edu.tr5Department of Medical ResearchChina Medical University HospitalChina Medical UniversityTaichungTaiwancmu.edu.cn202112620212021711202095202146202112620212021Copyright © 2021 M. Hafiz Uddin et al.This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

The closed-form wave solutions to the time-fractional Burgers’ equation have been investigated by the use of the two variables G/G,1/G-expansion, the extended tanh function, and the exp-function methods translating the nonlinear fractional differential equations (NLFDEs) into ordinary differential equations. In this article, we ascertain the solutions in terms of tanh, sech, sinh, rational function, hyperbolic rational function, exponential function, and their integration with parameters. Advanced and standard solutions can be found by setting definite values of the parameters in the general solutions. Mathematical analysis of the solutions confirms the existence of different soliton forms, namely, kink, single soliton, periodic soliton, singular kink soliton, and some other types of solitons which are shown in three-dimensional plots. The attained solutions may be functional to examine unidirectional propagation of weakly nonlinear acoustic waves, the memory effect of the wall friction through the boundary layer, bubbly liquids, etc. The methods suggested are direct, compatible, and speedy to simulate using algebraic computation schemes, such as Maple, and can be used to verify the accuracy of results.

1. Introduction

The nonlinear fractional evolution equations (NLFEEs) emerge frequently in diverse research field of science and applications of engineering. The fractional derivative has been happening in numerous physical problems, for example, recurrence subordinate damping conduct of materials, motion of an enormous meager plate in a Newtonian fluid, creep and relaxation functions for viscoelastic materials, and PIλDμ controller for the control of the dynamical system. Fractional-order differential equations describe the phenomena. The fractional-order differential equations are broadly used as generalizations of conventional differential equations with the integral order to explain different intricate phenomena in numerous fields including the diffusion of biological populations, electric circuit, fluid flow, chemical kinematics, control theory, signal processing, optical fiber, plasma physics, solid-state physics, and other areas . The concepts of dissipation, dispersion, diffusion, convection, and reaction are closely related to the abovestated phenomena, and nonlinear fractional partial differential equations (NLFPDEs) can be used to evaluate them exactly. Wave shape has an effect on sediment transport and beach morpho dynamics, while wave skewness has an impact on radar altimetry signals, and asymmetry has an impact on ship responses to wave impacts. Traveling wave solutions are a special class of analytical solutions for NLFEEs. Solitary waves are transmitted traveling waves with constant speeds and shapes that achieve asymptotically zero at distant locations. The appearance of solitary waves in nature is rather frequent in plasmas, fluids dynamics, solid-state physics, condensed matter physics, chemical kinematics, optical fibers, electrical circuits, bio-genetics, elastic media, etc. Consequently, it is important to search for the exact traveling wave solutions of NLFPDEs to understand the facts. Therefore, many researchers have been motivated on finding the exact solutions to nonlinear fractional-order differential equations, and significant progress has been made in analyzing the exact solutions of these types of equations. The major challenges, however, are that there is no unified numerical or analytical approach that can investigate all sorts of nonlinear fractional-order differential equations. Thus, several numerical and theoretical methods for finding solutions for NLFDEs have been established, for example, the differential transformation method [6, 7], the variational iteration method , the fractional subequation method , the Kudryashov  method, the homotopy perturbation method [13, 14], the homotopy analysis method , the exp-function method [16, 17], the G/G-expansion method and its various modification , the Chelyshkov polynomial method [23, 24], the multiple exp-function method , the finite difference method , the finite element method , the first integral method [28, 29], the modified simple equation method , the reproducing kernel method , the two variables G/G,1/G-expansion method [32, 33], and the Picard technique .

The time-fractional Burgers’ equation is crucial for modeling shallow water waves, weakly nonlinear acoustic waves propagating unidirectionally in gas-filled tubes, and bubbly liquids. Inc  studied the approximate and exact solutions to the time-fractional Burgers’ equation by the variational iteration method. Bekir and Guner  established the exact solution to the mentioned equation by using the G/G-expansion method. Bulut et al.  examined the analytical approximate solution to the suggested equation through the modified trial equation method. Recently, Saad and Al-Sharif  studied the exact and analytical solutions to this equation. As far as is known, the stated equation has not been investigated through the two variables G/G,1/G-expansion technique, exp-function strategy, and expanded tanh function method. Therefore, the aim of this study is to establish further general and some fresh solutions of the abovementioned equation using the suggested methods.

The residual segments of the article is schematized as follows: in Section 2, definition and preliminaries have been introduced; in Section 3, the two variables G/G,1/G-expansion method, the exp-function method, and the extended tanh function method have been described. In Section 4, the exact solutions to the suggested equation have established. In Section 5, physical interpretation and explanation of the extracted solutions are provided. In the lattermost part, the conclusions are given.

2. Definition and Preliminaries

Supposef:0, be a function. The α-order conformable derivative of f is interpreted as (1)Tαft=limε0ft+εt1αftε,for every t>0 and α0,1. If f is α-differentiable in some 0,a, a>0, and limt0+fαt exists; then, fα0=limt0+fαt. The following theorems point out few axioms that are satisfied conformable derivatives.

Theorem 1.

Consider α0,1 and let us suppose f and g be α-differentiable at a point t>0. Therefore,

Tαcf+dg=cTαf+dTαg, for all c and d

Tαtp=ptpα, for all p

Tαc=0, for all constant function ft=c

Tαfg=fTαg+gTαf

Tαf/g=gTαffTαg/g2

In addition, if f is differentiable, then Tαft=t1αdf/dt

Some more properties including the chain rule, Gronwall’s inequality, some integration techniques, Laplace transform, Tailor series expansion, and exponential function with respect to the conformable fractional derivative are explained in .

Theorem 2.

Let f be an α-differentiable function in conformable differentiable, and suppose that g is also differentiable and defined in the range of f. Then,(2)Tαfgt=t1αgtfgt.

The Caputo derivative is another important fractional derivative concept developed by Michele Caputo . This definition is particularly useful for finding numerical solutions. The definition of Riesz [40, 41] in relation to the fractional derivative, on the contrary, is also important for extracting numerical solutions. The two concepts are not discussed in depth here since the aim of this article is to establish exact solutions.

3. Outline of the Methods

In this part, we summarize the principal parts of the suggested methods to analyze exact traveling wave solutions to the NLFEEs. Assume the general NLFEE is of the form(3)Pu,Dtαu,Dxβu,DtαDtαu,DtαDxβu,DxβDxβ,=0,0<α1 and 0<β1,where u represents an unknown function, consisting the spatial derivative x and temporal derivative t, and P represents a polynomial of ux,t and its derivatives where the highest order of derivatives and nonlinear terms of the highest order are associated. Take into account the wave transformation(4)ξ=kxββ+ctαα,ux,t=uξ,where c and k are nonzero arbitrary constants.

By means of wave transformation (4), equation (3) can be rewritten as(5)Ru,u,u,u,=0,where the superscripts specify the ordinary derivative of u relating to ξ.

3.1. The Two Variables <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M54"><mml:mfenced open="(" close=")" separators="|"><mml:mrow><mml:mfenced open="(" close=")" separators="|"><mml:mrow><mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mi>G</mml:mi></mml:mrow><mml:mrow><mml:mo>′</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo>/</mml:mo><mml:mi>G</mml:mi></mml:mrow></mml:mrow></mml:mfenced><mml:mo>,</mml:mo><mml:mfenced open="(" close=")" separators="|"><mml:mrow><mml:mrow><mml:mn>1</mml:mn><mml:mo>/</mml:mo><mml:mi>G</mml:mi></mml:mrow></mml:mrow></mml:mfenced></mml:mrow></mml:mfenced></mml:math></inline-formula>-Expansion Method

Step 1: In this subsection, we apply the two variables G/G,1/G-expansion method to acquire the wave solutions of the NLFEEs. Take into account the second order ODEs(6)Gξ+λGξ=μ,

along with the following relations(7)ϕ=GG,ψ=1G.

In this manner, it gives(8)ϕ=ϕ2+μψλ,ψ=ϕψ.

The solutions to equation (6) depend on λ as λ<0,λ>0, and λ=0.

Case 1: when λ<0, the general solution to equation (6) is(9)Gξ=A1sinhλξ+A2coshλξ+μλ.

In view of that, we obtain(10)ψ2=λλ2σ+μ2ϕ22μψ+λ,

where σ=A12A22.

Case 2: ifλ>0, the solution to (6) is given as follows:(11)Gξ=A1sinλξ+A2cosλξ+μλ.

Therefore, we obtain(12)ψ2=λλ2σμ2ϕ22μψ+λ,

where σ=A12+A22.

Case 3: when λ=0, the solution of equation (6) is(13)Gξ=μ2ξ2+A1ξ+A2.

Therefore, we find(14)ψ2=1A122μA2ϕ22μψ.

where A1 and A2 are arbitrary constants.

Step 2: in agreement with two variables G/G,1/G-expansion scheme, the solution of (5) is presented as a polynomial of ϕ and ψ of the form(15)uξ=i=0Naiϕi+i=1Nbiϕi1ψ,

where ai and bi are arbitrary constants to be determined later.

Step 3: after balancing the maximum order of derivatives and nonlinear terms, which appear in equation (5), it can be fixed the positive integer N.

Step 4: setting (15) into (5) along with (8) and (10), this modifies to a polynomial in ϕ and ψ having the degree of ψ as one or less than one. If we compare the polynomial of similar terms to zero, then it will give a set of mathematical equations which can be unraveled by computational software and finally yield the values of ai, bi, μ, A1, A2, and λ, where λ<0; this condition provides solutions of the hyperbolic function.

Step 5: in a similar manner, we can examine the values of ai, bi, μ, A1, A2, and λ, and trigonometric and rational solutions can be established separately for the case of λ>0 andλ=0.

3.2. The Exp-Function Method

Within this section, the key components of the exp-function method are described for searching the traveling wave solution to the NLFDEs.

Step 1: the arrangement is to be communicated in the shape as indicated by the exp-function method:(16)uξ=n=cdpnexpnξm=pqqmexpmξ,

where c, d,p, and q are unknown positive integers, which can be evaluated later, and pn and qm are unidentified constants.

Step 2: the balancing principle between the highest-order linear and nonlinear terms presented in (5) and substituting (16) into (5) yieldc and p, and the balance of lowest-order linear and nonlinear terms yields the values of d and q.

Step 3: introducing (16) into (5) and setting the coefficient of expnξ to zero provides an arrangement of set of mathematical equations for pn, qm, c, and k. Then, unraveling the set with the aid of computer software, such as Maple, we attain the constants.

Step 4: substituting the values that showed up in step 3 into (16), we ascertain exact solutions to the NLFEEs in (3).

3.3. The Extended Tanh Function Method

In this section, the suggested extended tanh function method has been interpreted to obtain ample exact solutions to NLFEEs which was summarized by Wazwaz . The basic concept of this method is to present the solution as a polynomial of hyperbolic functions, and then, solving the coefficient of tanhμξ implies solving a system of algebraic equations. The core steps of the extended tanh function method for finding exact analytic solutions of nonlinear PDEs of the fractional order are as follows:

Step 1: we consider the wave solution as follows:(17)uξ=i=0NaiYi+i=1nbiYi,

wherein(18)Y=tanhμξ,

where μ is any arbitrary constant.

Step 2: taking uniform balance between the maximum order nonlinear term and the derivative of the maximum order appearing in equation (5) to determine the positive constant N.

Step 3: substitute solution (17) together with (18) into equation (5) with the value of n acquired in step 2, which yields the polynomials in Y. A set of algebraic equations for ai’s and bi’s are found by setting each the coefficient of the resulted polynomials to zero. With the help of symbolic computational software, namely, Maple, this set of equations for ai and bi can be solved.

Step 4: inserting the values that appeared in step 3 into equation (17) along with equation (18), we construct closed-form traveling wave solutions of nonlinear evolution equation (3).

4. Analysis of the Solutions

Here, we search further comprehensive exact analytic wave solutions for the stated time-fractional Burgers’ equation by means of the suggested methods. Let us consider the time-fractional Burgers’ equation as follows:(19)Dtα+puuxvuxx=0;t>0 and 0<α1,where p and v are arbitrary constants. The physical processes of unidirectional propagation of weakly nonlinear acoustic waves through a gas-filled pipe are described by the time-fractional Burgers’ equation. The fractional derivative results are obtained from the memory effect of the wall friction through the boundary layer. The similar formation can be found in several systems, namely, waves in bubbly liquids and shallow water waves. For equation (19), we recommend the subsequent wave transformation:(20)ξ=kxααctαα,ux,t=uξ,where c be the velocity of the traveling wave. For wave transformation (20), time-fractional Burgers’ equation (19) reduced to the ensuing integral order differential equation:(21)cu+kpuuk2vu=0.

Integrating equation (21) with zero constant, we obtain(22)cu+kpu22k2vu=0.

4.1. Solutions through Two Variables <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M132"><mml:mfenced open="(" close=")" separators="|"><mml:mrow><mml:mfenced open="(" close=")" separators="|"><mml:mrow><mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mi>G</mml:mi></mml:mrow><mml:mrow><mml:mo>′</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo>/</mml:mo><mml:mi>G</mml:mi></mml:mrow></mml:mrow></mml:mfenced><mml:mo>,</mml:mo><mml:mfenced open="(" close=")" separators="|"><mml:mrow><mml:mrow><mml:mn>1</mml:mn><mml:mo>/</mml:mo><mml:mi>G</mml:mi></mml:mrow></mml:mrow></mml:mfenced></mml:mrow></mml:mfenced></mml:math></inline-formula>-Expansion Method

Considering the homogeneous balance of the highest-order nonlinear term and highest-order derivative showing up in equation (22), the arrangements of equation (15) accept the shape(23)uξ=a0+a1ϕξ+b1ψξ,where a0,a1, and b1 are constants to be determined.

Case 1: for λ<0, embedding solution (23) into (22) along with equations (8) and (10) yields a set of algebraic equations, and by explaining these equations by computer algebra such as Maple, we achieve the following results:(24)a0=±λb11λ2σ+μ2,a1=±b1λλ2σμ2,c=±p2λb12λvλ2σ+μ2,k=±b1pλvλ2σ+μ2.

Inserting the top values into solution (23), we find the solution to equation (19) in the form(25)u11x,t=±λb11λ2σ+μ2±b1λλ2σμ2×λA1coshλξ+A2sinhλξA1sinhλξ+A2coshλξ+μ/λ±b1×1A1sinhλξ+A2coshλξ+μ/λ,

where σ=A12A22 and ξ=±b1pλ/vλ2σ+μ2xα/αp2λb12λ/vλ2σ+μ2tα/α.

Since A1 and A2 are basic constants, one might have picked self-assertively their values. If we take μ=0 and A1=0,A20 or A10,A2=0 in (25), we have(26)u12x,t=±b1σ±b1λσtanhλξ+b1sechλξ,(27)u13x,t=±b1σ±b1λσcothλξ+b1cosechλξ,

where ξ=±b1pλ/vλ2σ+μ2xα/αp2λb12λ/vλ2σ+μ2tα/α.

Case 2: in a comparative way, when λ>0, substituting (23) into (22) together with (8) and (12) yields an arrangement of algebraic equations for a0,a1, b1, and ω, and we acquire the following results by working out these equations:(28)a0=±λb11λ2σ+μ2,a1=±b1λλ2σμ2,c=±p2λb12λvλ2σ+μ2,k=±b1pλvλ2σ+μ2.

The substitution of these results into solution (23) possesses the following expression for the general solution of equation (19):(29)u14x,t=±λb11λ2σ+μ2±b1λλ2σμ2×λA1cosλξ+A2sinλξA1sinλξ+A2cosλξ+μ/λ±b1×1A1sinλξ+A2cosλξ+μ/λ,

where σ=A12+A22 and ξ=±b1pλ/vλ2σ+μ2xα/αp2λb12λ/vλ2σ+μ2tα/α.

If the unknown parameters are assigned as μ=0 and A1=0 and A20 or A10 and A2=0 in solution (29), it provides the next solitary wave solution:(30)u15x,t=±b1σ±b1λσ×tanλξ+b1×secλξ,(31)u16x,t=±b1σ±b1λσ×cotλξ+b1×cosecλξ,

where ξ=±b1pλ/vλ2σ+μ2xα/αp2λb12λ/vλ2σ+μ2tα/α.

Case 3: in the parallel algorithm when λ=0, using equations (22) and (21) along with (8) and (14), we achieve a set of mathematical equations whose solutions are

(32)a0=0,a1=±b1A122μA2,b1=b1,c=c,k=±pb1v1A122μA2.

Making use of these values into solution (23) produces the solution to equation (19) as(33)u17x,t=±b1A122μA2×μξ+A1μ/2ξ2+A1ξ+A2+b1μ/2ξ2+A1ξ+A2,where ξ=±pb1/v1/A122μA2xα/αctα/α.

It is substantial to observe that the traveling wave solutions u11u17 of the studied equation are inclusive and standard. The attained solutions have not been noted in the earlier study. These solutions are convenient to designate the physical processes of unidirectional propagation of weakly nonlinear acoustic waves via a gas-filled tube, shallow-water waves, and waves in bubbly liquids.

4.2. Solution by the Exp-Function Method

Considering the homogeneous balance, the solution of equation (16) takes the form(34)uξ=p1expξ+p0+p1expξq1expξ+q0+q1expξ.

Substituting equation (34) into (22) leads a equation in expnξ; here, n represents any whole number. Inserting each coefficient of this equation to zero yields a cluster of mathematical equations (for straightforwardness, here, we have discarded) for pis, qis, and ω. These mathematical equations are solved by computer algebra, namely, Maple, which gives the following outcomes:(35)Set 1: c=14p2p12q12v,k=12pp1q1v,p1=0,p0=0,p1=p1,q1=0,q0=q0,and q1=q1,Set 2: c=14p2p02q02v,k=12pp0q0v,p1=0,p0=p0,p1=0,q1=0,q0=q0,and q1=q1,Set 3: c=c,k=2cq1pp1,p1=0,p0=p0,p1=p1,q1=0,q0=p0q1p1,and q1=q1,Set 4: c=14p2p12q12v,k=12pp1q1v,p1=0,p0=p0,p1=p1,q1=q1,q0=p12q1+q1p02p0p1,and q1=q1,Set 5: c=18p2p12q12v,k=14pp1q1v,p1=p1,p0=0,p1=0,q1=q1,q0=0,and q1=q1,Set 6: c=c,k=2cq1pp1,p1=p1,p0=p0,p1=0,q1=q1,q0=p0q1p1,and q1=0,Set 7: c=14p2p12q12v,k=12pp1q1v,p1=p1,p0=p0,p1=0,q1=q1,q0=p12q1+q1p02p0p1,and q1=q1,Set 8: c=c,k=2cq0pp0,p1=p1,p0=p0,p1=p1,q1=p1q0p0,q0=q0,and q1=p1q0p0,Set 9: c=c,k=2cq1pp1,p1=p1,p0=0,p1=p1,q1=q1,q0=0,and q1=p1q1p1.

From the point of view of the above results, we achieve the following generalized solitary wave solutions:(36)u21x,t=p1exp1/2pp1/q1vxα/α+1/4p2p12/q12vtα/αq0+q1exp1/2pp1/q1vxα/α+1/4p2p12/q12vtα/α,(37)u22x.t=p0q0+q1exp1/2pp0/q0vxα/α+1/4p2p02/q02vtα/α,(38)u23x,t=p0+p1exp2cq1/pp1xα/αctα/αp0q1/p1+q1exp2cq1/pp1xα/αctα/α,(39)u24x,t=p0+p1exp1/2pp1/q1vxα/α+1/4p2p12/q12vtα/αq1exp1/2pp1/q1vxα/α+1/4p2p12/q12vtα/α+p12q1+q1p02/p0p1+q1exp1/2pp1/q1vxα/α+1/4p2p12/q12vtα/α,(40)u25x,t=p1exp1/4pp1/q1vxα/α1/8p2p12/q12vtα/αq1exp1/4pp1/q1vxα/α1/8p2p12/q12vtα/α+q1exp1/4pp1/q1vxα/α1/8p2p12/q12vtα/α,(41)u26x,t=p0+p1exp2cq1/pp1xα/αctα/αp0q1/p1+q1exp2cq1/pp1xα/αctα/α,(42)u27x,t=p0+p1exp1/2pp1/q1vxα/α+1/4p2p12/q12vtα/αq1exp1/2pp1/q1vxα/α+1/4p2p12/q12vtα/α+p12q1+q1p02/p0p1+q1exp1/2pp1/q1vxα/α+1/4p2p12/q12vtα/α,(43)u28x,t=p1exp2cq0/pp0xα/αctα/α+p0+p1exp2cq0/pp0xα/αctα/αp0q0/p0exp2cq0/pp0xα/αctα/α+q0+q1exp2cq0/pp0xα/αctα/α,(44)u29x,t=p1exp2cq1/pp1xα/αctα/α+p1exp2cq1/pp1xα/αctα/αp1q1/p1exp2cq1/pp1xα/αctα/α+p1q0/p0exp2cq1/pp1xα/αctα/α.

In particular, if p1=q1=p1 and p0=q0=1, solution (44) is simplified and offers the kink type solution of the form(45)u210x,t=tanh2cq1pp1xααctαα.

The choice of p1=q1=p1 and p0=q0=1 in (44) gives the singular kink solution:(46)u211x,t=coth2cq1pp1xααctαα.

It is significant to refer that the traveling wave solutions u21u211 of the considered Burgers’ equation are fresh and standard and were not established in the earlier investigations. It is deduced that physical systems should be assigned of unidirectional propagation of weakly nonlinear acoustic waves through a gas-filled tunnel and waves in bubbly fluids.

4.3. Solution Using the Extended Tanh Function Method

The homogeneous symmetry allows solution equation (17) as(47)uε=a0+a1Y+b1Y1.

Substituting (47) into (22) along with (18) makes the left hand side as a polynomial in Y. Setting each coefficient of this polynomial to zero, resulting a set of algebraic equations (for simplicity, we have omitted them to exhibition) for a0, a1, b1, k, and c. Computing the determined set of equations with the assistance of computer algebra, such as Maple, yields the succeeding results:(48)Set 1: c=2k2vμ,k=k,a0=2kvμp,a1=0,and b1=2kvμp,Set 2: c=2k2vμ,k=k,a0=2kvμp,a1=0,and b1=2kvμp,Set 3: c=2k2vμ,k=k,a0=2kvμp,a1=2kvμp,and b1=0,Set 4: c=4k2vμ,k=k,a0=4kvμp,a1=2kvμp,and b1=2kvμp,Set 5: c=2k2vμ,k=k,a0=2kvμp,a1=2kvμp,and b1=0,Set 6: c=4k2vμ,k=k,a0=4kvμp,a1=2kvμp,and b1=2kvμp.

Using the values of the parameters assembled above into solution (47) together with (18), we achieve the following solitary wave solutions:(49)u31x,t=2kvμp2kvμpcothμkxαα+2k2vμtαα,(50)u32x,t=2kvμp2kvμpcothμkxαα2k2vμtαα,(51)u33x,t=2kvμp2kvμptanhμkxαα+2k2vμtαα,(52)u34x,t=4kvμp2kvμptanhμkxαα+4k2vμtαα2kvμpcothμkxαα+4k2vμtαα,(53)u35x,t=2kvμp2kvμptanhμkxαα2k2vμtαα,(54)u36x,t=4kvμp2kvμptanhμkxαα4k2vμtαα2kvμpcothμkxαα4k2vμtαα.

The solutions established above by the extended tanh approach are advanced and progressive. These might be convenient to describe the relativistic electron and the physical processes of unidirectional propagation of weakly nonlinear acoustic waves via a gas-filled tube.

5. Physical Interpretation and Explanation

In this section, we mainly discuss about the physical interpretation of the determined solitary wave solutions, including kink, singular solitons, singular kink, and periodic wave of the NLFEEs. A graph is an effective approach for explaining mathematical concepts. It is capable of describing any circumstances in a straightforward and understandable manner. This segment explains the incidents by portraying 3D plots of some of the solutions that are found. The portraits are precedents of the solutions shown in Figures 16 using the computational software, namely, Mathematica.

3D plot of the kink type soliton of (26) sketched within the interval 0x50 and 0t50.

3D plot of the single soliton solution of (51) sketched within the interval 10x10 and 0t700.

3D plot of the periodic wave solution of (31) sketched within the interval 0x100 and 0t100.

3D plot of the singular kink type soliton of (46) sketched within the interval 0x and t10.

3D plot of the kink type soliton of (53) sketched within the interval 0x10 and 0t10.

3D plot of the singular kink type soliton of 0t10. (54) sketched within the interval 0x1000 and 0t1000.

The results of the time-fractional Burgers’ equation include the kink soliton, singular soliton, periodic soliton, and some general solitons which are displayed in Figures 16. Figure 1 is the kink shape soliton of solution (26) with the values of the parameters λ=1, v=1, p=1, b1=1, σ=2, and α=1/2 within the interval 0x50 and 0t50. The kink soliton is a soliton which rises or descends from one asymptotic state to another as ξ. Solution (51) represents the shape of the plane soliton characterized in Figure 2 for the values of parameters k=1, p=1, v=1, μ=1, and α=1/2 within the interval 10x10 and 0t700. Solution (31) represents the periodic wave solutions, plotted for λ=1, v=1, p=1, b1=1, σ=2, and α=1/2 within the interval 0x and t100 and labeled in Figure 3. When c=1, p=1, p1=2, q1=1, and α=1/2, solution (46) represents the singular kink type soliton characterized in Figure 4 within 0x and t10. On the contrary, for the values of k=1, p=1, v=1, μ=1, and α=1/2, solution (53) also represents the kink soliton illustrated in Figure 5 within the interval 0x and  t10. Finally, outcome (54) also represents the singular kink soliton for the values of parameters k=1, p=1, v=1, μ=1, and α=1/2 within the range 0x and t1000, which is labeled as Figure 6. The other figure of the solutions is analogous to the displayed figure; thus, for convenience, these are omitted here.

6. Conclusion

In this article, using three reliable approaches referring conformable the fractional derivative, we have established scores of advanced, further general, and wide-ranging solitary wave solutions to the time-fractional Burgers’ equation. The ascertained closed-form solutions of the considered equation include kink, single solitons, periodic solitons, singular kink, and some other kinds of solutions, including some free parameters. The obtained solutions are capable to analyze the phenomena of weakly nonlinear acoustic waves propagating unidirectionally in gas-filled tubes, shallow water waves, and bubbly liquids. The dynamics of solitary waves have been graphically depicted in terms of space and time coordinates which reveal the consistency of the techniques used. The accuracy of the results obtained in this study has been verified using the computational software Maple by placing them back into NLFPDEs and found correct. This study shows that all the methods implemented are reliable, effective, functional, and capable of uncovering nonlinear fractional differential equations arising in the field of nonlinear science and engineering. Therefore, we can firmly claim that the implemented methods can be used to compute exact wave solutions of other nonlinear fractional equations associated with real-world problems, and this is our next contrivance.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

MillerK. S.RossB.An Introduction to the Fractional Calculus and Fractional Differential Equations1993New York, NY, USAWileyOldhamK. B.SpanierJ.The Fractional Calculus1974New York, NY, USAAcademic PressPodlubnyI.Fractional Differential Equations1999San Diego, CA, USAAcademicHilferR.Applications of Fractional Calculus in Physics2000River Edge, NJ, USAWorld Scientific PublishingDiethelmK.The Analysis of Fractional Differential Equations2010Berlin, GermanySpringer-VerlagErturkV. S.MomaniS.OdibatZ.Application of generalized differential transform method to multi-order fractional differential equationsCommunications in Nonlinear Science and Numerical Simulation20081381642165410.1016/j.cnsns.2007.02.0062-s2.0-39149140685WangK.-L.WangK.-J.A modification of the reduced differential transform method for fractional calculusThermal Science20182241871187510.2298/tsci1804871w2-s2.0-85052437215JiJ.ZhangJ.DongY.The fractional variational iteration method improved with the Adomian seriesApplied Mathematics Letters201225122223222610.1016/j.aml.2012.06.0072-s2.0-84865455237IncM.The approximate and exact solutions of the space- and time-fractional Burgers equations with initial conditions by variational iteration methodJournal of Mathematical Analysis and Applications2008345147648410.1016/j.jmaa.2008.04.0072-s2.0-43949121726GuoS.MeiL.The fractional variational iteration method using He's polynomialsPhysics Letters A2011375330931310.1016/j.physleta.2010.11.0472-s2.0-78650881841GuoS.MeiL.LiY.SunY.The improved fractional sub-equation method and its applications to the space-time fractional differential equations in fluid mechanicsPhysics Letters A2012376440741110.1016/j.physleta.2011.10.0562-s2.0-84855193408AkbarM. A.AkinyemiL.YaoS. W.Soliton solutions to the Boussinesq equation through sine-Gordon method and Kudryashov methodResults Physics20212511010.1016/j.rinp.2021.104228GepreelK. A.The homotopy perturbation method applied to the nonlinear fractional Kolmogorov-Petrovskii-Piskunov equationsApplied Mathematics Letters20112481428143410.1016/j.aml.2011.03.0252-s2.0-79955470464GuptaP. K.SinghM.Homotopy perturbation method for fractional Fornberg-Whitham equationComputers & Mathematics with Applications201161225025410.1016/j.camwa.2010.10.0452-s2.0-78651257332ArafaA. A. M.RidaS. Z.MohamedH.Homotopy analysis method for solving biological population modelCommunications in Theoretical Physics201156579780010.1088/0253-6102/56/5/012-s2.0-81555196325BekirA.GünerÖ.CevikelA. C.Fractional complex transform and exp-function methods for fractional differential equationsAbstract and Applied Analysis201320131810.1155/2013/4264622-s2.0-84880174334AkbarM. A.AliN. H. M.New solitary and periodic solutions of nonlinear evolution equation by exp-function methodWorld Applied Science Journal2012171216031610WangM.LiX.ZhangJ.The ()-expansion method and travelling wave solutions of nonlinear evolution equations in mathematical physicsPhysics Letters A2008372441742310.1016/j.physleta.2007.07.0512-s2.0-37549033511BatoolF.AkramG.New solitary wave solutions of the time-fractional Cahn-Allen equation via the improved G/G-expansion methodEuropean Physical Journal - Plus2018133511110.1140/epjp/i2018-12025-y2-s2.0-85046475622AkbarM. A.AliN. H. M.ZayedE. M. E.A generalized and improved (G′/G)-Expansion method for nonlinear evolution equationsMathematical Problems in Engineering201220122245987910.1155/2012/4598792-s2.0-84859768511AkbarM. A.AliN. H. M.ZayedE. M. E.Abundant exact traveling wave solutions of generalized bretherton equation via improved (G′/G)-Expansion methodCommunications in Theoretical Physics201257217317810.1088/0253-6102/57/2/012-s2.0-84856834811IslamM. N.AkbarM. A.SpagnoloB.New exact wave solutions to the space-time fractional-coupled Burgers equations and the space-time fractional foam drainage equationCogent Physics201851142295710.1080/23311940.2017.1422957HamidM.UsmanM.HaqR. U.TianZ.A spectral approach to analyze the nonlinear oscillatory fractional order differential equationsChaos, Solitons and Fractals2021146110921HamidM.UsmanM.WangW.TianZ.A stable computational approach to analyze semi-relativistic behavior of fractional evolutionary problemsNumerical Methods for Partial Differential Equations2020115HamidM.UsmanM.ZubairT.HaqR. U.ShafeeA.An efficient analysis for N-soliton, Lump and lump-kink solutions of time-fractional (2+1)-Kadomtsev-Petviashvili equationPhysica A: Statistical Mechanics and Its Applications201952812132010.1016/j.physa.2019.1213202-s2.0-85065716616UsmanM.HamidM.LiuM.Novel operational matrices-based finite difference/spectral algorithm for a class of time-fractional Burger equation in multidimensionsChaos, Solitons & Fractals202114411070110.1016/j.chaos.2021.110701SeadawyA. R.The generalized nonlinear higher order of KdV equations from the higher order nonlinear Schrödinger equation and its solutionsOptik2017139314310.1016/j.ijleo.2017.03.0862-s2.0-85016519768LuB.The first integral method for some time fractional differential equationsJournal of Mathematical Analysis and Applications2012395268469310.1016/j.jmaa.2012.05.0662-s2.0-84864025292BekirA.GunerO.UnsalO.The first integral method for exact solutions of nonlinear fractional differential equationsJournal of Computational and Nonlinear Dynamics20151021610.1115/1.40280652-s2.0-84940489366AkbarM. A.HjN.AliM.WazwazA.Closed form traveling wave solutions of nonlinear fractional evolution equations through the modified simple equation methodThermal Science2018221341352AkgülA.BaleanuD.IncM.TchierF.On the solutions of electrohydrodynamic flow with fractional differential equations by reproducing kernel methodOpen Physics201614168568910.1515/phys-2016-00772-s2.0-85014642459Hafiz UddinM.AkbarM. A.KhanM. A.HaqueM. A.Close form solutions of the fractional generalized reaction Duffing model and the density dependent fractional diffusion reaction equationApplied and Computational Mathematics20176417718410.11648/j.acm.20170604.13LiL.-x.LiE.-q.WangM.-l.The (G′/G, 1/G)-expansion method and its application to travelling wave solutions of the Zakharov equationsApplied Mathematics-A Journal of Chinese Universities201025445446210.1007/s11766-010-2128-x2-s2.0-78650192468HamidM.UsmanM.WangW.TianZ.Hybrid fully spectral linearized scheme for time‐fractional evolutionary equationsMathematical Methods in the Applied Sciences20214453890391210.1002/mma.6996BekirA.GünerÖ.Exact solutions of nonlinear fractional differential equations by ( G′/G )-expansion methodChinese Physics B2013221111020210.1088/1674-1056/22/11/1102022-s2.0-84888632803BulutH.BaskonusH. M.PandirY.The modified trial equation method for fractional wave equation and time fractional generalized burgers equationAbstract and Applied Analysis201320131810.1155/2013/6368022-s2.0-84884864472SaadK.Al-SharifE. H.Analytical study for time and time-space fractional Burgers' equationAdvances in Difference Equations20172017110.1186/s13662-017-1358-02-s2.0-85029915984KhalilR.Al HoraniM.YousefA.SababhehM.A new definition of fractional derivativeJournal of Computational and Applied Mathematics2014264657010.1016/j.cam.2014.01.0022-s2.0-84893186929CaputoM.Linear models of dissipation whose Q is almost frequency independent--IIGeophysical Journal International196713552953910.1111/j.1365-246x.1967.tb02303.x2-s2.0-84977255207HamidM.UsmanM.ZubairT.HaqR. U.WangW.Innovative operational matrices based computational scheme for fractional diffusion problems with the Riesz derivativeThe European Physical Journal Plus20191341048410.1140/epjp/i2019-12871-y2-s2.0-85073186036UsmanM.HamidM.Ul HaqR.LiuM.Linearized novel operational matrices-based scheme for classes of nonlinear time-space fractional unsteady problems in 2DApplied Numerical Mathematics202116235137310.1016/j.apnum.2020.12.021WazwazA.-M.The extended tanh method for abundant solitary wave solutions of nonlinear wave equationsApplied Mathematics and Computation200718721131114210.1016/j.amc.2006.09.0132-s2.0-34247279667