Int. f. Engng Sci‘ Vol. 26, No. 4, pp. 313-323, 1988
~2~7225/88
$3.00+ 0.00
Copyright @ 1988 Pergamon Press plc
Printed in Great Britain. AU rights reserved
THE INFLUENCE OF VIBRATING PLATES ON POISEUILLE FLOW OF A BINARY MIXTURE Faculty of Mech~i~~
En~nee~ng, Department of Mechanics, Istanbul Technical Unive~ity, ~~rn~~su~~ Ist~bul, Turkey
Abstmct-A mixture of two incompressible Newtonian fluids is considered and the flow induced by steady vibrations of plates on poiseuille flow between two parallel plates, under a constant pressure gradient, is then investigated. Calculations are made for longitudinal vibrations of the plates. 1. INTRODUCTION In
recent years much work has been published on the continuum theory of mixtures. In the first theory, given by Truesdell and Toupin [l], equations of mass, momentum, moment of momentum and energy balance are basic postulates for each ingredient in the mixture. The second theory, given by Green and Naghdi [Z], uses only an equation of energy balance and entropy production inequality as basic postulates. Coincidentally, a paper by A, C. Eringen and 9. D. Ingram has been published which presents a unified approach for the derivation of basic conservation equations for a chemically reacting continuum [3]. Mixtures of fluids which are viscous in the sense that their constitutive equations contain gradients of the velocity fields of the constituents have been considered in the context of the continuum theory by Green and Naghdi [4], Craine et al. [53, Mills [6J, Bowen [7], Mtiller 181, Adkins [9], R. M. Bowen and Garcia [lo], Eringen and Ingram [ll] have investigated mixtures of chemically reacting fluids in addition there is a considerable literature in which general two-phase flows are viewed from a different ~e~oint, much of this work being concerned with suspensions of solid particles in fluids. Lastly we mention the time dependent problems which has been investigated by Craine 112, 131. In the work described here a binary mixture of incompressible, Newtonian fluids is considered. In Section 2 we use the governing field equations of the continuum theory of mixtures and the notations given by R. E. Craine [12, 131. The appropriate constitutive theory developed by the authors mentioned above is given in Section 2 [4-61. In Sections 3 and 4 of present paper, the flow induced by longitudinal vibrations of two infinite parallel plates (plane poiseuille flow), under a constant pressure gradient, is investigated and an exact solution is obtained. 2. BASIC THEORY In this section we summarize the theory, for more details the reader should consult Craine [12]. A mixture of two continua, in motion relative to each other, is considered. At an arbitrary time t it is assumed that at each place z in the region of space occupied by the mixture, two particles are situated, one belonging to each constituent. This is a basic assumption of the continuum theory of mixtures and its validity requires that the mixture appears homogeneous when viewed on the scale of the applied disturbance. The velocities of the l- and, 2-constituents are denoted by Vfl) and V@ respectively, and we define DC*) a _=_+.j_,d
Dr
dt
&j
42-l)
where a takes the values 1, 2. The governing kinematic equations and the energy equation of a binary mixture in which the constituents have a common temperature T and do not interact chemically are
(2.2)
314
The quantities pm> Ua, CT(~),FCC), e_3,P and q are in turn rhe delE&.yy,the fnicrmal energy pei? unit mass, the partial stress, the body force per unit mass of the ath c~nst~t~e~t~ the difksive force vector, the heat supply per unit vo~~~~, and the heat flux vecteor[4]. The partial stresses satisfy O@l4” Cr_i$= 0
<2.6)
where suffixes enclosed in square brackets refer to components of the ant-isymmetric part of a second order tensor and circular brackets have a similar meaning in relation to the symmetric ~~a~s~~s-~~he~ ~n~~ua~~~ymay be written [S] , part. If s, denotes the partial entropy then t
The formulation of constitutive equations f~sr a mixture of Newtonian fluids is now considered, attention being restricted to the case in which both constituents are incompsressible. At some time tOr SI and S, are separated and have constant densities plOl pzo respectiveiy. The mixture is formed by adding the element: of ~7~aa _XA, of mass dq and volume d& to ehe element of S, at YA, of mass $mz and volume dV,. There is assumed to be no change of total volume when the mixt~~ Is formed, so that a~ element of the mixture at xi of mass dmr t dmz oceq&s a volume de/, + d&2; Such a mixture was considered by Mills [6]. We now write
When S, and Sz are separated we have
Using (2.8) and (2.10) it follows that
For the densities of SIs,,,!&in the rni~t~r~ we write
So that the incompressibility
condition (2.11) now becomes
where y is a ~orn~os~t~~~ factor. The density p of the mixture is given by P = P1-k Pz and we have
The influence of vibrating plates on poiseuille flow of a binary mixture
315
(2.16)
(pzo - P)d~’ + (P - PlO)d~’ - ~iUj= 0
the quantities 5 and d(“) being defined by Ei
=
dZi”)
P,i
=
(2.17)
“&)
partial Helimholtz free energies A, are introduced,
For convenience
defined by the relations
A a = U, - TS,
(2.18)
As the basic constitutive postulate we assume that the quantities A,, S,, co, q and a(“) are arbitrary functions of p and T but depend linearly on 5, b, a, d’“) and I’, where J.j$;’ = #
bi = T,i,
JYik= j.@ - wp
Id+
a.I = vi(l)- $2)
(2.19)
Adopting an analogous method to that outlined in [12], it is possible to deduce from the Clausius-Duhem inequality that the constitutive equations have the reduced forms A, = A&
‘A =A@,
7%
(P$++
-(Pzo-P)
Pl’
T),
Oi 4i =
= -
s, = UP,
=
__dA dT
(2.20) (2.21)
P2=(P-P1o)(P2$$-")
aai - Agi + ybi
(2.22)
T(koai + k,bi)
(2.23)
Us’ = (-Pl+
n,d~)+ a3d~‘)6, + 2pldik
Us’ = (-Pz +
n,d~) + &df))6,
where the coefficients satisfy the inequalities
s
T),
+ Z!/igdik+ n,rik
(2.24)
+ 2/44dik+ 2cczdik - a5rik
(2.25)
E, ill, . . . , d5, y, ko, k3, pl, . . . , y, are functions
of p and T, and
Q 2 0,
k3 3 0,
,420,
,&~001
(Y,+!h)2s‘hh
(2.26)
2 ;1,+-/L,ao, 3
(2.27) (2.28)
3. FLOW
INDUCED
BY STEADY LONGITUDINAL TWO PARALLEL PLATES
VIBRATIONS
OF
In this section we consider a binary mixture of incompressible Newtonian fluids, to be contained between infinite parallel rigid plates normal to the y-axis of a rectangular Cartesian coordinate system X and to flow under the action of constant pressure gradients PI, P2 in the direction of the x-axis. Simultaneously the two plates are subjected to a sinusoidal vibration in their own planes, of equal amplitudes with angular frequency w, parallel to the x-direction. For convenience set x1 = x, x2 = y, x3 = z, and let the plates be located at y = T/z. We assume that body forces are absent, and that the temperature remains constant throughout the mixture (for discussion of this latter assumption see Craine [12]). Since the temperature may now be dropped from the list of variables in the problem, all the coefficients pl etc. appearing in the constitutive equations (2.24) and (2.25) depend only on the density of the mixture p. We assume a flow field of the form v(‘) = (u,(y,t),
0, 0) ,
vC2’ = (%(Y,
t), 0, 0)
P = P(Y,
t)
(3.1)
in which u,(y, t) and v,(y, t) are u,(y, t) = z+,(y) + L sin cot+ M cos ot
(3.2)
v,(y, t) = v,(y)
(3.3)
P = P(Y,
+ N sin ot + G cos wt t>
(3.4)
The iniluence of vibrating plates on poiseuille flow of a binary mixture
we have p # plo, p # pzo and d2(pA)/dp2 # 0. It must the mixture p is a constant. Since p has been proved have been assumed, all the coefficients in eqns (3.17) From (3.17), (3.18), (3.19) and (3.20) it follows constant A is
A=
_-=
a& dx
have dpldy = 0, from here the density of to be constant and isothermal conditions and (3.18) are constant. that A and fI(y) are constant, where
34 dx --=A
Pzo-P
317
P-P10
di2
(3.24)
3.x
On substituting eqns (3.2) and (3.3) into eqns (3.17) and (3.18) we find that (3.25) (3.26) (3.27) (3.28) (3.29) (3.30)
introducing the dimensionless
quantities,
(3.31)
where ,U is a typical viscosity dimensionless forms,
in the problem,
eqns (3.25)-(3.30)
take on the following
Replacing (3.33) by the sum of (3.32) and (3.33), the equations can be written
with the use of the boundary (3.26) gives,
conditions
G&F@ = 0, “Ls,(Th) = 0 on f = T&, integrating
Substituition of (3.40) into eqn (3.38) gives a linear differential linear equation, for solution &, and DOwe obtain
equation.
eqn
With the soMow
oi
where,
(3.43)
The last four of eqns (3.34)-(3.37) the notation
may be rewritten
Ri,=b-dti,
as two complex equations by introducing
R,=A+iC
(3.44)
Thus (3.45)
319
The influence of vibrating plates on poiseuille flow of a binary mixture
Replacing (3.46) by the sum of (3.45) and (3.46) the governing equations can be written. (3.47) [&
+ ,i&)D2 - i&n]&
where D denotes the differential operator
+ [(p2 + ,i&)D2 - ip2n]Z?, = 0
(3.48)
d/d7 and put
WI = [(fi2 + ,i&)D2 - i&n]Y,
Z&= [-(PI + ,&)D2 + @&I’
(3.49)
the eqn (3.48) is satisfied identically and (3.47) becomes (D4 + P1D2 + /?$I’ = 0
(3.50)
where
-(ih+P2+P3+rii4)-in
P1p2+~2~1+~~5) (
ie2
=el+
A=
(3.51)
PO
p =
-P&n2
+ in =e3+
2
ie,
Bo
seeking a solution of (3.50) of the form W = B exp AJ leads to a bi-quadratic with roots
equation for il
?L;= &+ itI = ; [ - ,& + (f?f - @$“] (3.52) a; = s, + it2= ; [--PI - (pf - 4p2)‘“] S, and ta being real constants. The complex numbers S, + it, have unique square roots with negative real part, say -1, + im,, where m, = -t,/21,
1, = 2-‘“[(S2, + @‘2 + S,]l”,
(3.53)
the solution for Y is Y = B exp(-1,
+ im,)jj + C exp(-l2
+ E exp[-(-1,
+ im2)jj+ D exp[-(-1,
+ im,)jj]
+ im2)y]
(3.54)
where B, C, D and E are complex constants, and it follows from (3.49) that l?, = B[(P2 + P&f - iP2n] * exp[(-1, + im,)j] + C[(p2 + @,)A,”- io2n] - exp[(-1, + D[(jQ2+ j13)A.f- ip2n] * exp[-(-lI + im&] + E[(j12 + ,~I,)dz- iP2n] - exp[-(-1, & = B[-(PI + ,k,)kf + i&n]exp[(-ll
+ im&] + C(-(PI + ,&)A$+ ip2n]exp[(-1,
f D[-(j& +&&If + ip,n]exp[-(-lI + E[-(pl
+ im,)y] + im2)9] -I-im2)Y]
+ im&j] + a4)Ag+ iP2n]exp[ - (-l2 + im2)y] (3.55)
Applying the boundary conditions (3.5) and (3.6) to eqns (3.55) gives jil = sin nZ,
R, = 1,
z= 1,
M=O,
on
J=FR
0, = sin nZ,
a2=1,
N=l,
G=O,
on
jj=T&
and
(3.56)
in (C - E)[exp(-
12
where
we must lame B = D,
Thus,
setting
It fQilOW§ frsm (3.55)
wkere
The influence of vibrating plates on poiseuilie flow of a binary mixture
i&,12”)
qm
iqf2l-o
321
q (gl
Fig. 1. Variation of fluid velocities. a = S/6, pz = 1.4, n = 1. sl,
s2, tl
and t2 being real constants, we calculate from eqns (3.51) and (3.52) that
[~(rl+r~)]liz]
,,=&+
S2=~[-e,-[~(r1+r,)lin]
r
IJZ
t*=;i-e2-t- &-r*+KJ) II t2=;[-ez
- [;
(-r* +
(3.64)
rq]
where rl =
ef - ef - 4e3
r2 =
e, . e2
-
2~
(3.65)
r3= (T?+ 4ry
G
’
,a 2
u,
(2lT)
“l,rl,
i7( 11, 2
Fig. 2. Variation of fluid velocities. p1 = 5/6, pz = 5, n = 1.
322
Fig. 3. Variation of fluid veiocities. p1 = 516, p2 = 1.4, PI= 100.
and it follows that
The numerical results are divided into two groups, The first group gives the variation of the velocities of the fluid with frequency at a fixed dimensionless time (f = mz/Zn, ~PX = I, 2, 3, 4; pz= 1, 100) and at a fixed distance from the plates. In the second set a fixed frequency of oscillation of the plates is assumed and the variation of the velocities with distance from the plates are found, again at a fixed dimensionless time. We use some numericai results which are E. Craine [I?]. With a varied selection of values for given to complete the discussion by
Fig. 4. Variation of fluid velocities. pI = 5/6, p2 = 5, n = 180.
The influence of vibrating plates on poiseuille flow of a binary mixture
323
these quantities computer programmes were run. Though in all the curve shown certain parameters are fixed (aI = 1, p3 + +1/2, ,& = +1/2, x5 = -l/2, PI = -1, pz = -1). Variation of fluid velocities is investigated for dimensionless plate vibration frequency ranging 12= 1 (Figs 1 and 2) and y1= 100 (Figs 3 and 4), taking viscosity coefficient ranging I.$= 1, 4 and 5, nS= mn/2 m(l,2, 3,4). Comparison of Figs 1 and 2 show that flow characteristics do not show a distinctive variation in Fig. 1 compared to Fig. 2. Comparing Figs 3 and 4, one concludes that figure exhibits nearly single fluid characteristics compared to Fig. 4.
REFERENCES [l] C. TRUESDELL and R. TOUPIN, The Classical Field Theories, Hundbuch der Physik, (Ed. S. FLOGGE), Band 111/l. Springer, Berlin (1960). [2] A. E. GREEN and P. M. NAGHDI, ht. J. EngrzgSci. 3,231 (1965). [3] A. C. ERINGEN and J. D. INGRAM, Znt. J. Engng Sci. 3, 197 (1965). [4] A. E. GREEN and P. M. NAGHDT, Q. J. Mech. uppl. Math. 22, 427 (1969). [5] R. E. CRAINE, A. E. GREEN and P. M. NAGHDI, Q. J. Mech. appl. Math. 23, 171 (1970). [6] N. MILLS, Znt. J. Engng Sci. 4, 97 (1966). [7] R. M: BOWEN, J. Chem. Phys. 49, 162.5 (1968). [S] I. MULLER, Arch. Rntl. Mech. Anal. 28, 1 (1968). [9] R. J. ATKIN, Q. J. Mech. appl. Math. 21, 171 (1968). [lo] R. M. BOWEN and D. J. GARCIA, Znt. J. Engng Sci. 8,63 (1970). [ll] J. D. INGRAM and A. C. ERINGEN, ht. J. Engng Sci. 5, 289 (1967). [12] R. E. CRAINE, Oscillations of a plate in a Binary mixture of incompressible Newtonian fluids. Znt. J. Engng. Sci. 9, 1177 (1971). [13] t9Fi)CRAINE
An oscillating plate in a binary mixture of hemihedral fluids. J. Appl. Math. Phys. (ZAMP)
[14] J. Y. IkAZAKIA and R. S. RIVLIN, Rheol. Acta 17,210-226 (1978). [15] J. Y. KAZAKIA and R. S. RIVLIN, Rheol. Acta 18, 224-255 (1979). (Received 21 April 1987)
24