Критический объем и плотность веществ их прогнозирование (стр. 3 из 3)

После 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;

= 2·(0,1483)1/3 cos(0,2950/3 + 2·3,14/3) = –0,6159;

= 2·(0,1483)1/3 cos(0,2950/3 + 4·3,14/3) = –0,4388;

= 1,0535 –(–2,3593/3) = 1,840 л/моль;

= –0,6159 –(–2,3593/3) = 0,171 л/моль;

= –0,4388 –(–2,3593/3) = 0,348 л/моль.

Наибольший корень

= 1,840 л/моль соответствует максимуму на изотерме 400 К и ограничивает метастабильные состояния пара слева. Корень , равный 0,171 л/моль, не имеет физического толкования, поскольку его значение меньше параметра b в уравнении Ван-дер-Ваальса. И, наконец, корень соответствует минимуму на изотерме 400 К и отделяет слева область пересыщенной жидкости от абсолютно неустойчивых состояний.

18. Давление в системе при соответствующем объеме пересыщенного пара (

) и пересыщенной жидкости ( ) находится из уравнения Ван-дер-Ваальса подстановкой в него требуемых значений температуры и объема.

= (0,08206·400)/(1,840–0,215)–38,72/1,8402 = 8,763 атм;

= (0,08206·400)/(0,348–0,215)–38,72/0,3482 = –72,928 атм.

19. Результаты расчета для прочих температур приведены в табл. 6.4.

Т, К	, атм	, л/моль	, атм	, л/моль
400	-72,928	0,348	8,763	1,840
500	-20,124	0,397	14,913	1,324
600	17,803	0,482	24,103	0,929
640	28,798	0,563	29,347	0,750