Курсовая работа: Термодинамический расчет заданной смеси
Федеральное агентство по образованию
Государственное образовательное учреждение высшего профессионального образования
Самарский государственный технический университет
Кафедра: «Технология органического и нефтехимического синтеза»
КУРСОВАЯ РАБОТА
ПО КУРСУ:
«Физико-химические свойства растворов»
Выполнил: студент III – ХТ – 2 Степанов В. Н.
Принял: к. х. н., доцент Саркисова В. С.
Cамара, 2005 г.
Исходные данные:
Tb, К |
Tc(эксп)К |
Pc(эксп),bar |
Vc(эксп), см3/моль |
W(эксп) |
Zc |
yi |
|
2.3-диметилбутан | 331.13 | 499.98 | 31.27 | 358 | 0.247 | 0.270 | 0.37 |
цис-1,2-диметилциклогексан | 402.90 | 606.00 | 29.3 | 461,73 | 0.22 | ||
метил-третбутиловый эфир | 328.30 | 497.10 | 34.3 | 334 | 0.1 | ||
индан | 451.10 | 684.90 | 39.5 | 389 | 0.31 |
Задание:
1. Для четырехкомпонентной смеси заданного состава рассчитать энтальпию, энтропию, теплоемкость в стандартном состоянии при заданной температуре.
2. Псевдокритические свойства: температуру, давление, объем, ацентрический фактор и коэффициент сжимаемости.
3. Плотность:
3.1. Ненасыщенной газовой и жидкой смеси при температуре, соответствующей приведенной температуре 0,95 в диапазоне приведенных давлений от 0,01 до 1. Построить график зависимости.
3.2. Плотность жидкой смеси на линии насыщения в диапазоне температур от 298К до критической с шагом 25К. Построить график зависимости.
4. Энтальпию, энтропию, теплоемкость смеси при заданной температуре в диапазоне приведенных давлений от 0,01 до 10. Построить график зависимости соответствующего свойства от давления.
5. Энтальпию испарения в стандартном состоянии при давлении отличном от 1 атм для диапазона температур для диапазона температур от 298К до критической с шагом 25К.
6. «Кажущуюся» стандартную энтальпию образования смеси в жидком состоянии .
7. Вязкость смеси при температуре 730 К и приведенном давлении 10.
8. Теплопроводность смеси при температуре 730К и приведенном давлении 10.
Решение:
1. Для расчета т/д характеристик смеси нужно знать параметры компонентов этой смеси.
Для 2,3-диметилбутана при вычислении используем следующую формулу:
где для диапазона температур :
и т.д.
Значения , и рассчитаны по табличным данным методом Бенсона.
При расчете используем уже найденные значения для диапазонов температур:
При расчете получаем: .
Для цис-1,2-диметилциклогексана при вычислении используем следующую формулу:
где для диапазона температур :
и т.д.
Значения , и рассчитаны по табличным данным методом Бенсона.
При расчете цис-1,2-диметилциклогексана используем уже найденные значения для диапазонов температур:
При расчете цис-1,2-диметилциклогексана получаем: .
Для метил-третбутилового эфира при вычислении используем следующую формулу:
где для диапазона температур :
и т.д.
Значения , и рассчитаны по табличным данным методом Бенсона.
При расчете используем уже найденные значения для диапазонов температур:
При расчете получаем: .
Для индана при вычислении используем следующую формулу:
где для диапазона температур :
и т.д.
Значения , и рассчитаны по табличным данным методом Бенсона.
При расчете используем уже найденные значения для диапазонов температур:
.
При расчете получаем: .
Для смесей в состоянии идеального газа энтальпия, энтропия и теплоемкость рассчитываются аддитивно с учетом переменной состава:
;
;
.
2. Псевдокритический объем смеси определяется, как функция от состава смеси и критических объемов чистых компонентов смеси: , при этом получена матрица:
1 | 2 | 3 | 4 | |
1 | 392.0816 | 265.4731 | 102.3749 | 342.5268 |
2 | 265.4731 | 178.7814 | 69.4187 | 231.5124 |
3 | 102.3749 | 69.4187 | 26.7200 | 89.4788 |
4 | 342.5268 | 231.5124 | 89.4788 | 299.0632 |
Псевдокритическая температура рассчитывается, как функция от критической температуры, критического объема и состава смеси:
При этом получена следующая матрица:
1 | 2 | 3 | 4 | |
1 | 196032.9584 | 146127.9295 | 51037.78151 | 200440 |
2 | 146127.9295 | 108341.532 | 38100.82371 | 149150.3 |
3 | 51037.78151 | 38100.82371 | 13282.512 | 52210.23 |
4 | 200440.0092 | 149150.2849 | 52210.23284 | 204828.4 |
Ацентрический фактор смеси рассчитывается аддитивно: .
Ацентрические факторы чистых компонентов смеси рассчитываются по уравнению:
,
, где: , Рс – критическое давление, атм.
Псевдокритическое давление смеси рассчитывается через уравнение: , где: – псевдокритический коэффициент сжимаемости, определяемый с помощью таблиц Ли-Кесслера и разложения Питцера.
3.
3.1. Плотность ненасыщенной газовой и жидкой смеси определяем по формуле:, где:
- молекулярная масса смеси, определяемая по правилу аддитивности.
- объем смеси, определяемый по уравнению Менделеева-Клапейрона:
,
- фактор сжимаемости смеси, определяемый с помощью разложения Питцера и таблиц Ли-Кесслера.
Для диапазона приведенных давлений от 0,01 до 1 имеем:
Pr |
P, атм |
z(0) |
z(1) |
z |
Vm, см3/моль |
ρm, г/см3 |
0.01 | 0.3316 | 0.9961 | -0.0012 | 0.9958 | 135736.7088 | 0.00075 |
0.05 | 1.6578 | 0.9803 | -0.0062 | 0.9787 | 26680.6560 | 0.00382 |
0.10 | 3.3156 | 0.96 | -0.0126 | 0.9567 | 13040.6132 | 0.00782 |
0.20 | 6.6312 | 0.9174 | -0.0262 | 0.9105 | 6205.5229 | 0.01644 |
0.40 | 13.2624 | 0.8206 | -0.0589 | 0.8051 | 2743.5055 | 0.03719 |
0.60 | 19.8936 | 0.6967 | -0.1110 | 0.6674 | 1516.3113 | 0.06728 |
0.80 | 26.5248 | 0.141 | -0.0540 | 0.1268 | 215.9884 | 0.47236 |
1.00 | 33.1559 | 0.1705 | -0.0607 | 0.1545 | 210.5945 | 0.48446 |
1.20 | 39.7871 | 0.1998 | -0.0678 | 0.1819 | 206.6516 | 0.49370 |
1.50 | 49.7339 | 0.2432 | -0.0788 | 0.2224 | 202.1250 | 0.50476 |
2.00 | 66.3119 | 0.3138 | -0.0967 | 0.2883 | 196.4949 | 0.51922 |
3.00 | 99.4678 | 0.4501 | -0.1310 | 0.4156 | 188.8181 | 0.54033 |
5.00 | 165.7797 | 0.7092 | -0.1943 | 0.6580 | 179.3773 | 0.56877 |
7.00 | 232.0916 | 0.9561 | -0.2526 | 0.8895 | 173.2122 | 0.58902 |
10.00 | 331.5594 | 1.3108 | -0.3339 | 1.2228 | 166.6762 | 0.61211 |
Зависимость имеет вид:
3.2. Плотность жидкой смеси на линии насыщения определяется по формуле: , где:
определяется по уравнению Ганна и Ямады: , где:
,
и Г – функции приведенной температуры.
Для диапазона температур :
;
Для диапазона :
;
Для диапазона :
; .
Для диапазона температур от 298К до 573К имеем:
Т,К |
Тr |
Г |
Vr0 |
Vs c |
Vm |
ρm, г/см3 |
298 | 0.5139 | 0.2368 | 0.3656 | 382.393653 | 131.0658 | 0.7784 |
323 | 0.5570 | 0.2307 | 0.3755 | 134.8551 | 0.7565 | |
348 | 0.6001 | 0.2244 | 0.3862 | 138.9467 | 0.7343 | |
373 | 0.6432 | 0.2179 | 0.3980 | 143.4386 | 0.7113 | |
398 | 0.6863 | 0.2112 | 0.4111 | 148.4628 | 0.6872 | |
423 | 0.7295 | 0.2043 | 0.4262 | 154.1862 | 0.6617 | |
448 | 0.7726 | 0.1973 | 0.4436 | 160.8103 | 0.6344 | |
473 | 0.8157 | 0.1901 | 0.4652 | 168.9805 | 0.6038 | |
498 | 0.8588 | 0.1827 | 0.4946 | 180.0206 | 0.5667 | |
523 | 0.9019 | 0.1751 | 0.5307 | 193.5710 | 0.5271 | |
548 | 0.9450 | 0.1674 | 0.5852 | 213.9129 | 0.4769 | |
573 | 0.9881 | 0.1594 | 0.7210 | 264.1337 | 0.3863 |
Зависимость имеет вид:
4. Энтальпию, энтропию, теплоемкость смеси при заданной температуре 730К в диапазоне приведенных давлений от 0,01 до 10 определяем с помощью таблиц Ли-Кесслера и разложения Питцера:
,
,
.
Для энтропии в диапазоне приведенных давлений от 0,01 до 10 получаем:
Pr |
P, атм |
Н(0) |
Н(1) |
Н, Дж/моль |
0.01 | 0.3316 | 0.005 | 0.0014 | -150668.93 |
0.05 | 1.6578 | 0.0252 | 0.0062 | -150772.42 |
0.10 | 3.3156 | 0.0508 | 0.0118 | -150902.96 |
0.20 | 6.6312 | 0.1022 | 0.0236 | -151165.77 |
0.40 | 13.2624 | 0.2072 | 0.0432 | -151696.91 |
0.60 | 19.8936 | 0.3146 | 0.0584 | -152234.03 |
0.80 | 26.5248 | 0.4252 | 0.0688 | -152780.47 |
1.00 | 33.1559 | 0.5382 | 0.0738 | -153331.62 |
1.20 | 39.7871 | 0.653 | 0.0736 | -153884.84 |
1.50 | 49.7339 | 0.8304 | 0.0624 | -154725.88 |
2.00 | 66.3119 | 1.1288 | 0.0232 | -156114.70 |
3.00 | 99.4678 | 1.6722 | -0.0292 | -158667.92 |
5.00 | 165.7797 | 2.2994 | 0.2868 | -162093.48 |
7.00 | 232.0916 | 2.5242 | 0.7384 | -163751.36 |
10.00 | 331.5594 | 2.5684 | 1.2668 | -164636.17 |
График зависимости имеет вид:
Для энтропии расчет ведем по формуле:
Для энтропий в диапазоне приведенных давлений от 0,01 до 10 получаем:
Pr |
P, атм |
S(0) |
S(1) |
S,Дж/(Кмоль) |
0.01 | 0.3316 | 0.0024 | 0.0014 | 638.31 |
0.05 | 1.6578 | 0.0122 | 0.0082 | 651.59 |
0.10 | 3.3156 | 0.0246 | 0.0160 | 657.23 |
0.20 | 6.6312 | 0.0500 | 0.0310 | 662.75 |
0.40 | 13.2624 | 0.1020 | 0.0606 | 668.02 |
0.60 | 19.8936 | 0.1554 | 0.0872 | 670.89 |
0.80 | 26.5248 | 0.2112 | 0.1114 | 672.76 |
1.00 | 33.1559 | 0.2688 | 0.1322 | 674.09 |
1.20 | 39.7871 | 0.3282 | 0.1492 | 675.08 |
1.50 | 49.7339 | 0.4198 | 0.1682 | 676.13 |
2.00 | 66.3119 | 0.5772 | 0.1876 | 677.17 |
3.00 | 99.4678 | 0.8622 | 0.2396 | 678.06 |
5.00 | 165.7797 | 1.1740 | 0.5794 | 678.97 |
7.00 | 232.0916 | 1.2792 | 0.9490 | 680.08 |
10.00 | 331.5594 | 1.3168 | 1.3508 | 681.85 |
График зависимости имеет вид:
Расчет теплоемкости ведем по формуле:
Для теплоемкости в диапазоне приведенных давлений от 0,01 до 10 получаем:
Pr |
P, атм |
C(0) |
C(1) |
Cp, Дж/(К∙моль) |
0.01 | 0.3316 | 0.0064 | 0.0068 | 297.50 |
0.05 | 1.6578 | 0.0318 | 0.0348 | 297.77 |
0.10 | 3.3156 | 0.0648 | 0.0698 | 298.12 |
0.20 | 6.6312 | 0.1328 | 0.1384 | 298.84 |
0.40 | 13.2624 | 0.2790 | 0.2726 | 300.34 |
0.60 | 19.8936 | 0.4402 | 0.3988 | 301.96 |
0.80 | 26.5248 | 0.6186 | 0.5118 | 303.69 |
1.00 | 33.1559 | 0.8152 | 0.6074 | 305.54 |
1.20 | 39.7871 | 1.0304 | 0.6812 | 307.49 |
1.50 | 49.7339 | 1.3866 | 0.7450 | 310.59 |
2.00 | 66.3119 | 2.0330 | 0.7616 | 316.00 |
3.00 | 99.4678 | 3.0174 | 1.2050 | 325.16 |
5.00 | 165.7797 | 3.0630 | 3.5034 | 330.57 |
7.00 | 232.0916 | 2.5654 | 4.4436 | 328.50 |
10.00 | 331.5594 | 2.0500 | 4.5874 | 324.53 |
График зависимости имеет вид:
5. Энтальпию испарения смеси в стандартном состоянии и при давлении, отличном от 1атм рассчитываем по уравнению Менделеева-Клапейрона:
, где:
при стандартном состоянии =1, при давлении, отличном от 1атм для определения используется выражение:
;
;
;
,
давление насыщенных паров определяется по корреляции Ли-Кесслера:
, где:
,. .
,
,.
Для диапазона температур от 298 до 573К:
Т, К |
Тr |
f(0) |
f(1) |
Pvpr |
∆Zv |
∆vH, Дж/моль |
∆vH0, Дж/моль |
298 | 0.5139 | -5.0752 | -6.2981 | 0.00119 | 0.9956 | 37254.7001 | 37418.78 |
323 | 0.5570 | -4.2589 | -5.0156 | 0.00377 | 0.9890 | 36071.81517 | 36471.57 |
348 | 0.6001 | -3.5658 | -3.9893 | 0.00988 | 0.9769 | 34737.17427 | 35559.11 |
373 | 0.6432 | -2.9702 | -3.1614 | 0.02229 | 0.9572 | 33212.48933 | 34697.05 |
398 | 0.6863 | -2.4528 | -2.4888 | 0.04465 | 0.9284 | 31478.11583 | 33906.18 |
423 | 0.7295 | -1.9985 | -1.9385 | 0.08130 | 0.8891 | 29531.02569 | 33213.72 |
448 | 0.7726 | -1.5956 | -1.4850 | 0.13709 | 0.8383 | 27373.36201 | 32654.6 |
473 | 0.8157 | -1.2346 | -1.1076 | 0.21726 | 0.7744 | 24991.39522 | 32273.04 |
498 | 0.8588 | -0.9077 | -0.7895 | 0.32763 | 0.6948 | 22319.37859 | 32124.32 |
523 | 0.9019 | -0.6083 | -0.5165 | 0.47499 | 0.5938 | 19164.33515 | 32276.68 |
548 | 0.9450 | -0.3306 | -0.2762 | 0.66803 | 0.4566 | 14981.04654 | 32813.54 |
573 | 0.9881 | -0.0696 | -0.0577 | 0.91869 | 0.2186 | 7396.108869 | 33835.83 |
Зависимость имеет вид:
6. «Кажущаяся» стандартная энтальпия образования смеси в жидком состоянии определяется по уравнению: , где:
- энтальпия парообразования смеси при стандартном давлении и температуре 298К,
- энтальпия образования газовой смеси при стандартном давлении и температуре 298К, определена по правилу аддитивности. Энтальпии образования чистых компонентов смеси определены методом Бенсона.
7. Вязкость смеси при атмосферном давлении рассчитываем методом Голубева:
т.к. , то:
,
где
При расчете вязкости для высоких давлений используем корреляцию для смесей неполярных газов, т.к. большинство газов в смеси – неполярные:
где: , приведенную плотность смеси определяем по формуле: , где V рассчитывается по уравнению Менделеева-Клапейрона: , а псевдокритический коэф. сжимаемости определяется из разложения Питцера с помощью таблиц Ли-Кесслера: . Условие применимости для : удовлетворяется.
8. Теплопроводность смеси при стандартном давлении рассчитаем с помощью корреляции Мисика-Тодоса:
где: - теплоемкость смеси в стандартном состоянии при температуре 730К,
.
Теплопроводность при высоком давлении и определяется по формуле: