While the disks are initially rotating eccentrically, the unsteady flow caused by their oscillations in their own planes and in the opposite directions is studied. The analytical solutions to the problem are obtained for both small and large times, and thus the velocity field is determined for every value of time. The variations of all the parameters on the flow are scrutinized by means of the graphical representations. In particular, the effect of the ratio of the frequency of oscillation to the angular velocity of the disks is analyzed. The dependence of the oscillations in both x- and y-directions on the flow is examined. The influence of the Reynolds number is also investigated.
1. Introduction
The instrument called the orthogonal rheometer consisting of two parallel disks rotating with the same angular velocity about noncoincident axes was originally developed by Maxwell and Chartoff [1]. Abbott and Walters [2] obtained an exact solution for the flow of a Newtonian fluid. They also found a solution by means of a perturbation analysis for a viscoelastic fluid. Later, Berker [3] proved that there is the existence of an infinite number of nontrivial solutions to the Navier-Stokes equations between eccentric rotating disks. Rajagopal [4] showed that this motion is one with constant stretch history. We refer the reader to the papers by Rajagopal [5] and Ersoy [6, 7] for a detailed list of references related to the flows between eccentric rotating disks.
Time-dependent flows between eccentric rotating disks have also attracted the attention of researchers [8–13]. In these studies, the unsteady flows are generated by the sudden motion of the disks rotating with the same angular velocity and the disks are not exposed to oscillation. Erdoğan [14, 15] was the first to study the unsteady motion induced by oscillations. He [14] studied the flow that the disks start to rotate eccentrically and the lower disk executes oscillations while the disks are initially rotating about a common axis. He [15] studied the flow that the disks start to rotate eccentrically and both the disks execute oscillations in the same direction while the disks are initially rotating about a common axis. It is clear that these flows are not symmetrical as indicated by Erdoğan [14, 15].
In addition, the reader may consult the references [16–24] for the analytical solutions about the unsteady flows induced by eccentric rotations of an oscillating disk and a fluid at infinity under various effects.
In this paper, the disks are initially rotating about noncoincident axes. For this reason, the initial condition is the solution obtained by Abbott and Walters [2]. The disks start to execute oscillations in their own planes and in the opposite directions, and thus the symmetrical condition is satisfied at all times. In order to obtain a more general solution, the oscillating disk velocity has two components. The problem is solved for both small and large times, and the two solutions are matched at a specific value of time. In other words, the velocity field is obtained for all times. The influences of the parameters acting on the flow are elucidated with the help of the figures.
2. Basic Equations
The flow field of the problem is bounded by two disks located at z=h and z=-h. Initially, the top and bottom disks are rotating about the z′- and z′′-axes with the same angular velocity Ω, respectively. The distance between the axes of rotation is shown in the y-direction by 2ℓ and the region between the disks is occupied by an incompressible Newtonian fluid. The motion of the fluid is examined after the disks start to execute oscillations in their own planes and in the opposite directions. The upper and lower disks oscillate in their own planes with the velocities U and -U, respectively, where U=(Uxsinnt,Uysinnt,0) and n is the frequency of the oscillation. It should be emphasized that the distance between the axes of rotation is fixed during the motion. Furthermore, a physical reality of the oscillations of the disks requires that the sine oscillation is more reasonable than the cosine oscillation in this paper. The geometry of the problem is shown in Figure 1.
Flow geometry.
Therefore, the initial and boundary conditions can be written in the following form:(2.1a)u=-Ωy+f^(z),ν=Ωx+g^(z)att=0for-h≤z≤h,(2.1b)u=-Ω(y-ℓ)+Uxsinnt,ν=Ωx+Uysinntatz=hfort≥0,(2.1c)u=-Ωy,ν=Ωxatz=0fort≥0,(2.1d)u=-Ω(y+ℓ)-Uxsinnt,ν=Ωx-Uysinntatz=-hfort≥0,where u and ν denote the velocity components in the x- and y-directions, respectively. The functions f^(z) and g^(z) obtained by Abbott and Walters [2] represent the eccentric symmetrical rotation for a Newtonian fluid and are given by
(2.2)f^(z)+ig^(z)=ΩℓsinhKzsinhKh,
where i=-1, K=Ω/(2ν)(1+i), and ν denotes the kinematic viscosity of the fluid. Equation (2.1c) reflects the symmetrical condition.
The velocity field for the flow under consideration is given by
(2.3)u=-Ωy+f(z,t),ν=Ωx+g(z,t).
We should note that this flow does not bring out a velocity component in the z-direction. Substituting (2.3) into the Navier-Stokes equations, one obtains
(2.4)ν∂2f∂z2-∂f∂t+Ωg=C1(t),ν∂2g∂z2-∂g∂t-Ωf=C2(t).
Using (2.1a)–(2.1d) and (2.3), we get(2.5a)f(z,0)=f^(z),g(z,0)=g^(z),(2.5b)f(±h,t)=±(Ωℓ+Uxsinnt),g(±h,t)=±Uysinnt,(2.5c)f(0,t)=0,g(0,t)=0.Introducing F(z,t)=f(z,t)+ig(z,t) and using (2.4), we have
(2.6)ν∂2F∂z2-∂F∂t-ΩiF=C(t).
The symmetrical condition gives C(t)=0, which implies the absence of the Poiseuille-type pressure gradient. The conditions for (2.6) become
(2.7)F(z,0)=f^(z)+ig^(z),F(±h,t)=±(Ωℓ+Uxsinnt)±iUysinnt,F(0,t)=0.
3. Solution for Small Times
Putting F(z,t)=H(z,t)e-iΩt, (2.6) takes the form
(3.1)ν∂2H∂z2=∂H∂t,
with the conditions
(3.2)H(z,0)=ΩℓsinhKzsinhKh,H(±h,t)=±[Ωℓ+(Ux+iUy)sinnt]eiΩt,H(0,t)=0.
The Laplace transform of H(z,t) is defined by the equation
(3.3)H-(z,s)=∫0∞H(z,t)e-stdt.
Taking the Laplace transform of (3.1) with the conditions (3.2), we have
(3.4)H-′′-sνH-=-ΩℓνsinhKhsinhKz,(3.5)H-(±h)=±Ωℓs-iΩ±(Ux+iUy)nn2+(s-iΩ)2,H-(0)=0,
where a prime denotes differentiation with respect to z. Applying the conditions (3.5), the solution of (3.4) is
(3.6)H-=Ωℓs-iΩsinhKzsinhKh+(Ux+iUy)nn2+(s-iΩ)2sinhs/νzsinhs/νh.
Letting φ=s/ν, (3.6) can be written as follows:
(3.7)H-=1sΩℓ1-iΩ/ssinhKzsinhKh-(Uy-iUx)2s[11+i(n-Ω)/s-11-i(n+Ω)/s]e-φ(h-z)-e-φ(h+z)1-e-2φh.
It is well known that the series ∑q=0∞Xq converges to (1-X)-1 for |X|<1. Using this binomial series, it is possible to obtain the solution for small times. Equation (3.7) can be written in the following form:
(3.8)H-=ΩℓsinhKzsinhKh∑q=0∞(iΩ)qsq+1-(Uy-iUx)2×∑p=0∞∑m=0∞{([i(Ω-n)]p-[i(n+Ω)]p)e-φ(h-z+2hm)-e-φ(h+z+2hm)sp+1}.
The inverse Laplace transform of (3.8) gives
(3.9)FΩℓ=sinhR/2(1+i)ζsinhR/2(1+i)-(Vy-iVx)2(cosτ-isinτ)×∑p=0∞ip(4τ)p[(1-k)p-(1+k)p]∑m=0∞[i2perfc1-ζ+2m2τ/R-i2perfc1+ζ+2m2τ/R],
or
(3.10)f-=P(1)P(ζ)+Q(1)Q(ζ)Δ+Vxsinτ-Vycosτ2(-α2T4+α4T8-α6T12+α8T16-⋯)-Vxcosτ+Vysinτ2(α1T2-α3T6+α5T10-α7T14+⋯),g-=P(1)Q(ζ)-Q(1)P(ζ)Δ+Vxsinτ-Vycosτ2(α1T2-α3T6+α5T10-α7T14+⋯)+Vxcosτ+Vysinτ2(-α2T4+α4T8-α6T12+α8T16-⋯),
where
(3.11)f-=fΩℓ,g-=gΩℓ,R=Ωh2ν,ζ=zh,τ=Ωt,Vx=UxΩℓ,Vy=UyΩℓ,k=nΩ,P(ζ)=sinhR2ζcosR2ζ,Q(ζ)=coshR2ζsinR2ζ,Δ=[P(1)]2+[Q(1)]2,αr=(4τ)r[(1-k)r-(1+k)r],Tr=∑m=0∞[irerfc1-ζ+2m2τ/R-irerfc1+ζ+2m2τ/R].
The properties of a function that characterizes the functions Tr defined here are given by Ersoy [25]. The series solutions shown by (3.10) converge rapidly for small times but converge slowly when τ increases. These solutions cannot be used for very large times. For this reason, it is necessary to introduce a different method for large times.
4. Solution for Large Times
For large times, we suggest a solution of the form
(4.1)F(z,t)=F0(z)+F1(z)cosnt+F2(z)sinnt,
where F0(z) corresponds to the case of n=0. Substituting (4.1) into (2.6), we get(4.2a)νF0′′-ΩiF0=0,(4.2b)νF1′′-nF2-ΩiF1=0,(4.2c)νF2′′+nF1-ΩiF2=0.The boundary conditions for (4.2a)–(4.2c) are
(4.3)F0(±h)=±Ωℓ,F1(±h)=0,F2(±h)=±(Ux+iUy),F0(0)=0,F1(0)=0,F2(0)=0.
From the solutions of (4.2a)–(4.2c) by the conditions (4.3), we have
(4.4)FΩℓ=sinhR/2(1+i)ζsinhR/2(1+i)-(Vy-iVx)2[sin(k-1)R/2(1+i)ζsin(k-1)R/2(1+i)-sinh(k+1)R/2(1+i)ζsinh(k+1)R/2(1+i)]coskτ+(Vx+iVy)2[sin(k-1)R/2(1+i)ζsin(k-1)R/2(1+i)+sinh(k+1)R/2(1+i)ζsinh(k+1)R/2(1+i)]sinkτ,
or
(4.5)f-=P(1)P(ζ)+Q(1)Q(ζ)Δ+{Vx[J(1)K(ζ)-K(1)J(ζ)2L-A(1)B(ζ)-B(1)A(ζ)2D]+Vy[J(1)J(ζ)+K(1)K(ζ)2L-A(1)A(ζ)+B(1)B(ζ)2D]}coskτ+{Vx[A(1)A(ζ)+B(1)B(ζ)2D+J(1)J(ζ)+K(1)K(ζ)2L]-Vy[A(1)B(ζ)-B(1)A(ζ)2D+J(1)K(ζ)-K(1)J(ζ)2L]}sinkτ,g-=P(1)Q(ζ)-Q(1)P(ζ)Δ+{Vx[A(1)A(ζ)+B(1)B(ζ)2D-J(1)J(ζ)+K(1)K(ζ)2L]+Vy[J(1)K(ζ)-K(1)J(ζ)2L-A(1)B(ζ)-B(1)A(ζ)2D]}coskτ+{Vx[A(1)B(ζ)-B(1)A(ζ)2D+J(1)K(ζ)-K(1)J(ζ)2L]+Vy[A(1)A(ζ)+B(1)B(ζ)2D+J(1)J(ζ)+K(1)K(ζ)2L]}sinkτ,
where
(4.6)A(ζ)=cosh(k-1)R2ζsin(k-1)R2ζ,B(ζ)=sinh(k-1)R2ζcos(k-1)R2ζ,J(ζ)=sinh(k+1)R2ζcos(k+1)R2ζ,K(ζ)=cosh(k+1)R2ζsin(k+1)R2ζ,D=[A(1)]2+[B(1)]2,L=[J(1)]2+[K(1)]2.
The solution for k=1 is
(4.7)f-=P(1)P(ζ)+Q(1)Q(ζ)Δ+[-Vyζ2+VxM(1)N(ζ)-N(1)M(ζ)2E+VyM(1)M(ζ)+N(1)N(ζ)2E]cosτ+[Vxζ2+VxM(1)M(ζ)+N(1)N(ζ)2E-VyM(1)N(ζ)-N(1)M(ζ)2E]sinτ,g-=P(1)Q(ζ)-Q(1)P(ζ)Δ+[Vxζ2-VxM(1)M(ζ)+N(1)N(ζ)2E+VyM(1)N(ζ)-N(1)M(ζ)2E]cosτ+[Vyζ2+VxM(1)N(ζ)-N(1)M(ζ)2E+VyM(1)M(ζ)+N(1)N(ζ)2E]sinτ,
where
(4.8)M(ζ)=sinhRζcosRζ,N(ζ)=coshRζsinRζ,E=[M(1)]2+[N(1)]2.
The solutions shown by (4.5) and (4.7) are not valid for small times. Table 1 compares the two different solutions of f- for the given specific values when the time varies.
Comparison of the two solutions of f- when τ increases (ζ=0.5, R=10, Vx=1, Vy=1, k=0.5).
Small-time solution (f-)
Large-time solution (f-)
τ=0.1
0.1799927555
0.0407381339
τ=0.5
0.1913381458
0.1385613370
τ=0.8
0.2231963485
0.2132298660
τ=1
0.2542444915
0.2626974023
τ=1.5
0.3528926478
0.3816911322
τ=2
0.4604354856
0.4881440727
τ=3
0.6281582038
0.6381439189
τ=4
0.6757990334
0.6759717470
τ=5
0.5936835522
0.5923659856
τ=6
0.4083508823
0.4077962408
τ=7
0.1674919804
0.1674516233
5. Results and Discussion
In this paper, the motion of the fluid between the disks executing oscillations in their own planes and in the opposite directions is studied while they are initially rotating noncoaxially. There is no flow perpendicular to the disks due to the fact that they are rotating with the same angular velocity at all times. The fluid layer in the plane z=0 rotates as if a rigid body about z-axis, which implies that the symmetrical condition is satisfied, since the disks oscillate in the opposite directions. The effects of all the parameters acting on the flow are revealed by means of Figures 2, 3, 4, 5, 6, and 7. Figures 2–4 display the effect of the frequency of oscillation. The influence of the oscillating direction of the disks on the flow is illustrated in Figures 5 and 6. Figure 7 shows how the flow depends on the Reynolds number. The main findings of the present analysis are pointed out below.
The solution of the problem is obtained for both small and large times. At some specific time, it is shown that the solutions for small times are in good agreement with those for large times. Thus, the velocity field is determined at all times.
It is shown that the velocity increases for the same small times when the frequency of oscillation increases. With the decrease of the frequency, the periodic motion takes place later.
When the oscillation takes place along the eccentricity direction, the y-component of the translational velocity is considerably affected, but the change in the x-component is almost imperceptible. On the other hand, a reverse effect is observed when the disks are forced to oscillate in the x-direction. In this case, the change in the x-component of the translational velocity is noticed clearly but the change in the y-component is insignificant.
Increasing the Reynolds number has the effect of decreasing the thickness of the boundary layer.
It is observed that the periodic motion presupposed in the solution for large times occurs.
As it is expected, there exists a phase lag between the flow velocity and the disk oscillation when the periodic motion occurs.
It is shown that a solution can be obtained even when the angular velocity of the disk is equal to the frequency of oscillation.
Variations of f- and g- with τ for k=0.5 (R=10, Vx=1, Vy=1).
Variations of f- and g- with τ for k=1 (R=10, Vx=1, Vy=1).
Variations of f- and g- with τ for k=1.25 (R=10, Vx=1, Vy=1).
Variations of f- and g- with Vx for Vy=0 (R=10, k=0.5, τ=0.5).
Variations of f- and g- with Vy for Vx=0 (R=10, k=0.5, τ=0.5).
Variations of f- and g- with R (k=0.5, τ=0.5, Vx=2, Vy=2).
Acknowledgments
The author would like to express his sincere thanks to the referees for their valuable comments and suggestions.
MaxwellB.ChartoffR. P.Studies of a polymer melt in an orthogonal rheometer196591415210.1122/1.548979AbbottT. N. G.WaltersK.Rheometrical flow systems—part 2: theory for the orthogonal rheometer, including an exact solution of the Navier-Stokes equations197040120521310.1017/S0022112070000125BerkerR.An exact solution of the Navier-Stokes equation: The vortex with curvilinear axis19822022172302-s2.0-001990003910.1016/0020-7225(82)90017-9ZBL0487.76039RajagopalK. R.On the flow of a simple fluid in an orthogonal rheometer198279139472-s2.0-002000282610.1007/BF02416565ZBL0513.76002RajagopalK. R.Flow of viscoelastic fluids between rotating disks1992341852062-s2.0-000039275410.1007/BF00417912ZBL0747.76018ErsoyH. V.An approximate solution for flow between two disks rotating about distinct axes at different speeds20072007162-s2.0-3424723771310.1155/2007/3671836718ZBL1141.76066V. ErsoyH.On the locus of stagnation points for a Maxwell fluid in an orthogonal rheometer2009356606642-s2.0-77952050835ErdoğanM. E.Unsteady viscous flow between eccentric rotating disks19953057117172-s2.0-001057092010.1016/0020-7462(95)00030-RErsoyH. V.Unsteady flow due to a sudden pull of eccentric rotating disks20013933433542-s2.0-003525411110.1016/S0020-7225(00)00040-9ErsoyH. V.Unsteady flow due to concentric rotation of eccentric rotating disks20033833253342-s2.0-003771698110.1023/A:1023374214783ZBL1032.76711MohyuddinM. R.Unsteady MHD flow due to eccentric rotating disks for suction and blowing20073131231352-s2.0-34548249631GuriaM.JanaR. N.GhoshS. K.Unsteady MHD flow between two disks with non-coincident parallel axes of rotation20073454254332-s2.0-3454840952710.1615/InterJFluidMechRes.v34.i5.30MajiS. L.GharaN.JanaR. N.DasS.Unsteady MHD flow between two eccentric rotating disks2009138796ErdoǧanM. E.Flow due to parallel disks rotating about non-coincident axis with one of them oscillating in its plane1999346101910302-s2.0-034231037110.1016/S0020-7462(98)00072-9ErdoganM. E.Unsteady flow between two eccentric rotating disks executing non-torsional oscillations20003546916992-s2.0-003422750710.1016/S0020-7462(99)00051-7ZBL1006.76103HayatT.AsgharS.SiddiquiA. M.Unsteady flow of an oscillating porous disk and a fluid at infinity19993442592652-s2.0-003330753010.1023/A:1004751322961ZBL0959.76089ErdoǧanM. E.Flow induced by non-coaxial rotation of a disk executing non-torsional oscillations and a fluid rotating at infinity20003821751962-s2.0-034231350010.1016/S0020-7225(99)00017-8HayatT.ZamuradM.AsgharS.SiddiquiA. M.Magnetohydrodynamic flow due to non-coaxial rotations of a porous oscillating disk and a fluid at infinity20034111117711962-s2.0-003741120310.1016/S0020-7225(03)00004-1HayatT.MumtazS.EllahiR.MHD unsteady flows due to non-coaxial rotations of a disk and a fluid at infinity20031932352402-s2.0-014171935010.1007/BF02484485HayatT.EllahiR.AsgharS.SiddiquiA. M.Flow induced by non-coaxial rotation of a porous disk executing non-torsional oscillations and a second grade fluid rotating at infinity20042865916052-s2.0-214281804910.1016/j.apm.2003.10.011HayatT.EllahiR.AsgharS.Unsteady periodic flows of a magnetohydrodynamic fluid due to noncoaxial rotations of a porous disk and a fluid at infinity2004401-21731792-s2.0-994423895910.1016/j.mcm.2003.09.035ZBL1096.76060HayatT.EllahiR.AsgharS.Unsteady magnetohydrodynamic non-Newtonian flow due to non-coaxial rotations of disk and a fluid at infinity2007194137492-s2.0-3374899953910.1080/00986440600642868GuriaM.DasB. K.JanaR. N.Oscillatory flow due to eccentrically rotating porous disk and a fluid at infinity20074254874932-s2.0-3514885872410.1007/s11012-007-9071-9ZBL1163.76441HayatT.EllahiR.AsgharS.Hall effects on unsteady flow due to non-coaxially rotating disk and a fluid at infinity200819589589762-s2.0-4314909079310.1080/00986440801906575ErsoyH. V.Examination of a special function defined by an integral201221616410.4236/ajcm.2012.21008