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

Серия 4. Химическое машиностроение и инженерная экология

2. Бахвалов Н.С., Жидков Н.П., Кобельков Г.М. Численные методы. - М.: Лаборатория Базовых Знаний, 2002 г. - 632 с.

3. Гданский Н.И. Геометрическое моделирование и машинная графика. - М.: МГУИЭ, 2003 г. - 236 с.

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

Артюшкин А.Ю., к.т.н. доц. Воронина П.В., к.т.н. доц. Татаринов A.B. ГБОУ ВПО «Московский государственный университет пищевых производств»

89037702868, atp@mgupp.ru

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

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

Введение

На современных химико-технологических производствах широко распространены аппараты периодического приципа действия, в которых реализуются биохимические и химические процессы, проводимые в условиях натуральной конвекции. Примерами могут служить биореакторы-ферментаторы, применяемые при производстве напитков брожения, дрожжевом производстве, производстве ферментных и медицинских препаратов и т.п. Такие технологические аппараты имеют обычно цилиндрическую или цилиндро-коническую форму и снабжены несколькими рубашками охлаждения (или нагревания, в зависимости от требуемых условий протекания процесса). В промышленности подобные аппараты принято называть танками. Обеспечение свободного конвективного движения жидкой среды в танке (термогравитационной конвекции или конвекции Рэлея) в течение всего периода процесса является зачастую обязательным условием протекания последнего. Так, в пищевой индустрии широкое распространение получили танки значительного объема (50-100 м3) с системой из трех-четырех охлаждающих рубашек, встроенных по вертикальной координате. Отметим, что время протекания (био)технологических процессов в танках обычно значительно и может составлять 10-15 суток [1].

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

Серия 4. Химическое машиностроение и инженерная экология

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

Будем использовать для математической постановки задачи принципы термодинамики при конечном времени [2], согласно которым между двумя взаимодействующими за заданный (и конечный) временной интервал термодинамическими подсистемами, имеющими разную температуру, возникает тепловой поток, зависящий от термодинамической движущей силы. Диссипацию энергии при теплообмене в этом случае можно определить количественно с помощью производства энтропии о, представляющую собой по определению неотрицательную величину и характеризующую непроизводительные потери работоспособной энергии [3]. Практический интерес представляет усредненное за временной интервал протекания процесса производство энтропии а , которое для совершенных с точки зрения термодинамики при конечном времени систем должно достигать (или быть максимально приближено) к предельной усредненной оценке аШ1П, то есть минимально возможному значению. Отметим, что указанная предельная оценка есть величина положительная и обращается в ноль только для идеализированных, обратимых процессов. Разность а - аШ1П определяет степень несовершенства системы.

Будем рассматривать задачу достижения минимально возможного а как оптимизационную с критерием в форме функционала

а = \ [X(0,©0)/(0,©0)^ ^ шт, (1)

1 ®0

где t - независимый аргумент (календарное время);

Т - заданное и конечное время протекания процесса,

0 - усредненная по пространственным координатам фазовая переменная, в рассматриваемом случае - температура жидкой среды в танке. Отметим, что такое усреднение оправдано самой формой критерия;

©0 - управление, в рассматриваемом случае - изменяющаяся во времени температура хладагента, находящегося в рубашке. На управление наложено автономное ограничение ©о е , то есть ©отш < ©о < &отах , гДе ©о™« и ©о^ - заданные и постоянные величины;

I(©, ©о) - тепловой поток, определяемый характером теплообмена;

X (©, ©о) - термодинамическая движущая сила, обусловливающая интенсивность теплового потока.

Отметим, что в соответствии с законом Гюи-Стодолы потери полезной работы, имеющие место в необратимых процессах, прямо пропорциональны приращению энтропии участвующих в процессах теплообмена сред [4]. Поэтому максимальный термодинамический коэффициент полезного действия соответствует минимальному производству энтропии, то есть наименьшему значению задачи (1).

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

I = «(0-©о), (2)

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

X = (1 ©о -1 ©). (3)

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

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

Серия 4. Химическое машиностроение и инженерная экология ^ (0,0О)dt = I = const, (4)

где I - заданная, постоянная и положительная величина

Во-вторых, в постановку задачи необходимо ввести определяющие связи, обусловливающие динамику процесса теплообмена, то есть выражения, описывающие математическую модель процесса. Отметим, что динамика рассматриваемых объектов корректно может быть описана интегро-дифференциальными уравнениями в частных производных, аргументами которых, помимо времени, являются пространственные координаты. Для цилиндрических и цилиндро-конических аппаратов принято применять цилиндрическую систему координат, аргументами которой в общем случае являются: высота слоя среды z , радиальная координата r и азимутальная координата (р. В инженерных приложениях принято пренебрегать азимутальной координатой и рассматривать три независимых аргумента: время, вертикальную и радиальную составляющую пространственных координат. Кроме того, для процессов, протекающих в условиях естественной термогравитационной конвекции, для плотности среды применяют приближение Буссинеска [5], согласно которому плотность среды считают функцией ее температуры, то есть

р = р° (l-Д0-00)), (5)

где: р - плотность среды в текущем сечении относительно координаты z,

р° - плотность среды при так называемой отсчетной температуре 0°; в качестве которой в нашем случае может быть выбрана температура среды в нижнем гоизонтальном сечении аппарата,

р - коэффициент объемного теплового расширения.

По Буссинеску, при малой величине /5 и малой вариации параметров среды (кинематической вязкости /и, температуропроводности % и др.) плотность и сами эти параметры считают постоянными во всех членах уравнений, кроме одного: вариация плотности сохраняется там, где она умножается на ускорение силы тяжести g [6].

Использование упрощающего приближения Буссинеска позволяет выписать математическую модель объекта в форме системы уравнений Навье-Стокса

dv V~ ~

— + v-Vv = - gfiQ + juAv, (6)

50

— + v V0 = ^A0, (7)

dt

Vv = 0. (8)

В выражениях (6) - (8) обозначено: v - v(vr, vz) - вектор скорости конвективного потока, v • Vv - конвективное ускорение,

0 = 0-0°и ~ = p -p0 - вариации температуры и давления соответственно. Символами V и А = V V обозначены операторы Гамильтона и Лапласа соответственно. Таким образом, постановка рассматриваемой оптимизационной задачи должна включать в себя критерий (l) и ограничения (4), (6) - (8).

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

Серия 4. Химическое машиностроение и инженерная экология

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

Откажемся от традиционной схемы решения поставленной задачи. Вместо этого будем использовать двухэтапный алгоритм. На первом этапе отбросим автономные ограничения на управляющую переменную и связи (б) - (8). Таким образом, действуя по логике В.Ф. Кротова, перейдем к расширенной задаче [7], включающей в себя критерий и интегральную связь (4). Аналитическое решение расширенной задачи, то есть оптимальное значение фазовой переменной ©*, будет являться оценочным для исходной постановки. На втором этапе будем осуществлять итерационный поиск наиболее близкого к оценочному физически реализуемого решения.

Для невырожденного случая расширенной задачи условие оптимальности примет вид

0о = а^шт Я(0*, 0о,/), (9)

где: Я - функция Лагранжа,

Л - множитель Лагранжа (постоянный коэффициент).

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

Я = 1 Т

С 2

а(&~&о) + л(а(&-&о) -1)

0©о

(Ю)

Необходимые условия оптимальности для (1о) в форме требования стационарности Я по ©о примут вид

дЯ _ 2(©-©о)0©о + (0-0о)20

= -^-^-^--аЯ = о. (11)

5©о (0©о)2 ( )

Отметим, что в поставленной задаче температуры представлены в шкале Кельвина, кроме того, нами принято, что © > ©о V t е [о;Т]. Учитывая эти обстоятельства, после

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

©

®*о =

лА+т©

(12)

Отметим, что в нашем случае Л < о.

Оценочное значение X может быть рассчитано из соотношения

1Т

где константу С определим из выражения

0®

С = [ (14)

•о М + М 4 '

0 + 0о

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

окажется близким к I*, то полученное допустимое решение (допустимое управление) является оптимальным и задача решена. Значительная разница между ^ и I свидетельствует об эффективности расширения.

Серия 4. Химическое машиностроение и инженерная экология Шаги алгоритма второго этапа решения.

1) Разбиваем время процесса на N интервалов Т, исходя из априорных сведений о динамике процесса.

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

1

©о =

л/1 -Л

(15)

где: ©0 = ©0/ ©, Л = Л © - безразмерные показатели. Таким образом, стартовое значение Л может быть выбрано, исходя из желамого (но не обязательно достижимого) соотношения между 0 и ©0.

1) Численно решаем с внутренним шагом At систему уравнений (6) - (8) и осуществляем операцию усреднения полученных значений © по координатам. При этом для каждого внутреннего шага перед процедурой численного решения пересчитываем ©0 и используем его для последующего внутреннего шага. Таким

образом, для первого заданного интервала Т получаем набор табулированных

значений © и ©0 в соответствии с внутренним шагом моделирования. Отметим,

что поскольку имеет место итеративный счет, целесообразно применить обычные инженерные допущения, упрощающие вид зависимостей (6) - (8).

2) Рассчитываем, используя накопленные в ходе решения модельных уравнений значения © и ©0 и константу С (14). Затем с помощью выражения (13) получаем

расчетное значение средней интенсиваности потока I.

5) Сравниваем заданное и полученное расчетное значение для средней интенсивности потока: 1I -11< , где - наперед заданное малое положительное число. Если условие не выполнено, то следует изменить значение Л и повторить шаги 2 ^ 4.

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

При переходе к следующему временному интервалу Тм полученное на предыдущем

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

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

Сделанный авторами вычислительный эксперимент, основанный на реальных данных процесса брожения пива в цилиндро-коническом танке (первый этап брожения), [8]

Серия 4. Химическое машиностроение и инженерная экология подтвердил работоспособность предлагаемого алгоритма.

Заключение

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

Литература

1. Кунце В. Технология солода и пива. - СПб: Издательство Профессия, 2008. - 1200 с.

2. Математические методы термодинамики при конечном времени/ В.А. Миронова, С.А. Амелькин, A.M. Цирлин. - М.: Химия, 2000. 384 с.

3. Цирлин A.M. Методы оптимизации в необратимой термодинамике и макроэкономике. -М.: ФИЗМАТЛИТ, 2003. - 416 с.

4. Техническая термодинамика. Учебник для вузов / под ред. В.И.Крутова. М.: Высшая школа. 1981. - 438 с.

5. В.К. Андреев, Ю.А. Гапоненко, О.Н.Гончарова, В.В. Пухначев. Современные математические модели конвекции. - М.: Физматлит, 2008. - 368 с.

6. Гетлинг A.B. Конвекция Рэлея-Бенара. Структуры и динамика. - М.: Эдиториал. УРСС, 1999. - 248 с.

7. Кротов В.Ф., Гурман В. И. Методы и задачи оптимального управления. - М.: Наука, 1973. -448 с.

8. Энергосберегающий алгоритм оптимального управления температурой брожения пива (термодинамический подход). / А.Ю. Артюшкин, В.И. Карпов, A.B. Татаринов // Известия вузов. Пищевая технология. - 2010. - № 4. - с. 103-106.

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

Веткин A.B., д.т.н. проф. Сурис А.Л.

Университет машиностроения 8(499) 267-12-10, avetkin@mail.ru

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

Ключевые слова: трубчатая печь, состав топлива, коэффициент ослабления луча, степень черноты, радиационный теплоперенос. В последнее время в трубчатых печах нефтеперерабатывающих заводов во многих случаях используется газообразное топливо, состав которого существенно отличается от состава природного газа. Это связано с необходимостью использования отходящих горючих газов от различных установок. Эти газы могут содержать большое количество тяжёлых алканов, непредельные углеводороды, азот. Кроме того, в некоторых случаях топливный газ трубчатых печей может содержать очень высокую концентрацию водорода (иногда - выше 90 %). В этой связи представляет интерес исследовать влияние изменения состава топлива на характеристики печей, спроектированных для работы на природном газе.

В настоящей работе представлены результаты исследований воздействия состава топ-