CC BY
15
УПРАВЛЕНИЕ, МОДЕЛИРОВАНИЕ, АВТОМАТИЗАЦИЯ
DOI: 10.24143/2072-9502-2017-2-125-135 УДК 574.6.663.1
Ю. Л. Гордеева, Ю. А. Комиссаров, Е. Л. Гордеева, А. Г. Бородкин
АЛГОРИТМ ИСПОЛЬЗОВАНИЯ УСЛОВИЙ ГУРВИЦА ДЛЯ ИССЛЕДОВАНИЯ НА УСТОЙЧИВОСТЬ СТАЦИОНАРНЫХ СОСТОЯНИЙ НЕПРЕРЫВНЫХ ПРОЦЕССОВ ФЕРМЕНТАЦИИ С НЕЛИНЕЙНОЙ КИНЕТИКОЙ
Приведен алгоритм оценки устойчивости стационарных состояний процесса микробиологического синтеза на основе математических моделей кинетики, базирующихся на неструктурированном подходе. Объект исследования - непрерывный биотехнологический процесс, в ходе которого помимо биомассы получается целевой продукт. Оценка устойчивости выполнена с использованием матрицы Гурвица. В процессе возможно использование сырья, содержащего компоненты для образования дополнительного количества основного субстрата. Алгоритм включает пять этапов: формирование математической модели, учитывающей биологические и технологические ограничения; определение стационарных состояний (одного или нескольких), для которых необходимо оценить условия устойчивости и осуществить расчет показателей процесса для принятых стационарных состояний в результате решения системы нелинейных алгебраических уравнений математической модели; запись системы уравнений первого приближения (система линейных дифференциальных уравнений первого порядка с постоянными коэффициентами) и оценка коэффициентов, отвечающих принятым стационарным состояниям; формирование матрицы Гурвица и вычисление ее элементов; расчет показателей необходимого и достаточного условия устойчивости на основании матрицы Гурвица. Приведены два варианта вычисления элементов матрицы Гурвица: преобразование системы четырех дифференциальных уравнений первого приближения в одно дифференциальное уравнение четвертого порядка, в ходе которого определяются расчетные соотношения для коэффициентов матрицы, и вычисление коэффициентов в результате решения векторного уравнения по системе первого приближения с введением собственных чисел. Проведена оценка устойчивости двух стационарных состояний (одно из которых оптимальное) процесса микробиологического синтеза получения молочной кислоты из сырья, содержащего компонент, воспроизводящий основной субстрат. Для численных значений констант (из литературных источников) оба стационарных состояния оказались устойчивыми.
Ключевые слова: микробиологический синтез, оценка устойчивости, матрица Гурвица.
Введение
Процесс непрерывного микробиологического синтеза технологически реализуется в стационарных условиях, т. е. показатели процесса не изменяются во времени (концентрация биомассы, субстрата, продукта и т. п.).
При реализации синтеза в процессе всегда имеют место возмущения (отклонения показателей процесса от стационарного значения). Если эти отклонения не приводят к нарушению технологического процесса и с течением времени значения показателей возвращаются к первоначальным, стационарное состояние полагается устойчивым. В противном случае для реализации процесса требуется внешнее управление.
Так как малых возмущений практически избежать не удается, то, если при малых возмущениях наблюдается нарушение технологического процесса, приводящее к невозможности обеспечить его функциональное назначение, процесс считается неустойчивым в малом.
Математический анализ такого рода устойчивости или неустойчивости в малом получил название «устойчивость первого приближения», оценка которой возможна по условиям Гурви-ца. Подробный математический анализ в доступной форме приведен в работах [1, 2].
Таким образом, настоящее сообщение будет базироваться на анализе устойчивости первого приближения. Ниже этот подход будет описан подробнее.
Объект исследования
Анализ устойчивости стационарных состояний технологического процесса возможен с использованием математической модели, воспроизводящей с достаточной точностью поведение показателей процесса.
Процессы микробиологического синтеза распространены достаточно широко. Однако здесь речь идет о промышленно значимых процессах получения целевых продуктов [3], в частности о процессах получения пищевых кислот [4].
Рассматриваемые процессы отвечают следующим условиям: в процессе синтеза производится биомасса X, получается целевой продукт P, расходуется основной субстрат S.
Сущность процесса определяется кинетикой роста биомассы, основным показателем которой является удельная скорость роста ц, т. к. биомасса является продуцентом образования продукта. Соотношения для удельной скорости роста сформированы на основе неструктурированного подхода.
Математически удельная скорость роста для рассматриваемого объекта является в общем случае нелинейной функцией X, S, P. Кроме того, в процессе ферментации удельная скорость роста в той или иной степени может быть ингибирована каждым из показателей X, S, P. В этом случае возможность ингибировании включается в математическое соотношение для ц.
Обращаясь к обзорам [5, 6], отметим, что для отдельных штаммов микроорганизмов при ферментации образуются, кроме целевого продукта, побочные, иногда представляющие практическую ценность. При этом таких продуктов образуется, как правило, незначительное количество [7].
С другой стороны, т. к. сырье является наиболее расходной статьей процесса, имеется тенденция [5, 6] использовать воспроизводимое сырье. Это отмечено в работах [8-10].
Особенность использования воспроизводимого сырья проявляется в том, что в этом сырье могут быть компоненты, которые в процессе организации технологии дают дополнительное количество основного субстрата S. Так, например, пшеничная мука (whole-wheat flour, wheat flour) в качестве сырья для ферментации содержит мальтозу, которая в процессе генерирует дополнительно основной субстрат - глюкозу, что повышает уровень эффективности использования сырьевых ресурсов. Схема такого процесса приведена на рисунке [10].
Схема процесса ферментации с получением дополнительного количества основного субстрата
Учет этой составляющей, а также составляющей побочного продукта позволяет записать обобщенную математическую модель рассматриваемого объекта в виде системы дифференциальных уравнений, представляющих материальный баланс по четырем компонентам - X, Р, М, где М есть концентрация побочного продукта или концентрация компонента, использование которого дает дополнительное количество основного субстрата.
Система уравнений записывается для нестационарных условий процесса, организованного по непрерывной схеме:
dS = F ( s, X, p, m) , (1)
dX = F2 (S, X, P, M) , (2)
dP
— = F3 (S, X, P, M) , (3)
dM = F4 (S, X, P, M) , (4)
где F1, F2, F3, F4 - нелинейные функции показателей процесса (S, X, P, M); S, X, P - концентрации субстрата, биомассы, продукта соответственно; M - концентрация побочного продукта или компонента сырья, деградирующего в основной субстрат (в последнем случае M в соотношениях (2) и (3) отсутствует, соотношение (4) значительно упрощается (будет показано ниже)).
Таким образом, объектом нашего исследования является процесс микробиологического синтеза, оформленный по непрерывному способу. Цель процесса - получение продукта в результате жизнедеятельности микроорганизмов. Математическая модель объекта представляет собой четыре уравнения материального баланса по компонентам: основному субстрату, биомассе, продукту, побочному продукту или компоненту сырья, восполняющему количество основного субстрата. Если в сырье отсутствует побочный продукт или компонент, дополнительно воспроизводящий основной субстрат, уравнения материального баланса записываются только по трем компонентам: S, X, P.
С учетом вышеизложенного система уравнений математической модели в общем виде представлена уравнениями (1)-(4).
Наиболее полно этому описанию объекта соответствует процесс получения молочной кислоты по непрерывному способу из пшеничной муки (wheat flour) [8].
Обращаясь к понятию «устойчивость первого приближения», запишем математические соотношения, определяющие это понятие.
Система уравнений математической модели (1)-(4) преобразуется в систему линейных дифференциальных уравнений первого порядка с постоянными коэффициентами, записанных для малых отклонений от стационарных значений.
Значения переменных представляются в следующем виде: S = Sct + 81; X=Хст + 52; P = Рст + 53; M = Мст + 84, где Sct, Хст, Рст, Мст - значения показателей для конкретного стационарного состояния; 81, 82, 83, 84 - малые отклонения от стационарных состояний.
Нелинейные функции Fi в правой части системы (1)-(4) представлены рядом Тейлора по малым отклонениям с сохранением только двух первых членов ряда (в этом смысл первого приближения). Так, для функции F1 (1) имеем
Л = Мст 62 + (^) 63 + №) 64,
' 1ст lasJ„ 1 Ux)„ 2 Up)„ 3 (дм)„ "
где частные производные вычисляются для конкретного стационарного состояния, т. е. они являются константами. В наших условиях F1ct = 0.
В итоге, обозначив:
Ч§).; Ч§ Ь 1=(f).; Ч!
и использовав эти обозначения для остальных уравнений системы (1)-(4), получим систему уравнений первого приближения:
= öjSJ + bx52 + q53 + dj54, (5)
dt
db
2 = Й25х + b2Ö2 + C2Ö3 + d2Ö4 , (6)
dt
= a3bl + b3b2 + C3S3 + d3b4, (7)
dt db
4 = a4b1 + b4b2 + c453 + d454. (8)
dt
В соответствии с [1, 2], если какое-либо решение системы устойчиво, устойчива и вся система.
Таким образом, система (5)-(8) служит основой для получения заключения об устойчивости первого приближения для любого рассматриваемого стационарного состояния.
Цель настоящей работы - получить алгоритм, выполнение которого позволит оценить условия устойчивости первого приближения для объекта микробиологического синтеза, описанного выше.
Описание алгоритма
Алгоритм включает пять последовательно выполняемых этапов, что обеспечивает решение поставленной задачи, т. е. получение возможности оценить устойчивость первого приближения принятого стационарного состояния.
Этап 1 является определяющим в формировании условий для математического обеспечения решения задачи.
Записываются система уравнений математической модели (1)-(4) и соотношение для удельной скорости роста ц с учетом возможного ингибирования по разным показателям процесса (уравнения получают с использованием источников [5, 6]).
Записываются численные значения констант.
Определяются пределы изменения входных показателей процесса «тт - «тах; Мтп - Мтах (по возможностям сырьевого обеспечения).
Определяются значения величины протока D : Dmln и Dmax.
Значение Dmax вычисляется по условию Dmax < Dпред, где Dпред есть значение протока, при котором исходный субстрат концентрацией «0 полностью вымывается из аппарата.
Так как по материальному балансу по X имеем уравнение для непрерывного процесса в виде D = ц, значение Dпред будет равно ц, в котором « = «0; X = 0; Р = 0; М = 0.
В результате имеется система уравнений математической модели, ограничения по входным показателям процесса, численные значения констант.
На этапе 2 определяются стационарные состояния (одно или несколько), для которых необходимо оценить условия устойчивости.
Стационарные состояния определяются заданием численных значений входных показателей £0, М0 и D, естественно, в заданных пределах.
Для принятых стационарных состояний вычисляются показатели процесса «ст, Хст, Рст, Мст по решению системы линейных алгебраических уравнений вида
^=0 с
F2 = 0 F3 = 0 '' ^ = 0,
На этапе 3 формируется система линейных дифференциальных уравнений первого порядка с постоянными коэффициентами на основе уравнений математической модели - система первого приближения.
Оцениваются значения коэффициентов системы для каждого принятого стационарного состояния.
Система уравнений первого приближения для рассматриваемого объекта имеет вид (5)-(8). Коэффициенты уравнений вычисляются разложением функций Fi в ряд Тейлора с использованием первых двух членов ряда, если иметь в виду, что в стационарном состоянии
-^ст F2ст F3ст F4ст 0-
Таким образом, коэффициенты а, Ъ., с, di вычисляются по следующим формулам:
дS)„• ' КдхА, ' (др' (дм
Значения частных производных получаем с использованием значений £ст, Хст, Рст, Мст для каждого принятого стационарного состояния. Эти значения получены на этапе 2.
В результате имеем систему линейных дифференциальных уравнений первого порядка с известными постоянными коэффициентами.
Реализация этапа 4 позволяет вычислить элементы матрицы Гурвица, которая является основой оценки условий устойчивости. Матрица Гурвица для рассматриваемого объекта имеет вид
А =
| Р1 Р0 0 0 1
Рз Р2 Р1 Р0
0 Р4 Рз Р2
V 0 0 0 Р4 J
(9)
Приведем два варианта вычисления элементов матрицы Гурвица Рь Р2, Р3, Р4. Первый вариант реализуется путем преобразования системы четырех дифференциальных уравнений первого порядка с постоянными коэффициентами в одно дифференциальное уравнение четвертого порядка, т. е. система (6)-(9) преобразуется в уравнение
^+р ^+Р2 ^+Рз +РА = 0. (10)
dt dt3 dt dt
Уравнение (10) записано относительно переменной 81. Отметим, что преобразование возможно по любой из переменных.
Элементы матрицы Гурвица Р. вычисляются по формулам, приведенным в табл. 1. Соотношения табл. 1 получены в последовательности преобразования (5)-(8) в уравнение четвертого порядка (последовательность преобразования не приводится).
Таблица 1
Соотношения для расчета коэффициентов Р0, Р1, Р2, Р3, Р4
Соотношение
а. ^ ; fдF¡ ^ ; ^ 1 ; fдF¡ ^ Ч^ 1/ '"(ах )„' с' "(ЗР JCI; ' _|5м JCI
М1 = Ъ1а2 + с1а3 + d1a4; М2 = Ь1Ь2 + с1Ъ3 + d1Ъ4; М3 = Ъ1с2 + с1с3 + d1c4; М4 = Ъ^2 + с1й3 + d1d4
N1 = М2а2 + М3а3 + М4а4; N2 = М2Ъ2 + М3Ъ3 + М4Ъ4; N3 = М2с2 + М3с3 + М4с4; N4 = М^2 + M3d3 + М^4
Ll = N2а2 + N3a3 + N4a4; L2 = ^Ъ2 + N3Ъ3 + N4Ъ4; L3 = N2c2 + N3c3 + N4c4; L4 = N2d2 + N3d3 + N4d4
К1 = М2 + Ъ1а1; К2 = М1Ъ1 - М2а1; К3 = М3Ъ1 - М2с1; К4 = М4Ъ1 - M2d1
Rl = М1Ъ1 + R2 = N1Ъ1 - N^1; R3 = N3Ъ1 - ^сь R4 = N4Ъ1 - N2d1
Sl = N1Ъ1 + L2; S2 = L1Ъ1 - L2a1; S3 = L3Ъ1 - L2c1; S4 = L4Ъ1 - L2d1
Gl = RзЪl + КзаА; G2 = К^1 - К^з; Gз = К^2 - K2Rз; G4 = KзR4 - К4Кз
Fl = Sзbl + МКзЪи F2 = KзSl - KlSз; Fз = KзS2 - K2Sз; F4 = KзS4 - K4Sз
Р0 = KзG4Ъl; Р1 = -(KзF4Ъl + KзG4alЪl); Р2 = -^1 - GlF4); Рз = -^2 - G2F4); Р4 = -^з - GзF4)
Р1 = Р1 / Р0; Р2 = Р2 / Р0; Рз = Рз / Р0; Р4 = Р4 / Р0; Р0 = 1
Второй вариант реализуется следующим образом.
Система уравнений (5)-(8) записывается в векторной форме:
8' = А ■ 5,
где А - матрица коэффициентов уравнений системы (5)-(8):
A =
' a b, c1 d1 Л
a2 b2 c2 d 2
a3 b3 c3 d 3
v a4 К c4 d 4 ,
Используя собственные числа X и матрицу А, записываем характеристическое уравнение в виде
с ^
a - х b
a2 b2 - x c2 d2
a
b b
C3 — X d 3 c4 d4 —X
= 0.
Раскрывая определитель, получаем уравнение вида
X4 + рх3 + P2X2 + P3X + P4 = 0.
На этом завершается этап 4.
Этап 5 является завершающим в описании алгоритма, который дает ответ на вопрос об устойчивости или неустойчивости рассматриваемых стационарных состояний.
Заполняется матрица Гурвица (10) и вычисляются следующие определители:
A = |P|; A2
P Po P3 P2
A3 =
Pi Po P3 P2
0 Pi
0 P4 P3
A4 =
P P0 0 0
P3 P2 P P0
0 P4 P3 P2
0 0 0 P4
(11)
Необходимое и достаточное условие устойчивости
А > о, i = 1,4.
Если какое-либо из А, < 0, стационарное состояние неустойчиво. Использование алгоритма приведено в следующем разделе.
Анализ устойчивости стационарных состояний процесса получения молочной кислоты из пшеничной муки
Исходные данные для моделирования в достаточно полном объеме представлены в работе [10]. Уравнения математической модели для нестационарных условий процесса по непрерывной технологии:
dS = —-1 |X + KMM + DS — S),
dt
dX dt
= |X — DX,
dP Y
— = ^ |X — DP
dt YS
(12)
(13)
a
3
S
ш _ -КМЫ + В(мо -М). (15)
Л
Уравнение для удельной скорости роста записано в виде
V К + £
' Р >
1 —
и
Р
тах у
(16)
Учитывая, что К, << £ (в результате оценок по эксперименту) соотношение (16) используется в виде
'1—р у
Р
тах у
В соотношениях (12)-(16) ц, цтах - удельная и максимальная скорость роста, ч-1; X, Р, М -концентрация компонентов, г/л; £0, М0 - концентрация компонентов во входном потоке, г/л; В - величина протока, ч-1; У5, УР - стехиометрические коэффициенты, безразм.; Кт, Ртах, п - константы.
Численные значения констант: цтах = 0,28 ч-1; Ртах = 98,6 г/л; К, = 0,5 г/л; = 0,053; УР = 0,52; Км = 0,035; п = 3,0.
Область задания £0, М0 и В: £0 = 50-120 г/л; М0 = 20-80 г/л; Впред = 0,28 ч-1 (определено по условию материального баланса по X).
К рассмотрению принято два стационарных состояния: первое - произвольное в допустимых пределах по £0, М0 и В; второе - для оптимальных условий, которые определены для максимальной продуктивности Qp = ВР = maxQP по величине протока В01й. На этом этап 1 завершен.
Стационарное состояние 1. Значения показателей получены решением уравнений (12)-(15) по условию
л,_лх _ лр_лм _ 0
Л Л Л Л Условия на входе: £0 = 100 г/л; М0 = 50 г/л; В = 0,15 ч-1.
Показатели стационарного состояния: £ст = 86,8731 г/л; Хст = 1,1971 г/л; Рст = 18,5205 г/л; Мст = 40,5405 г/л; Qp = 2,7781 г/(л • ч).
Стационарное состояние 2. Получено по решению задачи оптимизации для тех же исходных значений £ и М0. Оптимальное значение В0^ = 0,1181 ч-1.
Значения показателей в стационарных условиях: £ст = 81,3675 г/л; Хст = 1,5932 г/л; Рст = 24,65 г/л; Мст = 38,5714 г/л; maxQP = 2,9118 г/( л • ч). Завершается этап 2.
Далее вычисляются коэффициенты системы уравнений первого приближения (5)-(8) а,, Ь, с,, Л,. Для вычисления коэффициентов составляем табл. 2.
Стационарное состояние 1: а1 = -0,15; Ь1 = -2,8302; с1 = 0,1269; = 0,035; а2 = 0,0; Ь2 = 0,0; с2 = -0,006727; Л2 = 0,0; а3 = 0,0; Ь3 = 2,3207; с3 = -0,25407; Л3 = 0,0; а4 = 0,0; Ь4 = 0,0; с4 = 0,0; Л4 = -0,185. Стационарное состояние 2:
а1 = -0,118125; Ь1 = -2,2288; с1 = 0,14405; = 0,035; а2 = 0,0; Ь2 = 0,0; с2 = -0,0076347; Л2 = 0,0; а3 = 0,0; Ь3 = 1,8276; с3 = -0,2369; Л3 = 0,0; а4 = 0,0; Ь4 = 0,0; с4 = 0,0; Л4 = -0,1531.
Таблица 2
Соотношения для вычисления коэффициентов а;, Ь¡, с¡, й,
ai = -D *=- D 1 Y 1 S „ "DX„ C Ys (РШх - Р„) di = Km
a2 = 0,0 b2 = 0,0 „ nDXT C2 Ys (Р„ - Ра) d2 = 0,0
a3 = 0,0 Y b3 = 3 Y S c=-Y nDX- - D 3 Ys (Ршах - РС1 ) d3 = 0,0
a4 = 0,0 b4 = 0,0 C4 = 0,0 d4 = -( Km + D)
Завершен этап 3.
Для вычисления элементов матрицы Гурвица использует второй вариант расчета. Для системы уравнений первого приближения вида (5)-(8) матрица А с учетом численных значений констант аь Ьь с, di имеет вид
A =
Характеристическое уравнение:
0 0 0
( a, bi c1 d1 ^
0 0 C2 0
0 Ьз C3 0
V 0 0 0 d4 у
bi c1 d1
-X C2 0
Ьз C3 - 0
0 0 d 4 -
= 1 ; Р1 = - (a1 +
= 0.
или л + Р1л + Р2л + Р3л + Р4=0, где
Р3 = а1 (Ь3с2 - c3d4) + Ь3с2ё4; Р4 = -а-р3с1й4.
Таким образом, для первого стационарного состояния имеем: а1 = -0,15; с2 = -0,006727; Ь3 = 2,3207; с3 = -0,25407; d4 = -0,185; Р0 =1,0; Р1 = 0,58907; Р2 = 0,12807; Р3 = 0,01228; Р4= 0,0004332.
Для второго: а1 = -0,118125; с2 = -0,0076347; Ь3 = 1,8276; с3 = -0,2369; d4 = -0,1531; Р0 = 1,0; Р! = 0,508125; Р2 = 0,096291; Р3 = 0,008069; Р4 = 0,000252. Завершен этап 4.
Используя Р, записываем матрицу Гурвица:
(Р1 10 0 ^
Г =
Рз P2 P 1
0 Р4 Рз Р2 0 0 0 Р4
4 У
По матрице Гурвица вычисляются определители (11):
- для первого стационарного состояния: А1 = 0,58907; А2 = 0,0634; А3 = 0,000628;
А4 = 2,72-10;
- для второго стационарного состояния: А1 = 0,50125; А2 = 0,04086; А3 = 2,646-10-4; А4 = 6,6684-10-8.
Поскольку определители для первого и второго состояния положительны, по необходимому и достаточному условию [1, 2] оба стационарных состояния устойчивы по первому приближению.
Заключение
В заключение отметим, что разработанный алгоритм может использоваться вне зависимости от сложности математического представления кинетических соотношений и числа уравнений математической модели.
Единственным ограничением является то, что анализ проводится для условий первого приближения, т. е. в малой окрестности стационарного состояния, что в целом не позволяет делать заключение об устойчивости процесса при наличии больших возмущающих воздействий.
СПИСОК ЛИТЕРА ТУРЫ
1. Демидович Б. П. Лекции по математической теории устойчивости. М.: Изд-во Моск. ун-та, 1998. 480 с.
2. Жукова Г. С., Митрохин С. И., Дарсалия В. Ш. Дифференциальные уравнения. М.: РХТУ им. Д. И. Менделеева, 1999. 366 с.
3. Бирюков В. В. Основы промышленной биотехнологии. М.: КолосС; Химия, 2004. 296 с.
4. Смирнов В. А. Пищевые кислоты. М.: Легкая и пищ. пром-сть, 1983. 240 с.
5. Гордеев Л. С., Кознов А. В., Скичко А. С., Гордеева Ю. Л. Неструктурированные математические модели кинетики биосинтеза молочной кислоты. Обзор // Теоретические основы химической технологии. 2017. Т. 51, № 2. С. 8-25.
6. Гордеева Ю. Л., Рудаковская Е. Г., Гордеева Е. Л., Бородкин А. Г. Математическое моделирование биотехнологического процесса периодической ферментации получения молочной кислоты. Обзор // Теоретические основы химической технологии. 2017. Т. 51, № 3. С. 1-18.
7. Vazquez J. A., Murado M. A. Unstructured mathematical model for biomass, lactic and bacteriocin productions by lactic acid bacteria in batch fermentation // J. Chem. Technol. Biotechnol. 2008. Vol. 83. P. 91-96.
8. Akerberg C., Hofvendahl K., Zacchi G., Hahn-Hagerdal B. Modelling the influence of pH, temperature, glucose and lactic acid concentrations on the kinetics of lactic acid production by Lactococcus lactis ssp. lactis ATCC 19435 in whole-wheat flour // Appl. Microbiol. Biotechnol. 1998. Vol. 49, no. 6. P. 682-690.
9. Hofvendahl K., Hahn-Hagerdal B. L-lactic acid production from whole wheat flour hydrolysate using strains of Lactobacilli and Lactococci // Enzyme Microb. Technol. 1997. Vol. 20, no 4. P. 301-307.
10. Gonzales K., Tebbano S., Lapes F., Thorigne A., Givry S., Dumar D., Pareau D. Modeling the continuous lactic acid production process from wheat flour // Appl. Microbiol. Biotechnol. 2016. Vol. 100, no. 1. P. 147-159.
Статья поступила в редакцию 19.01.2017
ИНФОРМАЦИЯ ОБ АВТОРАХ
Гордеева Юлия Львовна - Россия, 109472, Москва; Московская государственная академия ветеринарной медицины и биотехнологии им. К. И. Скрябина; канд. техн. наук, доцент; доцент кафедры информационных технологий, математики и физики; l.s.gordeev@yandex.ru.
Комиссаров Юрий Алексеевич - Россия, 125047, Москва; Российский химико-технологический университет им. Д. И. Менделеева, д-р техн. наук, профессор; профессор кафедры электротехники и электроники; komiss@muctr.ru.
Гордеева Елена Львовна - Россия, 125047, Москва; Российский химико-технологический университет им. Д. И. Менделеева; канд. техн. наук, доцент; доцент кафедры высшей математики; l.s.gordeev@yandex.ru.
Бородкин Алексей Георгиевич - Россия, 125047, Москва; Российский химико-технологический университет им. Д. И. Менделеева; доцент кафедры процессов и аппаратов химической технологии, l.s.gordeev@yandex.ru.
Yu. L. Gordeeva, Yu. A. Komissarov, E. L. Gordeeva, A. G. Borodkin
ALGORITHM OF USING HURWITZ CONDITIONS
TO STUDY STABILITY OF STEADY STATES OF CONTINUOUS FERMENTATION PROCESSES WITH NONLINEAR KINETICS
Abstract. The authors of the article produce the algorithm to evaluate the stability of the steady states for microbiological synthesis using mathematical models for kinetics based on the unstructured approach. The object under study is a continuous biotechnological process aiming to obtain a target product, in addition to biomass. The estimation of stability is performed using Hurwitz matrix. In this process the use of raw materials containing components for additional producing of the primary substrate are possible. The algorithm includes five steps: formulation of the mathematical model that takes into account biological and technological limits; definition of steady states (one or several) for which it is necessary to evaluate conditions of stability and to calculate process parameters for the adopted steady states resulted from solving a system of nonlinear algebraic equations of the mathematical model; formulation of the equations of the first approximation (the system of first-order linear differential equations with constant coefficients) and evaluation of the coefficients corresponding to the adopted steady states; formation of Hurwitz matrix and calculation of its items; calculation of necessary and sufficient conditions for stability using Hurwitz matrix. Two methods of calculating Hurwitz matrix have been shown. The first method is implemented by transforming the system of four differential equations of the first approximation into one fourth-order ordinary differential equation. In the course of transformation there were determined calculated ratios for coefficients of the matrix. In the second method the coefficients are calculated by solving vector equation for the system of first approximation with the introduction of eigen values. The algorithm has been used for evaluating steady states of microbiological synthesis for lactic acid production from raw materials containing the component reproducing the main substrate. The estimation has been carried out for two steady states, one of which is optimal. For numerical values of the constants (from the literature) both steady states were stable.
Key words: microbiological synthesis, evaluation of stability, Hurwitz matrix.
REFERENCES
1. Demidovich B. P. Lektsii po matematicheskoi teorii ustoichivosti [Lectures on mathematical theory of stability]. Moscow, Izd-vo Moskovskogo universiteta, 1998. 480 p.
2. Zhukova G. S., Mitrokhin S. I., Darsaliia V. Sh. Differentsial'nye uravneniia [Differential equations]. Moscow, RKhTU imeni D. I. Mendeleeva, 1999. 366 p.
3. Biriukov V. V. Osnovy promyshlennoi biotekhnologii [Introduction into industrial biotechnology]. Moscow, KolosS; Khimiia Publ., 2004. 296 p.
4. Smirnov V. A. Pishchevye kisloty [Nutritional acids]. Moscow, Legkaia i pishchevaia promyshlennost' Publ., 1983. 240 p.
5. Gordeev L. S., Koznov A. V., Skichko A. S., Gordeeva Iu. L. Nestrukturirovannye matematicheskie mod-eli kinetiki biosinteza molochnoi kisloty. Obzor [unstructured mathematical models of kinetics of lactic acid biosynthesis. Review]. Teoreticheskie osnovy khimicheskoi tekhnologii, 2017, vol. 51, no. 2, pp. 8-25.
6. Gordeeva Iu. L., Rudakovskaia E. G., Gordeeva E. L., Borodkin A. G. Matematicheskoe modelirovanie biotekhnologicheskogo protsessa periodicheskoi fermentatsii polucheniia molochnoi kisloty. Obzor [Mathematical modelling of biotechnological process of periodic fermentation for obtaining lactic acid]. Teoreticheskie osnovy khimicheskoi tekhnologii, 2017, vol. 51, no. 3, pp. 1-18.
7. Vazquez J. A., Murado M. A. Unstructured mathematical model for biomass, lactic and bacteriocin productions by lactic acid bacteria in batch fermentation. J. Chem. Technol. Biotechnol., 2008, vol. 83, pp. 91-96.
8. Akerberg C., Hofvendahl K., Zacchi G., Hahn-Hagerdal B. Modelling the influence of pH, temperature, glucose and lactic acid concentrations on the kinetics of lactic acid production by Lactococcus lactis ssp. lactis ATCC 19435 in whole-wheat flour. Appl. Microbiol. Biotechnol, 1998, vol. 49, no. 6, pp. 682-690.
9. Hofvendahl K., Hahn-Hagerdal B. L-lactic acid production from whole wheat flour hydrolysate using strains of Lactobacilli and Lactococci. Enzyme Microb. Technol., 1997, vol. 20, no 4, pp. 301-307.
10. Gonzales K., Tebbano S., Lapes F., Thorigne A., Givry S., Dumar D., Pareau D. Modeling the continuous lactic acid production process from wheat flour. Appl. Microbiol. Biotechnol., 2016, vol. 100, no. 1, pp. 147-159.
The article submitted to the editors 19.01.2017
INFORMATION ABOUT THE AUTHORS
Gordeeva Yulia Lvovna - Russia, 109472, Moscow; Moscow State Academy of Veterinary Medicine and Biotechnology named after K. I. Skryabin; Candidate of Technical Sciences, Assistant Professor; Assistant Professor of the Department of Information Technologies, Mathematics and Physics; l.s.gordeev@yandex.ru.
Komissarov Yuriy Alekseevich - Russia, 125047, Moscow; Mendeleev Russian Chemical and Technological University; Doctor of Technical Sciences, Professor; Professor of the Department of Electrical Engineering and Electronics; komiss@muctr.ru.
Gordeeva Elena Lvovna - Russia, 125047, Moscow; Mendeleev Russian Chemical and Technological University; Candidate of Technical Sciences, Assistant Professor; Assistant Professor of the Department of Higher Mathematics; l.s.gordeev@yandex.ru.
Borodkin Aleksey Georgievich - Russia, 125047, Moscow; Mendeleev Russian Chemical-Technological University; Assistant Professor of the Department of Processes and Apparatus for Chemical Technology; l.s.gordeev@yandex.ru.
