Ufa, Russian Federation
student
Ufa, Russian Federation
Rassmatrivaetsya vopros vliyaniya neopredelennosti v parametrah matematicheskoy modeli na prakticheskie rezul'taty; proveden analiz OTR pri var'irovanii znacheniya energiy aktivacii; provodilsya vychislitel'nyy eksperiment pri raschete OTR dlya reakcii okisleniya monooksida ugleroda; polucheny eksperimental'nye dannye; postroen optimal'nyy temperaturnyy rezhim i issledovana zadacha analiza chuvstvitel'nosti optimal'nogo temperaturnogo rezhima dlya reakcii polucheniya okisi etilena v reaktore ideal'nogo vytesneniya s nepodvizhnym sloem katalizatora
modelirovanie, optimal'nyy temperaturnyy rezhim, vychislitel'nyy eksperiment, eksperimental'noe izmerenie, himicheskie reakcii, analiz chuvstvitel'nosti, intervaly, kineticheskie parametry, differencial'nye uravneniya
Всякое экспериментальное измерение включает в себя некоторую неопределенность, количественно выражаемую интервалом. Таким образом, значение измерения является объектом интервального анализа. Отсюда, если все вычисления, осуществляемые ранее для получения той или иной инварианты, повторить в интервальной арифметике, то результирующий интервал обладает двумя важными свойствами: гарантированности и независимости от закона распределения случайной величины.
В данной работе рассматривается вопрос влияния неопределенности в параметрах математической модели на практические результаты. Рассмотрение примеров приводит к следующему факту: знание энергий активации уже достаточно для прогноза оптимального температурного режима. В то же время следует заметить, что явное математические выражения в виде формул удаётся получить только для сравнительно простых случаях. В реальных ситуациях число дифференциальных уравнений математического описания может достигать десяти и более. Вид правых частей может также усложняться по сравнению с законом действующих масс. В этих ситуациях получить аналитические выражения для ОТР и условий качественной неизменности становится невозможным. И вопрос о чувствительности ОТР к вариации кинетических параметров становится задачей вычислительного эксперимента [1].
Вычислительный эксперимент проводился при расчете ОТР для реакции окисления монооксида углерода [2]. Обработка данных с целью решения обратной кинетической задачи проводилась на основе кинетического уравнения вида:
|
|
(1) |
х – концентрация этана или СО, n - порядки реакции по этану или СО.
Для таких функций найдены численные значения кинетических параметров: предэкспонентов, энергий активации и порядков реакции по ключевым компонентам.
Все экспериментальные данные получены в изотермических условиях при постоянном давлении. Математическое описание окисления этана и монооксида углерода в изотермическом реакторе идеального смешения представляется системой уравнений [3].
Входящие в правые части уравнений составляющие (Fi) равны:
|
F1 = - 0.5w1 – 3.5 w2; |
(2) |
|
F2 = - w1; |
(3) |
|
F3 = - w2; |
(4) |
|
F4 = w1 + 2 w2; |
(5) |
|
F5 = 3w2; |
(6) |
|
FN = - 0.5w1 + 0.5w2. |
(7) |
Система дифференциальных уравнений, построенная в соответствии с кинетической схемой, имеет вид:
|
|
(8) |
|
|
(9) |
|
|
(10) |
|
|
(11) |
|
|
(12) |
|
|
(13) |
|
|
(14) |
|
|
(15) |
|
|
(16) |
|
|
(17) |
|
|
(18) |
Система дифференциальных уравнений решалась с изменяемыми начальными данными.
В работе проведен анализ ОТР при варьировании значения энергий активации. Выбирается значение энергии активации, полученное для каждого из трех катализаторов. Это значение является базовым. Изменяя базовое значение в пределах ±25%, получаем 7 значений для каждой энергии активации. В результате расчета получены 63 графика: для 3 катализаторов, 3 наборов начальных данных и 7 значений энергий активации [4].
Вариация энергии активации Е1 в пределах 25% приводит к изменению температуры не более 20-25%. Иначе говоря, вариация энергии активации в некоторых пределах (в данном случае, до 25%) не меняет вид оптимального температурного режима. В данной работе такая ситуация наблюдается почти для всех наборов начальных данных и всех катализаторов.
Хотя на практике могут возникнуть такие ситуации, когда изменение энергии активации в небольших пределах приводит к достаточно большому изменению температуры. То есть, в определенных условиях может наблюдаться повышенная чувствительность к изменению кинетических параметров [5].
Анализ проведённого вычислительного эксперимента позволяет определить точность в кинетических константах, необходимую для надежного определения оптимальных характеристик технологического процесса.
Из всего вышесказанного можно сделать выводы, что в случае дробно-рациональной зависимости скорости реакции от концентраций получены явные математические выражения для оптимального температурного режима и чувствительности по кинетическим константам. Построен оптимальный температурный режим и исследована задача анализа чувствительности оптимального температурного режима для реакции получения окиси этилена в реакторе идеального вытеснения с неподвижным слоем катализатора. Модифицирован комплекс прикладных программ для расчета химико-технологических процессов с изменяющимися свойствами реакционной среды в реакторах с неподвижным слоем катализатора.
Уровень сервиса разработанного математического обеспечения позволяет использовать его при решении задачи расчёта и анализа оптимального температурного режима для широкого класса каталитических процессов. Разработанное математическое обеспечение использовано при расчете и анализе оптимального температурного режима для промышленно значимой реакции – очистки углеводородсодержащих газов от монооксида углерода в присутствии металлических катализаторов. Анализ вычислительного эксперимента по анализу чувствительности оптимального температурного режима к вариации кинетических констант позволяет оценить область вариации оптимального температурного режима при вариации кинетических констант в пределах, определяемых величиной погрешности измерений. Проведённые расчеты позволяют оценить надежность технологических решений, следующих из результатов математического моделирования соответствующих процессов.
1. Modelirovanie ob'ektov himicheskoy tehnologii // Iremadze E.O., Mustafina S.A., Spivak S.I. - Ufa, 2007.
2. Optimizaciya nestacionarnyh processov v usloviyah neopredelennosti kineticheskih parametrov // Mustafina S.A., Iremadze E.O., Spivak S.I. Trudy Srednevolzhskogo matematicheskogo obschestva. - 2007. - T. 9. - № 1. - S. 218-221.
3. Primenenie teorii chuvstvitel'nosti dlya optimizacii kataliticheskih processov i reaktorov // Mustafina S.A., Iremadze E.O., Spivak S.I. Obozrenie prikladnoy i promyshlennoy matematiki. - 2005. - T. 12. - № 2. - S. 447.
4. Chuvstvitel'nost' optimal'nogo temperaturnogo rezhima kataliticheskogo processa k variacii kineticheskih konstant // Iremadze E.O. Dissertaciya na soiskanie uchenoy stepeni kandidata himicheskih nauk / Ufa, 2002.
5. Neopredelennost' v kineticheskih konstantah i raschet optimal'noy temperatury / Iremadze E.O., Mustafina S.A., Spivak S.I. Matematicheskoe modelirovanie. - 2000. - T. 12. - № 3. - S. 21-22.
