Прогнозные оценки специалистов показывают, что удвоение в России внутреннего валового продукта в первое десятилетие XXI века вызывает при сохранении энергоемкости на уровне 2000 г. рост расчетной потребности в топливо-энергетических ресурсах в 1,6 — 1,7 раза.
Расчетная потребность в энергоресурсах в 2020 году достигает в этом случае 2,67 млрд. т у.т. Предполагаемые объемы добываемых энергоресурсов значительно меньше (~ 1 млрд. т у.т. в 2010 году и 1,265 млрд. т у.т. в 2020 году). Выходом из ситуации может быть активное проведение энергосберегающей политики [9].
Наиболее быстро возрастают объемы потребления газообразного топлива. По потреблению газа Россия занимает второе место в мире после США. Расход газа на экспорт и на внутреннюю потребность страны за последние годы составляет 700−800 млрд. м3. Его преимуществами являются удобство транспортирования, рентабельность и экологическая чистота. Этот вид энергоносителя используется во всех металлургических процессах и энергетических установках. Доля газа в топливно-энергетическом балансе отрасли составляет около 30% [23,36].
Проблема энергосбережения является актуальной и для рассматриваемых в работе плавильных процессов и установок в системах производства стекла, что обусловлено значительными масштабами потребления природного газа и существенным потенциалом энергосбережения, т. е. разностью между действительным удельным расходом топлива и его теоретическим минимумом [70].
Для стекловаренных установок характерна значительная доля тепловых потоков с газовыми отходами — Qro и теплопроводностью через ограждение рабочего пространства в окружающую средуQo c. Комплексная регенерация этих тепловых отходов может дать существенный энергосберегающий эффект. В данном исследовании рассмотрены энергоэффективные решения 6 по регенерации QT0 и Qoc, такие как регенеративный подогрев исходного материала, термохимическая регенерация (ТХР) и термическая регенерация Qoc посредством фильтруемой изоляции (ФИ) теплонапряженных участков ограждения.
Обзор литературы показал, что взаимодействие этих энергосберегающих аппаратов в составе стекловаренной высокотемпературной установки (ВТУ) при различных нагрузках и условиях работы проанализировано недостаточно подробно. Это приводит к реализации в промышленных условиях неоптимальных конструкций и невыгодных режимов работы ВТУ и, как следствие, к появлению больших резервов снижения энергоемкости в стекловаренных теплотехнологических системах. Сказанное обуславливает актуальность исследования процессов в этих аппаратах в целях разработки методики их расчета и выработки технических мероприятий, направленных на снижение расхода топлива в стекловаренных ВТУ.
Целью настоящей работы является изучение тепломассообменных процессов в аппаратах регенерации тепловых отходов стекловаренной установки, таких как термохимическая регенерация, подогрев исходного материала, термическая регенерация посредством фильтруемой изоляции теплонапряженных участков ограждения, и оценка роста энергетической эффективности установки вследствие применения данных энергосберегающих мероприятий.
Для достижения поставленной цели был сформирован ряд относительно самостоятельных, но взаимосвязанных задач, в частности:
— разработка математической модели тепломассообменных процессов и химического реагирования в термохимическом регенераторе — реакторе паровой конверсии (РПК);
— оценка пригодности разработанной математической модели путем сравнения с имеющимися экспериментальными данными;
— проведение численных исследований локальных характеристик полей концентраций газовых компонентов и температур в канале реакционного элемента (РЭ) термохимического регенератора посредством разработанной математической модели;
— численное исследование подогревателя исходного материала и фильтруемой изоляцииприменение разработанной модели для оценки энергетической эффективности и экономической целесообразности использования предложенных энергосберегающих мероприятий в составе стекловаренной установки с комплексной регенерацией тепловых отходов.
Научная новизна:
1. Разработана математическая модель реактора паровой конверсии, что позволяет определить эффективность термохимической регенерации при различных режимах работы РГЖ и использовать модель при проектной разработке РПК с учетом сложного комплекса теплои массообменных явлений в реакционном элементе.
2. Произведен расчет влияния режимных параметров РПК на эффективность ТХР в рамках тепловой схемы стекловаренной установки, работающей на природном газе.
Установлено, что главными факторами, определяющими скорость конверсии природного газа, являются температура, и содержание метана в смеси. Скорость конверсии возрастает на начальном участке реакционного элемента с увеличением температуры, и снижается близи выхода в связи с уменьшением концентрации метана в газовой смеси.
3. Проведенное исследование комплекса теплои массообменных явлений в реакционном элементе при изменении конструктивных, режимных и каталитических параметров показало, что при увеличении тепловой нагрузки или снижении расхода природного газа в реакционном элементе относительно режима, соответствующего состоянию химического равновесия наблюдалось существенное увеличение температуры газа и каталитической вставки реакционного элемента вблизи выходного коллектора. На коротком участке реакционного элемента температура газа и стенки реакционного элемента значительно возрастает до уровня температуры греющих дымовых газов, что является опасным для системы ТХР, так как может привести к пережогу стенки РЭ.
4. Впервые проведено численное исследование реактора паровой некаталитической конверсии, сконструированного в Газовом Технологическом Институте (GTI) [73], позволяющего получить конвертированный газ без использования катализатора.
5. Исследование совместного влияния рассмотренных энергосберегающих мероприятий показало, что комплексное использование ТХР совместно с ФИ позволяет снизить расход топлива на 50% по сравнению с исходной схемой.
Достоверность. Результаты численных исследований с использованием одномерных моделей, реализованных в среде MathCAD и двумерных моделей, созданных посредством пакета численного моделирования FLUENT, хорошо согласуются с результатами экспериментов и численных исследований других авторов.
Практическая ценность.
Практическая ценность работы состоит в том, что ее результаты позволяют расчетным путем определять энергетическую эффективность использования предложенных энергосберегающих мероприятий в различных условиях эксплуатации и при термохимической переработке других видов топлив, например жидких, а также при процессах отложения на поверхности теплообменных устройств.
Проведены расчеты по проектным данным высокотемпературной плавильной установки в системе производства стекловолокна с реактором паровой конверсии НПО «Техэнергохимпром», для реактора паровой конверсии природного газа, которые продемонстрировали возможность значительного повышения энергетической эффективности установки.
Основные положения, выносимые на защиту;
• Математическая модель и программа расчета процессов теплои массообмена в термохимическом регенераторе — реакторе паровой конверсии;
• Результаты численного исследования процессов паровой конверсии в реакционных элементах термохимического регенератора различной конструкции:
— реакторы каталитической и некаталитической паровой конверсии с прямоточными реакционными элементами;
— реакторы каталитической и некаталитической паровой конверсии с реакционными элементами, выполненными по типу трубок Фильда;
• Математическая модель регенеративного подогревателя исходного материала и фильтруемой изоляции теплонапряженных участков ограждения стекловаренной установки;
• Результаты расчета локальных характеристик (полей скорости и температуры) для регенеративного подогревателя технологического материала и фильтруемой изоляции.
• Результаты расчета энергетической эффективности комплексной регенерации тепловых отходов в стекловаренной установке.
Апробация работы. Основные положения работы, результаты теоретических и расчетных исследований докладывались и обсуждались на международных и всероссийских конференциях: Одиннадцатой, Двенадцатой и Тринадцатой Международных научно-технических конференциях студентов и аспирантов. Радиоэлектроника, электротехника и энергетика (Москва, 2005;2007 гг.) — Третьей Всероссийской школе-семинаре молодых ученых и специалистов «Энергосбережение — теория и практика». (Москва, 21−29 сентября 2006 г.) — Четвертой Российской национальной конференции по теплообмену (Москва, 23−27 октября 2006 г.).
Публикации. Материалы, отражающие содержание диссертационной работы и полученные в ходе ее выполнения, представлены 8 публикациями.
1. Крылов А. Н., Овчинников Е. В. Расчет местных характеристик объекта с учетом термовлажностных процессов // Одиннадцатая Международная науч.-тех. конф. студентов и аспирантов. Радиоэлектроника, электротехника и энергетика: Тез. докл. В 3-х т. — М., — 2005 — С. 387.
2. Крылов А. Н., Попов С. К., Сергиевский Э. Д. Моделирование процессов в регенеративном подогревателе технологического материала // Промышленная энергетика — 2005.-№ 1- С. 42−45.
3. Крылов А. Н., Сергиевский Э. Д. Численное исследование процессов теплообмена в регенеративных аппаратах // Двенадцатая Международная науч.-тех. конф. студентов и аспирантов. Радиоэлектроника, электротехника и энергетика: Тез. докл. В 3-х т. — М., — 2006. — С. 423−424.
4. Крылов А. Н., Сергиевский Э. Д. Моделирование тепломассообменных процессов при конверсии природного газа // Труды Четвертой национальной конференции по теплообмену: В 8 томах. — М., — 2006. — Т. 3 — С.262−265.
5. E.D. Sergievsky, А. N. Krylov Heat transfer processes in melting furnace elements while using transpiration insulation, fuel reforming and raw material preheating // 13th International Heat Transfer Conference, Sydney, Australia, 2006.
6. Сергиевский Э. Д. Крылов A.H. Математическая модель фильтруемой изоляции высокотемпературной печи при направленной подаче охладителя // Вестник МЭИ.- 2006.-№ 5. С. 115−120.
7. М. В. Исаев, И. А. Султангузин, Э. Д. Сергиевский, А. Н. Крылов, А. П. Яшин Моделирование горения газа в коксовой печи // Третья Всероссийская школа — семинар молодых ученых и специалистов «Энергосбережениетеория и практика»: Сб. трудов. — М., — 2006. — С. 184−186.
8. Сергиевский Э. Д., Крылов А. Н. Повышение энергоэффективности высокотемпературных установок // Сборник тез. докл. 4-й научной школыконф. «Актуальные вопросы теплофизики и физической гидрогазодинамики» — Алушта, — 2006 — С. 116−119.
Настоящая работа выполнена под руководством д.т.н., профессора Сергиевского Э. Д. в Московском Энергетическом Институте (Техническом Университете) на кафедре Тепломассообменных процессов и установок.
Автор выражает благодарность руководителю, Станиславу Константиновичу Попову и всем сотрудникам кафедры за помощь в проведенных исследованиях, а также за понимание и поддержку.
2.6. Выводы по главе.
Детальное изучение выполненных в последнее время работ показало, что моделирование процессов, происходящих в аппаратах термохимической рекуперации природного газа, выполнено на уровне моделей с сосредоточенными параметрами. Анализ и сравнение результатов, полученных с помощью моделей такого типа, показал, что они не отражают в полной мере процессы термодинамики и кинетики каталитической паровой конверсии. Это приводит к значительным погрешностям. Несмотря на конструктивную простоту, в реакционных элементах с каталитическими вставками имеет место сложный комплекс последовательно — параллельных теплои массообменных явлений.
Многообразие факторов результирующего процесса теплопередачи в реакторе, не определяемого аддитивно комплексом переносных и кинетических явлений, послужило основанием для составления модели, позволяющей получить локальные характеристики различных объектов реактора паровой конверсии.
В ряде работ, посвященных изучению процессов, происходящих в топливных ячейках, используются подробные математические модели, позволяющие находить локальные характеристики реагирующих веществ. По своей структуре эти модели близки к той, что представлена в данной работе. Однако, процессы в топливных элементах значительно отличаются от процессов, протекающих при конверсии природного газа. Поэтому полученные в этих работах результаты затруднительно использовать для исследования процессов конверсии.
На основании проведенного анализа сформулирована цель и конкретные задачи, решение которых позволит проанализировать особенности работы системы ТХР.
ГЛАВА 3. ТАРИРОВОЧНЫЕ РАСЧЕТЫ. ОЦЕНКА ПРИГОДНОСТИ ПРЕДЛАГАЕМЫХ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ.
На различных этапах создания расчетной методики было проведено сравнение полученных с ее помощью результатов с имеющимися литературными данными. Несмотря на ограниченность опубликованных данных, тарировка позволила подтвердить пригодность предлагаемой математической модели РПК, поскольку во всех случаях сравнение показало приемлемые результаты.
Результаты расчетов по кинетической модели и полученные данные по теплои массообмену последовательно оценивались следующим образом. Вначале было проведено сравнение температур газов, полученных с помощью предложенной одномерной модели, и полученных посредством двумерной модели FLUENT, затем была осуществлена параметрическая идентификация кинетической модели с использованием экспериментальных данных, полученных при исследовании каталитической конверсии.
3.1. Сравнение температур и состава смеси газов, полученных с помощью предложенной модели с распределенными параметрами, и полученных при помощи двумерной модели Fluent.
Расчеты с привлечением вычислительного комплекса FLUENT позволили определить локальные характеристики на входе в реакционные элементы, перепады давления и оценить изменения локальных величин при организации процесса взаимодействия смеси.
Распределение температуры газовой среды по длине реакционного элемента представлено на рис. 3.1.
Z, м.
Рис 3.1. Распределение температуры газовой среды по длине реакционного элемента.
1 — данные по модели с распределенными параметрами- 2 — данные по двумерной модели.
Значения температур, полученные по предложенной одномерной модели с распределенными параметрами, показаны на рис. 3.1. кривой 1.
Кривая 2 представляет собой среднемассовые значения температур нагреваемой среды, полученных с помощью FLUENT, в различных поперечных сечениях элемента. Характерный обратный изгиб кривой 2 на участке длиной около 1 м от входа и более низкие температуры газа объясняются влиянием начального участка на теплообмен. На оставшейся длине трубы течение носит установившийся характер, и динамика изменения температуры газа аналогична полученной в одномерной модели.
Выходная температура газа, рассчитанная по модели с сосредоточенными параметрами, составила 569 К, а по двумерной модели.
551 К. Расхождение результатов, полученных с помощью данных моделей, не превышает 12%.
Численное исследование позволило сделать вывод о том, что влияние начального участка и неравномерностей в потоке не оказывает существенного влияния на характер тепломассообменных процессов на большей части реакционного элемента рекуператора, следовательно для расчета процессов в реакционном элементе допустимо использовать одномерную модель.
3.2. Тарировка кинетической модели химических процессов.
Кинетическая модель, разработанная для одномерного расчетного алгоритма, позволяет рассчитать сложный механизм реакции конверсии. С ее помощью удается найти скорости образования и исчезновения отдельных газовых компонентов, результирующую скорость реакции и вызванное ей теплопотребление на каталитической вставке. Использование данной модели в составе общей математической модели реакционного элемента позволяет определить любые параметры газа в любой точке канала.
Параметрическая идентификация кинетической модели выполнена с использованием экспериментальных данных [59], полученных при исследовании паровой каталитической конверсии.
Рис. 3.2. демонстрирует изменение концентрации компонентов реагирующей газовой смеси по длине реакционного элемента при постоянной температуре конверсии 1073 К. Сплошными линиями показаны концентрации метана, водорода, водяного пара, углекислого газа и оксида углерода, полученные по разработанной кинетической модели, крестикамиэкспериментальные данные [59]. Концентрации газа непрерывно меняются вследствие протекания химической реакции, когда газ проходит через реактор. Из сравнения результатов видно, что расхождение между данными эксперимента и численного расчета не превышают 13%.
Удовлетворительное соответствие расчетных и опытных данных дает основание утверждать, что задача идентификации решена, и математическую модель тепломассообмена в реакторе паровой конверсии можно использовать для исследования и разработки конструкции этого важного элемента системы ТХР. концентрация.
Рис. 3.2. Изменение состава газовой смеси по длине реакционного элемента.
3.3. Тарировка расчета реакционного элемента с двойной циркуляцией газовой среды. Сравнение температуры смеси газов с результатами двумерной модели Fluent.
Для учета двойной циркуляции нагреваемой среды была использована усовершенствованная модель с распределенными параметрами, позволяющая также определить температуру парогазовой смеси во внутреннем канале.
Оценка пригодности результатов, полученных по этой модели, была проведена при помощи вычислительного комплекса FLUENT. С его помощью была составлена упрощенная двумерная осесимметричная модель реакционного элемента с двойной циркуляцией нагреваемой среды. Для простоты было принято, что реакция конверсии в реакционном элементе не происходит.
Данная модель позволила найти распределение температуры газовой среды по длине трубки Фильда. Температуры, полученные с использованием предложенной модели с распределенными параметрами и двумерной модели FLUENT, представлены на рис. 3.3.
Сравнение представленных на графике кривых 1,2 с кривыми 3,4 показало, что расхождение между результатами не превышает 15%. Следовательно, предложенная модель с распределенными параметрами может использоваться для расчета реакционных элементов с двойной циркуляцией.
Проведенные тарировочные расчеты позволили оценить достоверность результатов разработанных одномерных математических моделей с распределенными параметрами. Последовательное поэтапное сравнение данных по теплои массообмену в реакционном элементе с имеющимися данными экспериментов и результатами двумерной модели, позволило сделать вывод о том, что модель тепломассообмена в реакторе паровой конверсии можно использовать для исследования и разработки конструкции данного аппарата. длина реакционного элемента, м.
Рис. 3.3. Распределение температур газа в реакционном элементе по его длине.
1,2- температура парогазовой смеси во внутреннем канале, и реагирующего газа в кольцевом канале, найденная по предложенной модели- 3,4-температура парогазовой смеси во внутреннем канале, и реагирующего газа в кольцевом канале, найденная с использованием FLUENT- 5 — температура греющих дымовых газов.
ГЛАВА 4. ПАРАМЕТРИЧЕСКИЕ РАСЧЕТЫ СИСТЕМЫ ТЕРМОХИМИЧЕСКОЙ РЕГЕНЕРАЦИИ.
В данном разделе представлены результаты численного исследования теплои массообмена при паровой конверсии природного газа в реакционном элементе с цилиндрической каталитической вставкой при изменении конструктивных, режимных и каталитических параметров, и определение области их возможного варьирования.
Применение процесса паровой конверсии метана для термохимической регенерации тепла в данном технологическом цикле потребовало особого конструктивного оформления и определило специфичные режимные условия его организации. Были сформулированы [24,48] следующие основные требования к термохимическому регенератору:
— обеспечение условий реализации процесса при расходах водяного пара, существенно не превышающих стехиометрический;
— процесс не должен сопровождаться выпадением углерода;
— степень завершения процесса конверсии может быть ограничена на уровне Хсн4=0.8−0.9 от равновесной;
— малые сопротивления (порядка сотен мм водяного столба) тракта реагирующей смеси, что позволит использовать сетевой газ низкого давления и пар атмосферных параметров;
— надежность и простота конструкции, не требующей постоянного обслуживания при эксплуатации.
При проектно-конструкторской проработке реактора паровой конверсии появляется ряд вопросов по выбору теплотехнически и экономически оправданных режимно-конструктивных параметров. Это, во-первых, характеристики реакционного элемента:
— геометрические характеристики внутреннего и кольцевого канала;
— скорость парогазовой смеси во внутреннем и кольцевом канале;
— расход природного газа;
— удельный расход водяного пара;
— температура исходной парогазовой смеси;
Во-вторых, необходимы рекомендации по конструкциям реактора паровой конверсии в целом:
— взаимное движение дымовых газов и реагирующей парогазовой смеси;
— тип конструкции с точки зрения внешнего теплообмена;
— материал реакционных труб;
— диаметр, длина и количество труб при компоновке конкретного аппарата, определяемые как внутренней, так и внешней задачей.
В-третьих, требуется сформулировать некоторые общие требования, ограничения и рекомендации для проектирования системы ТХР, а именно:
— распределение регенерируемого тепла между энергосберегающими аппаратами;
— возможность конверсии части природного газа;
— целесообразность разбивки ТХР на ступени;
Целесообразным является исследование взаимодействия физических и химических процессов, протекающих в реакционном элементе ТХР. Это позволит не только качественно обосновать влияние конструктивных, режимных и каталитических параметров на процесс конверсии, но и определить пути повышения эффективности и надежности ТХР.
Таким образом, конечной целью является определение характеристик и закономерностей при изменении конструктивных, режимных и каталитических параметров, которые необходимы для проектирования ТХР.
4.1. Шероховатости в каналах реакционного элемента.
Длина участка стабилизации в процессе перехода к турбулентному течению зависит, как известно, от ряда факторов: уровня начальной турбулентности, шероховатости поверхности, неизотермичности, и др. Вследствие формирования отложений меняется шероховатость поверхности канала. При достижении некоторого критического значения высоты шероховатости возможна интенсификация процессов перестройки профилей скорости и температуры потока к характерным профилям развитого турбулентного течения, при этом уменьшается длина начального участка и интенсифицируется теплообмен. Возможна ситуация, когда интенсификация теплообмена на начальном участке вследствие роста шероховатости из-за отложений превосходит увеличение термосопротивления из-за роста толщины последних [22].
Наибольшее количество отложений образуется в местах поворотов, сужения потока, дефектов поверхности РЭ.
Для учета влияния шероховатостей поверхности на теплопередачу в реакционном элементе реактора паровой конверсии использована поправка для коэффициента конвективной теплоотдачи с поверхностей стенки РЭ и вставки. Согласно [12] в исследуемом диапазоне параметров течения газовой смеси шероховатости усиливают теплоотдачу приблизительно на 40%. Следовательно, для расчетов, в математической модели поправочный коэффициент принят равным =1,4. 12].
4.2. Повышение температуры стенки выше максимально допустимой по условию прочности.
Проведенное численное исследование системы термохимической регенерации показало, что при определенном сочетании режимных параметров наблюдалось повышение температуры стенки реакционного элемента РПК выше максимально допустимой по условию прочности. Установлено, что предпосылками к этому является снижение объемного расхода природного газа или увеличение температуры греющих дымовых газов.
Поясним причину возникновения данного явления. При снижении объемного расхода природного газа на входе увеличивается время пребывания газовой смеси в реакционном элементе. Это приводит как к более полному превращению парогазовой смеси в конвертированный газ, так и к более интенсивному нагреву газовой смеси. Кроме того, возрастает температура каталитической вставки.
Подобные результаты дает и увеличение температуры греющих дымовых газов. Это увеличение приводит к росту теплового потока от стенки реакционного элемента к вставке и газовой смеси. Данный рост тепловой нагрузки приводит к тому, что смесь газов имеет большую температуру, а на выходе степень завершения конверсии повышается. Динамика роста температуры стенок РЭ при повышении тепловой нагрузки противоточного реакционного элемента (увеличении температуры греющей среды) показана на рис. 4.1.
TW, K.
1200 1300 1400 1500 1600 1700.
Tfiue' ^.
Рис. 4.1. Зависимость температуры стенок РЭ Tw от средней температуры греющих дымовых газов Tjjue.
Чем ниже объемный расход природного газа на входе в ТХР или чем выше тепловая нагрузка, тем сильнее повышается температура конвертированного газа и тем выше степень конверсии на выходе из реакционного элемента. Рост степени конверсии продолжается до тех пор, пока кинетическая скорость образования продуктов конверсии в РЭ не станет равной скорости их разложения и образования исходных веществ, то есть пока не наступит состояние равновесия. В этом случае степень завершения реакции на выходе из РЭ близка к 100%, и концентрации исходных реагентов в газовой смеси малы (при данном уровне температур).
При дальнейшем повышении тепловой нагрузки или снижении объемного расхода природного газа в реакционном элементе наблюдалось увеличение температуры газа и каталитической вставки вблизи выхода из РЭ.
Данное явление в РЭ является опасным для системы ТХР, так как может привести к пережогу стенки РЭ. Значение температуры стенки Tw, и, следовательно, опасность пережога стенки оказались зависящими от особенностей протекания процесса конверсии в зоне кризиса.
Степень завершения реакции конверсии в области повышенной температуры достигает значительных величин, что приводит к значительному снижению скорости конверсии в этой области. Это происходит вследствие очень низкой концентрации исходных реагентов в газовой смеси. Прекращение химического реагирования приводит к тому, что из зоны реакции перестает поглощаться теплота химической реакции. На коротком участке реакционного элемента температура газа и реакционного элемента значительно возрастает до уровня температуры греющих дымовых газов.
Данное явление показано на рис. 4.2. Предельно допустимая температура стенки для стали 12Х18Н10Т, использованной для изготовления реакционных элементов, составляет 1250 °C или 1523 К. Кривые 1 и 2 характеризуют изменение температур реакционного элемента и греющих дымовых газов при нормальном режиме его работы. В данном случае температура стеки РЭ ниже предельно допустимой, показанной на графике прямой 5. Кривые 3 и 4 характеризуют изменение температур реакционного элемента и греющих дымовых газов при значительном увеличении тепловой нагрузки системы ТХР (при средней температуре греющей среды 1673 К). Из данного графика видно, что в этом случае на расстоянии около 1.5 м от входа в РЭ температура его стенок резко возрастает и становится выше предельно допустимой. Это приводит к опасности пережога стенок РЭ.
Рис. 4.2. Распределение температуры стенки РЭ и греющих при нормальном и аварийно опасном режиме работы 1,2- температуры стенки РЭ и дымовых газов при нормальном режиме работы- 3,4- температуры стенки РЭ и дымовых газов при аварийно опасном режиме работы- 5 — предельно допустимая температура стенки.
4.3. Оценка влияния химических и диффузионных процессов на теплообмен в реакционном элементе.
При определении коэффициента теплоотдачи с поверхности каталитической вставки требуется знать соотношение характерных времен химического превращения (тхим), пребывания газа в канале (тгидр), и турбулентной диффузии поперек канала {тдифф), которые определяют необходимость учета тех или иных физических процессов [22,25].
Для определения (тхим) используем результаты расчета по кинетической модели химических процессов в реакционном элементе. Характерное время пребывания определяется какгидр ~ jj, что позволяет учесть изменение времени пребывания вследствие изменения плотности реагирующей газовой смеси. Характерное время турбулентной диффузии в трубе тдифф.
PL определяется как тдифф ~ ?) > где DT — коэффициент турбулентной диффузии, D — гидравлический диаметр канала [22].
Скорость протекания химической реакции конверсии существенно меняется по длине реакционного элемента (см. рис. 4.3.).
Р> МОЛЬ кг-с 3 2 1.
О 0.5 1 1.5 2 2.5.
Длина реакционного элемента, м.
Рис. 4.3. Изменение скорости химической реакции конверсии по длине реакционного элемента.
Это происходит вследствие значительного изменения температуры РЭ (на величину порядка 300 К), и концентрации исходных веществ в реакционном элементе (для метана характерно изменение концентрации в 5 раз).
Для того, чтобы учесть различное влияние химической реакции на теплообмен на разных участках реакционного элемента, он был разделен по длине на три зоны равной длины (участки 1−3 в табл. 4.1). В первой зоне химическая реакция протекает наиболее интенсивно, во второй и третьей зоне скорость реакции значительно снижается.
Учет влияния химической реакции на теплообмен осуществлен путем введения поправки к при определении величины коэффициента теплоотдачи [35]. Величина к есть функция отношения характерных времен пребывания газа в канале в реакционном элементе и химического превращения: х ^ гидр W Т v хим J.
Характерные времена химических превращений, времена пребывания газа в канале, значения числа Дамкёллера Da, представляющего собой отношение времени пребывания ко времени химического превращения и поправка ki для расчета характеристик теплообмена представлены в табл. 4.1 и табл. 4.2:
ЗАКЛЮЧЕНИЕ
.
В данной диссертационной работе получены следующие основные результаты:
1. Разработана математическая модель реактора паровой конверсии, что позволяет определить эффективность термохимической регенерации при различных режимах работы РПК и использовать модель при проектной разработке РПК с учетом сложного комплекса теплои массообменных явлений в реакционном элементе.
2. Проведенная подробная поэтапная тарировка разработанной модели с использованием двумерных моделей и результатов экспериментов показала, что влияние начального участка и неравномерностей в потоке не оказывает существенного влияния на характер тепломассообменных процессов на большей части реакционного элемента ТХР. Для расчета процессов в реакционном элементе прямоточной конструкции и с двойной циркуляцией нагреваемой среды допустимо использовать одномерную модель.
3. Произведена оценка влияния режимных параметров на эффективность системы ТХР. Установлено, что система ТХР обеспечивает наиболее эффективную регенерацию теплоты газовых отходов при составе парогазовой смеси на входе с соотношением пара к газу, равным 2:1.
Установлено, что изменение внутреннего диаметра стенки реакционного элемента существенно влияет на интенсивность теплопередачи в реакционном элементе. Увеличение диаметра приводит к снижению температуры стенки и вставки на начальном участке реакционном элементе, и ее увеличению вблизи выхода.
Показано, что главными факторами, определяющими скорость конверсии природного газа, являются температура и содержание метана в смеси. Скорость конверсии возрастает на начальном участке с увеличением температуры и снижается близи выхода в связи с уменьшением концентрации метана в газовой смеси.
4. Проведенное исследование комплекса теплои массообменных явлений в реакционном элементе при изменении конструктивных, режимных и каталитических параметров показало, что при увеличении тепловой нагрузки или снижении расхода природного газа в реакционном элементе относительно режима, соответствующего состоянию химического равновесия наблюдается существенное увеличение температуры газа и каталитической вставки реакционного элемента вблизи выходного коллектора. На коротком участке реакционного элемента температура газа и реакционного элемента значительно возрастает до уровня температуры греющих дымовых газов, что является опасным для системы ТХР, так как может привести к пережогу стенки реакционного элемента.
5. Произведена оценка химического и диффузионного числа Дамкёллера для процесса конверсии в реакционном элементе. На основании результатов оценки расчет характеристик теплообмена в реакционном элементе произведен по зависимостям для нереагирующих потоков.
6. Впервые проведено численное исследование реактора паровой некаталитической конверсии, сконструированного в Газовом Технологическом Институте (GTI), позволяющего получить конвертированный газ без использования катализатора. В работе было установлено, что применение данного реактора приводит к значительному увеличению эффективности теплотехнического объекта, что подтверждается данными проведенного эксперимента.
7. Исследование совместного влияния рассмотренных энергосберегающих мероприятий показало, что применение ТХР теплоты газовых отходов для стекловаренной установки может обеспечить снижение расхода топлива на 25% относительно исходной схемы с регенеративным подогревом воздуха. Использование ТХР совместно с фильтруемой изоляцией в комплексной схеме позволяет снизить расход топлива на 50% ниже расхода в исходной схеме.