^{ПРОГНОЗИРОВАНИЕ} КРИТИЧЕСКОГО ОБЪЕМА

Для массовых расчетов может быть рекомендован метод Лидерсена, возможности и точность которого, однако, не следует переоценивать. По методу Лидерсена критический объем рассчитывается с использованием корреляции:

,(5.17)

где _v - парциальные вклады, значения которых, выраженные в кубических см³/моль, приведены в табл. 5.2. Расчет достаточно прост и не требует дополнительного комментария.

ПРОГНОЗИРОВАНИЕ АЦЕНТРИЧЕСКОГО ФАКТОРА

Фактор ацентричности был предложен в 1955 г. Питцером в качестве коррелирующего параметра, характеризующего ацентричность, или несферичность молекулы. Анализируя зависимость приведенного давления насыщенного пара различных веществ от приведенной температуры, Питцер с сотрудниками установили, что для аргона, криптона, ксенона, азота, кислорода, окиси углерода, метана и некоторых других веществ эта зависимость описывается практически одним уравнением. Однако расширение этого списка соединениями других классов дает серию практически прямых линий, наклоны которых различаются. Питцер и др. приняли приведенное давление насыщенного пара при определенной приведенной температуре в качестве характеристики вещества. При этих температурах приведенное давление инертных газов, выбранных в качестве простого вещества, составляет примерно 0,1. На основании этого наблюдения было сформулировано определение нового параметра - ацентрического фактора  как описывающего отклонение значения приведенного давления пара для определенного вещества от приведенного давления пара вещества сравнения в следующем виде:

(при T_r=0,7),(5.18)

где - давление насыщенного пара вещества при приведенной температуре T_r =0,7.

По определению Питцера ацентрический фактор является “мерой отклонения функций межмолекулярного потенциала от функций межмолекулярного потенциала сферических молекул вещества сравнения”. Значение  = 0 соответствует сферической симметрии в разреженном газе. Отклонения от поведения, характерного для простого вещества, очевидны, если  > 0. Для одноатомных газов ацентрический фактор близок к нулю. Для метана он еще очень мал. Однако для углеводородов с высокой молекулярной массой значение  возрастает и резко увеличивается с ростом полярности молекул.

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

Здесь следует заметить, что опыт нашей работы позволяет заключить, что приведенное выше ограничение является излишне категоричным. При соблюдении определенных условий корреляции с  могут использоваться и применительно к названным группам органических веществ.

Значения ацентрического фактора для многих веществ вычислены на основе лучших экспериментальных данных по упругостям паров, T_c и P_c соединений и содержатся в Приложении.

При отсутствии сведений об  для его прогнозирования могут использоваться:

• уравнение Эдмистера

;(5.19)

• уравнение Ли-Кеслера

,(5.20)

• уравнение Амброуза-Уолтона

,(5.21)

где - критическое давление, выраженное в физических атмосферах;

 = - приведенная нормальная температура кипения вещества;

- нормальная температура кипения вещества в градусах Кельвина;

- критическая температура в градусах Кельвина.

f⁽⁰⁾, f⁽¹⁾ – определены в описании метода Амброуза-Уолтона (раздел 7.3)

Завершая рассмотрение материала по критическим свойствам и критериям подобия, остановимся еще на одном важном и общем вопросе. Он касается критериев подобия. В настоящее время их предложено довольно много, мы познакомились с одним из них - ацентрическим фактором. В разд. 7 рассматривается еще один критерий подобия - и коэффициент Риделя. Оба критерия применяются весьма широко. Тем не менее универсальных подходов к выбору того или иного критерия подобия пока не создано, а значит, работы в этом направлении будут продолжены. Мы считаем целесообразным повторить те требования, которые перечислены Уэйлесом в его монографии [19] и относятся к дополнительным параметрам или критериям подобия:

• Эти параметры должны соотноситься с молекулярной структурой и электростатическими свойствами молекулы.

• Их можно определить при минимальном количестве экспериментальных данных.

• Критические свойства не должны оказывать непосредственное воздействие на их значения.

• При оценке этих параметров надо избегать использования данных о P-V-T , так как в противном случае теряется смысл приведенного уравнения.

Дополнительные параметры должны быть функцией температуры, предпочтительно приведенной.

Можно соглашаться или нет с перечисленными требованиями, но совершенно очевидно, что всему их комплексу не отвечает ни ацентрический фактор, ни критерий Риделя. Мало того, нам представляется ясным, что одной из причин успеха в их применении является именно согласованность их величин с критическими параметрами и P-T данными. В качестве носителя связи с P-T данными выступает температура кипения при одном из давлений, чаще при атмосферном.

Таким образом, развитие методов прогнозирования потребует, вероятно, и уточнения требований к критериям подобия.

6. ПРОГНОЗИРОВАНИЕ плотности газа и жидкости [6, 17-18]

Перед тем как перейти к прогнозированию, следует напомнить, что в зависимости от принятых температуры и давления вещество может находиться либо в насыщенном, либо в ненасыщенном состоянии. Давление над насыщенной жидкостью равно давлению ее насыщенного пара при данной температуре. Давление над ненасыщенной, переохлажденной или сжатой жидкостью больше давления ее насыщенного пара при избранной для расчета температуре. Для каждой из названных областей P-V-T пространства существуют самостоятельные подходы к прогнозированию плотности.

Прогнозирование плотности индивидуальных веществ с использованием коэффициента сжимаемости

Пример 6.1

Для изобутилбензола, имеющего критическую температуру 650 К, критическое давление 31 атм и ацентрический фактор 0,378, рассчитать с использованием таблиц Ли-Кеслера (табл. 4.6, 4.7):

• коэффициент сжимаемости при 500, 657 и 1170 К и давлении 1-300 атм,

• плотность при 500, 657 и 1170 К и давлении 1-300 атм;

дать графические зависимости:

• коэффициента сжимаемости от давления при указанных температурах,

• плотности от давления при указанных температурах.

Решение

Используем разложение Питцера (уравн. 4.34) и табл. 4.6, 4.7 для коэффициента сжимаемости.

1. Вычислим значения приведенных температур:

= 500/600 =0,769; = 657/650 =1,01; = 1170/650 =1,80.

1. Вычислим значения приведенных давлений:

= 1/31 =0,03226; = 300/31 =9,677.

Поскольку диапазон интересующих приведенных давлений совпадает с диапазоном , рассмотренных Ли-Кеслером, используем информацию о и для дискретных значений , представленных в табл. 4.6, 4.7.

Каждое из значений и получено линейной интерполяцией по температуре. Так, при 500 К ( = 0,769) и = 0,010 для имеем

(0,9935-0,9922)/(0,80-0,75)·(0,769-0,75)+0,9922 = 0,9927.

Прогнозирование плотности насыщенных жидкости и пара с использованием уравнений состояния вещества

Нахождение условий насыщения из уравнений состояния представляет собой достаточно сложную задачу, решение которой зачастую невозможно без привлечения вычислительной техники и специального программного обеспечения. Для простых уравнений состояния, таких как уравнение Ван-дер-Ваальса, эта задача может быть решена путем несложных вычислений. Однако необходимо помнить, что на практике при помощи уравнения Ван-дер-Ваальса можно лишь качественно оценить состояние насыщения. Для более точного представления насыщения разработаны другие уравнения состояния и специальные методы.

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

Пример 6.3

Для изобутилбензола при температурах 400, 500, 600 и 640 К, используя уравнение Ван-дер-Ваальса, рассчитать давление пара и объемы насыщения жидкости и пара. Определить также области метастабильных состояний пара и жидкости при указанных температурах. Критическая температура равна 650 К, критическое давление - 31 атм.

Решение

1. Запишем принцип Максвелла:

Площадь = .(6.1)

Выразим из уравнения Ван-дер-Ваальса значение давления и подставим его в подинтегральное выражение. Получим

.(6.2)

В данном случае имеется возможность найти аналитическое решение определенного интеграла

.(6.3)

Теперь задача сводится к отысканию значения P ^sat, при котором выражение 6.3 обратится в тождество. При его нахождении нам потребуется неоднократно определять значения объемов жидкости и пара для заданного P, т.е. находить решения (корни) кубического уравнения.

1. Перепишем уравнение Ван-дер-Ваальса в виде полинома по объему

.(6.4)

Корни данного уравнения можно найти, воспользовавшись формулами Кардано. Для этого перейдем к приведенному виду кубического уравнения, выполнив следующие преобразования. Обозначим коэффициенты в уравнении (6.4) через

; ;

и сделаем замену неизвестного V на Y:

;

тогда уравнение (6.4) примет приведенный вид

,(6.5)

где ; .

Число действительных решений кубического уравнения зависит от знака дискриминанта

.(6.6)

Если D > 0, то уравнение имеет одно действительное решение; если D < 0, то - три действительных решения; и если D = 0, то уравнение имеет либо два действительных решения, одно из которых двукратное, либо одно действительное трехкратное решение (последнее в случае p = q = 0).