Реферат: Критический объем и плотность веществ их прогнозирование

В пунктах 5-10 вычисления производились с округлением промежуточных величин на каждом шаге вычислений до значений, записанных в формулах. Далее приводятся результаты вычислений с точностью в 16 десятичных разрядов, и округление выполнено только при представлении окончательных величин.

11. Примем Psat = 3 атм. Повторим вычисления в пунктах 5-10. При 400 К и 3 атм объем пара составляет 9,878 л/моль, объем жидкости – 0,282 л/моль. Левая часть выражения (6.3) равна = 1,0515. Тождество не выполняется, но степень отклонения от него существенно уменьшилась.

12. Подбор давления насыщения следует продолжить. Теперь имеется два значения для левой части выражения (6.3) при соответствующих давлениях. Используя эти величины, можно оценить значение давления для следующего расчета путем линейной интерполяции.

= 1–(1–3)/(35,53–1,0515)·35,53 = 3,061 атм.

13. Повторим вычисления (пункты 5-12) для Psat = 3,061 атм. Получим:

= 9,658 л/моль; = 0,282 л/моль; = 0,473. Новое значение давления – 3,111 атм.

После 5 итераций, исключая расчет при Psat = 10 атм, имеем:

T = 400 K; P sat = 3,112 атм; = 9,480 л/моль; = 0,282 л/моль; = 8,7·10-5. Полученные значения давления и объемов жидкости и пара соответствуют условиям насыщения.

14. Результаты расчета для других температур приведены в табл. 6.3.

Таблица 6.3

Т , К	Psat , атм	, л/моль	, л/моль
400	3,112	0,282	9,480
500	9,888	0,322	3,235
600	22,328	0,410	1,322
640	29,127	0,515	0,850

15. Область метастабильных (пересыщенных) состояний пара и жидкости занимает пространство между бинодалью и спинодалью. Точки на изотермах, принадлежащие бинодали, определены выше, и их значения приведены в табл. 6.3.

Для определения конфигурации спинодали воспользуемся соотношением

,

т.е. условиями экстремальности для соответствующих точек изотермы. Далее продифференцируем уравнение Ван-дер-Ваальса по объему (при Т = const) и преобразуем полученное выражение к полиному по V. Получим кубическое уравнение (6.12), корни которого могут быть найдены изложенным выше способом (п.п. 5-9):

.(6.12)

16. Для 400 К имеем следующие значения коэффициентов уравнения (6.12):

= –[2·38,72/(0,08206·400)] = –2,3593;

= [4·38,72·0,2151/(0,08206·400)] = 1,0149;

= –[2·38,72·0,21512/(0,08206·400)] = –0,1092.

Коэффициенты приведенного кубического уравнения (6.5) соответственно равны:

= [3·1,0149 –(–2,3593)2]/3 = –0,8405;

= 2·(–2,3593)3/27 –(–2,3593·1,0149)/3 + (–0,1092) = –0,2838;

= (–0,8405/3)3 + (–0,2838/2)2 = –0,0019.

Значение D отрицательное, следовательно, уравнение имеет три действительных решения.

17. Найдем значения корней уравнения (6.12) при 400 К. Для этого выполним последовательно следующие вычисления:

= [–(–0,8405)3/27]1/2 = 0,1483;

= –(–0,2838)/(2·0,1483) = 0,9568;

= arccos (0,9568) = 0,2950 радиан;

= 2·(0,1483)1/3 cos(0,2950/3) = 1,0535;