Введение к работе
Актуальность темы исследования. Одной из важных задач физической химии является определение теоретического оптимального режима проведения реакции. Задача теоретической оптимизации состоит в возможности дать оценку предельным показателям реальных технологий. Одними из первых исследований по оптимизации каталитических процессов и аппаратов были работы академика Г.К. Борескова, чл.-корр. РАН М.Г. Слинько, Г.М. Островского, А.И. Бояринова, В.С. Бескова, А.В. Федотова, В.И. Быкова и др.
Следует отметить, что теоретическая оптимизация проводится на основе кинетической модели, точность которой существенно зависит от начальных значений входящих в нее параметров: констант скоростей реакций и энергий активации, определяемых из экспериментальных данных, то есть априори содержащих в себе некоторую погрешность. В работах С.И. Спивака, А.В. Балаева, А.В. Круглова, В.А. Холоднова, С.А. Мустафиной показано, что взаимосвязь между результатами оптимизации и неточностью кинетических параметров — задача недостаточно ясная даже для простых реакций. Актуальность ее решения определяется, прежде всего, запросами технологической практики, поскольку неопределенность в кинетических параметрах может привести не просто к некоторому отклонению от выбранного пути реализации технологии, а к принципиальному изменению всей стратегии ведения процесса. Важна и другая задача — необходимо оценить размеры участка неопределенности оптимального режима при отсутствии для кинетических параметров информации статистического характера. Отсюда возникает необходимость и актуальность проведения исследований по гарантированности прогноза этапа теоретической оптимизации в условиях интервальной неопределенности кинетических констант, построения соответствующих оценок изменения оптимальных решений.
При компьютерной обработке интервальных типов данных в существующих программах, например, в пакетах Matlab (модуль Intlab), Maple (модуль Intpakx) и Mathematica выявлены недостатки: несовместимость программного продукта с операционной системой, специфика языка программирования, увеличение времени расчета. В связи с этим возникает необходимость создания программного обеспечения для интервальных вычислений, обладающим удобным пользовательским интерфейсом и адаптированного к решению задач оптимизации в химической кинетике.
Цель работы. Разработка методики и математического обеспечения анализа влияния неопределенности в параметрах кинетической модели на результаты теоретической оптимизации на основе теории интервального анализа.
Задачи исследования.
-
Построение математических моделей химических реакций, учитывающих интервальную неопределенность начальных значений кинетических параметров.
-
Исследование области неопределенности с позиций качественной неизменности оптимального режима, выработка требований к предварительным исследованиям кинетики изучаемых химических реакций.
-
Количественная оценка границ изменения оптимальных условий при заданном изменении кинетических параметров математической модели.
-
Создание программно-методического обеспечения для анализа химических реакций с учетом интервальной неопределенности кинетических параметров и проведение вычислительного эксперимента поиска оптимальных температурных условий.
Личный вклад автора. Автором разработана методика вычислительного эксперимента по расчету оптимального температурного режима и чувствительности по кинетическим параметрам для систем реальной размерности, с произвольным видом нелинейности, отражающей зависимость скорости сложной реакции от концентрации и температуры.
Научная новизна. Построена методика анализа чувствительности оптимального температурного режима в условиях неопределенности кинетических параметров.
Показана необходимость проведения анализа неопределенности начальных данных реакций в два этапа: получения оптимальных условий постоянства стратегии ведения процесса и количественного определения границ температурных условий.
Путем вычислительного эксперимента найдены оптимальные температурные режимы для следующих процессов: окисления диоксида серы на платиновом катализаторе, пассивации никелевых катализаторов при упрощенной схеме, изомеризации ксилолов. Получены условия качественной неизменности оптимальных температурных режимов. Построены алгоритмы численного определения границ оптимальных температурных условий при вариации кинетических параметров внутри допустимой области их изменения.
Практическая значимость. Созданное программно-методическое обеспечение позволяет решать задачу анализа чувствительности оптимального температурного режима в условиях интервальной неопределенности начальных значений кинетических параметров.
Программный продукт зарегистрирован в объединенном фонде электронных ресурсов «Наука и образование» Института научной информации и мониторинга Российской академии образования (ОФЭРНиО ИНИМ РАО), Всероссийском научно-техническом информационном центре (ВНТИЦ).
Достоверность результатов обеспечивается использованием в качестве основы моделирования фундаментальных законов математики, химии, физики и выбором теоретически обоснованных численных методов, а также подтверждается удовлетворительным согласованием результатов проведенных расчетов с известными аналитическими решениями, экспериментальными данными и расчетами других исследователей.
Апробация работы. Основные положения работы докладывались и обсуждались на: XXII Международной научной конференции «Математические методы в технике и технологиях» (ММТТ-22, Псков, 2009); Международной научно-практической конференции «Роль классических университетов в формировании инновационной среды регионов. Фундаментальное естественно-научное образование — генерация знаний на базе научных исследований», проведенная в рамках празднования 100-летия Башкирского государственного университета (Уфа, 2009); Международной научной конференции «Дифференциальные уравнения и смежные проблемы» (г. Стерлитамак, 2008); IX Международной научной конференции «Дифференциальные уравнения и их приложения» (Саранск, 2010); Всероссийской научно- практической конференции «Обратные задачи в приложениях» (Бирск, 2008); Всероссийской молодежной конференции по квантовой и математической химии (Уфа, 2008); Всероссийской научно-практической конференции «Финансовая и актуарная математика» (Нефтекамск, 2009); XI, XII Всероссийском симпозиуме по прикладной и промышленной математике (Кисловодск, 2010; Дагомыс, 2010; Казань, 2011); Конкурсе научно-исследовательских работ аспирантов научных организаций и вузов Республики Башкортостан (диплом III степени)(Уфа, 2010); объединенном научном семинаре химического факультета и факультета математики и информационных технологий Башкирского государственного университета (руководители — профессор Спивак С.И., профессор Прочухан Ю.А., профессор Герчиков А.Я.); научных семинарах факультета математики и естественных наук Стерлитамакской государственной педагогической академии (руководители — профессор Мустафина С.А., профессор Кризский В.Н.).
Публикации. По теме диссертации опубликовано 15 статей и тезисов научных конференций, зарегистрирован 1 программный продукт. В совместных работах постановка задачи принадлежит профессору С.А. Мустафиной. Результаты, выносимые на защиту, принадлежат диссертанту.
Объем и структура работы. Диссертация состоит из введения, 4 глав, заключения. Полный объем составляет 129 страниц, включая 6 рисунков, 3 таблицы, библиографию.
