Cryocoolers - Part 1 - Fundamentals (1983, Springer US)
Cryocoolers - Part 1 - Fundamentals (1983, Springer US)
Part 1: Fundamenta ls
THE INTERNATIONAL CRYOGENICS MONOGRAPH SERIES
H. J. Goldsmid
Thermoelectric Refrigeration
G. T. Meaden
Electrical Resistance of Metals
E. S. R. Gopal
Specific Heats at Low Temperatures
M.G. Zabetakis
Safety with Cryogenic Fluids
D. H. Parkinson and B. E. Mulhall
The Generation of High Magnetic Fields
W. E. Keller
Helium-3 and Helium-4
A. J. Croft
Cryogenic Laboratory Equipment
A. U. Smith
Current Trends in Cryobiology
C. A. Bailey
Advanced Cryogenics
D. A. Wigley
Mechanical Properties of Materials
at Low Temperatures
C. M. Hurd
The Hall Effect in Metals and Alloys
E. M. Savitskii, V. V. Baron, Yu. V. Efimov,
M. I. Bychkova, and L. F. Myzenkova
Superconducting Materials
W. Frost
Heat Transfer at Low Temperature
I. Dietrich
Superconducting Electron-Optic Devices
V. A. Al'tov, V. B. Zenkevich, M. G. Kremlev, and
V. V. Sychev
Stabilization of Superconducting Magnetic Systems
G. Walker
Cryocoolers, Part 1: Fundamentals
Cryocoolers, Part 2: Applications
Cryocoolers
Part 1: Fundamentals
Graham Walker
The University of Calgary
Calgary, Alberta, Canada
ISBN 978-1-4899-5288-2
The rapidly expanding use of very low temperatures in research and high
technology during the last several decades and the concurrent high degree
of activity in cryogenic engineering have mutually supported each other,
each improvement in refrigeration technique making possible wider oppor-
tunities for research and each new scientific discovery creating a need for
a refrigerator with special features. In this book, Professor Walker has
provided us with an excellent exposition of the achievements of this
period, the fundamental principles involved, and a critical examination of
the many different cryogenic systems which have led to a new era of
low-level refrigeration.
I feel fortunate to have had a part in the developments discussed in
this book. During the early 1930s I constructed several rotary engines
using leather vanes. Their performance was not good, but I was able to
liquefy air. I had been impressed by the usefulness of leather cups in tire
pumps and in Claude-type engines for air liquefaction. I was trying to find
a way to avoid that part of the friction generated by a leather cup as a
result of the radial force of the working gas on the cylindrical part of the
cup. During the 1950s I built two efficient helium liquefiers in which
essentially leather pistons were used. A steel core was encased in a stack
of leather rings separated by thin steel rings of slightly smaller diameter,
the stack of alternate leather and steel rings being compressed to provide
a rigid piston which could be machined to fit the cylinder. The wearing
quality was excellent. One of these liquefiers provided thousands of liters
of liquid helium during about two years of use. Instead of becoming
smaller in diameter because of wear, the diameter of the pistons actually
increased by absorption of water. The necessity of controlling the water
content was a disadvantage.
I was also intrigued by Heylandt's crowned piston concept during the
1930s but ruled it out because of the sealing problem at the warm end of
the piston. I had seen some of them in action in oxygen plants. A slight
vii
viii Foreword
S. C. Collins
Preface
The kindest thing to be said about this book is that it is like the curate's
egg-good in parts. Many times in its composition I have felt I would be
better engaged in reading rather than writing a book on cryocoolers.
Perhaps my modest effort will stimulate others better qualified than I to
do the job.
Completion is due entirely to Klaus Timmerhaus, Associate Dean
of Engineering at the University of Colorado, and to Dr. Alan Clark of
the National Bureau of Standards, Boulder, Colorado, joint editors of the
International Cryogenics Monograph Series published by Plenum Press.
They saw a draft of my first book on Stirling-cycle machines and flattered
me with an invitation to write another on cryocoolers. With the starlight
in my eyes I signed the contract that has been their rod to beat me with
since. Now at the end I am most grateful for their unquenchable interest
and enthusiasm despite the discouragements with which I have confronted
them. They have proved the burr in my saddle that only a completed
manuscript will remove.
Luster has been added to my effort by the substantial contributions
of others. I am indebted first to Samuel Collins for the Foreword. Ray
Radebaugh of the Thermophysical Properties Division of the National
Bureau of Standards, Boulder, Colorado, contributed Chapter 11, dealing
with his interest and expertise in the fundamentals of cryocooling and
Chapter 12, dealing with the systems and techniques for achieving very
low temperatures.
Fred Chellis, senior Applications Engineer of Cryogenic Technology
Inc., Waltham, Massachusetts, has contributed Chapter 10, dealing with
the practical problems of cryocooler design and operation. Fred was highly
qualified to do this, having been around cryocoolers longer than he cared
to remember and really should have been the author of this book. He died
unexpectedly during the production of this volume; his death was a great
loss to the cryocooler community.
xi
xii Preface
come quickly to mind. More recently I have benefited much from mutual
exchanges with Bill Martini, Costa Rallis, and E. H. Cooke-Yarborough.
My colleagues at the University of Calgary withstand, good-humoredly,
my preoccupation with regenerative machines and my insistence on discuss-
ing with them matters I suspect are not their principal interest. In particular,
John Kentfield is a worthy foil off which I have bounced many ideas. The
Head of the Department, Dr. Peter Glockner, has helped in his encourage-
ment and advice.
I wrote this text as my principal activity on a six-month sabbatical
leave from the University of Calgary. I am most grateful to the University
for the opportunity to devote myself to the task unhampered by my normal
academic duties.
Many people helped in the production of the book. Bert Unterberger
and his assistants worked tirelessly on the diagrams and illustrations con-
tained herein. Karen Undseth, Edie Schulz, and Pamela Appleton did a
good job transforming my chicken-scratching info a readable text. Marlene
Stewart and Karen Undseth labored indefatigably in their customary fash-
ion on the index, corrections to the original text, and the galley proofs
thereby relieving me entirely of this onerous last lap. I am most grateful
for all their help.
My children Josephine and Christopher have been deprived of their
rightful allocation of my time but seem to have survived in good order.
Finally I owe my greatest thanks to my wife Ann. For over twenty
years now she has heard all about Stirling and other engines and still
manages a credible show of interest. With unfailing good humor she
somehow, in the end, makes it all worthwhile.
G. Walker
Calgary, Alberta
Contents of Part 1
Chapter 1 Introduction
Definition . . . . . . . . 1
Classification of Cryocoolers 1
Temperature Level . . 1
Refrigeration Capacity 1
Rotary Machines . . . 3
Reciprocating Machines 3
Mixed Units . 4
Large Systems 4
Small Systems 4
Heat Exchangers 5
Flow Regulation 5
Applications of this Text 6
Historical Development 6
Recent Development 10
Status Surveys 20
References . . . . . 26
Cooling . . . . 103
Air Cooling 104
Water Cooling 105
Spacecraft Radiative Cooling 105
Electrial and Electronic Systems 106
Drive Motors 106
Brushless de Motors 108
Electric Controls 109
References . . . . . . 110
Bibliography
General Reading 0 0 315
Government Reports 350
Relevant Conferences 367
Patents 367
Appendix I Glossary of Terms for Cryocoolers and List of
Organizations 0 0 0 0 0 0 0 0 0 o o o o o 375
Appendix II Organizations Having Substantial Interest in
Cryocoolers and Cryocooler Manufacturers 385
Organizations 385
Manufacturers 385
Introduction 0 0 0 0 387
Government Reports 387
NTISearches 0 0 388
Superintendent of Documents (SupDocs) 388
The Cryogenic Data Center 0 0 0 0 0 0 389
xxvi Contents of Part 2
Namelndex . . 395
Subject Index 357
Chapter 1
Introduction
DEFINITION
CLASSIFICATION OF CRYOCOOLERS
Temperature Level
LEVEL
6 5 4 3 2
700
;;;
~
~ 600 30
•
~
NOTE CHANGE
Of SCALE -........
z
'?
~
500
'""
\
~
~.
<r
"<r
~
\ \
w 400 20
<r
~
0
~
~
• 300
~
<r
w
~
0
1\ \
\
200 10
<r
"
~
0
<r
<
<r
100 Fig. 1.1. Theoretical ideal work to
generate a unit of refrigeration as a
• ~
0
""' 10 25 60
0
120
function of the refrigeration tem-
perature (assumed to be the Carnot
cycle with a sink temperature of
REFRIGERATION TEMPERATURE (K)
300K).
Refrigeration Capacity
For the purposes of classification, another important parameter is
refrigeration capacity. It is necessary to specify not only the refrigeration
capacity, but also the temperature at which it is associated. None of the
cryocoolers known at present can operate effectively over the whole cryo-
genic temperature range. Therefore, for classification and reference, bench-
mark temperature levels of 1 K, 4 K, 20 K, 80 K, and 120 K are proposed.
These correspond approximately to the liquefaction temperatures at normal
pressure of: helium (4 K), hydrogen (20 K), nitrogen (80 K), and methane
(120 K).
A possible classification scheme for cryocoolers is shown in Table 1.1.
In this classification five ranges of cooling capacity are prescribed: micro-
miniature, miniature, small, intermediate, and large.
The principal means of achieving cryogenic refrigeration in the level
1, 2, 3, 4 temperature range (120-4 K) is the expansion of compressed gas
in some form of engine. There are many varieties of cooling engines. They
may be broadly classified as
i. Rotary machines.
ii. Reciprocating machines.
Introduction 3
Cooling capacity
Rotary Machines
Rotary machines, turbines, and other forms of rotary expander, are
characterized by small size and weight, high speed, and high fluid flows at
relatively low compression ratios (maximum pressure/minimum pressure).
Rotary machines are relatively inefficient in small sizes but improve with
increases in size and speed. They are relatively inflexible for operation
away from the design condition and so are best suited for continuous,
steady-state operation. There are few bearings and no rubbing surfaces in
contact so that reliability and life expectancy can be very high. Rotary
machines are, inherently, in near perfect dynamic balance.
Reciprocating Machines
With reciprocating machines the effect of size is much less dramatic.
To a first approximation they are as efficient in small sizes as in large. Of
course the relative significance of the mechanical and thermofluid losses
tends to decrease as the size of the machine increases.
Reciprocating engines can operate over a very wide pressure range
and operate quite well even when conditions are reasonably removed from
4 Chapter 1
the design condition. They are relatively large, heavy machines and, to
minimize friction and wear, must operate at low speeds, often an order of
magnitude less than that of rotary systems. The increased number of moving
parts and various rubbing seals and sliding bearings involved decrease the
prospects for high reliability and long life so that relatively high maintenance
costs might be anticipated for reciprocating systems. The out-of-balance
reciprocating forces must be balanced in some way or absorbed in relatively
complicated foundation or mounting installations.
Mixed Units
A cryocooler generally comprises a compression unit and an expansion
unit. Both can be rotary or reciprocating or can be mixed. It is not
uncommon to find in intermediate and the smaller large-scale cooling
systems a combination of reciprocating compressor and turbine expander.
This is because the process of expanding a gas in a turbine is immeasurably
easier to achieve than compressing one. In a turbine expander the gas
readily decreases to a lower energy level. In a rotary compressor the gas
must be 'coaxed" by ingenious aerodynamic design to ascend in pressure
against all natural inclination. The "coaxing" to higher pressure can be
accomplished in more brutal fashion with a reciprocating compressor and
much greater compression ratios can be achieved.
Large Systems
Despite the above difficulties, the advantages of reduced size, weight,
and maintenance of rotary compressors compel their adoption in large
systems to an increasing degree. It is fundamentally a question of develop-
ment effort. Given the appropriate level of commercial application there
is little doubt that satisfactory rotary compressors could be developed for
all intermediate and large-scale applications. Substantial improvements
have been made and work continues with a variety of compression equip-
ment of the centrifugal, axial flow, and screw compressor variety.
Small Systems
Miniature, small, and intermediate cryogenic cooling systems almost
invariably utilize reciprocating systems. Valiant efforts have been made
with small turbosystems, but these remain of interest only in aircraft,
Introduction 5
Heat Exchangers
Heat exchangers are vital elements of all cryocoolers. Two principal
types of heat exchanger are used: (i) recuperative and (ii) regenerative. A
recuperative heat exchanger is a device where separate flow passages are
provided for the hot and cold fluids. The fluids are separated by a solid
wall and heat is transferred by conduction across the wall. The fluids may
flow continuously or periodically. A regenerative heat exchanger has a
single set of flow passages through which the hot and cold fluids flow
alternately and periodically. The regenerative matrix, often a porous finely
divided mass (of metal wire or spheres), may be thought of as a thermo-
dynamic sponge alternately accepting or rejecting heat as the hot or cold
fluid flows through it.
Both types of exchange may be used in counterflow or in parallel-flow
operation. In counterflow operation the hot and cold fluids flow in opposite
directions. In parallel flow operation the fluids flow in the same direction.
Counterflow operation is much more effective than parallel flow operation
and should always be used except in special circumstances.
The type of heat exchanger used, either regenerative or recuperative,
is a key feature in identifying the particular family of machine. Within this
family various machines may be further separated and classified by the
thermodynamic cycle on which they operate or by some other feature.
Cryocoolers with regenerative heat exchangers include Solvay and
Gifford-McMahon, Stirling, Ericsson, and Vuilleumier engines. Cryo-
coolers with recuperative heat exchangers include Linde, Hampson, Claude,
Collins, and Joule-Thomson engines.
Flow Regulation
Another key feature is the manner in which the flow of working fluid
is controlled, namely, by valves or by volume variations. All recuperative
engines have valves. Therefore, classification by means of flow regulation
is specific to regenerative engines. Here an arbitrary distinction is adopted:
i. Stirling engines where the flow is controlled by volume changes;
ii. Ericsson engines where the flow is controlled by valves.
This distinction between engines with and without valves is important and
has a profound impact on the design and operation of regenerative
machinery.
6 Chapter 1
HISTORICAL DEVELOPMENT
applications. Few details of these machines are known but Collins (1958)
describes them as "a commercial success and certainly matched the
efficiency of any low temperature reciprocating expander available at the
time."
An extensive theoretical treatment of refrigerative turbines was pub-
lished in 1939 by Kapitza. He concluded that a low-pressure turbine
liquefier was preferable to a high-pressure liquefier with reciprocating
expander. Not only was the unit thermodynamically superior but it was
also cheaper to build. Kapitza compared inward flow radial and axial flow
impulse turbines and concluded the former was the preferred machine for
gas liquefaction. Experimental results were presented for an inward flow
radial expander capable of handling 570 kg of air per hour.
The above account of the early history of cryocoolers is a much
abbreviated review of the historical survey given by Collins and Cannaday
(1958). This excellent book will be found in any reasonably equipped
library. Reference to it is strongly recommended.
RECENT DEVELOPMENT
HELIUM PRESSURE
TRANSFER
EXPANSION CHAMBER
PINCH TUBE
IHELIUM FILL]
L- CRANKCASE
Fig. 1.2. Small integral Stirling cooling engine for forward looking infrared applications.
tively large, heavy machines with separate compressor and expander units
characterized by exceptional durability and long life (Higa et al., 1977).
Gifford was a prolific writer and the Bibliography contains many entries
indicative of his exceptional activity.
Gifford-McMahon machines are manufactured by CTI, Waltham,
Massachusetts, and by Air Products and Chemicals Inc., Allentown,
Pennsylvania, (APCI). Ralph Longsworth of APCI, a former graduate
16 Chapter 1
student and favored protege of Professor Gifford, calls the APCI machine
a "modified Solvay engine." It is indistinguishable from the Gifford-
McMahon engine. Those wishing to clarify the distinction should refer to
Chapter 2-Elementary Thermodynamics of Cryocoolers.
The exceptional sensitivity of infrared detectors and superconducting
devices requires that the cold region of the cooling engine be isolated from
mechanical shock and vibration, thermal variation, and electromagnetic
interference. For systems mounted in gimbals to facilitate rotary or swivel
action it is important to keep the mounted mass as small as possible. Much
ingenuity has been exercised in efforts to satisfy these demands. Early
success was achieved by condensing liquid nitrogen on a cold cylinder head
and transferring the droplets along a tube to the cold region by the
evaporative boiling process known as "Leidenfrost boiling." However, this
technique appears to have been supplanted with the development of cooling
engines comprised of separate or "split" expander and compressor systems.
Gifford-McMahon cryocoolers were always of this type. They have
valves, two or four in number, with two connecting lines, one supply and
one return between the expander and compressor units. Miniature split-
Stirling cryocoolers are a comparatively recent development, although the
arrangement whereby a displacer operates in a separate cylinder which is
independently driven dates from the Stirling brothers' second engine in
1827. This was a standard beam engine intended for steam drive but
modified to operate as a hot-air engine. Many later hot-air engines had
separate cylinders and drive arrangements including the Robinson and
Heinrici engines (see Walker, 1980). The Kirk engines of 1860 were also
of this type, as was the free-displacer engine invented by Postle (1873)
and by Ringbom (1905). However, it does not appear to have been widely
used for cooling engines until quite recent times. Credit for the reinvention
of free displacer cooling engines belongs to Walter Higa of the Jet Propul-
sion Laboratory, Pasadena (Higa, 1965). Similar concepts appear to have
occurred to Stuart Horn at the Night Vision Laboratory and to Fred Chellis
at CTI about the same time.
The great advantage of the split-Stirling cryocooler is that it has no
valves and only one connecting tube. There are disadvantages, principally
in moving the displacer in properly synchronized relation to the compressor
piston. Despite this, the unit is sufficiently attractive to justify the present
concentrated development effort.
Another company active in Vuilleumier cryocooler development was
AiResearch Manufacturing Company of Phoenix, Arizona working on
contract for the NASA Goddard Space Flight Center, Greenbelt, Maryland.
Browning et al. (AiResearch, 1974) gave a comprehensive summary report
on this program.
Introduction 17
Glass Reinforced
Plastic Cylinder
Vacuum Space
Radiation Shields
Fig. 1.3. Diagrammatic cross section of
low-speed, low-cost, long-life cryogenic
cooling engine (after Zimmerman and
Radebaugh, 1977).
18 Chapter 1
STATUS SURVEYS
+t
10 2
- f-- .
--
--
f--+-
•
~ fi' "'
'""'- f-~ i
•• •
l>
0 r.· •0 '
~ '
1-lftj 0
I
NEW OLD -
f-- 0 ~ •·18·9K-
-I- l> 9 A - 10 • 30 K _
l>
O¢
-+-II
• · 30-90 K
Lunder Development _
Fig. 1.4. Efficiency of cryogenic cooling systems as a function of refrigerating capacity (after
Strobridge, 1974). (Efficiency=actual COP/Carnot COP.)
22 Chapter 1
the value of the constant from 6000 in 1974 to 10,000 in 1980, 17,000 in
1985, and 25,000 in 1990.
Figures 1.5 and 1.6 can be used to determine the volume and mass
of cryocoolers as a function of refrigeration temperature as shown in Fig.
1.8. In this figure, characteristics are sketched in for refrigerating capacities
of 1, 10, 100, and 1000 W.
Charts such as these have obvious attraction for planning purposes.
It must be recognized their value is rather limited. At first glance there
appears to be a surprisingly close similarity for all the systems. However,
the charts have logarithmic scales and small chart separations can be
differences of several times. For example, consider Fig. 1.4. At 1 W capacity
the range of variation in efficiency is from a low value of 0.4% to a high
value of 10%-a range of variation of 25 times. Similarly, at 1000 W
capacity the range of variation is between high and low values of 50% and
7%, a range of 7 times.
A most important factor not addressed in the Strobridge charts is the
useful life or the maintenance requirements of the different systems. A
cryocooler for aircraft or space use would likely be small, lightweight,
highly pressurized, fast acting, and consequently have a short life or operat-
ing period between maintenance. The same equation relating life, strain,
size, weight, and cost applies inexorably to cooling engines as to all other
machines or biological systems. Only the exact terms of the equation are
unknown and must be defined by experience.
Introduction 23
-gto'~t=~~~~~~~~~=1iit=~ftrt-~~~
u
101
INPUT POWER, kW
Fig. 1.7. Cost of cryogenic cooling system as a function of input power (after Strobridge, 1974).
24 Chapter 1
~ 10°r-~~~~~~~~~-+----~~
:::l
..J
~ 10- 1 r-~"d--~~----="""+-__:::,~----~~
'o-•L-------'-----L____::..______r:-_____L___j
10°r----,----,-----,----.----.-~
-,;
~
;
V>
,o•r-~~~~4-~~r---~~~4-~
"'~ IO'r-----t--~~--~~---f~--4-~
Fig. 1.8. Volume and mass of cryogenic
100 0~--~20,---~--~~---.~0~--~,0~0~ cooling systems as a function of refriger-
ation temperature (prepared from Figs. 1.5
TEMPERATURE (K) and 1.6) (after Strobridge, 1974).
Their findings are summarized on Figs. 1.9 and 1.10, reproduced from
the Jensen report. These show the specific weight (kilograms per watt) and
the coefficient of performance of a variety of cryogenic cooling systems as
a function of refrigerating capacity at 77 K.
The cryocoolers considered include the following:
i. A group of Ericsson cycle systems, Gifford-McMahon, Taconis,
and Solvay;
ii. Joule-Thomson systems;
iii. Brayton cycle systems;
iv. Vuilleumier cycle engines;
v. Stirling cycle engines.
In their report, Jensen et al. presented similar curves for cryocoolers
providing refrigeration at 20 and 4.2 K. They paid particular attention to
the problems of heat rejection by radiator, and assumed this to be included
in the total system.
A feature of Figs. 1.9 and 1.10 is the relatively low specific weight
and the high coefficient of performance of the Stirling cryocoolers. Another
interesting comparison is the order-of-magnitude difference between
Introduction 25
.."
"
0.
'"
..
(I)
(I)
:I;
10 100
CAPACITY (WATTS)
Fig. 1.9. Specific weight of cryogenic cooling systems as a function of capacity at 80 K (after
Jensen et al., 1971).
.."'
0
z
2
a:
...a:
0
"'
...
Cl.
0
,_
z
"'0
...... 001~----~~~~~~----------~~~~~~~~----1
"'0
0
CAPACITY (WATTS)
REFERENCES
Woodard, R. S., Welch, P. H., and Jansson, R. M. (1978). "Manufacturing Methods and
Technology for the Establishment of Production Techniques for a Split-Cycle Stirling
Cryogenic Cooler," Report No. 15181, Contract DAA B07-77-C-0631, U.S. Army
Elec. R. and D. Comm., Fort Monmouth, New Jersey (Martin-Marietta Corp., Orlando,
Florida).
Zimmerman, J. E., and Radebaugh, R. (1977). "Operation of a SQUID in a Very Low-Power
Cryocooler. App. of Closed Cycle Cryocoolers to Small Superconducting Devices," Proc.
of Conf. NBS, Boulder, October 1977, pp. 59-66 (Issued as NBS Spec. Pub. 508, April
1978).
Chapter 2
Elementary Thermodynamics of
Cryocoolers
INTRODUCTION
Matter can exist in three phases: gas, liquid, and solid. It can be
changed from one phase to another by the transfer of energy to or from
the system (the particular element of matter under consideration). It is a
matter of common experience that to melt ice (water in the solid phase)
we supply heat, whereupon the ice is converted to liquid. Further heating
causes the liquid to boil and evaporate to vapor or gas. The difference
between vapor and gas is that in one case (vapor) the temperature is below
and in the other case (gas) is above the critical temperature.
STATE PROPERTIES
namics and others of less interest for specialized applications. The combina-
tion of any two properties is sufficient to completely specify the state of a
pure (single component) substance. The six properties of general interest
are: pressure, volume, temperature, internal energy, enthalpy, and entropy.
Consider now the simple closed system shown in Fig. 2.1 consisting
of a piston and cylinder closed at one end. Assume the piston has no mass
and is free to slide without friction in the cylinder. We assume further
there are no leaks of the working fluid past the piston seal. Finally we
assume that some known mass of working fluid is contained in the cylinder,
say 1 kg. A weight of mass M placed on the piston will cause it to descend
until it comes to rest as shown in Fig. 2.1, a distance L from the bottom
of the cylinder.
Volume
For the system shown in Fig. 2.1 the state property volume, V, is the
total volume of the working fluid contained in the cylinder with diameter,
D, below the piston, i.e., V = (7r/4)D 2 L. We have unit mass of working
fluid in the cylinder so the volume measured is the specific volume.
Pressure
Pressure is the force exerted on the piston to keep it "floating" in the
cylinder a distance L from the bottom expressed as the force exerted per
unit area of the piston. In the system shown in Fig. 2.1 the total downward
force is the gravitational force acting on the mass M on the piston.
Elementary Thermodynamics of Cryocoolers 31
Temperature
The temperature, T, indicates the energy level of the working fluid.
Using the concepts of the kinetic theory of gases one can assume the
working fluid to be composed of molecules consisting of "billiard balls" of
matter either singly or in groups of two or three moving about in random
fashion as well as vibrating and rotating. The temperature is the degree of
activity of the billiard balls. Application of heat causes the billiard balls to
become more agitated and move faster with increased random collisions.
This is perceived as an increase in the temperature.
Temperature is measured on four scales: Centigrade, Fahrenheit,
Kelvin, and Rankine. Two useful temperatures are the ice-point and boiling
point of water at atmospheric pressures. These correspond to ooc and 100°C
for the Centigrade scale or 32°F and 212°F for the Fahrenheit scale.
The lowest temperature that can exist is -273.15°C. This is taken as
the starting point or zero datum of the Kelvin scale. The Centigrade unit
is the same as the Kelvin unit of temperature change. The absolute
minimum temperature (-273.15°C) corresponds to -459.67°F. This is
taken as the zero datum of the Rankine scale. The Fahrenheit unit is the
same as the Rankine unit of temperature change.
Internal Energy
Internal energy, U, is similar to temperature as a measure of the energy
level of the working fluid. According to Joule's law, the internal energy of
a perfect gas varies only with the temperature. In the special case of a
perfect gas the change in internal energy is the change in temperature
multiplied by the specific heat at constant volume, Cv, i.e.,
dU=CvdT (2.4)
32 Chapter 2
Enthalpy
The enthalpy, H, is property of state defined by three other state
properties in combination:
H=U+pV (2.5)
It is a useful measure to facilitate calculations for continuous flow systems
where combinations of the internal energy and the "flow work", p V, occur
frequently.
For a perfect gas the change in enthalpy is simply
dH=CpdT (2.6)
where CP is the specific heat at constant pressure.*
Entropy
Entropy, S, is a property of state useful in calculating the quantities
of energy with respect to temperature exchanged in changing from one
state to another. It is defined by the equation
dS = J(dQ)/T (2.7)
The change in entropy is the integral sum of the quantity of heat transferred,
Q, divided by the temperature at which it is transferred.
Books have been written about entropy and students become confused
by it. Everything that happens causes an increase in the entropy. It becomes
difficult to understand how this incessant aggregation does not result in
the whole world getting top-heavy with entropy. To the newcomer not of
philosophical bent, entropy is best regarded simply as a handy way to
calculate heat transfers.
is characteristic of the gas and equal to CR/ M), where R is the universal
gas constant and M the molecular weight of the fluid.
Many gases follow an equation of state:
pV=ZmRT (2.9)
Z is the coefficient of compressibility, a factor which varies widely as the
condition of the gas changes. It is determined experimentally. More compli-
cated equations of state exist, or the properties may be determined experi-
mentally and summarized in tables of thermodynamic data.
Properties of state for various fluids have been measured over many
levels of temperature and pressure. The results are best presented in tabular
form called thermodynamic tables. The tables are usually in two parts. The
first part contains properties of the saturated fluid, the second part contains
properties of the superheated fluid.
ISOTHERM
~ \
T
D C
DEW-POINT
BUBBLE -POINT
Fig. 2.2. Heating and evaporation process
displayed on the temperature-entropy
plane. s
34 Chapter 2
internal energy than the liquid. Let that state be represented by point C.
As more heat is added an increasing fraction of the working fluid will
evaporate. The temperature and pressure remain the same but the piston
will move a long way up the cylinder. The volume of the evaporated state
is much greater than the liquid state. The heat added to cause evaporation
of the working fluid is called the latent heat (or enthalpy) of vaporization.
In the process of vaporizing we can think of the fluid moving from
point B to point C. At some intermediate point, say D, we can talk about
the fluid having a "vapor fraction" x. This is the fraction of fluid at the
saturated vapor state C. It may be calculated simply as
. length ED Sv -SB
vapor fraction x = C (2.10)
length B Sc- SB
The fluid is never actually at any intermediate state for it can only exist at
either state B or state C. It is impossible for any particular element of the
fluid to be at state D. The intermediate state is simply a concept to help
visualize the fraction of fluid that has evaporated. The vapor fraction x is
simply the fraction of the total mass that has been evaporated from state
B to state C. Eventually all the liquid evaporates to vapor at state C called
saturated vapor. Further heat supplied to the cylinder causes the fluid to
increase in temperature along the line CE. It is then said to be superheated.
The line ABDCE is a line of constant pressure called an isobar. A
line of constant entropy, say FG, is called an isentrope and a line of constant
temperature, say HJ, is called an isotherm. If necessary, we could also show
lines of constant volume (isometric) or constant enthalpy (isenthalp) and
constant internal energy.
Point B is called the bubble point, where vapor first appears when
ascending the isobar from A to E. Point C is called the dew point where
liquid droplets first appear when descending the isobar from E to A.
If a larger mass were to be placed on the piston the pressure would
increase. If the above experiment were repeated, the same process of
evaporation would occur but at different values of temperature and entropy.
This is illustrated in Fig. 2.3. The new high-pressure isobar KLMN is shown
along with the original isobar ABDCE. Evaporation at the higher pressure
occurs at a higher temperature and starts, at point L, at an increased value
of the entropy. The increase in entropy during evaporation (SM -Sd is less
than the increase (Sc- SB) at the lower pressure.
Further isobars could be obtained for pressures greater and smaller
than the two considered so far. As the pressure increases, vaporization
occurs at higher temperatures. The loci of the bubble points BL and the
dew points CM may be drawn to obtain the envelope BLPMC. These loci
Elementary Thermodynamics of Cryocoolers 35
CRITICAL- POINT
~ GAS
7-----~A~R
:::--L~~;•_;.o-:~t--M
-------- .
T / •
.. E
. -------.------·.
.
A
e
D
are called the saturation curves. BLP is the saturated liquid line. CMP is
the saturated vapor line. The area encompassed by the saturation curves
is called the liquid-vapor envelope. Within the envelope, fluid exists partly
as saturated liquid and partly as saturated gas. Fluid at any state to the
right side of the envelope exists as a superheated vapor or gas. Fluid on
the left side of the envelope exists only as liquid.
The apex of the liquid-vapor envelope, point P, is called the critical
point. Properties of state at that condition are called the critical properties,
critical pressure, critical temperature, etc. At the critical point no distinction
can be made between a saturated liquid and a saturated vapor-they are
the same fluid.
Above the critical temperature, no liquid can exist whatever the
pressure exerted on it. The critical temperature is therefore taken as the
breakpoint distinguishing gas from vapor. The critical isotherm QPR is
shown on Fig. 2.3. Fluid above this temperature is said to be gaseous. Fluid
in the evaporated state below this temperature is said to be vapor.
which the system may change from state A to state B. Three such paths
are shown. The most direct ,route between A and B is the straight line
process ACE. An alternative involving some high temperatures is the path
ADB. Another alternative involving low temperatures is the process AEB.
Whatever the path followed the state properties at B will always be the
same. The final pressure, temperature, volume, or entropy is independent
of the process. However, the energy involved in accomplishing the change
depends greatly on the process.
It is useful to recall that entropy was defined as dS = JdQ/T. From
this it can be readily seen that rearrangement of the equation results in
dQ = JT dS, the heat transfer is equal to the integrated area of the T-S
diagram. Thus, the energy transfer of the process can be estimated by
computing the area enclosed on the T-S diagram. In Fig. 2.4 the energy
transfer associated with the process AEB is the shaded area enclosed by
AEBFG. This is the least energetic process of the three shown. Process
ACE would encompass a greater area and process ADB a still greater
area with increased energy flows.
p A B
T
and
dQ=dW+dU (2.18)
This is the mathematical statement of the first law for a nonflow system.
The first law is important to refrigerator operation. To produce
refrigeration, work must be supplied to drive the system. According to the
first law the amount of energy to be dissipated at the upper temperature
is the sum of the work supplied and the refrigeration generated.
The first law is also important in refrigeration where rubbing friction
takes place in the low-temperature region. The rubbing friction will require
input of additional work to keep the machine working and that work will
be converted (degraded) by frictional action into heat energy. Furthermore
the heat is produced in the very region one is trying to refrigerate. The
additional friction work degrades the refrigeration generated and increases
the heat to be dissipated.
The second law of thermodynamics is more subtle. The simplest
statement particularly apt for refrigerators is that "heat will not flow from
a low-temperature area to a high-temperature area without the application
of work." Enunciation of the second law in this form is attributed to
Clausius. It sounds so commonplace as to appear trivial: everyone knows
that heat will flow only from hot to cold and that to produce refrigeration,
work must be done; but one must never be deceived by the apparent
triviality of the second law.
Concepts for machines which contravene the two laws of thermody-
namics are called perpetual motion machines-perpetua mobilae. Those
which contravene the first law are called perpetual motion machines of the
first class. Those which contravene the second law are of the second class.
Angrist (1968) in a brilliant review of perpetual motion machines provides
examples of exceedingly subtle contraventions of the second law.
Elementary Thermodynamics of Cryocoolers 39
v s
:r\
PI:
II
~4I
·I I
I: I
v
Fig. 2.6. The Carnot cycle. The series of diagrams on the left side are the P-V and T-S
diagrams for operation with a gaseous working fluid. The diagrams on the right side are the
P-V and T-S diagrams for operation with a condensing working fluid. (a) Process 1-2,
isentropic compression; (b) process 2-3, isothermal compression; (c) process 3-4, isentropic
expansion; (d) process 2-3, isothermal expansion; (e) combined diagram for the cycle.
or water available for heat rejection from the cycle. In practice Tc would
be greater than the atmospheric ambient temperature to accomplish the
heat transfer. Temperature TE is the refrigeration temperature and is fixed
by the application. In practice TE would be less than the refrigeration
temperature actually required so as to induce heat to tranfer to the cylinder.
It is thermodynamically (and therefore economically) wasteful to operate
the engine at lower refrigeration or higher rejection temperatures than the
application requires or the ambient conditions allow.
Reflection on Fig. 2.6e will lead to the very important conclusion that
when operating between given temperature limits the cycle with the
maximum refrigerating effect per unit of work done is that with isothermal
heat addition and rejection to the cycle. Given upper and lower temperatures
of operation the best refrigeration performance always comes from the
cycle with isothermal heat transfers to and from the cycle.
COEFFICIENT OF PERFORMANCE
/ [\
50
v
,/EVAPORA T IDN
I
f.----
/
I/ ~
0
~
PHILIPS
' ~
'
~
STIRLING
''
u
"
z
0 30
~ 1-
20
c~
1--..... ~
"'"
Fig. 2.7. Comparison of the efficiency of JOUL THOMSON
domestic refrigeration, air conditioning, and the like. These are basically
Carnot systems with a condensing working fluid and are modified to replace
the expansion engine process.
The work per unit refrigeration increases dramatically as the refriger-
ation temperature decreases. Values calculated for a Carnot cycle system
with Tc of 300 K are quoted in Figure 1.1 (see Chapter 1). The work per
unit of refrigeration at 100 K is 2 watts per watt. At 10 K the ratio is 29
watts per watt and at 1 K it is 299 watts per watt.
Exactly the same cycle may be used as a heat pump to abstract heat
at ambient temperatures and to reject it from the system (say for heating
a building) at higher temperatures. It is customary to quote the coefficient
of performance of a heat pump as the ratio of heat rejected to work done:
(2.24)
This coefficient is greater, by unity, than the coefficient of performance
defined as a refrigerator. This is true even though the same machine may
be working over the same temperature limits, i.e.,
(COPheatpump = Oc/(Oc- QE) (2.25)
(COP) refrigerator= QEj(Qc- QE) (2.26)
so that
(COPheat pump= 1 +(COP) refrigerator (2.27)
It is fashionable to talk about refrigerators as heat pumps even though
they are used as refrigerators. In such unfortunate circumstances it is
necessary to be pedantic about the nomenclature to avoid the evident
problems that will arise without a precise definition of terminology.
44 Chapter 2
The coefficient of performance for the Carnot engine was found above
to be COP= TE/(Tc- TE). It is the maximum value that can be achieved
by a refrigerating engine operating between temperature limits TE and Tc.
However, it is less well known that other thermodynamic cycles can attain
the same performance. In fact, there are an infinite number of thermody-
namic cycles that can achieve the Carnot value. All these cycles incorporate
regenerative heat transfer processes in place of the adiabatic compression
and expansion of the Carnot cycle. The high-performance regenerative
cycles all embrace isothermal processes of heat addition and rejection from
the cycle. These are the necessary key features of high-performance ther-
modynamic cycles.
The most general form of idealized thermodynamic cycle, the Reitlin-
ger cycle, shown in Fig. 2.8, consists of two isothermal processes where
the heat is supplied to and drawn from the cycle and two polytropic
regenerative processes. It is essential in the ideal reversible engine that
external heat exchange occurs only at the upper and lower temperature
limits of the cycle. Heat reception and heat rejection during the polytropic
phase must be provided by thermal storage in a regenerator. This can be
conceived as a thermodynamic sponge, absorbing and releasing heat as
required.
The generalized regenerative cycle was first analyzed by Reitlinger
(1876) a century ago. It can exist in an infinite variety of forms distinguished
by the character of the polytropic regenerative phases. The Carnot cycle
is simply a special case of the Reitlinger cycle in which the regenerative
polytropic phases are replaced by adiabatic, isentropic phases. Two other
special cases have been given names: the Stirling and Ericsson cycles where
[]'
Tc
p T
TE 3 4
v s
Fig. 2.8. The generalized Reidinger isothermal cycle.
Elementary Thermodynamics of Cryocoolers 45
The Stirling cycle illustrated in Fig. 2.9 is similar to the Carnot cycle.
Operation is as follows. Consider a cylinder containing two opposed pistons,
with a regenerator between them. The regenerator may be thought of as
a thermodynamic sponge, alternatively releasing and absorbing heat. In
practice this is accomplished with a matrix of finely divided metal in the
~~.
·~ <:;J.
Tc 2
v (oJ s
REGENERATOR
(2)
(3)
(4)!
_!!J 1 DISPLACEMENT
"' / "' I
I
(I)
v
(2)
/ i/ (3)
~~ (4)
form of wires or strips. One of the two volumes between the regenerator
and the pistons is called the expansion space, and is maintained at low
temperature TE. The other volume is called the compression space, and is
maintained at ambient temperature T c- There is a temperature gradient
(Tc- TE) across the transverse faces of the regenerator and it is assumed
there is no thermal conduction in the longitudinal direction. The pistons
are assumed to move without friction or leakage of the working fluid
enclosed between them.
Assume that the compression-space piston is at the outer dead point,
and the expansion-space piston is at the inner dead point, close to the face
of the regenerator. All the working fluid is then in the compression space
at ambient temperature. The volume is a maximum, and the pressure and
temperature are represented by (1) on the P-V and T-S diagrams, shown
in Fig. 2.9. During compression (process (1-2), the compression piston
moves towards the inner dead point, and the expansion-space piston
remains stationary. The working fluid is compressed in the compression
space, and the pressure increases. The temperature is maintained constant
because heat Oc is abstracted from the compression-space cylinder.
In the transfer process (2-3), both pistons move simultaneously, the
compression piston towards, and the expansion piston away from, the
regenerator. The volume between them remains constant. The working
fluid is transferred through the porous metallic matrix of the regenerator
to the expansion space. In passing through the regenerator, the working
fluid is cooled from Tc to TE, by heat transfer to the matrix, and emerges
from the regenerator into the expansion space at temperature TE. The
decrease in temperature passing through the matrix at constant volume
causes a decrease in pressure.
In the expansion process (3-4 ), the expansion piston moves to the
outer dead point; the compression piston remains stationary at the inner
dead point, adjacent to the regenerator. As the expansion proceeds, the
pressure decreases as the volume increases. The temperature remains
constant because heat QE is abstracted from the cylinder surroundings. This
is the useful refrigeration of the cycle.
Finally both pistons move simultaneously to transfer the working fluid
(at constant volume) back through the regenerative matrix from the
expansion space to the compression space. In passing through the matrix
heat is transferred from the matrix so the working fluid increases in
temperature and emerges at Tc into the compression space.
The cycle is therefore composed of four heat-transfer processes:
Process (1-2)-Isothermal compression; heat transfer from the work-
ing fluid to the external sink at ambient temperature,
Tc.
Elementary Thermodynamics of Cryocoolers 47
s
Fig. 2.10. Stirling and Carnot cycles superimposed with common values for the minimum
and maximum temperatures, pressures, and volumes. Cross-hatched areas on the P- V plane
represent the increased work output of the Stirling cycle. Cross-hatched areas on the T-S
plane represent increased heat transfer of the Stirling cycle.
48 Chapter 2
Fig. 2.11. Ericsson and Carnot cycles superimposed with common values for the maximum
and minimum values of temperature, pressure, and volume. Cross-hatched areas represent
the increased work output and heat transfer of the Ericsson cycle.
TmoK
PRIME
MOVER~.
TE 3._?""-_ __.4
COOLING
ENGINE
s
Fig. 2.12. Stirling cycle as cooling engine and as prime mover. In both applications the
compression process occurs at ambient temperature Tc. The expansion process occurs at a
low temperature TE in the cooling engine and at a high temperature T max in the prime mover.
In both cases heat is absorbed during expansion and rejected during compression.
----u
Tatmos
3
2
by a diaphragm or piston from the working fluid in the Stirling cycle device,
the system can be classified as a Stirling cycle working as a prime mover.
In other instances, where the working fluid itself was the fluid to be
compressed and pumped, there are, invariably, valves or other flow-
controlling devices. According to the defintion given in Chapter 1 these
systems cannot be classified as Stirling engines but rather as Ericsson
engines. Nevertheless, when discussed in the literature they are frequently
called Stirling engines. In most cases, fluid is added to the cycle when the
pressure is low and withdrawn at a higher pressure. Heat is supplied at
high temperature and rejected at ambient temperature. Work during
expansion is greater than work during compression by an amount equivalent
to the "pump" work of the compressed fluid.
In the late 1950s when the Philips Stirling engine liquefiers were
introduced it was common for them to demonstrate the versatility of the
Stirling engine for operation as a prime mover, cooling engine, or a heat
pump. A simplified diagram of the engine used for the demonstration is
shown in Fig. 2.14. The engine was of a type known as the piston-displacer
machine rather than the opposed piston version used above to describe
the cycle.
Elementary Thermodynamics of Cryocoolers 51
EXPANSION
SPACE
FREEZER
REGENERATOR
DIS PLACER
COOLER
,_-,_,__COMPRESSION
SPACE
PISTON
In this engine the space above the dis placer, below the cylinder head,
is the expansion space. Between the underside of the displacer and the top
of the piston is the compression space. The two spaces were in communica-
tion through a water-cooled heat exchanger, regenerative heat exchanger,
and a finned freezer heat exchanger. The machine was mounted on a
bed-plate and coupled to an electric dynamometer capable of operating as
a motor or as a generator in either direction of rotation. Helium or hydrogen
gas was used as the working fluid. A thermocouple was attached to the
cylinder head and a large temperature scale was projected onto the wall
of the lecture room with a moving light spot to indicate the cylinder head
temperature.
The engine would be started initially as a refrigerator with the
dynamometer working as a motor driving the machine in the forward
direction. The cylinder head around the expansion space would cool rapidly
and on the temperature scale one could observe the decrease in tem-
perature. At about 80 K air condensed on the cylinder head and ran into
a collector at the base of the head. After a short time power to the motor
would be switched off and in less than a revolution the machine would
stop and self-start in the reverse direction.
It was now operating as a prime mover generating power. With the
reversal of rotation the expansion and compression spaces switched loca-
tions so the expansion space was between the displacer and piston. Heat
was drawn from the cooling water to drive the engine and rejected from
the cylinder head to the pool of liquid air serving as the engine coolant.
The reservoir of liquid air was soon dissipated and soon the temperature
of the cylinder head started increasing. The engine continued to run but
with progressively less vigor as the temperature increased until finally it
slowed to a very low speed.
At that point power would once more be supplied to the dynamometer
causing it to run again as a motor but preserving the reverse direction of
52 Chapter 2
rotation. With this power input the temperature of the cylinder head
continued to increase until eventually it attained and then exceeded the
ambient temperature. The machine was now operating as a heat pump
drawing heat from the "cooling" water and rejecting it at higher tem-
perature to the surroundings of the cylinder head now acting as the com-
pression space. The temperature of the head would increase rapidly,
whereupon the lighting would be lowered and the cylinder head could be
observed glowing cherry red. At that point Jan Kohler would break his
discourse, reach into his pocket, withdraw a cigar and casually lean down
to light it on the hot cylinder head. It was beautifully executed and done
so matter of factly. Following this, power to the motor would be switched
off and again, immediately, the engine would stop and starting running in
the opposite direction.
With this second reversal the direction of rotation was now as it had
been initially. The expansion and compression spaces were restored to the
configuration shown in Fig. 2.14. The engine now ran again as a prime
mover, withdrawing heat from the hot cylinder head, rejecting it to the
cooling water. After some time the cylinder head cooled and, eventually,
there was insufficient temperature difference between the expansion and
compression spaces and the engine stopped.
It was a remarkable demonstration, beyond the capability of any other
machine. The versatility for the same unit to operate as a prime mover,
refrigerator, and heat pump has practical implication for the future. A
Stirling refrigerator driven by a Stirling prime mover has already been
demonstrated by William Beale at Sun power Inc. Use of a common working
fluid has outstanding attractions to overcome seal problems. The combina-
tion may be a solar, isotope, or combustion heated cryocooler. More likely
in the near future it will be developed as a gas-fired air conditioning
refrigerator. The combination or duplex Stirling engine should not be
confused with the Vuilleumier cycle machines we shall consider later (see
Chapter 4).
Another concept to utilize the versatility of the Stirling engine is to
use it as a prime mover extracting heat from a thermal battery (lithium
fluoride) as an automotive propulsion engine. For downhill operation, or
when braking, the engine could be converted to function as a heat pump
feeding thermal energy into the thermal battery. In practice it is not
necessary to physically reverse the engine. The same effect can be gained
in some forms of engine by simply switching fluid connections (see Walker
1979). The combined regenerative braking/propulsion concept is attractive
in mining applications where trucks are used on loaded downhaul operation
and for underground mines. It will become attractive for road vehicles as
the fossil fuel supplies dwindle and the world moves into the "electric
Elementary Thermodynamics of Cryocoolers 53
economy." The Stirling engine with thermal storage will likely become the
preferred automotive propulsion unit. Commuters will plug in their car
for an electric (thermal) recharge on arrival at work or on arrival at home.
This consumer behavior prototype already exists in Canada, where, in
winter, it is customary when parking to "plug in" an engine block heater.
A. J. Oxley of the U.K. Rutherford Laboratory has explored use of
liquid nitrogen as the coolant for Stirling engines. Engines could be ener-
gized from an ambient temperature source to operate as a prime mover
for underground mining applications with no emission of noxious fumes.
Other possibilities for large-scale power generation exist using the tem-
perature difference between sea water and liquid natural gas now coming
into the United States, Japan, and Europe in very large quantiites.
REGENERATION
CONSTANT CONSTANT
EXPANSION COMPRESSION (1·2)
VOLUME (2~3) PRESSURE (2·2')
(3-4)
~~~-~
7
_ TE
Q IN AT
I' I I=I :I:::-+-:1I~r:·
REGENERATION
Q OUT AT
TE (REFRIGERATION) CONSTANT
VOLUME (4~ I)
Fig. 2.15. Conceptual arrangement of an engine for operating on the Rallis isothermal
regenerative cycle.
T-S planes for the cycle are shown in Fig. 2.16. The cycle is defined by
two isothermals, temperatures Tc and TE, bounded by constant volume
and constant pressure regenerative processes at the two extremes as shown
above. Tc is the ambient temperature at which cooling is available, TE is
the refrigeration temperature.
Processes (1-2)-Isothermal compression, heat is rejected from the
system to the surroundings.
Process (2-2')-Constant pressure regenerative cooling.
Process (2'-3)-Constant volume regenerative cooling.
Process (3-4) -Isothermal expansion, heat is extracted from the sur-
roundings, (this is the useful refrigerating effect of the
cycle).
Process (4-4')-Regenerative heating at constant pressure.
Process (4'-1)-Regenerative heating at constant volume.
For convenience in analysis we define the following nondimensional
parameters:
T = Tc/TE
rc = Vd V2 , volume compression ratio
re = V 4 / V 3 , expansion compression ratio
a= Tz·/Tz, regenerative cooling constant pressure ratio
b = T 3 /T2·, regenerative cooling constant volume ratio
d = T4·/T4, regenerative heating constant pressure ratio
e = TdT4·, regenerative heating constant volume ratio
e = (Ts- Tz)/(T3- Tz), regenerator effectiveness
COP= coefficient of performance defined as the ratio Qj W
W =specific external work required per cycle
Q = specific external heat extracted at the low temperature
(refrigeration)
Pm = W/swept volume, mean effective pressure
P = Pm!Pmax• mean effective pressure ratio of the cycle
Pmax = Pz = Pz·, maximum cycle pressure
Elementary Thermodynamics of Cryocoolers 55
Pmax
Pmin 1----+---"-4-----'4,
Vmin v Vmox
Now
and
(2.34)
and
a :::;; -1)
1 - e ( -'T' T -
RTc
= {(y -1) In re -[r -e(r -1)-1]}
r(y-1)
RTE
= - -{('}' - 1) In re - [ 'T - e ('T - 1) - 1]} (2.36)
( 'Y -1)
Elementary Thermodynamics of Cryocoolers 57
RTE
- 1 {( 'Y -1) In re- [T - e (r -1) -1]}
COP=-y~----------------------
RTE{r[1- a+ n (adre)]- (d -1 +In re)}
1 (y-1)lnre-[r-e(r-1)-1]
(2.37)
('Y - 1) {T [ 1 - a + In (adre)] - (d - 1 + In re)}
Pm = W/stroke volume
RTE{r[1- a+ In (adre)- (d -1 +In re)]}
v4,- v3
RTE{r[1-a +In (adre)- (d -1 +In re)]}
(2.38)
(dre -1) v3
P= Pm =RTE[r(1-a+ln(adre)-(d-1+1nre)]
Pmax (dre -1) v3 X Pz'
P= r[1-a+ln(adre)-(d-1+lnre)]
(2.39)
b(dre -1)
Similarly for case (b) where T 2, :s; Ts :s; T 2 the external refrigerating
effect is
=RTE In re -CvTE[ (D J
-1 -CpTc[ 1-eC~ 1 ) -a J
and if e = 1
1 TE
COP=--= the Carnot value (2.43)
r-1 Te-TE
2/2'
p
.~-
~
'
414 1
Fig. 2.17. Ideal Stirling cycle. v s
Elementary Thermodynamics of Cryocoolers 59
The mean effective pressure of the Stirling cycle is [from Eq. (2.38)]
_ T[ 1 - 1 + In (1 X 1 X r)] - ( 1 - 1 + In r)]
Pm=RTE (r- 1)V3
RTE (T -1) In r
(2.44)
v3 (r-1)
P=C~1)C~~) (2.45)
(y- 1) In r 1 -~-1 ) -
+ (- y[ T - s (T - 1) ~7 J
COP=----~----------~--------~
(y - 1){[1 - 1 +In (1 x 1 x r)]- (T - 1 +In r)}
(y -1) In r+y(T -1)(1-s)
( y - 1) (T - 1) In r
and when s = 1
.'~ 4 1/41
T
D
2'13
2
E
4
114'
212'
D
212' Tc 1/4'
~
3
4
1/4
'
T
3 4
Fig. 2.19. Cycle with constant-
volume regenerative cooling and
constant-pressure regenerative
v s heating.
Elementary Thermodynamics of Cryocoolers 61
(y-1)lnr. [ -(1-1)]
1-
T]
-y [ T-e(T-1)-;
COP=I 1 I
(y -1) { T[ 1-;+ In (r,/T)]- (1-1 +In r,)}
p -1 +In....:....;::'
= T___ (r./T) ---'-_-In____:r._ (2.54)
r. -1
p
2~
~I T
n
2 13~
1 Fig. 2.20. Cycle with constant-
414
414 1 pressure regenerative cooling
and constant-volume regenera-
v s tive heating.
62 Chapter 2
v s
Fig. 2.21. Rallis adiabatic regenerative cycle.
T2 = r~y-1)
therefore
fr~y- 1 ) =1
or
!=1/r~y- 1 )
therefore
g/r~y- 1 l = 1 or g = r~y-1)
T4' = dT1 = dT3 = (d/r)T1
but since T4' = Tde then 1/e = d/r or ed = r.
64 Chapter 2
RT1
=-
'Y -1
d
- [ (yrc<'Y - 1> -ya +a -1)- ( y--y-----=-rr--+-
d
T
1
Tr~'Y T
1)]
T
=
RT1)[
(y- 1 y(rc
(y-1) 1 1 ) ) +;(d-1)
y ( d-,~y-
-a)+(a-1)-; 1 J
=RT1
- - [ y(r/
< _0
-a)+(a-1)- 'Y ) (dr.'Y
(y- 1
< -1> 1
-1)+-(d-1) J
y-1 rr. T
= T~~\) [ Ty(r~y- 1 )- a)+ T(a -1)- r~'Y'Y_ 1 > (dr~y- 1 ) -1) + (d -1) J
(2.55)
= RTE [Ty(r~y- 1 )- a)+ T(a -1)- (y'Y_ 1 > (dr~y- 1 ) -1) + (d -1)]
(y-1) '·
(2.56)
- RTE
p =
m (y-1)(V4·- V3)
so that
- RTE
p =
m ( 'Y -1)(d -1) v3
P = b ('Y _ 1~(d _ 1) [ ry(r~y- 1 ) -a)+ r(a -1)- r~">''Y_ 1 l (dr~y- 1 ) -1) + (d -1)]
(2.58)
Thermal Regeneration
With perfect regeneration 0<7-4'l + 04' = 0(6-2') + 0(2'-3l so that Ts, the
temperature of the gas at the end of regenerative cooling, is Tc or TE.
With less than perfect regeneration, e < 1 (due to heat leaks or friction
effects) or with less than perfect external system cooling (so that T 6 > Tc
at the start of regenerative cooling) then T 5 , temperature at the end of
regeneration, will be greater than T 3 • The energy to reduce the temperature
of the gas from T 5 to T 3 will therefore have to be drawn from the
constant-pressure refrigeration process 04-7·
The heat extracted or available refrigeration will therefore be reduced
to
0 = 0(4-7)- 0(5-3)
1 -e (T- r-1)-
;-s, ~a
66 Chapter 2
Q=(y- 1
1 1 )-r+s(r-1)+;-
RTE ) { y [ 1-,~y- a] (a -r)}
-T-
COP=
( ~~~){
'
y[ 1- ,cY~ 1 l -r +s(r -1)+;]- (a ;r)}
e
Pseudo-Ericsson Cycle
Now consider the special case where the regenerative processes occur
only at constant pressure so that rc = re = r, b = e = 1, (1/ a)= d = T and
T 2 , s T 5 s T 6 only. The P- V and T-S planes are shown in Fig. 2.22. This
cycle may be called (after Rallis et al., 1977) the reversed pseudo-Ericsson
cycle. From Eq. (2.60)
Y[ 1- ,c!-1) -r+e(r-1)+;]-(;-1)
COP=------~--------------------------
ry ( r (y-1) - ;
1) +r (1; - 1) - ,cy-
'Y ( rr (y-1) - 1) + (r- 1)
1)
2- ( -1 -) - + e (T - 1)
,cy-1) T
(2.61)
(1-y)
and when s = 1, perfect regeneration,
(y-1) 1 )
COP= (,c:-1l (1-=- y) (2.62)
68 Chapter 2
1/4'
2'/3
v s
Fig. 2.22. Pseudo-Ericsson cycle.
(2.63)
Pseudo-Stirling Cycle
Another special case arises when the regenerative processes occur
wholly at constant volume. This may be thought of as the reversed pseudo-
Stirling cycle. In such a case, however, it is most unlikely that the heat
transfers to and from the system Q< 2 _ 6 l and Q< 4 _ 7 l would in fact occur at
constant pressure. Rather, it is more likely they would occur at constant
volume and therefore the equations derived above would not be applicable.
Elementary Thermodynamics of Cryocoolers 69
s
Fig. 2.23. Pseudo-Stirling cycle.
R
= - - [(T2- T1)- (T3- T4)]
(y -1)
RT3 [ (y-
= (y-
1) ( 1 )]
1) T(r -1)- 1- ,cy-1)
and since Pmax = T,Cy- 1 lP3 and RT3jp3 V3 = 1 the dimensionless mean
effective pressure ratio becomes
Pm (r(y- 1 ) -1)(T,Cy- 1 ) -1)
p = -- = --.,.27(-y------.-1')- - - - - (2.67)
Pmax Tr (y-1)(r-1)
70 Chapter 2
RTE [,<y- 1) -1
=(-y- 1) ,<-r 1) -y(r-1)(1-s)
J (2.68)
p T
3
4' 4
v (a) s 2
2'/3 5 6 2
p T 2'13
v (b) s 2
2
p T 3
v (c) s
Fig. 2.24. Three adiabatic cycles with various combinations of constant-pressure and constant-
volume regenerative processes.
PARAMETRIC EFFECTS
The above equations provide a basis for useful and interesting studies
of the effects of various parameters of the cycle. By way of example, results
computed for the ideal Stirling cycle are presented in Fig. 2.25. The
compression temperature Tc was assumed to be 320 K and the refrigeration
temperature, TE, 80 K. The figure shows the coefficient of performance as
a function of the volume compression ratio with different values of regen-
erator effectiveness. It is interesting to note the marked deterioration in
coefficient of performance with even the slightest departure from perfect
regeneration, particularly at low volume compression ratios. The volume
compression ratio in actual Stirling cryocoolers rarely exceeds 2. At this
value of volume compression ratio in Fig. 2.25 the COP with a regenerator
effectiveness of 0.9 is only 0.11, a third of the Carnot value (0.33 with
perfect regeneration).
72 Chapter 2
0.4
..,
REGENERlTOR
EFFECTIVENESS
..ob--
--- -
0.3
.ov -
--
1/ /
..oy
I . /
v
1/
I
1" /
~
f..--
I ; / v v
0 .I
I / 1/ ..
Fig. 2.25. Coefficient of performance
of the ideal Stirling cycle as a function
~ f..-- of the cycle volume compression ratio
4 6 8 10 12 with different values of regenerator
VOLUME COMPRESSION RATIO effectiveness.
Similar results calculated for the Rallis adiabatic cycle are presented
in Fig. 2.26. This shows the coefficient of performance as a function of the
volume compression ratio at different values of the regenerator effective-
ness. Temperatures of the working fluid at the beginning of compression
and expansion were assumed to be 320 K and 80 K as before.
An appreciation of the drastic consequence of replacing isothermal
processes by adiabatic processes may be gained by comparing Figs. 2.25
and 2.26. With adiabatic expansion the available refrigeration is substan-
tially consumed in completing the regenerative cooling process and the
effect is progressively cumulative as the volume compression ratio is
increased. It is sobering to realize that in practice, with cryocoolers operat-
ing at normal speeds (say 10 to 20Hz), the compression and expansion
are more likely to be adiabatic than isentropic processes.
In a cryocooler the maximum pressure is important for it determines
the strength requirement of the cylinder wall and the forces acting on the
drive mechanism. To a first approximation the maximum pressure is indica-
tive of the weight of the engine. The mean effective pressure is indicative
of the cycle power input per unit of piston displacement and hence, to a
first approximation, of the refrigeration capacity per unit volume of
machine. The ratio of mean effective pressure to maximum cycle pressure
is therefore an important guide to the refrigeration output per unit weight,
the specific capacity of the cryocooler.
Elementary Thermodynamics of Cryocoolers 73
I l
IDEA~OTHERMAL
I
STIRLING CYCLE
I
----
0.3 e--. 1
.."'
u
z
:IE
n
I
~
a: REGENERATOR
"'...a.. 0.2
EFFECTIVENESS
0 \' € '"1.0
....
z £•0.99
"'
u
...iL " \\ - <•0.975-
1- •. 0.95
~ _..... ~ v
1- <•0925
8
k------- l----
r:::g<;~~
<•0.9
v
<•085 -
0. I
17 <•0.8
I .--
-
1/
II
2
/ 4 6
./
8 10 12
VOLUME COMPRESSION RATIO
Fig. 2.26. Coefficient of performance of the pseudo-Stirling cycle as a function of the cycle
volume compression ratio with different values of regenerator effectiveness.
In Fig. 2.27 the pressure ratio (flm!Pmax) is compared for the adiabatic
and the isothermal reference cycles as a function of the volume compression
ratio. It can be seen from this figure there is a marked and progressive
decline in the output per unit weight for both cycles with increase in the
volume compression ratio. The isothermal cycle is superior to the adiabatic
cycle.
\
0.7
0.6
. 0.5
~
a
£
j
\~
\'
0.4
i 0. 3
~ ISOTHERMAL CYCLE
~~ ..........
-
-
2
~
Fig. 2.27. Comparison of the pressure ratio: 0 .I
ADIABATIC CYCLE
>-----
mean effective pressure/maximum cycle
pressure (flm/Pmaxl for the ideal and the
pseudo-Stirling cycles as a function of the 0
0 4 6 10 12
volume compression ratio. VOLUME COMPRESSION RATIO
74 Chapter 2
Wo Wo W
O.T
(b) Stirling Cycle
Fig. 2.28. Composite diagram of four different regenerative engines (after Radebaugh, 1977)
where 0 represents ambient temperature and h high temperature.
Elementary Thermodynamics of Cryocoolers 75
Stirling Cycle
On the Radebaugh diagram the Stirling cycle is shown in Fig. 2.28b.
It consists of the following:
i. A piston receiving work, W 0 , from an external source during the
piston upstroke and generating work, W, during the piston down-
stroke. When operating as a refrigerator the net cycle input work
is (W0 - W).
ii. A cooler rejecting heat 0 0 at temperature T 0 from the working
fluid to the surroundings.
iii. A displacer mechanically coupled to the piston but operating out
of phase with it so that fluid is displaced between the ambient
temperature compression space and the low-temperature
expansion space.
iv. A regenerative heat exchanger acting as a thermodynamic sponge
and receiving heat from the working fluid when the fluid is passing
from the ambient to low-temperature region. It releases heat to
the working fluid when the fluid is passing back from the low to
ambient temperature region.
v. A freezer receiving heat, Q, or abstracting it from the surroundings
at low temperature, T. This is the useful refrigeration product of
the system.
Vuilleumier Cycle
Another regenerative cycle related to the Stirling cycle is the Vuil-
leumier (pronounced Viamay) cycle shown in Fig. 2.28d. This is is some-
times described as the Stirling cycle with a thermal compressor rather than
a mechanical compressor.
Comparison of the diagrams for the Stirling and Vuilleumier cycles
shows the lower halves for both systems (the low-temperature region) to
be identical. The principal difference is that the Vuilleumier cycle has no
mechanical compressor but has another displacer and regenerator arrange-
ment with an additional heat exchanger called the heater.
76 Chapter 2
Solvay Cycle
The Solvay cycle engine shown in Fig. 2.28a consists of a compressor,
operating at ambient temperatures, compressing gas and receiving work
input W0 from an external source. The heat of compression Oo is rejected
from the system in the cooler at T 0 •
The compressed fluid is expanded in an expansion engine to produce
work output, W, and a low temperature of the fluid. The low temperatures
produced by expansion induce a flow of heat Q to the working fluid at
temperature T, in the freezer. A regenerative heat exchanger is interposed
between the ambient temperature compressor and the low-temperature
expander. The regenerator facilitates the attainment of very low tem-
peratures in the expansion space.
The Solvay engine is virtually identical to the opposed piston Stirling
engine shown in Fig. 2.9 used for the detail description of the Stirling cycle.
The only difference between that engine and the Solvay machine in Fig.
2.28a is the presence of valves to control flow in the Solvay machine. The
Elementary Thermodynamics of Cryocoolers 77
Gifford-McMahon Cycle
The Gifford-McMahon cycle system shown in Fig. 2.28c is similar to
the Solvay unit with the difference that a displacer is used instead of the
expansion engine of the Solvay system. The Gifford-McMahon engine is
very like the piston/displacer Stirling engine shown in Fig. 2.28b but has
valves isolating the compressor from the regenerator and displacer unit.
Use of valves requires the provision of two connecting ducts rather than
the single connecting duct of a Stirling engine. Separation of the compressor
unit and the use of valves to control the flow gives the Gifford-McMahon
engine the same advantage of independent compressor I expander operation
and high-pressure ratios as the Solvay engine.
The Gifford-McMahon engine has an advantage over the Solvay
engine. It uses a displacer rather than a piston. A displacer requires only
a low-pressure seal. A piston requires the use of a high-pressure seal
whereas this apparently trivial distinction is the biggest single reason why
Gifford-McMahon engines are used very widely for cryogenic cooling
engines whereas not a single Solvay engine is in production. Air Products
and Chemicals, Inc. sometimes refer to their cryocoolers as modified Solvay
engines. They are indistinguishable from Gifford-McMahon engines and
do not include the low-temperature expansion engine characteristic of
Solvay machines.
Parenthetically, it may be noted that all regenerative cryocoolers
include a recuperative heat exchanger for cooling and sometimes to improve
heat transfer at the expansion space.
78 Chapter 2
EXPANDER
COMPRESSOR
Fig. 2.29. General arrangement for a Joule-Brayton cycle refrigerating engine with recupera-
tive heat exchanger.
Elementary Thermodynamics of Cryocoolers 79
v
COMPRESSOR T
Pmax
v s
EXPANDER
v. A "freezer" heat exchanger wherein the cold gas leaving the engine
is heated in the constant-pressure heating process (5-6). The heat
is abstracted from the surroundings of the freezer and is the useful
refrigeration of the cycle.
v1. The low-pressure fluid returns through the low-pressure flow
passages of the counterflow recuperative heat exchanger and is
heated in the constant pressure process (6-1).
It is possible for the expander and cooler to be coupled so the output work
of the expander helps to drive the compressor. This reduces the work
required by the compressor. However, in many systems it is not convenient
for the expansion engine and the compression engine to be coupled. They
may be physically separated. The machines may even be of different types.
It is not unusual to find a high-speed expansion turbine supplied with fluid
compressed in low-speed reciprocating compressors. Furthermore, the
expansion work may be negligible compared with a compression work. It
is simply dissipated in a friction brake or used to drive a small generating
set or pump.
JOULE-THOMSON EXPANSION
s s
(a) GASEOUS PROCESS {b) LIQUEFACTION PROCESS
JOULE-BRAYTON CYCLE
erative heat exchanger. The term Joule-Brayton cycle used here should
therefore be understood as the abbreviated form of "Joule-Brayton cycle
with recuperative constant pressure heat exchange."
When the cycle is used with a regenerative heat exchanger it is the
cycle we have earlier called the pseudo-Ericsson cycle (see Part II, Chapter
2).
SIEMENS CYCLE
The Siemens cycle is a reference cycle similar to the Joule-Brayton
cycle but with an isothermal compression (process 1-3 in Fig. 2.30). The
cycle was named after Sir Charles Siemens who was granted a patent in
1857 for a refrigerator working approximately on the reference cycle. The
isothermal work of compression is substantially less than the isentropic
work. Rotary machines are always assumed to operate with isentropic
processes but reciprocating compressors run slower than rotary systems
and have better means of cooling. Sometimes it is more appropriate to
assume their ideal compression process to be isothermal rather than isen-
tropic. We could also assume the expansion process to be isothermal
[process (4-6) in Fig. 2.30] rather than adiabatic [process (4-5) in Fig.
2.30]. This cycle with isothermal compression and expansion and constant
pressure recuperative heat exchange has not been named. With regenerative
heat exchange the cycle is of course the ideal Ericsson cycle.
LINDE-HAMPSON CYCLE
A cryocooler utilizing a Joule-Thomson expansion is shown in Fig.
2.33. This is known as the Hampson, the Linde-Hampson, the Linde, or
~~~ Compressor
Compressor
0 0 .T 0 Valves
J.T. Valve
O.T
Fig. 2.33. Hampson liquefaction cycle (after Radebaugh, ja) Hampson Cycle
1977). (Joule- Thomson Cycle)
84 Chapter 2
CLAUDE CYCLE
I
Q,T
Fig. 2.34. Claude cycle refrigeration system with no liquefaction of
the working fluid (after Radebaugh, 1977). (b) Claude Cycle
FLOW~ CD
+----;
INTERNAL
ENERGY
/
~
KINETIC/
ENERGY
1----+
WORK
!
Fig. 2.35. Steady flow thermodynamic black box.
86 Chapter 2
m(Hl-Hz)=m(C~-Ci)+ W -0 (2.73)
The change in enthalpy from state 1 to state 2 is therefore the sum of the
change in kinetic energy, the work done, and the heat transferred.
Elementary analysis of recuperative cryocoolers can be approached
by considering enthalpy values at the different locations of the system.
Consider the Joule-Brayton cycle system shown in Fig. 2.30. The change
in enthalpy per unit mass flow (Hz- H1) during compression is the work
of compression that must be supplied to the cycle. This assumes the velocity
of the gas at 2 is about the same as at 1. The heat transfer 0 during
compression is zero because the process is assumed to be isentropic (i.e.,
A.S = A.O/T = 0).
Similarly, the heat transfer 0<3 -z> from the high-pressure fluid in the
water-cooled heat exchanger is (Hz- H3) and in the recuperative heat
exchanger 0<3-4) = (H3- H4). The work produced in the expander W4-s is
(H4 - Hs) and may or may not be applied to partially offset the work of
compression. Finally the useful refrigeration of the cycle is 0<6 -s> =
(H6 - H 5 ) and the heat transfer in the recuperative exchanger is 0<1-6) =
(Hl-H6).
The net work input to the cycle, Wno (assuming the expander work
output is used to reduce the compression work) is
(2.74)
(2.75)
Elementary Thermodynamics of Cryocoolers 87
(2.76)
Compressors
Figure 2.36 illustrates on the P- V and T-S planes, the different
processes for compression.
1. Process (1-2) is an isothermal compression represented by PV =
const.
ii. Process (1-4) is isentropic represented by p V, = const. y is the
ratio of specific heats of the gas, y = CP/ Cv.
m. Process (1-3) is polytropic between an isothermal and an isentropic
process. The compression can be represented by p V" = const,
where (1 :5 n :5 y ).
IV. Process (1-5) is polytropic represented by the equation pVm =
const, where m ~ y.
The ideal reference cycle for reciprocating compressors is the cycle with
isothermal compression. For rotary systems the isentropic compression
4
P maa
PV=C
I.
Pv"=c
T
p
I I
pvY.c
I,( I
P mon 7 PV •C
v s
Fig. 2.36. Compression processes on the P-V and T-S planes.
88 Chapter 2
cycle is the ideal reference cycle. In a rotary machine at high speed the
gas particles are moving so fast that there is simply no time for cooling to
take place. In the reciprocating compressor the gas passes through the
machine at a slower pace and so some cooling is possible. To achieve
isothermal compression would require either infinite rates of heat transfer
or an impracticably low speed of operation. In practice reciprocating
compressors follow the intermediate polytropic path, process (1-3) on Fig.
2.36. Rotary compressors follow the polytropic process (1-5).
The isothermal efficiency for a compressor (usually reciprocating) is
defined as
isothermal compression work area 1-2-6-7 Hz-Hz
T/(iso) =
actua l compression
.
war k area 1-3-6-7 H3-Hl
(2.77)
The isentropic efficiency for a compressor (usually rotary but sometimes
also used for reciprocating) is defined as
This usually ranges from 60% to 80% for reciprocating engines and up to
90% for rotary systems.
Reciprocating compressors have crankshafts, with main bearings, con-
necting rods with big and small end bearings, pistons with piston guide
rings and sealing rings. They may also be equipped with water cooling
requiring some pump work, with forced lubrication requiring oil pump
work and a fan for air cooling. The mechanical and hydraulic losses tend
to be appreciable compared with the pneumatic effects arising from the
relatively slow moving gas being compressed.
Elementary Thermodynamics of Cryocoolers 89
Rotary systems have only main shaft bearings and rotary shaft seals
which consume little power. The principal losses are pneumatic arising
from friction and windage effects of the high-speed rotation of the disks,
drums, and blades within the compression space.
Frequently the isentropic efficiency is quoted for reciprocating com-
pressors rather than isothermal efficiency. This was adopted by compressor
manufacturers to put their machines in a better light since the ratio of the
isentropic work to the actual shaft input is larger than the ratio of isothermal
work to actual shaft input. Both are arbitrary ratios and the actual work
required remains the same regardless of the efficiency used. When discussing
compressor efficiencies it is important to establish if the efficiency being
quoted is isentropic, isothermal or mechanical.
Yet another efficiency is much used, the volumetric efficiency, defined
as
volume of gas actually pumped
Tl(vol) = d'1sp1acement vo1ume of compressor (2.80)
The volume of gas passing through the compressor is always less than the
piston displacement volume. The reduction is due to a number of causes
illustrated by the work diagram shown in Fig. 2.37. This shows the ideal
work diagram 1-2-3-4 superimposed on an actual compressor indicator
diagram. The effects have been exaggerated for the sake of clarity. Due
to pressure losses in the intake duct (perhaps including a filter) the pressure
at the start of compression (5) is less than the theoretical value (1).
The mass of gas in the cylinder can be calculated as
PsVs
me= RTs (2.81)
where me is the mass of gas in the cylinder, ps, Vs, and T 5 are the pressure,
volume, and temperature at point (5), and R is the characteristic gas
7 6
Multistage Compression
When gases are compressed over a large pressure range it is advan-
tageous to perform the compression in two or more stages with intercooling
between stages and aftercooling the final stage. In Fig. 2.38 the P- V and
T-S diagrams for two-stage compression are superimposed on those for a
single stage.
The single-stage compression diagram is represented by process 1-11-
6-10. The two-stage compression process is represented by the low-
pressure stage, process 1-2-8-9, and the high-pressure stage, process
3-4-6-7. At the intermediate pressure the fluid is assumed to be cooled
to the original temperature. Process 1-3-5 is the isotherm for the tem-
perature at state 1.
Elementary Thermodynamics of Cryocoolers 91
12
I
I
I
I
I
I
I
\ I
T \ /~13
\2/ l
I
I
I
II Pmin
I
v
Fig. 2.38. Multistage compression.
Expanders
Processes for expansion are shown in Fig. 2.39. From an initial high
pressure the fluid may expand:
i. Isothermally, process (1-2) represented by pV =canst.
92 Chapter 2
LATE CUTOFF
Pmax
I-------,
""""6
PV•C ',
v s
Fig. 2.39. Expansion processes on the P-V and T-S diagrams.
with isothermal expansion. In situations where the heat leak from surround-
ings dominates the internal heat generation (usually the case with small
electronic or superconducting devices) it is advantageous to absorb heat
over a range of temperatures. In these situations a better comparison might
be to increase the temperature at the start of the adiabatic expansion to
generate the same quantity of refrigeration (at variable temperature) with
adiabatic expansion as with isothermal expansion. This modified adiabatic
expansion is represented by the process (8-9). Refrigeration is available
at variable temperatures in the process (9-10). The amount of refrigeration
is area 9-10-12-11 and has the same area as the isothermal process
3-5-7-6. With refrigeration at variable temperatures the duty and hence
the size of the recuperative heat exchanger is decreased since the tem-
perature is only reduced to (8) rather than (3).
REFERENCES
Stirling Cryocoolers
with double expansion was also developed which can achieve temperatures
of 20 K and has been used for research and hydrogen, neon, and helium
liquefaction.
Derivatives of the early engine have been incorporated in many cry-
ogenic systems involving cooling or gas liquefaction and reliquefaction.
These include standard cabinets with cold air recirculation for metallurgical
processing, freeze drying, food processing, and a whole range of research
applications requiring refrigeration at precisely controlled temperatures in
the range from ambient to 20 K. It is used for reliquefaction of "flash" gas
in liquid oxygen or other cryogenic reservoirs.
In addition to the basic Stirling cryocooler, much other work has been
carried out in miniature engines and in the intermediate size coolers shown
in Figs. 3.3 and 3.4. No significant application for the intermediate size
Cryoc oolers
Stirlil!g 99
100 Chapter 3
cooler has yet developed but will come with the growth of LNG transporta-
tion and storage systems and the eventual introduction of superconducting
systems.
Military interest in infrared thermal imaging equipment for night vision
and heat-seeking missile guidance focused research attention on miniature
cryocoolers in the 1960s. This is now sustained by continued military
interest, by spacecraft instrument applications, and by incipient applications
of superconducting electronic and electrical systems on a broad front .
The principal military interest in miniature cryocoolers was manifest
in the United States. Policies of the U.S. Government led Philips to
concentrate their development and production of miniature, militarily
supported cryocoolers at the associate company North American Philips.
The research facilities of the company are located in Briarcliff Manor, New
York and the production facilities were at Ashton, Rhode Island. In the
early 1970s, the production facilities were sold to a competitor, Cryogenic
Technology Incorporated (CTI), of Waltham, Massachusetts. CTI also
distributes in the United States the larger Philips cryogenic equipment
made in Eindhoven.
The research laboratory continues at Briarcliff Manor principally
engaged in research and prototype development of Stirling cryocoolers and
Fig. 3.5 . Miniature Stirling cryocooler with rhombic drive (after Lindale, 1978).
102 Chapter 3
Fig. 3.6. Cross section of miniature Stirling cryocooler with rhombic drive (after Lindale, 1978).
Stirling Cryocoolers 103
Fig. 3. 7. Ultracompact free-piston Stirling cryocooler with linear motor drive (Philips Research
Labs).
CLASSIFICATION
The use of valves for flow regulation has the advantage of increased
flexibility in flow control and timing with the possibility of virtually unre-
stricted pressure ratios (Pmax/ Pmin). Valves add to the mechanical com-
plexity of the system and provide sources of noise and additional points
of wear so the prospect for long life with low maintenance is reduced.
Stirling cryocoolers, machines operating on a regenerative cycle
without valves, can be classified into two principal groups, namely, machines
with mechanically induced compression and machines with thermally
induced compression. In this review, machines with thermal compression
are treated separately as Vuilleumier Engines in Chapter 4.
Machines with mechanical compressors may be further subdivided into
two families identified as single-acting and double-acting engines.
A- PISTON
9 - OISPLACER
C- EXPANSION SPACE
0-COMPRESSION SPACE
E- REGENERATOR
F- FREEZER
G-COOLER
other operating as free piston (or displacer). This has been called the
hybrid crank controlled (or disciplined) piston/free displacer engine. It is
widely used for split-Stirling coolers and can also be used for integral-
Stirling machines (see Walker, 1973).
A distinction also has to be made between machines with a separate
stationary regenerator and those where the regenerator is incorporated in
the displacer and moves with it.
Design variants for integral-Stirling engines are shown in Fig. 3.10
and for split-Stirling engines in Fig. 3.11. Design variants for two-piston
machines are shown in Fig. 3.12. The originator of a particular configuration
is given on these figures where it is known. In most cases, the system was
conceived for use as a prime mover and many of these cases have never
been made as cooling engines. This does not preclude future development
of new and different arrangements.
The spectrum of possibilities for two-piston arrangements is very wide.
It includes the conventional piston/cylinder systems shown in Fig. 3.12
and many others employing metal bellows and diaphragms, rotary assem-
blies and free-piston devices where motion is induced by fluidic forces
alone, or in combination with mechanical springs. The only two-piston
machine available commercially (from Philips) is the intermediate capacity
cooler shown in Fig. 3.4. This interesting machine utilizes a Joy crompressor
Stirling Cryocoolers 107
~ ~
KINEMATIC DRIVE
a) RHOMBIC
b) SCOTCH YOKE
c) CRANK/
CONNECTING ROO
{BEALE)
(WALKER)
crankcase and pistons to generate the hydraulic drive for the two opposed
pistons of the upper, central Stirling engine cooler. It is a difficult machine
to classify. The Stirling system pistons are activated by fluidic forces only,
pneumatic on one side, hydraulic on the other, to make it a free piston
SINGLE-ACTING 1
SPLIT STIRLING
r ENGINE
FREE PISTON
I
HYBRID CRANK DISCIPLINED PISTON
AND DISPLACER CONTROLLED PISTON AND DISPLACER
Yg
/FREE DISPLACER
KINEMATIC DRIVE
~ (HIGA)
I
PARALLEL-CYLINDER
T
PARALLEL-CYLINDER
I
'\ i
TWIN-EXPANSION CYLINDERS AT RIGHT
REGENERATIVE EXTERNAL CYUNDER ANGLES REGENERATIVE
DISPLACE A REGENERATOR
F
DISPLACER
~
(LAUBEREAU
SCHWARTZKOPH) (HEINRIC!)
{RAINBOW)
(ROBINSON)
PARALLEL CYLINDER
,-
~ (RIDER)
V CYLINDER
~
(FINKELSTEIN)
r--
OPPOSED PISTON
~
L___
engine. Yet, it has a kinematic vee cylinder crank connecting rod mechanism
through which the input work is supplied to the engine to generate the
hydraulic forces inducing the Stirling piston motion.
Regenerative Stirling engines using bellows are not uncommon in the
power systems field. The artificial heart engines of the McDonnel
Douglas/Richland Energy Laboratory and the Thermo-Electron Corpor-
ation incorporate a variety of bellows and diaphragms. The radio-
isotope/propane-fueled power generator developed by Cooke-Yarborough
at the U.K. Harwell Atomic Energy Establishment uses a diaphragm
(Walker, 1980).
act1ng I
I Double-
p1ston engmes J
I
Two- cylinder tw1n
system Compress1on
1n one cylmder and
expans1on 1n the
other for both
systems
(Franchot)
pressure and a large temperature difference between the top and bottom
faces.
The great advantage of double-acting engines is that the number of
reciprocating elements is half the number required in multiple arrange-
ments of single-acting engines. This can lead to major simplifications of
the kinematic drive arrangements and to reduced costs.
The principal disadvantage is the limited flexibility in design, and to
a lesser extent, operation. Moreover, for prototype development, it is
necessary to proceed with the entire multiple cylinder engine rather than
an experimental single-cylinder unit which can then be reproduced in
multiple cylinder versions. Nevertheless, the elimination of half the moving
parts is so significant an economy, that the major interest for prime mover
Stirling engines greater than 15 kW is concentrated on double-acting
arrangements.
Double-acting engines were invented in the 19th Century. Babcock
(1885) attributes the invention of the two-cylinder, twin system to the
Frenchman, Charles Louis Franchot in 1853. It was reinvented a century
later by Finkelstein (1959) for the prototype unit of a domestic refrigerator
and was later incorporated by him in the concept of a multiple-cylinder
free-piston unit for power generation.
The general form of double-acting engines with three or more cylinders
is ascribed by Babcock to Sir William Siemens who, in 1863, conceived
.....
.....
Q
Fig. 3.14. Four-cylinder double-acting Stirling engine invented by Sir William Siemens in 1863 (after Babcock, 1885).
(")
='
~
..;;-
~
Stirling Cryocoolers 111
the four-cylinder engine shown in Fig. 3.14. The engine had four cylinders
each containing a single piston and which were coupled in a square arrange-
ment through the associated heat exchangers and regenerators. It was
activated by a wobble-plate drive mechanism.
So far as is known, the Siemens engine was not taken beyond the
conceptual stage until it was reinvented, complete with wobble-plate, 80
years later by van Weenan in course of the Philips program. Following
some early success at Philips, the Siemens engine was abandoned in favor
of the single-acting "rhombic drive" engine. Later, in 1965, one of the
Philips licensees, General Motors, in the United States, resumed work on
double-acting engines and by 1970, Philips and the other licensees had
concentrated their efforts on double-acting engines for automotive engine
applications. The Philips unit, shown in Fig. 3.15, is a four-cylinder version
with a swash-plate drive* and exists in prototype 150- and 75-kW forms.
Until recently (1978), the engine was under development for automotive
application by Philips and the U.S. Ford Motor Co. in a 7-year $180 million
program funded by the U.S. Department of Energy. Ford abandoned its
development effort on Stirling engines in 1978.
An alternative kinematic mechanism has been adopted by the Philips
licensee United Stirling of Sweden using twin cranks and connecting rods
as shown in Fig. 3 .16. Another version with the cylinders arranged in a
shallow vee form was earlier favored. This used a single crankshaft but the
inclination of the cylinders increases both production and maintenance
costs. The parallel cylinder version is the unit intended for preproduction
evaluation in 40-, 75-, and 150-kW versions for automotive use, in mining
applications, for boats and stationary power generation (see Walker, 1980).
The United Stirling engine is the power unit for development of a Stirling
automotive unit in a program sponsored by the U.S. Department of Energy.
Mechanical Technology Inc. of Latham, New York and American Motors
are participating in the program.
Double-acting engines have not been used for cryogenic refrigerating
purposes. This is due entirely to the absence of need for cryogenic coolers
of intermediate or large capacity. Significant change is now at hand with
the development of relatively large-scale superconducting systems and
liquid natural gas transportation and reliquefaction systems. These develop-
ments will require new cryocooler concepts. Double-acting Stirling cryo-
coolers of intermediate and large-scale capacity will likely be prominent.
Already one major U.S. corporation is known to be actively developing a
Fig. 3.15 . Philips double-acting Stirling engine with swash-plate drive (after Walker, 1979).
Fig. 3.16. United Stirling double-acting Stirling engine with crank-connecting rod drive.
vehicle. The design was compact with a mass and volume half the conven-
tional vapor compression unit. Later on, other designs were evaluated for
railcar air conditioning systems and for cockpit cooling of fighter aircraft
and crop-spraying aircraft. In all cases, promising designs emerged but did
not proceed to prototype development.
The growth of interest and capability in single-acting free-piston
Stirling engines in the 1970s resulted in the development of concepts for
multiple-cylinder double-acting free-piston Stirling engine coolers and
prime movers. One engineering prototype development along the lines
proposed by Finkelstein (1978) is under way and further hardware develop-
ment is expected in the future.
One version of a three-cylinder free-piston Stirling cooling engine of
intermediate capacity is shown in Fig. 3.17. Three cylinders are used to
take advantage of the 120° phase difference in conventional3-phase electric
power supplies to energize the linear motor drive. Ingenious electrical
design is required to devise a motor capable of accepting high-frequency
(50 or 60 Hz) conventional power supplies and operating with a frequency
half to one third this value.
SPACE
"Ol'Ofl ORIVE
CONTAINED WITHIN
RECIPROCATOR
Fig. 3.17. D esign concept for three-cylinder free-pisto n double-acting Stirling cryocooler.
Stirling Cryocoolers 115
UHEAR G.AS
BEARING
SPRING
CAVITY
SPRING CAVITY
UNEAtt MOTOR
(b) (c)
(o)
Fig. 3.18. Alternative design features to accommodate the difference in swept volumes in the
expansion and compression spaces of double-acting Stirling cryocoolers.
116 Chapter 3
COMPR£SSION
S..ACE
COMPRESSION SPACE
n
II
II
II
- - - -----II
----~ IL~
I:Jl]·
L_
L___J
Fig. 3.20. Multielement cryogenic spot RAINBOW I CHELLIS (1970)
cooling.
Cylinder Head
Oispl a cer
Regenerator
Piston
(a )
(I)
Compte'SsiOn Regenero1 1'-'e E~~:panston ReQeneroltve
Cooli ng Heottng
(b )
Fig. 3.22. (a) Piston and displacer disposition, and (b) displacement-time diagrams for the
integral Stirling engine operating on the ideal Stirling cycle.
Fig. 3.23 . Piston and displacer disposition and the displacement-time diagrams for the
split-Stirling engine operating on the ideal Stirling cycle.
B c 0 A
~ (ol (b)
(I)
(c) (d)
A
~
EXPANSION
(2)
B
v
~
COMPRESSION (3)
c
v
with those given on Fig. 2.9, except for system A. Therefore, the four P- V
diagrams for the double-acting engine are also included in Fig. 3.24.
At station 1 the reciprocator in cylinder (a) is at the top dead center
position. The reciprocator in cylinder (b) is at the midstroke and descending.
In cylinder (c) the reciprocator is at the bottom dead center position and
in cylinder (d) the reciprocator is at midstroke and ascending. There is,
therefore, a phase difference of 90° between the reciprocator positions in
adjacent cylinders. The phase difference depends on the number of cylin-
ders, i.e., a= 360/n, where a is the phase difference and n is the number
of cylinders. Here n can be 3, 4, 5, or 6 for satisfactory operation.
At station 2 all the pistons have moved by one half stroke but the
phase relationship has been preserved. In cylinder (a) the reciprocator has
descended to midstroke, in cylinder (b) it has descended to the bottom
dead center position, in (c) it has ascended to the midstroke position, and
in (d) it has ascended to the top dead center position.
Consider what has happened in this process (1-2). System A, made
up of ( V c )a and ( V E )d, has decreased in volume to accomplish the isothermal
compression 1-2 shown on the P-V diagram for system A. Simultaneously
system B made up of (Vc)b and (VE)a maintained the total system volume
constant but pushed the working fluid from (Vc)b to (VE)a, that is the
constant volume regenerative cooling process (1-2) shown in the P- V
diagram for system B. While this was going on, system C was experiencing
an isothermal expansion and system D was experiencing a constant volume
regenerative heating process. If the pistons are now allowed to progress
another half stroke to the disposition shown in station 3, the analysis of
events will show the process (2-3), shown on the four P- V diagrams, will
occur. Similarly, we can continue through the processes (3-4) and the final
process (4-1).
To facilitate understanding, the time-displacement diagrams for the
four Stirling systems have been drawn in Fig. 3.25. These diagrams combine
the appropriate expansion and compression spaces constituting each system.
Thus, the upper diagram for system A, combines the motion of the
reciprocator in cylinder (d) (top half of the diagram) with the motion of
the reciprocator in cylinder (a) (bottom half of the diagram). The shaded
area in the top half represents the volume variation in the expansion space
of system A (located in cylinder (d) above the piston/displacer). The shaded
area in the lower half of the figure represents the volume variation of the
compression space of system A [located in cylinder (a) below the piston/dis-
placer]. At any instant, the sum of the two shaded areas represents the
total "live" working space in system A. The constant volume of the
connecting ducts containing the regenerator and associated heat exchangers
must be added to obtain the total working space. The same pattern of
Stirling Cryocoolers 123
So far it has been assumed the Stirling engine works on the ideal
Stirling cycle of two isothermal and two constant-volume regenerative
processes. Now it is time to consider the difficulties arising if attempts
are made to make cooling engines work on the ideal cycle. Having assessed
that problem, we can define a more realistic model with which to compare
the performance of actual engines.
124 Chapter 3
..,a:
=>
..,a:
(f)
(f)
a.
I
~ \""'
...
~
::l
...
II:
0.4 '---t--l
I
\
"'-...
r-
~
)
~
~
"" ~~
a. ozt---+---1
EXPANSION
I
IOMPRESSION TOTAL WORKING T • 2 , K = I 86
SPACE SPACE SPACE a;. 179. X•l
02 04 0 02 0.4 0.6 02 04 06 08 10
(V./VW) (Vc/VW) (V/Vwl
VOLUME
Fig. 3.27. Typical work diagrams for the expansion, compression, and total working spaces
of a Stirling cryocooler.
126 Chapter 3
Temperature-Entropy Diagrams
The introduction of continuous harmonic piston motion and finite void
volume in the heat exchangers, eliminates conditions inherent in the ideal
cycle that all the working fluid is instantaneously at the same state and
proceeds as a lumped mass around the cycle. With continuous harmonic
piston motion, the fluid is distributed about the machine in various spaces
at various temperatures.
Any particular element of mass inhabits a region of the engine and
when the engine is operating, the element will not move entirely from the
compression to expansion space, but will simply commute to and from,
say, the expansion space to the freezer. Another particle may inhabit the
regenerator and cooler and so on. It is not possible to draw meaningful
temperature-entropy diagrams for all the working fluid. Rather, individual
diagrams for particles at different regions of the engine can be constructed.
The individual envelopes drawn on the T-S plane can then be wholly
enclosed in boundaries that probably represent some kind of temperature-
entropy diagram.
Such representations have only limited use and are rarely drawn in
practice. Areas of these composite diagrams do not represent simple heat
transfer as in the idealized cycle.
Schmidt Cycle
The case discussed above with harmonic piston motion and isothermal
compression and expansion was analyzed by Gustav Schmidt (1871). It is
convenient to identify this case as the Schmidt cycle. The analysis is
presented in detail later in this chapter.
expansion spaces or the engine running at very low speeds. In real engines
running at realistic speed (15 to 50 Hz), conditions in the cylinders are
closer to adiabatic (no heat transfer) than isothermal (infinite heat transfer).
They may be close to isothermal in the freezer and cooler since these
components are designed to maximize heat transfer.
The departure from isothermal conditions in the compression and
expansion spaces causes a marked redistribution of working fluid in the
machine resulting in a decrease of refrigerating capacity and an increase
in the input work, thereby reducing the coefficient of performance. The
reduction may be as much as one half the ideal Carnot value.
Analysis of the Stirling engine with adiabatic processes is somewhat
more complicated than the isothermal Schmidt cycle. The theory has been
developed, and will be discussed later. For reference purposes, the adiabatic
case is called the Finkelstein adiabatic cycle.
VOLUME V8 VOLUME Vc
Fig. 3.28. Regenerator pressure drop and the effect on EXPANSION COMPRESSION
expansion cylinder work. SPACE SPACE
128 Chapter 3
Now include the heat exchanger thermal potentials of, say, 50 K in the
cooler and 20 K in the freezer. The thermodynamic temperature range of
operation is, therefore, increased to (290 +50) = 340 K and (80- 20) =
60 K. The coefficient of performance under these conditions is, therefore,
60
340-60 = 0 ' 214
Regenerator Contamination
The importance of the regenerator effectiveness to the ideal Stirling
cycle was demonstrated in Chapter 2 (see Fig. 2.21). Remarkably effective
regenerators (s > 0.95) are routinely used in Stirling cryogenic coolers with
matrices of finely divided materials in the form of wires or balls. This finely
divided material is also an excellent filter and effectively "strains out" any
contaminants contained in the working fluid passing through. The con-
taminants include lubricating oil leaking from the crank case, wear debris
from dry rubbing piston rings or seals, and condensed liquids or solids of
water and carbon dioxide from the working fluid. Blockage of regenerator
flow passages increases the fluid friction in the matrix (i.e., a larger pressure
drop) with all the unfortunate consequences to refrigerating capacity dis-
cussed above.
increases sharply at low temperatures. At the same time, the specific heat
of conventional matrix materials (copper and bronze wire, etc.) decreases
so that thermal saturation prevents effective operation of the matrix.
Change of the matrix material to lead produces an improvement. Lead
is not available in the form of very fine wires but can be obtained as fine
powder or "shot." Further improvement can be gained by the use of "rare
earth" materials (europium sulfide) with advantageous discontinuities in
the specific heat temperature characteristic. Eventually, though, thermal
saturation of the matrix prevails as the helium working fluid nears the
liquefaction temperature. The ultimate temperature achievable is rarely
less than 6 K. Zimmerman (see Chapter 1) has achieved 3.5 Kin a special
expansion machine with a "regenerative annulus," a pressure of less than
1 atmosphere of the helium working fluid and precooling with liquid helium
to 10 K.
Introduction
Theoretical analyses of Stirling engines acting as prime movers or
cooling engines have been developed with varying degrees of sophistication.
The most simple analysis is that for the ideal Stirling cycle, the thermo-
dynamic cycle, comprised of two isothermal and two constant-volume
regenerative processes. As we have seen, this involves such gross ideali-
zation as to be suitable only for the most elementary design calculations.
A more realistic analysis devised by Gustav Schmidt (1871) has become
the classical analysis of the cycle and is believed to give a reasonable
approximation of machine performance. The analysis is highly idealized,
and in practice, the indicated performance will likely be no better than
50% of the Schmidt cycle performance and often a good deal less. This is
because the Schmidt cycle assumes compression and expansion to be
isothermal. In practical engines running at 15Hz or more, the processes
are more nearly adiabatic. This apparently minor difference causes a redis-
tribution of the mass working fluid and has a critical effect on engine
performance.
Following introduction of the Schmidt analysis in 1871, nearly a
century was to elapse before Finkelstein (1960) devised the generalized
analysis of engines with processes of compression and expansion other than
isothermal. The Schmidt cycle with isothermal compression and expansion
then became a special case of the generalized Finkelstein analysis. Another
special case is the cycle with adiabatic processes in the engine cylinder (the
132 Chapter 3
Required Data.
1. Some reference temperature and pressure, or volume, say, condi-
tions, at state 1.
ii. Temperature ratio T = Tc/TE.
iii. Volume ratio r = V max/ V min·
For unit mass of working fluid, assumed to be a perfect gas, then V1 =
RTd Pt, from the characteristic gas equation.
The following data are required for each of the four cycle processes:
(a) Isothermal Compression Process (1-2). In this process, heat is
abstracted from the working fluid and rejected from the cycle at the
maximum cycle temperature Tc. Work is done on the working fluid equal
in magnitude to the heat rejected from the cycle. There is no change in
internal energy, and there is a decrease in entropy.
Heat transfer (Q) =work done (W) =P 1 V 1 ln (1/r) =RT1 ln (1/r). Change
in entropy (S 2 - S1 ) = R In (1/r ).
(b) Constant- Volume Regenerative Cooling Transfer Process (2-3).
In this process heat is transferred to the regenerative matrix from the
working fluid, decreasing the temperature from Tc toTE. No work is done,
and there is a decrease in the entropy, and internal energy, of the working
fluid.
= 1/('T -1)
This value corresponds to the Carnot value between the same temperature
limits.
If the instantaneous pressure is the same throughout the system, and equal
top, say, and if Te and Tc are constant at TE and Tc then, substituting for
the volumes, eliminating R, and rearranging
K/p = (Tc/TE)(l +cos c/J) +K[l +cos (c/J -a)]+ (2 V 0 Tc)/(VETv) (3.6)
If the temperature variation in the dead space is linear in the axial direction,
then the mean temperature
Stirling Cryocoolers 137
then y = .Jr 2 cos (<P -(3), where tan (3 = z/x, z = r sin (3, and x = r cos (3
since
.Jr 2 cos (<P - (3) = .Jr 2 (cos <P cos (3 +sin <P sin (3)
= r cos <P cos (3 + r sin <P sin (3
= x cos <P + z sin <P
r21T
=(h) Jo {Pmax(1- 8)/[1 + 8 COS (</J- O)]} d(<{J- (}) (3.12)
Q =P= JpdV
then
r21T
Q = -! Jo {Pmean V£[1- ~cos (</J- 0)] sin </J} d<{J
r21T
= -!Pmean VE Jo [sin <P- ~(cos <P cos(} sin <P +sin(} sin 2 <P)] d<{J
<P] 21T
=- !Pmean VE [ -~sin(} 2 0
which conforms to the required equation [Eq. (3.14)], therefore the heat
transferred in the expansion space, from Eq. (3.16), is given by
Q = 7rPmean VE8 sin 8/[1 + (1- 8 2)112 ] (3.17)
There is no need to invoke the approximate expression for p in Eq.
(3.15). Instead, Eq. (3.10b) is a suitable starting point, and can be combined
with (3.14) to calculate the heat transfer per cycle, Q, as follows:*
!
Q = ,( dV = ,( -pmin( 1 + 8 ) . VB sin <P d<P
jP T1+8cos(<{>-8) 2
= -VEPmin(1+8)I
2
where
I _ ,( sin <P d<P
j 1 +8 cos(</> -8)
let <P- 8 = t/1, then t/f/2 = t: then d<f> x dt/1 = 2dt/(1 + t 2) and
2 . 2 . ,( dt
= 7rSm8- sm8j (1 + 8 )+( 1 _ 8 )t2
. 2 sin 8 I 1 1 2 12
=27rsm8-( 1 _ 82 )172 tan- {t/[(1+8)/(1-8)] 1 } o/J~o
11"
*I am indebted to Dr. A. J. Oxley of the United Kingdom Rutherford Laboratory for the
alternative derivation given here-G. W.
140 Chapter 3
CRANK ANGLE
Fig. 3.29. Cyclic temperature variation in the compression and expansion space of a Stirling
cryocooler operating on the Finkelstein adiabatic cycle (after Lee, 1976).
to note the temperature of fluid in the expansion space was less than the
freezer temperature for much of the time. The mean temperature in the
expansion space was 90 K, a difference of 10 K below the freezer tem-
perature, TE. In the compression space, the mean temperature of the fluid
coincided more or less with the cooler temperature of 300 K.
Significant work using the adiabatic cycle approach was reported by
Qvale and Smith (1969) and by Rios and Smith (1969). They consider a
basic cycle with adiabatic compression and expansion and then separately
assess the effects of irreversibilities. This approach allows for independent
study of individual heat exchange components in a series of successive
approximations that can be extended to the required degree of complexity.
Little has been published about theoretical work on Stirling engines
at the Philips Company or their licensees. Mention is made in the various
papers about the extensive use of computers and the ability to predict
engine performance to one or two percent of test bed performance. Odd
straws in the wind gathered over the years have led to the understanding
that a family of computer programs have been developed by Philips. The
ensemble of programs are thought to include some overall cycle simulation
programs for the complete engine at various levels of sophistication and
complexity, whereas others are available for the detail simulation and
design of single components or subassemblies in the engine. The basis of
the Philips thermodynamic analysis simulation program is thought to be
similar to that developed by Smith, basically the Finkelstein adiabatic cycle
with successive approximations and corrections. The adherence by the
Philips Company over many years to similar design configurations has
Stirling Cryocoolers 145
Nodal Analysis
Nodal analysis of Stirling engines was pioneered by Finkelstein (1975)
but other nodal analysis programs were prepared independently. Develop-
ment was principally directed to prime movers, although the programs are
equally applicable to cooling engines. In nodal analysis programs, the
attempt was made to model the simultaneous energy and fluid flows
occurring in the engine and thus simulate exactly the engine cycle and
performance. This was achieved by writing and solving equations for the
conservation of mass, momentum, and energy for particular nodes, cells,
or elements of the engine. The equations were too complex for general
analytical solution and were solved numerically in terms of small incre-
mental time steps. The equations were invariably reduced to one-
dimensional form with additional simplifications. A comparative discussion
of the various programs has been given by Urieli (1979).
All the nodal analysis programs are basically similar in their general
approach. The design of engine to be simulated must be known in exact
detail to the extent that the mechanical arrangement, cylinder wall thickness
and material, heat exchanger tube diameter, fin dimensions or matrix pore
sizes, are all specified. This design is then broken down, as experience
dictates, into a number of nodes, elements, or control volumes. Some
operating conditions must be specified such as charge, pressure, and tem-
perature of the energy source and sink.
Differential equations for the conservation of mass, momentum, and
energy must be developed and generally converted to difference equations.
Empirical formulas for the aerodynamic friction and heat transfer effects
must be included as well as an equation of state for the working fluid. A
mathematically stable method must then be found for numerical solution
of the difference equations to resolve the pressure, temperature, and mass
distribution in the engine at the end of a particular time step, given the
conditions at the beginning of the step.
The usual procedure for solution is to start with some initial arbitrary
assumed conditions and then proceed through several engine cycles until
quasisteady state is achieved, when the instantaneous cyclic values of
pressure, temperature, and mass distribution are not significantly different
from the preceding cycle. The cyclic pressures and volumes are integrated
to calculate work transfer, heat flows are estimated, and the overall thermal
146 Chapter 3
w
a:
,_
:>
"a:~
::;
7: ~
c --
Fig. 3.30. Three-dimensional representation of the gas temperature versus node position and
cyclic angle in a Stirling engine prime mover (after Schock, 1978).
Stirlin2 Cryocoolers 147
Fig. 3.31. Three-dimensional representation of the gas temperatures versus node volumes
and cyclic angle in a Stirling engine prime mover (after Schock, 1978).
T=w 1
C = 1/K 1 p,
·~""7r-;::;:-
;; 7~ p,/v, t~~~ ,,..~.....,j/1
...--
~0
T=w 1
10
Crankcase
Fig. 3.33. Two-dimensional representation of Stirling cryocooler used for nodal analysis (after
Finkelstein eta/., 1971).
150 Chapter 3
water applications. Present status and potential applications for the com-
puter program are not known.
Dr. Allen Organ, formerly of Kings College, University of London,
now at the School of Engineering, University of Cambridge, has recently
published papers indicative of a powerful nodal analysis program in forma-
tion, but not yet fully developed. Readers are advised to contact Dr. Organ
directly for up-to-date information.
Several papers have hinted at the existence and use of nodal analysis
type programs at Philips and their licensees. No details have been given.
It is understood that for general design work, adiabatic cycle programs are
preferred over the more sophisticated nodal analysis programs.
Summary
Stirling engine simulation by nodal analysis is an expensive, time-
consuming activity to be reserved for those professionally engaged in
Stirling engine development or academics engaged in the training of
engineers at the graduate level.
At present, no comprehensive evaluation is available for the various
Stirling engine simulation programs described above. Martini (1978)
attempted the preparation of a design manual including a comparison of
the various theoretical procedures, but concentrated his efforts principally
on isothermal type analyses with only passing reference to the more sophisti-
cated nodal analyses. He later extended this work to higher-level simulation
programs (1980).
The program of Urieli and Finegold-Vanderbrug are completely listed
in the references cited above. However, neither of these programs has been
extensively validated by experimental data. The programs by Schock and
by Tew are neither listed nor have been validated experimentally. Neither
program is available commercially. Finkelstein's program is available com-
mercially on payment of moderate royalties, but no documentation of
validity with reference to practical engines is available. The Sunpower
program is available commercially and has been validated extensively
during the development of Beale-type free-piston engines. The applicability
of the program to engines with crank mechanisms (disciplined-piston
engines) has also been established recently.
For those wishing to obtain a simulation or optimization of a Stirlihg
engine concept, the choice is between Sunpower and Finkelstein. Another
possiblity is that Philips could offer consultancy services to evaluate new
concepts using their undoubted expertise and experience in the field.
It is simply out of the question to "dabble" in nodal analysis even
using the existing programs. The few companies and individuals wishing
152 Chapter 3
PRACTICAL DESIGN
Design Parameters
The principal independently chosen design parameters for Stirling
cooling engines are as follows:
1. The temperature ratio T = Tc/TE, the ratio of temperatures in the
compression and expansion spaces.
2. The swept volume ratio K = Vc!VE, the ratio of swept volumes
in the compression and expansion spaces.
3. The dead volume ratio X= Vv!VE, total internal volume of heat
exchangers and associated ducts expressed as a fraction of the
swept volume in the expansion space.
4. The phase angle a by which volume variations in the expansion
space lead those in the compression space.
Stirling Cryocoolers 153
0~1> 03 \
~£
g
w \WORK DON~
~ 02 \ --1-
lS
""'
0 \
HE~ LIFTED
-
v- ~
"" ~
0 ~I€£ 0
0 It---
........ --
/
v 100 K
EJ. I
200 K
EXPANSION SPACE TEMPERATURE
~ !'......
TE
300 K
Fig. 3.34. Effect of variation in expansion
space temperature TE on heat extracted,
input power and coefficient of performance
with constant values of K, X, and a.
Stirling Cryocoolers 155
~I,;
I 0
:: ~ 0 2 ,-----,-----,--...-----,
'"z
g
"'
<r
0
;!:
:5 011---fr----t--1-----j
~03
0]~e
-
1
0
!:~
'"0z
WO'K DONE
............
0
K
~ 02
;!:
"'
0 ~
--
..............
HEAT LIFTED
,"""-
~
r---
10
DEAD SPACE RATIO V0 / VE
20
156 Chapter 3
of the swept volume ratio Vel VE. Change in the temperature ratio 'T
produces a change in the refrigeration capacity and in the optimum swept
volume ratio. Increase in the value of T (decrease in the expansion space
temperature) produce a decrease in the heat extracted and an increase in
the optimum swept volume ratio.
Figure 3.36 shows the effect of variation in the dead space ratio, X.
Increase in the dead space causes a decrease in the heat lifted so the dead
space volume must be minimized for high refrigerating effects. In practical
engines, the situation is more complicated with increased dead space. A
positive improvement in refrigeration performance can result, provided
the increased dead space is effectively utilized as increased void volume
in the heat exchangers. This results in increased thermal capacity or bigger
flow passages with reduced fluid friction losses. It is vital that the dead
volume be limited to effective space and not wasted in oversized ducts,
large ports, and excessive clearance space.
Theoretical deterioration of refrigerating capacity with increase in
dead volume can be understood by considering the effect on the volume
compression ratio ( V max! V min). Increases in dead space decreases the
volume compression ratio and the pressure compression ratio (Pmax!Pmin).
Reduction in the pressure ratio results in decreased area of the expansion
space work diagram and a smaller refrigerating capacity.
Figure 3.37 shows the effect of variation in the phase angle a. The
refrigerating capacity is remarkably insensitive to changes in the phase
angle from 60° to 120° of crank rotation, but does exhibit an optimum
value around 90° at which the heat lifted is a maximum.
For any given combination of temperature ratio, 'T and dead space,
X, three-dimensional surfaces may be drawn as shown in Fig. 3.38. This
illustrates the heat lifted (OE! Pmax V T) as a function of swept volume ratio
K and the phase angle a. Different combinations of 'T and X result in the
Stirling Cryocoolers 157
generation of different surfaces, sitting one upon the other like onion skins.
The apex of the surface in Fig. 3.38 is the maximum refrigeration capacity
achievable with the given combination of r and X. The values of K and a
corresponding the apex of the surface are the optimum values necessary
to achieve the maximum refrigeration.
I #. v
#. '/ /
0·12
/-~ 0 / v
A~ ~.....-: Y!
8
~~~
v I
~ ....- ~
,/ P"
I
1\ i\ I
0·04
r-"
1-----
X""0•1,0·25, 0·5, 1·0, 1·5, 2·0
I I I I I
I \I
7
1/
II
6
5
I I I
I
0-.
" r--..
!'.:: .......
"~ 1'-.
I')._
~
~
:-
1 I T '::::::::: ::§
X"'0·1, 0·25, 0-5, 1·0, 1·5, 2·0
40 50 60 75 100 150 00 30
Fig. 3.39. Consolidated design chart for Stirl-
30
Expansion space temperatura (K) ing cryocoolers based on Schmidt theory.
To use the chart, one simply determines the expansion space tem-
perature at which refrigeration is required. Through this point, a vertical
line is drawn crossing the three diagrams. For different values of the dead
space, X, one can read off values of the maximum refrigeration capacity,
the optimum phase angle and swept volume ratio to achieve capacity. The
actual refrigerating capacity will be 30% to 50% of the indicated capacity.
Multiple-Expansion
The possiblity of an engine with multiple-expansion spaces coupled
to a single compressor was discussed above. The expansion process in all
spaces occurs simultaneously and at the same pressure. In electrical
engineering terms, this is "parallel expansion." There also exists the possi-
bility of "series expansion" at different temperature levels in a double
(triple or quadruple) expansion engine.
The principle of the double-expansion engine is illustrated in Fig. 3.40.
This shows the cross section of an integral Stirling single-acting piston-
displacer engine. The regenerator is contained within the displacer and
moves with it. The displacer-regenerator is arranged in two sections having
different diameters, with the top section substantially smaller than the lower
section. The stepped displacer operates in a cylinder which has correspond-
ing changes in diameter to fit the displacer.
Two expansion spaces are thereby created, namely, the customary
expansion space above the displacer, and an intermediate or "extra"
expansion space as shown in Fig. 3.40. The volume variations of the two
expansion spaces are in phase with each other. The two spaces are mutually
coupled through the porous regenerators and with the compression space
162 Chapter 3
20 K
60 K
300 K
so the pressure in all spaces is the same, apart from slight differences due
to regenerator pressure drop.
When the machine operates, it produces refrigeration at two tem-
perature levels, shown in Fig. 3.40 as 60 K in the intermediate expansion
space and 20 Kin the upper expansion space. The generation of intermedi-
ate temperature refrigeration creates a "thermal buffer" to absorb heat
leaks from ambient temperature sources and thereby protect the refriger-
ation generated at the lower temperature. The heat leak to the upper
expansion space is now much reduced for the temperature potential
between the buffer and low-temperature expansion space is only 60-20 K.
The drain of low-temperature refrigeration is correspondingly reduced.
Space limitations preclude a full discussion here of thermodynamic aspects
of multiple-expansion engines, but the same principle of cascade refriger-
ation at progressively lower temperatures is widely used in large-scale
recuperative gaseous liquefaction systems.
The pioneer work in double-expansion Stirling engines was carried
out at Philips in the early 1960s (Prast, 1963) using a modified air liquefier
of the type shown in fig. 3.1 with additional expansion space. With this
machine, temperatures of 20 K were readily achieved. It has been incorpor-
ated in a variety of hydrogen and helium liquefiers and other low-
temperature refrigeration systems.
Later on, Daniels and du Pre (1971) described a miniature triple-
expansion Stirling engine which attained a temperature of 7.8 K. A
diagrammatic cross section of the engine is shown in Fig. 3.41 and the
multi-element displacer in Fig. 3.42. Figure 3.43 is a view of the prototype
unit with the actual cylinder head. The tangle of tubes and wires around
the cylinder head are part of the apparatus for measuring the expansion
Stirling Cryocoolers 163
Free-Piston
This already long chapter could be greatly extended by proper con-
sideration of free-piston Stirling engines. In these machines, the pistons
Stirling Cryocoolers 165
Fig. 3.43. Laboratory prototype of miniature triple-expansion Stirling engine (after Daniels
and du Pre, 1971 ).
and displacers are not coupled through a kinematic chain (crank and
connecting rod) but are activated solely by fluidic forces acting on them.
Free-piston Stirling engines are in rapid development as prime movers
and, to a lesser extent, as cryocoolers (Haarhuis, 1978; de Jonge, 1979).
A comprehensive discourse on Stirling free-piston engines is contained in
166 Chapter 3
300
w
a: 160
.....,::>
a:
w
140 \' \
\
0. \
::.
~ 120
\
•oo I \
80
'
\ \ / - lsi STAGE
'
60
"---~---,
40
''
- - - 2nd STAGE
>O
- 3rd STAGE Fig. 3.44. Cooldown characteristics of
10 20 30 40 !)Q 60 70
a miniature triple-expansion Stirling
TIME (MIN) FROM START engine (after Daniels and du Pre, 1971).
:,,
tVe l ,..,. IVe)mo• tVelmln
I
Pmln IP~,
(V(lmu IVe:Jmln
.
IVtlmln
6 •
The mass ratio of piston to displacer and the area ratio of displacer
rod area to cylinder cross section area are critical to successful operation.
When these are in proper proportion, the displacer will be found to oscillate
in the cylinder as shown in Fig. 3.45, at the same frequency as the piston
but with a phase difference about 80° in advance of the piston. The motion
of the displacer in the cylinder causes the working fluid to move back and
forth between the expansion and compression spaces. The pressure is
decreasing when the displacer is closest to the piston. Thus, there is a
concentration of fluid in the expansion space and heat is abstracted from
the surroundings of the expansion space at low temperatures. Similarly,
the fluid is concentrated in the compression space when the pressure is
increasing; the displacer is widely separated from the piston and the piston
is rising on the compression stroke.
Figure 3.45 contains displacement-time curves for the piston and
displacer and pressure-time curves for the working space and bound space.
Typical pressure-volume diagrams are also included for the expansion,
compression, and total working spaces. On Fig. 3.45, the numbers 0 through
8 correspond to critical points of the cycle of operation in accordance with
the following key:
Fig. 3.46. Free displacer split-Stirling cryocooler, Model No. CM-3 80 K manufactured by
Cryogenic Technology Inc., Waltham, Massachusetts.
Assume first that the displacer has a low mass (therefore a low inertia)
and is completely unrestrained to motion in the cylinder within the confines
of the permissible stroke. The net fluidic force acting on the displacer will
cause the displacer to shuttle back and forth in the cylinder as shown in
the P-V plane, Fig. 3.47b, for both the compression and expansion spaces
in the cylinder. If the displacer had no mass, the motion would be the
straight line A-B-C. However, the displacer must have some mass and so
will follow the curved path A-E-B-F-C. It will return along the same path
and consequently there will be no refrigeration effect, for the area of the
work diagram on the return stroke will be exactly the same as on the
outward stroke. The net area of cycle work will be zero.
Now let the displacer be partially restrained by the sliding seals at the
warm end of the displacer, separting the working space and bounce space.
Static friction will prevent the displacer moving until a sufficient fluidic
Stirling Cryocoolers 171
(b)
·~·b_ v
(c)
v
·~~
Fig. 3.47. Diagrammatic cross sec-
tion of split-Stirling expander ele-
ment, including the pressure-time
and work diagrams for the com pre- {a) EXPANDER ELEMENT
sion and expansion spaces of the (d )
expander element. EXPA"lSION SPACE
force is acting to overcome the friction. The displacer will then move, but
dynamic friction will require some fluidic forces to be acting on the displacer.
Therefore, the line of motion followed on the P-V plane on the outstroke
will be different from that followed on the instroke and a definite network
diagram will be created. Diagrams arising from both a light and heavy seal
restraint are shown by Fig. 3.47d. Increase in the restraint increases the
separation of the outstroke and instroke paths, thereby increasing the cycle
work and the refrigerating effect.
This ingenious free-displacer split-Stirling concept was originated by
Walter Riga of the California Jet Propulsion Laboratory (Riga, 1965) and,
independently shortly thereafter, by Fred Chellis at Cryogenic Technology
Inc., Waltham, Massachusetts and Stuart Horn of the U.S. Army Night
Vision Laboratory, Fort Belvoir, Virginia.
The principal difficulty is the need to maintain the seal friction restraint
constant throughout the life of the unit. It is difficult to reproduce and
maintain constant the friction conditions of presently available dry rubbing
materials in a helium atmosphere. One solution (Walker, 1978) is to
incorporate a small auxiliary compressor feeding a separate duct to the
bounce space to provide a reproducible well-regulated fluidic drive force
on the regenerative displacer.
Horn et al. (1973) have presented an excellent theoretical analysis of
the free-displacer split-Stirling engine operation. No validation of the
analysis by comparison with experimental data has been published in the
open literature.
172 Chapter 3
The great majority of Stirling cooling engines in use are small systems
of low capacity used for local gas liquefaction or miniature electronic
instrument or infrared thermal sensors. Larger machines have been con-
structed, but presently there is no large demand for intermediate capacity
cooling systems. Future developments anticipated for the relatively near
term may change the situation. One development is the on-board
reliquefaction of vaporized liquid natural gas for marine LNG tankers. As
the "energy crisis" becomes more acute the marine transport of LNG is
increasing dramatically. It is possible to drive the propulsion engines of
the ship with the vaporized gas but it is difficult to arrange the rate of
boil-off to correspond exactly with the rate of consumption and some
reliquefaction may be required to eliminate the emission of combustible
vapors. Similar considerations apply to land-based peak-saving LNG
storage systems or tidewater terminal facilities.
Another area requiring cryocoolers of intermediate capacity is foreseen
for future superconducting applications. These include the superconducting
magnets of ship propulsion engines, magnetohydrodynamic power
generators, magnetic levitation of high-speed rapid transit systems, atomic
physics particle accelerators, and the magnetic plasma containment systems
of fusion power reactors. Superconducting electric cable transmission (likely
in conjunction with a liquid hydrogen/liquid natural gas cryogenic pipeline)
is another area where major requirements for intermediate cryogenic
cooling systems are required.
The future for these machines is, therefore, bright and it is timely to
look briefly at presently available machines and to speculate on future
development.
Fig. 3.48. Four-cylinder, single-acting integral-Stirling cryocooler, Philips Model No. PPG-
400.
and the mean pressure of the helium working fluid is 25 atmospheres. The
engine cool-down time is 10 min to attain rated capacity from a warm start
and a duty cycle of up to 300 hr of continuous operation is prescribed. The
unit requires 3 m 3 hr of cooling water at ambient temperatures. The weight
of the system including the bedplate and drive motor is 1000 kg, excluding
the concrete foundation. The installation envelope is 1m highx0.75 m
wide x 1.8 m long, again excluding the concrete foundation which adds
0.9 m to the height.
A machine of large capacity, again developed and manufactured by
Philips (Model No. PPG-2500) was shown earlier in Figs. 3.3 and 3.4. The
development of this machine has been described in detail by Dros (1965).
It is an interesting and unusual machine consisting of four single-acting
Stirling engine systems of the two-piston type in an opposed configuration.
174 Chapter 3
MODEL No PPG·400 .
SPEED
MEAN
- 1450rpm
~
0
PRESSURE - 25 otmos. u
COOLING l) 310• C
WATER 2)3m 3 /hr.
~
..
~
;!!;
3or----t----~----r----+----t----;
"'
.."'
;o
0
0
~
The unit embodies the hermetic rolling seal developed by Philips for their
Stirling prime-mover engines. It is driven by an electric motor supplying
shaft work to a vee-piston crank and connecting rod assembly (actually a
Joy air compressor unit). The pistons activate a hydraulic fluid connection
which operates on the compression and expansion space pistons causing
them to oscillate in the proper sequence. Massive external tie rods, easily
seen in Fig. 3.3, are used to contain the cylinder end caps. When the
machine is operating, these elastic members can be seen to be visibly
stretching and relaxing. It takes some getting used to this motion as normal
for machine operation.
The refrigeration characteristics for the engine, presented by Dros,
are reproduced in Fig. 3.50. The characteristics show the refrigeration
capacity (designated PE ), the input power (designated PM) and the coefficient
of performance relative to the Carnot value (designated as TilT/c) as a
function of the cold space temperature, TE. The discontinuity in the charac-
teristics, an increase in the slope of the refrigeration capacity curve and
the conversion to constant input power, PM, arises as a consequence of the
control system. When the power input to the drive motor approaches the
rated limit, 134 kW, working fluid is released from the cylinder to reduce
the pressure and so maintain the power input constant as the expansion
space temperature decreases further. Reduction in pressure of the working
Stirling Cryocoolers 175
0 .~ c
0
~
0 / ~~ 0 0 r- -!J!.nc ~
I ~~ ~~
30 0 ~60
20
VI}
/
/ ~
0 20
A
~ K' iO
0 / 1 / 0
/ H2 1/ He
0 0 0
0 50 100 150 200 250"1< 50 100 160 200 250"1<
--<>1£ -1£
ing the drive motor is 6500 kg excluding the concrete foundation. Overall
dimensions excluding the foundation are 1.95 m high, 1.6 m wide, and
2.97 m long.
*"The Werkspoor Cold Gas Refrigerator," Anon. VMF Review, 111(23), 125-131 (date
unknown).
Stirling Cryocoolers 177
DIS PLACER
EXPANSION
SPACE
COMPRESSION
SPACE
PI STON
.
o-
IZO
~
COEFF11CIENT OF PERFORMANCE
RELATIVE TO THE CARNOT ~
i
-
I
"''"
\\
o- 100 04 o-
i?
;!' -~ 0
o-
.",.
0 80
t'-.. 03~
'">
..
0
z
WERKSPOOR
MEAN
ENGINE
..~
/ ~ - 35 atmos
.'"
PRESSURE
>- 60
o- SPEED - 1!500rpm 02~
~ ~TEA
CYCLE INPUT
.
u COOLING - 25°C ~
.
~ 0
u ~
40
z
!0?
o-
r-- 0 ::'
~
0
~ REFRIGERATION CAPACITY o-
"''"
~ z
' 1.-----
20
'"
..~
----
u
::: Fig. 3.52. Refrigeration charac-
8 teristics of the single-cylinder
0 50 100 150 200 250 300 rhombic drive, integral Stirling
COLO SPACE TEMPERATURE (TEl(K) Werkspoor cryocooler.
120
y
I
~ 100
..
1-
:0
:;;
SHAFT INPUT
POWER P ~
:;:
lv
80
"~ I
0:
0
60
I
~
.
u
"u
"
~FULL LOAD
/
/
40
";::z
"
I
I
0:
"'~
k +
REFRIGERATION
20 /
~
"' / Fig. 3.53. Refrigeration capacity
and power input of the Werkspoor
-
0:
space and contaminating the regenerator, and egress of working fluid from
the working space). Replacement of the seals occupied four man-hours of
effort.
Vogelhuber and Parish presented data on installed and operating costs
of the Werkspoor engine liquefaction system and conventional Claude
cycle and cascade cycle systems. They concluded the cost for the Stirling
based system was generally favorable in liquefaction capacities up to
1.4 mmcfd. Such an installation would be a hefty liquefaction plant necessi-
tating duplicate four-cylinder Werkspoor engine with two-stage Joule-
Thompson expansion stages, precooling equipment, and all the associated
purification and storage facilities.
Nothing further about the Werkspoor engine was published. It is not
known if the engine is in production, if CVI still handles the engine in the
United States, and how many of the units were built.
The rhombic drive for Stirling engines was invented by Rolf Meijer
at Philips about 1952/53. To those familiar with Stirling engines, it is clear
from Fig. 3.51 that the Werkspoor engine owed much to contemporary
Philips engine prime movers. Werkspoor does much of the heavy machine
work and fitting of the Dros hydraulic engine and probably was familiar
with the earlier rhombic prime movers constructed by Philips.
Future Developments
There appear to be excellent future prospects for applications of
intermediate capacity engines as reviewed earlier. It is unlikely that any
180 Chapter 3
of the machines discussed above will enter volume production. They are
too complicated and expensive. There have to be better and easier ways
to accomplish the task.
There is every prospect that double-acting Stirling engines will come
to dominate the intermediate capacity field as inevitably as in the prime
mover automotive field. Double-acting engines have half the moving parts
of single-acting systems, so the economic motivation is overwhelming.
Double-acting engines can use a variety of kinematic drives, wobble-plates,
swash-plates, Scotch-yoke, crank and connecting rods. All these have been
discussed in more detail elsewhere (Walker, 1980). Already, multistage
four-cylinder wobble-plate helium compressors can be seen on prototype
development trials at Cryogenic Technology Inc., Waltham, Massachusetts.
With little imagination, one can foresee their adaptation to the Philips
swash-plate engine configuration shown in Fig. 3.15. Such developments
will relegate the Dros and Werkspoor engines to the status of historical
curiosities.
It is more than likely that double-acting free-piston Stirling engines
with hydraulic or linear electric power input will come to dominate the
field. In large capacities, these could be made so compact, so reliable
(elimination of seals, use of gas bearings, etc.), and so simple as to be
virtually unassailable as a commercial product. Already the literature
contains many references to free-piston Stirling engines as prime movers,
and Philips have in production a miniature double-acting free-piston Stirling
cooling engine (see Fig. 3.7). More significantly, a large U.S. manufacturer
of large-scale cryogenic equipment is presently working hard on a prototype
double-acting Stirling free-piston engine along the lines enumerated by
Finkelstein (1978).
REFERENCES
Babcock, G. H. (1885). "Substitutes for Steam." ASME Trans. 7, 680-741.
Beale, W. (1979). Private Communication. Sunpower Inc., Athens, Ohio.
Berchowitz, D. M., Rallis, C. J., and Urieli, I. (1977). "A New Mathematical Model for
Stirling Cycle Machine." Proc. 12th I.E.C.E.C., pp. 1522-1527, Washington, D.C.,
August 28-September 2.
Cairelli, J. E., and Thieme, L. G. (1977). "Initial Test Results with a Single Cylinder Rhombic
Drive Stirling Engine." Proc. ERDA, Highway Veh. Syst. Cont. Coord. Mtg., Dearborn,
Michigan, October 4-6.
Chellis, F., and Stewart, R. W. (1972). "Multiple Cold Finger Refrigeration System for
Cooling Infra-Red Detectors." AFFDL-TR-72-84, AD No. 903-939L (Cryogenic Tech-
nology Inc., Latham, Massachusetts).
Collins, S.C., and Cannaday, R. L. (1958). Expansion Machines for Low Temperature Processes.
Oxford University Press, Oxford.
Daniels, A., and du Pre, F. K. (1971). "Triple-Expansion Stirling Cycle Refrigerator." Adv.
Cryog. Eng. 16, 178-184.
de Jonge, A. K. (1979). "A Small Free-Piston Stirling Refrigerator." Paper No. 799245,
pp. 1136-1141, Proc. 14th I.E.C.E.C., Boston, Massachusetts.
Dros, A. A. (1965). "An Industrial Gas Refrigerating Machine with Hydraulic Piston Drive."
Philips Tech. Rev. 26(10), 297-306.
Finegold, J. G., and Vanderbrug, T. G. (1977). Stirling Engines for Undersea Vehicles."
Final Report No. 5030-63, J.P.L., March.
Finkelstein, T. (1959). "Air Engines." Engineer 207,492-497, 522-527, 568-571, 720-723.
Finkelstein, T. (1959). "Development and Testing of a Stirling Cycle Machine with Charac-
teristics Suitable for Domestic Refrigeration." English Electric Report W /M/3A.
Finkelstein, T. (1960). "Generalized Thermodynamic Analysis of Stirling Engines." S.A.E.
Paper No. 118B, January.
Finkelstein, T., Walker, G., and Joshi, T. (1970). "Design Optimization of Stirling-Cycle
Cryogenic Cooling Engines by Digital Simulation." Cryog. Eng. Conf., Paper K4, Boul-
der, Colorado, June.
Finkelstein, T. (1975). "Computer Analysis of Stirling Engines." Adv. Cryog. Eng. 20,
269-282.
Finkelstein, T. (1978). "Modeling Design and Optimization of Stirling Engines." Proc. 13th
I.E.C.E.C., San Diego, California, August 20-25.
Finkelstein, T. (1978). Pressure Compounding of Stirling Engines. Proc. 13th I.E.C.E.C.,
San Diego, California.
Gedeon, D. (1978). "Optimization of Stirling Cycle Machines." Proc. 13th I.E.C.E.C., San
Diego, California.
Goldwater, B., and Morrow, R. B. (1977). "Demonstration of a Free Piston Stirling Linear
Alternator Power Conversion System." Proc. 12th I.E.C.E.C., Paper No. 779249, pp.
1488-1495, Washington, D. C., August 28-September 2.
Haarhuiss, G. J. (1978). "The MC 80-A Magnetically Driven Stirling Refrigerator." Proc.
7th Int. Cryog. Eng. Conf., London, IPC Business Press, Guildford, U.K.
Higa, W. N. (1965). "A Practical Philips Cycle for Low Temperature Refrigeration." Cryog.
Techno!. 8, 203-209, July-August.
Horn, S. B., Lumpkin, M. L., Walters, B. T., and Acord, T. T. (1973). "Miniature Cryogenic
Cooler for TOW Night Sight." Proc. Closed Cycle Cryocooler Tech. and Appl., Vol. 1,
pp. 55-72, AFFDL-TR-73-149, WPAFB, Ohio, AD No. 918234.
Khan, M. I. (1962). "The Application of Computer Techniques to the General Analysis of
the Stirling Cycle." M.Sc. thesis, Univ. of Calgary.
Kirk, A. (1874). "On the Mechanical Production of Cold." Proc. Ins. Civil Eng. 37, 244-315,
London.
Kohler, J. W. L. (1965). The Stirling Refrigeration Cycle. Sci. Am. 212(4), 119-127.
182 Chapter 3
Lee, K. (1976). "The Stirling Cycle with Adiabatic Compression and Expansion." M.Sc.
thesis, Univ. of Calgary.
Lindale, E. (1978). "Stirling Cycle Refrigerators for Gamma Ray Detector." Report No.
PL-42-Cr78-0713, Johns Hopkins Univ., Applied Physics Lab., Laurel, Maryland (Philips
Laboratories).
Maki, E. R., and de Hart, A. 0. (1971). "A New Look at Swash-Plate Drive Mechanisms."
S.A.E. Trans. (Truck, Powerplant Fuels and Lubricants Mtg.), Vol. 80, Paper No. 710829,
St. Louis, Missouri, October.
Martini, W. (1978). Design Manual for Stirling Engines. NASA CR 135382, DOE/NASA
Contractor Report (NTIS, Springfield, Virginia). See also: Stirling Engine Design and
Feasibility for Automotive Use. (1979). (ed. M. J. Collie), Noyes Data Corp., Park Ridge,
New Jersey.
Prast, G. (1963). "A Philips Gas Refrigerating Machine for 20 K." Cryogenics 3, 156-160,
September.
Prast, G. (1965). "A Gas Refrigerating Machine for Temperatures down to 20 K and Lower."
Philips Tech. Rev. 26(1), 1-11, January.
Qvale, E. B. (1963). "An Analytical Model of Stirling-Type Engines." D.Sc. thesis, MIT.
Qvale, E. B., and Smith, J. L. Jr. (1969). "An Approximate Solution for the Thermal
Performance of a Stirling Engine Regenerator." Trans. A.S.M.E., I. Power 91, 109-112,
April.
Rios, P. A., Qvale, E. B., and Smith, J. L., Jr. (1968). "An Analysis of the Stirling Cycle
Refrigerator." Adv. Cryog. Eng. 13, 332-342.
Rios, P. A., and Smith, J. L. Jr. (1969). An Analytical and Experimental Evaluation of the
Pressure-Drop Losses in the Stirling Cycle. A.S.M.E. Paper No. 69-WA/Ener-8, Vol.
92, Ser. A., No.2, pp. 182-188, April.
Schmidt, G. (1861). "Theorie der Geschlossenen Calorischen Maschine Von Laubroy und
Schwartzkopff in Berlin." Z. Ver. Oster. Ing., p. 79.
Schmidt, G. (1871). "Theorie der Lehmannschen Calorischen Maschine." Z. Ver. dt. Ing.
15(1).
Schock, A. (1978). "Nodal Analysis of Stirling Cycle Devices." (a) Draft Report DOE; (b)
Proc. 13th I.E.C.E.C., San Diego, California.
Stoddart, D. (1960). "Generalized Thermodynamic Analysis of Stirling Engines." B.Sc. thesis,
Durham Univ.
Urieli, I. (1977). "A Computer Simulation of Stirling Cycle Machines." Ph.D. thesis, Univ.
of Witwatersrand, Johannesburg, South Africa, February.
Urieli, 1., Rallis, C. J., and Berchowitz, D. M. (1977). "Computer Simulation of Stirling Cycle
Machines." Proc. 12th I.E. C. E. C., Paper No. 779252, Washington, D.C., pp. 1512-1521,
August 28-September 2.
Urieli, I. (1979). A Review of Stirling Cycle Machine Analysis. Paper No. 799236, pp.
1086-1090, Proc. 14th I.E.C.E.C., Boston, Massachusetts.
Vogelhuber, W. W., and Parish, H. C. (1968). Compact LNG System Using Large Stirling
Cycle Cold Gas Refrigerator. Paper No. 28, Session 5, Proc. First Int. L.N.G. Conf.,
Chicago, Illinois (Institute of Gas Technology, Chicago, Illinois).
Walker, G. (1962). An Optimization of the Principal Design Parameters of Stirling Cycle
Machines. J. Mech. Eng. Sci. 4(3).
Walker, G., And Khan, M. (1965). "The Theoretical Performance of Stirling Cycle Machines."
S.A.E. Paper No. 949A, Int. Auto. Eng. Congress, Detroit, Michigan, January.
Walker, G. (1973). Stirling Cycle Machines. Oxford University Press, Oxford, England.
Walker, G., and Burn, K. (1976). "Exploratory Study of the Rainbow Variant Stirling Engine."
Proc. 11th I.E.C.E.C., Lake Tahoe, New York, September.
Walker, G. (1978). "Split-Stirling Cooling Engine with Fluidic Control of Displacer Motion."
Invention Disclosure, Calgary, Alberta, December.
Walker, G. (1980). Stirling Engines. Oxford University Press, Oxford, England.
West, C. (1971). "The Fluidyne Heat Engine." Report No. AERE-R-6776, Atomic Energy
Research Establishment, Harwell, Berks, England.
Stirling Cryocoolers 183
Zimmerman, J. E., and Radebaugh, R. (1977). "Operation of a SQUID in a Very Low Power
Cryocooler." Proc. of Conf. on Appl. of Closed Cycle Cryocoolers to Small Super-
conducting Devices, pp. 59-66, Boulder, Colorado, October (issued as NBS Spec. Pub.
508, April 1978-U.S. Govt. Printing Office).
Chapter 4
Vuilleumier Cryocoolers
INTRODUCTION
COLD
PRODUCTION
IW
77K
Fig. 4.1. Essential elements of a Vuilleumier cryocooler (after Pitcher and du Pre, 1970).
Advantages
The principal advantage claimed for Vuilleumier cryocoolers is the
potential for long life made possible by the low mechanical forces acting
on the drive mechanism, bearings, and seals. These low forces arise as a
consequence of the low pressure ratios (1.1 to 1.2) in Vuilleumier engines
compared with Stirling engines (2 to 2.5). The pressure ratios are low in
Vuilleumier engines for they are generated solely as the result of changes
in temperature, not both temperature and volume as in a Stirling engine.
Stirling cryocoolers are used where space and weight are at a premium.
The design is constrained to a machine operating at high pressure and
speed with maximum possible compression ratios to obtain the highest
refrigerating effect per unit mass or volume. The price to be paid for this
is a short operating life. It is often overlooked that it would be simple to
relax the output requirements imposed on Stirling engines and obtain the
same low speed, low compression ratio, low refrigeration output charac-
teristic of Vuilleumier engines. This would endow Stirling engines with the
same long life claimed as the principal attraction of Vuilleumier engines.
Disinterested comparative studies have shown the Stirling cryocooler
as more efficient and lower in cost, volume, and weight than the equivalent
Vuilleumier unit. The long-life characteristic of Vuilleumier engines has,
in fact, proved elusive despite a decade of intensive development. One
Vuilleumier Cryocoolers 187
particular difficulty has been the high wear rates of the hot rider, the sliding
guide ring, at the high-temperature end of the hot displacer.
The real advantage of the Vuilleumier cryocooler is the ability to use
a thermal energy input to produce refrigeration. This can be supplied by
concentrated solar energy, radioisotopes, combustion of fossil fuels or waste
materials, and finally, waste thermal streams of industrial plants or power
systems.
Solar and isotope energy made the Vuilleumier engine a strong candi-
date for cryogenic cooling in long-duration space exploration. This was
recognized by NASA, and exploratory developments have been undertaken
by the Goddard Spaceflight Center. Waste thermal streams or combustion
of waste products can be combined with Vuilleumier engines to produce
refrigeration at both cryogenic and higher temperatures for food processing
and air conditioning. Terrestrial applications of Vuilleumier engines for
air conditioning with solar energy are very attractive. In many locations,
air conditioning refrigeration is required only when the sun is shining. This
is a happier circumstance than the use of solar energy for power production,
where the main requirement for power arises at night.
In all these potential applications for Vuilleumier engines, similar
considerations apply to duplex Stirling systems, but so far as is known, no
comparisons have been made. Self-starting and complete independence of
an auxiliary electric power supply would weigh heavily in favor of duplex
Stirling systems. There are good prospects that an ingenious designer could
instill the same desirable characteristics in the Vuilleumier engine.
Thermodynamic Aspects
Most Vuilleumier cryocoolers have used electric resistance heating to
provide thermal input to the hot displacer cylinder and a motor to operate
the engine. Resistance heating for Vuilleumier engines is thermodynami-
cally inefficient. The second law of thermodynamics limits the fraction of
heat supplied that can be converted into work, the Carnot penalty. There-
fore, a high price is paid to produce electric power. Frequently, only 10%
to 30% of the energy supplied is transformed to electricity. The electric
power can be transformed in electric motors to useful work at relatively
high efficiency without payment of further thermodynamic penalties. The
work produced eventually becomes converted to low grade heat by way
of frictional dissipation.
The use of electricity in resistance heating produces no useful work.
All the energy is converted (degraded) to heat directly. If now that heat
is used in an engine, the Carnot penalty must be paid again. If the work
is used to drive a refrigerator, the refrigerator itself sustains a further
188 Chapter 4
Comparative Data
Chellis (1973, 1977) has given comparative data, reproduced in Tables
4.1 and 4.2 for a number of cryocoolers of different sizes and types. These
clearly show the weight, size, and efficiency penalties of the Vuilleumier
engines in comparison with equivalent Stirling machines.
HISTORICAL REVIEW
* I would be very glad to receive any information readers may have about the work, life, or
machines of Rudolph Vuilleumier, or about similar machines invented earlier in other parts
of the world.-G.W.
~
;=
i
Table 4.1. Comparison of Closed-Cycle Cryocoolersa
i:
--
Cooling Compressor Expander
capacity Input
Cryocooler power Size Wtb Size Wt MTBF
type Watts K (Watts) (in.) (Ib) (in.) (Ib) (hr) Application
Gifford-McMahon 3 20 2,100 28xl6x16 175 18x 10x6 22 20,000 Ground based
par amp
Gifford-McMahon 1 26 790 lOxlOxlO 20 11 x6x2.6 5 1,000 Airborne IR
Integral Stirling 1 26 475 14x 10x6 26 - - 250 Airborne IR
Integral Stirling 1 77 100 11 x4.5 x3 3.5 - - 250 Airborne IR
Split Stirling 1 77 135 10 X 5.5 X 2.5 5 3x2.5x2 0.3 1,000 Airborne IR
Split Vuilleumier 1 77 350 10x4x3 lOS 3 x2.5 x2 0.3 2,200 Airborne IR
a After Chellis (1977).
• Includes fan.
' Includes motor inverter and electronic temperature controller for hot end.
...
~
....
~
('l
=-
i
~
""
Vuilleumier Cryocoolers 191
l.,..__
ELECTRONICS
INTERFACE UNIT
fiiiFfiiiGlAATOfl
ATT~MUfT
('")
=-
Fig. 4.3 . Triple-expansion high capacity Vuilleumier cryocooler (after Renyer, 1973). i...
~
Vuilleumier Cryocoolers 195
Fig. 4.4. Philips Laboratories triple-expansion high capacity Vuilleumier cryocooler (after
Pitcher, 1975).
changes due to the ingress and egress of displacer rods, a bellows device
was incorporated in the crankcase system.
The rigidity of the rhombic drive mechanism permitted the use of
close tolerance regenerator seals (operating at ambient temperatures) in
both the hot and cold cylinders. Rubbing surfaces were entirely eliminated
in the "dry" helium-filled spaces of the engine to further enhance the
potential for extremely long lifetime.
On reading the Philips reports, one cannot help admiring the high
level of engineering competence displayed therein. At the same time, they
convey the impression they could have done a better job with a Stirling
cryocooler. Replacement of the Vuilleumier thermal compressor (heated
electrically) by a mechanical compressor (driven electrically) would have
converted the system to a Stirling cryocooler. There is little doubt it would
have been more efficient, lighter, smaller, and cheaper, with no sacrifice-in
fact, an enhanced potential-of long life.
Vuilleumier Cryocoolers 197
HEAT REJECTION, T0
HOT REGENERATOR
HEAT INPUT, Th
Support for this comes from the detailed failure analysis given by
Pitcher (1975). He concluded that the potential for a mechanical failure
rate was 6. 72 per 106 hours of operation and an electrical/ electronic failure
rate was 26 per 106 hours of operation. Of the total electrical/electronic
failures, 9.6 were attributed to the heater control. Elimination of the heater
would, therefore, reduce the electrical/electronic failure to 16 per 106
hours of operation. The use of a large motor for the increased work input
would reduce the prospect of motor controller failures (estimated at 9 per
106 hours of operation). At the same time, the mechanical failure rate
might increase to a balanced electrical and mechanical failure rate of 10
per 106 hours of operation.
Evidence of the more effective use of electric power in a Stirling engine
is contained in results presented by Pitcher (1975) as an appendix to the
main report. The results are reproduced in Table 4.3. It contains estimates
198 Chapter 4
Table 4.3. Estimate of Power Requirements for Multistage Stirling and Vuilleumier
Cooling Engines with Different Temperature Levels and Cooling Capacities"
Inlet
Cold Power coolant
Temperature (K) production (W) input temp.
Entry
no. 2 3 2 3 (W) (oC) Comments
INCH[$
0 I
L....L_J
SCALE
Fig. 4.7. Cross section of miniature Vuilleumier cryocooler (after Pitcher, 1970).
... selection of the Boeing material was based on wear rates, the predictability
of this material over extended period of operation and, to a lesser extent, a
low coefficient of friction. Flame-sprayed material was selected for the mating
surface because of wear rates, coefficient of friction and fabrication consider-
ations, primarily the ease of flame spraying the desired bearing mating
material.
Chrome carbide is used at the hot end of the refrigerator which may
operate at temperatures as high as 1200°F. Tungsten carbide is used for the
remainder of the bearings which operate at temperatures below 200°F.
Tungsten carbide cannot be used at the elevated temperatures of the hot end
due to the difference in thermal expansion between the tungsten carbide and
the base material Inconel 718.
MAGNETIC COUPLING
.:!Q
NNECTI NC.~ LINEAR BEARINGS HOT REGENERAT OR g
Q
FOR HOT DISPLACER
.."';-
COLD DIS PI A~FR
COLD REGENERATOR
COLD END HEAT
EXCHANGER
EVACUATED
SUPERINSULATION HOT D I SP LA~ F R
SUMP HEAT HOT END
COLD REGENERATOR SEAL
EXCHANGER EXCHANGER
LINEAR BEARINGS AMBIENT HEAT
FOR COLD DISPLACER PIPE ASSEMBLY
Fig. 4.8b. Cross section of AiResearch Vuilleumier cryocooler, 5 W capacity at 77 K, for NASA-Goddard Space Flight Center (after Browning
eta/., 1973).
.,.
(.>
=
204 Chapter 4
EGENERATOR
COLD EXPANSION VOLUME
Fig. 4.9 . Different freezer configurations for the AiResearch/ NASA Vuilleumier cryocooler
(after Browning eta/., 1973).
~~;:;~r:;,~~~~~~~~-::FPERFORATED PLATE
FLOW DISTRIBUTOR
COLD REGENERATOR
Fig. 4.10. Preferred freezer arrangement for the AiResearch/NASA Vuilleumier cryocooler
(after Browning eta/., 1973).
system might increase the mass of the engine by 0.5 to 3.6 kg. The saving
in mass would, therefore, be (19.9 -4.5) = 15.4 kg.
Ackermann (1971) made a theoretical analysis and parametric optimi-
zation of Vuilleumier cooling engines for the Night Vision Laboratory, Fort
Belvoir, Virginia. An analytical paper for steady state operation of an
idealized Vuilleumier cooling engine was contributed by Rule and Qvale
(1969). Qvale, then an assistant professor at Purdue University and now
at the Technical University of Copenhagen, had earlier worked with
Professor Smith at the Massachusetts Institute of Technology in the
development of a simulation program for Stirling engines.
EDITOR's NoTE: The following sections-Cycle of Operation, Inher-
ent Thermodynamic and Heat Transfer Losses, VM Cooler Variations,
and VM Accessories and Components-are reproduced, with minor
editorial changes, from the unclassified United States Air Force report
by Ronald White (1976). The report was prepared by White to satisfy
206 Chapter 4
HOT lt(C(N£AATOA.
SPAC[
PR£SSUM( Vf SSH
\I'AU.
Fig. 4.11. Cooler arrangement for the AiResearch/ NASA Vuilleumier cryocooler (after
Browning et al., 1973).
CYCLE OF OPERATION
low temperature and then rejects this heat at some intermediate tem-
perature. The heat input at the hot end of the machine provides the energy
required to produce the cooling effect at the cold end, while the heat from
both the hot and cold ends is rejected at the ambient (intermediate) end.
The VM refrigerator, in its simplest form, is composed of two cylinders
with displacers, two thermal regenerators, three heat exchange areas,
connecting passages, and a mechanism to drive the displacers in the proper
sequence.
The displacer is a long loose-fitting pistonlike device whose function
is to move (displace) gas from one end of the cylinder to the other (see
Fig. 4.14). The cylinder-displacer combination is arranged so that there is
a gas passage connecting the two ends of the cylinder. With this connection,
the effect of moving the displacer is to displace the gas from one end of
the cylinder to the other. Note that the total volume of the gas remains
constant (assuming an infinitely thin push rod). This is in contrast to a
piston-cylinder combination that would change the gas volume. Since the
COOLING COILS OR HEAT PIPES
FLOW DISTRIBUTION PASSAGE
ALUMINUM SUMP
SUMP PRESSURE
COOLING COLLAR
VESSEL WALL
PORT TO BACKSIDE
DISTRIBUTION OF HOT DISPLACER
FLOW PATH \-----_SUMP FILLER
BLOCK
Fig. 4.12. Heater arrangement for the AiResearch/NASA Vuilleumier cryocooler (after
Browning et al., 1973).
208 Chapter 4
75
70 /
~
/
w65
i"' / /
~ 60
"'...
/
/ TEST CONDITIONS
55
v PRESSURE 800 PSIG PEAK -
SUMP TEMPERATURE 14o" F
HOT END TEMPERATURE 1000°F Fig. 4.13. Performance characteristics
SPIEED 365 RiM of the prototype AiResearch/NASA
50
I Vuilleumier cryocooler (after Brown-
0 7
only VM that utilizes a piston is the special-case self driven VM, the term
"displacer" will be used throughout this report as a reminder that the total
gas volume of the system is not affected by the movement of the displacers.
The pressure on each end of the displacer is theoretically the same (in real
refrigerators it differs by the pressure drop through the regenerator and
connecting passages). Therefore, the energy needed to move the displacer
is very low, especially when compared to a refrigerator using pistons. Since
the two ends of the cylinder will be at significantly different tempera-
tures the displacer is long with thin walls and packed with insulation to
reduce the end-to-end heat conduction.
Next, a thermal regenerator is arranged to connect the two ends of
the cylinder so that gas passing from one end of the cylinder to the other
must pass through the regenerator. The thermal regenerator is an energy
storage device whose heat capacity greatly exceeds that of the gas. In
operation, gas passing through a regenerator is heated or cooled depending
on flow direction. For example, gas flowing from the ambient-temperature
end to the high-temperature end of a cylinder is heated by the hot regen-
erator. The energy added to the gas was stored in the regenerator packing,
DIS PLACER
GAS PASSAGE
CYLINDER
CYLINDER
REGENERATOR
MATRIX
DISPLACER
SEAL
efficient with no significant gas void spaces. Therefore the gas in the system
is distributed between the two (hot and ambient) active volumes. The term
"active volume" is used to describe the volume swept by the displacer.
When the displacer is in the top position the hot active volume is nearly
zero and all of the gas is in the ambient active volume. The ambient volume
is held at an ambient temperature (significantly below the hot temperature)
by the ambient heat exchanger. Using the perfect gas law, the gas is at low
temperature (volume is constant), therefore it is at a low pressure. When
the displacer is in the bottom position the ambient active volume is nearly
zero and all of the gas is in the hot active volume. Since this volume is
maintained at a high temperature (usually about 1200°F, 920 K) by addition
of heat through the hot end heat exchanger, the gas is hot and therefore
at high pressure. Since there are no valves in the system, all of the gas
throughout the system including any at the other end of the cylinder (or
in any dead volumes to be discussed later) is at the same pressure. Therefore,
it can be seen that by moving the hot displacer, the system pressure can
be varied from a maximum to a minimum in a cyclic way and that heat is
supplied to the hot end and rejected at the ambient end. The hot section
acts as a thermal compressor (whose effect on pressure is similar to a piston
in a dead ended cylinder). A similar thermal compressor was patented by
V. Bush (1939).
The manner in which the cold section, the second cylinder-displacer,
of the VM refrigerator uses the changing pressure to produce a cooling
effect is similar to that of the cold section of the Stirling cycle refrigerator,
since the method by which the cyclic pressure is produced (thermal
compressor or mechanical compressor) has little effect on the cold section.
Referring to Fig. 4.16, if the dis placer is in the top position and pressure
from the compressor is increasing, gas will flow into the ambient active
volume (the cold active volume is zero) and increase in pressure. This
compression causes the temperature of the ambient volume gas to increase,
but this heat is rejected by the ambient heat exchanger so that at the end
of the compression process the ambient active volume is filled with high-
pressure ambient temperature gas. The cold displacer is then moved to its
lower position and the gas is displaced from the ambient end, through the
cold regenerator to the cold end. As the gas passes through the cold
regenerator, it deposits its heat in the regenerator matrix and emerges at
the cold end as cold high-pressure gas (any decrease in volume of the gas
is made up by the compressor). After all the gas had been displaced to the
cold active volume, the compressor displacer is moved so that the pressure
decreases to the minimum. This decreases the pressure throughout the
refrigerator including the cold section. The gas in the cold active volume
undergoes an expansion (some of it goes back through the regenerator
Vuilleumier Cryocoolers 211
VOLUME
REGENERATOR
MATRIX
CYLINDER
toward the compressor) becoming colder low-pressure gas. This colder gas
then has the ability to absorb heat from the device to be cooled (refriger-
ation). After the gas has absorbed the heat, the cold displacer is moved to
the upper position, displacing the gas from the cold volume through the
regenerator where it picks up the heat that was deposited on the previous
half cycle and emerges at the ambient end where it will reject the heat
picked up from the refrigeration load. After this displacement is completed,
the compressor starts to increase the pressure and the cycle begins again.
From this description it can be seen that by combining two cylinder-
displacer-regenerator assemblies and phasing their motion properly, useful
refrigeration can be obtained.
To this point, nothing has been said about what sort of mechanism
drives the displacer push rods. In the previous discussion a stop-start square
wave motion was implied to simplify the discussion of the events in each
section. In almost all of the real VM refrigerators built to date, the displacers
have been driven harmonically. Harmonic drives (such as a crankshaft with
connecting rods) are easy to fabricate and avoid the problem of high
acceleration at the end of the stroke. However, harmonic drives do compli-
cate the discussion of the VM cycle. A schematic of a harmonically
(sinusoidal) driven VM refrigerator is shown in Fig. 4.17.
The cycle operates through the use of displacers moving the gas from
one section to another without the requirement to compress the gas in a
closed volume. Therefore, the pressures throughout the system are nearly
equal at any moment. The seals shown in the schematic are to force all of
the displaced gas through the regenerators. Since the pressure drop across
the regenerators in real refrigerators is only a few pounds per square inch,
the loading on the seals and displacer drive bearings is minimal, which
212 Chapter 4
GAS-FILLED
WORKING VOLUME HEAT EXCHANGER
(TYPICALLY HELIUM ABSORBS HEAT
AT HIGH PRESSURE) FROM LOAD AT LOW
TEMPERATURE
COOLING
CYLINDER
HEAT EXCHANGER
REJECTS HEAT
TO AMBIENT
(HOT)
POWER
Q DISPLACER
CRANKSHAFT
CRANK THROW
CRANKCASE
SOURCE OF
HEAT AT
HIGH TEMPERATURE 1 - - - - - '
(-1200 "F)
assumed that there are no pressure drops within the refrigerator in the
idealized model. The thermal regenerators are assumed to be perfect; i.e.,
no temperature difference is required between gas flows in each direction
in order to transfer heat, and therefore no heat flows through the regen-
erators over a complete cycle. The regenerators are assumed to have an
infinite heat capacity and therefore the temperature of the regenerator is
invariant with time. In the idealized refrigerator it is assumed that there
are no dead volumes, that is, all the volume inside the refrigerator is active
volume (VH, VA, and Vc). This means the gas volume in the regenerators,
heat exchangers, around the sides of the displacers, around crankcase parts,
at the ends of cylinders, etc. is assumed negligible (this will be modified
later). The heat conduction axially through the regenerator matrix, and
along the walls of the displacers and cylinders, is assumed negligible.
Attempts have been made to describe the cycle using a temperature-
entropy (T-S) diagram such as Fig. 4.18. It is apparent that the diagram
describes more than one unit mass of gas and the diagram is somewhat
artificial since a great many simultaneous operations occur. For an example,
let us look at three masses of gas that cycle between the ambient end and
the hot end. The mass of gas in the hot end nearest the hot regenerator
entrance will be the last mass in and (if no mixing occurs) the first mass
out. Its T-S diagram would look like Fig. 4.19a. The second mass (between
the first and last) would look like Fig. 4.19b. The third mass of gas is the
first one in the hot end and the last one out. It would have a T-S diagram
that looks like Fig. 4.19c. It can be seen that in order to use the T-S
diagram to describe this cycle a large number of individual T-S diagrams
for both ends of the refrigerator would need to be summed (mixing could
further complicate the analysis).
However, the pressure-volume (P- V) diagram for each volume
provides a clearer and more accurate picture of the operation than does
the T-S diagram.
The P- V or indicator diagrams for the three sections of a well-designed
V -M refrigerator (P 1 = P 2 = P max and P 3 = P 4 = P min) are shown in Figs.
4.20a-4.20c. Note that the ambient volume is the sum of the ambient end
volumes of the hot and cold cylinder.
214 Chapter 4
TH ) ; ,' 1,' I /
.... I 1 I
I I I I
W I I
~ //
=
1- I
I
~ I
~ TA / / I I
.... 'z ',
Fig. 4.19. Temperature-entropy diagram for the mass of gas in the hot space that is (a) last
to enter, first to leave the hot space; (b) arbitrarily intermediate the first and last elements
of gas to enter or leave the hot space; and (c) first to enter, last to leave the hot space (after
White, 1976).
The heat flow is equal to the work at that volume. There is, of course,
no net work resulting from the idealized refrigerator because no pressure
differences exist within the machine at any given time. Since the pressure
at each end of each displacer is equal and by elementary geometry the
change of volume at one end of a cylinder is exactly equal to the negative
of the volume change at the other end, we have
Since the displacer motion has been assumed harmonic, the pressure
changes are assumed to occur isothermally, and if the phase separation
angle is assumed to be 90°, the volume at each of the ends can be described
(a)PMAX
1~
12
~ I
I I
I!' I
ijl I
"'~ '
I
PMIN - - ,41 - 3
VOLUME-V ---+-
Fig. 4.20. Pressure-volume (indicator) diagram for the (a) cold volume; (b) hot volume, and
(c) ambient volume of a Vuilleumier cryocooler (after White, 1976).
Vuilleumier Cryocoolers 215
n= ~ [ VHM( 2 ~H + 2 ~J + Vc~ 2 ~c + 2 ~J
+ VH~ 2 ~A- 2 ~H) COS 8+ Vc~ 2 ~A- 2 ~J sin 8]
multiplying numerator and denominator by 2TA the equation becomes
n = 2 ~R [ VH~~;+1)+ VcM(~;+1)
+ V H~ 1 - ~;) cos () + V c~ 1 - ~;) sin 8J (4.5)
c=Vc~1-~;) (4.8)
Where a reflects the charge of gas in the refrigerator, b reflects the effect
of the movement of the hot end displacer on the total pressure, and c
reflects the effect of the cold-end displacer on the pressure.
216 Chapter 4
Oc f
= P dVc = nRTA V CM f2 1T
o a
b
cos 0
.
+ cos 0 + c sm 0
dO (4.11)
where QH is the heat input to the hot end and Oc is the heat input to the
cold end.
These equations can be integrated using formula 2.558-2. from Table
of Integrals Series and Products by I. S. Gradshteyn and I. M. Ryshik. Note
that a 2 >(b 2 +c 2) and that JdO/(a +b cos 0 +c sin 0) must be integrated
from 0 to 7T and from 7T to 27T with -'TT/2 :s; arctan 0 :s; 7T/2. The result is
OH
27TcnR
= b2+c2
[ a
VHMTA 1- (a2-b2-c2)1/2
J (4.12)
27TbnR [
Oc =- b2+c2 VcMTA 1
a
(a2-b2-c2)1/2
J (4.13)
These are the equations for the heat input to the hot end and the heat
absorbed by the cold end (per cycle) for an ideal VM refrigerator.
The relationship between the maximum hot active volume and the
maximum cold active volume is
VHM TH(TA- Tc)
(4.14)
VcM Tc(TH-TA)
This equation is derived from the ideal work and the ideal refrigeration
that takes place in the machine during a cycle, when the pressure at crank
position 1 is assumed equal to the pressure at crank position 2 as indicated
in Fig. 4.17.
The coefficient of performance of this ideal refrigerator is
(4.15)
Vuilleumier Cryocoolers 217
9- DEGREES
.
400:
..
300 ~
zoo ~
.
0
'00
cold volume, while the displacer side clearance void volume is usually about
0.17 times the active cold volume. The effect of void volumes is related to
the temperature of the void volume. The void volumes at the coldest
temperatures have the greatest effect on the refrigerator performance while
the void volumes at the hot end have the least effect on refrigerator
performance. The effect of void volumes is more pronounced in miniature
VM refrigerators, due in part to ordinary manufacturing tolerances. In
miniature VM refrigerators, decreasing the void volume is both essential
and expensive. Reduction of void volume requires close tolerances and
unusually shaped parts (especially in the crankcase). Void volume reduction
techniques used in current refrigerators even include the use of epoxy to
fill the screw slots in the internal screws.
The effect of void volume can be accounted for by adding terms to
Eq. (4.4) of the form:
where Vvl is a void volume and Tv 1 is the temperature of that void volume,
Vuilleumier Cryocoolers 219
The other equations (4.9)-(4.16) remain the same except the numerical
value of a has changed. Typical values of the reduced void volume ratio
(reduced void volume to cold end volume) for small VM refrigerators built
to date are in the range of 1.5 to 3.7, where reduced void volume is defined
as
Tc
Vv reduced = Vv actual X -T-- (4.18)
v actual
(4.24)
220 Chapter 4
Equations (4.9)-(4.13) and (4.16) are still valid, however, the values of a,
b, and c within these equations have changed slightly.
To summarize the theoretical VM section, if compressibility effects
and void volumes are to be accounted for, the equations to use are
(4.9)-(4.13), (4.16), and (4.18)-(4.24). Equation (4.13), Oc, is the gross
refrigeration produced by the refrigerator. To obtain the net refrigeration
produced by the refrigerator, all of the various cold end losses must be
subtracted from the gross refrigeration. Equation (4.12), QH, is the P- V
heat input to the hot cylinder. To obtain the actual heat (power) input
needed, the various hot end losses must be added to the P- V input.
Shuttle Loss
Shuttle loss is caused by the mismatch of thermal gradients between
the displacer and the cylinder. The cold cylinder wall is at ambient tem-
perature at one end and at cryogenic temperature at the other. It has a
fixed length and a gradient from warm to cold that is approximately linear.
The displacer is shorter than the cylinder wall by the length of the stroke,
however, it has the same temperature extremes: warm at one end and cold
at the other. When the displacer is at one extreme of its travel the tem-
perature gradients are somewhat mismatched. As it passes through its
stroke and reaches the other extreme the gradients are again mismatched,
but now they are mismatched in the opposite direction so that the displacer
picks up heat from the cylinder when it is at the warm end and it gives off
heat to the cylinder when it is at the cold end of its stroke. Hence, there
is a picking up of heat at the warm end, a shuttling of the displacer to the
cold end, where it drops off the energy to the cylinder. It is thus termed
a shuttle loss.
If the motion of the displacer is approximately harmonic and if the
thermal time lag of the cylinder and displacer materials is small compared
with the reciprocating time, the shuttle heat transfer can be computed from
Vuilleumier Cryocoolers 221
Pumping Loss
The pumping loss is due to the fixed clearance volume that must exist
between the displacer and the cylinder so that the displacer can move
without rubbing. This volume is bounded on one end (usually the ambient
end) by a displacer seal and is open on the other end. Because this volume
around the displacer allowed for running clearance is fixed, the mass of
gas in this volume at any one time is proportional to the (cycling) pressure.
At the minimum pressure point there is a minimum mass of gas in this
volume. As the pressure increases toward the maximum pressure point,
gas flows into this volume proportional to the increasing pressure. There-
fore, using the cold cylinder as an example, cold gas flows into this volume
and down into the warm areas where it picks up energy. Then as the gas
pressure falls from the maximum back to the minimum pressure, some
warm gas flows out of this volume into the cold regions. This causes a loss
of refrigeration called the pumping loss. This loss occurs on both the cold
end and the hot end of VM refrigerators when displacer seals are used
to force gas flow through the regenerators. The following equation for
pumping loss was derived by Leo (1971):
where Qpu is the pumping loss, De is the diameter of the cylinder, Ley is
the length of the cylinder' p max is the maximum pressure, p min is the
minimum pressure, N is the cycle speed, Cv is the specific heat of the gas,
T w is the temperature at the warm end, Tc is the temperature at the cool
end, S is the radial clearance between displacer and cylinder, Z is
the compressibility factor for gas, R is the gas constant, and kg is the
conductivity of the gas.
222 Chapter 4
Q = koAv (T _T ) (4.27)
D Lv w c
Q _7Tko(D~v-Dio)(Tw-Tc) (4.28)
v- 4Lo
where OeY is the conduction heat loss through the cylinder, key is the
conductivity of the cylinder material, Pmax is the maximum pressure, Dey
is the inside diameter of the cylinder, Ley is the length of the cylinder,
and u is the allowable stress.
Vuilleumier Cryocoolers 223
where F 1 is the friction force in the direction of travel and Y is the stroke.
It is rather difficult to evaluate this loss since much must be known about
the seal forces, rider loads, and the locations at which this loss is generated.
Regenerator Losses
The regenerators of the VM refrigerator must be very efficient if the
refrigerator is to be a practical device for producing refrigeration at tem-
peratures below 100 K. In fact, the overall performance of the machine
depends directly on the efficiency of the regenerators. The regenerator
losses can be divided into two types:
Heat Load Due to Friction in the Regenerator. This refrigeration loss
is caused by aerodynamic heating when the working medium flows through
the regenerator:
(4.31)
or
_ cPmax+Pmin) 3 V~MN 3[LR
Q RF- 3 3 3 2 2 (4.32)
2Z R T cgcAregP D.J
where QRF is the loss due to aerodynamic friction in the regenerator, Pmax
is the maximum system pressure, P min is the minimum system pressure,
VcM is the maximum displaced cold volume (if a hot regenerator is being
calculated this would be maximum displaced hot volume), Tc is the tem-
perature of the displaced volume (hot or cold), N is the cycle speed, f is
the friction factor, LR is the length of the regenerator, Z is the compressibil-
ity factor for the working fluid, R is the gas constant for the working fluid,
Areg is the cross-sectional area of the regenerator, p is the density of the
working fluid, D. is the equivalent hydraulic diameter, gc is the gravitational
conversion factor, and J is the conversion factor.
224 Chapter 4
Q =( 1 -N )C (T -T)(Pmax+Pmin)VcMN (4.33)
RH R P w c 2RTcZ
where
h'h"L~
(4.34)
Where QRH is the loss due to the limiting value of the film coefficient in
the regenerator, NR is the regenerator efficiency, Cv is the specific heat of
the working fluid, (Tw- Tc) is the temperature difference across the regen-
erator, VcM is the maximum displaced cold volume (if a hot regenerator
is being calculated this would be the maximum displaced hot volume), Tc
is the temperatare of the displaced volume (hot or cold), N is the cycle
speed, R is the gas constant of the working fluid, Z is the compressibility
factor, h' is the heat transfer coefficient-flow out, h" is the heat transfer
coefficient-flow in, LR is the length of the regenerator, G is the mass
velocity (mass flow rate per unit cross-sectional area), and rh is the hydraulic
radius. The value of NR for cold regenerators is usually in the range
from 0.995 to almost 1.0 and for hot regenerators is in the range of 0.9
to 1.0.
These equations for the regenerator losses are only approximate since
they are based on average mass flow rates rather than instantaneous mass
flow rates. Another way to compute the performance of the VM refrigerator
was partially outlined by Globe et al. (1971). It involves cutting the VM
refrigerator into a large number of control volumes and determining the
mass flow rate into and out of each of these control volumes as a function
of crank position. The losses, especially the regenerator and pumping losses,
can be calculated and their effect on the pressure can be iterated with the
mass flow equations so that a more accurate description of the refrigeration
at every crank position is obtained. This can be done for an existing
refrigerator design with the aid of a computer but is extremely difficult to
do when optimizing a new design since the control volumes themselves are
being changed during the optimization process. Therefore, the equations
presented above or similar equations are usually used for optimizing new
designs.
Vuilleumier Cryocoolers 225
Net Refrigeration
To obtain the net refrigeration, the cold end losses for each stage are
summed and subtracted from the gross refrigeration for that stage. In
multistage refrigerators some of the losses from colder stages appear as
increased refrigeration at warmer refrigeration stages. These include con-
duction losses, regeneration heat loss, shuttle loss, and pumping loss. These
losses from a colder stage are called interstage heat flow and should be
added to the gross refrigeration of the warmer stage. Other cold end losses
that were not discussed here but might be applicable depending on
refrigerator design are regenerator conduction, insulation, seal leakage
losses, and heat leaks down instrumentation leads.
VM COOLER VARIATIONS
Multistage VM Coolers
In many applications there is a need to produce useful refrigeration
at more than one temperature level at the same time. An example of this
is cooling an electronic device to a very low temperature while cooling a
Dewar heat shield surrounding the device to an intermediate temperature
to intercept the Dewar heat leaks. This saves considerable power since the
majority of the heat is removed at the higher temperature for much less
input power (and refrigerator size) than if all the heat had been removed
at the lower temperature.
Two mechanisms have been used to add additional cold stages to VM
refrigerators. The parallel cylinder system (Fig. 4.22) adds another smaller
diameter cold cylinder parallel to the first cold cylinder. In this configuration
the displacer seals are at the ambient temperature region of the displacers
(an advantage) but fabrication of the cylinder assembly is more difficult.
Keeping the two cylinders straight and parallel from brazing temperature
down through cryogenic temperature is difficult and expensive.
226 Chapter 4
Fig. 4.22. Multiple stage Vuilleumier cryocooler arrangement with parallel cylinders (after
White, 1976).
11 t STAGE
REGENERATOR
tlf STAGE
OISPLACER
Phase Angle
Phase angles other than 90° were investigated by E. B. Qvale and T.
T. Rule (1968) and by B. Leo (1971). Phase angles other than 90° complicate
the Eqs. (4.1)-(4.13) considerably and make the fabrication of parts more
difficult. The optimum phase angle is a function of the active swept volumes
of the hot and cold cylinders and of the thermal boundary conditions. The
investigations showed that for refrigerators with heat rejection tem-
peratures near room temperature, the optimum phase angle would be in
the range from 90° to 102° depending on cold end temperature (the lower
the temperature the greater the angle). Multistage refrigerators further
complicate these relationships and reduce the range of optimum phase
angles. Since the gross refrigeration is changing very slowly with respect
to phase angle near the optimum, most manufacturers are using a 90° phase
angle.
Similar Cycles
There are other heat powered refrigeration cycles similar to the VM
cycle. One by Bush (1938) is quite similar to the VM except that the two
ambient volumes (one at the ambient end of the cold cylinder and one at
the ambient end of the hot cylinder) are separated by a thermal regenerator
and reject heat to heat sinks at different temperatures. Another heat-
powered refrigerator was patented by Taconis (1951). It differs from the
VM in the timing of the movements (three instead of four motions) of the
displacers. Another heat-powered refrigerator was patented by Hogan
(1964). It produced cooling in the 10-20 K range, while the hot end
absorbed heat at room temperature and the heat rejection was at 77 K
(the heat was rejected to liquid nitrogen). A patent by Cowans (1968)
describes a modification to the VM refrigerator that allows it to drive its
own dis placers and produce useful shaft power. This is done by increasing
the cross-sectional area of either the hot displacer connecting rod or both
connecting rods, so that with the addition of connecting rod seals, and by
lowering the crankcase pressure below the minimum pressure in the VM
cycle, a net force can be created to drive the refrigerator. This has the
advantage that the small timing motor used on most VM refrigerators is
228 Chapter 4
not needed (but something must give it a shove to get it started). However,
this adds the life-limiting problem of dynamic connecting rod seals that
must be able to seal against the full cycle pressure (several hundred pounds
per square inch). This type of sealing problem is avoided by most VM
refrigerators since in a "pure" VM cycle the only dynamic seals in the
system are the displacer seals. Displacer seals usually experience very small
pressure differences of 5 to 15 psi which contributes to their very long life.
ing, slow on-off switching, and mechanical devices such as Curie point
switches and vapor bulb thermometers. On large VM refrigerators the
controller problem is reduced by calculating both the minimum power and
the maximum power required. The minimum power is then supplied by a
large heater with a simple on-off controller for malfunction protection
only. The difference between the minimum and maximum power is supplied
by a smaller heater with a proportional controller and necessary shielding.
This arrangement is more efficient and reduces the size and weight of the
controller and electronic filters.
Heaters
Two types of electric hot end heaters have been used in VM
refrigerators. The furnace type is a ribbon of heater wire wrapped on a
ceramic mandrel and held in place with cement. The ceramic furnace
surrounds the hot end of the hot cylinder and transfers heat by radiation.
The Calrod-type sheathed heater has been the most popular. Heat is
transferred either by radiation (Yoshikawa, 1970) or by brazing the heater
sheath directly to the hot cylinder. Since the watt density of the heaters
required by most VM designs is very high (for this type of heater) the
heater must be properly heat sunk to the hot cylinder or burnouts will
occur. Care also must be taken to be sure that the active (heat-producing)
portion of the heater terminates while still thermally connected to the hot
cylinder. The larger-diameter low-resistance lead-in wire that runs between
the heater wire and the terminal (inside the sheath) must be of a material
that does not embrittle or corrode when exposed to insulation or atmos-
pheric contaminants. Single-ended straight wire heaters have caused
numerous failures and have been largely abandoned in favor of two-ended
helically wound single wire heaters. Straight wire heaters are available in
smaller sheath diameters but must use smaller-diameter heater wire since
the total length of the heater wire is less. This, coupled with the possible
nonuniform reduction of heater wire diameter during the swaging of the
heater sheath and higher stresses during heater cycling have contributed
to numerous straight wire heater failures. The helical single-wire heaters
have larger heater sheaths (less convenient for the refrigerator designer)
but have a larger-diameter longer heater wire that does not change cross-
section (the helix angle changes) during the swaging of the heater sheath.
The heater wire is closer to the sheath (less temperature drop) and is less
sensitive to thermal cycling. These heaters are much more reliable than
straight wire heaters. In critical applications additional redundant heaters
are added to avoid scrapping an expensive hot cylinder assembly due to a
burned out heater.
230 Chapter 4
Motors
As mentioned earlier, ac induction motors with the rotor inside the
helium space and the stator with its windings outside the helium space
have been successfully used in a large number of applications. These motors
are either two or three phase and are purchased with a matching inverter.
Total efficiency for the motor and inverter is about 25%. In small VM
refrigerators, the motor power is a small fraction of the total power,
therefore cleanliness and reliability are more important than efficiency. In
a few applications an inverter could not be used (due to space or ambient
temperature problems) and since life was less critical a de brush-type motor
was used. Special brush materials were used along with special commutator
coatings. The motor windings were potted to reduce the generation of
contaminates. Very few data have been gathered on this motor, so its
limitations are still unknown. A brushless de motor is being used on one
refrigerator. The motor efficiency is expected to be at least 55%. It will
be in the oil-filled crankcase of the refrigerator and should present no
contamination problem.
Regenerators
The cold regenerators have been previously discussed and are discussed
in considerable detail in all of the references. The term "internal regen-
erator" is used for regenerators inside the displacer and the term "external
regenerator" is used for regenerators attached to or a part of the cylinder.
Cold regenerator matrix materials are usually screens of 100 mesh to 500
mesh in copper alloys or stainless steel and balls of monel or lead in sizes
down to 0.002 in. Lead balls are usually used for temperatures below 50 K
since lead is one of the few materials with appreciable specific heat at these
temperatures.
The hot regenerators are described as internal or external als0,
however, the forms of this matrix are more varied. Stacked screens, balls,
tubes, and the annulus have been used for the hot regenerator matrix. The
internal annular regenerator is composed of the walls of the displacer and
the cylinder. In this configuration the gas flows between the displacer and
the cylinder. A displacer seal is not used. This eliminates one wearing part
(the seal) and eliminates the pumping loss, but makes the radial location
of the displacer within the cylinder very critical, which in turn makes hot
rider ring wear extremely critical. As an example, if the regenerator is
designed with a 0.007-in. radial gap between the displacer and cylinder,
and if rider wear allows the displacer to be out of concentricity by 0.002 in.,
15% of the gas flow is on the narrow side of the regenerator while 85%
Vuilleumier Cryocoolers 231
of the flow is on the wide side. This causes the regenerator loss to be
doubled (Leo, 1971). Since lubrication and wear of the hot rider is a serious
problem the internal annular regenerator is rarely used anymore. The
external annular regenerator is composed of the walls of the cylinder and
one or more liner sleeves. The liner sleeves must be very thin and
concentric with the cylinder. A displacer seal riding on the sleeve assures
gas flow through the regenerator. Concentricity is a problem with this
regenerator also.
The tubular regenerator is composed of small-diameter thin wall tubes
constrained (usually in a single layer) between the cylinder (or displacer)
and a cylindrical liner with the tube axis in the direction of the cylinder
axis. The gas flows in the axial direction either through the tubes or in the
triangular-shaped spaces between the tubes and the liner (or cylinder). This
type of hot regenerator is easy to fabricate and is being used in several
VM refrigerators. The screen and ball regenerators are similar to the cold
regenerators except matrix elements are larger and of materials such as
monel and stainless steel.
Heat Rejection
Heat rejection has been accomplished in several ways. These include
rejection to forced ambient air (Globe et at., 1971), rejection to a pumped
liquid which in turn rejects the heat to air or a radiator (Doody, 1971),
and rejection by heat pipes (Yoshikawa, 1970; Browning and Potter, 1972;
AiResearch Mfg. Co., 1972, 197 4; Browning et at., 1972; Miller and Potter,
1973, 1974).
SPLIT-VUILLEUMIER CRYOCOOLER
iii. Shock and vibration resistance ratings higher than any other
existing coolers;
iv. Split configuration for remote cooling;
v. Low operating cost over a long in-service life;
vi. Simplified maintenance and service requirements.
Split coolers are attractive for applications where minimum vibration and
mass are important as in gimbaled, swiveling infrared sensors, or missile
sights. The compressor and expander units are coupled by a thin flexible
tube permitting independent motion of the expander cylinder. The void
volume of the connecting tube increases the dead volume and thereby
Vuilleumier Cryocoolers 233
Table 4.4. Specifications for CTI Model No. VM-1 Split-Vuilleumier Cooling
Enginesa
REFERENCES
Block, R. F. (1973). "The Design, Development and Test of a Heat Pipe Interface for the
Hot Cylinder of a V-M Refrigerator." Closed Cycle Cryogenic Cooler Technology and
Applications, pp. 271-288, AFFDL-TR-73-149, Vol. 1, WPAFB, Dayton, Ohio.
Browning, C. W., Miller, W. S., and Potter, V. L. (1972). "75 K Vuilleumier Cryogenic
Refrigerator Task IV, Model Fabrication." Report No. 72-8687, NASA GSFC,
Greenbelt, Maryland (AiResearch Manufacturing Co., Torrance, California).
Browning, C. W., and Potter, V. L. (1972). "75 K Vuilleumier Cryogenic Refrigerator. Task
II, Analytical and Test Program," Report No. 72-8497, NASA GSFC Greenbelt,
Maryland (AiResearch Manufacturing Co., Torrance, California).
Browning, C. W., and Potter, V. L. (1972). "75 K Vuilleumier Cryogenic Refrigerator, Final
Report for Task II-Analytical and Test Program." Contract NAS 5-21096, Goddard
Space Flight Center, Greenbelt, Maryland, August.
Browning, C. W., Miller, W. S., and Potter, V. L. (1972). "75 K Vuilleumier Cryogenic
Refrigerator-Final Report for Task IV Model Fabrication." Contract NAS 5-21096,
Goddard Space Flight Center, Greenbelt, Maryland, November.
Browning, C. W., Potter, V. L., Miller, W. S., and Gasser, M. G. (1973). "Developments
Toward Achievements of Long-Life Cryogenic Vuilleumier Refrigeration Systems."
Closed Cycle Cryogenic Cooler Technology and Applications, pp. 225-244, AFFDL-TR-
73-149, Vol. 1, WPAFB, Dayton, Ohio.
Bush, V. (1938). "Apparatus for Heat Transfer." U.S. Patent No.2, 127,286, August 16.
Bush, V. (1939). "Apparatus for Compressing Gases." U.S. Patent No. 2,157,229, May 9.
Chellis, F., and Hogan, W. H. (1964). "A Liquid-Nitrogen-Operated Refrigerator for
Temperatures Below 77 K." Adv. Cryog. Eng. 9, 545-551.
Chellis, F. (1973). "The Future of the Stirling Cycle Refrigerator in Airborne 1-R
Applications." Closed Cycle Cryogenic Cooler Technology and Applications, pp. 11-118,
AFFDL-TR-73-149, Vol. 1, WPAFB, Dayton, Ohio.
Chellis, F. (1977). "Design Compromises in the Selection of Closed Cycle Cryocoolers."
Proc. Appl. of Closed Cycle Cryocoolers to Small Superconducting Devices, Boulder,
Colorado, pp. 109-122, October (NBS Spec. Pub. 508).
Cohen, B., and Daniels, A. (1973). "Design and Development of a Free Displacer/Thermal
Compressor Cryogenic Refrigeration System." Final Report, Contract No. DAAK 02-71-
C-0421, N.V.L. USAECOM Fort Belvoir, Virginia, AD 920028 (Philips Laboratories).
Cowans, K. K. (1968). "Heat Powered Engine." U.S. Patent No. 3,379,026, April23.
Crouthamel, M. S., and Shelpuk, B. (1973). "A Combustion-Heated Thermally Activated
Vuilleumier Refrigerator." Adv. Cryog. Eng. 18, 339-351.
Daniels, A., and du Pre, F. K. (1971). "Miniature Refrigerators for Electronic Devices."
Philips Tech. Rev. 32(2), 49-56.
Doody, R. D. (1971). "Two-Stage Vuilleumier Cycle Cryogenic Refrigerator System for
Advanced Forward Looking Infrared (AFLIR) Applications." AFFDL-TR-71-17, DDC
Number AD886822, August.
Doody, R. D. (1975). "Long Life High Capacity Vuilleumier Refrigerator for Space Applica-
tions." AFFDL-TR-75-108, WPAFB, Dayton, Ohio (Hughes Aircraft Co., Culver City,
California).
Gasser, M.G., Yoshikawa, D. K., and Browning, C. W. (1972). "An Approach to Long-Life
VM Cryogenic Refrigerators for Space Applications." Appl. Cryog. Techno/. 4, 416-434.
Globe, J. B., Rhia, F. J. III, Gainey, R. T., Aske, H. D., and Doody, R. D. (1971). "Vuilleumier
Cycle Cryogenic Refrigeration System for Infrared Scanner Applications." AFFDL-TR-
71-18, WPAFB, Dayton, Ohio (Hughes Aircraft Co.).
Harkless, L. B. (1973). "Reliability Test Results on V-M Coolers." Closed Cycle Cryogenic
Cooler Technology and Applications," pp. 93-110, AFFDL-TR-73-149, Vol. 1, WPAFB,
Dayton, Ohio.
Hogan, W. H. (1964). "Closed-Cycle Cryogenic Refrigerator and Apparatus Embodying
Same." U.S. Patent No. 3,151,446, Oct. 6.
Horn, S. B., Cowan, K. C., and Berry, R. L. (1973). "Low Production Cost VM Coolers."
Closed Cycle Cryogenic Cooler Technology and Applications, pp. 73-82, AFFDL-TR-73-
149, Vol. 1, WPAFB, Dayton, Ohio.
Vuilleumier Cryocoolers 235
Horn, S. B., Lumpkin, M. L., Walters, B. T., and Acord, T. T. (1973). "Miniature Cryogenic
Cooler for TOW Night Sight." Closed Cycle Cryogenic Cooler Technology and Applica-
tions, pp. 52-70, AFFDL-TR-73-149, Vol. 1.
Lehrfeld, D., and Pitcher, G. K. (1973). "An Oil-Lubricated Triple Expansion VM Cooler
for Long Duration Space Missions." Closed Cycle Cryogenic Cooler Technology and
Applications, pp. 189-210, AFFDL-TR-73-149, Vol. 1, WPAFB, Dayton, Ohio.
Leo, B. (1970). "Designer's Handbook for Spaceborne Two-Stage Vuilleumier Cryogenic
Refrigerators." AFFDL-TR-70-54, Dayton, Ohio (Hughes Aircraft Co., Culver City,
California).
Leo, B. (1971). "Vuilleumier Cycle Cryogenic Refrigeration System Technology Report."
AFFDL-TR-71-85, WPAFB, Dayton, Ohio (Hughes Aircraft Co., Culver City,
California).
Magee, F. N., and Doering, R. D. (1968). "Vuilleumier Cycle Cryogenic Refrigerator
Development." AFFDL-TR-68-2, WPAFB (Hughes Aircraft Co., Ltd.).
Miller, W. S., and Potter, V. L. (1973). "Fractional Watt Vuilleumier Cryogenic Refrigerator-
Final Report for Task I Preliminary Design." Report No. 72-8846, NASA GSFC,
Greenbelt, Maryland (AiResearch Manufacturing Co., Torrance, California).
Miller, W. S., and Potter, V. L. (1974). "Fractional Watt Vuilleumier Cryogenic Refrigerator-
Final Report for Task II Analytical and Test Programs." Contract N AS 5-21715, Goddard
Space Flight Center, Greenbelt, Maryland, March.
Morgan, N. E. (1971). "Analysis and Preliminary Design of Airborne Air Liquefiers."
AFFDL-TR-71-171, WPAFB, Dayton, Ohio, AD 892500 (Hughes Aircraft Co., Culver
City, California).
Morgan, N. E. (1975). "Preliminary Design of a Cryocooler for Satellite Sensors (CCSS)."
AFFDL-TR-75-154, Vol. 1, WPAFB, Dayton, Ohio (Hughes Aircraft Co., Culver City,
California).
Pitcher, G. K. (1970). "Design and Development of Laboratory Model Vuilleumier
Cycle Cryogenic Refrigerator." Report No. 243, DA Prof. No. 663704P05201,
Contract No. DAAK02-69-C-0364, N.V.L. USAECOM, Fort Belvoir, Virginia (Philips
Laboratories).
Pitcher, G. K., aod du Pre, F. K. (1970). "Miniature Vuilleumier Cycle Refrigerator." Adv.
Cryog. Eng. 15, 447-451.
Pitcher, G. K. (1971). "Development of Spacecraft Vuilleumier Cryogenic Refrigerators-
Part !-Thermodynamic Design and Study of Life Limiting Components." AFFDL-TR-
71-147, Part 1, WPAFB (Philips Laboratories).
Pitcher, G. K. (1972). "Development of Spacecraft Vuilleumier Cryogenic Refrigerators Part
11-Development and Fabrication." AFFDL-TR-71-147, Part II, WPAFB (Philips
Laboratories).
Pitcher, G. K. (1973). "Mechanical Life of Space Cryocoolers." Closed Cycle Cryogenic Cooler
Technology and Applications," pp. 211-224, AFFDL-TR-73-149, Vol. 1, WPAFB,
Dayton, Ohio.
Pitcher, G. K. (1975). "Spacecraft Vuilleumier Cryogenic Refrigerator Development-Final
Report." AFFDL-TR-75-114, WPAFB (Philips Laboratories).
Pitcher, G. K. (1977). "Development of a Flight Design for an Oil Lubricated VM
Cryocooler." AFFDL-TR-77-95, WPAFB (Philips Laboratories).
Prast, G. (1963). "A Philips Gas Refrigerating Machine for 20 K." Cryogenics 3, 156-160,
September.
Renyer, B. L. (1973). "High Capacity Long-Life Vuilleumier Cycle Refrigerator." Closed
Cycle Cryogenic Cooler Technology and Applications, pp. 14 7-162, AFFDL-TR-73-149,
Vol. 1, WPAFB, Dayton, Ohio.
Rhia, F. J., Ill. (1971). "Development of Long Life High Capacity Vuilleumier Refrigeration
System for Space Applications." AFFDL-TR-71-92: Part I-AD 887399L. Part 11-
"Description of Preliminary Refrigerator Design and Results of Thermal Analyses and
Life Studies," AD 890601L. Part 111-"Refrigerator Design and Thermal Analyses,"
AD 901235L. Part IV-"Design and Mechanical Analyses of Refrigerator and Interface
Unit," AD 920 255.
236 Chapter 4
Richter, R., and Mahefkey, E. T. (1973). "The Applicability of a Solar Collector Thermal
Power System to a Vuilleumier Cooler." Closed Cycle Cryogenic Cooler Technology and
Applications, pp. 245-270, AFFDL-TR-73-149, Vol. 1, WPAFB, Dayton, Ohio.
Rule, T. T., and Qvale, E. B. (1969). "Steady State Operation of the Idealized Vuilleumier
Refrigerator." Adv. Cryog. Eng. 14, 343-352.
Russo, S. C. (1976). "Study of a Vuilleumier Cycle Cryogenic Refrigerator for Detector
Cooling in the Limb Scanning Infrared Radiometer." NASA CR 145078 AD/N77-1211
(Hughes Aircraft Co., Culver City, California).
Sherman, A. (1971). "Mathematical Analysis of a Vuilleumier Refrigerator." NASA TM-X-
65534, Ace. No. N-71-25812.
Sherman, A. (1972). "Selected Vuilleumier Refrigerator Performance Characteristics." Adv.
Cryog. Eng. 18, 352.
Sherman, A. (1978). "Cryogenic Cooling for Spacecraft Sensors, Instruments, and Experi-
ments." Astronaut. Aeronaut., November.
Sherman, A., Gasser, M., Goldowsky, M., Benson, G., and McCormick, J. (1979). "Progress
on the Development of a 3-5 year Lifetime Stirling Cycle Refrigerator for Space."
Cryogenic Eng. Conf., Madison, Wisconsin, August.
Taconis, K. W. (1951). "Process of and Apparatus for Heat Pumping." U.S. Patent No.
2,567 ,454, September 11.
Vuilleumier, Rudolph. (1918). "Method and Apparatus for Inducing Heat Changes." U.S.
Patent No. 1,275,507, August 13.
White, R. (1973). "Program for Predicting V-M Cooler Off-Design Performance. Closed
Cycle Cryogenic Cooler Technology and Applications," pp. 163-188, AFFDL-TR-73-
149, Vol. 1, WPAFB, Dayton, Ohio.
White, R. (1976). "Computer Program for Optimizing Three-Stage Vuilleumier Cycle
Cryogenic Refrigerators." AFFDL-TR-76-27, WPAFB, Dayton, Ohio.
White, R. (1976). "Vuilleumier Cycle Cryogenic Refrigeration." AFFDL-TR-76-17,
WPAFB, Dayton, Ohio, AD/A-027-055.
Yendall, E. F. (1960). "A Novel Refrigerating Machine." Adv. Cryog. Eng. 2, 188-196.
Yoshikawa, D. K. (1970). "75 K Miniature Vuilleumier Cryogenic Engine-Task 1." Report
No. 70-6854, NASA GSFC, Greenbelt, Maryland (AiResearch Manufacturing Co.,
Torrance, California).
Chapter 5
INTRODUCTION
motion. It is not known if actual engines were constructed and put into
service.*
The reinvention of the Solvay and Postle engines appears to have
occurred simultaneously. Some development of Solvay engines occurred
at A. D. Little Inc. but effort was concentrated principally on the Postle
displacer concept. Successful prototypes were developed and put into
production. They are widely known as Gifford-McMahon cryocoolers.
No Solvay engines are known to be in production anywhere. Cooling
engines made by Air Products and Chemicals Inc. of Allentown, Pennsyl-
vania are described as "modified Solvay" engines but have displacers rather
than expansion engines and are, therefore, Gifford-McMahon engines. The
lack of development effort given to Solvay engines is hard to understand.
It presumably arises from the difficult problem of a dynamic seal on the
piston to contain the working fluid in the system. The same problem exists
with Stirling and Claude cryocoolers and the seal technology developed
there could be applied to Solvay machines. Possibilities include the Philips
rolling seal, welded metal bellows, or the hydrodynamic lubricated, close
tolerance pistons used in Kaptiza, Collins, and Doll-Eder Claude cycle
expanders. There are many variants of Solvay machines. It is likely these
will receive development attention with the increasing need for intermedi-
ate-capacity cooling systems as superconducting applications materialize.
The Gifford-McMahon cryocooler was conceived at A. D. Little in
the mid-1950s and rapidly became a commercial product. The engine is
characterized as a large, heavy, slow-running system requiring minimal
maintenance. Higa et al. (1978) recently described an installation of 30
Gifford-McMahon cryocoolers with a cumulative total of one million hours
of operation over the preceding decade.
The durability of the Gifford-McMahon cryocooler makes it the system
of choice where size, weight, and efficiency are not of prime importance.
The coolers come in a variety of capacities and configurations. Table 5.1
summarizes the technical specification of Gifford-McMahon cryocoolers
available from Cryomech Inc., Jamesville, New York.
Initial production versions of Gifford-McMahon cryocoolers were
introduced by A. D. Little Inc. in the early 1960s and have remained in
production there and at the successor company Cryogenic Technology,
Inc., Waltham, Massachusetts. Figure 5.1 is a general view of a CTI Model
No. 0120 Gifford-McMahon cryocooler. This is a two-stage machine cap-
able of generating 1 W of refrigeration at 26 K.
The competing range of commercial cryocoolers was introduced by
Cryomech, Inc. of Jamesville, New York, a company founded by William
* The author would appreciate any information readers could supply on this point.
C'l
Q
=
I
.=-
,.,==
Table 5.1. Technical Specifications for Gifford-McMahon Cryocoolers"
10
==
::r
Q
Type of refrigerator l'l
riJ
Q
AL01 AL02 D101 GB02 GB05 GB12 ~
10
~
Output 10
" All require no cooling, are Gifford-McMahon cycles with heat exchangers, operate at 144 rpm, and are rated for 3000 hr endurance. Courtesy Cryomech, Inc.,
Jamesville, New York.
N
~
-=
240 Chapter 5
Fig. 5.1. Gifford-McMahon cryocooler CTI Model No. 0120 (courtesy Cryogenic
Technology Inc.).
COMPRESSOR
OISPLACER
OUTLET
VALVE
(sliding) gas seal is provided on the piston to stop leakage of gas from the
cylinder. The seal is located near the top of the piston and operates at
ambient temperatures. Location of the seal at the bottom of the piston
would reduce the available refrigeration generated because of thermal
degradation of the frictional work of the seal.
The regenerative heat exchanger consists of a porous matrix of finely
divided material, usually metallic, contained in a cylinder of low-conduc-
tivity material. The function of the regenerator is to act as a thermal sponge
abstracting heat from the gas when it is passing downwards through the
matrix and releasing heat to the gas when it is passing upwards through
the matrix.
Inlet and outlet valves control the flow of gas into and out of the
expansion engine system. The valves are coupled to the crankshaft so their
operation is synchronized with the piston position. The inlet and outlet
valves are connected through a loop containing a gas compressor and
adequate low-pressure and high-pressure reservoirs. A water cooler is
included in the high-pressure side delivery duct of the compressor to remove
the heat of compression.
The operation of the engine is as follows:
(1) Pressure Buildup. With the piston at the bottom of its stroke, the
inlet valve opens and high-pressure fluid enters, compressing the fluid
already in the engine from pressure p 1 to pressure p 2 •
(2) Intake Stroke. With inlet valve open, the piston is raised from the
bottom dead center position. Gas is admitted through the inlet valve and
regenerator into the bottom of the cylinder. Passing through the regen-
erator, the gas cools by heat transfer to the matrix material of the
242 Chapter 5
the mass of gas in the cylinder increases as fluid flows to the cylinder from
the regenerator. In Fig. 5.3, the expansion process is arbitrarily represented
by a straight line. The shaded area represents additional refrigeration
attributable to expansion in the cylinder of gas in the regenerator. This
partly compensates for the gas used to fill the regenerator void space which
appears to expand in the system to no useful purpose.
The gas required to fill the regenerator is not so great as might be
expected from a simple consideration of the void volume. The temperature
of the regenerator varies from the low refrigeration temperature at the
bottom end to ambient temperature at the top. The effective average value
of temperature may, therefore, be twice that of the expansion cylinder so
the mass required to fill the matrix is only half that required if it were all
at the expansion space temperature.
For best performance, the piston should have a discontinuous motion
with brief pauses at the end of each stroke. However, Gifford has indicated
that it is possible to use harmonic piston motion for simplicity without
significant loss of performance. Other constructional features are that the
piston may be quite loose fitting in the cylinder with up to 0.5 mm (0.02 in.)
clearance except for the gas tight seal at the warm end. To minimize heat
conduction, the cylinder must be made of thin-wall stainless steel tube and
the piston of low-conductivity, plastic material. The regenerator might be
a stack of dense mesh wire screen of 0.03-mm (0.001-in.) diameter copper
or phosphor bronze wire having 80 to 160 strands per em (200 to 400
strands per in.) contained in a thin-wall stainless steel tube.
A miniature Solvay cooling engine described by Gifford et al. (1960)
is shown in Fig. 5.4. The engine, intended for infrared detector applications,
was capable of attaining 55 Kin a single stage. The cylinder was a stainless
steel tube 5 mm (0.2 in.) in diameter with a wall thickness of 0.13 mm
(0.005 in.). The piston, made of linear base laminated phenolic, was desig-
ned to incorporate the regenerator consisting of fine wire mesh screen.
The inlet and outlet pressure of the helium working fluid was 1725 and
345 kPa (250 and 50 psi) with a consumption of 6 to 9 x 10-4 m 3 /sec
(0.2 = 0.3 fe /min).
The unit was equipped with a subassembly at ambient temperature
containing the inlet and exhaust valves, the piston seal, and cam shafts for
the valve and piston actuation. A 400-Hz, 115-V motor was used to drive
the crankshaft through a reduction gear at 500 rpm. The overall weight of
the cooling unit was 0.14 kg (10 oz) including the timing motor but excluding
the external gas compressor or compressed gas bottles.
The cooling characteristic for the device is shown in Fig. 5.5. It was
possible to cool to operating temperatures within 2-3 min and to maintain
temperatures in the 50-80 K range, indefinitely.
244 Chapter 5
o· 0
0
0 0
::
Fig. 5.4. Experimental miniature Solvay
cryocooler (after Gifford eta/., 1960).
130
120
t ItO
~ \
\
~ 100
a:
"'
0.
l
... 90
0
z
"'9
\\
80
8
70
~
60
........... .....__ Fig. 5.5 . Cooldown characteristic of the experi-
50
0 4 mental miniature Solvay cryocooler (after Gifford
ELAPSED RUNNING TIME - (mon ) eta/., 1960).
Gifford-McMahon, Solvay, and Poslle Cryocoolers 245
At (c) the inlet valve closes quickly and the exhaust valve opens slowly
and progressively. The pressure falls as fluid leaves the engine through the
exhaust valve. Fluid in volume 2 expands .as the pressure decreases and
the temperature of gas in volume 2 decreases. Heat transfer from the
cylinder walls increases the gas temperature to (d) in Fig. 5.7. Finally, the
displacer moves from the top to bottom of the cylinder causing fluid to
return back through the regenerator to volume 1 or exhaust. In passing
through the regenerator, the temperature of the gas increases to a value
slightly less than the mixture, following admission of high pressure gas. The
gas leaving the exhaust is, therefore, warmer than the intake feed gas. Heat
removed in this way is equal to the total refrigeration achieved. The net
result of investing the compressor work is that heat is pumped from a low
temperature to above the ambient level. The useful product is the refriger-
ation effect generated during expansion.
The work done by the gas on the displacer is shown by the work
diagram for volume 2. The work done by the displacer on the gas is shown
by the work diagram for volume 1. The pressures are the same in both
spaces so the areas of the work diagrams are identical (neglecting the
displacer rod area). For an ideal gas and adiabatic conditions, the area of
the diagram for volume 2 is equal to the change in enthalpy of the gas in
volume 2. It represents the net refrigeration available. Similarly, the area
of the diagram for volume 1 represents the change in enthalpy of a perfect
gas under adiabatic conditions.
REGENERATOR
6T
...
"'
~ 200
."'......"'
<(
1- 100
w
"'::>
!;;( 200
~
31-
and loosely fitting dis placer I cylinder sets with little mechanical stress. The
problems of construction are minimal. The additional stages operate exactly
like the single-stage device with gas supplied at the temperature of the
stage above. The expansion spaces describe P- V diagrams representing
the useful refrigerating output of the stage and the regenerative heat
exchangers prevent the loss of refrigeration to the stage above.
Multistaging is attractive because it provides a thermodynamically
efficient process to attain cryogenic refrigeration temperatures. For a given
refrigeration requirement, there is a decreased power requirement.
Expansion over given pressure and volume ratios produces the same
refrigeration regardless of temperature. However, the mass of gas involved
is inversely proportional to the temperature. A unit of refrigeration at 20 K
requires four times the gas and consequently four times the power, as at
80 K. Most cryocooler applications require the removal of heat at different
temperature levels, such as precooling before liquefaction, or the use of
cooled radiation shields to reduce heat leaks to the system at the lowest
temperatures. It is, therefore, thermodynamically advantageous for
refrigeration to be available at different temperature levels.
Furthermore, multistaging relieves the operating regimen of regen-
erators and heat exchangers. Consider the "refrigeration loss" Q,:
Q, =MCp fiT, (5.1)
where M is the mass flow of gas, CP is the heat capacity of the gas, and
fiT, is the temperature difference between the gas entering and leaving the
250 ChapterS
Fluidic-Driven Displacer
A refinement to the original Gifford-McMahon cryocooler occurred
with application of a technique known as "cryomatic gas balancing." The
displacer mechanical drive was eliminated and a free displacer used instead.
Movement of the displacer was caused by differential fluidic (gas pressure)
forces acting on the displacer arising from the inlet and exhaust of gas
through a simple rotary cored valve.
A schematic diagram of a single-stage Gifford-McMahon cryocooler
with fluidic-driven free displacer is shown in Fig. 5.12. It consists of two
connected coaxial cylinders of different length and diameter containing the
free displacer. The displacer is about three quarters the length of the larger
cylinder and has a coaxial extension to fit the smaller cylinder. This creates
three chambers of variable volume depending on the position of the
displacer. Chamber 1 is the volume between the top of the displacer
extension and the small cylinder head. Chamber 2 is the volume between
the top of the displacer and the large cylinder. Chamber 3 is the volume
between the bottom of the displacer and the large cylinder. Gas seals are
provided on both the displacer and displacer extension to prevent the
leakage of gas past the displacer from one chamber to another. Chambers
2 and 3 are interconnected by a thermal regenerator of finely divided
metallic material. The pressure in both is, therefore, the same, apart from
small differences due to the pressure drop in the regenerator matrix when
252 ChapterS
ROTARY VALVE
RECIPROCATI G
OISPLACER
gas is flowing through. Gas is admitted and exhausted from chamber 1 and
the regenerator conduit through a rotary cored valve. The same valve
regulates both the inlet of high-pressure gas from the compressor and the
exhaust of low-pressure gas.
The sequence of operations illustrated in Fig. 5.13 is as follows:
(1) Pressure Buildup. At the beginning of the cycle the displacer is at
the bottom of the cylinder so the volume of chamber 3 is zero and the
volumes of chambers 1 and 2 are at their maximum value. The rotary valve
is in such a position that chamber 1 is full of gas at a high pressure and
chambers 2 and 3 are connected to exhaust and at low pressure. This is
represented in Fig. 5.13a. When the valve has turned to connect chambers
2 and 3 to the inlet port the pressure increases as represented in diagram
(b). Fluid enters through the inlet port to fill chamber 2, the regenerator
and the clearance space in chamber 3.
(2) Intake Stroke. Further rotation of the valve opens chamber 1 to
exhaust and the pressure decreases. The upward force on the displacer due
to pressure in chamber 3 now exceeds the downward force due to pressure
in chambers 1 and 2, by an amount equal to l:l.p times the area of chamber
1. As a result, the displacer moves to the top of the cylindei (Fig. 5.12c),
displacing fluid from chamber 2 through the regenerator to chamber 3. In
passing through the regenerator, the fluid cools and additional fluid enters
to fill chamber 3 and the regenerator.
Gifford-McMahon, Solvay, and Postle Cryocoolers 253
ROTARY VALVE
ll
L VOLUME, CHAMBER 3
Fig. 5.13. Cycle of operation of the Gifford-McMahon cryocooler with free , fluidic-driven
displacer.
in enthalpy and, for adiabatic conditions with an ideal gas, the refrigeration
effect of the process per cycle. The diagram for chamber 2 is virtually
identical but reversed in direction to that for chamber 3. A temperature-
time history of gas in the engine would be the same as that presented in
Fig. 5.8.
In a description of the fluidic-drive displacer engine, Gifford (1966)
does not discuss valve operation. It is virtually certain that the same careful
attention to programmed valve opening would be required as with a
mechanically driven displacer. This gains improved effectiveness of the
regenerator resulting from steady flow during the pressure buildup and
expansion processes. A good deal of expertise is concealed by Gifford's
disarming comment:
... timing between pressurizing and depressurizing the chambers may be set
by proper selection of the components connecting the chambers with the
valve .. .
Fluidic-driven displacers have been used in multistage units, where
refrigeration is generated at different temperature levels. This requires an
additional displacer /cylinder unit and regenerator per stage. A schematic
digram of a two-stage cryomatic gas balancing refrigerator is shown in Fig.
5.14. An experimental unit of this type mentioned by Gifford (1966) had
..
0
>=
u UP
."'"'
u
~DOWN
90 180 270 360
c;
IN
seal of the Solvay engine has to withstand the full pressure difference
between the cylinder and the "crankcase." In the Gifford-McMahon
engine, the displacer seal has to sustain only the difference due to pressure
drop in the regenerator. Moreover, the difference reverses direction and
tends to cancel any leakage. A moderate leakage can be tolerated. The
displacer drive system is not required to sustain the large forces of a piston
drive and may be structurally simple and lightweight with fewer vibration
problems.
The principal disadvantage of the Gifford-McMahon cycle is the high
work input per unit of refrigeration (W/Q) characteristic. This arises from
the incomplete expansion and the mass of compressed gas consumed in
the regenerator void volume with little or no return of refrigeration.
The ideal work per unit of refrigeration for various cycles is shown in
Fig. 5.16. This was prepared by Gifford (1966) for refrigeration at 80 K
with helium as the working fluid, heat rejection at 300 K and a range of
pressure ratios p 2 /p 1 up to 5. The W/Q characteristic is drawn for the
Gifford-McMahon cycle, the Joule-Brayton cycle, and the ideal Stirling
and Carnot cycles. In discussing the figure, Gifford suggests that because
of the speed of operation, the reference cycle for Stirling cryocoolers should
be the Brayton cycle, the ideal reference cycle for small reciprocating and
turbine expanders. While allowing that the Gifford-McMahon system does
require more work for a given level of refrigeration, he suggests the
difference is not great at the small compression ratios of 1.5 to 2.5 used
in practice. Furthermore, the practical advantages outweigh the additional
work. The principal advantage is high reliability because the units are
simple, slow moving, and unstressed. The gas seal on a typical machine
100
90
/- GIFFORD McMAHON
80 / CYCLE
0 / .Y
BRAYTON CYCLE
(NO WORK RECOVERY)
/
04 v .---
~
0
...- BRAYTON CYCLE
4
oy /
0
TIRLING 8 CARNOT
0
may move at the rate of 0.05 to 0.15 m/s (2 to 5 in./s) compared with the
0.75 to 1.25 m/s (30 to 50 in./s) of small Stirling cryocoolers. Continuous
operating times of 5 to 10,000 hr are not uncommon and Gifford's cryo-
coolers are noted for their high reliability and extended maintenance-free
operation.
Heat-Balance Analysis
Ackerman and Gifford (1971) presented data on the experimental
performance of a Gifford-McMahon engine with a fluidic-drive free dis-
placer. The experimental refrigerator is shown in Fig. 5.17. The machine
was a single stage expansion engine with a displacer 19 mm (0.75 in.) in
diameter, 0.1 m (4 in.) long and with a stroke of 12 mm (0.5 in.). With
helium as the working fluid and pressures ranging from 689 to 1960 kPa
(100 to 285 psi) it was found that the best operating speed was 3Hz
(190 cycles/min).
CAPACITANCE
DISPLACEMEN
TRANSDUCER FLOW REVERSING
~ VALVE
~ • SEALS
~
c
...z
~
...
:ll
~
300f-----
:::.---- =---~
1
( 250osi~ ;;e--T•50°K
W• 5.8wath
INPUT CONDITIONS;
..,
z
Ci
\ P"• 3DO polg
PL • 9!5 pliQ
RPM•I90
"'"'
0::
\ wideal· 14. 5 wath
.,.,"':::>
0::
"'0::
Q.
0::
"'<.>
.,:::>z
0
:>-
1112.5 poig
1001------ -======-=-===~_jiOO paig
Fig. 5.18. Comparison of ideal and
1----- .SOOin.--------+1 actual work diagrams for the experi-
mental Gifford-McMahon cryocooler
DISPLACEMENT (after Ackermann et al., 1971).
Gifford-McMahon, Solvay, and Postle Cryocoolers 261
losses were found to be small and were ignored. No account was taken of
radiation losses.
Much of the large difference, 8 W, between the ideal and indicated
diagram refrigeration was attributed to the large void volume in the cold
region necessary to accommodate the instrumentation. In commercial
machines, the regenerative heat exchangers are frequently included in the
displacer to reduce the cold end voids and conduction heat loss because
the regenerator tube is eliminated. Motional heat leak then remains the
prime thermal loss.
DOUBLE -ACTING
- - - - - - - COMPRESSOR -CYLI DER___ -
COOLI G'
WATER
REGENERATIVE
DUCTS
from the valve seat moving inwards with the displacer. With the inlet valve
open, high-pressure fluid filled the cylinder and was cooled by the cooling
water heat exchanger on the cylinder end plate.
As the compressor piston reached the end of its stroke and started to
return, the pressure began to fall and fluid flowed from the dis placer cylinder
towards the compressor. The flow acting on the valve caused it to move
to the right, thereby closing the exhaust valve. Now the pressure in the
cylinder exceeded the decreasing pressure in the valve chamber, thereby
inducing the displacer to move from the inward to the outward position.
In this movement, fluid was transferred through the regenerative ducts to
the central expansion space.
Gifford-McMahon, Solvay, and Postle Cryocoolers 263
As the displacer approached the end of its outward stroke, the valve
stem impacted the cylinder cover plate causing the valve to separate from
the exhaust valve seat. Now with the exhaust valve open, fluid passed back
to the compressor cylinder all the while expanding in the expansion space
and generating refrigeration. When the compressor piston reached the end
of its expansion stroke and started to return, the flow of fluid reversed and
the operation of the system was repeated as described above.
The novel features of the self-acting valve arrangement of the Postle
engine endow the concept with the potential for development as a regenera-
tive cryocooler in the future. An alternative could be equipped with a gas
spring to operate at resonant frequency and be supplied with constant
pressure fluid from a reservoir.
REFERENCES
Ackermann, R. A., and Gifford, W. E. (1971). "A Heat Balance Analysis of a Gifford-
McMahon Cryorefrigerator." Adv. Cryog. Eng. 16, 221-29.
Collins, S.C., and Cannaday, R. L. (1958). Expansion Machines for Low Temperature Processes.
Oxford University Press, Oxford, England.
Gifford, W. E., and McMahon, H. 0. (1960). "A New Low Temperature Gas Expansion
Cycle-Part II." Adv. Cryog. Eng. 5, 368-372.
Gifford, W. E., and Hoffman, T. E. (1961). "A New Refrigeration System for 4.2 K." Adv.
Cryog. Eng. 6, 82-94.
Gifford, W. E. (1966). "The Gifford-McMahon Cycle." Adv. Cryog. Eng. 11, 152-159.
Higa, W. H., and Wiebe, E. (1978). "One Million Hours at 4.5 Kelvin." Proc. Appl. of
Closed-Cycle Cryocoolers to Small Superconducting Devices, pp. 99-108, Boulder,
Colorado, October 1977 (Published as NBS Spec. Pub. 508).
McMahon, H. 0., and Gifford, W. E. (1960). "A New Low Temperature Gas Expansion
Cycle-Part 1." Adv. Cryog. Eng. 5, 354-366.
Solvay, E. (1887). Deutsches Reichspatent No. 39,280.
Thomas, W. H. (1873). British Patent No. 709, "Producing Cold for Preserving Animal Food,
etc." (Available from Science Museum Library, London).
Chapter 6
INTRODUCTION
Isenthalpic Expansion
Consider the system shown in Fig. 6.1. It consists simply of a duct
containing a porous plug, a valve that may be partly "cracked" open, or
a metal plate with very fine holes drilled through it. Gas at high pressure
enters at station 1 and expands through the constriction to a low pressure
at station 2. Fluid passing through the plug experiences high fluid friction
which degrades the pressure but no work is done. For the purpose of
illustration it is assumed the duct is thermally insulated so that no heat is
transferred to or from the fluid. The fluid flows so slowly that the change
in kinetic energy in expansion is negligible. The duct is horizontal with no
change in potential energy.
Applying the steady flow energy equation to the system shown in Fig.
6.1 and equating changes in kinetic energy, potential energy, work, and
heat transfer to zero leaves
H1=H2
The actual flow in the valve is thermodynamically irreversible and the fluid
is not homogeneous. The fluid downstream at 2 can be reasonably assumed
to be mixed and homogeneous. It will be found to have the same enthalpy
as the inlet fluid.
Outlet conditions of temperature and pressure following isenthalpic
expansion from given inlet conditions are plotted on the temperature-
pressure plane shown in Fig. 6.2. The curves obtained are lines of constant
enthalpy. Starting at the right-hand, high-pressure, side of the figure and
®
Fig. 6.1. Duct with porous plug for Joule-Thomson expansion of gas.
Joule-Thomson Cooling Systems 267
moving along the constant enthalpy curve, A-B, the temperature initially
increases as the pressure decreases. The increase in temperature reaches
maximum at point C. Further isenthalpic expansion causes a decrease in
temperature.
Maximum inversion
Gas temperature, (K)
Above ambient
Carbon Dioxide 1500
Oxygen 760
Argon 722
Nitrogen 623
Air 602
Below ambient
Neon 250
Hydrogen 202
Helium 40
oxygen, argon, nitrogen, air, etc., and (b) those for which the maximum
inversion temperature is below ambient temperature, i.e., neon, hydrogen,
and helium. This difference in maximum inversion temperature, either
above or below the ambient value, explains why Joule and Thomson
observed the curious phenomena that on expansion through a porous plug
hydrogen experienced an increase in temperature whereas all the other
gases experienced a decrease.
Gases with low values of the maximum inversion temperature must
be precooled, by an expansion engine or other means, to a temperature
within the inversion envelope to ensure a decrease in temperature on
isenthalpic expansion.
Joule-Thomson Coefficient
The change in temperature for an isenthalpic pressure drop is defined
by the Joule-Thomson coefficient for a fluid: f.J.H = (aTjap)H the change in
temperature due to a change in pressure at constant enthalpy. The Joule-
Thomson coefficient is the slope of the isenthalpic curves shown in Fig.
6.2. The coefficient is zero along the inversion curve (at the peak of the
isenthalps) since the slopes of the isenthalps are zero. To the right of the
inversion curve the Joule-Thomson coefficient is negative and the tem-
perature increases on expansion. To the left the Joule-Thomson coefficient
(slope of the isenthalp) is positive and the temperature decreases on
expansion. Figure 6.3 shows the Joule-Thomson coefficient for air as a
function of temperature at several different pressures (reproduced from
Daunt, 1956).
Joule-Thomson Cooling Systems 269
08
06
~~~,
I .,., ~
04
I ~~
~~
); .------
-..;::
p•200 otm
(6.1)
(6.2)
and
c p
=(aH)
aT p
(6.4)
and
(aH)
ap
- v - T(~)
T aT p
(6.5)
(6.6)
*The following is a summary of the excellent treatment of this topic given by Barron (1966).
270 Chapter6
(6.8)
(6.10)
The first term in Eq. (6.1 0) is indicative of the departure of a real gas from
perfect gas behavior in terms of Joule's law. This states that the internal
energy of a gas is a function only of temperature. For real gases the first
term is always negative and contributes to a decrease in temperature on
expansion. This arises because the molecules of the gas become increasingly
separated as expansion proceeds and their microscopic potential energy is
increased. No external work or heat transfer occurs, so the increase in
potential energy must be offset by a decrease in the microscopic kinetic
energy of the gas molecules and is manifested by a decrease in temperature.
The second term in Eq. (6.10) is indicative of the departure of a real
gas from perfect gas behavior in terms of Boyle's law. This states that the
product of pressure and volume of a perfect gas is a function only of
temperature and is illustrated in Fig. 6.4 by the horizontal line. This shows
the product of pressure and volume (an isothermal) to be independent of
the pressure.
The behavior of real gases is also illustrated in Fig. 6.4. At low pressures
and at temperatures near saturation conditions gases are more compressible
than Boyle's law predicts since attractive forces are acting to condense the
gas. This means the second term in Eq. (6.10) is negative and contributes
to the production of a temperature decrease. At high pressures the
molecules are closer together and repulsive forces are brought into action.
The gas is therefore less compressible than Boyle's law predicts, and the
Joule-Thomson Cooling Systems 271
Linde-Hampson Systems
A schematic diagram of an elementary Linde-Hampson gas liquefac-
tion system is shown in Fig. 6.5. The system consists of a compressor, a
recuperative heat exchanger, a Joule-Thomson valve, and a liquid/gas
separator and liquid reservoir.
The system operates as follows:
(a) Compression (process 1-2"). Gas enters the compressor at condi-
tion (1) and is compressed to state (2). Work W is supplied and heat Q
rejected during compression of the gas. For the purpose of analysis the
compression may be assumed isothermal or isentropic. In practice it would
be a polytropic process (1-2") between the isothermal process 1-2 and the
isentropic process 1-2'.
272 Chapter6
trate the make-up/mixing process (7) is shown substantially less than the
compressor inlet (1).
(g) Make-up and mixing (process 8/7-1). As liquid gas is produced
and leaves at condition (5), new make-up gas must be admitted at the
compressor inlet to maintain constant mass flow and continuous operation.
The make-up gas is admitted at ambient temperature, (8), and mixes with
the return effluent (7) for admission to the compressors, (1).
For the purpose of elementary analysis it is assumed that:
i. the flows take place with no pressure drops,
n. there are no external heat leaks,
iii. the recuperative heat exchanger is perfect, and
iv. the compression process is isothermal.
Applying the steady-flow energy equation for mass, m, and enthalpy, h,
to the recuperative heat exchanger, the expansion valve and the separator
we obtain
where
(6.11)
m2 h?-hs
and in the idealized system h 7 = h 1 • The fraction of gas liquefied therefore
depends only on the pressure and temperature at (1) and the pressure ratio
p 2 /Pl· Many systems for air liquefaction use ambient pressure and tem-
perature for condition (1). This fixes h 1 and h 5 so the only variable is the
isothermal compression ratio p 2 /p 1 which determines h 2 •
The maximum fraction of gas will be liquefied when the optimum
pressure p 2 is that at which the enthalpy h 2 is a minimum. Maxima and
minima occur when the rates of change or slope of a function are zero. In
this case h 2 is a minimum when the rate of enthalpy change as a function
of pressure (at constant temperature) is zero, i.e.,
(ah) _0 (6.12)
ap T=Tt
It can be seen from Eq. (6.1) and Eq. (6.4) that Eq. (6.12) is equivalent to:
(ah)
ap T=Tt
= o= (p.,H · Cph=T, (6.13)
274 Chapter6
T T
s
Ia) (b) (c)
Fig. 6.6. Start-up and cool-down of a simple Linde-Hampson system where the gas has a
maximum inversion temperature above ambient.
be fluorine (normal boiling point 85 K), oxygen (90 K), air (78 K), methane
(113 K), argon (87 K) or nitrogen (77 K). In practice for reasons of economy
and safety nitrogen is always used except in unusual circumstances.
Helium has a maximum inversion temperature of only 40 K. Therefore
the only liquids suitable for helium precooling are hydrogen (normal boiling
point 20 K) or neon (27 K). On grounds of safety neon is the preferred
fluid but economic considerations dictate that hydrogen is almost always
used.
Cascade System
A helium liquefier operating on the precooled Linde-Hampson cycle
consists of a series of separate systems having different fluids at different
minimum temperature levels. This is known as a "cascade system." It was
used originally for the liquefaction of air. The schematic arrangement of
the compound system is shown in Fig. 6.9. The useful product, helium, is
available at 4 K. The helium is precooled by a hydrogen bath at 20 K. The
hydrogen is precooled in a nitrogen bath at 77 K.
A system of this type was described by Mann et al. (1960). The liquefier
was 0.5 m diameter and 2m high and was constructed at the Cryogenic
Engineering Laboratory of the National Bureau of Standards, Boulder,
Colorado. Mann reports the following performance of the helium liquefier:
i. Liquid helium: 2.4 kg per hour
ii. Helium flow rate: 12.6 kg per hour
iii. Liquid yield (a/b)= 0.19
iv. Liquid hydrogen evaporated: 1.2 kg per hour
v. Liquid nitrogen evaporated: 10 kg per hour.
Geist and Lashmet (1960) described a lower capacity (750 mW useful
refrigeration at 4.4 K) Joule-Thomson cascade helium liquefier using
nitrogen and hydrogen precooler circuits. The system was constructed by
Air Products and Chemicals Inc., Allentown, Pennsylvania as a prototype
for application with masers, cryotrons, superconducting magnets, and other
electronic devices requiring refrigeration at liquid temperatures. The unit
was notable for the use of Corblin* diaphragm compressors to eliminate
the need for fluid purification equipment. The nitrogen and hydrogen
circuit used two-stage compressors having inlet, intermediate and discharge
pressures of 0.1, 1.96 (1.69 for Hz) and 24.8 MPa [15, 285 (245 for Hz),
and 3600 psi], respectively. The helium compressor operated over the
pressure range 0.1 and 1.86 MPa (15 and 270 psi).
NITROGEN COMPRESSOR
The advantages of precooling and the cascade system are not confined
to gases with subambient maximum inversion temperatures (neon, hydro-
gen, and helium) but apply equally to gases of the other group of gases to
improve the liquid yield fraction. Figure 6.10, reproduced from Barron
(1966), shows the liquid yield fraction (y) (the fraction of total mass
expanded that is saturated liquid after expansion) as a function of the
temperature of compression (assumed to be isothermal). The figure refers
to a simple idealized Linde-Hampson system using nitrogen. It is evident
that the amount of liquid produced is a strong function of the compressor
temperature: the lower the compressor temperature the higher the fraction
of gas converted to saturated liquid after expansion.
A similar improvement can be obtained using the precooled Linde-
Hampson system. An auxiliary refrigerating system cools the compressed
gas between the compressor aftercooler and the main recuperative heat
Joule-Thomson Cooling Systems 279
"'
;;:
5"' 0.2f----l--+--\--++------r----l
::::;
;;
(f)
lL O.lf---------,f-------'rt---'<------"tc---1
0
Fig. 6.10. Fraction of gas liquefied in a simple
Linde-Hampson nitrogen liquefier as a func- ~
tion of the gas temperature at the end of
isothermal compression at several different
~
0:: 0 100
lL
compression ratios. COMPRESSION TEMPERATURE (K)
exchanger. The auxiliary refrigerator can use the same working fluid but
this is not essential since the auxiliary circuit is entirely independent. It is
often advantageous to use alternative fluids with more favorable thermody-
namic or thermofluid characteristics.
Barron (1966) considered a precooled Linde-Hampson cycle for
nitrogen liquefaction with Freon 12 as the refrigerant in the auxiliary
precooling circuit. For his study he assumed the nitrogen system to operate
between pressures of 0.1 and 20 MPa (15 and 3000 psi) and temperatures
of 77 and 294 K (138 and 530 R). The Freon 12 system was assumed to
operate between pressure of 103 and 590 kPa (15 and 86 psi) and tem-
peratures of 243 and 469 K (438 and 844 R). He found the addition of
the Freon 12 precooler increased the liquid yield by 40% compared with
the simple Linde-Hampson nitrogen liquefaction system. At the same time
the total work requirement increased only marginally so the work require-
ment per unit mass liquefied decreased by almost 40% compared with the
simple cycle. This improvement results from the closer approximation to
an ideal reversible system. It can be carried further with multiple refriger-
ation systems with continuing improvement but at a progressively diminish-
ing rate.
The cascade system was first applied for the liquefaction of air. One
system suggested by Keesom (1933) utilized auxiliary refrigeration systems
(in order of decreasing temperatures) of (a) ammonia used to liquefy, (b)
ethylene used to liquefy, (c) methane used to liquefy nitrogen. The system
is known as the "outer cascade cycle."
A later development, known as the "incorporated cascade cycle" uses
a compound working fluid having several components which are mixed
and processed together in a single compressor. The mixture circulates
through a series of heat exchangers and separators so that various fractions
280 Chapter6
of the gas become liquid and are separated at different temperature levels.
The various fractions are further expanded and subcooled to provide
refrigeration to cool and condense different components of the mixture.
Detailed consideration of this complex system is beyond the level of
this book. Suffice to say the incorporated cascade system is used commer-
cially in large installations for the liquefaction of natural gas. It is also used
in smaller plants of low cooling capacity. The principal advantage is the
high flexibility. Depending on the choice of components and the composi-
tion of the gas the system can be used to achieve refrigeration at any
temperature level between ambient and cryogenic temperatures. Carbonell
et a/. (1972) have described a development of the incorporated cascade
system using a complex working fluid and regenerative rather than
recuperative heat exchangers. The regenerators are switched periodically
to reverse the direction of fluid flow.
E~TRDPY
LOW PRESSURE
SEI'jiRATOR ~----•1
LIQUID PRODUCT
Fig. 6.11. Schematic diagram of gaseous liquefier using the Linde dual-pressure system.
Fig. 6.12. Ratio of the work required per FRACTION OF GAS RECIRCULATED
unit mass liquefied in the Linde dual- AT INTERMEDIATE PRESSURE
Prl!cool,ng
stations
Figure of Merit
To provide the means of comparison of gas liquefaction systems the
figure of merit is sometimes used. This is defined as the theoretical minimum
work required to liquefy unit mass of gas compared with the actual work
requirement of a system: FOM = Wid/ Wac· The figure of merit can range
from zero to unity and indicates how closely a system approaches ideal
operation.
CHARGE
ENTROPY
the order of 104 MPa for air), and far exceeds the capacity of any compress-
ion equipment existing or in prospect. Nevertheless the process provides
a useful basis on which to compare theoretical and practical liquefaction
systems.
Using the above concept for gas liquefaction Barron (1966) calculated
the ideal minimum work requirement for different fluids. His results are
reproduced in Table 6.2.
Ideal work of
Normal boiling liquefaction
Gas point (K) (kJ/kg)
Table 6.3. Performance of the Linde-Hampson System Using Different Fluids P1 =0.1 MPa, T 1 = 294 K, p 2 = 20 MPa, T2 =
294 K; Heat-Exchanger Effectiveness, E =1.00; Compressor Overall Efficiency Ylc =100% a
n
=-
i
~
~
Joule-Thomson Cooling Systems 287
Table 6.4. Comparison of Liquefaction Systems Using Air as the Working Fluid
T 1 = 294 K and p 1 =Atmospheric Pressure, Ylc =Overall Compressor Efficiency,
E = Heat Exchanger Effectiveness, i = Fraction of Mass Flow Recirculated a
the figure of merit to 0.201. This compared favorably with the cascade
system having the best figure of merit, 0.221. The simple Linde-Hampson
system therefore consumed about three times the work per unit mass
liquefied of the thermodynamically superior but very much more compli-
cated cascade or precooled, dual-pressure Linde cycle.
Hymatic Coolers
Miniature Joule-Thomson cooling systems have been developed for
electronic applications, particularly infrared thermal detectors. The
development was carried out in Great Britain by the Hymatic Engineering
Co. Ltd., Redditch, Worcester with Government sponsorship through the
Royal Radar Establishment, Malvern, Worcester. Accounts of this work
have been given by Mcinroy (1967), Parkinson (1959, 1967), and Nicholds
(1968, 1970).
Mcinroy described the construction and gave the performance data for
both liquid air and liquid hydrogen/liquid nitrogen Joule-Thomson
miniature liquefiers. Figure 6.15 is a cross section of the miniature air
liquefier described by Mcinroy as follows:
The miniature liquefier is enclosed by a representative cell for producing
temperatures of 79 K (air) or 77 K (N 2 ) in the region of the detecting element
of the cell. Liquid is produced in the collecting chamber by the Joule-Thomson
effect from high pressure air (60 to 400 atm) applied to the inlet.
The weight of the liquefier is 5.3 g and has a heat exchanger of 7.2 mm
diam x 25 mm long mounted on a hollow mandrel, or former, of Monel the
lower end of which is closed by a brass cap to which the nozzle assembly is
fitted. The cavity thus produced is filled with expanded perlite to reduce heat
leaks by convection. A sintered bronze filter of 2 micron pore size is placed
in the inlet cavity of the former to exclude harmful dust particles.
The original prototype liquefier built at the Royal Radar Establishment
had no nozzle assembly or inlet filter. Instead the end of the finned tube
heat exchanger was squeezed together to a fish-tail shape which restricted
the flow to the desired volume. This operated satisfactorily but required
very pure air for its operation as the fine slit became easily blocked.
Contaminant gases such as C0 2 at low vapor pressure sublimed at the cold
end to form minute particles which blocked the exist orifice. The addition
of an inlet filter and the porous nozzle effected a substantial improvement
in reliability, but the requirement for a very high degree of gas purity
remained. Indeed the problem remains today as perhaps the single greatest
impediment to the wider use of miniature Joule-Thomson coolers.
Joule-Thomson Cooling Systems 289
INLET
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0 FINNED
DEWAR 0
FLASK 0 >£AT
0 HAl«
0 COILS
0
0
0
0
0
,0
0
0
0
0
HOLLCW
Fig. 6.15. Cross section of miniature Joule- SCREW
>-.OULE
Thompson air or nitrogen liquefier (after
Mcinroy, 1967).
.
'9
/
. 3
2 •
\
2 7
....
u
/
~
w
I .....
....
(.)
v
~
·e ~·
"...z
v ~ !SO
0
\
/
:::;
~... g. 1-- L s•
z
... v
'
I
.
2
m
I
(.)
1'-
..
....0
" 0 /
!w
~
10 20 30 40 (,!) 0 10 20 30 40 10 20 30 40
Fig. 6.16. Performance of miniature Joule-Thomson nitrogen liquefier (after Mcinroy, 1967).
Joule-Thomson Cooling Systems 291
2 o r--c=F=:t::;::r=-r-, "20
o;;
160
!.,.:;:~ /
AIR
•
I
v
!
• ot CONSTANT 20MPo
1
~ I 51---t--+--t-f---;t'l ~1!5
J 120
I
0.10
HYO~OGEf>~V .!
a:
"'~10 I .
u
~
80
TIME TO H2
t IQUEFY
FROM TIME OF
PRODUCTION OF
LIQUID AIR
:~
",g Q. II \
31:
~ 0 51---t--+--t-+----l
"z
~0!5 0
1\_
0
..."' u
" "'0"'
0 04 08 12 ~ 0 04 0.8 12 0 04 08 12
Fig. 6.17. Performance of miniature Joule-Thomson hydrogen liquefier (after Mcinroy, 1967).
3 I"'-
"'-...
z
0
2 ~
~ ~
"'
:::>
0 ~
I
..... Fig. 6.18. Operating time for miniature
Joule-Thomson self-regulating cooler
(Hymatic Model No. MAC 222-200) as
260 280 300 320 340 360 380 a function of ambient temperature (after
OPERATING I STORAGE TEMPERATURE Stephens, 1970).
::
COOLER/CELL ASSEMBLY - . .
I
g-~-o '
I
AIR CJ SOLENOIDVAL\.~• CLEANER
I
AVERAGE FREE
DURING
CHANGING:
FREE AIR
THROUGHPUT
BETWEEN CJ GASCLEA~=WI.1. I
AIR DELIVERY
99 dm'/min
REACTIVATIONS
283m' ~ SELF SEALING COUPLING AIRCRAFT SHOE I MISSILE .I
(3.5ft'/min). (10,000ft') ~ ' PYLON
i.e., 7 BOTTLES i.e., 400
PEA HOUR BOTTLES
Fig. 6.19. Schematic arrangement for complete refrigerating system to provide miniature
Joule-Thomson cooler system for heat-seeking air-to-air missiles aboard aircraft (after
Stephens, 1970).
Joule-Thomson Cooling Systems 293
In all of the above we have explored various ways and means to utilize
Joule-Thomson expansion alone to produce low-temperature refrigeration
and gas liquefaction. Many other systems exist which incorporate both
cooling engines and a Joule-Thomson expansion system. One such combi-
nation mentioned in discussion of the Reitdijk expansion-ejector system
was described by Prast (1968). This combines a double-expansion Stirling
engine as the two-stage precooler of a helium liquefaction system. Haarhuis
(1967) described a similar system incorporating two Philips double-
expansion Stirling engines.
Chapter 5 includes the description of a helium liquefier built by A.
D. Little Inc. combining a three-stage Gifford-McMahon cooling engine
with a Joule-Thomson expansion process. Further development of that
concept was described by Stuart and Hogan (1965).
Chapter 7 describes many Claude cycle cryocoolers that include Joule-
Thomson expansion valves particularly where the fluid is to be liquefied.
It is customary for the fluid to be expanded partially in the engine and
partially in the Joule-Thomson valve. Expansion in the engine occurs
isentropically and achieves a greater temperature difference than isenthalpic
(JT) expansion over a given pressure ratio. However, expansion of two-
phase fluids in engines presents problems, so that the preferred combination
is partial expansion in an engine to approach the saturated vapor condition
and then Joule-Thomson expansion to achieve liquefaction.
REFERENCES
Baehr, H. D. (1963). "On the Thermodynamics of the Cold-Air Cycle with Throttling." Proc.
XI Int. Congr. of Ref./Prog. in Ref. Sci. Techno/. 1, 319-328, Munich, Germany (Pubs.
Pergamon Press, Oxford and Vertag Muller, Karlsruhe).
Barron, R. (1966). Cryogenic Systems. McGraw-Hill Book Co., Toronto.
Buller, J. S. (1971). "A Miniature Self-Regulating Rapid-Cooling Joule-Thompson Cryostat."
Adv. Cryog. Eng. 16, 205-213.
Carbonnell, E., Chovet, P., Johannes, C., Marinet, D., and Solente, P. (1972). "Refrigerator
without Moving Parts at Low Temperature Able to Cool Down to 90-100 K." Proc. 4th
Int. Cryo. Eng. Conf., Eindhoven, pp. 68-70.
Currie, R. B. (1967). "A Joule-Thomson Laboratory Expander." Adv. Cryo. Eng.12, 557-563.
Daunt, J. G. (1956). The Production of Low Temperatures down to Hydrogen Temperatures.
Handbuch der Physik Bd IV, Springer Verlag, Berlin.
Geist, J. M., and Lashmet, P. K. (1960). "Miniature Joule-Thomson Refrigeration Systems."
Adv. Cryog. Eng. 5.
Geist, J. M., and Lashmet, P. K. (1961). "Compact Joule-Thomson Refrigeration Systems
15-60 K." Adv. Cryog. Eng. 6.
Joule-Thomson Cooling Systems 295
Haarhuis, G. J. (1967). New Type Helium Liquefier." Proc. XII Int. Cong. of Ref./Prog., in
Ref. Sci. Techno/. 1, 121, Madrid (llR Paris).
Haisma, J., and Roozendaal, K. (1967). "Investigation of the Behaviour of an Expansion-
Ejector in the Low Temperature Region beyond the A- Transition of Helium." Proc. XII
Int. Cong. of Ref./Prog., in Ref. Sci. Techno/. 1, 111-210, Madrid (llR Paris).
Hampson, W. (1895). British Patent 10165.
Joule, J.P. (1852). Sci. Pap. 2, p. 216.
Joule, J.P., and Thomson, W. (1852). Phil. Mag. 4, 481.
Keesom, W. H. (1933). Commun. Phys. Lab. Univ. Leiden, suppl. 76a.
Kessler, G. (1970). "Joule-Thomson Hydrogen Refrigerator Target." Adv. Cryog. Eng. 15,
443-446.
Linde, H. V. (1897). German Patent 88824 [see also Z. ges Kiilteind 4, 23 (1897); and the
Engineer, Nov. 13 and 20 (1896)].
Mann, D. B., Bjorkland, W. R., Macinko, J., and Hiza, M. J. (1960). "Design, Construction
and Performance of a Laboratory Size Helium Liquefier." Adv. Cryog. Eng. 5, 346-353.
Mcinroy, J. (1967). Miniature Joule-Thomson Coolers." Proc. XII Int. Cong. of Ref./Prog.
in Ref. Sci. Techno/. 1, 59-68, Madrid (llR, Paris).
Nicholds, K. E. (1968). "Low Temperature Devices for Laboratory Operation." Proc. 2nd
Int. Cryo. Eng. Conf., pp. 65-66, Brighton Iliffe, Sci. and Tech. Pub!., Guildford, U.K.
Nicholds, K. E. (1970). (a) "Performance of Self-Regulating Joule-Thomson Minicoolers,"
pp. 277-282. (b) "Miniature Cryogenic Cooling Systems for an Upper Atmosphere
Infrared Research Programme," pp. 283-286, Proc. 3rd Int. Cryog. Eng. Conf., Berlin,
Iliffe Sci. and Tech. Pub!., Guildford, U.K.
Parkinson, D. H. (1959). "Some Problems in the Design of Helium Liquefiers Based on the
Joule-Thomson Effect." Proc. X Int. Cong. of Ref./Prog. in Ref. Sci. Techno/. 1, 53-57
(Pub!. Pergamon Press, Oxford).
Parkinson, D. H. (1967). "Miniature Refrigeration Systems-A Review." Proc. XII Int. Cong.
of Ref./Prog. in Ref. Sci. and Tech., Vol. 1, 69-77, Madrid (11R, Paris).
Prast, G. (1968). "A 3.5 K Refrigerator Based on the Three-Space Stirling Refrigerator."
Proc. 2nd Int. Cryo. Eng. Conf., Brighton, Iliffe Sci. and Tech. Pubs. Ltd., Guildford, U.K.
Rietdijk, J. A. (1966). "The Expansion-Ejector-A New Device for Liquefaction and
Refrigeration at 4 K and Lower." Included in Liquid Helium Technology, Bull. IIR,
Annexe 1966-65, p. 241.
Stephens, S. (1970). "A Self-Regulating Miniature Joule-Thomson Refrigerator." Appl.
Cryog. Techno/. 3.
Stuart, R. W., and Hogan, W. H. (1965). "A Small Helium Liquefier." Adv. Cryog. Eng. 10,
62-68.
Chapter 7
INTRODUCTION
T T
s s
(a) ISOTHERMAL {bl ISENTHALPIC
/
8
s s
(c) ISENTROPIC (d I POLYTROPIC
can be a liquid (gas previously liquefied) or wet vapor (gas cooled below
the critical temperature with liquid droplets entrained in it). The expansion
can approximate various idealized thermodynamic processes, namely,
isothermal, isenthalpic, isentropic and polytropic. These are represented
on the temperature-entropy diagram (T-S plane) as shown in Fig. 7.1.
Assume first the gas is at some high pressure and low temperature
TA. The gas expands to a low pressure and heat is supplied to bring the
temperature back to the initial temperature TA. An apparatus to carry out
such comparative work is shown in Fig. 7.2. It consists of a piston reciprocat-
CYLINDER HEAD
I VE
Fig. 7.1c. The cylinder head then becomes perfectly conducting and heat
is transferred from the cylinder surroundings to heat the gas to the final
temperature, TB. The shaded area E-B-C-D in Fig. 7.1c is larger than
the isenthalpic process and smaller than the isothermal process.
The polytropic process of Fig. 7.1d follows a similar pattern to the
isentropic with the difference that the temperature at the end of the
expansion Tk, is greater than at the end of isentropic expansion TE. This
difference may arise from two causes, limited heat transfer during the
expansion (usually the case with reciprocating engines but not turbines),
and irreversibilities (principally friction) which occur during expansion
(particularly in turbines) to increase entropy and, hence, the fluid tem-
peratures. The effect of limited heat transfer and irreversibilities during
expansion is to decrease the potential refrigerating effect compared with
the isentropic process. The potential refrigerating effect is the shaded area
E'-B-C' -D in Fig. 7 .1d, smaller than the area for the isentropic process
by area E-E'-C'-C.
Figure 7.1 shows clearly that in terms of refrigerating effect the order
of preference of the four processes is isothermal, isentropic, polytropic,
and isenthalpic.
In practice isothermal expansion cannot be achieved. It would require
infinite rates of heat transfer or the expansion to proceed infinitely slowly.
It is customary to adopt an isentropic process as the reference for expansion
engines. This cannot be achieved in practice and realistic cycle analyses
adopt an arbitrary value (60% to 85%) for the isentropic efficiency to
account for the polytropic process.
The same work is necessary to compress the gas regardless of the
subsequent expansion. Since the Joule-Thomson process produces so little
refrigerating effect compared with isentropic expansion it is to be expected
that Joule-Thomson cryocoolers will be less efficient than those incorporat-
ing expansion engines. The principal reason for using Joule-Thomson
expansion is simplicity. It is much easier and cheaper to provide a simple
throttle valve than an expansion engine with means to absorb the work
produced. All small domestic refrigeration systems operating on the vapor-
compression cycle have JT expansion valves.
In larger refrigerating and cryogenic cooling systems the additional
expense of expansion engines can be justified by the improved efficiencies.
Even there Joule-Thomson expansion systems are invariably used for the
final stage of expansion when there is partial liquefaction of the working
fluid. Two-phase working fluids are difficult to handle in expansion engines.
Daunt (1956) discussed isenthalpic and isentropic expansion in
more formal thermodynamic terms than the intuitive-descriptive method
followed above. The following is abstracted from Daunt's presentation:
Claude and Joule-Brayton Systems 301
as=
aT)
(ap T ( av)
s = C" aT "
(7.2)
(7.3)
whence
(7.4)
v
(o:H/as)=1-T
(aT)
a;; P (7.5)
Claude Cycle
A Claude cycle system for gas liquefaction, showninFig. 7.5, includes
(a) compressor equipped with an aftercooler, air or water cooled, (b)
recuperative heat exchangers, designated as high- (ambient), intermediate-,
and low-temperature exchangers, or alternatively upper-, intermediate-,
Claude and Joule-Brayton Systems 303
MAKE-UP
CHARGE
COMPRESSOR
------------
ENTROPY
Strictly speaking, the two separate streams: (a) process (3-9'), and (b)
process (3-4-5-6-8-9') cannot be represented on the T-s plane because
the diagram is constructed for unit mass of the working fluid. It is useful
to represent the processes as followed here, providing it is always recalled
that per unit mass of working fluid compressed, fraction (1- x) passes
Claude and Joule-Brayton Systems 305
through the expansion engine and only x is involved in the actual liquefac-
tion process.
y =x(H6-H7) (7.12)
Hs-H?
(9) The fraction of gas not liquefied but existing at state (8) is
(x -y) (7.13)
The work done by the expansion engine may be used to partially offset
the work of compression. In that case the net cycle work is
(We- WE)= (H;- Hl)-(1-x)(H3-H9) (7.14)
306 Chapter7
The work generated by the expansion engine is relatively small (about 10%
of the compression work). Accordingly the expansion engine and com-
pressor are generally not coupled. The expansion work can be absorbed
in useful ways to drive pumps or generators or simply dissipated in hydraulic
or electrical dynamometers.
30
8
~0 PRESSURE AT END OF
:::; COMPRESSION (MPa)
~
2.0
"'~
"'
J::
'
~
1-
,.[:5
"'"'50
1.0
"'"'
"'
0
a.
Fig. 7.6. Power requirement per unit mass
FLUID AIR
of gas liquefied as a function of the frac-
0 02 04 0.6 08 IO tion of gas expanded in the Joule-
FRACTION OF GAS EXPANDING IN Thomson valve (for air) (after Daunt,
JOULE-THOMSON VALVE 1956).
Claude and Joule-Brayton Systems 307
efficient operation, and (b) as the pressure increases the maximum efficiency
and the fraction to be expanded in the Joule-Thomson expansion increase.
The conditions of maximum efficiency for operation of an air liquefier
were provided by Davies (1949) in the data reproduced in Table 7.1. This
shows, for several different compressor pressures, the fraction of gas passing
through the expansion engine, i.e., (1- x ), the optimum temperature for
gas to enter the expansion engine, the theoretical power required in kW hr
per liter of liquid produced and, finally, the actual power required (kW hr
per liter of liquid produced). These data confirm, as in Fig. 7 .6, that the
optimum fraction of gas expanded in the engine decreases and the system
efficiency improves as the pressure increases. The temperature for
admission of the gas to the expansion engine increases as the pressure
increases. At a maximum cycle pressure of 200 MPa the temperature has
increased to 293 K, approximately the ambient temperature.
Heylandt Cycle
Heylandt was the first to appreciate that the optimum temperature
for admission to the expansion engine increased as the cycle pressure and
system efficiency increased. Heylandt introduced a modified Claude cycle
by eliminating the upper-stage recuperative heat exchanger of Fig. 7.5 and
elevated pressure level (for an air liquefier) from 4 MPa (used by Claude)
to 20 MPa. He further increased the fraction of mass flow passing to the
Joule-Thomson valve from 0.2 (at 4 MPa) to 0.4 (at 20 MPa).
As evident from Table 7.1 there is little difference in the theoretical
power requirement for the Heylandt-Claude cycle but there are substantial
practical advantages. The actual power required per unit mass liquefied
with the Heylandt system is only 75% of the power with the original Claude
system. Air is admitted to the expansion engine at ambient temperature
immediately following the aftercooling process on leaving the compressor.
About 60% of the air passes through the expansion engine and exits from
the engine at a temperature and pressure of about 150 K and 0.1 MPa
(condition 9' in Fig. 7.5). The relatively high temperature of operation of
the expansion engine has important practical advantages. The high-
temperature recuperator is eliminated and appreciable savings gained
thereby. The lubrication problems of the piston in the cylinder and the
thermal insulation or isolation problems are greatly minimized. The
efficiency of the Heylandt cycle is the highest of all simple Claude systems
and was much favored for air liquefaction systems particularly in the United
States.
Gardner and Smith (1960) have given a comprehensive summary of
the thermodynamics of the Heylandt cycle from the aspect of power
~
Table 7 .1. Conditions for Operation of Claude Cycle Air Liquefier at Maximum Efficiencya
Turbo
expansion Original
machine Claude Heylandt
(Kapitza) system cycle
n
=-
?.
~
....
Claude and Joule-Brayton Systems 309
A schematic of the flow system for the Collins air liquefier is shown
in Fig. 7.7. The system operates as follows:
(i) Inlet air (at atmospheric conditions) enters the compressor at condi-
tion (1).
(ii) The air is compressed to 13 MPa and passes to a watercooled
aftercooler.
(iii) The compressed fluid then passes through the reversing valves,
the high- and intermediate-temperature recuperators and the low tem-
perature reversing valves. At condition (6), downstream of the low-
temperature reversing valve, the pressure and temperature is 13 MPa and
125 K. By the time this temperature is achieved all the significant water
vapor, carbon dioxide, and other miscellaneous impurities in the inlet air
have precipitated on the walls of the recuperative heat exchangers. If
operated continuously the flow passages in the recuperators would soon
be blocked by accumulations of these precipitated solids. However, the
reversing valves operate at approximately 3-min intervals to cause the
return flow to pass through the passages containing the precipitated con-
taminants. As the return flow warms in passing through the recuperators
Claude and Joule-Brayton Systems 311
the contaminants reevaporate from the walls of the flow passages, are
entrained in the return flow and exhaust from the system at approximately
atmospheric conditions. In this system there is no recirculation of the return
effluent as in the other cases considered earlier.
(iv) At condition (6), pressure 13 MPa, temperature 125 K, the flow
separates into two streams. The minor stream enters the condenser/re-
cuperator where it condenses under pressure to the liquid phase, condition
(7), and hence passes to the liquid reservoir and subsequently to a
rectification column for separation into the liquid-nitrogen useful product
and gaseous oxygen for release to the atmosphere.
(v) The major stream at condition (6) passes to the expansion engine
where it expands to a final pressure and temperature of approximately
0.15 MPa and 87 K.
(vi) The low-pressure, low-temperature effluent enters the con-
denser/recuperator to accomplish the condensation of the liquid air in (d)
above and hence passes to exhaust through the upper stage recuperators.
In passing the effluent reevaporates and absorbs the contaminants precipi-
tated in the course of the initial cooling.
(vii) Prior to fluid actually entering the expansion engine the oppor-
tunity is taken to provide auxiliary cooling to the intermediate stage
recuperator. In passage to the expansion engine the fluid at condition (6),
pressure 1.3 MPa, temperature 125 K, is further subdivided. A small
amount of the gas is fed through a third set of flow passages in the
intermediate stage recuperator and reheated to about 200 K, condition
(12). The major stream is partially throttled to about 1.25 MPa, condition
(13), to provide sufficient pressure potential to induce flow through the
intermediate recuperator.
(viii) The two streams at condition (12) and (13) are mixed to provide
inlet pressure and temperature, condition (14), of 1.25 MPa and 146 K to
the expansion engine.
Auxiliary cooling in the intermediate recuperator is necessary to cool
the pressurized incoming stream to 125 K (at condition (5)), where all the
carbon dioxide will have precipitated out.
Collins (1956) conceived the use of reversing recuperative exchangers
as an alternative to Frankl regenerative exchangers for relatively low-
pressure systems. Application in systems with high-pressure streams would
pose formidable sealing problems in the low-temperature reversing valves.
On flow reversal the high-pressure fluid contained in the void volume of
the heat exchangers simply passes to exhaust and serves no useful purpose.
This so-called "carry-over loss" would incur a serious decline in the system
efficiency at high pressures.
312 Chapter7
MAKE UP
CHARGE
LIQUID
PRE COOLER
CONSTANT
PRESSURE . . __
ENTROPY S
Fig. 7.10. Schematic diagram of the dual-pressure Claude cycle cooling system.
Table 7.2. Comparison of Liquefaction Systems for Air as the Working Fluida
Work per
Liquid unit mass
yield, liquefied Figure of
y kJ/kg merit
"T1 = 294 K, p 1 = 0.1 MPa, 11< is the compressor overall efficiency, 11. is the expander overall efficiency, e
is the heat exchanger effectiveness (after Barron, 1966).
MULTIPLE-EXPANSION ENGINES
MAKE-UP
AFTER COOLER
UPPER-STAGE
RECUPERATOR
INTERMEDIATE
RECUPERATOR
RECUPERATOR/
CONDENSER
RESERVOIR
al PARALLEL b l SERIES
ARRANGEMENT ARRANGEMENT
unity and the pressure 3, in arbitrary units. The isotherm is drawn for a
perfect gas according to the equation PV = canst. The isotherm for T2
passes through the point where P = 3 and V = 0.5 units. This is equivalent
to specifying the temperature as T 2 is half the temperature at T 1 • Assume
an equal mass of gas expands in both the upper and lower expansion stages.
Then the volume of a given mass of gas at temperature T 2 will be half the
volume of the same mass of gas at temperature T 1 •
Diagrams (a) and (b) have been drawn for the above situation. In both
cases the expansion process was assumed to be an arbitrary polytropic
process so that the temperature is somewhat less at the end of expansion
than at the beginning. To a first approximation the work produced in the
upper stage engine is twice the work produced in the lower stage engine
(comparing the areas of the two diagrams).
Now consider the series arrangement represented in the diagrams (c)
and (d) of Fig. 7.14. In this case the sum of the two separate masses
experiences the expansion but in two stages. The isotherms for T 1 and T 2
have been drawn for the situation where double the mass of gas was present
for the diagrams (c) and (d) than was assumed for diagrams (a) and (b).
Claude and Joule-Brayton Systems 321
"'
SERIES ARRANGEMENT
0:
iil
<J)
"'
0:
Q_
2
(b) (d)
VOLUME
Fig. 7 .14. Work diagrams for the upper and lower stage expansion engines in parallel (diagrams
a and b) and series arrangements (diagrams c and d).
This simply means that the isotherms pass through points where, when
P = 3, the volume is 2 for T 1 and unity for T2 •
The arbitrary decision was made to select the expansion ratio to be
equal in both stages, i.e., Pr(upper) = PrOower) = .J3 = 1.73. The work diagrams
then correspond to those shown in diagrams (c) and (d) of Fig. 7.14. As
before an arbitrary polytropic expansion was assumed where the tem-
perature at the end of the expansion was somewhat less than the tem-
perature at the start.
The sums of the areas of diagrams (a) and (b) are greater by approxi-
mately 10% than the areas of the diagrams for (c) and (d). This indicates
that for this arbitrary idealized example the parallel arrangement will
produce more refrigeration (evidenced by the greater work output) than
the series arrangement. Confirmation of this view was provided by
Carbonnell (1974) in a discussion of the relative merits of parallel and
series arrangements. He confirmed the advantage of the parallel arrange-
ment in terms of a 10% increased refrigeration capacity for given constant
adiabatic efficiency and pressure expansion ratio.
Carbonn ell was concerned with large-scale refrigeration using turbines
and pointed out that the series arrangement was characterized by high flow
rates over low expansion ratios. These factors contribute to increased
turbine adiabatic efficiency great enough to favor series coupling for tur-
bines. Reciprocating engines, on the other hand, can handle high expansion
322 Chapter7
ratios at relatively low mass flow rates. This favors the parallel arrangement
where reciprocating expanders are used. Confirmation was provided by
the eminent Russian cryogenic engineer I. B. Danilov in the discussion to
a paper concerned with multiple-expansion helium liquefiers (Kapitza and
Danilov, 1968).
Thereafter was laid the foundation of the great French industrial syndicate
L' Air Liquide, Ltee.
Claude was much concerned about piston lubrication with metal rings
but eventually (1912) adopted the leather cup as a piston seal. It is surprising
this took so long, for leather had long been used on air expansion engines
for ordinary refrigeration and Claude must have been aware of this.
In 1920, Claude constructed hydrogen expansion engines to obtain
temperatures below the boiling point of liquid nitrogen. These engines
were used to separate hydrogen from coke oven gas. Collins (1958)
described hydrogen expansion engines "currently in use" constructed in
324 Chapter7
were given with estimates of the various losses due to heat leaks and ring
friction. Comparative data were presented for the predicted and actual
engine performance on test.
*The piston would not be naturally centered in the cylinder as one might expect. Instead
the high-pressure gas flowing past the piston from the high-pressure to the low-pressure
region causes the piston to move to one side of the cylinder and provide the flow path with
least restraint.
Claude and Joule-Bray ton Systems 327
I
Z4 ·
"*15I
I
i
I 13
I
I
i :::!8
37
38 /16
34
20cm
10 .._.._.,..-.ll
0,..H...._..._..,-ll....._,
Scale
Fig. 7.17. Sectional view of the Kapitza helium liquefier including an expansion
engine with hydrodyna mic lubrication (after Kapitza, 1934).
328 Chapter7
Description of Liquefier
A flow diagram of the actual cycle is given in Figure [7.19]. It differs
from the idealized cycle, not only because the heat exchangers and expansion
engines are necessarily imperfect, but also because practical considerations
have influenced the choice of apparatus and procedures. For the sake of
compactness of the liquefier and greater production of liquid from available
compressed helium, liquid nitrogen is employed to the extent of its utility.
The transfer of heat from gaseous helium to liquid nitrogen evaporating at
a constant temperature is irreversible, of course, and a net increase of entropy
is incurred. A final difference lies in the substitution of a throttle valve for
the fourth engine of the cycle.
In the flow diagram shown [Figure 7.19], compressed helium (about
12 atm) from the compressor is treated for entrained oil in an oil trap and
for vaporized oil in a refrigerated heat exchanger 3. Thereafter the stream
divides, about 8 percent going to heat exchanger 4, in which it is cooled to
80 K by means of liquid nitrogen and then expanded in engine E 1. the
remainder going to the principal heat exchanger 5. The temperature of the
gas in heat exchanger 5 ranges from room temperature at the upper end to
15 K at the bottom. At the zone of exchanger 5 where the temperature is
40-45 K, a second fraction (about 15 percent of the whole) of the compressed
helium is led off for expansion in the second engine, E 2 • At the lower end
w
CD
._::>
a::
w
u._
(f)
z
<l
._
a::
0
..J
w
I
(f)
z
._0
<l
0
<l
a::
4-STAGE
COMPRESSOR
Fig. 7 .19. Flow diagram of relatively large-scale Collins helium liquefier (after Collins, 1952).
Claude and Joule-Brayton Systems 331
Fig. 7.20. Collins-Claude cycle helium liquefier (CTI Model No. 1430).
( ·-'1,/
I. .I
b
-----'f---.>--./ HELIUM
IMPURE HELIUM GAS GAS BAG
FROM COLLECTING STATIONS (OPTIONAL)
were constructed by Cryovac Inc. of Columbus, Ohio for the U.S. Air
Force Arnold Engineering Development Center.
Zeitz (1963) described a closed cycle helium refrigerator for 2.5 K
constructed by Air Products and Chemicals Inc., Allentown, Pennsylvania.
This unit was a compound device consisting of a high-pressure (0.2-2 MPa)
double expander Joule-Brayton refrigerator generating refrigeration at
10 K to precool a low-pressure (0.009-0.5 MPa) Linde-Hampson system
generating refrigeration at 2.5 K. The expansion engines of the high-
pressure system are acknowledged to be "similar to the machines developed
by Collins and Kapitza."
Furihara et al. (1968) described what they whimsically called "a new
type of closed cycle helium refrigerator for laser applications." From the
description given it appeared to be a carbon copy of the Collins-Claude
cycle double expander with final-stage Joule-Thomson expansion. A
Collins-type expansion engine was also incorporated for helium expansion
in a compound helium-hydrogen-neon refrigerating system described by
Zeitz et al. (1965). It was constructed by Air Products and Chemicals Inc.
for the NASA Lewis Research Center to cool a cryogenic magnet used for
plasma heating and acceleration experiments.
Another large-scale (1.4 kW at 4.5 K) helium refrigerator incorporat-
ing two-stage Collins expansion engines was described by Winters and
Snow (1966). The unit was constructed by Air Products and Chemicals
Inc., Allentown, Pennsylvania for the NASA Langley Research Center.
The same company also constructed the lower-capacity Claude cycle system
using a single Collins expansion engine with no precooler described by
Meier and Currie (1968). These units were built to provide 1 W of useful
refrigeration at 4 K with additional refrigeration at 30 K. The units were
designed for long-term reliable automatic operation with redundant com-
pressor and expansion engines provided with automatic switches following
component failure. At the time of reporting (1968) 10 systems in field use
had accumulated 32,000 operating hours with only one machinery failure.
the piston surface. Work is transferred via the piston rod to the crank-gear
and hence to a dynamometer.
Gas flow and the resultant work diagram for the engine are illustrated
in Fig. 7.23 (after Kneuer et al., 1968). At condition 1, with the piston in
the bottom dead center position, the ring groove on the piston is aligned
with the inlet annular channel in the cylinder. High-pressure gas enters
the engine filling the clearance space, volume V1. principally the central
bore of the piston. As the piston rises the piston ring groove and the
cylinder inlet annular channel move out of alignment thereby closing off
admission of further fluid at high inlet pressure. The fluid contained in the
cylinder expands to a lower pressure and work is produced.
At condition (2) the lower edge of the piston clears the annular exhaust
channel in the cylinder releasing expanded fluid to exhaust. The pressure
drops rapidly to the exhaust pressure [condition (3)] as the piston ascends
to top dead center. Finally the piston descends closing off the exhaust
annular channel and compressing the gas retained in the cylinder to condi-
tion (4). Near the piston top dead center position the inlet port is again
uncovered and fresh charge is admitted to the cylinder to repeat the cycle.
There is a leakage of high-pressure gas to exhaust in the annulus
formed by the piston in the cylinder. The leakage is indicated, by the arrows
on Fig. 7.23, as a flow from (a) the inlet to the exhaust annular channel,
(b) the underside of the piston to the exhaust channel, and (c) the inlet
channel to the space above the piston. This leakage flow is the hydrodynamic
gas lubricant for the piston. Kneuer et al. (1968) identify the flow as a loss
that is unavoidable with valveless engines and give the following equation
for the magnitude of the flow loss:
v = (pt-p~)a 3 U/24p1L{3 (7.15)
where v is the volumetric rate of flow, P1 is the inlet pressure, P2 is the
336 Chapter?
,.,
1(P 1 )
... by means of a large number of grooves on the piston surface. The clearance
between successive grooves has a proper profile and achieves high stabilizing
forces.
et al. wrote that over 50 valveless engines had been incorporated in helium
liquefiers in the previous four years. Baldus and Sellmaier (1968) described
a complicated system to generate temperatures of 1.8 K in a helium II
refrigerator utilizing a Doll-Eder valveless expansion engine in a Claude
cycle system with a low-pressure Joule-Thomson expansion circuit generat-
ing 30 W of refrigeration at 1.8 K.
Publication of the original paper in 1964 stimulated considerable
research on valveless expanders. Leger et al. (1971) described, in fascinating
detail, the design and development of a Claude cycle helium liquefier with
two valveless expansion engines, in a parallel arrangement. One interesting
feature of the Leger engine was the use of different materials for the
cylinder and the piston. The piston was made of material having a lower
thermal expansion coefficient than the cylinder. On starting, the clearance
of the piston in the cylinder was as high as 11 f.1. and reduced, due to
differential contraction, to 2.5 f.1. at the normal running stage. The use of
different materials is contrary to conventional practice where the same
materials are used for piston and cylinder to maintain a constant clearance
during cooldown. In early engines with the piston and cylinder made of
the same steel, Leger found the piston cooled more slowly than the cylinder
and resulted in premature wear contributing to reduced performance at
low temperatures. In some cases the differential rate of cooling caused the
engine to seize and resulted in deformation of the cylinder.
Matsubara et al. (1968) have described a hybrid expander with the
inlet high pressure fluid entering through valveless ports as in the Doll-Eder
engine but exhausting through a crank-operated valve of conventional
design. This allows independent selection of exhaust valve timing but has
no other advantage over the Doll-Eder engine and introduces the com-
plexities of mechanically operated valves.
Another valveless, gas-lubricated expansion engine of exceptional
interest was described, albeit briefly, by Clarke et al. (1968). The engine
was a free-piston machine taking the form of a double-ended Doll-Eder
engine with a single reciprocating element, the piston, operating in both
cylinders. One end operated as a conventional Doll-Eder expansion engine
at low temperature. The other end operated as an auxiliary compressor
absorbing the work of the expansion engine to compress helium. The engine
had a bore and stroke of 31.75 mm (1~in.) and operated at a frequency
of 45Hz with an inlet pressure and temperature of 2.5 MPa and 24 K. It
was said that
Piston
Distance piece
30mm 10
31 mm OD
300 mm long
Sliding ring
Stroke
42 mm
Fig. 7 .24. Sketch of expansion engine with valves
operated by piston motion (after Claudet et a/.,
1972).
c)
Fig. 7.25. Piston and valve positions in Claudet-Verdier engine operation (after Claudet et
al., 1972).
spring force F holding the inlet valve open decreases and the inlet valve
closes.
(4) Expansion. The piston descends with both valves closed. The
pressure decreases as the expansion proceeds.
The performance of the prototype Claudet/Verdier engine was most
encouraging. With inlet-outlet pressures of 0.9 and 0.14 MPa they observed
isentropic efficiencies in the range 50% to 70%. The temperature difference
across the engine ranged from 5 K with an inlet temperature of 20 K to
26 K with an inlet temperature of 76 K.
Russian work on valveless piston expansion engines was described by
Dobrov and Orlov (1974). Their engine utilized a six-phase gas distribution
cycle in an attempt to overcome some of the deficiencies of the conventional
four-phase Doll-Eder eng1ne. The Dobrov-Orlov concept required mutual
but out-of-phase reciprocation of both the piston and cylinder, but it is
not clear how this was accomplished. One's feeling is that the added
complexity obviates the principal attraction of simplicity in the Doll-Eder
valveless engine.
Parulekar and Narayankhedkar (1972) and (1976) discussed the design
and performance of a valveless expansion engine with novel features
including a regenerative heat exchanger contained within the piston. The
piston was activated by alternating fluid pressure in the "warm" and cold,
expansion, part of the cylinder. Fluid pressure in the fluid drive region was
regulated by solenoid valves activated by the piston motion. There seems
every possibility this system could be made to operate at resonant frequency
as a free piston machine of the same family as that described by Clarke et
al. (1968). It could also be classified as an Ericsson engine of the types
described in Chapter 5.
340 Chapter7
~ HELIUM
~
~ BRASS
~
% STAINLESS STEEL
II P.T.F.E .
Fig. 7.26. Sectional view of formed bellows engine (after Long and Simon, 1954).
DISTANCE TUBE
I l ET
- - -- V ALV E TUBE
Fig. 7.27. Sectional view of welded bellows engine (after Smith, 1967).
~--CRANK CASE
LIOUIO NITROGEN
PRECOOLING
TIMING BELT
HEAT EXCHANGER
0 I STANCE TUBE
ELECTR IC
LIOU IO NITROGEN
MOTOR
RADIATION SHIELD
SURGE CHAMBERS
TIMING BELT
f
lr
Q,
=
"'
H[AO H(AT EXCHANG(R LINEAR ROlARY MOTOR ROTOR
....
STATOR ~
;r
EXPANSION CHAMBER
COMPRESSOR H[ AO I!=
EXPANO[R PISTON
i=
i"'
..
a
"'
-- -- --
GAS eEARING
Fig. 7.29. Diagrammatic view of rotary-stroking engine (after Schulte et al., 1965).
~
"""
344 Chapter7
Smith used bellows having an outer diameter of 3.17 em (1.25 in.) and
an inner diameter of 1.27 em (0.5 in.). Nineteen convolutions of 0.1 mm
(0.004 in.) thickness provided a stroke of 1.9 em (0.75 in.) and a displace-
ment of 1.1 x 10-4 m 3 (0.44 in?). The inlet and discharge pressures were
0. 51 and 0.11 MPa (7 5 and 16 psi). Two bellows to this design failed at
1.3 x 106 cycles (43 hrs) and 1.15 x 10 6 (38 hrs). The design was then
modified as detailed by Smith to a bellows of 35 convolutions and a
diaphragm thickness of 0.2 mm (0.008 in.) with the stress at the inner weld
always in compression. At the time of reporting this assembly had operated
for 7.7 x 106 cycles (258 hrs) without failure at temperatures of 28 and
18 Kat inlet and exhaust, respectively.
sqall resist the temptation. The subject is simply too big for this already
over-long chapter to contain. Instead we consign turbines to a separate
folio, companion to this volume, concerned primarily with large-scale
refrigeration, liquefaction, and gas processing. Such division is in the natural
order of things. Rotary machines, turbines, and compressors, are charac-
teristically best suited for high flow rate, low-pressure ratio systems as
compared with reciprocating devices which favor low-flow, high-pressure
ratio systems.
Much work has been done on miniature turbo expanders and compress-
ors for both Claude and Joule-Brayton cycle systems. It is appropriate to
briefly review that work below since our interest is limited to small- and
intermediate-capacity systems. However, a detailed consideration has been
left to the companion volume. Many papers on miniature turbine systems
may be found in the extensive Bibliography included here under the names
of Colyer and Sixsmith. Some of these are briefly discussed below.
MIXTURE
He&R-12
1034 k Pa
300K
ADSORBER
HELIUM
~ ~g ~Par---+-----f----4----,
~+-----,
648 k Pa
300K
flooded screw compressors and for the new mixed refrigerant cycle. The
Claude cycles included liquid-nitrogen precooled cycles in some cases. The
mixed refrigerant cycle emerged from this comparison indicating appreci-
able (15%) savings in the capital and operating cost.
No studies are known of the Terbot mixed refrigerant cycle in smaller
systems using reciprocating compressors or with auxiliary refrigerants other
than Rl2.
Reciprocating Compressors
Reciprocating compressors are massive, bulky, slow running (up to
50 Hz) machines having a high isentropic efficiency able to handle a virtually
unrestricted range of pressure ratios but comparatively small mass flow.
They are easy to start and to control. They are insensitive to sudden and
large changes in load. They have a high capital cost and, because they
require relatively frequent maintenance, have a high operating cost.
Their efficiency can be improved with multistage compression pro-
cesses having interstage cooling in air or water cooled heat exchangers.
Aftercoolers downstream of the compressor cool the compressed gas to
near ambient temperatures. Two stages of compression are commonly used
for the low-pressure Claude cycles (e.g. pressure ratios of up to 8 or 10).
Above this pressure ratio an increased number of stages may be used. It
is rarely worth going beyond four stages for the degree of improvement
occurs at a progressively diminishing rate.
Reciprocating compressors are found in a variety of shapes, sizes, and
arrangements. They may be double or single-acting crank-connecting rod
machines with automatic spring control or cam regulated valves. They may
be air-cooled or water-cooled, with a variety of drive systems, electric
motors, steam turbines, diesel engines. They exist in many mechanical
arrangements: inline multiple cylinder, vee-cylinder, radial cylinder, or
parallel cylinder with swash-plate drive.
An important aspect of compressors is the method of lubricating the
rubbing surfaces, the reciprocating piston rings, the rotating crankshaft
main bearings and the oscillatory big and little end bearings of the connect-
ing rods. Rubbing surfaces properly lubricated with oil have superior
Claude and Joule-Brayton Systems 349
Rotary Compressors
Rotary compressors may be characterized as small, light-weight, fast-
running machines having a relatively low isentropic efficiency able to handle
a limited range of pressure ratio with very high rates of mass flow. They
are easy to start and control but have a limited range of operation and are
acutely sensitive to sudden large changes in load. They are comparatively
inexpensive and require little maintenance with an excellent record of
reliability. There are three principal types of rotary compressor: axial flow,
centrifugal flow, and rotary screw.
Centrifugal flow compressors are the simplest type with a single rotating
disc carrying vanes. They are low cost and reliabile but have a relatively
low efficiency and are capable of handling only a limited pressure ratio.
Axial flow compressors comprise several stages of fixed and moving
blades in which the gas is successfully compressed. The compression ratio
per stage is very small but there may be many stages and a high overall
compression may be gained. Axial flow machines are expensive but can
have a higher efficiency than centrifugal machines. They are sensitive to
the physical characteristics of the fluid being compressed.
Claude and Joule-Brayton Systems 351
REFERENCES
Aberle, J. L., and Westbook, A. J. (1963). "Liquid Helium and Nitrogen Supply Systems for
Space Simulators." Adv. Cryog. Eng. 8, 190-198.
Anon.-CTI (1978). "Standard Helium Liquefier Model No. 1410." Tech. Mem. 91477,
CTI, Waltham, Massachusetts.
Baldus, W., and Sellmaier, A. (1968). "A Continuous Helium II Refrigerator." Adv. Cryog.
Eng. 18, 434-440.
Barron, R. (1966). Cryogenic Systems. McGraw-Hill Book Co., Toronto.
Becker, H., Doll, R., and Eder, F. X. (1968). "Measurement of the Efficiency of the Valveless
Doli-Eder Expansion Engine." Proc. Second Int. Cryog. Eng. Con[., Iliffe Sci. and Tech.
Pub., Guildford, U.K., pp. 9-11.
Breckenridge, R. W. (1969). "A 3.6 K Reciprocating Refrigerator." Adv. Cryog. Eng. 14,
387-393.
Campbell, D. N. (1967). "Miniature Closed Circuit Single Stage Claude Cycle Refrigerator
for 28 K." Proc. First Cryog. Eng. Con[., Tokyo, pp. 197-198, Heywood-Temple Ind.
Pub., London.
Carbonnell, E. (1974). "Large Scale Refrigerators with Turbines." Proc. Fifth Int. Cryog.
Eng. Con[., Kyoto, pp. 353-356, IPC Sci. and Tech., Guildford, U.K.
Clarke, M. E., and Gardner, J. B. (1968). "New Developments in Expansion Machinery for
Low Temperature Refrigerators." Proc. First Cryog. Eng. Con[., Tokyo, pp. 273-275,
Heywood-Temple Ind. Pub., London.
Class, C. R., Spero, R. P., and Mcintosh, G. E. (1960). "Efficient Utilization of Ortho-Para
Catalyst." Adv. Cryog. Eng. 3, 64-72.
Claude, G. (1913). Liquid Air Oxygen and Nitrogen. P. Blakiston Son and Co., Phila,
Pennsylvania.
Claudet, G., and Verdier, J. (1972). "Simplified Cryogenic Reciprocating Expansion Engine."
Proc. Fourth Int. Cryog. Con[., pp. 80-82, Eindhoven, IPC Sci. and Tech. Press Ltd.,
Guildford, U.K.
Collins, S. C. (1947). "A Helium Cryostat." Rev. Sci. Instrum. 18, 157.
Collins, S.C. (1952). "Helium Liquefiers." Science 116, 289.
Collins, S.C. (1956). Helium Liquefiers and Carriers. Handbuch der Physik Vol. XIV, 113-135.
Collins, S.C. (1966). "Helium Refrigerator and Liquefier." Adv. Cryog. Eng. 11, 11-15.
Collins, S.C., and Cannaday, R. L. (1958). Expansion Machines for Low Temperature Processes.
Oxford University Press, Oxford.
Collins, S.C., and Streeter, M. H. (1968). "Refrigerator for 1.8 K." Proc. First Cryog. Eng.
Con{., Tokyo, pp. 215-217, Heywood-Temple Ind. Pub., London.
352 Chapter7
Cowans, K. W., and Walsh, P. J. (1965). "Continuous Cryogenic Refrigeration for Three to
Five Micron Infrared Systems." Adv. Cryog. Eng. 10, 468-476.
Danilov, I. B., and Kovatchev, V. T. (1972). "Two-Stage Expansion Engine with Differential
Piston." Proc. Fourth Int. Cryog. Eng. Conf., pp. 87-89, IPC Sci. and Tech. Press Ltd.,
Guildford, U.K.
Daunt, J. G. (1956). The Production of Low Temperatures down to Hydrogen Temperatures.
Handbuch der Physik Vol. IV, Springer-Verlag, Berlin.
Davies, M. (1949). The Physical Principles of Gas Liquefaction and Low Temperature
Rectification. Longmans Green and Co. Ltd., London.
Dobrov, V. M., and Orlov, A. V. (1974). "The Valveless Piston Expansion Engine with
Six-Phase Gas Distribution Cycle." Proc. Fifth Cryog. Eng. Con{., Kyoto, IPC Sci. and
Tech. Press Ltd., Guildford, U.K.
Doll, R., and Eder, F. X. (1964). "Gas Lubricated Low Temperature Piston Expansion
Engine without Control Valves." Adv. Cryog. Eng. 9, 561-564.
Furihara, A., Kikuchi, K., and Nagano, H. (1968). "Closed Cycle Helium Refrigerator for
Laser Applications." Proc. First Int. Cryog. Eng. Con{., Tokyo, pp. 213-214, Heywood-
Temple Ind. Pub., London.
Gardner, J. B., and Smith, K. C. (1960). "Power Consumption and Thermodynamic Reversi-
bility in Low Temperature Refrigeration and Separation Processes." Adv. Cryog. Eng.
5, 32-46.
Gessner, R. L., and Colyer, D. B. (1968). "Miniature Claude and Reverse Brayton Cycle
Turbomachinery Refrigerators." Adv. Cryog. Eng. 13, 474-484. See also: Colyer,
D. B., and Gessner, R. L. (1968). "Miniature Cryogenic Refrigerator Turbomachinery."
Adv. Cryog. Eng. 13, 485-493.
Hood, C. B., Vogelhuber, W. W., and Barnes, C. B. (1964). "Helium Refrigerators for
Operation in the 10-30 K Range." Adv. Cryog. Eng. 9, 496-506.
Johnson, R. W., Collins, S. C., and Smith, J. L. (1971). "Hydraulically Operated Two-Phase
Helium Expansion Engine." Adv. Cryog. Eng. 16, 171-177.
Kapitza, P. (1934). "The Liquefaction of Helium by an Adiabatic Method." Proc. R. Soc.
London Ser. A 147, 189.
Kapitza, P. L., and Danilov, I. B. (1968). "Cascade Helium Liquefiers with Piston Type
Engines." Proc. First Int. Cryog. Eng. Con{., Tokyo, pp. 228-231, Heywood-Temple Ind.
Pub., London.
Kneuer, R., and Turnwald, E. (1968). "Measurements and Experience with Valveless Piston
Expansion Engines." Proc. Second Inst. Cryog. Eng. Con{., Brighton, U.K., pp. 12-15,
Iliffe Sci. and Tech. Pub. Ltd., Guildford, U.K.
Land, M. L. (1960). "Expansion Turbines and Engines for Low Temperature Processing."
Adv. Cryog. Eng. 5, 250-260.
Leger, P., Thomas, P., and Zemati, C. (1971). "Design of Valveless-Type Piston Expanders
with Variable Clearance for a 1.5 W, 4.5 K Helium Refrigerator." Proc. Third Int. Cryog.
Eng. Con{., Berlin, pp. 259-263, Iliffe Sci. and Tech. Pub., Guildford, U.K.
Long, H. M., and Simon, F. E. (1954). Appl. Sci. Res. 4, 237.
Maddocks, F. E. (1968). "Applications of Turbomachinery to Small-Capacity Closed Cycle
Cryogenic Systems." Adv. Cryog. Eng. 13, 463-473.
Malkov, M.P., and Danilov, I. B. (1968). "Cryogenic Equipment of the Institute for Physical
Problems." Proc. First Int. Cryog. Eng. Con{., Tokyo, pp. 225-228, Heywood-Temple
Ind. Pub., London.
Matsubara, Y., Ishizaki, Y., and Oshima, K. (1968). "New Type Expansion Engine for
Refrigeration." Proc. First Cryog. Eng. Con{., Tokyo, pp. 210-211, Heywood-Temple
Ind. Pub., London.
Meier, R. H., and Currie, R. B. (1968). "A 4 K Single Engine CyCle Helium Refrigerator."
Adv. Cryog. Eng. 13, 441-449.
Morain, W. A., and Homes, J. W. (1963). "An Analysis of the Performance of Large
Reciprocating Expansion Engines with the Aid of a Computer and Laboratory
Prototype." Adv. Cryog. Eng. 8, 228-235.
Claude and Joule-Brayton Systems 353
A. D. Little Company, [I] 10, 12, 14,237-240 Energy Research and Generation
AEG-Telefunken A. G., [I] 14 Incorporated (ERG), [I] 20
AiResearch Manufacturing Company, [I]
16, 199-202; [2] 387
Air Products and Chemical Incorporated
Fairchild Space and Electronics Company,
(APCI), [I] 15, 77, 238, 240, 293
[I] 150
American Motors, [I] Ill
Finkelstein, Theodore, [I] 109, 131, 142, 147
Flight Dynamics Laboratories, [I) 13
Ford Motor Company, [11 Ill
Baumann Institute of the Moscow High
Franchot, Charles Louis, [I] 109
Technological School, [2) 262
Beale, William, [I) 52, 166
British Oxygen Company, [I) 8, 19
Bush, Vannevar, [I) 188 General Motors, [I) Ill
Gifford, William, [I) 14, 237-263
Gorrie, John, [I] 8
Carnot, Sadi, [I] 39
Chellis, Fred, [I] 171
Claude, Georges, [I) 8, 9, 322, 323; [2] 375
Hampson, W, [I) 8
Collins, Samuel, [I] 9, 10, 326-329; [2] 376
Harwell Atomic Energy Establishment, [I]
Cowans, Ken, [I) 14
108
Crummett, Charles, [I) 9
Herschel, John, [I) 6, 95
Crummett, Orin, [I] 9
Heylandt, D., [I) 8
Cryogenic Data Center, U.S., [2] 389
Higa, Walter,[!) 16,171
Cryogenic Society of America, [2] 391
Horn. Stuart, [I] 171
Cryogenic Technology Incorporated (CTI),
Hughes Aircraft Company, [I] 14, 103,
[!) 10, 12, 14, 15, 100, 103, 238; [2] 387
191-193; [2) 387
Cryomech Incorporated, [I] 14, 238
Hughes Santa Barbara Research Center, [I]
14
Hymatic Engineering Limited, [I) 19,
Daniels, A., [I] 101
288-290
Daunt, John, [I] 14
Davis, Harvey, [I] 9
Defense Documentation Center (DDC),
U.S., [2] 388 International Institute of Refrigeration, [2]
du Pre, F. K., [I] 96, 101 391-392
3SS
356 Name Index
357
358 Subject Index
Oil flooding, [2] 380 Precooling, [I] 275-278, 312-314; [2] 380
Open-cycle refrigeration, [2] 130 Pressure, [I] 30-31
Optimization of design parameters, [I] Pressure drop, [I] 127-128; [2] 5, 24, 380
153-157, 255-257; [2] 266-272 Pressure excursion, [2] 380
Optimum recirculation fraction, [I] 306-307 Pressure generator, [I] 49-50
Organ nodal analysis, [I] 151 Pressure oscillation, [2] 305
Orthohydrogen, [2] 143 Pressure ratio, [I] 186; [2] 380
Oscillatory bearings, [2] 79, 80 Pressure-volume ( P- V) diagram, [I] 36-37
Oscillatory flow, [2] 4, 26-30, 41-43, 302-305 Prime mover, [I] 48-49; [2] 380
Oscillatory temperature, [2] 302-305 Publications, cryogenic engineering, [2]
Outer cascade cycle, [I] 279 262-264, 387-394
Overheating, [I] 228 Pulse-width modulation, [1]228
Pulsed refrigeration system, [2] 308-313
Pump work, [2] 5
Parahydrogen, [2] 143 Pumping loss, [1]221-222
Parallel flow, [2] 15
Parallel multiple-expansion engine
arrangement, [I] 318-322 Radiation, [I] 257
Paramagnetic materials, [2] 218-226 Rallis adiabatic regenerative cycle, [1]53,
Parametric effects, [I] 71-73 62-70
Peltier heat, [2] 146 Rallis isothermal regenerative cycle, [1]53-58
Pentaerythrityl fluoride, [2] 154 Rallis thermodynamic cycle, [I] 53; [2] 380
Perfect dynamic balance, [2] 77-79 RCA Vuilleumier cryocoolers, [I] 202-206
Perforated plate exchanger, [2] 4, 11-12 Reciprocating compressors, [1]347-349; [2]
Phase angle, [I] 156, 227; [2] 380 91
Phase equilibrium, [2] 136 Reciprocating cooling machines, [I] 2, 3, 4;
Philips nodal analysis, [I] 151 [2] 380
Philips Vuilleumier cryocoolers, [I] 195-199 Reciprocating expansion engines, [I]
Phonon drag effect, [2] 152-153 322-324 [2) 298-299
Phonons, [2] 137, 147, 152-153 Reciprocating masses, [2] 71-79
Photon cooling systems, [2] 168-171 Recuperative cycles, [I] 78-94
Piston crosshead system, [2] 94 Recuperative heat exchangers, [ 1]5, 78, 257,
Piston-displacer engine, [I] 50-53 346; [2] 1-30, 32, 380
Piston-displacer single-acting Stirling Recuperative system analysis, [1]85-87
engines, [I] I 06 Reduced length, [2] 39-41
Piston leakage, [I] 339 Reduced period, [2] 39-41
Piston motion, [I] 117-124, 243 Redundant units, [I] 347
Piston rings, [I] 324; [2] 92, 274 Refrigerant, [2] 139-143
Piston seals, [I] 169 Refrigeration, [I] 256
Piston side thrust, [2] 93-95 dilution refrigeration, [I] 19
Pistons, [I] I 05, 325-331; [2] 380 electrocaloric refrigeration, [1]19
Plate-fin exchangers, [2] 4, 9-11 magnetocaloric refrigeration, [1]19
Plated tube heat exchanger, [2] 23-24 Refrigeration capacity, [1]2, 128, 154,
Polytetrafluoroethylene (PTFE), [I] 749; [2] 157-159; [2] 380
87-88, 92 Refrigeration load, [2] 380
Polytropic process, [I] 87, 298-300 Refrigeration loss, [I] 250
Pomeranchuk refrigeration, [2] 144, 211-216 Refrigeration quality, [IJ I
Porosity, [2] 121, 380 Regenerative annulus, [I] 131; [2] 381
Postle cryocoolers, [I] 237-238, 261-263; [2] Regenerative cycles, [ 1]44-77; [2] 381
380 Regenerative displacer, [I] 233
Practical regenerative cycle, [I] 123-131 Regenerative heat exchangers, [I] 5, 6, 257;
Precooler heat exchanger, [I] 275-277 [2] I, 3, 30-54, 381
Subject Index 363
Void volumes, [I] 217-218; [2] 383 Vuilleumier cycle liquefiers, [I] 26
Volume, [1]30
Volume compression ratio, [2] 383
Volume variations, [1]5 Wall effect, [2] 18
Volumetric efficiency, [1]89 Water cooling, [2] I 05
Vuilleumier cooling engines, [1]5, 14, Wear, [2] 66, 93
185-233; [2] 26, 383 Weld joints, [2] 121
accessories and components, [1]228-229 Werkspoor cryocooler, [1]176-179
power input section, [I] 209-210 Wobble-plate, [I] Ill
refrigeration section, [I] 210-211 Work diagrams, [I] 124-126
Vuilleumier cycle, [ 1]75-76, 191, 206-220; [2] Working fluid, [1]175; [2]96-98, 383
116 Working space, [2] 383