tutorialETD Solns
tutorialETD Solns
1 import numpy as np
2 import matplotlib.pyplot as plt
3 import math
4 import scipy
5 import scipy.optimize as opt
6 from scipy import interpolate
Q1
Temperature of an object varies linearly as 𝑇 = 𝑎𝑋 + 𝑜𝑏 where 𝑋 is some property and 𝑎 and 𝑏 are constants.
The thermometer measures 10𝑜
when is 𝑋 10 and 35 when 𝑋 is halved. Determine the temperature when
𝑋 = 20 .
Click here for Soln.
Ans= -40
Q2
Determine the value of temperature corresponding to absolute zero using a constant volume gas thermometer.
Use freezing point and boiling point of water as and 50𝑜 100𝑜
respectively.
Click here for Soln.
Ans= -86.575
First Law
Q1
A rigid evacuated bottle having volume V is suddenly opened and atmospheric air ( 𝑃𝑜,𝑇𝑜 ) gradually flows in
and eventually the trapped air in the bottle reach in thermodynamic equilibrium with atmosphere. Determine the
net heat interaction during the entire filling process considering system as closed one.
Q2
Solve the above question considering system as open system
Q3
Determine the coefficient of thermal expansion for air at 20𝑜 𝐶 and 1 bar. Make suitable assumptions.
Click here for Soln.
localhost:8888/notebooks/Documents/academics/2019Odd/ETD/tutorialETD.ipynb# 1/17
11/20/2019 tutorialETD - Jupyter Notebook
Ans ≈ 0.003411𝐾 −1
Internal energy
Q4
Air at 300 K and 200 kPa is heated at constant pressure to 600 K. Determine the change in internal energy of
air.
Click here for Soln.
Δ𝑢 = 28.97 𝑘𝑔/𝑘𝑚𝑜𝑙
Using average value of specific heat at mean temperature (assuming linear variation with temperature)
𝑐𝑣,𝑎𝑣𝑔 = 𝑐𝑣@(𝑇1+𝑇2)/2 = 0.741 𝑘𝐽/𝑘𝑔𝐾 Δ𝑢 = 𝑐𝑣,𝑎𝑣𝑔 Δ(𝑇) = 222.37 𝑘𝐽/𝑘𝑔
In [ ]:
Q5
Determine the enthalpy of liquid water at 100𝑜 𝐶 and 15 Mpa.
Click here for Soln.
Using tables:
ℎ = 430.3 𝑘𝐽/𝑘𝑔
Using approximation @ saturation state (corresponding to temperature)
ℎ = 419 𝑘𝐽/𝑘𝑔
Click here for Soln.
Using tables:
Assuming the change of state from saturation state to high pressure at same temperature and finding enthalpy
change during pressure increase (using (𝑑ℎ = 𝑑𝑢 + 𝑝𝑑𝑣 + 𝑣𝑑𝑝) )
ℎ = 419 + 𝑣Δ𝑝
From table v @ saturation state = 0.001044𝑚3 /𝑘𝑔
localhost:8888/notebooks/Documents/academics/2019Odd/ETD/tutorialETD.ipynb# 2/17
11/20/2019 tutorialETD - Jupyter Notebook
In [ ]:
Q6
A stationary fluid system goes through a cycle comprising of following processes:
Q7
A working fluid of unit mass undergoes following cycle comprising of:
Given that 𝑃1 = 1 𝑏𝑎𝑟, 𝑃2 = 5 𝑏𝑎𝑟, 𝑣1 = 0.827 𝑚3 /𝑘𝑔 , 𝛾 = 1.4, calculate workdone in the cycle.
Click here for Soln.
localhost:8888/notebooks/Documents/academics/2019Odd/ETD/tutorialETD.ipynb# 3/17
11/20/2019 tutorialETD - Jupyter Notebook
In [ ]:
1 p1 = 125
2 p2 = 300
3 T1 = 325
4 T2 = 500
5 # pv^n = const.
6 temp = math.log(p1/p2)/math.log(T1/T2)
7 n = temp/(temp-1)
8 print('polytropic index:',n)
9 # specific work done (p1v1 -p2v2)/(n-1)
10 w = 0.287*(325 - 500)/(n -1)
11 print('specific work:',w)
Q9
A piston cylinder contains air at 600 kPa, 290 K and a volume of 0.01 𝑚3 . A constant pressure process gives 54
kJ of work out. Find the final volume and temperature of the air.
In [13]:
1 # W = PdV
2 p1 = 600 #kPa
3 T1 = 290 #K
4 V1 = 0.01 #m^3
5 W = 54 #kJ
6 #Constant pressure expansion
7 dV = W/p1
8 V2 = V1 + dV
9 print("final volume",V2,"m^3")
10 #assumption Ideal gas
11 T2 = V2/(V1)*T1
12 print("final Temperature",T2, "K")
Q10
A bicycle pump has a total stroke of 25 cm and is used to pump air into a tyre against a pressure of 3.5 bar.
Calculate the length of stroke necessary before the air enters the tyre when the piston is pushed in.
Assumptions:
during every stroke a fixed amount of air is compressed from atmospheric pressure (1 bar) to 3.5 bar (that
is when NRV opens and tyre gets filled with air)
air is ideal gas
What is known ?
what is required ?
𝑣2
soln 10
For the slow isothermal process 𝑝𝑣 = 𝑐
𝑃1 𝑣 1 = 𝑃2 𝑣 2
𝑣2 = 𝑣1 𝑃𝑃12
(25 − 𝑙𝑖 )𝐴 = 25𝐴 3.51
𝑙𝑖 = 17.85 𝑐𝑚
Q11
A helium gas is heated at constant volume from 100 kPa, 300 K to 500 K. A following process expands the gas
at constant pressure to three times the initial volume. What is the specific work in the combined process?
In [14]:
Q12
A soap bubble has a surface tension of 𝑠 = 3 × 10−4𝑁/𝑐𝑚
as it sits flat on a rigid ring of diameter 5 cm. You
now blow on the film to create a halfsphere surface of diameter 5 cm. How much work was done?
Click here for Soln.
𝑊 = 1.18 × 10−4𝐽
Q13
localhost:8888/notebooks/Documents/academics/2019Odd/ETD/tutorialETD.ipynb# 5/17
11/20/2019 tutorialETD - Jupyter Notebook
Assume that we fill a spherical balloon from a bottle of helium gas. The helium gas provides work to stretch the
balloon material and push back the atmosphere. Write the incremental balance for
𝑑𝑊ℎ𝑒𝑙𝑖𝑢𝑚 = 𝑑𝑊𝑠𝑡𝑟𝑒𝑡𝑐ℎ + 𝑑𝑊𝑎𝑡𝑚 to establish the connection between the helium pressure, the surface tension
𝑠 𝑃𝑜
, and as a function of the radius.
𝑃𝐻𝑒 = 𝑃𝑜 + 4 𝑠𝑟
Click here for Soln.
Q14
An initially deflated and now flat balloon is connected by a valve to a 12 𝑚3
storage tank containing helium gas
at 2 MPa and ambient temperature, 20𝑜 𝐶
. The valve is opened and the balloon is inflated at constant pressure,
𝑃𝑜 = 100 kPa, equal to ambient pressure, until it becomes spherical at 𝐷1
= 1 m. If the balloon is larger than
this, the balloon material is stretched, giving an inside pressure of
𝑃 = 𝑃𝑜 + 𝐶( 𝐷𝐷1 − 𝐷𝐷12 )
2
The balloon is inflated to a final diameter of 4 m, at which point the pressure inside is 400 kPa. The temperature
remains constant at 20𝑜 𝐶
. What is the maximum pressure inside the balloon at any time during the inflation
process? What is the pressure inside the helium storage tank at this time?
localhost:8888/notebooks/Documents/academics/2019Odd/ETD/tutorialETD.ipynb# 6/17
11/20/2019 tutorialETD - Jupyter Notebook
In [33]:
1 n = 30
2 x = np.linspace(1, 4, n)
3 C = 1600
4 P0 = 100
5 D1 = 1.6
6 P = P0 + C*(D1/x -D1**2/x**2)
7 #----------------------------
8 plt.plot(x, P, 'b--')
9 plt.grid(True)
10 #plt.show()
11 #-------------------
12 #x = np.linspace(-100,100,10000)
13 #y = np.sin(x)
14 tck = interpolate.splrep(x,P)
15 #tck
16 x0 = 2.0
17 y0 = interpolate.splev(x0,tck)
18 dydx = interpolate.splev(x0,tck,der=1)
19 tngnt = lambda x: dydx*x + (y0-dydx*x0)
20
21 plt.plot(x0,y0, "or")
22 plt.plot(x,tngnt(x), label="tangent")
23
24 #plt.show()
Out[33]:
[<matplotlib.lines.Line2D at 0xe90dfd0>]
In [31]:
1 def f(x):
2 return P0 + C*(D1/x -D1**2/x**2)
3 #return -2 * x**2 + 4 * x
4 max_x = opt.minimize_scalar(lambda x: -f(x), bounds=[1,4], method='bounded')
5 print(max_x)
fun: -499.9999999999804
message: 'Solution found.'
nfev: 10
status: 0
success: True
x: 3.2000007081104274
localhost:8888/notebooks/Documents/academics/2019Odd/ETD/tutorialETD.ipynb# 7/17
11/20/2019 tutorialETD - Jupyter Notebook
Any change in size of the control volume would require either a change in mass inside or a change in state
inside, neither of which is possible in a steady-state process.
In both processes a flow moves from a higher to a lower pressure. In the nozzle the pressure drop generates
kinetic energy, whereas that does not take place in the throttle process. The pressure drop in the throttle is due
to a flow restriction and represents a loss.
If you throttle a saturated liquid what happens to the fluid state? What if this is done to an ideal gas?
Click here for Soln.
The throttle process is approximated as a constant enthalpy process. Changing the state from saturated liquid
to a lower pressure with the same h gives a two-phase state so some of the liquid will vaporize and it becomes
colder.
If the same process happens in an ideal gas then same h gives the same temperature (h a function of T only) at
the lower pressure.
If you compress air the temperature goes up, why? When the hot air, high P flows in long pipes it eventually
cools to ambient T. How does that change the flow?
Click here for Soln.
Only the pressure can be felt upstream in a subsonic flow. In a supersonic flow no information can travel
upstream. The temperature information travels by conduction and even small velocities overpowers the
conduction with the convection of energy so the temperature at a given location is mainly given by the upstream
conditions and influenced very little by the downstream conditions.
Soln.
Use first law:
accordingly:
- temperature will be highest in Joule expansion (total loss, completely irreversible process)
- tempearture will be lowest in reversible expansion process
Important point to note here is that area under the curve is not work done for irreversible process
A constant mass goes through a process in which 100 J of heat transfer comes in and 100 J of work leaves.
Does the mass change state?
localhost:8888/notebooks/Documents/academics/2019Odd/ETD/tutorialETD.ipynb# 8/17
11/20/2019 tutorialETD - Jupyter Notebook
An ideal gas in a piston-cylinder is heated with 2 kJ during an isothermal process. How much work is involved?
A 25 kg steel tank initially at – 10𝑜 𝐶 is filled up with 100 kg of milk (assume properties as water) at 30𝑜 𝐶 . The
milk and the steel come to a uniform temperature of +5𝑜 𝐶 in a storage room. How much heat transfer is
needed for this process? Given: 𝐶𝑠𝑡𝑒𝑒𝑙 = 0.466 𝑘𝐽/𝑘𝑔.𝐾, 𝐶𝑚𝑖𝑙𝑘 = 4.18𝑘𝐽/𝑘𝑔.𝐾
Click here for Soln.
Soln.𝑉2 = 389𝑚/𝑠
Water at 180𝑜 𝐶 , 2000 kPa is throttled into a flash evaporator chamber having a pressure of 500 kPa. Neglect
any change in the kinetic energy. What is the fraction of liquid and vapor in the chamber?
Click here for Soln.
Second Law
1 T_h = np.linspace(40,80,100)
2 T_l = np.linspace(40,40,100)
3 dT = np.linspace(0,80,100)
4 # Initialize a Figure
5 fig = plt.figure()
6 # Add Axes to the Figure
7 ax = fig.add_axes([0,0,1,1])
8 plt.plot(T_h, 1-T_l/T_h)
9 T_h = np.linspace(40,40,100)
10 T_l = np.linspace(40,0,100)
11 plt.plot(T_l, 1-T_l/T_h)
12 T_h = np.linspace(80,80,100)
13 plt.plot(dT, dT/T_h)
14 ax.set(title = "$\eta$ V/S $T_{reservoir}$ for Heat Engine", xlabel = "$T_{L}$ or $T_{H
15 ax.legend(["$T_{L}$ = Const. (40)", '$T_{H}$ = Const.(40)','$\Delta T$', '$f_{xx}$', '$
16 #plt.show()
Out[2]:
<matplotlib.legend.Legend at 0x10a35350>
𝑇𝐿
variation of reverse carnot efficiency(COP) with 𝑇𝐻 and
localhost:8888/notebooks/Documents/academics/2019Odd/ETD/tutorialETD.ipynb# 10/17
11/20/2019 tutorialETD - Jupyter Notebook
In [3]:
1 T_h = np.linspace(45,85,100)
2 T_l = np.linspace(40,40,100)
3 dT = np.linspace(5,80,100)
4 fig = plt.figure()
5 ax = fig.add_axes([0,0,1,1])
6 plt.plot(T_h, T_h/(T_h-T_l))
7 T_h = np.linspace(45,45,100)
8 T_l = np.linspace(40,0,100)
9 plt.plot(T_l, T_h/(T_h-T_l))
10 plt.plot(dT, T_h/dT)
11 ax.set(title = "COP V/S $T_{reservoir}$ for Refrigerator", xlabel = "$T_{L}$ or $T_{H}$
12 ax.legend(["$T_{L}$ = Const. (40)", '$T_{H}$ = Const.(45)','$\Delta T$', '$f_{xx}$', '$
13
Out[3]:
<matplotlib.legend.Legend at 0x10b140b0>
localhost:8888/notebooks/Documents/academics/2019Odd/ETD/tutorialETD.ipynb# 11/17
11/20/2019 tutorialETD - Jupyter Notebook
In [10]:
Out[10]:
<matplotlib.legend.Legend at 0x11e649d0>
A refrigerator removes 1.5 kJ from the cold space using 1 kJ work input. How much energy goes into the
kitchen and what is its coefficient of performance?
Click here for Soln.
Soln. 𝑄𝐻 = 2.5 kJ
𝐶𝑂𝑃 = 1.5
A room is heated with a 1000 W electric heater. How much power can be saved if a heat pump with a COP of
2.0 is used instead?
Click here for Soln.
localhost:8888/notebooks/Documents/academics/2019Odd/ETD/tutorialETD.ipynb# 12/17
11/20/2019 tutorialETD - Jupyter Notebook
1 A farmer runs a heat pump with a 2 kW motor. It should keep a chicken hatchery at
$30^oC$, which loses energy at a rate of 10 kW to the colder ambient. What is the
minimum coefficient of performance that will be acceptable for the heat pump?
2 <div align="right">
3 <a href="#afr" class="btn btn-default" data-toggle="collapse">Click here for Soln.
</a>
4 </div>
5 <div id="afr" class="collapse">
6 <b> Soln.</b>
7 W = 2 kW
8 <br>
9 $Q_H = Q_{loss} = 10$ kW
10 <br>
11 COP = 10/2 = 5
12 </div>
For each of the cases below determine if the heat engine satisfies the first law (energy equation) and if it
violates the second law.
𝜂𝑠𝑡𝑒𝑎𝑚 𝑝𝑜𝑤𝑒𝑟 𝑝𝑙𝑎𝑛𝑡 = 0.42
Soln.
𝐶𝑂𝑃 = 1.38
Prove that a cyclic device that violates the Kelvin–Planck statement of the second law also violates the Clausius
statement of the second law or vice-versa.
A Carnot cycle heat engine has an efficiency of 40%. If the high temperature is raised 10% what is the new
efficiency keeping the same low temperature?
Click here for Soln.
localhost:8888/notebooks/Documents/academics/2019Odd/ETD/tutorialETD.ipynb# 13/17
11/20/2019 tutorialETD - Jupyter Notebook
In a cryogenic experiment you need to keep a container at −125°C although it gains 100 W due to heat transfer.
What is the smallest motor you would need for a refrigerator absorbing heat from the container and rejecting
heat to the room at 20°C?
Click here for Soln.
Two different fuels can be used in a heat engine operating between the fuel burning temperature and a low
temperature of 350 K. Fuel A burns at 2200 K delivering 30 000 kJ/kg and costs 105 Rs./kg. Fuel B burns at
1200 K, delivering 40 000 kJ/kg and costs 91 Rs./kg Which fuel will you buy and why?
Click here for Soln.
Soln. Select fuel B for more work per rupee though it has a lower thermal efficiency.
An ideal gas Carnot cycle with air in a piston cylinder has a high temperature of 1200 K and a heat rejection at
400 K. During the heat addition the volume triples. Find the two specific heat transfers (q) in the cycle and the
overall cycle efficiency.
Click here for Soln.
a) If the winter outside temperature is 0°C, what is the minimum power required to drive the heat pump?
b) For the same power as in part (a), what is the maximum outside summer temperature for which the
house can be maintained at 25°C?
Click here for Soln.
Entropy
Water at 100 kPa, 150𝑜 𝐶
receives 75 kJ/kg in a reversible process by heat transfer. Which process changes
entropy (s) the most: constant T, constant v or constant P?
Click here for Soln.
Soln. We know that "ds = dq/T" for reversible process, now from constant property lines on the T-s diagram.
The constant v line has a higher slope than the constant P line also at positive slope. Thus both the constant P
and v processes have an increase in T. As T goes up the change in s is smaller. Therefore the change in
entropy will be maximum in case of constant T (isothermal) process.
A substance has heat transfer out. Can you say anything about changes in s if the process is reversible? If it is
irreversible?
Click here for Soln.
Soln. If reversible: ds is less than zero, according to entropy balance eqn. If irreversible: can not say, it depends
upon the relative magnitude of dq/T and s_gen.
A substance is compressed adiabatically so P and T go up. Does that change entropy (s)?
localhost:8888/notebooks/Documents/academics/2019Odd/ETD/tutorialETD.ipynb# 14/17
11/20/2019 tutorialETD - Jupyter Notebook
Soln. If reversible: ds is zero, according to entropy balance eqn. If irreversible: If the process is irreversible then
s goes up as s_gen is always positive.
Soln. Use tables to calculate the missing properties: a) T = 64.9 C; x = 0.98; h = 2573.8 kJ/Kg b) x
undefined (superheated vapour); T = 682 𝑜𝐶
; s = 7.12 kJ/kgK c) x undefined (superheated vapour); P =
163.9 kPa
Soln. Use tables: state 1 is known given the temperature and wet state, 𝑣 = 𝑣𝑓 + 𝑥𝑣𝑓𝑔 𝑜= 0.223𝑚3 /𝑘𝑔 . Now
the state 2 has same specific volume (constant volume process) and Temperature = 140 𝐶 . Now using tables
find the other properties. x = 0.438; s = 4.015 kJ/kgK
Soln. Use tables: state 1 is known given the temperature and wet state, 𝑣 = 𝑣𝑓 + 𝑥𝑣𝑓𝑔 𝑜= 0.223𝑚3 /𝑘𝑔 Now
the state 2 has same specific volume (constant volume process) and Temperature = 140 𝐶 . Now using tables
find the other properties. x = 0.438; s = 4.015 kJ/kgK
Soln. Assume some specific heat capacity of steel say c = 0.46 kJ/kgK Both the blocks are finite bodies. so as
per the energy balance: 𝑇𝑓 = 137.5𝑜 𝐶 𝑠 = 0.1794𝑘𝐽/𝐾
and
Hint: Use tables and state 1 is known (given 2 properties) we know that throttling is constant enthalpy process.
So fix state 2 and determine the change in entropy using tables and interpolation. Δ𝑠 ≈ 0.419𝑘𝐽/𝑘𝑔𝐾
localhost:8888/notebooks/Documents/academics/2019Odd/ETD/tutorialETD.ipynb# 15/17
11/20/2019 tutorialETD - Jupyter Notebook
Hint: already discussed in the last class Entropy can be used to find the sense of flow in a well-insulated
duct given two distinct states.
localhost:8888/notebooks/Documents/academics/2019Odd/ETD/tutorialETD.ipynb# 16/17
11/20/2019 tutorialETD - Jupyter Notebook
Availibility
With keeping a source and sink temperature constant at 527𝑜 𝐶 27𝑜 𝐶
and respectively, 5000 kJ heat is
transfered from source to sink. Find out availability and unavailability for the reversible process.
Click here for Soln.
In [ ]:
localhost:8888/notebooks/Documents/academics/2019Odd/ETD/tutorialETD.ipynb# 17/17