 Open Access
 Total Downloads : 113
 Authors : Deepak Agarwal, Reshu Agarwal
 Paper ID : IJERTV5IS090366
 Volume & Issue : Volume 05, Issue 09 (September 2016)
 DOI : http://dx.doi.org/10.17577/IJERTV5IS090366
 Published (First Online): 20092016
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Flow of a NonNewtonian SecondOrder Fluid Over an Enclosed Torsionally Oscillating Disc in the Presence of the Magnetic Field
Deepak Agarwal* & Reshu Agarwal**
*Grd Girls Degree College, Dehradun
**University Of Petroleum & Energy Studies, Dehradun
Abstract – The problem of the flow of an incompressible nonNewtonian secondorder fluid over an enclosed torsionally oscillating discs in the presence of the magnetic field has been discussed. The obtained differential equations are highly nonlinear and contain upto fifth order derivatives of the flow functions. Hence exact or numerical solutions of the differential equations are not possible subject to the given natural boundary conditions, therefore the regular perturbation technique is applied. The flow functions H, G, L and M are expanded in the powers of the amplitude (taken small) of the oscillations. The behaviour of the radial, transverse and axial velocities at different values of Reynolds number, phase difference, magnetic field and secondorder parameters has been studied and shown graphically. The results obtained are compared with those for the infinite torsionally oscillating discs by taking the Reynolds number of outflow Rm and circulatory flow RL equal to zero. The transverse shearing stress and moment on the lower and upper discs have also been obtained.
Key Words: Flow; SecondOrder Fluid;Enclosed Torsionally Oscillating Disc, Magnetic Field.
1 INTRODUCTION
The phenomenon of flow of the fluid over an enclosed torsionally oscillating disc (enclosed in a cylindrical casing) has important engineering applications. The most common practical application of it is the domestic washing machine and blower of curd etc.. Soo1) has considered first the problem of laminar flow over an enclosed rotating disc in case of Newtonian fluid. The torsional oscillations of Newtonian fluids have been discussed by Rosenblat2). He has also discussed the case when the Newtonian fluid is confined between two infinite torsionally oscillating discs3). Sharma & Gupta4) have considered a general case of a secondorder fluid between two infinite torsionally oscillating discs. Thereafter Sharma & Singp) extended the same problem for the case of porous discs subjected to uniform suction and injection. Hayat6) has considered non Newtonian flows over an oscillating plate with variable suction. Chawla7) has considered flow past of a torsionally oscillating plane Riley & Wybrow8) have considered the flow induced by the torsional oscillations of an elliptic cylinder. Bluckburn9) has considered a study of two dimensional flow past of an oscillating cylinder. Sadhna Kahre10) studied the steady flow between a rotating and porous stationary disc in the presence of transverse
magnetic field. B. B. Singh and Anil Kumar11) have considered the flow of a secondorder fluid due to the rotation of an infinite porous disc near a stationary parallel porous disc. Present paper is extended work of Reshu Agarwal13) who has considered Flow of a NonNewtonian SecondOrder Fluid over an Enclosed Torsionally Oscillating Disc.
Due to complexity of the differential equations and tedious calculations of the solutions, no one has tried to solve the most practical problems of enclosed torsionally oscillating discs so far. The authors have considered the present problem of flow of a nonNewtonian secondorder fluid over an enclosed torsionally oscillating disc in the presence of the magnetic field and calculated successfully the steady and unsteady part both of the flow functions. The flow functions are expanded in the powers of the amplitude (assumed to be small) of the oscillations of the disc. The nonNewtonian effects are exhibited through two dimensionless parameters 1(=n2/1) and 2(=n3/1), where 1, 2, 3 are coefficient of Newtonian viscosity, elasticoviscosity and crossviscosity. n being the uniform frequency of the oscillation. The variation of radial, transverse and axial velocities with elasticoviscous parameter 1, crossviscous parameter 2, Reynolds number R, magnetic field m1 at different phase difference is shown graphically.
2 FORMULATION OF THE PROBLEM
The constitutive equation of an incompressible second order fluid as suggested by Colemann and Noll12) can be written as:
ij = – pij + 21dij + 22eij + 43cij ————— ( 1 ) where
d ij = Â½ (ui,j + uj,i),
,i
,i
e ij = Â½(ai,j + aj,i) + um um,j,
j.
j.
c ij = dim dm ———– ( 2 )
p is the hydrostatic pressure, ij is the stresstensor, u i and a i are the velocity and acceleration vector.
The equation (1) together with the momentum equation for
no extraneous force
i,
i,
( ui/t + um ui,m) = tm m ———————–( 3 )
and the equation of continuity for incompressible fluid
ui
ui
,i = 0 ————————————— ( 4 )
where is the density of the fluid and comma (,) represents covariant differentiation, form the set of governing equations.
In the three dimensional cylindrical set of co ordinates (r, , z) the system consists of a finite disc of radius rs (coinciding with the plane z = 0) performing rotatory oscillations on the type rCos of small amplitude
, about the perpendicular axis r = 0 with angular velocity
in an incompressible secondorder fluid forming the part of a cylindrical casing or housing. The top of the casing (coinciding with the plane z = z0 rs) may be considered as a stationary disc (stator) placed parallel to and at a distance equal to gap length z0 from the oscillating disc. The symmetrical radial steady outflow has a small mass rate m of radial outflow ( m for net radial inflow). The inlet condition is taken as a simple radial source flow along z axis starting from radius r0. A constant magnetic field B0 is applied normal to the plane of the oscillating disc. The induced magnetic field is neglected.
The conditions on M on the boundaries are obtainable form the eq.(7) for m as follows:
M(1,) – M(0,) = 1, ————————— ( 9 )
which on choosing the discs as streamlines reduces to M(1,) = 1, M(0,) = 0
————————————————————— ( 10 )
Using eqs.(1) and expression (6)in equation (3) and neglecting the squares & higher powers of Rm/Rz (assumed small), we have the following equations in dimensionless form:
– (1/z0)(p/) = – nz0{H (Rm/Rz) (M/)}+ 2z0
(H2 2HHG2)
+2z0(Rm/Rz)(2HM/)
2z0(RL/Rz)(2LG/)+(1/z0){H (Rm/Rz)(M/)}(22/z0)[(n/2)
{( Rm/Rz)(M/) H}
+ 2(H2 HHiv) + (Rm/Rz)( 2/) (HM+HM+HM+ HMiv) ( RL/Rz) (22/)(LG + LG)](432/z0) {(Rm/Rz)(1/2) (HM +HM +HM)
( RL/Rz)(1/2)(2LG + LG) + (/4)(H2 G2 2HH)}+(B 2z /){H+ (R /R )(M/)} — — ( 11 )
0 0 m z
Assuming (u, v, w) as the velocity components along the cylindrical system of axes (r, , z)the boundary conditions of the problem are:
z = 0, u = 0, v = r ei(Real part), w = 0,
z = z0, u = 0, v = 0, w = 0, ——————–( 5 )
where the gap z0 is assumed small in comparison with the disc radius rs. The velocity components for the axisymmetric flow compatible with the continuity criterion can be taken as 4).
U = – H(,) + (Rm /Rz) M(,)/
V = G(,) + (R /R )L(,)/,
0 = nz0{G + (RL/Rz)(L/)}22z0) (HGHG)2z0 (Rm/Rz)(2MG/) 2z0(RL/Rz)(2HL/)+ (1/z0){G + (RL/Rz)(L/)}+ (22/z0)(n/2) {G+( RL/Rz)(L/)}+ (RL/Rz)( 2/) (HL +HL+ HL
2
2
+HL) +(2) (HG HG)+( Rm/Rz)(22/) (MG+MG)]+(232/z0){(HG HG)+( RL/Rz)(1/)(HL + HL + HL)+( Rm/Rz)(1/)(2MG + MG) – (B0 z0/){G+(RL/Rz)(L/)} —( 12 )
(1/z0)(p/) = 2nz0H+42z0HH 21H/z0 (22/z0){nH+ 222 (HH+GG)+2(22HH + 2HH) ( Rm/Rz)22 (HM+HM) + (RL/Rz )
22(LG + LG)}(232/z0) {2(HH + GG) +
L z 14HH (Rm/Rz) (HM + HM) + (RL/Rz)(LG +
W = 2H(,). ————————————— (6)
where U = u /z0, V = v /z0, W = w /z0, = r /z0, ,
are dimensionless quantities and H(,), G(,), L(,), M(,) are dimensionless function of the dimensionless
LG)} —— ( 13 )
where B0 and are intensity of the magnetic field and conductivity of the fluid considered. R (=nz02/1) is the Reynolds number, 1(=n2/1), 2(=n3/1) and (=/n) are
the dimensionless parameter, m2 = B 2z 2/ is the
variables = z/z0 and = nt. Rm(=m/2z 0 1), RL (=L/2z01) are dimensionless number to be called the
dimensionless magnetic field.
0 0 1
Reynolds number of net outflow and circulatory flow
respectively. Rz(=z02/1) be the flow Reynolds number. The small mass rate m of the radial outflow is represented
Differentiating (11) w.r.t. and (13 ) w.r.t. and then
eliminating 2p/. from the equation thus obtained. We get
nz0{H(Rm/Rz)M/}22z0 (HH+GG)
by
+(R
/R )(22z /) (HM+HM)(R /R )(22z /)
z0 m z 0
L z 0
m = 2 r u dz ————– ( 7 ) 0
Using expression (6), the boundary condition (5) transform for G, L & H into the following form:
G(0,) = Real(ei), G(1,) = 0,
L(0,) = 0, L(1,) = 0,
H(0,) =0, H(1,) = 0,
H(0,) = 0, H(1,) = 0.
————————— ( 8 )
(LG+LG)(1/z0){(Rm/Rz)(Miv/)Hiv} (22/z0)
[(n/2){(Rm/Rz)(Miv/) Hiv}2(2HH+HHiv+HHv+4GG)+ (Rm/Rz)( 2/) (2HM+
HivM+2HM+2HMiv+HMv)(RL/Rz)(22/)
(2LG+LG+LG)] (232/z0){(Rm/Rz)(1/)(HivM+2HM+ 2HM+
HMiv) (RL/Rz)(1/) (3LG+2LG+LG)(HHiv + 3GG+2HH)} + (B02z0/)
{H+(Rm/Rz)(M/)}= 0 ————————— ( 14 )
On equating the coefficients of and 1/ from the equation

& (14), we get the following equations:
G = RG + 2R(HGHG)1G21(HGHG)
22(HGHG) +m2G —————–(15) L = RL+2R(MG + HL)1L 21(HL+HL+HL+HL+2MG+
2MG)22(HL+HL+HL+2MG +MG) +m2L
———— ( 16 )
Hiv = RH+2R(HH+GG)1Hiv
21(HHiv+HHV+2HH+4GG) 22 (HHiv
+2HH+3GG) +m2H —————————– ( 17 ) Miv = RM+2R(HM+HMLGLG)1Miv 21(2HM+HivM +
2HM +2HMiv+HMv4LG2LG2LG)
22(HivM+2HM+
2HM +HMiv3LG2LGLG) +m2M —– ( 18 )

SOLUTION OF THE PROBLEM Assuming the relationship m2 = m12, equations (15) (18) becomes
G = RG + 2R(HGHG)1G21(HGHG)22(HG HG)
1
1
+ m 2G ——————(19)
L = RL+2R(MG + HL)1L 21(HL+HL+HL+HL+2MG+ 2MG)
1
1
22(HL+HL+HL+2MG +MG) +m 2L — ( 20 )
Hiv = RH+2R(HH+GG)1Hiv 21(HHiv+HHv+2HH+4GG) 22 (HHiv +2HH+3GG) +
1
1
m 2H — —————– ( 21 )
Using (23) and (32), the boundary conditions (8) & (10) for N = 0, 1 transforms to
G0s(0) =0, G0t(0) = 1, G1s(0) =0, G1t(0) = 0,
G0s(1) =0, G0t(1) = 0, G1s(1) =0, G1t(1) = 0,
H0s(0) =0, H0t(0) = 0, H1s(0) =0, H1t(0) = 0,
H0s(1) =0, H0t(1) = 0, H1s(1) =0, H1t(1) = 0,
H0s(0) =0, H0t(0) = 0, H1s(0) =0, H1t(0) = 0,
H0s(1) =0, H0t(1) = 0, H1s(1) =0, H1t(1) = 0,
L0s(0) =0, L0t(0) = 0, L1s(0) =0, L1t(0) = 0,
L0s(1) =0, L0t(1) = 0, L1s(1) =0, L1t(1) = 0,
M0s(0) =0, M0t(0) = 0, M1s(0) =0,
M1t(0) = 0,
M0s(1) =0, M0t(1) = 0, M1s(1) =0,
M1t(1) = 0,
M0s(0) =0, M0t(0) = 0, M1s(0) =0, M1t(0) =
0,
M0s(1) =1, M0t(1) = 0, M1s(1) =0, M1t(1) = 0.
—— ( 33 )
Applying (32) & (33) in eqs. (24) to (31), we get G0s() = G1s() = 0,
G0t() = [sinh {d(1)}]/sinh d,
[R{(1+1 1 1 [R{(1+1 1 12)1/2}/ {2(1+ 2)}]1/2 + i
2)1/2}/ {2(1+ 2)}]1/2 + i
where d = {iR/(1+i1)}1/2 = [R{1+(1+1 1 2)1/2 }/{2(1+ 2)}]1/2 = A+ iB,
G0(,) = Real{ei G0t()},
= [cos.{cosh{(2)A}.cosB – coshA.cos{(2)B}} sin
.{sinhA. sin {(2)B} sinh{(2)A}.sinB}] /(cospAcospB),
1
1
G1t() = [(m 2Sinh d)/{2d(1+i1) Sinpd}] [{m12 Cosh d(1)}/
{2d(1+i1) Sinh d}],
G1(,) = Real{ei G1t()},
= (m 2/2){( – ).Cos – ( – ).Sin },
Miv = RM+2R(HM+HMLGLG)1Miv21(2HM+HivM
1 17 19
18 20
+
2HM +2HMiv+HMv4LG2LG2LG)22(HivM+2HM+
1
1
2HM+HMiv3LG2LGLG) + m 2M ———— ( 22 )
Substituting the expressions
G(,) = N GN (,)
L(,) = N LN (,)
H(,) = N HN (,)
M(,) = N MN (,) —————————— ( 23 )
into (19) to (22) neglecting the terms with coefficient of 2 (assumed negligible small) and equating the terms independent of and coefficient of , we get the following equations:
G0 = R G0/ – 1 G0/ —– ( 24 )
2G — ( 25 )
2G — ( 25 )
G1 = R G1/2R(H0G0H0G0)1G1/21(H0G0H0G0) 22(H0G0 H0G0) + m1 0
L = R L / – L / —— ( 26 )
where, 1 = CosB.SinhA,
2 = SinB.CoshA,
3 = (Cos2B.CospA1)/2,
4 = (Sin2B.SinpA)/2,
5 = CosB(1).CoshA(1),
6 = SinB(1).SinhA(1),
7 = CosB.SinhA,
8 = SinB.CoshA,
9 = A3 – B4,
10 = B3 + A4,
11 = A7 – B8,
12 = B7 + A8,
13 = 9 – 110,
14 = 10 + 19,
15 = 11 – 112,
16 = 12 + 111,
14
14
17 = (113 + 214)/(132+ 2),
2 + 2),
2 + 2),
18 = (213 – 114)/(13 14
0 0 1 0
L = RL /2R(M G + H L ) L / 19 = (515 + 616)/(152 + 162)
1 1 0 0 0 0 1 1
2 (H L +H L +H L 20 = (615 – 516)/(152 + 162),
1 0 0 0 0 0 0
+H0L0+2M0G0+ 2M0G0)22(H0L0+H0L0+H0L0+2M0G0
+M0G0) + m1 0
+M0G0) + m1 0
2L ————( 27 )
H iv = RH / – H iv/ — ( 28)
G(,) = G0(,) + G1(,). L0s() = L0t() = L1s() = 0,
L1t() = {sinh d /sinh d}[{(A1A3)/2}(1/2d3 1/6d)+ A2/2d2]
+ cosh{d(1)} [{(A1 A3)/2}{/2d3 – 2/2d + 3/3d} + (A2/2d 2}] +
0 0 1 0
H iv = RH /+2R(H H +G G ) H iv/2 (H H iv sinh{d(1)}
1 1 0 0
0 0 1 1
1 0 0
+H0H0 0 0 0 0
+H0H0 0 0 0 0
v+2H H +4G G )
0 1
0 1
22(3G0G0+H0H iv+2H0H0) + m 2H0 ———— ( 29 )
M0 0 1 0
M0 0 1 0
iv = R M / – M iv/
——————– ( 30 )
[{(A1 A3)/2}{(2)/2d2}+ A2(2)/2d],where, A1 = 12R/{(1+i1)sinh d}, A2 = 24(1+2)d/{(1+i1)Sinh d},
A3= 6d2(41+22)/ {(1+i1)Sinh d}, L1(,) = Real{ei L1t()},
M iv = RM / + 2R(H M + H M L G L G ) –
1M1 – 21 0 0 0 0 0 0 0 0
1M1 – 21 0 0 0 0 0 0 0 0
1 1 0 0
0 0 0
0 0 0
= (N7+N9N5).cos – (N8+N10N6).sin ,
v+ 2H M iv)2
v+ 2HM iv)2
(2H0M0+H0 0 0 0 0 0 0 0 0 0
(2H0M0+H0 0 0 0 0 0 0 0 0 0
iv/ (2H M +H ivM + 2H M – 4L G – 2L0G02L0G0+H0M0 0 0 2
ivM +2H M – 3L G 2L G – L G
+ H0M0 1 0
+ H0M0 1 0
iv) + m 2M —————— ( 31 )
where, N1 = [{6R(121+62)(A2B2)}(cosB.sinhA1sinB.coshA) 2AB(121+62) (1.cosB.sinhA+sinB.coshA)]
/[(cosB.sinhA sinB.coshA)2 + ( .cosB.sinhA+sinB.coshA)2],
1 1
Taking Gn(,)= Gns()+ eiGnt() N = [{6R(12 +6 )(A2B2)}( cosB.sinhA+sinB.coshA)
2 1 2 1
Ln(,)= Lns()+ eiLnt() 2AB(12 +6 ) (cosB.sinhA – .sinB.coshA)]
1 2 1
Hn(,)= Hns()+ e2iHnt()
Mn(,)= Mns()+ e2iMnt()—————— ( 32 )
/[(cosB.sinhA1sinB.coshA) 2 + (1.cosB.sinhA+sinB.coshA)2],
N3 = [24A(1+2)(cosB.sinhA1sinB.coshA)+24B(1+2) (1.cosB.sinhA+sinB.coshA)] /[(cosB.sinhA1sinB.coshA) 2 + (1.cosB.sinhA+sinB.coshA)2],
N4 = [24B(1+2)(cosB.sinhA1sinB.coshA)24A(1+2) (1.cosB.sinhA+sinB.coshA)] /[(cosB.sinhA1sinB.coshA)2 + (1.cosB.sinhA+sinB.coshA)2],
I = (cosB.sinhA.cosB.sinhA +sinB.coshA.sinB.coshA)/
W7 = cos2B.sinpA, W8 = sin2B.cospA,
W9 = 2A(cos2B.cospA1) 2Bsin2B.sinpA,
W10 = 2B(cos2B.cospA1)+ 2Asin2B.sinpA,
W11 = 1+CeCcosD,
1 2 2 W12 = DeCsinD,
(sin B+sinh A),
I2 = (sinB.coshA.cosB.sinhAcosB.sinhA.sinB.coshA)/ (sin2B+sinpA),
W13 = W3(W5W7)W4(W6W8) W9W11+W10W12,
W = W (W W )+W (W W )
I3 = [3(A2+B2)( A33AB2)A{(A33AB2)2+ (3A2BB3)2}] / [6(A2+B2)
14 4 5 7
3 6 8
{(A33AB2)2+ (3A2BB3)2}],
I4 = [B(A33AB2)2 +(3A2BB3)2} 3(A2+B2) (3A2BB3) ] / [6(A2+B2)
{(A33AB2)2+ (3A2BB3)2}],
I = {(A2B2)N +2N AB}/ [2{( A2B2)2 + 4A2B2}],
W10W11W9W12, W15 = W1W13W2W14, W16 = W2W13+W1W14,
W17 = 4A24B2C2+D2,
5 3 4
2 2
2 2 2 2 2
W18 = 8AB2CD,
I6 = {(A B )N42N3AB}/ [2{( A B ) + 4A B }],
I7 = [3( A2+B2)( A33AB2)(3A22A3){(A33AB2)2+ (3A2BB3)2}] / [6( A2+B2){ (A33AB2)2 +(3AB2B3)2}],
I8 = [3( A2+B2)(3A2BB3)+(3B22B3){(A33AB2)2+ (3A2BB3)2}] / [6( A2+B2){ (A33AB2)2 +(3AB2B3)2}],
I = [{N (A2B2)+2ABN }] / [2{( A2B2) 2+4A2B2}],
W19 = W17(8A28B2)16ABW18, W20 = W18(8A28B2)+16ABW17, W21 = 4eC{(C+2).cosD + D.sinD}
+eC{(C2).cosDDsinD}, W22 = eC{(C2).sinD + D.cosD}
C
9 3 4
e {DcosD(C+2).sinD},
I10 = [{N4(A2B2)2ABN3}] / [2{( A2B2) 2+4A2B2}],
W = W W W W ,
I = (A2B2)/[(A2B2)2+4A2B2],
23 19 21 20 22
11
I12 = 2AB/[(A2B2)2+4A2B2],
W24 = W20W21+W19W22,
Q = (W W +W W )/(W 2+W 2),
2 2 3
15 23
16 24
23 24
I13 = (AN3+BN4)/( A +B ), I14 = (AN4BN3)/( A2+B2),
Q4 = (W16W23W15W24)/(W232+W242),
X = Q (1eCcosD)Q eCsinD,
3 3 4
N5 = (N1I1I3+N2I1I4+I1I5N1I2I4+N2I2I3I2I6),
N6 = (N1I2I3+N2I2I4+I2I5+N1I1I4N2I1I3+ I1I6),
N7 = (N117+N2I8+I9)cosB(1).coshA(1)(N1I8N2I7+I10) sinB(1
).sinhA(1),
N8 = (N1I8N2I7+I10)cosB(1).coshA(1)+(N1I7+N2I8+I9) sinB(1
).sinhA(1),
N9 = [(N1I11N2I12+I13)cosB(1) .sinhA(1)+(N2I11+N1I12I14) sinB(1
).coshA(1)](2/2),
N10 = [(N1I11N2I12+I13)sinB(1).coshA(1)(N2I11+N1I12I14)
cosB(1).sinhA(1)](2/2),
L(,) = L0(,) + L1(,) = L1(,).
X4 = Q4(1eCcosD)+Q3eCsinD, X5 = (1eCcosD),
X6 = eCsinD,
6
6
5 6
5 6
X7 = (X3X5+X4X6)/(X52+X 2); X8 = (X4X5X3X6)/(X 2+X 2); X9 = W1W9W2W10,
2),
2),
X10 = W2W9+W1W10, X11 = W19X5W20X6, X12 = W20X5+W19X6,
12
12
X13 = (X9X11+X10X12)/(X112+X12 X14 = (X10X11X9X12)/(X112+X 2),
Y3 = {C(Q3Q5)+D(Q4Q6)}/(C2+D2),
Y = {C(Q Q )D(Q Q )}/(C2+D2),
4 4 6 3 5
H0s() = H0t() = H1t() = 0,
H1t() = c1ef / f 2 + c2ef / f 2 + A4sinh 2d(1)/{8d2(4d2f2)} + c3 + c4, Where,
c1 = [{c2(1ef )/(1ef )} + A4f(cosh 2d1}]/{4d(4d2f2)(1ef ),
c2 = A4f [f(1 ef )(2d.cospdsinh 2d) 2d(1+fef )(cosh 2d 1)]/{8d2 (4d2 f2) (42ef – 2ef – fef + fef ),
c3 = (1/f)[c2 c1 + {fA4cosh 2d /{4d(4d2f2)}],
c4 = (1/f 2)[c1 + c2 + {A4f 2sinh 2d/ {8d2(4d2f2)}], A4 = [(81 + 62)d3 2Rd]/{(1+2i1) sinpd},
= [2R{21+(1+41 1
= [2R{21+(1+41 1
f = (2iR/(1+2i1}1/2
2) 1/2}/{2(1+4 2)}]1/2
Y5 = Q1cos2B.cospAQ2sin2B.sinpA, Y6 = Q2cos2B.cospA+Q1sin2B.sinpA, Y7 = 4AW174BW18,
Y8 = 4BW17+4AW18,
8
8
8
8
Y9 = (Y5Y7+Y6Y8)/(Y72+Y 2), Y10 = (Y6Y7Y5Y8)/(Y72+Y 2), Q5 = X7+X13,
Q6 = X8+X14, Q7 = Y3+Y9, Q8 = Y4+Y10,
Y11 = [(Q3+Q5)(C2D2)+2CD(Q4+Q6])/ [(C2D2)2+4C2D2],
Y = [(Q +Q )(C2D2)2CD(Q +Q ])/
2 2 Â½
12 4 6 3 5
+ [2R{(1+41 )21}/ {2(1+41 )}]
= C + iD,
H1(,) = Real{e2i H1t()},
= (J1+J3+J5+J7 ).cos 2 – (J2+ J4+J6+J8).sin 2,
where, X1 = (81+62)(A33AB2)2RA,
X2 = (81+62)(3A2BB3)2RB,
Y1 = (cos2B.cospA121sin2B
.sinpA)/2.
2+Y 2),
2+Y 2),
Y2 = {21(cos2B.cospA1)+ sin2B.sinpA}/2.
Q1 = (X1Y1+X2Y2)/(Y1 2
1 2
1 2
Q2 = (X2Y1X1Y2)/(Y 2+Y 2), W1 = Q1CQ2D,
W2 = Q2C+Q1D,
W3 = C(1eCcosD)+DeCsinD, W4 = D(1eCcosD)CeCsinD, W5 = 2Acos2B.cospA
2Bsin2B.sinpA, W6 = 2Bcos2B.cospA+
2Asin2B.sinpA,
[(C2D2)2+4C2D2],19
19
Y13 = [W19(Q1W7Q2W8)+W20 (Q2W7+Q1W8)] /(W 2+W202),
Y14 = [W20(Q1W7Q2W8)+W19 (Q2W7+Q1W8)] /(W192+W202),
Q9 = (Y11+Y13),
Q10 = (Y12+Y14),
O1 = eC(Q5cosDQ6sinD), O2 = eC(Q6cosD+Q5sinD), O3 = eC(Q3cosD+Q4sinD), O4 = eC(Q4cosDQ3sinD), J1 = [O1(C2D2)+2CDO2]/[(C2
D2)2+4C2D2],
J2 = [O2(C2D2)2CDO1]/[(C2D2) 2+ 4C2D2],
J3 = [O3(C2D2)+2CDO4]/[(C2D2)2+ 4C2D2],
J4 = [O4(C2D2)2CDO3]/[(C2D2)2+ 4C2D2],
J5 = Q7+Q9, J6 = Q8+Q10,
Y15 = cos{2B(1)}.sinh{2A(1)},
Y16 = sin{2B(1)}.cosh{2A(1)}, J7 = [W19(Q1Y15Q2Y16)+W20(Q2Y15+
19 20
19 20
Q1Y16)]/ {W 2+W 2},
2+W 2},
2+W 2},
J8 = [W20(Q1Y15Q2Y16)+W19(Q2Y15+ Q1Y16)]/{W19 20
H(,) = H0(,) + H1(,) = H1(,).
.M0(,) = M0s() = – 23 + 32, M0t() = M1t() = 0,
1
1
M1(,) = M1s() = – (m 2/20)(2554+432),
M(,) = M0(,) + M1(,).

RESULTS AND CONCLUSION
For Rm = RL = 0, the results are in good agreement with those obtained by Sharma and Gupta4) (for s1 = s2 = 1 and B0 = 0) and Sharma and Singp) (for s1= s2=1 and suction parameter = 0).
The variation of the radial velocity with for different values of elasticoviscous parameter 1 = 0.5, – 0.4, 0.3; when crossviscous parameter 2 = 2, m1 = 2, R =
0.5, Rm = 0.05, RL = 0.049, Rz = 2, = 0.02, = 5 at phase
difference = 2/3, /3 and 0 is shown in fig(1), fig(4) and fig(7) respectively. In all these figures, the behaviour of the radial velocity is being represented through the approximate parabolic curve with vertex upward. It is evident from fig (1) that the radial velocity decreases with an increase in 1 before the common point of intersection which is just before the line = 0.5 while it increases with an increase in 1 after this common point of intersection and from fig (4) & fig (7), the radial velocity increases with an increase in 1 before the common point of intersection and decreases with an increase in 1 after this common point of intersection. The value of the radial velocity is same approximately in the middle of the gaplength for all the values of elasticoviscous and phase difference parameters. In fig (1) the point of maxima of radial velocity is little beyond the half of the gaplength whenever in fig(4) it is in the middle of the gaplength and in fig(7) it is being shifted a little towardsthe oscillating disc from the middle of the gaplength.
The variation of the transverse velocity with for different values of elasticoviscosity parameter 1 = 0.5,
0.4, 0.3; when crossviscous parameter 2 = 2, m1 = 2, R
= 0.5, Rm = 0.05, RL = 0.049, Rz = 2, = 0.02, = 5 at
phase difference = 2/3, /3 and 0 is shown in fig(2), fig(5) and fig(8) respectively. In fig (2), the transverse velocity increases linearly and in fig(5) & fig(8), it decreases linearly throughout the gaplength. It is also evident from fig (2) and fig (5) that the transverse velocity increases with increase in 1 and from fig(8), the transverse velocity decreases with increase in 1 throughout the gap length.
The variation of the axial velocity with for different values of elasticoviscous parameter 1 = 0.5, – 0.4, 0.3; when crossviscous parameter 2 = 2, m1 = 2, R =
0.5, Rm = 0.05, RL = 0.049, Rz = 2, = 0.02, = 5 at phase
difference = 2/3, /3 and 0 is shown in fig(3), fig(6) and fig(9) respectively. In fig(3), the axial velocity increases in the first half and then decreases in the second half of the gaplength and forming the bell shape curve with point of maxima at = 0.5. It is also evident that the axial velocity
increases with increase in 1 throughout the gaplength. In fig (6) the behaviour of the axial velocity is being represented through the bell shape curve with point of maxima at = 0.5 for 1 = 0.5 & 0.4 whenever point of minima at = 0.5 for 1 = 0.3. It is also evident that the axial velocity decreases with increase in 1 and values of the axial velocity remains negative throughout the gap length for 1 = 0.3. In fig(9), the behaviour of the axial velocity is being represented through the bell shape curve with of minima at = 0.5. It is also evident that the axial velocity decreases with an increase in 1 and remains negative throughout the gaplength.
The variation of the radial velocity with for different values of crossviscous parameter 2 = 30, 20, 6; when elasticoviscous parameter 1 = 2, m1 = 2, R = 0.5,Rm = 0.05, RL = 0.049, Rz = 2, = 0.02, = 5 at phase
difference = 2/3, /3 and 0 is shown in fig(10), fig(13) and fig(16) respectively. It is evident from fig (10) and fig(16) that the radial velocity increases with increase in 2 before the common point of intersection and it decreases with an increase in 2 after the common point of intersection while in fig (13), the radial velocity decreases with increase in 2 before the common point of intersection (lying approximately the middle of the gaplength) and increases with an increase in 2 after the common point of intersection. The value of the radial velocity is same approximately in the middle of the gaplength for all the values of crossviscous and phase difference parameters. In fig (10) the point of maxima of the radial velocity is approximately in the middle of the gaplength whenever in fig(13) it is little beyond the half of the gaplength and in fig(16) it is being shifted a little towards the oscillating disc from the middle of the gaplength.
The variation of the transverse velocity with for different values of crossviscous parameter 2 = 30, 20, 6; when crossviscous parameter 1 = 2, m1 = 2, R = 0.5, Rm
= 0.05, RL = 0.049, Rz = 2, = 0.02, = 5 at phase
difference = 2/3, /3 and 0 is shown in fig(11), fig(14) and fig(17) respectively. In fig(11), the transverse velocity increases linearly and in fig(14) & fig(17), it decreases linearly throughout the gaplength. It is also evident from fig(11) & fig(14) that the transverse velocity slightly increases with increase in 2 in the first half and decreases with increase in 2 in the second half of the gaplength whenever in fig(17) the transverse velocity slightly decreases with increase in 2 in the first half and increases with increase in 2 in the second half of the gaplength. Because of their slightly increment or decrement, these figures are being overlapped.
The variation of the axial velocity with at different values of crossviscous parameter 2 = 30, 20, 6; when elasticoviscous parameter 1 = 2, m1 = 2, R = 0.5, Rm = 0.05, RL = 0.049, Rz = 2, = 0.02, = 5, phase
difference = 2/3, /3 and 0 is shown in fig(12), fig(15) and fig(18) respectively. In above three figures, axial velocity increases in the first half and then decreases in the second half throughout the gaplength. It is also evident from figure (12), the axial velocity decreases with increase
in 2 and from fig (15) & fig (18) the axial velocity increases with increase in 2. the behaviour of the axial velocity is being represented through the normal curve symmetric about = 0.5 line.
The variation of the radial velocity, transverse velocity and axial velocity with at different Reynolds number R = 1, 0.5, 0.1; when crossviscous parameter 2 = 10, elasticoviscous parameter 1 = 2, m1 = 2, Rm = 0.05, RL = 0.049, Rz = 2, = 0.02, = 5, phase difference =
/3 is shown in fig(19), fig(20) and fig(21) respectively. It is evident from fig (19) that the radial velocity decreases with an increase in Reynolds number R before the common point of intersection while increases with an increase in Reynolds number R after the common point of intersection. In fig (20), transverse velocity decreases linearly throughout the gaplength. It is also evident from this
0.009
0.008
0.007
0.006
0.005
0.004
0.003
0.002
0.001
0
At Xaxis At Yaxis U
1 = 0.5 – – – – –
1 = 0.4
1 = 0.3 _ _ _ _
figure that transverse velocity increases with increase in Reynolds number R. In fig (21), the axial velocity increases in the first half and decreases in the second half of the gap length. It is also evident from this figure that axial velocity increases with increase in Reynolds number R.
The variation of the radial velocity and transverse velocity with at different magnetic field m1 = 1, 30, 40; when crossviscous parameter 2 = 10, elasticoviscous parameter 1 = 2, Reynolds number R = 0.5, Rm = 0.05, RL
= 0.049, Rz = 2, = 0.02, = 5, phase difference = /3 is shown in fig(22) and fig(23) respectively. It is seen from this figure that the radial velocity increases with increase in m1 upto = 0.28, then it decreases with increase in m1 upto
= 0.75 and then it increases with increase in m1 upto = 1.
0 0.5 1 1.5
Fig(1) variation of radial velocity U with for different elasticoviscous parameter 1 at = 2/3
0
0 0.5 1 1.5
0.5
1
In fig (23), the transverse velocity decreases linearly for m1 = 1, it increases upto = 0.28 and decreases thereafter for m1 = 30, 40. it is also seen from this figure that transverse velocity increases with increase in m1.
There is no magnetic field term in axial velocity so variation of axial velocity at different m1 is not possible.
The transverse shearing stress on the lower and upper discs respectively is obtained as ;
(z) z = 0 = 1[(/z 0){1/(cospA cos2B)}{(A sinpA+ Bsin2B).cos +
(A sinpB – Bsin2A).sin} + (RL/Rz) (1/z 0) {N7+N9 N5}=0. cos
– (N8+N10N6) =0.sin]
1.5
2
2.5
3
At Xaxis At Yaxis V
1 = 0.5 – – – – –
1 = 0.4
1 = 0.3 _ _ _ _
and
(z) z = z0 = 1[(2/z 0){1/(cospA cos2B)}{(A sinhA.cosB+ BsinB coshA).cos + (A coshA.sinB BsinhA cosB).sin} + (RL/Rz)(1/z 0) {N7+N9
N5}=1. cos – (N8+N10N6) =1.sin].
Fig(2) variation of transverse velocity V with for different elastico viscous parameter 1 at = 2/3
0.012
0.01
At Xaxis At Yaxis W
1 = 0.5 – – – – –
3
2.5
At Xaxis
At Yaxis V
0.008 1
2
= 0.4
1 = 0.5 – – – – –
1 = 0.4
0.006
0.004
0.002
0
1 = 0.3
_ _ _ _
1.5
1
0
1 = 0.3 _ _ _ _
0 0.5 1 1.5
0 0.5 1 1.5
Fig(3) variation of axial velocity W with for different elasticoviscous parameter 1 at = 2/3
0.008
Fig(5) variation of transverse velocity V with for different elasticoviscous parameter 1 at = /3
0.0035
0.007
0.006
0.005
At Xaxis At Yaxis U
1 = 0.5 – – – – –
1 = 0.4
0.003
0.0025
0.002
0.0015
At Xaxis At Yaxis W
1 = 0.5 – – – – –
1 = 0.4
1 = 0.3 _ _ _ _
0.004
0.003
0.002
0.001
1 = 0.3 _ _ _ _
0.001
0.0005
0
0.0005
0.001
0 0.5 1 1.5
0
0 0.5 1 1.5
Fig(4) variation of radial velocity U with for different elasticoviscous parameter 1 at = /3
0.0015
Fig(6) variation of axial velocity W with for different elasticoviscous parameter 1 at = /3
0.009
0.008
0.007
0.006
0.005
At Xaxis At Yaxis U
1 = 0.5 – – – – –
1 = 0.4
0
0.001 0 0.5 1 1.5
0.002
0.003
0.004
0.004
0.003
0.002
1 = 0.3 _ _ _ _
0.005
0.006
0.007
At Xaxis At Yaxis W
1 = 0.5 – – – – –
1 = 0.4
0.001
0
0 0.5 1 1.5
0.008
0.009
0.01
1 = 0.3 _ _ _ _
Fig(7) variation of radial velocity U with for different elastico viscous parameter 1 at = 0
Fig(9) variation of axial velocity W with for different elasticoviscous parameter 1 at = 0
6
At Xaxis
5 At Yaxis V
4 1 = 0.5 – – – – –
1 = 0.4
0.009
0.008
0.007
0.006
3 1 = 0.3
2
1
_ _ _ _
0.005
0.004
0.003
0.002
At Xaxis At Yaxis U
2 = 30 – – – – –
2 = 20
0
0 0.5 1 1.5
0.001
0
2 = 6
_ _ _ _
Fig(8) variation of transverse velocity V with for different elastico viscous parameter 1 at = 0
0 0.5 1 1.5
Fig(10) variation of radial velocity U with for different crossviscous parameter 2 at = 2/3
0
0.5
0 0.5 1 1.5
0.012
0.01
0.008
At Xaxis At Yaxis U
2 = 30 – – – – –
2 = 20
1
1.5
2
2.5
At Xaxis At Yaxis V
2 = 30 – – – – –
2 = 20
0.006
0.004
0.002
2 = 6 _ _ _ _
2 = 6
3
_ _ _ _
0
0.002
0 0.2 0.4 0.6 0.8 1 1.2
Fig(11) variation of transverse velocity V with for different crossviscous paramet2er at = 2/3
Fig(13) variation of radial velocity U with for different crossviscous parameter 2 at = /3
0
0.001
0.002
0.003
0 0.5 1 1.5
3
2.5
2
At Xaxis At Yaxis V
2 = 30 – – – – –
2 = 20
0.004
0.005
0.006
0.007
At Xaxis At Yaxis W
2 = 30 – – – – –
2 = 20
1.5
1
2 = 6 _ _ _ _
0.008
0.009
2 = 6
_ _ _ _
0.5
0
Fig(12) variation of axial velocity W with for different crossviscous parameter 2 at = 2/3
0 0.2 0.4 0.6 0.8 1 1.2
Fig(14) variation of transverse velocity V with for different cross viscous parameter 2 at = /3
0.03
0.025
0.02
At Xaxis 6
At Yaxis W
5
2 = 30 – – – – –
2 = 20 4
2 = 6 _ _ _ _
At Xaxis At Yaxis V
2 = 30 – – – – –
2 = 20
0.015 3
2 = 6
_ _ _ _
0.01 2
0.005 1
0 0
0 0.5 1 1.50 0.5 1 1.5
Fig(15) variation of axial velocity W with for different crossviscous parameter 2 at = /3
Fig(17) variation of transverse velocity V with for different cross viscous parameter 2 at = 0
0
0.009
0.008
0.007
At Xaxis At Yaxis U
2 = 30 – – – – –
2 = 20
0 0.5 1 1.5
0.002
0.004
0.006
At Xaxis
0.006
0.005
0.004
0.003
0.002
2 = 6
_ _ _ _
0.008
0.01
0.012
0.014
0.016
At Yaxis W
2 = 30 – – – – –
2 = 20
0.001
0
0.018
2 = 6 _ _ _ _
0 0.5 1 1.5
Fig(16) variation of radial velocity U with for different crossviscous parameter 2 at = 0
Fig(18) variation of axial velocity W with for different crossviscous parameter 2 at = 0
0.009
0.008
0.007
0.006
0.005
0.004
0.003
0.002
0.001
0
At Xaxis At Yaxis U
R = 1 ——— R = 0.5
R=0.1 _ _ _ _
0.016
0.014
0.012
0.01
0.008
0.006
0.004
0.002
0
At Xaxis At Yaxis W
R = 1 ——— R = 0.5
R=0.1 _ _ _ _
0 0.5 1 1.5
0 0.5 1 1.5
Fig(21) variation of axial velocity W with for different Reynolds number R at = /3
Fig(19) variation of radial velocity U with for different Reynolds number R at = /3
3
0.009
0.008
0.007
2.5
2
1.5
1
At Xaxis At Yaxis V
R = 1 ——— R = 0.5
R=0.1 _ _ _ _
0.006
0.005
0.004
0.003
0.002
0.001
At Xaxis At Yaxis U
m1= 1 ———
m1= 30 m1= 40 _ _ _ _
0.5
0
0 0.5 1 1.5
Fig(20) variation of transverse velocity V with for different Reynolds number R at = /3
0
0 0.5 1 1.5
Fig(22) variation of radial velocity U with for different m1 at = /3
4.5
4
3.5
3
2.5
At Xaxis At Yaxis V
m1= 1 ———
m1= 30 m1= 40 _ _ _ _
2
1.5
1
0.5
0
0 0.5 1 1.5
Fig(23) variation of transverse velocity V with for different m1 at = /3

REFERENCES


Soo, S. L., Laminar flow over an enclosed rotating disc, Trans. A.
S. M. E. Jour. Basic Engg., 80 (1958) 287.

Rosenblat, S., Torsion oscillation of a plane in a viscous fluid, J, Fluid Mech., 6 (1959)206.

Rosenblat, S., Flow between torsionally oscillating discs, J, Fluid Mech., 8 (1960) 388.

Sharma, H. G. and D. S. Gupta, Flow of a nonNewtonian secondorder fluid between two infinite torsionally oscillating discs, Indian J. pure appl. Math., 14 (1983) 1091

Sharma, H.G. and K. R. Singh, flow of a secondorder fluid between two torsionally oscillating discs with uniform suction and injection,

Hayat, T (Pak QAZAM) (with Abbas, Q, (Pak QAZAM) Khan, Masood (Pak QAZAM) Siddiqui A. M. (IPASY), Non Newtonian flows over an oscillating plate with variable suction. Arch. Mech.55(2003)

Chawla, S.S., Flow past a torsionally oscillation plane, Z Anglew Math Phy 45 (1994)no6, 986990 (76 D)

N. Riley and M. F. Wybrow, The flow induced by torsional oscillation of an elliptic cylinder, J. Fluid Mech (1996) 279

H. M. Bluck Burn, A study of two dimensional flow past of an oscillating cylinder. J.Fluid Mech (1999) 255.

Khare Sadhna, Indian J. Pure & appl. Math. 8 (1977), 808815.

B. B. Singh & Anil Kumar, Indian J. Pure & appl. Math., 20(9), 931943,(1989)

Coleman, B.D. and W. Noll, An approximation theorem for functions with application in continuum mechanics, Arch. Ratl. Mech. Anal., 6(1960)355.

Reshu Agarwal, International Journal of Engineering Technology, Management and Applied Sciences, February 2016, Volume 4, Issue 2.