Advances in Mathematical Physics
Volume 2016 (2016), Article ID 6021462, 11 pages
http://dx.doi.org/10.1155/2016/6021462

Research Article

Iterative Methods for Solving the Fractional Form of Unsteady Axisymmetric Squeezing Fluid Flow with Slip and No-Slip Boundaries

A. A. Hemeda and E. E. Eladdad

Department of Mathematics, Faculty of Science, Tanta University, Tanta 31527, Egypt

Received 15 January 2016; Accepted 11 April 2016

Academic Editor: Giorgio Kaniadakis

Copyright © 2016 A. A. Hemeda and E. E. Eladdad. 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.

Abstract.

An unsteady axisymmetric flow of nonconducting, Newtonian fluid squeezed between two circular plates is proposed with slip and no-slip boundaries. Using similarity transformation, the system of nonlinear partial differential equations of motion is reduced to a single fourth-order nonlinear ordinary differential equation. By using the basic definitions of fractional calculus, we introduced the fractional order form of the fourth-order nonlinear ordinary differential equation. The resulting boundary value fractional problems are solved by the new iterative and Picard methods. Convergence of the considered methods is confirmed by obtaining absolute residual errors for approximate solutions for various Reynolds number. The comparisons of the solutions for various Reynolds number and various values of the fractional order confirm that the two methods are identical and therefore are suitable for solving this kind of problems. Finally, the effects of various Reynolds number on the solution are also studied graphically.

1. Introduction

The squeezing of an incompressible viscous fluid between two parallel plates is a fundamental type of flow that is frequently observed in many hydrodynamical tools and machines. Compression and injection molding, polymer processing, and modeling of lubrication systems are some practical examples of squeezing flows where their usage is found. The modeling and analysis of squeezing flow has been started in the nineteenth century and continues to receive significant attention due to its vast applications areas in biophysical and physical sciences. The first work in squeezing flows was laid down by Stefan [1] who developed an ad hoc asymptotic solution of Newtonian fluids. An explicit solution of the squeeze flow, considering inertial terms, has been established by Thorpe and Shaw [2]. However, P. S. Gupta and A. S. Gupta proved that the solution given in [2] fails to satisfy boundary conditions [3]. Verma [4] and Singh et al. [5] have established numerical solutions of the squeezing flow between parallel plates. Leider and Byron Bird performed theoretical analysis of power-law fluid between parallel disks [6]. Qayyum et al. present in [7] analysis of unsteady axisymmetric flow of nonconducting, Newtonian fluid squeezed between two circular plates with slip and no-slip boundaries using OHAM and in [8] the authors model and analyse the unsteady axisymmetric flow of nonconducting Newtonian fluid squeezed between two circular plates passing through porous medium channel with slip boundary condition using HPM. Furthermore, in [9] analytical solutions for squeeze flow with partial wall slip are introduced by Laun et al., while in [10] Ullah et al. present approximation of first-grade MHD squeezing fluid flow with slip boundary condition using DTM and OHAM. In [11] analytical solution of squeezing flow between two circular plates is presented by Rashidi et al. Islam et al. [12] studied Newtonian squeezing fluid flow in a porous medium channel, Siddiqui et al. [13] investigated the unsteady squeezing flow of viscous fluid with magnetic field, Marinca et al. [14] applied an optimal homotopy asymptotic method to the steady flow of a fourth-grade fluid past a porous plate, and Idrees et al. [15] applied the optimal homotopy asymptotic method to squeezing flow. For more details see [16, 17].

A large number of perturbation methods, which can solve nonlinear boundary value problems analytically, are discussed in the literature. These methods, however, have a limitation of assuming small or large parameters. Recently, Daftardar-Gejji and Jafari have proposed a new technique, called new iterative method (NIM), for solving linear and nonlinear boundary value problems. The method minimize the limitations that are usually associated with perturbation methods and at the same time it takes full advantages of the traditional perturbation methods. This method has proven useful for solving a variety of nonlinear equations such as algebraic equations, integral equations, ordinary and partial differential equations of integer and fractional order, and systems of equations as well. The NIM is used by Hemeda in [18] for solving fractional partial differential equations and in [19] for solving the th-order integrodifferential equations, but in [20] he used the NIM for solving fractional physical differential equations. Also, this method is used by Daftardar-Gejji and Jafari in [21] for solving nonlinear fractional equations and by Bhalekar and Daftardar-Gejji in [22] it is used for solving partial differential equations while in [23] it is used for solving evolution equations. In the late of the nineteenth century, Emile Picard proposed a good method called Picard iteration method (or shortly) Picard method (PM), which is used widely by researches for solving linear and nonlinear boundary value problems. This method has not any of the difficulties presented in perturbation methods and takes full advantages of the traditional perturbation methods as well as NIM. The PM is used by Youssef and El-Arabawy in [24] for initial value problems and by Ibijola and Adegboyegun in [25] for solving nonlinear differential equations. Also, this method is used by Hemeda in [26] for solving the fractional gas dynamics and coupled Burgers’ equations. In this method, the solution takes the form of a rapidly convergent series with easily computable components.

The importance of the study of the fractional forms of the differential equations is due to their wide appearing in many of the mathematical, physical, and chemical problems. So the aim of this work is to continue in this study by preparing and using the NIM and PM where there is not any of the above-mentioned difficulties in the perturbation methods for solving the fractional order form of an unsteady axisymmetric squeezing fluid flow between two circular plates with slip and no-slip boundaries. Also the effects of various of Reynolds number and fractional order on the solution are studied tabularly and graphically.

2. Formulation of the Problem

In this section, the unsteady axisymmetric squeezing flow of incompressible Newtonian fluid with density , viscosity , and kinematic viscosity , squeezed between two circular plates having speed , is considered with a fractional form. It is assumed that, at any time , the distance between the two circular plates is . Also it is assumed that -axis is the central axis of the channel while -axis is taken as normal to it. Plates move symmetrically with respect to the central axis while the flow is axisymmetric about . The longitudinal and normal velocity components in radial and axial directions are and , respectively. For more physical explanation and details, see [7, 8].

The equations of motion are

where is the vorticity function and is the generalized pressure.

The boundary conditions on and are

where is the velocity of the plates. The boundary conditions in (4) are due to symmetry at and no-slip at the upper plate when . If we introduce the dimensionless parameter

(1), (2), and (3) transforms to

The boundary conditions on and are

By eliminating the generalized pressure between (7) and (8), we obtain

where is the Laplacian operator.

Defining velocity components as [3]

we see that (6) is identically satisfied and (10) becomes

where

Here both and are functions of but we consider and as constants for similarity solution. Since , integrate first equation of (13), and we get

where and are constants. The plates move away from each other symmetrically with respect to when and . Also the plates approach each other and squeezing flow exists with similar velocity profiles when , , and . From (13) and (14), it follows that . Then, (12) becomes

Using (9) and (11), we determine the boundary conditions in case of no-slip and slip at the upper plate as follows:

3. Fractional Calculus

In this section, we mention some basic definitions of fractional calculus which are used in the present work.

Definition 1. The Riemann-Liouville fractional integral operator of order , of a function and , is defined as [27]

For the Riemann-Liouville fractional integral operator, , we obtain

Definition 2. The fractional derivative of in the Caputo sense is defined as [28]

For the Caputo fractional derivative operator, , we obtain

For the Riemann-Liouville fractional integral and Caputo fractional derivative operator of order , we have the following relation:

Remark 3. According to the previous fractional calculus, (15) can be rewritten in the following fractional order form:

4. Analysis of the Considered Methods

In this section, we discuss the considered methods with preparing them for solving any fractional differential equation.

4.1. New Iterative Method (NIM)

To illustrate the basic idea of this method, we consider the following general functional equation [18–23]:

where is a nonlinear operator from a Banach space and is a known function (element) of a Banach space . We are looking for a solution of (23) having the series form:

The nonlinear operator can be decomposed as

From (24) and (25), (23) is equivalent to

The required solution for (23) can be obtained recurrently from the recurrence relation:

Then,

The r-term approximate solution of (23) is given by

4.2. Solving General Fractional Differential Equation by NIM

To solve any fractional differential equation of arbitrary order , we consider the following general fractional differential equation of order :

subject to the initial values

where is a linear operator, is a nonlinear operator, is a nonhomogeneous term, and is the fractional differential operator of order . In view of the fractional integral operators, the initial value fractional problem (30a) and (30b) is equivalent to the fractional integral equation:

where , , and is the inverse of . The required solution for (31) and hence for (30a) and (30b) can be obtained recurrently from the recurrence relation (27).

4.3. Picard Method (PM)

To illustrate the basic idea of this method, we consider the following general fractional differential equation of arbitrary order [24–26]:

where is the fractional differential operator of order . In view of the fractional integral operators, the initial value fractional problem (32a) and (32b) is equivalent to the fractional integral equation:

where , , and is the inverse of . The required solution for (33) which is also the solution for (32a) and (32b) can be obtained as the limit of a sequence of functions generated by the recurrence relation:

where .

5. Applications

In this section, we illustrate the application of the two considered methods to solve the nonlinear fractional order ordinary differential equation (22) subject to the boundary conditions (16a) and (16b).

5.1. NIM

Using (22), (16a), and (16b), the initial value fractional order problem

according to (31), is equivalent to the fractional integral equation:

where and are constants. Let . Therefore, according to (27), we can obtain the following first few components of the new iterative solution for (35)

and so on. In the same manner, the rest of components can be obtained. The 4-term solution for (35) in series form is given by

In the special case, , (38) becomes

Using the boundary conditions in (16a) with the initial conditions in (35), the unknowns and for fixed values of in (39) can be easily determined. In case of no-slip boundary, then and . For , the 4-term solution, obtained by the NIM in (39), is therefore

Similarly, using the boundary conditions in (16b) with the initial conditions in (35), in case of slip boundary, then and . For , the 4-term solution becomes

5.2. PM

Also the initial value fractional problem (35), according to (33) and (34), is equivalent to the fractional integral equation:

Therefore, according to (34), we can obtain the following first few components of Picard solution for (35):

and so on. In the same manner the rest of components can be obtained. The 4-order term solution for (35), in series form, is given by

In the special case, , (44) becomes

Also, in case of no-slip boundary where , the 4-order term solution is

and, in case of slip boundary, the 4-order term solution becomes

From the previous results for (35), obtained by the two considered methods, it is clear that the approximate solution for (35) obtained by PM in (44), (45), (46), and (47) is the same approximate solution as obtained by NIM in (38), (39), (40), and (41). Therefore, the two methods are identical in solving this problem and hence the two methods are suitable for solving this kind of problems.

The residual error of the problem is

where is the 4-term approximate solution in (38) or (44) for (35).

If , then will be the exact solution. However, this usually does not occur in nonlinear problems.

It is clear from the obtained results that the above-considered methods minimize the limitations of the ordinary perturbation methods. In the same time, these methods take full advantages of the traditional perturbation methods. Therefore, these methods are powerful methods for solving the nonlinear fractional order differential equations.

6. Numerical Results and Discussion

In this work, an unsteady axisymmetric flow of nonconducting, incompressible Newtonian fluid squeezed between two circular plates is considered. The resulting nonlinear fractional order boundary value problem is solved analytically in cases of no-slip and slip boundaries using NIM and PM.

Tables 1, 2, and 3 present the solutions in case of no-slip boundary for different values of the fractional order along with absolute residual errors at for fixed value of the Reynolds number in each table, while Tables 4, 5, and 6 are in case of slip boundary and slip parameter . The results in the six tables indicate that, as and , . Therefore, the approximate solutions converge to the exact solution for the considered problem.

Table 1: Solutions for different values of and at , in case of no-slip boundary.

0.0
0.0
0.0
0.0
0.0
0.0
0.0

0.1
0.149500
0.149500
0.149500
0.149500
0.149500
5.63438⁻¹⁵

0.2
0.296005
0.296004
0.296003
0.296002
0.296001
7.21923⁻¹³

0.3
0.436536
0.436527
0.436520
0.436515
0.436511
1.14405⁻¹¹

0.4
0.568140
0.568108
0.568083
0.568064
0.568049
5.97444⁻¹¹

0.5
0.687900
0.687816
0.687750
0.687697
0.687654
2.94232⁻¹¹

0.6
0.792945
0.792761
0.792611
0.792490
0.792392
2.46997⁻⁹

0.7
0.880455
0.880099
0.879805
0.879563
0.879364
2.03243⁻⁸

0.8
0.947673
0.947045
0.946519
0.946080
0.945714
1.09632⁻⁷

0.9
0.991913
0.990880
0.990004
0.989263
0.988637
4.67357⁻⁷

1.0
1.010560
1.008960
1.00758
1.006400
1.005390
1.69343⁻⁶

Table 2: Solutions for different values of and at , in case of no-slip boundary.

0.0
0.0
0.0
0.0
0.0
0.0
0.0

0.1
0.149501
0.149500
0.149500
0.149500
0.149500
4.55608⁻¹³

0.2
0.296016
0.296012
0.296008
0.296006
0.296004
5.84710⁻¹¹

0.3
0.436608
0.436581
0.436561
0.436546
0.436534
9.26826⁻¹⁰

0.4
0.568418
0.568324
0.568250
0.568193
0.568148
4.85819⁻⁹

0.5
0.688694
0.688445
0.688247
0.688088
0.687963
1.90307⁻⁹

0.6
0.794814
0.794270
0.793827
0.793467
0.793175
1.93785⁻⁷

0.7
0.884309
0.883263
0.882394
0.881676
0.881085
1.59180⁻⁶

0.8
0.954885
0.953051
0.951504
0.950206
0.949120
8.52601⁻⁶

0.9
1.004440
1.001450
0.998887
0.996708
0.994859
3.59984⁻⁵

1.0
1.03108
1.026470
1.022470
1.019020
1.01606
1.28930⁻⁴

Table 3: Solutions for different values of and at , in case of no-slip boundary.

0.0
0.0
0.0
0.0
0.0
0.0
0.0

0.1
0.149501
0.149501
0.149500
0.149500
0.149500
3.51580⁻¹²

0.2
0.296027
0.296019
0.296014
0.296010
0.296007
4.51141⁻¹⁰

0.3
0.436680
0.436635
0.436602
0.436576
0.436557
7.15258⁻⁹

0.4
0.568695
0.568538
0.568416
0.568321
0.568247
3.76300⁻⁸

0.5
0.689480
0.689070
0.688741
0.688479
0.688270
1.10362⁻⁸

0.6
0.796656
0.795763
0.795033
0.794438
0.793954
1.44776⁻⁶

0.7
0.888089
0.886381
0.884956
0.883773
0.882795
1.18733⁻⁵

0.8
0.961917
0.958942
0.956417
0.954287
0.952498
6.31480⁻⁵

0.9
1.016570
1.011760
1.007600
1.004050
1.001010
2.64053⁻⁴

1.0
1.050790
1.043440
1.037000
1.031420
1.026570
9.34686⁻⁴

Table 4: Solutions for different values of and at and in case of slip boundary.

0.0
0.0
0.0
0.0
0.0
0.0
0.0

0.1
0.075250
0.075250
0.075250
0.075250
0.075250
7.27474⁻¹³

0.2
0.151997
0.151997
0.151998
0.151999
0.151999
9.13051⁻¹¹

0.3
0.231727
0.231732
0.231737
0.231740
0.231743
1.50939⁻⁹

0.4
0.315910
0.315930
0.315946
0.315959
0.315968
1.07916⁻⁸

0.5
0.405994
0.406047
0.406090
0.406124
0.406151
4.84153⁻⁸

0.6
0.503398
0.503514
0.503609
0.503686
0.503749
1.60795⁻⁷

0.7
0.609512
0.609735
0.609921
0.610074
0.610200
4.31515⁻⁷

0.8
0.725690
0.726081
0.726410
0.726685
0.726916
9.85311⁻⁷

0.9
0.853250
0.853886
0.854429
0.854889
0.855278
1.97750⁻⁶

1.0
0.993468
0.994447
0.995291
0.996016
0.996637
3.56318⁻⁶

Table 5: Solutions for different values of and at and in case of slip boundary.

0.0
0.0
0.0
0.0
0.0
0.0
0.0

0.1
0.075250
0.075250
0.075250
0.075250
0.075250
5.89262⁻¹¹

0.2
0.151990
0.151992
0.151995
0.151996
0.151997
7.39595⁻⁹

0.3
0.231680
0.231697
0.231711
0.231720
0.231728
1.22281⁻⁷

0.4
0.315732
0.315792
0.315839
0.315876
0.315905
8.74585⁻⁷

0.5
0.405489
0.405646
0.405772
0.405873
0.405953
3.92681⁻⁶

0.6
0.502219
0.502559
0.502838
0.503065
0.503250
1.30608⁻⁵

0.7
0.607104
0.607750
0.608290
0.608738
0.609110
3.51392⁻⁵

0.8
0.721230
0.722348
0.723298
0.724099
0.724775
8.05617⁻⁵

0.9
0.845586
0.847384
0.848935
0.850263
0.851400
1.62693⁻⁴

1.0
0.981059
0.983786
0.986174
0.988246
0.990047
2.95870⁻⁴

Table 6: Solutions for different values of and at and in case of slip boundary.

0.0
0.0
0.0
0.0
0.0
0.0
0.0

0.1
0.075250
0.075250
0.075250
0.075250
0.075250
4.54678⁻¹⁰

0.2
0.151983
0.151987
0.151991
0.151994
0.151995
5.70694⁻⁸

0.3
0.231634
0.231663
0.231684
0.231701
0.231713
9.43682⁻⁷

0.4
0.315556
0.315655
0.315733
0.315794
0.315841
6.75191⁻⁶

0.5
0.404995
0.405251
0.405458
0.405624
0.405756
3.03391⁻⁶

0.6
0.501075
0.501625
0.502079
0.502451
0.502756
1.01056⁻⁴

0.7
0.604788
0.605823
0.606696
0.607426
0.608035
2.72554⁻⁴

0.8
0.716984
0.718754
0.720276
0.721570
0.722671
6.27313⁻⁴

0.9
0.838373
0.841182
0.843641
0.845764
0.847608
1.27436⁻³

1.0
0.969523
0.973721
0.977461
0.980735
0.983639
2.33762⁻³

Figures 1 and 2 indicate the approximate solutions for different values of at in cases of no-slip and slip boundaries, respectively, with . Figures 3 and 4 show the residual errors at for different values of in cases of no-slip and slip boundaries, respectively, with . In addition to the above figures, the effect of Reynolds number on velocity profiles in case of no-slip boundary is shown in Figure 5. In these profiles, we varied as and observed that the normal velocity is increased with increasing (Figure 5(a)). It is also noted that the normal velocity monotonically increases from to for fixed value of at a given time. Figure 5(b) describes the effect of on the longitudinal velocity in case of no-slip boundary. It is clear that this component of velocity decreases near the wall but increases near the central axis of the channel.

Figure 1: Solution for at in case of no-slip boundary.

Figure 2: Solution for at , in case of slip boundary.

Figure 3: Residual for (continuous curve), 0.3 (dashed curve), and 0.5 (dotted curve) at in case of no-slip boundary.

Figure 4: Residual for (continuous curve), 0.3 (dashed curve), and 0.5 (dotted curve) at , in case of slip boundary.

(a)

(b)

Figure 5: (a) Velocity profiles (normal velocity) for (continuous curve), 0.3 (dashed curve), and 0.5 (dotted curve) at in case of no-slip boundary. (b) Velocity profiles (longitudinal velocity) for (continuous curve), 0.3 (dashed curve), and 0.5 (dotted curve) at in case of no-slip boundary.

The effect of Reynolds number on velocity profiles in case of slip boundary is depicted in Figure 6. In these profiles, we varied as for fixed slip parameter . It is noted that the normal velocity decreases as increases (Figure 6(a)). It is also noted that the longitudinal velocity decreases near the central axis of the channel but increases near the walls when increases (Figure 6(b)).

(a)

(b)

Figure 6: (a) Velocity profiles (normal velocity) for (continuous curve), 0.3 (dashed curve), and 0.5 (dotted curve) at , in case of slip boundary. (b) Velocity profiles (longitudinal velocity) for (continuous curve), 0.3 (dashed curve), and 0.5 (dotted curve) at , in case of slip boundary.

7. Conclusion

In this work, an analytical solution for an unsteady axisymmetric squeezing flow of incompressible Newtonian fluid between two circular plates in fractional form is obtained using the NIM and PM in cases of no-slip and slip boundaries. Analysis of the residual errors confirms that the NIM and PM are identical and efficient schemes. Convergence of the considered methods is confirmed by absolute residual errors for different values of the Reynolds number . Therefore, we concluded that the considered methods can be effectively used in various fields of science and engineering as they give better results in terms of accuracy.

Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.

References

M. J. Stefan, “Versuch über die scheinbare adhäsion,” Sitzungsberichte/Österreichische Akademie der Wissenschaften in Wien Mathematisch-Naturwissenschaftliche Klasse, vol. 69, pp. 713–721, 1874.

J. F. Thorpe and W. A. Shaw, Developments in Theoretical and Applied Mechanics, Pergamon Press, Oxford, UK, 1976.

P. S. Gupta and A. S. Gupta, “Squeezing flow between parallel plates,” Wear, vol. 45, no. 2, pp. 177–185, 1977.

R. L. Verma, “A numerical solution for squeezing flow between parallel channels,” Wear, vol. 72, no. 1, pp. 89–95, 1981.

P. Singh, V. Radhakrishnan, and K. A. Narayan, “Squeezing flow between parallel plates,” Ingenieur-Archiv, vol. 60, no. 4, pp. 274–281, 1990.

P. J. Leider and R. Byron Bird, “Squeezing flow between parallel disks. I. Theoretical analysis,” Industrial and Engineering Chemistry Fundamentals, vol. 13, no. 4, pp. 336–341, 1974.

M. Qayyum, H. Khan, M. T. Rahim, and I. Ullah, “Analysis of unsteady axisymmetric squeezing fluid flow with slip and no-slip boundaries using OHAM,” Mathematical Problems in Engineering, vol. 2015, Article ID 860857, 11 pages, 2015.

M. Qayyum, H. Khan, M. T. Rahim, and I. Ullah, “Modeling and analysis of unsteady axisymmetric squeezing fluid flow through porous medium channel with slip boundary,” PLoS ONE, vol. 10, no. 3, article e0117368, 2015.

H. M. Laun, M. Rady, and O. Hassager, “Analytical solutions for squeeze flow with partial wall slip,” Journal of Non-Newtonian Fluid Mechanics, vol. 81, no. 1-2, pp. 1–15, 1999.

I. Ullah, H. Khan, and M. T. Rahim, “Approximation of first grade MHD squeezing fluid flow with slip boundary condition using DTM and OHAM,” Mathematical Problems in Engineering, vol. 2013, Article ID 816262, 9 pages, 2013.

M. M. Rashidi, A. M. Siddiqui, and M. T. Rastegari, “Analytical solution of squeezing flow between two circular plates,” International Journal for Computational Methods in Engineering Science and Mechanics, vol. 13, no. 5, pp. 342–349, 2012.

S. Islam, H. Khan, I. A. Shah, and G. Zaman, “An axisymmetric squeezing fluid flow between the two infinite parallel plates in a porous medium channel,” Mathematical Problems in Engineering, vol. 2011, Article ID 349803, 10 pages, 2011.

A. M. Siddiqui, S. Irum, and A. R. Ansari, “Unsteady squeezing flow of a viscous MHD fluid between parallel plates, a solution using the homotopy perturbation method,” Mathematical Modelling and Analysis, vol. 13, no. 4, pp. 565–576, 2008.

V. Marinca, N. Herişanu, C. Bota, and B. Marinca, “An optimal homotopy asymptotic method applied to the steady flow of a fourth-grade fluid past a porous plate,” Applied Mathematics Letters, vol. 22, no. 2, pp. 245–251, 2009.

M. Idrees, S. Islam, S. Haq, and S. Islam, “Application of the optimal homotopy asymptotic method to squeezing flow,” Computers & Mathematics with Applications, vol. 59, no. 12, pp. 3858–3866, 2010.

R. J. Grimm, “Squeezing flows of Newtonian liquid films an analysis including fluid inertia,” Applied Scientific Research, vol. 32, no. 2, pp. 149–166, 1976.

H. Khan, S. Islam, J. Ali, and I. Ali Shah, “Comparison of different analytic solutions to axisymmetric squeezing fluid flow between two infinite parallel plates with slip boundary conditions,” Abstract and Applied Analysis, vol. 2012, Article ID 835268, 18 pages, 2012.

A. A. Hemeda, “Solution of fractional partial differential equations in fluid mechanics by extension of some iterative method,” Abstract and Applied Analysis, vol. 2013, Article ID 717540, 9 pages, 2013.

A. A. Hemeda, “New iterative method: application to nth-order integro-differential equations,” INFORMATION, Japan, vol. 16, no. 6, pp. 3841–3852, 2011.

A. A. Hemeda, “New iterative method: an application for solving fractional physical differential equations,” Abstract and Applied Analysis, vol. 2013, Article ID 617010, 9 pages, 2013.

V. Daftardar-Gejji and H. Jafari, “An iterative method for solving nonlinear functional equations,” Journal of Mathematical Analysis and Applications, vol. 316, no. 2, pp. 753–763, 2006.

S. Bhalekar and V. Daftardar-Gejji, “New iterative method: application to partial differential equations,” Applied Mathematics and Computation, vol. 203, no. 2, pp. 778–783, 2008.

S. Bhalekar and V. Daftardar-Gejji, “Solving evolution equations using a new iterative method,” Numerical Methods for Partial Differential Equations, vol. 26, no. 4, pp. 906–916, 2010.

I. K. Youssef and H. A. El-Arabawy, “Picard iteration algorithm combined with Gauss-Seidel technique for initial value problems,” Applied Mathematics and Computation, vol. 190, no. 1, pp. 345–355, 2007.

E. A. Ibijola and B. J. Adegboyegun, “A comparison of Adomian's decomposition method and Picard iteration method in solving nonlinear differential equations,” Global Journal of Science Frontier Research F: Mathematics and Decision Sciences, vol. 12, Version 1.0, no. 7, 2012.

A. A. Hemeda, “Iterative methods for solving fractional gas dynamics and coupled Burgers' equations,” International Journal of Mathematical Modelling and Numerical Optimisation, vol. 6, no. 4, pp. 282–298, 2015.

I. Podlubny, Fractional Differential Equations, vol. 198 of Mathematics in Science and Engineering, Academic Press, San Diego, Calif, USA, 1999.

M. Caputo, Elasticita E Dissipazione, Zani-Chelli, Bologna, Italy, 1969.