Scipediacontent: Created page with "==Abstract== A relatively new method called q-Homotopy Analysis Method (q-HAM) is adopted in this paper to obtain an analytical solution of the time fractional Rosenau–Hyma..."

2017-04-12T08:09:32Z

Created page with "==Abstract== A relatively new method called q-Homotopy Analysis Method (q-HAM) is adopted in this paper to obtain an analytical solution of the time fractional Rosenau–Hyma..."

New page

==Abstract==

A relatively new method called q-Homotopy Analysis Method (q-HAM) is adopted in this paper to obtain an analytical solution of the time fractional Rosenau–Hyman equation in series form. Our analysis shows the simplicity nature of the application of q-HAM to nonlinear fractional differential equations. The convergence rate of the method used is faster in the sense that just very few terms of the series solution are needed for a good approximation due to the presence of the auxiliary parameter ''h'' comparable to exact solutions. Numerical solution obtained by this method is compared with the exact solution and solutions obtained by other analytical methods of the equation under various conditions. The numerical results are obtained using Mathematica 9 and MATLAB R2012b.

==Keywords==

Rosenau–Hyman equation; Fractional derivative; q-Homotopy Analysis Method

==1. Introduction==

The frequent use of analytical methods in solving nonlinear differential equations has the restriction that the variables involved in such methods are discretized and therefore the numerical techniques are always rounded off, see [[#b0005|[1]]]. For some other analytical methods for solving nonlinear differential equations see [[#b0010|[2]]] and [[#b0015|[3]]].

The Rosenau–Hyman equation was found as a simplified model for the study of the role which nonlinear dispersion plays in pattern formation in liquid drops and it has found diverse applications in modelling of various problems in physics and engineering.

Modern advances in fractional differential equations are motivated by new examples of its applications in visco-elasticity, fluid mechanics, electro-chemistry, mathematical biology, and physics. Examples of such applications include the use of fractional derivatives in the model of nonlinear oscillation of earthquake [[#b0020|[4]]] and in fluid-dynamic traffic model to eliminate the deficiency arising from the assumption of continuum traffic flow [[#b0025|[5]]]. Experimental data obtained for seepage flow in porous media suggest that the differential equations associated with them are fractional and in recent times fractional differential equations have proven to be valuable tools for the modelling of many physical phenomena [[#b0030|[6]]].

Many analytical methods have been successfully put to use to obtain solutions of the Rosenau–Hyman equation, such as the semi-analytical methods, Dissipative perturbation methods, Variational Iteration Method and Homotopy Perturbation Method, see [[#b0035|[7]]], [[#b0040|[8]]] and [[#b0045|[9]]]. In recent, a modified HAM called q-Homotopy Analysis Method was introduced in [[#b0050|[10]]], see also [[#b0055|[11]]], [[#b0060|[12]]], [[#b0065|[13]]], [[#b0070|[14]]], [[#b0075|[15]]], [[#b0080|[16]]] and [[#b0085|[17]]]. The reliability of this method in terms of the rate of convergence was proven from the fact that the algorithm contains fraction factor which aids the convergence. The conclusion was reached that the method is better than the usual HAM in terms of fast convergence.

No attempt has been made regarding analytical solutions of the time-fractional Rosenau–Hyman equation using q-Homotopy Analysis Method to the best of our knowledge. In this paper, we consider this equation subject to some appropriate initial condition. Comparative analysis of our results is carried out with exact solution and with other analytical results for this problem. The effects of the auxiliary parameter ''h '' and the fractional order <math display="inline">\alpha </math> on the solution are also shown. The numerical results of the problems are presented graphically.

==2. Preliminaries==

We give some definitions and some known results in this section. Caputo’s fractional derivative is adopted in this work.

==Definition 2.1. ==

The Riemann–Liouville’s (RL) fractional integral operator of order <math display="inline">\alpha \geqslant 0</math>, of a function <math display="inline">f\in L^1(a\mbox{,}b)</math> is given as

{| class="formulaSCP" style="width: 100%; text-align: center;"
|-
|
{| style="text-align: center; margin:auto;"
|-
| <math>I^{\alpha }f(t)=\frac{1}{\Gamma (\alpha )}{\int }_0^t{\left(t-\tau \right)}^{\alpha -1}f(\tau )d\tau \mbox{,}\quad t>0\mbox{,}\quad \alpha >0\mbox{,}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (1)
|}

where <math display="inline">\Gamma </math> is the Gamma function and <math display="inline">I^0f(t)=f(t)</math>.

==Definition 2.2. ==

The fractional derivative in the Caputo’s sense is defined as [[#b0030|[6]]],

{| class="formulaSCP" style="width: 100%; text-align: center;"
|-
|
{| style="text-align: center; margin:auto;"
|-
| <math>D^{\alpha }f(t)=I^{n-\alpha }D^nf(t)=\frac{1}{\Gamma (n-\alpha )}{\int }_0^t{\left(t-\tau \right)}^{n-\alpha -1}f^{\left(n\right)}(\tau )d\tau \mbox{,}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (2)
|}

where <math display="inline">n-1<\alpha \leqslant n\mbox{,}\quad n\in N\mbox{,}\quad t>0</math>.

==Lemma 2.1. ==

''Let ''<math display="inline">t\in (a\mbox{,}b]</math>''. Then''

{| class="formulaSCP" style="width: 100%; text-align: center;"
|-
|
{| style="text-align: center; margin:auto;"
|-
| <math>\left[I^{\alpha }{\left(t-a\right)}^{\beta }\right](t)=</math><math>\frac{\Gamma (\beta +1)}{\Gamma (\beta +\alpha +1)}{\left(t-a\right)}^{\beta +\alpha }\mbox{,}\quad \alpha \geqslant 0\mbox{,}\quad \beta >0\mbox{.}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3)
|}

==3. q-Homotopy Analysis Method (q-HAM)==

Differential equation of the form

{| class="formulaSCP" style="width: 100%; text-align: center;"
|-
|
{| style="text-align: center; margin:auto;"
|-
| <math>N\left[D_t^{\alpha }u(x\mbox{,}t)\right]-f(x\mbox{,}t)=</math><math>0</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (4)
|}

is considered, where ''N '' is a nonlinear operator, <math display="inline">D_t^{\alpha }</math> denote the Caputo fractional derivative, <math display="inline">\left(x\mbox{,}t\right)</math> are independent variables, ''f'' is a known function and ''u'' is an unknown function. To generalize the original homotopy method, the zeroth-order deformation equation is constructed as


{| class="formulaSCP" style="width: 100%; text-align: center;"
|-
|
{| style="text-align: center; margin:auto;"
|-
| <math>(1-nq)L\left(\phi (x\mbox{,}t\mbox{;}q)-u_0(x\mbox{,}t)\right)=</math><math>qhH(x\mbox{,}t)\left(N\left[D_t^{\alpha }\phi (x\mbox{,}t\mbox{;}q)\right]-\right. </math><math>\left. f(x\mbox{,}t)\right)\mbox{,}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (5)
|}

where <math display="inline">n\geqslant 1\mbox{,}q\in \left[0\mbox{,}\frac{1}{n}\right]</math> denotes the so-called embedded parameter, ''L '' is an auxiliary linear operator, <math display="inline">h\quad \not =\quad 0</math> is an auxiliary parameter, and <math display="inline">H(x\mbox{,}t)</math> is a non-zero auxiliary function.

It is clearly seen that when <math display="inline">q=0</math> and <math display="inline">q=\frac{1}{n}</math>, Eq. [[#e0025|(5)]] becomes

{| class="formulaSCP" style="width: 100%; text-align: center;"
|-
|
{| style="text-align: center; margin:auto;"
|-
| <math>\phi (x\mbox{,}t\mbox{;}0)=u_0(x\mbox{,}t)\quad \mbox{and}\quad \phi \left(x\mbox{,}t\mbox{;}\frac{1}{n}\right)=</math><math>u(x\mbox{,}t)</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (6)
|}

respectively. So, as ''q '' increases from 0 to <math display="inline">\frac{1}{n}</math>, the solution <math display="inline">\phi (x\mbox{,}t\mbox{;}q)</math> varies from the initial guess <math display="inline">u_0(x\mbox{,}t)</math> to the solution <math display="inline">u(x\mbox{,}t)</math>.

If <math display="inline">u_0(x\mbox{,}t)\mbox{,}L\mbox{,}h\mbox{,}H(x\mbox{,}t)</math> are chosen appropriately, solution <math display="inline">\phi (x\mbox{,}t\mbox{;}q)</math> of Eq. [[#e0025|(5)]] exists for <math display="inline">q\in \left[0\mbox{,}\frac{1}{n}\right]</math>.

Expansion of <math display="inline">\phi (x\mbox{,}t\mbox{;}q)</math> in Taylor series gives


{| class="formulaSCP" style="width: 100%; text-align: center;"
|-
|
{| style="text-align: center; margin:auto;"
|-
| <math>\phi (x\mbox{,}t\mbox{;}r)=u_0(x\mbox{,}t)+\sum_{m=1}^{\infty }u_m(x\mbox{,}t)q^m\mbox{.}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (7)
|}

where

{| class="formulaSCP" style="width: 100%; text-align: center;"
|-
|
{| style="text-align: center; margin:auto;"
|-
| <math>u_m(x\mbox{,}t)={\frac{1}{m!}\frac{{\partial }^m\phi (x\mbox{,}t\mbox{;}q)}{\partial q^m}}_{q=0}\mbox{.}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (8)
|}

Assume that the auxiliary linear operator ''L '', the initial guess <math display="inline">u_0</math>, the auxiliary parameter ''h '' and <math display="inline">H(x\mbox{,}t)</math> are properly chosen such that the series [[#e0035|(7)]] converges at <math display="inline">q=\frac{1}{n}</math>, then we have

{| class="formulaSCP" style="width: 100%; text-align: center;"
|-
|
{| style="text-align: center; margin:auto;"
|-
| <math>u(x\mbox{,}t)=u_0(x\mbox{,}t)+\sum_{m=1}^{\infty }u_m(x\mbox{,}t){\left(\frac{1}{n}\right)}^m\mbox{.}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (9)
|}

Let the vector <math display="inline">u_n</math> be defined as follows:

{| class="formulaSCP" style="width: 100%; text-align: center;"
|-
|
{| style="text-align: center; margin:auto;"
|-
| <math>{\overset{\rightarrow}{u}}_n=\left\{u_0(x\mbox{,}t)\mbox{,}u_1(x\mbox{,}t)\mbox{,}\ldots \mbox{,}u_n(x\mbox{,}t)\right\}\mbox{.}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (10)
|}

Differentiating Eq. [[#e0025|(5)]]''m''-times with respect to the (embedding) parameter ''q '', then evaluating at <math display="inline">q=0</math> and finally dividing them by <math display="inline">m!</math>, we have what is known as the <math display="inline">m</math>th-order deformation equation [[#b0090|[18]]] as

{| class="formulaSCP" style="width: 100%; text-align: center;"
|-
|
{| style="text-align: center; margin:auto;"
|-
| <math>L\left[u_m(x\mbox{,}t)-{\chi }_m^{{_\ast}}u_{m-1}(x\mbox{,}t)\right]=</math><math>hH(x\mbox{,}t)R_m\left({\overset{\rightarrow}{u}}_{m-1}\right)\mbox{.}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (11)
|}

with initial conditions

{| class="formulaSCP" style="width: 100%; text-align: center;"
|-
|
{| style="text-align: center; margin:auto;"
|-
| <math>u_m^{\left(k\right)}(x\mbox{,}0)=0\mbox{,}\quad k=</math><math>0\mbox{,}1\mbox{,}2\mbox{,}\ldots \mbox{,}m-1\mbox{.}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (12)
|}

where

{| class="formulaSCP" style="width: 100%; text-align: center;"
|-
|
{| style="text-align: center; margin:auto;"
|-
| <math>R_m\left({\overset{\rightarrow}{u}}_{m-1}\right)={\frac{1}{(m-1)!}\frac{{\partial }^{m-1}\left(N[D_t^{\alpha }\phi (x\mbox{,}t\mbox{;}q)]-f(x\mbox{,}t)\right)}{\partial q^{m-1}}}_{q=0}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (13)
|}

and


{| class="formulaSCP" style="width: 100%; text-align: center;"
|-
|
{| style="text-align: center; margin:auto;"
|-
| <math>{\chi }_m^{{_\ast}}=\begin{array}{ll}
0 & m\leqslant 1\\
n & \mbox{otherwise,}
\end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (14)
|}

==4. Fractional Rosenau–Hyman equation==

We consider the time-fractional Rosenau–Hyman equation.

Let


{| class="formulaSCP" style="width: 100%; text-align: center;"
|-
|
{| style="text-align: center; margin:auto;"
|-
| <math>\frac{{\partial }^{\alpha }u}{\partial t^{\alpha }}=</math><math>u\frac{{\partial }^3u}{\partial x^3}+u\frac{\partial u}{\partial x}+</math><math>3\frac{\partial u}{\partial x}\frac{{\partial }^2u}{\partial x^2}\mbox{,}\quad t>0\mbox{,}\quad 0<\alpha \leqslant 1</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (15)
|}

subject to the initial condition

{| class="formulaSCP" style="width: 100%; text-align: center;"
|-
|
{| style="text-align: center; margin:auto;"
|-
| <math>u(x\mbox{,}0)=-\frac{8c}{3}{cos}^2\left(\frac{x}{4}\right)\mbox{.}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (16)
|}

The exact solution to this problem for <math display="inline">\alpha =1</math> and ''c'' arbitrary constant, is given [[#b0095|[19] as follows:]]

{| class="formulaSCP" style="width: 100%; text-align: center;"
|-
|
{| style="text-align: center; margin:auto;"
|-
| <math>u(x\mbox{,}t)=-\frac{8c}{3}{cos}^2\left(\frac{x-ct}{4}\right)\mbox{.}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (17)
|}

Analytical solution of this problem is considered in [[#b0035|[7]]] using VIM and HPM.

===4.1. Application of q-HAM===

To solve the problem given in [[#e0075|(15)]] by q-HAM, we choose the linear operator

{| class="formulaSCP" style="width: 100%; text-align: center;"
|-
|
{| style="text-align: center; margin:auto;"
|-
| <math>L[\phi (x\mbox{,}t\mbox{;}q)]=D_t^{\alpha }\phi (x\mbox{,}t\mbox{;}q)</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (18)
|}

with property that <math display="inline">L[d_1]=0\mbox{,}d_1</math> is constant.

Using <math display="inline">u_0(x\mbox{,}t)=-\frac{8c}{3}{cos}^2\left(\frac{x}{4}\right)</math> as initial approximation, we define the non-linear operator as

{| class="formulaSCP" style="width: 100%; text-align: center;"
|-
|
{| style="text-align: center; margin:auto;"
|-
| <math>N[\phi (x\mbox{,}t\mbox{;}q)]=D_t^{\alpha }\phi (x\mbox{,}t\mbox{;}q)-</math><math>\phi (x\mbox{,}t\mbox{;}q){\phi }_{xxx}(x\mbox{,}t\mbox{;}q)-</math><math>\phi (x\mbox{,}t\mbox{;}q){\phi }_x(x\mbox{,}t\mbox{;}q)-</math><math>3{\phi }_x(x\mbox{,}t\mbox{;}q){\phi }_{xx}(x\mbox{,}t\mbox{;}q)\mbox{.}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" |
|}

We construct the zeroth order deformation equation

{| class="formulaSCP" style="width: 100%; text-align: center;"
|-
|
{| style="text-align: center; margin:auto;"
|-
| <math>(1-nq)L\left[\phi (x\mbox{,}t\mbox{;}q)-u_0(x\mbox{,}t)\right]=</math><math>qhH(x\mbox{,}t)N\left[D_t^{\alpha }\phi (x\mbox{,}t\mbox{;}q)\right]\mbox{.}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (19)
|}

Choosing <math display="inline">H(x\mbox{,}t)=1</math>, the ''m''th-order deformation equation is

{| class="formulaSCP" style="width: 100%; text-align: center;"
|-
|
{| style="text-align: center; margin:auto;"
|-
| <math>L\left[u_m(x\mbox{,}t)-{\chi }_m^{{_\ast}}u_{m-1}(x\mbox{,}t)\right]=</math><math>hR_m\left({\overset{\rightarrow}{u}}_{m-1}\right)\mbox{,}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (20)
|}

with initial condition for <math display="inline">m\geqslant 1\mbox{,}u_m(x\mbox{,}0)=0\mbox{,}{\chi }_m^{{_\ast}}</math> is as defined in [[#e0070|(14)]] and

{| class="formulaSCP" style="width: 100%; text-align: center;"
|-
|
{| style="text-align: center; margin:auto;"
|-
| <math>R_m\left({\overset{\rightarrow}{u}}_{m-1}\right)=D_t^{\alpha }u_{m-1}-</math><math>\sum_{k=0}^{m-1}u_k{\left(u_{m-1-k}\right)}_{xxx}-</math><math>\sum_{k=0}^{m-1}u_k{\left(u_{m-1-k}\right)}_x-3\sum_{k=0}^{m-1}{\left(u_k\right)}_x{\left(u_{m-1-k}\right)}_{xx}\mbox{.}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (21)
|}

So, the solution to Eq. [[#e0075|(15)]] for <math display="inline">m\geqslant 1</math> becomes

{| class="formulaSCP" style="width: 100%; text-align: center;"
|-
|
{| style="text-align: center; margin:auto;"
|-
| <math>u_m(x\mbox{,}t)={\chi }_m^{{_\ast}}u_{m-1}+{hI}_t^{\alpha }\left[R_m\left({\overset{\rightarrow}{u}}_{m-1}\right)\right]\mbox{.}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (22)
|}

Therefore, we obtain components of the solution using q-HAM successively as follows:

{| class="formulaSCP" style="width: 100%; text-align: center;"
|-
|
{| style="text-align: center; margin:auto;"
|-
| <math>u_1(x\mbox{,}t)={\chi }_1^{{_\ast}}u_0+{hI}_t^{\alpha }\left[D_t^{\alpha }u_0-\right. </math><math>\left. u_0{\left(u_0\right)}_{xxx}-3{\left(u_0\right)}_x{\left(u_0\right)}_{xx}-\right. </math><math>\left. u_0{\left(u_0\right)}_x\right]</math>
|-
|<math>=\frac{2{hc}^2}{3}sin\left(\frac{x}{2}\right)\frac{t^{\alpha }}{\Gamma (1+\alpha )}</math>
|-
|<math>u_2(x\mbox{,}t)={\chi }_2^{{_\ast}}u_1+{hI}_t^{\alpha }\left[D_t^{\alpha }u_1-\right. </math><math>\left. u_0{\left(u_1\right)}_{xxx}-u_1{\left(u_0\right)}_{xxx}-\right. </math><math>\left. 3\left({\left(u_0\right)}_x{\left(u_1\right)}_{xx}+\right. \right. </math><math>\left. \left. {\left(u_1\right)}_x{\left(u_0\right)}_{xx}\right)\right]</math>
|-
|<math>-{hI}_t^{\alpha }\left[u_0{\left(u_1\right)}_x+u_1{\left(u_0\right)}_x\right]</math>
|-
|<math>=\frac{2(n+h){hc}^2}{3}sin\left(\frac{x}{2}\right)\frac{t^{\alpha }}{\Gamma (1+\alpha )}+</math><math>\frac{h^2c^3}{3}cos\left(\frac{x}{2}\right)\frac{t^{2\alpha }}{\Gamma (1+2\alpha )}</math>
|-
|<math>+\frac{4h^2c^3}{9}(cos(x)-1)\frac{t^{2\alpha }}{\Gamma (1+2\alpha )}\mbox{.}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (23)
|}

Following the same approach, <math display="inline">u_m(x\mbox{,}t)</math> for <math display="inline">m=3\mbox{,}4\mbox{,}5\mbox{,}\ldots </math> can be obtained using Mathematica.

Then the series solution expression by q-HAM can be written in the form of


{| class="formulaSCP" style="width: 100%; text-align: center;"
|-
|
{| style="text-align: center; margin:auto;"
|-
| <math>u(x\mbox{,}t\mbox{;}n\mbox{;}h)=-\frac{8c}{3}{cos}^2\left(\frac{x}{4}\right)+</math><math>\sum_{i=1}^{\infty }u_i(x\mbox{,}t\mbox{;}n\mbox{;}h){\left(\frac{1}{n}\right)}^i\mbox{.}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (24)
|}

Eq. [[#e0120|(24)]] is an appropriate solution to the problem [[#e0075|(15)]] in terms of convergence parameter ''h'' and ''n''.

==5. Numerical results==

In this section, we give some numerical results using series solution obtained above. Comparison is made with exact solution for a special case using the 2-term series solution. We also show the graph displaying the best choice of h for fast convergence and the effects of different fractional order <math display="inline">\alpha </math> on the solution obtained.

===5.1. Approximate solution vs exact solution===

Exact solution is known in the case of <math display="inline">\alpha =1</math> and so we present the numerical result (2-term series solution) obtained by the q-Homotopy Analysis Method and the exact solution of Eq. [[#e0075|(15)]] under some conditions.

==Remark 5.1. ==

It should be noted that we have used only 2-term of the series solution obtained by our method to make [[#f0005|Fig. 1]] as against 5-term series solution obtained by the variational iterative method and homotopy perturbation method used in [[#b0035|[7]]]. [[#f0005|Fig. 1]] shows a perfect match with exact solution. This shows the effectiveness of the homotopy analysis method over other analytical methods due to the ability to control or choose appropriately the auxiliary parameter ''h''.



{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;"
|-
|

[[Image:draft_Content_854382069-1-s2.0-S1110016816000673-gr1.jpg|center|368px|q-HAM solution plot of u for h=-1.02,n=1,c=1 and α=1 against exact solution.]]

|-
| 

Figure 1.

q-HAM solution plot of ''u '' for <math display="inline">h=-1.02\mbox{,}n=1\mbox{,}c=1</math> and <math display="inline">\alpha =1</math> against exact solution.


|}

===5.2. The ''h''-curve===

The question that comes to mind when following this method of solution is how does one choose the auxiliary parameter ''h'' to get a good approximate solution. The answer is in the ''h''-curve. Apparently, our choice in the plots can be seen directly from the graph, the range of which is by drawing a horizontal line on the curve parallel to ''x''-axis. [[#f0010|Fig. 2]] is made with <math display="inline">c=1\mbox{,}n=1</math> and <math display="inline">\alpha =1</math>.



{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;"
|-
|

[[Image:draft_Content_854382069-1-s2.0-S1110016816000673-gr2.jpg|center|369px|q-HAM solution plot against h with fixed x=π13 and t=0.01.]]

|-
| 

Figure 2.

q-HAM solution plot against ''h '' with fixed <math display="inline">x=\frac{\pi }{13}</math> and <math display="inline">t=0.01</math>.


|}

===5.3. Solution plots with different values of <math display="inline">\alpha </math>===

Here, we give the solution plots of the 2-term series solution given by our method. This shows the effect of the fractional order <math display="inline">\alpha </math> on the obtained solution (see [[#f0015|Figure 3]] and [[#f0020|Figure 4]]).



{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;"
|-
|

[[Image:draft_Content_854382069-1-s2.0-S1110016816000673-gr3.jpg|center|366px|q-HAM solution plot of Eq. (15) for different values of α with fixed ...]]

|-
| 

Figure 3.

q-HAM solution plot of Eq. [[#e0075|(15)]] for different values of <math display="inline">\alpha </math> with fixed <math display="inline">x=\pi /13\mbox{,}c=1\mbox{,}h=-1.02</math> and <math display="inline">n=1</math>.


|}



{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;"
|-
|

[[Image:draft_Content_854382069-1-s2.0-S1110016816000673-gr4.jpg|center|369px|q-HAM solution plot of Eq. (15) for different values of α with fixed ...]]

|-
| 

Figure 4.

q-HAM solution plot of Eq. [[#e0075|(15)]] for different values of <math display="inline">\alpha </math> with fixed <math display="inline">t=0.5\mbox{,}c=1\mbox{,}h=-1.02</math> and <math display="inline">n=1</math>.


|}

==6. Conclusion==

The main achievement of this paper is the demonstration of the successful application of the q-Homotopy Analysis Method (q-HAM) in obtaining an analytical solution of the time fractional Rosenau–Hyman equation. Our results confirm that the method is very effective for handling solutions of a class of non-linear partial differential equations of fractional order system.

The comparisons made with the exact solution and other analytical methods such as a variational iterative method (VIM) and homotopy perturbation method (HPM), see [[#b0035|[7]]], expose the accuracy nature of q-HAM in the sense that just two terms of the series solution are needed to obtain better approximation using the auxiliary parameter ''h''. This method is a potential analytical method for further works in strongly nonlinear fractional differential equations.

==References==

<ol style='list-style-type: none;margin-left: 0px;'><li>
[[#b0005|[1]]] H.-K. Pang, H.-W. Sun; Multigrid method for fractional diffusion equations; J. Comput. Phys. (2011) [http://dx.doi.org/10.1016/j.jcp.2011.10.005 http://dx.doi.org/10.1016/j.jcp.2011.10.005]</li>
<li>
[[#b0010|[2]]] Md. N. Alam, M.G. Hafez, M.A. Akbar; H-Or- Roshid; exact traveling wave solutions to the <math display="inline">\left(3+1\right)</math>-dimensional mKdV–ZK and the <math display="inline">\left(2+1\right)</math>-dimensional Burgers equations via exp<math display="inline">\left(-\Phi (\eta )\right)</math>-expansion method; Alexandria Eng. J., 54 (2015), pp. 635–644</li>
<li>
[[#b0015|[3]]] S. Gupta, D. Kumar, J. Singh; Analytical solutions of convection-diffusion problems by combining Laplace transform method and homotopy perturbation method; Alexandria Eng. J., 54 (2015), pp. 645–651</li>
<li>
[[#b0020|[4]]] J.H. He; Semi-inverse method of establishing generalized variational principles for fluid mechanics with emphasis on turbomachinery aerodynamics; Int. J. Turbo Jet Engines, 14 (1) (1997), pp. 23–28</li>
<li>
[[#b0025|[5]]] J.-H. He; Some applications of nonlinear fractional differential equations and their approximations; Bull. Sci. Technol. Soc., 15 (2) (1999), pp. 86–90</li>
<li>
[[#b0030|[6]]] I. Podlubny; Fractional differential equations; Mathematics in Science and Engineering, vol. 198, Academic Press, San Diego, Calif, USA (1999)</li>
<li>
[[#b0035|[7]]] R. Yulita, M.S.M. Nooorani; Solving the fractional Rosenau–Hyman equation via variational iteration method and homotopy perturbation method; Int. J. Diff. Eq. (2012), p. 14</li>
<li>
[[#b0040|[8]]] J. Garralón, F.R. Villatoro; Dissipative perturbations for the K(n,n) Rosenau–Hyman equation; Commun. Nonlinear Sci. Numer. Simulat., 17 (2012), pp. 4642–4648</li>
<li>
[[#b0045|[9]]] M. Dehghan, J. Manafian, A. Saadatmandi; Application of semi-analytical methods for solving the Rosenau–Hyman equation arising in the pattern formation in liquid drops; Int. J. Numer. Meth. Heat Fluid Flow, 22 (6) (2012), pp. 777–790</li>
<li>
[[#b0050|[10]]] M.A. El-Tawil, S.N. Huseen; The Q-homotopy analysis method (Q-HAM); Int. J. Appl. Math. Mech., 8 (15) (2012), pp. 51–75</li>
<li>
[[#b0055|[11]]] O.S. Iyiola, F.D. Zaman; A fractional diffusion equation model for cancer tumor; AIP Adv., 4 (2014), p. 107121</li>
<li>
[[#b0060|[12]]] O.S. Iyiola, G.O. Ojo; Analytical solutions of time-fractional models for homogeneous gadner equation and non-homogeneous differential equations; Ain Shams Eng. J., 5 (2014), pp. 999–1004</li>
<li>
[[#b0065|[13]]] O.S. Iyiola, M.E. Soh, C.D. Enyi; Generalised homotopy analysis method (q-HAM) for solving foam drainage equation of time fractional type; Math. Eng. Sci. Aerospace, 4 (4) (2013), pp. 105–116</li>
<li>
[[#b0070|[14]]] O.S. Iyiola, G.O. Ojo; On the analytical solution of Fornberg–Whitham equation with the new fractional derivative; Pramana J. Phys., 85 (4) (2015), pp. 567–575</li>
<li>
[[#b0075|[15]]] M.S. Mohamed, Y.S. Hamed; Solving the convection-diffusion equation by means of the optimal q-homotopy analysis method (Oq-HAM); Results Phys., 6 (2016), pp. 20–25</li>
<li>
[[#b0080|[16]]] O.S. Iyiola; On the solutions of non-linear time-fractional gas dynamic equations: an analytical approach; Int. J. Pure Appl. Math., 98 (4) (2015), pp. 491–502</li>
<li>
[[#b0085|[17]]] O.S. Iyiola; Exact and approximate solutions of fractional diffusion equations with fractional reaction terms; Prog. Fract. Diff. Appl., 2 (1) (2016), pp. 21–30</li>
<li>
[[#b0090|[18]]] S.-J. Liao; An approximate solution technique not depending on small parameters: a special example; Int. J. Non-linear Mech., 30 (3) (1995), pp. 371–380</li>
<li>
[[#b0095|[19]]] P.A. Clarkson, E.L. Mansfield, T.J. Priestley; Symmetries of a class of nonlinear third-order partial differential equations; Math. Comput. Modell., 25 (8–9) (1997), pp. 195–212</li>
</ol>

Samuel-Iyiola et al 2016a - Revision history

Scipediacontent: Scipediacontent moved page Draft Content 854382069 to Samuel-Iyiola et al 2016a

Scipediacontent: Created page with "==Abstract== A relatively new method called q-Homotopy Analysis Method (q-HAM) is adopted in this paper to obtain an analytical solution of the time fractional Rosenau–Hyma..."

← Older revision	Revision as of 08:48, 12 April 2017
(No difference)