aa r X i v : . [ phy s i c s . c l a ss - ph ] J u l Multiple solutions of steady MHD flow of dilatantfluids Z AKIA H AMMOUCH ∗ LAMFA, CNRS UMR 6140, Universit´e de Picardie Jules VerneFacult´e de Math´ematiques et d’Informatique, 33 rue Saint-Leu 80039 Amiens, France
Abstract.
In this paper we consider the problem of a steady MHD flow of a non-Newtonian power-law andelectrically conducting fluid in presence of an applied magnetic field. The boundary layer equa-tions are solved in similarity form via the Lyapunov energy method, we show that this problem hasan infinite number of positive global solutions.
Keywords:
Asymptotic solution; Boundary-layer; Degenerate differential equation; MHD flow;Power-law fluid; Similarity solution.
MSC:
1. Introduction
The study of non-Newtonian fluid flows has considerable interests, this is primarily because of thenumerous applications in several engineering fields. Such processes are wire drawing, glass fiberand paper production, crystal growing, drawing of plastic sheets etc. For more details about thebehavior in both steady and unsteady flow situations, together with mathematical models, we referthe reader to the books [1] by Astarita and Marucci, [2] by Bohme and the references therein. Oneparticular non-Newtonian model which has been widely studied is the Ostwald-de Wael power-lawmodel [3][4], which relies the shear stress to the strain rate u y by the expression τ xy = k | u y | n − u y , (1.1)where k is a positive constant, and n > is called the power-law index . The case n < is referredto pseudo-plastic or shear-thinning fluid, the case n > is known as dilatant or shear-thickeningfluid. The Newtonian fluid is a special case where the power-law index n is equal to one. In thepresent work we shall restrict our study to the case n > .The magnetohydrodynamics (MHD) flow problems find also applications in a large variety ofphysical, geophysical and industrial fields [5]. It is also interesting to study the flow of non-Newtonian fluids with externally imposed magnetic fields. To the author knowledge MHD flowof non-Newtonian fluids was first studied by Sarpkaya [6]. In [7] Sapunkov derived the equationsdescribing the similarity solutions for the non-Newtonian flow when the external applied magneticfield varies as x m − , in presence of a pressure gradient, he used the method of series expansion.Later, Djuvic [8] employed a Crocco’s variables to study the unsteady flow with exponentially ex-ternal velocity (in time). Recently, Liao [9] introduced a powerful technique (homotopy analysis)to give analytic solutions of MHD viscous flows of non-Newtonian fluids over a stretching sheet. ∗ Email adress: [email protected]
1n this paper, we reconsider the steady two-dimensional laminar flow of an incompressible viscouselectrically conducting dilatant fluid over a stretching flat plate with a power-law velocity distribu-tion in the presence of a perpendicular magnetic field. Our interest in this work has been motivatedby the work of Chiam [10], who have considered the flow over an impermeable flat plate, for whichsimilarity solutions were found via the Crocco transformation.
2. Derivation of the model
Consider a steady two-dimensional laminar flow of an incompressible dilatant and electrically con-ducting fluid of density ρ , past a semi-infinite flat plate. Let ( x, y ) be the Cartesian coordinates ofany point in the flow domain, where x − axis is along the plate and y − axis is normal to it. Assumethat a magnetic field H ( x ) , is applied normally to the plate.The continuity and momentum equations can be simplified, within the boundary-layer approxima-tion, into the following equations (see[7][10]) u x + v y = 0 , (2.1) uu x + vu y = ν ( | u y | n − u y ) y + u e u ex + σµ H ρ ( u e − u ) . (2.2)Accompanied by the boundary conditions u ( x,
0) = U w ( x ) , v ( x,
0) = V w ( x ) and u ( x, y ) → u e ( x ) as y → ∞ . (2.3)Where the functions u and v are the velocity components in the x and y directions respectively, u e ( x ) = U ∞ x m is the free-stream velocity. The parameters ν, n, µ, σ and H are the kinematicviscosity, the flow behavior index, the magnetic permeability, the electrical conductivity of thefluid, and the magnetic field intensity respectively. The functions U w ( x ) = u w x m ( u w > and V w ( x ) = v w x m (2 n − − nn +1 are the stretching and the suction/injection velocities respectively.In term of the stream-function ( ψ which satisfied u ( x, y ) = ψ y ( x, y ) and v ( x, y ) = − ψ x ( x, y ) ),equations (2.1),(2.2) can be reduced to the single equation ψ y ψ xy − ψ x ψ yy = ν ( | ψ yy | n − ψ yy ) y + u e u ex + σµ H ρ ( u e − ψ y ) , (2.4)subject to ψ y ( x,
0) = u w x m , ψ x ( x,
0) = − v w x m (2 n − − nn +1 and ψ y ( x, y ) = U ∞ x m as y → ∞ . (2.5)According to Sapunkov [7], similarity solutions for problem (2.4),(2.5) exist only if the magneticfield has the following form H ( x ) ∼ x m − .To look for similarity solutions we define the following η := Ayx − a and ψ ( x, y ) := Bx b f ( η ) , (2.6)where f is the transformed dimensionless stream function and η is the similarity variable. Thanksto (2.6), the function f satisfies the new boundary value problem ( | f ′′ | n − f ′′ ) ′ + af f ′′ + m (1 − f ′ ) + M (1 − f ′ ) = 0 ,f (0) = α, f ′ (0) = δ, f ′ ( ∞ ) = 1 , (2.7)2f and only if a = 1 + m (2 n − n + 1 , b = 1 + m ( n − n + 1 , a − b = m, and the parameters A and B satisfy AB = u ∞ and νB n − A n − = 1 . Where the primes denote differentiation with respect to η , the function f ′ ( η ) denotes the normal-ized velocity and the parameters M = σµ H ( n + 1) u ∞ ρ , α = − ( n + 1) v w ( m + 1)( νu n − ∞ ) n +1 and δ = u w u ∞ , are respectively: The Hartmann number, the suction/injection and the stretching parameters.Such problems have been investigated by several authors for example, Anderson et al. [11], Zhanget al. [12] and Kumari and Nath [13].In the same context, Chiam [10] studied Problem (2.1)-(2.3). To look for similarity solutions, hesolved the following boundary value problem n | f ′′ | n − f ′′′ + f f ′′ + β (1 − f ′ ) + M (1 − f ′ ) = 0 ,f (0) = 0 , f ′ (0) = 0 , f ′ ( ∞ ) = 0 . (2.8)Where β = m ( n +1)(2 n − m +1 . We aim here to stress that for n = 1 , equation (2.7) can be degenerate atsome point η s for which f ′′ ( η s ) = 0 (for more details see [15]) and then any solution of (2.7) is notnecessarily of C (0 , ∞ ) . Hence equations (2.7) and (2.8) are not equivalent.Let us notice that for the Newtonian case ( n = 1) , problem (2.7) reduces to the Falkner-Skan flowin Magnetohydrodynamics, which has been studied by Hildyard [17], Aly et al. [18] and Hoernel[19]. The case m = M = 0 leads to the generalized Blasius problem (see [20]). While the case m = − M , by a suitable scaling, is referred to the mixed convection of a non-Newtonian fluid in aporous medium (see for example [21]). We note also that in absence of the magnetic field, problem(2.7) is simplified to the Falkner-Skan flow for non-Newtonian fluids. A complete study on thissubject is given in [22] by Denier and Dabrowski.Very recently, Aly et al. [18] reported a theoretical and numerical investigations on the existenceof solutions to problem (2.7) for Newtonian fluids ( n = 1) , say f ′′′ + m +12 f f ′′ + m (1 − f ′ ) + M (1 − f ′ ) = 0 ,f (0) = α ≥ , f ′ (0) = δ, f ′′ (0) = γ. (2.9)They showed that problem (2.9) has multiple solutions for any δ ∈ (0 , Γ) and γ ∈ R satisfying γ ≤ m δ + M δ − M + m ) δ, (2.10)where Γ = − M m " r m M ( m + M ) > . In the present work, we aim to extend theirresults to the non-Newtonian dilatant fluids ( n > , by using a condition on γ which is differentfrom (2.10) and without any restriction on the parameter δ .3 . Non-uniqueness of solutions Guided by the analysis of [14],[15] and [16], we aim to prove the existence of solutions to problem(2.7), for related values of the parameters m, M, n, α, δ and γ . This result will be established bymean of the so-called shooting method, the boundary value problem (2.7) is then converted intothe following initial value problem ( | f ′′ | n − f ′′ ) ′ + af f ′′ + m (1 − f ′ ) + M (1 − f ′ ) = 0 ,f (0) = α, f ′ (0) = δ, f ′′ (0) = γ. (3.1)Where the real number γ is the shooting parameter.The initial value problem (3.1) can be transformed into the equivalent first order ordinary differen-tial system f ′ = g,g ′ = | h | − nn h,h ′ = − af | h | − m (1 − g ) − M (1 − g ) , (3.2)with the conditions f (0) = α, g (0) = δ, h (0) = | γ | n − γ. (3.3)By the classical theory of ordinary differential equations, problem (3.2),(3.3) has a unique local(maximal) solution for every γ = 0 . Let f γ denotes this solution and (0 , η γ ) , η γ ≤ ∞ , denotes itsmaximal interval of existence. The main task now is to show how existence of solutions dependson γ .The local solution f γ satisfies the following | f ′′ γ | n − f ′′ γ + af ′ γ f γ − M ( f γ + α ) = | γ | n − γ + aαδ − ( M + m ) η + ( a + m ) Z η f ′ γ ( τ ) dτ. (3.4)Equation (3.4) will be used for proving the main results.D EFINITION f γ is said to be a solution to (3.1) if f ∈ C (0 , ∞ ) , | f ′′ γ | n − f ′′ γ ∈ C (0 , ∞ ) and satisfies lim η →∞ f ′ γ ( η ) = 1 (i) and lim η →∞ f ′′ γ ( η ) = 0 (ii) α ∈ R ) T HEOREM α ∈ R , δ > , M > , n > and − n < m < − M . For any γ satisfying | γ | n − γ > − aαδ ( ⋆ ) , problem(3.1)admitsaglobalunboundedsolution. Proof.
From a physical point of view, it is more convenient to prove the result for the cases α ≥ (suction) α < (injection) separately.We have to show that f γ is a positive monotonic increasing function on (0 , η γ ) , globally definedand going to infinity with η . For this sake we define the Lyapunov Energy function by V ( η ) = 1 n + 1 | f ′′ | n +1 − m f ′ − M f ′ + ( M + m ) f ′ . (3.5)4hich satisfies V ′ ( η ) = − af f ′′ . Then V is monotonic decreasing on (0 , η γ ) . On the other hand, from equation (3.4) and condition ( ⋆ ) we see that f γ ′ and f γ are positive on (0 , η γ ) as long as f γ exists. Using the Lyapunov function V we see that f γ ′ and f γ ′′ are bounded, since V is bounded from below by M +4 m . If f γ were alsobounded, say f → η ∞ L with L ∈ (0 , ∞ ) (since f γ is positive). Then f γ ′ ( η ) → η ∞ which impliesthat f ′′ γ ( η k ) → k ∞ , where ( η k ) k ≥ is a sequence tending to infinity with k . Using again (3.4) todeduce | f ′′ γ ( η k ) | n − f ′′ γ ( η k ) + af ′ γ ( η k ) f γ ( η k ) = − M ( f γ ( η k ) + α ) + | γ | n − γ + aαδ − ( M + m ) η k + ( a + m ) Z η k f ′ γ ( τ ) dτ. Letting k → ∞ , the right hand side goes to zero while the left hand side goes to minus infinity,which is impossible. then f γ is a global unbounded solution to (3.1).From the above f ′ γ and f ′′ γ are bounded and f ′ γ is monotonic increasing on ( η , ∞ ) , for η largeenough. Then there exists l > such that lim η →∞ f ′ γ ( η ) = l, and there exits a sequence ( ζ k ) k ,tending to infinity with k such that lim k →∞ f ′′ γ ( ζ k ) = 0 . Making recourse to the Lyapunov function V we get lim η →∞ f ′′ γ ( ζ k ) = 0 .Assume now that f ′′ γ is not monotonic on any interval [ η , ∞ ) . Then, there exists a sequence ( τ k ) k going to infinity with k such that: • ( | f ′′ γ | n − f ′′ γ ) ′ ( τ k ) = 0 , • | f ′′ γ | n − f ′′ γ ( τ k ) is a local minimum , • | f ′′ γ | n − f ′′ γ ( τ k +1 ) is a local maximum . From (3.1) , we have f ′′ γ ( τ k ) = − m (1 − f ′ γ ( τ k )) + M (1 − f ′ γ ( τ k )) af γ ( τ k ) . Since f ′ γ is bounded and f γ goes to infinity with η k , we get easily from the above that f ′′ γ goes tozero with η .Now we show that f γ satisfies (i). Recall that f ′ γ is a positive bounded function then f ′ γ → l with l ∈ ( α, ∞ ) . At infinity we have f γ ∼ ηl and from identity (3.4) we get | f ′′ γ | n − f ′′ γ ∼ η [ ml + M l − ( M + m )] + o (1) as η approaches infinity, this leaves only the possibility that l is either or − Mm − , thanks to thepositivity of f ′ γ we deduce that l = 1 .To finish we show the result for α < . In such case, the function f γ is negative on a smallneighborhood of zero. According to (3.4) f γ cannot have a local maximum, then two possibilitiesarise: • Either f γ < ∀ η ∈ (0 , η γ ) • Or ∃ η ⋆ such that f γ < on (0 , η ⋆ ) , f γ ( η ⋆ ) = 0 and f γ > ∀ η > η ⋆ . Assume that the first assertion holds, then α > f γ ( ∞ ) ≤ and f ′ γ ( ∞ ) = 0 , f ′ being positive weuse again (3.4) to get that f ′′ γ is positive. A contradiction. Then, f γ has exactly one zero η ⋆ . Wedefine the shifted function h by : η h ( η ) = f γ ( η + η ⋆ ) , which satisfies h (0) = 0 , h ′ (0) = δ and h ′′ (0) > , and we use the above analysis to conclude that h is an unbounded global solution to (3.1). (cid:4) .2. Reversed flows ( δ < Now we pay attention to the case of reversed flows ( δ < . First, we show that the shootingparameter has to be positive.P
ROPOSITION f γ be a solutionto (3.1) with m ∈ ( − n , − M ) , α < , δ < and γ ≤ ,thenthecondition(i)is failed. Proof.
Let δ < , if γ ≤ then f γ ′′ is negative on some (0 , η ) , for η small, and equation (3.1)can be written as ( f γ ′′ e F ) ′ = − e F n | f ′′ γ | − n (cid:2) m (1 − f ′ γ ) + M (1 − f ′ γ ) (cid:3) , where F ( η ) = an Z η f γ | f γ ′′ | − n dτ . From this we see that η f ′′ e F decreases and then f ′′ γ ( η ) ≤ for all η ∈ (0 , η γ ) . It follows that f ′ γ is decreasing on (0 , η γ ) and then the condition (i) could notbe satisfied. (cid:4) T HEOREM δ < , α > and m ∈ ( − n , − M ) . Forany γ > satisfying αγ n − δ γ n − + aα δ − M α > ⋆⋆ ) , problem(3.1)hasa globalunboundedsolution. Proof.
Let f γ be the local solution of (3.1), define the auxiliary function G ( η ) = f γ f γ ′′ | f ′′ γ | n − − f ′ γ | f ′′ γ | n − + af γ f ′ γ − M f γ , (3.6)which satisfies G ′ ( η ) = − ( m + M ) f γ + (cid:20) a + m + ( n − a n (cid:21) f ′ γ f γ + n − n f ′ γ f ′′ γ − h m (1 − f ′ γ ) + M (1 − f ′ γ ) i , (3.7)and G (0) > .Since f ′ γ < , the function f γ is negative on a small neighborhood of zero. Assume that there exists η ∈ (0 , ∞ ) such that f γ ( η ) > , f ′ γ < ∀ η ∈ [0 , η ) and f γ ( η ) = 0 . Hence G is a monotonic nondecreasing function on (0 , η ) and then G ( η ) ≤ . Then G ( η ) ≤ for all η ∈ (0 , η ) , in particular G (0) ≤ , which is a contradiction with ( ⋆⋆ ) . Therefore we have : • Either f γ > and f ′ γ ≤ ∀ η ≥ • Or ∃ η > f γ > , f ′ γ < ∀ η ∈ (0 , η ) , f ′ γ ( η ) = 0 and f γ ( η ) is a local maximum . Assume that the first assertion holds, then f γ has a finit limit at infinity, say L ∈ (0 , ∞ ) and thereexists a sequence ( χ k ) k ≥ tending to infinity with k such that f ′ γ ( χ k ) goes to zero at infinity. If f ′ γ is monotonic (resp. non-monotonic on any interval ( η, ∞ ) ) we get f ′ γ goes to zero at infinity andthen f ′′ γ ( δ k ) goes to zero at infinity for a sequence ( δ k ) k ≥ going to infinity with k (resp. f ′′ γ ( δ k ) = 0 and f ′ γ ( δ k ) goes to zero at infinity). Because G (0) < G ( δ k ) , we obtain a contradiction by takingthe limit as k goes to infinity.Now, we claim that the function f γ cannot have a local maximum after η . Actually, assume thereexists η > η such that f γ ( η ) is a local maximum. At this point the function G takes a negativevalue and satisfies G ( η ) ≥ G (0) a contradiction. Since f γ is monotonic increasing after η wededuce as the above that is a global solution.Next, we argue as in the proof of Theorem. 3.1 to show that f γ is unbounded at infinity and satisfies(i) and (ii). (cid:4) .3. Flow with large initial velocity ( δ ≫ ) In this subsection, we construct asymptotic solutions to problem (3.1) when the real δ is very large.Adopting the method used in [23] by Aly et al., we assume that such solutions can be written underthe following form f ( η ) = η + ξ r g ( t ) , where t = ξ s η, ξ = δ − and r, s ∈ R . Then problem (3.1) reads ξ ( r +2 s )( n − ( | g ′′ | n − g ′′ ) ′ + aηξ − r g ′′ + agg ′′ − (2 m + M ) ξ − ( r + s ) − mg ′ = 0 ,g (0) = αξ − r , g ′ (0) = ξ − ( r + s ) , g ′ ( ∞ ) = 0 . (3.8)Setting r = n − n +1 and s = − nn +1 , ensures that the highest derivative remains present in the resultingproblem.As ξ goes to infinity, we deduce ( | g ′′ | n − g ′′ ) ′ + agg ′′ + mg ′ = 0 ,g (0) = 0 , g ′ (0) = 1 , g ′ ( ∞ ) = 0 . (3.9)Problem (3.9) describes the steady free convection flow of a non-Newtonian power-law fluid overa stretching flat plate embedded in a porous medium. In [15], it was shown that for m ∈ ( − n , any local solution g , whith positive values of τ ( τ = g ′′ (0)) , is global and satisfies the followingasymptotic behaviour g ( t ) ∼ t m (2 n − m ( n − , as t → ∞ . Consequently, a solution f for positive γ and large δ (if it exists), may have the following large η -behaviour f ( η ) ∼ η h δ − m ( n − η m ( n − m ( n − i .
4. Concluding remarks
Based on the similarity transformation approach, the boundary layer equations for flows of purelyviscous non-Newtonian dilatant and electrically conducting fluids are investigated. Using a shoot-ing argument, it is shown that the relevant problem admits an infinite number of solutions ([24][25][26]and [27]), this is due to the arbitrariness of the shooting parameter γ . From a physical point of view,we underline that γ = f ′′ (0) originates from the local skin friction coefficient C f x , and the localReynolds number Re x = ( u w x m ) − n x n νk via the the formula C f x Re x n +1 = 2 (cid:16) an (cid:17) n +1 | γ | n − γ. In conclusion, we may expect that the solutions determined above are physically acceptable. How-ever, only experiments are able to prove their physical existence.
Acknowledgements
I would like to thank the anonymous referee for his constructive suggestions, which have improvedthe earlier version of this work. 7 eferences [1] A
STARITA
G., and M
ARRUCCI
G., Principles of Non–Newtonian Fluid Mechanics,McGraw-Hill, 1974.[2] B
OHME
H., Non–Newtonian Fluid Mechanics, North–Holland Series in Applied Mathemat-ics and Mechanics, 1987.[3] A
CRIVOS
A., S
HAH
M.J., P
ETERSEN
E.E, Momentum and heat transfer in laminarboundary-layer flows Equations of non-Newtonian fluids past external surfaces, A.I.Ch.EJ., 6 312 (1960).[4] E CE M.C. , B ¨ UY ¨ UK E., Similarity solutions for free convection to power-law fluids from aheated vertical plate, App. Math. Lett., 15 1-5 (2002).[5] P
AVLOV
K. B., Magnetohydrodynamic flow of an incompressible viscous fluid caused bydeformation of a surface, Magnitnaya Gidrodinamika, 4 146-147 (1974).[6] S
ARPKAYA
T., Flow of non-Newtonian fluids in a magnetic field, AIChE J. 7 324-328 (1961).[7] S
APUNKOV
S.Y., Self-similar solutions of non-Newtonian fluid boundary in MHD,Mekhanika Zhidkosti I Gaza 2 77-82 (1967).[8] D
JUVIC
D. S., Hiemenz magnetic flow of power-law fluids, ASME J. Appl. Mech., 41 822-823 (1974).[9] L
IAO
S. J., On the analytic solution of magnetohydrodynamic flows of non-Newtonian fluidsover a stretching sheet, J. Fluid Mech., 488 189-212 (2003).[10] C
HIAM
T. C., Solutions for the flow of a conducting power-law fluid in a transverse magneticfield and with a pressure gradient using Crocco variables, Acta Mech., Vol. 137 225 - 235(1999).[11] A
NDERSON
H. I., B
ACH
K.H., D
ANDAPAT
B.S., Magnetohydrodynamic flow of a power-law fluid over a stretching sheet, Int. J. Non-Linear Mech., 27 929-939 (1992).[12] Z
HANG
Z., W
ANG
J., S HI W., A boundary layer problem arising in gravity driven lami-nar film of power-law fluids along vertical walls, ZAMP. (J. Appl. Math. Phy.), 55 769-780(2004).[13] K
UMARI
M., N
ATH
G., MHD boundary layer flow of a non-Newtonian fluid over a contin-uously moving surface with a parallel free stream, Acta. Mech., 146 139-150 (2001).[14] B
RIGHI
B., On a similarity boundary layer equation, ZAA(Zeitschrift fr Analysis und ihreAnwendungen), 21 931-948 (2002).[15] H
AMMOUCH
Z., ´Etude math´ematique et num´erique de quelques probl`emes issus de la dy-namique des fluides, PhD Dissertation, Universit´e de Picardie Jules Verne, Octobre 2006.[16] G
UEDDA
M., Similarity solutions of differential equations for boundary layer approxima-tions in porous media, ZAMP. (J. Appl. Math. Phy.), 56 749-762 (2005).[17] H
ILDYARD
L., Falkner-Skan Problems in Magnetohydrodynamics, The Phys. of Fluids,Vol.15 782–793 (1972). 818] A LY E. H., B
ENLAHSEN
M., G
UEDDA
M., Similarity solutions of a MHD boundary-layerflow past a continuous moving surface, Int. J. Engng. Sci. 45 486-503 (2007).[19] H
OERNEL
J-D., On the similarity solutions for a steady MHD equation, Com. Nonlin. Sci.Num. Sim., 13 1353-1360 (2008)[20] N
ACHMAN
A., T
ALIAFERRO
S., Mass transfer into boundary-layers for power-law fluids,Proc. R. Soc. Lon. A, 365 313-326 (1979).[21] P OP I., Mixed convection to power-law type non-Newtonian fluids from a vetical wall, Pol.Plast. Tech and Eng., 30 47-65 (1991).[22] D
ENIER
J.P., D
ABROWSKI
P.P., On the boundary-layer equations for power-law fluids, Proc.R. Soc. Lond. A, 460 3143-3158 (2004).[23] A LY E.H, E
LLIOTT
L., I
NGHAM
D.B., Mixed convection boundary-layer flow over a verti-cal surface embedded in a porous medium, Eur. J. Mech. B Fluids, 22 529-543 (2003).[24] M
AGYARI
E., A LY E.H., Mechanical and thermal characteristics of a mixed convectionboundary-layer flow in a saturated porous mediuM Int. J. Heat Mass Trans., 49 3855-3865(2006).[25] L
IAO
S.J., A new branch of solutions of boundary layer flows over a permeable stretchingplate Int. J. Nonlinear. Mech., 42 819–830 (2007).[26] R
ILEY
N., W
EIDMAN
P. D., Multiple Solutions of the Falkner-Skan Equation for Flow Pasta Stretching Boundary, SIAM J. Appl. Math., 49 1350-1358 (1989).[27] R