Исходные р –v и Т – s диаграммы.
рис.1
0 – к – политропное сжатие в компрессоре;
к – а – изобарный отвод тепла в теплообменный аппарат;
а – с – политропное сжатие в цилиндре двигателя;
с – z – изохорный подвод тепла в цилиндре двигателя;
z – z – изобарный подвод тепла в цилиндре двигателя;
z – b – политропное расширение в цилиндре;
b – 1 – изохорный отвод тепла из цилиндра двигателя;
1 – 2 – изобарный подвод тепла перед турбиной;
2 – 3 – политропное расширение газов в турбине;
3 – 0 – изобарный отвод тепла из цикла.
Расчёт параметров в основных точках цикла.
Точка "0"
Задано р0= 0,103 МПа
t0 =17° С
Т0 = t0+273,15 =17+273,15= 290,15 (K)
По уравнению состояния определяем:
v0= = = 0,8085
Точка "к"
Задано:
=3
= 0,805
= рк= р0 = 0,1033 = 0,309 (МПа)
=
s=0 в обратимом адиабатическом процессе в изоэнтропийном процессе
s=s
s
s определяем методом интерполяции используя таб. 1 стр. 56 [1]
s = 6,6081+ (290,15-273,15) = 6,6655
s =6,6655+0,287*ln3=6,9808
Находим обратной интерполяцией Т
Т - температура в изоэнтропийном процессе.
Т =373,15+ (6,9808 – 6,9218) = 396,32 (К)
Находим интерполяцией h
h = 373,8 + (396,32 – 373,15) = 397,2
h0 = 273,2+ (290,15 – 273,15) = 290,3
hк = h0+ =290,3+ =423,1
Т =373,15+
Находим показатель политропы:
0,3407
n – 1= 0,3407n
n =
Находим v в т. к
v =v =0,8085
Проверка v = = 0,3918
Относительная погрешность
vк= = 100=0,03
Точка а
Задано
Степень охлаждения воздуха = =1,15
Та= = = 366,84 (К)
Для изобарного процесса:
=
va = = = 0,3408
ра= рк =0,309 (МПа)
Проверка v = = 0,3407
Относительная погрешность
va= = 100=0,03
4. Точка с
Задано: n1=1,368
= =12,5
vc = = = 0,0273
pc = pa = pa = 0,30912,51,368 = 9,784 (МПа)
Тc = Тa = Ta = 366,8412,51,368-1 929,26 (К)
Проверка v = = 0,0273
Относительная погрешность
vс= = 100=0,00
Точка z
Задано: рz = pz = 13,0 МПа
vz = vz = 0,0273
=
Тz = = = 1234,71 (К)
Проверка Т = = 1236,59 (К)
Относительная погрешность
Tz= = 100=0,15
Точка z
Задано: рz = 13,0 МПа
Степень предварительного расширения: = =1,65
vz = vz = 0,02731,65 = 0,0450
= =
Tz = Tz = 1234,711,65 = 2037,27 (K)
Проверка Т = = 2038,33 (К)
Относительная погрешность
Tz= = 100=0,05
7. Точка b
Задано: n2=1,248
z = =0,81
vb = v1 = vaz = 0,34080,81 = 0,2760
pb = pz = pzz = 13,00,811,248 = 1,352 (МПа)
Тb = Тz = 2037,27 1299,32 (К)
Проверка Т = = 1300,18 (К)
Относительная погрешность
Tb= = 100=0,07
8. Точка 1
Задано: v1 = vb = 0, 2760
w = = =1,14
p1 = = 0,2711 (МПа)
Процесс изохорный:
Т1 = Tb = 1299,32 = 260,54 (K)
Проверка Т = = 260,71 (К)
Относительная погрешность
T1= = 100=0,07
9. Точка 2
Задано: р2 = р1 = 0,2711 МПа
Из условия расчета цикла:
q2 = q1
Процесс b –1 – изобарный
q2 = ub -1
Процесс 1 – 2 – изохорный
q1 = h1 - 2
u1 = 159,1 + (260,54 – 223,15) = 188,0
u2 = 1021,0 + (1299,32 – 1298,15) = 1022,1
h1 = 223,1 + (260,54 – 223,15) = 260,5
h2 = 159,1 + (260,54 – 223,15) = 188,0
q1 = h1 - 2 = h2 - h1 = q2
h2 = q2+h1 = 834,1+260,5 = 1094,6
Обратной интерполяцией находим Тz:
Tz = 1023,15 + (1094,6 - 1072,4) = 1042,49 (К)
Процесс изобарный:
v2 = v1 = 0,2760 = 1,1043
Проверка v = = 1,1036
Относительная погрешность
v2= = 100=0,06
Точка "3"
Задано:
р3 = р0 = 0,103 МПа
= 0,825
= = 0,825
l = lрасп. = -h = h2 – h3
l = h2 - h
h - теоретическое значение энтальпии при изоэнтропном процессе.
h3 - действительное значение энтальпии.
s=0 в обратимом адиабатическом процессе в изоэнтропийном процессе
=
s = s - s - Rln
s
S определяем методом интерполяции используя таб. 1 стр. 56 [1]
S = 7,9939 + (1042,49 - 1023,15) = 8,0151
S = 8,0151+0,287*ln =7,7374
Находим обратной интерполяцией Т
Т - температура в изоэнтропийном процессе.
Т =773,15+ (7,7374 – 7,6808) = 815,48 (К)
Находим интерполяцией h
h = 792,4 + (815,48 – 773,15) = 839,1
h3 = h2+(h2 - h ) =1094,6+(1094,6-839,1)0,825 = 883,8
Т3 = 773,15 +
Находим показатель политропы:
0,2037
n – 1= 0,2037n
n =
Находим v3
v3=v2 =1,1043
Проверка v = = 2,3850
Относительная погрешность:
v3= = 100=0,08
Расчет удельной внутренней энергии, удельной энтальпии и удельной энтропии в основных точках термодинамического цикла.
1. Точка 0
T0 = 290,15 K
u0 = 194,8 + (290,15 – 273,15) = 207,0
h0 = 290,3
s = 6,6655
2. Точка k
Tк = 421,87 K
uк = 266,7 + (421,87 – 373,15) = 302,1
hк = 423,1
s = 6,9218 + (421,87 – 373,15) = 7,0458
3. Точка а
Tа = 366,84 K
uа = 230,7 + (366,84 – 323,15) = 262,2
hа = 323,4 + (366,84 – 323,15) = 367,4
s = 6,7769 + (366,84 – 323,15) = 6,9035
4. Точка с
Tс = 929,26 K
uс = 693,9 + (929,26 – 923,15) = 699,0
hс = 958,8 + (929,26 – 923,15) = 965,7
s = 7,8772 + (929,26 – 923,15) = 7,9500
5. Точка z
Tz = 1236,59 K
uz = 953,8 + (1236,59 – 1228,15) = 965,8
hz = 1304,9 + (1236,59 – 1228,15) = 1320,7
s = 8,2014 + (1236,59 – 1228,15) = 8,2142
6. Точка z
Tz = 2037,27 K
uz = 1652,2 + (2037,27 – 1973,15) = 1714,3
hz = 2218,5 + (2037,27 – 1973,15) = 2299,0
s = 8,7827 + (2037,27 – 1973,15) = 8,8230
7. Точка b
Tb = 1299,32 K
ub = 1022,1
hb = 1393,6 + (1299,32 – 1298,15) = 1394,0
s = 8,2718 + (1299,32 – 1298,15) = 8,2729
8. Точка 1
T1 = 260,54 K
u1 = 188,0
h1 = 260,5
s = 6,4005 + (260,54 – 223,15) = 6,5556
9. Точка 2
T2 = 1042,49 K
u2 = 778,7 + (1042,49 – 1023,15) = 795,4
h2 = 1094,6
s = 8,0151
10. Точка 3
T3 = 855,94 K
u3 = 570,5 + (855,94 – 773,15) = 965,8
h3 = 883,8
s = 7,6806 + (855,94 – 773,15) = 7,7917
Полученные данные сводим в таблицу 1.
Таб.1 Параметры в основных точках термодинамического цикла |
s0 |
кДж/кгК |
6,6655 |
7,0458 |
6,9035 |
7,8845 |
8,2142 |
8,8230 |
8,2729 |
6,5556 |
8,0151 |
7,7917 |
h |
кДж/кг |
290,3 |
423,1 |
367,4 |
965,7 |
1320,7 |
2299,0 |
1394,0 |
260,5 |
1094,6 |
883,8 |
|
u |
кДж/кг |
207,0 |
302,1 |
262,2 |
699,0 |
965,8 |
1714,3 |
1022,1 |
188,0 |
795,4 |
638,1 |
|
T |
К |
290,15 |
421,87 |
366,84 |
929,26 |
1234,71 |
2037,27 |
1299,32 |
260,54 |
1042,49 |
855,94 |
|
v |
м3/кг |
0,8085 |
0,3919 |
0,3408 |
0,0273 |
0,0273 |
0,0450 |
0,2760 |
0,2760 |
1,1043 |
2,3869 |
|
Р |
МПа |
0,103 |
0,309 |
0,309 |
9,784 |
13,0 |
13,0 |
1,352 |
0,2711 |
0,2711 |
0,103 |
|
Параметры
Точки |
0 |
k |
a |
c |
z |
z |
b |
1 |
2 |
3 |
Расчет термодинамических процессов цикла.
1. Процесс 0 – к – политропный процесс сжатия.
u = u0 – uk = 302,1 – 207,0 = 95,1
h = h0 – hk = 423,1 – 290,3 = 132,8
s = s – s - Rln = 7,0458 – 6,6655 – 0,287ln3 = 0,0650
Теплота:
q = cnT = cvm (Tk – T0)
k = = = 1,3964
cvm = = = 0,7220
q = 0,7220 (421,87 – 290,15) = 22,2
Работа:
l = = = -73,1
Проверка по первому закону термодинамики:
q* = u + l = 95,1 – 73,1 = 22
Относительная ошибка:
q = 100 = 100 = 0,9
Проверка по второму закону термодинамики:
q** = Тсиs = 351,910,0650 = 22,9
Тси = = = 351,91(К)
Относительная ошибка:
q = 100 = 100 = 3,15
Изменение эксергии:
ех = h – T0s = 132,8 – 290,150,0650 = 113,9
Эксергетический КПД:
ех = = = = 0,858
Потери эксергии:
d = h - ех = 132,8 – 113,9 = 18,9
2. Процесс к – а – изобарный процесс охлаждения.
u = uк – uа = 262,2 – 302,1 = - 39,9
h = hк – hа = 367,4 – 423,1 = -55,7
s = s - s - Rln = 6,9035– 7,0458 – 0 = - 0,1423
Теплота:
q = h = -55,7
Работа:
l = рv = 0,309103 (0,3408 – 0,3919) = -15,8
Проверка по первому закону термодинамики:
q* = u + l = -39,9 – 15,8 = - 55,7
Относительная ошибка:
q = 100 = 100 = 0,00
Проверка по второму закону термодинамики:
q** = Тсиs = 393,72(-0,1423) = - 56,0
Тси = = = 393,72 (К)
Относительная ошибка:
q = 100 = 100 = 0,54
Изменение эксергии:
ех = h – T0s = -55,7 – 290,15(-0,1423) = -14,41
Эксергетический КПД:
ех = 0
Потери эксергии:
d = ех = 14,41
3. Процесс а – с – политропный процесс сжатия в цилиндре.
n1 = 1,368
u = uс – uа = 699,0 – 262,2 = 436,8
h = hс – hа = 965,7 – 367,4 = 598,3
s = s – s - Rln = 7,8845 – 6,9035 – 0,287ln = - 0,0106
Теплота:
q = cnT = cvm ( Tс – Tа)
k = = = 1,3697
cvm = = = 0,7766
q = 0,7766 (929,26 – 366,84) = - 2,0
Работа:
l = = = - 438,6
Проверка по первому закону термодинамики:
q* = u + l = 436,8 – 438,6 = - 1,8
Относительная ошибка:
q = 100 = 100 = 10
Проверка по второму закону термодинамики:
q** = Тсиs = 605,10(-0,0106) = - 6,4
Тси = = = 605,10 (К)
Относительная ошибка:
q = 100 = 100 = 220
Изменение эксергии:
ех = u – T0s + p0v = 436,8 – 290,15(- 0,0106) + 0,103103(0,0273 – - 0,3408) = 407,6
Эксергетический КПД:
ех = = = 0,929
Потери эксергии:
d = l - ех = 438,6 – 407,6 = 31,0
4. Процесс с – z – изохорный процесс.
u = uz – uс = 965,8 – 699,0 = 266,8
h = h z –hс = 1320,7– 965,7 = 355,0
s = s – s - Rln = 8,2142 – 7,8845 – 0,287ln = 0,2481
Теплота:
q = u = 266,8
Работа:
l = 0
Проверка по второму закону термодинамики:
q** = Тсиs = 1074,760,2481 = 266,6
Тси = = = 1074,76 (К)
Относительная ошибка:
q = 100 = 100 = 0,07
Изменение эксергии:
ех = u – T0s + p0v = 266,8 – 290,150,2481 + 0= 194,8
Эксергетический КПД:
ех = = = 0,730
Потери эксергии:
d = q - ех = 266,8 – 194,8 = 72,0
5. Процесс z – z – изобарный процесс охлаждения.
u = uz – uz = 1714,3 – 965,8 = 748,5
h = hz – h z = 2299,0 – 1320,7 = 978,3
s = s - s - Rln = 8,8230– 8,2142 – 0 = 0,6088
Теплота:
q = h = 834,1
Работа:
l = рv = 13,0103 (0,0450 – 0,0273) = 224,6
Проверка по первому закону термодинамики:
q* = u + l = 607,4 – 224,6 = 832,0
Относительная ошибка:
q = 100 = 100 = 0,25
Проверка по второму закону термодинамики:
q** = Тсиs = 1602,640,6088 = 975,7
Тси = = = 1602,64 (К)
Относительная ошибка:
q = 100 = 100 = 0,27
Изменение эксергии:
ех = u – T0s + p0v = 748,5 – 290,150,6088 + 0,103103(0,0450 – - 0,2730) = 590,1
Эксергетический КПД:
ех = = = 0,838
Потери эксергии:
d = q – (ех + l ) =978,3-(590,1+230,1) 158,1
6. Процесс z – b – политропный процесс расширения в цилиндре.
n2 = 1,248
u = ub– uz = 1022,1 – 1714,3 = - 692,2
h = hb – hz = 1394,0 – 2299,0 = - 905,0
s = s – s - Rln = 8,2729 – 8,8230 – 0,287ln = 0,0995
Теплота:
q = cnT = cvm ( Tb – Tz)
k = = = 1,3074
cvm = = = 0,9380
q = 0,9380 (1299,32 – 2037,27) = 165,8
Работа:
l = = = 854,0
Проверка по первому закону термодинамики:
q* = u + l = - 692,2 + 854,0 = 161,8
Относительная ошибка:
q = 100 = 100 = 2,41
Проверка по второму закону термодинамики:
q** = Тсиs = 1640,730,0995 = 163,4
Тси = = = 1640,73 (К)
Относительная ошибка:
q = 100 = 100 = 1,45
Изменение эксергии:
ех = u – T0s + p0v = -692,2 – 290,150,0995 + 0,103103(0,2760 – - 0,0450) = - 697,3
Эксергетический КПД:
ех = = = 0,9895
Потери эксергии:
d =( ех +q) - l = (-697,3 + 165,8) – 854,0 = 9,1
7. Процесс b – 1 – изохорный процесс.
u = u1 – ub = 188,0 – 1022,1 = - 834,21
h = h z –hс = 260,5– 1394,0 = - 1133,5
s = s – s - Rln = 6,5556 – 8,2729 – 0,287ln = -1,2561
Теплота:
q = u = - 834,1
Работа:
l = 0
Проверка по второму закону термодинамики:
q** = Тсиs = 646,47(-1,2561) = - 812,0
Тси = = = 646,47 (К)
Относительная ошибка:
q = 100 = 100 = 2,65
Изменение эксергии:
ех = h – T0s = -1133,5 – 290,15(-1,2561) = - 469,6
Эксергетический КПД:
ех = = = 0,979
Эксергия теплоты:
exq = q = (-834,1) = - 459,7
Потери эксергии:
d = ех-exq = -469,6 – -459,7 = 9,9
8. Процесс 1 – 2 – изобарный процесс.
u = u2 – u1 = 795,4 – 188,0 = 607,4
h = h2 – h 1 = 1094,6 – 260,5 = 834,1
s = s - s - Rln = 8,0151 – 6,5556 – 0 = 1,4595
Теплота:
q = h = 834,1
Работа:
l = рv = 0,2711103(1,1043 – 0,2760) = 230,1
Проверка по первому закону термодинамики:
q* = u + l = 748,5 – 230,1 = 978,6
Относительная ошибка:
q = 100 = 100 = 0,03
Проверка по второму закону термодинамики:
q** = Тсиs = 563,931,4595 = 823,1
Тси = = = 563,93 (К)
Относительная ошибка:
q = 100 = 100 = 1,32
Изменение эксергии:
ех = h – T0s = 834,1 – 290,151,4595 = 410,6
Эксергетический КПД:
ех = = = 1,569
Эксергия теплоты:
exq = q = 834,1 = 404,9
Потери эксергии:
d =exq – (ех + l ) = 404,9 - (410,6 +224,6) = -230,3
9. Процесс 2 – 3 – политропный процесс.
n = 1,2258
u = u3– u2 = 638,1 – 795,4 = - 157,3
h = h3 – h2= 883,8 – 1094,6 = - 210,8
s = s – s - Rln = 7,7917 – 8,0151 – 0,287ln = 0,0543
Теплота:
q = cnT = cvm ( T3 – T2)
k = = = 1,3401
cvm = = = 0,8432
q = 0,8432 (855,94 – 1042,49) = 52,0
Работа:
l = = = 209,3
Проверка по первому закону термодинамики:
q* = u + l = -157,3 + 209,3 = 52,0
Относительная ошибка:
q = 100 = 100 = 0,00
Проверка по второму закону термодинамики:
q** = Тсиs = 946,150,0543 = 51,4
Тси = = = 946,15 (К)
Относительная ошибка:
q = 100 = 100 = 1,15
Изменение эксергии:
ех = h – T0s = - 210,8 – 290,150,0543 = - 266,6
Эксергетический КПД:
ех = = = 0,863
Эксергия теплоты:
exq = q = 52,0 = 16,0
Потери эксергии:
d =( ех +q) - l = (-226,6 + 16,0) – 209,3 = 33,3
10. Процесс 3 – 0 – изобарный процесс.
u = u0 – u3 = 207,0 – 638,1 = - 431,1
h = h0 – h 3 =290,3 – 883,8 = 593,5
s = s - s - Rln = 6,6655 – 7,7917 – 0 = -1,1262
Теплота:
q = h = -593,5
Работа:
l = рv = 0,103103 (0,8085 – 2,3869) = - 162,6
Проверка по первому закону термодинамики:
q* = u + l = – 431,1 - 162,2 = -593,7
Относительная ошибка:
q = 100 = 100 = 0,03
Проверка по второму закону термодинамики:
q** = Тсиs = 523,01(-1,1262) = -589,0
Тси = = = 523,01 (К)
Относительная ошибка:
q = 100 = 100 = 0,76
Изменение эксергии:
ех = h – T0s = -593,5 – 290,15(-1,1262) = -266,7
Эксергетический КПД:
ех = 0
Потери эксергии:
d = ех =266,7
Значения, полученные в результате расчета, сводим в таблицу 2.
Термодинамический КПД цикла:
t = 0,587
Эксергетический КПД цикла:
ex = t = 0,587
Среднее давление цикла:
Рt = 2,641(МПа)
Таб.2 Параметры отдельных процессов и всего цикла. |
ex |
|
0,858 |
0 |
0,929 |
0,730 |
0,838 |
0,9895 |
0,979 |
1,569 |
0,863 |
0 |
|
d |
кДж/кг |
18,9 |
14,4 |
31,0 |
72 |
158,1 |
9,1 |
9,9 |
-230,3 |
33,3 |
266,7 |
383,1 |
|
ex |
кДж/кг |
113,9 |
-14,41 |
407,6 |
194,8 |
590,1 |
-697,3 |
-469,6 |
410,6 |
-226,6 |
-266,7 |
42,39 |
|
l |
кДж/кг |
-73,1 |
-15,8 |
-438,6 |
0 |
230,1 |
854,0 |
0 |
224,6 |
209,3 |
-162,6 |
827,9 |
|
q |
кДж/кг |
22,2 |
-55,7 |
-2,0 |
266,8 |
978,3 |
165,8 |
-834,1 |
834,1 |
52,0 |
-593,5 |
828,2 |
|
s |
кДж/кгК |
0,0650 |
-0,1423 |
-0,0106 |
0,2481 |
0,6088 |
0,0995 |
-1,2561 |
1,4595 |
0,0543 |
-1,1262 |
0,0000 |
|
h |
кДж/кг |
132,8 |
-55,7 |
598,3 |
355,0 |
978,3 |
-905,0 |
-1133,5 |
834,1 |
210,8 |
-593,5 |
0,0 |
|
u |
кДж/кг |
95,1 |
-39,9 |
436,8 |
266,8 |
748,5 |
-692,2 |
-834,1 |
607,4 |
-157,3 |
-431,1 |
0,0 |
|
Параметры
Процессы |
0 – k |
k – a |
a – c |
c – z |
z – z |
z – b |
b –1 |
1 – 2 |
2 – 3 |
3 – 0 |
Итого |
Построение p – v и T – s диаграмм.