JAM Journal of Applied Mathematics 1687-0042 1110-757X Hindawi Publishing Corporation 259608 10.1155/2012/259608 259608 Research Article Boundary Layer Equations and Lie Group Analysis of a Sisko Fluid http://orcid.org/0000-0002-0046-9090 Sarı Gözde 1, 2 Pakdemirli Mehmet 1, 2 Hayat Tasawar 3, 4 Aksoy Yiğit 1 Mahomed Fazal M. 1 Department of Mechanical Engineering Celal Bayar University Muradiye 45140 Manisa Turkey bayar.edu.tr 2 Applied Mathematics and Computation Center Celal Bayar University Muradiye 45140 Manisa Turkey bayar.edu.tr 3 Department of Mathematics Quaid-i-Azam University 45320 Islamabad 44000 Pakistan qau.edu.pk 4 Department of Mathematics King Saud University P.O. Box 2455 Riyadh 11451 Saudi Arabia ksu.edu.sa 2012 20 12 2012 2012 27 09 2012 05 12 2012 2012 Copyright © 2012 Gözde Sarı et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Boundary layer equations are derived for the Sisko fluid. Using Lie group theory, a symmetry analysis of the equations is performed. A partial differential system is transferred to an ordinary differential system via symmetries. Resulting equations are numerically solved. Effects of non-Newtonian parameters on the solutions are discussed.

1. Introduction

Due to the inadequacy of Newtonian fluid model which predicts a linear relationship between the shear stress and velocity gradient, many non-Newtonian fluid models were proposed to explain the complex behavior. Usually, the stress constitutive relations of such models inherit complexities which lead to highly nonlinear equations of motion with many terms. To simplify the extremely complex equations, one alternative is to use boundary layer theory which is known to effectively reduce the complexity of Navier-Stokes equations and reduce drastically the computational time. Since there are many non-Newtonian models and new models are being proposed continuously, boundary layer theory for each proposed model also appears in the literature. It is beyond the scope of this work to review vast literature on the boundary layers of non-Newtonian fluids. A limited work on the topic can be found in .

In this work, boundary layer theory is developed for the Sisko fluid, a non-Newtonian fluid model which combines the features of viscous and power law models. A complete symmetry analysis of the boundary layer equations is presented. Using one of the symmetries, the partial differential system is transformed into an ordinary differential system. The resulting ordinary differential system is numerically solved by a finite difference algorithm. Effect of non-Newtonian parameters on the velocity profiles is shown in the graphs.

As a non-Newtonian model, Sisko  proposed a new model which is later called sisko fluid. Some of the recent work on Sisko fluids is as follows: an analytical solution was presented using homotopy analysis method for the flow of a magnetohydrodynamic Sisko fluid through a porous medium . Wire coating analysis was investigated by withdrawal from a bath of Sisko fluid . The equation modeling thin film flow of a Sisko fluid on a moving belt qualitatively was analyzed and a series solution was found using homotopy analysis method . Unsteady unidirectional flow of a Sisko fluid bounded by a suddenly moved plate was studied . Unsteady flow of a Sisko fluid over a moving wall with suction or blowing was considered . A numerical study for the unsteady flow of a magnetohydrodynamic Sisko fluid in annular pipe was presented . Heat transfer was also incorporated into the annular pipe problem . To the best of authors’ knowledge, boundary layer analysis and symmetry reductions of Sisko fluids do not exist in the literature.

2. Boundary Layer Equations

The Cauchy stress tensor for Sisko fluid is (2.1)T=-pI+S, where (2.2)S=[a+b|12tr(A12)|n-1]A1,A1=L+LT,L=grad  V.V is the velocity vector and A1 is the first Rivlin Ericksen tensor. a and b are the material constants. The model is a combination of viscous and power-law models. For a=0, the model exhibits power-law behavior whereas, for b=0, the flow is Newtonian. Steady-state, two-dimensional, incompressible equations of motion including mass conservation can be written as (2.3)u*x*+v*y*=0ρ(u*u*x*+v*u*y*)=-p*x*+Sxxx*+Sxyy*ρ(u*v*x*+v*v*y*)=-p*y*+Sxyx*+Syyy*, where x* is the spatial coordinate along the surface, y* is vertical to it, u* and v* are the velocity components in the x* and y* coordinates. The shear stress components are inserted into the equations of motion and the usual boundary layer assumptions are made, that is, x*~O(1), y*~O(δ), u*~O(1), v*~O(δ). The highest order terms are retained and the momentum equations become (2.4)ρ(u*u*x*+v*u*y*)=-p*x*+a2u*y*2+bn(u*y*)n-12u*y*20=p*y* from which dependence of pressure on y* is eliminated. Dimensionless variables and parameters are defined as follows: (2.5)x=x*L,y=y*L,u=u*V,v=v*V,p=p*ρV2ε1=aρVL,ε2=bnVn-2ρLn, where L is a characteristic length and V is a reference velocity. Expressing the pressure in terms of outer velocity, the final dimensionless boundary layer equations become (2.6)ux+vy=0uux+vuy=UdUdx+ε12uy2+ε2(uy)n-12uy2.

The classical boundary conditions for the problem are (2.7)u(x,0)=0,v(x,0)=0,u(x,)=U(x). For ε2=0 or n=1, the equations reduce qualitatively to those of Newtonian fluid.

3. Lie Group Theory and Symmetry Reductions

Lie group Theory is employed in search of symmetries of the equations. Details of the theory can be found in [28, 29]. The infinitesimal generator for the problem is (3.1)X=ξ1(x,y,u,v)x+ξ2(x,y,u,v)y+η1(x,y,u,v)u+η2(x,y,u,v)v. A straightforward and tedious calculation (see [28, 29] for details) yields (3.2)ξ1=3ax+b,ξ2=ay+c(x),η1=au,η2=c(x)u-av. The classifying relation for the outer velocity is (3.3)(3ax+b)ddx(UU)+a(UU)=0. Symbolic packages developed to calculate symmetries fail to produce the above results due to the arbitrary outer velocity function, and hand calculation becomes inevitable for precise results. There are two finite parameter Lie group symmetries represented by parameters “a” and “b,” the former corresponding to scaling symmetry and the latter to translational symmetry in the x coordinate. There is an additional infinite parameter Lie group symmetry represented by c(x) function. For a general stress tensor which is an arbitrary function of the velocity gradient, the symmetries were calculated previously . Our results for the special case of Sisko fluid confirm the previous calculations (see principal Lie Algebra presented in ).

Usually boundary conditions put much restriction on the symmetries which may lead to removal of all the symmetries. In our case, however, some of the symmetries remain stable after imposing the boundary conditions. For nonlinear equations, the generators should be applied to the boundaries and boundary conditions also . Applying the generator to the boundary y=0 yields c(x)=0. Applying the generator to the boundary conditions does not impose further restrictions and hence the symmetries valid for the equations and boundary conditions reduce to (3.4)ξ1=3ax+b,ξ2=ay,η1=au,η2=-av with the classifying relation (3.3) remaining unchanged.

Selecting parameter “a” in the symmetries, the associated equations which define similarity variables are (3.5)dx3x=dyy=duu=dv-v. Solving the system yields the similarity variables (3.6)ξ=yx1/3,u=x1/3f(ξ),v=g(ξ)x1/3. From (3.3) with b=0, U(x)=x1/3. Substituting all into the boundary layer equations yields the ordinary differential system (3.7)f-ξf+3g=0f2-ξff+3gf=1+3ε1f′′+3ε2  (f)n-1f′′. The boundary conditions also transform to (3.8)f(0)=0,g(0)=0,f()=1.

4. Numerical Results

Equation (3.7) with the boundary conditions (3.8) are numerically integrated using a finite difference code. The finite difference code implements the three-stage Lobatto IIIA formula that provides continuous solutions of fourth order accuracy in the problem domain. Description of the numerical method is given by Shampine et al.  and Kierzenka and Shampine . In Figure 1, f function and in Figure 2, g function related to the x and y components of the velocities are drawn for different ε1 parameters. Boundary layer becomes thicker for an increase in ε1.

Effect of ε1 parameter on the similarity function f related to the x component of velocity (ε2=1, n=2).

Effect of ε1 parameter on the similarity function g related to the y component of velocity (ε2=1, n=2).

There is an increase in y component of velocity which is negative over the whole domain excluding the surface for an increase in ε1 as can be seen from Figure 2. A similar trend is observed for parameter ε2, although not as influential as in the previous case (see Figures 3 and 4).

Effect of ε2 parameter on the similarity function f related to the x component of velocity (ε1=1, n=2).

Effect of ε2 parameter on the similarity function g related to the y component of velocity (ε1=1, n=2).

On the contrary, a reverse effect is observed for power index n, that is, as n increases, boundary layer becomes narrower eventually converging to a fixed thickness (see Figure 5). For the term related to y component of velocity, higher values of power index n cause a decrease which eventually converges to a solution as can be seen from Figure 6.

Effect of power index n on the similarity function f related to the x component of velocity (ε1=1, ε2=1).

Effect of power index n on the similarity function g related to the y component of velocity (ε1=1, ε2=1).

5. Concluding Remarks

Boundary layer equations of Sisko fluid derived for the first time. Lie group theory is applied to the equations. Equations admit two finite parameter Lie group transformations and an infinite parameter Lie group transformation. The infinite parameter Lie group transformation is not stable with respect to usual boundary layer conditions. Using the scaling symmetry which is one of the finite parameter transformations, the partial differential system is transferred into an ordinary differential system. Resulting ordinary differential system is solved numerically using a finite difference scheme. Effects of non-Newtonian parameters on the boundary layers are discussed in detail. It is found that an increase in the parameters ε1 and ε2 causes the boundary layer to thicken whereas a reverse effect is observed for power index n.

Acknowledgments

This work is completed during mutual short visits of T. Hayat to Turkey and M. Pakdemirli to Pakistan. Funding supports of TUBITAK of Turkey and HEC of Pakistan are highly appreciated.

Acrivos A. Shahand M. J. Petersen E. E. Momentum and heat transfer in laminar boundary layer flows of non-Newtonian fluids past external surfaces AIChE Journal 1960 6 312 317 Schowalter W. R. The application of boundary-layer theory to power-law pseudoplastic fluids: similarity solutions AIChE Journal 1960 6 25 28 Bizzell G. D. Slattery J. C. Non-newtonian boundary-layer flow Chemical Engineering Science 1962 17 10 777 782 2-s2.0-0010414563 Srivastava A. C. The flow of a non-Newtonian liquid near a stagnation point Zeitschrift für Angewandte Mathematik und Physik 1958 9 1 80 84 2-s2.0-34250611717 10.1007/BF01596862 Beard D. W. Walters K. Elastico-viscous boundary-layer flows Proceedings of Cambridge Philosophical Society 1964 60 667 674 MR0171475 ZBL0123.41601 Astin J. Jones R. S. Lockyer P. Boundary layers in non-Newtonian fluids Journal Mechanique 1973 12 3 527 539 2-s2.0-0015661901 Rajagopal K. R. Gupta A. S. Wineman A. S. On a boundary layer theory for non-Newtonian fluids International Journal of Engineering Science 1980 18 6 875 883 2-s2.0-0019221447 Rajagopal K. R. Gupta A. S. Na T. Y. A note on the Falkner-Skan flows of a non-Newtonian fluid International Journal of Non-Linear Mechanics 1983 18 4 313 320 10.1016/0020-7462(83)90028-8 MR718753 ZBL0527.76010 Garg V. K. Rajagopal K. R. Flow of a non-Newtonian fluid past a wedge Acta Mechanica 1991 88 1-2 113 123 10.1007/BF01170596 MR1111096 Massoudi M. Ramezan M. Effect of injection or suction on the Falkner-Skan flows of second grade fluids International Journal of Non-Linear Mechanics 1989 24 3 221 227 2-s2.0-0024912212 Garg V. K. Rajagopal K. R. Stagnation point flow of a non-Newtonian fluid Mechanics Research Communications 1990 17 6 415 421 10.1016/0093-6413(90)90059-L MR1078763 ZBL0727.76003 Pakdemirli M. Şuhubi E. S. Boundary layer theory for second order fluids International Journal of Engineering Science 1992 30 4 523 532 10.1016/0020-7225(92)90042-F MR1150129 ZBL0747.76013 Pakdemirli M. The boundary layer equations of third-grade fluids International Journal of Non-Linear Mechanics 1992 27 5 785 793 2-s2.0-0026923728 Pakdemirli M. Boundary layer flow of power-law fluids past arbitrary profiles IMA Journal of Applied Mathematics 1993 50 2 133 148 10.1093/imamat/50.2.133 MR1220765 ZBL0773.76004 Bhatnagar R. K. Gupta G. Rajagopal K. R. Flow of an Oldroyd-B fluid due to a stretching sheet in the presence of a free stream velocity International Journal of Non-Linear Mechanics 1995 30 3 391 405 2-s2.0-58149365641 Hagen T. Renardy M. Boundary layer analysis of the Phan-Thien-Tanner and Giesekus model in high Weissenberg number flow Journal of Non-Newtonian Fluid Mechanics 1997 73 1-2 181 189 2-s2.0-0031280060 Yürüsoy M. Pakdemirli M. Symmetry reductions of unsteady three-dimensional boundary layers of some non-Newtonian fluids International Journal of Engineering Science 1997 35 8 731 740 10.1016/S0020-7225(96)00115-2 MR1464855 ZBL0904.76006 Yürüsoy M. Pakdemirli M. Group classification of a non-Newtonian fluid model using classical approach and equivalence transformations International Journal of Non-Linear Mechanics 1999 34 2 341 346 10.1016/S0020-7462(98)00037-7 MR1658926 Aksoy Y. Pakdemirli M. Khalique C. M. Boundary layer equations and stretching sheet solutions for the modified second grade fluid International Journal of Engineering Science 2007 45 10 829 841 10.1016/j.ijengsci.2007.05.006 MR2354835 ZBL1213.76060 Sisko A. W. The flow of lubricating greases Industrial and Engineering Chemistry 1958 50 12 1789 1792 Khan M. Abbas Z. Hayat T. Analytic solution for flow of Sisko fluid through a porous medium Transport in Porous Media 2008 71 1 23 37 10.1007/s11242-007-9109-4 MR2377316 Sajid M. Hayat T. Wire coating analysis by withdrawal from a bath of Sisko fluid Applied Mathematics and Computation 2008 199 1 13 22 10.1016/j.amc.2007.09.032 MR2415798 ZBL1137.76009 Akyıldız F. T. Vajravelu K. Mohapatra R. N. Sweet E. Van Gorder R. A. Implicit differential equation arising in the steady flow of a Sisko fluid Applied Mathematics and Computation 2009 210 1 189 196 10.1016/j.amc.2008.12.073 MR2504135 ZBL1160.76002 Molati M. Hayat T. Mahomed F. Rayleigh problem for a MHD Sisko fluid Nonlinear Analysis. Real World Applications 2009 10 6 3428 3434 10.1016/j.nonrwa.2008.09.023 MR2561355 Hayat T. Moitsheki R. J. Abelman S. Stokes' first problem for Sisko fluid over a porous wall Applied Mathematics and Computation 2010 217 2 622 628 2-s2.0-77955659964 10.1016/j.amc.2010.05.099 Khan M. Abbas Q. Duru K. Magnetohydrodynamic flow of a Sisko fluid in annular pipe: a numerical study International Journal for Numerical Methods in Fluids 2010 62 10 1169 1180 10.1002/fld.2068 MR2640689 Khan M. Munawar S. Abbasbandy S. Steady flow and heat transfer of a Sisko fluid in annular pipe International Journal of Heat and Mass Transfer 2010 53 7-8 1290 1297 2-s2.0-74549135885 10.1016/j.ijheatmasstransfer.2009.12.037 Bluman G. W. Kumei S. Symmetries and Differential Equations 1989 81 New York, NY, USA Springer xiv+412 MR1006433 Stephani H. Differential Equations: Their Solution Using Symmetries 1989 New York, NY, USA Cambridge University Press xii+260 MR1041800 Shampine L. F. Gladwell I. Thompson S. Solving ODEs with MATLAB 2003 chapter 3 Cambridge, UK Cambridge University Press viii+263 10.1017/CBO9780511615542 MR1985643 Kierzenka J. Shampine L. F. A BVP solver based on residual control and the MATLAB PSE Transactions on Mathematical Software 2001 27 3 299 316 10.1145/502800.502801 MR1918120