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

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

КУРЧЕВ Андрей Олегович

МОДЕЛИРОВАНИЕ ТЕПЛОВЫХ ПРОЦЕССОВ В РЕГЕНЕРАТИВНЫХ УТИЛИЗАТОРАХ ТЕПЛОТЫ С ФАЗОВЫМИ ПЕРЕХОДАМИ

В НАСАДКЕ

Автореферат

диссертации на соискание ученой степени кандидата технических наук по специальности 05.14.04 - Промышленная теплоэнергетика

1 О ИЮН 2010

Иваново 2010

004603505

Работа выполнена в Государственном образовательном учреждении высшего профессионального образования «Ивановский государственный архитектурно-строительный университет».

Научный руководитель: доктор технических наук, профессор

ЕЛИН Николай Николаевич

Официальные оппоненты: доктор технических наук, профессор

АВДЮНИН Евгений Геннадьевич кандидат технических наук, доцент КУЛАГИН Станислав Михайлович

Ведущая организация: ЗАО «ИВЭНЕГОСЕРВИС»

Защита состоится 25 июня 2010 г. в // часов на заседании совета но защите докторских и кандидатских диссертаций Д 212.064.01 при ГОУ ВПО «Ивановский государственный энергетический университет имени В.И. Ленина» по адресу: 153003. г. Иваново, ул. Рабфаковская 34, ауд. Б-237.

Отзывы в двух экземплярах, заверенные печатью организации, просим направлять по адресу: 153003, г. Иваново, ул. Рабфаковская 34. Ученый совет ИГЭУ.

Тел. (4932) 38-57-12, факс (4932) 38-57-01.

E-mail: uch_sovet@ispu.ru

С диссертацией можно ознакомиться в библиотеке ГОУ ВПО «Ивановский государственный энергетический университет имени В.И. Ленина» (ИГЭУ). Автореферат диссертации размещен на сайте ИГЭУ www.ispu.ru

Автореферат разослан «¿0» мая 2010 г.

Ученый секретарь совета доктор технических наук, профессор

Мошкарин А.В.

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

Актуальность темы диссертации. Повышение тепловой эффективности производственных процессов - потребителей тепловой энергии, может быть достигнуто путем утилизации их вторичных энергоресурсов (ВЭР). В промышленности строительных материалов потери теплоты иногда достигают 40...50% от подводимой теплоты, а общий потенциал ВЭР оценивается в несколько миллионов тонн условного топлива. Другим примером являются компрессорные станции систем магистрального транспорта газа, оборудованные газотурбинными приводами, на которых имеется большое количество среднепотенциаль-ного тепла уходящих газов.

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

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

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

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

Разработка таких моделей является актуальной научной и технологической задачей, что и определило цель настоящей работы, которая выполнялась в рамках ФЦП «Интеграция» (2.1 - AI 18 Математическое моделирование ресурсосберегающих и экологически безопасных технологий) и планов НИР ИГАСУ, а также частично в рамках международного договора о научно-техническом сотрудничестве между кафедрой прикладной математики ИГЭУ и Горным институтом г. Алби, Франция.

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

Научная новизна результатов работы заключается в следующем:

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

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

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

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

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

Практическая ценность результатов состоит в следующем:

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

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

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

4. Разработанные математические модели, инженерные методы расчета и оценки эффективности, а также средства компьютерной поддержки моделирования и расчета нашли практическое применение в практике исследовательских и проектных работ в НТЦ «Промышленная энергетика».

Автор защищает:

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

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

3. Методику расчета теплофизических параметров парожидкостной смеси, заключенной в обогреваемой/охлаждаемой ячейке постоянного объема.

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

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

Основные положения диссертации были доложены, обсуждены и получили одобрение на XXI и XXII Международной научной конференции «Теоретические основы создания, оптимизации и управления энерго- и ресурсосберегающими процессами и оборудованием», Саратов - 2008; Псков - 2009; XIII и XV Международной конференции «Информационная среда вуза», Иваново, ИГАСУ, 2006, 2008, а также на научных семинарах кафедры гидравлики, водоснабжения и водоотведения ИГАСУ и кафедры прикладной математики ИГЭУ (2006-2010 гг.).

Публикации. Основное содержание диссертации опубликовано в 8 печатных работах, в том числе в 2-х изданиях, предусмотренных перечнем ВАК.

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

Содержание работы

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

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

Физическая сущность рассматриваемого теплового процесса состоит в теплообмене массивного протяженного тела с обтекающим его потоком газа. Известные интегральные методы его расчета, в которых вместо локальных распределенных значений температур теплоносителя и насадки используются значения, осредненные по длине канала и по времени цикла нагрева или охлаждения насадки, не могут обеспечить требуемую точность расчетных прогнозов. Тепловые потоки между газом и насадкой, рассчитываемые по этим температурам, оказываются не адекватными реальным. Более точное решение задачи могло бы быть выполнено на основе дифференциальных уравнений процесса теплопереноса в газе и насадке, но в этом случае модель сводится в двум сопряженным через граничные условия дифференциальным уравнениям, как минимум одно из которых является уравнением в частных производных с источ-никовым членом, описывающим фазовый переход. Это исключает возможность их аналитического решения, а численные методы трудно адаптируются к широкому спектру разнообразных условий процесса. В последнее время в работах H. Berthiaux, В.Е. Мизонова, C.B. Федосова и ряда других авторов для решения подобных задач были эффективно использованы ячеечные модели и связанный с ними математический аппарат теории цепей Маркова. Несомненным достоин-

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

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

В заключение главы сформулированы детализированные задачи исследования.

Вторая глава посвящена разработке ячеечной модели прогрева насадки с фазовым переходом типа «плавление» обтекающим ее потоком горячего газа - первой стадии процесса работы насадки в регенеративном цикле. В основу модели положено описание процесса в элементарной ячейке насадки, все тепловые параметры которой могут считаться равномерно распределенными по ее объему. Расчетная схема процесса в ячейке показана на рис. 1. Будем рассматривать теплофизическое состояние процесса через малые промежутки времени Ат. Тогда текущее время процесса может быть выражено как тк=(к-1)Дт, где к=1,2,... - номер перехода из одного состояния в другое, то есть к может рассматриваться как дискретный целочисленный аналог времени. Изменение температуры ячейки при подводе к ней теплоты при отсутствии фазовых переходов описывается рекуррентным равенством

Рис.1. Структура тепловых потоков в ячейке (а) и среды в ней (б).

1к+1=1к+а,(е-0-а2(1к-и,), (!)

а.Б.Дт аЯ,Ат

где а, , (2) а2=^—, (3)

срУ срУ

- приведенные безразмерные коэффициенты теплоотдачи (индекс к при

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

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

При достижении в ячейке температуры плавления материала 1те вступают в силу другие соотношения, описание которых сопровождается логическими операторами. Если происходит нагрев ячейки, то есть в соответствии с (1) по-

1с +1 1с и+|

лучается г > I, но если г превосходит температуру плавления 1те, то эта температура принимается равной температуре плавления 1к+'= ^ причем это справедливо до тех пор, пока т5к>0, то есть пока в ячейке присутствует твердая фаза. При этом вся разница подведенной от источника и отданной среде теплоты идет на исчезновение твердой (прирост жидкой) фазы, масса которой в конце временного перехода составит

ш8к+1 = - (дс^ - ддошк)/г„ (4)

где г - удельная теплота плавления, Д(58к и Д0О1|,к - количества теплоты, переданной в течение перехода от источника и в окружающую среду, рассчитываемые по соответствующим разностям температур.

Прогрев материала до температуры плавления рассчитывается по (1) с теплоемкостью и плотностью твердой фазы (с=с5, р=р5, а|=а]5, а2=а28). Как только плавление полностью заканчивается (т5к+|=0), начинается прогрев однородной жидкой фазы, который снова рассчитывается по (1) с с=с|, р=рь а1=ац, а2=а;|. В процессе плавления, когда в ячейке присутствуют и твердая, и жидкая фазы, для точного расчета необходимо использовать теплоемкость и плотность (а далее - и теплопроводность) их смеси. Чтобы не загромождать этими соотношениями расчетные зависимости, приближенно будем считать, что эти свойства одинаковы в твердом и жидком состоянии, а, следовательно, и в состоянии смеси.

к+1 к

Если происходит охлаждение ячейки, то есть I < I, и оказывается, что 1к+|<1те, то принимается, что 1к+'= ^ а разница подведенной от источника и отданной среде теплоты идет на формирование твердой (убыть жидкой) фазы.

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

Расчетная схема процесса теплового взаимодействия насадки с обтекаю-

к ' \ '

ш8 \

А, , /

Рис.2. Влияние переменной температуры теплового источника (жирная линия) на кинетику прогрева и формирования расплава в ячейке (а|=0,5,1=1, 32=0).

I "¡> ¡(Ш)1-

'' И '

а)

'$»^80111

Ах

Н

Рис.3. Расчетная схема процесса (а) и ее ячеечное представление (б)

щим ее газом показана на рис.3. Из регенератора выделен элемент, включающий в себя половину ширины канала для прохода газа и половину ширины насадки. В силу симметрии считается, что поперечные тепловые потоки за пределы этого элемента отсутствуют. Длина элемента разбита на ш элементарных ячеек длиной Ах = Ь/т. В к-й момент времени тепловое состояние газа или насадки во всей цепи может быть представлено векторами-столбцами

[V 1

1к 182

о"' 1к' . 6т.

и т.д.

(5)

При переходе от одного состояния к другому с векторами состояния происходят следующие изменения:

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

АО^а^/Дх.Ч^-О.'Дг + а^ЧЛх.Ч^-^Лг, (6)

где первое слагаемое соответствует конвективному, а второе - радиационному теплообмену;

- перенос массы газа вдоль ячеек и пополнение первой ячейки за счет внешней подачи газа (вектор источников массы т,,/*)

твьи=М1к*(т,к+тЕГк); (7)

- пополнение запаса теплоты в первой ячейке за счет внешней подачи газа (вектор источников теплоты <38гк), ее потеря в ячейках за счет передачи теплоты в насадке и перенос вдоль цепи ячеек (вектор стоков теплоты А(}к)

0К"+'=М8к*(08Чс>?гк-Л0к), (8)

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

мв =

у61 о

о

о О

о о

1 — v

83

О

о о о

1 - у„

(9)

8(т-1) ' '8т_

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

V =у =

81 61 Ау та '

(Ю)

где У8| - размерная скорость газа в ячейке, рассчитываемая по уравнению неразрывности с учетом уравнения состояния, так как плотность газа меняется по длине канала;

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

к+1_

д^ЧсЛ'тЛ,

(Ц)

где сг - вектор теплоемкостей газа в ячейках; в приведенных выше формулах операторы .* и ./ означают поэлементное умножение и деление векторов;

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

<5гк+,=Рок*(Огк+А<Л,

1йк+1= <Злк+1./(сс.* те), если т,,к=0 или тик=1,

(12) (13а)

-и, если ^ и 0< тц<1; т,; - т,," - АО /г. ¡=1,...,т, (136)

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

Ро =

1-с1 (1 0 0 0

с! 1-2<1 (1 0 0

0 д 0 0

0 0 0 0 й

0 0 0 0 й 1-е!

(14)

где

, Ат Ах ....

а = —= аг —, (15)

сгрг Ау" Ау

- доля теплоты, переносимая путем теплопроводности в соседние ячейки из данной за один временной переход (А^— коэффициент теплопроводности материала насадки, аг- коэффициент ее температуропроводности).

Равенства (6)-(15) вместе с уравнением состояния газа являются полной ячеечной моделью конвективно-радиационного теплообмена газа с насадкой с фазовым переходом типа плавление.

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

У8=0,2 у8=0,5 У8=0,8

Рис.4. Кинетика распределения температур газа и насадки и формирования жидкой фазы в насадке при различных скоростях движения газа (г=2; а2=0,1; ^=0,5)

За базовую кинетику принято значение у=0,5, которому соответствует средний столбец графиков. За 200 временных переходов достигается полное плавление материала насадки. При у=0,2 (то есть в 2,5 раза меньшем расходе газа) за то же время достигается полное плавление материала лишь в первых по

ходу газа 4-х ячейках из 11-ти, а в последних ячейках плавление даже не начинается. Наоборот, при v=0,8 процесс протекает гораздо быстрее, после 150-ти переходов насадка практически полностью расплавлена и прогрета до температуры газа. При этом газ перестает передавать теплоту насадке, и регенеративный процесс теряет смысл.

В третьей главе рассмотрены вопросы моделирования циклов нагрева и охлаждения насадки, базирующиеся на разработанной выше модели. При общей продолжительности цикла N временных переходов на нагрев отводится Ng переходов, на охлаждение - Na (здесь и далее индекс g относится к горячему газу, а индекс а - к нагреваемому воздуху) Соотношение Ng и Na определяется соотношением объемов (масс) горячего газа, который необходимо охладить, и холодного воздуха, который необходимо нагреть. Охлаждение насадки нагреваемым воздухом может быть описано в точности той же моделью, что и ее нагрев горячим газом. Естественно, что параметры газа, его скорость и коэффициенты теплоотдачи будут разными на стадии нагрева и охлаждения. Единственное принципиальное различие может возникнуть, если горячий газ и нагреваемый воздух движутся противотоком, что практически всегда имеет место в регенераторах с твердой насадкой. В этом случае элементы переходной матрицы для газа перемещаются из-под главной диагонали наверх, а ненулевые элементы вектора подачи газа относятся к последней ячейке цепи.

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

AQz = Z AQk = sum(AQk), (16)

i

где оператор sum означает суммирование всех элементов вектора. Считается, что AQx положительно при нагреве насадки и отрицательно при ее охлаждении.

Полная теплота, запасенная насадкой за к переходов, рассчитывается как

QL(k) = ¿AQk = cumsum(sum( AQk)), (17)

Ы1

где оператор cumsum означает кумулятивное (накопительное) суммирование.

Очевидно, что при установившемся режиме количество теплоты в насадке в конце каждого цикла «нагрев-охлаждение» такое же, как в его начале.

По аналогичным соотношениям рассчитывается полное содержание жидкой фазы в частично расплавленной насадке

m

ma00 = Xm!i/ m = sum(mll')/m. (18)

i=I

Если расчет начинается с холодной насадки, то ее тепловое состояние в

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

На рис.5 показано влияние продолжительности цикла нагрева-охлаждения на изменение температуры газа и воздуха на выходе и эволюцию распределения параметров газа и насадки по длине регенератора, когда продолжительность нагрева и охлаждения одинаковы. В случае (а) с продолжительностью полного цикла N=250 температуры газа и воздуха на выходе остаются практически постоянными, а температура конечной ячейки насадки остается постоянной за весь период цикла. Из графика распределений для этого случая следует, что постоянна не только температура последней ячейки, но и всех ячеек насадки, то есть во всех ячейках материал находится в частично расплавленном состоянии, а не в полностью жидком или полностью твердом. Очевидно, что в этом случае полностью реализуются все преимущества насадки с фазовым переходом, обеспечивающей практически постоянную температуру уходящего газа при нагреве насадки и постоянную температуру подогреваемого воздуха пои ее охлаждении.

10 « 20 к/20

Рис.5. Изменение параметров газа и насадки в цикле нагрева и охлаждения насадки в установившемся цикле при различной его продолжительности (а|=0,1: а2=0,0!; у=0,6; <3=0,2; г=0,5; 1^=0,5)

В случае (б) при N=500 температура конечной ячейки насадки также остается постоянной, но у температуры газа имеется небольшое повышение в конце цикла нагрева, а у температуры воздуха - понижение в конце цикла охлаждения. Это объясняется тем, что в начальной зоне насадки в конце каждого из циклов фазовый переход полностью завершается и происходит прогрев или охлаждение индивидуальной фазы. Эти эффекты еще более выражены в случае (в), когда N=1000. Здесь зона частичного расплава существует лишь в части цикла нагрева и охлаждения, и температура насадки и газа меняется значительно в течение цикла. Таким образом, очевидно, что с точки зрения поддержания стабильной температуры обоих газов преимущественным является установившийся режим, при котором насадка всегда находится в двухфазном состоянии, когда ее температура постоянна. При этом рациональная продолжительность цикла зависит от удельной теплоты плавления, поскольку, чем она больше, тем дольше при прочих равных условиях длится фазовый переход.

На рис.6 показан пример кинетики накопления теплоты насадкой без и с

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

А\

/ у ■ \ \

/ у /7 \ \ \ \ \

]А

Рис.6. Накопление и отдача насадкой теплоты при наличии (жирн.) и отсутствии (тонк.) фазового перехода. (в]=0,2; а2Ю,01;с1=0,2; у=0,9; г=1; и=0,5)

мерной моделью.

Если при насадке в виде твердого тела сК</с1()т =1/с, где с - теплоемкость материала, то при насадке с рассмотренным выше фазовым переходом «плавление-отвердевание», где цикл нагрева и охлаждения организован так, что он не завершается полностью, температура насадки практически постоянна и равна температуре плавления, а <к|/сК)£ = 0. С точки зрения поддержания высокого температурного напора этот случай наиболее благоприятный. Однако далеко не всегда удается подобрать материал или сплав с подходящей температурой плавления, а цикл нагрева и охлаждения не всегда легко организовать без полного завершения фазовых переходов. Другим вариантом насадки с фазовым переходом является насадка с кипением жидкости в постоянном объеме насадки.

Здесь ск^рЕ > 0, но температура с ростом возрастет медленнее, чем в твердотельной насадке, то есть ее тепловые свойства занимают промежуточное положение между твердотельной насадкой и насадкой с фазовым переходом в виде плавления. Преимуществом такой насадки являются большие возможности регулирования ее характеристик, например, путем изменения давления в ее объеме. Однако для расчета процесса с такой насадкой необходимо знать количественно характеристику ^¡/¿С^.

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

Пример результатов моделирования показан на рис.7,8. Рис.7 иллюстрирует установление цикла нагрева-охлаждения при различных соотношениях Б-Ы^^+Ма) продолжительности нагрева к общей продолжительности цикла. На рис.8 показано изменение температуры газа и воздуха в течение цикла, которое менее стабильно, чем при насадке с плавлением, но более стабильно, чем при твердой насадке.

Рис.7. Формирование стационарного цикла Рис.8. Изменение температуры го-

работы насадки по X] (У8о=1м/с, у=500°С, рячего и холодного газа на выходе

4,0=0,5м/с, 1„и=20 "С, ае-5ит/м град, при различной длине регенератора 1 - 5=0,4; 2 - 0,45; 3 - 0,5). (1 _ 1=10м, 2 - 5м, 3 - 2,5м).

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

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

Проверка адекватности математической модели производилась:

1. Для твердотельной насадки путем сравнения результатов расчетов температур продуктов сгорания, воздуха и насадки в выходном сечении с результатами теплотехнического обследования регенеративного подогревателя воздуха вагранки СМ-5232, выполненного из шамотного кирпича.

2. Для насадки с фазовым переходом путем сравнения результатов расчетов времени, необходимого для осуществления фазового перехода в аккумуляторе фазового перехода (АФТ), используемого в системах кондиционирования и вентиляции, заполненного этилен-карбонатом с температурой 36,4 °С, с фактическим временем.

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

Задача конструктивного расчета регенератора заключается в определении его геометрических размеров: длины проточной части, в которой размещается насадка, и размеров каналов для прохода теплоносителей. Исходными данными при этом являются расходы теплоносителей, их теплофизические свойства, температуры на входе в регенератор и желаемые температуры на выходе из него, теплофизические свойства материала насадки, а также располагаемый перепад давления, который можно использовать для создания потока теплоносителей через проточную часть регенератора. Используя зависимости аэродинамического сопротивления насадки от эквивалентного диаметра проходных сечений её каналов, можно получить взаимосвязь приведенной скорости теплоносителя \уо и предельно допустимой длины проточной части насадки Ь. Для вывода зависимости длины проточной части насадки от эквивалентного диаметра проходных сечений её каналов для заданных параметров потоков теплоносителей и располагаемого перепада давления, используется уравнение теплового баланса, в котором коэффициент теплоотдачи от теплоносителя к насадке ач отнесен к объему канала, в котором расположена насадка. Для этого величину «обычного» коэффициента теплоотдачи а, отнесенного к площади поверхности насадки, омываемой теплоносителем, следует умножить на величину удельной поверхности насадки (отношение площади поверхности теплообмена между газом и насадкой к её объему). Выполненные численные эксперименты показали, что при увеличении эквивалентного диаметра проходных сечений каналов насадки величина а быстро уменьшается, достигая минимального значения при ~ 3 см, а затем увеличивается за счет увеличения «радиационной» составляющей, а величина ау монотонно уменьшается.

Для расчета взаимосвязи между эквивалентным диаметром проходных сечений каналов насадки и длиной её проточной части Ь для заданных пара-

метров потока газа Др, с8, рц, 181П, 1£0Ш и температуры фазового перехода материала насадки ^ получено уравнение

сер,

(19)

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

- увеличение температуры горячего (уменьшение температуры холодного) теплоносителя на входе в регенератор;

- уменьшение температуры горячего (увеличение температуры холодного) теплоносителя на выходе из регенератора;

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

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

С целью оценки эффективности применения насадки с фазовым переходом вместо твердотельной были выполнены расчеты для системы регенеративного подогрева воздуха вагранки СМ-5232 для производства стекловаты, при расходе воздуха 12000 м3/час. Использование насадки с фазовым переходом позволяет обеспечить стабильные температуры теплоносителей на выходе из регенератора в течение всего цикла. Динамика изменения этих температур показана на рис. 9.

80 т/т..,%

......................Г...... ... ■> — .....-

\ Ч б).

N

ч

V " - !' V ^ ^

** «ч

6» 80 т/т.«,

Рис.9. Результаты расчета динамики температур дымовых газов (а) и воздуха (б) на выходе из регенератора: 1 - для твердотельной насадки, 2 - для насадки с фазовым переходом типа плавления.

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

Для рассматриваемого объекта экономия тепловой энергии составит около 5,800 ГДж/час или 1,384 Гкал/час или 4,76%, а требуемая поверхность теплообмена уменьшится в в 2,4 - 3,13 раза.

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

Разработанная математическая модель и ее программно-алгоритмическое обеспечение внедрены в НТЦ «Промышленная энергетика», где также приняты к внедрению рекомендации по расчету тепловой и экономической эффективности регенеративного подогрева направляемого на сгорание воздуха.

Основные результаты диссертации

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

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

3. Показано, что наиболее эффективным является применение насадки с фазовым переходом «плавление-отвердевание», а наиболее рациональным циклом ее нагрева и охлаждения является цикл, когда она в течение всего цикла пребывает в частично расплавленном состоянии.

4. Предложена методика расчета и получены эмпирические зависимости изменения температуры и фазового состояния парожидкостной смеси в замкнутой обогреваемой/охлаждаемой ячейке насадки.

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

6. Разработанная математическая модель и ее программно-алгоритмическое обеспечение внедрены в НТЦ «Промышленная энергетика», где также приняты к внедрению рекомендации по расчету тепловой и экономической эффективности регенеративного подогрева направляемого на сгорание воздуха.

Основное содержание диссертации опубликовано в следующих печатных работах автора

В изданиях, рекомендованных ВАК:

1. Курчев А.О. Математическое моделирование теплового процесса в регенераторе с насадкой с фазовыми переходами [Текст] / Курчев А.О., Елин H.H., Мизонов В.Е., Субботин В.И. // Промышленная энергетика. №12, 2008, с.33-36.

2. Курчев А.О. Моделирование регенеративного теплообмена с фазовыми переходами в насадке [Текст] / Курчев А.О., Елин H.H., Мизонов В.Е. // Вестник ИГЭУ, вып.З, 2008, с.44-46.

В других изданиях:

3. Мизонов В.Е. Математическая модель теплового процесса в регенераторе с насадкой, допускающей фазовые переходы [Текст] / Мизонов В.Е., Курчев А.О., Елин H.H., Лезнов B.C. // Труды XXI МНК «Математические методы в технике и технологиях - ММТТ21», т.5, Саратов, 2008, с.118-119.

4. Мизонов В.Е. Сравнение тепловых характеристик насадок регене-рационных теплообменников [Текст] / В.Е. Мизонов, А. О. Курчев, H.H. Елин, В.Б. Медведев // Информационная среда вуза: Материалы XV Междунар. конф. Иваново, ИГАСУ. 2008, с.782-783.

5. Мизонов В.Е. Математическое моделирование тепломассообмена в контактных утилизаторах теплоты продуктов сгорания газового топлива [Текст] / В.Е. Мизонов, А.О. Курчев, H.H. Елин.// Информационная среда вуза: Материалы XV Междунар. конф. Иваново, ИГАСУ. 2008, с.779-781.

6. Медведев, В.Б. Тепловые характеристики различных типов насадки в регенеративных теплообменниках [Текст]/ В.Б.Медведев, А.О. Курчев, В.А. Огурцов, H.H. Елин.// Ученые записки инженерно-строительного факультета / ИГАСУ - Иваново, 2008. - Вып. 4. - с.226-231

7. Мизонов В.Е. Конструктивный расчет проточной части регенеративных теплообменников [Текст] / В.Е. Мизонов, А.О. Курчев, H.H. Елин.// Труды XXII МНК "Математические методы в технике и технологиях -ММТТ22", т.9, Псков, 2009, с.99-100.

8. Курчев, А.О. Моделирование процессов с управляемой газотурбинной энергетической установкой [Текст]/ А.О. Курчев, В.М. Ткачев.// Информационная среда вуза: Материалы XIII Междунар. конф. Иваново, ИГАСУ. -2006, с.577-581.

КУРЧЕВ Андрей Олегович

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

АВТОРЕФЕРАТ диссертации на сонекяние ученой степени доктора технических наук

Подписано в печать 18.05.2010. Формат 60x84 1/16. Печать плоская. Усл. печ. л. 1,39 Тираж 100 экз. Заказ № 119. ГОУ ВПО «Ивановский государственный энергетический университет им. В.И. Ленина» 153003, Иваново, ул. Рабфаковская, 34.

Отпечатано в УИУНЛ ИГЭУ.

Условные обозначения.

Введение.

Глава 1. Использование регенеративного подогрева для утилизации теплоты продуктов сгорания теплоиспользующих установок.

1.1. Анализ существующей практики утилизации вторичных энергетических ресурсов теплоиспользующих установок и основных направлений их использования.

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

1.3. Оценка эффективности эксплуатации регенеративных теплообменных аппаратов.

1.4. Методы расчета регенеративных теплообменников.

1.5. Математическое моделирование взаимодействия газа с насадкой. Ячеечные модели прогрева одно- и двухмерных объектов.

1.6. Постановка задачи исследования.

Глава 2. Моделирование тепловых процессов в обтекаемой газом насадке с фазовыми переходами.

2.1. Описание кинетики нагрева-охлаждения и плавления-отвердевания в элементарной ячейке.

2.2. Уровни декомпозиции моделируемого объекта и его расчетные схемы.

2.3. Описание кинетики нагрева и плавления в одномерном объекте.

2.4. Модель сопряженного теплообмена газа с насадкой.

2.5. Выводы по главе 2.

Глава 3. Моделирование тепловых циклов нагрева и охлаждения насадки в регенераторе.

3.1. Цикл нагрева и охлаждения и установление его параметров.

3.2. Расчетное исследование циклов работы регенератора с плавящейся насадкой.

3.3. Процессы и циклы в насадке с фазовым переходом «кипение-конденсация».

3.3.1. Моделирование состояния парожидкостной смеси в насадке.

3.3.2. Расчетное исследование циклов работы насадки с фазовым переходом «кипение-конденсация».

3.4. Выводы по главе 3.

Глава 4. Инженерная методика расчета регенеративных теплообменников с насадками с фазовыми переходами.

4.1. Тепло физические свойства заполнителей замкнутых ячеек насадок регенераторов.

4.1.1. Особенности теплофизических свойств жидких металлов и сплавов.

4.1.2. Особенности теплофизических свойств парожидкостных смесей.

4.2. Расчет теплообмена в регенеративном подогревателе.

4.3. Расчет геометрических параметров насадки.

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

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

4.6. Выводы по главе 4.

Основные результаты диссертации.

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

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

Одним из возможных направлений увеличения глубины утилизации ВЭР при приемлемых габаритах регенератора является использование насадки, внутри которой могут происходить фазовые переходы, имеющие значительную удельную теплоту. Это позволяет увеличить полную теплоемкость насадки и тепловые потоки между газом и насадкой и насадкой и воздухом, поскольку при фазовом переходе изменение ее температуры значительно меньше, чем при обычном прогреве, и сохраняется высокая разность температур, обеспечивающая повышенные тепловые потоки. Кроме того, более стабильная во времени температура насадки позволяет иметь более стабильные температуры газа и воздуха на выходе из регенератора и обойтись без используемых для этого регенераторов с вращающейся насадкой типа «Юнгстрем», имеющих значительную стоимость и габариты. Вместе с тем, преимущества регенераторов с насадкой с фазовыми переходами могут быть реализованы только при соблюдении надлежащих температурных режимов и конструктивных параметров регенератора. Задачи их проектирования не могут быть решены с помощью известных методов расчета, в которых вместо актуальных локальных значений температур теплоносителей используются значения, осреднен-ные по длине канала и времени цикла, а также ряд других далеко идущих допущений.

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

Разработка таких моделей является актуальной научной и технологической задачей, что и определило цель настоящей работы, которая выполнялась в рамках ФЦП «Интеграция» (2.1 - All8 Математическое моделирование ресурсосберегающих и экологически безопасных технологий) и планов НИР ИГ АСУ, а также частично в рамках международного договора о научно-техническом сотрудничестве между кафедрой прикладной математики ИГЭУ и Горным институтом г. Алби, Франция.

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

Научная новизна результатов работы заключается в следующем:

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

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

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

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

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

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

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

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

4. Разработанные математические модели, инженерные методы расчета и оценки эффективности, а также средства компьютерной поддержки моделирования и расчета нашли практическое применение в практике исследовательских и проектных работ в ЗАО «ИВЭНЕРГОСЕРВИС».

Автор защищает:

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

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

3. Методику расчета теплофизических параметров парожидкостной смеси, заключенной в обогреваемой/охлаждаемой ячейке постоянного объема.

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

Основные положения диссертации были доложены, обсуждены и получили одобрение на Международной научной конференции «Теоретические основы создания, оптимизации и управления энерго- и ресурсосберегающими процессами и оборудованием», Иваново, ИГХТУ, 2007; XIV Международной конференции «Информационная среда вуза», Иваново, ИГ АСУ, 2007, а также на научных семинарах кафедры гидравлики, водоснабжения и водоотведе-ния ИГАСУ и кафедры прикладной математики ИГЭУ (2006-2008гг.).

Публикации. Основное содержание диссертации опубликовано в 7 печатных работах, в том числе в 2-х изданиях, предусмотренных перечнем ВАК.

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

Основные результаты диссертации

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

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

3. Показано, что наиболее эффективным является применение насадки с фазовым переходом «плавление-отвердевание», а наиболее рациональным циклом ее нагрева и охлаждения является цикл, когда она в течение всего цикла пребывает в частично расплавленном состоянии.

4. Предложена методика расчета и получены эмпирические зависимости изменения температуры и фазового состояния парожидкостной смеси в замкнутой обогреваемой/охлаждаемой ячейке насадки.

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

6. Разработанная математическая модель и ее программно-алгоритмическое обеспечение внедрены в НТЦ «Промышленная энергетика», где также приняты к внедрению рекомендации по расчету тепловой и экономической эффективности регенеративного подогрева направляемого на сгорание воздуха.

1. Андрющенко А.И. Основы термодинамики циклов теплоэнергетических установок / А. И. Андрющенко. М.: Высшая школа, 1977. - 280 с.

2. Бакластов, A.M. Проектирование, монтаж и эксплуатация теп-ломассообменных установок/ A.M. Бакластов, В.А. Горбенко, П.Г. Удыма. -М.: Энергоиздат, 1981. 336 с.

3. Бажан, П.И. Справочник по теплообменным аппаратам / П. И. Ба-жан, Г. М. Каневец, В. М. Селиверстов. М.: Машиностроение, 1989. - 366 с.

4. Баруча-Рид, А. Т. Элементы теории Марковских процессов и их приложения / А. Т. Баруча-Рид. М.: Наука, 1969. - 511 с.

5. Боришанский, В.М. Жидкометаллические теплоносители./В.М. Боришанский, С.С. Кутателадзе, Новиков И.И. и др. — М.: Атомиздат, 1967. 442 с.

6. Бродянский, В. М. Эксергетический метод термодинамического анализа / В.М. Бродянский; М.: Энергия, 1973. - 217 с.

7. Варгафтик, Н.Б. Теплофизические свойства веществ: справоч-ник/Н.Б. Варгафтик М.: Машиностроение, 1972.- 720 с.

8. Высокотемпературные теплотехнические процессы и установки / И.И. Перелетов, JT.A. Бровкин, Ю.И. Розенгарт и др., Под ред. А.Д. Ключникова. — М.: Энергоатомиздат, 1989.- 336 с.

9. Дацковский, В.М. О расчете вращающегося регенератора/В.М. Дацковский В.М. // Теплоэнергетика. 1965. - №8. - с. 93-95.

10. Жуков В. П. Применение теории цепей Маркова к динамическому моделированию теплообменных аппаратов / В. П. Жуков, Е. В. Ба-рочкин, В. Е. Мизонов, Г. В. Ледуховский // Изв.ВУЗов, Химия и химическая технология. 2005. - Т. 48. - Вып. 4. - С. 87-89.

11. Исаченко, В. П. Теплопередача: учебник для вузов / В. П. Исаченко, В. А. Осипова, А. С. Сукомел; 4-е изд., перераб. и доп. - М.: Энерго-атомиздат, 1981. - 416 с.

12. Ионин А.А. Газоснабжение: учебник для вузов /А.А. Ионин М.: Стройиздат, 1989. -439 с.

13. Кафаров, В. В. Оптимизация тепломассообменных процессов и систем / В. В. Кафаров, В. П. Мешалкин, Л. В. Гурьева. М.: Энергоатомиз-дат, 1988.- 192 с.

14. Кейс, В.М. Компактные теплообменники/В.М. Кейс, А.Л. Лондон М.: Энерия, 1967. - 472 с.

15. Келоглу Ю.П. Краткий металловедческий справочник./ Ю.П. Ке-логлу, К.М. Захариевич, М.И. Карташевская — Кишинев: Картя Молдовеняска, 1969.-167 с.

16. Кемени, Дж. Счетные цепи Маркова: Пер. с англ. / Дж. Кемени, Дж. Снелл, А. Кнепп. М.: Наука, Гл. ред. физ.-мат. лит., 1987. - 416 с.

17. Кириллов, П.Л. Справочник по теплогидравлическим расче-там/П.Л. Кириллов, Ю.С. Юрьев, В.П. Бобков М.: Энергоатомиздат, 1984. -296 с.

18. Кирсанов, Ю.А. Двухмерная теплопроводность в твердом теле при циклических четырехпериодных граничных условиях третьего рода/ Ю.А. Кирсанов// Известия АН СССР Энергетика, 1966.- №2,- с. 69-74.

19. Кирсанов, Ю.А. Математическое моделирование тепловых процессов в регенеративном воздухоподогревателе/ Ю.А. Кирсанов // Теплоэнергетика. 1999. № 1. - с. 51-54.

20. Кирсанов, Ю.А. Выбор уравнения для замыкания сопряженной задачи циклического теплообмена твердого тела с холодным и горячим теплоносителями/ Ю.А. Кирсанов //Изв. Вузов. Проблемы энергетики. 2003. -№1-2.

21. Кирсанов, Ю.А. Тепловой расчет регенеративного воздухоподогревателя/ Ю.А. Кирсанов // Изв. Вузов. Авиац. Техника. 1999. - №1. -с.32-35.

22. Кирсанов, Ю.А. Влияние нестационарности на теплоотдачу в регенеративном воздухоподогревателе/ Ю.А. Кирсанов //Изв. Вузов. Авиац. Техника. 2003. - №1. - с.22-25.

23. Кирсанов, Ю.А. Циклический сопряженный теплообмен потоков теплоносителей с твердым телом/ Ю.А. Кирсанов// Известия АН СССР -Энергетика, 1998.-№5.-с. 113-119.

24. Кошкин, В. В. Нестационарный теплообмен / В. В. Кошкин, Г. А. Дрейцер, С.А. Ярхо; М.: Машиностроение, 1973. - 347 с.

25. Краткий справочник физико-химических величин./под. ред. К.П. Мищенко, А.А. Равделя. М.: Химия, 1972. - 200 с.

26. Кулинченко В.Р. Справочник по теплообменным расчетам/В.Р. Кулинченко Киев: Техника, 1990. - 165 с.

27. Курчев, А.О. Моделирование процессов с управляемой газотурбинной энергетической установкой./В.М. Ткачев.// Сб. ст. XIII Междунар. конф. «Информационная среда вуза», Иваново, ИГАСУ. 2006, с.577-581.

28. Курчев, А.О. Моделирование регенеративного теплообмена с фазовыми переходами в насадке./ А.О. Курчев, Н.Н. Елин, В.Е. Мизонов. // Вестник ИГЭУ. 2008. - вып. 3. - с. 44-46.

29. Курчев, А.О. Сравнение тепловых характеристик насадок реге-нерационных теплообменников./ А. О. Курчев, В.Б. Медведев, Н.Н. Елин. // Сб. ст. XV Междунар. конф. «Информационная среда вуза», Иваново, ИГА-СУ. 2008. - с.782-783.

30. Курчев, А.О. Математическое моделирование теплового процесса в регенераторе с насадкой с фазовыми переходами./А.О. Курчев, Н.Н. Елин, В.Е. Мизонов, В.И. Субботин // Промышленная энергетика 2008. -№12.- с.33-36.

31. Ли, Ц. Оценивание параметров марковских моделей по агрегированным временным рядам / Ц. Ли, Д. Джадж, А. Зельнер. М.: Статистика, 1977.-355 с.

32. Лыков, А. В. Теория тепло- и массопереноса / А. В. Лыков, Ю. А. Михайлов. М. - Л., Госэнергоиздат, 1963. - 536 с.

33. Лыков, А. В. Тепломассообмен. Справочник. /А.В. Лыков М.: Энергия. 1972, 560 с.

34. Медведев, В.Б. Тепловые характеристики различных типов насадки в регенеративных теплообменниках./ В.Б.Медведев, А.О. Курчев, В.А. Огурцов, Н.Н. Елин.// Ученые записки инженерно-строительного факультета / ИГ АСУ Иваново, 2008. - Вып. 4. - с.226-231.

35. Новожилов, Ю.Н. Схема рециркуляции дымовых газов в котлах с вращающимися регенеративными воздухоподогревателями/Ю.Н. Новожилов // Промышленная энергетика, — 2002. № 7. — с. 26-30.

36. Сазанов Б.В. Теплоэнергетические системы промышленных предприятий. /Б.В. Сазанов, В.И. Ситас. М.: Энергоатомиздат, 1990.- 304 с.

37. Семененко, Н.А. Вторичные энергоресурсы промышленности и энерготехнологическое комбинирование/Н.А. Семененко М.: Энергия, 1968.-296 с.

38. Систер, Г. В. Принципы повышения эффективности тепломассо-обменных процессов / В. Г. Систер, Ю. В. Мартынов. — Калуга: Издательство Н. Бочкаревой. 1998. - 508 с. ISBN 5 - 89552 - 036 - 7.

39. Строительные машины. Справочник. В 2-х т. Под ред. В.А. Баумана. т.2. Оборудование для производства строительных материалов и изделий. М.: Машиностроение, 1977. - 496 с.

40. Теория тепломассообмена. Учеб. для вузов./ Под ред. А. И. Леонтьева. М., Высш. шк., 1979. 495 с.

41. Телегин, А.С. Теплотехника и нагревательные устройства/А.С. Телегин, В.Г. Авдеева М.: Машиностроение, 1985. - 247 с.

42. Тепло- и массообмен. Теплотехнический эксперимент: Справочник / Е.В. Аметистов, В. А. Григорьев, Б.Т. Емцев и др.; под общ. ред. В. А. Григорьева, В. М. Зорина. М.: Энергоиздат, 1982. - 512 с.

43. Теплотехническое оборудование и теплоснабжение промышленных предприятий: Учебник для техникумов / Голубков Б.Н., Данилов О.Н., Зосимовский J1.B. и др.; Под ред. Б.Н. Голубкова. 2-е изд., перераб. - М.: Энергия, 1979.-554 с.

44. Троянкин, Ю.В. Реконструкция регенераторов промышленных печей с целью экономии топлива / Ю.В. Троянкин// Промышленная энергетика. 2004. - № 5. - с. 22-23.

45. Уонг, X. Основные формулы и данные по теплообмену для инженеров. Справочник / X. Уонг. М.: Атомиздат, 1979. - 265 с.

46. Цой, П. В. Системные методы расчета краевых задач тепломассо-переноса / П. В. Цой; М.: Издательство МЭИ, 2005. 568 с.

47. Фан-Лянь-Цэнь. Дискретный принцип максимума / Фан-Лянь-Цэнь, Валь Чу-сен. М.: Мир, 1967. - 183 с.

48. Федосов, С. В. Моделирование тепловых процессов в регенеративных утилизаторах теплоты уходящих газов промышленных печей. / С. В. Федосов, В.Е. Мизонов, Н. Н. Елин, С. В. Хавер // Строительные материалы. -2007.-№9.-С. 14-16.

49. Фраас, А. Расчет и конструирование теплообменников/А. Фраас, М. Оцисик. Пер. с англ. - М.: Атомиздат. 1971. - 326 с.

50. Хаузен, X. Теплопередача при противотоке, прямотоке и перекрестном токе/Х. Хаузен М.: Энергоиздат, 1981.-384 с.

51. Хейлигенштедт, В. Теплотехнические расчеты/В. Хейлиген-штедт М.: ОНТИ, 1937. - 432 с.

52. Хоблер, Т. Теплопередача и теплообменники/Т. Хоблер, М. Тадеуш Пер. с польского А.В. Плиса. Под ред. П.Г.Романкова.- JL: Госхимиз-дат,, 1961.-820 с.

53. Хавер С.В. Моделирование прогрева и охлаждения насадки регенеративного теплообменника/С.В. Хавер, Н.Н. Елин, В.Е. Мизонов, А.Б. Иванов// Изв. вузов «Химия и хим. технология», т. 50, вып. 12, 2007, с. 105107.

54. Шак, А. Промышленная теплопередача/ А. Шак М.: Металлург-издат, 1961.-528с.

55. Шорин, С.Н. Теплопередача: Учебник для вузов./С.Н. Шорин -М.:Высш. школа, 1964. 490 с.

56. Эксергетический метод и его приложения / под ред. В. М. Бродян-ского. -М.: Мир, 1967. -172 с.

57. Юдаев, Б.Н. Теплопередача: Учебник для вузов./Б.Н. Юдаев -М.:Высш. школа, 1973. 360 с.

58. Application of the theory of markovian chains to simulation and analysis of processes with granular media / V. Mizonov, H. Berthiaux, V. Zhukov. // With contributions by K. Marikh, E. Barantseva, D. Ponomarev. Albi -2002.

59. Berthiaux, H. Analysis of Grinding Processes by Markov Chains / H. Berthiaux // Chemical Engineering Science, 55 (2000), P. 4117-4127.

60. Berthiaux, H. Applications of Markov Chains in Particulate Process Engineering: A Review/ H. Berthiaux, V. Mizonov // The Canadian Journal of Chemical Engineering. V.85, No.6, 2004, pp.1143-1168.

61. Berthiaux, H. Application of the theory of Markov chains to model different processes in particle technology/ H. Berthiaux, V. Mizonov, V. Zhukov //. Powder Technology 157 (2005) 128-137.

62. Djerroud, D. A non-linear cell model of the drying process in a paddle dryer/D. Djerroud, V. Mizonov, P. Arlabosse, H. Berthiaux // Proc. of Int. Conf. "Science and Technology of Particles", May 23-25 2007, Albi, France. CD edition, Paper 137.

63. Fan, L.S. A Stochastic Model for Particle Disintegration / L.S. Fan, , R.C. Srivastava // Grinding Mechanism, Chemical Engineering Science, 36 (1981), P. 1091-1096.

64. Marikh, K. Flow Analysis and Markov Chain Modelling to Quantify the Agitation Effect in a Continuous Mixer/K. Marikh, H. Berthiaux, V. Mizonov, E. Barantseva, D. Ponomarev //. Chemical Engineering Research and Design. 2006, 84(A11), pp.1059-1074.

65. Marikh, K. Algorithme de construction de modeles markoviens multidimensionnels pour le melange des poudre / K. Marikh, V. Mizonov, H. Berthiaux, E. Barantzeva, V. Zhukov // Congms Francophone de Gi-mie des Procedes GP 2001-Nancy 17-19 october 2001.

66. Mizonov, V. Application of Multi-Dimensional Markov Chains to Model kinetics of Grinding with Internal Classification / V. E. Mizonov, et al Proc. of the 10-th symposium on Comminution Heidelberg, 2002, 14 p.

67. Mizonov, V. Application of the Theory of Markovian Chains to Processes Analysis and Simulation/ V. Mizonov, H. Berthiaux, K. Marikh, V. Zhu-kov // Ecole des Mines d'Albi, 2000, 61 p.

68. Mizonov, V. Application of the Theory of Markov Chains to Simulation and Analysis of Processes with Granular Materials/ V. Mizonov, H. Berthiaux., V. Zhukov // Ecole des Mines d'Albi, 2002, 64p.

69. Mizonov, V. On Possible Instability of Throughputs in Complex Milling Circuits / V. Mizonov, et al Proc. of the 4th International Conference for Conveying and Handling of Particulate Solids, v.l. Budapest, Hungary, May 2003, pp. 8.23-8.26.

70. Tamir, A. Applications of Markov Chains in Chemical Engineering / A. Tamir. Amsterdam: Isevier publishers, 1998, - 604 p.