Ce - Module 2
Ce - Module 2
CRYOGENIC ENGINEERING
−𝑊ሶ −𝑊ሶ
= 𝑦
𝑚ሶ 𝑚ሶ 𝑓
28-May-23 Dr. Kishor Kumar V. V., MED, GCEK 2
GAS LIQUEFACTION SYSTEMS
𝑊ሶ 𝑖 −𝑊ሶ 𝑖 ൗ𝑚ሶ 𝑓
𝐹𝑂𝑀 = =
𝑊 −𝑊ሶ ൗ𝑚ሶ 𝑓
ሶ
TRACE KTU
▪ It is highly impractical to attain this pressure in
liquefaction systems
▪ Thus it is not an ideal process for a practical
system
ሶ
𝑣2 𝑣 2ሶ
𝑄ሶ 𝑛𝑒𝑡 − 𝑊ሶ 𝑛𝑒𝑡 = 𝑚ሶ ℎ + + 𝑔𝑧 − 𝑚ሶ ℎ + + 𝑔𝑧
2 2
𝑜𝑢𝑡𝑙𝑒𝑡𝑠 𝑖𝑛𝑙𝑒𝑡𝑠
TRACE KTU
In liquefaction systems, the kinetic and potential energy changes are much
smaller than the enthalpy changes and these energy terms are neglected
ሶ
𝑄ሶ 𝑛𝑒𝑡 − 𝑊ሶ 𝑛𝑒𝑡 = 𝑚ℎ
ሶ − 𝑚ℎ
ሶ
𝑜𝑢𝑡𝑙𝑒𝑡𝑠 𝑖𝑛𝑙𝑒𝑡𝑠
TRACE KTU
TRACE KTU
−𝑊ሶ 𝑖
𝑚ሶ TRACE KTU
ሶ 𝑠 − ℎ − ℎ =−𝑊ሶ 𝑖
= 𝑇1 𝑠1 − 𝑓 1 𝑓 ሶ 𝑚𝑓
TRACE KTU
▪ The transportation of these gases across the world is done in liquid
state
▪ Gas refrigeration can be used as a precoolant for liquefiers and also in
applications where low temperature gases are required
TRACE KTU
▪ This process is called as refrigeration
𝑄𝐿
𝐶𝑂𝑃 =
𝑄𝐻 − 𝑄𝐿
TRACE KTU
▪ COP represents Watt of cooling effect obtained per Watt
of power input at a particular temperature
▪ Best performance is delivered by a refrigeration system
operating on reversed Carnot cycle
28-May-23 Dr. Kishor Kumar V. V., MED, GCEK 13
Refrigeration-Basics
TRACE KTU
▪ 𝑊𝑝 𝑎𝑛𝑑 𝑊𝑐 𝑎𝑟𝑒 𝑡ℎ𝑒 𝑤𝑜𝑟𝑘 𝑖𝑛𝑝𝑢𝑡 𝑎𝑛𝑑 𝑐𝑜𝑜𝑙𝑖𝑛𝑔 𝑒𝑓𝑓𝑒𝑐𝑡 𝑟𝑒𝑠𝑝𝑒𝑐𝑡𝑖𝑣𝑒𝑙𝑦 𝑖𝑛 𝑊𝑎𝑡𝑡𝑠
▪ 𝐴𝑠 𝑇𝐿 𝑑𝑒𝑐𝑟𝑒𝑎𝑠𝑒𝑠, 𝑡ℎ𝑒 𝐶𝑎𝑟𝑛𝑜𝑡 𝐶𝑂𝑃 𝑑𝑒𝑐𝑟𝑒𝑎𝑠𝑒𝑠
28-May-23 Dr. Kishor Kumar V. V., MED, GCEK 14
Methods of production of low temperature
TRACE KTU
1. Throttling method 2. Heat Exchanger
TRACE KTU
3. Compression/Expansion
28-May-23 Dr. Kishor Kumar V. V., MED, GCEK 16
Refrigerator
▪ Operates in a closed thermodynamic cycle
▪ The rate of mass flow is constant at any point inside
the system
TRACE KTU
▪ The heat is exchanged between the cold end and the
object to be cooled
▪ This cold end heat exchanger can also be used to
liquefy gases
TRACE KTU
in an open thermodynamic cycle
▪ The mass deficit occurring due to the loss of
working fluid is replenished by a make up gas
connection
Joule-Thomson effect
𝑣𝑓 2 𝑣𝑖 2
𝑄𝑛𝑒𝑡 − 𝑊𝑛𝑒𝑡 = 𝑚ሶ 𝑓 ℎ𝑓 + + 𝑔𝑧𝑓 − 𝑚ሶ 𝑖 ℎ𝑖 + + 𝑔𝑧𝑖
2𝑔 2𝑔
neglected
TRACE KTU
▪ The changes in velocities and datum levels are very small and can be
ℎ𝑓 = ℎ𝑖
TRACE KTU
𝜇𝐽𝑇
TRACE KTU
=
𝜕𝑇
𝜕𝑝 ℎ
𝜇𝐽𝑇
>0
TRACE KTU
Effect
Cooling
<0 Heating
=0 No Effect Isenthalpic expansion of a real
gas
𝜕ℎ
𝜕𝑝
𝜕𝑝
𝜕𝑇
TRACE KTU
𝜕𝑇
𝜕ℎ
= −1
𝑑𝑠 =
𝜕𝑠
𝜕𝑇 𝑝
𝑑𝑇 +
𝜕𝑠
𝜕𝑝 𝑇
𝑑𝑝
𝑇 ℎ 𝑝
𝜕𝑠 𝜕𝑠
𝑇𝑑𝑠 = 𝑇 𝑑𝑇 + 𝑇 𝑑𝑝
𝜕𝑇 𝑝
𝜕𝑝 𝑇
𝜕𝑠 𝜕𝑣
𝑇 = 𝑐𝑝 𝑇𝑑𝑠 = 𝑐𝑝 𝑑𝑇 − 𝑇 𝑑𝑝
𝜕𝑇 𝜕𝑇 𝑝
TRACE KTU
𝑝
𝜕𝑠 𝜕𝑣 𝜕𝑣
=− 𝑑ℎ = 𝑐𝑝 𝑑𝑇 − 𝑇 − 𝑣 𝑑𝑝
𝜕𝑝 𝜕𝑇 𝜕𝑇 𝑝
𝑇 𝑝
𝜕𝑣 𝜕ℎ 𝜕𝑣
𝑑ℎ = 𝑐𝑝 𝑑𝑇 − 𝑇 − 𝑣 𝑑𝑝 =− 𝑇 −𝑣
𝜕𝑇 𝑝 𝜕𝑝 𝜕𝑇
𝑇 𝑝
TRACE KTU
ℎ = 𝑓 𝑝, 𝑇
𝜂𝐽𝑇 =
𝜕𝑇
𝜕𝑝
=−
𝜕ℎ
𝜕𝑝
𝜕𝑇
𝜕ℎ
𝜕ℎ 𝜕ℎ ℎ 𝑇 𝑝
𝑑ℎ = 𝑑𝑝 + 𝑑𝑇
𝜕𝑝 𝑇
𝜕𝑇 𝑝
1 𝜕𝑣
𝜂𝐽𝑇 = 𝑇 −𝑣
𝑐𝑝 𝜕𝑇 𝑝
28-May-23 Dr. Kishor Kumar V. V., MED, GCEK 29
1 𝜕𝑣 Joule-Thomson effect
𝜂𝐽𝑇 = 𝑇 −𝑣
𝑐𝑝 𝜕𝑇 𝑝
➢ An ideal gas will not experience
𝑅𝑇
𝐹𝑜𝑟 𝑎𝑛 𝑖𝑑𝑒𝑎𝑙 𝑔𝑎𝑠, 𝑣 = temperature change upon expansion
𝑝
through an expansion valve
𝜕𝑣
𝜕𝑇
TRACE KTU
=
𝑅 𝑣
=
𝑃 𝑇 ➢ Gases are imperfect at low enough
𝑝
𝑚ℎ
ሶ 2 = 𝑚ሶ − 𝑚ሶ 𝑓 ℎ1 + 𝑚ሶ 𝑓 ℎ𝑓
TRACE KTU
Fraction of mass flow that is liquefied
𝑚ሶ 𝑓 ℎ1 − ℎ2
=𝑦=
𝑚ሶ ℎ1 − ℎ𝑓
TRACE KTU
For an isothermal compression, 𝑄𝑅 = 𝑚𝑇
𝑚𝑇
ሶ 1 (𝑠2 − 𝑠1 ) − 𝑊 = 𝑚(ℎ
ሶ 1 (𝑠2 − 𝑠1 )
ሶ 2 − ℎ1 )
TRACE KTU 𝑚ሶ 𝑓
𝑚ሶ
=𝑦=
ℎ1 − ℎ2
ℎ1 − ℎ𝑓
Work required per unit mass liquefied is
−𝑊 −𝑊 ℎ1 − ℎ𝑓
= = 𝑇1 𝑠1 − 𝑠2 − ℎ1 − ℎ2
𝑚ሶ 𝑓 𝑚𝑦
ሶ ℎ1 − ℎ2
28-May-23 Dr. Kishor Kumar V. V., MED, GCEK 37
Simple Linde-Hampson system 𝑦=
ℎ1 − ℎ2
ℎ1 − ℎ𝑓
TRACE KTU𝜂𝐽𝑇 =
𝜕𝑇
𝜕𝑝 ℎ
=−
𝜕ℎ
𝜕𝑝 𝑇
𝜕𝑇
𝜕ℎ 𝑝
𝜕ℎ
𝜕𝑇 𝑝
= 𝑐𝑝
𝜕ℎ
= 0 = 𝜂𝐽𝑇 𝑐𝑝
𝜕𝑝 𝑇=𝑇
1
For optimum performance (maximum liquid yield),
state 2 should lie on the inversion curve
28-May-23 Dr. Kishor Kumar V. V., MED, GCEK 39
Simple Linde-Hampson
system
TRACE KTU
ℎ2 at room temperature and hence y is negative
TRACE KTU
➢ A separate refrigeration system using a fluid such as carbon dioxide, ammonia,
or a Freon compound is used to cool the main gas stream
TRACE KTU
Simple Linde-Hampson system Precooled Linde-Hampson system
28-May-23 Dr. Kishor Kumar V. V., MED, GCEK 45
Precooled Linde-Hampson system
TRACE KTU
▪ Compression process in precooling system is
assumed to be adiabatic
▪ R134a, NH3 and CO2 are the common
refrigerants in the precooling system
TRACE KTU
TRACE KTU
𝑚ℎ
ሶ 2 + 𝑚ሶ 𝑟 ℎ𝑑 = 𝑚ሶ − 𝑚ሶ 𝑓 ℎ1 + 𝑚ሶ 𝑟 ℎ𝑎 + 𝑚ሶ 𝑓 ℎ𝑓
TRACE KTU
𝑚ሶ 𝑓
𝑚ሶ
=𝑦=
ℎ1 − ℎ2
ℎ1 − ℎ𝑓
+𝑟
ℎ𝑎 − ℎ𝑑
ℎ1 − ℎ𝑓
𝑟=
𝑚ሶ 𝑟
𝑚ሶ
𝐸𝑖𝑛 = 𝐸𝑜𝑢𝑡
𝑚ℎ
ሶ 1 − 𝑊𝑐1 = 𝑚ℎ
ሶ 2 − 𝑄𝑅
TRACE KTU
For an isothermal compression, QR = mT
𝑚ℎ
ሶ 1 − 𝑊𝑐1 = 𝑚ℎ
ሶ 2 − mT
ሶ 1 s2 − s1
ሶ 1 s2 − s1
−𝑊𝑐1 = mT
ሶ 1 s1 − s2 − 𝑚ሶ ℎ1 − ℎ2
𝐸𝑖𝑛 = 𝐸𝑜𝑢𝑡
−𝑊𝑐2 = 𝑚ሶ 𝑟 (ℎ𝑏 − ℎ𝑎 )
𝑊𝑐 = 𝑊𝑐1 + 𝑊𝑐2
−𝑊𝑐1 = mT
ሶ 1 s1 − s2 − 𝑚ሶ ℎ1 − ℎ2
ሶ 1 s1 − s2 − 𝑚ሶ ℎ1 − ℎ2 +𝑚ሶ 𝑟 (ℎ𝑏 − ℎ𝑎 )
−𝑊𝑐 = mT
−𝑊𝑐 ሶ 𝑠2 − ℎ1 − ℎ2 + 𝑟 ℎ𝑏 − ℎ𝑎 𝑚ሶ 𝑟
= 𝑇1 𝑠1 − 𝑟=
𝑚ሶ 𝑚ሶ
gas is 𝑚𝑅𝑇
ሶ 1 𝑙𝑛
TRACE KTU
𝑝2
𝑝1
• Two compressors
• A 3-fluid heat exchanger
TRACE KTU • Two J-T expansion devices
• Two liquid containers
TRACE KTU 𝑚ℎ
ሶ 3 = 𝑚ሶ 𝑖 ℎ2 + 𝑚ሶ − 𝑚ሶ 𝑓 − 𝑚ሶ 𝑖 ℎ1 + 𝑚ሶ 𝑓 ℎ𝑓
𝑚ሶ 𝑓 ℎ1 − ℎ𝑓 = 𝑚ሶ ℎ1 − ℎ3 − 𝑚ሶ 𝑖 ℎ1 − ℎ2
𝑚ሶ 𝑓 ℎ1 − ℎ3 ℎ1 − ℎ2 𝑚ሶ 𝑖
= −𝑖 𝑖=
𝑚ሶ ℎ1 − ℎ𝑓 ℎ1 − ℎ𝑓 𝑚ሶ
ℎ1 − ℎ3 ℎ1 − ℎ2
TRACE KTU 𝑦=
ℎ1 − ℎ𝑓
−𝑖
ℎ1 − ℎ𝑓
𝑚ሶ 𝑖
Where 𝑖 =
𝑚ሶ
𝐸𝑖𝑛 = 𝐸𝑜𝑢𝑡
𝑚ሶ − 𝑚ሶ 𝑖 ℎ1 − 𝑊𝑐1 = 𝑚ሶ − 𝑚ሶ 𝑖 ℎ2 − 𝑄𝑅1
−𝑊𝑐1 = 𝑚ሶ − 𝑚ሶ 𝑖 ℎ2 − 𝑚ሶ − 𝑚ሶ 𝑖 ℎ1 − 𝑚ሶ − 𝑚ሶ 𝑖 𝑇1 𝑠2 − 𝑠1
−𝑊𝑐1 = 𝑚ሶ − 𝑚ሶ 𝑖 𝑇1 𝑠1 − 𝑠2 − ℎ1 − ℎ2
28-May-23 Dr. Kishor Kumar V. V., MED, GCEK 58
Linde dual-pressure system
Applying the first law
𝐸𝑖𝑛 = 𝐸𝑜𝑢𝑡
𝑚ℎ
ሶ 2 − 𝑊𝑐2 = 𝑚ℎ
ሶ 3 − 𝑄𝑅2
𝑚ℎ
ሶ 2 − 𝑊𝑐2 = 𝑚ℎ
ሶ 3 − 𝑚𝑇
ሶ 1 𝑠3 − 𝑠2
−𝑊𝑐2 = 𝑚ሶ 𝑇1 𝑠2 − 𝑠3 − ℎ2 − ℎ3
𝑊𝑐 = 𝑊𝑐1 + 𝑊𝑐2
−𝑊𝑐1 = 𝑚ሶ − 𝑚ሶ 𝑖 𝑇1 𝑠1 − 𝑠2 − ℎ1 − ℎ2
TRACE KTU
−𝑊𝑐2 = 𝑚ሶ 𝑇1 𝑠2 − 𝑠3 − ℎ2 − ℎ3
−𝑊𝑐
= 𝑇1 𝑠1 − 𝑠3 − ℎ1 − ℎ3 − 𝑚ሶ 𝑖
𝑚ሶ 𝑖=
𝑖 𝑇1 𝑠1 − 𝑠2 − ℎ1 − ℎ2 𝑚ሶ
TRACE KTU
➢ If the expansion engine is reversible and adiabatic, the expansion
process is isentropic and a much lower temperature is attained than
for an isenthalpic expansion
TRACE KTU
𝑚ℎ
ሶ 2 = 𝑊𝑒 + 𝑚ሶ − 𝑚ሶ 𝑓 ℎ1 + 𝑚ሶ 𝑓 ℎ𝑓
TRACE KTU 𝑊𝑒 = 𝑚ሶ 𝑒 ℎ3 − 𝑚ሶ 𝑒 ℎ𝑒
𝑚ℎ
ሶ 2 = 𝑚ሶ 𝑒 ℎ3 − 𝑚ሶ 𝑒 ℎ𝑒 + 𝑚ሶ − 𝑚ሶ 𝑓 ℎ1 + 𝑚ሶ 𝑓 ℎ𝑓
𝑚ሶ 𝑓 ℎ1 − ℎ𝑓 = 𝑚ሶ ℎ1 − ℎ2 + 𝑚ሶ 𝑒 ℎ3 − ℎ𝑒
𝑚ሶ 𝑓 ℎ1 − ℎ𝑓 = 𝑚ሶ ℎ1 − ℎ2 + 𝑚ሶ 𝑒 ℎ3 − ℎ𝑒
TRACE KTU 𝑦=
𝑚ሶ 𝑓
𝑚ሶ
=
ℎ1 −ℎ2
ℎ1 −ℎ𝑓
+𝑥
ℎ3 −ℎ𝑒
ℎ1 −ℎ𝑓
𝑚ሶ 𝑒
Where expander mass flow rate, x =
𝑚ሶ
−𝑊𝑐 + 𝑚ℎ
ሶ 1 = −𝑄𝑅 + 𝑚ℎ
ሶ 2
TRACE KTU 𝑄𝑅 = 𝑚𝑇
ሶ 1 𝑠2 − 𝑠1
−𝑊𝑐 + 𝑚ℎ
ሶ 1 = −𝑚𝑇
ሶ 1 𝑠2 − 𝑠1 + 𝑚ℎ
ሶ 2
−𝑊𝑐 = 𝑚ሶ 𝑇1 𝑠1 − 𝑠2 − ℎ1 − ℎ2
𝑚ሶ 𝑒 ℎ3 = 𝑊𝑒 + 𝑚ሶ 𝑒 ℎ𝑒
𝑊𝑒 = 𝑚ሶ 𝑒 ℎ3 − ℎ𝑒
−𝑊𝑐 = 𝑚ሶ 𝑇1 𝑠1 − 𝑠2 − ℎ1 − ℎ2
𝑊𝑒 = 𝑚ሶ 𝑒 ℎ3 − ℎ𝑒
−𝑊𝑛𝑒𝑡
= 𝑇1 𝑠1 − 𝑠2 − ℎ1 − ℎ2 − 𝑥 ℎ3 − ℎ𝑒
𝑚ሶ
➢ With precooled system, the temperature of the gas entering the basic
Linde-Hampson part of the liquefier can be lowered below ambient
temperature by choosing the correct fluid to precool the system
28-May-23 Dr. Kishor Kumar V. V., MED, GCEK 72
LIQUEFACTION SYSTEM FOR NEON AND
HYDROGEN
➢ Any fluid that has a triple point temperature below that of the
maximum inversion temperature of neon or hydrogen can be used as
a precoolant
𝑚ሶ 𝑁2 ℎ2 − ℎ1 ℎ1 − ℎ𝑓
𝑧= = +𝑦
𝑚ሶ ℎ𝑐 − ℎ 𝑎 ℎ𝑐 − ℎ𝑎
𝑚ሶ 𝑁2 𝑧
=
𝑚ሶ 𝑓 𝑦
ℎ1 − ℎ2 ∆ℎ1 ∆ℎ2
𝑦= + 𝑥1 + 𝑥2
ℎ1 − ℎ𝑓 ℎ1 − ℎ𝑓 ℎ1 − ℎ𝑓
𝑚ሶ 𝑒1 𝑚ሶ 𝑒2
𝑥1 = 𝑥2 =
𝑚ሶ 𝑚ሶ
𝑊𝑒1 = 𝑚ሶ 𝑒1 ∆ℎ1
𝑊𝑒2 = 𝑚ሶ 𝑒2 ∆ℎ2
𝑄ሶ 𝑅 − 𝑊ሶ 𝑖 = 𝑚ሶ ℎ𝑓 − ℎ1 = −𝑚ሶ ℎ1 − ℎ𝑓
Heat transfer process is reversible and isothermal
From II law,
𝑄ሶ 𝑅 = 𝑚𝑇
ሶ 1 𝑠2 − 𝑠1 = −𝑚𝑇
ሶ 1 𝑠1 − 𝑠𝑓
−𝑊ሶ 𝑖 ሶ 𝑠 − ℎ − ℎ =−𝑊ሶ 𝑖
= 𝑇1 𝑠1 − 𝑓 1 𝑓
𝑚ሶ ሶ 𝑚𝑓
𝑄ሶ 𝑅 − 𝑊ሶ 𝑖 = 𝑚ሶ ℎ𝑓 − ℎ1 = −𝑚ሶ ℎ1 − ℎ𝑓
𝑄ሶ 𝑅 = 𝑚𝑇
ሶ 1 𝑠2 − 𝑠1 = −𝑚𝑇
ሶ 1 𝑠1 − 𝑠𝑓
ሶ 1 𝑠1 − 𝑠𝑓 − 𝑊ሶ 𝑖 = −𝑚ሶ ℎ1 − ℎ𝑓
−𝑚𝑇
−𝑊ሶ 𝑖 ሶ −𝑊ሶ 𝑖
= 𝑇1 𝑠1 − 𝑠𝑓 − ℎ1 − ℎ𝑓 = ሶ
𝑚ሶ 𝑚𝑓