Скачать .docx |
Реферат: Термодинамический расчет смеси
Федеральное агентство по образованию
Государственное образовательное учреждение высшего профессионального образования
Самарский государственный технический университет
Кафедра: «Технология органического и нефтехимического синтеза»
КУРСОВАЯ РАБОТА ПО КУРСУ:
«Физико-химические свойства растворов»
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 |
Ацентрический фактор смеси рассчитывается аддитивно: .
Ацентрические факторы чистых компонентов смеси рассчитываются по уравнению:
,
,
где: , Рс – критическое давление, атм.
Псевдокритическое давление смеси рассчитывается через уравнение: , где: – псевдокритический коэффициент сжимаемости, определяемый с помощью таблиц Ли-Кесслера и разложения Питцера.
2.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 |
Зависимость имеет вид:
2.2. Плотность жидкой смеси на линии насыщения определяется по формуле: , где:
определяется по уравнению Ганна и Ямады: , где:
,
и Г – функции приведенной температуры.
Для диапазона температур :
;
Для диапазона :
;
Для диапазона :
; .
Для диапазона температур от 298К до 573К имеем:
Т,К |
Тr |
Г |
Vr 0 |
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 |
Зависимость имеет вид:
3. Энтальпию, энтропию, теплоемкость смеси при заданной температуре 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 |
График зависимости имеет вид:
4. Энтальпию испарения смеси в стандартном состоянии и при давлении, отличном от 1атм рассчитываем по уравнению Менделеева-Клапейрона:
, где:
при стандартном состоянии =1, при давлении, отличном от 1атм для определения используется выражение:
;
;
;
,
давление насыщенных паров определяется по корреляции Ли-Кесслера:
, где:
,. .
,
,.
Для диапазона температур от 298 до 573К:
Т, К |
Тr |
f(0) |
f(1) |
Pvpr |
∆Zv |
∆v H, Дж/моль |
∆v H0 , Дж/моль |
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 |
Зависимость имеет вид:
5. «Кажущаяся» стандартная энтальпия образования смеси в жидком состоянии определяется по уравнению: , где:
- энтальпия парообразования смеси при стандартном давлении и температуре 298К,
- энтальпия образования газовой смеси при стандартном давлении и температуре 298К, определена по правилу аддитивности. Энтальпии образования чистых компонентов смеси определены методом Бенсона.
6. Вязкость смеси при атмосферном давлении рассчитываем методом Голубева:
т.к. , то:
,
где
При расчете вязкости для высоких давлений используем корреляцию для смесей неполярных газов, т.к. большинство газов в смеси – неполярные:
где: , приведенную плотность смеси определяем по формуле: , где V рассчитывается по уравнению Менделеева-Клапейрона: , а псевдокритический коэф. сжимаемости определяется из разложения Питцера с помощью таблиц Ли-Кесслера: . Условие применимости для : удовлетворяется.
8. Теплопроводность смеси при стандартном давлении рассчитаем с помощью корреляции Мисика-Тодоса:
где: - теплоемкость смеси в стандартном состоянии при температуре 730К,
.
Теплопроводность при высоком давлении и определяется по формуле: