Инструкция моделирование химических и фазовых равновесий при высоких температурах, страница 17

Каждая таблица промежуточных результатов  содержит  следующие данные: номер итерации, максимальную относительную погрешность вычисление (eps), подготовки  матрицы  коэффициентов (tau1), время решения системы линейных уравнений (tau2), время формирования начального приближения для следующей итерации (tau3), давление (p, МПа), температуру (T, K), удельный объем (v, m**3/кг), энтропия (S, кДж/(кг*К),  энтальпия (I,кДж/кг), внутреннюю энергию (U, КДж/кг), значения неопределенных множителей Лагранжа (L[эл]), (числа молей вещества в каждой из фаз системы (моль/кг), содержание  компонентов в газовой фазе и конденсированных растворах (моль/кг).

Для p, T, v, S, I, U, а также для мольного состава фаз на печать выводится по два значения.  Первое - результат итерации, второй  -  величина,  выбранная  в  программе  в  качестве начального приближения для следующего шага.

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

Для контроля предусмотрена возможность вывода -   печать  тех значений термодинамические свойств индивидуальных веществ, которые использовались в текущем расчете.  Чтобы реализовать ее,  в список директив  следует  включить директивы <Contr.  Она вызывает печать таблицы значений (p)i,  (»)i,  (I)i, (S)i после завершения каждого цикла итераций.

Парциальные давления  компонентов,  выведенные  в размерности физ. тм,  позволяют проверить значение любой константы  равновесия (Kp).  Приведенная  свободная энергия Гиббса и энтальпия компонентов,  включают   в   себя   величину   энтальпии   образования: (I)i=(dH)i[298] + (H)i[T] - (H)i[298]. Размерность термодинамических свойств определяется директивой ввода <PrSI или <Prte  и может быть либо Дж-моль-К, либо кал-моль-К .

Пример:            <inte<prte<prdom<prp<contr>

-ДЕК.90,

DATА  ДЕК.90, VAR = CONTROL,

v = 0.01, T = 1000,1500,

(100%AlCl3);

< Prdom >

p=2, i=0, ox=(O2[0]), fu=(H2[0]),                                      

alpha=1.2;

Исходными данными предусмотрено проведение серии расчетов с вариацией исходного состава:

(N0%O2[0])

(1%H2[0])

N0= .9524E+01

Приняты к исполнению директивы:  <INSI   <PRSI   <LIST   <PRDOM  <ION <PRM

Вариант расчета:   000                    *                    (1)

Заданы параметры:

P = 20000+01

I = 00000-00

Состав рабочего тела  моль/кг     O= 56564+02   H= 94273+02

**   АСТРА.4/pc   **

Московский Государственный Технический университет им. Баумана

23 фев.96

Xapaктepиcтики paвновесия – СИ

P= 2OOOO+O1    T= 34636+O4  V= 842O8+OO  S= 15297+O2  I= OOOOO-OO  U=-16842+O4

M= 5848O+O2   Cp= 29633+O1  k= 11963+O1  Cp"= 1O959+O k"= 11822+O1 A= 13763+O4

Mu= 94O57-O4  Lt= 4258O+OO  Lt"= 23833+O1  

MM= 171OO+O2 Cp.г= 29651+O1  k.г= 11963+O1

MM.г= 171OO+O2  R.г= 48624+O3    Z= OOOOO-OO   Пл= OOOOO-OO   Bm= 61833+OO

Coдepжaние компонeнтoв - моль/кг

O  17261+O1       -1*O  19496-O6       O2  5O494+O1   +1*O2 1O185-O7 

-1*O2  6O868-O7         O3  26O84-O6        H  2O634+O1   -1*H  16157-O7  

H2  49193+O1         OH  7O57O+O1    -1*OH  98518-O6    HO2  15478-O1

H2O  37647+O2     +1*H2O  32393-O7      H2O2 15998-O2 +1*H3O  29316-O5  эл.газ  17143-O5

МГТУ * АСТРА.4/pc 1996

рис. 7.

6.  ПРИМЕРЫ ПРОБЛЕМ, РАЗРЕШИМЫХ С ПОМОЩЬЮ

ПРОГРАММНОГО КОМПЛЕКСА  Астра.4/рс а)  Расчет состава продуктов диссоциации, термодинамических и теплофизических свойств газовой смеси,  состоящей из 20% (CO2) и 80% (O2) при давлении 0.01 МПа и в диапазоне температуры 2000-4000 К с шагом 500 К.

Результаты: