In this article, time fractional Fornberg-Whitham equation of He’s fractional derivative is studied. To transform the fractional model into its equivalent differential equation, the fractional complex transform is used and He’s homotopy perturbation method is implemented to get the approximate analytical solutions of the fractional-order problems. The graphs are plotted to analysis the fractional-order mathematical modeling.
Time fractional Fornberg-Whitham equationfractional complex transformHe’s homotopy perturbation methodIntroduction
Over the last few years, the study of the fractional calculus and applications in the area of life science, physics and the engineering has been paid a great attention. The fractional calculus are also used in many other fields, such as optics, solitary waves, control theory of dynamical systems, and so on, which can be derived by linear or nonlinear fractional order differential equations. Recently, the studies of nonlinear problems and their effects are of widely significance. Many analytical and approximation methods have been presented to solve nonlinear fractional differential equations such as [1–11]. The homotopy perturbation method (HPM) [12–14] is widely applied to various science and engineering problems. This method was first proposed by He [12]. Fractional complex transform was suggested also by He et al. [15–19], which converts the fractional differential equation into its equivalent differential equation, so that the HPM can be effectively used. Now it is considered as a powerful method to find the approximation solutions of nonlinear fractional order differential equations.
In this paper, we study the time fractional Fornberg-Whitham (FW) equation [20] as follows:
here β is the fractal dimensions of the fractal medium and ∂β∂tβ is He’s fractional derivative defined [21–23]:
∂βu∂tβ=1Γ(n−β)dndtn∫t0t(s−t)n−β−1[u0(s)−u(s)]ds
where u0(x,t) is the solution of its starting point of the nonlinear fractional order model. When β=1, Eq. (1) becomes to be the original Fornberg-Whitham (FW) model which is a significant mathematical equation in mathematical physics. This equation was presented by Parkes et al. [24–26] in 1978 which describes nonlinear water waves with peakon solutions. The peakon solution is a special solitary wave solution which is peaked in the limiting case. The FW equation has been found to require peakon results as a simulation for limiting wave heights as well as the frequency of wave breaks. Now many scholars have researched the FW model for the fractional-order derivative because of its fractional calculus applications.
HPM Procedure
To illustrate the basic ideas of this method [27–32], we consider the following nonlinear functional equation:
A(u)−f(r)=0,r∈Ω
with the following boundary conditions:
B(u−∂u∂n)=0,r∈Λ
where A is a general functional operator, B is a boundary operator, f(r) is a known analytical function, and Λ is the boundary of the domain Ω. The operator A can be decomposed into two operators L and N, where L is linear, and N is nonlinear operator. Eq. (4) can be written as follows:
L(u)+N(u)−f(r)=0
Using the homotopy technique, we construct a following homotopy:
where the homotopy parameter p is considered as a small parameter.
Obviously, we can get:
H(ξ,0)=L(ξ)−L(ξ0)
H(ξ,1)=A(ξ)−f(r)
The HPM uses embedding parameter p as an expanding parameter, and basic assumption is that the solution of Eq. (6) can be specified as a power series in p:
ξ=ξ0+pξ1+p2ξ2+p3ξ3+…
When p→1, the approximate solution of Eq. (1):
u=limp→1ξ=ξ0+ξ1+ξ2+ξ3+…
HPM Implementation
According to the HPM, the first step is to transform the fractional model equation into its equivalent differential equation:
y=tβΓ(1+β)
Eq. (1) can be written into its equivalent differential equation form:
∂u∂y−∂3u∂x2∂y+∂u∂x−u∂3u∂x3+u∂u∂x−3∂u∂x∂2u∂x2=0
Applying the HPM process, substituting Eq. (11) into Eq. (7), we have a set of equations that is to be simultaneously solved:
In Figs. 1 and 2, the approximate solution obtained by fractional complex transform and HPM and actual solution of Eq. (1) are plotted. It is observed that HPM solutions are in closed contact with the exact solution. In Figs. 3 and 4, the solutions of Eq. (1) at various fractional-order of the derivatives are plotted which showns that the balancing phenomena between dispersion and nonlinearity is valid.
Approximate solution of Eq. (1) for β=1
Exact solution of Eq. (1) for β=1
Approximate solution of Eq. (1) for β=0.8
Approximate solution of Eq. (1) for β=0.6
In Figs. 5 and 6, the the approximate solutions of Eq. (1) at different fractional order of the derivatives are plotted. This behavior illustrated that there is directly relationship between both the width and the height of the solitary wave and the value of β. We can change the shape of solitary wave with the help of adjusting the value of fractal order β are plotted. In Figs. 7 and 8, the numerical results illustrated which confirms the high accuracy of the method.
The relationship between β and amplitude of Eq. (1) for t=0.3
The relationship between β and amplitude of Eq. (1) for t=0.5
The relationship between exact solution (blue) and approximate solution (red) of Eq. (1) for β=1, t=0.005
The relationship between exact solution and approximate solution (red) of Eq. (1) for β=1, t=0.015
Conclusion
In this manuscript, a mathematical technology is used to find the solution of time fractional Fornberg-Whitham equation. The fractional-derivatives are discussed within He’s fractional derivative. The solutions are determined for fractional-order problems which shows the high accuracy and efficiency. This is easy and can be extended to other nonlinear differential equations with fractal derivatives in science and engineering.
Funding Statement: The work was supported by the National Natural Science Foundation of China under Grant No. 11561051.
Conflicts of Interest: The authors declare that we have no conflicts of interest to report regarding the present study.
ReferencesHe, J. H., Ji, F. Y. (2019). Taylor series solution for Lane-Emden equation. ,57(8),1932–1934. DOI 10.1007/s10910-019-01048-7.He, J. H. (2019). The simplest approach to nonlinear oscillators. ,15(1),102546. DOI 10.1016/j.rinp.2019.102546.Wang, Y., An, Y., Wang, X. Q. (2019). A variational formulation for anisotropic wave traveling in a porous medium. ,27(4),1950047. DOI 10.1142/S0218348X19500476.Wang, K. L., He, C. H. (2019). A remark on Wang’s fractal variational principle. ,27(8),1950134. DOI 10.1142/S0218348X19501342.He, C. H., Shen, Y., Ji, F. Y., He, J. H. (2019). Taylor series solution for fractal bratu-type equation arising in electrospinning process. ,28(1),2050011.He, J. H. (2020). Taylor series solution for a third order boundary value problem arising in architectural engineering. ,11(4),1411–1414. DOI 10.1016/j.asej.2020.01.016.Ahmad, H. (2019). Variational iteration algorithmi with an auxiliary parameter for wave-like vibration equations. ,38(3–4),1113–1124. DOI 10.1177/1461348418823126.Ahmad, H., Seadawy, A. R., Khan, T. A., Thounthong, P. (2020). Analytic approximate solutions for some nonlinear parabolic dynamical wave equations. ,14(1),346–358. DOI 10.1080/16583655.2020.1741943.He, J. H. (2021). A modified Li-He’s variational principle for plasma. (ahead-of-print).He, J. H. (2019). Lagrange crisis and generalized variational principle for 3-D unsteady flow. ,30(3),1189–1196. DOI 10.1108/HFF-07-2019-0577.He, J. H., Sun, C. (2019). A variational principle for a thin film equation. ,57(5),2075–2081. DOI 10.1007/s10910-019-01063-8.He, J. H. (1999). Homotopy perturbation technique. ,178(3–4),257–262. DOI 10.1016/S0045-7825(99)00018-3.He, J. H. (2000). A coupling method of a homotopy technique and a perturbation technique for non-linear problems. ,35(1),37–43. DOI 10.1016/S0020-7462(98)00085-7.He, J. H. (2005). Application of homotopy perturbation method to non-linear wave equation. ,26(3),695–700. DOI 10.1016/j.chaos.2005.03.006.He, J. H., Li, Z. B. (2012). Converting fractional differential equations into partial differential equations. ,1(2),331–334. DOI 10.2298/TSCI110503068H.Li, Z. B., He, J. H. (2010). Fractional complex transform for fractional differential equations. ,15(5),970–973. DOI 10.3390/mca15050970.Li, Z. B., Zhu, W. H., He, J. H. (2012). Exact solutions of time fractional heat conduction equation by the fractional complex transform. ,16(2),335–338. DOI 10.2298/TSCI110503069L.Ain, Q. T., He, J. H. (2019). On two-scale dimension and its applications. ,23(3B),1707–1712.He, J. H., Ji, F. Y. (2019). Two-scale mathematics and fractional calculus for thermodynamics. ,23(4),2131–2133. DOI 10.2298/TSCI1904131H.Singh, J., Kumar, D., Kumar, S. (2013). New treatment of fractional Fornberg-Whitham equation via laplace transform. ,4(3),557–562. DOI 10.1016/j.asej.2012.11.009.Wang, Q. L., Shi, X., He, J. H., Li, Z. B. (2018). Fractal calculus and its application explanation of biomechanism of polar bear hairs. ,26(6),1850086. DOI 10.1142/S0218348X1850086X.He, J. H. (2014). A tutorial review on fractal space time and fractional calculus. ,3(11),3698–3371.He, J. H. (2011). A new fractal derivation. ,15(1),S145–S147. DOI 10.2298/TSCI11S1145H.Parkes, E. J., Vakhnenko, V. O. (2005). Explicit solutions of the Camassa-Holm equation. ,26(5),1309–1316. DOI 10.1016/j.chaos.2005.03.011.Whitham, G. B. (1967). Variational methods and applications to water waves. ,299(1456),6–25.Fornberg, B., Whitham, G. B. (1978). A numerical and theoretical study of certain nonlinear wave phenomena. ,289(1361),373–404.He, J. H. (2005). Periodic solutions and bifurcations of delay dierential equations. ,347(4–6),228–230. DOI 10.1016/j.physleta.2005.08.014.He, J. H. (2005). Application of homotopy perturbation method to nonlinear wave equations. ,26(3),695–700. DOI 10.1016/j.chaos.2005.03.006.He, J. H. (2005). Limit cycle and bifurcation of nonlinear problems. ,26(3),827–833. DOI 10.1016/j.chaos.2005.03.007.Das, S., Gupta, P. K. (2010). An approximate analytical solution of the fractional diusion equation with absorbent term and external force by Homotopy Pertur bation Method. ,65(3),182–190. DOI 10.1515/zna-2010-0305.Yildirim, A., Kocak, H. (2009). Homotopy perturbation method for solving the space-time fractional advection dispersion equation. ,32(12),1711–1716. DOI 10.1016/j.advwatres.2009.09.003.He, J. H. (2006). Homotopy perturbation method for solving boundary value problems. ,350(1–2),87–88. DOI 10.1016/j.physleta.2005.10.005.