 Open Access
 Total Downloads : 16
 Authors : Ateka Pathan , M. G. Timol
 Paper ID : IJERTCONV3IS23022
 Volume & Issue : NCETRASECT – 2015 (Volume 3 – Issue 23)
 Published (First Online): 24042018
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Similarity Solution of MHD Boundary Layer Flow of PrandtlEyring Fluids
Special Issue – 2015
International Journal of Engineering Research & Technology (IJERT)
ISSN: 22780181
NCETRASECT2015 Conference Proceedings
Ateka Pathan
Department of Mathematics,
Shree J. P. Arts and Science College, Bharuch, Gujarat, INDIA
M. G. Timol
Department of Mathematics,
Veer Narmad South Gujarat University, Surat7, Gujarat, INDIA
Abstract This paper investigates the magneto hydrodynamics (MHD) boundary layer flow of nonNewtonian PrandtlEyring fluids past a right angle wedge. The use of Lie group transformation method socalled scaling group symmetry analysis is made to derive possible similarity transformations of the present flow problem. The important conclusion drawn from the present analysis is that for all those nonNewtonian fluids whose shearing stress is composite function of rate of strain, the similarity solutions exist only for the flows past 90 degree wedge. The solution of similarity equation, which is nonlinear ordinary differential, is obtained by KellerBox method for the various values of flow parameters. Effect of imposed magnetic field on the flow of nonNewtonian PrandtlEyring fluids is also studied.
KeywordsNonNewtonian Fluids; PrandtlEyring Fluids; MHD Boundary layer flow; similarity solution; Scaling group symmetry; NonNewtonian Fluids; PrandtlEyring Fluid;

INTRODUCTION
The partial differential equations governing the motion of fluid flow problems are usually nonlinear in nature and hence cannot be solved easily. Whenever possible these differential equations are reduced to ordinary differential equations by employing transformations to obtain similar solutions. Ames [1], BlumanKumai[2],Hansen[3],Seshadri Na[4] and Stephani[5] have discussed application of groups and symmetries to partial differential equations arising from natural phenomena and technological problems. Symmetry groups are invariant transformations which do not alter the structural form of the equation under investigation. The advantage of the symmetry method is that it can be applied successfully to nonlinear partial differential equations governing the motion of fluid. Sophus Lie developed a transformation, currently known as Lie group of transformation, which maps a given differential equation to itself. The differential equations remain invariant under some continuous group of transformations usually known as symmetries of a differential equation.
The classical theory of Newtonian fluid depends upon the hypothesis of linear relationship between stress tensor and rate of strain tensor. The fluids which do not follow such a linear relationship are called nonNewtonian fluids. Non Newtonian fluids are generally divided in to two categories like isotropic and homogeneous and when they are subjected to a shear the resultant stress depends only on the rate of shear. However, such types of fluids show diverse behavior in response to applied stress. Numbers of rheological models
have been proposed to explain such a diverse behavior. Some of this models are ; Powerlaw fluids, Sisko fluids, Ellis fluids, Prandtl fluids Williamson fluids, Sutterby fluids ReinerRivlin fluids, Bingham plastic , PrandtlEyring fluids , PowellEyring fluids , ReinerPhilippoff etc.
When we consider electrically conducting nonNewtonian fluids flowing under the influence of external magnetic field, the study becomes interesting .This is because in such situation magnetic forces produced in it could influence the motion of the fluids in significant way and hence such interaction problems have great practical applications. The problem of twodimensional magneto hydrodynamic boundary layer equation for laminar incompressible flow past flat plate has been investigated by Rossow [6] and Greenspan et al [7].Rossow[6] has considered transverse magnetic field where as Greenspan et al [7] have considered longitudinal magnetic fields on the velocity and temperature distributions. Timol et al [8] have investigated threedimensional magneto hydrodynamic boundary layer flow with pressure gradient and fluid injection. Similarity transformation for both steady and unsteady threedimensional MHD boundary layer flow of purely viscous nonNewtonian fluid has been derived by Manisha et al[9].They have also derived Similarity Analysis in MHD Heat and Mass Transfer of NonNewtonian Power Law Fluids Past a Semiinfinite Flat Plate [10].
To investigate the nonNewtonian effects, the class of solutions known as similarity solutions place an important role. This is because that is the only class of the exact solution for the governing equations which are usually non linear partial differential equations (PDEs) of the boundary layer type. Further this also serves as a reference to check approximate solutions.
It is well known that similarity solutions for the PDEs governing the flow of Newtonian and nonNewtonian fluids exist only for limited classes of main stream velocities at the edge of the boundary layer. For example, for two dimensional laminar boundary layer flow of Newtonian fluids, similarity solutions are limited to the well known FalknerSkan solution Rajagopal et al [11]. Most of the generalization of the FalknerSkan solutions and approximate solutions in the literature are limited to the power law fluids; this is because they are mathematically the easiest to be treated among most of the nonNewtonian fluids.
Special Issue – 2015
International Journal of Engineering Research & Technology (IJERT)
ISSN: 22780181
NCETRASECT2015 Conference Proceedings
Hansen and Na [12] were probably first to derive similarity analysis of laminar incompressible boundary layer equations of the nonNewtonian fluids whose stress and rate of strain are related by arbitrary continuous function, while extending the work of Lee and Ames [13].They have used linear group of transformations to derive their similarity equations which includes many of the nonNewtonian viscoinelastic fluids. Timol and Kalthia [14] are probably first to derive similarity solution for the class of three dimensional boundary layer flow of viscoinelastic non Newtonian fluids.
For the derivation of the constitutive equations governing the motion of nonNewtonian fluids, the mathematical structure of stressstrain relationship, which is non linear, is important in when found in functional form. This relationship may be implicit or explicit. In the present paper we have consider such relationship in the form of general arbitrary continuous function of the type:
Here is the shearing stress and u/y is the rate of the strain of the fluids.
For the similarity analysis many techniques are available, among them the similarity methods which invoke the invariance under the group of transformations are known as group theoretic methods. These methods are more recent and are mathematically elegant and hence they are widely used in different fields. The group theoretic methods involve mainly two different types of groups of transformations, namely, assumed group of transformations, spiral group transformations are the assumed group of transformations and are mainly due to Birkhoff [15] and Morgan [16].
So Motivated by these , we represent in the present paper Using Scaling group transformation technique ,class of similarity solution for steady ,two dimensional laminar MHD boundary layer flows of incompressible nonNewtonian is derived .From the present analysis it is interesting to observe that for nonNewtonian viscoinelastic fluids of ay model , which is characterized by the property that its stress and the rate of strain can be related by arbitrary continuous function given by equation (1), the similarity solutions exist only for the flows past 900 wedge , as shown in Fig.1.
Figure 1: Schematic Diagram of flow past 90Â° wedge

GOVERNING EQUATIONS
The equation of motion for incompressible electrically conducting nonNewtonian fluid can be written as
With stressstrain relationship is given by,
Together with boundary conditions,
Where Reynolds number and magnetic parameter
Substitute these quantities in equation (1) to (5) and dropping the asterisk, for simplicity
We
With stressstrain relationship is given by,
Special Issue – 2015
International Journal of Engineering Research & Technology (IJERT)
ISSN: 22780181
NCETRASECT2015 Conference Proceedings
Introducing stream function such that,
By solving above equations with , we get
Equation of continuity (7) gets satisfied identically, Equation (8) and (9) becomes,
With boundary conditions,
Introducing equation (21) into equation (15) result in
With the boundary conditions, equation (13) becomes
Introducing equations (22) in equation (11)(13), we get following similarity equation
ButU is independent of y, G() must be constant
.Therefore G() assume Unity
i.e. G() 1, G() 0 And also assume that S0 1
With the boundary conditions, And the StressStrain functional relationship is given by,

SIMILARITY SOLUTION OF THE PROBLEM By using scaling linear group transformation
for the dependent and independent variables. From equation

one obtains


NUMERICAL SOLUTION OF PROBLEM
The transformed highly nonlinear ordinary differential equation (25) subject to the boundary conditions (26) is solved numerically by using KellerBox method (Cebeci and Bradshaw [17]; Keller [18]). This method is second order accurate and allows uniform and nonuniform grid size. The numerical algorithm to solve the twopoint nonlinear boundary value problem is as explained below.
The boundary value problem of (24) in f is reduced to a first order system of three simultaneous ordinary differential equations:
(15), into the Eqs. (1112) results in
equation
dfi d (28)
Fi , f1, f2,… fn ;
i 1, 2, 3 n
Where,
f1 f ,
f2 f ',
f3 f "
The boundary conditions (24) become,
f1 0 0, f2 0 0, f2 U0
Then, after choosing
, the numerical infinity, a grid for
And
The differential equation are completely invariant to the proposed linear transformation, the following coupled algebraic equations are obtained
the closed interval 0, is chosen and the above system of first order equations are transformed into a system of finite difference equations (FDEs) by replacing the differential terms by forward difference approximation and the non differential terms by the average of two adjacent grid points. The numerical method gives approximate values of f , f , f at all the grid points. By adding the boundary
conditions (28) to the system of FDEs, we obtain a nonlinear system of algebraic equations in which the number of equations and unknowns are the same. Subsequently, the
Special Issue – 2015
International Journal of Engineering Research & Technology (IJERT)
ISSN: 22780181
NCETRASECT2015 Conference Proceedings
linearization of these FDEs was done by Newton's method. The resulting systems of linear equations were solved by a
block tridiagonal solver. The step size in and the
position of the edge of the boundary layer in are to be adjusted for different values of the parameters to maintain accuracy. For brevity, further details on the solution process are not presented here. It is worth mentioning that a uniform grid
of 0.01was found to be satisfactory for a convergence
criterion of 106 in all most all the cases.

RESULTS AND DISCUSSION
Using the algorithm discussed above (which is unconditionally stable), the numerical solutions are obtained for several series of values of the magnetic parameter Mn, the flow parameters, and the numerical values are plotted in
Figure 3: Comparison of wedge velocity for various
Figs. 27. Without loss of generality U0
is taken as unity.
These figures depict the influence of flow parameters under effect of magnetic field on the velocity component along the wedge of surface, socalled the wedge velocity (Figs. 2, 3) and on the variation of local shearstress, hence skin friction at the wedge (Figs. 4, 5). Fig. 6, 7 shows the magnetic effect on wedge velocity profile and shearstress at wedge. Figs. 2,3 are the graphical representation of the velocity profiles f along the wedge for different values of and . It shows
the effect of magnetic field on velocity profile. These figures distinguish the velocity with and without magnetic fields. It is worth to observe that the velocity along the wedge increases under the effect of magnetic field for the given values of the flow parameters.
Figure 2: Comparison of wedge velocity for various
Figure 4: Comparison of local shearstress for various
Figure 5: Comparison of local shearstress for various
Special Issue – 2015
International Journal of Engineering Research & Technology (IJERT)
ISSN: 22780181
NCETRASECT2015 Conference Proceedings

It is interesting to note that the deductive group theoretic method based on general group of transformation is applied to derive proper similarity transformations for the nonlinear partial differential equation with the stressstrain functional relationship condition, governing the flow under consideration.

It is to be observing that similarity solutions for all non
Newtonian fluids exist only for the flow past 90 wedge.

The present similarity equation is solved by KellerBox method.

The numerical solutions have produced and compare for fluid flow with and without magnetic field.

It is worth to note that all solutions have derived for non dimensional quantities and hence these results are applicable for all types of under considered non Newtonian fluids.
Figure 6: Influence of magnetic field on wedge velocity
Figure 7: Effect of magnetic field on local shearstress
Figs. 4,5 are represent the comparison of shear stress
f 0
REFERENCES

W.F. Ames, Nonlinear partial differential equations in engineering, Academic press, New York 1969.

C.E. Bluman and J.D. Cole, Similarity methods for differential equations, SpringerVerlag, New York 1974.

A.G. Hansen, Similarity analysis of boundary value problems in engineering, Prentice Hall Inc., N. J. 1964.

R. Seshadri and T.Y. Na, Group Invariance in Engineering Boundary Value Problems, SpringerVerlag, N. Y. 1985.

H. Stephani, Differential Equations: Their Solution using Symmetries, Univ. Cambridge Press, England 1989.

V.J. Rossow, On the flow of electrically conducting fluids over a flat plate in the presence of transverse magnetic field, NASA NT 3971,1957.

H.P. Greenspan and G.F. Carrier, Magneto fluid dynamic flow past a flat plate, JFM, vol. 6, 7782, 1959.

M. G. Timol and A. G. Timol, Three dimensional magnetofluid dynamic flow with pressure gradient and fluid injection , Ind. J. pure and App. Maths, vol. 19, no. 2, 208212, 1988.

M. Patel and M.G. Timol, Similarity solutions of three dimensional MHD boundary layer flow of nonNewtonian fluids, ISH Journal of Hydraulic Engineering, Special Issue, Vol.14, no.3. pp. 8293, 2008.

M. Patel and .G. Timol, Similarity Analysis in MHD Heat and Mass Transfer of NonNewtonian Power Law Fluids Past a Semiinfinite Flate Plate, International Journal of Applied Mechanics and Engineering, vol.15, pp 12051214, 2010.

K. R. Rajagopal, A. S. Gupta and T. Y. Na, A note on the Falkner Skan flows of a nonNewtonian fluid, Int. J. NonLinear Mechanics ,
at wedge with and without magnetic field. It suggests that
under the influence of magnetic field the shear stress at wedge decreases and hence the local skin friction coefficient C f decreases for given values of the flow parameters. These figures are also depicting the effect of flow parameters on local skin friction.
Fig. 6,7 depicts the magnetic effect on fluid flow. It is
interesting to notice that as the magnetic field strength increase there is sharp fall down in the shear stress at wedge and hence the local skin friction decreases sharply. This shows the essential effect of magnetic field on non Newtonian PrandtlEyring fluids.


CONCLUDING REMARKS:

The similarity solutions for MHD laminar incompressible boundary layer equations of all non Newtonian PrandtlEyring fluids are derived.
Vol. 18 (4) (1983), 313320.

A. G. Hansen and T. Y. Na, Similarity solutions of laminar incompressible boundary layer equations of nonNewtonian fluids ASME J. basic eng. (1968), 7174.

S. Y. Lee and W. F. Ames, Similarity solutions for nonNewtonian fluids, A.I. Ch.E.J. 12 (1966), 700.

M. G. Timol and N.L Kalthia., Similarity solutions of three dimensional boundary layer equations of non Newtonian fluids, Int. J. NonLinear Mechanics, 21(6) (1986), 475481.

G. Birkhoff, Hydrodynamics, Princeton University Press, Princeton, NJ. (1960).

A. J. A. Morgan, The reduction by one of the number of independent variables in some systems of partial differential equations, Quart. J. Math.3 (1952), 250259.

T. Cebeci, P. Bradshaw, Physical and Computational Aspects of Convective Heat Transfer, SpringerVerlag, New York, 1984.

H.B. Keller, Numerical Methods for Twopoint Boundary Value Problems, Dover Publ., New York, 1992.