автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.18, диссертация на тему:Методы, алгоритмы и программы моделирования кинетики химических и биохимических процессов с использованием интервального анализа

005018184

На правах рукописи

ХаибШУСр

ХАИДАРОВ АНДРЕЙ ГЕННАДЬЕВИЧ

МЕТОДЫ, АЛГОРИТМЫ И ПРОГРАММЫ МОДЕЛИРОВАНИЯ КИНЕТИКИ ХИМИЧЕСКИХ И БИОХИМИЧЕСКИХ ПРОЦЕССОВ С ИСПОЛЬЗОВАНИЕМ

ИНТЕРВАЛЬНОГО АНАЛИЗА

05.13.18 - Математическое моделирование, численные методы и комплексы программ

АВТОРЕФЕРАТ ДИССЕРТАЦИИ НА СОИСКАНИЕ УЧЕНОЙ СТЕПЕНИ КАНДИДАТА ТЕХНИЧЕСКИХ НАУК

1 9 ДПР 2012

005018184

На правах рукописи

ХАЙДАРОВ АНДРЕЙ ГЕННАДЬЕВИЧ

МЕТОДЫ, АЛГОРИТМЫ И ПРОГРАММЫ МОДЕЛИРОВАНИЯ КИНЕТИКИ ХИМИЧЕСКИХ И БИОХИМИЧЕСКИХ ПРОЦЕССОВ С ИСПОЛЬЗОВАНИЕМ

ИНТЕРВАЛЬНОГО АНАЛИЗА

05.13.18 - Математическое моделирование, численные методы и комплексы программ

АВТОРЕФЕРАТ ДИССЕРТАЦИИ НА СОИСКАНИЕ УЧЕНОЙ СТЕПЕНИ КАНДИДАТА ТЕХНИЧЕСКИХ НАУК

Работа выполнена в федеральном государственном бюджетном образовательном учреждении высшего профессионального образования «Санкт-Петербургский государственный технологический институт (технический университет)»

Научный руководитель -заслуженный работник высшей школы РФ, доктор технических наук, профессор, заведующий кафедры математического моделирования и оптимизации химико-технологических процессов, СПбГТИ(ТУ),

Холодное

Владислав Алексеевич

Официальные оппоненты -заслуженный деятель науки РФ, доктор технических наук, профессор, заведующий кафедры машины и аппараты химических производств, СПбГТИ(ТУ),

Веригин

Александр Николаевич

Кандидат технических наук, доцент кафедры процессов и аппаратов химической технологии, Санкт-Петербургская государственная химико-фармацевтическая академия,

Маркова

Алла Валентиновна

Ведущая организация - Институт информатики Академии Наук Республики Татарстан, г. Казань.

Защита состоится « 2£» апрглЛ г. в /4 час., ауд. ££ на заседании диссертационного совета Д 212.230.03 при федеральном государственном бюджетном образовательном учреждении высшего профессионального образования «Санкт-Петербургский государственный технологический институт (технический университет)»

L ?

С диссертацией можно ознакомиться в фундаментальной библиотеке института"' Г

Отзывы на автореферат в двух экземплярах, заверенные печатью, просим направлять по адресу: 190013, Санкт-Петербург, Московский пр. д.26, Санкт-Петербургский Государственный технологический институт (технический университет), Ученый совет-тел. 494-93-75; факс: 712-77-91; e-mail: dissovet@technolog.edu.ru.

Автореферат разослан « 23» миыт 2012 г.

Ученый секретарь /'

диссертационного совета y^Jl^ В.И. Халимон

Общая характеристика работы

Актуальность работы

Компьютерное моделирование химических реакторов доказало свою актуальность и перспективность, поскольку позволяет выполнять расчет материальных, тепловых балансов и осуществлять поиск наилучших режимов функционирования.

Решающее значение при выборе условий проведения химико-технологических процессов имеют вопросы определения скорости химических превращений, изучаемые химической кинетикой. Кинетические уравнения, содержащие необходимую информацию об основных закономерностях химических превращений, являются первоосновой математической модели химического реактора.

В настоящее время накоплен значительный опыт по решению традиционных для химической кинетики прямых и обратных задач.

Следует отметить большой вклад в разработку программных комплексов для исследования кинетических моделей Ю.В. Шарикова, В.И. Коробова, В.Ф. Очкова и др. В этих работах рассматривается традиционный подход к моделированию кинетики химических реакций с помощью компьютерных технологий и сетевых расчетов на Mathcad Calculation Server.

При решении прямых и обратных задач для сложных химических и биохимических реакций используются различные пакеты оптимизации GAMS, MATLAB, Mosaic и др., и применяются методы онлайн оптимизации: AMPL - сайт, BARON , Network-Enabled Optimization System Server.

В этих и ранних работах не рассматриваются задачи, связанные с неопределенностью кинетических параметров. Моделирование сложных химических и биохимических процессов основано на номинальных значениях кинетических параметров.

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

В научной литературе опубликованы работы, в которых рассмотрены вопросы влияния неопределенности исходной информации на результаты моделирования технических объектов (Grossmann I.E., Островский Г.М., Зиятдинов Н.Н, Дворецкий С.И., Егоров А.Ф., Холодное В.А. и др.). В работах Вощинина А.П., Нариньяни А.С., Шарого С.П., Шокина Ю.И., Левина В.И. и других при моделировании и оптимизации технических объектов в условиях неопределенности развиваются интервальные методы анализа.

Автором впервые при решении прямых и обратных задач химической кинетики для сложных химических и биохимических реакций используются интервальные методы.

В диссертационном исследовании речь идет об интервальной неопределенности кинетических параметров. Интервальная неопределенность предполагает естественную для задач химической технологии ограниченность кинетических параметров и их представление диапазоном возможных значений переменных.

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

Актуальность диссертационного исследования определяется тем, что в нем приведены разработанные автором методики, алгоритмы и программы для моделирования кинетики химических, биохимических реакций в химических реакторах в условиях интервальной неопределенности кинетических параметров с использованием интервальных методов. В работе впервые предложен интервальный метод вычисления чувствительности значений концентраций веществ к изменению констант скоростей реакций.

Применение интервальных методов, алгоритмов и программ позволит на новом уровне подойти к проблеме моделирования сложных химических и биохимических процессов в химических реакторах.

Цель диссертационной работы. Разработка совокупности взаимодополняющих методов, алгоритмов и программ для определения кинетических параметров сложных химических и биохимических процессов на основе экспериментальных данных с использованием интервальных методов.

Для достижения поставленной цели решены следующие задачи:

• проведен анализ существующих алгоритмов и программ для сложных химических и биохимических процессов на основе экспериментальных данных;

• обоснован выбор и разработан алгоритм адаптации метода огибающих для решения систем дифференциальных уравнений при решении прямой и обратной задач химической кинетики в интервальной постановке;

• адаптирован метод глобальной оптимизации в пакете программ Mathematica для решения обратных задач химической кинетики в интервальной постановке;

• создан пользовательский интерфейс для решения прямых и обратных задач кинетики химических и биохимических процессов в интервальной постановке;

• разработаны методы, алгоритмы и программы для решения прямых и обратных задач химической кинетики в интервальной постановке. Протестированы предложенные методы, алгоритмы и программы с использованием вычислительного и натурного экспериментов;

• предложена методика моделирования кинетики химических и биохимических процессов на основе разработанных алгоритмов;

• Проведен анализ чувствительности концентраций веществ к изменению кинетических параметров.

Объект исследования. Химические и биохимические процессы в химических реакторах, рассматриваемые в диссертационном исследовании.

Предмет исследования. Кинетические модели химических и биохимических процессов в условиях интервальной неопределенности кинетических параметров.

Методы исследования. В диссертации использованы методы: системного анализа, химической кинетики, математического моделирования и оптимизации химико-технологических процессов, компьютерной и интервальной математики.

Инструменты исследования. Интерактивные вычислительные системы Wolfram Mathematica, программирование пользовательского интерфейса.

Обоснованность научных результатов обеспечивается применением современных математических методов моделирования кинетики химических, биохимических реакторов, тестированием разработанных алгоритмов и программ и совпадение результатов вычислений с известными результатами натурных экспериментов.

Достоверность теоретических разработок подтверждена совпадением результатов вычислительного эксперимента на ПК с литературными и экспериментальными данными, что позволяет сделать вывод об эффективности разработанных алгоритмов и программ с использованием интервальных методов.

Научная новизна работы:

Проведено комплексное исследование кинетики важнейших химических и биохимических процессов в химических реакторах на основе интервального анализа с применением компьютерной технологии математического моделирования, натурного и вычислительного экспериментов:

1. Интервальные модели важнейших химико-технологических процессов:

а) тонкого органического синтеза: получение бета нафталин сульфокисло-ты, альфа нафталин сульфокислоты, алкилирования фенилацетонитрила в микрореакторе в интервальной постановке;

б) биохимического биокаталитического восстановления этилового эфира циклогексанонкарбоновой кислоты в этиловый эфир (111,28)-цис-2-гидроксициклогексанкарбоновой кислоты в присутствии Басскаготусея

в микроструктурном реакторе в интервальной постановке.

2. Численные методы:

а) интервальный метод для определения кинетических параметров химических и биохимических процессов по экспериментальным данным;

б) метод и алгоритм автоматизированного задания схемы химических реакций;

в) предложенный метод анализа чувствительности параметров математических моделей на основе интервального подхода;

г) модификация интервального метода огибающих для решения систем дифференциальных уравнений, описывающих химическую кинетику.

Основные положения, выносимые на защиту:

1. Методы и алгоритмы:

а) решения прямой и обратной задачи кинетики химических и биохимических реакций в интервальной постановке;

б) вычисления чувствительности концентраций веществ к изменению кинетических параметров;

в) алгоритм автоматического генерирования системы дифференциальных уравнений на основе различных подходов к составлению схемы химических реакций;

г) модификация алгоритма огибающих для решения систем дифференциальных уравнений с учетом погрешности экспериментальных данных.

2. Математические модели в интервальной постановке:

а) для исследования процессов тонкого органического синтеза;

б) для описания кинетики биохимических процессов в микрореакторах;

3. Разработанные программные комплексы:

а) для решения прямых, обратных задач химической кинетики, в том числе и в интервальной постановке;

б) для исследования чувствительности концентраций веществ к изменению констант скоростей реакций.

4. Результаты математического моделирования исследованных процессов.

Практическая значимость и реализация результатов работы. На основе теоретических результатов работы предложены и разработаны алгоритмы и программы для моделирования химических, биохимических реакций в реакторах. Для иллюстрации работоспособности предлагаемых методов и алгоритмов решены задачи по комплексному исследованию различных химических и биохимических процессов в реакторах.

Разработанные методы, алгоритмы и программы используются в учебном процессе в Санкт-Петербургском государственном технологическом институте (техническом университете), Казанском химико-технологическом университете, Березников-ском филиале Пермского государственного технического университета и в Тамбовском государственном техническом университете.

Работа была выполнена в рамках следующих проектов при поддержке Правительства РФ:

• Проект «Фундаментальные исследования закономерностей взаимосвязи основных параметров каталитических процессов в микроструктурных реакторах в условиях неопределённости исходной информации», реализуемый в рамках Аналитической ведомственной целевой программы "Развитие научного потенциала высшей школы" 2009-2010 гг.

• Проект «Интервальный подход к решению задач по математическому моделированию и оптимизации химических и биохимических процессов в микроструктурных реакторах в условиях неопределенности исходной информации», реализуемый в рамках Аналитической ведомственной целевой программы "Развитие научного потенциала высшей школы" 2009-2010 гг.

• Проект «Математические модели и методы интервального анализа для исследования и оптимизации биокаталитических процессов в микрореакгоре», реализуемый в рамках федеральной целевой программы «Научные и научно-педагогические кадры инновационной России» на 2009-2013 гг.

Апробация работы. Материалы диссертации докладывались и обсуждались на Всероссийских и международных научных конференциях: «Системы компьютерной математики и их приложения», Смоленск, СмолГУ, 2011 г.; «Математические методы в технике и технологиях ММТТ-24», Киев, КПИ, 2011 г.; «Технологии высокопроизводительных вычислений и компьютерного моделирования», Санкт-Петербург, СПбГТИ(ТУ), 2009 г.; «Актуальные проблемы химико-технологического образования»' РХТУ им. Д.И.Менделеева, 2011г. Получен сертификат участника ММТТ-24.

Публикации. По материалам диссертационной работы опубликовано 9 печатных работ, из которых 4 входит в рекомендуемый ВАК перечень.

Структура и объем диссертации. Диссертация состоит из введения, пяти глав

заключения, списка литературы. Работа изложена на 149 страницах основного текста '

содержит 55 рисунков, 17 таблиц, библиографический список литературы включает 121 наименований.

б

Основное содержание работы

Во введении обоснована актуальность решаемой научной проблемы, показана её практическая значимость, сформулированы цель и задачи исследования, приведены положения, выносимые на защиту, показана их научная новизна и практическая значимость.

1 Аналитический обзор научно-технической литературы

Моделирование химического реактора является важнейшей задачей при моделировании всего процесса в целом, Основой моделирования реактора является химическая кинетика. Кинетические уравнения являются первоосновой математической модели химического реактора и получаются при решении прямых и обратных задач химической кинетики.

Моделирование и оптимизация процессов на основе средних значений параметров математического описания не позволяет гарантировать режим функционирования, который может возникать в процессе эксплуатации. В связи с этим в последнее время используется перспективное и быстро развивающееся направление последних лет - математическое моделирование и оптимизация на основе интервальных данных, когда параметры математического описания задаются не числами, а интервалами.

На основе литературного обзора были поставлены задачи диссертационного исследования.

2 Алгоритмы численных методов и комплексов программ с использованием интервальных методов

Интервальный подход обусловлен погрешностью измерения значений концентраций при проведении кинетических исследований. В последнее время для учета погрешностей при обработке экспериментальных данных получил распространение интервальный анализ.

В главе 2 подробно рассмотрены существующие математические методы, на основании которых был разработан и далее реализован в среде МаШетайса модуль для интервального решения задач химической кинетики.

Рассмотрены следующие вопросы: основы интервальной арифметики, методы решения систем обыкновенных дифференциальных уравнений в интервальной постановке, алгоритмы нахождения интервалов констант скоростей реакций при решении об-1 ратной задачи.

На основе рассмотренных алгоритмов был выбран метод огибающих, который при сравнительно небольших вычислительных затратах позволяет построить сходящееся двустороннее решение с приемлемой точностью.

2.1 Алгоритм метода огибающих для решения прямой кинетической задачи

Для построения двустороннего решения системы дифференциальных уравнений используется модифицированный автором метод:

1. Задание исходных данных;

2. Вычисление множества решений для интервальных параметров для заданной системе дифференциальных кинетических уравнений;

3. Выбор А - шага разбиения исходного временного интервала;

4. Реализация последовательного цикла по исходному времени пребывания для его отдельных участков [ть Г, + /г]:

а) если временной участок требует построения огибающей, то огибающая рассчитывается для заданного участка;

б) если временной участок не требует построения огибающей, то в качестве границ двустороннего решения выбираются соответствующие частные решения;

5. Уменьшение шага Н и возврат к п.З в случае, если точность полученного двустороннего решения системы дифференциальных уравнений ниже заданной;

6. Вывод полученных результатов в табличном или графическом виде.

2.2 Постановка интервальной обратной кинетической задачи

Задачу интервального оценивания можно сформулировать так: по экспериментальным данным построить числовой интервал, относительно которого можно сказать, что внутри этого интервала находится оцениваемый параметр. Интервальное оценивание особенно необходимо при малом числе наблюдений, когда точечная оценка мало надежна.

Интервальная модель изменения концентраций веществ задается следующими соотношениями;

АС, - _ _

-¡г = /у №. ],[*,, к,],..., [к^, к„У,С, (г), С2 (г)..., С„ (г))

(1)

С,(0) = с; (У = 1,2.....п)

где С° - заданные начальные концентрации веществ;

п - количество участвующих в реакции веществ;

М, - верхняя и нижняя граница значений констант скоростей реакции для

''= 12.....т-Если этот интервал велик, то оценка, произведенная с его помощью, мало

пригодна для практики.

2.3 Алгоритм для определения интервалов констант скоростей реакции и проверка адекватности математической модели:

1. Задание верхних и нижних границ для расчета интервалов констант скоростей реакции. г

2. Задание целевой функции - суммы квадратов отклонений расчетных значений концентраций от экспериментальных данных с учетом весовых коэффициентов.

РЪ,к3,...*1)=£а,рс;-с>)\ (2)

где целевая функция, с;, С> - соответственно экспериментальные

и расчетные концентрации, - весовой коэффициент, /-количество определяемых констант, п - количество экспериментальных точек в одном наборе, т - количество наборов экспериментальных точек.

3. Задание ограничений по концентрациям веществ.

Задачу минимизации целевой функции можно сформулировать следующим образом: найти минимум целевой функции Zпри ограничениях Су(г) >0, j = 1,2,...п. При этом минимум целевой функции отыскивается внутри области ограничений с

использованием штрафной функции P(kt ,к2.....к,).

Z = F(kuk2.....к,) + Р(к„к2,...,к,), (3)

4. Определение средних значений констант к, в заданных интервалах, при которых целевая функция минимальна. Эти константы являются промежуточными для дальнейшего поиска.

5. Задание ограничений по константам скоростей реакций:

к,-е, <к' ^k,+ei (4)

где Л, - промежуточные значения констант скоростей реакции, st - отклонения от

промежуточного значения для г-ой константы, к'- значения констант;

6. Поиск значений е,, при которых обеспечивается попадание некоторого числа экспериментальных точек (N) внутрь расчетных интервалов концентраций для полученных значений интервальных констант:

N-A-M, (5)

где N- сумма всех экспериментальных точек, лежащих внутри расчетных интервалов для всех концентраций, А - коэффициент, задающий соотношение между количеством попадающих в интервал точек и шириной интервала, М-общее количество экспериментальных точек.

Так как для одного значения А существует бесчисленное множество наборов значений s,, нужно выбрать оптимальный набор, при котором ширина интервалов концентраций будет минимальна

7. Проверка адекватности по известным методикам с использованием критерия Фишера и скорректированного коэффициента детерминации R^ :

^ = (6) N-т-\

где N- число опытов, т - число подбираемых параметров математического описания, R - коэффициент детерминации.

На рисунке 1 представлена блок-схема алгоритма решения обратной задачи химической кинетики в интервальной постановке.

3 Практическое использование комплекса программ для моделирования кинетики химических процессов

В данной главе предлагается программный комплекс, который позволяет решать прямые и обратные задачи в традиционной и интервальной постановке с помощью Wolfram Mathematica. Алгоритмы для данного комплекса созданы с применением интервальных методов.

Комплекс алгоритмов и программ протестирован на процессе получения ^-нафталинсульфокислоты. Для удобства пользователей создан удобный интерфейс, а янные для расчета сохраняются в обычном текстовом файле.

Рисунок 1 - Блок- схема алгоритма для решения обратной задачи химической кинетики

в интервальной постановке

3.1 Исследование математической модели синтеза моносульфокислот нафталина с использованием интервального анализа

Математическая модель была построена по литературным данным и собственным исследованиям процесса синтеза сульфокислот нафталина.

Экспериментальные данные были получены на лабораторной установке периодического действия в изотермических условиях при различных температурах ведения процесса.

Математическое описание реакции сульфирования нафталина в периодическом реакторе в изотермических условиях для различных гипотез можно представить в виде:

S = S0 - (А + В + 2D) W = WI1 + A + B + 2D

где А, В, D- соответственно концентрации а- и /?- нафталин сульфокислот и ди-сульфокислоты, кмоль/м3;

N, S, W- концентрация нафталина, сульфирующего агента и воды, кмоль/м3;

индекс 0 относит соответствующую концентрацию к начальному моменту времени;

ki, кь кз, к4, к5- интервальные константы скоростей реакций (нижняя черта -нижняя граница, верхняя черта - верхняя граница) сульфирования нафталина в aw.fi-положения и реакции изомеризации а-нафталинсульфокислоты в /3-нафталинсульфокислоту, м3/(кмолъч).

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

С использованием предложенного во второй главе алгоритма по разработанной программе Interval Kinetic впервые были определены интервальные оценки кинетических параметров для данного процесса с учетом различной погрешности определения концентрация альфа-, бета- и дисульфокислот нафталина. Относительная погрешность определения концентраций альфа-, бета- сульфокислот 10%, а относительная погрешность измерения дисульфокислот 15%, что принималось во внимание в алгоритме решения обратной задачи и соответствует числу значащих цифр в значениях кинетических параметров в таблице 1.

Адекватность полученных моделей проверялась на основе вычисления критерия Фишера и на основе скорректированного коэффициента детерминации.

Полученные результаты приведены в таблице 1.

'¿А _____

Ц- = [k^NP-lk^lAtW + С) - [кМАР at ~~ — —

Ц- = - ft.Wlf + С) - [к5Л]ВР

at — — —

at — —

N = N0-(A + B + D)

(f\fr + Cy' для SO, Р = • W''(W + С)'1 для HSO, W-'iW + Cy1 для S2Ot

(7)

И

Таблица 1 - Рассчитанные интервальные значения логарифмов предэкспоненциальных

Параметр Реагент

SO3 HS03+ s2o6

Inkoi [16.94,17.011 Г13.59, 13.651 [12.94, 12.99]

In ko2 [13.98,14.041 [13.81, 13.87] [12.91, 12.96]

In коз Г28.26,28.371 [22.16,22.25] [20.20, 20.281

In ko4 [19.49, 19.561 [23.17,23.271 [25.77, 25.871

In k05 [15.09,15.151 [12.37,12.421 [12.07, 12.121

In к0б [15.26, 15.321 [12.03, 12.08] [11.18, 11.231

Ei, ккал/молъ [11088, 111331 [6893, 6921] [5864, 5887]

E-2, ккал/молъ [13136, 131891 [13137, 131891 [13137, 13189]

E3, ккал/молъ [20679, 207621 [14186,14243] [11389, 11434]

E4, ккал/молъ [20379, 204611 [23776, 23871] [27672, 277831

Es, ккал/молъ [13136, 131891 [8991, 9027] [7592, 76231

Ee, ккал/молъ [11788, 118351 [7592,76231 [6234, 62591

-------------"I----™ w 4..1L* w IjiUI w llu^lDWpyR^k-na IHliUl^a U

том, что сульфирующим агентом является S03.

Аналогичные исследования были проведены для сульфирования нафталина в среде дихлорэтана до a-нафталин сульфокислот. Результаты моделирования и адекватность полученных результатов реальному процессу представлены в диссертации в разделе 3.5.2 и подтверждают работоспособность разработанного программного комплекса.

3.2 Исследование математической модели алкилирования фенилацетонит-рила в микроструктурном реакторе на основе модели идеального вытеснения с использованием интервального анализа

Прямая задача

Реакция алкилирования фенилацетонитрила при условии раздела фаз проходит на границе контакта водной и органической фаз. Применение микрореактора обеспечивает эффективное перемешивание фаз и смещение реагентов к поверхности раздела фаз, что позволяет получить максимальный выхода целевого продукта.

Для решения прямой задачи химической кинетики процесса алкилирования с помощью разработанного комплекса программ Interval Kinetic сначала необходимо ввести исходные данные: схему реакций, начальные приближения для констант скоростей. На рисунках 2 и 3 представлены формы ввода начальных данных для процесса алкилирования:

—

CI(QH7N)= |3.013 I * 13.013 |

C2(NaOH)= |3.31 1 + 13-31 |

Cj(CeH6N)= 1° 1 + |0 I

Q(Na)= 1° + |0 1

Cs(H20)= |o J + |0 1

Cs(CjH5B0= I + 13.24 |

C7(CioHuN)= 1° 1 + lo "I

ce(Bf)= |0 1 + |o I

C9(CioH10N)= |o 1 + lo I

CwiCuHuN)- |o 1 - lo I

веществ в реакции алкилирования

2. 1

3. 1

]свН6м + [Г ]свН6м + ГГ

|№ + [1 1 С2Н5Вг

|н2о.

№ЮН

кЗ

СюНцМ + ^

Вг

4. 1

СюНцМ+у.

lBrJil.iT-

СвНбМ+|1

С2Н5Вг

5. 1

|с10н11м + 11

lNaOH-ii.IT

I С10Н10М + 11

N3+1

6- 1

|Са0Н10Ы+|1

N3+1

|н2о-

7. |1 | С10Н10М +11 | С2Н5Вг -

к а

-»11_|с10ним + [1_

| С12Н15М + ¡1 1 Вг

Jн2o

ЫаОН

8.1

I С12н15ы +¡1

Вг-

I с10н10м + II

С2Н5Вг

Рисунок 3 - Матрица стехеометрических коэффициентов реакции при решении прямой задачи кинетики для реакции алкилирования

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

Сг'М = -Оч ЙИ с2М) + с2'М = -(к1саис2М) + с3'м = к1с1[4с2м-к2

с* М = ка С1М с2м - к2 сз'М - ^ С! и с2и - к2

Се'М = -<к3 СзИ СеМ) +

с7'и = к3 с3и СеИ - к5 Св'и = к3 с3м СеМ - кз Сэ'[1] = к5 С2[1] С7[1] + кв Сю'И = -(И-з СюМ СвМ)

к2с3мс4мс5м

кг с3М С,и с5и - к5 с2р] с7М -1% С4М С5И Сэр] СзМ С4Й С5М - к3 СзМ СеМ + к, с7М СвИ СзМ с,И с5м ■+ к5 с2м с7м - ^ С4И с5м сэи сзМ С4И С5Ю + к5 сги С7М - ко С4М С5М сэи кв СаоМ СвМ + с7М СвИ - к7 СеМ Сэ[1] ОгМ с7[1] - к, С7[1] СзМ + кб С$р] С5[1] Сэр] ОюМ СвМ - к4 С7И СеМ + к7 СеМ СЭМ

Сю[Ч Св[Ч - кб С,[1] С5[1] Сэ[1] - к7 Сб[1] С9[1] + к7 СеМ сэи

Рисунок 4 - Система дифференциальных уравнений для реакции алкилирования Обратная задача

Для решения обратной задачи в программе необходимо задать эксперименталь-

Время С1 СЗ С9 СЮ

1 1.371 0.171 0.012 0.021

2 1.077 0.112 0.009 0.056

3 0.835 0.072 0.01 0.06

4 0.677 0.048 0.007 0.1

Рисунок 5 - Исходные значения концентраций для решения обратной задачи алкилирования фенилацетонитрила

В результате решения обратной задачи были получены интервальные значения констант скоростей реакции.

В таблице 2 представлены экспериментальные и расчетные значения концентрации при известных экспериментальных значениях концентраций для четырех веществ.

Таблица 2 - Результаты решения обратной задачи химической кинетики для процесса алкилирования фенилацетонитрила с учетом относительной погрешности

экспериментальных данных

Обозначения веществ Моменты времени, мин

1 2 3 4 1 2 3 4

Экспериментальные значения концентраций, кмоль/м3 Расчетные значения концентраций (верхняя и нижняя границы), кмоль/м3

с, 1.371 1.077 0.835 0.677 [1.362, 1.479] [0.975, 1.096] [0.793, 0.911] [0.682, 0.799]

с2 - - - - [1.624, 1.745] [1.210, 1.332] [1.010, 1.128] [0.888, 0.999]

С3 0.171 0.112 0.072 0.048 [0.172, 0.208] [0.094, 0.111] [0.063, 0.075] [0.048, 0.056]

с4 - - - - [1.565, 1.686] [1.978, 2.100] [2.182, 2.300] [2.310, 2.422]

С5 - - - - [1.565, 1.686] [1.978, 2.100] [2.182, 2.300] [2.310, 2.422]

Сб - - - - [1.748, 1.881] [1.248, 1.378] [1.015, 1.137] [0.875, 0.991]

С7 - - - - [1.302, 1.437] [1.741, 1.886] [1.937, 2.087] [2.049, 2.203]

С8 - - - - [1-359, 1.492] [1.862, 1.992] [2.103, 2.225] [2.249, 2.365]

с9 0.012 0.008 0.010 0.007 [0.009, 0.012] [0.008, 0.012] [0.007, 0.010] [0.006, 0.009]

Сю 0.021 0.056 0.060 0.100 [0.020, 0.026] [0.044, 0.060] [0.059, 0.084] [0.071, 0.103]

4 Особенности использования алгоритмов и программ для математического моделирования биохимических процессов в микрореакторах

В данной главе диссертации рассматривается интервальная методика математического моделирования для процесса биокаталитического восстановления этилового эфира циклогексанонкарбоновой кислоты в этиловый эфир (Ш,28)-цис-2-гидроксициклогексан карбоновой кислоты в присутствии 5асскаготусм сгге\та проводимого в микрореакторе и реакторе смешения.

Предложено три математические модели описания процесса в микроструктурном реакторе. В автореферате приводится математическая модель, описывающая кинетику для одной из гипотез протекания биохимического процесса. Результаты для двух других моделей приведены в диссертации. Полученные результаты позволяют судить о химизме биохимического процесса в микроструктурном реакторе.

dr — '

где X, S, Pl, P2 ,N— соответственно концентрации биомассы (фермент), субстрата, продуктов и этанола (моль/л);

Мтах - максимальная удельная скорость роста микроорганизмов, (ч'1);

K¡ - коэффициент полунасыщения (константа Михаэлиса), константа, равная концентрации субстрата, при которой скорость роста равна У2 максимальной (концентрация полунасыщения), (моль/л);

Yx/s - константа выхода биомассы для субстрата, (моль/моль);

Ypi,s - константа выхода продукта Р1 для субстрата, (моль/моль);

УР2/.s - константа выхода продукта Р2 для субстрата, (моль/моль);

m¡ - константа поддержания субстрата, (моль/моль/ч);

К*' - константа образования этанола;

kd - константа «отмирания», (ч'1).

На рисунке 6 представлены результаты расчетов по трем моделям для одной из выборок, где сплошной линией обозначаются кинетические кривые, а крестиками обозначаются экспериментальные данные.

Модель Моно 1.3 Модель Мово 2

Рисунок 6 - Интервальные зависимости концентраций субстрата А, продуктов Р1, Р2 и этанола N от времени для моделей 1, 2, 3 (выборка 1)

Интервальные значения констант для разных моделей отличаются незначительно. Выборки 1 и 3 описываются близкими наборами констант. В выборке 2 константа образования этанола KN почти в два раза меньше, чем в выборках 1 и 3, что может быть связано с расположением экспериментальных точек выборки в начале временного интервала. Все модели Моно удовлетворительно описывают экспериментальные данные для реактора смешения и микрореактора. При близких значениях /w, константы Ks, Кл, Ух/s, Yx/a,Ypi/a, YP2/a, kd, KN одинакового порядка. Для обоих субстратов скорости их поддержания стремятся к нулю.

5 Алгоритмы и программы для исследования чувствительности кинетики

химических реакций

Одна из задач исследования химических реакций - это оценивание влияния значений констант скоростей на ход протекания реакции. По рекомендациям ИЮПАК стадию, определяющую скорость, устанавливают вычислением контролирующего фактора (CF). Для реакции i, константа скорости которой kh CF, при постоянстве остальных k¡ и

констант равновесия KJt равен: CF= . Стадия, с существенно большим по

сравнению с другими CF, является определяющей.

В работе рассматриваются два возможных подхода для оценки влияния констант скоростей реакции к изменению концентраций веществ. Первый подход - традиционный анализ чувствительности с помощью дополнительной системы дифференциальных уравнений, отражающих чувствительность. Второй предлагаемый подход - оригинальный интервальный метод оценивания влияния концентраций к изменению констант скоростей химической реакции в различные моменты времени.

5.1 Постановка задачи

Модель химической кинетики записывается в виде системы дифференциальных уравнений первого порядка: dC.

-jp = /;(*,.*„...,*.;С1РС,...,С,) U = 1,2,...,и)

с заданными начальными условиями СЛ0 = С% где С) - концентрации исходных веществ, kq - константы скоростей отдельных стадий реакции, q = ¡,...,m. Решением этой системы дифференциальных уравнений, для заданных начальных условий, является набор концентраций q от времени, которые зависят от параметров - констант скоростей Kqt

5.2 Оценивание чувствительности значений констант скоростей Дифференциальный метод анализа чувствительности

Рассмотрим чувствительность на примере набора элементарных реакций газофазного окисления метана при температуре Т=1000К:

СЯ4 + 02 -й-> CH¡ + НОг НСО + 02 СО + Н02

СН3 + 02 —Ctfз оо сн4+но2-^->н2о2+снг

СН3 оо > снг о+он сн2о + яо2 —я202 + НСО

CHi+0H-^H20 + 0H со + он ш >со2 + я

сн2о+он >Я20+СЯ3 сн4 + н- »■ >ся3 +нг

CH20+02 -1Í-+ нсо+но2 сн2о + н~^нсо + н

По исходной схеме реакций по предложенному алгоритму программой автоматически генерируется система дифференциальных уравнений, характеризующая

изменение концентраций веществ. В результате дифференцирования \ исход-

ам, ^ а ;

ной системы уравнений по кд получаем т - систем дифференциальных уравнений, каждая из которых содержит и - уравнений (всего получается т • п уравнений). Полученные системы уравнений следует решать совместно с исходной системой при соответствующих начальных условиях. При этом получаются не только значения концентраций Су , но и т ■ п функций для каждого момента времени, которые характеризуют оценку чувствительности концентраций С, к изменению различных констант скоростей реакции кд. Для оценки чувствительности воспользуемся нормированным критерием: к,дС^,к) ЗС,(/,*) , , „

ок_

Этот критерий чувствительности позволяет проводить оценку влияния кя на концентрацию Су в момент времени I. Таким образом, становится возможным определить степень влияния констант на различные концентрации, то есть выяснить какие из констант скоростей кч и на каких стадиях реакции являются определяющими, а какие являются мало влияющими.

Для оценки степени влияния на концентрации веществ необходимо к этой исходной системе дифференциальных уравнений добавить систему уравнений, полученную дифференцированием исходной системы уравнений по к,.

Для определения чувствительности относительно т-констант скоростей реакции

8 (дС,

дк„

81

. Также

необходимо решить систему из п(т+1) уравнений для функций Су и

возможно совместно решать только по две системы: исходную и полученную дифференцированием по одной константе, после чего - для следующей константы и т.д. для всех кч.

На рисунке 7 представлена логарифмическая зависимость критериев чувствительности от времени по концентрации формальдегида для различных констант скоростей кд (д = 4, 6, 8). Анализ графиков позволяет увидеть, что значения критериев чувствительности могут быть положительными, отрицательными или нулевыми. В ходе реакции значение критерия чувствительности может существенно меняться, как это видно из рисунка 7.

-о.оою

/

\

0.0011 ; \ \ / д\ > 1 ' \ __ гс^/гьйц

0.0010

0.0001 __ь ; \ \ // \ -- < 1 , гс^/йп!^

0.0091

0.0005 1

Рисунок 7 - Зависимость чувствительности концентрации формальдегида (С6) от времени (Г, с) к изменениям кч (д = 4, б, 8)

Как следует из полученных результатов, концентрация формальдегида наиболее чувствительна к константам к4, к6 и кв в момент времени /я = 310 секунды, а также она чувствительна к константам к;нк9ъ тот же момент времени. Значения критериев чувствительности по к4, к6, к8 положительны в момент времени Гя, а значения критериев по к5 и ко отрицательны. Чувствительность концентрации формальдегида к изменению к,, къ кз, к7, кц относительно мала вблизи 1т, т.к. находится в области нуля из-за того, что происходит изменение знака критерия чувствительности.

Предлагаемый интервальный метод анализа чувствительности

Для определения чувствительности концентраций к изменению констант скоростей реакции нами предлагается использовать другой подход, суть которого заключается в применении интервальных вычислений при решении системы исходных дифференциальных уравнений.

Система дифференциальных уравнений решается для констант скоростей реакции, которые представляют собой интервальные оценки с разбросом относительно среднего значения.

В результате решения системы дифференциальных уравнений получаются функции, которые характеризуются не одним значением концентрации в каждый момент времени, а интервалом. Этот интервал имеет нижнюю и верхнюю границу значений концентрации. Ширина полученных интервалов концентраций характеризует степень влияния конкретной константы скорости реакции на конкретную концентрацию. Чем шире рассчитанный интервал, тем большее влияние оказывает константа на эту концентрацию. Чем уже интервал, тем меньшее влияние оказывает данная константа.

К сожалению, при этом подходе невозможно определить положительное или отрицательное влияние оказывает данная константа на концентрацию. Определяется только степень влияния. В этом состоит недостаток этого метода.

Ширина интервала на протяжении времени меняется. В качестве критерия оценки влияния константы скорости реакции кч (критерия чувствительности) на концентрацию Су предлагается вычислить значение площади между кривой верхней и нижней оценки концентрации для всех моментов времени. Площадь будет тем больше, чем больше ширина интервалов между верхней и нижней оценкой концентрации.

Несмотря на то, что числовые значения чувствительности для дифференциального метода различаются на несколько порядков от значений чувствительности для интервального метода, соотношения чувствительностей для различных веществ для этих методов совпадают, что говорит о тождественности этих подходов.

Поскольку значения констант скоростей различаются, то различаются и их интервальные значения. Очевидно, что если константа имеет самую большую интервальную оценку, то ее варьирование оказывает самое незначительное влияние на процесс в сравнении с влиянием других констант, которые имеют меньшие диапазоны интервальных оценок. Сделаны выводы о чувствительности концентраций веществ к изменениям констант скоростей исследуемых химических реакций.

Считаю своим долгом поблагодарить директора института технической химии г. Дрездена профессора Решетиловского В.П., доцента, к.т.н. Чепикову ВН и к.т.н. Боровинскую Е.С. за оказанные консультации по моделированию процессов в микрореакторах.

Выводы

1. На основе интервального анализа проведено комплексное исследование кинетики важнейших химических и биохимических процессов в химических реакторах и микро-

18

структурных реакторах: сульфирование нафталина, окисления метана, алкилирование фенилацетонитрила, биокаталитическое восстановление этилового эфира циклогекса-нонкарбоновой кислоты в этиловый эфир (111,25)-цис-2-гидроксициклогексан карбоно-вой кислоты в присутствии Saccharomyces cerevisia.

2. Метод огибающих в совокупности с методом глобальной оптимизации случайного поиска позволяет в разработанном комплексе программ Interval Kinetic эффективно реа-[изовывать численные методы решения прямой и обратной задач химической кинетики (ля различных процессов.

i. Предложенный способ ввода исходных данных позволяет значительно упростить адание схемы реакции и генерирования системы дифференциальных уравнений. ■. В диссертации с использованием предложенного комплекса программ впервые бы-ги получены интервальные оценки кинетических параметров для важнейших процессов :имической технологии.

!. Предлагаемый интервальный метод оценки чувствительности позволяет оценивать тепень влияния значения констант скоростей химической реакции на концентрации ¡еществ.

>. Предложенная впервые интервальная модель кинетики сульфирования нафталина 'бедительно подтвердила представление о том, что сульфирующим агентом при сульфировании является диоксид серы.

Впервые составлена интервальная модель кинетики процесса для реакции алкилиро-;ания фенилацетонитрила, которая адекватно описывает качественные и количествен-[ые свойства этого процесса в микроструктурном реакторе.

!. Впервые составлена интервальная математическая модель для описания биокатали-ического восстановления этилового эфира циклогексанонкарбоновой кислоты в этило-1ЫЙ эфир (111,28)-цис-2-гидроксициклогексан карбоновой кислоты в присутствии iaccharomyces cerevisia в микрореакторе, которая позволяет подойти к проблеме интен-ификации исследуемого биохимического процесса.

». Предложены комплексы программ, реализованные в виде модулей, зарегистриро-¡анных в объединенном фонде электронных ресурсов "Наука и образование" (ОФЭР-1иО).

Список основных публикаций по теме диссертации

. Хайдаров, А.Г. Использование системы Wolfram Mathematica для оценивания чув-твительности значений констант скоростей в системах уравнений химической кинети-м на примере окисления метана / А.Г. Хайдаров, В.А. Холодное, Е.С. Боровинская, Ш. Решетиловский // Известия Смоленского государственного университета -2010 М,-С. 99-109.

!. Хайдаров, А.Г. Интервальный метод оценки чувствительности констант скоростей фи решении системы уравнений химической кинетики окисления метана / А.Г. Хайда-юв, В.А. Холодное, Е.С. Боровинская, В.П. Решетиловский // Информационные систе-ш и технологии. -2010, № 6. -С. 61-65.

i. Хайдаров, А.Г. Программный комплекс для исследования кинетики химических реакций с использованием интервальных методов / А.Г. Хайдаров, В.А. Холоднов, Е.С. ¡оровинская, В.П. Решетиловский // Материалы конф. Тринадцатая межвузовская учеб-ю-методическая конференция «Актуальные проблемы химико-технологического обра-ования». -М: РХТУ им. Д.И. Менделеева, -2011. - С. 83- 84.

г. Хайдаров, А.Г. Исследование кинетики химических реакций с использованием ин-ервальных методов / А.Г. Хайдаров, В.А. Холоднов В.А., Е.С. Боровинская, В.П. Реше-

тиловский // Материалы XII Международная конференция «Системы компьютерной математики и их приложения», Смоленск: СмолГУ, -2011. Вып. 12. - С. 67-69.

5. Хайдаров, А.Г. Вычислительный эксперимент для сложной системы химических реакций с использованием Mathcad / В.А. Холоднов, Е.С. Боровинская., A.B. Гайков, А.Г. Хайдаров // Материалы XI Международная конференция «Системы компьютерной математики и их приложения», Смоленск: СмолГУ, -2010. Вып. 11,- С. 74-76.

6. Хайдаров, А.Г. «Исследование чувствительности кинетики химических реакций с помощью интервальной математики» / А.Г. Хайдаров, В.А. Холоднов В.А., Е.С. Боровинская, В.П. Решетиловский // Сб. трудов XXIV Международная научная конференция «Математические методы в технике и технологиях ММТТ-24» Т.З. Киев: г. Саратов «ИППОЛиТ-XXI век» -2011.- С. 125-127.

7. Хайдаров, А.Г. Задачи применения интервального метода в системном анализе химических технологий// Материалы VI Всероссийской межвузовской конференции молодых ученых. Сб. трудов конференции молодых ученых, Выпуск 7. Тезисы докладов научной школы «Технологии высокопроизводительных вычислений и компьютерного моделирования» . - СПб: СПбГУ ИТМО. -2009. - С. 36.

8. Хайдаров, А.Г. Свидетельство о регистрации электронного ресурса «Интервальная кинетика». //№17896, зарегистр. в ОФЭРНиО 10.02.2012 г.

9. Хайдаров, А.Г. Свидетельство о регистрации электронного ресурса «Анализ чувствительности при исследовании кинетики химических реакций». // № 17895, зарегистр. в ОФЭРНиО 10.02.2012 г.

Отпечатано с оригинал-макета. Формат 60x90'/i6 Печ.л. 1,25 .Тираж экз. 90. Зак. Ks 54

Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Санкт-Петербургский государственный технологический институт (технический университет)»

190013, Санкт-Петербург, Московский пр., 26 Типография издательства СПбГТИ(ТУ), тел. 49-49-365, e-mail: publ@technolog.edu.]

61 12-5/3539

Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Санкт-Петербургский государственный технологический институт (технический университет)»

МЕТОДЫ, АЛГОРИТМЫ И ПРОГРАММЫ МОДЕЛИРОВАНИЯ КИНЕТИКИ ХИМИЧЕСКИХ И БИОХИМИЧЕСКИХ ПРОЦЕССОВ С ИСПОЛЬЗОВАНИЕМ

ИНТЕРВАЛЬНОГО АНАЛИЗА

05.13.18 - Математическое моделирование, численные методы

и комплексы программ

На правах рукописи

ХАЙДАРОВ АНДРЕЙ ГЕННАДЬЕВИЧ

ДИССЕРТАЦИЯ НА СОИСКАНИЕ УЧЕНОЙ СТЕПЕНИ КАНДИДАТА ТЕХНИЧЕСКИХ НАУК

Научный руководитель -

доктор технических наук, профессор Холоднов В. А.

Санкт-Петербург

2012 г

Основное содержание работы

Введение.........................................................................................................................................4

1 Литературный обзор................................................................................................................11

1.1 Основные положения моделирования процессов в химических реакторах................11

1.2 Использование интервального анализа при моделировании химико-технологических процессов и систем..................................................................................................................13

1.3 Применение методов математического моделирования и численных методов для решения обратных задач.........................................................................................................15

1.3.1 Постановка обратной кинетической задачи......................... ....................................15

1.3.2 Построение функционалов и методы их минимизации при решении обратной кинетической задачи..........................................................................................................16

1.3.3 Особенности решения обратных задач....................................................................17

1.3.4 Задача нахождения кинетических параметров сложного процесса......................18

1.4 Глобальная оптимизация..................................................................................................21

1.5 Методология интервального анализа и ее применение в химической кинетике........23

1.6 Применение методов оптимизации для решения задач по построению кинетических моделей.....................................................................................................................................25

1.7 Обзор существующих комплексов программ.................................................................26

1.7.1 Анализ прикладных программных комплексов для решения задач химической кинетики...............................................................................................................................26

1.7.2 Программные системы, разработанные для использования через Интернет.......28

Выводы.....................................................................................................................................30

2 Алгоритмы математического моделирования численных методов и комплексов программ для интервальных методов.......................................................................................31

2.1 Интервальная арифметика................................................................................................31

2.2 Интервальные методы решения систем обыкновенных дифференциальных уравнений.................................................................................................................................33

2.2.1 Теория интервального решения обыкновенных дифференциальных уравнений 33

2.2.2 Методы решения систем дифференциальных уравнений в интервальной постановке............................................................................................................................38

Выводы.....................................................................................................................................55

3 Использование комплекса программ для решения практических задач химической кинетики с помощью интервальных методов...........................................................................57

3.1 Объекты исследования и постановка задачи интервального метода...........................57

3.2 Требования, предъявляемые к комплексу программ.....................................:...............59

3.3 Структура и назначение комплекса программ Interval Kinetic.....................................59

3.4 Разработка и применение комплекса программ Interval Kinetic...................................61

3.4.1 Алгоритм нахождения верхних и нижних границ концентраций.........................62

3.4.2 Алгоритм нахождения интервалов констант скоростей реакции..........................63

3.4.3 Разработка пользовательского интерфейса.............................................................68

3.4.4 Системные требования к программному продукту................................................68

3.5 Решение практических задач химической кинетики.....................................................69

3.5.1 Процесс получения (3-нафталинсульфокислоты.....................................................69

3.5.2 Применение комплекса программ для моделирования кинетики а - нафталин сульфокислоты....................................................................................................................79

3.5.3 Применение комплекса программ для моделирования кинетики алкилирования фенилацетонитрила.............................................................................................................84

Выводы.....................................................................................................................................96

4 Особенности использования алгоритмов и программ для математического моделирования биохимических процессов в реакторе смешения и микрореакторе............97

4.1 Микроструктурные реакторы...........................................................................;...............97

4.2 Использование интервального анализа для исследования биокаталитического процесса в микрореакторе......................................................................................................98

4.2.1 Модель биохимического процесса...........................................................................99

4.2.2 Интервальная неопределенность в микрореакторе...............................................101

4.3 Исследование влияния интервальной неопределенности в реакторе смешения......106

Выводы...................................................................................................................................112

5 Алгоритмы и программы для исследования чувствительности кинетики химических реакций.......................................................................................................................................113

5.1 Постановка задачи...........................................................................................................113

5.2 Оценивание чувствительности значений констант скоростей на примере окисления метана.....................................................................................................................................114

5.2.1 Дифференциальный метод анализа чувствительности.......................................114

5.2.2 Интервальный метод анализа чувствительности..................................................121

5.3 Применение комплекса для оценивания чувствительности значений констант скоростей на примере алкилирования фенилацетонитрила..............................................123

5.4 Интервальный метод анализа чувствительности.........................................................129

Выводы...................................................................................................................................137

Заключение................................................................................................................................139

Литература.................................................................................................................................141

Введение

Актуальность работы

Компьютерное моделирование химических реакторов доказало свою актуальность и перспективность, поскольку позволяет выполнять расчет материальных, тепловых балансов и осуществлять поиск наилучших режимов функционирования.

Решающее значение при выборе условий проведения химико-технологических процессов и химико-технологических систем имеют вопросы скорости химических превращений, изучаемые химической кинетикой. Кинетические уравнения, содержащие необходимую информацию об основных закономерностях химических превращений, являются первоосновой математической модели химического реактора и получаются при решении, так называемых прямых и обратных задач химической кинетики.

В настоящее время накоплен значительный опыт по решению традиционных для химической кинетики прямых и обратных задач. Следует отметить большой вклад в развитие построения кинетических моделей и программных комплексов работ д.т.н., проф. Ю.В. Шарикова, доц. В.И. Коробова, д.т.н., проф. В.Ф. Очкова и др. В отличие от традиционного подхода в работах этих авторов решение рассматривается с помощью современных систем компьютерной математики Mathcad, Maple, а также организации сетевых расчетов на Mathcad Calculation Server. Но в данных работах игнорируется вопрос существования неопределенности кинетических параметров и моделирование процессов, решения обратных задач основаны на номинальных значениях кинетических параметров.

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

подсистема химического превращения вносит существенный вклад в структуру себестоимости получения любого конечного продукта.

В последнее время в научной литературе опубликованы работы, в которых рассмотрены вопросы влияния неопределенности исходной информации на результаты моделирования технических объектов (prof. Grossmann I.E., д.т.н., проф. Островский Г.М. и др.). В работах А.П. Вощинина, С .П. Шарого, Ю.И. Шокина и других при моделировании и оптимизации технических объектов в условиях неопределенности используются интервальные методы анализа.

В настоящей работе впервые при решении прямых и обратных задач для сложных химических и биохимических реакций используется интервальный анализ. В работе речь идет об интервальной неопределенности

кинетических параметров.

В диссертационной работе при решении прямых и обратных задач используются как традиционные, так и интервальные методы с использованием систем компьютерной математики Mathematica.

Предложенные автором методы, алгоритмы и программы для моделирования кинетики химических и биохимических реакций с использованием интервального метода опробованы для исследования

сложных процессов.

Актуальность диссертационного исследования определяется тем, что в нем приведены разработанные методики, алгоритмы и программы для моделирования кинетики химических, биохимических реакций в химических и микроструктурных реакторах в условиях интервальной неопределенности кинетических параметров с использованием интервальных методов. Впервые предложен оригинальный метод определения чувствительности значений концентрации веществ к изменению констант скоростей реакций.

Разработка этих методов, алгоритмов и программ позволит на новом уровне подойти к проблеме моделирования процессов в химических и микроструктурных реакторах.

Цель диссертационной работы. Разработка совокупности взаимодополняющих методов, алгоритмов и программ для определения кинетических параметров сложных химических и биохимических процессов на основе экспериментальных данных с использованием интервальных методов.

Для достижения поставленной цели решены следующие задачи:

• Проведен анализ существующих алгоритмов и программ для сложных химических и биохимических процессов на основе экспериментальных данных с использованием интервальных методов;

• Обоснование выбора и адаптация метода огибающих для решения систем дифференциальных уравнений при решении прямой и обратной задач химической кинетики в интервальной постановке;

• Адаптация глобального метода оптимизации пакета программ МаЛетайса для решения обратных задач в интервальной постановке химической кинетики;

• Создание пользовательского интерфейса для решения прямых и обратных задач кинетики химических и биохимических процессов в интервальной постановке;

• Разработка методов, алгоритмов и программ для решения прямых

0 и

и обратных задач химическои кинетики в интервальной постановке. Тестирование предложенных методов, алгоритмов и программ с использованием вычислительного и натурного экспериментов;

• Методика моделирования процессов на основе разработанных алгоритмов;

• Анализ чувствительности концентраций веществ к изменению

кинетических параметров.

Методы исследования. В диссертации использованы методы: системного анализа, химической кинетики, математического моделирования и оптимизации химико-технологических процессов, компьютерной и интервальной математики.

Инструменты исследования. Интерактивные вычислительные системы Wolfram Mathematica, программирование пользовательского интерфейса.

Обоснованность научных результатов обеспечивается применением современных математических методов моделирования кинетики химических, биохимических реакторов и микрореакторов. Тестированием разработанных алгоритмов и программ на контрольных примерах и совпадение результатов вычислений с результатами натурных экспериментов.

Достоверность теоретических разработок подтверждена совпадением результатов вычислительного эксперимента на ЭВМ с литературными данными, что позволяет сделать вывод об эффективности разработанных алгоритмов и программ с использованием интервальных методов.

Основные положения, выносимые на защиту:

1. Методы и алгоритмы:

а) решения прямой и обратной задачи кинетики химических и биохимических реакций в интервальной постановке;

б) вычисления чувствительности концентраций веществ к изменению кинетических параметров;

2. Алгоритм автоматического генерирования системы дифференциальных уравнений на основе различных подходов к составлению схемы химических реакций;

3. Модификация алгоритма огибающих для решения систем дифференциальных уравнений с учетом погрешности экспериментальных данных;

4. Программные комплексы:

а) для решения прямых, обратных задач кинетики в интервальной постановке;

б) для исследования чувствительности концентраций веществ к изменению констант скоростей реакций;

5. Математические модели в интервальной постановке:

а) для исследования процессов тонкого органического синтеза;

б) для описания кинетики биохимических процессов в микрореакторах;

6. Результаты математического моделирования исследованных процессов.

Научная новизна работы:

Проведено комплексное исследование кинетики важнейших химических и биохимических процессов в химических реакторах и микроструктурных реакторах в условиях интервальной неопределенности кинетических параметров с применением современной технологии математического моделирования, натурного и вычислительного эксперимента:

• Применение интервального метода для определения кинетических параметров химических и биохимических процессов в микрореакторах по экспериментальным данным;

• Модели исследованных химико-технологических процессов:

а) тонкого органического синтеза: получение бета нафталин сульфокислоты, альфа нафталин сульфокислоты, алкилирования фенилацетонитрила в микрореакторе в интервальной постановке

б) биохимического биокаталитического восстановления этилового эфира циклогексанонкарбоновой кислоты в этиловый эфир (111,28)-цис-2-гидроксициклогексанкарбоновой кислоты в присутствии БасскаготусеБ сегеугягае в микроструктурном реакторе в интервальной постановке.

• Способ и алгоритм автоматизированного задания схемы химических реакций;

• Модификация метода анализа чувствительности параметров математических моделей на основе интервального метода;

• Модификация интервального метода огибающих для решения систем дифференциальных уравнений описывающих химическую кинетику.

Практическая значимость и реализация результатов работы. На

основе теоретических результатов работы предложены и разработаны алгоритмы и программы для решения моделирования химических, биохимических реакций в реакторах и микрореакторах. Для иллюстрации работоспособности предлагаемых методов и алгоритмов решены задачи по комплексному исследованию различных химических процессов в реакторах.

Разработанные методы, алгоритмы и программы используются в учебном процессе в Санкт-Петербургском государственном технологическом институте, Казанском национальном исследовательском технологическом университете, Березниковском филиале Пермского государственного технического университета и в Тамбовском государственном техническом университете.

Работа была выполнена в рамках следующих проектов при поддержке Правительства РФ:

• Проект «Фундаментальные исследования закономерностей взаимосвязи основных параметров каталитических процессов в микроструктурных реакторах в условиях неопределённости исходной информации», реализуемый в рамках Аналитической ведомственной целевой программы "Развитие научного потенциала высшей школы" 2009-2010 гг.

• Проект «Интервальный подход к решению задач по математическому моделированию и оптимизации химических и биохимических процессов в микроструктурных реакторах в условиях неопределенности исходной информации», реализуемый в рамках Аналитической ведомственной целевой программы "Развитие научного потенциала высшей школы" 2009-2010 гг.

• Проект «Математические модели и методы интервального анализа для исследования и оптимизации биокаталитических процессов в микрореакторе», реализуемый в рамках федеральной целевой программы «Научные и научно-педагогические кадры инновационной России» на 20092013 гг.

Апробация работы.

Материалы диссертации докладывались и обсуждались на Всероссийских и международных научных конференциях: «Системы компьютерной математики и их приложения», Смоленск, СмолГУ, 2011 г.; «Математические методы в технике и техн�