Abstract

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 [119].

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 [20] 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 [21]. Wire coating analysis was investigated by withdrawal from a bath of Sisko fluid [22]. 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 [23]. Unsteady unidirectional flow of a Sisko fluid bounded by a suddenly moved plate was studied [24]. Unsteady flow of a Sisko fluid over a moving wall with suction or blowing was considered [25]. A numerical study for the unsteady flow of a magnetohydrodynamic Sisko fluid in annular pipe was presented [26]. Heat transfer was also incorporated into the annular pipe problem [27]. 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 where is the velocity vector and is the first Rivlin Ericksen tensor. and are the material constants. The model is a combination of viscous and power-law models. For , the model exhibits power-law behavior whereas, for , the flow is Newtonian. Steady-state, two-dimensional, incompressible equations of motion including mass conservation can be written as where is the spatial coordinate along the surface, is vertical to it, and are the velocity components in the and coordinates. The shear stress components are inserted into the equations of motion and the usual boundary layer assumptions are made, that is, , , , . The highest order terms are retained and the momentum equations become from which dependence of pressure on is eliminated. Dimensionless variables and parameters are defined as follows: where is a characteristic length and is a reference velocity. Expressing the pressure in terms of outer velocity, the final dimensionless boundary layer equations become

The classical boundary conditions for the problem are For or , 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 A straightforward and tedious calculation (see [28, 29] for details) yields The classifying relation for the outer velocity is 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 “” and “,” the former corresponding to scaling symmetry and the latter to translational symmetry in the coordinate. There is an additional infinite parameter Lie group symmetry represented by function. For a general stress tensor which is an arbitrary function of the velocity gradient, the symmetries were calculated previously [18]. Our results for the special case of Sisko fluid confirm the previous calculations (see principal Lie Algebra presented in [18]).

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 [28]. Applying the generator to the boundary yields . 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 with the classifying relation (3.3) remaining unchanged.

Selecting parameter “” in the symmetries, the associated equations which define similarity variables are Solving the system yields the similarity variables From (3.3) with , . Substituting all into the boundary layer equations yields the ordinary differential system The boundary conditions also transform to

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. [30] and Kierzenka and Shampine [31]. In Figure 1, function and in Figure 2, function related to the and components of the velocities are drawn for different parameters. Boundary layer becomes thicker for an increase in .

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

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

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 and causes the boundary layer to thicken whereas a reverse effect is observed for power index .

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.