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

.(6.4)

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

; ;

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

;

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

,(6.5)

где ; .

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

.(6.6)

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

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

, .(6.7)

Тогда решениями приведенного уравнения (6.5) будут

;(6.8)

;(6.9)

,(6.10)

от которых заменой

(6.11)

снова можно перейти к решениям кубического уравнения (6.4).

3. Вычислим характеристические константы уравнения Ван-дер-Ваальса. Для удобства вычислений примем следующие единицы измерения: V - л/моль , P - атм, Т - К. Тогда R = 0,08206 л·атм/(моль·К);

 a = 27·0,082062·6502/(64·31)=38,72 л·атм;

 b = 0,08206·650/(8·31)=0,2151 л.

4. Давление насыщения находится методом последовательных приближений. В качестве первого приближения при Т = 400 К примем давление насыщения равным 10 атм.

5. Рассчитаем значения коэффициентов уравнения (6.4):

 = –(0,2151+0,08206·400/10) = – 3,4975;

 = 38,72/10 = 3,872;

 = – (38,72·0,2151/10) = – 0,8329.

6. Далее вычислим коэффициенты приведенного кубического уравнения (6.5) и значение дискриминанта D:

 = [3·3,872–(–3,4975)2]/3 = – 0,2055;

 = 2·(–3,4975)3/27–(–3,4975·3,872)/3+(–0,8329)=0,5121;

 = (–0,2055/3)3+(0,5121/2)2 = 0,0652.

Значение дискриминанта (D) получилось положительным, что говорит о единственном действительном решении уравнения (6.5). Следовательно, значение давления выбрано неверно.

7. Предположим, что давление насыщения равно 1 атм. Повторим вычисления в пунктах 5 и 6.

 = –(0,2151+0,08206·400/1) = –33,04;

 = 38,72/1 = 38,72;

 = –(38,72·0,2151/1) = –8,329;

 =[3·38,72 –(–33,04)2]/3 = –325,2;

 = 2·(–33,04)3/27 –(–33,04·38,72)/3+(–8,329) = –2254;

 = (–325,2/3)3+(–2254/2)2 = –3632.

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

8. Найдем эти решения, но прежде вычислим вспомогательные величины и

 = [–(–325,2)3/27]1/2 = 1129;

= –(–2254)/(2·1129) = 0,9982;

= arccos (0,9982) = 0,0600 радиан;

 = 2·(1129)1/3·cos(0,0600/3) = 20,82;

= 2·(1129)1/3 cos(0,0600/3 + 2·3,14/3) = –10,75;

= 2·(1129)1/3 cos (0,0600/3 + 4·3,14/3) = –10,09.


Информация о работе «Критический объем и плотность веществ, их прогнозирование»
Раздел: Химия
Количество знаков с пробелами: 16804
Количество таблиц: 2
Количество изображений: 0

Похожие работы

Скачать
9601
2
1

... -7 -47,5 Двойная связь между атомами углерода 7 – Тройная связь между атомами углерода 14 –   Неаддитивный метод Тина и Каллуса Величина мольного объема жидкости при нормальной температуре кипения представлена в качестве функции критического объема: ,(6.13) где  и  выражены в см3/моль. Это простое соотношение хорошо прогнозирует  для органических чистых жидкостей, погрешность ...

Скачать
36551
48
10

... Выпишем парциальные вклады для температуры, давления и объема: Группа к-во ΔT ΔP СН2 (цикл) 4 0,04 -0,0112 Сумма 4 0,04 -0,0112 Критическая температура.   Критическое давление. ; о-Ксилол, 1,2-диметилбензол Выпишем парциальные вклады для температуры, давления и объема: Группа кол-во ΔT ΔP CН3 2 0,0282 -0,0024 -СН=(цикл) 4 0, ...

Скачать
87973
33
5

... степени, 10% —второй. Таким образом, наибольшую опасность представляют бескар­касные здания без фундамента из местных материалов, жители которых могут серьезно пострадать. Прогнозирование и оценка обстановки при чрезвычайных ситуациях природного характера осуществляется с использованием специальных методик и рекомендаций , часть некоторых из них изложена в данной методичке и в учебном пособии « ...

Скачать
34451
46
2

... Задание №5 Для первого соединения рассчитать плотность вещества при температуре 730 К и давлении 100 бар. Определить фазовое состояние компонента. Для определения плотности вещества воспользуемся методом прогнозирования плотности индивидуальных веществ с использованием коэффициента сжимаемости. где -плотность вещества; М- молярная масса; V-объем. Для данного вещества найдем коэффициент ...

0 комментариев


Наверх