Cryogenic Engineering
Cryogenic Engineering
ENGINEERING
Second Edition
Revised and Expanded
Thomas M. Flynn
CRYOCO, Inc.
Louisville, Colorado, U.S.A.
ISBN: 0-8247-5367-4
Headquarters
Marcel Dekker, 270 Madison Avenue, New York, NY 10016, U.S.A.
tel: 212-696-9000; fax: 212-685-4540
The publisher offers discounts on this book when ordered in bulk quantities. For more infor-
mation, write to Special Sales=Professional Marketing at the headquarters address above.
10 9 8 7 6 5 4 3 2 1
To Rita
Dr. Flynn’s Cryogenic Engineering 2nd Edition was written for a specific audience,
namely, the professional engineer or physicist who needs to know some cryogenics
to get his or her job done, but not necessarily make a career of it. The 2nd Edition
was written to follow closely the cryogenic engineering professional course given
annually by Dr. Flynn for 25 years, and accordingly has been thoroughly tested
to be very practical. This 2nd (and last) edition includes over 125 new literature
citations, and features more than 130 new graphs and tables of data, which may
no longer be available elsewhere.
This book is for the engineer and scientist who work for a living and who need cryo-
genics to get a job done. It is a deskbook, containing hundreds of tables and chart of
cryogenic data that are very hard to come by. Examples and sample calculations of
how to use the data are included. It is not a textbook. Instead, it assumes that the
reader already has basic engineering and science skills. It is practical, using the mea-
surement units of trade—SI, U.S. customary, and hybrid systems—just as they are
commonly used in practice. It is not a design text, but it does contain many useful
design guidelines for selecting the right system, either through procurement or
in-house construction. In short, it is the cryogenics book I would like to have on
my desk.
This book was written to gather into one source much of the technology devel-
oped at the National Bureau of Standards (NBS) Cryogenic Engineering Laboratory
in Boulder, COL, over the last 40 years.
In the early 1950s, there was a need for the rapid development of a liquid
hydrogen technology, and the major responsibility for the progress of this new engi-
neering specialty was entrusted to the Cryogenics Section of the Heat and Power
Division of NBS in Washington, D.C. Russell Scott led the work as chief of that sec-
tion. Scott soon became the individual immediately in charge of the design and con-
struction in Boulder (in March 1952) of the first large-scale liquefier for hydrogen
ever built. This was the beginning of the Cryogenic Engineering Laboratory of NBS.
Again, in the late 1950s—when the nation was striving to regain world leader-
ship in the exploration of space—the NBS Cryogenics Engineering Laboratory,
which had by then matured under Scott’s leadership, assumed a pre-eminent role
in the solution of problems important to the nation. Scott, having had the foresight
to establish a Cryogenic Data Center within the laboratory, was able to provide a
focal point for information on many aspects of cryogenic engineering.
As a result of all this pioneering in the field of low-temperature engineering, a
considerable amount of valuable technology was developed that in the course of nor-
mal events might have been lost. Scott recognized this, and the result was his book
Cryogenic Engineering, an important first in its field. Its quality, authority, and com-
pleteness constitute a lasting tribute to him.
This present book is a mere shadow of Scott’s work but is intended once more
to up-date and preserve some of the cryogenics developed at NBS. There is only one
author’s name on the cover. Nonetheless, this book is a product of the collaboration
of hundreds of good men and women of the NBS Cryogenic Engineering Labora-
tory. It is written to preserve the technology they developed.
vii
2. Basic Principles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
1. Introduction . . . . 23
2. Thermodynamics . . . . 27
3. Heat Transfer . . . . 55
4. Momentum Transfer . . . . 68
5. Heat Leak and Pressure Drop in Cryogenic Transfer Lines . . . . 69
6. Cooldown . . . . 73
7. Summary . . . . 76
3. Cryogenic Fluids . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
1. PVT Behavior of a Pure Substance . . . . 77
2. Temperature–Enthalpy and Temperature–Entropy Diagrams of Pure
Substances . . . . 81
3. Properties and Uses of Cryogenic Fluids . . . . 83
ix
7. Insulation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 445
1. Introduction . . . . 445
2. Heat Transfer . . . . 446
3. Vacuum Insulation . . . . 448
4. Evacuated Porous Insulation . . . . 468
5. Gas-Filled Powders and Fibrous Materials . . . . 475
6. Solid Foams . . . . 477
7. Multilayer Insulation . . . . 483
8. Vapor Barriers . . . . 510
9. Protective Enclosures . . . . 513
10. Liquid and Vapor Shields . . . . 514
11. Composite Insulations . . . . 518
12. Other Materials . . . . 525
13. Placement of Insulation Systems . . . . 525
14. Adhesives . . . . 526
References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 873
1
Cryogenic Engineering Connections
1. FOREWARNED
In his masterpiece, Civilization, Sir Kenneth Clark writes not a foreword but rather a
fore-warned. Clark writes to the effect that he does not know if his recounting is the
actual history of civilization or not. He only knows that his story is his own view of
the history of civilization. Having said that, Clark found much relief from the
tedium of accuracy and was able to tell a more entertaining, and still highly accurate,
story.
I wish to do the same. I do not know if the recounting to follow in this chapter
is accurate or not. I only know that it is my own view of how it happened. I hope it is
accurate—I did not deliberately make any of it up—but who knows?
I am an American and unabashedly proud of that fact. Therefore, let me begin
with the American who may have gotten it all going, John Gorrie (1803–1855).
2. THE ENTREPRENEURS
John Gorrie was born in Charleston South Carolina in 1803 and graduated from
the College of Physicians and Surgeons in New York City in 1833. In the spring of
1833 Gorrie moved to Apalachicola, a small coastal town in Florida situated on
the Gulf of Mexico at the mouth of the river after which the town is named. By
the time he arrived it was already a thriving cotton port, where ships from the
north-east arrived to unload cargoes of supplies and load up with cotton for the
northern factories. Within a year Gorrie was involved with town affairs. He served
as mayor, city treasurer, council member, bank director, and founder of Trinity
Church.
In 1834, he was made postmaster and in 1836 president of the local branch of
the Pensacola Bank. In the same year, the Apalachicola Company asked him to
report on the effects of the climate on the population, with a view to possible expan-
sion of the town. Gorrie recommended drainage of the marshy, low-lying areas that
surrounded the town on the grounds that these places gave off a miasma com-
pounded by heat, damp and rotting vegetation which, according to the Spallanzani
theory with which every doctor was intimate, carried disease. He suggested that only
brick buildings be erected. In 1837, the area enjoyed a cotton boom, and the town
population rose to 1500. Cotton bales lined the streets, and in four months 148 ships
arrived to unload bricks from Baltimore, granite from Massachusetts, house framing
from New York. Gorrie saw that the town was likely to grow as commerce increased
1
2 Chapter 1
and suggested that there was a need for a hospital. There was already a small medical
unit in operation under the auspices of the US Government, and Gorrie was
employed there on a part-time basis. Most of his patients were sailors and waterside
workers, and most of them had fever, which was endemic in Apalachicola every
summer.
Gorrie became obsessed with finding a cure to the disease. As early as 1836, he
came close to the answer, over 60 years ahead of the rest of the world. In that year he
wrote: ‘‘Gauze curtains, though chiefly used to prevent annoyance and suffering
from mosquitoes, are thought also to be sifters of the atmosphere and interceptors
and decomposers of malaria.’’ The suggestion that the mosquito was the disease
carrier was not to be made until 1881, many years after Gorrie’s death, and for
the moment he presumed that it came in some form of volatile oil, rising from
the swamps and marshes.
By 1838, Gorrie had noticed that malaria seemed to be connected with hot,
humid weather, and he set about finding ways to lower the temperature of his
patients in summer. He began by hanging bowls full of ice in the wards and circu-
lating the cool air above them by means of a fan. The trouble was that in Apala-
chicola ice was hard to come by. Ever since a Massachusetts merchant named
Frederic Tudor had hit on the idea of cutting ice from ponds and rivers in winter
and storing it in thick-walled warehouses for export to hot countries, regular ice
shipments had left the port of Boston for destinations as far away as Calcutta.
But Apalachicola was only a small port which the ships often missed altogether:
if the ice crop was poor, the price rose to the exorbitant rate of $1.25 a pound.
In 1844, Gorrie found the answer to the problem. It was well known that com-
pressed gases which are rapidly allowed to expand absorb heat from their surround-
ings, so Gorrie constructed a steam engine to drive a piston back and forward in a
cylinder. His machine compressed air, causing it to heat, and then the air through
radiant coils where it decompressed and cooled, absorbing heat from a bath of
brine. On the next cycle the air remained cool, since the brine had given up most
of its heat. This air was then pumped out of the cylinder and allowed to circulate
in the ward. Gorrie had invented air-conditioning. By bringing the cold brine into
contact with water, Gorrie was then able to draw heat from the water to a point
where it froze. Gorrie’s application of compressed and decompressed gas as a cool-
ant in radiant coils remains the common method for cooling air in modern refrig-
eration systems. His first public announcement of this development was made on 14
July 1850 in the Mansion House Hotel, where M. Rosan, the French Consul in
Apalachicola, was celebrating Bastille Day with champagne. No ice ship had
arrived, so the champagne was to be served warm. At the moment of the toast
to the French Republic four servants entered, each carrying a silver tray on which
was a block of ice the size of a house-brick, to chill the wine, as one guest put it,
‘‘by American genius’’.
In May of the following year Gorrie obtained a patent for the first ice-making
machine (see Fig. 1.1), the first patent ever issued for a refrigeration machine. The
patent specified that the water container should be placed in the cylinder, for faster
freezing. Gorrie was convinced his idea would be a success. The New York Times
thought differently: ‘‘There is a crank,’’ it said, ‘‘down in Apalachicola, Florida,
who claims that he can make ice as good as God Almighty!’’ In spite of this, Gorrie
advertised his invention as ‘‘the first commercial machine to work for ice making and
refrigeration.’’ He must have aroused some interest, for later he was in New Orleans,
selling the idea that ‘‘a ton of ice can be made on any part of the Earth for less than
$2.00.’’ But he was unable to find adequate backing, and in 1855 he died, a broken
and dispirited man. A statue of John Corrie now stands in Statuary Hall of the US
Capitol building, a tribute from the State of Florida to his genius and his importance
to the welfare of mankind.
Three years after his death a Frenchman, Ferdinand Carré, produced a
compression ice-making system and claimed it for his own, to the world’s acclaim.
Carré was a close friend of M. Rosan, whose champagne had been chilled by
Gorrie’s machine eight years before.
Just before he died, Gorrie wrote an article in which he said: ‘‘The system is
equally applicable to ships as well as buildings . . . and might be instrumental in pre-
serving organic matter an indefinite period of time.’’ The words were prophetic,
because 12 years later Dr. Henry P. Howard, a native of San Antonio, used the
air-chilling system aboard the steamship Agnes to transport a consignment of frozen
beef from Indianola, Texas, along the Gulf of Mexico to the very city where Gorrie
had tried and failed to get financial backing for his idea. On the morning of Saturday
10 June 1869, the Agnes arrived in New Orleans with her frozen cargo. There it was
served in hospitals and at celebratory banquets in hotels and restaurants. The New
Orleans Times Picayune wrote: ‘‘[The apparatus] virtually annihilates space and
laughs at the lapse of time; for the Boston merchant may have a fresh juicy beefsteak
from the rich pastures of Texas for dinner, and for dessert feast on the delicate,
luscious but perishable fruits of the Indies.’’
At the same time that Howard was putting his cooling equipment into the
Agnes, committees in England were advising the government that mass starvation
was likely in Britain because for the first time the country could no longer feed itself.
Between 1860 and 1870, consumption of food increased by a staggering 25%. As the
population went on rising, desperate speeches were made about the end of
democracy and nationwide anarchy if the Australians did not begin
4 Chapter 1
immediately to find a way of sending their sheep in the form of meat instead of
tallow and wool.
Thomas Malthus (1766–1834), the first demographer, published in An Essay
on the Principle of Population (1798). According to Malthus, population tends to
increase faster than the supply of food available for its needs. Whenever a relative
gain occurs in food production over population growth, a higher rate of population
increase is stimulated; on the other hand, if population grows too much faster than
food production, the growth is checked by famine, disease, and war. Malthus’s the-
ory contradicted the optimistic belief prevailing in the early 19th century, that a
society’s fertility would lead to economic progress. Malthus’s theory won supporters
and was often used as an argument against efforts to better the condition of the
poor. (the poor should die quietly). Those (the vast majority) who did not read
the complete essay, assumed that mass starvation and imminent and inevitable, lead-
ing to the so-called Malthusian revolution in England. There were food riots in the
streets. Food storage and transportation was seen as a critical global issue facing
mankind.
Partially for these reasons, two Britons, Thomas Mort and James Harrison,
emigrated to Australia and set up systems to refrigerate meat. In 1873 Harrison
gave a public banquet of meat that had been frozen by his ice factory, to cele-
brate the departure of the S.S. Norfolk for England. On board were 20 tons of
mutton and beef kept cold by a mixture of ice and salt. On the way, the system
developed a leak and the cargo was ruined. Harrison left the freezing business.
Mort then tried a different system, using ammonia as the coolant. He too gave
a frozen meat lunch, in 1875, to mark the departure for England of the S.S.
Northam. Another leak ruined this second cargo, and Mort retired from the busi-
ness. But both men had left behind them working refrigeration plants in Austra-
lia. The only problem was to find the right system to survive the long voyage to
London. Eventually, the shippers went back to Gorrie’s ‘‘dry air’’ system. Aboard
ship, it was much easier to replace leaking air than it was to replace leaking
ammonia. Even though ammonia was ‘‘more efficient’’, it was sadly lacking while
air was ubiquitous. This NH3=air substitution is one early example of a lack of
‘‘systems engineering’’. There are often unintended consequences of technology,
and the working together of the system as a whole must be considered. NH3
refrigerators were thermodynamically more efficient, but could not be relished
with NH3 readily.
3. THE BUTCHERS
Figure 1.2 The S. S. Strathleven, carrying the first successful consignment of chilled beef
from Australia to England. Note the cautious mixture of steam and sail, which was to
continue into the 20th century.
4. THE BREWERS
Figure 1.3 Unloading frozen meat from Sydney, Australia, at the South West India Dock,
London. This shows the hold of the Catania, which left port in August 1881 with 120 tons of
meat from the same exporters who had filled the Strathleven.
6 Chapter 1
produced by a method using a yeast which ferments on the surface of the beer vat
over a period of 5–7 days, when the ideal ambient temperature range is from 60 F
to 70 F. Beer brewed in this way suffers less from temperature changes while it is
being stored, and besides, Britain rarely experiences wide fluctuations in summer
temperatures. But in Germany beer is produced by a yeast which ferments on the
bottom of the vat. This type of fermentation may have been introduced by monks
in Bavaria as early as 1420 and initially was an activity limited to the winter months,
since bottom fermentation takes place over a period of up to 12 weeks, in an ideal
ambient temperature of just above freezing point. During this time the beer was
stored in cold cellars, and from this practice came the name of the beer: lager, from
the German verb lagern (to store). From the beginning there had been legislation in
Germany to prevent the brewing of beer in the summer months, since the higher tem-
peratures were likely to cause the production of bad beer. By the middle of the 19th
century every medium-sized Bavarian brewery was using steam power, and when the
use of the piston to compress gas and cool it became generally known, the president
of the German Brewer’s Union, Gabriel Sedlmayr of the Munich Spätenbrau brew-
ery, asked a friend of his called Carl Von Linde if he could develop a refrigerating
system to keep the beer cool enough to permit brewing all the year round. Von Linde
solved Sedlmayr’s problem and gave the world affordable mechanical refrigeration,
an invention that today is found in almost every kitchen.
Von Linde used ammonia instead of air as his coolant, because ammonia lique-
fied under pressure, and when the pressure was released it returned to gaseous form,
and in so doing drew heat from its surroundings. In order to compress and release
the ammonia, he used Gorrie’s system of a piston in a cylinder. Von Linde did
not invent the ammonia refrigeration system, but he was the first to make it work.
In 1879 he set up laboratories in Wiesbaden to continue research and to convert
his industrial refrigeration unit into one for the domestic market. By 1891, he had
put 12,000 domestic refrigerators into German and American homes. The modern
fridge uses essentially the same system as the one with which Von Linde chilled
the Spätenbrau cellars.
5. THE INDUSTRIALISTS
Interest in refrigeration spread to other industries. The use of limelight, for instance,
demanded large amounts of oxygen, which could be more easily handled and trans-
ported in liquid form. The new Bessemer steel-making process used oxygen. It may
be no coincidence that an ironmaster was involved in the first successful attempt to
liquefy the gas. His name was Louis Paul Cailletet, and together with a Swiss
engineer, Raoul Pictet, he produced a small amount of liquid oxygen in 1877.
At the meeting of the Académie des Sciences in Paris on 24 December 1887,
two announcements were made which may be recognized as the origins of cryogenics
as we know it today. The secretary to the Académie spoke of two communications he
had received from M. Cailletet working in Paris and from Professor Pictet in Geneva
in which both claimed the liquefaction of oxygen, one of the permanent gases.
The term ‘‘permanent’’ had arisen from the experimentally determined fact
that such gases could not be liquefied by pressure alone at ambient temperature,
in contrast to the nonpermanent or condensable gases like chlorine, nitrous oxide
and carbon dioxide, which could be liquefied at quite modest pressures of
30–50 atm. During the previous 50 years or so, in extremely dangerous experiments,
8 Chapter 1
pursue the liquefaction of oxygen (he made a name developing ice-skating rinks), his
use of the cascade system inspired others like Kamerlingh Onnes and Dewar.
In the early 1880s one of the first low-temperature physics laboratories, the
Cracow University Laboratory in Poland, was established by Szygmunt von
Wroblewski and K. Olszewski. They obtained liquid oxygen ‘‘boiling quietly in a test
tube’’ in sufficient quantity to study properties in April 1883. A few days later, they
also liquefied nitrogen. Having succeeded in obtaining oxygen and nitrogen as true
liquids (not just a fog of liquid droplets), Wroblewski and Olszewski, now working
separately at Cracow, attempted to liquefy hydrogen by Cailletet’s expansion tech-
nique. By first cooling hydrogen in a capillary tube to liquid-oxygen temperatures
and expanding suddenly from 100 to 1 atm, Wroblewski obtained a fog of liquid-
hydrogen droplets in 1884, but he was not able to obtain hydrogen in the completely
liquid form.
The Polish scientists at the Cracow University Laboratory were primarily inter-
ested in determining the physical properties of liquefied gases. The ever-present prob-
lem of heat transfer from ambient plagued these early investigators because the
cryogenic fluids could be retained only for a short time before the liquids boiled
away. To improve this situation, an ingenious experimental technique was developed
at Cracow. The experimental test tube containing a cryogenic fluid was surrounded
by a series of concentric tubes, closed at one end. The cold vapor arising from the
liquid flowed through the annular spaces between the tubes and intercepted some
of the heat traveling toward the cold test tube. This concept of vapor shielding is used
today in conjunction with high-performance insulations for the long-term storage of
liquid helium in bulk quantities.
All over Europe scientists worked to produce a system that would operate to
make liquid gas on an industrial scale. The major problem in all this was to prevent
the material from drawing heat from its surroundings. In 1882, a French physicist
called Jules Violle wrote to the French Academy to say that he had worked out a
way of isolating the liquid gas from its surroundings through the use of a vacuum.
It had been known for some time that vacua would not transmit heat, and Violle’s
arrangement was to use a double-walled glass vessel with a vacuum in the space
between the walls. Violle has been forgotten, his place taken by a Scotsman who
was to do the same thing, much more efficiently, eight years later. His name was
Sir James Dewar, and he added to the vessel by silvering it both inside and out
(Violle had only silvered the exterior), in order to prevent radiation of heat either
into or out of the vessel.
6. THE SCIENTISTS
10 Chapter 1
Figure 1.6 Sir James Dewar lecturing at the Royal Institution. Although Violle preceded
Dewar in the development of the vacuum flask, there is no evidence that Dewar knew of
his work when he presented the details of his new container to the Royal Institution in 1890.
liquids. Meanwhile his work on the absorptive capacity of charcoal at low tempera-
tures paved the way towards the development of all-metal, double-walled vacuum
vessels.
Another first for Dewar was his use of mixtures of gases to enhance J-T
cooling, a topic revisited in cryogenics in 1995.
Historically, the expansion of a mixture of hydrogen and nitrogen was
employed by Dewar in 1894 in attempts to liquefy hydrogen at that time. Dewar
wrote: ‘‘Expansion into air at one atmosphere pressure of a mixture of 10% nitrogen
in hydrogen yielded a much lower temperature than anything that has been recorded
up to the present time.’’
Because his flasks (unlike those of Jules Violle) were silvered inside and out,
Dewar’s flasks could equally well retain heat as cold. Dewar’s vessel became known
in scientific circles as the Dewar flask; with it, he was able to use already liquid gases
to enhance the chill during the liquefaction of gases whose liquefaction temperature
was lower than that of the surrounding liquid. In this way, in 1891, he succeeded for
the first time in liquefying hydrogen.
Dewar had considerable difficulty in finding competent glass makers willing to
undertake the construction of his double-walled vessels and had been forced to get
them made in Germany.
By 1902 a German called Reinhold Burger, whom Dewar had met when visit-
ing Germany to get his vessels made, was marketing them under the name of
Thermos. The manufacture of such vessels developed into an important industry
in Germany as the Thermos flask, and this monopoly was maintained up to 1914.
Dewar never patented his silvered vacuum flask and therefore never benefited finan-
cially from his invention.
The word ‘‘Cryogenics’’ was slow in coming.
The word cryogenics is a product of the 20th Century and comes from the
Greek—kroB—frost and—ginomai—to produce, engender. Etymologically, cryo-
genics means the science and art of producing cold and this was how Kamerlingh
Onnes first used the word in 1894. Looking through the papers and publications
of Dewar and Claude, it appears that neither of them ever used the word; indeed
a summary of Dewar’s achievements by Armstrong in 1916 contains no mention
but introduces yet another term, ‘‘the abasement of temperature.’’
In 1882, Kamerlingh Onnes (1853–1926) was appointed to the Chair of Experi-
mental Physics at the University of Leiden in the Netherlands and embarked on
building up a low-temperature physics laboratory in the Physics Department. The
inspiration for his laboratory was provided by the systematic work of Van der Waals
at Amsterdam, and subsequently at Leiden, on the properties of gases and liquids. In
1866, Van der Waals had published his first paper on ‘‘the continuity of liquid and
gaseous states’’ from which the physical understanding of the critical state, and of
liquefaction and evaporation was to grow.
He operated an open-door policy encouraging visitors from many countries to
visit, learn, and discuss. He published all the experimental results of his laboratory,
and full details of the experimental apparatus and techniques developed, by introdu-
cing in 1885 a new journal ‘‘Communications from the Physical Laboratory at the
University of Leiden.’’
As a result, he developed a wide range of contacts and a growing track record
of success, so that from the turn of the century his Leiden laboratory held a leading
position in cryogenics for almost 50 years, certainly until the mid-1930s. In 1908, for
example, he won the race with Dewar and others to liquefy helium and went on to
discover superconductivity in 1911.
Onnes’ first liquefaction of helium in 1908 was a tribute both to his experimen-
tal skill and to his careful planning. He had only 360 L of gaseous helium obtained
by heating monazite sand from India. More than 60 cm3 of liquid helium was pro-
duced by Onnes in his first attempt. Onnes was able to attain a temperature of
1.04 K in an unsuccessful attempt to solidify helium by lowering the pressure above
a container of liquid helium in 1910.
It is interesting to compare Dewar’s approach with that of his rival,
Kamerlingh Onnes. Having successfully liquefied hydrogen in 1898, Dewar had been
able to monopolize the study of the properties of liquid hydrogen and published
many papers on this subject. His attempts in 1901 to liquefy helium in a Cailletet
tube cooled with liquid hydrogen at 20.5 K, using a single isentropic expansion from
pressures up to 100 atm, led at first to a mist, being clearly visible. Dewar was
12 Chapter 1
7. THE ENGINEERS
Dr. Hampson was medical officer in charge of the Electrical and X-ray Departments,
Queens and St. John’s Hospitals, Leicester Square, London. He was a product of the
Victorian age, with a classics degree at Oxford in 1878, having subsequently acquired
his science as an art and a living; but he possessed an extraordinary mechanical flair.
He was completely overshadowed by Dewar at The Royal Institution, Ramsay at
University College and Linde at Wiesbaden, each with their considerable laboratory
facilities. Indeed, Dewar seems to have been unable to accept Hampson as a fellow
experimentalist or to acknowledge his contribution to cryogenics.
And yet, Hampson with his limited facilities was able to invent and develop a
compact air liquefier which had a mechanical elegance and simplicity which made
Dewar’s efforts crude and clumsy in comparison. Indeed, Hampson’s design of heat
exchanger was so successful that it is still acknowledged today. The Hampson type
coiled tube heat exchanger is widely used today. In 1895, when Hampson and Linde
independently took out patents on their designs of air liquefier, the Joule–Thomson
effect was known and the principle of recuperative cooling by so-called self-
intensification had been put forward by Siemens as early as 1857. The step
forward in both patents was to break away from the cascade system of cooling
and to rely entirely on Joule–Thomson cooling together with efficient heat exchanger
designs.
Figure 1.7 Linde two-stage compressor and 3 L=hr air liquefier (1895).
14 Chapter 1
Figure 1.8 Tripler’s air liquefier with steam-driven compressor and tube-in-shell heat
exchangers (1898).
caused the gas to cool isentropically rather than isenthalpically, as is the case in
Joule–Thomson expansion. Claude’s first engines were reciprocating engines using
leather seals (actually, the engines were simply modified steam engines). During
the same year, Claude established l’Air Liquide to develop and produce his systems.
The increase in cooling effect over the Joule–Thomson expansion of the
Linde=Hampson=Tripler designs was so large as to constitute another technical
breakthrough. Claude went on to develop air liquefiers with piston expanders in
the newly formed Société L’Air Liquide.
Although cryogenic engineering is considered a relatively new field in the US, it
must be remembered that the use of liquefied gases in US industry began in the early
1900s. Linde installed the first air-liquefaction plant in the United States in 1907, and
the first American-made air-liquefaction plant was completed in 1912. The first com-
mercial argon production was put into operation in 1916 by the Linde company in
Cleveland, Ohio. In 1917 three experimental plants were built by the Bureau of
Mines, with the cooperation of the Linde Company, Air Reduction Company,
and the Jefferies-Norton Corporation, to extract helium from natural gas of Clay
County, Texas. The helium was intended for use in airships for World War I.
Commercial production of neon began in the United States in 1922, although Claude
had produced neon in quantity in France since 1920.
Around 1947 Dr. Samuel C. Collins of the department of mechanical engineer-
ing at Massachusetts Institute of Technology developed an efficient liquid-helium
laboratory facility. This event marked the beginning of the period in which liquid-
helium temperatures became feasible and fairly economical. The Collins helium
cryostat, marketed by Arthur D. Little, Inc., was a complete system for the safe,
economical liquefaction of helium and could be used also to maintain temperatures
at any level between ambient temperature and approximately 2 K.
The first buildings for the National Bureau of Standards Cryogenic Engineering
Laboratory were completed in 1952. This laboratory was established to provide engi-
neering data on materials of construction, to produce large quantities of liquid hydro-
gen for the Atomic Energy Commission, and to develop improved processes and
equipment for the fast-growing cryogenic field. Annual conferences in cryogenic engi-
neering have been sponsored by the National Bureau of Standards (sometimes spon-
sored jointly with various universities) from 1954 (with the exception of 1955) to 1973.
At the 1972 conference at Georgia Tech in Atlanta, the Conference Board voted to
change to a biennial schedule alternating with the Applied Superconductivity Confer-
ence. The NBS Cryogenic Engineering Laboratory is now part of history, as is the
name ‘‘NBS’’, now the National Institutes for Science and Technology (NIST).
The impact of cryogenics was wide and varied. The Dewar flask changed the social
habits of the Edwardian well-to-do: picnics became fashionable because of it. In time
it changed the working-man’s lunch break and accompanied expeditions to the tropics
and to the poles, carrying sustenance for the explorers and returning with hot or cold
specimens. Later, it saved thousands of lives by keeping insulin and other drugs from
going bad. Perhaps its most spectacular impact was made, however, by two men whose
work went largely ignored, and by a third who did his work in a way that could not be
ignored. The first was a Russian called Konstantin Tsiolkovsky, whose early use of
liquid gas at the beginning of the 20th century was to lie buried under governmental
lack of interest for decades. The second was an American called Robert Goddard,
who did most of his experiments on his aunt’s farm in Massachusetts, and whose only
reward was lukewarm interest from the weather bureau.
On 16 March 1926, Dr. Robert H. Goddard conducted the world’s first suc-
cessful flight of a rocket powered by liquid-oxygen–gasoline propellant on a farm
near Auburn, Massachusetts. This first flight lasted only 212 sec, and the rocket
reached a maximum speed of only 22 m=s (50 mph). Dr. Goddard continued his
work during the 1930s, and by 1941 he had brought his cryogenic rockets to a fairly
high degree of perfection. In fact, many of the devices used in Dr. Goddard’s rocket
systems were used later in German V-2 weapons systems (Fig 1.9).
The third was a German, Herman Oberth, and his work was noticed because it
aimed at the destruction of London.
His liquid gases were contained in a machine that became known as Vengeance
Weapon 2, or V-2, and by the end of the Second World War it had killed or injured
thousands of military and civilian personnel. All three men had realized that certain
gases burn explosively, in particular hydrogen and oxygen, and that, since in their
liquid form they occupy less space than as a gas (hydrogen does this by a factor
of 790) they were an ideal fuel. Thanks to the principles of the Dewar flask, they
could be stored indefinitely, transported without loss, and contained in a launch
vehicle that was essentially a vacuum flask with pumps, navigation systems, a
combustion chamber and a warhead.
(One of Oberth’s most brilliant assistants was a young man called Werner von
Braun, and it was he who brought the use of liquid fuel to its most spectacular
expression when his brainchild, the Saturn V, lifted off at Cape Canaveral on 16 July
1969, carrying Armstrong, Aldrin, and Collins to their historic landing on the
moon.)
16 Chapter 1
Figure 1.9 The V-2 liquid-fueled rocket used a mixture of oxygen and kerosene. Originally
developed at the experimental rocket base in Peenemunde, on the Baltic, the first V-2 landed in
London in 1944. When the war ended, German engineers were working on a V-3 capable of
reaching New York.
significance of his discovery, Kamerlingh Onnes was totally frustrated by the lack of
any practical application because he found that quite small magnetic fields, applied
either externally or as self-fields from internal electric currents, destroyed the
superconducting state. The material suddenly quenched and acquired a finite electri-
cal resistance.
For the next 50 years, very little progress was made in applying superconduc-
tivity, although there was a growing realization that a mixed state of alternate lami-
nae of superconducting and normal phase had a higher critical field which depended
on the metallurgical history of the sample being studied. In fact, an impasse had
arisen by the late 1950s, and little progress was being achieved on the application
of superconductivity.
On the other hand, great progress had been achieved by this time in the UK,
USSR, and US in the theoretical description of superconductivity. This progress
stemmed largely from two sets of experimental evidence; firstly, demonstrations of
the isotope effect, indicating that lattice vibrations must play a central role in the
interaction leading to superconductivity; secondly, the accumulation of evidence that
an energy gap exists in the spectrum of energy states available to the conduction elec-
trons in a superconductor. In 1956, Bardeen, Cooper, and Schrieffer (see Table 1.1.
for a list of the many Nobel Laureates indebted to cryogenics) proposed a successful
theory of superconductivity in which conduction electrons are correlated in pairs
with the same center of mass momentum via interactions with the lattice. The theory
made predictions in remarkably good agreement with experimental data and
provided the basis for later developments such as quantum mechanical tunneling,
magnetic flux quantization, and the concept of coherence length.
The break out of the impasse in the application of superconductivity came from
systematic work at the Bell Telephone Laboratories, led by Matthias and Hulm. In
1961, they published their findings on the brittle compound Nb3Sn, made by the
high-temperature treatment of tin powder contained in a Niobium capillary tube. This
compound retained its superconductivity in a field of 80 KG (8 Tesla), the highest field
available to them, while carrying a current equivalent to a density of 100,000 A=cm2.
The critical field at 4.2 K, the boiling point of helium, was much greater than 8 T and
was therefore at least 100 times higher than that of any previous superconductor.
This finding was followed by the discovery of a range of compounds and alloys
with high critical fields, including NbZr, NbTi, and V3Ga; the race was on to man-
ufacture long lengths of wire or tape and develop superconducting magnets for
commercial applications. At first, the early wires were unstable and unreliable and
it soon became clear that a major research effort was needed.
18 Chapter 1
(Continued)
20 Chapter 1
(Continued)
1973 B. D. Josephson, I. Giaever, and L. Esaki awarded the Nobel prize for discovery
of the Josephson Junction (tunneling supercurrents)
1978 Peter Kapitza receive the Nobel for the characterization of HeII as a superfluid
1987 J. G. Bednorz and K. A. Mueller awarded the Nobel Prize for discovering
high-temperature superconductors
Y–BA–Cu–O ceramic superconductors found
1996 D. Lee, D. D. Osheroff, and R. C. Richardson receive the Nobel prize of the
discovery of superfluidity in helium-3
1997 S. Chu and Claude Cohen-Tannoudji awarded the Nobel for discovering methods
to cool and trap atoms with lasers
1998 R. B. Laughlin, H. L. Stormer, and D. C. Tsui receive the Nobel for discovering a
new form of quantum fluid excitations at extremely low temperatures
2001 E. Cornell, W. Ketterle, and C. Wieman awarded the Nobel for achieving the
Bose–Einstein condensate at near absolute zero
Many Nobel laureates received their honors either because of work directly in cryo-
genics, or because of work in which cryogenic fluids were indispensable refrigerants.
Table 1.1 provides an imposing list of such Nobel Prize winners.
Two examples from this list of Nobel laureates are considered to illustrate the
application of cryogenics to education and basic research.
The first example involves Donald Glaser who invented the bubble cham-
ber. Glaser recognized the principal limitation of the Wilson cloud chamber,
namely that the low-density particles in the cloud could not intercept enough
high-energy particles that were speeding through it from the beams of powerful
accelerators. Glaser’s first bubble chamber operated near room temperature and
used liquid diethylether. To avoid the technical difficulties presented by such a
complex target as diethylether, Prof. Luis Alvarez, of the University of Califor-
nia, devised a bubble chamber charged with liquid hydrogen. Since hydrogen is
the simplest atom, consisting only of a proton and an electron, it interferes only
slightly with the high-energy processes being studied with the giant accelerators,
although it is readily ionized and serves quite well as the detector in the bubble
chamber.
The first really large liquid-hydrogen bubble chamber was installed at the
Lawrence Radiation Laboratory of the University of California and became opera-
tive in March 1959. A remark was made that this 72-in. liquid-hydrogen bubble
chamber would be the equivalent of a Wilson cloud chamber one-half mile long.
They could not really be equal, however, because the cloud chamber does not have
the great advantage of utilizing the simplest molecules as the detector.
The second example chosen from the list of Nobel laureates that serves to
link cryogenics to basic science is John Bardeen, who was cited for his work in
22 Chapter 1
2
Basic Principles
1. INTRODUCTION
In the last half of the 19th century, roughly from 1850 to the turn of the century, the
principal export of the United States was cotton. The second most important export
by dollar volume was ice. Americans developed methods for harvesting, storing, and
transporting natural ice on a global scale. In 1810, a Maryland farmer, Thomas
Moore, developed an ice-box to carry butter to market and to keep it hard until sold.
This early refrigerator was an oval cedar tub with an inner sheet of metal serving as a
butter container that could be surrounded on four sides by ice. A rabbit skin pro-
vided the insulated cover. Moore also developed an insulated box for home use. It
featured an ice container attached to the lid and a 6 cu ft storage space below. Ice
harvesting was revolutionized in 1825 when Nathaniel J. Wyeth invented a horse-
drawn ice cutter (see Fig. 2.1.) It was Wyeth’s method of cutting blocks of similar
size quickly and cheaply that made the bountiful supplies of natural ice resources
of the United States available for food preservation. A steam-powered endless chain
was developed in 1855 that could haul 600 tons of ice per hour to be stored. Wyeth
and Tudor patented a means to prevent the ice blocks from freezing together by pla-
cing saw-dust between the layers. Uniform blocks reduced waste, facilitated trans-
portation, and introduced ice to the consumer level (see Figs. 2.2 and 2.3).
The amount of ice harvested each year was staggering. In the winter of 1879–
1880, about 1,300,000 tons of ice was harvested in Maine alone, and in the winter of
1889–1890, the Maine harvest was 3 million tons. The Hudson River supplied about
2 million tons of ice per year to New York City during that period. Even that was
not enough, for during those same years, New York City imported about 15,000 tons
from Canada and 18,000 tons from Norway annually.
Ice changed the American diet. Fresh food preserved with ice was now pre-
ferred to food preserved by smoking or salting. Fresh milk could be widely distrib-
uted in the cities. By now, the people were accustomed to having ice on demand, and
great disruptions of the marketplace occurred when ice was not available. Ice was,
after all, a natural product, and a ‘‘bad’’ winter (a warm one) was disastrous to
the marketplace and to the diet. It was time for the invention of machines to make
ice on demand and free the market from dependence on the weather.
Machines employing air as a refrigerant, called compressed air or cold air
machines, were appearing. They were based on the reverse of the phenomenon of
heating that occurred when air was compressed, namely that air cooled as it
expanded against resistance. This phenomenon of nature had been observed as early
as the middle of the 18th century. Richard Trevithick, who lived until 1833 in
23
24 Chapter 2
Figure 2.1 Ice harvesting and storage in the 1850s. Accumulated snow was first removed
with a horse-drawn ice plane. Then a single straight groove was made with a hand tool.
The ice was marked off from this line into squares 22 in. on a side. Actual cutting was begun
with Wyeth’s invention and finished with hand saws. The uniform blocks of ice were pulled
across the now open water to the lift at the icehouse. Here, they were stored under an insula-
tion of straw or sawdust. The availability of a relatively dependable supply of ice all year,
combined with uniformity of size of blocks, did much to create the market for ice.
Figure 2.2 Delivering ice in New York, 1884. In this neighborhood, apparently, icemen did
not enter the kitchen but sold their wares at the curb.
Basic Principles 25
Figure 2.3 Commercial ice delivery in New York, 1884. We see delivery to a commercial
establishment, possibly a saloon. Though ice is available, the butcher a few doors down the
street still allows his meat to hang unrefrigerated in the open air.
Cornwall, England, constructed engines in which expanding air was used to convert
water to ice. In 1846, an American, John Dutton, obtained a patent for making ice
by the expansion of air. The real development of the cold air machine, however,
began with the one developed by John Gorrie of Florida and patented in England
in 1850 and in the United States in 1851.
Later, we shall see how Gorrie’s machine evolved into true cryogenics, not
merely ice production. For now, however, in the late 1800s, it was possible to buy
an ice-making machine for home use like that shown in Fig. 2.4. This primitive
refrigerator had the same components as any refrigerator of today. It also followed
exactly the same thermodynamic process as all of today’s refrigerators.
All refrigerators involve exchanging energy (work, as the young lady in Fig. 2.4
is doing) to compress a working fluid. The working fluid is later expanded against a
resistance, and the fluid cools (usually). Engineers have made it their job to deter-
mine how much cooling is produced, how much work is required for that cooling,
and what kind of compressor or other equipment is required. More than that, we
have made it our business to optimize such a cycle by trading energy for capital until
a minimum cost is found. We can make any process more ‘‘efficient’’, that is, less
energy-demanding, by increasing the size, complexity, and cost of the capital equip-
ment employed. For instance, an isothermal compressor requires less energy to
operate than any other kind of compressor. An isothermal compressor is also the
most expensive compressor to build, because an infinity of intercooling stages are
required for true isothermal compression.
As another example of trading energy costs for capital costs, consider the
Carnot cycle (described later), which is the most thermodynamically efficient heat
engine cycle possible. An essential feature of the Carnot engine is that all heat
26 Chapter 2
Figure 2.4 Home ice machine. Machines such as this were available from catalog centers in
the late 1880s. The figure illustrates that energy (work) and capital equipment are at the heart
of all cryogenic processes.
The science that deals with these trade-offs between energy and capital is ther-
modynamics. Thermodynamics is about money. We need to predict how much a
refrigerator will cost to run (the energy bill) and how much it will cost to buy (the
capital equipment bill). Accordingly, our first job is to define heat and work, the
energy that we must buy. Work is fairly easily defined as the product of a force
andR a displacement. Its symbolic representation takes many different forms, such
as P dV, but it is always the product of an external force and the displacement that
goes with that force. Heat is more troublesome to calculate. The driving force for
heat, DT, seems obvious and analogous to the force or P term of the work equation.
But what is the equivalent displacement for heat? What gets moved, if there is such
an analogy?
Basic Principles 27
Carnot reasoned, by the flow of water, that there must be a flow of heat. He
also reasoned that in the case of heat, something must be conserved, as energy itself
is conserved (here he was right). Carnot discovered that the quantity that was con-
served, and that was displaced, was not heat but rather the quantity Q=T. This led
Carnot to the concept of entropy. In this book, entropy (S) is used to calculate heat
(Q). That calculation is analogous to the one for work. Just as dW ¼ P dV, likewise
dQ ¼ T dS. There are certainly more sophisticated implications of entropy, but for
us, entropy is just a way to calculate how much heat energy is involved.
We need these ways to calculate heat and work because it is heat and work
(energy) that we pay for to run the equipment. You will never see a work meter
or a heat meter on any piece of cryogenic equipment. Instead it is a fact of nature
that these two economically important quantities must be calculated from the ther-
modynamically important quantities, the things that we can measure: pressure (P),
temperature (T), and density (actually, specific volume). Hence, we need two more
tools.
One essential tool is an equation of state to tie pressure, volume, and tempera-
ture together.
The second essential tool to help us figure out the cost is the thermodynamic
network that links not only P, V, and T together but also links all of the other ther-
modynamic properties (H, U, S, G, A, . . . —you know what I mean) together in a
concise and consistent way. In real life, we usually measure pressure and temperature
because it is possible to buy pressure gauges and thermometers and it is not possible
to buy entropy or enthalpy gauges. Hence, we need a way to get from P and T to H
and S. (The beautiful thermodynamic network that gets us around from P to T to
heat and work is shown later in Fig. 2.7.) This thermodynamic network is the only
reason we delve into partial differential equations, the Maxwell relations, and the
like. The payoff is that we can go from easily measured quantities, pressure and tem-
perature, to the hard-to-measure, but ultimately desired, quantities of heat and work.
You might want to look at Fig. 2.7 now to discover a reason to plow through all the
arcane stuff between here and there. Once we do get our hands on the thermody-
namic network, we can calculate heat and work for any process, for any fluid, for
all time.
We will start our cryogenic analysis with some definitions. Admittedly, this is a
boring approach. However, it will ensure that we are all singing from the same sheet
of music. We begin with nothing less than the definition of thermodynamics itself.
2. THERMODYNAMICS
Thermodynamics is built upon four great postulates (guesses), which are called,
rather stuffily, the four laws of thermodynamics.
The second law was discovered first; the first law was discovered second; the
third law is called the zeroth law; and the fourth law is called the third law. I am
sorry to tell you that things do not get much better.
The zeroth law of thermodynamics just says that the idea of temperature makes
sense.
The first law of thermodynamics is simply the conservation of energy principle.
Inelegantly stated, the first law says that what goes into a system must either come
out or accumulate. Undergraduates often call this the ‘‘checkbook law’’. If you can
balance your checkbook, you can do thermodynamics.
28 Chapter 2
The second law is the entropy principle. The second law is also a conservation
equation, stating that in the ideal engine, entropy (S) is conserved. All real engines
make entropy.
The third law just says that there is a temperature so low that it can never be
reached.
These colloquial expressions of the four great principles of thermodynamics
have a calming effect. We will not be able to quantify our analyses unless we use
much more precise definitions and, expressly, the mathematical statements of these
principles. That is the task we must now turn to.
Thermodynamics is concerned with energy and its transformations into various
forms. The laws of thermodynamics are our concepts of the restrictions that nature
imposes on such transformations. These laws are primitive; they cannot be derived
from anything more basic. Unfortunately, the expressions of these laws use words
that are also primitive; that is, these words are in general use (e.g., energy) and have
no precise definition until we assign one. Accordingly, we begin our discussion of the
fundamental concepts of thermodynamics with a few definitions.
Thermodynamics: Thermodynamics is concerned with the interaction between a
system and its surroundings, the effect of that interaction on the properties of the
system, and the flow of heat and work between the system and its surroundings.
System: Any portion of the material universe set apart by arbitrarily chosen but
specific boundaries. It is essential that the system be clearly defined in any thermo-
dynamic analysis.
Surroundings: All parts of the material universe not included in the system.
The definitions of system and surroundings are coupled. The system is any
quantity of matter or region mentally set apart from the rest of the universe, which
then becomes the surroundings. The imaginary envelope that distinguishes the
system from the surroundings is called the boundary of the system.
Property of the system, or state variable: Any observable characteristic of the system,
such as temperature, pressure, specific volume, entropy, or any other distinguishing
characteristic.
State of the system: Any specific combination of all the properties of the system such
as temperature, pressure, and specific volume. Fixing any two of these three proper-
ties automatically fixes all the other properties of a homogeneous pure substance
and, therefore, determines the condition or state of that substance. In this respect,
thermodynamics is a lot like choosing the proper size dress shirt for a gentleman.
If you specify the neck size and sleeve length, the shirt size is defined. The size of
a man’s dress shirt is determined by two variables: neck size and sleeve length. A
single-phase homogeneous thermodynamic system is defined by fixing any two
thermodynamic variables. If you can pick out shirt size correctly, you can do
thermodynamics.
Energy: A very general term used to represent heat, work, or the capacity of a system
to do heat or work.
Heat: Energy in transition between a system and its surroundings that is caused to
flow by a temperature difference.
Work: Energy in transition between a system and its surroundings. Work is always
measured as the product of a force external to the system and a displacement of the
system.
If there is no system displacement, there is no work. If there is no energy in
transition, there is no work. Consider the following conundrum: It has been raining
steadily for the last hour at the rate of 1 in. of rain per hour. How much rain is on the
Basic Principles 29
ground? Answer: none. Rain is water in transition between heaven and earth. Once
the water hits the earth, it is no longer rain. It may be a puddle or a lake, but it is not
rain.
Likewise, heat and work are energy in transition between system and surround-
ings. Heat and work exist only while in transition. Once transferred, heat and work
become internal energy, or increased velocity, or increased kinetic energy, but they
are no longer heat and work. Heat and work cannot reside in a system and are
not properties of a system.
Because heat and work are in transition, they have direction as well as magni-
tude, and there is a sign convention to tell us which way they are moving.
The convention among chemists and chemical engineers is:
Thus, the expression Q W is the algebraic sum of all the energy flowing from
the surroundings into the system. Q W is the net gain of energy of the system.
Process: The method or path by which the properties of a system change from one
set of values in an initial state to another set of values in a final state. For example, a
process may take place at constant temperature, at constant volume, or by any other
specified method. A process for which Q ¼ 0 is called an adiabatic process, for
instance.
Cyclic process: A process in which the initial state and final state are identical. The
Carnot cycle, mentioned earlier, is such a cyclic process. We will now examine the
Carnot cycle in the light of the preceeding definitions.
In 1800, there was an extraordinary Frenchman who was a statesman and gov-
ernment minister. His work was so important that he was known as the ‘‘organizer
of victory of the French Revolution’’. Lazare Nicholas Marguerite Carnot (1753–
1823) was also an outstanding engineer and scientist. A century later, this glorious
family tradition was carried on by his grandson, Marie François Sadi Carnot
(1837–1894), another gentleman engineer, who served as the President of the
Republic of France from 1887 to 1894, when he was assassinated.
The Carnot family member we are interested in is neither of these, but a cousin
who lived between them. Nicholas Leonard Sadi Carnot (1796–1832), called ‘‘Sadi’’
by everyone, was very interested in steam engines. He wanted to build the most effi-
cient steam engine of all to assist in military campaigns. Sadi needed the steam
engine that required the least fuel for a given amount of work, thereby reducing
the logistical nightmare of fueling the war engines.
He began by considering the waterwheel. He reasoned that it should be possi-
ble to hook two waterwheels together, one driven by falling water in the conven-
tional way. The other, joined to the first by a common shaft, should be able to
pick up the water spent by the first and carry it up to the headgate of the first. In
this perfect, frictionless engine, one waterwheel should be able to drive the other for-
ever. We know, as Sadi knew, that this is not possible because of the losses (friction)
30 Chapter 2
in the system. However, the idea led Sadi to the conclusion that the most efficient
engine would be a reversible one, an engine that could run equally well in either
direction. The perfect waterwheel could be driven by falling water in one direction
or simply run in reverse to carry (pump) the water to the top again. A pair of such
reversible engines could thus run forever. Anything less than a reversible engine
would eventually grind to a halt. Hardly anyone paid any attention to Carnot and
his seminal concept, and indeed, Carnot’s work was further obscured by his early
death.
Carnot’s perfect abstraction was noticed, however, by the German physicist
Rudolph Clausius and the Glaswegian professor of natural philosophy William
Thomson. The principle of the conservation of energy was well established by this
time. Crudely stated: the energy that goes into a system must either come out or
accumulate. In the steady state, exactly the energy that goes into the engine must
come out, regardless of the efficiency of the engine. Clausius and Thomson both
saw in Carnot’s work that a perfect engine would conserve energy like any other
engine, but in addition, the quantity Q=T would also be conserved. In the perfect
(reversible) engine, Qin=Tin ¼ Qout=Tout. Not only is energy conserved (the first
law) but also the quantity Q=T is conserved in the ideal engine.
This astonishing fact leads us to the concept of entropy, S. Entropy was defined
then as now as S ¼ Qrev=T because of the remarkable discovery that the quantity
Q=T, entropy, is conserved in the ideal (reversible) engine. We can summarize the
work of Clausius, Thomson, and Carnot by saying that in an ideal engine, Sin ¼ Sout;
in any real engine, Sout Sin. We now have a concise description of the most efficient
engine possible, Sin ¼ Sout. All real engines, regardless of their intended product,
produce entropy: SoutSin for any real engine.
This is the reason that the concept of reversibility is so important, and it leads
to the following definition:
Reversible process: A process in which there are no unbalanced driving forces. The
system proceeds from the initial state to the final state only because the driving forces
in that direction exceed forces opposing the change by an infinitesimal amount. In
short, no forces are wasted.
Basic Principles 31
the equation
dS ¼ dQrev =T ð2:1Þ
where T is the absolute temperature of the system.
The first law of thermodynamics cannot be properly formulated without the
recognition of internal energy U as a property. Internal energy is a concept. Like-
wise, the second law of thermodynamics can have no complete description without
the recognition of entropy as a property. Entropy is also a concept, just as internal
energy is a concept, but is usually much harder to grasp.
Abstract quantities such as internal energy U and entropy S are essential to the
thermodynamic solution of problems. Accordingly, a necessary first step in thermo-
dynamic analysis is to translate the problem into the terminology of thermody-
namics. One of the real systems most often encountered in cryogenic engineering
is the one described by pressure, temperature, and specific volume. These systems
are called PVT systems and are composed of fluids, either liquids or gases or both.
For such systems, yet another postulate is required.
Postulate 4. The macroscopic properties of homogeneous PVT systems in
equilibrium states can be expressed as functions of pressure, temperature, and
composition only.
This postulate is the basis for all the thermodynamic equations that will follow
and has enormous utility. It is the basis for assuming the existence of an equation of
state relating P, V, and T. Note that it neglects the effects of electric, magnetic, and
gravitational fields and also neglects viscous shear and surface effects.
As mentioned, these postulates are primitive; they cannot be derived from
anything more basic.
32 Chapter 2
written as
mgDZ mDv2
DðPE + KEÞ ¼ þ ð2:3Þ
gc 2gc
where m is the mass of the system, DZ is the elevation of the system above
a reference plane, v is the velocity of the system, g is the local acceleration
of gravity, and gc is the gravitational constant.
Example 2.1. Water is flowing at a velocity of 100 ft=s through a pipe elevated 250 ft
above a reference plane. What are the kinetic energy (KE) and the potential energy
(PE) of 1 lbm of the water?
The work done on a body to increase its velocity leads to the kinetic energy
concept. From the first law, the work done on the body to accelerate it is trans-
formed into the kinetic energy of the body.
Work is defined as force times displacement, or
W ¼F l
and
dW ¼ F dl for a constant force F
The force may be obtained from Newton’s second law of motion,
ma
F¼
gc
Therefore,
ma
dW ¼ dl
gc
By definition, acceleration is
dv
a¼
dt
or
m dl
dW ¼ dv
gc dt
Basic Principles 33
This equation may now be integrated for a finite change in velocity from v1 to v2:
Z
m v2 m v2 2 v1 2
W¼ v dv ¼
gc v 1 gc 2 2
or
The quantity mv2=2gc was defined as the kinetic energy by Lord Kelvin in 1856.
Thus,
mv2
KE ¼ ð2:4Þ
2gc
Similarly, the work done on a body to change its elevation leads to the concept
of potential energy.
If a body of mass m is raised from an initial elevation z1 to a final elevation z2,
an upward force at least equal to the weight of the body must be exerted on it, and
this force must move through the distance z2 z1. Since the weight of the body is the
force of gravity on it, the minimum force required is given by Newton’s second law
of motion:
1 1
F¼ ma ¼ mg
gc gc
where g is the local acceleration of gravity and gc is as before. The minimum work
required to raise the body is the product of this force and the change in elevation:
g
W ¼ F ðz2 z1 Þ ¼ m ðz2 z1 Þ
gc
or
g g mzg
W ¼ mz2 mz1 ¼ D
gc gc gc
34 Chapter 2
The work done on a body in elevating it is said to produce a change in its poten-
tial energy (PE), or
mzg
W ¼ DPE ¼ D
gc
Thus, potential energy is defined as
g
PE ¼ mz ð2:5Þ
gc
This term was first proposed in 1853 by the Scottish engineer William Rankine
(1820–1872).
For the problem statement above,
ft 1 s2 lbf
Potential energy ¼ ð1 lbm Þð250 ftÞ 32:174 2
s 32:174 ft lbm
¼ 250 ft lbf
Potential energy and kinetic energy are obvious mechanical terms. The energy
put into a system may be stored in elevation change or in velocity change as well as in
internal energy U by way of property changes.
The internal energy U represents energy stored in the system by changes in its
properties T, P, V, etc. DU is a composite effect of changes in the motion and
configuration of all of the ultimate particles of a system.
It is essential to note that the left-hand side of Eq. (2.2) describes a system–
surroundings interaction, while the right-hand side describes properties of the
system alone. This is a very valuable aspect of the first law equation, for it permits
the calculation of system–surroundings interactions (heat and work) from changes in
system properties only. Some of the consequences of this attribute of the first law
equation are discussed below.
As stated earlier, work is a system–surroundings interaction—energy in transit
between the system and surroundings, a precise kind of energy that crosses the
system boundary. Therefore, work must always have two components: an external
(surroundings) force and an internal (system) displacement. Work is always
Fexternal Xsystem
where X is the system displacement.
For fluids, work is always
Pexternal Vsystem
Only in the reversible case is Pexternal ¼ Pinternal. Then, and only then, may the P of
the system be substituted for Pexternal and equations of state be used to relate P
and V.
Because the left-hand side of Eq. (2.2) involves a system–surroundings interac-
tion, questions of path come into play (e.g., adiabatic, reversible, constant external
pressure, etc.). Indeed, the path must always be specified for calculations involving
the left-hand side of Eq. (2.2), i.e., heat and work. In thermodynamics, it is essential
to distinguish between quantities that depend on path and those that do not. Indeed,
we commonly make this distinction in ordinary affairs. For instance, the straight-line
distance between Denver and Los Angeles is a fixed quantity, a property. On the
other hand, other quantities for an auto-mobile trip between Denver and Los
Basic Principles 35
Angeles—miles traveled, miles per gallon consumed, etc.—depend very much on the
path.
As stated, the right-hand side of Eq. (2.2) describes attributes of the system
only and contains properties and potentials of the system only. Thus, it can always
be evaluated from system properties or potentials alone without considering the
surroundings. The question of path is simply irrelevant, because there is no
system–surroundings interaction to be found on the right-hand side of the first law
statement.
It is a constant goal of thermodynamics to evaluate system–surroundings inter-
actions (e.g., heat and work) in terms of properties and potentials of the system
alone. The essence of the first law concept is (a) that it is a conservation equation
and (b) that a system–surroundings interaction can be expressed as a function of
system parameters only.
The first law is valid for all systems under all cases, as written in Eq. (2.2).
There are numerous special cases, and some of the more important will now be
considered.
36 Chapter 2
Dv2 g
Q Ws ¼ H 2 H 1 þ þ DZ
2gc gc
or
Dv2 g
Q Ws ¼ DH þ þ DZ
2gc gc
If velocity and elevation changes are negligible, then Q Ws ¼ DH. This is the
flow system analog of the first law for closed systems.
Consideration of other special cases of the first law leads to other auxiliary
functions, discussed below.
Basic Principles 37
QP ¼ P0 ðV2 V1 Þ þ DU
¼ P0 V2 P0 V1 þ U2 U1
¼ ðP0 V2 þ U2 Þ ðP0 V1 þ U1 Þ
and since P0 ¼ P2 ¼ P1 ¼ constant, this may be rewritten as
QP ¼ ðP2 V2 þ U2 Þ ðP1 V1 ¼ U1 Þ
which, from the definition of enthalpy, H U þ PV
QP ¼ H 2 H 1 ¼ DH
Thus, at constant pressure,
Q ¼ ðDHÞP
Differentiating this expression with respect to temperature suggests the function
dQ dH
¼ CP
dT P dT P
which is the definition of heat capacity at constant pressure.
38 Chapter 2
Q has a quality as well as a quantity, and that quality depends on the temperature.
Also, the effect of adding Q to a system is always to change the properties of
the system, so there must be a property change that describes the addition of Q.
We postulate a property in thermodynamics to describe these observations that
meets the following qualitative requirements:
1. The addition of Q causes a property change in the system.
2. The way (path) in which Q is added fixes the magnitude of the property
change.
3. Adding Q at a lower temperature causes a greater change than adding the
same amount of Q at a higher temperature.
Qrev ¼ TDS
Since Q is not a property, its path must be specified, and Qrev is chosen to complete
the definition.
Thus, we have developed qualitatively Postulate 3, which was stated earlier:
there exists a property called entropy S, which is an intrinsic property of a system
that is functionally related to the measurable coordinates that characterize the
system. For a reversible process, changes in this property are given by
dQrev
dS ¼
T
The definition of entropy, Qrev ¼ TDS, can always be written in differential
form because neither reversibility nor any other path has any significance in differ-
ential notation:
dQ ¼ T dS
U ¼QW
becomes in differential form, for all cases and for all substances,
dU ¼ T dS P dV ð2:7Þ
This is a universal relationship, just like the first law. One can always write the first
law in this differential form for all cases. Although derived for a reversible process,
Eq. (2.7) relates properties only and is valid for any change between equilibrium
states in a closed system. The question of reversibility=irreversibility has no signifi-
cance to Eq. (2.7) because (a) the concept of path has no bearing on any equation
written in differential form and (b) the equation dU ¼ T dS P dV contains proper-
ties of the system only. Path is a concept relevant to system–surroundings interac-
tions, and ‘‘surroundings’’ are completely absent from this description of the
system alone.
Basic Principles 39
or
Z
T dS ¼ Q þ lost work
The lost work is energy that might have done useful work but is dissipated instead,
causing unnecessary randomness and disorder within the system.
The second law can be stated formally as follows.
Second law of thermodynamics. The entropy change of any system and its surround-
ings, considered together, is positive and approaches zero for any process that
approaches reversibility.
The second law of thermodynamics is a conservation law only for reversible
processes, which are unknown in nature. All natural processes result in an increase
in total entropy. The mathematical expression of the second law is simply
DStotal 0
where the label ‘‘total’’ indicates that both the system and its surroundings are
included. The equality applies only to the limiting case of a reversible process.
Since the energy change of the system for the complete cycle must be zero,
DUcycle ¼ 0 ¼ Q1 Q2 W
40 Chapter 2
or
Q1 Q2 ¼ W
where W is the net work, which is the difference between the work produced in the
expansion process (b–c) and the work required in the compression process (d–a).
Also, the absolute values of Q1 and Q2 are used, temporarily abandoning the usual
sign convention.
The efficiency of this reversible cycle is
W Q1 Q2
Z¼ ¼
Q1 Q1
since Q1 ¼ T1DS and Q2 ¼ T2DS
The efficiency of the Carnot engine can be expressed in terms of the tempera-
ture of the heat reservoirs alone:
T1 T2
Z¼
T1
Carnot stated (and proved) the following principles, which are known today as
the Carnot principles:
1. No engine can be more efficient than a reversible engine operating between
the same high-temperature and low-temperature reservoirs. Here the term
‘‘heat reservoir’’ is taken to mean either a heat source or a heat receiver or
sink.
2. The efficiencies of all reversible engines operating between the same
constant-temperature reservoirs are the same.
3. The efficiency of a reversible engine depends only on the temperatures of
the heat source and heat receiver.
The two equations for the efficiency of the Carnot cycle can be combined to
yield
W Q1 Q2 T1 T2 T2
Z¼ ¼ ¼ ¼ 1
Q1 Q1 T1 T1
or
T2
W ¼ Q1 1
T1
Basic Principles 41
But
W ¼ Q1 Q2
Therefore,
T2
Q1 Q2 ¼ Q1 1
T1
which reduces to
Q1 Q2
¼
T1 T2
or S1 ¼ S2, which is to repeat that entropy is conserved in a reversible engine.
In 1848, Lord Kelvin (William Thomson of Glasgow) derived a temperature
scale based on the above equation, independent of the thermodynamic substance
employed in the cycle. He defined the scale such that
T 1 Q1
¼
T 2 Q2
42 Chapter 2
Basic Principles 43
as exact differentials. Thus, changes in these functions depend only on the initial and
final states of the system and not on the method or path by which the change takes
place. These functions are
U ¼QW Internal energy
H ¼ U þ PV Enthalpy
G ¼ H TS Gibbs free energy
A ¼ U TS Helmholtz free energy
The four state functions written in differential form are
dE ¼ T dS P dV ðEq: 2:8Þ
dH ¼ T dS þ V dP ðEq: 2:9Þ
dG ¼ S dT þ V dP ðEq: 2:11Þ
dA ¼ S dT P dV ðEq: 2:13Þ
It would be imprecise to stop with the impression that these four state func-
tions are derived from considering the maximum possible work connected with a
reversible process. Such a consideration only suggests the function. New thermody-
namic functions such as H, G, and A cannot in general be defined by the random
combination of variables. As a minimum, dimensional consistency is required.
Instead, a standard mathematical method exists for the systematic definition of
the desired functions; this is the Legendre transformation. The functions for G, H,
and A are in fact Legendre transformations of the fundamental property relation
of a closed PVT system, dU ¼ TdS PdV.
A consequence of the Legendre transformation is that the resulting expressions
for G, H, and A are exact differential equations with the following special property:
If z ¼ f(x,y) and dz ¼ Mdx þ Ndy, then the property of exactness is that
@M @N
¼
@y x @x y
Applying this property of exact differential equations to the four equations above
yields the following Maxwell relations:
@T @P
¼ ð2:14Þ
@V S @S V
@T @V
¼ ð2:15Þ
@P S @S P
@S @V
¼ ð2:16Þ
@P T @T P
@S @P
¼ ð2:17Þ
@V T @T V
The last two are among the most useful, because they relate entropy changes to an
equation of state.
The development of many of these relations is shown in Fig. 2.7, as proposed
by Dr. A. J. Kidnay of the Colorado School of Mines.
Basic Principles 45
The molecules are in random motion. Their velocity, a vector, can be resolved into
three components. Take one component normal to the wall. Since the motion is ran-
dom, on the average, one-third of the velocities are normal to the wall and half of
those, one-sixth of the total, are moving toward the wall at any time (the other half
are moving away from the wall).
The total number of molecules in the small box shown in the figure is the number
per unit volume, n=V t, times the volume A dL. One-sixth of that number is moving
toward the wall and will strike the wall in the time dt that it takes to travel the distance
dL at the average velocity, u ¼ dL=dt. Thus, the rate at which molecules strike the wall is
1 n 1 1 n u 1 n
A dL ¼ A dL ¼ Au
6 Vt dt 6 Vt dL 6 Vt
46 Chapter 2
When molecules strike the wall, their velocity is reversed. So the total change in velocity
is 2u. The rate of change of momentum is the rate of striking the wall times the mass per
molecule (i.e., molecular weight times the change in velocity, 2u). Thus,
1 n 1 Mn
F¼ AuðMÞð2uÞ ¼ Au2
6 Vt 3 Vt
1 Mu2 3
Mu2 ¼ RT or ¼ RT
3 2 2
where Mu2=2 is the average kinetic energy per mole. Thus, we have shown that for an
ideal gas the kinetic energy is directly proportional to the absolute temperature. Or, put
another way, temperature is a direct measure of kinetic energy. Notice also that
pressure is equal to kinetic energy per unit volume.
This short derivation also provides some insight into the universal gas
constant R. R is the chemist’s version of Boltzmann’s constant k. R is in units per
mole, and k is in units per molecule. Both R and k are scaling factors between tem-
perature and energy. For physicists, E ¼ kT. For chemists, it is PV ¼ RT, where the
product PV is energy. Several values of R are given in Table 2.1.
We have already defined the heat capacities Cv and Cp. The ideal gas derivation
can also breathe some life into these two important quantities. Heat capacity
describes the ability of some substance to ‘‘soak up’’ energy.
Consider the molecules of a monatomic ideal gas. The only way a monatomic
gas can soak up energy is to increase in kinetic energy. That is, the only way a mona-
tomic ideal gas can ‘‘contain’’ energy (short of energies large enough to affect elec-
tronic or nuclear changes) is by the motion of its molecules—translational kinetic
energy. Thus, if energy is added to a box of such molecules at constant volume,
the change in their internal energy will be
dH ¼ CP dT ¼ dU þ P dV ¼ CV dT þ R dT ¼ ðCV þ RÞ dT
Thus, CP ¼ CV þ R. This is true for any ideal gas (monatomic or not, because the
work expression is always the same).
For polyatomic molecules, the internal energy includes the translational kinetic
energy, but the molecules also can ‘‘soak up’’ energy in the form of rotation and
vibration about their bonds. Polyatomic molecules thus have greater CV (and CP)
Basic Principles 47
Table 2.1 Values of the Gas Constant R for 1 mol of Ideal Gas
Temp. unit Pressure unit Volume unit Energy unit R per gram mol
K Calorie 1.9872
K atm cm3 82.057
K atm L 0.082054
K mmHg L 62.361
K bar L 0.08314
K kg=cm2 L 0.08478
R per lb mol
R Btu (IT) 1.986
R hp hr 0.0007805
R kW hr 0.0005819
R atm cu ft 0.7302
R in. Hg cu ft 21.85
R mmHg cu ft 555.0
R psia cu ft 10.73
psfa cu ft ft lb 1545.0
K atm cu ft 1.314
K mmHg cu ft 998.9
48 Chapter 2
PV
¼ 1 þ Br þ Cr2 þ ð2:20Þ
RT
2. As an expansion in pressure:
PV
¼ 1 þ B 0 P þ C 0 P2 þ ð2:21Þ
RT
The virial coefficients (B, C, . . . ) can be calculated from first principles (two-body
interactions, three-body interactions, etc.) and are functions of temperature only.
Generalized or corresponding state forms of the equation of state are also useful:
PV
¼ Z; Z ¼ ZðTr ; Pr Þ ð2:22Þ
RT
where Z is the compressibility factor and Tr and Pr are the reduced temperature and
pressure, e.g., Tr ¼ T=Tc. T is the system temperature, and Tc is the critical tempera-
ture of the system component.
Another useful generalized equation of state is the three-parameter corre-
sponding states form employing the acentric factor o,
Z ¼ f ðTr ; Pr ; oÞ ð2:23Þ
S ¼ f ðP; TÞ
@S @S ð2:24Þ
dS ¼ dT þ dP
@T P @P T
Step 3. Evaluate the partial differential coefficients from the basic Maxwell relations.
For example, evaluate the first coefficient, (@S=@T)P. Because this coefficient
involves an entropy term (S) that changes with T at constant P, it is likely that H
will be involved.
From the definition of dH,
dH ¼ T dS þ V dP
Basic Principles 49
we have
@H @S
¼T þ0
@T P @T P
Since at constant pressure dP ¼ 0, and since
@H
CP
@T P
@S
CP ¼ T
@T P
therefore,
@S 1
¼ CP ð2:25Þ
@T P T
So far, there are no restrictions on this equation.
50 Chapter 2
Step 4. Now evaluate the second coefficient, (@S=@P)T. From the relationship
dG ¼ S dT þ V dP
and the property of exactness described before, we have
@S @V
¼ ð2:26Þ
@P T @T P
Step 5. The total differential, dS, is now complete as
CP @V
dS ¼ dT dP ð2:27Þ
T @T P
and
dT dP
dS ¼ CP R
T P
We now have the total derivative of S as a function of properties only for any
process. The integration to yield DS is straightforward. This equation is now restricted
to the ideal gas state, since that was the equation of state chosen to evaluate @V=@TP.
This is a perfectly general method for any thermodynamic function: DS, DH,
DU, DG, etc. Since such equations contain only properties of the system, they are
independent of path and can be tabulated for various fluids. Such tabulations form
the basis for the various thermodynamic charts: T–S, H–S, etc.
Example 2.2. Calculate internal energy (U) and entropy (S) as a function of T
and V.
It is sometimes more convenient to take T and V as independent variables
rather than T and P. The method described above is perfectly general and can be
applied here. Writing U and S as total differentials:
@U @U
dU ¼ dT þ dV ð2:28Þ
@T V @V T
@S @S
dS ¼ dT þ dV ð2:29Þ
@T V @V T
Recall that
@U
CV ¼ ð2:30Þ
@T V
Basic Principles 51
and
@U @S
¼T P ð2:32Þ
@V T @V T
52 Chapter 2
Process W Q DU DH
Isobaric (constant P)
Real P(V2 V1) nCP(T2 T1) QW nCP(T2 T1)
Ideal gas nR(T2 T1) nCV(T2 T1)
Isometric (constant V)
Real 0 nCV(T2 T1) Q DU þ D(PV)
Ideal gas nCV(T2 T1) nCP(T2 T1)
Isothermal (constant T) Q DU DU þ W QW DU þ D(PV)
Real nRT ln(V2=V1) nRT ln(P1=P2) 0 0
Ideal gas nRT ln(P1=P2) nRT ln(V2=V1)
Isentropic (constant S)
Real (nonflow) DU 0 W DU þ D(PV)
Ideal gas (nonflow) (P2V2 P1V1)=(1 g) nCV(T2 T1) nCP(T2 T1)
nCV(T2 T1)
Real (flow) DH 0 DH D(PV) W
Ideal (flow) g(P2V2 P1V1)=(1 g)
Polytropic
Real (nonflow) Q DU DU þ W QW DU þ D(PV)
(P2V2 P1V1)=(1 n)
Ideal (nonflow) Q nCV(T2 T1) nCV(T2 T1) nCP(T2 T1)
Real (flow) Q DH DH þ W DH D(PV) QW
Ideal (flow) n(P2V2 P1V1)=(1 n)
Cyclic Qnet Wnet 0 0
The choice of these parameters will generally be made on the basis of experience or
convenience.
Application of the first law to the ideal work of isothermal compression in a
refrigerator cycle with appropriate substitution for the heat summation reduces to
Ws =m
_ ¼ T1 ðS1 S2 Þ ðH1 H2 Þ ð2:38Þ
where ṁ is the mass rate of flow through the compressor, S is the specific entropy, H
is the specific enthalpy, and the subscripts 1 and 2 refer to the inlet and exit streams,
respectively, of the isothermal compressor. In a real system, an overall compressor
efficiency factor will have to be incorporated in this equation to obtain meaningful
results.
Evaluation of the refrigeration duty in a refrigerator cycle is dependent on the
process chosen. For example, if the refrigeration process is isothermal, the refrigera-
tion duty can be evaluated using the relation
Q¼m
_ TðS2 S1 Þ ð2:39Þ
where T is the temperature at which the refrigeration takes place and the subscripts 1
and 2 refer to the inlet and exit streams of the process. On the other hand, if the
refrigeration process involves only the vaporization of the saturated liquid refriger-
ant at constant pressure, evaluation of the duty reduces to
Q¼m
_ hfg ð2:40Þ
where hfg is the specific heat of vaporization of the refrigerant. If the process is
performed at constant pressure with the absorption of only sensible heat by the
Basic Principles 53
refrigerant (e.g., warming a cold gas), the duty can be evaluated from
_ ðH2 H1 Þ
Q¼m ð2:41Þ
Similar first law balances could be written for other units of equipment used at low
temperatures.
54 Chapter 2
Example 2.3. Show that the Joule–Thomson coefficient for a perfect gas is always
zero.
The Joule–Thomson coefficient has been defined as
@T
m¼ ð2:42Þ
@P H
The expression (@H=@T)P has previously been identified as the heat capacity, Cp. The
expression (@H=@P)T can be found by differentiating Eq. (2.9 ) and combining the
result with Eq. (2.16 ) to yield
@H @V
¼ V T
@T P @T P
becomes
@V
dH ¼ CP dT þ V T dP ð2:44Þ
@T P
Setting dH ¼ 0 (constant enthalpy), Eq. (2.45 ) gives the Joule–Thomson coeffi-
cient in terms of measurable thermodynamic properties as
1 @V
mJT ¼ T V ð2:45Þ
CP @T P
For a perfect gas, V ¼ RT=P, and
@V R V
¼ ¼
@T P P T
Basic Principles 55
3. HEAT TRANSFER
All cryogenic systems operate, that is to say complete their required tasks, at very
low temperatures. This fact requires a reasonably strong understanding of heat
transfer and all associated technologies. As stated earlier, the best that one can hope
for in some cases (e.g., cryogenic fluid storage) is to severely limit the energy transfer
across boundaries. Since heat always flows from hot to cold and cannot be stopped,
how it is controlled is the design aspect relating to heat transfer. In other cases (like
gasification of cryogenic fluids) the ability to maximize the amount of heat intro-
duced into the fluid is the desired design state. Due to the nature of the temperature
difference between the everyday world we live in and the extremes of operating ther-
mal conditions of cryogenic systems, a number of unique problems present them-
selves. These include, but not limited to: (1) variability of material properties, (2)
thermal insulation, (3) fluid property differences for near-critical-point convection,
(4) effects of thermal radiation, (5) equipment design characteristics and operation,
(6) selection of instrumentation and control sensors, (7) transfer of cryogenic fluids,
and (8) disposal and=or control of cryogenic boil-off. Each of these problems will be
addressed in the remaining chapters of this book. Heat transfer at low temperatures
is governed by the same three mechanisms present at ambient and elevated tempera-
tures: conduction, convection, and radiation. It is not surprising, therefore, that all
the general equations are appropriate for low-temperature applications as long as
they are adjusted for the property changes in both materials and fluids.
56 Chapter 2
determined. The values of KH and KC are also useful when their difference is divided
by the difference of TH and TC and becomes the apparent thermal conductivity of
the material, expressed as
KM ¼ ðKH KC Þ=ðTH TC Þ
Figures 2.9 and 2.10 are graphical depictions of thermal conductivity integrals
for some metals and insulators, respectively.
The conduction shape factor S is derived from the Fourier equation and simply
places all the physical parameters of the shape of the material into a single integral.
This allows for various geometric shapes to be evaluated and conduction shape fac-
tors to be determined based on the desired physical dimension of interest, such as
wall thickness or sphere diameters. Table 2.5 gives some representative conduction
shape factors for various geometries of interest. These geometries are pictured in
Fig. 2.11.
With this simple expression it is possible to get a good engineering estimate
on the heat balance around various cryogenic components under steady-state
conditions
dT
Q ¼ kA ð2:47Þ
dx
where Q is the heat transferred, k the material thermal conductivity (normally a
function of temperature), A the area of the material through which the heat is trans-
ferred and dT=dx the change in temperature observed in the material in the x direc-
tion due to the heat transfer.
If Q, k, and A are constant, the relation becomes
T2 T1
Q ¼ kA ð2:48Þ
Dx
This relation is useful when determining the heat transfer through a layer of insula-
tion or a series of insulation materials. If the conduction is through insulation
wrapped around a pipe, the relationship is modified to
dT T1 T2
Q ¼ kA ¼ 2 p Lk ð2:49Þ
dr lnðr2 =r1 Þ
where L is the length of the pipe, r1, r2 the inside and outside radii of the insulation
and T1, T2 the corresponding temperatures at these locations.
For nonsteady-state conductive heat transfer, a heat balance may be made
around a volume element of material of dimensions dx, dy, and dz. In the x direc-
tion, the heat transferred into the volume element is given by
@T
dQx ¼ k dy dz ð2:50Þ
@x
The heat leaving the volume element is
@dQx
dQ0x ¼ dQx þ dx ð2:51Þ
@x
If k is constant, then
Basic Principles 57
Figure 2.9 Thermal conductivity integrals of metals. Data provided from S. Courts,
Lakeshore Cryotronics.
@2T
dQ0x dQx ¼ k dy dz dx ð2:52Þ
@x2
Similar expressions are valid for the other two directions. If heat is stored in the
volume element, then
58 Chapter 2
Figure 2.10 Thermal conductivity integrals of insulating materials. Data provided from
S. Courts, Lakeshore Cryotronics.
Basic Principles 59
so
2
@T @2T @2T @ T
¼a þ þ ð2:55Þ
@y @x2 @y2 @z2
where y is the time, CP and r are the heat capacity and density of the material,
60 Chapter 2
Q ¼ hc AðTs Tb Þ ð2:56Þ
where hc is the individual heat transfer coefficient, Ts is the surface temperature, and
Tb is the bulk fluid temperature. The individual heat transfer coefficient is critical in
the solution of Eq. (2.56) and must generally be obtained with the aid of some
empirical relationship that has been found to correlate the convective heat transfer
reasonably well for a specific geometry or flow pattern. The empirical relationships
used to obtain the individual heat transfer coefficients at low temperatures are very
similar to the relationships used at ambient temperatures.
Tables 2.6–2.8 provide some useful heat transfer relationships for
low-temperature heat exchangers for a variety of operating conditions. The units
for the relationships in these tables are in consistent American engineering units.
Conversion factors to other systems of measurement are found in the Appendix.
If heat is transferred from one fluid to another separated by a pipe as in a heat
exchanger, then Eq. (2.56 ) can be modified to include the individual heat transfer
coefficients for both fluids in an overall heat transfer coefficient as
Q ¼ Uo Ao DToverall ð2:57Þ
where Uo is based on the outside surface area Ao. If the inside surface area is used,
the relationship can be modified to
Q ¼ Ui Ai DToverall ð2:58Þ
1
Ui ¼ ð2:59Þ
1=hi þ Drp Ai =kp Ap;lm þ Ai =ho Ao
where hi is the inside heat transfer coefficient, Drp is the thickness of the pipe, Ai is
the inside surface area, kp is the thermal conductivity of the pipe, Ap,lm is the log
mean area of the pipe given by (Ao Ai)[ln (Ao=Ai)], Ao is the outside surface area,
and ho is the outside heat transfer coefficient. A similar expression can be written for
Uo. For a flat wall, Eq. (2.59 ) reduces to
1
Ui ¼ Uo ¼ ð2:60Þ
1=hi þ Dx=k þ 1=ho
Empirical pressure
Type Flow conditions Empirical heat transfer equation drop=length equation
" #
0:0668ðDe =LÞRe Pr k DP 32 G 2
Straight tubular pipe Flow inside, Re < 2100; h ¼ 3:658 þ 2=3 D
¼
1 þ 0:04ðDe =LÞRe Pr e DL Re gc De r
no phase change (gas)
" 2=3 # 0:14
De CP mf DP 0:158 G2
DP 0:952 G2
Collins tubing Flow in annulus, h ¼ 0.118 CPG Re0.2 Pr2=3 ¼ 0:2 max
DL Re gc De r
400 < Re < 10,000;
no phase change (gas)
For tubular pipe: De ¼ inside pipe diameter; for noncircular tubes, De ¼ 4(AcL=Aw), where Ac ¼ inside free-flow cross-sectional area, L ¼ length of pipe, and Aw ¼ heat transfer
or wetted tube surface area; for annulus, De ¼ D2 D1, where D2 ¼ inside diameter of outer pipe and D1 ¼ outside diameter of inner pipe; Dh ¼ diameter of helix.
For Collins tubing: De ¼ 4(Va=Aw), where Va ¼ active flow volume in annulus, Aw ¼ wetted area in annulus, including fins.
Evaluate fluid properties for all relations at the mean film temperature, Tm ¼ 0.5(Tw þ Tb), where Tw ¼ tube wall temperature and Tb ¼ bulk temperature of fluid. Constants in
these equations are for use with 51 units.
61
© 2005 by Marcel Dekker
62
Table 2.7 Empirical Heat Transfer and Corresponding Pressure Drop Relationships for Flows Normal to Outside of Tubes on Shell Side of Coiled Tube
Heat Exchangers
Re ¼ DoGmax=m, where Do ¼ outside diameter of tubes; Gmax ¼ m_ =Amin, with Amin ¼ minimum open free-flow area between tubes; XT ¼ transverse pitch=tube outside diameter;
N ¼ total number of tubes in line of flow; XL ¼ longitudinal pitch=tube outside diameter; n ¼ 0.43 þ (1.13=XL). Fluid properties for all relations are evaluated at the mean film
temperature, Tm ¼ 0.5(Tw þ Tb), where Tw ¼ tube wall temperature and Tb ¼ bulk temperature of fluid.
Chapter 2
© 2005 by Marcel Dekker
5367-4 Flynn Ch02 R2 090804
Basic Principles 63
Table 2.8 Empirical Heat Transfer and Corresponding Pressure Drop Correlations for
Brazed Aluminum Plate-Fin Heat Exchangersa
64 Chapter 2
Figure 2.12 Pooling boiling depictions: (a) saturated conditions generate vapor and raises
the bulk fluid temperature while (b) subcooled conditions shown allows for no vapor genera-
tion but does raise in the bulk fluid temperature, as well.
boiling. There is a further subdivision of two primary boiling types (nucleate and film)
based upon whether the boiling fluid is under saturated or sub-cooled conditions. The
nucleate boiling conditions are depicted in Fig. 2.12 and 2.13, while 2.14 illustrates
the filming boiling process. Most experimental work at low temperatures has been
conducted on pool boiling. Such studies have shown that the boiling action at these
temperatures is also largely dependent on the temperature difference between the sub-
merged surface and the boiling liquid. Unfortunately, the available information on
boiling liquids at low temperatures is limited and also shows wide scatter in the
experimental data. Figure 2.15 shows typical experimental nucleate and film pool
boiling data for nitrogen at 1atm compared with predictive correlations. Figures
2.16–2.18 give similar data for oxygen, hydrogen, and helium. See ‘‘Suggested Read-
ing’’ at the end of the book for the original data sources called out on Figs. 2.15–2.18.
Figure 2.13 Forced convection boiling: (a) subcooled conditions again allow for no vapor
generation while (b) saturated or bulk boiling generates vapor and any number of different
vapor-liquid flow conditions.
Basic Principles 65
Figure 2.14 Film boiling occurs when a vapor covers the super-heated bottom surface: (a) if
conditions are subcooled, then no vapor will be generated, and (b) at saturated conditions
large amounts of vapor can be formed.
For some materials, for example glass or plastic, this can be written as
aþtþr¼1 ð2:63Þ
Figure 2.15 Experimental and predictive nucleate and film boiling of nitrogen at 1 atm.
66 Chapter 2
Figure 2.16 Experimental and predictive nucleate and film boiling of oxygen at 1 atm.
where Eb is the emissive power of a blackbody at the same temperature. The emissive
power for a blackbody is given by
Eb ¼ sT 4 ð2:65Þ
Figure 2.17 Experimental and predictive nucleate and film boiling of hydrogen at 1 atm.
Basic Principles 67
Figure 2.18 Experimental and predictive nucleate and film boiling of helium at 1 atm.
E ¼ EsT 4 ð2:66Þ
When two bodies within visual range of each other are separated by a medium that
does not absorb radiation, an energy exchange occurs as a result of reciprocal pro-
cesses of emission and absorption. The net rate of heat transfer from one surface at
T1 to another surface at T2 can be calculated from
where A1 is the surface area of one of the bodies, FA is a shape and orientation factor
for the two bodies relative to area A1, and FE is the emission and absorptance factor
for the two bodies. If the surface of one body is small or enclosed by the surface of
the other body, FA ¼ 1; otherwise, FA must be evaluated. For a surface A1 radiating
to a large area or space, FE ¼ E1. For two large parallel plates with a small space
between them,
1
FE ¼ ð2:68Þ
1=E1 þ 1=E2 1
1
FE ¼ ð2:69Þ
1=E1 þ ðr1 =r2 Þ2 ½ð1=E2 Þ 1
68 Chapter 2
For the case where FA ¼ 1 and FE ¼ E1, we can define a radiative heat transfer coeffi-
cient as
s T14 T24
hr ¼ ð2:71Þ
T1 T1 T2
This can then be used to determine the heat transferred by radiation from the
relation
Q ¼ hr A1 ðT1 T2 Þ ð2:72Þ
4. MOMENTUM TRANSFER
Cryogenic fluids are often handled near their boiling point. Thus, flow is usually compli-
cated because vaporization occurs, resulting in two-phase flow. The transition from liquid
to vapor depends on the heat influx and pressure changes encountered in the system.
For single-phase flow, the most important point to recognize is the distinction
between laminar and turbulent flow. For low temperatures, the same concepts that
have been developed for ambient flow apply. For flows with Reynolds numbers less
than 2100, laminar flow is assumed to occur; and for flows with Reynolds numbers
greater than 3000, turbulent flow is assumed to occur. The range in between 2100
and 3000 is normally considered the transition region between laminar and turbulent
flow. Thus, the empirical relationships that have been used successfully in calculat-
ing, for example, the pressure drop in a pipe at ambient conditions are also valid
for systems using cryogenic fluids.
Tables 2.5–2.7 provide some useful relationships to determine the empirical heat
transfer coefficient and pressure drop per unit length of many common heat exchangers.
In cryogenic systems, two-phase flow cannot be avoided, particularly during
the cooldown of transfer lines. If a cooldown situation is encountered frequently,
complete system performance reliability requires design for this condition. It takes
only a slight localized heat influx or pressure drop to cause the evaporation of part
of the fluid. Vapor in such lines can significantly reduce the liquid transfer capacity.
This transition from liquid to vapor is usually an unequilibrated condition and
presents an undefined problem for design purposes.
The Martinelli–Nelson correlation has been found to be useful for estimating pres-
sure drop for a single-phase stream under steady-state flow conditions. However, this cor-
relation, developed for ambient conditions, has a tendency to underestimate the pressure
drop by as much as 10–30% because of the high tendency of cryogenic fluids to vaporize.
Another fluid problem encountered rather frequently is that of cavitation in
pumps. Cavitation begins at a pressure slightly below where vaporization just begins.
The difference in pressure between cavitation and normal vaporization is equal to the
net positive suction head minus the drop in head pressure in the pump. Normally, this
must be determined experimentally and is a characteristic of the pump and the fluid.
A critical flow can be achieved in the transfer of a cryogenic fluid. In such a
situation, the flow rate of the fluid cannot be altered by lowering the downstream
Basic Principles 69
pressure as long as the supply pressure remains constant. Such a situation is due to
the increased down-stream vaporization, which reduces the cross-sectional area
available for liquid transfer.
Cooldown of cryogenic systems involves a period of surging flow and pressure
fluctuations with accompanying sound effects. The primary factors that affect surging
and maximum pressure include inlet pressure and temperature, pipe length and
diameter, pipe precooling, throttling, and liquid properties. Other factors such as
heat capacity, external heat leak, and pipe wall thickness have a minor effect on
surging, because the pipe temperature undergoes little change during such a surge.
Cryogenic fluid transfer presents a number of challenges to any designer but two
problem areas, in particular, have wide ranging impacts on the overall system design:
(1) heat leak in cryogenic transfer, and (2) pressure drop during cryogenic transfer of
fluids.
The transfer piping design should accommodate the desired heat transfer
requirements as well as the fluid transfer conditions necessary to accommodate the
overall system design. A comparison of the heat leak into a bare pipe vs. those with
various insulations techniques is given in Table 2.9 for liquid nitrogen. The effect of
improved insulation techniques, such as vacuum jacketed pipe (VJP), is obvious and
can result in reduced heat loads of several orders of magnitude over bare pipes. VJP
comes in rigid and flexible configurations to accommodate the designer’s need for
making bends, etc. The flexible configuration requires a more complex and slightly
less efficient insulation scheme and therefore do not offer the same performance as
the rigid versions but is still very good vs. the nonvacuum jacketed alternatives. Both
rigid and flexible transfer lines come in various sizes and a number of commercial
vendors are available to supply these types of units.
Table 2.10 gives estimated heat leak for liquid nitrogen, oxygen, and hydrogen
with various line sizes. An additional consideration of VJP is the required cooldown
heat removal necessary to allow the line to efficiently transfer cryogenic fluids. Note
that flow is not a parameter. This is the case because fully developed turbulent flow is
assumed. For laminar flow to exist, the pipe diameter would have to be extremely
large, resulting in excessive costs and unnecessarily high heat leak into the system.
At turbulent flow conditions, the controlling insulating parameter to heat leak into
the fluid is the vacuum jacket around the transfer pipe. Solid thermal conductivity
through the pipe walls, as well as the heat transfer coefficient to turbulent cryogenic
70 Chapter 2
Table 2.10 Estimate Heat Leak for Vacuum Jacketed Piping, Both Rigid and Flexible.
(Courtesy of PHPK Technologies, Inc.)
flow are both negligible to that of the radiation heat load through the vacuum space
per length of pipe. This radiation heat load is a first order effect relating to surface
areas, the larger the surface area, the larger the heat load. The internal and external
surface areas of pipes are driven by their respective diameters.
Table 2.10 also shows very little increase in heat leak in changing the fluid
transferred. All three fluids have heat leak conditions within about 10% of each other
even though the temperature changed from 90 to 20 K. From simple radiation heat
transfer calculations alone, one might expect a heat leak increase by a factor of over
400 (¼(90=20)4). However, the actual case is that the hydrogen temperatures do a
much better job of cryopumping than do oxygen or nitrogen cooled surfaces. Thus
at 90 K, there is still some small residual gas conduction, which is not present at
20 K. Therefore, the apparent thermal conductivities are not that much different.
Table 2.11 illustrates using vacuum jacket valves the cooldown loss in Btu for
various sizes at two temperatures, liquid nitrogen and liquid hydrogen, and then
compares the cooldown loss to the steady-state heat leaks of these units. The heat
leak into vacuum jacketed (VJ) valves is expected to be higher due to the (1) increase
of the penetrations into the VJ around the valve body, (2) the required increased
Table 2.11 Estimated Performance Data for Different Size Vacuum Jacketed Valves.
(Courtesy of PHPK Technologies, Inc.
Valve Size 80 K 20 K 80 K 20 K
1
2" 81 86 10.5 (3.1) 11.6 (3.4)
3
4" 81 86 10.5 (3.1) 11.6 (3.4)
1" 99 105 10.5 (3.1) 11.6 (3.4)
112" 189 200 13.7 (4.0) 15.0 (4.4)
2" 324 342 18.1 (5.3) 20.0 (5s.9)
3" 990 1045 39.7 (11.6) 44.1 (12.9)
4" 1260 1330 47.9 (14.0) 53.1 (15.6)
Basic Principles 71
support for the internal valve components, (3) heat from the power source to drive
the valve, (4) instrumentation required to operate and monitor the valve, etc. Addi-
tional heat leak considerations should be included for the mechanical connections
between VJP segments and the components included in the fluid systems. Each indi-
vidual connection can be roughly estimated to be about 1.24 times a one-foot seg-
ment of the same size line.
The second design problem is pressure drop in cryogenic transfer lines. Two
figures (Figs. 2.19 and 2.20) are provided to show the expected pressure drop per foot
of pipe vs. flow in gallons per minute for nitrogen and hydrogen. Inside pipe dia-
meter is the parameter for the data lines. Additionally, Fig. 2.21 is provided and
illustrates the pressure drop differences of rigid pipe vs. flexible pipe for flowing
liquid nitrogen. The increase in pressure drop can be directly attributed to the inter-
nal design use of welded bellows to form the interior pipe walls. This is true in all, not
just cryogenic, fluid transfer involving flexible lines.
The values shown in Table 2.10 and the graphs of Figs. 2.19 and 2.20 can be
expected from currently available, commercially fabricated transfer lines using
current vacuum-jacketing technology. These values have been confirmed by the two
largest US manufacturers of this equipment. At their request, no credit is given, as
they have no control over how the data will be used and accordingly wish to assume
no responsibility for potential misuse by others.
72 Chapter 2
Example 2.4. Given the following data, calculate the heat leak and pressure drop
for liquid hydrogen in a commercially available vacuum-insulated transfer line:
Liquid hydrogen transfer rate: 250 gpm
Maximum allowable heat leak: 800 Btu=hr
Acceptable pressure drop: 20–30 psi
Desirable pipe size: 2 in. i.d.
Heat leak. Table 2.10 shows that a 2 in. pipe will have a heat leak of 0.85 Btu=hr per
foot of pipe in which liquid hydrogen is flowing. The various fittings (spacers, valves,
bayonet connectors, etc.) will each amount to 1.24 equivalent feet of pipe.
Therefore, if 800 Btu=hr can be accommodated, the maximum run of equiva-
lent pipe is 800=0.85 ¼ 941 ft of pipe. If there were 24 components (a large number),
their equivalent length would be 24 1.24 ¼ 41 ft. Thus, the pipe run could comprise
900 ft of pipe and 24 components.
Pressure drop. Figure 2.20 for hydrogen flow shows a pressure drop of 0.034 psi per
foot of 2 in. pipe.
If the tolerated pressure drop is 20 psi, then the maximum length of pipe is
20=0.034 ¼ 588 ft of pipe.
If the acceptable pressure drop is 30 psi, then the maximum length of pipe is
30=0.034 ¼ 882 ft of pipe.
Basic Principles 73
Figure 2.21 Pressure drop estimates for rigid and flexible transfer lines with liquid nitrogen.
(Courtesy of PHPK Technologies, Inc.).
Since heat transfer and pressure drop follow Rayleigh’s equivalence, a safe
assumption is that the pressure drop caused by a component is about 1.24 times that
of an equivalent length of pipe. Let us pick a pipe run of 500 ft of actual pipe with 24
components in it. The 24 components represent 24 1.24 ¼ 30 equivalent feet of
pipe. The heat leak for this case is then
ð530 ft of equivalent pipeÞ ½0:85Btu=hr ftÞ ¼ 450:5 Btu total heat leak
The pressure drop for this case is
ð530 ft equivalent pipeÞ ð0:034 psi=ftÞ ¼ 18 psi
Summary. For this sample case, 500 ft of 2 in. I.D. pipe carrying liquid hydrogen with
24 fittings will have a heat leak of 450 Btu and a pressure drop of 18 psi, which is well
within our stated parameters. Other cases are readily calculated in a similar manner.
6. COOLDOWN
The amount of liquid required for cooldown can be estimated rapidly by equating
the sensible heat of a pipeline and its insulation to the heat of vaporization of the
fluid. This, of course, does not consider heat leaking into the system, frictional losses,
74 Chapter 2
or changing temperature of the gas discharge from the system. Should any one of
these items contribute a sizable heat load to the system, then the amount of liquid
required for cooldown would obviously be low. In the final analysis, economics will
normally be the most important factor in the selection of a piping system. Such
economics will depend largely on the duration of the transfer and the time
between transfers. These times must be evaluated rather closely if the economic
evaluation is to be valid.
In Chapter 4, we will see that the heat capacity of all crystalline materials goes
to zero as the absolute temperature goes to zero. Accordingly, the amount of cryo-
genic fluid needed for cooldown also goes to zero as the initial temperature of the
solid goes to zero. In other words, prechilling a liquid hydrogen system with liquid
nitrogen will save a lot of liquid hydrogen.
Figures 2.22–2.25 give cooldown values for nitrogen, oxygen, hydrogen, and
helium. They are based on the heat capacity of aluminum (Al); stainless steel (SS),
typically a 304 series; and copper. They all show that the amount of coolant per mass
of metal decreases rapidly as the initial cooldown temperature decreases. For each
metal, there is a pair of MAX and MIM curves. The MAX curve is the maximum
amount of coolant required if only the heat of vaporization of the liquid is available
to cool the metal. This is the worst case. The MIN curve is the minimum amount of
coolant required if both the heat of vaporization of the respective liquid and all of the
sensible heat capacity of the liquid turned gas are available to cool the metal. This is
the best case and will occur only if all of the boil-off gas is somehow in perfect heat
exchange with the metal being cooled, such that the boil-off gas vents at the same
temperature as the warm metal. The vent gases must move slowly past the metal
in order to achieve this ‘‘best case’’ cooling. In other words, all of the possible refrig-
eration is captured; none is wasted. Only by very careful and expensive design, such
as in the cooling of some space satellites, is the minimum ever achieved.
Figure 2.22 Cooldown to liquid nitrogen temperatures. Mass of liquid nitrogen required per
mass of metal vs. starting temperature.
Basic Principles 75
Figure 2.23 Cooldown to liquid oxygen temperatures. Mass of liquid oxygen required per
mass of metal vs. starting temperature.
It is instructive to see just how much precooling a liquid hydrogen system with
liquid nitrogen can help. Referring to Fig. 2.2 for hydrogen, we see that starting at
540 R (80 F), it will take (MAX case) about 0.2 1b of liquid hydrogen to cool 1 lb of
stainless steel to 36 R (20 K). Now suppose that the stainless steel system is pre-
cooled with liquid nitrogen to 139 R (77 K). The same MAX curve now predicts that
only 0.01 1b of liquid hydrogen will be needed per pound of steel to cool it the rest of
the way to 36 R (77 K). The liquid hydrogen use has been cut by a factor of 20. The
price we pay (see Fig. 2.19) is that 0.4 1b of liquid nitrogen is used per pound of metal
Figure 2.24 Cooldown to liquid hydrogen temperatures. Mass of liquid hydrogen required
per mass of metal vs. starting temperature.
76 Chapter 2
Figure 2.25 Cooldown to liquid helium temperatures. Mass of liquid helium required per
mass of metal vs. starting temperature.
for precooling from 540 R (300 K) to 139 R (77 K). Such a trade is very economical
indeed.
7. SUMMARY
As one would expect, the basic principles of thermodynamics, heat transfer, and
momentum transfer apply equally as well to cryogenic engineering as they do to
all other fields of engineering. Appropriate use must be made of the sometimes unu-
sual low-temperature properties of fluids and materials (see Chapter 3 and 4).
Simple engineering correlations, usually based on the behavior of a large
number of near-room-temperature fluids, must be used with caution. Pressure drop
correlations, for instance, are usually based on the flow of single-phase streams
under steady-state conditions. Thus they may underestimate the pressure drop in
cryogenic systems by 10–30% due to the high vaporization tendency of cryogenic
fluids.
Finally, the usual modeling assumptions for normal engineering work may not
apply. Cryogenic fluids are nearly always at their boiling point, are seldom adiabatic,
and are often composed of two phases. Radiative heat transfer can seldom be
neglected.
This is not to say that cryogenic systems cannot be rigorously modeled, for
they can be, if one takes into account all of the peculiar properties of cryogenic
materials and properties. Rather, correlations must be used with caution, and a
return to basics is usually the safest approach.
3
Cryogenic Fluids
The relationship of specific or molar volume to temperature and pressure for a pure
substance in an equilibrium state can be represented by a surface in three
dimensions, as shown in Fig. 3.1. The shaded surfaces marked S, L, and G in
Fig. 3.1 represent, respectively, the solid, liquid, and gas regions of the diagram.
The unshaded surfaces are the regions of co-existence of two phases in equilibrium,
and there are three such regions: solid–gas (S–G), solid–liquid (S–L), and liquid–gas
(L–G).
Heavy lines separate the various regions and form boundaries of the surfaces
representing the individual phases. Because of the difficulty of visualizing and draw-
ing these surfaces, it is customary to depict the data on projections such as the PT
and PV planes. Figure 3.2 is an example of a PT plane for water; other materials
have similar diagrams.
The heavy line passing through points A and B marks the intersections of the
two-phase regions and is the three-phase line, along which solid, liquid, and gas
phases exist in three-phase equilibrium. According to the phase rule, such three-
phase, single-component systems have zero degrees of freedom; they exist for a given
pure substance at only one temperature and one pressure. For this reason, the pro-
jection of this line on the PT plane of Fig. 3.2 (shown to the left of the main diagram)
is a point. This is known as a triple point.
The triple point is merely the point of intersection of the sublimation and
vaporization curves. It must be understood that only on a P–T diagram is the triple
point represented by a point. (On a P–V diagram it is a line.)
The phase rule also requires that systems made up of two phases in equilibrium
have just one degree of freedom, and therefore the two-phase regions must project as
lines on the PT plane, forming a P–T diagram of three lines—fusion (or melting),
sublimation, and vaporization—which meet at the triple point.
The P–T diagram is divided into areas by three lines: the solid–vapor, solid–
liquid, and liquid–vapor lines. The triple point is a region in which liquid, solid,
and vapor can exist in equilibrium simultaneously.
The P–T projection of Fig. 3.2 is shown to a larger scale in Fig. 3.3.
Upon projecting the surface on the PT plane, the whole solid–vapor region
projects into the sublimation curve, the whole liquid–vapor region projects into
the vaporization curve, the whole solid–liquid region projects into the fusion curve,
and finally the ‘‘triple point line’’ projects into the triple point.
77
78 Chapter 3
The solid lines of Fig. 3.3 clearly represent phase boundaries. The fusion curve
(line 2–3) normally has a positive slope, but for a few substances (water is the best
known) it has a negative slope. This line is believed to continue upward indefinitely.
The two curves 1–2 and 2–C represent the vapor pressures of solid and liquid,
respectively.
The terminal point C is the critical point, which represents the highest pressure
and highest temperature at which liquid and gas can coexist in equilibrium. The
termination of the vapor pressure curve at point C means that at higher temperatures
and pressures no clear distinction can be drawn between what is called liquid and
what is called gas. At the critical point, there is no distinction between the liquid
and vapor phases.
Cryogenic Fluids 79
Homogeneous fluids are normally divided into two classes, liquids and gases.
However, the distinction cannot always be sharply drawn, because the two phases
become indistinguishable at the critical point.
Thus, there is a region extending indefinitely upward from and indefinitely to
the right of the critical point that is called simply the fluid region. The fluid region,
existing at higher temperatures and pressures, is marked off by dashed lines that do
not represent phase changes but depend on arbitrary definitions of what constitutes
liquid and gas phases.
The region designated as liquid in Fig. 3.3 lies above the vaporization curve.
Thus, a liquid can always be vaporized by a sufficient reduction in pressure at
constant temperature. The gas region of Fig. 3.3 lies to the right of the sublimation
and vaporization curves. Thus, a gas can always be condensed by a sufficient
reduction in temperature at constant pressure. Because the fluid region fits neither
of these definitions, it can be called neither a gas nor a liquid.
The gas region is sometimes considered to be divided in two parts, as shown by
the dash-dot vertical line of Fig. 3.3. A gas to the left of this line, which can be
condensed either by compression at constant temperature or by cooling at constant
pressure, is frequently called a vapor. A vapor is a gas existing at temperatures below
Tc, and it may therefore be condensed either by reduction of temperature at constant
P or by increase of pressure at constant T.
The P–T projection of Fig. 3.3 provides no information about the volumes of the
systems represented. This is, however, given explicitly by the projection of Fig. 3.4,
where all surfaces of the three-dimensional diagram show up as areas on the PV plane.
The liquid and gas regions of the PV projection are shown in more detail by
Fig. 3.5.
Since the PV plane is perpendicular to the PT plane in the P,V,T, coordinate
system, the solid–vapor, liquid–vapor, and solid–liquid areas of the PV plane appear
as lines on the PT plane or phase diagram.
The primary curves of Fig. 3.5 give the pressure–volume relations for saturated
liquid (A–C) and for saturated vapor (C–B). The area lying below the curve ACB
80 Chapter 3
Figure 3.4 Projection of the P-V-T surface of a pure fluid on PV and PT planes.
represents the two-phase region where saturated liquid and saturated vapor coexist at
equilibrium. Point C is the critical point, for which the coordinates are Pc and Vc.
A number of constant-temperature lines or isotherms are included on Fig. 3.5.
The critical isotherm at temperature Tc passes through the critical point. Isotherms
Cryogenic Fluids 81
for higher temperatures T1 and T2 lie above the critical isotherm. Such isotherms do
not cross a phase boundary and are therefore smooth.
Lower temperature isotherms T3 and T4 lie below the critical isotherm and are
made up of three sections. The middle section traverses the two-phase region and is
horizontal, because equilibrium mixtures of liquid and vapor at a fixed temperature
have a fixed pressure, the vapor pressure, regardless of the proportions of liquid and
vapor present. Points along this horizontal line represent different proportions of
liquid and vapor, ranging from all liquid at the left end to all vapor at the right.
The locus of these end points is represented by the dome-shaped curve labeled
ACB, the left half of which (from A to C) represents saturated liquid, and the right
half (from C to B), saturated vapor.
The horizontal line segments become progressively shorter with increasing tem-
perature, until at the critical point the isotherm exhibits a horizontal inflection. As
this point is approached, the liquid and vapor phases become more and more nearly
alike, eventually becoming unable to sustain a meniscus between them and ultimately
becoming indistinguishable.
The sections of isotherms to the left of line AC traverse the liquid region and
are very steep, because liquids are not very compressible; that is, it takes large
changes in pressure to effect a small volume change. A liquid that is not saturated
(not at its boiling point, along line AC) is often called subcooled liquid or
compressed liquid.
The sections of isotherms to the right of line CB lie in the vapor region. Vapor
in this region is called superheated vapor to distinguish it from saturated vapor as
represented by the line CB.
82 Chapter 3
Cryogenic Fluids 83
Figure 3.8 T–S diagram phase boundaries near the triple point.
Property Oxygen Nitrogen Neon e-Hydrogena Helium-4 Air Fluorine Argon Methane
Normal boiling point (K) 90.18 77.347 27.09 20.268 4.224 78.9 85.24 87.28 111.7
Density (kg=m3) 1141 808.9 1204 70.78 124.96 874 1506.8 1403 425.0
Heat of vaporization (kJ=kg) 212.9 198.3 86.6 445.6 20.73 205.1 166.3 161.6 511.5
Specific heat (kJ=(kg K)) 1.70 2.04 1.84 9.78 4.56 1.97 1.536 1.14 3.45
Viscosity (kg=(m s) 106) 188.0 157.9 124.0 13.06 3.57 168 244.7 252.1 118.6
Thermal conductivity(mW=(m K)) 151.4 139.6 113 118.5 27.2 141 148.0 123.2 193.1
Dielectric constant 1.4837 1.434 1.188 1.226 1.0492 1.445 1.43 1.52 1.6758
Critical temperature (K) 154.576 126.20 44.4 32.976 5.201 133.3 144.0 150.7 190.7
Critical pressure (MPa) 5.04 3.399 2.71 1.293 0.227 3.90 5.57 4.87 4.63
Temperature at triple point (K) 54.35 63.148 24.56 13.803 – – 53.5 83.8 88.7
Pressure at triple point (MPa 103) 0.151 12.53 43.0 7.042 – – 0.22 68.6 10.1
a
For a discussion of forms of hydrogen, see Sec. 3.7.
85
© 2005 by Marcel Dekker
Table 3.2 Properties of Cryogenic Liquids at Their NBP, Critical and Triple Points
86
Property Helium Hydrogen Neon Nitrogen Air Fluorine Argon Oxygen Methane
Molecular weight 4.0026 2.01594 20.1790 28.0134 28.975 37.9968 39.9480 31.9988 16.0420
Normal boiling point, K 4.2 20.4 27.1 77.3 78.9 84.95 87.3 90.2 111.7
Normal boiling point, R 7.57 36.7 48.8 139.2 141.8 152.9 157.1 162.4 201.1
NBP Liquid density, kg=m3 124.9 70.9 1204 810.8 874.0 1502 1403 1134 425.0
NBP Liquid density, lbm=ft3 7.80 4.43 75.17 50.61 54.56 93.7 87.56 70.8 26.53
NBP Vapor density, kg=m3 16.845 1.3395 9.5924 4.624 4.487 5.629 5.7765 4.4668 1.8176
NBP Vapor density, lbm=ft3 1.0516 0.08362 0.59883 0.28866 0.28011 0.3513 0.36062 0.27885 0.11346
NBP Heat of vaporization, kJ=kg 20.7 446.3 85.71 198.8 205.1 157.4 161.6 213.1 510.33
NBP Heat of vaporization, Btu=lbm 8.92 191.9 36.849 85.32 88.2 67.8 69.5 91.63 219.2
Critical temperature, K 5.2 33.2 44.4 126.12 133.3 144.31 150.9 154.6 190.7
Critical temperature, R 9.36 59.74 79.9 227.0 240 259.8 271.2 278.3 343.3
Critical pressure, MPa 0.23 1.31 2.71 3.38 3.90 5.22 4.87 5.06 4.63
Critical pressure, Atm 2.26 12.98 26.84 33.5 38.7 51.47 48.3 50.1 45.8
Triple point, K 2.177(l) 13.9 24.56 63.2 60.6 53.5 83.8 54.4 88.7
Triple point, R 3.92(l) 25.1 44.2 113.7 109.1 96.37 150.8 98.0 159.7
Pressure at triple point (MPa)103 4.73(l) 7.2 43.0 12.8 7.035 0.25 68.6 0.15 10.1
Pressure at triple point, atm 0.0469(l) 0.0711 0.426 0.1268 0.0697 0.00249 0.679 0.00150 0.099
Gas density NTP 1 atm, 70 F
r lb=f3 0.010343 0.0052089 0.052148 0.072445 0.07493 0.09803 0.10355 0.082787 0.041552
r kg=m3 (g=L) 0.16569 0.083438 0.83532 1.1604 1.2000 1.5704 1.6555 1.3261 0.66559
Volume ratio
GAS NTP; 1 atm 70 F 754.2 848.5 1433.7 695.3 728.1 956.7 842.0 860.6 634.7
Liquid NBP
Ft3 GAS NTP; 1 atm 70 F 26.63 29.96 50.9 24.55 25.71 33.78 29.70 30.39 22.41
Chapter 3
Liter Liq. NBP
Cryogenic Fluids 87
Table 3.3 Enthalpy of Phase Changes and Density of Five Common Cryogenic Fluids
Fluid He P-H2 N2 O2 C1
Molecular weight 4.0026 2.01594 28.0134 31.9988 16.0420
Figure 3.11 shows the vapor pressure of several cryogenic liquids. Again, the
unusual liquid range of oxygen is shown.
Figure 3.12 gives the unusual correlation between the NBP and critical tem-
perature for several cryogenic fluids. In the absence of data, a rule of thumb is that
Tc is about 1.6 times TNBP or that the NBP is about 60% of the critical temperature.
Helium, of course, marches to its own drummer. It is surprising that the other quan-
tum fluids, hydrogen and neon, are so well behaved.
Density
MW K
F
C atm psi mm lb=ft3 g=L lb=ft3 g=L lb=ft3 g=L
Density
MW K
F
C lb=ft3 g=L lb=ft3 g=L BTU=lb J=g
Chapter 3
O2 31.9988 90.191 297.33 182.76 71.240 1141.2 0.27885 4.4668 91.597 213.05
CH4 16.0420 111.69 258.63 161.46 26.372 422.44 0.11347 1.8176 219.40 510.33
MW K
F
C atm psia lbs=ft3 g=l
Critical point
He 4.0026 5.1953 450.32 267.95 2.2453 32.996 4.3476 69.641
p-H2 2.01594 32.938 400.38 240.21 12.670 186.20 1.9577 31.360
N2 28.0134 126.19 232.52 146.96 33.534 492.81 19.547 313.11
O2 31.9988 154.58 181.42 118.57 49.771 731.43 27.228 436.14
CH4 16.0420 190.55 116.68 82.6 45.391 667.06 10.154 162.65
89
© 2005 by Marcel Dekker
5367-4 Flynn Ch03 R2 090804
90 Chapter 3
Table 3.5 Weight Volume Equivalents for Air, Argon, Nitrogen, Oxygen, Helium and
Hydrogen
Weight lb Gallons L Cu ft Cu M L
Volume Equivalents
1 atm, 70 F
Gas, 1 atm 70 F Ft3 Gas 1 atm, 70 F
3
MW lb=ft g=L Liquid NBP Liter liquid NBP
Cryogenic Fluids 91
Weight of liquid or gas Volume of liquid or normal boiling point Volume of gas at 70 F
lb kg cu ft l quarts gallons cu ft cu m
Oxygen
1.000 0.454 0.0140 0.396 0.419 0.105 12.08 0.342
2.205 1.000 0.0308 0.873 0.923 0.231 26.63 0.754
71.5 32.43 1.000 28.32 29.92 7.481 863.5 24.45
2.524 1.145 0.0353 1.000 1.056 0.264 30.48 0.863
2.388 1.083 0.0334 0.946 1.000 0.250 28.84 0.817
9.560 4.336 0.134 3.786 4.000 1.000 115.5 3.271
8.28 3.756 0.116 3.285 3.472 0.868 100.0 2.832
2.92 1.324 0.0408 1.155 1.221 0.305 35.31 1.000
Air
1.000 0.454 0.0183 0.518 0.548 0.137 13.35 0.379
2.205 1.000 0.0404 1.144 1.209 0.302 29.44 0.834
54.57 24.77 1.000 28.32 29.92 7.481 728.5 20.63
1.927 0.875 0.0353 1.000 1.056 0.264 25.73 0.729
1.824 0.828 0.0334 0.946 1.000 0.250 24.35 0.690
7.294 3.311 0.134 3.786 4.000 1.000 97.37 2.758
7.493 3.402 0.137 3.880 4.099 1.025 100.0 2.832
2.646 1.201 0.0485 1.374 1.449 0.363 35.31 1.000
Argon
1.000 0.454 0.0114 0.323 0.341 0.0853 9.681 0.274
2.205 1.000 0.0252 0.714 0.755 0.189 21.35 0.605
87.4 39.6 1.000 28.32 29.92 7.481 846.1 23.96
3.085 1.399 0.0353 1.000 1.056 0.264 29.86 0.846
2.919 1.324 0.0334 0.946 1.000 0.250 28.26 0.800
11.71 5.310 0.134 3.786 4.000 1.000 113.4 3.212
10.33 4.686 0.118 3.342 3.533 0.883 100.0 2.832
3.648 1.655 0.0417 1.181 1.248 0.312 35.31 1.000
Nitrogen
1.000 0.454 0.0198 0.561 0.593 0.148 13.79 0.391
2.205 1.000 0.0437 1.239 1.310 0.327 30.41 0.861
50.4 22.86 1.000 28.32 29.92 7.481 695.2 19.69
1.78 0.807 0.0353 1.000 1.056 0.264 24.55 0.695
1.68 0.762 0.0334 0.946 1.000 0.250 23.17 0.656
6.75 3.062 0.134 3.786 4.000 1.000 93.10 2.637
7.25 3.289 0.144 4.078 4.310 1.077 100.0 2.832
2.56 1.161 0.0508 1.439 1.521 0.380 35.31 1.000
Helium
1.000 0.454 0.128 3.625 3.832 0.958 97.08 2.749
2.205 1.000 0.283 8.00 8.468 2.117 214.1 6.063
7.80 3.541 1.000 28.32 29.92 7.481 757.3 21.45
0.2753 0.125 0.0353 1.000 1.056 0.264 26.73 0.757
0.2605 0.118 0.0334 0.946 1.000 0.250 25.29 0.716
1.045 0.474 0.134 3.786 4.000 1.000 101.5 2.875
1.033 0.469 0.132 3.738 3.948 0.987 100.0 2.832
0.3637 0.165 0.0466 1.320 1.396 0.349 35.31 1.000
(Continued)
92 Chapter 3
Weight of liquid or gas Volume of liquid or normal boiling point Volume of gas at 70 F
lb kg cu ft l quarts gallons cu ft cu m
Hydrogen
1.000 0.454 0.226 6.400 6.764 1.691 191.6 5.426
2.205 1.000 0.498 14.10 14.90 3.726 422.4 11.96
4.43 2.011 1.000 28.32 29.92 7.481 848.7 24.04
0.156 0.0708 0.0353 1.000 1.056 0.264 29.89 0.846
0.148 0.0672 0.0334 0.946 1.000 0.250 28.35 0.803
0.594 0.270 0.134 3.786 4.000 1.000 113.8 3.223
0.522 0.237 0.118 3.342 3.532 0.883 100.0 2.832
0.184 0.0835 0.0416 1.178 1.244 0.311 35.31 1.000
3.1. Oxygen
Oxygen was the base used for chemical atomic weights, being assigned the atomic
weight 16.000, until 1961, when the International Union of Pure and Applied Chem-
istry (IUPAC) adopted carbon-12 as the new basis. Oxygen has eight isotopes. Natu-
rally occurring oxygen consists of three stable isotopes of atomic mass numbers 16,
17, and 18, having abundances in the proportion 10,000:4:20. Oxygen condenses into
a light blue liquid whose density, 1134.2 kg=m3 (70.8 lbm=ft3) at its boiling tempera-
ture, is slightly greater than that of water at room temperature. At 0.987 bar (1 atm)
pressure, liquid oxygen boils at 90.2 K (162.4 R) and freezes at 54.4 K (97.9 R).
In the liquid state, there is thought to be some weak transient association of
oxygen molecules forming O4, which is said to be responsible for the blue color of
both liquid and solid oxygen. Ozone (O3), a highly active allotropic from of oxygen,
is formed by the action of an electrical discharge or ultraviolet light on oxygen.
Ozone’s presence in the atmosphere (amounting to the equivalent of a layer 3 mm
thick at ordinary pressures and temperatures) is of vital importance in preventing
harmful ultraviolet rays of the sun from reaching the earth’s surface.
Table 3.7 lists the temperature, pressure, and specific volume of each of the
various phases of oxygen at several important points. Tables 3.8–3.10 give properties
of oxygen in metric and English units at three of these fixed points.
Figures 3.14 and 3.15 and T–S diagrams for oxygen. A Mollier chart (pressure
vs. heat content) is given in Fig. 3.16, and the Joule–Thomson inversion chart in Fig.
3.17. Table 3.11 presents data for the J–T curve.
94 Chapter 3
Figure 3.11 Vapor pressure curves. CP ¼ critical point; TP ¼ triple point; l ¼ lambda point
(see Sec. 3.8.1).
Figure 3.12 Ratio of the critical temperature to the NBP for various cryogens. For C terms,
see Fig. 3.11.
Cryogenic Fluids 95
Critical point
T ¼ 154.576 0.010 K T ¼ 154.581 278.237 R
P ¼ 49.76 0.02 atm (731.4 psia; 5.043 MPa)
V ¼ 73.37 0.10 cm3=mol (0.03673 ft3=lbm; 436.2 kg=m3)
Normal boiling point
T ¼ 90.180 0.01 K T ¼ 90.188 162.324 R
P ¼ 1 atm (14.696 psia; 101.3 kPa)
V (liquid) ¼ 28.05 0.028 cm3=mol (0.01404 ft3=lbm; 1141 kg=m3)
V (vapor) ¼ 7150 7 cm3=mol (3.579 ft3=lbm; 4.48 kg=m3)
Normal melting point
T ¼ 54.362 0.001 K T ¼ 54.372 97.852 R
P ¼ 1 atm (14.696 psia; 101.3 kPa)
V (liquid) ¼ 24.49 0.024 cm3=mol (0.01226 ft3=lbm; 1307 kg=m3)
Triple point
T ¼ 54.351 0.001 K T ¼ 54.361 97.832 R
P ¼ (1.50 0.06) 103 atm (0.0220 psia; 151.7 Pa)
V (solid) ¼ 23.55 0.03 cm3=mol (0.01179 ft3=lbm; 1359 kg=m3)
V (liquid) ¼ 24.49 0.024 cm3=mol (0.01226 ft3=lbm; 1307 kg=m3)
V (vapor) ¼ 2.97 0.01 106 cm3=mol (1487 ft3=lbm; 10.8 g=m3)
Solid-solid transitions. There are two such transitions at atmospheric pressure:
at 43.8 K and at 23.9 K
96 Chapter 3
Oxygen
Chemical Symbol: O2
CAS Registry Number: 7782-44-7
DOT Classification: Nonflammable gas
DOT Label: Oxidizer
Transport Canada Classification: 2.2 (5.1)
UN Number: UN 1072 (compressed gas); UN 1073 (refrigerated liquid)
US Units SI Units
Physical constants
Chemical formula O2 O2
Molecular weight 31.9988 31.9988
Density of the gas at 70 F (21.1 C) and 1 atm 0.08279 lb=ft3 1.326 kg=m3
Specific gravity of the gas at 70 F (21.1 C) 1.105 1.105
and 1 atm (air ¼ 1)
Specific volume of the gas at 70 F (21.1 C) 12.08 ft3=lb 0.7541 m3=kg
and 1 atm
Boiling point at 1 atm 297.33 F –182.96 C
Freezing point at 1 atm 361.80 F 218.78 C
Critical temperature 181.43 F 118.57 C
Critical pressure 731.4 psia 5043 kPa abs
Critical density 27.22 lb=ft3 436.1 kg=m3
Triple point 361.82 F at 218.79 C at
0.02147 psia 0.1480 kPa abs
Latent heat of vaporization at boiling point 91.7 Btu=lb 213 kJ=kg
Latent heat of fusion at 361.1 F (218.4 C) 5.959 Btu=lb 13.86 kJ=kg
melting point
Specific heat of the gas at 70 F (21.1 C) and
1 atm
CP 0.2197 Btu=(lb)( F) 0.9191 kJ=(kg)( C)
CV 0.1572 Btu=(lb)( F) 0.6578 kJ=(kg)( C)
Ratio of specific heats (CP=Cv) 1.40 1.40
Solubility in water, vol=vol at 32 F (0 C) 0.0491 0.0491
Density of the liquid at boiling point 9.52 lb=gal or 1141 kg=m3
71.23 lb=ft3
Density of the gas at boiling point 0.2799 lb=ft3 4.483 kg=m3
Gas=liquid ratio (gas at 70 F (21.1 C) and 860.5 860.5
1 atm, liquid at boiling point, vol=vol)
Cryogenic Fluids
Triple point Normal boiling point Standard conditions
Property Solid Liquid Vapor Liquid Vapor Critical point STP(0 C) NTP(20 C)
Temperature (K) 54.351 54.351 54.351 90.180 90.180 154.576 273.15 293.15
Pressure (mmHg) 1.138 1.138 1.138 760 760 37,823 760 760
Density (mol=cm3) 103 42.46 40.83 0.000336 35.65 0.1399 13.63 0.04466 0.04160
Specific volume (cm3=mol) 103 0.02355 0.02449 2975 0.028047 7.1501 0.07337 22.392 24.038
Compressibility factor, Z ¼ PV=RT — 0.0000082 0.9986 0.00379 0.9662 0.2879 0.9990 0.9992
Heats of fusion and vaporization (J=mol) 444.8 — 7761.4 6812.3 — 0 — —
Specific heat (J=(mol K))
Cg, at saturation 46.07 53.313 –108.7 54.14 –53.2 (very large) — —
CP, at constant pressure — 53.27 29.13 54.28 30.77 (very large) 29.33 29.40
CV at constant volume — 35.65 20.81 29.64 21.28 (38.7) 20.96 21.04
Specific heat ratio, g ¼ CP=CV — 1.494 1.400 1.832 1.446 (large) 1.40 1.40
Enthalpy (J=mol) 6634.46 6189.6 1571.8 4270.3 2542.0 1032.2 7937.8 8525.1
Internal energy (J=mol) 6634.4 6189.6 1120.0 4273.1 1817.5 662.3 5668.9 6089.5
Entropy (J=(mol K)) 58.92 67.11 209.54 94.17 169.68 134.42 202.4 204.5
Velocity of sound (m=s) — 1159 141 903 178 164 315 326
Viscosity
(N sec=m3) 103 — 0.6194 0.003914 0.1958 0.00685 (0.031) 0.01924 0.02036
Centipoises — 0.6194 0.003914 0.1958 0.00685 (0.031) 0.01924 0.02036
Thermal conductivity, k(mW=(cm K)) — 1.929 0.04826 1.515 0.08544 —b 0.2428 0.2575
Prandtl number, Npr ¼ uCp=k — 5.344 0.7392 2.193 0.7714 — 0.7259 0.7265
Dielectric constant, e pffiffiffiffi (1.614) 1.5687 1.000004 1.4870 1.00166 1.17082 1.00053 1.00049
Index of refraction, n ¼ ea (1.271) 1.2525 1.000002 1.219 1.00083 1.0820 1.00027 1.00025
Surface tension (N=m) 103 — 22.65 — 13.20 — 0 — —
Equiv. vol.=vol. liquid at NBP 0.8397 0.8732 106,068 1 254.9 2.616 798.4 857.1
a
Long wavelengths.
b
Anomalously large.
Gas constant: R ¼ 62,365.4 cm3 mmHg=(mol K).
97
Values in parentheses are estimates based on the NBS-1955 temperature scale using 90.180 K as a fixed point for the normal boiling temperature of oxygen. (The IPTS-1968
temperature scale uses 90.188 K as the normal boiling temperature of oxygen, but the reported data used in these tables have not yet been converted to the new temperature
scale.) 273.15 K ¼ 0 C ¼ 32 F ¼ 491.67 R base point (zero values) for enthalpy, internal energy, and entropy are 0 K for the ideal gas at 1 atm pressure.
98
Property Solid Liquid Vapor Liquid Vapor Critical point STP(32 F) NTP(68 F)
Chapter 3
temperature scale used 90.188 K as the normal boiling temperature of oxygen, but the reported data used in these tables have not yet been converted to the new temperature
scale.) 273.15 K ¼ 0 C ¼ 32 F ¼ 491.67 R. Calculated from property values given in this table. Base point (zero values) for enthalpy, internal energy, and entropy are 0 K for the
ideal gas at 1 atm pressure.
Cryogenic Fluids 99
100 Chapter 3
101
© 2005 by Marcel Dekker
5367-4 Flynn Ch03 R2 090804
102 Chapter 3
Range of values:
ln P ¼ A1 þ A2 T þ A3 T 2 þ A4 T 3 þ A5 T 4 þ A6 T 5 þ A7 T 6 þ A8 T 7 ð3:1Þ
where
104 Chapter 3
The estimated uncertainty is 0.02% from the normal boiling point to the critical
point. Below the normal boiling point the uncertainty increases, reaching about 3%
at the triple point. The equivalent uncertainty is 0.02 K in the range of 60–154.58 K.
Figures 3.19 and 3.20 are graphs of the vapor pressure of oxygen.
f. Virial Coefficients. The virial coefficients (Table 3.12) are usually defined
from the virial equation in density,
e1
¼ Ar þ Br2 þ Cr3 ð3:3Þ
eþ2
where A ¼ 0.12361, B ¼ 3.2 104, C ¼ 1.21 103, and r is in g=cm3. The equa-
tion is valid over the range 54.35–340 K, 0.2–340 atm, and will yield reasonable
values upon extrapolation.
Figure 3.19 Vapor pressure of oxygen in atm and mmHg vs. temperature in K.
The uncertainty in the term e 1 varies from 0.15% at low densities to less than
0.05% at high densities.
h. Other Properties. Tables 3.13–3.15 give values for the viscosity of gaseous
and liquid oxygen as well as the density of saturated liquid oxygen at various tem-
peratures. Figures 3.21–3.29 were compiled by the National Bureau of Standards,
Cryogenic Engineering Laboratory for the US Air Force, Wright Air Development
Division in WADD Technical Report 60–56, V. J. Johnson ed., October 1960. They
are commonly referred to as the WADD data. In Figs. 3.21–3.29 you will find best
values of the melting curve, dielectric constant, density, specific heats, and viscosity,
for oxygen.
3.1.2. Uses of Oxygen
Commercial oxygen demand in the United States is about 45 109 kg (50 million
tons) per year, at an average growth rate of 5.1% annually. For the most part,
106 Chapter 3
tonnage oxygen (that delivered via pipeline) is consumed in steel plants. Present con-
sumption in other industries has seldom reached levels requiring individual oxygen
plants; those industries are usually supplied with liquid ‘‘merchant oxygen’’, which
is the name given to liquid oxygen delivered to the consumer by truck. The tonnage
oxygen plants currently fill about 80% of the oxygen demand.
The first continuous machine for liquefying air was built by Karl von Linde of
Munich, Germany, in 1895. The first liquid air installation of commercial size in the
United States was built for the American Linde Co., now Linde Division of Union
Carbide, in Niagara Falls, New York, in 1907. It had a capacity of 21,240 m3
(750,000 cu ft) per month of oxygen, equivalent to about 907 kg (1 ton) per day.
There is no substitute for oxygen in any of its uses, and it is not recycled or
reclaimed.
a. Steel. About 60–70% of the oxygen produced in the United States is
consumed in the manufacture of iron and steel. Concurrently with the development
of pelletizing and other improvements in blast furnace practice, low-cost tonnage
oxygen became available in the 1950s and was instrumental in the development of
the basic oxygen furnace (BOF), which is a converter using oxygen instead of air.
b. Chemical. The chemical industry uses about 12% of the total oxygen
produced. Major uses are in partial oxidation of methane to produce acetylene; in
the oxidation of ethylene (the chief hydrocarbon produced in the United States) to
Temp. (K) Viscosity (cP) 103 Temp. (K) Viscosity (cP) 103
108 Chapter 3
The T–S diagram for oxygen (Fig. 3.30) contains much more information than
first meets the eye. The answers to the above questions are readily found on the
oxygen T–S chart.
First, find the 150 psia isobar for oxygen and follow it to the two–phase region.
The isobars and isotherms are coincident horizontal lines under the two–phase
dome. Therefore, the boiling point of oxygen is read directly from the temperature
scale to be about 216 R. The published value is 216.38 R at 150 psia.
Likewise, the heat of vaporization can be found in two independent ways.
First, it is only necessary to read the enthalpy values for saturated vapor (Hg) and
saturated liquid (Hf). Along the 150 psia isobar we find Hg ¼ 40 Btu=lbm and
Hf ¼ –34 Btu=lbm. The difference is Hfg ¼ Hg Hf ¼ 40 (34) ¼ 74 Btu=lbm. The
published value is 74.39 Btu=lbm.
Heat of vaporization is also TDS. From the TS chart we find Sg ¼
1.167 Btu=(lbm R) and Sf ¼ 0.824 Btu=(lbm R). Accordingly, DS ¼ 1.167 0.824 ¼
0.343 Btu=(lbm R). TDS in appropriate units is (216.4 R) [0.343 Btu=
(lbm R)] ¼ 74.2 Btu=lbm, which compares well with the published value of
74.4 Btu=lbm.
3.2. Nitrogen
Nitrogen (mol wt. 28.0134) has two stable isotopes of mass numbers 14 and 15 with a
relative abundance of 10,000:38. Liquid nitrogen is of considerable importance to the
cryogenic engineer because it is a safe refrigerant. Because it is rather inactive
110 Chapter 3
chemically and is neither explosive nor toxic, liquid nitrogen is commonly used in
hydrogen and helium liquefaction cycles as a precoolant.
Nitrogen is the major constituent of air (78.09% by volume or 75.45% by
weight). The atmosphere of Mars, by comparison, is 2.6% nitrogen. The estimated
amount of this element in out atmosphere is more than 4000 billion tons. From this
inexhaustible source, it can be obtained by liquefaction and fractional distillation.
Liquid nitrogen is a clear, colorless fluid that resembles water in appearance.
At 1 bar pressure, liquid nitrogen boils at 77.3 K (139.2 R) and freezes at 63.2 K
(113.8 R). Saturated liquid nitrogen at the normal boiling point has a density of
808.9 kg=m3 (50.5 lbm=ft3) in comparison with water at 520 R (60 F), which has a
density of 998.0 kg=m3 (62.3 lbm=ft3). One of the significant differences between
the properties of liquid nitrogen and those of water (apart from the differences in
normal boiling points) is that the heat of vaporization of nitrogen (Table 3.16) is
more than an order of magnitude less than that of water. At its normal boiling point,
liquid nitrogen has a heat of vaporization of 198.3 kJ=kg (85.32 Btu=lbm), while
water has a heat of vaporization of 2255 kJ=kg (970.3 Btu=lbm).
Table 3.17 gives the physical constants of nitrogen.
112 Chapter 3
Table 3.18 gives some of the important fixed points for nitrogen.
Figures 3.31–3.33 are TS charts for nitrogen. Figure 3.34 is the Joule–
Thomson inversion curve.
Nitrogen has a potential safety hazard in that a bare pipe of liquid nitrogen at
77 K will condense an air mixture containing approximately 50% liquid oxygen.
Thus, one must exercise extreme caution that so-called inert liquid nitrogen does not
in fact become an unsuspected source of liquid oxygen, with the attendant risks cited
earlier. Several explosions and deaths have been attributed to the phenomenon of
oxygen enrichment of the atmosphere in the presence of liquid nitrogen cooled
surfaces. See Chapter 10.
359:093
log p ¼ a ð3:4Þ
T
334:64
log p ¼ þ 7:577 0:00476 T ð3:5Þ
T
114 Chapter 3
Figure 3.27 Specific heat (CP) of solid oxygen below 4 K (NASA SP3071).
where p is in mmHg and T is in K. (This equation is also valid for liquid nitrogen
from 63.14 to 77.35 K; 96.4–760 mmHg.)
For liquid nitrogen,
255:821
log p ¼ 6:49594 ð3:6Þ
T 6:600
where p is in mm Hg and T is in K for the range 64–78 K
or
316:824
log p ¼ þ 4:47582 0:0071701T þ 2:940 105 T 2 ð3:7Þ
T
where, p is in atm, T is in K and 1–32 atm.
116 Chapter 3
Heat of vaporization, hv
Solid nitrogen
62.00 6775.0
62.0018 6787.4
62.0172 6762.4
Liquid nitrogen
67.9588 5901.6
67.9620 5899.0
68.00 5899.0
73.0913 5739.1
73.0887 5732.1
73.10 5735.2
77.395 5735.2
78.00 5579.4
78.0147 5563.1
78.0153 5571.8
80 1313
85 1266
90 1213
95 1155
100 1086
105 1010
110 918
115 803
120 643
125 328
126.1a 0
a
Critical point.
Source: Furukawa and McCoskey (1953), Millar and Sullivan
(1928).
where e is the dielectric constant, r is the density, and p is the specific polarization, a
property of the substance having dimensions of specific volume:
p ¼ An þ Bn r þ Cn r2 ð3:9Þ
d. Other Properties. Tables 3.19–3.23 present data on the solid and liquid
vapor pressure, density, viscosity of nitrogen. Figures 3.35–3.53, taken from the
WADD report mentioned before, show the heat of vaporization, vapor pressure,
melting curve, density, thermal expansivity, thermal conductivity, specific heats,
dielectric constant, surface tension, and viscosity of nitrogen.
118 Chapter 3
Nitrogen
Chemical Symbol: N2
Synonyms: LIN (liquid only)
CAS Registry Number: 7727-37-9
DOT Classification: Nonflammable gas
DOT Label: Nonflammable gas
Transport Canada Classification: 2.2
UN Number: UN 1066 (compressed gas); UN 1977 (refrigerated liquid)
US Units SI Units
Physical Constants
Chemical formula N2 N2
Molecular weight 28.01 28.01
Density of the gas at 70 F (21.1 C) and 1 atm 0.072 lb=ft3 1.153 kg=m3
Specific gravity of the gas at 70 F (21.1 C) 0.967 0.967
and 1 atm (air ¼ 1)
Specific volume of the gas at 70 F (21.1 C) 13.89 ft3=lb 0.867 m3=kg
and 1 atm
Density of the liquid at boiling point and 1 atm 50.47 lb=ft3 808.5 kg=m3
Boiling point at 1 atm 320.4 F 195.8 C
Melting point at 1 atm 345.8 F 209.9 C
Critical temperature 232.4 F 146.9 C
Critical pressure 493 psia 3399 kPa abs
Critical density 19.60 lb=ft3 314.9 kg=m3
Triple point at 1.81 psia (12.5 kPa abs) 346.0 F 210.0 C
Latent heat of vaporization at boiling point 85.6 Btu=lb 199.1 kJ=kg
Latent heat of fusion at melting point 11.1 Btu=lb 25.1 kJ=kg
Specific heat of the gas at 70 F (21.1 C)
and 1 atm
CP 0.249 Btu=(lb)( F) 1.04 kJ=(kg)( C)
CV 0.177 Btu=(lb)( F) 0.741 kJ=(kg)( C)
Ratio of specific heats, CP=CV 1.41 1.41
Solubility in water vol=vol at 32 F (0 C) 0.023 0.023
Weight of liquid at boiling point 6.747 lb=gal 808.5 kg=m3
Gas=liquid ratio (gas at 70 F (21.1 C) and 696.5 696.5
1 atm, liquid at boiling point, vol=vol)
Critical point
T ¼ 126.20 K (227.16 R)
P ¼ 33.55 atm (493.0 psia;3.399 MPa)
r ¼ 0.01121 mol=cm3 (19.60 lbm=ft3; 314 kg=m3)
Normal boiling point
T ¼ 77.347 K (139.225 R)
P ¼ 1 atm (14.696 psia; 101.3 kPa)
rgas ¼ 0.0001647 g mol=cm3 (0.2880 lbm=ft3; 4.614 kg=m3)
rliquid ¼ 0.02887 g mol=cm3 (50.49 lbm=ft3; 808.9 kg=m3)
Triple point
T ¼ 63.148 K (113.666 R)
P ¼ 0.1237 atm (1.818 psia; 12.53 kPa)
rgas ¼ 2.41 105 g mol=cm3 (4.21 102 lbm=ft3; 674 g=m3)
rliquid ¼ 0.03098 gmol=cm3 (54.18 lbm=ft3; 867.96 kg=m3)
120 Chapter 3
reverse problem at many plants; they are looking for new oxygen uses to balance the
demand for nitrogen.
Example 3–2 Nitrogen Sample Problem. Nitrogen gas at 138 K and 220 atm
pressure is expanded through a throttling valve. What is the final temperature?
What is the fraction liquefied?
This problem statement is typical of the conditions found in a Joule–Thomson
nitrogen liquefier. The nitrogen gas has been compressed to 200 atm and precooled
to 138 K in a countercurrent heat exchanger just prior to expansion through the J–T
122 Chapter 3
valve. We know the starting conditions, 138 K, 200 atm. The thermodynamic process
is expansion through a J–T valve, which is isenthalpic.
The nitrogen T–S diagram of Fig. 3.54 contains all the information we need to
solve this problem.
The starting point is 138 K, 200 atm, which also fixes the enthalpy at 160 J=gK.
Since the process is isenthalpic, we follow the 160 J=gK isenthalp to the final pressure
of 1 atm. Note that temperature and pressure lines are coincident and horizontal in
the two-phase liquid–vapor region. The final temperature is easily found, namely the
NBP of nitrogen, 77 K. The fraction liquefied is given by the lines of constant quality
(X) in the two-phase region. Our final quality falls between X ¼ 60% vapor and
X ¼ 70% vapor. We estimate the final quality to be about 63% vapor.
124 Chapter 3
63.14a
64.73 0.8622
77.31 0.8084
77.32b
77.5
78.00 0.8043
90.58 0.7433
99.36 0.6922
111.89 0.6071
119.44 0.5332
125.01 0.4314
125.96c 0.31096
3.3. Air
Dry air (avg mol wt. approximately 29.00) is a mixture consisting principally of
nitrogen, oxygen, and argon with traces of other gases as shown in Table 3.25.
When air is liquefied, the carbon dioxide is usually removed; thus, for practical
purposes, liquid air can be considered to consist of 78% nitrogen, 21% oxygen, and
1% argon, the other constituents being present in negligible amounts. Sometimes the
presence of argon is ignored and liquid air is considered to be a binary mixture of
21% oxygen and 79% nitrogen. Since argon has a vapor pressure between those
of oxygen and nitrogen, this assumption is a rather good approximation for some
purposes.
Liquid air has a density of 874 g=L at its normal boiling point of 78.9 K.
Figure 3.55 is a T–S diagram for air.
63.9 0.292
64.3 0.290
64.8 0.284
69.1 0.231
69.25 0.228
71.4 0.209
76.1 0.165
77.33 0.158
111.7 0.074
1 1.778
2 1.780
5 1.784
10 1.793
20 1.809
30 1.827
40 1.846
50 1.867
60 1.890
70 1.916
80 1.946a
90 1.978a
100 2.015a
a
Extrapolated.
Source: Kestin and Wang (1958).
126 Chapter 3
b. Vapor Pressure. It should be noted that there are separate curves for the
vapor pressure of saturated liquid (bubble) and saturated vapor (dew). Air is a multi-
component mixture, and at a given pressure the saturated liquid is at a different tem-
perature than the saturated vapor (see Table 3.28 and Figure 3.56).
Figure 3.36 Vapor pressure of solid and liquid nitrogen (entire pressure range).
Bubble curve:
290:70
log P ¼ 3:5713 þ 0:001494T ð3:10Þ
T
128 Chapter 3
Figure 3.37 Vapor pressure of solid and liquid nitrogen (low pressure range).
Dew curve:
333:88
log P ¼ 4:0816 ð3:11Þ
T
where P is in atm and T is kelvin.
130 Chapter 3
Figure 3.40 Density of gaseous nitrogen (saturated vapor) to the critical point.
132 Chapter 3
Figure 3.41 Density of gaseous nitrogen (saturated vapor) triple point to the NBP.
phase (xi) and the vapor phase (yi) at equilibrium is given by the equilibrium
constant (Ki) :
Ki ¼ yi =xi ð3:12Þ
If both the liquid and vapor phases are ideal mixtures, it can be shown from
Raoult’s and Dalton’s laws that
Ki ¼ Pi =P ð3:13Þ
where Pi is the vapor pressure of component i and P are total system pressure.
134 Chapter 3
For real systems, values of K are determined experimentally and are often
tabulated as the distribution coefficient function, ln (Kp=po), where po is some refer-
ence pressure, such as 101.3 kPa (1 atm). Values for the distribution coefficient func-
tion for nitrogen, oxygen, and argon are given in Table 3.31.
136 Chapter 3
138 Chapter 3
140 Chapter 3
Figure 3.49 Dielectric constant for solid and liquid nitrogen at 1 atm.
142 Chapter 3
144 Chapter 3
has led to the gradual disuse of air. As it is essentially a mixture of nitrogen and oxy-
gen, air does not have the advantage of an invariant boiling temperature. A mixture
containing 80% nitrogen begins boiling at about 79 K, but since the initial vapor has
a composition richer than 80% nitrogen, the remaining liquid gradually becomes
richer in oxygen as boiling proceeds. This increases both the boiling temperature
of the remaining liquid and the hazards associated with oxygen-rich liquid.
This has caused explosions in the vacuum pumps used to reduce the pressure above
liquid air.
Most liquid air is now produced as an intermediate step in the production of
oxygen and nitrogen by distillation. The principal interest in liquid air is in the
preparation of pure nitrogen, oxygen, and rare gases.
3.4. Argon
Argon (at wt. 39.948) has three stable isotopes of mass numbers 36, 38, and 40,
which occur in a relative abundance in the atmosphere of 338:63:100,000.
Liquid argon is a clear, colorless fluid with properties similar to those of liquid
nitrogen. At 1 bar pressure, liquid argon boils at 87.3 K (157.1 R) and freezes
at 83.8 K (151.0 R). Saturated liquid argon at 0.987 bar (1 atm) is more dense than
oxygen, as one would expect, since argon has a larger molecular weight than oxygen
[argon density ¼ 1403 kg=m3 (87.5 lbm=ft3) for saturated liquid at 0.987 bar (1 atm) ].
See Fig. 3.57 for a T–S diagram of argon.
Argon is present in atmospheric air in a concentration of 0.934% by volume
or 1.25% by weight. The atmosphere of Mars contains 1.6% 40Ar and 5 ppm 36Ar,
for comparison. Since the boiling point of argon lies between that of liquid oxygen
and that of liquid nitrogen (slightly closer to that of liquid oxygen), a crude grade
of argon (90–95% pure) can be obtained by adding a small auxiliary argon recovery
column in an air separation plant. About 900 kg (1 ton) of argon is recovered for
every 36 103 kg (40 tons) of oxygen. This is a little less than 50% of the argon
passing through the oxygen column and represents the point of present maximum
146 Chapter 3
Air
Synonyms: Compressed air, atmospheric air, the atmosphere (of the earth)
CAS Registry Number: None (for nitrogen, 7727–37-9; for oxygen, 7782–44-7)
DOT Classification: Nonflammable gas
DOT Label: Nonflammable gas
Transport Canada Classification: 2.2
UN Number: UN 1002 (compressed gas); UN 1003 (refrigerated liquid)
Physical constants
Chemical name Air Air
Molecular weight 28.975 28.975
Density of the gas at 70 F (21.1 C) 0.07493 lb=ft3 1.2000 kg=m3
and 1 atm
Specific gravity of the gas at 70 F (21.1 C) 1.00 1.00
and 1 atm (air ¼ 1)
Specific volume of the gas at 70 F (21.1 C) 13.346 ft3=lb 0.8333 m3=kg
and 1 atm
Boiling point at 1 atm 317.8 F 194.3 C
Freezing point at 1 atm 357.2 F 216.2 C
Critical temperature 221.1 F 140.6 C
Critical pressure 547 psia 3771 kPa abs
Critical density 21.9 lb=ft3 351 kg=m3
Latent heat of vaporization at 88.2 Btu=lb 205 kJ=kg
normal boiling
point
Specific heat of gas at 70 F (21.1 C)
and 1 atm
CP 0.241 Btu=(lb)( F) 1.01 kJ=(kg)( C)
CV 0.172 Btu=(lb)( F) 0.720 kJ=(kg)( C)
Ratio of specific heats (CP=CV) 1.40 1.40
Solubility in water, vol=vol at 32 F (0 C) 0.0292 0.0292
Weight of liquid at normal boiling point 7.29 lb=gal 874 kg=m3
Density of liquid at boiling point 54.56 lb=ft3 874.0 kg=m3
and 1 atm
Gas=liquid ratio (liquid at boiling point, 728.1 728.1
gas at 70 F and 1 atm), vol=vol
Thermal conductivity
at 148 F (100 C) 0.0095 Btu=(hr) 0.0164 W=(m)( C)
(ft)( F=ft)
at 32 F (0 C) 0.0140 Btu=(hr) 0.0242 W=(m)( C)
(ft) ( F=ft)
at 212 F (100 C) 0.0183 Btu=(hr) 0.0317 W=(m)( C)
(ft) ( F=ft)
economic recovery. Some oxygen plants do not recover argon, and US argon
production is thus about a quarter of the theoretical maximum.
The first commercial production, amounting to about 227 m3 (8000 ft3), or
363 kg (800 lb), was carried out by the American Cyanamid Co., Niagara Falls,
Ontario, Canada, in 1914–1915. The argon was produced, not by the air liquefaction
Nitrogen 28 78.09
Oxygen 32 20.95
Argon 40 0.93
Carbon dioxide 44 0.02–0.04
Neon 20.2 18 10-4
Helium 4 5.3 104
Krypton 83.7 1.1 104
Hydrogen 2 0.5 104
Xenon 131.3 0.08 104
Ozone 48 0.02 104
Radon 222 7 1018
148 Chapter 3
route that is now universal, but by chemical reaction of the other constituents of the
air with solids. Production of argon by liquefaction of air was begun by the Amer-
ican Linde Co. (now Linde Division of Union Carbide Corp.) in 1915. Production
was small until about 1943. Before this time, argon was mixed with nitrogen and
used chiefly for filling incandescent lamps; later, similar mixtures were used in fluor-
escent lamps. After 1943, demand for argon as a shield gas in arc welding grew to
exceed all other uses.
Table 3.27 The Heat of Vaporization of Oxygen–Nitrogen Mixtures at 101.3 kPa (1 atm)
Note: There is a maximum at 82% oxygen. Argon is usually considered part of the nitrogen. The heat of a
21% oxygen mixture is 49.04 cal=g.
K 83.8 87.3
kJ=kg 27.78 161.6
150 Chapter 3
is constant within the accuracy of the data. e is the dielectric constant and r is the
density.
For liquid argon, values for the dielectric constant are given in Table 3.33
For gaseous argon at STP, according to Jelatis (1948),
e ¼ 1:000554
d. Density. The data tabulated below are from the work of Baly and Donnan
(1902). These results may be expressed by the formula
Density
3
Pressure (atm) Temp. (K) Volume (cm =mol) g=cm3 lb=ft3
Temp. (K) 101.3 kPa 202.6 kPa 506.6 kPa 101.3 kPa 202.6 kPa 506.6 kPa 101.3 kPa 202.6 kPa 506.6 kPa
Chapter 3
© 2005 by Marcel Dekker
5367-4 Flynn Ch03 R2 090804
84.0 1.4233
84.5 1.4201
85.0 1.4169
85.5 1.4136
86.0 1.4104
86.5 1.4072
87.0 1.4039
87.5 1.4007
88.0 1.3975
88.5 1.3942
89.0 1.3910
89.5 1.3878
90.0 1.3845
e. Other Properties. Figures 3.58–3.66 give the density, specific heats, melting
curve, surface tension, and viscosity of argon, taken from the WADD report men-
tioned before (Johnson, 1960)
154 Chapter 3
Argon
Chemical Symbol: Ar
Synonym: LAR (liquid only)
CAS Registry Number: 7440–37–1
DOT Classification: Nonflammable gas
DOT Label: Nonflammable gas
Transport Canada Classification: 2.2
UN Number: UN 1006 (compressed gas); UN 1951 (refrigerated liquid)
US Units SI Units
Physical constants
Chemical formula Ar Ar
Molecular weight 39.95 39.95
Density of the gas at 70 F (21.1 C) 0.103 lb=ft3 1.650 kg=m3
and 1 atm
Specific gravity of the gas at 70 F 1.38 1.38
(21.1 C) and 1 atm
Specific volume of the gas at 70 F 9.71 ft3=lb 0.606 m3=kg
(21.1 C) and 1 atm
Density of the liquid at boiling point 87.02 lb=ft3 1394 kg=m3
and 1 atm
Boiling point at 1 atm 302.6 F 185.9 C
Melting point at 1 atm 308.6 F 189.2 C
Critical temperature 188.1 F 122.3 C
Critical pressure 711.5 psia 4905 kPa abs
Critical density 33.44 lb=ft3 535.6 kg=m3
Triple point 308.8 F at 9.99 psia 199.3 C at 68.9 kPa abs
Latent heat of vaporization at boiling 69.8 Btu=lb 162.3 kJ=kg
point and 1 atm
Latent heat of fusion at triple point 12.8 Btu=lb 29.6 kJ=kg
Specific heat of the gas at 70 F (21.1 C)
and 1 atm
CP 0.125 Btu=(lb)( F) 0.523 kJ=(kg)( C)
CV 0.075 Btu=(lb)( F) 0.314 kJ=(kg)( C)
Ratio of specific heats 1.67 1.67
Solubility in water at 32 F (0 C) vol=vol 0.056 0.056
Weight of the liquid at boiling point 11.63 lb=gal 1394.0 kg=m3
156 Chapter 3
f. Other Uses. About 1% of the demand for argon results from its use in
lamps. Infrared and photoflood lamps are filled with a mixture of about 88% argon
and 12% nitrogen. The high density and low reactivity of argon help retard the rate
of evaporation of the tungsten filament, thus prolonging the life of the lamp; the
nitrogen prevents arcing.
3.5. Neon
Neon (at wt. 20.183) has three stable isotopes of mass numbers 20, 21, and 22, which
occur in atmosphere air in the relative abundance 10,000:28:971. Five unstable iso-
topes are also known.
Neon is a rare gas element present in the atmosphere to the extent of 1 part in
65,000 of air. It is obtained by liquefaction of air and seperated from the other gases
by fractional distillation. It is a very insert element; however, it is said to form a com-
pound with fluorine. Of all the rare gases, the discharge of neon is the most intense at
ordinary voltages are currents.
Temperature entropy diagrams for neon are shown in Fig 3.67 and 3.68.
158 Chapter 3
160 Chapter 3
162 Chapter 3
c. Melting Point Temperature of Neon. Mills and Grilley (1955) give the melt-
ing curve of neon as
p ¼ a þ bT c ð3:20Þ
where p is the pressure in kg=cm2, T is the temperature in K, a ¼ 1057.99,
b ¼ 6.289415, c ¼ 1.599916.
d. Vapor Pressure of Solid Neon. Between 15 and 20.5 K, the following
equation may be used:
111:76
log np ¼ þ 6:0424 ð3:21Þ
T
where vp is in cmHg and T is in K.
e. Vapor Pressure of Liquid Neon. Verschaffelt (1929) gives the following
vapor pressure equation for liquid neon between 25 and 35 K:
T log np ¼28:100 þ 36:00ðT 35Þ 0:003333ðT 35Þ2
ð3:22Þ
þ 0:000400ðT 35Þ3 þ 0:00002667ðT 35Þ4
where vp is in atm and T is in K on the Leiden scale (0 C ¼ 273.09 K).
f. Density. The density of solid neon (Johnson, 1960) is in the range:
4.3 1443
24.57a 1444
a
Normal melting point value.
164 Chapter 3
Chemical Symbol: Kr Ne Xe
CAS Registry Number: 7439-90-9 7440-01-9 7440-63-3
DOT Classification: Nonflammable gas Nonflammable gas Nonflammable gas
DOT Label: Nonflammable gas Nonflammable gas Nonflammable gas
Transport Canada 2.2 2.2 2.2
Classification:
UN Number:
(compressed gas): UN 1056 UN 1065 UN 2036
(refrigerated liquid): UN 1970 UN 1913 UN 2591
US Units SI Units
Krypton
Chemical formula Kr Kr
Molecular weight 83.80 83.80
Density of the gas at 70 F 0.2172 lb=ft3 3.479 kg=m3
(21.1 C) and 1 atm
Specific gravity of the gas at 2.899 2.899
70 F (21.1 C) and 1 atm
(air ¼ 1)
Specific volume of the gas at 4.604 ft3=lb 0.287 m3=kg
70 F (21.1 C) and 1 atm
Boiling point at 1 atm 244.0 F 153.4 C
Melting point at 1 atm 251 F 157 C
Critical temperature 82.8 F 63.8 C
Critical pressure 798.0 psia 5502 kPa abs
Critical density 56.7 lb=ft3 908 kg=m3
Triple point 251.3 F at 10.6 psia 157.4 C at
73.2 kPa abs
Latent heat of vaporization 46.2 Btu=lb 107.5 kJ=kg
at boiling point
Latent heat of fusion at 8.41 Btu=lb 19.57 kJ=kg
triple point
Specific heat of the gas at 70 F
(21.1 C) and 1 atm
CP 0.060 Btu=(lb)( F) 0.251 kJ=(kg)( C)
CV 0.035 Btu=(lb)( F) 0.146 kJ=(kg)( C)
Ratio of specific heats, CP=CV 1.69 1.69
Solubility in water, vol=vol at 0.0594 0.0594
68 F (20 C)
Weight of the liquid at boiling point 20.15 lb=gal 2415 kg=m3
Density of the liquid at boiling point 150.6 lb=ft3 2412.38 kg=m3
Gas=liquid volume ratio (liquid 693.4 693.4
at boiling point, gas at 70 F
and 1 atm, vol=vol)
Neon
Physical Constants
Chemical formula Ne Ne
Molecular weight 20.183 20.183
Density of the gas at 70 F (21.1 C) 0.05215 lb=ft3 0.83536 kg=m3
and 1 atm
Specific gravity of the gas at 70 F (21.1 C) 0.696 0.696
and 1 atm (air ¼ 1)
(Continued)
166 Chapter 3
US Units SI Units
Specific volume of the gas at 70 F (21.1 C) 19.18 ft3=lb 1.197 m3=kg
and 1 atm
Boiling point at 1 atm 410.9 F 246.0 C
Melting point at 1 atm 415.6 F 248.7 C
Critical temperature 379.8 F 228.8 C
Critical pressure 384.9 psia 2654 kPa abs
Critical density 30.15 lb=ft3 483 kg=m3
Triple point 415.4 F at 6.29 psia 248.6 C at 43.4 kPa abs
Latent heat of vaporization at boiling point 37.08 Btu=lb 86.3 kJ=kg
Latent heat of fusion at triple point 7.14 Btu=lb 16.6 kJ=kg
Specific heat of the gas at 70 F (21.1 C)
and 1 atm
CP 0.25 Btu=(lb)( F) 1.05 kJ=(kg)( C)
CV 0.152 Btu=(lb)( F) 0.636 kJ=(kg)( C)
Ratio of specific heats, CP=CV 1.64 1.64
Solubility in water, vol=vol at 68 F (20 C) 0.0105 0.0105
Weight of the liquid at boiling point 10.07 lb=gal 1207 kg=m3
Density of saturated vapor at 1 atm 0.5862 lb=ft3 9.390 kg=m3
Density of the gas at boiling point 0.6068 lb=ft3 9.7200 kg=m3
Density of the liquid at boiling point 75.35 lb=ft3 1207 kg=m3
Gas=liquid volume ratio (liquid at boiling point, 1445 1445
gas at 70 F and 1 atm, vol=vol)
Xenon
Physical Constants
Chemical formula Xe Xe
Molecular weight 131.3 131.3
Density of the gas at 70 F (21.1 C) and 1 atm 0.3416 lb=ft3 5.472 kg=m3
Specific gravity of the gas at 70 F and 4.560 4.560
1 atm (air ¼ 1)
Specific volume of the gas at 70 F (21.1 C) 2.927 ft3=lb 0.183 m3=kg
and 1 atm
Boiling point at 1 atm 162.6 F 108.2 C
Melting point at 1 atm 168 F 111 C
Critical temperature 61.9 F 16.6 C
Critical pressure 847.0 psia 5840 kPa abs
Critical density 68.67 lb=ft3 1100 kg=m3
Triple point 169.2 F at 11.84 psia 111.8 C at 81.6 kPa abs
Latent heat of vaporization at boiling point 41.4 Btu=lb 96.3 kJ=kg
Latent heat of fusion at triple point 7.57 Btu=lb 17.6 kJ=kg
Specific heat of the gas at 70 F (21.1 C)
and 1 atm
CP 0.038 Btu=(lb)( F) 0.269 kJ=(kg)( C)
CV 0.023 Btu=(lb)( F) 0.096 kJ=(kg)( C)
Ratio of specific heats, CP=CV 1.667 1.667
Solubility in water, vol=vol at 68 F (20 C) 0.108 0.108
Weight of the liquid at boiling point 25.51 lb=gal 3057 kg=m3
Density of the liquid at boiling point 190.8 lb=ft3 3057 kg=m3
Gas=liquid volume ratio (liquid at 558.5 558.5
boiling point, gas at 70 F and 1 atm,
vol=vol)
Pressure Temperature
Property kPa atm psia K R
24.57a 1248
25.17 1238
27.17b 1204
30.13 1149
36.05 1017
39.08 928
43.02 749
44.38c 483
a
Triple point temperature.
b
Normal boiling temperature.
c
Critical temperature.
168 Chapter 3
g. Dielectric Constant of Gaseous Neon. One value has been published at 70 F
and 1 atm. The Clausius–Mosotti function may be used to estimate values at other
conditions with e ¼ 1.0001274.
h. Other Properties. Figures 3.69–3.76 give the density, specific heats, thermal
conductivity, vapor pressure, and melting curve of neon, taken from the WADD
report mentioned before (Johnson 1960).
170 Chapter 3
172 Chapter 3
174 Chapter 3
3.6. Fluorine
Liquid fluorine (at wt. 18.9984) is a light yellow liquid having a normal boiling point
of 85.3 K (153.6 R). At 53.5 K (96.5 R) and 0.987 bar (1 atm), liquid fluorine freezes
as a yellow solid, but upon subcooling to 45.6 K (82 R), it transforms to a white
solid. Liquid fluorine is one of the most dense cryogenic liquids [density at normal
boiling point ¼ 1504.3 kg=m3 (93.9 lbm=ft3)].
Fluorine will react with almost all substances. If it comes in contact with
hydrocarbons, it will react hypergolically with a high heat of reaction, which is some-
times sufficiently high that the metal container for the fluorine is ignited. Satisfactory
handling systems for liquid fluorine have been operated, but this substance presents
even greater difficulty than oxygen because of its higher reactivity. The flow of liquid
fluorine through a system is preceded by a process of ‘‘passivation’’ in which fluorine
gas is admitted slowly to the system in order to build up a layer of passive fluorides
on its surfaces. This layer then presents a chemically inert barrier to the large quan-
tities of fluorine that follow.
Fluorine is highly toxic. The fatal concentration range for animals is
200 ppm hr; i.e., for an exposure time of l hr, 200 ppm of fluorine is fatal; for an expo-
sure time of 15 min, 800 ppm is fatal; and for an exposure time of 4 hr, 50 ppm is
fatal. The maximum allowable concentration for human exposure is usually consid-
ered to be approximately 1 ppm hr. The presence of fluorine in air may be detected by
its sharp, pungent odor for concentrations as low as 1–3 ppm.
Fluorine was first isolated in 1886 by the French chemist Henri Moissan after
nearly 75 years of unsuccessful attempts by several others. For many years after its
isolation, fluorine remained little more than a scientific curiosity, to be handled with
extreme caution because of its toxicity. Commercial production of fluorine began
during World War II when large quantities were required in the fluorination of ura-
nium tetrafluoride (UF4) to produce uranium hexafluoride (UF6) for the isotopic
separation of uranium-235 by gaseous diffusion in the development of the atomic
bomb. Today, commercial production methods are essentially variations of the
Moisson process, and safe techniques have been developed for the bulk handling
of liquid fluorine.
The major use of fluorine is in the form of fluorspar. Although not a major
commodity in terms of total quantity produced, fluorspar is a critical raw material
for the aluminum, chemical, and steel industries of the world. There is at present
no significant use of cryogenic liquid fluorine per se.
Table 3.37 gives the physical constants for fluorine.
176 Chapter 3
Chemical Symbol: F2
CAS Registry Number: 7782-41-4
DOT Classification: Nonflammable gas
DOT Label: Poison and Oxidizer
Transport Canada Classification: 2.3 (5.1)
UN Number: UN 1045
U.S.Units SI Units
Physical constants
Chemical formula F2 F2
Molecular weight 38.00 38.00
Density of the gas
at 32 F (0.0 C) and 1 atm 0.106 lb=ft3 1.70 kg=m3
at 70 F (21.1 C) and 1 atm 0.098 lb=ft3 1.57 kg=m3
Specific gravity of the gas
at 70 F (21.1 C) and 1 atm 1.312 1.312
Specific volume of the gas
at 70 F (21.1 C) and 1 atm 10.17 0.635
Density of the liquid
at 306.8 F (188.2 C) 94.2 lb=ft3 1509 kg=m3
at 320.4 F (195.8 C) 97.9 lb=ft3 1568 kg=m3
Boiling point at 1 atm 306.8 F 188.2 C
Melting point at 1 atm 363.4 F 219.6 C
Critical temperature 199.9 F 128.8 C
Critical pressure 756.4 psia 5215 kPa abs
Critical density 35.8 lb=ft3 573.6 kg=m3
Triple point 363.4 F at 219.6 C at
0.324 psia 0.223 kPa abs
Latent heat of vaporization
at 306.8 F (188.2 C) 74.8 Btu=lb 173.5 kJ=kg
Latent heat of fusion
at 363.4 F (219.6 C) 5.8 Btu=lb 13.4 kJ=kg
Weight of liquid per gallon
at 306.8 F (188.1 C) 12.6 lb=gal 1509.8 kg=m3
at 320.4 F (195.8 C) 13.1 lb=gal 1569.7 kg=m3
Heat capacity of the gas, CP
at 32 F (0 C) 0.198 Btu=(lb)( F) 0.828 kJ=(kg)( C)
at 70 F (21.1 C) 0.197 Btu=(lb)( F) 0.825 kJ=(kg)( C)
Heat capacity of the gas, CV
at 70 F (21.1 C) 0.146 Btu=(lb)( F) 0.610 kJ=(kg)( C)
Ratio of specific heats, CP=CV
at 70 F (21.1 C) 1.353 1.353
Thermal conductivity of gas
at 32 F (0 C) and 1 atm 0.172 Btu.in.=hr F 0.248 W=m.K
Viscosity
of gas at 32 F (0 C) and 1 atm 0.0218 cP 0.0218 mPa s
of liquid at 314 F (192.2 C) 0.275 cP 0.275 mPa s
85.25 166.2
90 160.9
95 154.9
100 148.6
105 142.0
110 134.6
115 126.6
120 117.7
125 107.6
130 95.65
135 80.81
140 59.12
144 14.85
178 Chapter 3
90.0 0.767
169.3 1.424
200.0 1.680
289.1 2.345
327.1 2.547
3.6.5. Viscosity
Tables 3.40 and 3.41 present data on the viscosity of gaseous and liquid fluorine,
respectively.
100 0.0862
150 0.134
200 0.183
250 0.228
273 0.247
300 0.269
350 0.308
180 Chapter 3
3.7. Hydrogen
Natural hydrogen (at wt. 1.008) is a mixture of two stable isotopes; hydrogen (1H),
of atomic mass 1, and deuterium (2H), of atomic mass 2. The abundance ratio is
about 6400:1, although slight variations occur, depending on the source of the
hydrogen. Actually, since molecular hydrogen is diatomic, nearly all the deuterium
atoms in natural hydrogen exist in combination with hydrogen atoms. Accordingly,
ordinary hydrogen is a mixture of H2 and HD (1H2H) molecules in the ratio 3200:1.
A very rare radioactive isotope of hydrogen of atomic mass 3 (3H), called tritium,
also exists. Tritium has a half-life of about 12.5 years.
Hydrogen is the most abundant of all the elements in the universe, and it is
thought that the heavier elements were, and still are, being built from hydrogen
and helium. It has been estimated that hydrogen makes up more than 90% of all
the atoms or three-fourths of the mass of the universe. Production of hydrogen in
the United States alone now amounts to 85 million m3 (3 billion ft3). The physical
constants for hydrogen are shown in Table 3.43.
182 Chapter 3
184 Chapter 3
Hydrogen
Chemical Symbol: H2
CAS Registry Number: 1333-74-0
DOT Classification: Flammable gas
DOT Label: Flammable gas
Transport Canada Classification: 2.1
UN Number: UN 1049 (compressed gas); UN 1966 (refrigerated liquid)
Normal Hydrogen
Chemical formula H2 H2
Molecular weight 2.016 2.016
Density of the gas at 70 F (21.1 C) 0.00521 lb=ft3 0.08342 kg=m3
and 1 atm
Specific gravity of the gas at 32 F 0.06960 0.06960
(0 C) and 1 atm (air ¼ 1)
Specific volume of the gas at 70 F 192.0 ft3=lb 11.99 m3=kg
(21.1 C) and 1 atm
Boiling point at 1 atm 423.0 F 252.8 C
Melting point at 1 atm 434.55 F 259.2 C
Critical temperature 399.93 F 239.96 C
Critical pressure 190.8 psia 1315 kPa abs
Critical density 1.88 lb=ft3 30.12 kg=m3
Triple point 434.55 F at 259.2 C at
1.045 psia 7.205 kPa abs
Latent heat of vaporization at 191.7 Btu=lb 446.0 kJ=kg
boiling point
Latent heat of fusion at triple point 24.97 Btu=lb 58.09 kJ=kg
Specific heat of the gas at 70 F
(21.1 C) and 1 atm
CP 3.425 Btu=(lb)( F) 14.34 kJ=(kg)( C)
CV 2.418 Btu=(lb)( F) 10.12 kJ=(kg)( C)
Ratio of specific heats
(CP=CV) 1.42 1.42
Solubility in water, vol=vol at 60 F 0.019 0.019
(15.6 C)
Density of the gas at boiling point 0.083 lb=ft3 1.331 kg=m3
and 1 atm
Density of the liquid at boiling point 4.43 lb=ft3 (0.5922 lb=gal) 70.96 kg=m3
and 1 atm
Gas=liquid ratio (liquid at boiling point, 850.3 850.3
gas at 70 F (21.1 C) and 1 atm), vol=vol
Heat of combustion at 70 F (21.1 C)
and 1 atm
Gross 318.1 Btu=ft3 11,852 kJ=m3
Net 268.6 Btu=ft3 10,009 kJ=m3
Para Hydrogen
(The following constants are different from normal hydrogen)
Boiling point at 1 atm 423.2 F 252.9 C
Melting point 434.8 F 259.3 C
(Continued)
186 Chapter 3
metric tons=day), Mama Bear (3.5 US tons=day or 3.2 metric tons=day), and Papa
Bear (30 tons=day or 27 metric tons=day) were designed and built by Air Products
and Chemicals, Inc. Papa Bear, which went on-stream in West Palm Beach, Florida,
in February 1959, was the world’s first truly large-scale liquid hydrogen plant.
From 1960 to the present, several additional large-scale commercial liquid
hydrogen plants were built in North America. For example, in 1964 the Linde Divi-
sion of Union Carbide built the largest single plant to date, a 60 US ton=day (54
metric tons=day) plant in Sacramento, California. Between 1965 and 1976, Air Pro-
ducts constructed a 66 US tons=day (60 metric tons=day) two-plant liquid hydrogen
production complex in New Orleans, Louisiana, and between 1982 and 1990 three
liquid hydrogen plants were built in Canada by Air Products, Liquid Air, and BOC.
Today, nine large-scale liquid hydrogen plants are in operation in the United
States and Canada with a combined capacity of about 190 US tons=day (172 metric
tons=day). Some older plants like Papa Bear and the Linde Sacramento plant no
longer exist. Table 3.44 lists some key milestones in the creation of large-scale com-
mercial liquid hydrogen plants, and Fig. 3.83 shows the location, size, and ownership
of liquid hydrogen plants in North America in 1990. The aerospace program, which
created the initial demand for large quantities of liquid hydrogen, accounted for only
20% of liquid hydrogen demand in 1990. A variety of commercial applications
account for the remaining demand.
The diversity of commercial liquid hydrogen containers for both stationary sto-
rage and transportation are listed below. Three insulation technologies are
employed. In increasing order of thermal efficiency and cost, they are evacuated per-
lite insulation, evacuated multilayer insulation, and liquid nitrogen shielded multi-
layer insulation. Large stationary storage containers are reasonably efficient, using
just an evacuated annulus filled with perlite, due to their relatively low surface-to-
volume ratio. Smaller stationary containers and all transportation containers require
multilayer insulation for acceptable efficiency (i.e., low vent rates). Liquid nitrogen
shielded multilayer insulation is widely used for liquid helium storage and transport,
Date Event
but for liquid hydrogen it is typically limited to portable containers that must be
shipped long distances (often overseas) without venting for several days.
Boiloff from product liquid hydrogen storage tanks at the liquid hydrogen
plant site is usually redirected back to the liquefier plant for reliquefaction.
Transportation
Over-the-road tankers 13,000–15,000
Seagoing containers 11,000
Inland waterway barges 250,000
Railcars 28,000–34,000
Storage
Large tanks
Spheres at Kennedy Space Center 830,000
Spheres at Air Products’ New Orleans plant 500,000
Horizontal cylinders at Air Products’
Sacramento plant 70,000
Customer tanks
Horizontal=vertical cylinder 1500–9000
Horizontal cylinder 9000–20,000
a
Insulation types include perlite, multilayer, and liquid nitrogen shielded multilayer.
b
1000 US gal ¼ 3.785 L.
188 Chapter 3
Figure 3.83 Liquid hydrogen plants in the United States and Canada, 1990.
Figure 3.84 illustrates the two states of the hydrogen molecule: the higher
energy ortho state in which the protons in the two nuclei spin in the same direction,
and the lower energy para state in which the protons spin in opposite directions.
Temperature determines the relative abundance of each state, as shown in Table 3.45
and by the ortho–para thermodynamic equilibrium diagram in Fig. 3.85. At ambient
and higher temperatures, hydrogen consists of 25% para and 75% ortho and is
referred to as ‘‘normal’’ hydrogen (n-H2). At the normal boiling point of liquid
hydrogen (20.4 K), equilibrium shifts to almost 100% para. Conversion from ortho
to para is exothermic, with the heat of conversion of 302 Btu=lb (702 J=g) substan-
tially exceeding the latent heat of vaporization of normal boiling point hydrogen
of 195 Btu=lb (454 J=g).
The equilibrium mixture of o-H2 and p-H2 at any given temperature is called
equilibrium hydrogen (e-H2). At the normal boiling point of hydrogen, 20.4 K
(36.7 R), equilibrium hydrogen has a composition of 0.21% o-H2 and 99.79%
p-H2, or practically all parahydrogen.
Two deuterium atoms can combine to form a deuterium molecule, often
denoted D2 and sometimes called heavy hydrogen. As in the case of ‘‘light’’ hydro-
gen, the nuclear spins of the atoms in a deuterium molecule can assume different
Figure 3.84 Ortho- and parahydrogen differ in the manner in which their two protons spin
in the nucleus of the hydrogen atom.
spatial orientations with respect to each other, although the orientations are a bit
more complicated. Ortho- and para deuterium have different sets of possible spatial
orientations, and the ortho–para equilibrium concentrations are temperature-
dependent, as in the case of hydrogen. The nucleus of the deuterium atom consists
190 Chapter 3
of one proton and one neutron, so that the high-temperature composition (compo-
sition of normal deuterium) is two-thirds orthodeuterium and one-third paradeuter-
ium. In the case of deuterium, p-D2 converts to o-D2 as the temperature is
decreased, in contrast to hydrogen, in which o-H2 converts to p-H2 upon a decrease
in temperature. (The prefix ortho- simply means common and does not inherently
refer to the more energetic form.) The hydrogen deuteride molecule does not have
the symmetry that hydrogen and deuterium molecules possess; therefore, HD exists
in only one form.
Temperature–entropy diagrams for parahydrogen and deuterium are presented
in Figs. 3.86–3.91
Table 3.46 lists the heat of conversion for normal hydrogen to parahydrogen.
Figure 3.92 presents similar data graphically.
Table 3.47 gives the fixed points for hydrogen, and Tables 3.48 and 3.49 sum-
marize the fixed point properties of normal hydrogen and parahydrogen.
A catalyst is normally required to convert between the ortho and para forms of
hydrogen. Without a catalyst the ortho to para reaction proceeds only very slowly in
the liquid phase and not at all in the vapor phase. In the liquid phase, unconverted
normal liquid hydrogen will slowly convert to p-H2, releasing sufficient heat to
vaporize most of the liquid. Thus, ortho-to-para conversion must be incorporated
into the hydrogen liquefier so that the product liquid will consist almost entirely
of stable para-hydrogen.
To minimize the impact of heat removal due to ortho-to-para conversion on
the thermodynamic efficiency of the liquefier cycle, it is important to accomplish
ortho-to-para conversion at the highest possible temperature levels consistent with
thermodynamic equilibrium. Continuous conversion would be the most thermody-
namically efficient. In continuous conversion, ortho-to-para conversion occurs
simultaneously with cooling so that the hydrogen cooling path closely follows the
equilibrium curve. (Some offset from full equilibrium is necessary to provide some
driving force for conversion). Continuous conversion, however, requires a complex
integration of catalyst and heat exchange surface that is difficult to design into prac-
tical systems. A common design approach is to approximate continuous conversion
by using several simple adiabatic catalyst converters interdispersed with heat
exchange. Starting at around liquid nitrogen temperatures, hydrogen is partially
cooled, removed from the heat exchanger, and sent through an adiabatic converter
where the hydrogen partially rewarms due to the heat of conversion. Conversion to
parahydrogen is accomplished to the extent permitted by approaching near equili-
brium at the converter outlet temperature. The partially converted hydrogen is
returned to the heat exchanger, and the stepwise heat exchange-conversion process
is repeated until the product hydrogen from the final converter is at least 95% para-
hydrogen.
Two types of ortho–para catalysts are commercially available: an iron oxide
catalyst and a nickel silicate catalyst. The nickel silicate catalyst is the more active
of the two. Anywhere from 50 to a few hundred pounds of catalyst would be used
for each ton per day of liquefier plant capacity.
Most of the physical properties of hydrogen and deuterium, such as vapor
pressure, density of the liquid, and triple point temperature and pressure, are mildly
dependent upon the ortho–para composition. The transition from one form to the
other involves a substantial energy change, as shown earlier in Table 3.45.
Unfortunately, the uncatalyzed conversion of orthohydrogen to parahydrogen
is relatively slow, which leads to an important consideration in the liquefaction and
storage of hydrogen that does not exist for other gases. As mentioned, equilibrium
hydrogen at room temperature consists of 75% orthohydrogen and 25% parahydro-
gen. At 20.39 K, however, the approximate temperature at which hydrogen boils at
atmospheric pressure, the equilibrium concentration, is 99.8% parahydrogen. If one
starts with equilibrium or normal hydrogen at room temperature and liquefies it, the
liquid orthohydrogen will slowly convert to parahydrogen. Since parahydrogen is
192 Chapter 3
the lower energy form, heat is liberated in the conversion process, and this is suffi-
cient to cause evaporation of nearly 70% of the hydrogen originally liquefied. Lique-
faction cycles, therefore, usually incorporate a catalyst through which the newly
liquefied hydrogen must pass. This causes the conversion to take place in the lique-
fier, where the heat can be absorbed conveniently, and eliminates the otherwise
gradual evaporation of the stored liquid from this cause.
194 Chapter 3
dq dx
¼ Hm ¼ Hkmx2 ð3:27Þ
dt dt
where H is heat generated when converting unit mass from ortho- to parahydrogen.
196 Chapter 3
dm 1 dq Hk
¼ ¼ mx2 ð3:28Þ
dt L dt L
where L is the latent heat of vaporization per unit mass.
Heat of conversion
10 527.139 253.9865
20 527.140 253.987
20.39 527.138 253.986
30 527.138 253.986
33.1 527.138 253.986
40 527.117 253.976
50 526.845 253.845
60 525.531 253.212
70 521.770 251.400
80 513.932 247.623
90 500.757 241.275
100 481.671 232.079
120 427.248 205.857
150 322.495 155.385
200 163.774 78.91
250 70.524 33.98
298.16 28.558 13.76
300 27.562 13.28
198 Chapter 3
Critical point
T ¼ 32.976 0.05 K (59.357 R) T ¼ 33.19 K (59.742 R)
P ¼ 12.759 atm (1.2928 MPa) P ¼ 12.98 atm (1.315 MPa)
r ¼ 15.59 mol=L (31.43 kg=m3) r ¼ 14.94 mol=L (30.12 kg=m3)
Normal boiling point
T ¼ 20.268 K (36.482 R) T ¼ 20.39 K (36.702 R)
>P ¼ 1 atm (0.101325 MPa) P ¼ 1 atm (0.101325 MPa)
rliq ¼ 35.11 mol=L (70.78 kg=m3) rliq ¼ 35.2 mol=L (71.0 kg=m3)
rvap ¼ 0.6636 mol=L (1.338 kg=m3) rvap ¼ 0.6604 mol=L (1.331 kg=m3)
Triple point
T ¼ 13.803 K (24.845 R) T ¼ 13.957 K (25.123 R)
P ¼ 0.0695 atm (0.00704 MPa) P ¼ 0.0711 atm (0.00720 MPa)
rsolid ¼ 42.91 mol=L (86.50 kg=m3) rsolid ¼ 43.01 mol=L (86.71 kg=m3)
rliq ¼ 38.21 mol=L (77.03 kg=m3) rliq ¼ 38.3 mol=L (77.2 kg=m3)
rvap ¼ 0.0623 mol=L (0.126 kg=m3) rvap ¼ 0.0644 mol=L (0.130 kg=m3)
Equation (3.28) gives the rate of variation of m with t. To solve this equation, it
is convenient to first express x as an explicit function of t, i.e., solve Eq. (3.25). The
general solution is 1=x ¼ kt þ A, where A is the constant of integration.
Starting with normal liquid in the dewar, x ¼ 0.75 when t ¼ 0, from which
A ¼ 1.33. then
1
¼ kt þ 1:33
x
or
1
x¼ ð3:29Þ
kt þ 1:33
Substitution of this value for x into Eq. (3.28) gives
!
dm kH m
¼
dt L ðkt þ 1:33Þ2
!
dm kH dt
¼
m L ðkt þ 1:33Þ2
kH 1
ln m ¼ þ In B ð3:30Þ
L kð1:33 þ ktÞ
where ln B is a constant of integration. Now write
m H
ln ¼ ð3:31Þ
B Lð1:33 þ ktÞ
or
m H
¼ exp ð3:31aÞ
B Lð1:33 þ ktÞ
Cryogenic Fluids
Triple point Normal boiling point Standard conditions
Property Solid Liquid Vapor Liquid Vapor Critical pointa STP(0 C) NTP(20 C)
Temperature (K) 13.957b 13.957b 13.957b 20.390b 20.390b 33.19b 273.15 293.15
Pressure (mmHg) 54.04 54.04 54.04 760 760 9865 760 760
Density (mol=cm3) 103 43.01 38.30 0.0644 35.20 0.6604 14.94 0.04460 0.04155
Specific volume (cm3=mol) 103 0.02325 0.026108 15.519 0.028409 1.5143 0.066949 22.423 24.066
Compressibility factor, Z ¼ PV=RT – 0.001621 0.9635 0.01698 0.9051 0.3191 1.00042 1.00049
Heats of fusion and vaporizationc (J=mol) 117.1 911.3 – 899.1 – 0 – –
Specific heat (J=(mol K)
Cs, at saturation 5.73 13.85 –46.94 18.91 –33.28 (Very large) – –
CP, at constant pressure – 13.23 21.22 19.70 24.60 (Very large) 28.59 28.89
CV, at constant volume – 9.53 12.52 11.60 13.2 (19.7) 20.30 20.40
Specific heat ratiod, g ¼ CP=CV – 1.388 1.695 1.698 1.863 (Large) 1.408 1.416
Enthalpye (J=mol) 321.6 438.7 1350.0 548.3 1447.4 1164 7749.2 8324.1
Internal energye (J=mol) 317.9d 435.0 1234.8 545.7 1294.0 – 5477.1 5885.4
Entropye (J=mol k) 20.3 28.7 93.6 34.92 78.94 54.57 139.59 141.62
Velocity of sound (m=s) – 1282 307 1101 357 – 1246 1294
Viscosity, m – 0.026 0.00074 0.0132 0.0011 (0.0035) 0.00839 0.00881
N s=m2 103
centipoise – 0.026 0.00074 0.0132 0.0011 (0.0035) 0.00839 0.00881
Thermal conductivity, k(mW(cm K)) 9.0 0.73 0.124 0.99 0.169 –f 1.740 1.838
Prandtl number, Npr ¼ mCp=k – 2.24 0.623 1.29 0.809 – 0.6849 0.6848
Dielectric constant, e 1.287 1.253 1.00039 1.231 1.0040 1.0937 1.000271 1.000253
pffiffi
Index of refractiong,h, n ¼ e 1.134 1.119 1.000196 1.1093 1.0020 1.0458 1.000136 1.000126
Surface tension (N=m) 103 – 3.00 – 1.94 – 0 – –
Equiv. vol.=vol. liquid at NBP 0.8184 0.9190 546.3 1 53.30 2.357 789.3 847.1
Gas constant: R ¼ 8.31434 J(mol K); molecular weight ¼ 2.01594; mole ¼ gram mole.
a
Values in parentheses are estimates.
b
These temperatures are based on the IPTS-1968 temperature scale. To compare with the corresponding temperatures based on the NBS-1955 temperature scale, 0.01 K should
be subtracted; i.e., use 13.947, 20.380, and 33.18, respectively, for the triple point, normal boiling point, and critical point when determining property differences between
199
normal and parahydrogen.
c
Heats of fusion and vaporization calculated from enthalpy differences.
d
Calculated from property values given in this table.
e
Base point (zero values) for enthalpy, internal energy, and entropy are 0 K for the ideal gas at 0.101325 Mpa (1 atm) pressure.
f
Anomalously large.
g
Long wavelengths.
h
Index of refraction calculated from dielectric constant data.
200
Triple point Normal boiling point Standard conditions
Property Solid Liquid Vapor Liquid Vapor Critical pointa STP(0 C) NTP(20 C)
Temperature (K) 13.803 13.803 13.803 20.268 20.268 32.976 273.15 293.15
Pressure (mmHg) 52.82 52.82 52.85 760 760 9696.8 760 760
Density (mol=cm3) 103 42.91 38.207 0.0623 35.11 0.6636 15.59 0.5459 0.04155
Specific volume (cm3=mol) 103 0.02330 0.026173 16.057 0.028482 1.5069 0.064144 22.425 24.069
Compressibility factor, Z ¼ PV=RT – 0.001606 0.9850 0.01712 0.9061 0.3025 1.0005 1.0006
Heats of fusion and vaporizationb (J=mol) 117.5 905.5 – 898.3 – 0 – –
Specific heat (J=(mol K)
Cs at saturation 5.73 13.85 46.94 18.91 33.28 (Very large) – –
Cp, at constant pressure – 13.13 21.20 19.53 24.50 (Very large) 30.35 30.02
CV, at constant volume – 9.50 12.52 11.57 13.11 19.7 21.87 21.70
Specific heat ratioc g ¼ CP=CV – 1.382 1.693 1.688 1.869 (Large) 1.388 1.383
Enthalpyd (J=mol) 740.2 622.7 282.8 516.6 381.7 77.6 7656.6 8260.6
Internal energyd (J=mol) 740.4c 622.9 169.8 519.5 229.0 5.7 5384.5 5822.0
Entropyd (J=(mol K) 1.49 10.00 75.63 16.08 60.41 35.4 127.77 129.90
Velocity of sound (m=s) – 1273 305 1093 355 350 1246 1294
Viscosity,m – 0.026 0.00074 0.0132 0.0011 0.0035 0.00839 0.00881
N s=m2 103
centipoise – 0.026 0.00074 0.0132 0.0011 0.0035 0.00839 0.00881
Thermal conductivity,k(mW(cm K)) 9.0 0.73 0.124 0.99 0.169 –e 1.841 1.914
Prandtl number, Npr ¼ mCp=k – 2.24 0.623 1.29 0.809 0.6866 0.6855
Dielectric constant, e 1.286 1.252 1.00038 1.230 1.0040 1.098 1.00027 1.00026
pffiffi
Index of refraction f,g, n ¼ e 1.134 1.119 1.00019 1.109 1.0020 1.048 1.00013 1.00012
Surface tension (N=m) 103 – 2.99 – 1.93 – 0 – –
Equiv. vol.=vol. liquid at NBP 0.8181 0.9190 563.8 1 52.91 2.252 787.4 845.1
Chapter 3
Gas constant: R ¼ 8.31434 J(mol K); molecular weight ¼ 2.01594; mole ¼ gram mole.
a
Values in parentheses are estimates.
b
Heats of vaporization calculated from enthalpy differences.
c
Calculated from property values given in this table.
d
Base point (zero values) for enthalpy, internal energy, and entropy are 0 K for the ideal gas at 0.101325 MPa (1 atm) pressure.
e
Anomalously large.
f
Long wavelengths.
g
Index of refraction calculated from dielectric constant data.
when t ¼ 0, m ¼ mo, so
mo 3H mo
¼ exp or B¼ ð3:32Þ
B 4L expð3H=4LÞ
The solution then becomes
m 3H H
exp ¼ exp ð3:33Þ
mo 4L Lð1:33 þ ktÞ
or
m H 3H
¼ exp ð3:33aÞ
mo Lð1:33 þ ktÞ 4L
L ¼ 906.09 J=mol (216.56 cal=mol) for n-H2 at 20 K (36 R) and 894.87 J=mol
(213.88 cal=mol) for 99.8% p-H2 at 20 K. 890 J=mol (215 cal=mol) is a reasonable
average value. The value for H is 1416.19 J=mol (338.65 cal=mol), and the value
for the rate constant k is 0.0114=hr1.
These numerical substitutions reduce Eq. (3.33) to the form
m 1:57
¼ exp 1:18 ð3:34Þ
mo ð1:33 þ 0:0114tÞ
or
m 1:57
ln ¼ 1:18 ð3:34aÞ
mo 1:33 þ 0:0114t
If the original composition of the liquid is not 0.75 orthohydrogen at t ¼ 0, then
a new constant of integration based on the starting composition can be evaluated for
Eq. (3.29). The new expression for m=mo then follows directly as outlined above, and
Fig. 3.93 shows the results.
Triple point 10.1325 Mpa (100 atm) 30.3975 MPa (300 atm)
Uncertainty is estimated to be 0.1% for the liquid phase and 0.5% for the solid phase.
b. Heat of Vaporization. Range of values (J=g):
202 Chapter 3
Figure 3.93 Fraction of liquid hydrogen evaporated due to ortho–para conversion as a func-
tion of time.
Uncertainty is estimated to be 1.24 J=g for p-H2 and 2.5 J=g for n-H2.
The heat of vaporization for parahydrogen and normal hydrogen is tabulated
in Table 3.51 and graphed in Figs. 3.94 and 3.95.
c. Vapor Pressure. The vapor pressure is the pressure, as a function of
temperature, of a liquid in equilibrium with its own vapor.
20.39 20.39
Normal Equilibrium Normal Equilibrium
hydrogen, hydrogen, deuterium, deuterium, Hydrogen
75% o-H2 0.21% o-H2 66.67% o-D2 97.8% o-D2 deuteride
Triple point
temp., K 13.96 13.81 18.72 18.69 16.60
Triple point pressure
atm 0.07105 0.0695 0.1691 0.1691 0.122
kPa 0.07197 0.0704 0.1713 0.1713 0.124
Normal boiling
point, K 20.39 20.27 23.57 23.53 22.13
Critical temp., K 33.19 33.1 38.3 35.9
Critical pressure
atm 12.98 12.8 16.3 14.6
kPa 13.15 13.1 16.5 14.8
Critical volume,
cm3=mol 66.95 65.5 60.3 62.8
204 Chapter 3
Parahydrogen:
ðnpÞ
ln ¼ n1 x þ n2 x2 þ n3 x3 þ n4 ð1 xÞn5 ð3:35Þ
ðnpt Þ
1 Tt =T
x¼
1 Tt =Tc
with
Tt ¼ 13:8 K; TC ¼ 32:938 K
Figures 3.96 and 3.97 present curves for the vapor pressure of p-H2 below and
above p ¼ 1 atm, respectively.
Normal hydrogen:
where
A ¼ 1:55447; B ¼ 31:875; C ¼ 2:39188
and T is in kelvins.
Table 3.52 presents vapor pressure values for the isotopic forms of hydrogen.
The vapor pressure equations of Table 3.52 agree with the experimental data
for hydrogen to within the experimental error over most of the liquid range to the
critical point. The values for deuterium and hydrogen deuteride should not be used
at pressures much above 1 atm.
Figures 3.98 and 3.99 present vapor pressure curves for n-H2.
d. Virial Coefficients. The virial coefficients are commonly defined in two
ways as follows:
and
PV ¼ RT þ B0 ðTÞP þ C 0 ðTÞP2
The first two terms of these series may be inverted by the following:
CðTÞ BðTÞ2
BðTÞ ¼ B0 ðTÞ and C 0 ðTÞ ¼ ð3:37Þ
RT
206 Chapter 3
where P is the pressure, T is the temperature, r is the density, V is the 1=r, and R is
the gas constant. B(T ), B0 (T ) andC(T ), C0 (T ) are virial coefficients of a power series
expansion in density and pressure, respectively.
Both series are theoretically infinite in length; however, the coefficients beyond
the first two [B(T ) and C(T ) ] are of less interest because of their complexity. The
temperature at which B(T) ¼ 0 is called the Boyle point. Either of the above two
Material State A B C
equations is adequate to describe the PVT surface for densities up to about one-half
the critical density.
Range of values:
Figure 3.98 Vapor pressure of normal hydrogen below 0.101325 MPa (1 atm).
208 Chapter 3
Figure 3.99 Vapor pressure of normal hydrogen above 0.101325 MPa (1 atm).
for the J–T curve for parahydrogen are given in Table 3.55. The uncertainty in the
data is estimated to be 5%.
Figure 3.100 shows the J–T inversion curve for parahydrogen.
f. Dielectric Constant. The dielectric constant for fluid parahydrogen may be
calculated from
e1
¼ Ar þ Br2 þ Cr3 ð3:39Þ
eþ2
where A ¼ 0.99575, B ¼ 0.09069, C ¼ 1.1227, and e is the dielectric constant. The den-
sity, r, must be in units of grams per cubic centimeter. The equation is valid over the
range of the tables and may be extrapolated in the fluid phase with reasonable results.
Range of values (dimensionless) :
210 Chapter 3
Specific volume
26 0.207 4.84
26.5 0.211 4.75
27 0.215 4.66
27.5 0.219 4.58
28 0.222 4.49
28.5 0.226 4.42
29 0.230 4.34
29.5 0.234 4.27
30 0.238 4.20
31 0.246 4.06
32 0.254 3.93
33 0.262 3.81
34 0.270 3.70
35 0.278 3.60
36 0.286 3.50
37 0.294 3.40
38 0.302 3.31
39 0.310 3.23
40 0.318 3.15
42 0.334 3.00
44 0.350 2.86
46 0.365 2.74
48 0.381 2.62
50 0.397 2.52
55 0.437 2.29
60 0.477 2.10
65 0.516 1.94
70 0.556 1.80
Above 0 K, however, all matter has thermal energy in the form of rapid
random motion of its atoms, and the balance of forces among the particles becomes
dynamic rather than static. In cooling matter slowly to 0 K, then, classical theory
would predict the loss of thermal energy of the material through the loss of the
kinetic energy of its atoms until at 0 K there would exist perfect motionless order.
Quantum theory, on the other hand, shows that each atom has an irreducible
minimum of kinetic energy amounting to (1=2) hv, where h is Planck’s constant,
6.6 1027 erg sec, and v is the frequency of oscillation of the atom. Even at 0 K,
when a substance has lost all its thermal energy, it will still have this zero-point
energy.
This amount is not large, and in most cases it is effectively inundated by the
thermal energy of matter at higher temperatures. At very low temperatures, however,
the zero-point energy becomes a significant fraction of the total energy of some
substances.
Solid hydrogen, for example, has a zero-point energy of about 200 cal=mol,
which counteracts about 50% of its computed lattice energy of 400 cal=mol.
K
R MPa atm psia mol=cm3 103 lbm=ft3
212 Chapter 3
0.005 1.01515
0.010 1.03046
0.015 1.04594
0.020 1.06158
0.025 1.07739
0.030 1.09336
0.035 1.10950
0.040 1.12580
0.045 1.14226
0.050 1.15889
0.055 1.17569
0.060 1.19265
0.065 1.20977
0.070 1.22705
0.075 1.24449
0.080 1.26210
Another effect of the zero—point energy is to reduce the depth of the potential
well in Fig. 3.101. The ‘‘well’’ region below the horizontal axis actually represents the
liquid phase of a fluid. It is only here that the fluid can exist in two different densities
(interatomic separations) at the same potential energy (temperature). The potential
well therefore represents the liquid–vapor equilibrium region of a fluid. Reducing
the depth of the potential well reduces the temperature range over which vapor–
liquid equilibrium is possible. Therefore, the depth of the potential well is propor-
tional to the heat of vaporization (or sublimation) of the fluid. Hence, the low heat
of vaporization of liquid hydrogen is a macroscopic manifestation of the zero-point
energy of liquid hydrogen.
We shall see shortly that the zero-point energy is manifested even more promi-
nently in liquid helium. It should be noted that neon also exhibits strong quantum
mechanical effects. That is why these three fluids—hydrogen, helium, and neon—
are referred to as the quantum fluids. This is also the reason the PVT properties
of these three fluids are omitted from the usual generalized correlations of state.
These equations depend upon the similarity among fluids, especially the rule that
fluids equally distant from their critical points behave similarly. That is to say, at
the same reduced temperature, TR ¼ T=Tc, and the same reduced pressure,
PR ¼ P=PC, the fluids will have the same compressibility factor, Z ¼ PV=RT.
Helium, hydrogen, and neon do not follow the rules of normal fluids and do not
fit well the generalized equations of state. They are quantum fluids instead.
214 Chapter 3
216 Chapter 3
Figure 3.107 Melting line for parahydrogen, critical point and above.
218 Chapter 3
Figure 3.108 Melting line from triple point to critical point pressure for parahydrogen.
hydrogen is used as a reducing atmosphere and for brazing of metals and silicon
compounds for solid-state components, vacuum tubes, light bulbs, and crystal grow-
ing.
g. Other Uses. Hydrogen use in other areas is difficult to quantify. These use
areas include glass manufacture via the float glass process, oxygen-hydrogen glass
cutting, artificial gem production, refrigeration, pharmaceuticals, nuclear fuel pro-
cessing, and nuclear research using liquid hydrogen bubble chambers.
3.8. Helium
Evidence of the existence of helium (He) was first obtained by Janssen during the
solar eclipse of 1868 when he detected a new line in the solar spectrum; Lockyer
and Frankland suggested the name helium for the new element. In 1895, Ramsay
discovered helium in the uranium mineral clevite, and it was independently
discovered in clevite by the Swedish chemists Cleve and Langlet about the same
time.
Onnes first liquefied helium in 1908 in his laboratory at Leiden, Netherlands,
by using liquid hydrogen precooling in a Joule–Thomson liquefier. Liquid helium
was not available outside of Leiden until 1923.
In 1920, Aston discovered that, in addition to the common isotope helium-4,
there existed a very rare isotope helium-3. The helium-4 atom consists of two elec-
trons orbiting a nucleus of two protons plus two neutrons. The helium-4 nucleus
is the alpha particle associated with radioactive decay and other atomic processes.
The helium-3 atom consists of two electrons orbiting a nucleus of two protons plus
one neutron.
220 Chapter 3
Helium-4 is by far the more common of the two isotopes. Ordinary helium gas
contains about 1.3 104% helium-3, so that when we speak of helium or liquid
helium, we normally are referring to helium-4 (molecular weight 4.0026). Liquid
helium-4 has a normal boiling point of 4.216 K (7.57 R) and a density at the normal
boiling point of 124.96 kg=m3 (7.8 lbm=ft3), or about one-eighth that of water.
The physical constants for helium are shown in Table 3.58.
Figure 3.111 Vapor pressure of normal and para-hydrogen below the triple point.
the entropy difference between the liquid and gaseous states (equal to the heat of
vaporization divided by the temperature) is about 80 cal=(mol K) at the normal boil-
ing point (Trouton’s rule). In liquid helium, the influence of zero-point energy causes
this number to drop to about 20 cal=mol. This zero-point energy accounts for the
very low heat of vaporization of liquid helium. It is as if the zero-point energy
had already started the pot to boil.
The large zero-point energy of liquid helium has a profound and totally unex-
pected effect on its phase diagram.
Figure 3.114 shows the phase diagram of a normal fluid on the left and the
phase diagram of helium on the right. The normal diagram shows that there are
three phases, solid, liquid, and vapor; that these are separated by the melting, vapor
pressure, and sublimation curves; that there is a critical point where the vapor and
liquid phases are no longer distinguishable; and that there is a solid–liquid–vapor
coexistence point, the triple point.
222 Chapter 3
Figure 3.112 Difference in the specific heat of para- and normal hydrogen.
Figure 3.113 Difference in the thermal conductivity of para- and normal hydrogen.
Helium
Chemical Symbol: He
CAS Registry Number: 7440-59-7
DOT Classification: Nonflammable gas
DOT Label: Nonflammable gas
Transport Canada Classification: 2.2
UN Number: UN 1046 (compressed gas); UN 1963 (refrigerated liquid)
US Units SI Units
Physical constants
Chemical formula: He He
Molecular weight 4.00 4.00
Density of the gas at 70 F (21.1 C) and 1 atm 0.0103 lb=ft3 0.165 kg=m3
Specific gravity of the gas at 70 F (21.1 C) 0.138 0.138
and 1 atm
Specific volume of the gas at 70 F (21.1 C) 97.09 ft3=lb 6.061 m3=kg
and 1 atm
Density of the liquid at boiling point and 1 atm 7.802 lb=ft3 124.98 kg=m3
Boiling point at 1 atm –452.1 F –268.9 C
Melting point at 1 atm None None
Critical temperature –450.3 F –267.9 C
Critical pressure 33.0 psia 227 kPa abs
Critical density 4.347 lb=ft3 69.64 kg=m3
Triple point None None
Latent heat of vaporization at boiling point 8.72 Btu=lb 20.28 kJ=kg
and 1 atm
Latent heat of fusion at triple point No T.P. No T.P.
Specific heat of the gas at 70 F (21.1 C)
and 1 atm
CP 1.24 Btu=(lb)( F) 5.19 kJ=(kg)( C)
CV 0.745 Btu=(lb)( F) 3.121 kJ=(kg)( C)
Ratio of specific heats (CP=CV) 1.66 1.66
Solubility in water, vol=vol at 32 F (0 C) 0.0094 0.0094
Weight of the liquid at boiling point 1.043 lb=gal 125.0 kg=m3
The helium phase diagram (Fig. 3.114) exhibits many of the characteristics dis-
played by normal fluids: the existence of a critical point on the liquid–vapor curve,
for instance. One difference between helium and all other substances is the shape of
the solid–liquid (melting) curve. For other substances, this line intersects the liquid–
vapor line at some point, the triple point. For helium, there is no temperature
and pressure at which the solid, liquid, and vapor phases of helium can exist in
equilibrium. Helium has no triple point. Linear extrapolation of the solid–liquid line
would place a triple point at about 1 K. However, from about 2 K down, this extra-
polated line deviates from its higher temperature linear behavior and approaches 0 K
horizontally at a value of 25 atm. This means that helium can never be solidified at
atmospheric pressure merely by reducing its temperature, as is the case with all other
substances. In addition, it is necessary to subject the helium to a pressure of 25 atm.
This is a direct result of the quantum mechanical zero-point energy of helium.
224 Chapter 3
Figure 3.114 Vapor pressure curve for ordinary liquids and helium.
Consider the gedanken experiment of sliding down the helium vapor pressure
curve, starting at the critical point, 5.2 K and 2.24 atm. As the pressure is reduced,
the vapor pressure is reduced correspondingly, reaching the normal boiling point
of 4.2 K and 1 atm. Our progress down the vapor pressure curve continues unevent-
fully until at about 2.17 K there is a totally unexpected event. A totally new liquid
phase of helium-4 is formed that has the properties of zero viscosity, anomalously
large thermal conductivity, zero absolute temperature, and zero viscosity. This is
called the superfluid.
Classical theory predicts a state of perfect order, and hence zero entropy, at
absolute zero. Liquids become solids (freeze) as the temperature is decreased, to
achieve the most orderly, crystalline state. The large zero-point energy of helium,
however, keeps the interatomic separation too large for the atoms to fall into a lat-
tice, that is, to freeze into a solid. Nevertheless, helium must still approach zero
entropy, so it does the next best thing. Instead of forming a solid under its own vapor
pressure, helium creates a second liquid phase of zero entropy, called He II. This
superfluid is helium’s answer to remaining a liquid (the atoms are too far apart
for a solid) and at the same time heading for perfect order. Helium II plays some-
thing of the role that solids do for normal fluids.
Liquid helium-4 is odorless and colorless and somewhat difficult to see in a
container because its index of refraction is so near that of the gas (nr ¼ 1.02 for liquid
He). The heat of vaporization of liquid He at the normal boiling point is 20.73 kJ=kg
(8.92 Btu=lbm), which is only 1=110 that of water.
Table 3.59 gives some of the important fixed points of helium.
Figures 3.115–3.117 are T–S charts for helium-4. The Joule–Thomson
inversion curve is given in Fig. 3.118 and data for the J-T curve are listed in
Table 3.60.
Critical point
T ¼ 5.2014 K (9.3625 R)
P ¼ 2.245 atm (32.99 psia; 2.275 kPa)
r ¼ 0.017399 mol=cm3 (4.348 lbm=ft3)
Normal boiling point
T ¼ 4.224 K (7.604 R)
P ¼ 1 atm (14.696 psia; 1.01325 kPa)
rgas ¼ 0.004220 mol=cm3 (1.054 lbm=ft3)
rliq ¼ 0.03122 mol=cm3 (7.802 lbm=ft3)
Lower lambda point
T ¼ 2.177 K (3.919 R)
P ¼ 0.0497 atm (0.730 psia; 5036 Pa)
rliq ¼ 0.03653 mol=cm3 (9.127 lbm=ft3
Upper lambda point
T ¼ 1.763 K (3.174 R)
P ¼ 29.74 atm (437.1 psia; 3013 M Pa
rliq ¼ 0.04507 mol=cm3 (11.26 lbm=ft3)
226 Chapter 3
usually considered a transition of the ‘‘second order,’’ i.e., the Gibbs energy has a
discontinuity in the second derivative.
As mentioned earlier, He II is often referred to as a superfluid. Immediately
below the lambda point, the flow of liquid through narrow slits or channels becomes
very rapid. Figure 3.121 shows the viscosity of liquid helium measured by an oscillat-
ing disk. Helium I has a viscosity of about 3 106 (Pa s), whereas this experiment
showed that He II has a viscosity of about 1 107 Pa s at 1.3 K. However, Kapitza
and Allen and Misener showed that the viscosity for flow through thin channels
(104–105 cm) was about 1012 Pa s and was independent of the pressure drop or
length of channel and dependent only on the temperature. Thus, the viscosity of
helium below the lambda point is different for bulk flow and for flow through very
thin channels.
To explain these viscosity effects, a ‘‘two-fluid’’ model is used. The liquid is
considered to be composed of two fluids where the total density r is made up of a
normal density rn and a superfluid density component rs such that r ¼ rn þ rs.
228 Chapter 3
The former is termed the normal fluid, while the latter has been designated the
superfluid because under certain conditions the fluid acts as if it had no viscosity.
Figure 3.122 shows the temperature dependence of rn and rs below the lambda
transition. This picture is very useful in explaining many of the physical observations
of He II. However, the detailed theoretical description of the helium problem is
much more complicated, with the two-fluid model being only a first approximation.
To resolve the viscosity paradox, assume that at the lambda point all the fluid
is normal fluid with a normal viscosity and that at absolute zero all the fluid is super-
fluid with zero viscosity. In the thin-channel experiment described above, only the
superfluid atoms, which have zero entropy and do not interact, can flow through
K
R atm psia mol=L lb=ft3
230 Chapter 3
the slit. On the other hand, the oscillating disk is damped by the normal fluid and
thus accounts for the shape of the viscosity curve below the lambda point.
The flow of He II through very thin channels is accompanied by two very inter-
esting thermal effects called the thermomechanical effect and the mechanocaloric
effect.
The thermomechanical effect was discovered in 1938 and is illustrated in
Fig. 3.123. When heat is applied to the fluid in the inner container, the superfluid
component (being cold) tends to move toward a region of higher temperature.
The superfluid can flow very rapidly through the narrow channel, whereas the
flow of the normal component is inhibited by the channel resistance. Thus, a
Figure 3.121 The total viscosity and normal viscosity of helium II.
232 Chapter 3
being zero at the lambda point and becoming nearly temperature-independent below
1.5 K.
Sound can be propagated in liquid He II by at least three mechanisms. First, or
ordinary, sound is the transfer of energy by a pressure wave. Second sound is a tem-
perature wave caused by out-of-phase oscillations of superfluid and normal compo-
nents of the He II. The velocity of second sound is zero at the lambda point and rises
to a value of about 20 m=s between 1 and 2 K.
Third sound is a wave motion present in helium films in which the superfluid
component oscillates but the normal component remains fixed to the walls. The wave
motion appears as an oscillation in the thickness of the film. The velocity if propaga-
tion of third sound is approximately 0.5 m=s.
Yet another form of wave propagation has been postulated to exist in helium-3
at very low temperatures; this form has been termed zero sound. The attenuation of
ordinary sound waves (first sound) in helium-3 becomes infinite as the temperature
goes to zero. In this limit, zero sound replaces first sound. Zero sound is described as
an oscillation in the
pffiffishape
ffi of the Fermi surface. In the limit of T ¼ 0, the velocity of
zero sound will be 3 times the velocity of first sound (at finite temperature). It is not
at all obvious what type of experiments should be tried to detect zero sound. The
situation is analogous to that which occurred after second sound was predicted to
exist in helium-4: Experiments to detect it failed until it was realized that thermal
rather than acoustical techniques were required.
Other Properties. Figures 3.126–3.142 give the density, thermal conductivity,
specific heats, heat of fusion, heat of vaporization, vapore pressure, surface tension,
viscosity, of helium, taken from the WADD report mentioned before (Johnson,
1960).
234 Chapter 3
3.8.2. Helium-3
The isotope helium-3 is very rare and is difficult to isolate from the isotope helium-4.
Interest in helium-3 is two fold. First, its use in low-temperature refrigerators extends
the minimum attainable temperature (outside of using magnetic cooling techniques)
to about 0.3 K. Second, its properties are of fundamental interest in relation to the-
ories of quantum statistical mechanics.
The first liquefaction of helium-3 was reported by Sydoriak, Grilly, and
Hammel in 1948. Although they had only a total of 20 cm3 of gas at STP to work
236 Chapter 3
with, they were able to ascertain the critical point (Tc ¼ 3.35 K, Pc ¼ 875 mmHg) and
measure the saturated vapor pessure from the critical point down to 1.2 K.
Previous to this, it had been seriously questioned whether helium-3 could be
liquefied at all under its saturated vapor pressure. The zero-point energy of
helium-3 should be four-thirds as large as that of helium-4, which is in itself very
large. It was felt that this would make the zero-point energy higher than the potential
energy associated with van der Waals attraction at all interatomic spacings and,
therefore, the gas might not liquefy at any temperature without large external pres-
sures.
The behavior of the light isotope 3He differs from that of 4He. This is partly
due to its greater zero-point energy. Since 3He has only three-fourths the mass of
4
He, its frequency of vibration is higher, and consequently the quantity (1=2) hv is
greater. This zero-point energy effectively counteracts about 95% of the cohesive
energy of 3He so that it very nearly remains a gas all the way to 0 K. Both
liquid helium-3 and liquid helium-4, in fact, display many different attributes more
characteristic of gases than of other liquid, and this is directly traceable to their high
zero-point energy. This is not the whole story, however. Helium-3 obeys a different
set of quantum statistics than helium-4 owing to its antisymmetrical nuclear struc-
ture. For our purposes, the high zero-point energy of these quantum fluids is suffi-
cient explanation.
Liquid helium-3 is a clear, colorless substance having a normal boiling point of
3.19 K (5.74 R) and a density at the normal boiling point of 58.95 kg=m3
(3.68 lbm=ft3). The heat of vaporization of liquid helium-3 at the normal boiling point
is only 8.48 kJ=kg (3.65 Btu=lbm) —so small that there was some doubt in the minds of
early investigators that helium-3 could be liquefied at normal atmospheric pressure.
As in the case of liquid helium-4, liquid helium-3 remains in the liquid state under
its own vapor pressure all the way down to absolute zero. Helium-3 must be
compressed to 28.9 bar (29.3 atm) at 0.32 K (0.576 R) before it will solidify.
An approximate pressure–temperature diagram for helium-3 is shown in Fig.
3.140 These data were used by Arp and Kropschot (1962) to prepare Fig. 3.141
From this figure the practical interest in helium-3 becomes clear. Liquid helium-3
238 Chapter 3
240 Chapter 3
242 Chapter 3
244 Chapter 3
in equilibrium with its vapor at a given pressure is significantly colder than liquid
helium-4 at the same pressure. In an ordinary dewar system, it is difficult to reduce
the helium vapor pressure to less than a corresponding temperature of about 0.8 K.
With a modest (closed cycle) pumping system, it is not difficult to reach 0.3 K by
246 Chapter 3
using helium-3. This decrease of about 0.5 K may seem small, but in reality it is more
like a gain of a factor of about 3 in temperature.
Mixtures of helium-3 and helium-4 are not completely miscible at very low
temperatures, as shown in Fig. 3.142. In fact, as absolute zero is approached,
helium-3 appears to be completely insoluble in helium-4. At other temperatures,
a phase separation occurs in any mixture whose average concentration, at the given
temperature, falls beneath this curve. The concentrations of the separated phases are
given by the two intersections of this curve with the (horizontal) line of constant
temperature. This separation into two liquid phases, and the difference in vapor
pressures, forms the basis for the helium-3=helium-4 dilution refrigerator used
to obtain temperatures close to absolute zero.
248 Chapter 3
can be used, where p is the specific polarizability in cm3=g and r is the density in
g=cm3. The uncertainty of tabulated values of dielectric constant shown in Table 3.61
is estimated to be 0.01%.
e. Melting Curve. The melting pressure of helium can be expressed as
Cryogenic Fluids
Spec. heat (J=(g K)) Thermal
Pressure Densitya Enthalpya Entropya conductivity Viscositya Dielectric Prandtl
Temp. (K) (atm) (g=cm3) 100 (J=g) (J=(g .K)) CV CP (mW=cm. K) (g=(cm s) 106 constant e number Npr
249
(Continued)
3.40 0.4100 13.72 6.393 2.706 2.07 2.95 0.186 35.5 1.05157 0.562
0.7235 29.32 9.463 3.29 6.98 0.0784 9.30 1.00268 0.828
3.50 0.4637 13.60 6.722 2.790 2.11 3.10 0.188 35.1 1.05113 0.580
0.8085 29.52 9.314 3.30 7.18 0.0813 9.65 1.00300 0.852
3.60 0.5220 13.48 7.068 2.875 2.14 3.28 0.190 34.7 1.05066 0.599
0.9008 29.69 9.169 3.30 7.40 0.0843 10.0 1.00334 0.878
3.70 0.5849 13.35 7.432 2.962 2.18 3.47 0.191 34.2 1.05016 0.621
1.001 29.84 9.025 3.31 7.66 0.0874 10.4 1.00371 0.908
3.80 0.6528 13.21 7.816 3.050 2.21 3.68 0.192 33.8 1.04963 0.646
1.109 29.96 8.884 3.32 7.94 0.0906 10.7 1.00411 0.942
3.90 0.7257 13.06 8.221 3.141 2.25 3.92 0.194 33.3 1.04906 0.674
1.228 30.06 8.743 3.33 8.27 0.0939 11.1 1.00455 0.980
4.00 0.8040 12.90 8.650 3.234 2.28 4.19 0.195 32.8 1.04845 0.707
1.356 30.12 8.603 3.33 8.65 0.0974 11.5 1.00503 1.02
4.10 0.8878 12.73 9.105 3.330 2.32 4.51 0.196 32.3 1.04780 0.745
1.496 30.15 8.463 3.34 9.10 0.101 11.9 1.00555 1.07
4.20 0.9772 12.54 9.508 3.429 2.36 4.88 0.196 31.8 1.04710 0.791
1.649 30.14 8.322 3.35 9.63 0.105 12.4 1.00612 1.13
4.224 1.000 12.50 9.711 3.454 2.37 4.98 0.196 31.7 1.04692 0.884
1.689 30.13 8.287 3.35 9.78 0.106 12.5 1.00626 1.15
4.30 1.073 12.35 10.10 3.532 2.40 5.32 0.197 31.3 1.04635 0.847
1.818 30.08 8.177 3.35 10.3 0.110 12.8 1.00674 1.20
4.40 1.174 12.13 10.66 3.640 2.44 5.86 0.197 30.7 1.04553 0.916
2.005 29.98 8.029 3.36 11.1 0.114 13.3 1.00744 1.28
Chapter 3
4.50 1.282 11.89 11.25 3.753 2.49 6.55 0.199 30.1 1.04463 0.994
2.213 29.81 7.876 3.37 12.1 0.128 13.8 1.00821 1.39
4.60 1.397 11.63 11.90 3.873 2.54 7.44 0.200 29.5 1.04364 1.10
251
© 2005 by Marcel Dekker
5367-4 Flynn Ch03 R2 090804
252 Chapter 3
Cryogenic Fluids
T (K) 0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09
0.5 0.016342 0.022745 0.031287 0.042561 0.057292 0.076356 0.10081 0.13190 0.17112 0.22021
0.28121 0.35649 0.44877 0.56118 0.69729 0.86116 1.0574 1.2911 1.5682 1.8949
2.2787 2.7272 3.2494 3.8549 4.5543 5.3591 6.2820 7.3365 8.5376 9.9013
11.445 13.187 15.147 17.348 19.811 22.561 25.624 29.027 32.800 36.974
41.581 46.656 52.234 58.355 85.059 72.386 80.382 89.093 98.567 108.853
1.0 120.000 132.070 145.116 159.198 174.375 190.711 208.274 227.132 247.350 269.006
292.169 316.923 343.341 371.512 401.514 433.437 467.365 503.396 541.617 582.129
625.025 670.411 718.386 769.057 822.527 878.916 938.330 1000.87 1066.67 1135.85
1208.51 1284.81 1364.83 1448.73 1536.61 1628.62 1724.91 1825.58 1930.79 2040.67
2155.35 2274.99 2399.73 2529.72 2665.08 2805.99 2952.60 3105.04 3263.48 3428.07
1.5 3598.97 3776.32 3960.32 4151.07 4348.79 4553.58 4765.68 4985.18 5212.26 5447.11
5689.88 5940.76 6199.90 6467.42 6743.57 7028.47 7322.31 7625.21 7937.40 8259.02
8590.22 8931.18 9282.06 9643.02 10014.3 10395.9 10788.2 11191.2 11605.1 12030.1
12466.1 12913.7 13372.8 13843.6 14326.1 14820.7 15327.3 15846.3 16377.7 16921.7
17478.2 18047.7 18630.1 19225.5 19834.1 20455.9 21091.1 21739.7 22402.0 23077.9
2.0 23767.4 24470.9 25188.1 25919.2 26664.2 27423.3 28196.3 28983.2 29784.2 30599.1
31428.1 32271.1 33128.0 33998.6 34882.8 35780.3 36690.9 37614.3 38550.2 39500.3
40465.6 41446.6 42443.5 43456.5 44485.7 45531.3 46593.5 47672.5 48768.6 49881.8
51012.3 52160.2 53325.8 54509.2 55710.5 56930.0 58167.8 59423.8 60698.8 61992.0
63304.3 64635.2 65985.4 67354.8 68743.5 70152.0 71580.2 73028.1 74496.0 75984.2
2.5 77493.1 79022.2 80572.2 82142.9 83734.6 85347.2 86981.2 88636.7 90313.8 92012.6
93733.1 95476.0 97240.8 99028.2 100838 102669 104525 106403 108304 110228
112175 114145 116139 118156 120198 122263 124353 126465 128603 130765
132952 135164 137401 139663 141949 144260 146597 148961 151349 153763
156204 158671 161164 163684 166230 168802 171402 174028 176682 179364
3.0 182073 184810 187574 190366 193187 196037 198914 201820 204755 207719
210711 213732 216783 219864 222975 226115 229285 232484 235714 238974
(Continued)
253
© 2005 by Marcel Dekker
254
Table 3.64 (Continued )
T (K) 0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09
242266 245587 248939 252322 255736 259182 262658 266166 269706 273278
(Continued)
276880 280516 284183 287883 291615 295380 299178 303008 306871 310768
314697 318659 322654 326684 330747 334845 338976 343141 374341 351575
3.5 355844 360147 364485 368860 373269 377714 382194 386710 391262 395849
400471 405130 409825 414556 419324 424128 428968 433846 438760 443713
448702 453729 458794 463897 469038 474218 479435 484691 489985 495317
500688 506098 511547 517036 522564 528132 533739 539387 545075 550805
Chapter 3
© 2005 by Marcel Dekker
5367-4 Flynn Ch03 R2 090804
Temp. (K) 3 atm 5 atm 6 atm 10 atm 15 atm 30 atm 50 atm 70 atm
The only helium plant in the western world outside the United States is near Swift
Current, Saskatchewan.
Helium is a low-density inert gas with a combination of unique properties that
make its availability essential for a wide range of industrial and scientific endeavors.
These include applications not only in cryogenics, but also in special welding
technology, space exploration, chromatography, heat transfer, and controlled
atmosphere.
a. Pressurizing and Purging. Most of the helium for pressurizing and purging
is used in the space program. Helium is used to pressurize liquid oxygen, liquid
hydrogen, kerosene, and the hypergolic propellants. Helium is used also to pressurize
pneumatic systems on boosters and spacecraft. A Saturn booster such as that used
on the Apollo lunar missions requires about 13 106 ft3 of helium for a firing plus
more for checkouts. At the height of the Apollo program, almost one-third of all
domestic helium production was used for this purpose.
Temp. (K) 3 atm 5 atm 10 atm 20 atm 40 atm 60 atm 80 atm 100 atm
256 Chapter 3
4
Mechanical Properties of Solids
1. INTRODUCTION
A knowledge of the properties and behavior of materials used in any cryogenic sys-
tem is essential for proper design considerations. Often the choice of materials for
the construction of cryogenic equipment will be dictated by consideration of
mechanical and physical properties such as thermal conductivity (heat transfer along
a structural member), thermal expansivity (expansion and contraction during cycling
between ambient and low temperatures), and density (weight of the system relative to
its volume). Since properties at low temperatures are often significantly different
from those at ambient temperature, there is no substitute for test data, and fortu-
nately there are now several excellent data compilations (McClintock and Gibbons,
1960; Durham et al., 1962; Campbell, 1974; Johnson, 1960). To help make sense out
of all of the data that do exist, and to help estimate properties when no data are
present, it is useful to have certain general rules in mind. Such is the purpose of
the following discussion.
In any branch of engineering design, the choice of material is dictated by ques-
tions of safety and economy. The range of choice in cryogenic design is perhaps lim-
ited by the issue of low-temperature embrittlement. The purpose of this section is to
show which materials that are ductile at ordinary temperatures retain their ductility
at low temperatures. This determination requires an explanation of why ductility
exists at all in some materials and not others and a knowledge of the temperature-
related mechanisms that govern ductility. Common indices of strength of materials
are also explored along with the temperature dependence of this important property.
When a bar of a structural solid is subjected to an elongating force, the first response
is a slight stretching that is directly proportional to the force producing it. If the
elongating force is released during this initial period, the material will return to very
nearly its original length. This is called elastic behavior, and all solid materials exhi-
bit it in some measure. At some higher value of stress, however, the material will no
longer behave elastically. It may break without further deformation (brittle beha-
vior), or it may take on a permanent deformation (ductile behavior).
The ductile or brittle behavior of structural materials is usually determined by
the familiar stress–strain test shown in Fig. 4.1. The figure represents a steel bar of
cross-section a, on which a gauge length l has been marked. Let P be the value of any
257
258 Chapter 4
axial tensile load that gradually increases from zero value until the steel bar breaks.
As the value of P increases, values of the elongation e are taken.
Stress:
S ¼ P=A
Strain:
E ¼ e=l
The relationship between s and E is the stress–strain curve found experimentally and
is shown in Fig. 4.2.
Another class of materials (nonmetals) are capable of extreme great elastic
deformations. These are called elastomeric materials.
Usually, materials exhibiting plastic deformation under stress are the more
desirable for structures. Ductility is desirable so that accidental stresses beyond
design values can be redistributed to safer levels by means of plastic flow.
Brittle materials have no such mechanism to protect them from excessive stress.
When a local overstress occurs in a brittle material, the result is often failure of the
piece rather than a deformation. It is useful to know, therefore, which structural
materials that are acceptable (ductile) at normal temperatures will remain ductile
at low temperatures. This requires that we know the mechanisms of ductility in solids
and how those mechanisms change with temperature. To begin, let us examine the
structure of solids.
Figure 4.2 The stress–strain relationship of a ductile material (a) and a brittle material (b).
The theory of elasticity treats solids as continuous elastic media. Solids are not
continuous media. They are composed of atoms bound together in more or less
regular arrays. Essentially, three different types of solid structures exist, represented
by glasses, plastics, and metals.
3.1. Glasses
The structure of an idealized glass is shown in Fig. 4.3 (Jones, 1956). The chemi-
cal composition of this glass is G2O3, where G represents a metallic element and
O represents oxygen. The principal characteristic of the glass to be noted is the lack
of any spatial order, even though the stoichiometric relationship between G and O is
observed. Liquids also show such lack of order, and glasses can be described to a first
approximation as extremely viscous liquids.
3.2. Plastics
Plastics, or polymers, are composed of giant long-chain molecules. Polymer mole-
cules may have from tens to thousands of atoms each and are essentially linear. They
lie in tangled disarray; a segment of a typical polymer molecule is shown in Fig. 4.4.
The intermolecular force that binds the polymer molecules to one another is the
rather weak van der Waals force. Some cross-linking among the polymer molecules
may also be present. Cross-linking is actual chemical bonding between adjacent
molecules and is much stronger than the van der Waals forces. The more highly
cross-linked (thermoset) polymers form much more rigid solids than those that are
only sparingly cross-linked (thermoplastics).
3.3. Metals
Compared to glasses and polymers, metals have a highly ordered structure. Metal
atoms are arranged in symmetrical crystal lattices. The three most common metal
crystal lattices are shown in Fig. 4.5 (Quarrel, 1959). The face-centered cubic (fcc)
260 Chapter 4
lattice consists of a cube that has an atom at each of its eight corners and one at the
center of each of the six crystal faces. The hexagonal close-packed (hcp) lattice con-
sists of a right hexagonal prism with an atom at each of its 12 vertices, an atom at the
center of each of its two regular hexagonal ends, and three more atoms located mid-
way between the ends of the prism, each on a line parallel with the vertical axis of the
prism and running through the center of the equilateral triangle formed by three
neighboring atoms in the end faces. The body-centered cubic (bcc) lattice is a cube
with an atom at each of its eight corners and an atom at the center of the body of
the cube. Other more complicated metal structures exist, but only these three basic
types will be considered here.
The drawings of Fig. 4.5 are useful for visualizing the slip planes that exist in
various metal crystal systems. For example, from the fcc lattice shown in Fig. 4.5a, it
is easy to see how planes of atoms might slide rather easily over one another from
upper left to lower right. This lattice also has several other such slip planes. Although
it is not readily apparent from the figure, a three-dimensional model of these crystal
systems would show that the bcc lattice offers the least number of slip planes of the
three shown, and the hcp system falls in between the fcc and bcc systems.
Real crystals do not have the perfect geometric arrangement of atoms shown in
Fig. 4.5 (Van Buren, 1960). For example, there are always some impurity atoms
present. These impurity atoms cause irregularities in the crystal lattice. Real
crystals may also have locations where atoms are simply out of place. An extra atom
in the interstices between regular crystal sites is called an interstitial atom. A
missing atom is called a vacancy. Interstitial atoms and vacancies are examples of
‘‘point defects’’ in the crystal.
Extended imperfections can also occur. An extra plane of atoms can be
present in a crystal lattice. Extended imperfections of this type are called edge
dislocations.
Another type of extended imperfection is common in crystals: the screw dislo-
cation. This type of dislocation occurs when part of a crystal has slipped one atomic
distance relative to its adjacent part.
Both dislocations and point defects are nearly always present to some degree in
real metallic crystals (Cottrell, 1953). The concept of a geometrically perfect crystal is
useful in explaining certain properties of crystals, but there are other properties that
can be adequately explained only by consideration of the defects. The concept of
dislocations, for example, helps explain the yield stress and ductility of metals.
The regular arrangement of atoms in the crystal (and not so much its defects) helps
explain the modulus of elasticity.
4. DUCTILITY
As mentioned, the role of dislocations is useful in explaining the yield strength and
ductility of certain crystalline solids.
One could account for elasticity by imagining a regular crystal structure with-
out dislocations such as those of Fig. 4.5. If a certain horizontal shearing stress were
applied to the top of the crystal lattice, then the structure would bend, but it would
snap back to its original shape when the stress was removed. Above a certain stress,
however, plastic shear would suddenly take place along one of the planes of atoms,
resulting in a permanent deformation.
The ‘‘perfect crystal’’ was the model originally used to explain the plastic defor-
mation of metals. However, when the stress necessary to shear a perfect crystal was
calculated on the basis of interatomic forces, the results were a thousand times
greater than the experimentally observed stresses. Taking the presence of disloca-
tions into account gives better results, because the stress required to force a disloca-
tion through a crystal is of the required magnitude (Barrett, 1957). The movement of
one plane of atoms over another in a crystal actually occurs progressively as the
dislocation moves through the crystal. The resulting crystal deformation is the same
as if the process had taken place suddenly. Shearing a perfect crystal lattice structure
is like moving a heavy rug all at once; forcing a dislocation through the crystal is
like moving a wrinkle through the rug. The displacement of the rug is the same in
both cases.
The ability of crystalline materials to deform plastically (their ductility)
depends primarily on the mobility of dislocations within the crystal (Fisher, 1957).
The mobility of dislocations depends on temperature, as one might expect, and on
several other things as well. The number of slip systems available in a crystal influ-
ences dislocation mobility. Slip systems are directions within the crystal in which the
planes can slip easily over one another. Dislocations can move most easily through a
crystal in the direction of slip planes.
The movement of dislocations can also be impeded in many ways. Impurity
atoms at lattice sites can lock dislocations in place. Vacancies and interstitial atoms
also provide obstacles to the movement of dislocations. The existence of hard foreign
particles in the crystal can also impede the mobility of dislocations.
Grain boundaries, which are discontinuities in the orientation of neighboring
crystals, also are obstacles. A metal having a small grain size will have a higher yield
stress than one of large grain size. One reason annealed metals are more ductile than
work-hardened ones is that work hardening reduces the size of grains. Finally,
dislocations themselves provide obstacles to the progress of other dislocations.
Glasses have no slip systems and are therefore very brittle. When stress is
applied to a piece of glass, the atomic bonds rupture at some point on the surface.
The crack then propagates, causing a fracture. The final stage of rupture of a metal
262 Chapter 4
4.1.1. Metals
Since the number of slip systems is not usually a function of temperature, the ducti-
lity of fcc metals is relatively insensitive to a decrease in temperature (Brick, 1953).
Metals of other crystal lattice types tend to become brittle at low temperatures. Crys-
tal structure and ductility are related because the fcc lattice has more slip systems
than the other crystal structures. In addition, the slip planes of bcc and hcp crystals
tend to change at low temperature, which is not the case for fcc metals. Therefore,
copper, nickel, all of the copper–nickel alloys, aluminum and its alloys, and the
Copper Iron
Nickel Carbon and low alloy steels
All copper–nickel alloys
Aluminum and all its alloys Molybdenum
Austenitic stainless steels containing Niobium
more than 7% nickel Zinc
Most bcc metals
Zirconium Most plastics
Titanium
Most fcc metals
Polytetrafluoroethylene
austenitic stainless steels that contain more than approximately 7% nickel, which are
all face-centered cubic, remain ductile down to low temperatures if they are ductile at
room temperature. Iron, carbon and low-alloy steels, molybdenum, and niobium,
which are all body-centered cubic, become brittle at low temperatures. The hexago-
nal close-packed metals occupy an intermediate place between fcc and bcc behavior.
Zinc undergoes a transition to brittle behavior in tension, whereas zirconium and
pure titanium remain ductile (Corruccini, 1957).
The thermal vibration of atoms in the crystal lattice is strongly temperature-
dependent and is less effective in assisting dislocation motion at low temperatures.
The interaction of dislocations with thermal vibrations is complicated, but it is non-
etheless satisfying to find that ductility usually decreases somewhat with a decrease
in temperature.
The complete situation regarding brittleness in metals at low temperature
depends on more than just crystal structure. For instance, bcc potassium and beta
brass remain ductile down to 4.2 K, and the bcc metals lithium and sodium show
no signs of brittleness down to 4.2 K. In general, however, brittleness will not
occur in face-centered cubic metals in which dislocations cannot be firmly locked by
impurity atoms.
4.1.2. Elastomers
As mentioned, the tendency to show a ductile–brittle transition is correlated with the
lattice type. Thus, the fcc metals show only a few cases of this effect and for struc-
tural purposes may be regarded as almost uniformly well behaved. These include
copper, nickel, aluminum, the solid solution alloys of each of these, and the austeni-
tic stainless steels. In contrast, bcc metals for the most part show brittle behavior
(though the transition zones of some can be depressed to low temperatures). The fer-
ritic steels are by far the most prominent of these. Prominent among the established
structural metals with a different lattice are the hexagonal close-packed metals mag-
nesium and titanium. The impact strength of magnesium is low at all subambient
temperatures because its brittle transition zone is above room temperature. Tests
on commercial titaniums indicate that ductility is retained in tension to low tempera-
tures if the amounts of the interstitial solutes carbon, oxygen, nitrogen, and hydro-
gen are small. However, notched-bar impact tests show a transition above ambient
temperature.
For temperatures much below 200 K, it is common practice to use the
face-centered cubic metals almost exclusively, especially where shock and vibration
are encountered. However, less expensive steels can be used in many less critical
applications, especially for temperatures above 150–200 K.
All but one plastic or elastomer become brittle at low temperatures. Polytetra-
fluoroethylene (PTFE) is unique in that it can still be deformed plastically to a small
degree at 4 K.
Plastics and elastomers do not respond to stress as metals do. The less cross-
linked elastomers yield by uncoiling their long-chain molecules and by sliding over
one another. The thermal energy of the material at room temperature facilitates this
motion.
At low temperatures, however, the attractive intermolecular forces are more
effective than the thermal energy ‘‘lubricant,’’ and the material deforms less readily.
This effect is especially pronounced through its ‘‘glass transition’’ in the temperature
range at the onset of brittleness.
264 Chapter 4
266 Chapter 4
Figure 4.9 The curve of ductility vs. temperature for several important structural materials.
(1) 2024 T4 aluminum; (2) beryllium copper; (3) K Monel; (4) titanium; (5) 304 stainless steel;
(6) C1020 carbon steel; (7) 9% nickel steel.
temperature, a material may exhibit considerable ductility in the tensile test but not
in the impact test.
Figure 4.10 shows this effect for an ordinary carbon steel. This figure illustrates
that the ductility of a material is affected by the type of stress system and the rate of
application of this stress system. Between T2 and T3, for example, the carbon steel
displays ductile behavior in a simple uniaxial stress system (tensile test) and displays
brittle characteristics at high rates of loading (impact test). Increasing either the
strain rate or the complexity of the stress system moves the curves to the right. This
amounts to an increase in the brittle transition temperature. Uniformly ductile or
brittle behavior is observed above T4 and below T1.
The tensile elongation and the impact energy vs. temperature of annealed oxy-
gen-free copper are shown in Fig. 4.11. The tensile elongation of copper is quite high
Figure 4.11 Tensile elongation and impact energy of annealed oxygen-free copper.
Figure 4.12 Tensile elongation and impact energy of a chromium–nickel stainless steel (AISI
316, annealed).
268 Chapter 4
Figure 4.13 Tensile elongation and impact energy of a chromium stainless steel (AISI 430,
annealed).
In general, the ultimate tensile strength of a solid material is greater at low tempera-
tures than it is at ordinary temperatures. This is true for both crystalline and
noncrystalline solids and for many heterogeneous materials as well (e.g., glass-
reinforced plastics). For metals, the strength at 4 K may be 2–5 times that at room
temperature. For plastics, the strength at 76 K may be 1.5–8 times greater than
the room-temperature value. Glasses show less change in strength at low tempera-
tures; at 76 K, glasses have between 1.5 and 2 times their room-temperature strength.
Reduction of the thermal energy of the metal lattice at lower temperatures is
responsible for part of the increase in the strength of metals. The decreased thermal
vibration at low temperatures also results in stronger plastics, because it becomes
more difficult for the long-chain molecules to slide over one another.
The change in strength of glass at low temperatures is not as great as that of
metals and plastics. The important cohesive forces of glasses are the chemical bonds
between atoms (rather than the intermolecular forces of plastics). Such chemical
bond forces are less affected by temperature.
There are exceptions to the above general rule. The thermal and mechanical
history of a material, various solid-state phase changes, the form of stress system,
and the rate of application of this stress system can also affect mechanical properties.
Figure 4.14 Tensile elongation and impact energy of plain carbon steels.
Figure 4.15 The Charpy impact strength at low temperatures for several common structural
materials. (1) 2024 T4 aluminum; (2) beryllium copper; (3) K Monel; (4) titanium; (5) 304
stainless steel; (6) C1020 carbon steel; (7) 9% nickel steel.
Some of the exceptions to the general rule of increasing strength with decreasing
temperature can be explained qualitatively by phase changes during deformation.
Some austenitic stainless steels undergo a partial transformation under strain from
a face-centered cubic to an intermediate hexagonal close-packed phase and then to
a body-centered cubic lattice (martensite). The martensite form is stronger than
the parent structure and contributes to the high strength of these materials. The ten-
sile strength of austenitic steels shows the usual increase with decreasing temperature
to some low temperature and then a decrease, reflecting a maximum in the phase
transformation curve.
5.1. Metals
The yield strength of ductile metals increases with decreasing temperature. The yield
strength at 20 K is between one and three times the room-temperature yield strength
for fcc and hcp metals. An even greater tendency toward increases in yield strength
at low temperatures is shown for bcc metals. In fact, for body-centered cubic steels,
the yield strength increases so rapidly with decreasing temperature that it becomes
greater than the fracture strength. Thus, below a certain temperature, these materials
fracture before they reach their yield strengths.
The yield and ultimate tensile strengths of annealed copper are shown in
Fig. 4.16. This figure shows that the yield strength of a material can be insensitive
to temperature while its tensile strength is increasing by more than a factor of 2.
The material has ‘‘work hardened.’’ That is, the material has hardened itself by
generating obstacles to dislocation motion.
The behavior of a plain carbon (bcc) steel is shown in Fig. 4.17 for contrast.
This behavior is typical of the metals that become brittle at low temperatures. The
yield and tensile strength curves approach one another as temperature is reduced.
When the two curves converge, the material is brittle in the tensile test.
270 Chapter 4
5.2. Plastics
Only a few plastics have been tested at temperatures below 200 K. Of these, only
Teflon showed ductility down to the lowest test temperature, which was 4 K. How-
ever, reinforced plastics such as the glass fiber laminates can have good properties,
the tensile strength parallel to the laminations increasing at low temperatures and
the modulus being approximately constant. Mylar breaks with fragmentation in a
tensile test and so is obviously brittle. Yet in films about 2.540 105 m (0.001 in.)
or less in thickness, Mylar shows remarkable flexibility in bending tests at tempera-
tures as low as 20 K.
Figure 4.18 shows the tensile strength of several plastics. The increase in
strength as the temperature is depressed is accompanied by a rapid decrease in elon-
gation and impact resistance. The glass-reinforced plastics are the only plastics that
retain appreciable impact resistance as the temperature is lowered. Figure 4.19 shows
the tensile strength and impact energy of a glass fiber-reinforced epoxy resin sample.
These glass-reinforced epoxies are very useful for cryogenic applications because of
their high strength-to-weight ratios and high strength-to-thermal conductivity ratios
(Weitzel et al., 1960).
5.3. Glass
The strength of glass at room temperature varies inversely with load duration, is
sensitive to atmospheric water vapor, and is also sensitive to rather minute surface
Figure 4.17 Tensile and yield strength of a plain carbon steel (SAE 1010 m, cold-rolled).
Figure 4.18 Tensile strength of selected plastics. (1) Polyethylene terephthalate (Mylar); (2)
polytetrafluoroethylene (Teflon); (3) polytrifluoromonochloroethylene (Kel-F); (4) polyvinyl
chloride; (5) nylon.
Figure 4.19 Tensile strength and impact energy of a glass fiber-reinforced epoxy resin.
272 Chapter 4
Table 4.2 Breaking Stress of a Borosilicate Crown Optical Glass—Median Values from
Probability Plots
Breaking stress
Condition psi=s MPa=s psi Mpa psi MPa psi MPa psi MPa
Abraded 800 5.52 7500 51.7 9500 65.5 10,400 71.71 10,400 71.71
Abraded 10 0.0689 5500 37.9 7500 51.71 10,400 71.17 10,600 73.08
Abraded 1 0.00689 5000 34.5 6400 44.13 10,400 71.17 10,200 70.33
Unabraded 800 10,400 71.7 18,000 124.11
Figure 4.20 Ultimate stress for several engineering materials. (1) 2024 T4 aluminum; (2) ber-
yllium copper; (3) K Monel; (4) titanium; (5) 304 stainless steel; (6) C1020 carbon steel; (7) 9%
nickel steel; (8) Teflon.
Figure 4.21 Yield stress for the same engineering materials as in Fig. 4.20.
c. Monel
d. Beryllium copper
e. Teflon
274 Chapter 4
6. ELASTIC CONSTANTS
Elastic constants are physical properties, which relate stress to strain, or force per
unit area to relative length change. Either stress or strain can be induced by any elas-
tically coupled force—mechanical, thermal, magnetic, or electrical. Since solids resist
both volume change and shape change, they have at least two independent elastic
constants. Most materials considered in this chapter are quasi-isotropic; their macro-
scopic elastic constants do not depend on direction. They are characterized elasti-
cally by two independent elastic constants. Composites and textured aggregates
are anisotropic rather than quasi-isotropic. They are characterized by five indepen-
dent elastic constants for the transverse-isotropic case and nine constants for the
orthotropic case. This section defines and describes the elastic constants used most
often to characterize polycrystals.
6.1. Compressibility, KT
All matter—gas, liquid, or solid—responds to pressure change and exhibits at least
one elastic constant, the compressibility. It is given by
KT ¼ ð1=V Þð@V =@PÞT ð4:1Þ
where P, T, and V denote pressure, temperature, and volume. Figure 4.22 illustrates
the mechanical deformations with which compressibility is associated. KT is always
positive and has units of reciprocal pressure.
kT ¼ ðA=‘Þðd‘=dLÞT ð4:2Þ
where A denotes the cross-sectional area of a rod-shaped, hydrostatically loaded spe-
cimen, l denotes length, and L denotes load. For isotropic solids, linear and bulk
compressibilities relate according to
3k ¼ KT ð4:3Þ
B has pressure units, and, herein, elastic constants with pressure units are called
elastic moduli. For describing solids, hydrostatic stress, s, is often preferred over
pressure, which is a negative stress. Thus,
BT ¼ V ð@s=@V ÞT ð4:5Þ
A high bulk modulus usually implies strong interatomic forces, high cohesive energy,
high melting point, and high elastic-stiffness moduli.
276 Chapter 4
Figure 4.24 Mechanical deformations with which shear modulus and torsional modulus are
associated.
Figure 4.24 illustrates the mechanical deformations with which the shear
modulus is associated. A cube sheared on one face has a shear strain related to
the length change of the face diagonal by
g ¼ 2D‘=‘ ð4:9Þ
In torsion (see Fig. 4.24), t is the torsional stress. For isotropic materials, simple-
shear and torsional G’s are identical. Both shear and torsion conserve volume but
change shape. The bulk modulus, as described above, describes volume change
without shape change.
Figure 4.25 Mechanical deformation and strains with which Poisson’s ratio is associated.
mean that a material tends to maintain constant volume during uniaxial deforma-
tion. Liquids exhibit constant volume during deformation, and some solids (such
as rubber) tend to.
The four elastic constants just described—B, E, G, and n—are sometimes
referred to as the bulk, engineering, and macroscopic, polycrystalline, practical, or
technological elastic constants. These four constants are those used most frequently
to describe materials such as polycrystalline aggregates.
B E G n
3BE 1 E
B, E – –
9B E 2 6B
9BG 1 3B 2G
B, G – –
3B þ G 2 3B þ G
3B 1 2v
B, v – 3B (12v) –
2 1þv
GE E
E, G – – 1
3ð3G EÞ 2G
E E
E, v – –
3ð1 2vÞ 2ð1 þ vÞ
2G 1 þ v
G, v 2G (1 þ v) – –
3 1 2v
278 Chapter 4
T (K) Ni Fe Cu Ti Al Mg
technological alloys. These data are taken from "Materials at Low Temperatures,"
Am. Soc. Metals (1983).
Table 4.8 gives ambient-temperature elastic constants for 24 cubic elements for
which reliable data exist. This table also includes mass density, r.
Table 4.9 gives zero-temperature values extrapolated from 4 K data. Since the
values shown in Tables 4.4–4.9 come from many different sources, some small differ-
ences in values will occur among the tables.
The data in Tables 4.4–4.9 are all taken from ‘‘Materials at Low Tempera-
tures,’’ Am. Soc. Metals, 1983.
Table 4.5 Young’s Modulus of Six Metals at Selected Temperatures in Units of GPa
T (K) Ni Fe Cu Ti Al Mg
Table 4.6 Shear Modules of Six Metals at Selected Temperatures in Units of GPa
T (K) Ni Fe Cu Ti Al Mg
7. MODULUS OF ELASTICITY
The constant of proportionality relating tensile stress and strain in the range of elas-
tic response of a material is the modulus of elasticity, or Young’s modulus, E. Two
other commonly used elastic moduli are (1) the shear modulus G, which is the rate of
change of shear stress with respect to shear strain at constant temperature in the
elastic region, and (2) the bulk modulus B, which is the rate of change of pressure
(corresponding to a uniform three-dimensional stress) with respect to volumetric
strain (change in volume per unit volume) at constant temperature.
T (K) Ni Fe Cu Ti Al Mg
280 Chapter 4
E
B¼ ð4:14Þ
3ð1 2mÞ
(in some texts, the symbol K is used for the bulk modulus),
E
G¼ ð4:15Þ
2ð1 þ mÞ
2Gð1 þ mÞ
E¼ ð4:16Þ
1þm
where
E 2G
m¼
2G
vt
r B E G vm vl (cm ms1)
(g=cm3) A (GPa) (GPa) (GPa) v (cm ms1) (cm ms1) y (K)
Ag 10.634 2.9912 108.72 87.02 31.84 0.3666 0.1950 0.3770 0.1730 226.46
Al 2.732 1.2073 79.38 78.24 29.29 0.3357 0.3674 0.6583 0.3274 430.56
Au 19.478 2.8436 180.32 83.19 29.23 0.4231 0.1391 0.3355 0.1225 161.74
C 3.516 1.2109 442.00 1139.37 532.23 0.0704 1.3416 1.8099 1.2304 2239.62
Cr 7.231 0.6865 201.00 300.58 120.16 0.2508 0.4526 0.7068 0.4076 589.57
Cu 9.024 3.1904 142.03 132.51 49.28 0.3445 0.2625 0.4798 0.2337 344.40
Fe 7.922 2.3431 169.00 213.17 82.64 0.2898 0.3603 0.5936 0.3230 472.44
Ge 5.340 1.6651 76.52 132.92 54.90 0.2105 0.3544 0.5295 0.3207 373.36
Ir 22.659 1.5698 366.67 557.12 223.43 0.2468 0.3485 0.5416 0.3140 429.62
K 0.910 7.6267 3.66 3.10 1.14 0.3587 0.1261 0.2386 0.1120 90.52
Li 0.5375 9.2400 13.30 10.70 3.90 0.3650 0.3049 0.5868 0.2706 326.74
Mo 10.223 0.9024 265.29 335.72 130.22 0.2891 0.3481 0.6552 0.3569 474.53
Na 1.005 8.5000 7.24 5.98 2.20 0.3623 0.1664 0.3181 0.1478 147.43
Nb 8.616 0.5183 173.03 110.87 39.79 0.3932 0.2430 0.5122 0.2149 276.60
Ni 8.972 2.3859 187.60 233.52 90.33 0.2925 0.3541 0.5860 0.3173 475.98
Pb 11.593 3.8379 48.79 30.34 10.86 0.3964 0.1095 0.2336 0.0968 105.33
Pd 12.092 2.4552 195.43 134.72 48.63 0.3851 0.2266 0.4640 0.2005 275.92
Pt 21.578 1.4828 288.40 183.54 65.83 0.3939 0.1976 0.4175 0.1747 238.44
Rb 1.629 8.0657 3.21 2.39 0.87 0.3761 0.0824 0.1638 0.0730 55.31
Si 2.331 1.5640 99.22 163.21 66.57 0.2258 0.5916 0.8980 0.5344 648.87
Ta 16.754 1.6154 194.21 191.15 71.54 0.3360 0.2319 0.4158 0.2066 263.77
Th 11.888 3.4545 58.10 75.58 29.45 0.2832 0.1754 0.2862 0.1574 163.83
V 6.050 0.8130 157.04 135.23 49.85 0.3565 0.3230 0.6078 0.2870 399.19
W 19.313 0.9959 314.15 417.76 163.40 0.2784 0.3240 0.5248 0.2909 384.39
or
3K E
m¼
6K
It has been found that Poisson’s ratio for isotropic materials does not change
appreciably with change in temperature in the cryogenic range. Therefore, all three
of the elastic moduli above vary in the same manner with temperature. Accordingly,
in the discussion that follows, ‘‘modulus’’ shall mean Young’s modulus E, which is
the slope of the initial portion of the stress–strain curve of Fig. 4.1.
We expect materials to get stiffer at low temperatures, and that is the usual
case. However, the change in modulus of elasticity with temperature is not very great
in metals compared to the change in strength over the same temperature range. Some
general statements can be made regarding the temperature dependence of the elastic
modulus.
The elastic moduli of most polycrystalline metals increase by about 10% in
going from room temperature to about 20 K, and below this there is little change
(Campbell, 1974). For some alloys, the modulus of elasticity may either decrease
or increase in value or even exhibit maxima or minima as the temperature is lowered,
depending upon their composition. The moduli of noncrystalline plastics are 2–20
times their room temperature values at 4 K, whereas the moduli of glasses, also
282 Chapter 4
8. FATIGUE STRENGTH
A simple reverse bending test is the usual method used to measure fatigue strength.
Fatigue strength is defined as the stress required to cause failure after a certain num-
ber of bending cycles and is given as sf.
Figure 4.27 Young’s modulus of several metals at low temperatures. (1) 2024 T4 aluminum;
(2) beryllium copper; (3) K Monel; (4) titanium; (5) 304 stainless steel; (6) C1020 carbon steel;
(7) 9% nickel steel; (8) Teflon.
Metals exhibiting very small change in Metals exhibiting very large change in
modulus of elasticity from 300 to 20 K modulus of elasticity from 300 to 20 K
Aluminum Aluminum
356 1100
7075 2024
Cobalt
Elgiloy
Copper
Berylco 25
Pure coppera
Iron Iron
Vascojet 1000a 17–7 pH
2800 (9% Ni) 304, 310
4340a 321, 347
Nickel Nickel
Inconel K Monel
Inconel Xa Pure nickel
Titanium Titanium
A-110-ATa C-210-AV
B-120-VCAa
a
No significant temperature dependence.
284 Chapter 4
Figure 4.28 Fatigue strength of several materials at 106 cycles. (1) 2024 T4 aluminum; (2)
beryllium copper; (3) K Monel; (4) titanium; (5) 304 stainless steel; (6) C1020 carbon steel;
(7) 9% nickel steel; (8) Teflon.
Some materials, such as carbon steel, have the property that fatigue failure will
not occur if the stress is maintained below a certain value, called the endurance limit
se no matter how many cycles have elapsed.
Fatigue strength data at low temperatures are not as common as data on yield
strength and ultimate tensile strength, because the fatigue tests are more time-
consuming and hence more expensive to perform. As one would expect, however,
the data that have been reported all show that fatigue strength increases as tempera-
ture decreases, in the same manner as the yield and ultimate tensile strength.
Fortunately, for aluminum alloys at least, it has been found that the ratio of
fatigue strength to ultimate strength remains fairly constant as the temperature is
lowered. Therefore, the fatigue strength of these alloys varies with temperature in
the same manner that the ultimate strength varies with temperature. This fact may
be used in estimating the fatigue strength for nonferrous materials at cryogenic
temperatures if no fatigue data are available.
Representative data of fatigue strength of several materials at low temperature
are shown in Fig. 4.28.
Young’s modulus is 5–20% greater at the lower temperatures, and fatigue properties
are also improved at the lower temperatures (with the exception of 2024 T4 alumi-
num). Annealing of these metals and alloys can affect both the ultimate and yield
strengths.
The body-centered cubic (bcc) metals and alloys are usually undesirable for
low-temperature construction. This class of metals includes iron, the martensitic
steels (low-carbon and the 400 series of stainless steels), molybdenum, and niobium.
If not already brittle at room temperature, these materials exhibit a ductile-to-brittle
transition at low temperatures. Working of some steels, in particular, can induce the
austenite-to-martensite transition.
The hexagonal close-packed (hcp) metals exhibit structural properties inter-
mediate between those of the fcc and bcc metals. For example, zinc suffers a duc-
tile-to-brittle transition, whereas zirconium and pure titanium do not. Titanium
and its alloys, having an hcp structure, remain reasonably ductile at low tempera-
tures and have been used for many applications where weight reduction and reduced
heat leakage through the material have been important. Small impurities of O, N, H,
and C can have a detrimental effect on the low-temperature ductility properties of Ti
and its alloys.
Plastics also increase in strength as the temperature is decreased, but this is
accompanied by a rapid decrease in elongation in a tensile test and a decrease in
impact resistance as the temperature is lowered. The glass-reinforced plastics also
have high strength-to-weight and strength-to-normal conductivity ratios. Conver-
sely, all elastomers become brittle at low temperatures. Nevertheless, many of these
materials, including rubber, Mylar, and nylon, can be used for static seal gaskets
provided they are highly compressed at room temperature prior to cooling.
The strength of glass under constant loading also increases with decreases in
temperature. Since failure occurs at a lower stress when the glass has surface defects,
the strength can be improved by tempering the surface.
See Figs. 4.29–4.32 for the temperature dependence of several important
module: E.
s ¼ pr=t ð4:17Þ
where, p is the internal pressure, r the radius and, t the thickness of cylinder.
286 Chapter 4
Figure 4.29 Temperature variation of (E) Young’s modulus for six metals.
Figure 4.30 Temperature variation of bulk modulus (B)(reciprocal compressibility) for six
metals.
Figure 4.31 Temperature variation of shear modulus (G) for six metals.
Figure 4.32 Temperature variation of Poisson’s ratio (v) for six metals.
288 Chapter 4
s ¼ Pr=2t
The ASME Code that transforms this equation for practical use requires that
the minimum thickness for spherical shells, hemispherical heads, elliptical heads, or
ASME torispherical heads be determined from the equation
PDK PDo K
th ¼ ¼ ð4:19Þ
2sa ew 0:2P 2sa ew þ 2PðK 0:1Þ
where D is the inside diameter of the spherical vessel or hemispherical head, the
inside major diameter for an elliptical head, or 2 (crown radius) for the ASME
torispherical head. The value of the constant K is given by
1
K¼
6½2 þ ðD=D1 Þ2
where D1 is the minor diameter of the elliptical head. For the ASME torispherical
head, K ¼ 0.885.
The outer shell, which has a high vacuum on one side and atmospheric pressure
on the other, is subject to the failure mode of external pressure collapse. A
thin-walled external shell will collapse under an externally applied pressure at stress
values much lower than the yield strength for the material.
2Eðt=Do Þ3
Pc ¼ ð4:21Þ
1 m2
where E is the Young’s modulus of shell material, t the shell thickness Do the outside
diameter of shell and m the Poisson’s ratio for shell material.
A long cylinder in this analysis has an L=D such that
2:42Eðt=Do Þ5=2
Pc ¼ ð4:23Þ
ð1 m2 Þ3=4 ½L=Do 0:45ðt=Do Þ1=2
where L is the distance between stiffening rings for the outer shell.
The heads that enclose the ends of the outer shell cylinder are also subject to
collapsing pressure. The critical pressure for a hemispherical, elliptical, or torisphe-
rical head (or for a spherical vessel) is given by
0:5Eðth =Ro Þ2
Pc ¼ ð4:24Þ
½3ð1 m2 Þ1=2
where th is the thickness of the head and Ro is the outside radius of the spherical
head, or the equivalent radius for the elliptical head, or the crown radius of the
torispherical head. The equivalent radius for elliptical heads is given by Ro ¼ K1D,
where D is the major diameter and K1 is the shape factor.
Stresses on the piping needed to fill and empty cryogenic containers must also
be taken into account.
The minimum wall thickness for piping under internal pressure is given by the
ASA Code as
PDo
t¼
2sa þ 0:8P
290 Chapter 4
where P is the design pressure, Do the outside diameter of pipe, and sa the allowable
stress of pipe material.
For piping subjected to external pressure, Eq. (4.21) may be used.
In addition to pressure-induced stress, thermal stresses on the piping are often
a serious consideration. When the deformation arising from change of temperature is
prevented, temperature stresses arise that are proportional to the amount of defor-
mation that is prevented. Let a be the coefficient of expansion per degree of tempera-
ture, l1 the length of bar at temperature T1, and l2 the length at temperature T2. Then
l2 ¼ l1 ½1 þ aðT2 T1 Þ
If, subsequently, the structured member is cooled to a temperature T1, the propor-
tionate deformation is s ¼ a(T2 T1), and the corresponding unit stress is
s ¼ E a xðT2 T1 Þ:
Stress:
s ¼ EE ¼ bDTE
¼ ð1:656 107 Þð255:6 77:8Þð1:86 1013 Þ ð4:26Þ
¼ 5:48 108 Pa
a. Stress < ultimate strength of 1.17 109 Pa. Therefore, material will not
rupture.
b. Stress > yield strength of 2.76 108 Pa. Therefore, permanent set will
occur. Not satisfactory.
c. Working stress ¼ (1=5)(1.17 109 Pa) ¼ 2.34 108 Pa.
Stress 5.48 108 Pa > 2.34 108 Pa.
Stress exceeds working stress. Not satisfactory.
d. No. Aluminum and brass are both fcc metals and should serve equally
well, all other things being the same.
292 Chapter 4
Allowable Stress
Carbon steel (for outer shell only) SA-285 Grade C 94.8 13,750
SA-299 129.2 18,750
SA-442 Grade 55 94.8 13,750
SA-516 Grade 60 103.4 15,000
Low-alloy steel SA-202 Grade B 146.5 21,250
SA-353-B (9% Ni) 163.7 23,750
SA-203 Grade E 120.6 17,500
SA-410 103.4 15,000
Stainless steel SA-240 (304) 129.2 18,750
SA-240 (304L) 120.6 17,500
SA-240 (316) 129.2 18,750
SA-240 (410) 112.0 16,250
Aluminum SB-209 (1100-0) 16.2 2,350
SB-209 (3004-0) 37.9 5,500
SB-209 (5083-0) 68.9 10,000
SB-209 (6061-T4) 41.4 6,000
Copper SB-11 46.2 6,700
SB-169 (annealed) 86.2 12,500
Nickel alloys (annealed) SB-127 (Monel) 128.2 18,600
SB-168 137.9 20,000
One reason annealed metals are more ductile than work-hardened ones is that work
hardening increases the grain size.
For large storage vessels, 9% nickel steels are commonly used with higher boil-
ing cryogens (T > 75 K), whereas aluminum alloys and austenitic steels can generally
be used over the entire range of liquid temperatures. To speed up cooldown time and
lower the cost of the cooldown of the container, the thickness of the inner shell is
kept as small as possible and is only designed to withstand the maximum internal
pressure and the bending forces. Section VIII of the ASME Boiler and Pressure
Vessel Code provides all the necessary design equations and should be adhered to
for safety reasons when designing vessels of capacities greater than 0.2 m3.
Thermal
Young’s Shear Thermal expansion Specific
Temperature Density modulus modulus Poisson’s conductivity (mean) heat
( F) (lb=in3) (104 psi) (106 psi) ratio (Btu=ft=h= F) (106= F) (101Btu=lb= F)
3.5Ni
70 0.284 29.6 11.5 0.282 20 6.6 108
150 – 30.5 11.9 0.281 17 5.7 84
5Ni
70 0.282 28.7 11.2 0.283 18 6.6 108
260 – 30.1 11.8 0.277 12 5.2 60
320 – 30.3 11.8 0.277 9.2 4.9 36
9Ni
70 0.283 28.3 10.7 0.286 16 6.6 108
260 – 29.6 11.2 0.281 10 5.2 60
320 – 29.8 11.3 0.280 7.5 4.9 36
36Ni (Invar)
70 0.282 22.1 8.08 0.284 8.0 0.65
320 – 20.4 7.37 0.307 3.6 1.00
452 – 20.6 7.33 0.332 0.18 0.77
294 Chapter 4
relative values for the properties of nickel alloys at both room and cryogenic
temperatures.
Commercial nickel and nickel alloys are available in a wide range of wrought
and cast grades; however, considerably fewer casting grades are available. Wrought
alloys tend to be better known by trade names such as Monel, Hastelloy, Inconel,
and Incoloy. There is also a class of superalloys that are nickel-based, strengthened
by intermetallic compound precipitation in a face-centered cubic matrix.
Most wrought nickel alloys can be hot- and cold-worked, machined, and
welded successfully. The casting alloys can be machined or ground, and many can
be welded or brazed. Nearly any shape that can be forged in steel can be forged
in nickel and nickel alloys. However, because nickel work-hardens easily, severe
cold-forming operations require frequent intermediate annealing to restore soft
temper. Annealed cold-rolled sheet is best for spinning and other manual work.
296 Chapter 4
304 154 87 90
304L 118 67 67
310 142 89 86
347 120 66 57
320 F, and 425 F are shown in Table 4.14. The toughness of these four steels
decreases somewhat as the temperature is decreased from room temperature to
320 F, but on the further decrease to 425 F the toughness remains about the
same. It should be noted that the toughness of these steels is excellent at 425 F
and that the toughness of Types 304 and 310 is significantly higher than those of
Types 304L and 347. Shown in Table 4.15 are the results of impact tests on the same
four steels at 320 and 425 F for various thickness of plate. These results demon-
strate that the toughness of the austenitic stainless steels is not markedly affected by
plate thickness.
Table 4.16 shows the results of impact tests done on Type 304 stainless steel
after up to 1 year of exposure at 320 F. These results indicate that Type 304 is very
0 85.0
30 min 72.5
6 months 73.0
12 months 75.0
stable and does not exhibit any marked degradation of toughness due to long-time
exposure at 320 F.
298 Chapter 4
5
Transport Properties of Solids
1. THERMAL PROPERTIES
302 Chapter 5
later. In addition, for solids (and liquids), the difference between CP and CV is
usually less than 5% at room temperature. Accordingly, we shall limit the present
discussion to CV, to lattice effects only, and to electrical insulators. The specific heat
contribution of free electrons in a metal will be taken up after the lattice effects have
been explored.
(a) Models of Lattice Specific Heat. The various models of lattice specific heat
will be discussed in the order of their historical development, which is also more or
less the order of their exactness.
Dulong and Petit. In 1911, Dulong and Petit observed that the heat capacity
of many solids, both metals and nonmetals, was very nearly independent of tempera-
ture near room temperature, and furthermore that it had a value of 6 cal=(mol K )
regardless of the substance. This experimentally observed fact can be explained if
it is assumed that each lattice point can absorb energy in the same way as every other
lattice point (the equipartition of energy). In this model, the atoms are assumed to
oscillate as independent classical particles. Thus the total thermal energy of 1 mol
of material would be
kT
U ¼ ð3Þð2ÞNo ¼ 3No kT ð5:2Þ
2
where 3 is the number of degrees of freedom, 2 arises from the fact that the bound
oscillator has equal amounts of both kinetic and potential energy, No is Avogadro’s
number, the number of particles per mole, and (kT=2) is the thermal energy of each
degree of freedom.
From Eq. (5.1), then,
@U
CV ¼ ¼ 3No k ¼ 6 cal=ðmol KÞ
@T
or
2 !
yE eyE =T
CV 1 ¼ 3R ð5:4Þ
T ðeyE =T 1Þ2
where Nok has been replaced by its equivalent R, the general gas constant, and where
yE ¼ hn=k, which has the dimensions of temperature and is the Einstein characteristic
temperature. The characteristic temperature yE is the parameter that permits varia-
tion from one material to another. It permits plotting all materials on one curve by
replacing T by T=yE as the abscissa. Expression (5.4) fits experimental data quite well
for all materials except that at low temperatures it drops below the experimental
data. At intermediate and high temperatures, the fit is good, and Eq. (5.4) provides
the approach to a limiting Dulong–Petit value of 3R at high T. Although now super-
seded by more exact models, the Einstein model included the most important effect,
the quantization.
Nernst and Lindemann. The Einstein model was a considerable improvement
over that of Dulong and Petit, but it still did not have the correct low-temperature
shape. Real atoms in the lattice do not vibrate at only one frequency, but are coupled
to one another in a complex way that depends upon both the lattice type and the
nature of the interatomic forces.
Nernst and Lindemann improved upon the Einstein model by assuming that
the atoms vibrate at two frequencies, one equal to the Einstein frequency n, and
one at half that value, n=2. This treatment, while an improvement, was still a
considerable oversimplification and was on shaky theoretical grounds because it
required the introduction of a half quantum of energy.
Debye. Debye made a major advance in the theory of heat capacity at low
temperatures by treating a solid as an infinite elastic continuum and considering
the excitement of all possible standing waves in the material. These range from
the acoustic vibrations up to a limiting frequency, n m. Instead of a system of simple
304 Chapter 5
harmonic oscillators all vibrating at only one frequency (Einstein) or two frequencies
(Nernst and Lindemann), Debye derived a parabolic frequency distribution for the
atoms vibrating in the lattice. The Debye model gives the following expression for
the lattice heat capacity per mole (CV1):
3 Z yD =T 4 x 3
T x e dx T T
CV 1 ¼ 9R 2
¼ 3R D ð5:5Þ
yD 0 ðex 1Þ yD yD
where R is the universal gas constant per mole, T is the absolute temperature, yD is
the Debye characteristic temperature, and x is a dimensionless variable defined by
the expression x ¼ hn=kT. In this last expression, h is Planck’s constant, n is the
frequency of vibration, and k is Boltzmann’s constant.
In the abbreviated form of the equation, D(T=yD) stands for the expression in
the integral, the Debye function. The Debye expression for CV vs. T=yD is plotted in
Fig. 5.1.
The important thing to notice about the Debye function is that for a given
substance, the lattice heat capacity is dependent only on a mathematical function
of the ratio of the Debye characteristic temperature yD to the absolute temperature.
This mathematical function is the same for all materials, and only yD changes from
material to material. yD thus appears as a constant characteristic of the material at
hand.
The value of yD is an adjustable parameter for the temperature scale, similar to
yE for the Einstein model. Its value usually lies between 200 and 400 K, although
values occur far outside this range. Values of yD are given in Table 5.1.
The Debye function predicts heat capacity surprisingly well, especially consid-
ering that Debye merely assumed a parabolic form for the vibration spectrum of
every solid without regard for the individual differences to be expected among
different solids.
At high temperatures (T > 2yD), the specific heat given by Eq. (5.5) approaches
a constant value of 3R, which is the Dulong–Petit value, as it should be. At low tem-
peratures (T < yD=12), the Debye function approaches a constant value of
D(0) ¼ 4p4=5; thus, the lattice specific heat at temperatures less than yD=12 may be
expressed as
or
CV ¼ AT 3 ð5:7Þ
where A is a constant.
From Eq. (5.6), we see that the lattice contribution to the specific heat of solids
varies as the third power of the absolute temperature at very low temperatures. This
behavior is useful in separating the effects of electrons from those of the lattice for
electrical conductors, and this will be used later.
The Debye temperature, yD, is not only important to estimating heat capacity
but is also a useful property in itself. yD can be calculated directly from first
principles by the theoretical expression
hn a 3N 1=3
yD ¼ ð5:8Þ
k 4pV
305
© 2005 by Marcel Dekker
5367-4 Flynn Ch05 R2 090904
306 Chapter 5
where h is the Planck’s constant; n a the speed of sound in the solid; k the
Boltzmann’s constant and N=V is the number of atoms per unit volume for the solid.
In practice, the Debye temperature is usually determined by selecting the value
that makes the theoretical specific heat curve fit the experimental specific heat curve
as closely as possible. In such cases, yD has lost its physical significance as a measure
of limiting frequency and has become an ‘‘effective’’ characteristic temperature.
A modification of the Debye theory is sometimes used to calculate heat
capacity in a semiempirical way. Essentially, yD is treated as a function of tempera-
ture rather than as a characteristic constant. A few heat capacities are measured, and
yD is allowed to vary with temperature to get the best fit between the Debye function
and experiment. Usually, the resulting variation of yD with temperature is not too
great for metals. Plastics, however, are characterized by large changes in their intera-
tomic force constants with temperature (equivalent to large changes in their elastic
moduli). As a result, yD for these materials increases greatly with decreasing
temperature.
The situation for compounds and alloys is more complicated. Near room
temperature, the specific heat of a compound or alloy can be approximated quite
well by linear combination of the specific heats of the constituent elements (the
Kopp-Joule rule). This is because the temperature is sufficiently high that all modes
of lattice vibration are fully excited (except for a few hard substances of low atomic
weight such as diamond, BeO, B2O3), and the gram atomic heat capacity is about the
same for all the constituent elements (the Dulong-Petit value).
This procedure is not proper at low temperatures because there the specific heat
is determined by the arrangement and bonding of the lattice, and for a compound
this bears no relation to the lattice characteristics of the component elements. This
is seen at once by comparing the solid properties of an alkali halide with those of
the alkali metal and the halogen as separate entities. Nevertheless the Debye function
is found to represent simple compounds that form lattices with only one type of
bonding (such as alkali halides, other binary salts, and some oxides). The yD value
in this case bears no relation to the yD value of the constituent elements.
Other Models. Although the model of Debye gives good agreement with
observation in many cases, deviations are observed. The cause seems to be that
the detailed crystal structure is not introduced. Much further work has been done,
beginning with Born and von Karman and Blackman and continued by many others.
These models attempt to express the detailed situation of actual crystal lattices. This
leads immediately to dispersion; that is, the velocity of sound depends upon the
frequency. In all these models, the phonon (quantized sound wave or lattice vibra-
tion) picture is retained, and it has become one of the well-established particles.
The frequency spectra differ from model to model, and each crystal leads to a differ-
ent frequency spectrum. These other models are very important to achieving a
detailed understanding of the crystal lattice. For process design purposes, however,
the Debye function and yD yield adequate estimations. As usual, though, good
experimental data are much to be preferred.
The specific heats of several materials used in low-temperature construction are
shown in Fig. 5.2 as a function of temperature. Values for other materials are tabu-
lated in several excellent references (Corruccini and Gniewek, 1960; Johnson, 1960;
Landolt-Bornstein, 1961), and a few values are abstracted in Table 5.2 (Gopal, 1966).
Example 5.1. Determine the lattice specific heat of chromium at 200 K as given by
the Debye function.
Temperature (K) Al Mg Cu Ni a-Mn a-Fe g-Fe Cr 18–l8 SSa Monel Fused silica Pyrex Teflon
Chapter 5
© 2005 by Marcel Dekker
5367-4 Flynn Ch05 R2 090904
different. The additive rule will fail for alloys having crystal structures and elastic
constants widely different from those of either of the major components.
It was noted that the Debye and other functions give CV, whereas one desires
CP for most purposes. The difference between CP and CV is given by several
thermodynamic expressions, one of which is
CP CV ¼ a2 VT=bT ð5:10Þ
in which V is the molar volume, a is (1=V )(dV=dT )P, which is equal to three times
the coefficient of linear expansion, and bT is –(1=V )(dV=dP)T. However, these
coefficients are not always known. Since the difference between the two specific heats
is small, 1–10% of CP, it can often be neglected or calculated approximately from the
expression
CP CV ¼ ACV 2 T ð5:11Þ
in which A ¼ a2V=bTCV 3 is regarded as a constant and is obtained from data at
normal temperatures.
If, to go one step further, the data are not available for calculating A at any
temperature from the above formula, A can be obtained to a less satisfactory degree
of accuracy from a formula by Nernst and Lindemann
A ¼ 0:0214=Tm ð5:12Þ
for specific heat in cal=(mol K) and Tm in K. Finally, it may be admitted that if CV
itself has to be estimated by an indirect method, the accuracy of the resulting heat
capacities will probably not be such as to justify taking CP – CV into account at all.
Often the property actually desired (as in heat balances) is the enthalpy differ-
ence over a large temperature interval, and this difference is fortunately much less
sensitive to estimation errors than is the specific heat.
310 Chapter 5
so that at 26 K, where T ¼ yD=12, the lattice specific heat is 0.275 cal=(g-atom) and
the electronic specific heat is only 0.0046 cal=(g-atom). At 4 K, the two contributions
to the specific heat are 0.001 and 0.007 cal=(g-atom), respectively, and at 1 K they are
0.0000155 and 0.000178 cal=(g-atom). Thus the electronic contribution, which is of
considerable theoretical interest, can most accurately be determined by measure-
ments at temperatures below about 4 K. This is a typical example of the contribu-
tions to physical theory that can be made as a result of low-temperature research.
The coefficient g in Eq. (5.13 ) is given by
where a is the number of valence electrons per atom; me the electron mass; M the
atomic weight of material; R the specific gas constant for material; T the
absolute temperature; h the Planck’s constant; No the Avogadro’s number and
N=V is the number of free electrons per unit volume.
Example 5.2. Determine the ratio of the electronic contribution to the specific heat,
Cve, to the lattice specific heat (as determined by the Debye function) for copper at
(a) 5.6 K, (b) 0.56 K, and (c) 0.056 K.
a. At 5.6 K:
J
CV e ¼ ð1:13 102 Þ ð5:6 KÞ ¼ 0:0645 J=ðkg KÞ
kg K2
T=yD ¼ 5:6=309:4 ¼ 0:0181
3
T 3 cal
CV ¼ 233:8 R ¼ 2:75 10 ð4184Þ ¼ 11:5 J=ðkg KÞ
yD gK
CV e =CV ¼ 5:61 103
b. At 0.56 K, by similar calculations,
CV e =CV ¼ 0:561
c. At 0.056 K,
CV e =CV ¼ 56:1
312 Chapter 5
For the present, we are concerned with structural materials, so this overview is
limited to dielectrics (insulators) and metals. Furthermore, since dielectrics have only
one heat transport mechanism (phonons) and metals have two (phonons plus
electrons), we shall begin by examining heat transport in dielectrics.
one-third of the phonons are traveling in the desired direction, they travel at a
velocity U, and they travel a distance L before being diverted. To arrive at a knowl-
edge of how the thermal conductivity should vary with temperature, it is necessary to
consider the various mechanisms by which phonon scattering takes place and the
dependence upon temperature of each mechanism. This amounts to examining the
temperature dependence CV, U, and L (Eq. (5.16)).
The specific heat CV is the lattice specific heat of the solid, as already discussed.
Its variation with T is plotted in Fig. 5.3b. The average velocity U of the phonons is
the mean of the velocity of sound. U varies only slightly with temperature as plotted
in Fig. 5.3c. The phonon gas differs from a real gas in that the number of particles
varies with the temperature, increasing in number as temperature is increased.
At high temperatures, the large number of phonons leads to collisions between
phonons. Then, as T increases, L decreases, as shown in Fig. 5.3d.
At low temperatures, phonon–phonon collisions disappear and L is determined
by the distance between imperfections in the crystal. The phonons collide with
imperfections such as impurity atoms, dislocations, intercrystalline boundaries, or
finally, the specimen boundaries. The distance between these does not depend on
temperature, so L becomes constant at low T.
The general shape of the curve of K vs. T can now be readily expressed in terms
of the other curves. U is almost constant for all T, so it has negligible effect on the
shape. At low T, L is constant. Thus the K curve has the shape of the one remaining
variable, CV. At high T, CV is constant, so the K curve has the shape of the L curve,
which drops off with temperature. At intermediate temperature, a broad maximum
occurs in the transition between the high- and low-temperature behavior.
These considerations lead to the general temperature dependence of Kt shown
in Fig. 5.3a, which is experimentally observed also. For all materials, it starts at zero
for absolute zero temperature, rises to a maximum, and then falls asymptotically
toward zero as the temperature becomes very high.
Materials differ in the height of the maximum and its position in temperature,
as can be seen in Fig. 5.4, which shows the thermal conductivity curves for several
samples. The variability of thermal conductivity with impurity and cold work and
the suppression of the maximum in Kt is shown in Fig. 5.3. Note especially the ther-
mal conductivity of the dielectric, corundum. Its peak conductivity, around 37 K, is
almost 15 times that of copper at room temperature, making it an excellent thermal
314 Chapter 5
This expression holds when T is equal to or greater than yD and gives a 1=T
dependence on temperature.
When T is considerably less than yD, i.e., below about yD=10, then
3
T yD
Kt ¼ K0 exp
yD bT
path increases until it reaches the order of the dimensions of the crystal. Phonons are
then reflected from its surfaces. With L and U now constant, Kt depends on the
dimensions of the crystal and varies as CV with temperature.
We have already seen that the heat capacity of a dielectric varies as T 3 at very
low temperatures, and so we also expect a T 3 temperature dependence for Kt at the
lowest temperatures. This dependence has in fact been observed experimentally.
The type of thermal conductivity curve suggested by all of these considerations
is exemplified by that of corundum in Fig. 5.4. Above temperatures corresponding to
about yD=10, phonon–phonon collisions govern. These are decreased at lower
temperatures to cause a maximum near yD=20, and below this the size of the crystal
limits the mean free path and the T 3 dependence causes a drop to zero at 0 K.
The phonon interaction and boundary scattering process establish the general
shape of the conductivity curve for a perfect crystalline dielectric. Further scattering
may occur at internal lattice imperfections such as vacancies, interstitial atoms,
impurity atoms (among which may be included isotopes of the primary constituents),
and dislocations. These added resistances have the effect of lowering and flattening
the conductivity curve and, in the limit of a glassy or disordered solid, lead to a curve
without a maximum.
Kt ¼ Ke þ Kg ð5:18Þ
Ke ¼ ð1=3ÞCe Ue Le ð5:19Þ
where Ce is the electronic heat capacity; Ue the mean velocity of the electron and Le is
the mean free path of electrons between collisions with obstacles.
Although Ce is very small, Ue is the speed of light, not the speed of sound, as is
the case with phonons. Furthermore, electrons effectively scatter phonons, reducing
the contribution of Kg still further. In fact, the striking difference between the very
high phonon conduction of dielectrics and its very low value in metals is due to
the effectiveness of the free electrons in scattering phonons. Accordingly, Ke is much
greater than Kg for most pure metals.
The details of thermal conductivity in metals can be seen by examining the
temperature dependence of the electron heat capacity Ce and the lattice heat capacity
Cl. The former is given by
k2
Ce ¼ p2 ZNo T ð5:20Þ
m v2
316 Chapter 5
where m is the effective electron mass: n is the electron velocity at the energy of the
Fermi face: Z the number of free electrons per atom and No is Avogadro’s number
(number of atoms per mole). The coefficient of T in the preceding equation is a
constant of the material, and it is convenient to give it the symbol g. Then
Ce ¼ gT ð5:21Þ
that is, the electron specific heat varies linearly with T. Ordinarily, gT is small com-
pared to the lattice specific heat. (See Sec. 1.1.1.) However, it becomes appreciable at
very high T, where the lattice CV remains constant at a value 3R. The electron CV
also becomes appreciable at low T, where the Debye model predicts that the lattice
CV varies as (T=yD)3. Then the total CV is given by
3
T
CV ¼ gT þ 234R ð5:22Þ
yD
For copper, at 1 K, the first term is roughly four times the magnitude of the
second term in this equation.
In dealing with the thermal conductivity due to the electrons, CV0 , the specific
heat per unit volume, is needed. This differs from CV, the specific heat per mole,
by having No replaced by n, the number of electrons per unit volume.
For metals, the thermal conductivity is almost entirely the result of energy
carried by the electrons rather than by the phonons. The two heat fluxes are additive,
but the phonon conductivity becomes negligible because collisions of phonons with
electrons drastically reduce the phonon mean free path in metals.
In treating the thermal conductivity due to electrons, the solid is considered to
contain an electron gas. The gas thermal conductivity expression may again be
applied, that is
K ¼ ð1=3ÞCV0 UL ð5:23Þ
in much the same way as for phonons. At high T, the electrons collide with the
phonons, so L decreases as the number of phonons increases. At low temperature,
the number of phonons is small and the electron mean free path is limited by dis-
tances between crystal imperfections. These are impurities, dislocations, intercrystal-
line boundaries, or, finally, the specimen walls. Looking again at the product of CV0 ,
V, and L, the quantity V is constant, so K depends on CV0 and L. At low T, L is
constant, so K follows CV0 and varies linearly with T. At high T, CV0 still varies linearly
with T, but L begins to fall off so rapidly that the product begins to drop off. At
intermediate T, a broad maximum occurs.
Figures 5.6 and 5.7 show the thermal conductivity curves for several solids,
with the highest curve having the highest purity. It will be noted that the metals
of high purity exhibit a maximum of conductivity at low temperatures that in some
cases is many times the room temperature value. Moreover, the conductivities of
these metals approach a room temperature value that is almost temperature-
independent. The height (or existence) of the maximum is quite sensitive to certain
slight impurities, although it will be noted that some of the commercial coppers
exhibit maxima.
Table 5.3 gives the thermal conductivity and thermal conductivity integral for a
few common materials of construction (Johnson, 1960).
The variability of thermal conductivity with alloying and cold working and the
suppression of the maximum are shown clearly in Figs. 5.5 and 5.6. Inconel, Monel,
318 Chapter 5
and stainless steel are structural alloys that also exhibit these properties and are thus
useful in cryogenic service, which requires low thermal conductivity over the entire
temperature range.
Alloys and highly cold-worked metals represent cases in which scattering from
imperfections in the lattice is so great that this effect predominates over phonon
scattering even at room temperature. The resulting flattening of the thermal conduc-
tivity curve is illustrated in Fig. 5.4 for nickel–chromium stainless steel and brass.
Furthermore, pure metals having a small number of conduction electrons per atom
(e.g., antimony and bismuth) and semiconductors (e.g., germanium) are all metals
for which phonon thermal conduction can no longer be safely neglected.
Experimental thermal conductivity values for most materials used in
low-temperature design are readily available in the literature (Powell, 1963;
Touloukian, 1964; Powell et al., 1966).
Example 5.3. In support members for cryogenic storage vessels and other equip-
ment, it is desirable to minimize the heat transfer through the members but maximize
Temperature RT RT RT RT RT RT
(K) Kt 4:2 K Kt dT Kt 4:2 K Kt dT Kt 4:2 K Kt dT Kt 4:2 K Kt dT Kt 4:2 K Kt dT Kt 4:2 K Kt dT
319
© 2005 by Marcel Dekker
5367-4 Flynn Ch05 R2 090904
320 Chapter 5
the strength of the members at the same time. An important parameter in this situa-
tion is the strength-to-conductivity ratio. Determine this ratio, Sy=Kt [in (N K)=(W
m)] at 22.2 K for
a. 2024 T4 aluminum
b. 304 stainless steel
c. Monel
d. Beryllium copper
e. Teflon
Using Fig. 4.21 for Sy values, and Fig. 5.7 for Kt values:
a. 2024 T4 aluminum
Sy ¼ 7.24 108 N=m2
Kt ¼ 17.3 W=(m K)
Sy=Kt ¼ 4.18 107 N K=(W m)
b. 304 stainless steel
Sy ¼ 1.32 109 N=m2
Kt ¼ 2.77 W=(m K)
Sy=Kt ¼ 4.77 108 N K=(W m)
c. Monel
Sy ¼ 1.03 109 N=m2
Kt ¼ 7.61 W=(m K)
Sy=Kt ¼ 1.35 108 N K=(W m)
d. Beryllium copper
Sy ¼ 8.96 108 N=m2
Kt ¼ 9.52 W=(m K)
Sy=Kt ¼ 9.41 107 N K=(W m)
e. Teflon
Sy ¼ 1.38 108 N=m2
Kt ¼ 0.138 W=(m K)
Sy=Kt ¼ 1.0 109 N K=(W m)
all) in the copper. Thus they are localized at grain boundaries. This condition is also
responsible for the free-machining characteristic.
The extreme sensitivity of the low-temperature conductivity of nearly pure
metals to impurities and cold working is in strong contrast to the specific heat and
expansivity, both of which are somewhat sensitive to the type of lattice structure
and to the interatomic forces but relatively insensitive to local imperfections of the
lattice. The variability shown by coppers is demonstrated by other metals. This
serves to illustrate that where conductivity is concerned, there are no ‘‘pure’’ metals
but only alloys of various degrees of dilution. Fortunately for the cryogenics
designer, the only other widely used cryogenic metal (aluminum) is available in only
one nominally pure commercial grade.
As noted earlier, the lattice conduction of pure metals and dilute metal alloys is
generally small, so ke can be taken as the total conduction. This principle has been
used as the basis for a generalized low temperature equation for thermal conductivity
as represented in Fig. 5.9, where k is equal to k=km, T is equal to T=Tm, and Tm is
the temperature at which k is a maximum, km. In analogy to reduced equations of
state for thermodynamic properties of fluids, this correlation is referred to as the
principle of corresponding states. The reducing parameters chosen are the tempera-
ture and thermal conductivity at the maximum in conductivity. The results of this
correlation are illustrated in Fig. 5.9. The figure represents 1002 data points for a
total of 22 metals (83 specimens) in varying states of impurity and imperfection
content. The standard deviation of the data is only 3.2%.
322 Chapter 5
Figure 5.9 Reduced thermal conductivity vs. reduced temperature of 22 metallic elements.
The k is k=km where km is the maximum value of k and T is T=Tm where Tm is the tempera-
ture corresponding to km.
were the case, there would be no thermal expansion. Instead, the intermolecular
potential energy curve (forces of attraction and repulsion) is not symmetrical.
The general shape of the intermolecular potential, as illustrated in Fig. 5.10, is
asymmetric. This results from a repulsive force that increases very rapidly as the
two atoms approach each other closely (one atom is considered at rest at the origin
in Fig. 5.10) plus an attractive force that is not such a strong function of distance.
Their algebraic sum provides a curve having a minimum at r0, the classical
equilibrium separation of the atoms at 0 K. Increasing the temperature from 0 K
corresponds to raising the energy of the system slightly, and now thermal vibrations
take place between two higher values of r. But because of the asymmetry of the
potential well, the average distance no longer corresponds to r0; it is slightly greater
and results in the observed thermal expansion of materials. The mean spacing of the
atoms increases with temperature as the energy (temperature) of the material
increases. Thus, the coefficient of thermal expansion increases as temperature is
increased.
Both the specific heat and the coefficient of thermal expansion arise from the
intermolecular potential, and accordingly these two properties are related. Our
present understanding of this relation is largely due to Gruneisen (1926), who
developed an equation of state for solids based on the lattice dynamics of Debye.
He showed that the dimensionless ratio g should be a universal constant:
g ¼ aV =bT CV ð5:24Þ
or
g ¼ aV =bS CV ð5:25Þ
where
1 @V
a¼ volume coefficient of thermal expansion
V @T P
1 @V
bT ¼ isothermal compressibility
V @P T
1 @V
bS ¼ adiabatic compressibility
V @T S
324 Chapter 5
Since both CV and a depend on the intermolecular potential of the lattice, they
are related to a first approximation by the Debye temperature yD
V dyD
g¼ ð5:27Þ
yD dV
yD ¼ 309; MW ¼ 63:54
T=yD ¼ 22=309 ¼ 0:0712 < 1=12
3
T
CV ¼ 233:8R
yD
1:987
CV ¼ ð233:8Þ ð0:0712Þ3 ¼ 0:00264 cal=ðg KÞ ¼ 11:04 J=ðkg KÞ
63:54
Material g Material g
Li 1.17 Mo 1.57
Na 1.25 W 1.62
K 1.34 Fe 1.60
Mg 1.51 Ni 1.88
Al 2.17 Co 1.87
Sn 2.14 Pd 2.23
Pb 2.73 Pt 2.54
Sb 0.92 NaCl 1.63
Bi 1.14 KCl 1.60
Cu 1.96 KBr 1.68
Au 3.03 FeS2 1.47
Zn 2.01 PbS 1.94
Cd 2.19 Ta 1.75
Figure 5.11 Coefficient of linear expansion for several metals as a function of temperature
in K.
326 Chapter 5
by the addition of filler material of low expansivity such as glass fiber, silica,
alumina, and asbestos.
As in the case of heat capacity, anomalies can also occur in the thermal expan-
sivity of materials. Various transformations that occur in the solid state can alter the
interatomic forces of the material and lead to changes in its overall dimensions. For
example, the ferromagnetic Curie point of dysprosium at low temperatures leads to
an irregularity in its expansivity. Other examples are KH2PO4 with a ferroelectric
Curie point at 122 K, and the ‘‘glass’’ transitions of soft rubbers. In fact, the mea-
surement of thermal expansion is a useful method by which solid-state transitions
can be investigated. Tables 5.5–5.7 give values of thermal expansion for some
structural materials at low temperatures.
Example 5.5. An aluminum distillation column 20 m high is placed in an air separa-
tion service. How much will the column shrink as it is cooled from 300 K to an aver-
age temperature of 89 K?
Table 5.5 Mean Linear Thermal Expansion of Various Metals at Selected Temperaturesa
1020
Tempera- Low-carbon Stainless Yellow
ture (K) Ni Mg Zn Cu Al steel steel Inconel Monel brass
0 0 0 0 0 0 0 0 0 0 0
20 0 1 1 1 0 0 0 0 0 1
40 1 5 9 4 3 1 0 1 1 4
60 4 12 28 12 9 4 4 5 6 16
80 12 29 57 26 24 10 13 12 15 34
100 23 55 93 45 46 20 32 24 28 58
120 38 87 133 67 71 32 49 38 45 85
140 55 124 176 92 104 47 72 55 64 115
160 74 164 221 120 138 63 97 74 84 147
180 95 208 267 149 175 81 122 95 107 180
200 117 254 314 179 214 101 150 117 130 215
220 140 303 363 210 254 121 179 140 154 249
240 164 353 413 241 295 142 208 163 179 283
260 188 403 465 273 339 164 240 187 207 320
280 213 453 518 305 384 187 272 212 233 357
300 239 503 572 337 418 210 304 238 261 397
a
Net expansion from 0 K to any other temperature in the table, DL=L0 105. DL is the change in length
from 0 K to T. L0 is the length at 0 K.
Source: Corruccini and Gniewek (1961).
been reported to exhibit this characteristic. (This effect for Pyrex glass between 40 C
and 0 C can be seen in Table 5.6.) Some anisotropic materials such as zinc show a
negative expansivity in one crystallographic direction although they have an average
positive expansivity, but for the materials just mentioned, it is the average linear
expansivity that has been reported negative. In 304 stainless steel, this occurrence
is ascribed to the partial transformation of the material at low temperatures to a sec-
ond phase, martensite.
Figure 5.13 shows the range of Thermal expansion values for various
categories of materials.
Estimating radiant heat transfer across an insulating vacuum requires data on the
emissivity, absorptivity, and reflectivity of the surfaces involved. The common defini-
tions for these terms are as follows.
Total normal emissivity (En) is the rate of radiant energy emission normal to the
surface, divided by the corresponding rate from a blackbody.
Total hemispherical emissivity (eh) is the rate of radiant energy emission into a
hemisphere (centered on the normal to the emitting surface), divided by
the corresponding rate from a blackbody.
Total absorptivity (a) is the fraction of energy incident upon a surface that is
absorbed. Quantities an and ah analogous to en and eh may be distin-
guished.
Reflectivity (r) is the fraction of incident energy that is reflected; r ¼ 1 a. The
value of the reflectivity also depends on the angle of incidence.
0 0 0 0 0 0 0 0 0
20 1 30 3 14 10 10 21 27
Chapter 5
© 2005 by Marcel Dekker
Transport Properties of Solids
Table 5.7 Thermal Expansion of Selected Metalsa
Aluminum Aluminum Aluminum Titanium Copper ETP Brass Monel K Stainless 2800 Steel
2024 5083 7039 6A1-4V & OFHC 70=30 &S steel 304 (9% Ni)
Temperature
( R) A B A B A B A B A B A B A B A B A B
329
(Continued)
330
Aluminum Aluminum Aluminum Titanium Copper ETP Brass Monel K Stainless 2800 Steel
2024 5083 7039 6A1-4V & OFHC 70=30 &S steel 304 (9% Ni)
Temperature
( R) A B A B A B A B A B A B A B A B A B
360 228 10.56 242 11.06 244 11.10 99 4.55 179 8.44 201 9.46 146 6.67 169 7.78 128 5.81
380 207 10.79 219 11.31 222 11.37 90 4.64 162 8.56 182 9.60 133 6.82 153 7.92 116 5.96
400 185 10.98 196 11.55 199 11.62 81 4.73 145 8.68 162 9.72 119 6.97 137 8.07 104 6.09
420 163 11.15 173 11.75 175 11.85 71 4.81 127 8.78 143 9.84 105 7.10 121 8.19 92 6.21
440 141 11.30 149 11.94 151 12.06 61 4.89 110 8.87 123 9.95 91 7.22 104 8.32 79 6.32
460 118 11.47 125 12.13 127 12.27 51 4.97 92 8.96 103 10.05 76 7.33 87 8.45 66 6.42
Chapter 5
© 2005 by Marcel Dekker
5367-4 Flynn Ch05 R2 090904
Figure 5.13 Range of thermal expansion values for various categories of materials.
332 Chapter 5
At a fixed temperature, the emissivity must equal the absorptivity. (If these
differed, there could be a net transfer of heat between two surfaces at the same
temperature—a violation of the second law of thermodynamics.) Helpful tabulations
of emissivity values have been prepared by Fulk (1959), and the emissivities of a few
selected materials are shown in Table 5.8.
3. ELECTRICAL PROPERTIES