Adomian Decomposition Method Applied to Study. Nonlinear Equations Arising in non-newtonian flows

Transcription

1 Applied Mathematical Sciences, Vol. 6,, no. 98, Adomian Decomposition Method Applied to Study Nonlinear Equations Arising in non-newtonian flows A. M. Siddiqui (a), M. Hameed (b) (a) Department of Mathematics, Pennsylvania State University, York campus, York, PA 74 98, USA. (b) Division of Mathematics & Computer Science, USCUPSTATE 8 University Way, Spartanburg, 9, USA B. M. Siddiqui (c), B. S. Babcock (a) (c) Department of Computer Science, COMSATS Institute of Information Technology, Park Road, Chak Shahzad, Islamabad, Pakistan. Abstract Adomian s decomposition method (ADM) is employed to solve nonlinear differential equations which arise in non-newtonian fluid dynamics. We study basic pipe flow problems of a third grade and 6-constant Oldroyd non-newtonian fluids. The technique of Adomian decomposition is applied with elegant results. The solutions obtained show the reliability and efficiency of this analytical method. Numerical solutions are also obtained by solving nonlinear ordinary differential equations using Chebyshev spectral method. We present a comparative study of the analytical and numerical solutions. The analytical results agree well with the numerical results, which reveal the effectiveness and convenience of the Adomian decomposition method. Keywords: Decomposition Method; Spectral method; Third grade fluid; Oldroyd fluid; Couette flow; Poiseuille flow

2 489 A. M. Siddiqui et al. Introduction In recent years, there has been great interest in the application of a recently developed decomposition method, due to Adomian [, ], in solving nonlinear ordinary and partial differential equations. A considerable amount of research work has been invested in applying this method to a wide class of linear and nonlinear ordinary differential equations, partial differential equations and integral equations. The ADM has been successfully applied to solve nonlinear equations in studying many interesting problems arising in applied sciences and engineering [, 4, 6, 7, 9,, ], and is usually characterized by its higher degree of accuracy. It provides analytical solution in the form of an infinite series in which each term can be easily determined. The ADM has been proved to be an effective and reliable method for handling differential equations, linear or nonlinear. Unlike the traditional methods, the ADM needs no discretization, linearization, spatial transformation or perturbation, thus has some significant advantages over the other traditional analytical techniques. The ADM provides an analytical solution in the form of an infinite convergent power series. The rapid convergence of the series solution obtained by ADM has been discussed by Cherrualt et al. [5] and provides insight into the characteristics and behavior of the solution as in the case with the closed form solution. Non-Newtonian fluids are of great interest to researchers because of their practical importance in engineering and industry. The classical Navier-Stokes equations have been proved inadequate to describe and capture the characteristics of complex rheological fluids as well as polymer solutions []. These kinds of fluids are generally known as non-newtonian fluids. Most of the biological and industrial fluids are non-newtonian in nature. Few examples of such fluids are blood, tomato ketchup, honey, mud, plastics and polymer solutions. The inadequacy of the classical theories to describe these complex fluids has led to the development of different new theories to study non-newtonian fluids. There are different models which have been proposed to describe the non-newtonian flow behavior. Among these, the fluids of differential type [, ] have received considerable attention. Fluid of a third grade is a subclass of fluids of differential type, which has been studied successfully in various types of flow situations [4, 5] and is known to capture the non-newtonian affects such as shear thinning or shear thinking as well as normal stresses. The constitutive equations for the third grade fluid, as in the case of most non-newtonian fluids, give rise to complicated and highly nonlinear model equations. This intrinsic nonlinearity due to the presence of normal stresses makes the problem difficult to solve analytically even in simple geometries. In the present work, we successfully apply for the first time the Adomian decomposition method to solve the problems modeling the flow of a third grade fluid

3 Adomian decomposition method 489 and 6-constant Oldroyd fluid in an axisymmetric tube. The decomposition method is seen to provide a direct scheme for solving these nonlinear problems without any need for linearization or any restrictive assumptions. Moreover, the method is found to greatly reduce the size of computational work while maintaining high accuracy. The formulas for Adomian polynomials for many forms of nonlinear terms in differential equations show their effectiveness in our problems as they do in a variety of applications in other fields. We have also investigated this problem numerically and obtained solutions using higher order Chebyshev spectral collocation method. Spectral methods have been successfully used in finding the numerical solution of various problems and in studying computational fluid dynamics problems. Spectral methods are proved to offer a superior intrinsic accuracy for derivative calculations [6]. We present numerical solution for both cases and compare analytical solutions obtained by Adomian decomposition method with numerical solutions. We present a comparative study and find that ADM solutions compare well with the numerical solutions therefore prove that the ADM is powerful and effective which can easily handle wide class of nonlinear problems.. Basic Equations The basic equations governing the flow of an incompressible fluid neglecting the thermal effects are: v =, () % % dv ρ % = ρ f + T. () dt % % % Here ρ denotes the fluid density, v the velocity vector and f the body force vector % % per unit mass. The operator d dt denotes the material time derivative and T % is the stress tensor which for a third grade fluid is given by T = pi + S % % %, () where I is the unit tensor and extra stress S is defined as % % S = μa+ αa + αa + βa + β( AA + AA) + β( tra ) A, (4) % % % % % % % % % % % μ is the coefficient of the shear viscosity, α, α, β, β and β material constants and p represents the pressure. The tensors A, A and A are given by % % % A = (grad v) + (grad v) T, (5) % % % d A = A+ A(grad v) + (grad v) T A, (6) % dt % % % % d A = A + A(grad v) + (grad v) T A. (7) % dt % % % %

4 489 A. M. Siddiqui et al The constitutive equation of a 6-constant Oldroyd fluid is given by T = pi + S % % %, (8) wherein the extra stress S can be written as % DS λ λ5 DA S + λ % + ( SA+ AS ) + [ tr( S) A] A = μ( A+ λ % + λ4a ), (9) % Dt %% % % % % % % Dt % where μ, λ, λ, λ, λ 4 and λ 5 are material constants. The convected derivative D is Dt defined as DS ds T % = % LS SL, Dt dt %% %% and the first Rivlin Ericksen tensor A can also be written as % where T = +, A L L % % % L = (grad v). % % () () (). Problem Formulation We consider steady, isothermal, laminar, fully developed flow of a non- Newtonian incompressible fluid down the axis of a long pipe of circular cross section with diameter R. Choosing cylindrical coordinates (, r θ, z), we take the axis of symmetry the pipe along the z-axis. For a unidirectional axisymmetric fully developed flow we choose the velocity filed and the extra stress, respectively, as v= (,, u( r)), () % S = S() r. (4) % % Thus the equation of continuity () is satisfied identically whereas the equation of motion (), in the absence of body forces, yields p d = + ( rsrr ), (5) r r dr p d = + ( rsrz ). (6) z r dr On differentiating (5) with respect to z and (6) with respect to r, we obtain p d d = + ( rsrr ) z r r dr, (7) p d d = + ( rsrz ) r z dr r dr. (8) From equations (7), (8) and (4) we conclude that

5 Adomian decomposition method 489 d ( rsrz ) = constant. (9) rdr Now if we substitute (9) into (6) we arrive at p dp = = constant. () z Substituting () into (6) and integrating we conclude that rdp C S = rz + r, () where C is a constant of Integration and dp is the constant pressure gradient. In order to determine the value of the constant C we use the symmetry condition S rz () =, () which yields C =. () Thus, equation () reduces to rdp Srz =. (4) In the next section, we derive the resulting differential equations for the pipe flow under consideration for both fluids. Namely, the third grade fluid and the 6-constant Oldroyd fluid. Third grade Fluid On making use of assumptions () and (4) in the constitutive equation for third grade fluid (4), we obtain the expression for stress component S rz, which after substituting into (4), yields a nonlinear governing equation given as du β du r dp + =, (5) dr μ dr μ where β = β + β. (6) The boundary condition for the pipe flow problem is the no slip condition at the walls ur () = at r = R. (7) Here we would like to point out that the β represents the non-newtonian contribution to the flow. If we set β = in (5), the equation for the famous Hagen- Poiseuille flow is recovered.

6 4894 A. M. Siddiqui et al 6-Constant Oldroyd Fluid As before, we use () and (4) in the constitutive equation (9) for 6-constant Oldroyd fluid to obtain stress component S rz. Then the expression for S rz is combined with the momentum equation (4) to yield [8] the first order non-linear differential equation du du r dp du + α r dp α =, dr dr μ dr μ (8) where α = λλ 4 ( λ4 λ)( λ+ λ5), (9) α = λλ ( λ λ)( λ+ λ5). () The associated no slip boundary condition for the pipe flow is u = at r = R. () 4. Adomian's Decomposition Method To give a brief resume of the decomposition method, we consider the following differential equation Lu+ Nu = g, () where L is a linear operator, N a non-linear operator, g the source term and u the function of independent variable y alone. If we apply L to (), the result can be written as u = f + L g L Nu, () where the function f represents the term arising from the solution of the homogeneous equation Lu = (4) which involves the constants of integration. The decomposition method consists of representing the solution u(y) by the decomposition series, u( y) = u ( y). (5) The non-linear term Nu(y) from the series of Adomian polynomials is given by n= n= Here, the Adomian polynomials A n are given by A = f( u ), n Nu( y) = A. (6) n

7 Adomian decomposition method 4895 etc. d A = u f ( u ), du d u d A = u f ( u ) + f ( u ), du! du d d u d A = u f( u ) + uu f( u ) + f( u ), (7) du du! du If we define u = f + L g, (8) and make use of (5) - (6) in (), then the later equation can be written as which gives etc. un( y) = u L An n= n=, (9) u u u u = L A, =, (4) =, =, L A L A 4 L A Finally, we summarize the result as follows: u = u, n n= where u, u, u,..., are evaluated according to (4). We point out that the convergence of this method has been discussed in [5]. In the sequel, we apply the Adomian decomposition method to two boundary value problems: The first one studies the pressure driven flow of a third grade fluid in a circular pipe described by (5) and (7) and the second one analyses the pressure driven flow of a 6-constant Oldroyd fluid in a circular pipe described by (8) and (). 4(a) Solution for a Third Grade Fluid The resulting differential equation (5) along with the boundary condition (7) for this fluid may be written as

8 4896 A. M. Siddiqui et al r dp β du Lu = μ μ dr, ur () = at r = R, (4) where is a linear invertible operator. d L =, (4) dr Applying the inverse operator L on both sides of the eqn. (4), we obtain r dp β du L Lu() r L L =. μ μ dr Equivalently, we write r dp β du ur () C L L = +, (4) μ μ dr where C is the constant of integration. We decompose u (r) and the non-linear du term Nu() r =, respectively, as follows: dr ur () = u(), r n n= Nu() r = A. The first few Adomian polynomials of the non-linear term Nu(r) are given by du = n= A, dr du du A =, dr dr du du du du A = +, dr dr dr dr du du du du du du A = + 6 +, dr dr dr dr dr dr n (44) (45) A du du4 du du du du du du du 4 = dr dr dr dr dr dr dr dr dr

9 Adomian decomposition method 4897 Now, we use equation (44) in (4) and write the sum as r dp β un() r = C+ L L An. n= μ μ (46) n= We identify the zero th component as () r dp u r = C+ L μ, (47) and the remaining components as the recurrence relation, β un+ = L An, n. (48) μ n= From equation (48), we write the first few components as follows: β u() r = L [ A], μ β u() r = L [ A], (49) μ β u() r = L [ A], μ β u4() r = L [ A] μ M The corresponding boundary condition (no slip condition), after making use of (44), becomes un ( r) = at r = R, n. (5) Consequently, we obtain dp r R u =, (5) μ 4 4 β dp r R u() r =, μ μ 4 4 (5) β dp r R u() r =, μ μ 6 6 (5) β dp r R u() r = R, μ μ 8 8 (54) 4 9 β dp r R u4() r = 55. μ μ (55) Hence the approximate solution u(r) up to order four in a series form is given by

10 4898 A. M. Siddiqui et al u() r = u() r + u() r + u() r + u() r + u4() r +L and summarizing the result, we write the fourth order approximate solution for the velocity field u(r) in the form 4 4 dp r R β dp r R ur () = μ μ μ β dp r R β dp r R + μ μ 6 6 μ μ β dp r R (56) μ μ The other quantities of practical interest in a circular pipe flow of radius R can be computed easily. The volume flow rate Q is given by the formula R Q = π u( r) rdr. (57) By substituting the velocity field from equation (56) into (57) and performing the integration, we obtain the following expression for Q, 4 dp π R Q = μ β dp π R β dp π R + (58) μ μ 6 μ μ β dp πr β dp πr L. μ μ μ μ The Hagen Poiseuille equation [6] follows from equation (5) by setting the parameter β = and (56) gives the corresponding solution, namely dp r R ur () =. μ The average velocity ur () for the circular geometry under consideration is given by the formula Q ur () = (6) π R and after making use of equation (58) for Q, the above expression takes the form, 4 dp R β dp R ur () = + μ 4 μ μ 6

11 Adomian decomposition method β dp R β dp R + (6) μ μ 8 μ μ 4 9 β dp R 55 + L. μ μ On setting the parameter β =, one can deduce the Newtonian average velocity dp R ur () =. (6) μ 4 For the sake of completeness, we list here the normal stress components of the constitutive equation for third grade fluid (4) by making use of the assumptions () and (4). These are du() r rr = ( α+ α) S dr, (6) du() r Szz = α dr. (64) The expression for normal stress difference turns out to be du() r Srr Szz = α dr, 4 α dp β dp β dp 4 Srr Szz = r r r + μ μ μ μ μ β dp 6 β dp 8 r + 55 r + L. (65) μ μ μ μ By setting the parameter α equal to zero, the difference becomes S S =, (66) which agrees with that of the Newtonian case. rr zz 4(b) Solution for 6-Constant Oldroyd Fluid The governing equation (8) for pipe of 6-constant Oldroyd fluid may be written in operator form as where αr α r dp du dp du Lu = + μ dr μ dr, (67)

12 49 A. M. Siddiqui et al d L =. dr The corresponding boundary condition associated with pipe flow is the no slip condition u = and r = R. (68) As before, we apply the inverse operator L - on both sides, to get r dp du α dp du L L u() r = L αl + L r, μ dr μ dr which may be written as r dp du α dp du ur () = C+ L αl + L r, (69) μ dr μ dr where C is the constant of integration. Substituting the decomposition series ur () = u() r, n n= () = n + n n= n=, (7) Nu r A B into both sides of (69), we obtain r dp α dp un() r = C+ L αl An + L r Bn n= μ n= μ n=, (7) where du du An =, Bn = dr dr. (7) We have already generated the Adomian polynomials for A n in our previous problem. The Adomian polynomials for B n, obtained with the help of eqn. (7), are du = B dr, du du B =, dr dr du du du B = +, (7) dr dr dr du du du du B = +, dr dr dr dr

13 Adomian decomposition method 49 du du4 du du du 4 = + +, B dr dr dr dr dr M The rest of the polynomials, if desired, can be constructed in a similar manner. We identify the zero th component by r dp u () r = C+ L μ, (74) and the remaining components u n+ (r) by the recurrence relation α dp un+ () r = αl An + L r Bn, n= μ n= n (75) From (75), the first few u n+ (r) are α [ ] dp u() r = αl A + L r B, μ α [ ] dp u() r = αl A + L r B, μ α [ ] dp u() r = αl A + L r B, μ α [ ] dp u4() r = αl A + L r B, μ (76) M Substituting the decomposition series (7), into the boundary condition (44), we get u () n r = at r = R, n. (77) Using the boundary condition (77), we obtain from (74) and (75) dp r R u() r = μ, (78) 4 4 dp r R u() r = ( α α), μ dp r R u() r = ( α α)( α α), μ dp r R u() r = ( α α)( α α)( 6α 5α) μ 8 8, 5 7

14 49 A. M. Siddiqui et al ( α α ) ( α α )( α α )( α α ) ( α α ) u () r = dp r R ( α α)( 6α α) + α( α α). μ Finally, the fourth order approximate solution for the velocity profile is given by 4 4 dp r R ( α α) dp r R u = μ μ dp r R + ( α α)( α α) μ dp r R ( α α)( α α)( 6α 5α) μ ( α α ) ( α α )( 6α 5α )( α α ) ( α α ) + ( α α )( α α ) α ( α α ) 7 9 dp r R (79) + 6. μ In view of (79), expressions for practically important quantities, namely the volume flow rate Q, the average velocity urand () normal stress difference Srr Szz can be computed as in the preceding problem Numerical Method In this section, we present numerical methods used to find the solution of the problem of the fully developed pipe flow for the third grade non-newtonian fluid as well as for 6-constant Oldroyd fluid. We will solve non-dimensionalize governing equations for these problems. Our interest is to compare the analytical solutions obtained by Adomian decomposition method with numerical solutions. We find numerical solutions using Chebyshev spectral collocation method. Chebyshev polynomials have been successfully employed in finding numerical solution of various boundary value problems and in the study of computational fluid dynamics [6]. In solving ordinary differential equations, linear or nonlinear, to high accuracy on a simple domain, and if the data defining the problem are smooth, then spectral methods are usually the best tool. They can often achieve high order of accuracy (ten digits) of accuracy as compared with the finite difference or finite element methods. We consider the flow of a non-newtonian third grade fluid in a circular pipe. The governing equation was given in (5) as

15 Adomian decomposition method 49 du β du r dp + =. (8) dr μ dr μ When finding the analytical solution, we solved the above first order equation subject to the no slip boundary condition at the walls of the pipe. To facilitate the numerical solution, we differentiate the above equation with respect to r, keeping in mind that the pressure does not vary with r, and obtain a second order differential equation and solve it subject to the no slip at the walls and symmetry condition at the centerline. The system we solve numerically is given below: d u 6β du d u dp + = (8) dr μ dr dr μ with boundary conditions given as u'( r) = at r =, (8) ur () = at r= R, where the prime denotes the derivative with respect to r. First, we non-dimensionalize the equation and the boundary conditions by introducing the following non-dimensional parameters: u* u r z U Rp, r*, z*, *, p* U R β λ R μ R μ U. (8) The non-dimensional equation (after dropping *) together with the boundary conditions yield d + 6λ =, dr dr dr (84) u'( r) = at r =, ur () = at r=. (85) To solve equation (84) subject to boundary conditions (85), we employ spectral collocation method of Chebyshev type. We look for an approximate solutionφ N which is a global Chebyshev polynomial of degree N defined on the interval [-, ] by TN ( y) = cos Nθ, θ = cos y. (86) We discretize the interval by using collocation points to define the Chebyshev nodes in [-, ], namely jπ yj = cos, N j =,,... N. The function φ ( y) is approximated by an interpolating polynomial which is constructed in terms of the values of φ at each of the collocation points by employing a truncated Chebyshev series of the form

16 494 A. M. Siddiqui et al N φ( y) = φ ( y ) = % φ T ( y ), j =,,..., N, (87) N j k k j k = where % φ k represents the series coefficients. The derivatives of the functions at the collocation points are given by N φn ( yj) = Djkφk, dy % (88) k = where D represents the derivative matrix [7] and is given by j + k c ( ) j D = j k; k =,,..., N. (89) jk c y y k j k The higher order derivatives are computed as simply multiple powers of D, i.e., i N φn i ( yj) = jk k, i D φ dy k = k =,,..., N, (9) where i is the order of the derivative. As described above, the Chebyshev polynomials are defined on the finite interval [-,]. Therefore to apply Chebyshev spectral method to our equation, we make a suitable linear transformation and transform the physical domain [,] to a computational domain [-, ]. We sample the unknown function u at the Chebyshev T points to obtain the data vector u = [ u( y), u( y),... u( y N )]. The next step is to find a Chebyshev polynomial φ of degree N that interpolates the data, i.e., φ ( yj) = uj, j =,,... N, and obtain the spectral derivative vector u '' by differentiating φ and evaluating at the grid points, i.e., u'' j = φ ''( tj), for j =,,..., N. This transforms the nonlinear differential equation into a nonlinear algebraic equations which can be solved by Newton s iterative method. Using the numerical method described in this section, we first consider flow of a third grade fluid in a pipe, which is modeled by the non-dimensional equation (84) with boundary conditions given in (85). A comparison of numerical solutions and the ADM solution is presented in the graphs to show the validity of the analytical solutions. The results are presented for different values of non-newtonian parameter and imposed pressure gradient. The numerical results are found be in good agreement with the analytical results obtained by the Adomian s decomposition method.

17 Adomian decomposition method 495 Fig. Radial velocity field u for third grade pipe flow, (a) Comparison numerical and ADM solution for λ =., dp (b) Comparison numerical and ADM solution for λ =. dp =-.5. dx =-.5, dx Fig. (a) Effect of non-newtonian parameter λ on radial velocity for third grade pipe flow for fixed pressure gradient. (b) Effect of imposed pressure gradient dp/dx on radial velocity for third grade pipe flow for fixed λ. We now consider the flow of a non-newtonian 6-constant Oldroyd fluid in a circular pipe of radius R. Model equation for this problem which is derived in earlier section and is given below du + α du α r dp du r dp =. (9) dr dr μ dr μ When finding the analytical solution, we solved the above first order highly nonlinear equation subject to no-slip boundary condition at the walls of the pipe. As in the third grade case, we differentiate the above equation to facilitate the numerical solution, and get a second order differential equation. We solve the resulting equation subject to the no slip condition at the walls and symmetry condition at the centerline. The equation we solve numerically is give as

18 496 A. M. Siddiqui et al α α α r d u du d u dp du dp du d u dp + dr dr dr μ dr μ =, (9) dr dr μ with boundary conditions given by u'( r) = at r =, (9) ur () = at r= R. We non-dimensionalize the equation and the boundary conditions by introducing the following non-dimensional parameters: u r z U α U α Rp u* =, r* =, z* =, λ* =, λ * =,p* = U R R μ R μ R μu Then, the non-dimensional equation (after dropping *) with boundary conditions becomes λ λ λ r d u du d u dp du dp du d u dp + dr dr dr =, (94) dr dr dr u'( r) = at r =, ur () = at r=. (95) Fig. Radial velocity field u for 6 constant Oldroyd pipe flow (a) Comparison numerical and ADM solution for λ =., λ =. and fixed pressure gradient dp =-, (b) Comparison numerical and ADM solution for λ =., λ =.4 and same fixed pressure gradient. dx

19 Adomian decomposition method 497 Fig. 4 Radial velocity field u for 6 constant Oldroyd fluid in a pipe (a) Effect of material constant λ, with keeping λ =., and fixed pressure gradient (b) Effect of material constant λ, with keeping λ =., and fixed pressure gradient. 6 Discussion & Conclusion Numerical and analytical solutions are obtained for the incompressible flow of a non-newtonian third order fluid and 6-constant Oldroyd fluid through an axisymmetric tube. Analytical solutions are obtained by the Adomian decomposition method and the numerical solutions are obtained using highly accurate Chebyshev spectral methods. To the best of our knowledge, this is the first attempt to apply this robust and highly effective analytical technique as well as highly accurate Chebyshev spectral method to study the non-newtonian flows in this geometry. The proposed numerical method offers a superior intrinsic accuracy for the derivative calculations. The numerical results indicate the usefulness of the spectral methods in obtaining accurate solutions to the nonlinear problems arising in non-newtonian fluid mechanics. As compared to other numerical techniques, such as finite differences, the nonlinearity is not a major complication for spectral methods. It is found that the analytical method used in this work (ADM) is more powerful and easy to apply for obtaining analytic approximate solutions. It provides analytic, verifiable, and rapidly convergent approximations which allow further insight into the characteristics and behavior of the solution as is the case with closed form solutions. Moreover, this method introduces the solutions, without any need for linearization or discretization which produces a rapidly converging series solution, components of which are elegantly computed. We would like to point out that even first few terms of the solution series give a pretty accurate approximation which compares well with the numerical results and proves the validity of this method.

20 498 A. M. Siddiqui et al The first set of results represents the solutions for the fully developed Poiseuille flow in a pipe. Fig. a,b shows the comparison of numerical solutions with the analytical solutions for the velocity field in the case of a third order fluid. The solutions are compared for different values of non-newtonian parameter and imposed pressure gradient and we notice a close agreement between the two solutions. Fig. (a) shows the influence of third grade parameter on the velocity profile. It is observed that for a given pressure gradient, the increase in non- Newtonian parameter λ results in the increase of velocity. Fig. (b) depicts the velocity profiles for the fixed non-newtonian parameter λ =. and for different values of the imposed pressure gradients. Fig., 4 shows the results for a different non-newtonian fluid model known as the Oldroyd 6-constant fluid. In our formulation of the pipe flow only two of the material constants λ and λ appear. We again observe a close agreement between the two solutions. In fig. 4, we show the effects of these parameters on the flow profile. We find that for a fixed valueλ =., the increase in the λ results in the increase of the radial velocity and vice versa. On the other hand for a fixed λ =., the increase in the value of λ results in the decrease of the radial velocity. Analytical expressions for the other physical quantities of interest such as volume flux Q, average velocity u, shear stress distribution and normal stress are also computed and presented for both third grade fluid as well as 6-constant Oldroyd fluid. The velocity filed and other physical quantities for the classical viscous fluid can be recovered from the given expressions as a special case in appropriate limits. References [] G. Adomian, A review of the decomposition method and some recent results for nonlinear equation, Math. Comput. Model. (7) (99) 7-4. [] G. Adomian, Solving Frontier Problems of Physics: The Decomposition Method, Kluwer, Dordrecht, (994). [] E. Babolian, A. R. Vahidi, G. H. Asadi Cordshooli, Solving differential equations by decomposition method, Appl. Math. Comput. 67 () (5) [4] H. Bayramoglu, N. Bildik, The solution of two dimensional nonlinear differential equation by the Adomian decomposition method, Appl. Math. Comput. 6 () (5) [5] G. Adomian, Y. Cherruault, Decomposition method: a new proof of convergence, Math. Comput. Model. 8() (99) -6.

21 Adomian decomposition method 499 [6] M. Dehghan, Application of the Adomian decomposition method for twodimensional parabolic equation subject to nonstandard boundary specifications, Appl. Math. Comput., 57(4) [7] S. Pamuk, Solution of the porous media equation by Adomian s decomposition method, Phys. Lett. A 4 (5) [8] M. Ahmed, Q. K. Ghori, S. Islam, A. M. Siddiqui, Homotopy analysis of Couette and Poisuille flows for a fourth grade fluid, Acta Mechanica, 8(5), no. -4. [9] A. M. Wazwaz, A reliable technique for solving the weakly singular secondkind Volterra-type integral equations, Appl. Math. Comput. 8 (996) [] Xing-Guo Luo, Qing-Biao Wu, Bing-Quan Zhang, Revisit on partial solutions in the Adomian decomposition method: Solving heat and wave equations, J. Math. Anal. Appl. () (6) 5-6. [] S. Nadeem, N. Akbar, Effects of heat transfer on the peristaltic transport of MHD Newtonian fluid with variable viscosity: Application of Adomian decomposition method, Commun Nonlinear Sci Numer Simulat. 4 () (9) [] J. E. Dunn, K. R. Rajagopal, Fluids of differential type: Critical review and thermodynamic analysis, Internat. J. Eng. Sci. (5) (995) [] C. Truesdell, W. Noll, S. Antman, The non-linear fields theories of Mechanics ( Ed.), Springer, 4. [4] K. R. Rajagopal, On the stability of third-grade fluids, Arch. Mech. (6) (98) 867. [5] R. L. Fosdick, K. R. Rajagopal, Thermodynamics and stability of fluids of third grade, Proc. R. Soc. Lond. A 69 (98) [6] M. Hussaini at. el., Spectral Methods in Fluid Dynamics Springer (988). Received: March,