Моделирование и оптимизация технологических систем с распределенными параметрами Текст научной статьи по специальности «Математика»

УДК 62-501.6

Вестник СибГАУ 2014. № 3(55). С. 55-62

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

Н. Д. Демиденко1, Л. В. Кулагина2

1 Специальное конструкторско-технологическое бюро «Наука» Красноярского научного центра Сибирского отделения Российской академии наук Российская Федерация, 660036, г. Красноярск, Академгородок, 50 2Сибирский федеральный университет Российская Федерация, 660041, г. Красноярск, просп. Свободный, 79 Е-mail: [email protected]

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

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

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

Ключевые слова: математическое моделирование, системы с распределенными параметрами, оптимальное управление, тепломассообмен.

Vestnik SibGAU 2014, No. 3(55), P. 55-62

SIMULATION AND OPTIMISATION IN TECHNOLOGICAL SYSTEMS WITH DISTRIBUTED PARAMETRES

N. D. Demidenko1, L. V. Kulagina2

Special Design and Technological Bureau "Nauka" Krasnoyarsk Scientific Centre of Siberian Branch Russian Academy of Sciences, 50, Akademgorodok, 660036, Krasnoyarsk, Russian Federation

2Siberian Federal University 79, Svobodny prosp., 660041, Krasnoyarsk, Russian Federation Е-mail: [email protected]

The article presents a method of simulation and optimal control for rectification plants, containing technological furnaces and rectification columns. This complicated system is shown as a system with distributed parameters, as for static and dynamic modes simulation a mathematical formalism of differential equations in partial differential was applied. As a rule scientific literature deals with operation modes outcomes for single plants. Due to complicated mathematical formalism and computation investigations dealing with a number of objects with distributed parameters are infrequently. The optimal control tasks solutions for complicated plants are poorly shown in scientific literature A phe-nomenological approach was applied for mathematical model of heat processes in technological furnaces and separation processes in multicomponent mixtures for rectification columns. The model contains partial differential equitations for heat- and mass-exchange processes and hydrodynamics for flue gas flows in technological furnace and massexchange equitations for rectification column. Described processes occurs with recirculation flow interaction and heat-exchange processes' boundary conditions for rectification columns are described in ordinary derivatives. The conditions are given at the opposite ends of the plant. The optimal control modes solutions were formulated an optimal control task for desired product and necessary conditions for optimality. Due control functions are included in the main boundary derivatives, control variations in the range and on the boundary are dependent, this works on dual problem structure. Arguments of variational calculus was applied. Necessary conditions for optimality are containing original boundary-value problem conjugated relatively to the Lagrangians. Original and boundary tasks solution determines optimal control and technological processes parameters. The article gives numerical outcomes with control of flow of raw materials in debutanizer plant for sulfuric acid alkylation of isobutane with butenes. A numeric algorithm was worked out. The peculiarity of the algorithm is solving a boundary conditions tasks for some parameters at the opposite ends of the area. Moreover Lagrangians for conjugated system are given at a finite time. Optimal characteristics of heat-exchange process of desired product in the dephlagmator and cube were obtained.

Keywords: mathematical modeling, systems with distributed parameters, optimal control, heat-and-mass exchange.

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

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

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

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

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

Уравнения процессов в технологической печи. При исследовании процесса горения капель жидкого топлива в воздухе в основном представляет интерес распределение концентраций компонентов, плотности, температуры сырья, температуры и скорости дымовых газов в печи при статических и динамических режимах работы. Исходя из одномерности движения потока, математическая модель нестационарного горения может быть представлена следующими уравнениями [1]:

1. Уравнение неразрывности

^ + = 0, (1)

дt д£

где р - массовая плотность смеси, и - скорость движения смеси.

Для покомпонентной модели процесса горения уравнение (1) можно записать в виде

д(РХ) + д(рхи) =_ Р£ (2)

дг 31 т , )

где I - линейный размер; х - концентрация горючего вещества в смеси (0< х <1); т - время сгорания. 2. Уравнение движения в виде

' ди ди Л дР р I — + и— 1+— = 0,

1 дг д1) д1

где Р - давление.

3. Уравнение сохранения энергии

р^!+и I Л=Рт' - ? т)+

+К1 (т1 - т-)+к (т2 - т-),

(3)

(4)

где ' - теплота сгорания топлива; Тп - температура дымовых газов; б(ТП) - потери на излучение; Т.1,

Т 2

температура нагреваемого сырья, движущегося

в нижнем и верхнем направлении; К. - коэффициент теплопередачи; - энтропия, причем

5 = с, 1пР ( = 1 -1,4). РУ

4. Уравнение теплообмена между нагреваемым сырьем, движущимся в низ и верх печи, и дымовыми газами можно записать в следующем виде:

дТ1 дТ„ —— - w

дг д£

^ = К2 (Тп - Т;1)- б(Тп);

дТ;

2

дг

- + w-

дТ;

2

д£

= К2 (Тп - Т;2 )- б(Тп),

(5)

(6)

К2 - коэффициент теплопередачи; w - скорость течения сырья.

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

а б

Рис. 1. Схема ректификационной установки: а - технологическая печь; б - ректификационная колонна; 1 - камера радиации (топочная камера); 2 - камера конвекции; 3 - дымовая труба; 4 - конвекционные трубы; 5 - радиантные трубы; 6 - вход сырья; 7 - выход сырья; 8 - дымовые газы

Начальные условия:

р(£,0) = р0, х(£,0) = х0, и(£,0) = и 0, тп(^,0) = тп0,Тс1 (£,0) = т^,Т2 (^,0) = Гс20. Граничные условия:

р(0, г) = а1, х(0, г) = а 2, и (0, г) = а 3, Тп (0, г) = а4, Тс1 (Ь, г) = а5,Т? (0, г) = Т,1 (0,г) = <

(7)

(8)

где Ь - длина печи.

На базе системы (1)-(8) можно сформулировать и решить задачи оптимального управления с различными управляющими параметрами [2-6].

Расчет динамических режимов трубчатых печей. Рассмотрим следующую тепломассообменную задачу для процессов в трубчатой печи. Для этого приведем систему (1)-(6) к следующему виду:

др др ди — = _и — _ р—,

дг д£ д£

дх дх х

дг ~ д£ т,

ди = ди ^ дТП ЯТП др

дг " д£ д£ р д£,

дТп ¡л \гг ди дТп

—п = (1_ У )Тп--и—п +

дг У ' п д£ д£

х С т

6 (Тп, Тс1, Тс2

)

С-Р

+К (Тс1 _ Тп)+к (Тс2 _ Тп ),

(9)

дТ1 дТ

= м-

+ К (Тп _ТС1 )-

дг д£ _6 (Тп), г е [0,Т], I е [0,Ь],

дтс

дТс2

дг

_6 (Тп ),

д£ г е[0, т]

+ К (Тп _ТС2)-

£ е[0, Ь1

где Ь1 - точка вывода сырья из печи. Таким образом, сформулирована краевая задача (7)-(9). Здесь температура сырья Т.1 задается в точке I = Ь, так как сырье подается сверху в печь, и таким образом мы имеем противоточный технологический процесс.

Кривые разгона на выходе печи получены при возмущении на ± 20 % с шагом 5 % на входе печи по температуре сырья и по температуре дымовых газов (рис. 2-5).

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

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

200 400 600 N00 1000 '.С

Рис. 2. Кривые разгона по плотности потока в зависимости от температуры потока (от 424 до 636 °С)

Рис. 3. Кривые разгона по скорости потока дымовых газов в зависимости от температуры сырья (от 216 до 324 °С)

Рис. 4. Кривые разгона по температуре дымовых газов в зависимости от температуры сырья (от 216 до 324 °С)

Щх^ - = k (y - У( x))+f № W xf,

) d(Vy) , * ч

V y '--^ = k (y (x) -y), 0 < t < T, 0 < £ < L.

dt d£ v '

Краевые условия: при € = 0

d ( H xk Xk )

(10)

dt

= L (0, t) x (0, t) - V (0, t)y (0, t) - Wxk (t), (11)

y (^t) = a [y* (xk ) - xk (t)] + xk (t), при t = L d (Hxx )

dt

■ = Vd (t) yd (t)-( (t) + D (t))xd (t),

Vdyd (t) - V (L, t)y (L, t) = Ldxd (t) - L (L, t) x(L, t), (12)

yd (t) = у (L, t) + Ed [y* (xd) - x (L,t)].

Начальные условия:

x(£,0) = x0 (£), y (i,0) = y0 (£), 0 < £ < 1

xd (0) = xd0, xk (0) = xt0, 0 < a, Ed < 1.

(13)

ности в жидкости и паре; ха - концентрация целевого продукта в дефлегматоре; хк - концентрация целевого продукта в кубе. Возмущающее воздействие по концентрации целевого продукта в сырье Е2 (рис. 6, а); у* - равновесная концентрация целевого продукта в паровой фазе; Б, Ш - отбор целевого продукта вверху и внизу колонны; - орошение.

Величины потоков при этом удовлетворяют условиям

W(t) + D (t) = F(t), D (t) + Ld (t) = V (1,t), W(t) + V (0,t) = L(0,t), Vd (t) = V (1,t).

(14)

Рис. 5. Кривые разгона по температуре сырья в зависимости от температуры потока газов (от 424 до 636 °С)

Для получения необходимых условий оптимальности (условий стационарности) используются методы классического вариационного исчисления [5-15]. Управления предполагаются кусочно-непрерывными, а соответствующие им решения - непрерывными и кусочно-гладкими.

Здесь рассматривается следующая модель процесса [16]:

Предполагается, что удерживающие способности Hx, Hy постоянны, V не зависит от t.

Анализ условий (14) показывает, что только два из четырех потоков (W, D, L, V) являются независимыми. Выбор тех или иных двух независимых потоков в качестве управлений определяет соответствующие задачи. Эти задачи рассматриваются как задачи оптимального управления в классе кусочно-непрерывных управлений с критерием качества:

TL 2

I = Ц(у(4/)-е*(I,/)) )t, (15)

0 0

где е* - заданное значение концентрации целевого продукта.

Потоки L, D фиксированы, а V, Vd, Ld исключаются согласно (14). Управление F выбирается в классе кусочно-непрерывных функций и принимает значение в промежутке

F < F(t) < F

1 mm — J W — J max*

Переходя к нормальной форме дифференциальных уравнений (10)-(12), получаем следующую задачу. Рассматривается процесс ректификации, описываемый уравнениями

1

xt =

H.

yt =

Hy

(L + L )1 x + k (y -/)

+ Fxf Фх X, x 'г = Z(1), (D + L)Z(2) + k (( - y)) 7

(16)

y '/= Z

= Z(2)

при краевых условиях

x't— [((+F)x-Vy-Wk > Xk, Hk

y=a[y* - xk ]+xk, £= 0, 0 <t<T; x'dt ~ [(D+L) -(L+D)xd >Xd,

(17)

£=L, 0<t<T,

xd

Здесь величина Е(г) - поток сырья в жидкой фазе, подводимый в колонну и являющийся управлением; х (Ь, г), у (Ь, г) - концентрации целевого продукта

в жидкой и паровой фазах; Ь(1,г) - поток жидкости, У(1, Н) - поток пара; Нх, Ну - удерживающие способ-

(Б+Ь)( - У) - Ь( - х) = 0, уа - у - Ей (у* - у) = 0,

при начальных условиях (13) и ограничениях на управления

(Р - )(^тах - Е)- и2 = 0, (18)

где и - вспомогательное управление.

Задача состоит в том, чтобы в множестве кусочно-непрерывных функций F, удовлетворяющих условию (18), найти такую, что соответствующее ей решение задачи (16)—(18), (13) дает минимум интегралу (15). Применяя известную процедуру вариационного исчисления, получаем необходимое условие оптимальности:

Ьг -

L + L

Hx

V

Ь/ = k (у* )'

Пг +ТГ П = к

Hy

при I = о, о < г < T

ах к1) _ х к1)

н.

1)

ь ( - В ) + х <?>( -1 - а (у* )),

хк

х (1)(т)_ о,

X ^(Ь + F )-Ь (ь + £ )_ 0.

V

_п__

КНУ НХку

- X к2)_ о;

при г = Т, о < I < Ь, Ь = о, п = 0; при I = ь, о < г < т

ах

(1)

^ ^х^-ьх^), х<?>(Т)_ о,

н

хй

V .х(}) + х(2V+х(3)_ о,

НX,

( , \ х (2) Ь х й -

н.

Ь —

хй \

н

-ь - х<4 (у*) _ о,

хй

V НУ

-- х

V - х(й3)(1 - Ей)_ о;

при о < г < т

н _/-

н

л

г д2 Ь _ дЬ г (1)

-X + Ф х xF--С

^дЕ х др у

м -

х1

(1)

н

-(Хк -X) + у(тт + Fmax -2F)_ о,

Хк

ум _ о,

(19)

(2о)

(21)

(22)

3) далее полагаем ^п+1) = ^п) - аН;

4) предельные значения Е(и+1) при п ^ да дают оптимальные управления.

На рис. 6 приведены результаты расчетов по оптимальному управлению для промышленной колонны К-34 [17-2о] установки сернокислотного алкилирова-ния изобутана бутиленами (разделяемая многокомпонентная смесь сведена к бинарной). Основные параметры: В = 27,об кмоль/ч, Ж = 76,59 кмоль/ч, Ьй = 45,21 кмоль/ч, Нхй = 5о кмоль, Нхк = 3о кмоль.

. кол . Золи

о,г

Ща,

0,8Щ

ОМ ■Якг

0,28

0,20

0,12

5 ¡60

80

, чол 9м ц йг 1 _1 1 0

•

, пАЛ Золи ■ $¿2 / А2

, 1 , !

у «т

-— 1 1 . 1

ГО

го

где Ь, П, х(1)й, х(2)й, х(3)й, х(1)к, х(2)к - множители Лагранжа.

Алгоритм решения задачи оптимального управления содержит следующие этапы:

1) задается начальное приближение управляющей функции F0(t);

2) решается система уравнений (1о)-(13) и (19)-(21);

Рис. 6. Графики изменения концентрации бутана в сырье (а), дефлегматоре (б) и кубе (в) при управлении потоком сырья (г)

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

а

б

в

г

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

Библиографические ссылки

I. Демиденко Н. Д. Моделирование статических и динамических режимов в трубчатых печах // Управление, вычислительная техника и информатика : Вестн. Томского гос. ун-та. 2012. № 3(20). С. 13-21.

2. Демиденко Н. Д., Кулагин В. А., Шокин Ю. И. Моделирование и вычислительные технологии распределенных систем. Новосибирск : Наука. 2012. 424 с.

3. Демиденко Н. Д. Оптимальное управление режимами технологических печей // Вестник СибГАУ. 2013. № 3 (49). С. 182-187.

4. Демиденко Н. Д., Кулагина Л. В. Численный метод исследования стационарных режимов в технологических печах // Engineering & Technologies : J of Siberian Federal University. 2014. 1(7). 55-61.

5. Демиденко Н. Д., Кулагина Л. В., Мельник И. Н. Анализ нестационарных режимов в системах контроля и управления распределенными параметрами // Информатика и системы управления : межвуз. сб. науч. тр. Вып. 9. Красноярск: ГУ НИИ информатики и процессов управления. 2004. С. 135-144.

6. Демиденко Н. Д., Кулагина Л. В. Оптимальное управление режимами работы технологических печей в нефте-перерабатывающей промышленности // Фундаментальные исследования. 2005. № 2. С. 43-44.

7. Демиденко Н. Д. Моделирование и оптимизация тепломассообменных процессов в химической технологии. М. : Наука. 1991. 240 с.

8. Демиденко Н. Д., Кулагина Л. В. Методы и средства оптимального управления теплотехнологи-ческими процессами в трубчатых печах // Химическое и нефтегазовое Машиностроение. 2006. № 3. С. 8-9.

9. Демиденко Н. Д., Потапов В. И., Шокин Ю. И. Моделирование и оптимизация систем с распределенными параметрами. Новосибирск : Наука, 2006. 551 с.

10. Дьяконов С. Г., Елизаров В. И., Лаптев А. Г. Математические основы и моделирование процессов разделения веществ. Казань : Изд-во Казан. ун-та, 1993. 437 с.

II. Шахтахтинский Т. Н., Келбалиев Г. И. Аналогия и подобие нестационарных процессов конвективного тепломассопереноса // Докл. РАН. 1990. Т. 315, № 6. С. 1427-1429.

12. Кафаров В. В., Громов В. Ю., Матвейкин В. Г. Задачи управления объектами химической технологии при наличии нечеткости // Докл. РАН. 1994. Т. 337, № 5. С. 628-630.

13. Кафаров В. В., Громов В. Ю., Матвейкин В. Г. Математическое моделирование не полностью наблюдаемых химико-технологических объектов // Докл. РАН. 1994. Т. 337, № 1. С. 68-69.

14. Самарский А. А., Михайлов А. П. Математическое моделирование: Идеи. Методы. Примеры. М. : Физматлит, 2001. 320 с.

15. Дворецкий С. И., Егоров А. Ф., Дворецкий Д. С. Компьютерное моделирование и оптимизация технологических процессов и оборудования : учеб. пособие. Тамбов : Изд-во Тамб. гос. техн. ун-та, 2003. 224 с.

16. Демиденко Н. Д., Кулагина Л. В. Особенности сжигания топлива в технологических установках // Химическая техника. 2006. № 7. С. 35-37.

17. Кулагина Л. В. Проблема ингерентности и ликвидности перспективных энергогенерирующих систем // Социальные проблемы инженерной экологии, природопользования и ресурсосбережения : Материалы Всерос. науч.-практ. конф. Вып. IX. Красноярск : Краевое НТО, 2003. С. 45-56.

18. Демиденко Н. Д. Управляемые распределенные системы. Новосибирск : Наука. 1999. 392 с.

19. Демиденко Н. Д., Кулагина Л. В. Повышение эффективности и экологической безопасности сжигания топлива в технологических установках // Энергоэффективность систем жизнеобеспечения города : Материалы VI Всерос. НПК. Красноярск : ИПЦ КГТУ. 2005. С. 276-281.

20. Кулагина Л. В. Анализ теплофизических и гидродинамических процессов при сжигании различных видов топлив в энергетическом котле КВ-ТМ-180-150-25 // Социальные проблемы инженерной экологии, природопользования и ресурсосбережения : Материалы Всерос. НПК. Вып. X. Красноярск : Краевое НТО. 2004. С. 29-45.

References

1. Demidenko N. D. [Modelirovanie staticheskih i dinamicheskih rejimov v trubchatih pechah]. Upravlenie, vichislitel'naya technika i informática. Vestnik Tomskogo gosudarstvennogo universiteta. 2012, no. 3(20), p. 13-21. (In Russ.)

2. Demidenko N. D., Kulagin V. A., Shokin Y. I.

Modelirovanie i vichislitel'nie tehnologii raspredelenih sistem. [Simulation and Computing Technologies for Distributed Systems]. Novosibirsk, Nauka Publ., 2012, 424 p.

3. Demidenko N. D. [Optimal'noe upravlenie rejimami technologicheskih pechei]. Vestnik SibGAU. 2013, no. 3 (49), p. 182-187. (In Russ.)

4. Demidenko N. D., Kulagina L. V. [Chisleniy metod isledovaniya stacionarnih rejimov v technologicheskih pechah]. Journal of Siberian Federal University. Engineering & Technologies. 2014, Vol. 7, Issue 1, p. 55-61. (In Russ.)

5. Demidenko N. D., Kulagina L. V., Mel'nik I. N. [Non-stationary modes analysis in control and management systems with distributed parameters]. Informatika I sistemi upravleniya. 2004, no. 5, p. 135-144. (In Russ.)

6. Demidenko N. D., Kulagina L. V. [Optimal control of the operation modes in technological furnaces for refining industry]. Fundamental'nie issledovaniya. 2005, no. 2, p. 43-44. (In Russ.)

7. Demidenko N. D. Modelirovanie i optimizaciya teplo-masoobmehhih procesov v chimicheskoy technologiy. [Simulation and Optimization for Heat-

Mass-exchange Processes in the Chemical Engineering] Moscow, Nauka Publ., 1991, 240 p.

8. Demidenko N. D., Kulagina L. V. [Means and facilities for optimal control of the technological processes in tube furnaces]. Himicheskoe i neftegezovoe machinostroenie. 2006, no. 3, p. 8-9. (In Russ.)

9. Demidenko N. D., Potapov V. I., Shokin U. I. Modelirovanie i optimizaciya sistem s raspredelennimi parametrami [Simulation and optimization for systems with distributed parameters]. Novosibirsk, Nauka Publ., 2006, 551 p.

10. D'yakonov S. G., Elizarov V. I., Laptev A. G. Matematicheskie osnovi i modelirovanie procesov razdeleniya veshestv [Mathematical foundations and simulation of processes of separation of substances]. 1993, Kazan, Kazan University Publ., 437 p.

11. Shahtahtinskiy T. N., Kelbaliev G. I. [Analogie and similarity of non-stationary processes of convective heat-and mass-transfer]. Doklady RAN. 1990, Vol. 315, no. 6, p. 1427-1429. (In Russ.)

12. Kafarov V. V., Gromov V. U., Matveikin V. G. [Tasks of objects' management in chemical engineering while fuzziness] Doklady RAN. 1994, Vol. 337, no. 5, p. 628-630. (In Russ.)

13. Kafarov V. V., Gromov V. U., Matveikin V. G. [Mathematical simulation for not fully observed chemical and technological properties]. Doklady RAN. 1994, Vol. 337, no. 1, p. 68-69. (In Russ.)

14. Samarskiy A. A., Mihailov A. P. Matematicheskoe modelirovanie: idei, metodi, primeri [Mathematical simulation: ideas, methods, examples]. Moscow, Fizmatlit Publ., 2001, 320 p. (In Russ.)

15. Dvoreckiy S. I., Egorov A. F., Dvoreckiy D. S. Komp'uternoe modelirovanie i optimizaciya tehnologicheskih procesov I oborudovaniya [Computer simulation and optimization of industrial processes and hardware]. Tambov, Tambov. State. tehn. University Publ., 2003, 224 p.

16. Demidenko N. D., Kulagina L. V. [Features of fuel burning in the technological furnaces]. Chimicheskaya technika. 2006, no. 7, p. 35-37. (In Russ.)

17. Kulagina L. V. [Ingerente and liquidity tasks in perspective power generation systems]. Materialy Vseros. nauchyu-praktyu konf. "Sotsial'nyye problemy inzhenernoy ekologii, prirodopol'zovaniya i resursosberezheniya". Vyp. IX. [Proceedings of All-Russia. Teach-practical conference "Social problems of engineering ecologie, environmental management and resource saving"]. Issue. IX, Krasnoyarsk, Krayevoye NTO Publ., 2003, p. 45-56. (In Russ.)

18. Demidenko N. D. Upravlyaemie raspredelennie sistemi [Controlled distributed systems]. Novosibirsk, Nauka Publ., 1999, 392 p. (In Russ.)

19. Demidenko N. D., Kulagina L. V. [Improving the efficiency and environmental security of fuel combustion in industrial plants]. Materialy VI Vseros. nauch.-prakt konf. "Energoefectivnost' sistem jizneobespecheniya goroda" [Proceedings of the VI All-Russia. scientific-practical conference "Energy efficiency of city lifesupport systems"]. Krasnoyarsk, IPTS KGTU Publ., 2005, p. 276-281. (In Russ.)

20. Kulagina L. V. [Analysis of thermal and hydrodynamic processes during combustion of different fuels in the power boiler KV-TM-180-150-25] Materiali Vseros. nauch.-prakt konf. "Social'nie problemi injenernoy ekologii, prorodopol'zovaniya iresursosberejeniya" [Materials Vseros. scientific-practical conference "Social problems of environmental engineering, environmental management and resource saving"]. Issue. X. Krasnoyarsk, 2004, p. 29-45. (In Russ.)

© fteMHgemo H. ft, Kyrarma H. B., 2014