Super-Laminiir Flow Beaills I I S Seils
Super-Laminiir Flow Beaills I I S Seils
Super-Laminiir Flow Beaills I I S Seils
NO
GEBniBfiTrf ;•• % ^ I
Neither the United Statess nor the Commission, nor any person acting on
behalf of the Commission:
B. Assumes any liabilities with respect to the use of, or for damages
resulting from the use of any information, apparatus, method, or
process disclosed in this report.
« G - C 5 ^ . MN , / ^
Proceedings of the Technical Meeting on:
&5SIRACTS
11 S\JCXiB^- SCI^^
E. B. Arwas, Editor
December, 1966
MTI-2349
1
k TABLE OF CONTENTS
Page
FOREWORD iii
6.1
4
6.1 - Experiments with Hydrodynamic Journal
Bearings of Various Materials and Designs
in Sodium at Temperatures to 800 F 6.1.1
F. T. Schuller, W. J. Anderson and
Z. Nemeth, National Aeronautics and
Space Administration
4
f
Ill
FOREWORD
The technical meeting of November 1 and 2, 1966 was held to provide for an ex-
change of Information between the various groups active in research, development
and application programs in this area. During the course of three technical
sessions, seventeen presentations were made by representatives of U.S. Government
Agencies and commercial contractors, on recent progress, current effort and anti-
cipated future requirements in process-fluid lubrication and its applications
to nuclear and aerospace turbomachinery.
This report contains all these presentations, reproduced from the manuscripts
furnished by the authors to the editor. For continuity and ease of reference,
these presentations have been assembled here according to their principal topic
into the sections listed below:
Section 1 - Summary
Section 2 - Introductory Address
Section 3 - Background Review
Section 4 - Technology I - Fundamentals
Section 5 - Technology II - Seals
Section 6 - Technology III - Bearings and Materials
Section 7 - Technology IV - Rotor-Bearing Dynamics
Section 8 - Application to Large Turbomachines
by
Review
Technology
Applications
The following represents a brief summary of the meeting. More detailed discussions
will be found in the other papers which make up this report.
2. Superlaminar theory neglecting inertia effects has been developed for the
calculation of load, frictional losses, stiffness and damping for several
bearing types. This theory is in excellent agreement with practice.
by
Nicholas Grossman
Chief, Special Technology Branch
Division of Reactor Development and Technology
USAEC
Many of you are aware of the long standing, well organized and highly
successful program in gas lubrication coordinated by the Fluid Mechanics Branch
of the Office of Naval Research. Each of us here today is a direct beneficiary
of the farsighted and effectively administered ONR program. There are a number
of enthusiastic supporters who materially contributed to this undertaking, but I
feel we should specifically mention two of our friends, whose unfailing devotion
stands out as an example for us all: Mr, Stanley Doroff of ONR, who is directly
responsible for procuring the support,defending the program, and dispensing the
funds with Solomon-like impartiality and wisdom -- and Professor Dudley Fuller
of Columbia University, who in the role of Technical Secretary of the Coordinating
group, has been the program chairman, recording secretary, conciliator and tower
of strength.
I recall that during the 1965 summer meeting of the ONR Gas Bearing Co-
ordinating Group, Columbia University was our obliging host. We were fortunate
to have Dean Robert Dunning as our keynote speaker. As you know. Dean Dunning
was one of the pioneers who helped move atomic energy from the physics laboratory
into the realm of engineering.
During his talk to the ONR Gas Bearing Group in 1965, Dean Dunning mentioned
how he and his engineering associates working on the design of the gaseous
diffusion plant some twenty years earlier, recognized the practical advantages
of gas lubrication and made record notes accordingly! Thus if we wish to
establish a date to connect gas lubrication and atomic energy, we can go back to
the wartime Manhattan Engineering District and state that engineers working with
atomic energy on an industrial scale recognized the natural affinity between
process-fluid lubrication and atomic energy at the very outset. In my opinion,
had rigid security classification not been necessary during the war, the develop-
ment of process-fluid lubrication would have progressed at a much faster rate,
and certainly would have been a practical industry application at a much earlier
date.
It is well known that the idea of gas-lubricated machinery was not incorp-
orated in the gaseous diffusion plant built during the war. It is at this
junction that we must appreciate the foresight of the Fluid Mechanics Branch of
ONR for recognizing the industrial potential of gas lubrication and organizing
a sound research program to cover this vast field in an orderly way. Naturally,
a clever idea like gas lubrication cannot be kept as a monopoly of any one group,
and just as it was recognized in America, similar development efforts were pursued
in Europe.
One of the obvious consequences of the war in Europe was the acute shortage
of modern central electric power stations. Consequently, there was greater
impetus to develop central power stations -- and their construction was pushed
much more vigorously in Europe than in the United States. It is understandable.
therefore, that industrial applications of gas bearings in gas blowers appeared
in Europe a few years ago.
Coming back to the Atomic Energy Commission -- its role in the promotion
of industrial application is specified in Public Law 83-703, popularly known as
"The Atomic Energy Act of 1954." Section 3 states in part: "It is the purpose
of this Act to effectuate ....... a program of conducting, assisting, and
fostering research and development in order to encourage maximum scientific and
industrial progress." How successful the Commission has been in carrying out its
charter can be readily assessed in terms of nuclear power stations built or being
constructed, the use of radioisotopes and other byproducts now in common use,
and nuclear propulsion for naval purposes. What is amazing about the success
of atomic energy as an established segment of American industry is not that it
came into existence, but rather the phenomenal speed of that development which
surprised even its most optimistic advocates.
Today you will hear about turbine-driven circulators and power conversion
systems using water and liquid metal bearings and seals in a turbuluent regime.
Future developments will lead to the engineering application of two-phase
lubricated systems in turbomachinery Since the AEG is developing two-phase
systems, the Commission's interest in process-fluid lubrication is therefore
apparent. In feet, most of the advanced process-fluid lubricated machines
2.4
What then is the goal of the AEC for future development? I will restrict
my remarks to the program of the Division of Reactor Development and Technology.
Current commercial nuclear power plants use a uranium "burning" reactor operating
in the thermal neutron spectrum which converts a portion of the uranium into a new
fissionable material: plutonium. There is thus a newly created material that can
be used for power generation. However, the amount of new material produced is
small compared with the uranium that was used or "burned" in the fission process.
Therefore, looking at this process from the point of view of most effective
utilization and conservation of our national reserves, thermal reactors are con-
sumers of fissionable material. One of the amazing aspects of nuclear fission --
at least to us engineers -- is that in a suitable fast neutron spectrum more new
fissionable material can be "created" than is being "used". This concept is
popularly known as "breeder reactor". The orderly development of this system
is a high priority goal of civilian nuclear power development programs.
Before closing I want to state that the Space Nuclear Systems Division of the
AEC, as well as other government agencies, are also pursuing the development of
liquid metal bearings and seals, where their application offers attractive benefits
to their specific requirements. It is one of the purposes of this s5nnposium to
afford us all an opportunity to become better acquainted with the efforts of
other groups, and offer the benefit of our experience -- to one another. I am
very pleased with the progress to date and I hope you will share this pride.
Thank you.
3.1
BACKGROUND REVIEW
'-'11
V ' PROCESS-FLUID LUBRICATION OF TURBOMACHINERY BEARINGS
by
E. B. Arwas
Mechanical Technology Inc.
Latham, N. Y.
ABSTRACT
INTRODUCTION
The papers presented in this technical meeting have ranged over much of the recent
^
work as well as the current and future needs of this technology. Thus, one of the
presentations (Section 8.1) discussed the bearing and seal requirements for projected,
large, nuclear, liquid sodium facilites. Other presentations (Sections 6.2, 7.1,
and 8.3) reviewed recent and current turbomachinery developments involving water
or liquid metal lubricated bearings and seals operating in the turbulent flow
regime. In still other presentations (Sections 4 through 6 ) , some of the most
recent theoretical approaches, empiricisms and experimental analyses of turbulent
flow lubrication with low viscosity fluids are described. The important point is
that all this effort, comprising research, practical hardware development and
definition of specification for future applications, is proceeding simultaneously.
This is a healthy condition for a technology. It is the hope behind technical
meetings such as this one to provide forums where the hardware development engineers
and the research investigators are kept cognizant of one another's progress and
changing requirements
(b) prepare numerical procedures and computer programs for calculating the
steady-state and dynamic characteristics of various bearing geometries,
from the turbulent lubrication theory,
(c) test the validity of the theory by means of bearing and rotor-bearing
dynamics tests with different bearing geometries and rotor arrangements
and,
(d) apply the verified theoretical analysis and the computation procedures
in the design of bearings and seals for hardware developments.
3.3
The theoretical analysis and its experimental verification have been published in
contract reports and in the open literature (References 1 through 5). The results
obtained to date and the questions still remaining are noted in Section 4.1,
together with our current conclusions regarding the regions of bearing operation
where vortex flows, developed turbulence and fluid inertia effects either singly
or in combination, will govern performance. The information gathered in this
research is being applied to practical bearing designs, such as the ones for
HTGR circulator, which are described in Section 8.3.
Currently, under continuing AEC sponsorship this work on bearing and rotor-bearing
dynamic performance is being extended to operation with a liquid metal lubricant
(NaK). This test program is currently underway, the facility and the proposed
test program are briefly described in an Appendix to this paper.
Since the specifics of the effort at M.T.I, are covered here in the presentation
of other M.T.I, engineers, the balance of this presentation will be limited to a
general review of:
3.4
TABLE 1
CHARACTERISTIC PROPERTIES OF HIGH SPEED. OIL LUBRICATED
SLIDER BEARINGS
ADVANTAGES:
While load capacity is not generally a problem in high speed machines, the
relatively high viscosities of many oils, combined with their excellent boundary
lubrication properties offer a very desirable margin of safety in event of mal-
function or abnormal conditions. In event of a failure, such as the loss of
blades in a turbine and the consequent very high dynamic loads, the high ultimate
load capacity and the boundary lubrication properties may, in some cases,
reduce the severity of the damage to the rotor.
Oils are generally non-corrosive, which increases the choice of bearing materials.
Finally, the low surface tension and good wettability improve the start up and
we believe, on the basis of some fairly preliminary investigations (References 6
and 7) , that these properties also tend to reduce the likelihood of cavitation
erosion in most applications. This is because they tend to produce a steady
striation type of film rupture in the regions of negative static pressure in the
film, instead of small, gas or vapor filled bubbles which can then collapse as they
move to positive static pressure regions.
Despite these advantages,and they are important ones, oil lubrication is difficult
in many of the advanced, closed cycle systems required in nuclear, aerospace and
undersea application. Thus, oil fouling and contamination of some of these closed
systems cannot be tolerated and it is extremely difficult to prevent if oil
lubricated bearings are used. The problem in the way of development of dynamic,
zero-leakage seals are, we believe, far greater than those of process-fluid
lubrication.
Lubricating oils are also limited in their temperature capability and cannot be
used at temperatures exceeding say 600 to 700 F, or at very low temperatures
approaching the cryogenic range. Thus, complex bearing temperature control provisions
are needed, as well as use of heat dams or other undesirable design comprises.
Because of high viscosity, friction losses in high speed oil lubricated slider
bearings are large and represent a sizeable penalty on efficiency.
Hydrocarbon oils cannot sustain radio-activity for long periods without breaking
down, which rules them out in some nuclear applications.
The complexity of the auxiliary lube-oil system cannot be minimized. The complex
of lube pumps, coolers, filters, controls and other items is complicated and
bulky. Moreover, while its component parts are generally well developed they have
nonetheless proved to be a major source of unreliability. Finally, of course, the
lube-oil system sometimes constitutes a fire hazard.
TABLE 2
ADVANTAGES OF FULL-FILM PROCESS-FLUID LUBRICATION
IN HIGH SPEED. HIGH TEMPERATURE.
CLOSED CYCLE SYSTEMS
4, the bearings can operate at the level of radio activity of the system
Thus 9 in the case of gas or vapor lubrication, the self-generated bearing load
capacity is extremely small (of the order of a few PSI), due to the very low
absolute viscosities of gases and vapors. Accordingly, externally pressurized
bearings have to be used to support large steady state or dynamic loads.
The fluid film damping in gas or vapor lubricated bearings is relatively small,
and accurate rotor response analysis should be made to insure that system and
component resonant frequencies are outside the operating range. The compress-
ibility of gases and vapors complicates the analysis of gas bearings, introducing
an additional parameter (the compressibility number A) and making the lubrication
differential equation non-linear. An important practical consequence of the
non-linearity is that the stiffness and damping of gas-bearing films are highly
frequency dependent. In particular, gas bearing damping tends to diminish rapidly
when the film is subjected to high frequency excitations. Since excitations at
multiples of the running speed are not unusual in some turbomachines (e.g. alternator
rotors), they can produce large vibrations if they excite a system resonance,
due to the reduced damping of the gas film. It is thus particularly important to
analyse accurately the dynamic behavior of high speed, gas bearing supported
systems,
The low load capacity ®f self-acting gas bearings makes them relatively intolerant
of thermal distortions, which can adversely affect the film geometry or produce
misaligning moments and edge loads. This can be aggravated by the fact that the
thermal capacity of the gas flow through self-acting bearings is small, so that
the heat generated in the bearing is removed by conduction through the bearing
elements, resulting in thermal gradients and corresponding distortions. In past
applications J the problem of thermal distortions has been controlled by the use
of heat dams and heat shunts, by locating the bearings in isothermal regions of
the rotor and by external cooling. Currently, effort is being devoted to develop-
ment of flexure mounts and conformable surfaces to permit bearing operation in the
presence of large thermal gradients.
Gases provide little or no boundary lubrication so that in all cases (in parti-
cularly for self-acting gas bearings which do not use hydrostatic jacking and
where a large number of starts and stops are anticipated over the life of the
machine) the materials must be selected to sustain frequent rubs without damage.
Currently, experience with vapor lubrication is much smaller than that with gas
lubrication. However, since lubrication of the bearings located at the turbine
end of a rotor with the vapor supplied to the turbine (e.g. steam, metal vapors
or others) can potentially result in major simplification and increased relia-
bility of the turbomachine5 an increased level of effort in this area of
lubrication technology is anticipated. Here, the effects of phase change, bearing
erosion and, in the case of externally pressurized bearings, restrictor erosion
have to be considered (Ref. 6) .
With liquid lubricants, the problems, at least from the standpoint of load
capacity, fluid film damping and stability, are somewhat less critical than
with gas lubrication, due to the higher absolute viscosities or liquids and
the absence or compressibility effects.
Many of these liquids, however, have very low kinematic viscosities, leading to
onset of superlaminar flows in high speed bearing films. Table 4 on Page 3.10
TABLE 4 -- VISCOSITIES OF SOME PROCESS FLUIDS
ABSOLUTE KINEMATIC
ABSOLUTE VISCOSITY, n VISCOSITY KINEMATIC VISCOSITY,v VISCOSITY
FLUID (Ib-sec/in. ) RATIO* ...__{ in,, i/secl RATIO*
lists the absolute and kinematic viscosities of a number of fluids and compares
them with those of an SAE 10 oil. This illustrates the relative load capacities
and the turbulence onset speeds, since the absolute viscosity of the lubricant
is a measure of the potential bearing load carrying capability (for self-acting
bearings), while the speed at onset of turbulence is directly proportional to
the kinematic viscosity of the fluid.
All of these changes are important in high speed bearings so that, for rational
bearing and rotor-bearing dynamics design, it is necessary to know whether the
bearing film is laminar or turbulent and when fluid inertia affects are sizeable.
The appropriate theory should then be used to compute the steady-state and dynamic
performance characteristics of the bearing.
A recent study conducted jointly by the Army Engineers Reactor Group and MTI
(Ref. 7 ) , illustrated some of the potential advantages of process fluids
lubrication. In this study, the oil lubrication systems and the bearing power
losses of three closed loop gas turbines developed in connection with the U.S. Army
program on compact, mobile nuclear power conversion systems were studied. They
were then compared with a gas lubricated system (CSG-1) that has been proposed and
for which a dynamic simulator was built and evaluated on externally pressurized
gas bearings. Table 5, (Page3.13) reproduced from Ref. 7, lists the power losses
in each case. It is noted that here even for machines in the 2000 shaft H.P.
range, there is a significant reduction in power loss associated with process-
fluid lubrication. With machines having a smaller power output, the percentage
reduction in power loss will be still greater.
The dynamic simulator of the CSGl was built and tested at MTI The rotor weight
and design point speed of the simulator were 90 lbs. and 28,000 RPM respectively.
The unit was operated at up to 30,000 RPM. Figure 4 is a photograph of the simulate
parts. It should beroted, however, that full development and operational experience
with gas lubricated, closed loop gas turbines in the 2000 HP range has not been
undertaken, so that the potential long term operation of such units and their
tolerance to off-design conditions are not yet known.
To date, the practical experience with gas bearing turbomachinery has been princi-
pally with motor driven gas bearing compressors used mostly in nuclear and chemo-
nuclear installations. There are a number of such units operating in the United
States and in Europe (Ref. 8). The several machines operating in the United States
have accumulated over 40,000 hours of highly successful field operation, the
longest time on any one machine being 28,000 hours.
For fewer turbine driven gas bearing machines have been developed todate. Figure 5
shows the first such unit to have been developed and tested at its design speed
(24,000 RPM) and its design point turbine inlet temperature (1300 F). This unit
was developed and operated for the U.S. Bureau of Mines and the U.S. Atomic Energy
TABLE 5
COMPARISON OF BEARING SYSTEM PARASITIC LOSSES FOR OIL- AND GAS-LUBRICATED TURBOCOMPRESSORS
Notes: 1. Numbers marked with an asterisk are measured values. All other numbers are calculated
values.
2. Bearing loss data for the CSG-1 machine is based on the CSG-1 rotor-bearing simulator
described in Reference 2.
3. The above losses are for the turbocompressor bearing system only. Losses associated
with reduction gearing and generator bearings are not included.
HI I-2
3.14
Commission by MTI, This unit, which has a rated shaft turbine power of 84.8 HP,
is representative (in terms of gas circulation alone) of 100 to 300 HP circulators.
In terms of closed Brayton Cycle power generation, the turbocompressor is re-
presentative of the gas generator section for a 30 to 100 KW(e) gas-turbine power
plant. Over 125 hours of test operation have been accumulated with this machine
at turbine inlet temperatures ranging from 1000 to 1400 F.
Gas bearing supported compressors and expanders are also required for cryogenic
systems where very long, maintenance free operating life is required. Recently a
two stage regenerative helium compressor was operated at MTI at speeds up to
114,000 RPM.
This particular unit, shown in Figure 6, has a rotor weighing 1.5 lbs. which
comprises the two regenerative compressor stages, an internal cooling fan and
thrust plate and a 3 phase synchronous, hysteresis type drive motor. The rotor
is supported in two 3/4" diameter tilting pad bearings, designed for stable
operation at up to 150,000 RPM. The unit is cited here, because in the course of
developing it, some preliminary experimental investigations were made of the
influence of electro-magnetic forces on such small, high speed self-acting gas
bearings. When a clean, sinusoidal AC voltage was applied, the unit operated at
up to 114,000 RPM without difficulty. When, however, an approximately square
wave output was used, the system harmonics resulted in large, unbalanced magnetic
forces on the bearings which exceeded their capabilities. Figure 7 illustrates
orbits obtained in the two cases. These show the need for careful electrical
design, to prevent large, unbalanced magnetic forces in designing motors and
alternators which are to be supported in self-acting gas bearings.
Water lubricated bearings have been used extensively and successfully in many
applications. Water lubricated pumps are, for example, commonplace in marine
and other service. Water lubricated canned-rotor pumpsj using graphitar bearings
are standard in nuclear submarine service.
Water lubricated bearings are currently being developed for the HTGR helium
circulators, because of their advantages in this application^ as discussed in
Section 8.3 of this report. Figure 8 is a drawing of the test rig for this
machine.
3.15
Liquid metals have also been used as bearing lubricants. Thus, mercury is the
lubricant used in the SNAP 2 rotating unit discussed in Section 7.1 of this
report and in other,earlier units (Ref. 9 ) .
TABLE 6
SOME DESIGN CONSIDERATIONS WITH SELF-ACTING GAS BEARING MACHINES
As part of the current task, tables of the steady-state and dynamic characteristics
of a number of tilting pad bearings (four-pad bearings with L/D = 1/4, 1/2 and 1
and six-pad bearing with L/D - 1/4) have been computed for Reynolds numbers up to
60,000.
In the parallel, experimental work under the task currently in progress, the
steady-state and dynamic load capacity, as well as the stiffness and damping
coefficients are being accurately determined for a four pad, L/D = 1 bearing,
for comparison with the theoretical data. The lubricant used is NaK 78 and
the measurements made are of the journal location and orbit size and geometry,
for a range of values of steady-state loads, dynamic loads and rotational speed.
The characteristics of the test apparatus are:
The test program calls for determination of the bearing steady-state and dynamic
load capacities, as well as of its stiffness and dynamic coefficients, from
measurements of journal center displacement and orbits over the above noted
load and speed ranges. It is the Intent of the test program to obtain this data
on the steady-state and dynamic characteristics of the bearings and rotor-
bearings system with very high accuracy, in order to establish the degree of
confidence with which the existing turbulent lubrication theory may be used to
predict alkali metal lubricated bearing performance. Specifically, any deviations
between calculated and measured performance due, for example, to fluid inertia
effects at pad entrances and exits (velocity lead losses) are to be noted, as
functions of Reynolds number and eccentricity ratio . Tilting pad bearings were
selected because their high stability makes them, prime candidates for high speed,
process-fluid lubricated turbomachinery.
The test rig is mounted in a dry box under nitrogen cover gas. The NaK is fur-
nished from a 50 GPM, 40 psi ffeK circulation assembly and gas purifier. The
assembled dry box and NaK circulation assembly are shown in Figure 9. Figure 10
is a photograph of the test rotor and Figure 11 shows one of the pillow blocks
with the bearing pads and seal rings. Figure 12 shows the assembled test rig,
but with the upper half of the cylindrical housing and the loader bearing omitted.
The hydrostatic loader bearing is shown in Figure 13. The rotor is driven through
a splined quill coupling by a variable speed (0 to 36000RPM) motor and MG set.
There are six measurement planes, one at each end of the two liquid metal lubri-
cated bearing; immediately outside the seal rings, and one at each end of the
nitrogen loader bearing. Two capacitance probes installed at 90 degree to one
another are located in each measurement plane to measure the journal displacements
and orbits.
The tests were started in October 1966 and the unit was operated first in an
"easy to use" fluid for check out purposes and later in NaK. To date, three
test sequences have been conducted in NaK. The first sequence was at speeds
up to 12,000 RPM and was made for the purpose of check out of the system, in-
cluding the NaK circulation loop, feed and drain systems, instrument performance
and others. In the second two sequences, which were conducted at speeds up to
26,400 RPM, test data was obtained first with a balanced shaft and then with an
unbalance of 0.29 oz.in., corresponding (at the highest speed of the test)to a
dynamic load of 180 lbs. per bearing (i.e. 45 psi on projected area). In both
sets of tests the static load was the 57 lb. weight of the rotor, which is
equally supported on the two test bearings. Figure 14 shows photographs of some
of the orbits at different speeds. Currently, preparations are in progress for
conducting other test runs with larger values of steady-state and dynamic loads.
REFERENCES
2. H.G. Elrodj Jr., C.W. Ng and C.H.T. Pan, "A Theory for Turbulent Films and
Its Application to Bearings/' Topical Report under AEC Contract AT(30-1)-3363,
AEC Report No. NYO-3363-2 5 MTI Report 65TR9, March 1965 (also published
under the same title as ASME Paper 66Lubl2, June 1966).
3. J.W. Lund et al: "Rotor Bearing Dynamics Technology," Final Reports under
USAF Contract No. AF 33(615)-1895^ U.S. Air Force Reports AFAPL-TR-65-45
Parts III and V^ MTI Report 64TR14 and 65TR15, May 1965.
5. (a) F.K. Orcutt, "The Steady-State and Dynamic Properties of the Tilting
Pad Bearing in the Laminar and Turbulent Flow Regimes," Topical Report
under NASA Contract NAS-w-1021, MTI Report 65TR32, June 1965 (condensed
version published under the same title, as ASME Paper 66-LUB-19, June 1966).
(b) F.K. Orcutt and C.W. Ng, "Steady-State and Dynamic Properties of the
Floating Ring Bearing in the Laminar and Turbulent Flow Regimes,"
Topical Report under NASA Contract NAS-w-1021, MTI Report 65TR33,
June 1965.
7. P. W. Curwen, G.B. Manning, R.A, Harmon, "A Comparison of Oil and Gas
Lubrication Systems for Closed-Loop Gas-Tarbine Machinery," ASME Paper
accepted for presentation at 12th Annual Gas Turbine Conference and
Products Show, Houstin, Texas, March 1967.
10. M.J. Wallace, "Summary Report of Potential Liquid Metal Bearing Materials for
SNAP 50/SFUR Pumps," Pl€AC-^68, November 1965 (prepared under Contract
AT(30-l)-2789. TID-4500, CAT.UC-25) .
R.Gs Frank (Editor), "Materials for Potassium Lubricated Journal Bearings,
NASA CR 11011, General Electric Co. Progress Report under NASA Contract
NAS 3-2534.
F.K. Orcutt and C.H.T. Pan, "An Experimental Study of Film Rupture in
Journal Bearings with Low Kinematic Viscosity Lubricants," Topical Report
under AEC Contract AT(30-1)-3363, AEC Report No. NYO-3363-4, MTI Report
No. 65TR13, March 1965.
TURBINE
BEAR
Fig. 1 Schematic Diagram of CSN = 2 Lube-Oil and Seal Gas System MTI-1604 oo
(p>-{><H (P>H>^
LOW DP
][H> DETECTOR
EXTERNAL
VI IX}-^ NITROGEN
SUPPLY
Fig. 4 Rotor and Gas Bearing Components for CSG-1 Full Scale Simulator MTI-1625
t.
Fig. 5 Bu-Mines/AEC, 24,000 RPM Gas-Bearing Turbocompressor
00
Fig. 6 Rotor, Housing and Bearing Parts for 150,000 RPM, Two-Stage
Regenerative Compressor
i.
r
MTI-2341
o
Wi
r
Fig. 9 Dry Box and NaK Circulation Loop for Liquid Metal Bearings L!
Test Program (Gas Purifier Not Shown) i^
MTI-2350
F i g . 10 Test Rotor
MTI-2344
1
«•
Fig. 12 Assembled Test Rig (Upper Half of Housing and N^ Loader Bearing
Not Shown)
% .
r
F i g . 13 N„ Loader Bearing CO
U1
MTI-2347
N = 5000 RPM N = 13,000 N = 21,000 RPM
Rotor Balanced Rotor Balanced Rotor Balanced
Scale: 1cm = 0.8 mils Scale: 1cm = 0.8 mils Scale: 1cm = 0.8 mils
SECTION 4
TECHNOLOGY I - FUNDAMENTALS
4,
by
ABSTRACT
I. INTRODUCTION
Conventional fluid film bearings operate in a manner that the viscous shear
stress predominates in the fluid-dynamical process, which obeys the classical
lubrication theory originated by Reynolds (ref. 1). Reynolds' lubrication
theory is predicated on the propositions that a state of laminar flow pre-
vails and that the fluid film thickness is considerably smaller than the
characteristic dimension of the bearing surface. In recent years, trends
toward process fluid lubricated bearings bring to light the need to consider
lubrication films of fluids, which, because of a low kinematic viscosity,
cannot be adequately described by Reynolds' lubrication theory. In table I,
several typical process fluids and their kinematic viscosities are listed
for comparison with SAE No. 10 oil. Depending on the hydrodynamic and the
geometrical parameters in effect, one or more super-laminar phenomena may
dominate the operation of a process fluid lubricated bearing. The purpose
of this work is to discuss the relevance of these phenomena to lubrication
problems, to review the progresses thus far achieved, and to project future
needs of research. The fluid-dynamic point of view will be adopted, there-
for this work is as much directed to fluid seals as it is to fluid film
bearings.
CLASSIFICATION OF SUPER-LAMINAR FLOWS
3U
--J- - 0 (1)
3x.
3
(a^j) (2)
^ 1 at - "j 3xj ^ij 3x
j
(4)
(5)
(6)
for i = 1,3
(7)
4.1.
(8)
^l3^ (9)
^Xj ^""i ^^"2
for i »» 1,3
so that the left hand side of eq. (2), which accounts for the Inertia forces
in the fluid film, can be neglected for i « 1,3 ; and p can be regarded as
independent of X2. In typical oil bearings
^ ^ % 100 (11)
U
f :^ io~3 (12)
LI
and
AT ^ I (13)
Thus the inequalities of eq. (10) are readily justified for an oil film.
B. Flow Instabilities
The significance of the inertia effect is more than the simple prevalence
of transient and/or convective forces. Two types of breakdown can happen to a
laminar flow when the kinematic viscosity is sufficiently small.
One of these is caused by the centrifugal force of a curved layer of
flow. Its classical form exists in the annulus between an inner rotating
cylinder and an outer stationary cylinder [ref. 2]. When the rotational
speed exceeds a critical value, depending on the kinematic viscosity and
the radius ratio, a system of toroidal secondary vortices will appear as
illustrated in Fig. 1. When the radii of the inner and outer cylinders
are nearly equal, the critical speed for secondary vortices is
(14)
and the axial spacing of the vortices is approximately same as the annular
gap AR. At their first occurrence, these vortices are time independent.
Until this critical speed is reached the torque required to rotate the
inner cylinder is directly proportional to the rotational speed. Above
the critical speed, the torque will exceed this condition of being proportion-
al to the speed as shown in Fig. 2. The pertinence of the secondary vortices
in the operation of process fluid journal bearings and shaft seals is self-
evident; here the annular gap AR is the nominal radial clearance C. Secondary
vortices do not occur if the outer cylinder instead of the inner one rotates.
au.
^ . 0 ^i6>
ox.
MTI-2355
^ ir
'\ .rr ^
+ "3 i ^
.- !-^^ --'-^
„^ "i'-'d'^ <^''
where barred quantities are temporal averages and primed quantities are
the fluctuating parts. The last term in eq. (17) represents turbulent
momentum transport which causes the mean stress to be higher than the
laminar viscous stress. The criteria for neglecting the left hand side
of eq. (17) are
1 C > T 1 C ? 1 fia\
and
P V C , 103
y
Since C^ itself is a function of the Reynolds number, above criteria can be
pVC C C C
expressed in terms of , —, r-—, and —. A composite flow regimes map with
DVC C C C
-— as the ordinate, and either —, or 777;fr, or • as the abscissa is shown in
y L VAi K
Fig. 3. The line marking the importance of inertia effects is based on
I- £ _ 1 L_ = 0 25
C^ L C. VAT
C =-^
and for ^— > 2000, C is calculated according to [ref, 4], V is the velocity
f y pVC 1
y
of the moving surface.
Worthy of special notice is the fact that typical process fluid bearings
are designed with — between 10~3 and 3 x 10~3. In this range, the regime of
flow with secondary vortices but without turbulence occupies a very small por-
tion of the map. Also worthy of notice is the fact that the line marking the
importance of inertia effects is almost vertical for — — > 2000. When
inertia effects are to be considered, they include transient, convective,
centrifugal, and Coriolis accelerations.
III. FLOW WITH SECONDARY VORTICES
From the flow regime chart shown in Fig. 3, one would expect that the first
transition from laminar flow to occur in bearings with C/R > 10" would be to a
form of vortex flow rather than to turbulence. It is, therefore, of considerable
practical importance to know as precisely as possible at what speeds transition
to vortex flow will occur. The vortex transition boundary shown in Fig. 3 is
for the classical case of concentric cylinders with no pressure gradients in the
flow. However, in loaded journal bearings, both axial and circumferential pressure
gradients exist. Also, in loaded journal bearings, the radial clearance varies
circumferentially. Obviously, then, the problem of the stability of flows in
journal bearings to development of Taylor vortices is very much more complex than
for the classical case of concentric cylinders.
where W is a mean axial flow velocity determined by dividing the axial volume flow
rate by the annular cross-section area for the flow.
One can note in Fig. 4 that for (N„ ) . - = 0, flow stability increases with
increasing eccentricity ratio. This is due to the net effect of circumferential
pressure flows induced by rotation of the eccentric cylinder. This effect is
discussed below.
DiPrima (Ref. 9) showed theoretically that for the case of concentric rotating
cylinders, a uniform pressure flow in the direction of rotation tends to make the
flow less stable to onset of Taylor vortices while a negative pressure flow tends
to make the flow more stable. When a cylinder rotates within a stationary, non-
concentric outer cylinder, a negative pressure flow is induced in the region of
maximum clearance while a positive pressure flow is induced in the region of
minimum clearance. To determine the net effect of these pressure flows on the
stability of the flow, DiPrima applied his abovementioned analysis to the local
flow at every circumferential point around non-concentric rotating cylinders
(Ref. 10). The well known Sommerfeld (Ref. 11) solution for journal bearings
was used to calculate the magnitude of the local pressure flows.
H In Fig. 6 are shown various experimental data for critical speed for onset
of vortices in flow between non-concentric cylinders. For these data, the point
of onset of vortices was determined by torque measurements with silicone fluids
as the test fluids (Refs. 12 and 13). As can be seen, the measured transition
speeds for onset of vortices all are higher than predicted by the minimum transi-
tion speed curve of DiPrima. Also, one can note that transition speeds become
greater as clearance ratio increases, particularly at large values of e.
For the typical clearance ratios found in bearings, the range of operating
conditions in which vortex flow will occur without turbulence is fairly narrow
(see Fig. 3). As we can see later, once fully developed turbulence sets in, the
transport mechanism in the flow will soon be dominated by turbulent fluctuations
and the effect of the secondary vortex flow which is present soon becomes negligi-
ble. Nonetheless, it is still of practical interest to develop a "vortex theory
of lubrication'" which can essentially "bridge the gap" between the operating range
in which laminar theory applies and the range in which turbulent theory applies.
Such a vortex theory of lubrication has been developed (Ref. 13) based on an
analysis by DiPrima (Ref. 14). In his analysis, DiPrima developed the following
theoretical relationships for flow between concentric cylinders in the vortex
regime at speeds just above Taylor critical speed.
V T (Q)
-^ = 1/2 + Q/6 + 1 --S- F(Q) G(Q) (21)
Vi
T (Q)
= 1 - Q + 1 -- ^ H(Q) (22)
pV] /h
For T > T
c
where
Q 3 V /V = 3 ^ — —
•^ p^ c -^ [ 12yRi 36 Vi
V
p mean circumferential velocity due to pressure gradient
V
c mean circumferential velocity due to rotation = Vi/2
Vl
surface velocity of inner cylinder
Rl
radius of inner cylinder
R2
radius of outer cylinder
V
fluid viscosity
Vih
T local Taylor number = 4
Rl + R2
h local radial clearance between inner and outer cylinder
« C (1 + e cos 8)
Relations (21) and (22) are derived for the case where both clearance and
pressure gradient do not vary circumferentially. In the case of non-concentric
rotating cylinders, both clearance and pressure gradient do vary around the
cylinders. Q, therefore, varies circumferentially. To obtain an approximate
calculation of the effect of vortex motion on shear stress and pressure gradient
in the flow between non-concentric cylinders, one can apply relations (21) and
(22) locally at each circumferential point using the appropriate local values
of Q and sum these local contributions. The procedure for doing this is describ-
ed in Ref. 13. From this procedure, one can calculate the effect of vortex motion
on the overall viscous drag associated with rotating non-concentric cylinders.
Typical results for calculated viscous drag are shown in Fig. 7 along with
some experimental measurements taken from Ref. 13. The ordinate used in this
figure is G/G , the ratio of rotational torque in the vortex regime to the torque
4.1.11
that would occur if flow remained laminar. The solid curve shown in Fig. 7 is the
theoretical curve calculated using values of critical Taylor number, T (Q), which were
determined from DiPrima's analysis for C/R^ -> 0 (see Fig, 5 ) . However, as was shown
earlier, the critical speed for first onset of vortices was found experimentally to
increase with C/Rj. To allow for this, the theoretical curve for C/R^ = 0 is shifted
horizontally to the right in Fig. 7 so as to align the point of calculated increase
in torque with the measured transition point for onset of vortices. Two such shifted
curves are shown: one for C/Ri = 0.0104 (dashed curve 1) and one for C/Rj = 0.099
(dashed curve 2). For C/R^ = 0.0104, agreement between measured values of torque
and the shifted theoretical curve are fairly good at speeds just above critical
speed. For the case of C/R-[ = 0.099, measured values of torque increase more rapid-
ly at speed above Taylor transition speed than would be predicted by the corresponding
shifted theoretical cuinre.
Once the inertia effects are separated out from pressure profiles measured in
the vortex or turbulent flow regimes, one can assess, directly, the influence of
vortex motion or turbulence on the magnitude and shape of the profile. In Fig. 9
are plotted three experimental profiles measured, respectively, under conditions
of laminar flow, vortex flow, and vortex flow with perhaps the beginning of turbu-
lence. Inertia effects have been separated out from the profiles shown. The
ordinate is the dimensionless pressure PC^/6pViRi. If the flow remained laminar, all
of the profiles would follow the same curve, since this dimensionless pressure
correctly accounts for the effect of speed on the laminar profile. The fact that
the profiles measured at /i/2T* - 83 and /TJTF = 169 have larger magnitude than
that measured at /1/2T = 13.6 is a direct indication of the effect of vortex mo-
tion or turbulence. One can note that at i/l/2T = 169, superlaminar effects have
increased the peak magnitude of the experimental profile by approximately 40%.
The development of secondary vortex flow has a more pronounced effect on the
magnitude of pressure profiles than upon their shape. In Fig. 10 is plotted the
ratio Pmax / ( PLTmax
) V S Reynolds
•' number for non-concentric cylinders
•' at various
eccentricity ratios. Here P is the peak magnitude of the measured profile
while (P.f ) is the theoretical peak pressure that would be obtained if flow
remained laminar. The point of onset of vortices in the flow is characterized by
the ratio of Pmax to (P^) L max becoming greater than unity. Theoretical curves
based on both vortex theory and the turbulent theory of Ng and Pan (Ref. 4) are
compared with the experimental data in Fig. 10. Data shown are for C/Ri = 0.0104.
The significant thing to note in this figure is that for N ** > 2400, the turbulent
theory begins to provide a quite accurate calculation of the peak pressure, even
though it is known from visual observations that a vortex secondary flow still per-
sists in the flow. Apparently, as soon as turbulence becomes fully developed in the
flow, the turbulent transport mechanism dominates over that due to the vortex flow.
Re V ^
4.1.13
The effect of vortex flow and turbulence on the shape of the pressure profiles
around non-concentric cylinders is shown in Fig. 11 where various experimental pro-
files are plotted normalized to a peak value of 1.0. As noted earlier, flow at the con-
dition /T/XI - 169 may have contained some turbulence as well as vortices, since the
Reynolds number corresponding to the flow condition was N « 1660. The turbulent
theory of Ng and Pan (Ref. 4) does predict the shape of the provile at •/TfTf = 169
more accurately than does vortex theory. However, referring back to Fig. 10 we see
that turbulent theory does not predict the magnitude of the peak pressure at this
condition (Np « 1660) as accurately as does vortex theory. Essentially, at this
condition, neither vortex theory alone nor turbulent theory alone can be expected
to adequately describe the flow.
For the case in which an outer cylinder rotates about a stationary inner
cylinder, vortex flow will not develop. In this situation, the friction factor,
A, for the flow will be a function only of the Reynolds number. A is defined as
X = 2T/PVI^ where T is the shear stress on the surface of the rotating cylinder.
In the case of rotating, inner cylinder, where vortices will develop, the
friction factor becomes a function of both Reynolds number and clearance ratio.
The dependence of A on these parameters is expressed quite accurately by the
following empirical formulas developed by Wendt* (Ref. 15),
(R2 - Rl) R2
X «
*The range of Reynolds numbers over which these formulae are to be used has
been modified slightly by the writer.
The dependence of friction factor on clearance ratio can be seen directly in
Fig. 13. Here are plotted curves of A vs Re for C/Ri = 0.01, 0.002, 0.001 and
0.0005, calculated from Wendt's empirical formulae. The Taylor critical speeds
for predicted onset of vortices at different clearance ratios are indicated along
the line for laminar friction factor. Also shown is the friction factor curve based
on the turbulent theory of Ng and Pan. This latter theory neglects effects of
clearance ratio since it neglects consideration of vortex flow.
It should be noted that Wendt's formulae are based on data taken with test
cylinders having clearance ratios greater than 0.095. Application of Wendt's
formulae to clearance ratios on the order of 10~3 is therefore questionable.
Hence, the discrepancy at high Reynolds numbers between the turbulent theory
curve and the curves calculated from Wendt's formulae may not necessarily in-
dicate an inadequacy of the turbulent theory.
4.1.15
S (23)
dx. 'ij dx, + 2 ^ ^ + ^ (Ph) = 0
G.. is a diagonal matrix (G.. = 0 when i 5^ j) when one of the coordinates (say,
•^ J -^ J
x,) is chosen to be parallel to the surface velocity (V„ = 0 ) . In particular,
1
G.. = TTT 5.. for the laminar film, where 6.. is the Kronecker delta. For tur-
ij 12 ij Vih h3 Sp iJ h3 Sp
bulent films, G. . depends on -77-, -TT T^—•, and —- -^ . Results computed in
' ij '^ V ' |j,V dx^' laV dx ^
ref. 19 for G . are reproduced in Figs. 14 through 17. In Figs. 14 and 15 these are
iJ Vih !h3 5p . , op
respectively plotted against —-— for various values of —- -5— with ^ = 0;
f -^ ^ °„ V I [j,V dx, dx^
in Figs. 16 and 17, — -^ is the variable with -5^ = 0. When neither
laV d x dx^
dn dp
^ nor ^ is zero, the following interpolation formula was found to be
dx-^ dx^
satisfactory:
V^h h Sp
G. . h ^P K
[iV Sj
29
G. .
/V^h
h Bp
a^ + b^, 01+
20 V (24)
'' M-V 3xj^ 3t IJ |j.V d x „
where
-1 b/a
9 tan
is the angle between the direction of the pressure gradient and the surface
velocity.
From Figs. 14 through 17, it is seen the effects of pressure gradients are
small so long as
V .2
,^, ^ 1 ^ 1 < 0.25 X 10-2
Vh jiv 3x'
In this range
which was treated in ref. 4 and the results are reproduced in Fig. 18. Using the
approximation of eq. (25) in eq. (23), one clearly sees that the equation becomes
linear in p, this is referred to as the "linearized turbulent lubrication theory."
It is quite evident that the linearized turbulent lubrication theory would be quite
accurate for nearly constant film thickness conditions. It is rather surprising that
3D
the non-linear effects of -r^— seem to remain small for a full cylindrical journal
bearing even when the eccentricity is very high as shown in Fig. 19. It appears
that the linearized turbulent lubrication theory should be adequate for all smooth-
film self-acting bearings.
The linearized turbulent lubrication theory has also been used to calculate
bearing stiffness and damping which are needed in the prediction of dynamic rotor
response. Such a prediction was corroborated by test for a rotor supported by
four-shoe tilting-pad journal bearings. Correlation between test and theory for
this bearing is shown in Fig. 20 for both static and dynamic conditions. The agree-
ment is quite gratifying.
V. INERTIA EFFECTS
Experimental evidence of the convective inertia in the vortex flow regime was
indicated in Fig. 8 for a full cylinder. The over-all consequence is primarily a
shift of the center of pressure such that the attitude angle would be a little
larger than 90° and, to a lesser extent, the bearing force is somewhat increased.
It should be recalled that the clearance-to-radius ratio is considerably higher
than common lubrication practice.
An inertia effect is also manifested along the edge at which flow enters the
fluid film. Some data on this phenomenon was reported in ref. 16. Wall pressure
and velocity profile at the inlet edge of a cylindrical partial-arc pad were
measured for several combinations of inlet film thickness, rotational speed, and
inlet flow. This data was correlated by the formula
p = p + C ^ (26)
*^o ^a e 2
where
pU 2
p = p +
*^o "^e
well with the linearized turbulent lubrication theory as illustrated in Fig. 24,
However, discrepancies suggesting a ram effect are still evident in both the load-
deflection and attitude angle curves, for instance, for C/R = 4 x 10-3 Re ^ 5820
and e = 0.27. Indications of increased disparity between the linearized turbulent
lubrication theory and test data for a larger C/R are also apparent in the dynamic
orbits shown in Fig. 25. It is somewhat surprising that similar discrepancy was
not apparent in the data of ref. 21. This may in part be due to a smaller clear-
ance-to-radius ratio, and also in part due to the close spacing between the con-
secutive pads.
In ref, 22, the inertia effect at an internal edge of a step pad was also
studied experimentally. It was found that a steep pressure drop occurs here. This
is probably due to the sudden acceleration of flow as it enters the step. For the
particular conditions studied, which include a fairly large gap, the land region of
the step pad practically does not support any load.
VI. CONCLUSIONS
The experience cited above can be examined in terms of the flow-regime map.
The mean
aean film thickness
thickness, C , is used to calculate the parameters which form the
coordinates of the map. C mean is defined as follows: *^
mean
Bearing Type C
mean
Full Cylindrical Mean Radial Clearance
Pocket Bearing
Land Region Mean Radial Clearance
Pocket Region Mean Radial Clearance plus Pocket Depth
Partial Arc Mean Radial Clearance
Pivoted Pad Concentric Clearance at Pivot
NOMENCLATURE
C^ coefficient of friction x/ -^ p V ^
r z m
C mean film thickness
mean
D diameter
G rotational torque
<«Re>axial V^^
p pressure
P max P in film
max
Q 3V /V
^ P c
Re VC/v
i2
S Sommerfeld nximber - , L^
NyLD
t time
4.1.21
Vl h 2
t e m p e r a t u r e , l o c a l Taylor ntimber 4 —^
V Rl + R2
load
Cartesian coordinates
temporal average of ( )
fluctuating component of ( )
arc angle
Kronecker delta
R2 - Rl
eccentricity ratio
kinematic viscosity
fluid density
attitude angle
4.1.22
•I
Subscripts
2 normal direction
REFERENCES
Reynolds, 0., "On the Theory of Lubrication and Its Application to Mr. Beauchamp
Tower's Experiments," Phil. Trans. Roy. S o c , London, vol. 177, part 1, 1886.
Ng, C.W. and Pan, C.H.T., "A Linearized Turbulent Lubrication Theory," Journal
of Basic Engineering, Trans., A.S.M.E., Series D, Vol. 87, 1965, pp. 675-688.
Krueger, E.R., and DiPrima, R.C., "On the Stability of Spiral Flow between
Rotating Cylinders," Journal of Fluid Mechanics, vol. 19, part 4, pp. 528-538,
1964.
DiPrima, R.C., "The Stability of a Viscous Fluid between Rotating Cylinders with
an Axial Flow," Journal of Fluid Mechanics, Vol. 9, Part 4, pp. 621-631, 1960.
Donnelly, R.J., and Fultz, D., "Experiments on the Stability of Spiral Flow
between Rotating Cylinders," Proc Nat. Acad., Sci. 4i6, 1960
Snyder, H.A., "Experiments on the Stability of Spiral Flow at Low Axial Reynolds
Numbers, "Proceedings of the Royal Society of London, Series A, Vol. 265, 1962,
pp. 198-214.
DiPrima, R.C., "The Stability of a Viscous Flow between Rotating Cylinders with
a Pressure Gradient Acting Round the Cylinders," J. Fluid Mech. 6^, 462-468 (1959).
DiPrima, R.C., "A Note on the Stability of Flow in Loaded Journal Bearings"
ASLE Transactions, 6, pp. 249-253 (1963).
Chou, Y.T., and Saibel, E.,M "The Effect of Turbulence on Slider Bearing Lubri-
cation," Trans. A.S.M.E., Ser. E, vol. 81, March, 1959.
Elrod, Jr., H.G., and Ng, C.W., "A Theory for Turbulent Fluid Films and Its
Application to Bearings," ASME Paper No. 66-Lub-12.
Burton, R.A.,and Carper, H.J., "An Experimental Study of Annular Flows with
Applications in Turbulent Film Lubrication," ASME Paper Np. 66'-Lub-14.
Orcutt, F.K., and Arwas, E.B., "The Steady-State and Dynamic Characteristics
of a Full Circular Bearing and a Partial Arc Bearing in the Laminar and
Turbulent Flow Regimes," ASME Paper No. 66-LubS-4.
TABLE I
TABLE OF FUNCTIONS OF Q
llOr
lOOi
90=
5 80'
€=0
*2
Rj = Iil2 IN.
R2 = 2.00 m
i
UJ
MTI-2313
4.1.29
O
>
MTI-2305
4.1.30
0.
10 20 30 40 50 60 70 80
^NRi) AXIAL
MTI-2314
4.1.31
JOURNAL CENTER
ECCENTRICITY, e
Figure 5 - The Taylor Vortex Transition Value of V C jC as a Function of e for Various Values of
4.1.32
100
^ C/R, = 0.099 1
•
o C/Rp 0.0104 r
90 —^THEORY
i 1ii_\/i« 1
C/Ri
\jf i « i
-^ 0\j 1
(DiPRIMA, REF.IO) 1
1
1
80 1
1 s
:
V
c^70 /
1
7
/o
Qi
/ y 0
/ /
^ 0 /
o 0
>" 60 i /
/
/
/ /
» /
> /
50 /
/
9/ 4>
1^ O
y" /
y i /
/
.^ #w
Ox
.^<T
y^ /
X**" •
, > • * )^^^
40 t? N
5
•
/
• , ^ ^
y
30
0 0.2 04 0.6 0.8 1.0
ECCENTRICITY RATIO, €
Figure 6 - Measurements of Transition Speed for Onset of Vortices
in Flow Between Eccentric Cylinders
MTI-2304
1.8
O RUN
A KUN Z
C/R, = 0.0104
o RUN 3
1.6 €=0.51 „ V RUN 4
C/R, =0.099
6 = 0.50 ^
.X*
14- ^ - ' ^
NOTE•• CURVES 1 AND 2 ^ "
ARE THEORETICAL CURVE SHIFTED VORTEX THEORY, C / R , — 0 ^ ^ ^
^ - - ^
FOR COiPARISON WITH iEASURED "
DATA ^ ^ o
^ 1.2 1 I i . E 1 i
o
^ 7*. e^ •
0 __ m
^ ^ .—^^ ^ 4
THEORETICAL TR•ANSITION iEASURED ' iEASURED TRANSITION
r\r\ 1 HIT- r\ i r\
ruiiMi , t / r TRANSITION
0.6
0 10 20 30 40 50 60 70 80 90
2C
V,C/i
y R1+R2
1.8
>/" "X .CURVE 1
1.4 I ^y^N\^A/CURVE 2
f
/ 1
r
®
1.0 / 1
f
> u
0.6 4/ \\
\\
^/
0. / u
V
—f-1
>
o> M
CD """"**"- ^ i ' ^
^ \ ^ /
•0.
€ = Q68,V T T = I69 , NRe=l660'll
^
JL- ^
f
/
//
•0, o PROFILE WITHOUT u y
INERTIA EFFECT ll A
^
^ PROFILE INCLUDING 11
INERTIA EFFECT /^
^
A INERTIA CORRECTION 11 \
\ \
J
{CURVE 2 - C U R V E 1) y/
\
1 1 1 1 1 1 . ^
1.8,
0 40 80 YtJO \i50 ^ ^ >2' 10 2«JO 3J?0 36 0
ANGULAR POSITION, DEGREES
f
//
9 ^ V - "LAilNAR THEORY (SOiiERFELD^
//
0.8 A
/.w
ir——»
kit
4 / \
/
/
.^ .(tmV
04 ^
^yf
cr. a^
^
0
to EXPERliENTAL
1 K L^
X ^y-
#^
A Y P = 69,NRe=l660
i
y
/
/
-1.6
0 40 80 120 160 200 A 240 280 320 360
ANGULAR POSITION, DEGREES
4>
30
• €=0.68
o €=051
EXPERliENTAL
2.6 A c=0.35
A e=0.20
2.2
.r
TURBULENT THEORY, €=0.68
o
^ E 1.8
TAYLOR VORTEX
S 14 THEORY, €=0.50
0.6
0.35/ >¥
CRITICAL REYNOLDS
NUiBERS
02
200 800 1600 2400 3200
REYNOLDS NUMBER, NRE = V|C/F
MTl-2306
1.0
TURBULENT THEORY EXPERliENTAL
^
0.8 ^ L i / i T = l69 • .^1660
ND. .
VORTIXTHEORY;^ \V|TM3.6O iRg=l34
o
E 0.6 ^
Q_
^ \
Q_ 0.4
cc
• ^ l A i l N A R THEORY
k
^ \
0.2 ^F =
c/>
CO
LiJ
IE
0^
Q_
Q \ ^ ^
yj -
•0.2
—J
f = 0.68
« ^ - 0.4 \
^
o _ 0.6 \ d
/y'-
0
0.8 \
\ /%
^
H.O,
0 40 80 120 160 200 240 280 320 360
ANGULAR POSITION, DEGREES
Tfe» t
%
EXPERliENTAL
o C/R| =0.099
\
WkJ
p^
C ^
• C/R =00104 (CORRECTED)
!^QaOA#%
•"'-*='«*4^
0 ^ 3
0.001
100 1000 10,000
ViC
REYNOLDS' NUiBER NR^
.^ X
\ X. .•SilTH AND FULLER
\
X K C/R| = 0.00293
V -" ^ ^
" ^ w
! 1 1
r--" X,
\ .
^
' ^ - - . . ^ C/R, = 0.01
^
I :RIT CAL
^ •'*-x__'''''*'*'*^'-oi_J
^, 1 ^''-^^^..^^^^^ r-^ : ^
SPEEDS C/R|=0.()005 -'-^
II
\ —
^^ 1 ^*-^
. ; : =
>
^/Ri = 0.0005-
^""
€=0
0.0001 »»«j
1x10"
—
^
- ^
— ^ ^
_.
^ \
1x10^ - ^
)^
[iV Ox
V
^ ^
—
o
0.01
X
N
TxLO^
"~-
lxl08__
0.1
1x10"^ p
4 p
1x10
^**'*'*^
^ \
llxlO^
\
X.
^
O
0.01
^
3
^p-Lix >n6
N^
P -t-. J
1x10'
^
ramiisr**®*""
X
P
, .8
1x10 —
1 L_ MM 1 1 ™UJ_L 1 1 „JJJJ.
10 10- 10 10'
REYNOLDS NUMBER = —
MXI-2295
4.1.42
0.1 _ _ _ ^ _ _
-
-
:24
^
^^^k.
Si^
0.01
XV
^
o
— 1x10^
VX
——.^JV
IxlO''
\
———.^^
5x10^
REYNOLDS NUMBER = —
I
0.1
__
— ,
Z^Io^
LxlO^
"™"
h^ Sp 1x10^
=
3
V.
\
O
0.01
__.
1x10^
X
" ^
\
.
1 xlO'7
__ • ^
5x10'' .
MTI-2299
^11' ^33
o o
I
o 11 1.11^ 11
00 / 7
^
00
\
o
1 // /^
o h- /
p
a. o
h-'
O < <
J/ /
\
o
p
<
p
/ .
V
/
^
1/V / /
y 7
t.
5x10
MTI-2303
46
02 04 06 08
ECCENTRCTY RJTIO igi ml ECC'IJTRCTY RATO ..11 ml
Comparison between measured a n d calcuSoted s t e a d y state Comparison between measured a n d calculated steady-state load
load capacity — m = 0 city - m = 0 5
"
u ^ --.» ^ CEiTERLliE
BEARIN6
^ / e
^
3£ r ^ o
-«^ CENTERLINE
3 1 / o
SHAFT
1
1 -j—i /
'A
« 1f
s THEORETICAL CURVES AND
EXPERIMENTAL DATA POINTS
Cleardiueiatu (Cj,//0 > X 1<> ' m/ii
om 05 W 200LB
SleiideineNsiatio {I ID) - 1 (I
» m 0 W 200 LB
iVic length (8) 80deg
- - » m 0, W ZOLB
PlV( t pOMtl 111 (9p/(3) ~ OOT
0 I f 0 12
SHAFT SPEED RPM « l O '
r
OO
--
y /
J
10
3700
10
""
—
y J4
/
/
/
/
11,100 ^ A /
_.
— L/ /
r-
1375
/
- _ ™ ^
/
1.0 1.0
0.4 0,5 0.6 0.7 08. 0.9 1.0 0.4 0.5 0.6 0.7 0.8 0.9 1.0
Eccentricity Ratio Eccentricity Ratio
C/R = 3.9 X lO"^ C/R = 2.9 X lO"-^
Fig. 21 60-Degree Partial Arc Bearing -
Load Versus Eccentricity
90^
Linearized
80*^ Turbulent
70' 0-2 Theory
— - Re = 1375
60^ Re = 14800
0.4
EXPERIMENTAL DATA
0.S5 C.8.0IL
50^ L/D = i/a p-so®
\ ">< /Ml / / — ®
I
20'
40
LB LOAD
2.9 » 10 *
0.8 ® '^°^
40^ + 20 "1 LB LOAD
t @0 C/R=
1 0 * I20J 3.9 t 10"^
o,^^\Ki:^
N--^^ \^^
TO
C
1-1
(0
ho
O
O
H
rt
H-
OJ
I—"
S>
M
O
to
*J^ ; C ? ^ |
* / A
y
"I^ /y
^ V ^ __>w^
^^^^-^ ^2L_I
%f
THEORETICAL RESULTS > - " ••--^^ 7"^ ' — _
I3,3M
^ \ Xl665
^
a.
r —Vi
»|z . ^^ - \V ^ *<r
^ . ^ ^ « .^ ——
^^ —
^^^^ _^^^^
.4^ >-*^
/ Re C/R
/
/ j^
..ai^ X 1665 2 « 10 5
• 5820 2 « 10''
/ 0 5820 4 « lO'
/ A 13304 4 « lO'
^ f
01 02 03 04 05 06 07 0£ 09 10
ECCENTRICITY RATIO, «
10^ _ — . ~-™i ^ ™„
. .^^^
-y
-- - /
— ^ - __^_,
— - -r 1 "~^"~^
-4__- .„^^„__--^
J — ^
:
\ — .
^
'? •;z=^
BIZ 1_ ^ „u ^
' ^ lOO* PARTIAL ARC
a ^ o C/R 2 «10 5 INCH/INCH
^^,„.H^"*^^!NAR FLOW 50 cs OIL EXPERIMENTAL
OATA
^
03 04 05 06 07 08 09 10
ECCENTRICirr RATIO, ,
>r Y
1 nA
.24
•i ^ 1 111 1 16
S =0.120 S=0.046 .20-
S = 0.0906 14
^ "\ s 14 *-<^
.20
-—THEORETICAL THEORETICAL .16 THEORETICAL
«^^' > — THEORETICAL •I2
--' __ X • ^
EXPERIMENTAL 16
EXPERIMENTAL EXPERIMENTAL / ./ ^^- - .-- "" r-r- ^ k EXPERIMENTAL ....--' .— ' T ^ •^.
12 s
— • .12- / / ^4'^ r^ .X" 08 v
^ A' \
j»=. M I' /' ^M
h^
.08 • " " ^ 06 06
^ ^ ' ^/
/
c
•a" .04 > \
/-r
^0 04 A> 04 ^
1
A /^ / /'
1
•1^
/ /^
1 04^
/ V _ }i f 4/
/ 02
i //
02-
// „
-2 8 -24 -2 0 -16
1
1 I-.0 8 -.04
.P^y 1.1 5 .20 .24 .28 "-.28 -24 720 -.16 i / ' / 0 8 -04 M ^ St
1
15 20 .24 .28" X-
-2
° /7- - 4 -12 -10 -08 -06 -04 -02 02 04 06 08 10 12 u
14 / • ' 18 20"
V
"-22 - 2 0 / ^4 -16 -14 -12 -10 -08 -06 -04 -02 02 04 06 08 10 12 14 //{% 2o"
--«o . - s
•^
-.08
^
/ INCHES X10''
i.^ 0 ^
-.08
INCHES X10"' // -04
<^. - y
/. / INCHES X I C ' //
/ ^
-04 i INCHES X W '
• - _ / ^ r
r
-06
-.12- -.12 ^
.M.. 1 -^
-.16 -.16- \ x? ^^ . / " -^ -+•;/ ^
-.20-
\ '^~ • - .
_- —
—o— — "f ^ \ -J .,- -•\l ^ --
\N y \ ,^
_0>i
V, ..«-^
X ' --> ~r -14
^-
Y " " INCHES X10'' ^ INCHES X 10'' T, -16 -^-
r" Y NCHES X 10"' Y NCHES X I C
Y 'r j 1 1 1 1
14
S = 0 0372
.20 -—THEORETICAL 16
S =0.0338 THEORETICAL 12 i„™ a— Re =8314
EXPERIMENTAL EXPERIMENTAL 10
—-"- - r — ^•s
— THEORETiCAL .16 .12- ^ \ C/R=2x 10"^ IN /m
EXPERIMENTAL
.12^
— -- - 0 ~ ^ V
.08- i - "' u • —
*——~.~~\-^
1
y
^' ^
1 l\
B.S'"*' f
i ^ 04
V
^ .-**" < 1 .rf7
*>'**''* <^ /»
.*^ ^ ^ -f ^ --^ 02- i /
x.^4
-.2 8 -14 -.20x^6 -.12 / -.04 k /
^'
^l^
../ /'
^''^16 .20 .24 .28
1 ) 1
- V
-20 -16 ^
•
/m»
,^.0 4
1
-rtft
r^"
** .12 .1B .20 .24
1 1
INCHES X10"'
yy
"-24 -22 - 2 0 . / ' ' -16 / - a
/ 1/
y ^
-10 0 8 - 0 6 -04 -02 02
"
04 06 08
^PTJ-f
y y
'O^yf^U^.' 18 20 22
INCHES X ID"'
Y INCHES X I C '
C/R = 4x 10"^ IN./IN.
MTI-2298
4.1.51
10^
Actual and
Contemplated Bearing Experimental
10== Designs Pump Bearings
Test Refer-
Fluid Remarks ence No.
\
\ _@®LAN. 1. 3 KW HG RANKINE Hydrostatic Pads 6. Partial Arc 60° Pad Silicone Data suggests 23
Oil inlet ram effect
\ 2. SNAP 8 4 Shoe Tilting Pad 7. Partial Arc 100° Pad
' ^ ^ ' ^ ^ ^ ^ ' ' ^ Good agreement 25
with linearized
30% pad loading (a) C/R = 2 x 10"^ Silicone
turbulent lubri-
Oil
cation theory
(b) C/R = 4 X 10~3 Silicone Data suggests
3. General Atomics Oil inlet ram effect
(a) Lower Journal Hybrid Pocket
(b) Upper Journal 8. 4 Shoe Tilting Pad
(a) Unloaded Silicone Good agreement 21
4. Liquid Hydrogen Pump Oil with linearized
Hybrid Pocket
(b) 50% Pad Loading Silicone turbulent lubri-
cation theory
Oil
5. Liquid Oxygen Pump Hybrid Pocket
10 ^ 10"
Figure 26 - Flow Regimes of Various Tests and Design
MTI-2296
4.
by
Ralph Burton
Southwest Research I n s t i t u t e
San A n t o n i o , T e x a s
ABSTRACT
INTRODUCTION
This apparatus was based upon a 6-ft diameter wheel, 20 in .long. About this
surface, which could be turned at speeds above 100 rpm, various stationary bear-
ing configurations could be supported, including tilted pads, stepped pads,
spiral grooves and the full bearing. End leakage could be restricted by plastic
curtains, thus making possible the simulation of long or short bearings, depend-
ing upon the presence or absence of the curtains. Reynolds number in excess
of hu /2v = R = 2500 were possible, where h is film thickness, u is journal
surface speed and V is kinematic viscosity. This is well within the turbulent
regime.
Figure 2 shows pressure variations for two tilt angles, and end leakage fully
restricted. Again the leading edge pressure plays a dominant role.
When entrance region velocity measurements were made for different leading edge
pressures in tilted pads and stepped pads, it was found that the correlation
presented in Figure 3 prevailed. Here the dimensionless pressure is related
to the velocity in such a way that results can be correlated in terms of a
dimensionless "impact" pressure according to the equation shown. If the entering
flow had been carried at wheel speed (u , where u = 2u, ) then brought to a halt
S S DC
in a frictionless process, as on a pitot impact tube, the following relationship
would prevail:
4.
2 2 2
p/pu, = u /2u, (1)
be be
or
2
2 + - V = 2
o 2 (2)
pu, 2u,
be be
Stepped Pads
Pressure distributions for stepped pads are shown in Figures 4 and 5. The
"unshielded" case represents unrestricted end leakage and the "shielded" one
represents zero end leakage. Of special interest is the large pressure drop
at the step where film thickness changes from 0.5 in. to 0.25 in. For the pad
with leakage the pressure drops to a sizeable negative value at this point. This
2
dimensionless pressure would correspond to -14.7 psig when (14.7)(144) = pu, (0,1)
or pu, - 21,200. For a fluid with the density of water this would occur at
be •'
approximately u = 200 ft/see. This is not an excessive speed for projected high-
s
performance applications. Thusj it would not be surprising to see cavitation under
a pad in the turbulent operating region.
Spiral Grooves
Figure 6 shows pressure tap loeatioas for a spiral groove configuration in the
large-scale appratus. Figure 7 gives pressure measurements corresponding to |
the tap locations indicated in Figure 6. The spiral angle is 45 , and the ratio
of groove-clearance to land clearance is 2:1. Note that negative dimensionless
pressures occur over the land and at one station in the groove. The maximum
negative magnitude is very near that measured for the stepped pad.
Short Bearing
Conclusion
The above examples serve to illustrate the rather sizeable interial effects in
turbulent lubricant film. These are: (1) leading edge "ram effect" which moves
center of pressure forward and otherwise significantly affects pressure distri-
bution under pads, (2) convergence effect at a step which causes a strong pres-
sure drop and may lead to cavitation, (3) negative pressure over a land of a
spiral groove configuration which can, likewise, lead to cavitation, and (4)
tendency of pressure to be high in thick-film regions and low in thin-film
regions (as in the Bernoulli effect), thus giving rise to a possible negative-
spring effect in turbulent, short, full bearings.
References
(1) Carper, H.J., and Burton, R.A. "A Technique for Micromanometry," AIAA
Journal, Vol. 4, No. 3, 1966.
(2) Burton, R.A., and Carper, H.J., "An Experimental Study of the Annular
Flows with Applications in Turbulent Film Lubrication," ASME Paper 66-Lub-14.
.24
LIMITED END LEAKAGE
A
P .12
X
£
Fig. 2 Pressure variation for two tilt angles, for a pad without
end leakage
4.2.7
o
o
.3 o\
—
.2 _. _
X O
o \^
P 1 OBB \ . «•
«.j2 \y
P%c o\ 0
0 _
\D
—
P
-.1 —
Ho \
-.2
!
-J L 1 f
0 .1 .2_ .3 .4 .5
U2
2Ub2,
Fig. 3 Variation of pad-entry pressure with inlet velocity
MII-2111
1 1 r—
UNSHIELDED
0.540 -0.280
A R»2400
O R«f200
SHIELDED -
L.E. T.E.
.6 0.540 0.265"
- A 0.540 0.280 _
R«2400
' .4 — -
Y-"^
K
• P \ .^y^
-
I -^
I " b .2 4 .6 .8 1.0
• MTI-2113
4.2.10
•I
i
t
i
6"K
Fig. 6 Layout of experimental groove and land arrangements,
showing location of pressure taps MTI-2114
SEALED
0.021 OPEN
-0.043
-0.059
-0.091
-0.080
-0.080
-0.027 1/4"
-0.043 CLEARANCE
0.289
-0.016
0.289
0.257
0.278
0.203
0.I3I
1/2"
CLEARANCE
^ O R = 2400 O x
^
""yC^\$4^r7C\JKW
Wf ^=0^
^
7?^Xg
^
\
\
\M\yy^-^
\
\
\
\
\
>
r X x
^
::2^-^w^
x^^^^ \ J 2 ! z\^rjj\^
::3j\\^/ y^
C///
X//
1
x< 9 QO /
/ 1
^
/
> / ^
i \ '
0^
•I
4.3.1
by
ABSTRACT
DISCUSSION
The test apparatus for the visual study consisted basically of a stationary
glass outer cylinder having a 4 inch internal diameter with an aluminum, rotating
inner cylinder, 3.64 inches in diameter. The inner cylinder could be positioned
at arbitrary eccentricites within the outer cylinder by means of eccentric bush-
ings at the ends of the cylinders. A sketch of the test apparatus is shown in
Fig. 1.
The test fluids used in the experiments were silicone fluids of various
viscosities. Flow visualization was accomplished by means of fine aluminum part-
icles, suspended in the flow.
The movie presented contained scenes showing both Taylor vortex flow and
turbulent flow. Taylor vortices are a form of flow instability which arise due
to centrifugal forces. The instability takes the form of toroidal vortices spaced
periodically in the axial direction. The characteristic flow pattern in a pair
of Taylor vortices is shown in Fig. 2. In Fig. 3 is shown a photograph of non-
concentric rotating cylinders with fully developed vortex flow. Vortices are
absent from the ends of the cylinder due to the axial pressure flow there which
stabilizes the flow against development of vortices.
When the cylinders are concentric, the vortex activity in the flow is
symmetrical around the cylinders. When the cylinders are eccentric, however, the
intensity of vortex activity varies considerably around the cylinder although the
axial wavelengths of the vortices remain constant. This circumferential varia-
tion in vortex intensity is shown in Fig. 4 where a set of 18 photographs taken
at 20° intervals around the test rig are assembled into a single "unwrapped" view
of the flow around the cylinders. The eccentricity ratio for this figure is
0.475. The widest part of the annulus between the cylinders is at 6 = 0° and the
narrowest portion is at 9 = 180°. Rotation of the inner cylinder is in the dir-
ection of increasing 6. The dark vertical lines at 6 = 107°, 227° and 347° in
Fig. 4 are the tie rods connecting the upper and lower flanges of the test rig.
When the cylinders of the test rig are positioned at fairly large eccentri-
city ratios (e _> 0.5) the sequence of flow patterns that occurs as rotational
speed is increased is quxte different from that occurring when the cylinders are
concentric. When the rotational speed of the inner cylinder is increased above
the Taylor critical speed, wave instability occurs first at a position just down-
stream of the point of maximum annular clearance. Because of the fact that flow
conditions and mean flow velocity change rapidly in the circumferential direc-
tion, wave instability for non-concentric cylinders never demonstrates the clean-
ly defined sequence of modes that occurs with concentric cylinders. Instead, the
wavy instability in the wide portion of the annulus becomes steadily stronger and
more irregular as speed is increased until the flow breaks down into a pattern of
large irregular eddies which completely obscure the vortex pattern. This break-
down occurs when the Reynolds number based on the mean clearance is of the order
of 500.
In the bottom three pictures in Fig. 6 is shown the flow at the same loca--
tions as in the top pictures but at a much higher Reynolds number (N = 6500).
At this condition the flow in the narrow portion of the channel appears to be
turbulent Judging from its 'grainy'" appearance. In the picture to the bottom
left in Fig. 6, showing the wide portion of the channel, one can note the re-
appearance of an axially periodic, vortex structure in the turbulent flow.
The various flow patterns described above for both concentric and non-con-
centric cylinders are typical of those which are shown in the motion picture
taken of the flow. The motion picture is a 16mm one, approximately 500 feet in
length, with no sound track. The film was prepared by MTI for National Aero-
nautics and Space Administration under contract NASw- 1021. A copy of the film
may be obtained from
A technical report on the experimental project in which the motion film was
made is listed as reference 2.
REFERENCES
To motor
Support
iw~m
Drive Quill
k<rtm m<^>
D
Suspending
Test Section
Eccentric Bushing
for Varying Eccentricity
of Cylinders
Outer Cylinder
2
Precision Bore
4.001 + 0.002 in I.D.
Gas Pime.
24 in. in Length Outer Cylinder
Counter Balance
Weight for Torque
Measurement
Altminum Inner Cylinder
3.641 + 0.001 in. C D .
Support -zJ~
Jluflgf
i»
MTI-2104
4.3.8
•I
MTl-2105
I 4.3.9
240 ~ aO
HTE-ZiaS
M
O
y^
a §
as
ON 03
4S>
o
o
00
O
t. -
4,3.11
•^<»^-"r ^"-^/v^
by
ABSTRACT
DISCUSSION
The location of the film-cavity interface and the pressure jump occurring across
the meniscus of this interface are of interest to lubrication engineers, since
such information bears on the proper boundary conditions for Reynolds' pressure
equation. Various boundary conditions have been proposed in the literature. The
Swift-Stieber condition which is most widely used ends the film where dp/d0 equals
zero. Birkhoff and Hays in a recent paper propose a condition of incipient counter-
flow for film separation in a diverging clearance. Other older conditions which
eliminate the negative pressures or terminate the film at minimum clearance are
still in common use.
This earlier work points out the need for determining correct boundary conditions
at the edge of a lubricating film, and has provided the motivation for the analysis
of free surface separation given in this paper. Here we focus attention on the
liquid-vapor interface of a film which is attached to a stationary surface and
* The work described in this presentation and its references was conducted by the
author at Columbia University, as part of his doctoral candidacy requirements. The
work was carried out under the guicJance of Professor H.G.Elrod,Jr, Columbia Univ.
and was sponsored by Mechanical Technology Inc.
is being swept along by an opposing moving surface. The physical situation thus ^P
described is, depicted in Fig. 1. It is one of fundamental interest, and appears
in many applications besides cavitation in bearings.
In the theory we consider a smoothly separating film, i.e. fluid velocities only
in the x-y plane. The object of the analysis was to find the shape of the liquid-
gas interface from the point of attachment out over the lower moving surface.
The fluid was assumed to be Newtonian in character with constant density, viscosityj
and surface tension. The following conditions are satisfied exactly in the analysis.
On the interface:
Some typical solutions of film height and surface velocity are shown in Fig, 2
for different values of the parameter N. N is defined as:
N - (3UU/T)^^-^
For any given solution, one can recover the constants in the quadratic velocity
i ^
4.4,3
representation, and thereby construct the streamlines of the flow. This has been
done for the case of (N) ~ 0.5 in Fig. 3. The separating streamline (x = 1,0),
which divides the recirculating flow from the flow being swept away on the moving
plate, is shown. A further consequence of the assumptions is that the fluid is
stagnant along a line normal to the boundary at 0 = jt/2 . The solution therefore,
corresponds to the situation depicted in Fig. 3. When 0 is close to Jt/2, the y
coordinate samples the flow far upstream. For fully developed upstream flow, the
assumed quadratic velocity distribution was proved to be perfectly consistent with
Reynolds' equation. This consistence was one of the objects of the analysis and
was exploited in developing film edge boundary conditions for Reynolds' equation.
Boundary conditions for ase with Reynolds' equation can now be developed. First
the pressure jump across the meniscus due to surface tension establishes the
pressure in the film just upstream of separation. Second, the flow carried under
the meniscus must equal the upstream Reynolds' flow. This establishes the pressure
gradient at film separation. The pressure and pressure gradient boundary conditions
for use with Reynolds' equation are shown in Fig. 5, in non-dimensional form, as
a function of parameter N. The pressure gradient boundary conditions according
to Swift-Stieber and Birkhoff-Hays, which are independent of the parameter N, are
also shown. As an example of the application of the boundary condition at film
separation of Fig. 5 the pressure profile for an infinitely-wide, parabolically
varying clearance was calculated and then plotted in Fig. 6. These pressure pro-
files are for the case ofa fully flooded inlet, that is, the film is assumed to
begin infinitely far upstream. The abscissa 0 of Fig. 6 is the x-coordinate,
distorted such that 0 = 90 corresponds to x - oo and 0 = 0 corresponds to the
point of minimum film thickness. Note the relatively large sub-ambient pressure
region associated with the pressure profiles for small (N) . Wl^
The object of the experimental portion of the investigation was to test the
validity of the boundary conditions derived from the theory. The experiment
was of simple design requiring only measurements which could be made easily
and accurately. The basic set-up consisted of a cylindrical slider bearing
riding in a fluid film. Measurements were made of the load on the slider and of
the location of the inlet and exit film boundaries. For comparison, Reynolds'
equation was solved for the geometry of the slider using the boundary conditions
of the theory as well as those of Swift-Stieber and Birkhoff-Hays. The re-
lationship between the load and the location of film separation was computed for
nearly flooded inlet conditions. The experiments were run for very light loads
and anall values of (N). Under these conditions, surface tension effects become
important and the present theory shows the greates t disagreement with others in the
literature. However, it should be emphasized that under many conditions of load
and N, the various results do not differ substantially.
Figure 7 shows a schematic sketch of the experimental apparatus. The slider was
a one inch length cut for a cylindrical glass lens of two inch radius. This
section was then cut in half. The moving surface was a glass disc placed over a
phonograph turntable. The lubricant was Dow Corning silicone oil. The slider was
attached to the end of a phonograph tone arm by means of a dial indicator attachment.
The tone arm had a point pivot which allowed free motion of the arm in all planes.
The turntable speed was first set to give a desired value of N. The ball and socket
joint of the dial indicator attachment was used to set the slider at any arbitrary
value of the angle J. The angle 7 was measured optically making use of a mirror
attached to the upper surface of the slider. It should be noted that the angle 7
which is measured optically bears a geometric relationship to the angular distance
between the point of minimum film thickness and the edge of the glass slider.
Therefore, when film separation coincides with the edge of the glass slider, the
measurement of 7 is tantamount to measuring the location of film separation. The
sliding weight on the tone arm was used to load the slider to the point where the
film was just on the verge of receding from the slider edge. Fine adjustment of
load was made by lowering a very light chain onto the slider surface. The location
of film inlet was measured by means of small scale attached to the side of the fQI
slider. This film inlet location measurement does not require too much accuracy
4.4.5
since small deviations from the fully flooded inlet condition result In small
differences in the load carrying capacity of the slider.
As load was applied to the slider, a point was reached where the film began to
recede upstream from the slider edge. This recession was always accompanied by
some flow in the lateral directions giving the familiar periodic striated flow
pattern observed by other investigators. The edge of the lens-slider therefore
served two purposes: one to stabilize the separating film and therefore force
compliance with the planar-flow of the thereotical model; and two to provide
a means of locating the film edge. Figure 8 shows the recession of the film
from the slider edge and the start of the striations for the case of N = ,5.
With 1.6 grams on the slider the film is separating from the slider edge. With
an additional .2 grams, the film begins to recede upstream and the striated flow
pattern is initiated. With 3.6 grams on the slider, the striated flow pattern
is fully developed.
The experimental data for N = 0.5 are plotted in Fig. 9 along with the predictions
of the present theory, the Swift-Stieber and the Birkhoff-Hays theories. Along
side of each data point is recorded the observed location of film inlet 0,. (9 equals
o
90 corresponds to a fully flooded inlet). In all theoretical predictions, the
pressure at inlet is assumed to be zero, and no correction has been made for the
weight of the fluid. The data are in best agreement with the present theory. However,
agreement becomes gradually poorer when the load on the slider is reduced. For
such lightly loaded conditions, the theory underestimates the load carrying
capacity of the slider. In Fig. 10 the same experimental data are compared with
the theory where the inlet pressure has been corrected for the hydrostatic head
Pgh.I a small correction has also been applied to account for the weight of the
fluid. The agreement between theory and experiment is excellent. Experimental
results for the case (N) =1,0 also gave excellent agreement with the theory.
For more detailed description of the experiment, discussion of experimental errors
and the calculation of load carrying capacity based on the boundary conditions of
the present theory, Ref. 2,
2. For (N) smaller than approximately 1,34, a stagnation point appears on the
liquid-vapor interface at a height of three asymptotic film thicknesses,
marking the onset of flow recirculation just upstream of the separation point.
References:
2. J.C„ Coyne and H.G. Elrod, Jr., "Slider-Bearing Test of a Theory For Film
Rupture or Separation" MTI Report 66TR345 June, 1966.
^^>/^^^ R^TE
\
27 - XBRETHERTON'S
\ DATA
A PRESENT THEORY
15- 0 MODIFIED TAPE THEORY
9 - V
6 - TAYLOR'S\
DATA \ .
^^^"-\ COX'S
3 - ^^5^___ DATA
« ^ A
2 - o
1 1 I I I 1
0.2 •0.5 -1.0 -1.5 -2.0 -5.0
Fig. 4 Comparison of Bubble-Experiment Data with the Present Theory and also
with the Tape Theory which is Modified by Retaining the Demoninator
in the Curvature Term
.001 190
R='ii
APPROACHES-1
AS N-*0
R=10
5
APPROACHES -^TIS
AS N — 0 "^^
N =-0.2
-60" -30 30 60
Fig. 6 Pressure Profiles for an Infinitely-Wide, Parabolically-Varying
Clearance (h = h + T^IIT) where 2r/h = 900. The Pressure Profiles
were Calculated using Boundary-Conditions in Accordance vith the
Present Theory.
r
BALANCE ABOUT
X AXIS
4:^
Fig. 7 Sketch of the Apparatus
U3
MTI-2166
4.4.14
•I
Fig. 10 Gamma (7) versus Load Data for N = -0.5, i-iU = -0.76 dynes/cm, compared
with the corrected theory (solid Imes: ¥, = -89, -88, -87, -85, and
-80 Deg,) The first digit nSii of 9. alongside the data has been omitted.
t.
5.1
SECTION 5
TECHNOLOGY II - SEALS
5.1.1
INTRODUCTION
The viscoseal is a special case of a general class of parallel groove
geometries which have many potential applications for use as p r e s s u r e
generation devices. For example, the parallel groove has been applied
in spiral groove thrust bearings (ref. 1), spherical and conical spiral
groove bearings (ref, 1), helical groove pumps (ref. 2), as well as heli-
cal groove seals (viscoseaL ref. 3).
The viscoseal is potentially useful for sealing liquid metals, such as
mercury, potassium, and sodium, in space electric-power generation
systems that require 1 to 3 years of unattended operation and near zero
leakage (refs, 4 and 5). As pointed out in references 5 and 6, the visco-
seal has inherent reliability and long life because of the absence of solid
surfaces in rubbing contact.
To date, analyses of parallel groove geometry have been limited to
creeping flow solutions where convective inertia effects have been ne-
glected, Convective inertia is the fluid mass acceleration due to a spa-
tial velocity change, for example, a change in velocity due to an obstacle
in the flow path. This is contrasted to an inertia effect caused by a tem-
poral velocity change, for example, a variation in velocity at a point in a
bearing due to shaft runout (which is a function of time). The creeping
flow analysis is restricted to flows where the modified Reynolds number
(see DISCUSSION) is much less than 1; that is,
Re = R e ^ l^\ « 1
or
Re"< 0.01
where
Re, ^ U L ^ p U L
£ < i
L
Creeping flow analyses for pumps, screw extruders, seals, and
thrust and journal bearings have been published by many authors ( e . g . ,
refs. 1 to 26). Appendix C contains an outline and discussion of the pub-
lished creeping flow solutions. The main criticism of the creeping flow
solution (generally used in the laminar flow regime) is that its validity is
restricted to Reynolds numbers below most engineering applications.
A number of semiempirical and turbulent flow formulations have been
devised (refs. 4, 18, 27, and 28) for operation in a regime designated as
turbulent, but they neglect the convective inertia t e r m s which are impor-
tant physically. The main criticism of these semiempirical analyses and
turbulent flow model is that the existence of the turbulent regime is a s -
sumed. Turbulence, however, may not be necessary to explain reported
experimental results and may not even exist. Even if turbulence does
exist, convective inertia effects a r e thought to be far more significant.
Kettleborough (ref. 29) numerically analyzed the slider bearing with in-
ertia, turbulent, and viscous t e r m s considered. When inertia only was
considered, the results were in qualitative agreement with published
(turbulent-attributed) slider bearing experimental results. Kettleborough
concluded that the turbulence term did not appear to greatly affect the
operation of the slider bearing,
A bearing number
AG empirical sealing parameter (defined in ref. 5)
II absolute viscosity of fluid
p kinematic viscosity of fluid
p* kinematic viscosity of fluid at reference state
I coordinate in direction of plate relative velocity vector
p density
p* density at reference state
ff shape parameter
# dissipation function
(p angle traversed in polar plane
V vector differential operator Del or Nabla
2
V Laplacian operator
Subscripts:
c based on clearance or film thickness
g groove
L based on characteristic length
M moving surface
n normal direction
r ridge or land
S stationary surface
t tangential direction
z across ridge-groove
e plate relative velocity direction
f] normal to plate velocity direction (axial length)
APPARATUS AND PROCEDURE
Experimental Study of P r e s s u r e Patterns
Figure 1 shows a schematic drawing of the experimental apparatus
and viscoseal assembly used to study end effects and pressure patterns.
The rotor of the viscoseal is attached to the power input shaft which is
supported by externally pressurized gas journal bearings. Axial thrust
due to sealed fluid cavity p r e s s u r e against the rotor is resisted by the
externally pressurized gas thrust bearing. The power input shaft is
driven by a variable-speed electric drive and step-up transmission. A
magnetic pickup monitors the shaft speed. The viscoseal housing is at-
tached to a support shaft and roller bearing assembly, which permits
axial adjustment. A pump pressurizes the seal housing cavity and circu-
lates a cooling fluid flow from the reservoir, through a heat exchanger, to
the test cavity and then back to the reservoir. This cooling flow is neces-
sary to control fluid film temperature in the viscoseal. The steel rotors
have a centerline average surface finish of 20x10" inch and a r e mounted
within 0. 0002 inch of the total indicator reading. The viscoseal housing
is constructed from acrylic plastic and contains thermocouples and p r e s -
sure taps as shown in figure 2. P r e s s u r e taps, in a line parallel to the
axis, a r e located at groove edges and midland and midgroove positions;
pressure taps a r e also located near both edges along one groove. Two
sets of four pressure taps (90^ apart) were used to aline the housing with
respect to the rotor within 0. 0004 inch as determined by calibration.
Thermocouples were located flush with the bore and groove root in a line
parallel to the axis and at each midgroove and midland position. Figure 2
also gives the groove dimensions used in this evaluation. The housing had
1° helix angle. Groove widths of 0. 19 inch and
a 2. 001-inch bore and a 14-^
land widths of 0. 13 inch were provided by using five helix grooves (five
starts). These geometric proportions, which are based on the optimum
relations given by reference 19, included a 0. 015 inch groove depth and a
rotor radial clearance of 0.005 inch.
The kinematic viscosity-temperature relations for the mineral oil used
were found to be the same before and after use in the test. These values
can be found in table I.
Experimental Study of Gas Ingestion and Sealing Capacity
In studies on sealing water (fig. 3), the water p r e s s u r e and flow to the
pressurized cavity (at the viscoseal high-pressure end) were controlled by
throttling valves on the inlet and outlet water lines. This arrangement
provided temperature control of the sealed fluid. A transparent housing
permitted visual observations on the rate of gas flow out of the pressurized
cavity; and, in studies on grooved rotors, a stroboscope permitted obser-
vations of the gas ingestion process in the rotating grooves.
Figure 4 is a diagrammatic sketch of the experimental apparatus used
for evaluating viscoseals in sodium. The temperature was maintained by
a furnace (resistance heaters) surrounding the viscoseal housing. The
temperature was monitored by thermocouples attached to the seal housing
outer diameter and by one thermocouple submerged in the sodium in the
pressurized cavity.
The enclosure surrounding the viscoseal assembly is pressurized
with argon (after prior evacuation) to slightly above ambient to insure
exclusion of air. P r i o r to the introduction of liquid sodium, the input
shaft is set at some fixed speed. Sodium is introduced by pressurizing
the reservoir and venting the leak detector tank (see sodium supply system
in fig. 5). The sodium first fills the pressurized cavity (viscoseal high-
pressure end) and then is allowed to fill the leak detector tank to a p r e -
determined level. The leak detector tank float displacement is monitored
by recording the differential transformer output that is produced by the
transformer core attached to the top of the float. Thermocouples placed
at different levels in the tank provide calibration during the filling process.
h
rotor arrangement. In a plane orthogonal to the rotor centerline, the
p r e s s u r e increases a c r o s s the groove (from point A to B) in the direction
of rotor rotation and decreases across the land (from point B to C). The
increase and decrease in p r e s s u r e repeat for each groove-land pair and
produce a saw-tooth p r e s s u r e profile in the orthogonal plane. The p r e s -
sure increases from the low-pressure end (plane 3) to the high-pressure
end (plane 1), and the p r e s s u r e pattern has a helical twist corresponding
to the helical grooves. (For a rotating grooved shaft this pressure pattern
would be rotating.)
The pressure gradients in the planes orthogonal to the rotor center-
line a r e essentially linear for axial positions not influenced by the end
effects. Figure 7 shows typical experimental results. These pressure
profiles a r e the same a s those for the orthogonal planes shown in figure 6
except that a single groove-land pair is unwrapped to form a plane figure.
(Neglecting curvature is not a significant e r r o r since the ratio of rotor
radius to groove depth is 65 to 1.)
In these orthogonal planes, the ratio of pressure to cavity pressure
P/P increases linearly over the groove and then decreases linearly over
the land. This linear relation was found to hold to the maximum test speed
of 6000 rpm, which produced a modified Reynolds number (Re*) of 0. 25.
End Effect
"" When the helical grooves connect directly with the pressurized
cavity, the p r e s s u r e developed along the land leading edge decays sharply
near the high-pressure end of the seal. Typical results a r e shown in fig-
^^ ure 8 for 1000 and 6000 rpm. The pressure pattern at the end was essen-
tially the same for the speed range investigated (1000 to 6000 rpm) except
that the results for 5000 rpm (not shown) and 6000 rpm start to show a
slight shortening of the end-effect length. This shortening is probably due
to convective inertia effects. In all cases, the end-effect axial length was
approximately equal to one-half the axial groove width. This suggests that
end-effect length could be taken to one-half groove width without significant
error. The pressure decay within the groove is due to the equalization of
the land leading- and trailing-edge pressure, since a common p r e s s u r e
must exist in the cavity.
Gas Ingestion and Sealing Capacity Observations
When Sealing Water
Grooved housing. - The series of photographs in figure 9 shows the
gas ingestion process when sealing water with a grooved housing and
smooth outer diameter rotor. At a modified Reynolds number (Re*) of
1. 75 (2000 rpm), no gas ingestion is evident, the water film is clear, the
interface is stable, and the scavenging length is nonwetted.
At Re* = 2. 58 (3000 rpm) the film had large gas pockets, gas was
ingested as evidenced by air bubbles passing out of the transparent vent
line, and the sealing capacity (parameter, AG) was lower than that at
Re* = 1. 75 (2000 rpm). Various degrees of gas ingestion and film rupture
a r e evident for Re* -^ 3. 50 (4000 rpm) to 6. 10 (7000 rpm). At Re* = 7. 00
(8000 rpm), no detectable gas ingestion rate was evident, and gas ingestion
was not detectable from Re* =• 7. 00 (8000 rpm) to the highest speed evalu-
ated, which gave Re* ^ 10. 84 (12 000 rpm). At Re* = 7.00 and above,
the seal wetted length had two distinct regions (fig. 10(f)): a gas-liquid
region, and a region of homogeneous liquid at the high-pressure end.
5.1.15
h
Between Re* = 3 . 5 and 6. 10, it was observed that the bubble size
progressively decreased when moving from the low-pressure end to the
high-pressure end. It is probably that with longer seal lengths the c o r r e -
sponding higher p r e s s u r e s would reduce the bubble size sufficiently to al-
low escape over the lands, and thus gas ingestion would be eliminated. It
was observed that concentricity of rotor and housing markedly affect gas
ingestion occurrence.
Grooved rotor. - In a series of visual observations using a grooved
rotor, the gas ingestion started at Re* = 1, 9 and increased in rate with
increasing Reynolds number (maximum investigated, Re* = 10.84). Below
Re* = 1 . 9 the liquid-to-gas inteiiace was stable and no gas ingestion was
observed.
Gas ingestion mechanism. - These observed gas ingestion character-
istics of the grooved rotor and the grooved housing are illustrated in fig-
ure 10. The probable mechanism producing the difference between the gas
ingestion process of the grooved rotor as compared with that of the grooved
housing is illustrated in figure 11. For the grooved housing, when the cen-
trifuge action on the liquid becomes great enough, the gas bubbles are d i s -
placed to the rotor surface and pass over the lands. Therefore, the bub-
bles a r e not pumped to the high-pressure end. For the grooved rotor, the
centrifugal action on the liquid forces the gas bubbles into the grooves;
therefore, the bubbles are pumped to the high-pressure end of the seal.
When a viscoseal is ingesting gas, the axial pressure gradient is non-
linear and increases as the liquid-to-gas volume ratio increases in the di-
dj^ rection toward the high-pressure end. For no gas ingestion, the axial
pressure gradients a r e linear as shown in figure 8, which contains data
5.1.16
increases to the rotor outer diameter, the gas flows back out the annulus
formed by the rotor and housing. This causes a decrease in bubble size
and a sharp drop in the leak detector readout trace. If the interface is far
enough away from the seal low-pressure end, the liquid will be scavenged
back within the nonwetted length, and no loss of liquid will occur. The in-
gestion process will begin again and the net result is a saw-toothed trace,
each sharp drop indicating gas blowback. If the interface is sufficiently
close to the seal end, liquid loss will occur, and the leak detector readout
will show a saw-toothed trace with a downward trend, which represents
a liquid loss at each blowout (fig. 14(c)). This latter action may be the
"seal breakdown" phenomenon reported in reference 21 and the "secondary
leakage" phenomenon reported in reference 4. (In sealing water, the r e -
peating bubble growth and subsequent seal blowout was visually observed
by means of a transparent housing.)
starts (see table II) with and without secondary grooves on the lands. In
both cases, a higher sealing coefficient was obtained through the use of
secondary grooves. Secondary grooves also operate at lower power ab-
sorption levels for equal sealing capacity, as was determined by the
amount of heating power required to maintain operating temperature.
This lower power absorption is due to less land area as compared with
the land area of the seals without secondary grooves.
Sodium was sealed-for p r e s s u r e s ranging from 2 to 100 psig with
negligible leakage loss. The operating range was within Reynolds numbers
(Re*) of 8. 66 to 50. 5, and the fluid temperature was in the range 300°
to 625° F. The maximum length of any one run was 8 hours and accumu-
lated testing time amounted to 50 hours for 10 runs. The oxygen content
of the sodium for the 10 runs varied between 50 to 95 parts per million,
as determined by the analytical method described in reference 31.
THEORETICAL ANALYSIS OF THE VISCOSEAL
Basic Model
y-direction:
R,J^^^,^\ (3)
Dt ^
Dt W2 3y2 3z2y
where # is the dissipation function
.2 ,^,.2 ,^...21 + /9Z + /aw ^ avA2
*-p^©^ey lax
9UA2
ay/ \ay az/
F x - ^y- Fz == 0
This modified Reynolds number gives the relative magnitude of the inertia
forces to the viscous forces; that is,
Re* = I^srtia forces
Viscous forces
Thus, inertia forces can be neglected only if Re* < 0. 01.
After the formal ordering procedure (see appendix A), the three-
dimensional flow field equations are as follows:
5.1.23
h
Conservation of momentum:
,V au + w
, au m+p^ ap = constant
where — (9)
ay az p ax ay 2 ax
•I
if s 0 (10)
ay
Conservation of m a s s :
i l +l ^ = 0 (12)
9y az
The boundary conditions for an arbitrary shaped groove-ridge geome-
try (see fig. 20) are:
(1)
v= 0
u = U cos a /at y = 0
w = -U sin a.
(2)
u =v =w =0 at y = h(z)
For example, when h(z) is a rectangular groove-ridge pair the boundary
condition is:
u =v =w=0 at c < y < c + h (on side of groove)
u = v = w = 0 at y = c (on ridge)
u = v = w = 0 at y = c + h (on groove root)
(3) Periodicity of groove-ridge pairs at the clearance interface.
For 0 < y < h(z)
v(0, y) =- v(L sin a,y)
w(0, y) ^ w(L sin a, y)
P(0,y) = P ( L s i n « , y ) - (AP^+g^z
where (AP J ^ = const
5.1.25
W
All velocity derivatives must also be the same to satisfy the periodicity
boundary conditions
ap N(AP^_)
r+g'z
ax 27rR cos a
These equations (9), (11), and (12) can be solved for the velocity com-
ponents. The velocity components can be obtained for arbitrary (AP „) .
Then the axial pressure gradient dF/dV can be found from
ap ^Paxial ^N)(AP^ )
= = _2— cot a
dh' V L
where L ' is the axial length (see fig. 21). Qualitative pressure distribu-
tion from experimental data in the x - and z-directions is shown in figures
22 and 23.
By examining the analysis, boundary conditions, and the derived equa-
tions, the axial pressure gradient dF/dV is a function of (a, U, p, i^),
where cr is a shape parameter that is a function of (b/h , c / h , a/b, a).
5.1.26
4
Since convective inertia effects are important, new geometric parame-
ters appear which are not implied in creeping flow theory. Most notable of
these is the aspect ratio b/h , which implies the degree of convective ef-
fect. When the b/h and a/b ratios are known, N, the number of helix
or groove starts, can be found. (Direct effect of N is seen from boundary
condition (4).) Thus the number of groove or helix starts also is an impor-
tant parameter. Since inertia is important, a high density is desired for
pressure generation.
The equations (9) to (12) and boundary conditions are the minimum
number of terms that are necessary to mathematically describe the physics
of the internal flow. This is the simplified form of the Navier Stokes equa-
tions which should be used for an analytical solution; however, a mathe-
matical function to describe the groove-ridge shape may be too complex
for a closed-form analytical solution. Thus the equations are solvable by
using high speed digital computer numerical methods. Since numerical
schemes can accommodate all two independent variable terms with equal
ease, it is suggested that the complete set of two-dimensional equations be
solved numerically, i. e.,
v i I + w l Z = - l l P + . f l ! z + l?Z^ (13)
ay az p ay ^^^2 3^2j
vi^ + w i l = - l i P + . / l ! ^ + i M (11)
ay az p az \^^2 ^^2)
11 + ^ ^ = 0 (12)
ay az
5.1
^!l,^!l.J.[^.fl],JL^,^ (15)
ay az Pr^gy2 g^V PSC pgC
Comparing the above equation with the x-direction momentum equation
it is seen that the mathematical form of both equations is identical.
Substituting 6 =pu
and
4
The equations a r e of a form that might tempt one to reduce the num-
ber of independent variables by one by trying to find the admissable flow
field which would result in a similarity transform solution. This tech-
nique fails because of the characteristic lengths that a r e present in this
type of confined flow.
The equations a r e solvable by using high-speed digital computer nu-
merical methods. A mathematical function to describe the groove-ridge
shape may be too complex for a closed-form analytical solution. It i
should be interesting to find out if the step is still the optimum shape
with convective inertia effects included. For creeping flow, Lord Rayleigh
found in 1912 that the step was the optimum shape by using calculus of vari-
ations.
The equations and ordering procedure should be reexamined carefully
before extending these results beyond the assumptions used in this analysis.
Discussion of Modified Reynolds Number
A check was made on the significance of the modified or reduced
Reynolds number by calculating Re* from the available published data at
the claimed onset of turbulence (breakpoint). In table III it can be noted
that the departure in all cases from a constant sealing coefficient occurred
when Re* was in the range of 0. 26 to 2. 2, which suggests that the so-
called beginning of turbulent operation is really the beginning of significant
convective inertia effects. This would be in agreement with Kettleborough
(ref. 29), who found qualitative agreement between his inertia solution and
experimental turbulent results.
For seals 5, 6, and 7 (see ref. 27 and table III), only the groove to |n
ridge ratio has changed, As the ratio increases, the breakpoint decreases.
5.1
(10)
ay
aw w 9w _ 1 ap ^ J^a^w ^ a V (11)
V +
ay az p az 'Vay2 az2
Conservation of mass:
3v ^ aw _ Q (12)
ay az
With the following boundary conditions
(1) At the moving smooth plate surface,
u = U cos a V=0 w = -V sin a
(2) At the fixed parallel groove geometry surface,
u -= V = w = 0
5.1.35
(4)
ap ^ N(APr^g)^
ax 27rR cos a
From the development and analysis of the equations, the following
results and conclusions can be made:
1. A modified or reduced Reynolds number evolved which gives the
relative magnitude of the convective inertia forces to the viscous forces:
pfu^+vll.w^ZV-^+Mf^ +^ +^ l (A2)
ax ay az/ ay \^j ^^2 9^2^
L cos ff c L sin a
U COS a V U sin a
The characteristic "dimensions c and L can be seen in figure 17.
It is noted that a trigonometric relation exists between x and z charac-
teristic lengths and that c/L < 1. The characteristic length was selected
to be the length of the groove-ridge pair in the direction of the smooth
plate velocity vector. This is the principal length over which the driving
force acts. The flow behavior over each groove-ridge pair is the same.
Let p* and v* be the fluid density and viscosity corresponding to a
reference state (P*, T*). Therefore,
17=ii=l
and
p = £^=.i
P*
for incompressible fluids with constant therm ©physical properties. Now
to find the relative magnitude of V so that all t e r m s of the incompressible
continuity equation will be of the same order, the incompressible continuity
equation is nondimensionalized:
5.1.39
b
To have all t e r m s of the equation the same order of magnitude, —— must
_ , U c
be of the order of 1. Therefore, V = U — , and thus v = — —.
L U c
Next the equations of motion are nondimensionalized. The pressure
term must first be nondimensionalized; even though the pressure is not
the driving force, it is the function of primary importance. Thus the char-
acteristic p r e s s u r e is not extremely important. Since the pressure gradient
is related to the rotor speed, it appears to be logical to nondimensionalize
the p r e s s u r e with respect to the dynamic pressure pU (also called the
dynamic head):
PU2
It can be seen that the characteristic dynamic pressure pU has r e -
duced the number of nondimensional groups by one (i. e., if P * (a refer-
ence pressure) had been used, another nondimensional group would be
present in the equations). Also, since the incompressible case is being
analyzed, the p r e s s u r e appears in the dynamic role only, not in the thermo-
dynamic role.
Nondimensionalization of x-Direction Momentum Equation
When the x-direction momentum equation is nondimensionalized,
equation (1) results in
u _31 au + w----==
- + v-— du 1 ap
—
ax ay az „^„2 ax
•^ cos a
I7* 1 d\ ^ (V^ d\ ^ 1 d\
+ —
UL
Next a Reynolds number is defined based on the characteristic length L:
Re, =^
where
p = p*'v = V*
and thus
17 au . „ au , =. 31 1 ap"
U +v .+ W — = - —
3x a? dz „^„2 ax
' cos a
+ ^ I ^ a ^ ^ fhy 2 -d^
^ ^ 1 a\ (A6)
Re^ ^2„. .~2
— \c J ^--2
— 2 —2
cos^ff ax^ ^"^ ' ay"" sin a 3z
The convective acceleration and pressure t e r m s are of unit order; thus,
the viscous term must be of unit order also. The largest term of the
22 -.2^
2— ,.. , ,xv2
viscous forces is (^] •^-^, since (—] < 1; therefore, the _ i — ( _
\c / p-»2 VLJ Re, \ c ,
ar ^"' ""L '^ ^
should be of unit order. The modified Reynolds number is defined by
Re* = R e L l ^ I • This is the identical "reduced" or "effective" Reynolds
number that is used in the bearing theory to show the true relative magnl
tude of the inertia and viscous forces; therefore,
Rp* = Ifl^rtia forces
Viscous forces
(e. g., see ref. 32). Now for curved geometries, a Reynolds number
based on clearance has appeared. This is especially true in viscoseal
analysis where
5,1.41
nL
which leads to
2
R e * = RCj^
(i) - M I )
It does not matter which Reynolds number is used to find the relative
inertia force to viscous force ratio provided that c/L < 1, which c o r r e -
sponds to the conditions
e2 »
a2
3y2 3z2
and
32 32
»
3y^ 3x2
1. i i ~ H + 13 0 0 0 l ^ + m ^
^«L V 3x2 9^2 9^2
h
1 d\ ^ fh\ d\ ^_ /c 1 d\
Re. (i)z: 2 —2
COS « 9x
c / 3^2
ay
L /' sm
. 2a az
5-2 __
l i both equations it is seen that d\/dy is the largest term. Now the
modified or reduced Reynolds number is
2
Re*-=ReL^|] =Re^^|^=L75
l ^ + f 9 v ^ ^ 3 v _ / L X2 - .3^ P ^ ^ 1 3^
ax 3f az \c / 3y Re cos^a dx^ ^^ ' 3y^ sin^a 3z^J
(A8)
Again using the definition of Re* and examining the terms of unit order
u^, v ^ , w l Z , 1 8 ^ mi)
ax' ay' az' Re* g™2
thus
2 n^
0(1)
\c / ay
Since
ay
therefore
^ s 0 (A9)
ay
l ^ +^ i l +w ^ + 4 f ^ f / l V ^ B i n « w + cos«l)2=^/Lf al
ax ay 3z R \ C / \ R / \ c / ay
b Now
2
±. fL\ f±\ (-sin aw + cos au)^ -' O'H)
Comparing the largest viscous term, pressure and centrifugal inertia term
Now use a typical Re, (as in example) = 103 and apply restriction
- < 0.01
R-
Again
9P^0
ay
if
^ < 0.01
R~~
which substantiates the claim that the centrifugal effects should be small.
In this model Coriolis forces were also thought to be negligible. The
above centrifugal effect is for the case when the grooves a r e on the sta-
tionary housing. A slight modification is required for the groove-ridge
pattern rotating. Note that as c/R — 1 the centrifugal force effects are
no longer negligible.
Nondimensionalization of z-Direction Momentum Equation
In a like manner the z-direction momentum equation is nondimension-
alized and appears in a similar form as the nondimensional x-direction
momentum equation (A6)
U ! H + v ^ ^ .+
-^w^ ^^-^
' ^ -= 1 3P_^^ V* 3^,, (hf ^+_L 3^
Sx 3y" 3z sm^a 9^ UL 2 —2 3y^ siE^ff 3z 2
cos ff ax
(A 10)
Examining the viscous t e r m s
.2 . ^
1 1 )V
Re, ^^ ^ 3y2 sin2« a^2_
cos 2 « 3 x 2
The largest term is (L/c)^ ( 3 ^ / 3 y ^ ) and it would appear that the 3 ^ / 3 1 ^
term could be neglected; however, this term must be retained due to a
physical argument. Because of the "edge" there is a sharp gradient in the
z-direction and the nature of the flow is such that the 3 w/az term must
be retained. Thus the z-direction momentum equation in dimensional form is
1 ap _^ „/3^w ^ d \
u i Z + v^E + w ^ (All)
3x ay az p dz
.3y 3z
This was pointed out by Dr. A. Mager of Aerospace Corporation, Los
Angeles, California (personal communication, October 5, 1966).
Nondimensionalization of the Energy Equation
The energy equation for steady temperature variation for an incom -
pressible fluid with constant thermophysical properties and negligible
compression work is
3T 3T „. 3T k /3^T 3^T 3^T\ M * Q*
(A12)
U —=- + V -^-^ + w
3x 3y 3z pgC \ 9^2 pgC pgC
3y az
Nondimensionalize as before and nondimensionalize the temperature by
letting
e= and # ' = —
'M
5.1.47
k Thus
_ 30 - 30 . ^ 30 3^0
u +V +w
ax ay 3z IpgC Uc' _ „ 22 a^ ax
VL] cos ,^2
—2 ,^2
-^
3y U i .^J
sin^a 3z^
#^
Ql
LpgC(TM V PgC(T^-"Tg)
now
/ k J i L\ J ^ _ J _ / L\ ^ 1 ^ J___
I pgC y Uc cy P r Re^ I c / PrRe* P*
or
^ aa
ax
- 30
u— + Y — + w — = —^
3y
•'
— 30
az
1
P*
e
(tP 3^0 ^ 3^0 + i £
cos 2 a 3x
„~2 3^2
3y ^^
1
sin^a 3z^
3^0
#' Q*
LpgC(T^-Tg) pgC(T^^Tg)
..
V— 30 +
^w,„ 30
— k ^ll + ^]+JL.^ + ^
3y 3z pgC \^gy2 3z 2 / pgC pgC
5.1.48
4
P r = CM
k
therefore
tana=I^ (Bl)
L
Now going the same distance along the x-axis
cot a = 2 ^ ^ ^ (B2)
L'
Equating (Bl) and (B2)
tan a cot a
or
(p = 2v cos^a (B3)
The arc distance (path length) in the x-direction is:
S = YCR^)^ + (L')2
Upon substitution for L' and (p
S = 27rR cos a (B4)
The pressure differential at a from traversing the z-direction from
0 is
<^Po-a ^t-^Pr+g'zW
then
3P __ N(^Pr+g)z , N ( A P „ g ) , ^^^^
3x S 27rR cos a
It is seen that aP/ax = constant which follows from fully developed
flow along the groove-ridge direction and verified by experiment in refer-
ence 5.
APPENDIX C
SUMMARY OF CREEPING FLOW ANALYSIS
The optimization from creeping flow analysis has been formulated by
considering a flat-plate model (unwrapped cylinders, see fig. 27) and the
following:
(1) General incompressible plane Couette flow is assumed in the
^-direction (Poiseuille flow and simple Couette flow). The |-direction
momentum equation is
9 \ 1 ap
9^2 II d^
^ h 2/i d | V h/
(2) Poiseuille flow is assumed in the TJ-direction. The 17-direction
momentum equation is
^ \ „ 1 3P
g^2 II df]
w =A(l£V(?~h)
^ 2iU V3i7 /
Now the condition Q = 0 is imposed in the control volume; that is,
Q - Q^ + Q^ = f f \ d7] d? + ffw^ d? d | = 0
which means that there is no net end leakage or closed channel flow. This
is a restriction that limits the optimization and p r e s s u r e gradient predic-
tion to a seal only.
Solving for p r e s s u r e gradient that satisfies the Reynolds equation
(ref. 34) results in a sealing parameter
5.I.51
AP ^ 6jUUG
or
G == <=^ ' ^ p
G
GliVW / A
where A is the bearing number (ref. 34) and G is a geometry factor.
See references 1 to 26 for various expressions and values for G.
The creeping flow analysis has some paradoxes. It can best be illus-
trated by referring to figures 17 and 28 and by observing the resolution of
the drag force into components along the groove (x-direction) and normal
to the groove (z-direction). In creeping flow, the drag force is in equilib-
rium with the pressure force. Thus along the groove (x-direction) the
drag force is indeed in equilibrium with the pressure force. In the normal
groove direction (z-direction), however, the pressure and drag forces are
in the same direction, which is also the direction of a very large pressure
gradient (see figs. 22 and 23). Therefore, an unbalance of forces exists
in this direction. Physically, the Q = 0 restriction is not possible. Also,
this point can be illustrated from ^~ and rj-direction momentum equations,
respectively: ^2
^ \ _ 1 ap
9^2 M ai
' \ „ 1 ap
9^2 M a?]
To resolve the paradox of having the drag and pressure forces un-
balanced in the z-direction, a convective inertia force is in equilibrium
with the pressure and drag forces (a nonlinear effect). The step is causing
a convective change that results in generation of an axial p r e s s u r e gradient.
The role of the step is now described.
If a concentric rotor is rotating about a stationary smooth sleeve, no
axial pressure gradient is generated. This is the classic axisymmetric
Couette flow (see ref. 36). Furthermore, if an axial pressure gradient is
imposed ( e . g . , by having the rotor translate in the axial direction) an
"uncoupled" helical flow will result. Subsequently, of course, there will
be net flow out.
In the helical groove seal (viscoseal), the step or edge h is the
pressure generation mechanism. Without the step, that is, the limiting
case when h — o (rotor and housing smooth), there is no axial p r e s s u r e
5.1
V • VP = /iV^(V . V) = 0
or
V^P - 0
which is Laplace's equation, since pressure is a scalar quantity. Conse-
quently, the multitude of mathematical analog methods of solving linear
equations and specifically Laplace's equation can be used. Thus, solu-
tions for end effect can readily be found (refs. 1 and 34).
Muijderman (ref. 1) used conformal mapping to analyze the spiral
grooved bearing. Again, this technique can be readily applied because
pressure is a potential function.
REFERENCES
1. Muijderman, E. A.: Spiral Groove Bearings. Philips Res. Rept.
Suppl. 1964, no. 2 . : See also Scientific Lubr., vol. 17, no. 1,
Jan. 1965, pp. 12-17.
5.1.54
4
2. Rowell, H. S.; and Finlayson, D . : Screw Viscosity Pumps. Engi-
neering, vol. 114, Nov. 1922, pp. 606-607.
3. Boon, E. F . ; and Tal, S. E . : Hydrodynamic Seal for Rotating Shafts.
DEG. Inf. Ser. 13, United Kingdom Atomic Energy Authority, 1961.
4. King, Alan E . : Screw Type Shaft Seals for Potassium Lubricated
Generators. IEEE Trans, on Aerospace, vol. AS-3, Supplement,
June 1965, pp. 471-479.
5. Ludwig, Lawrence P . ; Strom, Thomas N . ; and Allen, Gordon P . :
Gas Ingestion and Sealing Capacity of Helical Groove Fluid Film
Seal (Viscoseal) Using Sodium and Water as Sealed Fluids. NASA
TND-3348, 1966.
6. Ludwig, Lawrence P . ; Strom, Thomas N . ; and Allen, Gordon P . :
Experimental Study of End Effect and P r e s s u r e Patterns in Helical
Groove Fluid Film Seal (Viscoseal). NASA TN D-3096, 1965.
7. Vohr, J. H.; and Chow, C. Y.: Characteristics of Herringbone-
Grooved Gas-Lubricated Journal Bearings. J. Basic Eng., vol. 87,
no. 3, Sept. 1965, pp. 568-578.
8. Billet, A. B . : Hydraulic Sealing in Space Environments. Proceedings
of the Second International Conference on Fluid Sealing, B. S. Nau,
H. S. Stephens, and D. E. TurnbuU, e d s . , British Hydromechanics
Research Association, Harlow, Essex, England, 1964, pp. C2-17 -
C2-36.
9. Rowell, H. S.; and Finlayson, D. : Screw Viscosity Pumps. Engi-
neering, vol. 126, Aug. 1928, pp. 249-250.
5.1.55
h
10. Rowell, H. S.; and Finlayson, D . : Screw Viscosity Pumps. Engi-
neering, vol. 126, Sept. 1928, pp. 385-387.
11. Rogowsky, Z . : Mechanical Principles of the Screw Extrusion Ma-
chine. Engineering, vol. 162, no. 4213, Oct. 11, 1946, pp. 358-
360.
12. Strub, R. A . : Spindle Drag Pump. Machine Design, vol. 25, July
1953, pp. 149-151.
13. Pigott, W. T . : P r e s s u r e s Developed by Viscous Materials in the
Screw Extrusion Machine. ASME T r a n s . , vol. 73, Oct. 1951,
pp. 947-955.
14. Anon.: Theory of Extrusion, Ind. Eng. Chem., vol. 45, oo. 5,
May 1953, pp. 969-993.
15. Eccher, Silvio; and Valentinotti, Aldo: Experimental Determination
of Velocity Profiles in an Extruder Screw. Ind. Eng. Chem., vol.
50, no. 5, May 1958, pp. 829-836.
16. Griffith, R. M . : Fully Developed Flow in Screw Extruders. Ind.
Eng. Chem. Fundamentals, vol. 1, no. 3, Aug. 1962, pp. 180-187.
17. Squires, P . H.: Screw Extrusion - Flow Patterns and Recent Theo-
retical Developments. SPE T r a n s . , vol. 4, no. 1, Jan. 1964, pp.
7-16.
18. McGrew, J. M.; and McHugh, J. D.: Analysis and Test of the Screw
Seal in Laminar and Turbulent Operation. J. Basic Eng,, vol. 87,
no. 1, Mar. 1965, pp. 153-162.
19. Zotov, V. A . : Research on Helical Groove Seals. Russ. Eng. J . ,
vol, 10, Oct. 1959, pp. 3-7.
Asanuma, T . : Studies on the Sealing Action of Viscous Fluids.
Paper No. A3 presented at the F i r s t Mernational Conference on
Fluid Sealing. Cranfield, England, British Hydromechanics Re-
search Association, Harlow, Essex, England, Apr. 1961.
Lessley, R. L . ; and Hodgson, J. N . : Low-Leakage Dynamic Seal-
to~Space. Paper 65-GTP-14, ASME, Feb. 1965.
Whipple, R. T. P . : The Inclined Groove Bearing. Rep. AERE-T/
R-622 (rev.), Research Group, Atomic Energy Research Establish-
ment, United Kingdom Atomic Energy Authority, 1958.
Gruntz, Robert D . ; and Rackley, Ray A.: Snap 50/Spur Power Con-
version System-Objectives,, Current Status and Lunar Applications.
Paper No. 650321, SAE, May 1965.
Golubiev, A. L : Studies on Seals for Rotating Shafts of High-Pressure
Pumps. Wear, vol. 8, no. 4, July/Aug. 1965, pp. 270-288.
Holan, Karel: Sealing in Engineering. Proceedings of the Second
International Conference on Fluid Sealings, B. S. Nau, H. S.
Stephens, and D. E. Turnbull, e d s . , British Hydromechanics Re-
search Association, Harlow, Essex, England, 1964, pp. E5-73 -
E5-88.
Stair, William K.: Analysis of the Visco Seal. Rep. No. ME 65-587-2
University of Tennessee, Jan. 18, 1965.
Stair, William K.: Theoretical and Experimental Studies of Visco-
Type Shaft Seals. Rep. No. ME 66-587-5, University of Tennessee,
Apr. 28, 1966.
5.1.57
•I
TABLE I. - VISCOSITY
OF MINERAL OIL
Temperature Kinematic
.Op
' viscosity,
°C cs
W//////A /////A/////A
0 0
( 0 I 10 ir ) (O 1 0 (0 0
0 0 6 ^ en -, 0 Oo
V////////////////////A
0 5^
( jeo-"! (O ] 0
S'^C
f-^—-— ——•—i
V//////////////////////A
0 =(6 0 0 0 010
(0 10 I n ( 0 0^ in ) (0 Oi (0 00--
V////////////////////////i
O 01'-
(0 1 ) (0 0
se''0'"jar
(ircoves
wmm^^'^m^
0 4 0 0 0 0
( 0 ] iO i n ) I (O L30 (O a 0 4
Usec! a s b a b e v
TABLE m . - CRITICAL MODIFIED REYNOLDS NUMBER FOH ONSET OF TURBULENCE OF VAMOUS EXPERIMENTERS
Experimenter Reference Seal Diameter, Radial Ratio of Ridge Groove Helix C h a r a c t e r i s t i c C l e a r a n c e Modified
and fluid in. clearance, radial width, width, angle, length, Reynolds Reynolds
c, c l e a r a n c e to a, b, a, L, number, number.
in. r o t o r r a d i u s , in. in. deg in. R«c Re*
c/R
Stair 27 ai 1.2430 0.0042 0.0068 0.0934 0.1596 14.5 1.01 400 1.7
Water 2 1-. 2465 .00235 .0038 .1176 .0494 9.67 .98 (b) (b)
2B 1.2420 .0047 .0076 .1176 .0494 9.67 .98 (b) (b)
3 1.2461 .00265 .0043 .0828 .0842 9.67 .98 350 .95
3B 1.2408 .0053 .0085 .0828 .0842 9.67 .98 (b) (b)
4 1.2461 .00265 .0043 .0514 .1156 9.67 .98 200 .95
4B 1.2408 .0053 .0085 .0514 .1156 9.67 .98 400 2.2
5 1.2455 . 00295 .0047 .1585 .0607 5.81 1.83 600 .97
6 1.2455 .00295 .0047 .1070 .1083 5.81 1.83 200 .32
7 1.2455 .00295 .0047 .0609 .1600 5.81 1.83 160 .26
King 4 and 27 3C 2 0.0029 0.0029 0.1063 0.1063 3.86 3.2 550 0.50
Oil, IB 2 .0020 .0020 .1063 . 1063 3.86 3.2 400 .26
water, po- 2E 2 .0026 .0020 .1000 . 1000 7.26 1.54 650 1.1
tassium 5D 2 .0015 .0015 .1000 .1000 7.26 1.54 650 .63
McGrew and 18 S F 9 6 - 5 oil 1 0.00325 0.0065 (c) (c) 14.5 3.14 400 0.41
McHugh Oil
Lessley 28 Configura- 2 0.0032 0.0032 0.075 0.122 14.5 0.79 350 1.4
Oil tion B
Water Configura- 2 .0032 .0032 .083 .136 4 3.13 500 . 51
Mercurs^ tion C
D
O
I
xl Liquid under
-Liquid-vapor ,-^.,Externally pressurized gas journal bearings
pressure -o
interface
Roller
bearings T-
Support shaft-
,, I \ •- Viscoseal rotor
^Viscoseal housing
^r-Groove depth, hg
/ 0.015 in.
05
CO
to
I
H
Housing twre,
2.001 in.
Section A-A
(a) Axial locations.
r A ^ ^ Pressure taps
m Circumferential locatfons.
Figure 2. - Pressure tap and thermcKOuple focations.
5.1.63
k Water out
-Groovewidth, b'
Groove depth, h^
1
T
^ - Input
Support shaft Rotation shaft
CD
to
I
I Radial clearance, c
Argon gas
atmosphere
t
^ ^ probe / :
-Heated shroud
w -Resistance heater wires
^Resist -Helix angle, a
I Line to reservoir
Differentia
to transformer-|
to
I
Leak rate
monitoring
tank Air-operated
Argon gas valve
pressure
Hand valve
Float
Sodium level
'777777777TT77777777I272Z777777777/jy^/7y.
•• Plane orthogonal to
sealcenterline
Pressurized
cavity
Characterfstfc
length - ABC
Housing, m •* Rotor
nal grooves- rotation
Rotor
centerline
Plane
Pressure decreas-
ing over land
Rotor
Pressure increas rotation
ing over groove
Sealcenterline
MII-2066
Pressure tap
n Pressure decreasing
r i ci5ur c iMci ycjbmy uver groovy nvpr land
r 1^
-• Axial distance
f r n m r7i\i\^i
^-^
. ^ ^ ^ - - \ \ \
.-^
X
^^^
.--- - ^ \ ^ 3- '
^ \ ^
(I) Rotational speed, 2000 rpm; average temperature of oil film, 46.2° C 1115° F); temperature
change along film, 0.8°C (1.5°F); cavity pressure, 30.4 Newtons per square centimeter
(44.1 Wsq in. absolute).
(b) Rotational speed, 4000 rpm; average temperature of oil film, 138° F; temperature change along
film, 4° F; cavity pressure, 60.2 Ib/sq in. absolute.
1.1.
End 1 1
-effect * i 0.030 in. from land
len( th leading edge
.„„
0.030 in. from land
1.0
\ \ tr ailing edge
N
\
.9
\ \
•^ \
\\ \
.8 \
\
s
s
7 \
.1 .2 .3 .4 .5
Axial distance
High-
pressure I
end | r ^ Liquid-gas
interface
t Rotation
^Rotor
High-
pressure
end
ging or •
nonwetted
-« Ingestion direction ''"^tti
MTI-2071
I 5,1.69
I*
High-
pressure
end
-Ingestion
direction
(c) Modified Reynolds number, 3 50(4000 rpm)
High-
pressure
end
Figure 9 Continued
> •
MTI-2073
5.1.70 I
•I
High-
pressure
end
wetted
-Ingestion direction length
(e) Modified Reynolds number, 5 25 (6000 rpm)
High-
pressure
end
Figure 9 Concluded
••
MTI-2072
H
CO
to
I
2 3 4 5 6 7 8
Modified Reynolds number
(a) Grooved housing. Bubble motion pri- (W Grooved rotor. Bubble motion has axial
marily cirtumferential. component due to helical grooves.
Figure 11. - Gas ingestion mechanisms of grooved housing compared with grooved rotor viscoseal.
72
^^^^K^'
High-
pressure
eno
^
itiJ^l
1 ^
End of scavenging
i 3h
MII-2067
5.1.73
Figure 13. - Gas ingestion rate when sealing liquid sodium at 400° F.
i_ To leak detector
pSealed liquid
Direction of
gas blowout
"^Liquid-gas interface
, '-Gas bubbles ingested towards high-pressure end
t \ . "-Gas bubble increasing in size
^—To reservoir
(see fig. 5)
(a) Mechanism involved in gas ingestion into pressurized cavity.
(b) Schematic liquid level oscillograph trace for repeated cycles of gas ingestion, subsequent
gas blowout, and no liquid leakage at blowout.
-Gas blowout
Gas ingestion-, Zero gas ingestion;
I
Leakage at blowout
Time
(c) Schematic liquid level oscillograph trace for repeated cycles of gas ingestion, subsequent
gas blowout, and liquid leakage at blowout.
^
^
•r^ ^ y -
._^¥ ^V
ZZ_B: /
y /
p <- '
^Z23
^I
U
1 1 1 M M 1
—o— Smooth-bore housing,
E
to (geometry 1)
to D- Smooth rotor, five grooves
in housing (geometry 5)
.1 i 1 1
103 10^ 105
Reynolds number based on clearance
Figure 15. - Comparison of sealing parameters for grooved rotor and housing
operating in sodium. Sodium temperature, 165° to 335° C (329° to 635° F).
P0)
dl
r^sr
~'a"
«
r/" '
CH Z J
V
r-ija
.•^L / X,
.B o t
i
r ,{/
c/^
I 1
—c) S moot i-b<)re hous ng, five groov eson •otor (ge 3m ;try 1)
— c ] — - S moot i-b<)re hous ng, five groov esand secondary and
gro ove >o 1 rotor (geometry 2)
,1 1 1
(a) Five helix starts.
1— [ —
J / ^'
55
y
^f'X}
j ^
M> r 7Sr^
r .
/ o (Y X)
, 1 ^ ^
C
'^ cf'
& o
1 1 M ' 1 1 1 1 1
— O — Smooth -bor e hous ng, 10 grooves on rotor (geometry . *
—n—Smooth -bor e hous ng, 10 grooves and secondary land
groo ves on roto r (geometry 4)
1 1 L_L
M 1 , 1 1 1 11
103 10^ 105
Reynolds number based on clearance
CD
^0
_ .jmferential
I unwrapped cylinder
Figure 17. - Mathematical model of parallel groove-ndge geometry of infinite horizontal extent
•»—Scavenging or
nonwetted region
Cavity
pressure
Average pres-
sure ridge
leading edge
Averse pres-
Pressure sure ridge
trailing edge
,,-Ambient pressure
Axial distance
High pres-
sure end
C-66-2693
Figure 18. - Region of analysis validity (homogeneous fluid) of a grooved housing viscoseal, Reynolds
number based on clearance or film thici<ness, 1700; modified or reduced Reynolds number, 19;
400-microinch shaft rotational movement during photographic exposure; shaft speed, a)Oo'rpm;
sealed fluid water
High pres-
sure end
mw-'^'
Figure 19. - Viscoseal with grooved housing operating in a regime of water-air mix-
ture (gas ingestion). Analyses does not apply for this discontinuous fluid film
case; 250-microinch shaft rotational movement during photographic exposure;
shaft speed, 5000 rpm, sealed fluid, water
U sin a
c , (3) 1 (3)| y
^ ^
B»
2irR
H
cn
to
to
Circumferential
Leading^
edge-/" y;/^f ^
CD-8606
Figure 22 Pressure distribution over parallel groove geometry (Qualitative
graphical representation from experimental data )
5.1.80
nLcosa
Coordinate along ridge-groove, x
Reference
pressure-
nL sin a -
f^a -^
Ridge'"
Groove
Coordinate across ridge-groove, z
Figure 23. - Pressure profiles (qualitative graphical representation from experimental data)
E-3691
- -o- -TT
"Breakpoint:
-c
-&- rg3:
i^tt :»is^
fe
60 80 m 160 200 400 600 80) low 2(M0 40)0
Clearance Reynolds number
Figure 24, - Experimental sealing parameter for seals 5, 6, and 7 from reference 27.
1 1
Significant /
convective /
inertia y
H>-Actual region effects /
of creeping
flow validity yf Generally accepted by
y// current researchers
Repre-
sentation
of top view
.^^
Schematic
representation
of end view
to where Pj > PQ
Po ^
to
I Section A-A of groove-ridge pattern
U
gg????S^S
^^m^^^^^^
x ^ ^ ^ ^mm^^^M^^^^^^^.^ ^
Figure 27 /Vlodel used m creeping flow solution
(High pressure)
Ux-0.96IJ,
Uj-0.26 0
(Low pressure)
TzMzimzmzzzzzm^^^^^^zzzzzzmzzzzzmzim.
m>.
VWZ^^^^W7Z^7777777'^7Z^77?7Zy7Z^^^
1000 \ or 1000 c
Figure 29 - Example where the creeping flow mathematical model can be considered a valid
physical model Modified Reynolds number (Re ) much less than 1 since clearance c is
much less than characteristic length L
ABSTRACT
Experimental study of the visco seal reveals that the laminar seal performance
is in excellent agreement with the analysis of Boon and TaI. In the turbulent
region, however, the sealing performance exceeds the theoretical laminar pre-
diction and the optimum laminar seal geometry is not optimum for turbulent
operation. A theoretical equation has been obtained for the sealing coefficient
of a visco seal which is suitable for both laminar and turbulent flow. Experimental
results from 13 different seal geometries confirm the utility of the analysis.
Air ingestion has been observed during turbulent operation of each experimental
seal. This phenomenon, however is not fully understood, may become a major
problem in the application of the visco seal.
I. INTRODUCTI ON
The viscoseal, viscosity pump, screw seal, spiral groove seal or bearing are various
names given the device whose working principle is based on the pressure generation
in a viscous fluid, enclosed in a narrow annulus or slit, by means of grooves on a
rotating shaft or plate. The viscoseal, shown basically in Figure I, is an old but
little used device which functions without sliding contact and, unlike other
clearance-type seals, holds promise of zero leakage. The lack of contact, long
reliable life, and the possibility of zero leakage has recently created a renewed
interest in the viscoseal as a sealing element for use in the critical applications
associated with the nuclear energy and space programs (I).
I I. BACKGROUND
The performance of the viscoseal may be expressed in terms of a sealing coefficient
defined as:
A = ^
c Z^
^ t T ( 1 - ?- ) ( s ^ - i )
which is based on the work of Boon and Tal (7) . Figure 2 presents a
comparison of several laminar analyses with data for t e s t seal No. 1.
2. At high Reynolds numbers the sealing coefficient improves (/V
decreases) and is a function of Reynolds number as well as seal geometry.
3. There is a smooth systematic transition from laminar to turbulent
operation and the Reynolds number at which the transition begins is a
function of seal geometry. As T is increased in a seal t e s t s e r i e s ,
having the same a and B, the onset of turbulence begins at lower
Reynolds numbers. These observations suggest that whatever the
theoretical expression for the turbulent A might b e , it should become
equivalent to equation (3) at low values of Reynolds number.
The a n a l y s e s of the turbulent v i s c o s e a l are very limited and are
largely empirical. McGrew and McHugh (5) made use of the Prandtl
mixing length concept and deduced the form of the turbulent sealing
coefficient a s :
Ao=«=l*<=2Re^)-l (4)
Since the experimental constants were determined for a single seal only,
equation 4 is of limited utility in design.
TURBULENT SEALING COEFFICIENT
Making use of the same physical model employed by Boon and
Tal (7) for the laminar a n a l y s i s , a theoretical sealing coefficient for the
turbulent v i s c o s e a l has been derived (2). The procedure utilized in the
turbulent a n a l y s i s was to modify the equations of motion and continuity
used in the laminar a n a l y s i s by replacing instantaneous velocity
components and pressure with mean components and turbulent fluctua-
ting components in order to represent turbulent flow. From these basic
equations the Reynolds equations for turbulent flow were formed.
Following a method similar to that employed by Pai, and representing
the velocity profiles for Poiseuille and Couette flow as power s e r i e s ,
the equations for the Poiseuille and Couette velocity components were
obtained. The flow components in the seal were the same as those
identified in the laminar a n a l y s i s (9). Setting the sum of the flow
components to zero, the turbulent sealing coefficient was obtained a s :
h^h r 11
3
A^ = h + K,.
5 I.
(5)
L 4
where
3 9
(7)
ig - p r t ,
^3 = PV [ y + P ^ I - t )] (8)
and
y + ^ (1 - y ) + g(p - 1)
^4 = ^ [l- ^ + ¥ (9)
I + 6^1 - t )
The terms I , I , I , and I . are functions of the seal geometric parameters ^
X Zi sj T:
P
In the turbulent range the coefficient of resistance for flow normal to
the grooves was expected to be different from the coefficient for flow
parallel to the grooves. In Figure 3 two lines were selected to
represent estimates of the r e s i s t a n c e coefficients for these two flow
directions. The lower l i n e , defined by the equation
F.P.= ^ . (23)
c
5.2.11
NOMENCLATURE
a Axial land width, in.
b ' Axial groove width, i n .
c Radial c l e a r a n c e , in.
D Seal diameter, in.
4jr4
P^p^ _ Re
___E_ Friction parameter
'c
f Resistance coefficient, Figure 3
h Groove depth, i n .
L Active seal lengthy in,
1 Axial threaded length of s e a l , in,
n ^ Number of thread starts
P pressure, Ibf./in.^
q Power l o s s , i n . I b f . / s e c .
Uc/'
l^g _ —j^— Reynolds number based on clearance
t = tan A Tangent of the helix angle
U Surface velocity of the s-eal, i n . / s e c .
X X coordinate in the direction of motion
y y coordinate along shaft axis
z z coordinate in the radial direction
a Helix a n g l e , degrees
h + c
p__ Dlmensionless parameter
^ b
" = ——- Dlmensionless parameter
71 Tl coordinate
A Sealing coefficient a s presented in this work
/^ Absolute v i s c o s i t y , Ibf. s e c . / i n ,
s t coordinate
P Density, Ibf. s e c . ^ / i n . ^
^ Wall shearing s t r e s s , I b f . / i n . 2
Y Dissipation function
Subscripts 1
! Denoting ^ direction
5.2.19
NOMENCLATURE
a Axial land width, in.
b ' Axial groove width, in.
c Radial c l e a r a n c e , in.
D Seal diameter, in.
4ff A
A3
P^p^ _ Re
_—_c_ Friction parameter
c
f Resistance coefficient, Figure 3
h Groove depth, i n .
L Active seal length, i n .
1 Axial threaded length of s e a l , i n .
n Number of thread starts
P pressure, Ibf./in.^
q P o w e r l e s s , in. I b f . / s e c .
Uc^
l^g _ ___i— Reynolds number based on clearance
c yM
t = tan A Tangent of the helix angle
U Surface velocity of the s e a l , i n , / s e c .
X X coordinate in the direction of motion
y y coordinate along shaft axis
z z coordinate in the radial direction
0 Helix a n g l e , degrees
h + c
p _ — c. — Dimensionless parameter
b Dimensionless parameter
a + b
71 TX coordinate
A Sealing coefficient as presented in this work
/^ Absolute v i s c o s i t y , Ibf. s e c . / i n .
S t coordinate
P Density, Ibf. s e c . 2 / i n . ' *
^ Wall shearing s t r e s s , I b f . / i n . ^
Y Dissipation function
5.2.20
REFERENCES
D c h a b a
Test Groove Land Groove Angle
Seal Diameter Clearance Depth Width Width Angle (h+c), b
No. Inches Inches Inches Inches Inches Degrees G (a+b)
T e s t S e a l Number 2B 3B 4B
A ^ (Laminar) 11.73 13.09 12.39 12,86 17.54 14.29 17.20 14.08 12.44 14.10
A £ (Laminar) 10,98 12.86 12.25 11.62 17.90 16.02 14.40 14.92 11.33 11.75
A p
0.93 0.98 0.99 0.90 1.02 1.12 0.84 1.06 0,91 0,83
Ax
(Laminar)
din A E
0.698 0.731 0.695 0.676 0.709 0.709 0.709 0.414 0.430 0.438
dlnRe
(Turb.)
'A
7 L h c ^c
'////A//////A^//////. /
Ui
seal
^"T^ interface a=(|«-y)frD t a n a / n
'D tan «/rH
b= y i r D t a n a/n
^ = ambient pressure
N5
Fig. 1 Basic Elements of a Visco Seal
MTI-2079
A Sealing Coefficient
1-1.
OQ
_ _
.^_^ T ' " T
n
1 1 1 •
o
as
n
en
o 1
O 1
s
j
o 1
i-h
1 Q^
o (
H
n> 1 1
en
ta
3 X
W
s f
00 < m >^
m
H
1 o
c
Q
0
-»•
O <
f\\ ^ m
5 1
0 3- w Q
? 3 J2,. o n
P Q
rt 1 rf" 3
i-h CO
O ^ M »
H 9 \
o - n2 11
O
3
m \ 3"
as
O
en
o
5
O 1 fi Q
W pP J- S
O t
T W 00 Oi 3
o 1 1
5*
1
&. 1
9Z-Z-S
ISJ
Fig. 3 Resistance Coefficient versus Reynolds Number for Pipe Flow
.28
36
32 / / /
/ /
H 2 8 ^ a = 210.15° U 57 Q P,7°J 5.81°/
z
y
^24
LL
Li.
7
/ / /
LLJ
O20 'T"' "
u 1
/
O 16
y
<
UJ 12
uo
\J><^ ^^
r
< i
Rec = 100
D = 1.250"
1
8 1
1
c = 0.003"
7 = 0.5
0
0 4 6 8 10 12 14 16
P, CLEARANCE RATIO
"^"^^
- —
: > - - ^
^yf^
CD >i '• II
O O O r-
O UJ Q
o = ^
11
l\3
P
\ (0
\(J}
\b)
\ ,
-(j|-Va
\
\
V \
CO
CM CD
in
^Y^\°p\ CO
o
o
1 VA-J
CM
N
A
O
II M-l
8 CM oa
\ \ O CO
IW
O DC o
s
o
4J
LLJ o
U c
3
•
1
PM
o
o ©I to o
o
en tn
cd
II
<
9o o
j j
d o
•S^
1
LJJ
(D o
•r-l TJ
^60 'IJ O M-(
0)
o
CO
«
II II II 4) U Ul
Q uK ^ , / /
^
60O
d
•r4 11
r-l
r
CO K
— ^ 0)
r-l d
CO . H
CM O
•H i n
4 J CM
<u .
U .-4
o
<u II
C D ( D ^ ^ O ( D ( D ^ C\J
O
O
1 N 3 D I 3 3 3 0 3 9NnV3S V 00
•r-l
Ul
10
Fig. 8 Theoretical Sealing Coefficient versus Re for a 5.81° and 20.15° with
7 of 0.3 and 0.7 ^
Ul
Fig, 9 Visco Seal Test Facilities
Oil In /^-Bearing feds
Thrust
Plate
®
-Test Spind le /
® ® ® ®
S^bnt
7 Inlet
^-The rmocouptes Support
Block
Torque Arm •Oil Drain
K80 V
^-o- L-75 S
m
1-70 E
•Intercept Pressure
•Supply Pressure
-0A5~*
L
, r 3
Effective Seal Length (in.)
Spindles 1-4B
3/16-^ |*-b
•Spindles 5^6^7
Ex p e r i m e n t a l C u r v e by Screw N 0.1 c = 0 . 0 0 4 2 ()
M c G r e w & McHugh-y ,— O € = 0.1 0=14.50''
-- — • € = 0.6 0= 3.38
y= 0.631
^ u rP d rs""^ 1 i
10 /
^Laminar
^ - ^
Theory ^ ^ " ^
h%J^<^bx.
' ^T^i^^ i
1
Mctjrew & Mcnugh \f
k'**
\
^
\
45-a. i
\ ^ X
10 10^ 10
^
o
OQ
WJvNO
r°°Ti
D
CO i^
n
*w CO
fD t
O
H 1
fD
J ^ Si
O
JSb
b- 'ch
m XI XI
3 © ©
CL
><!
(D !S
l-i 4i». kdJ
P
rt O O O
P
D
'O b^
O
KO
O
^
f i P
as P ,
o o o ]
CO b b ^
n " o o b
hO NO
l-h CO • • ft
/ •
o o D
J • ^
fD ^/l
p ^w
rf =^n
en
I-h
O
M
CO
fD
1^
^
p
Jp> D
P
a.
1 /
^1
8e'2-g
80 1 _ — « _ ^ — 1
Screw a 0 c
D 5 5.8 V' 6.96 0.277 0.00295"
A 6 5.81^ 7.00 0.502 0.00295"
• 7 5,81 ^ ZOO 0.724 000295"
The o r y
5 nn QCfam
7
6 1^:€^ % ^ ^
\ r u , I
10
^^^i ^
^ . ^ ^ f j bo.
- * T_j
r1^^^
^ \
10^ 10
Re.
MTI-2093
O
Re.
Fig. 16 Theoretical and Experimental Sealing Coefficients for Seals 2 and 2B.
A Sealing Coefficient
O
n>
Co
p
m
xs
fD
H
g-
(D
3
03
p
OP
o
o
m
Ml
[2
ro
P
rt
m
i-ti
O
OS
P
O.
UJ
lf7*Z'e
Fig. 18 Theoretical and Experimental Sealing Coefficients for Seals 4 and 4B.
MTL|209
5
Experiment by
\
Smith & M c G r e w & McH ugh
Fuller -J ^ 1
/ V
ocrew a 0 Y c
2 X D \
o 1 14.50^ 3.38 0.631 0.00420'
.? D 2 9.67*" 5.94 0.296 0.00235"
4 A 3 9.67*" 6.16 0.504 0.00265"
0.10
• \
s^ ^X^^ _. 4 967° 5.84 0.692 0.00265"
^ ^ ^
K ng —J
^^^v^S
1^ Smith & Fuller
0.01
\
1 McGrew &
McHugh
^ ^
10' 10
Re,
Fig. 20 Theoretical and Experimental Friction Parameters for Seals 2B, SB, and 4B.
F i g . 21 T h e o r e t i c a l and Experimental F r i c t i o n P a r a m e t e r s for Seals 5 , 6, and 7.
"^
\ .
SSas—•ss^
\ Of
^
V
1 —
^ ^ ^
1 i 1
\ . — 2E
%^ »
r"^^^
__
1
M n as?s exDerimenia
. . 1
curves IB "X^ X
-A,
~ Scir e w ' a ^ ' y c 5D
h%Sv ^ ^
\J
SECTION 6
bottom support bearing to dissipate heat and prevent excessive lower sup-
port bearing temperatures due to soak back from the high temperature
sodium in the test vessel. The test journal was mounted and keyed to the
bottom end of the test shaft. The test bearing was mounted in a housing
in the test vessel as shown schematically in figure 1(b).
Liquid sodium at 400° F was introduced to the test vessel and heated
to the desired test temperature by means of an induction heater. The
induction heater coil around the test vessel (fig. 1(b)) does not come into
physical contact with the outer surface of the test vessel and therefore
does not inhibit its free swinging motion.
A drain was provided at the bottom of the test vessel to facilitate
draining of contaminated sodium.
Sodium Supply System
The sodium supply system was a noncirculating, once-through system
composed of the following major components: the supply tank, the supply
line, the sodium filter, the filter by-pass line, and the control valve.
The supply tank was sized to contain 20 gallons of sodium and was
equipped with a fill valve, dual thermowells, diffusional cold trap, p r e s -
sure transmitter, vapor trap, and a sodium supply line. The liquid metal
was supplied to the bearing test rig by means of a differential pressure be
tween the tank and the test rig. A diffusional cold trap was provided to
control the sodium oxide content during operation of the sodium system.
The cover gas for pressurizing and venting the supply tank was passed
through a wire-mesh demister (vapor trap) which prevented any sodium
vapors from getting into the cover gas system.
A micrometallic filter was provided in the supply line to the bearing
rig permitting continuous filtering of all supply sodium. The porous
micrometallic filter unit was so designed as to permit replacing the filter-
ing element.
Bearings and Journals
Hydrodynamic bearings of four configurations were tested fully im-
mersed in liquid sodium. Bearings with two axial grooves, three axial
grooves, a herringbone groove journal with a plain bearing, and tilting
pad bearings with three pads were evaluated. One of the three groove
bearings was run with an axial flow pump attached to the test shaft, p r e s -
sure feeding the test bearing through a hole in the test journal.
The bore and length of the bearings in all cases were nominally
% inches. The journal outside diameter and bearing inside diameter
were machined to a 4 to 8 microinch finish, r m s .
A three pad configuration was chosen for the tilting pad bearings
because it affords greater load capacity than a configuration of more than
three pads. Also, it is much easier to accurately aline all pivot points
on the common pivot circle center with a three pad configuration than one
with more than three pads. The load was applied symmetrically between
two support points (pivots).
The bearing and journal materials were: Stellite Star J, Mo-0. 5Ti,
Hastelloy X, Inconel, and titanium carbide (K184B). The composition
and hardnesses of these materials a r e given in table I.
Instrumentation
Dual chromel alumel thermocouples were attached to the test bearing
back and into the liquid sodium in the test vessel. The induction heater
was controlled by one of the dual thermocouples in the sodium bath.
Two capacitance probes, which measured the movement of the test
vessel during a test run, were mounted outside of the test vessel on the
vessel cover, 90° from each other. The signal from the probes was fed
through displacement m e t e r s to an x-y display in an oscilloscope where
the actual pattern of motion of the test vessel could be observed. The
orbital frequency of the test vessel was measured by means of a f r e -
quency counter.
Shorting probes were used in the test vessel to indicate the level of
the sodium. These probes would short out when sodium came into con-
tact with them, thereby either lighting a level indicator light or closing
the main sodium supply valve automatically.
Test shaft speed was measured with a magnetic pickup head mounted
in close proximity to a six-toothed gear on the test shaft. The signal from
the pickup was displayed on a four-channel frequency counter.
PROCEDURE
Pretest Preparation for Two and Three Axial Groove Bearings
P r i o r to each test run, the test bearing was pressed into its housing
with a slight interference fit. The bearing was then machined in place to
a predetermined inside diameter at room temperature that would result in
the desired bore size at test temperature. Nine bore gage readings, each
accurate to within 0. 001 inch, were averaged and used as a measure of the
bearing bore. The outside diameter of the mating journal was then ground
to a size that would result in the desired clearance for the test bearing.
To insure a minimum amount of runout, the journal outside diameter was
machined to within 0.0002 inch concentricity with its inside diameter.
The test bearing housing was then assembled into the test vessel
which was carefully raised into position around the test journal on the
shaft by means of the lower air cylinder.
Pretest Preparation for the Tilting Pad Bearings
The radii of the three pads of a bearing assembly were checked after
delivery from the vendor to insure an accurate geometry. The pads were
assembled into an annular housing by means of a threaded pivot and nut
arrangement. Desired preload was obtained by adjustment of the threaded
pivots until the bearing surfaces of the pads made intimate contact around
a presized set-up plug. The plug was then removed and the tilting pad
housing was assembled in a manner like that for the axial groove bearings.
General Pretest Preparation
After the bearing and journal were assembled and the test vessel
raised to its run position with the upper and lower gas bearings turned
on, the test vessel was filled with alcohol and drained as a final cold
cleaning procedure. The test vessel was then purged with argon and a
cover gas of argon was supplied to the test vessel throughout the test.
The test vessel was preheated to 500° F and liquid sodium at about 400° F
was introduced into the test vessel through the inlet port under about
5 pounds per square inch p r e s s u r e from the 20-gallon supply system.
Sodium flowed into the test vessel until it made contact with the liquid
level probe which automatically closed the main sodium supply valve p r e -
venting overfilling.
T - i^ A ^ L N
^ 120 C^
and the theoretical turbulent torque calculated from the equation given by
Smith and Fuller (ref. 6) where
T^= T^ (0.039 Re^-^"^)
The theoretical torques a r e for a full circular bearing with zero eccen-
tricity, so they a r e only rough approximations for a tilting pad bearing,
and only for zero preload. Bearing torque values indicated that turbulent
flow conditions prevailed over the greater part of the speed range, due
primarily to the low sodium viscosity. The experimental transitional
speed, where the bearing passes from the laminar to the turbulent regime,
occurred at a higher speed than that predicted by theory (critical t r a n s i -
tional speed, Nrp). This may be due to the fact that this transition occurs
over a range of speeds rather than at one definite speed. Similar results
with the transition speed were reported in reference 7.
No abrupt increase in torque readings was noticed when turbulent flow
conditions were approached. The rather gradual increase in torque with
speed made it difficult to specify exactly the speed at which full turbulence
was attained.
Bearing Instability
The bearing instability of principal concern here is half-frequency
whirl. Figure 6 shows oscilloscope traces of bearing motion obtained with
a two axial groove bearing in sodium at 500° F with a 10-pound radial load.
At 4100 rpm the trace indicated stable bearing operation (fig. 6(a)), When
the speed was increased to 5000 rpm, however, the increase in attitude
angle was sufficient to sustain half-frequency whirl. The whirl pattern
observed on the oscilloscope screen is shown in figure 6(b). K the bearing
is allowed to operate unstably, the supporting film between the bearing and
journal soon breaks down, and the bearing eventually fails.
One of the most undesirable characteristics of the two and three groove
bearings was their instability. The result of such instability is graphically
shown in figure 7. Test bearing J - 2 is shown in this figure after 257 min-
utes of operation at 11 psi unit bearing load in 500° F sodium under half-
frequency whirl conditions. The excessive wear shown is the result of
unstable bearing operation.
Of the 14 two and three groove cylindrical bearings tested, five showed
excessive wear due to half-frequency whirl and one seized because of this
instability.
One of the three groove bearings (bearing number M-9, table II) was
loaded sufficiently at each speed throughout its evaluation to keep it run-
ning stably. After 290 minutes of running time in 500° F sodium, the
bearing was removed and no measurable wear was present on either the
journal or the bearing. The maximum load on this bearing was 26, 7 psi
at a maximum speed of 10 000 rpm indicating that a three groove cylin-
drical bearing will run successfully in sodium if properly loaded to sup-
p r e s s half-frequency whirL
Theory indicates that herringbone groove bearings operate at consid-
erably lower attitude angles than do smooth bearings resulting in more
favorable stability characteristics. This type of bearing assembly was
indeed more stable than the two and three groove bearings and plain
journal assemblies. However, the herringbone groove bearing assembly
did show evidence of half-frequency whirl at low load conditions. A more
judicious design of the herringbone groove journal might lead to a bearing
that would be stable even at zero load conditions. Such a bearing has been
run at zero load in air up to 60 000 rpm without any evidence of half-
frequency whirl (ref. 8).
Figure 8 shows the relative stability of the four different cylindrical
bearing configurations tested. The two and three axial groove bearing
configurations were the least stable of the four since they required the
highest load at any specific speed to maintain stable operation. The plain
bearing with a herringbone groove journal was the most stable of the four
since it required the lowest load at a given speed to keep it running stably.
Axial grooved bearings appeared to require a linear increase in load with
speed to maintain stable operation whereas the herringbone groove b e a r -
ing was stable at 25-pound load at speeds of 7000 to 10 000 rpm.
Experimental data on the threshold of instability of two and three axial
groove bearings correlated well with the theoretical curves reported in
reference 9. Figure 9 shows theoretical curves of Sommerfeld number
plotted against the dimensionless critical rotor mass for a 100 partial
bearing and full circular journal bearing. The data points for the two
axial groove and three axial groove bearings generally fall between the
curves indicating good correlation.
The procedure for determining the threshold speed is to calculate the
dimensionless rotor mass
/iDL
Enter figure 9 with this value and determine the corresponding Sommerfeld
number using the appropriate bearing curve. The speed corresponding to
this Sommerfeld number is the threshold speed, that is, the rotor speed at
onset of instabilriy.
The tilting pad bearings were the most stable of the four configurations
tested. However, half-frequency whirl was observed and measured with a
frequency counter on bearing number T-3 at 12 000 rpm and zero load.
Tilting pad bearings T-3 and T-2A1 were run up to 11 000 rpm at zero load
without exhibiting any half-frequency whirl instability.
Although not tested at zero load, the remaining three bearings of the
tilting pad group showed good stability al light loads of 4. 5 psi to speeds
of 11 000 rpm.
Material Compatibility
Table V lists the bearing and journal material combinations that had
good wear and seizure properties in sodium to 800° F.
Stellite Star J material mated with Hastelloy X, titanium carbide
(K184B), or Inconel showed the best wear and seizure properties. Also
titanium carbme |K184B' mated with Mo-0. 5Ti, showed excellent promise.
Materials such as Has.eli'»y X and Inconel having high nickel content,
were prone to catastrophic seizure wh^^n paired in a bearing and journal
combination.
Figure 10 shows the rehults ol a seizure of a three axial groove b e a r -
ing due to an incompatible bectni.g and journal material combination. The
bearing material was Mo-0 5Ti and the journal, Hastelloy X. The surface
of the journal shows a gallmg typical ol this type of failure.
Another example oi the r^sui s of pairing 'ncompatible materials is
shown in figure I L The bearing material was Inconel and the journal
material was Hastelloy X, boih high in mckel content. Seizure resulted
in galling of the pads and iournah The beaiing radial clearance at the
pivots of 0. 0003 inch would not allow many particles to pass through the
bearing without initia+ing suriace damage. With a poor combination of
materials, galling, severe surface damage, and possible seizure quickly
follow the initial surface- damage
Figure 7 shows a*" example ^t manng materials with good compatibil-
ity. The combnid+ic^n ol Stt-'llne Siar J and Hastelloy X material showed
excellent seizure resistance S'nre «his bearing and journal combination
did not seize even aH.er ptol -^ngrd operuion with whirl present which pro-
duced the exce&sivf' w-'ar shov^n in *he figure Figure 12 shows a photo-
graph of the pivot arrangement used on the tilting pad bearings: a sphere
against a flat. Figure 13 shows a typical example of surface damage to
the pivot that was observed in some tests, even after very short runs at
light load. This damage occurred with both the titanium carbide (K162B)
and the Stellite Star J pivot materials which were always mated against
themselves. As shown in table IV, bearing number T-1 ran for 175 min-
utes with a maximum load of 15 pounds on the pivot producing a Hertz
s t r e s s of only 52 600 psi which was sufficient, however, to cause slight
pivot surface damage.
SUMMARY OF RESULTS
A series of hydrodynamic journal bearing experiments was run in
sodium at 500° and 800° F at speeds to 12 000 rpm and loads to 70 pounds.
Four different configurations were tested; cylindrical bearings with two
and three axial grooves, plam cylindrical bearings with herringbone groove
journals, and tilting pad bearings with three pads. The bearing bore in all
cases was 1. 5 inches and all bearings had a length to diameter ratio of 1.
The following results were obtained:
1. The tilting pad bearings were the most stable. Following in order
were: (a) plain cylindrical bearing with a herringbone groove journal, (b)
three axial groove bearing, pressure fed from an axial shaft pump through
a hole in the journal, and 'c'" three and two,axial groove bearings.
2. Stellite Star J material mated with Hastelloy X titanium carbide
(K184B), or Inconel showed the best wear and seizure properties. Also
titanium carbide i'K184B) mated with Mo-0. 5Ti showed excellent promise.
Materials having high n-ickel ronten*-, such as Hastelloy X and Inconel,
6.1,17
41.1 V
TTDC^
MN'DL / R V
Sommerfeld number.
w vc^
T^ Petroff's torque for laminar flow, zero load: T, = M7r^D^LN/120C^
T^ turbulent torque: T^ = T^(0. 039 Re°- ^'^) (Smith and Fuller, ref. 6)
W bearing load, lb
W load due to rotor mass, lb (W^ = Mg)
Jc3m
-1— dimensionless critical rotor mass
o
p-
m
TABLE I. - NOMINAL COMPOSITION AND HARDNESS OF BEARING AND JOURNAL MATERIALS
Stellite Star J C-62 .. 2.5 32 40.5 3 .... 2.5 .... 17 ... .. .... 2.5
Mo-0.5Ti 8-87 „ »-« -- . . . . - .... .... .... .... 0.5 .. 99.5 ...
Inconel B-75 toB-95 -- 0.04 15 .... 7 0.35 78 0.20 .... . . . -- .... ...
Titanium C-67 3 3 .... 40 .... . . . 50 4 "'"
Carbide {K-184B)
(nickel bonded)
E-3570
M-1 M o - 0 5Tl K-6 TiC 2 Groove 500 0.0010 3000 to 8000 4 5 to 2 0 . 0 470 At v a r i o u s speed B e a r m g and j o u r n a l showed m o d e r a t e w e a r
and load conditions due to half frequency w h i r l Increasmg
load stabilized b e a r m g . B e a r m g did not
s e i z e Ran d e l i b e r a t e half-frequency-whirl
operation at m o m e n t a r y i n t e r v a l s
M-7 Mo-0 5Ti K-B TiC 2 Groove 800 0.0010 4000 to 11 000 4 5 to 2 0 . 0 703 None p r e s e n t B e a r m g showed m o d e r a t e w e a r at both e n d s
m d i c a t m g the t e s t v e s s e l had cocked. Bear-
ing did not s e i z e . Sodium had m i g r a t e d to
top g a s b e a r m g of t e s t v e s s e l .
J-1 S t e l h t e Star J HX-B Hastelloy X 3 Groove 500 0 0013 5000 22.0 5 Not s e t up t o Contammant p a r t i c l e s c o r e b e a r m g a n d
observe caused seizure almost immediately after
s t a r t of t e s t
J-2 Stellite s t a r J HX-C Hastelloy X 3 Groove 500 0 0012 5000 to 7000 11 0 257 At 7000 r p m and B e a r m g and j o u r n a l showed s e v e r e w e a r but
11 0 p s i load no s e i z u r e r e s u l t e d . Held b e a r i n g at 7000
r p m and 1 1 0 p s i d e l i b e r a t e l y for 90 m m u t e s
to o b s e r v e r e s u l t s of h a l f - f r e q u e n c y - w h i r l
operation
J-3 Stelhte Star J HX-D Hastelloy X 3 Groove 500 0.0011 5000 to 7000 9.0 260 Not set up to B e a r m g and j o u r n a l showed heavy w e a r
observe probably due to u n s t a b l e operation at light
load B e a r m g did not s e i z e .
HX-B Hastelloy X J-5 Stelhte 3 Groove 500 O.OOU 5000 11.0 30 Not set up to B e a r m g s e i z e d a f t e r 30 m m u t e s a p p a r e n t l y
Star J observe due to c o n t a m m a n t p a r t i c l e .
HX-D Hastelloy X J-6 Stelhte 3 Groove 500 0.0015 3000 to 7000 4. 5 to 8. 9 460 Not s e t up to B e a r m g and j o u r n a l showed heavy w e a r but
Star J observe did not s e i z e Light l o a d s probably c a u s e d
the b e a r m g to r u n unstably.
HX-5 Hastelloy X J-3 Stelhte 3 Groove 800 0.0008 5000 Not r e c o r d e d Not Too s h o r t a run B e a r m g s e i z e d i m m e d i a t e l y , probably due
Star J recorded to o b s e r v e to c o n t a m m a n t p a r t i c l e . B e a r i n g showed a
s c o r e m a r k completely a r o u n d on one end.
^By observation of oscilloscope pattern as the test vessel was moved by hand, to the limits of its outward motion, in a circular orbit about the stationary journal.
E-3570
Bearmg Bearmg Journal Journal Bearmg Test Measured'' Journal speed Ohitload Total test Observed Remarks
number material number material type temper- radial range tested, range tested, time, mstability
ature, clearance rpm psi mm
°F at t e s t
temperature,
m.
K-3 TiC J-7 Stelhte 3 Groove 500 0.0017 5000 4.5 315 Not set up to B e a r m g and j o u r n a l showed-heavy w e a r but
Star J observe did not s e i z e , B e a r m g apparently w a s
o p e r a t m g unstably b e c a u s e of light load.
M-9 Mo-0.5Ti J-U Stelhte 3 Groove 600 0.0010 5000 to 10 000 1 1 . 0 t o 2 6 . 7 290 At v a r i o u s speed B e a r m g and j o u r n a l showed discoloration
Star J and load conditions but no m e a s u r a b l e w e a r , B e a r m g w a s
purposely loaded sufficient to keep it rvm-
nmg stably. Also r a n d e l i b e r a t e half-
f r e q u e n c y - w h i r l operation at m o m e n t a r y
mtervals.
M-12 M o - 0 5Tl J-1 Stelhte 3 Groove 800 O.OOU 5000 to 10 000 13. 3 to 22. 2 355 None p r e s e n t B e a r m g w a s loaded sufficient to keep it
Star J runnmg stably, however, b e a r i n g s e i z e d at
8000 r p m during shutdown 1/8 mch gall
m a r k a r o u n d b e a r m g possibly c a u s e d by
c o n t a m m a n t p a r t i c l e , o t h e r w i s e light w e a r .
M-10 M o - 0 5Ti J-2 Stelhte 3 Groove 800 0.0010 4000 to 7000 15. 5 t o 1 7 . 8 242 None p r e s e n t B e a r m g s e i z e d at 7000 r p m and 17.8 psi
Star J load. Two s m a l l w e a r a r e a s on each end of
b e a r m g m d i c a t m g the t e s t v e s s e l had cocked.
M-11 Mo-O.STl HX-7 Hastelloy X 3 Groove 800 0.0013 Not r e c o r d e d 0 Not Present B e a r m g s e i z e d i m m e d i a t e l y at 0 load with
recorded immediately defmite mdication of half-frequency w h i r l .
W e a r on both ends of b e a r m g w a s o b s e r v e d .
HX-7 Hastelloy X J-8 Stelhte 3 Groove 500 0.0012 4000 to 9000 1 3 . 3 to 22. 2 647 At v a r i o u s speed A shaft s c r e w pump forced sodium mto the
Star J a n d load conditions b e a r m g through a hole m the j o u r n a l . B e a r -
mg showed s m a l l w e a r a r e a m loaded zone
probably due to d e l i b e r a t e half-frequency-
w h i r l operation at m o m e n t a r y m t e r v a l s .
By observation of oscilloscope pattern as the test vessel was moved by hand, to the limits of its outward motion, in a circular orbit about the stationary journal.
E-3570
TABLE in. - RESULTS OF PLAIN CYLINDRICAL BEARING AND HERRINGBONE GROOVE JOURNAL TESTS IN SODIUM
Bearing Bearing Journal Journal Bearing type Test Measured^ Journal speed Unit load Total Groove angle Number Width Depth Observed Remarks
number material number material temper- radial range tested, range tested, test (measured of of of instability
ature, clearance rpm psi time, from a per- grooves grooves grooves,
°F at test min pendicular and ^ and in.
temper- to the journal lands lands,
ature, axis) ((!) in.
°F
MP-3 Mo-0. 5Ti K-14 TiC Plain bearmg. 500 0.0013 4000 to 12 000 ''O to 17. 8 490 33° 20 0,064 0.0014 At various Light wear on bearing
Herringbone (on lands) speed and and journal, due to
groove journal load con- deliberate half-frequency-
ditions whirl operation at mo-
mentary intervals.
MP-4 Mo-0.5Ti K-A TiC Plain bearing. 800 0. 0008 5000 to 7000 ''O to 20. 0 90 33° 20 0.064 0.0014 At 5000 Bearing seized due to
Herringbone (on lands) rpm and overload caused by
groove journal 0 load heater expansion pushing
test vessel increasing
radial load.
By observation of oscilloscope pattern as the test vessel was moved by hand, to the limits of its outward motion, in a circular orbit about the stationary Journal.
Zero load only momentary since ha If-frequency-whirl occurred.
#
E-3570
Bearing Bearing Journal Journal Bearing Test Measured^ P r e l o a d J o u r n a l speed Unit load Total Observed Remarks
number material number material type temper- radial coeffi- r a n g e t e s t e d , r a n g e t e s t e d , t e s t instability
ature, clearance cient rpm psi time,
°F at t e s t min
temperature,
°F
T-1 Inconel J-10 Stellite Star J 3 Tilting 500 " 0 ^ = 0.0021 0.19 5000 t o J l 000 4 . 5 to 6 . 7 175 None p r e s e n t Very little w e a r o b s e r v e d with m o s t w e a r on
pad unloaded pad. Pivot s u r f a c e s showed slight
''C = 0.0017 surface damage.
T-2 Inconel J-12 Stellite S t a r J 3 Tilting 500 Cp = 0 . 0 0 1 1 0.36 5000 t o 9000 4 . 5 to 3 1 . 1 311 None p r e s e n t B e a r i n g s e i z e d a t 9000 r p m a t 3 1 p s i b e c a u s e
pad of overloading. Loaded p a d showed m o s t w e a r .
Cp = 0.0007
T-2A1 Inconel J-13 Stellite Star J 3 Tilting 500 C = 0.0028 0.50 3000 to 11 000 0 to 8 . 9 590 None p r e s e n t Unloaded pad showed m o s t w e a r indicating
pad possibly insufficient p r e l o a d . J o u r n a l showed
Cp = 0.0014 light w e a r . P i v o t s showed mating s u r f a c e
damage.
T-3 Hastelloy X J - 1 4 Stellite Star J 3 Tilting 600 C = 0.0036 0.72 2000 t o 12 000 0 to 1 7 . 8 1013 At 12 000 r p m V e r y little w e a r on a l l s h o e s . P i v o t s showed
pad and and 0 load mating surface d a m a g e . C l e a r a n c e a t 800° F
800 Cp = 0.0010
w a s not appreciably different from 500° F v a l u e .
T-IA Inconel HX-8 Hastelloy X 3 Tilting 800 C = 0.0010 0.70 5000 to 8000 4.6 175 None p r e s e n t B e a r i n g s e i z e d when higher load than 4 . 5 p s i
pad w a s a t t e m p t e d , probably due t o tight c l e a r a n c e
C = 0.0003 and incompatible m a t e r i a l s . W e a r evenly d i s -
tributed between p a d s .
By observation of oscilloscope pattern a s the test vessel was moved by hand, to the limits of its outward motion, in a circular orbit about the stationary journal.
C = bearing radial clearance before preload.
*^C = bearing radial clearance at pivot location after preload.
TABLE V. - BEARING AND JOURNAL MATERIAL
Inconel
Titanium carbide
(K-184B)
— '••• ——.•..—»••.»••.•.
CS-41286
(b) Detailed view of test-bearing installation.
Figure 1. - Concluded.
6.1.28 I
•I
C-66-2144
CS-41282
C-66-2149
CS-41284
•P
Figure 3. - Plain bearing (MP-3) versus Herringbone groove journal (K-14) after test.
I 6.1.29
I* JOURNAL
0 INCH 1
1.1II 11111
PAD A PAD B PAD C
C-65-2416
CS-41283
Figure 4. - Failure due to overloading. Tilting pad bearing (T-2) versus journal (J-12).
61—
Tf [0.039 ReO-57]
4 —
ITH AND FULLER)
2 —
UTT^D^LN
(PETROFF)
120 Cr
BEARING •*
TORQUE, .08
( I N . LB)
NSITIONAL SPEED
.06
TTD C ^ V C ,
04
TRANSITIONAL
02 SPEED ( N j )
(TAYLOR
CRITERION)
01
2000 4000 6000 1 0 , 0 0 0 20,000
JOURNAL SPEED, RPM CS-412J
(a) Zero load.
Figure 5. - Comparison of experimental and theoretical friction
torques for a 3-pad tilting pad bearing. Lubricant, 500° F
^ sodium; bearing material, HastelloyX; Journal material,
Stellite Star J; nominal diameter, 1.5 inches; nominal
length, 1.5 inches; radial clearance, 0.0036 inch.
6.1.30
.4
.039 R e ° - " ]
.2
ND FULLER)
.1
.08
BEARING
TORQUE, 06
( I N . LB) • (PETROFF)
.04
ONAL SPEED
V D
o «. V c w «
02
Oil
1000 2000 4000 6000 10,000 20,000
JOURNAL SPEED, RPM
CS-41293
(b) 10 pound load.
Figure 5. - Concluded.
60
3-GROOVE BEA
50
2-GROOVE BEARING
r-3-GROOVE BEARING
40
/ PRESSURE LUBRICATED
STABLE
BEARING OPERATION
LOAD, 30 ,-HERRINGBONE
LB GROOVE JOURNAL
20
I/ UNSTABLE
'/ OPERATION
10
1
.8 THEORETICAL FULL
CIRCULAR BEARING
.6
THEORETICAL
100 DEG BEARING
2 —
1 ± 111
02 .04.06 .1 .2 4 .6 2 4
HN'DL
S =
W {§;) CS-41287
C-66-2148
CS-41281
Figure 10. - Seizure due to incompatible material combination. IV\o -0.5Ti bearing.
•P
(M-11) versus Hastelloy " X " journal (HX-7).
MTI-2070
I 6.1.33
JOURNAL
0 INCH 1
lllllllll
C-66-2151
CS-41279 PAD A PAD C
Figure 11. - Tilting pad bearing seizure due to incompatible material combination. Inconel
bearing (T-IA) versus Hastelloy " X " journal (HX-8).
0 1
lllllllll
.—SPHERICAL PIVOT
INCH
CS-41285
Figure 12. - Spfierical pivot and flat mating surface in pad. (Titanium
carbide (K162B)) versus (Titanium carbide (K162B)).
< • 5-41278
SURFACE DAMAGE TO SPHERI
CAL PIVOT. (TITANIUM CARBIDE
(K162B)
Screw seals, rotating channel and slinger squeeze seals were investigated in the
dynamic seal program. Interface instability was identified as one of the major
problems. It was resolved for the case of the slinger squeeze seal, which was
successfully tested in potassium.
DISCUSSION
Alkali metal bearings and seals comprise two of the basic components necessary
for reliable operation of Rankine cycle turbomachinery for space. As a result
of the high temperatures and radiation fields in nuclear space power systems,
and the general complexity of multi-fluid systems, present turbogenerator
bearings and seals employ the primary thermodynamic fluid of the system in
their operation. SPPS is engaged in the development of alkali metal bearings
and seals for incorporation in a full turbogenerator system.
A typical 350 KWe turbogenerator for space power conversion - General Electric
concept - is shown in Figure 1. It consists of a six-stage turbine driving
an axial gap generator. Each component is supported on two radial bearings
and a thrust bearing. The two components are connected by an elastic coupling.
The design has been optimized for minimum rotor weight and distribution of
this weight over four equal size pivoted pad bearings. Each component,
turbine as well as generator, has its own thrust bearing primarily to fix the
rotor position accurately, in relation to its stator. Thrust loads are mini-
mized in the turbine by proper fluid design, and in the generator by carefully
positioning the rotor in the middle of the stator. The generator cavity is
open to space and is based on the availability of zero leakage seals. The
Within SPPS, bearing and seal development is being carried out under
four different programs, three of which are basic, while the last one will
combine the knowledge gained in the three previous ones. Of the three basic
programs, the first one is devoted to the study of bearing stability which
of bearing rotor response behavior. This program is carried out under NASA
o
contract and uses water at 120 F as a lubricant (Appendix A and B).
suitable bearing materials. Under this program (Contract No. NAS 3-2534,
Appendix C), accurate physical properties have been established for these
is being carried out in two highly sophisticated test rigs, one operating in
a vacuum and the other one in a potassium atmosphere. This program favors
extremely hard and stable refractory materials primarily Tungsten and Titanium
Carbides using special binders. This program will not be further discussed
here because of the limitation in time. It covered three years and is close
time.
The third program was carried out under Air Force contract and was
sealant (Appendix D). Its main application will be in the electric generator
for space.
and has as its objective the simulation of the bearing and seal arrangement
in a space generator. This presents the more difficult problem than the
space turbine due to rotor and containment materials involved. This program
vacuum (Appendix E ) .
The first portion of the Bearing Stability Program was carried out in
1961 through 1963 under NASA Contract No. NAS 3-2111. The test rig in this
time was designed and manufactured by SPPS. Testing and test evaluation,
however, was carried out in the Bearing and Lubricant Center in Schenectady
by Gerry Fox and his staff. The program was continued in 1965 under a new
regime. Out of this program, the pivoted pad and the three-lobe bearing
emerged as the outstanding candidates for stable operation up to high
speeds. These two bearing types are being further investigated under the
the four-pad pivoted pad bearing by running tests with one or two pads
fixed.
motions are monitored by four Bently gages per plane, arranged in four
separate planes.
have been made to measure bearing loads and to establish the load vector
by size and angle during testing up to 500 cycles per second. The Bearing
and Lubricant Center has had good success in making such measurements by
using strain gaged buttons for shaft speeds up to 7,000 rpm. For high
selected due to their high spring constant and their extreme capability of
data handling system for data taking. In high frequency testing, the normal
comes to the evaluation of test results, both of these methods are inadequate
This is done by using peak-to-peak detectors for forces and displacement and
average level detectors for equilibrium position of the shaft. The other
inputs, such as speed, torque, flows, temperatures and pressures, are straight
shaft supported by two test bearings approximately 12 inches apart. The shaft
power supply through a quill shaft to speeds of 30,000 rpm. Side loads are
operated. The lower bearing assembly can be adjusted both transversely and
motion and force transmitted to the bearings are made in four planes, located
Figure 4, The gages have Teflon tips to protect them against water, as
gage readings, the shaft is silver plated in the zones of gage measurement
(Figure 5 ) . A measuring accuracy within + 50 micro-inch is an objective
Figures 7 and 8 show the test rig under testing conditions at room
(closed).
five compartments, one handling the lube flow, pressure and temperature
control, another for power and frequency remote control of a variable fre-
quency power supply arranged in another building. The third panel holds a
angle and force measurements. The fifth panel holds the equipment required
present test program. It has four pads with gimbaling points at 55 percent
partial arc. Figure 11 shows the calculated load carrying capacity of this
bearing with oil and with liquid metal as lubricants. This graph brings out
the weakness of the low viscosity bearing resulting in a load carrying capa-
generated by G,E,, Bearing and Lubricant Center, Another one was published
by MTI under Air Force contract, SPPS has used both of these programs and
has generated an additional program which combines VAST and LSE , offering
test results with those predicted will be possible. Figure 12 shows a typical
in turn is used to calculate fluid film forces as shown in Figure 13, Figures
two different levels of shaft unbalance in the easy-fluid test rig, calculated
by using the MTI computer programs. In these curves, the occurrence of critical
speeds can be identified with the high amplitude zones, and it can be seen
that for the larger value of unbalance, the rotor amplitude actually exceeds
and varying support structure stiffness, both of which are difficult to pre-
duct analytically.
seal, Figure 16. Considerable testing in water was done to identify the
vapor pressure, the molecular diffusion from a free liquid surface to the
the liquid. Any migration of liquid along stationary walls toward the
vacuum side results in re-ingestion of the liquid into the seal before it
was accomplished in the liquid metal bearing and seal facility in Figure 21,
diameters of interest.
tank during 173 hours of testing the dynamic seal. Most of this was accumu-
The liquid metal bearing test rig, shown in Figure 25, is similar to the
previously discussed water test rig, in that two test bearings support a
6.
vertical shaft driven by an electric motor. The rotor weighs 70 pounds and
Along with the very accurate manufacturing requirements of the bearing, the
design requirement, since the rotor of the generator must be made from magnetic
steel, the casing from stainless steel and the bearing is fabricated from
radial springs, radial faced teeth, and conical faces as shown in Figure 27,
the latter approach having been selected as the most practical. The bearing
detail is shown in Figure 28, and comprises four segmented pads, retained by
The pads are weight balanced so that their centers of gravity coincide with
the pivot contact points to minimize fretting motions. The entire bearing
potassium. A side load of 200 to 300 pounds can be applied to the rotor
by generating a magnetic field with an induction coil. The test rig has
other new features such as a high temperature induction gage measuring the
oxide ring which is part of the rotor. The lower shaft end is exposed to
vacuum and sealed off by a zero leakage seal modified for a low power require-
startup. Figures 29 and 30 show rotor and stator hardware. The test rig is
presently being assembled and will undergo its easy-fluid checkout in
November 1966. Liquid metal testing is scheduled for the first quarter
of 1967.
Conclusions
even the problems involved were poorly identified, much knowledge and much
available and will eventually reach close agreement with actual test condi-
bearings and by developing seals which will combine the zero leakage feature
bearing not because of its resistance to partial frequency whirl, but also
1967.
APPENDIX A
The test rig shall be designed with the facility to alter and to control
accurately the following input variables:
1. Shaft speed
2, Bearing loads
Unit directional
Rotating
3. Inlet lubricant pressure
4. Inlet lubricant temperature
5, Lubricant flow rate
1. Bearing load
2. Inlet lubricant pressure
3. Inlet and outlet lubricant temperature
4. Lubricant flow rate
5. Shaft speed
6. Friction torque
7. Shaft position in the bearing with respect to time
8. Eccentricity
APPENDIX A - (Cont'd)
The test lubricant for this test phase shall be distilled water.
APPENDIX B
Specific Requirements
( Length .
Diameter 1.0
Variables Range
The Contractor shall provide a friction and wear test facilities capable
of maintaining a vacuum in the range of 10~^torr pressure or better. This
facility shall be designed to alter, control, and read out accurately the
following variables:
APPENDIX C - (Cont'd)
3. Induced torque,
4. Test environment
pressure, torr
potassium level.
Objectives
Requirements
Evaluation
1. Preliminary experiments
3. Thermal-cycling t e s t
Test Conditions
The test conditions are potassium at 1400 F and 15 psi on one side of
the seal, a pressure of 10" mm Hg on the other side, and a shaft speed of
24,000 rpm.
APPENDIX E
Objectives
Specifications
Temperature Level i n B e a r i n g s
and S e a l s 600°F
n -^
III
MTI-2129
6.2.24 I
•I
ElKltlC D i m MOTOI
If,
TOIOUE ftCIUPS ouiit mm
, •• /
UMIAIANCI MiCS
%i TEST SHAfT
EOiCI lUIIOItS
r^
I
lOAOfft l E A i m S
INDUCTANCE 6AUCE DISTANCE OETECTOfiS
I -.jum^-
I
AOJUSIAIIE SUIASSEMIIY
HI »Al
TEST IEAtlN6
» i •» All IN tPROIES)
</
wATti o n ^
Displacement
Sensor
Holder (4)
iP i (
•I ®
(4)
strumentation / T \
Locations \U
Side Load
Probe Lead
i-cickin;.; Cij] 'i!.>t;
jBjMlgfeiJ^
^^H^^<*^
•.
^^^^^?<
-
^^^^^^^^^
^^BP**- -^ Purj^f;
^ ^" , ^' Air In
Air Qui: (Typ . )
Teflon CsTi Pro tec t: i vc Curt a i n
\ \^ s ^ ^.
(HI-213Z
I ().1.11
Y-J-^fwifl^W*- . .*tT,r?;^'-
Test Shaft
Test Journal Diameter (Lower End)
"titiirrriiiir"-'"f;r '1
;;v:jrfTrjrn|nijii,|r,;[M,j,,pp|if^rf;Mi,n
["•:K.. ° '"'" ^ ^ 3 * 5
! •
MTI-2134
00
ML.
Fig. 7. Test Rig in Room Temperature Test. F i g . 8. Test Rig Temperature C o n t r o l l e d 120 F.
MTI-2135 MTI-2136
% , t. -
SPEED POMER AND
VIBRATIONS PHASE ANRLE SHAFT MOTIONS
TORQUE FREOUENCY
FORCE MEAS. I TEMPERATURES CONTROL
1.875"
\-r-\
2000 „
, y ^ Shaft Dlaseter » 1.25 In.
"^ C Q . 3 mils
Mm _
800 - _
600 Unidirectional
» ^ Bearing
Stiffness, 1(
400 ^ tfe./ln.
.
200 „ -
100 ™
Radial Displacement » , .^^^^ I
80 Radial Clearance ' ' H ^ ^
60
B B
> / ^ ^ Oj^/^""'""^ I
40
" 10^
y*^ ^ ^ Potassium X. L
I.CbslO* 1.4 1.8 2.2 i.6 3.0
20 -
X totatlve Speed, r^~
10 __l i 1 1 i I 1
4000 MOO 12,000 16,000 20,000 24,M0 28,000
Rotative "Sneed, rpm Fig. 12. Bearing Stiffness Vs. Speed and Side Load.
MTI-2139 MTI-2140 •
200 1
T
Plvoted-Pad Bearing
Shaft Diameter^ 1.25 in.
Lubricant^ Water 120®P
Radial Clearance^ 1.5 mils
Preload Coefficient, 0.4
150 Unbalance5 Upper Bearing 0.50 gr.
Lower Bearing 0,50 gr.
Side Load Per Bearing, 25 lbs.
Maximum
Force, 100
lb.
50
0 10 15 20
R o t a t i v e Speedy rpm x 10"3
"
1.?
1.0 • ""^^^^==rf========,_^^ ;
I.';
irtJfIf^rt^'•
0.6 -
Shdf. DiflBteter. 1.2'j in.
0.4 Ratlial Clearance, 1.5 mil
Unbalance* 0.2S qr. i n .
Side load, 25 lbs.
0.2 .
MTI-2142
Radial
oisplaces»ntt
nils Squdlbrlijfl
OHplacesent
£ 10 14 }6 18 20
Rotative Speed, <
MTI-2143
OBJECTIVE
SPECIFICATIONS
SEAL FLUID POTASSIUM AT 600°F
SHAFT DIAMETER 1"
SHAFT SPEED 24,000 RPM
EXTERNAL: LOW VACUUM (lO'S TORR)
LEAKAGE: 1 TO 10 LBS. IN 10.000 HOURS
Vapor
Pressure
(nmHg)
Rotating Dlse
Squeeze Seal
MTI-2145 MTI-2146
ON
0^
i'iacu;.;^^
^.-.iif."
r^OLK^ijifiiS P'jMr-
MTI-2147
r -
t
"N
X
Vsxr^^
0^
U3
MTI-2148
0^
00
t. .
MTI-2149
6..2.39
Sealing
Differential
Pressure,
psi
1 i i — 1 — T " — 1 —
100 ^ J
^1a. » 6 In.
f
80 »- J
/
Conswned
60 ™
/
/
"4
40 —
/
J
20
0
_.
^:+l-
y _ j _
^
_^
Ma. = 4 In.
to
4>
O
(OsfOOi r••?i,^-^
s 01)-^
is SH ^ C/U!
i- i. ' i *. 5-. •.
MII-2153
ON
4>
m
SI
C
s
TIC^ IO%Cb
s
C A ^ B0LOV B07
3
a- KENNAMlTAL 7170
! •
..mmmi marff
CYLINDRICAL ^M^^i^^
© SPRiNe
^ STRAmMT
f^ TAPemo
<^z^zizz:^\
,^ CONICAL I CONTACT ^rm.
^oo mi
'zzzinz:^
•I
K'-^'Hi'.
r-iVOT PIN
i>iVCT PiiiTE
iOUSNAl
BfA?2iHG PAD
'•,
MTI-2156
r
4>
MTI-2157
MOTOR SUPPORT CASING
TOP CASING
MID CASING
BEARING HOUSING
SUPPORT CASING
ROnOM CASING
VACUUM ATTACHMINT
CASING
S. Frank Murray
ABSTRACT
J ^
The obvious problems of selecting materials for high temperature corrosion
resistance and strength have masked the fact that many of these liquid-metal
lubricated bearing and seal surfaces may fail by seizure long before corrosion
becomes a significant factor. Some of the hard, refractory materials, such as
the cemented carbides, appear to offer a means for achieving good sliding be-
havior, corrosion resistance and high temperature strength all in one package.
Actually, these materials are pacifiers. Design and fabrication problems gen-
erally outweigh the advantages that the cermets offer. In addition, the use
of these materials may limit the size of machine components."— ^
There are materials available with good sliding behavior and corrosion resist-
ance in NaK, but their physical characteristics impose many limitations on
design and fabrication. The purpose of this discussion is to review briefly
the status of the work on bearing materials selection and to recommend an
approach to the problem of achieving good sliding characteristics without
penalizing the bearing or seal design.
Since that time, there have been a number of other investigations made to select
bearing material combinations for various liquid metal environments, e.g., Refs.
4-9. In general, the conclusions have been much the same. The cemented car-
bides and tungsten or molybdenum are the most compatible sliding combinations
6
for liquid metal-lubricated bearings. Stellite Star J has also been added to
the list of promising materials. Certainly, these studies were useful since
they broadened our knowledge of the behavior of these materials in various
liquid metals. However, from the standpoint of the design engineer, who is
concerned with building practical machinery, this approach of trying to use
solid carbides or refractory metals appears to be the wrong way to attack this
problem.
For these reasons, MTI has leaned heavily on the use of coatings to provide good
bearing surfaces for unconventional fluid-lubricated systems.
EXPERIENCE ON COATINGS FOR BEARING APPLICATIONS
Most of the work which has been done at MTI on using coatings for bearings has
been directed toward gas-lubricated bearing systems. However, this work does
apply to the problem of liquid metal lubrication. In most of the gas bearing
work, the gases were inert and no contaminating oxide films were present to
provide surface protection. In addition, these non-reactive gases provided no
boundary lubrication. These two conditions apply equally well to lubrication
by NaK or by inert gases, although NaK is certainly a more severe case.
Basically, our approach has been to select one base alloy for the shaft, bear-
ings, and any other critical structural members. This alloy is chosen for
strength, density, corrosion resistance and any other essential attributes.
Then, plasma or flame-sprayed oxide or cermet coatings are used on the journal
and bearing surfaces to provide the necessary resistance to sliding damage.
These coatings are applied, then ground back to a standard thickness of .003
inches. In this way, problems of matching thermal expansion are minimized.
Experience has also shown that these thin coatings are much more capable of
withstanding thermal or mechanical shock than are the solid cermets or ceramics
Even among these so-called wear resistant coatings, there are wide gradations
in effectiveness. Under NASA contract, MTI has evaluated some promising gas
bearing material combinations for use in an inert environment (Ref. 10). In
those tests, a single hydrodynamic tilting pad test bearing was dead weight
loaded down against a 1.5 inch diameter test shaft. A photograph of the par-
tially assembled test rig is shown as Figure 2. Capacitance probes, mounted
on each side of the pad, were used to monitor lift-off speed, running film
thickness and wear. The purpose of those tests was to measure changes in bear-
ing performance as a result of sliding contacts.
oxide coated bearing against a chrome oxide coated shaft. The contacting sur-
faces were lightly polished and the capacitance probes indicated less than
forty microinches of wear. The third bearing, a hardened M-50 tool steel (60Rc)
running against a cobalt-bonded tungsten carbide coated shaft was effective for
about 350 starts and stops. At this point, the capacitance probe signals in-
dicated that the bearing performance was beginning to degrade. The test was
continued for the full thousand starts and stops, but there the bearing was not
lifting effectively at the end of the test. Examination of the bearing surfaces
after the test showed that the tool steel had transferred and welded to the
carbide coated shaft and that this transferred material was scoring the pad
surface in one area.
These same bearing combinations were also subjected to a series of high speed
rubs by bringing the shaft up to a stable speed of 60,000 rpm and then impacting
the pad against the shaft under various shock loads. Time in contact was on the
order of .02 - .05 seconds. The appearance of the bearings after those high
speed rub tests is shown in Figure 4. The chrome oxide coated pad running
against the chrome oxide coated shaft was still in excellent condition after
ninety rubs. The capacitance probe signals indicated that the bearing was
operating on a larger film thickness at the end of the test than it was in the
beginning. The contacting surfaces were highly polished and smooth.
As was the case in the start-stop tests, the M-50 tool steel pad running against
the cobalt bonded tungsten carbide coated shaft was a very poor combination.
This test had to be stopped before the full sequence of shock loads could be run
because the pad was no longer lifting. The pad was found to be severely cold-
worked and the shaft was covered with transferred metal.
The tungsten carbide coated pad running against an Al 0 coated shaft showed a
gradual deterioration in running performance throughout the test. After ninety
rubs, the probe signals indicated that intermittent contact was taking place
between the pad and the shaft. Examination of the test specimens showed that a
smooth film of the binder metal was smeared on the Al 0„ coated shaft. The
carbide-coated pad was severely abraded, either by the transferred metal or by
carbide grains being torn out of the pad surface.
Similar results have also been obtained with nickel and cobalt bonded tungsten
carbide coated pads running against the same coatings on the shafts.
From the results of these tests the self-mated chrome oxide coating was selected
as being the most promising material for good sliding behavior. This coating has
since been used successfully in a number of gas and water-lubricated machinery
applications.
a) Chrome oxide
b) Chrome oxide, undercoated with molybdenum to improve the adherence of
the oxide.
c) Cobalt-bonded tungsten carbide.
The coatings were sprayed on one end of a hardened, cylindrical, 410 stainless
steel test specimen which was part of a dumbbell tensile test bar. The geometry
and dimensions of the specimens are shown in Figure 5. The coatings were then
ground back to a standard thickness of .003 inches. This specimen configuration
was used so that some measure of the bond strengths of the coatings could be
obtained before and after exposure to NaK. To determine the coating bond strength,
the coated face of the specimen was cemented to the other half of the tensile
test bar with an epoxy resin. The bar was then pulled on a tensile test machine.
Normally, the break will occur at the interface between the coating and the sub-
strate since this bond is rarely as strong as the epoxy joint. Figure 6 is a
6
photograph of some typical test specimens showing the type of coating failure
which should normally occur. By comparing the bond strength of control speci-
mens with the bond strength of specimens which have been exposed to NaK, a
measure of the resistance of the coatings to attack by the liquid metal can
thus be obtained. This bond strength test, coupled with visual observation of
the condition of the coating on the surface, is a reasonable measure of the
resistance of the coating to corrosive attack.
The first test was a 300 hour static test with the specimens immersed in NaK at
200F. The specimens were then removed, examined and subjected to the tensile
test. Table 3 gives the bond strength values for the control specimens (which
had not been exposed to NaK), and for the specimens which had been exposed to
200F NaK. Two discrepancies were noted. First, the bond strength of the chrome
oxide coating, which was undercoated with molybdenum, actually increased by a
significant amount after exposure to the 200F NaK. At this time, we have no
explanation for this behavior.
The second discrepancy was the fact that the joints broke at the epoxy when the
tungsten carbide coated specimens, which had been exposed to NaK, were tested.
At that time, it was felt that this was simply the result of a poor epoxy joint.
However, repeat tests gave essentially the same results and in no case were we
able to get good adhesion between the epoxy and any of the carbide coated speci-
mens after the carbide coatings had been exposed to NaK. It is possible that
the NaK had reacted with the carbide surface or the metal binder to form a film
which did not make a good bond with epoxy. This introduced some uncertainty
into the test results but time did not permit us to resolve this.
The second part of the screening study consisted of a series of fifty hour tests
in 600F NaK. After each fifty hour cycle, the specimens were removed, washed
and examined. The reason for removing the specimens every fifty hours was to
see if the formation of caustic, as the result of exposure to air, would cause
more corrosion than the continuous immersion in NaK.
After the first fifty hour test at 600F, the chrome oxide coatings were found to
be completely disintegrated. Not only had the coating separated from the surface,
the particles of chrome oxide had also separated from each other and had ended
up in the bottom of the basket as fine debris. The molybdenum undercoat was
still intact on the specimen which had been undercoated. A photograph of one of
the chrome oxide coated specimens after this test is shown in Figure 7.
The cobalt-bonded tungsten carbide specimens had a dark film on the surface, but
the coatings were still intact. These carbide coated specimens were then put
back for another fifty hours. After this second run, there was a faint suspicion
of pitting starting to appear on the surfaces. This could have been due to the
fact that the NaK had washed out some grinding debris which had been embedded in
the surface. In retrospect, this explanation seems reasonable since the pitting
never became any worse.
The tests were suspended temporarily at this point while a new set of specimens
was obtained. This time, 410 stainless specimens with a nickel-bonded carbide
coating (25% WC, 7% Ni + mixed W-Cr carbides) were prepared. This nickel-bonded
carbide is known to be more corrosion resistant than the cobalt-bonded material.
In addition, a tilting pad bearing with this nickel-bonded carbide coating on a
416 stainless substrate was also included. This pad bearing had been ground and
lapped to a 1-3 rms finish for some previous tests. It was included here be-
cause the highly polished surface would make it easier to see if the coating was
actually pitting.
Finally, a type 316 stainless flat with an Al 0 coating was also evaluated just
to see if NaK would attack the Al 0 . The Al 0 coating was as-sprayed.
The specimens described above, plus the two cobalt-bonded tungsten carbide speci-
mens which had already been subjected to the two fifty hour cycles in 600F NaK,
were then run through three more fifty hour cycles at 600F in NaK.
At the end, the specimens were washed and examined. The appearance of the sur-
faces is summarized in Table 4. Figure 8 is a photograph of the polished carbide
6
coated specimen and the Al 0„ coated flat after the test. The carbide coated
pad appeared to be in very good condition, perhaps a little duller than before
the test but certainly not pitted. The Al 0„ coated flat had discolored badly
and looked much rougher. However, the coating was still intact. This test
should be repeated with an Al 0„ coating which has been ground back to a smooth
finish. The result does not prove that Al 0„ is unsuitable, in fact, Al 0„ has
been reported to be resistant to NaK at much higher temperatures but it does
indicate that there maybe some questions as to the ability of the coating to
withstand corrosion by NaKo
The shiny cylindrical specimen, with the threaded end, which is also shown in
Figure 8, is a test sample of 6061 T-6 aluminum alloy that had been plated with
0002 - .0003" of electroless nickel. This aluminum specimen had gone through
the five 50 hour cycle tests at 600F with no apparent corrosion taking place.
Table 5 summarizes the results of the bond strength tests on the carbide-coated
specimens before and after exposure to 600F NaK. Once again, poor epoxy bonds
were obtained on the specimens which had been exposed to NaK, particulary with
the cobalt-bonded carbide. However, the results do indicate that the losses in
bond strength, if any, were not large.
Based on the results of these static corrosion tests, the nickel-bonded tungsten
(W-Cr) carbide coating was selected as being the best choice for the bearing and
journal surfaces in the MTI NaK loop. As far as sliding behavior was concerned,
tbis coating was known to be satisfactory for start-stop operation. Its ability
to withstand high speed rubs was not outstanding, but seemed to be adequate for
several lightly-loaded momentary contacts. The corrosion resistance of the car-
bide coating appeared to be suitable at temperatures up to at least 600F.
This conclusion of the suitability of the carbide coatings for use in 600F NaK
was also shown by Wallace (Ref, 6) who had evaluated cobalt-bonded tungsten car-
bide coatings for their sliding behavior in NaK at 600F, and had found that these
coatings had good wear resistance.
Figure 9 is a photograph of the shaft and pads from the MTI NaK test rig. The
journals and pad surfaces have the nickel-bonded carbide coating. Although the
number of operating hours on this facility is still low, these bearings have
been run in water, alcohol and NaK. As far as sliding behavior is concerned,
the coatings appear to be very satisfactory.
One point which should be emphasized here is the capability of these coatings
for protecting larger areas such as face seal and thrust bearing surfaces. There
is no limit to the diameter of a part as far as the application of coatings is
concerned.
There are other problems in selecting the best coating for a given application
and in establishing reasonable temperature limitations for the coatings but this
information can be developed in a straightforward manner.
6.3
CONCLUSIONS
Flame or plasma-sprayed hard coatings offer many significant advantages in sim-
plifying the design and fabrication of process fluid lubricated bearing and seal
components. This Is particularly true for large diameter parts where mismatches
in the thermal expansion coefficients of materials may otherwise pose insur-
mountable difficulties. The effectiveness of certain of these coatings in
protecting the surfaces against sliding damage in inert environments has been
demonstrated in a number of practical applications. All of the test results
which have been obtained to date show that these coatings can also be used
effectively for liquid metal-lubricated systems.
The major drawback to the use of coatings is a lack of adequate quality control.
Either we live with the solid cermets and refractory metals and accept the design
problems and limitations which these materials impose, particularly for larger
bearings and seals, or we learn to use coatings and develop the necessary quality
control procedures which can ensure a reliable product.
REFERENCES
1. Vail, D. B. "Compatibility of Materials in Liquid Metal". USAEC Report
KAPL-589. August 18, 1951.
2. Vail, D. B, "Compatibility of Materials in Liquid Metal; Second Report".
USAEC Report KAPL-1021. January 5, 1954.
CONCLUSIONS
Flame or plasma-sprayed hard coatings offer many significant advantages in sim-
plifying the design and fabrication of process fluid lubricated bearing and seal
components. This is particularly true for large diameter parts where mismatches
in the thermal expansion coefficients of materials may otherwise pose insur-
mountable difficulties. The effectiveness of certain of these coatings in
protecting the surfaces against sliding damage in inert environments has been
demonstrated in a number of practical applications. All of the test results
which have been obtained to date show that these coatings can also be used
effectively for liquid metal-lubricated systems.
The major drawback to the use of coatings is a lack of adequate quality control.
Either we live with the solid cermets and refractory metals and accept the design
problems and limitations which these materials impose, particularly for larger
bearings and seals, or we learn to use coatings and develop the necessary quality
control procedures which can ensure a reliable product.
REFERENCES
1. Vail, D, B. "Compatibility of Materials in Liquid Metal". USAEC Report
KAPL-589. August 18, 1951.
2. Vail, D. B. "Compatibility of Materials in Liquid Metal; Second Report".
USAEC Report KAPL-1021. January 5, 1954.
STRUCTURAL MATERIALS
347 Stainless 9.3 0.29 Good
410 Stainless 5.7 0.28 Good
Columbium - iZr 3.82 0.31 Good
Haynes 25 9.4 0.33 Good
TABLE 3
RESULTS OF BOND STRENGTH TESTS BEFORE AND AFTER EXK)SURE TO 200 F NAK FOR 300 HOURS
ALL TESTS IN DUPLICATE
Coating After One After Two After Three More 50 Hour Tests
50 Hour Test 50 Hour Tests
Cobalt-Bonded Tungsten Dark film, Faint suspicion Dark film, easily washed off.
Carbide on 410 easily washed of pitting Same as after first two 50
Stainless Steel off hour tests
Cobalt-Bonded Tungsten
3490 psi About 40%
Carbide on 410 Stainless
Mri-2120
ON
TUNGSTEN CARBIDE
Fig. 3 Tilting Pad Bearings After 1000 Starts and Stops in Argon Atmosphere at 4 psi
H^^Bks
r
TUNGSTEN CARBIDE
Fig. 4 Tilting Pad Bearings After High Speed Rubs at 60,000 rpm in Argon Atmosphere
MTl-1963
6.3.22
TITLE
^ ^ ^ 'r" 1. / -^ /
/ a ^ ^:^£y^£/7^X C^A^y^-^^d^Z
• ,S0€> O/A.
L
7Zy/j scf^/'/9££ ra
J £>^ r^i^^j
APPHOVCO //.Jo>vej
CHECKED
^//(•/i C j '"»
SCALE
sy/i ^^^3
P R A W N ^ ^ J^^^y^^.«>iiC^^ £^^^
MTI-2122
CZTJ-IIH
in
J009 5B 3BN oa aansodxa ZB^JV snaratoads asax jo aotiBJBaddv g -STJ
CO
•10 Huini'oo o^'TV
S u i a c o o DpTqai?3 (:tD-jvi) u a n s B u n x
popuoq "[3^^TN 'poddr."!
Kv^',^5U^^'%^V'-''^~;,''
UJ
Fig. 9 Shaft and Tilting Pad Bearings from MTI NaK Test
Showing Coatings on Journal and Pad Surfaces
M|[|||24
7.1
SECTION 7
by
R. Kasuba
TRW^Inc.
Cleveland, Ohio
SUMMARY
The CRU V's were designed following a combined experimental and analytical
approach. With the past and present general scarcity of rigorous, analytical
methods in deriving a successful mercury lubricated rotor-bearing system
operating beyond the laminar regime, the initially selected analytical tech-
niques were progressively extended or empirically adjusted as needed.
The CRU V bearing system exhibited the critical damping factors of about .25
to .5 for the considered geometric and operational parameters. As the conse-
quence the shaft orbits and bearing reactions were effectively limited for the
entire speed range and the predicted critical speed could not be experimentally
observed.
The three-pad bearing provided sufficiently large stability regions by select-
ing appropriate pressures and geometric parameters. For example, at the CRU V
operating speed of 36,000 rpm the lubricant supply pressures of 225 psia were
sufficient to prevent fractional whirl instabilities for all variations of
bearing clearances and unidirectional and rotating load combinations.
Finally, the results described here were not only dependent upon development
of improved analytical techniques, but also reflect advanced design and
dimensional control of critical system elements through fabrication and
assembly procedures.
7.1.3
INmODPCTIOM
This presentation s e t s as i t s primary task to discuss b r i e f l y some of the r o t o r -
bearing dynamical q u a n t i t i e s and responses of the Mercury Hankine Power Systems
(SNAP 2) turbo-altornator u n i t s which are being developed for the United States
Atomic Energy Comission by TW. "Hie above turbo-alternator u n i t i , nore popu-
l a r l y named as ttie Combined Rotating Units (CRU's) use a liquid mercury l u b r i -
cated bearing system. A degree of success in developing the CRU's may be
i l l u s t r a t e d by the accumulated operational times. For exanple, the COT ¥ type
r o t o r s have accuinjlated almost 22,000 hours. Individual CRU ? u n i t s show
operational times up to ii700 hours. An additional 20,000 hours were accrued
by Mercury Rankine turbo-Mchineiy u n i t s of e a r l i e r designs.
Design and development of the CRU's was not an ea^y task. The rotor-bearing
lubricated bearings. The usual budgetary and engineering schedules for develop-
the basic CRU rotor bearing system's design was pursued along the p a r a l l e l
preselected bearings.
7.1.4 fl|
adjusted as needed^
t h e e l e c t r i c a l power r e q u i r e m e n t s j t h u s l e a v i n g the r e s p o n s i b i l i t y f o r o p t i m i z a -
follows:
Load c a r r y i n g c a p a c i t y
Power l o s s
Bearing f i M c o e f f i c i e n t s t o p r o v i d e s a t i s f a c t o r y range of b e a r i n g -
r o t o r responses
Stability
Resistance to cavitation-erosion
C o m p a t i b i l i t y of r e q u i r e d m a t a r i a l s
Easy CRU i n t e g r a t i o n
7.1.6
above mentioned bearings, the three-pad journal bearing and the hydrostatic
The selected hybrid three-pad bearing shown in Figure ii has t h r e e supply pads
equally spaced around the circumference and three equally sized coitip en sating
orifices. The three-pad bearing can be also described as geometrically
i d e n t i c a l to the three-sector bearing with the exception of the a x i a l drain
grooves.
The main experimental portion of the previously mentioned combined CRU design
CRP V Unit
The CRU ¥ and the associated test setup at the system's test level are illustrated
a single location far outboard of the turbine bearing. Limited direct instru~
mentation and structural inaccessibility made the actual CRU ? somewhat less
suitable for the fundamental and convenient research of the bearings. For
7.1.7
implementation of the shaft orbit data taken at a single axial location, several
accelerometers were installed on the CRU ? housing which was in turn supported
cross-checking element of the calculated shaft orbits for U.e entire rotor
tion levels with the predicted counterpart values was quite satisfactory and
The basic instrumentation and display systems used in monitoring the rotor
oscilloscopes, and magnetic tape recording and playback systems for studies
of startups and shutdowns. The capacitance probe signals were displayed both
and dynamic phase angles were referenced to the scribe marks in the shaft at
The basic test rig shown in Fi^nire 3 has capabilities for operation at elevated
temperatures to 500®F« The test rotor which is dynamically similar to the CRU
either upward or downward direction by means of the hydrostatic f^as type load
pads and the desired initial mode of unbalance can be selected by inserting
The basic instrumentation was similar to that used in the actual CRU's; however,
capacitajxe probes xn sets of two located 90 degrees apart a t two axial loca-
tions. Techniques used in displaying the capacitance probe data a t the system's
level were also employed at t h i s t e s t i n g l e v e l with the sarrple data presented
in Figure 5.
7.1.9
d e f i n i t i o n of l o a d s , a r e a s of s t a b l e ( n o n - f r a c t i o n a l frequency w h i r l ) o p e r a t i o n
i n v o l v i n g mercury l u b r i c a t e d t h r e e - p a d b e a r i n g s o p e r a t i n g in a non-laminar
regime.
r o t o r s a r e s u b j e c t e d t o b o t h r o t a t i n g and u n i d i r e c t i o n a l l o a d s of v a r i o u s
primary l o a d i n g mode.
of t h e r o t a t i n g magnetic l o a d s were e x p e r i m e n t a l l y e s t a b l i s h e d w h i l e t h e
The approximate range of the CRU ? load magnitudes can be given ass
Rotating magnetic unbalance loads - 2^0 to 7.0 pounds
or o n - s i t e balancing.
The multi-plane and two-plane balancing methods were considered and used a t
and the number of balancing planes can be selected by examining the computed
The force transmitted into the bearing film or bearing reactions can be
expressed as:
inertia effects.
The above effective film coefficients were obtained following hagg and
Sankey from the steady-state solutions of the following basic s e t of equa-
tions of motion;
M X + C^ X + K^ X = F cos (wt + 0)
F and M are the applied unbalance force and the applicable rotor mass per
to the major and minor axes of the journal locus. Employing tie measured
and
Y = S^i 0 = 0Y for X = 0
"cos^X . __...2
+ m U;
h • -^X
^sin ^X
c.
X - -J^
^sin ^Y
c.
Y " -JTW
The required displacements and phase angles for solving the above r e l a t i o n s h i p s
For the case where the predominant loads are r o t a t i o n a l loads such as i s
the case in CRU's the measured o r b i t s are c i r c u l a r and thus the effective
K^ =• Ky and C^ = C^
which then arbitrarily could be referred to as the radial bearing film stiff-
C
7.1.13
graphically in Figures 6 and 7, which show the effects of supply pressures and
clearances which will be designated here as the standard and large (1.3 x
standard) clearances. The basic trends of the hot bearing test results (ijOO^F)
A. For the considered range of loads which represent the eccentricity ratios
below 0.5 the bearing film coefficients of a three-pad bearing can be taken
B. The bearing film stiffness curves suggest again two distinctive regions
the inception of the high speed region somewhere above 25,OOU rpm range
ture from 70° to ijOO^F. VJithin the same temperature range the
7.1.14
Stability Threshold
Two stability threshold patterns of bearing supply pressure versus shaft speed
with no load, and unidirectional load versus bearing flow at a given speed are
experimentally for the vertical and horizontal shaft operation and can be
considered as the limiting patterns for the system since they represent the
The instabilities observed vrith the three-pad bearings for several orifice
sizes did not exhibit an abrupt transition between stable and unstable opera-
tion. The whirl orbit of these bearings gradually increased at the onset of
orbit plateau when the bearing was operated viell into the unstable region.
gas or oil lubricated bearings the whirl orbits typically increase from a
stable small orbit quite rapidly to a whirl orbit approximating the bearing
clearance orbit. The observed fractional frequency whirl was at .[48 to .Ii9
analyzers and also by four notches in the Lissajous pattern or two notches
- \
7.1.15
In summary, a t the CRU V operating speed of 36,000 rpm the lubricant supply
pressures of 225 p&ia were sufficient to prevent fractional whirl i n s t a b i l -
i t i e s for a l l operational v a r i a t i o n s of bearing clearances and unidirectional
and r o t a t i n g load combinations.
7.1.16
i n s t a l l a t i o n of t h e c a p a c i t a n c e p r o b e s .
of d i s c u s s i o n w i l l be s e t by e s t a b l i s h i n g t h e t h e o r e t i c a l r e l a t i o n s h i p between
depending on b e a r i n g c l e a r a n c e s and o o e r a t i o n a l c o n d i t i o n s .
ranges:
speed plateaus
bearing c h a r a c t e r i s t i c s .
7.1.17
With orbit measure! ^nts taken at a single axial location far outboard of the
turbine bearing, thr; maximum orbit at that location for some load configura-
tions may not determine the actual critical soeed. The phase angles between
for the equivalent rigid maps type members - aisplacement relationships are
not directly transferable to the flexible tjne lumbers under various load
about 20^ or luwer for the entire CRU speed range, the small changes in
the component level test data which indicatestinefirst critical speed for
the dynamicu.iiy equivalent CRU V shaft in the predicted critical speed ranges
changes of orbits and phase angles between the aoplied shaft forces and
journal displacement vectors over the entire speed range from 0 to 145,000 rpm
or above.
The behavior of rotor whirl orbits in the entire operational CRU V speed
range is shown in Figure 11, This presentation contains two types of datat
Eccentric3 t i e ? ;
e^ = .0001 in-
ep « ^0002 in«
e« = .0001 in*
(me)^ and (me)p - balancing weights needed for a close counterbalance of the
at 36,000 rpm
The presented curves show the corresponding onaft ortoits and phase angles
between the total rotor rotating force vector and the maximum displacement
the effects of the vectorial load combination as shoi^n in Figure 12^ with
the Curve II representing the response when only the rotating magnetic load
was rotated by 90 degrees from the original positiono Curve III represents
the load condition as Curve I with the load (me)„ removed. This is the load
condition which may be identified with some of the later CRU V units. As it
may be expectedj, the change in localized loads at the flexible rotor section
The CRU V housing and i t s own suspension system within the t e s t r i g s form an
housing-rotor responses*
of poor correlation of the shaft extension orbits in the case of CRU unit ID
related to the ootained satisfactory correlations between the measured and pre-
butions.
At the beginning of the CRU V testi T nroprnm tie desired or required instru-
mentation was not available and some elements had to be improvised and designed,
such as the caDacitance orobe extension stems, probe cooling systems, seals,
etc. VJith time the CRU ¥ test instrums^ntation was progressively improved to
of the gap, long duration and thermal drifts, etc. Acceleroneters on the other
hand were not affected by most of tie listed factors and correlations between
The directly measured rotor orbits were comprised by tie actual shaft extension
orbit and the physical shaft runout at tie nrobe location. The shaft exten-
sions were ground to less tlan 100 micro-incles (TTH) referred simultaneously
7.1.21
to both journal surfaces. However, in most cases even such a small runout
orbits. Thus, for obtaining a true shaft extension orbit from wl ich the journal
ship between the recorded runout and directly measured orbit at the capacitance
probe locations.
environment.
units, the capacitance probfs were removed leaving the housing accelerometers
iP^^
CAPACITANCE
PROBE
-ALTERNATC«
BEARING
^
THRUST
BEARING
TURBINE
BEARING
^
r "
FIGURE 2 CRU V TEST SET-UP
N3
FIGURE 3 BEARING TEST RIG
mKIMIWmBESlWSM MASS mmimi mmm
EAOHSMioNirrmsr BYlNSmLMIOHOFSCRliS
mffiOPTlOML
. - 1. _ s \ • » II
;i.-.:.;.
I..
•^^^•^•?sms^
1^ ::M^^m^m:^m^:r^
'l^rfii^v
:\W.\ _»..?«apiv>
VI* ^ti* %
LUBRICAN'l -MERCURY
TEST CONDITIONS:
S U P P L Y P R E S S U R E Z50 AND 400 P S I
ROTATING L O A D S :
.OO^^S AND . 0 0 7 Z S I N - O Z / B R G
T E M P E R A T U R E : ROOM
u
SO
m BEARING F I L M DAxM i n N G
J ^ ^ - -
40
O o ^ ^ r^
Pi
z
< 30
^
Q
20
180
160 /
CO BEARING F I L M S T I l F N E S S
CO
W
Z 140 /
H 120
A[
en
az
(^ en
100
p = 4 0 0 PSI
v#
r ^
rTl^^^
80
60
J^zz^ C ^ '^
p = 250 PSI
40
20
10 20 50 40 50
SHAFT SPEED - KRPM
7.1.28
50
U 40
30
o
'z Q
« 5 20
0-
160
140
en
W
Z 120
h
H
m 7 100
o
d X
80
o
«
w
pq
40
20
10 20 30 40
SHAFT SPEED - KRPM
7.1,29
Three-Pod Bearing
No Rotating Load
Lubricant Supply ^ 100°F
300
Drain Pressure- 6 psio
1 .3 X Standard Bearing Clearance
20 30
SHAFT SPEED (RPM X 10 3)
FIGURE 9 STABILITY THRESHOLD WITH UNIDIRECTIONAL LOAD
4 5 6 7 8
BEARING FLOW (LB/MIN)
FIGURE 10 CRU V CRITICAL SPEEDS VERSUS BEARING FILM STIFFNESS CURVES
7
80 ~
70 ~
1 SECOND MODE
CRITICAL SPEED
/
2 /
DS 60 -
^^_^^,^,,*«*'«*'***°FIRST MODE
o / ^t^^"^^^"^ CRITICAL SPEED
V)
D
Z /
<
CO 60 /
U' X
O /
2
^
/
C 40
J
J ^
/
/ / VA ^ y OPERATING SPEED
ii
L«»«««(
a. """^7^ S6, 000 RPM
/
/
<::: 30 , /
/ /
en
I1 1
A 1 / BEARING F I L M S T I F F N E S S CURVES
/ L STANDARD C L E A R A N C E , p= 400 P S I
2. STANDARD C L E A R A N C E , p= 250 P S I
20
3 . L 3 X STANDARD C L E A R A N C E p = 400 P S I
/ i 4 . L 3 X STANDARD C L E A R A N C E , p = 250 P S I
10
!
10-^ 5 6 7 8 9 10' 10° 10'
BEARING F I L M S T I F F N E S S (K). L B / I N
FIGURE I I SHAFT ORBITS DURING TRANSIENT STATE
OPERATION (SPINUPS AND SPINDOWNS)
ENVELOPE OF MEASURED
SHAFT EXTENSION ORBITS
FOR VARIOUS CRU V UNITS
C A t C U ^
CALCULATED
JOURNAL «BITS FOR
TYPICAL CRU V LOADS
LZZZZH L _ ^ ^
\
——J—-
r
— '
L
H
/ •?•
© © (8
©©
400
0| 191°
1 //
-i
f
300 _ II
200
01 113° ^.--..^^
100 ^ ^ ^ ^ 1 2 8 °
0|,l 140° ^i '''°^^^xJL^
0,1, 163°
© © ®
LOCATIONS LOAD DISTRIBUTIONS
1 ALTERNATOR BEARING (me) I BALANCING PLANE EFFECTS
2 BALANCING PLANE (me)3 BALANCING PLANE EFFECTS
3 BALANCING PLANE Fm ROTATING MAGNETIC LOAD
4 TURBINE WHEEL (me)6 THERMALLY INDUCED UNBALANCE
5 TURBINE WHEEL (me)/ LOCALIZED ASSEMBLY OR
6 TURBINE BEARING THERMALLY INDUCED
7 SHAFT EXTENSION UNBALANCE
8 CAPACITANCE PROBE ei,e2, DISTRIBUTED UNBALANCES
LOCATION
«3
7.1.34
Z
0 12 400
zu
^Z
X
-
I
3D 1C
i
4D
CRU V UNITS
i
2E ID
O
Z
5C
.60
4700 HRS.
2500 HRS,
Z> =40
- y
^Z
go 998 HRS
> < 20 HRS 683 HRS.
760 HRS. ^
.20
eg r-^
u
3D
i
IC 4D
i 1
2E ID 5C
CRU V UNITS MEASURED
^ PREDICTED
7.2.1
by
ABSTRACT
The test rotor is a uniform, flexible shaft with heavy wheels mounted
at the ends and in the middle. The rotor is supported in two silicone fluid
lubricated, tilting-pad bearings. The rotor amplitude caused by an induced
unbalance has been measured over a speed range of 3000 to 24,000 RPM for three
different rotor configurations, obtained by removing one or both end wheels.
This speed range extends to or through the third ciritcal speed for each of
the rotor configurations. The results are compared with the theoretical values
and, in general, the agreement is found to be good.
#
INTRODUCTION
To date, the shaft diameters of high-speed, process-fluid lubricated machines
have generally been selected with the intent of maintaining high shaft rigidity,
so that all flexural critical speeds are well above the highest operating speed.
Application of this design approach, however, becomes increasingly difficult
as operating speeds are raised, in addition to introducing needlessly high parasitic
losses in the bearings and seals. Accordingly, it is anticipated that as the
level of experience with high speed, process-fluid lubricated turbomachinery
increases, there will be greater acceptance of the flexible rotor design approach.
Parasitic bearing power losses may be minimized by adopting the flexible rotor
approach to the design of turbomachinery. This is particularly the case with
units that operate at high speeds and are lubricated with low viscosity cycle
fluids, such as water or liquid metals, so that the bearings operate beyond the
transition from laminar flow. In the superlaminar flow regimes, but where the
4
Reynolds Number is less than 10 , the bearing and seal losses vary as the 3.5
power of the diameter. At higher values of Reynolds Number, the exponent is
still larger. Moreover, the use of flexible rotors permits greater machinery
compactness and it may also allow design simplification. For these reasons, it
has long since been generally adopted for high speed, oil lubricated machines
such as steam and gas turbines.
In units that operate at high speeds and where the steady-state load on the radial
bearings is small, reductions in bearing diameter are generally permissible from
the standpoint of load capacity. The effect of increased flexibility of the
rotor will, however, cause highly significant changes in the dynamic response of
the rotor-bearings system. Accordingly, an accurate and detailed investigation
of rotor response to unbalanced loads is necessary during the design phase, in
order to arrive at arrangements that are not subject to high vibration amplitude
within operating speed ranges. Such an analysis must couple the dynamic chara-
cteristics of the bearings and those of the rotor, as well as the pedestals
or other flexible supports.
7.
Analytical procedures for computing the dynamic properties of fluid film bearings
in turbulent as well as in laminar flow regimes are currently available, from
earlier studies of turbulent fluid films 1,2 . These dynamic properties of
fluid film bearings, expressed in the form of stiffness and damping coefficients,
are then coupled with the dynamic properties of the rotor and those of the
pedestals or other flexible supports, to analyze the dynamic response of the
overall rotor-bearings-supports system, using procedures which are also currently
available 3 .
Before proceeding with the review of the analysis and the test data, it is appro-
priate to emphasize here the need for investigating the response of high-speed,
rotor-bearing systems over the entire speed range of interest. It is not generall
sufficient to merely compute the critical speeds, considering the rotor and bear-
ing flexibilities and then adjust the rotor-bearings design so that these critical
speeds lie outside the operating speed range. While critical speed calculations
are ©Wiously valuable, they provide only a part of the information which the
desigaer needs. In particular, the following information is not provided by
critical speed calculations, but may be obtained from calculation of the unbalance
response of the rotor-bearings system, including the effects of bearing damping:
2. Depending on the rotor and bearing designs, there will be great diff-
erences in the amplitudes of vibration which are encountered at the
critical speeds. In most cases, bearing film damping is the dominant
controlling Influence on rotor vibration at the critical speeds. Bear-
ing damping can be so effective that, in many machines, it is very diffi-
cult to identify the critical speeds without special instrumentation.
On the other hand, if the bearings ate not effective in controlling the
rotor vibration either because of their location or because they have
insufficient damping, there may be very large amplitudes at the critical
speeds. Thus, unbalance response calculations can be used to determine
the importance of the critical speeds and to evaluate the adequacy of the
proposed rotor and bearing designs.
ANALYSIS *
At any given speed and in the absence of rotor mass unbalance, the axis of the
shaft is in a stationary equilibrium with a shape determined by the static de-
flection of the rotor and the eccentricity of the journals in the bearings.
However, when the rotor contains an unbalance, the centrifugal forces cause the
rotor to whirl and deform such that each point of the shaft axis whirls in a
closed orbit around the corresponding equilibrium position. The whirl orbits
are closely elliptical and their size, orientation and phase angle vary along
the length of the rotor. It is the purpose of the analysis to set up a method
for calculating these whirl orbits for an arbitrary rotor subjected to specified
mass unbalances.
<i)
The four spring coefficients K ,K ,K and K , and the four damping coefficients
BXX ,Bxy ,Byx and Byy are determined directly
J
from the rpartial derivatives of F^ and
F , evaluated at the journal center equilibrium position corresponding to the
particular rotor speed. They may be obtained by perturbing and solving the lubri-
cation equation (Reynolds Equation) [3,65889,101. Thus, the coefficients depend
on the bearing geometry, the lubricant viscosity, the static bearing reaction and
the rotor speed. They can be expressed In dimensionless form as functions of the
Although Eq. (1), strictly speaking, is valid only for small amplitudes, the
expressions prove to represent the actual fluid film force also for large
amplitudes with surprisingly good accuracy as long as the minimum film thickness
during the whirling motion does not become less than approximately 25 per cent
of the radial clearance 13 I . In an actual application the rotor amplitude is
rarely allowed to be that large and the chosen linearized representation of the
fluid film forces is, therefore, adequate for most practical purposes. This
seems also confirmed by the present experimental investigation as discussed later.
For analysis purposes the rotor itself can be represented by a series of stations
connected by shaft sections of uniform cross-section as indicated in Fig. 1.
The stations are chosen to coincide with any concentrated rotor masses such as
wheels, impellers, etc. and, furthermore, stations are intro<luceci at each location
where the shaft diameter changes, at each bearing location and at both rotor ends.
Thus, each station can be assigned a mass, a transverse mass n»ment of inertia, a
polar mass moment of inertia, an unbalance and a set of 8 bearing spring and damping
coefficients. For any particular station, several or all of these quantities may
be equal to zero depending on which station is considered. Specifically, it should
be noted that the mass of the shaft is treated separately and is not included at
the rotor stations. Then applying a force balance and a moment balance, the jump
in shear force and bending moment across a rotor station can be expressed in terms
of the rotor amplitude and slope at the station. The equations are derived in
Appendix A where they are given as Eqs. (9) to (12).
For each shaft section the length, the cross-sectional area, the cross-sectional
moment of Inertia and the shape factor for shear deformation are known. The shaft
material properties: Youngs modulus, shear modulus and mass density, are also
given. Then, the relationships for the amplitude, the slope, the bending moment
and the shear force between the two ends of the shaft section can be established
7.2.
The equations for the rotor stations together with the equations for the shaft
sections establish a set of recurrence formulas by which a step-by-step calcula-
tion of the rotor can be performed. Assuming the rotor ends to be free, a succ-
essive application of the recurrence formulas allows coaaputing the amplitude, the
slope, the bending moment and the shear force along the rotor in terms of the
amplitude and the slope at the first station. Adding the contribution from the
unbalance forces and equating the bending moment and the shear force at the last
station to sero, makes it possible to determine the unknown amplitude and slope at
the first station. By back-substitution, the amplitude at each station can finally
be calculated. The procedure is given in more detail in Appendix A and is readily
progranmed for a digital computer. Such a program has been used to obtain the re-
sults for comparison with the experimental data as discussed later.
EXPERIMENTAL APPABATUS
The test rotor is a cylindrical steel shaft with overall length of 41 inches and a
diameter of 2.5 inches except for a centrally located integral 4isk of 6 inches
diameter and 6 inches length. Detachable disks of 6 inches diameter, 3 inches
length can be mounted on the ends of the rotor by heavy interference fit and lock
nuts which clamp them against shoulders on the shaft. The total weight of the
rotor, exclusive of the detachable end masses, is 88 lbs. The center disk, not
counting the inner 2,5 inch diameter section, weights 36 lbs. and the end disks
each weigh 18 lbs.
The rotor is supported by two tilting-pad journal bearings whose centerlines are
12.5 inches on either side of the rotor center plane. These are four-pad bearings
with slenderness ratio (L/D) of 1.0 and bearing diameter of 2.5 inches. Other
design characteristics include:
-3
Clearance ratio (C/R) - 3 x 10 in/in (based on machined clearance)
Pad arc length - 80 degrees
Pivot position - 0.55 of the arc length (44 degrees) n^asured from the
leading edge.
Geometrical preload coefficient - 0.5
The bearings are rigidly mounted in housings machined from a single block of
aluminum jig plate. The housings are bolted and keyed to a massive structural f|
4
steel base whose moment of inertia about the horizontal axis exceeds 9,000 in .
There were no detectable traces of flexure or vibration of the base during the
experiments.
The rotor is driven by an electric motor through a crowned spline coupling. The
coupling can accommodate misalignment, between the motor and rotor axes, of up to
0.030 inch without restraint to the shaft. The motor is bolted and keyed to the
same base surface which supports the bearing housings. Motor and test rotor
shafts were aligned with less than 0.002 inch total runout measured radially and
0.0015 inch total runout measured axially on a 2 inch radius. A variable-frequency
motor-generator set supplying power to the motor provides for variable speed
operation within a range from about 3000 to 24,000 rpm.
The journal bearings are lubricated by a very low viscosity silicone fluid, 0.65
cs at 77F and 0,51 cs at 130F. Lubricant at controlled temperature is pumped into
the bearing housing so the bearings are operated submerged in the lubricant. There
are clearance seals and scavenging rings on either side of the bearing housings.
The motion of the rotor is measured in both horizontal and vertical directions at
each of five measurement planes along its length. There are measurements planes
at each end, at the center, and at the inboard sealing ring of each bearing housing
(8.65 in, from the rotor center plane). These measurements are obtained with eddy-
current, non-contacting proximity sensors. Readout is by oscilloscope with the
signals connected to the horizontal and vertical deflection plates of the cathode
ray tube. With this arrangement, the rotor center axis at the measurement plane is
represented by a spot on the oscilloscope screen and motions of the rotor axis are
reproduced by the spot on the screen.
Accurate measurement of the response of the rotor to imposed unbalance load requires
care in minimizing other causes of real or apparent shaft axis motion. The principal
such causes are shaft surface out-of-roundness or non-concentricity and residual
rotor unbalance. Rigid specifications were imposed on rotor surface roundness and
concentricity with the result that subsequent gaging showed all surfaces to be
7.2
round and concentric with the bearing journals within less than 0,00012 inch total
runout. The rotor was balanced first in a balancing machine and then in place in
the apparatus. In-place balancing was done by a trial and error procedure. The
single, centered, disk rotor was balanced so that the maximum orbit amplitude at
any measurement plane over the entire speed.range was 0.00014 inch radius. For
the two-disk system, which is more sensitive to unbalance, the maximum orbit
radius was 0.00021 inch at the light end and all other positions were always below
0.00015 inch radius. The three-disk system was balanced at speeds up to 10,800
rpm. At this speed, the maximum whirl orbit radius was 0.00035 inch. At all
speeds below 10,000 rpm, the maximiun orbit radius was below 0.00013 inch.
In order to minimize the effects of runouts and residual unbalance on the test
results, it is desirable to Impose large deliberate unbalances to obtain sizeable
response orbits. On the other hand, the assumption of linear bearing stiffness
and damping about the steady-state rotor position will be strained if the orbits
are too large. For most operating conditions, a good balance between these con-
siderations was considered to be: the larger of the orbit radii at the two
stations adjacent to the bearings should be between 0.20 and 0.40 times the bearing
radial clearance measured at the pivots, or, between about 0.0004 and 0.0008 Inch
orbit radius. This limit was not applied to the three-disk rotor for reasons which
will become clear. During the experiments, the amount of unbalance weight was
adjusted as the speed was changed In order to remain within the limits given above.
Since the orbit amplitudes at the center and end positions were usually larger than
those at the bearing positions, the orbits with Imposed unbalances were usually at
least 4 times larger than the corresponding orbits with no Imposed unbalance.
A computer program based on the preceeding analysis was used to calculate the
unbalance response of each of the three rotor configurations, that Is, with one,
two or three disks. Dynamic stiffness and damping values for the test bearings
were obtained from the data given in Ref. 5- The stiffness and damping vary with
the shaft speed and steady-state load and, since the test bearings operate in the
turbulent flow regime, with the mean film Reynolds Number. In the case of the
four-pad bearing with steady-state load line midway between the pivots, synmetry
results in elimination of the cross-coupling terms and equality of stiffness and
damping in x and y directions. Sample dimensional stiffness and damping values
are given in Table 1 to illustrate the values obtained.
The computer program output included the calculated vibration amplitude at each
of the planes where experimental measurements were made. Because of test bearing
symmetry, the calculated vibration is a circular orbit whose radius is given In
mlls/oz-in. total unbalance.
Flexibility of the pad pivots was Introduced into the calculations as a stiffness
in series with the bearing. The calculated Hertzian contact stiffness of the
combined pivots at 100 lb. bearing loa<
load is 5 x 10 lb/in. and this is the value
used in the calculations at all speeds,
The undamped critical speeds of the three rotor-bearing systems were calculated
using the Prohl method (Ref. 11) together with the theoretical bearing stiffness
data. This was done for comparison with the calculated and measured unbalance
response data. The comparisons are significant since it is common practice to
use calculated critical speeds in rotor-bearing design analysisjbut much less
common to perform unbalance response calculations. The results are tabulated
below for each of the three rotor configurations.
Stiffness, lb/in.
(1) One-Disk Two-Disk Rotor Three-Disk
Speed,rpm Rotor Light-End Heavy-End Rotor
Damping, lb.sec/in.
The theoretical curves show the one-disk rotor system to be well behaved through
the operating speed range. The orbit amplitudes level off about midway between
the first and second calculated critical speeds and remain about the same or
gradually decline through the remainder of the operating range. The experimental
data show the same features and the absolute levels of vibration amplitude agree
rather well with the calculated levels. While there is some experimental scatter,
the measurements show no marked differences between the coupling and free ends
or between the two angular positions of the unbalance weights.
According to theory, the response orbits should always be circular. The actual
measured orbits were not always so, although their ellipticity was seldom marked.
There were no recorded instances when the ratio of major to minor axes was greater
than 3 and very few when it exceeded 2. When the orbits were elliptical, the major
and minor axes were measured and the mean between them was taken as the experimental
data point. There are a number of possible reasons for non-circular orbits of which
the most significant is believed to be asymmetry of the bearings, because of
slightly different preload between the two opposing pairs of pads.
Calculated curves and experimental data for the two-disk rotor with the unbalance
weights in-line are given in Figs. 5 through 9. This configuration experiences
7.2.15
Results for the two-disk rotor when the unbalance weights at the two disks are 180
degrees out-of-line are given in Figs. 10 through 14. This results in a sharply
defined conical mode vibration peak at about 11,000 rpm with noticeably larger
amplitude than occurred when the unbalance weights were in-line. Again, agree-
ment between experimental and theoretical results is good with respect to both
the trend and the whirl orbit amplitudes.
At the third, free-free mode, critical speed of the three-disks rotor system, the
bearings are very close to the rotor nodal points. The result is a sharp peak
with very large whirl orbit amplitudes at the calculated third critical speed.
Whirl orbit amplitudes for the out-of-line arrangement of unbalance weights are
shown in Figs. 15 through 17.
The total unbalance weight used when passing through the response peak was about
0.18 in-oz which results in very small orbits and comparatively poor agreement at
higher speeds when the response amplitudes are much lower. When passing through
the third critical speed with this small unbalance, the or-bits at the stations
adjacent to the bearings were appreciably larger than the bearing clearance.
Because of concern over the well-being of the apparatus and its operators, the
rotor was accelerated through the response peak and there was no steady-state
operation at that point. Peak amplitude data were obtained by using a camera
attachment for the oscilloscope with the shutter held open during the passage
through the peak. The largest recorded orbit amplitude is shown as the data
point at the peak in Figs. 15 through 17,
The modes of rotor vibration during the experiments were determined from obser-
vations of the phase angles between the direction of the unbalance force and the
response at the different measurement planes. The measurement was accomplished
by using a sensor which responds with a sharp voltage peak as a mark on the shaft
passes. The amplified voltage peak is used to momentarily dim the beam of the
oscilloscope so that there is a short interruption of the orbit on the screen
which identifies the location of the rotor axis on its orbit at the time the mark
on the shaft passes the sensor. The phase angle is measured on the oscilloscope
as the angle between the break in the orbit and the known angular position of the
unbalance weights at the time the mark passes the sensor. The accuracy of the
measurement is adversely affected if the orbit is not truly circular. For this
reason and because calculated phase angles are of comparatively little significance
in design analysis, the emphasis has been placed on the response amplitude results.
However, a typical sample of the agreement between measured and calculated phase
angles is given in Table 2 for the two-disk rotor with the unbalance weights
in line.
I
TABLE 2
Calculated
Speed, rpm Light End Light End Center Heavy End Heavy End
Bearing Bearing
3000 30 32 32 34 195
6000 80 89 90 90 220
9000 124 136 146 170 280
10,000 131 146 158 195 288
11,000 137 154 163 260 318
12,000 143 158 163 340 340
15,000 156 168 158 126 33
18,000 163 174 164 134 60
Measured
l»
DISCUSSION
Comparison of the measured response of the test rotor-bearing system with the re-
sults calculated from the combination of the rotor response analysis described in
this paper and theoretical bearing stiffness and damping has shown good agreement.
In particular, there is agreement in the following respects
(1) For each of three rotor configurations, the trend of the response to
unbalance as the speed was varied over a range extending to or through
the third undamped critical speed was substantially as calculated,
(2) With few exceptions, the measured and calculated amplitudes of vibration
are comparable. Examination of either measured or calculated results for
the purpose of determining the relative sensitivity to unbalance of the
three different rotor configurations clearly leads to the same conclusions.
critical speeds. Results for the two-disk system with unbalance weights in-line
are especially Interesting since they show two well defined peaks whose mode
shapes correspond to those of the first two critical speeds. These vibration
peaks occur at 1.7 and 1.4 times the first and second critical speeds respectively.
These differences between critical speed and rotor response results are caused
by the effects of bearing damping which is not included in the critical speed
analysis.
The experimental portions of this paper are taken from research supported by
the Atomic Energy Commission under Contract Number AT(30-1)-3363. The
authors wish to thank Mr. N. Grossman of the AEC and Mr. E.B. Arwas of M.T.I,
for their support and guidance in carrying out the experimental work.
7.2.21
NOMENCLATURE
I
Nomenclature (cont'd) •I
V Shear force, lbs.
'^ Phase angle for whirl orbit, see E q s . (19) and (20)
Indices
X x-component
\A y-component
•I
7.2.23 \
I
Nomenclature (cont'd)
Superscripts
••
REFERENCES
2. H.G. Elrod, Jr., C.W.Ng and C.H.T. Pan, "A Theory of Turbulent Films and Its
Application to Bearings," Topical Report under AEC Contract AT(30-I)-3363,
AEC Report No. NYO-3363-2, MTI Report 65TR9, March 1965 (Also published
under the same title as ASME Paper 66Lubl2, June 1966).
3. J.W. Lund et al: "Rotor Bearing Dynamics Technology," Final Reports under
USAF Contract No. AF 33(615)-1895, U.S. Air Force Reports AFAPL-TR-65-45
Parts III and V, MTI Report 65TR14 and 65TR15, May 1965.
5. (a) F.K. Orcutt, "The Steady State and Dynamic Properties of the Tilting Pad
Bearing in the Laminar and Turbulent Flow Regimes," Topical Report under
NASA Contract NAS-w-1021, MTI Report 65TR32, June 1965 (condensed version
published under the same title, as ASME Paper 66 Lub 19, June 1966).
(b) F.K. Orcutt and C.W. Ng, "Steady State and Dynamic Properties of the
Floating Ring Bearing in the Laminar and Turbulent Flow Regimes,"
Topical Report under NASA Contract NAS-w-1021, MTI Report 65TR33,
June 1965,
6. B. Sternlicht, "Elastic and Damping Properties of Cylindrical Journal
Bearings," JOURNAL OF BASIC ENGINEERING, TRANS. ASME, Series D, Vol. 81,
1959, pp. 101-108.
7. A.C. Hagg and G.O. Sankey, "Some Dynamic Properties of Oil-Film Journal
Bearings with Reference to the Unbalance Vibration of Rotors," JOURNAL
OF APPLIED MECHANICS, Vol. 23, TRANS. ASME. Vol. 78, 1956, pp. 302-306.
8. A.C. Hagg and G.O. Sankey, "Elastic and Damping Properties of Oil-Film
Journal Bearings for Application to Unbalance Vibration Calculations,"
JOURNAL OF APPLIED MECHANICS, Vol. 25, TRANS. ASME. Vol, 80, 1958, p. 141.
10. P.C. Warner and R.J. Thoman, "The Effect of the 150-Degree Partial Bearing M
on Rotor-Unbalance Vibration," JOURNAL OF BASIC ENGINEERING, TRANS. ASME, ^
Series, D, Vol. 86, 1964, pp. 337-347.
11. J.W. Lund and B. Sternlicht, "Rotor-Bearing Dynamics with Emphasis on
Attenuation," JOURNAL OF BASIC ENGINEERING, TRANS. ASME, Vol. 84, Series D,
1962, pp. 491-502.
12. J.W. Lund, "Spring and Damping Coefficients for the Tilting Pad Journal
Bearing," ASLE Transactions, Vol. 7, 1964, pp. 342-352.
14. M.A. Prohl, "A General Method for Calculating Critical Speeds of Flexible
Rotors," JOURNAL OF APPLIED MECHANICS, Vol. 67, 1945, pp. A-142 to A-148.
7.2.26
APPENDIX A
This appendix describes the analysis for calculating the whirl amplitudes of an
unbalanced rotor. The analysis takes into account the anisotropic stiffness
and damping characteristics of the bearings which cause the rotor whirl orbit to
be elliptical. Furthermore, the gyroscopic moments of the rotor wheels are in-
cluded. The analysis is an extension of the Prohl Method 14 and is readily
programmed for a digital computer.
To define the rotor motion, introduce a cartesian coordinate system with the
X-axis vertical downwards, the y-axis horizontal and the z-axis along the rotor
axis. The origin of the x-y-system coincides at each axial location with the
static deflection of the rotor, corrected at any given speed for the eccentricity
of the journals in the bearings. Thus, the rotor amplitudes are x and y, the
corresponding slope components are 6 and w , and the bending moment and the
shear force in the rotor are denoted as M and V, respectively. They have the
components M and M , and V and V .
X y X y
(3)
Considering the n'th rotor station as shown in Fig. 1, a force equation and a
moment equation yields:
IT» 5^+<"!,. f ' Mj,-M,„ (6)
The rotor Is forced to vibrate by the unbalance forces with the frequency CO.
The vibration Is harmonic whereby the amplitude, the slope, the bending moment,
the shear force and the unbalance can be expressed In complex notation:
where: vC-U|+t U ^ ^ X sXc4t>f| j^'tc. ^ and only the real parts apply.
Zxyf-Kx'^i'^^)'^ (8)
and analogously for Z)fjj ^Zfejf and Zjj^ , With these definitions, and dropping the
bar notation, Eqs. (4) to (7) become:
.28
The n'th station is connected with the (n+1)'th station by a shaft section of
length- C|^ with a cross-sectional area An , a cross-sectional moment of Inertia
Ih and a cross-section shape factor rt^ for shear deformation. The shaft
material has a mass density ^^ , Youngs modulus C.^, and shear modulus Q^ « Then,
as shown in Appendix B, it is possible to express ^niii^^+i j H ^ i|.^Eid Vxin-^l
as linear combinations of X^fO^^ ^ tiHn ^^^ ^ (the relationships for the
y-components are identical). The equations are given as Eqs. (40) to (43) In
Appendix B. The coefficients in the equations are functions of the shaft pro-
perties and the speed of the rotor. For simplicity, It is assumed that unbalance^
rotary Inertia and gyroscopic moments in the shaft can be Ignored, and that
these effects instead are included at the rotor stations.
Equations (9) to (12) together with Eqs. (40) to (43) are a set of recurrence re-
lationships from which the rotor amplitudes can be computed. Assuming the rotor
to have free ends, the bending moments and shear forces at the ends are zero:
(13)
1;I-W:H: (14)
M)r|'
n
= {^H!:+ =0 (15)
^nJ
where {<^ls a 2 by 4 matrix and d j l s a 4 by 4 matrix. The matrix elements C^l,
and dy are complex. They are found simplest by performing a total of 5 rotor
calculations. In the first calculation, set Vj-sj and Mj - @i = 0^} - Uj, - (?
whereby C^n X^^i ^du^dti ^ d^j and cf^j are obtained. Next, set ^j-' and
y^"^&ir ^j~ U.^~Q and determine the second columns of the matrices. Proceeding
In this way, all the coefficients are obtained. Solving Eq. (15) for ^ij^ij®i
and (Ij allows computing Xj, and u^ from Eq. (14) for all the rotor stations.
Noting that:
XH- CHK.-^iK
"Sw
and similarly for U^ the semiaxes and the orientation of the ellptlcal whirl
orbit are calculated from
*/g
(17)
(18)
(19)
30
where a^ Is the major semlaxls, b^ the minor semlaxls, If^ the angle from the
x-axls to the major semlaxls In the direction of rotor rotation, and '^ is the
phase angle. The definition of the phase angle is such that If the x-y-coordinate
system Is rotated the angle >!, into an x'-y'-system (i.e. x' is along the major
semlaxls), then the rotor motion can be expressed as:
(20)
If the value for the minor semlaxls Is negative, the rotor precesses backwards.
For those bearing types where £.^y^~ c.u^ and ^y«"'"^«fy It Is seen from
Eqs. (9) to (12) that:
Hence, the whirl orbit becomes circular and the outlined calculation procedure
can be considerably simplified. This condition applies to the four shoe tilting
pad bearing used In the present tests.
7.2.31
APPENDIX B
In this appendix a uniform shaft section Is considered and the relationships are
established for the amplitude, slope, bending moment and shear force at one end
(22)
® " o ^ ^ "^ar
and similarly for the y-component, where Q Is the modulus of shear and «4 Is
a cross-sectional shape factor ( d ^ 0.75 for a circular cross section).
The bending equation Is given by:
H.^ElIf (23)
jA | p = - 1^ +WW (24)
and analogously for the y-dlrectlon. Finally, if per unit length the shaft has
a transverse mass moment of Inertia |-- and a polar mass moment of Inertia la
a moment equation results In:
= c/fu,+ . » ^ i k + i l k i i ^ ^ l l i k ] (27)
and similarly for the slope, the bending moment and the shear force. Only the
real part applies. Furthermore, setting:
With these assumptions, Eqs. (29) and (30) become Identical and It is only
necessary to consider the first equation. The following definitions are Introduced:
^ = (If^f
/= gl
(31)
(32)
Dropping the bar notation and setting *A-J.J-* jp-0^ Eq(29) becomes:
p* EI dz
Substitution from Eq. (36) allows determining "^x^W ^^-^^ ^ <• Next, the
four constants can be evaluated by setting:
(41)
2 I s
(42)
(43)
where:
(51)
L- 2 - ^ (52)
^^^--m
. ^
''''
X,-^j4 K^f>tL (54)
z •A l (55)
It should be noted that AJU / "tn and n^*, are the static Influence coefficients
neglecting shear deformation. Furthermore, for sufficiently low speeds and suff-
iciently many mass stations, the shear effect can be Ignored and pS'O, Then,
the coefficients ^j„ to ^^„ become equal to 1 and the shaft mass can be lumped
with good approximation at the ends of the shaft section, setting Mt^-O in the
above equations. In most practical cases, this will be accurate enough.
UNBALANCE
BEARING n n^ I BEARING
• L . ^
2 ""3 ^"4
•^ — z»zijs-^q
" ^
IfttAj^.an.EntGj,,^
^0, .^^ntl
I.
I
STATION 0 STATION ( n ^ l )
Sign Conftntion for Amplitydei Slope, Binding Moment end Shear Force
0.5 . '0.5
A
A
A
/
n
1•
K
o
'"'"**'**%_
iL
i
i
A
1
u - ^
A
/
^x. ***»«»»». _—...»«« """"^
^ . j t f ^ ^
/A o
f A
/ O
m
K • 1
V
I I •
/
/
/
/
r1
i r » I F™ ,f\ » . ™ , < - ~ , - i - , , p-\ff~rf^s*fi •yrf^
mtUKtflLAL Kt5ULI5
EXPERIi ENTAL DATA
o ® COUPLING END
AA FREE END
0 i\ £\ YiI 16 2D 24 28
SPEED, RPM X 10^
MTI-2171
0.5 0.5
1 A
A
r-.„.
1 1
"T~~
BEARING POSITIONS
<
$
1 \
©
^
>
4
( L„
/ A W^^ -^^^
" / • ^t^
'x.
/X
^ X•
n
/
r r —
« J f © •
THEORETICAL RESULTS
EXPERI MENTAL DATA
tUUP Liib tNU titm m
A A FREE END 1BEAR
m
0 4 8 12 16 20 24 28
SPEED, RPi X 10^
Fig. 3. Unbalance Response of the One-Disk Rotor, Bearing Positions
7.2.39
m
»
/ o N o
»
J
/
f o
m
ps ^ "
m
.^^^^
I »
o c; •
m
/
^ i.O i ®
/
1
/< )
/
O
I
7
CO
CO
o i
J
m
sn
go. /
/
/
®
CENTERLINE
n i
fi. 1i
.05
*- • '^
J
CO i — ^l ^ ^
oa.
CO
—JH A ; ' 1
Ti ir-rsDi-Tir' h1 n r - o i ii •r n ^ \
0.01
8 12 16 20 24 28
SPEED, RPM X 10^
/ m
5 "*7" ^ - ^ n^
/ o — — •
o / °
c
? LO / >
r.
J
O J
I
CO
1 I€AVY EN[)
0.!5 i 0.5
/ 1
-~-w~~] k
Q /
i rlEORETICAL RESULTS
m / o . E^ PERIiEiTAL DATA
0.
. 11
o
Q.
€0
0.01
0 I\ i1 12 H5 20 24 28
SPEED, RPi X 10^
Fig. 5. Unbalance Response of the Two-Disk Rotor, Weights In-Line, Heavy End Positions
7.2.41
®
m
>
/ " o o o
1.0 ® »
o D
/ ®
m
1
M " /
O
• /
€/5 /
® i
r
CO
Q
OS
0.01
8 12 16 20 24 28
Fig. 6. Unbalance Response of the Two-Disk Rotor, Weights In-Line, Heavy End Bearing
Position
MII-2175
7.2.42
.0
X- \ \ <>
__JO__ i) \
m/
X• • \
o o \ ^^^-^ ^ -
/ / ^ •
/° \ / ^ 9
11
- — _
O I.
r
_J A 1•
m • /
KJ
o
I
CO
• / f
CO 0. rFMTFRI IMF
Q
1o 0.5 1 0.5 1
a:
m
tt: ft n
o
ijj
CO XUrAorTII^AI D C C I I I TO
o i ncuni-1 iwHL ncouLij
Q_ o . EXPERIiENTAL DATA
CO
0.01
0 8 16 20 24 28
SPEED, RPi X 10^
1 1 ®
® (»
o
m 3
o o
J< ' \ L
f
O
^
o ^
»».»_
V
/ \
A
1
>
o ®l
•/
/
/
r o
1
o
1 J
LIGHT END BEARING
/ 1 0.5 i 0.51
/
/ I — "
_
k
_
1 1 '• L_J
/
/ THEORETICAL RESULTS
I o ® EXPERIMENTAL DATA
0 I\ £\ 12 16 20 24 28
SPEED,RPi X 10'
Fig. 8. Unbalance Response of the Two-Disk Rotor, Weights In-Line, Light End Bearing
Position
MTI-2177
l.l.kk
<•
A
. \
#
l \ i _ —
o
o
/o V
•
\
v^
7
o
c
D
/ 1 iriHT Fwn 1
1
" / 0.5 I 0.5 i
/
. /
ft A
o
/
• / THEORETICAL RESULTS
r o • EXPERIiENTAL DATA
8 12 16 20 24 28
SPEED, RPi X lO'
Fig. 9. Unbalance Response of the Two-Disk Rotor, Weights-In-Line, Light End Position
MII-2178
7.2.45
jf%. O.b
m '>
-^r^^l
/
r
r< »
m
\
\
V
0.5
®/
>
1itAV r tm
4
m
s;7
i /
/
/
^
1
1
if
o J
®/
THEORETICAL RESULTS
o» EXPERIiENTAL DATA
1 1 1 1 ^ J
0 8 12 16 20 24 28
SPEED, R P i X 10'
10. Unbalance Response of the Two-Disk Rotor, Weights Out-of-Line, Heavy End
Position
MTI-2179
l.l.hb
©
t
~v o
•
©
<>
V •
!y f
^ * * 1 t _
^^ •^-Xs,^^ 3
7 ^ - N
s^
m ^ ^ x^
/
© -^
• ''^'*"^^ w. b
/f « I
- ^
/
• —-H
<
(/
o
05
•/
OJ
• / 1
/ 11L=H¥ 1 l»ilL/ UL.mi^ii'^v 1
/
/
/
/
'
— THEORETICAL RESULTS
o • EXPERIiENTAL DATA
1 1^1 1 1 J
8 12 16 20 24 28
SPEED, RPi X 10^
Fig. 11. Unbalance Response of the Two-Disk Rotor, Weights Out-of-Line, Heavy End
Bearing Position
MTI-2180
7.2.47
„1
@
f'\
1 I ®
O
/ \"• ^ «1
1.0 • 1
m
/ V
f V.^
"J 1
r*®»»-.«
r— ^
/
/
o
I /
0
CO J/ CENTFRLINE.
o
©
7 05
•
0.5
® /
0.1
CO
i
CO
ixJ
oc
THEORETICAL RESULTS
o• EXPERIiENTAL DATA
1 I 1 1 1
0.0I|
3 ^i i1 i;2 16 20 24 28
SPEED, RPi X 10'
Fig. 12. Unbalance Response of the Two-Disk Rotor, Weights Out-of-Line, Center Positi on
Mri-2181
7.2.48
, ri »
•© • • 1
•
\
yj /
o 1.0 /o
f
\ .
^ r •
s^\ j t »
J >
ffi
z>
7 __o P—»^
---^^
/
O t
/
1
f
ft 1
CO
•/ 0.5
•
CO J^X^^^^Sl
?BIT RADIU
O i T
# /
' 0.5
p
— THEORETICAL RESULTS
o» EXPERIiENTAL DATA
1 I 1 ...,..±™^ L^ L
8 12 16 20 24 28
SPEED, RPi X 10^
Fig, 13. Unbalance Response of the Two-Disk Rotor, Weights Out-of-Line, Light End
Bearing Position
MTI-2182
1.2.h9
m
>
PV
k
o
» m
e o
<^1)
1
<>
o
\
®
X ^
1 ^ S^ 1 >**^
J
/
/
<
/
f
J
1
%
•
0.5
n
1, n
m
0.5
/
/ L GHT END
/
/
\ m
'
THEORETICAL RESULTS
o ® EXPERIiENTAL DATA
1 1 1 1 1
8 12 16 -20 24 28
SPEED, RPi 110^
Fig. 14. Unbalance Response of the Two-Disk Rotor, Weights Out-of-Line, Light End
Position
MTI-2183
7.2.50
A
1
END POSITIONS
0
I i 0.33
0 t 0.33
h I ' •'
iO.33
\
\
A
(
i\
-
T \
ml
A
V
o
<^
THEORETICAL RESULTS
/
• EXPERIiENTAL DATA
2 o » COUPLING END
/
A A FREE END
/
i
<
J ^\
J
A
0 i\ iI 12 16 20 24 28
SPEED, R P i X 10'
1 1 1 1
A
BEARING POSITIONS
1 1 1 1
i, 0.33
f 0.33
" A - A
.
1( 133
1
&
f THLORtTICAL RESULTS
EXPERIMENTAL DATA
o® COUPLING END BEARING
1
4
©
•
/ \
J X•\
M
^
r1 n V
V
f b.
\ \
1 /
I
^
^
'/I
/
/ -'"""—
0 4 £\ 12 HD 2D 24 28
SPEED, RPi t 10^
®
1 0.33 t 0.33
K A
o
I 0.33
j
-™ THEORETICAL RESULTS
L_ ^ o ® EXPERIMENTAL DATA
,
(f
ml \
V®
\ I[s
i
/i
>
\ .
1 X.
J
I — _
1
1'
0 A
J i^ 12
SPEED, R P i K 10'
16 20 24 28
>
by
Robert W. Dickenson
Atomics International
Canoga Park, California
SUMMARY
INTRODUCTION
Fast breeder reactors of power levels of 1000 Mwe or greater are presently
considered to be necessary for the continued conservation of nuclear
resources in an expanding nuclear power economy. Component requirements
in the large future systems transcend any experience now available; the
largest sodium cooled reactor built to date has a design power of
300 Mwt, with components sized to meet this power generation capability.
Therefore, an order of magnitude extrapolation in component size for
power generation at about 2500 Mwt is a near term objective.
Sodium, because of its high chemical reactivity, and its rather poor
lubricating qualities is not suitable for anti-friction bearings. Conse-
quently, sodium systems have demanded bearings lubricated by the
*This work was supported by the U.S. Atomic Energy Commission under
Contract AT(ll-l)-GEN-8
8.1.2
working fluid and rotating seals to retain inert gas over the sodium wher- y
ever mechanical shafting had to be introduced into the heat transfer
systems, almost since the inception of sodium cooled reactors. Sodium
also becomes highly radioactive in the primary circuit. While this does
not affect the bearing performance, it does introduce stringent require-
ments for reliability, accessibility and repairability. In small plants,
the requirements could sometimes be dodged by devices such as frozen
sodium shaft seals which permitted conventional bearings to be used out-
side the sodium, with varying degrees of success. However, with the
advent of even the moderately large plants constructed to date, it has
become apparent that at least for pumps the impeller and shaft must be
stabilized by sodium lubricated bearings. This automatically requires
the use of a rotating shaft seal to retain inert gas over the sodium in the
pump casing; we are faced again with a large extrapolation from current
practice,
CURRENT APPLICATIONS
Bearings
The largest sodium pumps built to date have been installed in the
Enrico Fermi Fast Breeder Reactor, with a capacity of approximately
11,000 gpm. These have been operated quite satisfactorily, indicating
that the principle of sodium lubricated bearings is a sound one. The
bearing dimensions for these pumps are 12-in. diameter and 14-in.
length. The six pumps installed at the Hallam Nuclear Power Facility,
of 7500 gpm capacity and identical bearing size, have operated equally
satisfactorily for a total of 70,000 hours on all pumps.
Seals
An inert gas blanket must be maintained over the sodium pool within
the free surface pumps. This gas is normally argon, in fast reactors,
which is maintained at 1 to 5 psi over atmospheric pressure, to permit
outleakage of argon rather than inleakage of air in case of seal failure.
In addition to minimum inert gas leakage, the rotating seal must not per-
mit any contaminants to enter the sodium stream. Oil lubricated face
type seals must be carefully designed so that the lubricant is trapped
before it can leak down the pump shaft, and the oil volume should be
arranged such that the traps can contain any conceivable amount of
leakage. Carbon is quite deleterious to sodium systems, carburizing
and embrittling thin stainless steel sections such as fuel cladding, and
its introduction into a sodium system must therefore be avoided. It is
possible to find non-carbon bearing lubricants; the effects of their
decomposition products in sodium has not, to my knowledge, been com-
pletely explored. The simplest course is therefore to assure that
potential contaminants cannot escape from the seal into the sodium.
Bearings
Seals
by
Donald R. Nixon
Westinghouse Electric Corporation
Atomic Equipment Division
Cheswick, Pa.
SUMMARY
A review of sodium pump and sodium cooled reactor technology was made by
Westinghouse in 1963 under U.S. AEC Contract No. AT(30-1)-3123*. A purpose
of this review was to establish criteria and limits for the sodium pumps which
will be required for large (1000-1500 MWe) sodium cooled power reactors. This
study indicated possible pumping requirements of 5000-10,000 hydraulic horse-
power per pump.
DISCUSSION
The conceptual design shown Figure 1 depicts a 6000 HP electric motor coupled
to a vertical pump shaft. The shaft extends down through a five foot thick
shielding plug into the pump tank, A free surface of sodium is maintained in
the pump tank by gravity or cover gas pressure balancing. By providing an inert
gas cover above the sodium level in the tank, ingress of air is prevented and
the requirements of the pump shaft seal (located just above the floor level)
are reduced to that of sealing 5-50 psig gas pressure and sodium vapor. Radial
baffling within the pump tank is provided to prevent vortexing of sodium around
the shaft. A labyrinth or series of conical baffles is shown above the free
surface to limit the passage of sodium vapor and to reduce convective heat transfer
in the cover gas area. A sodium lubricated radial bearing is shown just above the
pump impeller at the lower end of the shaft. The upper radial and thrust bearings
are located above floor level and above the gas shaft seal. They are shown as
conventional oil lubricated bearings.
The hydraulic or pump end shown is a double suction pump in which the sodium is
drawn into the top and bottom of the impeller and discharged through a diffuser
into the pump casing or volute. The pump internals (impeller, diffuser, lower
The leading particulars for the pump hydraulic design shown are:
The overall length of this unit is approximately 45 feet and the diamter of the
pump is five feet.
Westinghouse is presently performing parametric and design studies under U.S. AEC
Contract No. AT(30-1)-3762 directed towards development of Engineering data for
large mechanical centrifugal sodium pumps.
Due to the size and speed requirements envisioned for the large sodium pumps of
the future, development of reliable turbulent flow sodium bearings and gas seals
will be necessary.
H
OQ
©
m
o
IM^^ri: 3 S' ,'
n
CD
8.3.1
by
J. Yampolsky
General Atomics Division of General Dynamics Corp.
San Diego, California
and
D. F. Wilcock
Mechanical Technology Inc.
Latham, New York
ABSTRACT
The development of a new gas cooled power reactor, has required the engineering
of a new steam-driven, high-pressure, helium circulator. Because of the
simplifications possible in seals and lubricant systems, water lubricated
bearings were selected.
A shrouded step design was chosen for both journal and thrust bearings. These
designs offer the advantages of a large tolerance to wear, no moving parts,
and the flexibility of running as hybrid or even hydrostatic bearings if desired.
DISCUSSION
The 330 MWe Power Plant being constructed by General Atomic for the Public Service
Corporation of Colorado has the entire primary circuit consisting of the reactor
core, steam generators and circulators, contained entirely within a prestressed
concrete pressure vessel. The output power of the plant is produced by a con-
ventional 2400 psi 1000 F/IOOO F reheat steam turbogenerator which is located
together with the remainder of the associated power plant equipment in a separate
W turbine building. The only penetrations of the PCRV required are those for
power production; feed water in, high pressure steam out, cold reheat steam in,
and hot reheat steam out. The concept of the series steam turbine driven helium
circulator was conceived by General Atomic to best provide the means for primary
helium flow in the reactor power plant. Each of the circulators consist of a
single stage axial flow steam turbine that is in series or tandem with the power
producing steam turbogenerator. These driver turbines are located in the cold
reheat line between the discharge from the high pressure element of the turbo-
generator and the reheater. Figure 1 shows the series steam turbine driven
circulator installed in the bottom head of the PCRV. The machine is oriented
vertically and its housing forms the closure to the primary gas circuit. The
arrangement of the elements comprising the circulator is shown in Figure 2. The
compressor and turbine rotors are overhung from a center bearing housing on
opposite ends of the shaft. The bearing housing which contains the bearings and
seals is supported by the conical section that forms the closure to the primary
system.
The choice of the bearing type and lubrication system is strongly influenced by
the operational and system requirements. The loads imposed on the radial bearings
are those resulting from unbalance since the orientation of the machine is vertica
while the thrust bearing loads result from the pressures acting at each end of the
machine and the axial reactions of the turbo machine elements. It should be borne
in mind that although nuclear plants are, in general, base loaded, they require
starting and stopping and operation under emergency conditions. The thrust loads
for all of the required conditions are identified in Figure 3, where the speed
(plant load) is the abscissa and the thrust loads and pressure are the ordinates.
the lubricant also leaks. The lubricant and buffer gas are then separated in
the accessory system and returned to their respective systems. Moreover,
since the pressures in the two process fluids vary considerably from each other
at part load, large flows of buffer gas would be required or each of the bearings
would have to be pressure referenced to the process it is contiguous to, with
resulting leakage flows in the lubrication system between the two pressure systems.
If, however, one of the process fluids is used in the lubricamt, a single pressure
system referenced to the other process fluid is feasible with a considerable sav-
ings of equipment and complication as compared to the two-pressure lubricant system
and a two buffer gas system. The choice of water as the lubricant allows for this
simplification.
As was pointed out earlier, hydrostatic bearings are required by the operating
conditions so that the boundary lubrication proeprties of oil are not necessary
for startup or shutdown, and it is possible to always have large clearances
separating the possible rubbing surfaces. Another contributing factor to the
choice of water as the lubricant, is that redundancy of the lubrication supply
is easily achieved by the use of the feed water system as an emergency lubricant
supply.
This paper describes the preliminary design for the shrouded step thrust bear-
ing. It is discussed here because it is the critical bearing in the design
because of power loss and thrust capacity considerations.
BEARING DESCRIPTION
Approximate Analysis
Several assumptions were made which permitted a closed form solution to be obtaine
In addition to the usual assumptions involved in establishing Reynold's equation,
the bearing was assumed to be a linear strip, thus ignoring the curvature of the
boundaries and the difference in length between inner and outer edges. Further-
more, pressure across the step (in the radial directiion) was assumed to be con-
stant, so that a one dimensional solution could be obtained. Finally, the corner
areas beyond the step (Bxe) were ignored both for their contribution to flow and
their contribution to load capacity.
The basic differential equaticns were the same by summing flows into and out of the
strip of length dx running across the bearing. Both shear flow and pressure flow
were considered in the direction of motion, and end flow out of the ends of the
strip was assumed to be due to pressure flow only.
Because these bearings were to operate in the turbulent regime, with Reynolds'
numbers as high as 50,000, provision was made to introduce an effective viscosity,
corrected for the influence of turbulent on viscosity.
^ - M — ^ ^ = 0 (1)
dx^ ^^ ^^z ^2
where
|a is the effective viscosity in the x direction over the step and |i
X-ii z
is the effective viscosity in the z direction over the lands. u^
•' xl
is the effective viscosity in the x direction over the lands.
8.3.5
Equation (1) applies along the step. Assuming a linear pressure drop over the
lands, Equation (1) maybe integrated and with the appropriate boundary conditions,
giving the following relationship for the load:
1 e . ,
W = + 2 i s^" h ^
(2)
^ + ^ sinh A ^ + ^ cosh
where
3 -, 1/2
2E! / ^ I \ I \
^ = bL (3)
H = (4)
Equation (2) is conservative in that it does not include the load carried by
pressure in the feed groove itself or over the land between the drain groove and
the feed groove.
Calculated Performance
It was desired to have this bearing fully operable over the full speed range from
30 percent to 100 percent of full speed. Bearing proportions were analyzed there-
fore at the worst load condition at 30 percent speed.
The effect of the width of the side shroud is shown in Figure 5, holding the total
radial width of the bearing constant. While the maximum load is obtained at a land
width of about 0.3 inches, the fluid flow through the bearing increases rapidly as
the land width is decreased. The dotted curve in Figure 5 shows the load per unit
flow, perhaps a better measure of overall bearing performance. Therefore, a value
of land width of 0.475 was selected as a compromise between load capacity and flow.
The step depth is also a critical design parameter. If the film thickness, h^, is
held constant the step depth is gradually increased, the load carried by the bear- •I
ing increases rapidly at first then reaches a optimum value at approximately 5 times
the film thickness. The optimum point is a function of feed pressure, film thick-
ness and the other design variables. A value of 5 was selected for this design,
corresponding to the optimum at the 30 percent speed point with the desired minimum
film thickness of .002 inches.
The length of the land beyond the step, e, cannot be optimized with these approximate
equations. These equations indicate an increasing load capacity as small e is
increased. However, due to the approximations made in eliminating the corner flows,
it is very doubtful whether this is the case. For this reason, small e was chosen
at an arbitrary value that would provide reasonable flow rates, and was roughly
comparable with the side land width.
The pressure generated at the step is shown in Figure 6 for two supply pressure.
Here the step depth is kept constant at 0.012 inches and the h /h^ ratio is varied
by varying the film thickness h^ . It is seen that the pressure at the step in-
creases rapidly with h„/h^ for the low supply pressure. This calculation, made
at 110 percent speed, shows one evidence of the hybrid or partically hydrostatically
nature of this bearing. At small values of h„/h^ the step pressure is very little
greater with a high supply pressure than with a low supply pressure, suggesting that
pressure drop along the step is a strong factor when a high supply pressure is used
in the feed groove.
In Figure 7 the pressure at the step is shown as a function of speed, with both the
step depth and the bearing film thickness held constant. With a low supply pressure
the pressure at the step increases regularly with speed at a rate faster than linear.
For a high supply pressure, 250 psi, the pressure at the step is nearly equal to
the feed pressure at 0 speed drops to a minimum of about 1/2 feed pressure at 4000 rpm
and then increases again. This again is evidence of the quasi-restrictor action
of the step towards lubricant flow in the bearing.
the first part of the paper, because in the course of design the load curve has
shifted. The solid curve for film thickness at an inlet pressure of 250 psi
shows that at no time does the bearing go below the desired minimum film thick-
ness of .002 inches. The solid curve marked flow in gallons per minute gives
the flow through a single pad. Note that the flow is very high at full speed and
again is very high at very low speeds.
Because it was noted that high flows occurred where the film thickness was much
greater than the minimum required, the bearing was redesigned to include a flow
controlled restrictor in series with the bearing. Using a higher, but available,
supply pressure of 600 psi plus a restrictor, the film thickness and flow rate
curves shown by the dash lines were obtained. By this means, the flow rate was
kept nearly constant over the entire speed range, and the film thickness was
maintained above the desired minimum value from 30 percent speed to full speed.
A single restrictor in series with the six bearing pads was used rather than an
individual restrictor per pad. Thus the restrictor served primarily as a constant
flow control device, reducing the speed pressure at the bearing in the speed
ranges where large flows would otherwise be obtained.
While the approximate solution discussed above takes account of turbulence effects
in the bearing, it ignores the influence of other inertia affects. These inertia
affects can be expected to manifest themselves in three ways. These are Bernoulli
pressure drops at steps or entrance areas, ram effects at entrance, and centrifugal
effects. The scale of these effects can be large in bearings of this size. For
example, the Bernoulli pressure drop from an entrance groove into an area with an
average velocity equal to half the runner velocity, is more than 100 psi. For
these reasons 5 a more exact analysis is underway in which these inertia effects
are accounted for in every point in the film. The solution is realized on a
digital computer in two dimensions (R,p) using an iterative technique.
The inlet, exit and transition boundary conditions have been carefully considered.
At a transition such as the step, conservation of total pressure has been main-
tained. At the exit boundaries, the static pressure in the film is made equal
to the external ambient pressure, on the assumption that there is no sensible
recovery of the velocity head possible. At the inlet, no velocity head correction
is made on the assumption that the fluid exiting from the land between the drain
and the feed groove into the feed groove forms a boundary layer which increases
in thickness due to energy imparted by the runner as it traverses the feed groove.
In fact, one can hypothesize everything from a Bernoulli pressure to a full ram
pressure rise at the entrance to the step. Until experimental information is
available, the compromise assumption of no pressure change is made.
A few results from this program are now available. Figure 9 shows a pressure pro-
file along the circumferential mid plane of the bearing for a Reynold's number of
2500 and for a step depth equal to the film thickness. In this representation,
the feed groove is at Station 3 and the step at Station 8. The drain groove is
at Station 0, 10. The lower, solid, curve is the solution with no inertia terms
but including the turbulent effects. The upper curves show the solution with the
inertia correction, including the sharp drop at the step. It can also be that the
values obtained on the first iteration are quite close to the final converged
solution. The centrifugal affect is shown in Figure 10. Here under the same
conditions as in Figure 9, one can compare the inner half profile with the outer
profile. Perhaps surprisingly, the outer half profile is lower than the inner
half profile, indicating that the influence of the longer leakage path overrides
the centrifugal effect. The influence of Reynold's numbers is shown in Figure 8,
in which the dimensionless load is plotted as a function of Reynold's number.
Note that while the dimensionless load appears to decrease as Reynold's number
increases, the true bearing load will increase with Reynold's number since it is
proportional to the square of the angular velocity. Figure 11 shows the curve
with no inertia as the lower solid line^ and the first term correction as the
upper dash line. The middle curve showing the conversion solution has been ob-
tained up to a Reynold's number of 2500, and appears to be reasonably extrapolatable
beyond that point. Thus the first term correction may serve as a guide at Reynold's
numbers for which full conversion cannot be obtained.
8.3.9
PLANNED PROGRAM
It is planned to carry out the experimental program on a test rig in which the
inertia of the circulator rotor is simulated. In fact, the actual shaft and discs
will be used. This experimental work will provide information on the bearings,
seals and rotor dynamics. Figure 13 shows the contemplated measurements. There
will be five planes in which the radial displacement of the shaft relative to the
stator and the axial displacement of the thrust runner relative to the thrust
bearing will be measured. In addition, the static pressure profiles in one of the
pockets of both the thrust and journal bearing will be measured. In addition,
the flows, temperatures and pressures in the seal and lubrication system are to
be observed.
CONCLUSIONS
The step bearing appears to have particular advantages for this installation which
requires a bearing capable of long life with no attention for 10 years or more
while operating in a high pressure water environment. The step bearing offers
the decided advantage of excellent hydrodynamic performance in addition to capa-
bility as a hydrostatic bearing. Because of the relatively large step depth,
and the small rate of variation in load capacity with step depth when operating
near the optimum point on the step depth in the load versus step depth curve,
the step bearing can be expected to be very insensitive to wear. Thus, in the
bearing under discussion here, wear of one to two mils can be tolerated without
seriously affecting the bearing performance. Although inertia effects can play
a major role in the load capacity of this bearing when operating at high Reynold's
numbers and the nature of the boundary condition at the inlet of the step is an
unknown area which can benefit from experimental investigation, it is possible to
meet most any exigency that might occur through a modification of the hydrostatic
component of this type of hybrid bearing.
SERIES STEAM TURBINE CIRCULATOR
PCRV INSTALLATION
.DIFFUSER
COMPRESSOR
STATOR•
.COMPRESSOR
ROTOR
PCRV LINERv
STEAM
TURBINE
CASING-
-PRIMARY
CLOSURE
WATER
INLETS
PRIMARY CLOSURE
INLET STEAM /
8000 — 800
PRESS.
\
INLET He \
y^
6000
PRESS. \
> -
^\ .——
7—z
y^ — 600
^ ^
O /
7^
START-UP
THRUST ^OUTLET
4000 STEAM 400
CIRCULATOR PRESS.
DOWNWARD
PRESSURE
AXIAL \ EMERGENCY (PSIA)
THRUST DEPRESSURIZEXL 2 0 0
(LBS) 2000
DOWN y \ THRUST \
THRU ST
\ \
0 \ ^>
\
11 AXIAL THRUST
DURING NORMAL
^A
UP \
THRUST OPERATION
2000 200
\ -
20 40 60 80 100
% SPEED
Drain
Tied
P
400
300
W
(LBS)
200
0
0 05 1.0
b
400
PS!
300 Pi == 250/
A •20
200
/
A
/
/
100
/
^
Q
0 2 3 4 5 6
h2/h|
300
Pi = 250
^ /
PSI
200
"^ Pi = 20
100
i 1 f 1 l i l t 1 J I ! .
0,
0 5000 10,000 15,000
RPi
Figure 7: Pressure at Step Versus Speed
8.3.1
40 10
5x93/4 x6PAD THRUST
h,p;=250
A = 0.012 WATER
w
KIPS
30
h
MILS
Q
GPi
20 0
10
_ - . - » » — ^ ' ^ ^ ^ p . = 6 0 0 + RESTRICTOR
0
0 5000 10000 15000
RPM
U3
Re-2500, R/C-IOOO, Pg-0.0844
A / C - l , R-0.5125
0.25 (APPROXIiATELY AT THE MIDDLE
OF THE POCKET)
0.20
3
ICL
4 6
NUMBER OF STATION
0.12 \
\
3
^
ICL • • - « - . ^
0.08
NO'INERTIA^v X .
__. (BOTH HALVES) O c 1
^ ^ ,
0.04
^TfP
0 1 cr
Ota
•• LAN!)J \ \
,
\°
\
0
0 I 2 3
NUMBER OF STATION
V X ^X
^
V
\
\ ^ N_
\ \ _. laiiTi 1 1 s
S»w.
\ »«.. V^ %!, wn n 1ST 1 trtivi turtrttt 1 luiM
/
NO INERTIA^ "X, X.*^' ^**^^^»^^^^ /
^^
^"•^Z^ " " " ' ! . ' - - ^»» « "" «»
i|^2275
00
Figure 12: Journal and Thrust Bearing w
N3
MTl-2264
SNSTRUMENTATiON BEARING TEST RIG
EXTERNAL
RESTRICTOR
^ - T A P E RECORMR
I CH. 8 4.9
- -OSCILLOSCOPE
THRUST BEARING
PRESSURE. TAPPING
-STEAM TUHBINC
INERTIA DISC
— TAPE RECOWJER
CH. 10 & 11
i/JP\ 1___JT^
V i ^ 2 4 PT. V l y 12 PT.
RACK RACK
MTO MTO
NOTES
(, TWO DISPLACEMENT PROBES ARE INSTAU-EO
REPORT DISTRIBUTION
•
DISTRIBUTION LIST
Rocketdyne Division
Nucleonics Subdivision
6633 Canoga Avenue
Canoga Park, California 91304
Atten: Robert S. Siegler (1)
•
ADDITIONAL DISTRIBUTION LIST FOR REPORT NYO 3363-6
Clarence E. Miller
Reactor Engineer
Division of Reactor Develop. & Tech.
Components Branch, U.S.Atomic Energy Coram.
Washington, D.C. 20545
Robert Spies J.R. Wetch
Supervisor Atomics International
Rocketdyne P.O. Box 309
6633 Canoga Ave., Canoga Park, California 91304 (1)
Canoga Park, California 91304 (1)
J. Yampolsky
Professor W.K. Stair Staff Member
Dept. of Mechanical & Aerospace Enginee General Atomic
The University of Tennessee P.O. Box 608
Knoxvtlle, Tennessee (1) San Diego, California (1)