International Journal of Fluids Engineering. ISSN 974-8 Volume 5, Number (), pp. 57-76 International Research Publication House http://www.irphouse.com Unsteady Mhd Flow of a Non-Newtonian Fluid Down and Open Inclined Channel with Naturally Permeable Bed K. Bhagya Lakshmi, G. S. S. Raju and N. V. R. V. Prasad Department of Mathematics, CMRTechnical Campus, kandlakoya, Hyderabad, A.P., India. Department of Mathematics, J N T U A College of Engg. Pulivendula, Y.S.R Dist., A.P., India. Department of Mathematics, S.V.G.S. Junior College, Nellore, A.P., India. E-mail: rajugss@yahoo.com Abstract: This paper deals with unsteady MHD flow of a Walters fluid (Model B ) an open inclined channel of width a and depth d under gravity, with naturally permeable bed, the walls of channel being normal to the surface of the bottom, under the influence of a uniform transverse magnetic field. The free surface is exposed to atmospheric pressure. A uniform tangential stress is applied at the free surface in the direction of flow. The naturally permeable bottom of the channel is taken at an angle β with the horizontal. Flow of fluid both in porous medium and in free fluid region is studied with the same pressure gradient. The exact solution of velocity distribution has been obtained by using Laplace transform and finite Fourier sine transform techniques. We have evaluated the velocity distribution and the flux of the fluid in different cases of time dependent pressure gradient g(t) viz., i) constant, ii) exponentially decreasing function of time and iii) cosine function of time. The effects of magnetic parameter M, viscoelastic parameter K, permeable parameter K, Reynolds number R and time t on velocity distribution u and flux C f in three different cases are investigated. Key words: Non-Newtonian fluid, Open inclined channel, Porous medium, Magnetic field.
58 K. Bhagya Lakshmi et al INTRODUCTION The flow of a liquid in an open inclined channel with a free surface has wide applications in the designs of drainage, irrigation canals, flood discharge channels and coating to paper rolls etc. Vanoni [] has evaluated velocity distribution in open channels. Johnson [6 ] has studied the rectangular roughness in the open channel. SatyaPrakash [7] considered analytically a viscous flow down an open inclined channel with plane bottom and vertical walls under the action of gravity. The free surface was exposed to atmospheric pressure and bottom was taken as impermeable. Bakhmeteff [], Henderson [5] and Chow [] have discussed many types of open channel flows. Gupta etal [4] have studied the flow of a viscous fluid through a porous medium down an open inclined channel. VenkataRamana&Bathaiah [] have studied the flow of a hydro magnetic viscous fluid down an open inclined channel with naturally permeable bed under the influence of a uniform transverse magnetic field. Unsteady laminar flow of an incompressible viscous fluid between porous and parallel flat plates have been investigated by Singh [8], taking (i) both plates are at rest (ii) Generalized plane coutte flow. Rheology is the science of deformation and flow of matter. The aim of rheology is to predict the deformation on flow resulting from the application of a given force system to a body or vice versa. The subject of rheology is of technological importance as in many branches of industry, the problem arises of designing apparatus to transport or process substances, which can't governed by the classical stress, strain velocity relations. In the manufacture of rayon, nylon or other textiles fibers, viscoelastic effects are encountered when the spinning solutions are transported or forced through spinnerets and in the manufacture of lubricating greases and rubber. Further viscoelastic fluid occurs in the food industry, e.g. emulsions, pastes and condensed milk. Non-Newtonian fluids have wide importance in the present day technology and industries. The Walters fluid is one of such fluid. The constitutive equations governing motion of Walters fluid (Model B') are P = pg + P. P = η e K e. The equation of motion and continuity are ρ ν = P + P. ν.4 wherep ik is the stress tensor, is the density, p is the pressure, g ik is the metric tensor of a fixed coordinate system x', i is the velocity vector, in the contra variant form is e = + ν e ν e ν e.5 It is converted derivative of the deformation rate tensor e ik defined by e = ν + ν.6
Unsteady Mhd Flow of a Non-Newtonian Fluid Down 59 Here, h is the limiting viscosity at small rates of shear which is given by η = N(τ)dτ.7 K = τn(τ)dτ.8 N(t) being the relaxation spectrum as introduced by Walters [, 4]. This idealized model is a valid approximation of Walters s fluid (Model B') taking very short memory into account so that terms involving τ N(τ)dτ, n.9 Have been neglected. In this paper an attempt has been made to study the unsteady MHD flow of a Walters fluid (Model B') down an open inclined channel under gravity of width 'a' and depth 'd', with naturally permeable bed, the walls of the channel being normal to the surface of the bottom, under the influence of a uniform transverse magnetic filed. A uniform tangential stress is applied at the free surface in the direction of flow. The naturally permeable bottom of the channel is taken at an angle 'b' with the horizontal. We have evaluated the velocity distribution by using Laplace transform and finite Fourier Sine transform techniques. We have evaluated the velocity distribution and flux of the fluid in different cases of time dependent pressure gradient g (t), viz., (i) The fluid flows in the steady state for t (ii) Unsteady state occurs at t > and (iii) Unsteady motion is influenced by time-dependent pressure gradient. The velocity distribution and flux have been obtained in some particular cases i.e., when (i) g(t) = C (ii) g(t) = C e -bt (iii) g(t) = C Cos bt; where C and b are constants. The effects of magnetic parameter 'M', viscoelastic parameter K, porosity parameter K, Reynolds number 'R' and time 't' are investigated on the velocity distribution and the flux of the fluid.. FORMULATION & SOLUTION OF THE PROBLEM The unsteady MHD flow of a Walters fluid (Model B') down an open inclined channel of width 'a' and depth 'd' under gravity, with naturally permeable bed, the walls of channel being normal to the surface of the bottom, under the influence of a uniform transverse magnetic field. A uniform tangential stress 'S' is applied at the free surface. The free surface is exposed to atmospheric pressure. The naturally permeable bottom of the channel is taken at an angle 'b' (o<bp/) with the horizontal. The x-axis is taken along central line in the direction of the flow at the free surface, y-axis along the depth of the channel and the z-axis along the width of the channel. A uniform magnetic field of intensity 'H o ' is introduced in the y-direction. Therefore, the velocity and magnetic fields are given by and. The fluid being slightly conducting, the
6 K. Bhagya Lakshmi et al magnetic Reynolds number much less than unity, so that the induced magnetic field can be neglected in comparison with the applied magnetic field [Sparrow &Cess()]. Fluid flow in porous medium is governed by Darcy's law and fluid in free flow region is governed by Naveir-stoke's equations. In the absence of any input electrical field the equations of continuity and motion of the unsteady MHD Walters fluid ( Model B') flowing down an open inclined channel of t> are =. ρ K σμ H u. + ρg cos β =..4 The Darcy's equation for the flow in the porous medium is Q = + ρg sin β.5 Where ρis the density of the fluid. g is the acceleration due to gravity p is the pressure μis the coefficient of viscosity K is the non-newtonian parameter K is the permeability of the medium Q is the velocity in the porous medium μ is the magnetic permeability σis the electrical conductivity of the fluid H is the magnetic field The boundary conditions are t, u = u t >, z = ±a, u = y =, μ = S y = d, u = u = (u Q).6 Where u is the initial velocity u is the slip velocity s is the dimensionless constant depending on the porous material S is the uniform tangential stress
Unsteady Mhd Flow of a Non-Newtonian Fluid Down 6 We introducing the following non-dimensional quantities. u = u = Q = K = t = y = z = S = Q = K = p = In the view of the above non-dimensional quantities equations.&.5reduces to ( dropping * ) = R + sin β + + K + Mu.8 Q = KR.9 Where R = (Reynolds number) F = (Froude number) M = (Magnetic parameter) The non-dimensional boundary conditions are t, u = u t >, z = ±L =, u = y =, = SR y =, u = u = sin β KR u. Assuming R + sin β = g(t) att > = Patt. Substituting z = Lin equation.8, we get = g(t) + + K The boundary conditions. reduces to t, u = u t >, ξ =, πandu = y =, = SR y =, u = u + Mu.
6 K. Bhagya Lakshmi et al = [Kg(t) u ]. Since u is the initial velocity i.e.,, by taking g(t)=p in equation., u is given by u = u + Where sh = sinh ch = cosh th = tanh u = B = + αthα q = α = q + M () sin nξ.4 sin nξ.5 Now to solve equation., we take Laplace transform to equation. w.r.t. 't' defined as u(y, ξ, s) = u(y, ξ, t) e dt, s >.6 We get + + () g (s) = g(t) u = + g (s) u u K α.6 e dt.7 The transformed boundary conditions are y =, = y =, u = u = u.8 On taking the finite Fourier sine transform of equation.7 w.r.t. defined as u (y, N, s) = u(y, ξ, t) sin Nξ dξ.9 We get u = + + ( )( ) + () ( )
Unsteady Mhd Flow of a Non-Newtonian Fluid Down 6 + cos Nπ N s P α ch αy ch α + SR α g (s) + H ( K s) ( )( ) sh α( y) + P K q ch α α H K s. The transformed boundary conditions are y =, = y =, u = u = u. On inverting finite Fourier sine transform as given by Sneddon [9] u(y, ξ, s) = u (y, N, s) sin Nξ. In equation. we get u = + () + + () ( ) sin Nξ. On inverting the Laplace transform as defined by Snedden [9] u(y, ξ, t) = u(y, ξ, s) e dt.4 In equation. we obtain u = u + + h(u) Where a = (r + ) A = h(u) = () () () () () g( u)du () () d = sin Nξ.5
64 K. Bhagya Lakshmi et al If we take the limit K, K and M in equation.5 then we get the velocity distribution in the case of viscous non-magnetic and impermeable bed. In this case the velocity distribution is u = R + sin β Sin Nξ.6 Which is in agreement with SatyaPrakash [7]. Now we evaluated the velocity distribution and flux in three different cases viz., Case When g(t) = C, using in equation.5, we get u = u + + + () () () () () () C = u dy dz =u + + () () + Case When g(t) = C e -bt, b>, using in equation.5 we get u = u + () () + () sin Nξ.7.8 () () () sin Nξ.9 () C = u + + () + () Case When g(t) = C Cos bt, using in equation.5 we get +. ch αy u = u ch α + ( ) (r + )π u d e cos a y a + α
Unsteady Mhd Flow of a Non-Newtonian Fluid Down 65 + cos Nπ 4( ) P e cos a y π N (r + )π (a + α SR ) α () sh α( y) ch α + sin Nξ. + () + 4L π cos Nπ N () C = u th α α + u d e a + α 8 P e (r + ) π (a + α ) SR α ch α ch α.. CONCLUSIONS We discuss velocity distribution and flux in three distinct time dependent pressure gradients viz., (i) g(t) = C (ii) g(t) = C e -bt (iii) g(t) = C Cos bt Figures to, 4 to 6 and 7 to 9 are drawn to investigate the effects of velocity distribution 'u' against 'y' for different values of viscoelastic parameter K magnetic parameter 'M' and time 't' respectively in three different cases (i), (ii) & (iii). We noticed that in all the three cases the velocity distribution increases with increase in M or t, where as it decreases with increase in K. Further we observed that the velocity distribution decrease upto y =.8 and then increase with increase in 'y' in figures to. Also we noticed that the velocity distribution 'u' decrease with increase in 'y' in figures 4 to 9. Figures to are prepared to bring out the effects of Reynolds numbers R, viscoelastic parameter K on the flux C f of the fluid. We observed that C f decreases with increase in R or K respectively in all the three cases. In figures to 5, C f is drawn against t for different values of K in cases (i), (ii) & (iii). We noticed that C f decreases with increase in K, where as it increases with increase in 't' in all three cases. From figure 6, one can conclude that C f increases with increase in K or M in case (i). From figures 7 & 8 we conclude that C f decreases with increase in 'M', where as it increases with increase in K respectively in cases (ii) & (iii).from tables to, we see that the velocity distribution is plotted against 'y' fordifferent values of 'K' respectively in three different cases. We noticed that 'u' increase in 'K' where as decreases with increase in 'y' respectively in three different cases. From table 4 we conclude that the slip velocity u B decreases with increase in M, where as it increases
66 K. Bhagya Lakshmi et al with increase in 'K'. Case..8 K =. K =.5 K =.5.6 u.4...4.6.8 y Fig. : u against y for different K Case.4. K =. K =.5 K =.5.8 u.6.4...4.6.8 y Fig. : u against y for different K
Unsteady Mhd Flow of a Non-Newtonian Fluid Down 67 Case.4. K =. K =.5 K =.5.8 u.6.4. 8 7 6 5..4.6.8 Fig.: u against for different K Case 4 y M = M = M = 4 4 u..4.6.8 Fig.4 : u against y for different M y
68 K. Bhagya Lakshmi et al 8 7 6 5 4 Case 5 M = 4 M = M = u..4.6.8 y Fig.5: u against y for different M Case 6 8 7 6 M = 4 M = M = 5 4 u..4.6.8 y Fig.6 : u against y for different M
Unsteady Mhd Flow of a Non-Newtonian Fluid Down 69 5 Case 7 t =. t =.5 t =. 5 u 5..4.6.8 y 5 Fig.7:u against y for different t Case 8 t =. t =.5 t =. u 5 5..4.6.8 y Fig. 8 : u against y for different t
7 K. Bhagya Lakshmi et al Case 9 5 t =. t =.5 t =. 5 u 5..4.6.8 Fig.9 : u against y for different t Case y C f.8.6.4..8.6.4. K = K = K = 4 5 Fig.: C f against R for Different K R
Unsteady Mhd Flow of a Non-Newtonian Fluid Down 7 Case.5 C f.5.5 K = k = k = 4 5 R Fig.: C f against R for different K Case.5.5 C f.5 K = K = K = 4 5 R Fig. : C f against R for different K
7 K. Bhagya Lakshmi et al 5 4 Case K = K = K = 5 C f.5.5.75.5.5 t 8 7 6 Fig; : C f against t for different K Case 4 K = K = K = 5 5 4 C f.5.5.75.5.5 Fig.4 : C f against t for different K t
Unsteady Mhd Flow of a Non-Newtonian Fluid Down 7 8 7 6 5 4 Case 5 K = K = K = 5 C f.5.5.75.5.5 t Fig.5 : C f against t for different K Case 6.5.5 C f.5 M = 4 M = 8 M =.5..5..5 K Fig.6 : C f against K for different M
74 K. Bhagya Lakshmi et al Case 7.5 C f.5.5 M = M = 6 M =.5..5..5 K Fig.7 : C f against K for different M Case 8.5.5 C f.5 M = M = 6 M =.5..5..5 K Fig.8 : C f against K for different M
Unsteady Mhd Flow of a Non-Newtonian Fluid Down 75 Table. Case u against y for different K y K=.5 K =. K =.5 K =. K =.5..45647.4577.45756.4578.4584..999.46.4.47.447.4.44.45.45.48.5.6.964.97.97.9.9.8.557.566.566.59.544..45.59.59.5684.585 Table. Case u against y for different K y K=.5 K =. K =.5 K =. K =.5..497.467.46.44.454..967.9697.9647.96464.96484.4.7.79.7.76.78.6.979.986.99.999.995.8.757.764.768.777.778..4468.4944.59.56.586 Table. Case u against y for different K y K=.5 K =. K =.5 K =. K =.5..489.559.599.66.647..989.996.9.99.949.4.689.686.684.6847.68447.6.7648.778.7757.7784.785.8.744.754.757.757.759..447.4947.565.565.5866 Table.4 u B against K for different M y M = 5 M = M = 5 M = M = 5.5.49.48.474.465.456..4966.4954.494.499.497.5.595.596.59.597.589..5566.556.55585.55569.5555.5.5766.5769.57575.57575.57558
76 K. Bhagya Lakshmi et al REFERENCES [] Bakhmeteff B.A. 9 Hydraulics of Open channels Newyork, McGraw - Hill Book Company. [] Beavers G.S. 967 Boundary Conditions at a naturally Joseph D.D. permeable wall. J. Fluid Mech.,, 97. [] Chow V.T. 959 Open Channel Hydraulics McGraw - Hill, Inc., Newyork. [4] Gupta P.C. 98 Indian Journal. Theory.Phys., 8, 65. Chowdary J.S. Sharma R.G. [5] Henderson F.M. 966 Open channel flow, the Macmillan Company, Newyork. [6] Johnson J.W. 944 Trans. Am. Geophys. Un. 5, 96. [7] SathyaPrakash 97 Indian J. pure appl. Math.,. [8] Singh R. 99 Ph.D. Theasis, Rajasthan University, Jaipur. [9] Sneddon I.N. 974 The use of integral transforms, Tata McGraw - Hill Book Co. Ltd., [] Sparrow E.M. 96 Trans, ASME, J. Appl. Mech. Cess R.D. 9,, 8. [] Vanoni V.A. 94 Velocity distribution in open channels, Civil Engg., II, 56-57. [] VenkataRamana S. 97 Hydro magnetic viscous flow down in an Bathaiah D. open inclined channel with natural permeable bed. Ganit (J.Bangladesh Math Soc., ) Vol., No., 94-. [] Walters K. 96 Quart. J.Math. Appl. Mech., 444 [4] Walters K. 96 J.Mech., 479