Flow of viscoelastic fluids between rotating disks

Few boundary-value problems in fluid mechanics can match the attention that has been accorded to the flow of fluids, Newtonian and non-Newtonian, between parallel rotating disks rotating about a common axis or about distinct axes. An interesting feature which has been recently observed is the existence of solutions that are not axially symmetric even in the case of flow due to the rotation of disks about a common axis. In this article we review the recent efforts that have been expended in the study of both symmetric and asymmetric solutions in the case of both the classical linearly viscous fluid and viscoelastic fluids.The support of the Air Force Office of Scientific Research is gratefully acknowledged.     

On finite deformation of viscoelastic rotating disks

The governing equations for a rotating disk under large elastic or viscoelastic deformations, are formulated in terms of two coupled first-order ordinary differential-integral equations with explicit derivatives. The constitutive equations developed by Christensen are used for the formulation. A recurrent formula for the hereditary integrals is obtained. This recurrent formula makes it possible for the same governing equations and numerical methods to be applied to both elastic and viscoelastic problems.     

Numerical investigations of flow of a viscoelastic fluid between rotating coaxial disks

Summary A numerical investigation is made of the problem of flow of a viscoelastic fluid of the Rivlin-Ericksen type between rotating coaxial disks. The finite-difference analogues of the governing nonlinear ordinary differential equations are written and the resulting equations are solved using point successive overrelaxation method (SOR) under the appropriate boundary conditions. Two cases of interest are treated, namely, when one of the disk is at rest while the other rotates with a constant angular velocity and when both the disks rotate with constant angular velocities but in the opposite directions. Typical examples of Reynolds numberR in the range 10 R 1000 are described for various values of the non-Newtonian parameters and results are compared with those for a classical viscous fluid.

---

Zusammenfassung Die Strömung einer viskoelastischen Flüssigkeit vom Rivlin-Ericksen-Typ zwischen zwei rotierenden koaxialen Kreisscheiben wird numerisch untersucht. Die das Problem beschreibenden nichtlinearen Differentialgleichungen werden in analogen Gleichungen für finite Differenzen umgeschrieben und unter angepaßten Randbedingungen mit einem punktweisen Über-Relaxations-Verfahren (SOR) gelöst. Es werden zwei interessierende Fälle behandelt, zuerst der Fall, bei dem eine Platte ruht und die andere mit konstanter Winkelgeschwindigkeit rotiert, als zweiter derjenige, bei dem beide Platten mit gleichen konstanten Winkelgeschwindigkeiten, jedoch entgegengesetztem Drehsinn rotieren. Typische Beispiele für Reynoldszahlen zwischen 10 und 1000 werden für verschiedene Werte der nicht-newtonschen Parameter beschrieben und die Ergebnisse mit denen für newtonsche Flüssigkeiten verglichen.

---

---

With 9 figures     

Mathematical and numerical modelling of large creep deformations for annular rotating disks

A simulation model is presented for the creep process of the rotating disks under the radial pressure in the presence of body forces. The finite strain theory is applied. The material is described by the Norton-Bailey law generalized for true stresses and logarithmic strains. A mathematical model is formulated in the form of a set of four partial differential equations with respect to the radial coordinate and time. Necessary initial and boundary conditions are also given. To make the model complete, a numerical procedure is proposed. The given example shows the effectiveness of this procedure. The results show that the classical finite element method cannot be used here because both the geometry and the loading (body forces) change with the time in the creep process, and the finite elements need to be redefined at each time step.     

Optoelectronic-strain-measurement system for rotating disks

A recently developed optical technique is presented in the following paper for measuring in-plane deformations and strains of rotating hardware. The technique is fully described. Data are presented which were obtained in preliminary tests using the system on a rotating disk. The experiments show that the system is capable of making radial-deformation measurements to within 9.5 μm over a 10-mm range.     

Theoretical and experimental investigation of the flow of a viscoelastic fluid in the gap between two rotating disks

A theory of the nonlinear viscoelastic behavior of polymer fluids has been constructed in [1]. The theory was used in [2] to investigate the motion of a nonlinear viscoelastic medium under steady and unsteady deformation rates in simple shear flow, and a comparison was made with experiment. The experiments in [2], which were performed on a cone-plate Weissenberg rheogoniometer, indicate that this arrangement is unsuitable for measurements of normal stresses under unsteady conditions in fluids with a fairly high viscosity. Below, we will show the suitability of using a disk-disk Weissenberg rheogoniometer to measure normal stresses in this case for unsteady conditions (transition to steady flow and stress relaxation). In this regard, a theoretical study of the flow of a viscoelastic fluid in the gap between rotating disks is needed. Note that in this case new information will be obtained from a comparison with simple uniform shear flow, since in the flow of a polymer between two disks all three normal stress components contribute to the axial force, while in the gap between a cone and a plate only the first normal stress difference contributes to the normal force.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 25–30, March–April, 1976.     

Fracture-mechanics investigations of cracks in rotating disks

Stress-intensity factorsK were determined both analytically and by using a photoelastic method for the simple case of a rotating solid disk containing radial cracks. Good agreement is found not only between the calculated and the experimentalK factors, but also between the static and the dynamic toughness values determined in ASTM tension tests and spin-burst tests. This confirms the applicability of linear-elastic fracture mechanics and the validity of the brittle-fracture criterion. In addition, the use of the simple superposition procedure is justified as a basis for the analysis. The possibilities of and the limitations on applying these results to practical situations are considered.     

Flow between two unsteadily rotating disks

Summary This is a theoretical investigation of the unsteady laminar flow of a viscous incompressible fluid between two infinite parallel disks, which are rotating with angular velocities varying with time. The solution is obtained in the form of a series expansion about the quasi-steady state. The deviation of the actual instantaneous state of the flow from the quasi-steady state is determined.     

Stretching surface in rotating viscoelastic fluid

The boundary layer flow over a stretching surface in a rotating viscoelastic fluid is considered. By applying a similarity transformation, the governing partial differ- ential equations are converted into a system of nonlinear ordinary differential equations before being solved numerically by the Keller-box method. The effects of the viscoelastic and rotation parameters on the skin friction coefficients and the velocity profiles are thor- oughly examined. The analysis reveals that the skin friction coefficients and the velocity in the x-direction increase as the viscoelastic parameter and the rotation parameter in- crease. Moreover, the velocity in the y-direction decreases as the viscoelastic parameter and the rotation parameter increase.     

Frictional contact of two rotating elastic disks

We study the problem of constrained uniform rotation of two precompressed elastic disks made of different materials with friction forces in the contact region taken into account. The exact solution of the problem is obtained by the Wiener-Hopf method.An important stage in the study of rolling of elastic bodies is the Hertz theory [1] of contact interaction of elastic bodies with smoothly varying curvature in the contact region under normal compression. Friction in the contact region is assumed to be negligible. If there are tangential forces and the friction in the contact region is taken into account, then the picture of contact interaction of elastic bodies changes significantly. Although the normal contact stress distribution strictly follows the Hertz theory for bodies with identical elastic properties and apparently slightly differs from the Hertz diagram for bodies made of different materials, the presence of tangential stresses results in the splitting of the contact region into the adhesion region and the slip region. This phenomenon was first established by Reynolds [2], who experimentally discovered slip regions near points of material entry in and exit from the contact region under constrained rolling of an aluminum cylinder on a rubber base. The theoretical justification of the partial slip phenomenon in the contact region, discovered by Reynolds [2], can be found in Carter [3] and Fromm [4]. Moreover, Fromm presents a complete solution of the problem of constrained uniform rotation of two identical disks. Apparently, Fromm was the first to consider the so-called “clamped” strain and postulated that slip is absent at the point at which the disk materials enter the contact region.Ishlinskii [5, 6] gave an engineering solution of the problem on slip in the contact region under rolling friction. Considering the problem on a rigid disk rolling on an elastic half-plane, we model this problem by an infinite set of elastic vertical rods using Winkler-Zimmermann type hypotheses. Numerous papers of other authors are surveyed in Johnson’s monograph [7].The exact solution of the problem on the constrained uniform rotation of precompressed rigid and elastic disks under the assumptions of Fromm’s theory is contained in the papers [8, 9]. In the present paper, we generalize the solution obtained in [8, 9] to the case of two elastic disks made of different materials.     

Slow viscoelastic radial flow between parallel disks

A perturbation analysis is presented for the steady-state radial flow of a third-order fluid between two parallel disks. The results include previous perturbation analyses in which various other rheological models were used. The pressure drop needed to maintain the radial flow is less than that for the Newtonian creeping-flow solution because of fluid inertia and shear-thinning viscosity, whereas the normal stresses have the opposite effect. Possible use of the “radial flow viscometer” for experimental evaluation of second-order constants is also discussed. Finally, molecular stretching in the flow system is examined using the elastic dumbbell model for a polymer molecule.     

高速转动圆盘流动数值模拟研究

重点研究高速离心压气机叶轮与机匣间的间隙流动及其温度分布。研究将离心压气机简化为高速转动圆盘,搭建了相关实验平台,并开展了相应的数值模拟研究。通过改变转动圆盘的转速和轴向进入的冷却流的流量,研究了转速和流量对于间隙内温度和速度分布的影响。结果显示,转速是影响温度变化的最主要因素,转速越大,温度越高;同等幅度的流量变化对温度的影响则较小。研究发现,在实验和模拟对应的大雷诺数条件下,无量纲的速度分布基本不受到圆盘转速、冷却流量和温度场的影响。     

Attractors of a rotating viscoelastic beam

We investigate the non-linear oscillations of a rotating viscoelastic beam with variable pitch angle. The governing equations of motion are two coupled partial differential equations for the longitudinal and transversal displacements. Using a perturbation technique and Galerkin's projection, we reduce the equations of motion to a non-autonomous ordinary differential equation. Our regular perturbation technique is based on the expansion of longitudinal displacement and the amplitude of first transversal mode in terms of a small parameter. We numerically generate the Poincaré maps of the reduced equations and reveal that the system exhibits regular and chaotic attractors. The regular attractors are stable limit-cycles that are relevant to stable, short-period oscillations of the beam. A bifurcation analysis has also been performed when the pitch angle is constant.