Unsteady Flow Produced by Oscillations of Eccentric Rotating Disks

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 xand y-directions on the flow is examined. The influence of the Reynolds number is also investigated.


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 gan 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 gan 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.

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 U x sin nt, U y sin nt, 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.
Therefore, the initial and boundary conditions can be written in the following form: 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 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 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.5c
Introducing F z, t f z, t ig z, t and using 2.4 , we have 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

Solution for Small Times
Putting F z, t H z, t e −iΩt , 2.6 takes the form H 0, t 0.

3.2
The Laplace transform of H z, t is defined by the equation Taking the Laplace transform of 3.1 with the conditions 3.2 , we have where a prime denotes differentiation with respect to z. Applying the conditions 3.5 , the solution of 3.4 is Letting ϕ s/ν, 3.6 can be written as follows:

3.7
It is well known that the series ∞ q 0 X q 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
The inverse Laplace transform of 3.8 gives where

3.11
The properties of a function that characterizes the functions T r 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.

Solution for Large Times
For large times, we suggest a solution of the form where F 0 z corresponds to the case of n 0. Substituting 4.1 into 2.6 , we get The boundary conditions for 4.2a -4.2c are From the solutions of 4.2a -4.2c by the conditions 4.3 , we have Mathematical Problems in Engineering or

4.5
where The solution for k 1 is

4.8
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.

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.i 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.
ii 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.
iii When the oscillation takes place along the eccentricity direction, the y-component of the translational velocity is considerably affected, but the change in the xcomponent 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 vi As it is expected, there exists a phase lag between the flow velocity and the disk oscillation when the periodic motion occurs.vii It is shown that a solution can be obtained even when the angular velocity of the disk is equal to the frequency of oscillation.

5 −Figure 2 :
Figure 2: Variations of f and g with τ for k 0.5 R 10, V x 1, V y 1 .

Figure 3 :
Figure 3: Variations of f and g with τ for k 1 R 10, V x 1, V y 1 .

Figure 7 :
Figure 7: Variations of f and g with R k 0.5, τ 0.5, V x 2, V y 2 .

Table 1 :
Comparison of the two solutions of f when τ increases ζ 0.5, R 10, V x 1, V y 1, k 0.5 .