Алгоритм решения этой системы уравнений был реализован в виде
компьютерной программы
,
позволяющей исследовать системы про
-
извольного химического состава благодаря подключению базы дан
-
ных
,
содержащей термодинамические и термохимические свойства
индивидуальных веществ
,
а также параметры потенциальной функции
Леннард
–
Джонса
[12].
Программа позволяет исследовать как адиаба
-
тически изолированные
,
так и изотермические системы
,
находящиеся
при постоянном объеме или давлении
.
Установление степени достоверности результатов моделирования
с использованием предлагаемого подхода возможно только путем со
-
поставления их с известными экспериментальными данными
.
Однако
среди опубликованных сведений о непосредственных измерениях кон
-
центраций веществ в ходе химических превращений практически не
-
возможно найти примера
,
когда замеры концентраций одновременно
проводились для достаточного числа соединений
,
чтобы можно было
удостовериться в соблюдении баланса массы химических элементов
.
Кроме того
,
опубликованные сведения практически всегда охватывают
только часть временн
´
ого интервала химических превращений
.
По этой причине в настоящей работе сделана попытка определить
степень практической пригодности разработанного метода путем срав
-
нения с результатами расчетов
,
полученными по какой
-
либо признан
-
ной и проверенной кинетической модели
.
В качестве такой была выбра
-
на модель горения углеводородных топлив
GRI-Mech,
опубликованная
в сети
Internet
на странице университета Беркли
[16].
Там же приведена
информация о ее экспериментальной проверке
.
Необходимые расчеты
по этой кинетической модели проводились с использованием компью
-
терной программы
Kintecus [17].
На рис
. 1
представлены результаты моделирования кинетического
процесса химических превращений в системе
0,01%CH
4
+0,02%O
2
+
+99,97%Ar
для условий
T
= 2000
K
и
p
= 0
,
0984
МПа по описанно
-
му выше методу
,
а на рис
. 2
приведены концентрационные кривые
,
вычисленные с использованием формальной кинетической модели
GRI-Mech.
Сопоставление результатов по абсолютным значениям кон
-
центраций
(
моль
/
см
3
)
и по характеру их последовательного образова
-
ния позволяет сделать вывод о принципиальной работоспособности
метода
.
При этом
,
конечно
,
необходимо принять во внимание тот факт
,
что в описываемой термодинамической модели время можно рассма
-
тривать только как условный параметр
,
однонаправленный с реальным
временем
.
В процессе моделирования использовались термодинамические и
термохимические свойства индивидуальных веществ
,
опубликован
-
ные в работах
[4, 13].
Параметры потенциальной функции Леннард
–
ISSN 1812-3368.
Вестник МГТУ им
.
Н
.
Э
.
Баумана
.
Сер
. “
Естественные науки
”. 2005.
№
3
33
