Авторефераты по всем темам  >>  Авторефераты по техническим специальностям  

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

Кулешов Олег Юрьевич

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

Специальность 05.14.04 - Промышленная теплоэнергетика

Автореферат

диссертации на соискание ученой степени

доктора технических наук

Саратов 2012

Работа выполнена в Федеральном государственном бюджетном

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

Саратовский государственный технический университет

имени Гагарина Ю.А.

Научный консультант:	доктор технических наук, профессор Седелкин Валентин Михайлович
Официальные оппоненты:	Кузнецов Валерий Алексеевич доктор технических наук, профессор, ФГБОУ ВПО Белгородский технологический университет им. Шухова В.Г., профессор кафедры Энергетика теплотехнологий Парамонов Александр Михайлович доктор технических наук, профессор, ФГБОУ ВПО Омский государственный технический университет, профессор кафедры Теплоэнергетика Щёлоков Анатолий Иванович доктор технических наук, профессор, ФГБОУ ВПО Самарский государственный технический университет, заведующий кафедрой Промышленная теплоэнергетика
Ведущая организация:	ФГБОУ ВПО Национальный исследовательский университет МЭИ (г. Москва)

Защита состоится л 23 октября 2012 г. в 14.00 часов на заседании диссертационного совета Д212.242.07 при ФГБОУ ВПО Саратовский государственный технический университет имени Гагарина Ю.А. по адресу: 410054, г. Саратов, ул. Политехническая,77, корпус 1, ауд. 319.

С диссертацией можно ознакомиться в научно-технической библиотеке ФГБОУ ВПО Саратовский государственный технический университет имени Гагарина Ю.А.

Автореферат разослан л____ __________ 2012 г.

Ученый секретарь

диссертационного совета                                        Ларин Е.А.

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

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

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

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

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

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

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

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

Предметом исследований является сложный, сопряженный теплообмен в промышленных печах.

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

Задачи исследования:

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

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

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

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

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

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

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

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

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

10.аТехнико-экономический анализ эффективности предложенных решений по совершенствованию тепловых режимов и конструкций промышленных печей.

Научная новизна:

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

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

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

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

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

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

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

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

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

Практическая значимость и реализация работы:

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

2.аРезультаты математического моделирования и экспериментальных исследований зонального и локального теплообмена в нагревательных и реакционных трубчатых печах, а также разработанные на их основе рекомендации использованы при наладке и оптимизации режимов работы печи пиролиза бензина типа SRT-II этиленовой установки ЭП-450 ОАО Нижнекамскнефтехим, при реконструкции печей типа ББ1 Норильского ГРС с целью увеличения их тепловой мощности, при совершенствовании режимов работы трубчатой печи каталитической паровой конверсии метана ППР-600 установки производства аммиака ОАО Воскресенские минеральные удобрения, при анализе схем отопления и режимов эксплуатации печи ЦД4 установки первичной перегонки нефти АВТ-6 Новополоцкого НПЗ.

3.аПолученные результаты использованы в учебном процессе СГТУ имени Гагарина Ю.А, при подготовки инженеров по специальностям Энергетика теплотехнологий, Промышленная теплоэнергетика, Машины и аппараты химических производств, Машины и аппараты пищевых производств.

4.аПредложенные методики расчета среднезональных и локальных характеристик теплообмена могут послужить базой для разработки ведомственных руководящих технических материалов (ВРТМ) и стандартов предприятия (СТП) по проектированию трубчатых печей, подогревателей газа и нефти, хлебопекарных печей, факельных установок для сжигания сбросных газов.

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

Положения, выносимые на защиту:

1. Зональная математическая модель сложного теплообмена в камерных печах, а также методология деления рабочего пространства на расчетные зоны.

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

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

4. Математические модели и методики расчета сопряженного теплообмена в технологических печах различных производств: реакционных трубчатых печах пиролиза и каталитической паровой конверсии углеводородов; конвейерных хлебопекарных печах.

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

6. Результаты численных исследований теплообмена в печах и факельных установках в зависимости от конструктивных и режимных параметров.

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

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

Апробация работы. Основные научные и прикладные результаты диссертационной работы докладывались на: Седьмой Всесоюзной конференции по радиационному теплообмену (Ташкент, 1991); Международном совещании-семинаре Теплофизические проблемы промышленного производства (Тамбов, 1992); Минском международном форуме по тепломассообмену - ММФ-92 (Минск, 1992); 1-ой Российской национальной конференции по теплообмену (Москва,1994); 12-th Intern. Cong. of Chem. and Process Eng. CHISAТ96 (Praha, 1996); The First European Congress on Chemical Engineering- ECCE-1 (Milano, 1997); Международной конференции Проблемы эффективного использования энергоносителей и низкосортных топлив в промышленности (Саратов, 1998); 2-ой Российской национальной конференции по теплообмену (Москва, 1998); 13-th International Congress of Chemical and Process Engineering CHISAТ98 (Praha, 1998); 2-ой Международной научной конференции Повышение эффективности тепломассообменных процессов и систем (Вологда, 2000); 4-ом Минском Международном форуме по тепломассообмену ММФ-2000 (Минск, 2000); 14-th International Congress of Cemical and Process Engineering CHISAТ2000 (Praha, 2000); Всероссийской научной конференции Тепло- и массообмен в химической технологии (Казань, 2000); 2-ой Всероссийской научно-технической конференции Прогрессивные процессы и оборудование металлургического производства (Череповец, 2001); Международной научно-технической. конференции Энергосбережение в теплоэнергетических системах (Вологда, 2001); 4-ой Международной Теплофизической школе Теплофизические измерения в начале XXI века (Тамбов, 2001); Юбилейной Международной научно-практической конференции Пищевые продукты XXI века (Москва, 2001); Международной научно-технической конференции Моделирование, оптимизация и интенсификация производственных процессов и систем (Вологда, 2001); 3-ей Международной научно-технической конференции Повышение эффективности теплообменных процессов и систем (Вологда, 2002); Научно-методической конференции Современные научные и информационные технологии (Энгельс, 2003); Научно-технической конференции, посвященной 50-летию ЭТИ СГТУ (Энгельс, 2006); Шестой Международной теплофизической школе Теплофизика в энергосбережении и управлении качеством (Тамбов, 2007); XXI Международной научной конференции Математические методы в технике и технологиях - ММТТ-21 (Саратов, 2008); Научно-практической конференции Синтез инноваций: направления и перспективы (Саратов, 2009); Пятой Российской национальной конференции по теплообмену (Москва,2010); ХХV Международной научной конференции Математические методы в технике и технологиях - ММТТ-25 (Волгоград, 2012); XIV Минском международном форуме по тепломассообмену - ММФ- XIV (Минск, 2012).

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

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

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

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

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

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

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

Более универсальным и эффективным представляется зональный метод, который основан на разбиении расчетной области на однородные объемные и поверхностные зоны и записи системы нелинейных алгебраических уравнений зональных тепловых балансов, коэффициенты переноса в которых вычисляются на основе моделирования процессов теплообмена. Современные зональные методы в сочетании с имитационным моделированием излучения на основе метода Монте-Карло обладают повышенной точностью (принципиально не ограниченной) в отношении расчета радиационной составляющей теплообмена в пространственно неоднородных излучающих и поглощающих системах со сложной геометрией границ. Поэтому зональные методы получили интенсивное развитие и применение для расчета реальных печей и топок. Основы прикладного зонального метода расчета радиационного и сложного теплообмена были заложены в работах Г.Л. Поляка, Ю.А. Суринова, А.С.аНевского, В.Н. Адрианова, Х. Хоттеля. Дальнейшее развитие методы расчета топок получили в работах В.М. Седелкина, В.Г. Лисиенко, Ю.А. Журавлева и др.

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

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

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

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

В свернутой форме система уравнений зональных тепловых балансов имеет вид:

;                                        (2.1)

где N - общее число зон в расчетной области - зональной геометрической модели печи; M - число зон, непосредственно контактирующих с j-ой зоной; Ti - абсолютная температура i-ой зоны; Pij - коэффициент радиационного обмена между зонами i и j; Ωij - коэффициент конвективно-турбулентного обмена между зонами i и j; Cj = f(Qj) - свободный член уравнения, включающий в себя внутренний тепловой источник или сток в j-ой зоне; i, j - зоны источник и приемник теплоты.

Коэффициенты радиационного обмена в (2.1) записываются в виде:

                                       (2.2)

где ; - приведенные разрешающие обобщенные угловые коэффициенты излучения в полосе спектра k; N′ - число поверхностных зон в расчетной области; (N-N′) - число объемных зон в расчетной области; Fi - площадь поверхностной зоны; Vi - объем газовой зоны; σ0 - постоянная Стефана-Больцмана; Z - число рассматриваемых полос спектра в модели излучения продуктов сгорания; bi,k - доля излучения АЧТ в k-ой полосе спектра при температуре i-ой зоны; χi,k - коэффициент поглощения продуктов сгорания в k-ой полосе спектра при физических условиях в i-ой зоне; εi - степень черноты i-ой поверхностной зоны-источника излучения; εj - степень черноты j-ой поверхностной зоны-приемника излучения; δij - дельта-символ Кронекера.

Коэффициенты конвективно-турбулентного обмена в уравнениях (2.1) записываются в виде

               (2.3)

где αij - коэффициент конвективной теплоотдачи между объемной и поверхностной зонами i и j; ср,i - теплоемкость продуктов сгорания в i-ой объемной зоне; ρu+ρuт - массовый поток газов между объемными зонами i и j с учетом турбулентной составляющей.

Первый член уравнений (2.1) учитывает радиационный перенос, второй - конвективно-турбулентный перенос, третий - тепловой источник или сток в зоне.

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

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

На основе анализа конструкций камерных печей выделены возможные пространственные геометрические типы объемных и поверхностных зон для печей коробчатого и цилиндрического типа. Зоны также классифицированы по оптическому типу в зависимости от их свойств и положения в зональной геометрической модели печи. Для печей коробчатого типа всего выделено 19 геометрических и 7 оптических типов зон, для печей цилиндрического типа - 15 геометрических и 7 оптических типов зон.

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

Методика статистических испытаний (реализация метода Монте-Карло) основана на прослеживании пути лучей (или их пучков) в зональной геометрической модели и включает в себя: 1) выбор излучающей точки в пределах зоны; 2) выбор направления единичного луча; 3) определение координат и длины луча и его поглощения при пересечении объемной зоны; 4) определение взаимодействия луча с поверхностной зоной, на которую он падает (для обычных поверхностей происходит полное поглощение и на этом единичное испытание заканчивается; для поверхностей теплообмена, в частности - трубных экранов, происходит частичное поглощение; для зеркальных поверхностей - отражение по законам геометрической оптики); 5) статистическую обработку результатов моделирования излучения с данной зоны на другие зоны системы, т.е. вычисление первичных ОУК. Методика учитывает селективность излучения на основе модели широких полос спектра.

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

Предложены принципы построения зональных геометрических моделей печей различного типа. Разработаны алгоритмы и программы для расчета ОУК в печах коробчатого и цилиндрического типов.

Реализация метода имитационного моделирования излучения составляет основную часть всех зональных расчетов. Полный расчетный алгоритм зонального метода представлен на рис. 3.1. (сплошные линии). Все расчетные процедуры включены в главный итерационный процесс решения нелинейной тепловой задачи, в том числе и вычисление ОУК.

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

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

Соотношения для коррекции ОУК имеют вид:

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

или ;                                                (3.1)

для объемной зоны-приемника излучения

или ;                                (3.2)

где ψij,k - обобщенный угловой коэффициент излучения между зонами i и j в k-ой полосе спектра; εj,k - степень черноты j-ой объемной зоны; Dij,k - пропускательная способность среды; верхний индекс л0 означает базовый режим.

Среднегеометрическая пропускательная способность среды между зонами i и j имеет выражение:

;                                                                (3.3)

где Lij - среднегеометрическое расстояние между зонами i и j ; χk - коэффициент поглощения среды в k-ой полосе спектра.

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

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

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

Получено уточненное выражение для степени черноты газовых зон в рамках зонального подхода путем разложения экспоненты в ряд Маклорена:

,                                                (3.4)

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

;                                                                (3.5)

где liа=а4Viа/аFiоб - среднегеометрическая длина пути луча в объемной i-ой зоне; Vi и Fiоб - объем зоны и площадь ограничивающих плоскостей для газовой зоны i, соответственно; χi,k - коэффициент поглощения продуктов сгорания в объемной зоне i в k-ой полосе спектра излучения; N0 - число учитываемых членов ряда Маклорена.

С учетом выражения (3.4) предложено уточненное выражение для РОУК объемной зоны-источника с учетом корреляции с интегральной по спектру степенью черноты газовой среды

,                                        (3.6)

,                                        (3.7)

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

Относительный коэффициент поглощения среды в полосе спектра в выражении (3.6) задается в виде:

δχi,kа=аχi,kа/-χi  ,                                                                        (3.8)

,  ,

где - средний коэффициент поглощения среды в ИК-области спектра; χi,k - действительный коэффициент поглощения среды в полосе спектра; bi,k - доля излучения АЧТ в полосе спектра.

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

Выражение для локального удельного теплового потока к элементарной площадке dFM с центральной точкой М имеет вид:

;                                        (4.1)

где N - общее число зон в расчетной области; εj - степень черноты поверхности нагрева в зоне j; ТМ - абсолютная температура поверхности трубы в точке М; σ0а - постоянная Стефана-Больцмана; αiM - локальный коэффициент конвективной теплоотдачи к элементарной площадке dFM от контактирующей объемной зоны i; PiM - локальный коэффициент радиационного обмена.

окальные коэффициенты радиационного обмена имеют выражение:

                               (4.2)

где N′ - число объемных зон в расчетной области; (N - N′) - число поверхностных зон в расчетной области; Z - число рассматриваемых полос спектра излучения; bi,k - доля излучения АЧТ в k-ой полосе спектра при температуре Ti ; χi,kа - коэффициент поглощения продуктов сгорания в k-ой полосе спектра в объемной зоне i; εi - степень черноты поверхностной зоны i; Vi - объем газовой зоны i; Fi - площадь поверхностной зоны i; ΨiM,k - локальные РОУК между зоной i и элементарной площадкой dFM в k-ой полосе спектра излучения; Kεаi - поправка на нелинейность степени черноты объемной зоны i для k-ой полосы спектра.

окальные РОУК рассчитываются по выражению:

;                                                        (4.3)

где ψiM,k - локальный ОУК между зоной i и элементарной площадкой dFM в k-ой полосе спектра излучения; ψnM,k - локальные ОУК между зонами n расчетной области и элементарной площадкой dFM в k-ой полосе спектра излучения; Ψin,kа - среднезональные РОУК между зоной i и зонами n расчетной области в k-ой полосе спектра излучения; Rn - коэффициент диффузного отражения для поверхностной зоны n или коэффициент изотропного рассеяния для объемной зоны n.

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

Для реализации метода впервые получены соотношения взаимности для ОУК как для поверхностных, так и для объемных зон:

для поверхностной зоны,

,                                                                (4.4)

для объемной зоны,

.                                                (4.5)

Выражение локальных ОУК через местные ОУК имеет вид:

                                       (4.6)

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

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

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

Реакционные трубчатые печи характеризуются узкой топочной камерой коробчатого типа, одно-и двухрядными двухсветными экранами и системами обогрева от настенных излучающих горелок или настильных факелов (рис.а5.1). Это обеспечивает высокую интенсивность радиационного теплообмена, но в то же время предъявляет жесткие требования к условиям теплоподвода. Наиболее жесткими режимами нагрева характеризуются реакционные трубчатые печи пиролиза и каталитической конверсии углеводородов, в которых реагирующий поток нагревается до температуры 850-900С при теплонапряженности поверхности труб до 100 кВт/м2 . Реакционные трубы выполняются из высоколегированных сталей, работающих на пределе своих прочностных возможностей.

а	б	в

Рис. 5.1. Схемы реакционных трубчатых печей

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

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

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

Описание сложного теплообмена в радиантной камере печи дается в рамках зонального подхода системой нелинейных алгебраических уравнений зональных тепловых балансов (2.1).

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

• для теплового потока

;                                (5.1)

• для температуры

,                                                                        (5.2)

где nj и lj - число и длина труб (или их участков) входящих в j-ю экранную зону; dтр - диаметр трубы; Тj - средняя температура экранной поверхности (наружной стенки труб) в пределах j-ой зоны; Тж - локальная температура внутреннего потока в трубах; k - локальный коэффициент теплопередачи через стенку труб, который вычисляется по формуле

;                                                                (5.3)

αвн - локальный внутренний коэффициент теплоотдачи; ∑(δ/λ)ст - локальная сумма термических сопротивлений стенки трубы с учетом слоя внутренних отложений.

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

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

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

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

Система уравнений продольного тепломассопереноса в ядре потока при δ→0:

уравнение теплопереноса

;                                                (5.4)

уравнение движения

;                                                        (5.5)

уравнение массопереноса

;                        (5.6)

уравнение состояния

.                                                                        (5.7)

Система уравнений радиального тепломассопереноса в пристеночном ламинарном слое δ:

уравнение массопереноса свободных радикалов

;                        (5.8)

уравнение теплопереноса

;                                        (5.9)

;                                                                        (5.10)

Расчетная область:

;                                        (5.11)

.

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

;

;                                                (5.12)

;

.                                                                        (5.13)

Обозначения: n - общее количество компонентов реагирующей смеси в кинетической модели процесса пиролиза; m - долгоживущие молекулярные компоненты реагирующей смеси; (n-m) - короткоживущие свободные радикалы; Т - температура; Р - давление; u - скорость; ср - теплоемкость; ρ - плотность; М - число компонентов реагирующей смеси; сi - концентрация i-го компонента реагирующей смеси; Hi - теплота i-ой реакции; Wi - скорость i-ой реакции; Di - коэффициент молекулярной диффузии i-го компонента; ji - поток массы i-го компанента; λ - коэффициент теплопроводности; qw - тепловой поток на стенке; αвн - коэффициент внутренней теплоотдачи; Rw - радиус трубы; δ - толщина пристеночного ламинарного слоя; ξтр и ζм - коэффициенты трения линейный и местный, соответственно; индекс f соответствует ядру потока; индекс w соответствует внутренней поверхности стенки трубы.

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

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

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

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

Исследованы характеристики сопряженного теплообмена в серийной печи пиролиза углеводородов типа SRT-II этиленовой установки ЭП-450 Нижнекамского нефтехимического комбината.

Конструктивная схема печи показана на рис 5.1а. Печь узкокамерная с центральным вертикальным однорядным экраном. На боковых стенах топки расположены излучающие плоскопламенные горелки фирмы John Zinc. Сырье пиролиза - прямогонный бензин. Топливо - метано-водородная фракция. В случае пиролиза бензина (основное сырье) змеевик выполняют разветвленным комбинированным из труб различного диаметра. Всего в топочной камере расположено 4 змеевика.

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

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

Рис. 5.2. Температура газов на расстоянии 10 см (1) и 50 см (2) от излучающей стенки: --  - расчет; , ♦ - эксперимент	Рис. 5.3. Плотность потока излучения, падающего на экран: --  - расчет; - эксперимент

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

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

Проведено расчетное исследование режима работы печи при перераспределении 20 % расхода топлива с последнего 5-го вертикального ряда горелок на 3-й ряд. Результаты расчета отражены в графической форме на рис.5.4, 5.5. За счет перераспределения теплоподвода происходит изменение теплонапляженности труб во 2-ой и 3-ей секциях змеевика на ≈10 кВт, что способствует росту теплоподвода в основную зону реакции на 8% и снижению теплоподвода на последнюю 10-ю трубу на 12%. В результате это смещает равновесие реакции пиролиза в сторону увеличения выхода этилена на ≈1%. При этом максимальная температура выходной трубы снижается на 15С (рис. 5.4).

Распределение удельных тепловых потоков по периметру труб №6 и №10 (второй и третьей секций) до и после регулирования показано на рис.5.5. Снижение максимальной теплонапряженности и температуры стенки выходной трубы способствует увеличению срока ее службы, замедлению коксообразования и увеличению рабочего пробега печи с 70 до 100 суток, т.е. на ≈40%.

Рис. 5.4. Распределение параметров по длине трубчатого реактора: --  - расчетные параметры до изменения режима; ----  - то же после; , × - экспериментальные данные	Рис. 5.5. Локальная теплонапряженность по окружности труб: --  - до изменения режима; -----  - после

Для серийной нагревательной трубчатой печи ЗР2 проведено сравнительное исследование двух возможных схем отопления: 1) от вертикальных рядов излучающих беспламенных горелок и 2) от подощелевых настильных факельных горелок.

Конструктивная схема печи ЗР2 представлена на рис. 5.1в.

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

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

Рис. 5.7. Зависимость интегрального тепловосприятия экрана от длины факела	Рис. 5.8. Зависимость коэффициента неравномерности тепловосприятия экрана от длины факела

Результаты исследования закономерностей теплоотдачи настильных факелов применены для совершенствования режимов работа серийной трубчатой печи каталитической паровой конверсии природного газа ППР-600 (печь первичного риформинга установки производства аммиака производительностью 600 т/сутки).

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

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

Рис. 5.9. Изменение параметров вдоль реакционных труб по ходу движения потока:

q - теплонапряженность; tст - температура стенки трубы; tпот - тепмпература потока; x - степень конверсии сырья; P - давление

-- - параметры до изменения режима; ---- - то же после изменения;

Рекомендовано увеличить длину факелов, для чего проведено сравнительное расчетное исследование теплообмена в топочной камере печи при двух вариантах относительного количества вторичного воздуха на горение: номинальном 40% и измененном 100%, соответственно. Из рис.а5.9 видно, что при большем количестве вторичного воздуха (диффузионное горение) наблюдается сглаживание максимумов теплонапряжения и их смещение в сторону начального и среднего участков реакционных труб. Этому способствуют и наклонные настильные стены. При этом реализуются более мягкие термические условия работы реакционных труб, более благоприятные условия протекания химических процессов и как результат - рост степени конверсии сырья на ~10%. На установках производства аммиака в трубчатых печах проводится первичная конверсия сырья до 60-70%. Поэтому производительность печи по сырью может быть увеличина до достижения требуемой степени конверсии сырья соответственно на 10%.

Режимы теплообмена в нагревательных трубчатых печах были исследованы на примере печи газовой промышленности ББ1 коробчатого типа и печи нефтеперерабатывающей промышленности ЦД4 цилиндрического типа (рис.а5.10).

Возможности увеличения производительности печи были исследованы на примере серийной нагревательной трубчатой печи ББ1. Конструктивная схема печи представлена на рис. 5.10а. Печь узкокамерная с настенными экранами и фронтальными рядами панельных горелок.

Разработана зональная геометрическая модель топки печи и проведено математическое моделирование теплообмена.

а	б

Рис. 5.10. Схемы нагревательных трубчатых печей

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

Рис. 5.11. Средняя теплонапряженность труб перевального и подового экранов: № 1-12 - номера труб на схеме печи; -- - расчет; • - эксперимент	Рис. 5.12. Локальная теплонапряженность экранных труб по окружности: № 1, 5, 8, 11 - номера труб на схеме печи; линии - расчет; • - эксперимент (труба № 5)

Для печи ББ1 Норильской ГРС имеются подробные данные натурного обследования, включая локальные характеристики теплообмена по окружности труб. Обследование выполнено институтом ВНИПИгаздобыча.

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

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

Рис. 5.13. Распределение теплонапряженности по высоте экрана при тепловой нагрузке на секцию печи 10 МВт: -- - расход первичного воздуха 100%; сдняя теплонапряженность 42,44 кВт/м2 ------ - расход первичного воздуха 60%, вторичного воздуха 40%; среднее теплонапряженность 41,15 кВт/м2	Рис. 5.14. Распределение теплонапряженности по высоте экрана при тепловой нагрузке на секцию печи 12 МВт: -- - расход первичного воздуха 100%; среднее теплонапряжение 54,44 кВт/м2 ------ - расход первичного воздуха 60%, вторичного воздуха 40%; сднее теплонапряжение 50,95 кВт/м2

Конструктивная схема печи показана на рис. 5.10б. Печь ЦД4 используется в нефтеперерабатывающей промышленности и состоит из четырех секций. Она имеет центральный призматический рассекатель-распределитель с четырьмя вогнутыми гранями для настильного сжигания газообразного топлива в диффузионном факеле. Через отверстия в рассекателе на уровне 1/2 высоты топки производится подвод вторичного воздуха. Первичный воздух поступает через туннели диффузионных горелок за счет естественной тяги.

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

При необходимости производительность печи может быть увеличена на 20% с сохранением первоначальной теплонапряженности печных труб (рис.а5.14).

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

В настоящее время преимущественное распространение получили печи с канальным обогревом и рециркуляцией продуктов сгорания в топку печи для создания нужной температуры теплоносителя ~500С (Рис 6.1). Каналы обычно располагаются сверху и снизу печного конвейера в виде плоской или трубчатой поверхности нагрева. Высота надконвейерного пространства - до 50 см. Температура греющих каналов до 300-400С. Среда пекарной камеры - паровоздушная смесь с высоким содержанием водяного пара (60-70% об.) в области греющих каналов. Температура среды - до 250С. Газовая динамика выражена слабо. В этих условиях преобладающим является радиационный теплообмен - до 80%.

Рис. 6.1. Конструктивная схема хлебопекарной печи с канальным обогревом:

1,2,3,4 - греющие трубные каналы с соответствующими порядковыми номерами; 5 - топка; 6 - посадочно-выгрузочное окно; 7 - цепной люлечный конвейер; 8 - увлажняющее устройство

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

Зональное описание внешнего теплообмена в пекарной камере дается системой уравнений (2.1).

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

• для теплового потока ;                (6.1)
• для температуры  ,                                                (6.2)

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

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

уравнение переноса теплоты в твердой фазе (скелете хлеба)

               (6.3)

уравнение переноса жидкой влаги в твердой фазе (скелете хлеба)

                               (6.4)

уравнение переноса теплоты в газовой фазе (порах хлеба)

                                       (5.5)

уравнение фильтрации пара в газовой фазе (порах хлеба)

                                                                       (6.6)

;

уравнение состояния пара

.                                                                                (6.7)

Расчетная область задачи:

X=x/L; L=f(Fo).                                (6.8)

Начальные условия: τ=0: t=tнач ; u=uнач .                                        (6.9)

Граничные условия:                (6.10)

1 период выпечки  X=1: q=qконд ; X=0: q=qт ;

2, 3 периоды выпечки X=1: q=qконв + qлуч ; X=0: q=qт ;

где а2 - влагопроводность материала скелета выпекаемого изделия; с - теплоемкость; Fo - число Фурье; L - высота изделия; qт , qконв , qлуч , qконд - плотность теплового потока, обусловленная теплопроводностью, конвекцией, излучением и конденсацией пара на поверхности выпекаемого изделия; u - влажность материала; uг - гигроскопическая влажность; X=x/L - безразмерная координата; λ - коэффициент теплопроводности; ρ - плотность; τ - время; τвып - время выпечки; φ - относительная влажность среды пекарной камеры; εа - пористость; γ - удельная поверхность; αконв - коэффициент теплоотдачи в порах; j - поток вещества; ΔМфаз - источник массы вещества при фазовом переходе; Qфаз=ΔМфазrρ1,2 - источник теплоты при фазовом переходе ; r - теплота парообразования. Индексы: 0, 1, 2 - скелет хлеба, пар и жидкость соответственно; нач. - начальное значение; w - поверхность.

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

Исследован теплообмен в современной печи Р3-ХПА с тупиковой технологической камерой. Конструкция печи показана на рис. 6.1.

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

Рис. 6.2. Изменение температуры: tтр1 , tтр3 - 1 и 3 трубные каналы; tс - среда пекарной камеры; tп - поверхность выпекаемых изделий	Рис. 6.3. Изменение плотности тепловых потоков на поверхности выпекаемых изделий: qтрр,пад - падающий лучистый от трубных каналов; qср,пад - падающий лучистый от среды пекарной камеры; q∑р,рез - результирующий лучистый; qк - конвективный; q∑ - суммарный

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

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

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

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

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

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

Зональная геометрическая модель факела представляет собой усеченный конус с углом раскрытия, зависящим от условий истечения (для свободного прямоструйного факела 12,5), и высотой равной длине факела. По длине факел разбивается на зоны.

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

Длина произвольного свободного факела определяется на основе длины свободного диффузионного факела с поправочными множителями:

,                                                                (7.1)

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

Длина свободного диффузионного факела рассчитывается по выражению:

;                                                        (7.2)

где - приведенное число Архимеда; dг - диаметр сопла горелки; Lсх - стехиометрическая длина факела, т.е та его часть, на которой подсасывается теоретически необходимое количество воздуха; g - ускорение свободного падения; ρг и ρв - плотности газа и воздуха на срезе горелки; uг - скорость газа на срезе горелки.

Зональное описание сложного теплообмена в факеле имеет вид:

система уравнений зонального теплообмена

;                                (7.3)

где N - число факельных зон.

коэффициенты радиационного обмена

;                                                (7.4)

где ψij,k - ОУК излучения между факельными зонами; δij Цдельта символ Кронекера.

коэффициенты конвективного обмена

;                                                                        (7.5)

;                                                                                (7.6)

где Vпсj-1 и Vпсj - объемные расходы продуктов сгорания в сечениях факела для предыдущей и текущей факельной зоны; ρ0j-1 , ρ0j и сj-1 , cj - плотности и теплоемкости продуктов сгорания в предыдущей и текущей факельной зоне.

свободный член

;                                                                (7.7)

где Vг - расход горячего газа через горелку; Qнс - теплота сгорания газа; xj-1 и xjа - степень выгорания горючего газа в предыдущей и текущей факельных зонах.

степень черноты факельных зон

.                                                                (7.8)

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

Суммарный поток падающего излучения от всех зон факела к площадке поверхности грунта dFМ :

;                (7.9)

где εфi - степень черноты i-ой факельной зоны; Fiа - площадь поверхности i-ой факельной зоны; Ti - абсолютная температура i-ой факельной зоны; σ0 - постоянная Стефана-Больцмана; φiM - угловой коэффициент с i-ой факельной зоны на элементарную площадку dFМ .

Для вычисления угловых коэффициентов между факельной зоной и элементарной площадкой грунта в диссертации предложен численный метод, который предполагает разбиение поверхности факельной зоны на малые площадки ΔFф (рис. 7.1), вычисление элементарных угловых коэффициентов и их суммирование для определения общего углового коэффициента с факельной зоны. Для вертикальных факельных установок учитывается отклонение факела от вертикального положения за счет воздействия ветра.

Угловой коэффициент с факельной зоны на площадку dFM рассчитывается по формуле:

.                                        (7.10)

где cosφΔF и cosφM - косинусы углов между линией r, соединяющей центры площадок, и нормалями к площадке факела ΔF и элементарной площадке поверхности грунта dFM

На основании аналитической геометрии в пространстве получены выражения для cosφΔF , cosφM и расстояния r для чего осуществлен перевод определяющих геометрических параметров конической поверхности факельной зоны из локальной системы координат {o′; x′, y′, z′} в глобальную {o;x,y,z} , связанную с поверхностью грунта (рис 7.1).

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

Полученные результаты расчета в виде распределения локальной плотности потока падающего излучения на поверхности грунта для двух типов факельных установок приведены на рис. 7.2 и 7.3.

Рис. 7.2. Распределение локальной плотности тепловых потоков q,акВт/м2 для вертикальной факельной установки (изолинии): ⊗ - ствол; ↑ - направление ветра	Рис. 7.3. Распределение локальной плотности тепловых потоков q,акВт/м2 для горизонтальной факельной установки (изолинии и значения в узлах сетки): ⇒ - горизонтальный факельный ствол

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

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

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

Выводы по работе

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

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

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

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

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

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

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

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

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

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

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

Публикации в изданиях, рекомендованных  ВАК Минобрнауки России

1. Кулешов, О.Ю. Зональная математическая модель и методика расчета сложного теплообмена в промышленных хлебопекарных печах / О.Ю. Кулешов, В.М. Седелкин // Вестник Саратовского государственного технического университета. - 2010. - № 3(46). - С. 136-143.
2. Кулешов, О.Ю. Исследование режимов работы трубчатых печей с настильными факелами на основе математического моделирования / О.Ю. Кулешов, В.М. Седелкин // Промышленная энергетика. - 2011. - №4. - С. 33-36.
3. Кулешов, О.Ю. Зональная математическая модель и методика расчета сопряженного теплообмена в радиантной секции трубчатых печей / О.Ю. Кулешов, В.М. Седелкин // Вестник Саратовского государственного технического университета. - 2011.- №1 (48). - С. 181-187.
4. Кулешов, О.Ю. Методика расчета сопряженного теплообмена в технологических трубчатых печах в рамках зонального подхода / О.Ю. Кулешов, В.М. Седелкин // Известия вузов. Проблемы энергетики. - 2011.- №5-6. - С. 47-54.
5. Кулешов, О.Ю. Метод расчета теплообмена в рабочей зоне открытых факельных установок / О.Ю. Кулешов, В.М. Седелкин // Химическое и нефтегазовое машиностроение. - 2011. - №9. - С. 12-14.
6. Кулешов, О.Ю. Исследование локального результирующего теплообмена в экранированных топках трубчатых печей на основе математического моделирования / О.Ю. Кулешов, В.М. Седелкин // Промышленная энергетика. - 2011. - №11. - С. 34-37.
7. Кулешов, О.Ю. Коррекционный зональный метод расчета радиационного и сложного теплообмена в высокотемпературных установках / О.Ю. Кулешов, В.М. Седелкин // Вестник Саратовского государственного технического университета. - 2011. - №4 (60). - С. 157-161.
8. Кулешов, О.Ю. Анализ характеристик сложного теплообмена в промышленных хлебопекарных печах в зависимости от конструктивных и режимных параметров / О.Ю. Кулешов, В.М. Седелкин // Вестник Саратовского государственного технического университета. - 2011. - №4 (60). - С. 161-165.
9. Кулешов, О.Ю. Метод расчета локальных характеристик сложного теплообмена в экранированных топках в рамках зонального подхода / О.Ю. Кулешов, В.М. Седелкин // Научный журнал Труды Академэнерго. - 2012. - №1. - С. 32-43.
10. Кулешов, О.Ю. Математическое моделирование локального результирующего теплообмена в экранированных топках / О.Ю. Кулешов, В.М. Седелкин // Тепловые процессы в технике. - 2012. - Т.4, №3. - С. 118-124.
11. Кулешов, О.Ю. Анализ эффективности применения различных систем сжигания газообразного топлива в реакционных трубчатых печах / О.Ю. Кулешов, В.М. Седелкин // Химическое и нефтегазовое машиностроение. - 2012. - №4. - С. 9-12.
12. Кулешов, О.Ю. Расчетный анализ локальной теплонапряженности экранных труб в реакционных трубчатых печах / О.Ю. Кулешов, В.М. Седелкин // Химическое и нефтегазовое машиностроение. - 2012. - №5. - С. 15-18.
13. Кулешов, О.Ю. Уточнение зонального метода расчета сложного теплообмена в огнетехнических установках в части учета радиационных свойств продуктов сгорания / О.Ю. Кулешов, В.М. Седелкин // Известия вузов. Проблемы энергетики. - 2012. - №3-4. - С. 20-27.
14. Кулешов, О.Ю. Новый подход к анализу тепловых режимов промышленных печей с использованием метода коррекции зональных оптико-геометрических характеристик излучения / О.Ю. Кулешов, В.М. Седелкин // Промышленная энергетика. - 2012. - №6. - С. 39-43.
15. Кулешов, О.Ю. Математическое моделирование процессов горения и теплообмена в открытых факельных установках / О.Ю. Кулешов, В.М. Седелкин // Вестник Самарского государственного технического университета. Серия Технические науки. - 2012. - №1 (33). - С. 196-202.
16. Кулешов, О.Ю. Исследование режимов радиационно-конвективного теплообмена в промышленных хлебопекарных печах на основе математического моделирования / О.Ю. Кулешов, В.М. Седелкин // Промышленная энергетика. - 2012. - №7. - С. 28-33.
17. Кулешов, О.Ю. Сравнительный анализ результирующего теплообмена в реакционных трубчатых печах при различных схемах отопления / О.Ю. Кулешов, В.М. Седелкин // Промышленная энергетика. - 2012. - №8. - С. 26-31.

Публикации в зарубежных научных периодических изданиях:

18. Кулешов, О.Ю. Методика расчета теплового излучения в зоне действия открытых факелов / В.М. Седелкин, А.В. Паимов, О.А. Толоконникова, О.Ю. Кулешов //Инженерно-физический журнал (НАН Беларуси). - 1993. - Т.64, №3. - С. 297-299.
19. Kuleshov, O.Yu. The Zone-Element Method in Application to the Combined Heat Transfer Problems/ O.Yu. Kuleshov, V.M. Sedelkin // Heat transfer research. - 2002. - vol.33, №7,8. - С. 496-501. (ISSN: 1064-2285 Print)

Публикации в других изданиях:

20. Математические модели и пакет программ для расчетов сложного теплообмена в огнетехнических установках / В.М. Седелкин, А.В. Паимов, М.С. Угольников, О.Ю. Кулешов, О.П. Полонская // Седьмая Всесоюз. конф. по радиационному теплообмену: Материалы конф., Ташкент, 21-23 окт. 1991. - Ташкент, 1991. - С. 109-110.
21. Кулешов, О.Ю. Разработка математического обеспечения САПР огнетехнических установок / В.М. Седелкин, О.Ю. Кулешов, А.В. Паимов // Теплофиз. проблемы пром. производства: Материалы Междунар. совещания семинара. - Тамбов: ТТИ, 1992. - С. 102.
22. Определение локальных характеристик радиационного теплообмена в факельных системах / А.В. Паимов, В.М. Седелкин, О.Ю. Кулешов, М.С. Угольников // Тепломассообмен- ММФ-92. Радиационный и комбинированный теплообмен. Т.2. - Минск: ИТМО, 1992. - С. 43-45.
23. Кулешов, О.Ю. Методика расчета на ПЭВМ радиационного теплообмена в открытых факельных системах / В.М. Седелкин, О.Ю. Кулешов // Радиационный и сложный теплообмен: Тр. 1-ой Российской нац. конф. по теплообмену. - М.: Изд-во МЭИ, 1994. - Т.9. - С. 196-199.
24. Кулешов, О.Ю. Математическая модель и методика расчета сопряженного теплообмена в трубчатых реакторах / В.М. Седелкин, О.Ю. Кулешов // Радиационный и сложный теплообмен: Тр. 1-ой Российской нац. конф. по теплообмену. - М.: Изд-во МЭИ, 1994. - Т.9. - С. 190-195.
25. Кулешов, О.Ю. Разработка математической модели сопряженного теплообмена в трубчатых реакторах / В.М. Седелкин, О.Ю. Кулешов // Материалы 1-ой Науч. техн. конф. (апрель 1993 г.): Сборник / Сарат. гос. техн. ун-т. Технол. ин-т. - Энгельс, 1995. - С. 226-243. - Деп. ВИНИТИ 10.03.95, № 660-В 95.
26. Kuleshov, O.Yu. Mathematical support CAD/CAM of process of gasous waste open jet-combustion / V.M. Sedelkin, O.Yu. Kuleshov // 12-th Intern. Cong. of Chem. and Process Eng. CHISAТ96: Summaries topics 8. (Praha, august 25-30 1996). - Praha, 1996. - P. 120.
27. Kuleshov, O.Yu. Mathematical support CAD/CAM of process of gasous waste open jet-combustion / V.M. Sedelkin, O.Yu. Kuleshov // 12-th Intern. Cong. of Chem. and Process Eng. CHISAТ96: Full text of the paper, Praha, Czech. Republik, 25-30 August, 1996. - Praha, 1996. - № P.9.64.-4 p.
28. Kuleshov, O.Yu. Mathematical Model of Conjugated heat transfer in Tubular reactors with outside heating by combustion products / V.M. Sedelkin, O.Yu. Kuleshov // 12-th Intern. Cong. of Chem. and Process Eng. CHISAТ96: Summaries topics 5. (Praha, august 25-30 1996). - Praha, 1996. - P. 61.
29. Kuleshov, O.Yu. Mathematical Model of Conjugated heat transfer in Tubular reactors with outside heating by combustion products / V.M. Sedelkin, O.Yu. Kuleshov // 12-th Intern. Cong. of Chem. and Process Eng. CHISAТ96: Full text of the paper, Praha, Czech. Republik, 25-30 August, 1996. - Praha, 1996. - № P.1.18.-4 p.
30. Kuleshov, O.Yu. The mathematical modeling of processes into tubular pyrolysis reactor with account of inner radiative heat transfer effect / O.Yu. Kuleshov, V.M. Sedelkin // The First European Congress on Chemical Engineering - ECCE-1: Proceedings. - Milano, 1997. - Vol. 3. - P. 1815-1817.
31. Kuleshov, O.Yu. Method for numerical investigation of complex problem of burning and outer heat radiating of open gasous jets / V.M. Sedelkin, O.Yu. Kuleshov // The First European Congress on Chemical Engineering - ECCE-1: Proceedings. - Milano, 1997. - Vol. 3. - P. 1819-1821.
32. Kuleshov, O.Yu. New zone-node method and technique for numerical solving of conjugate heat exchange problem in tubular furnace / V.M. Sedelkin, O.Yu. Kuleshov // The First European Congress on Chemical Engineering - ECCE-1: Proceedings. - Milano, 1997. - Vol. 3. - P. 1897-1899.
33. Кулешов, О.Ю. Современные математические модели и методы расчета теплообмена в трубчатых печах / В.М. Седелкин, О.Ю. Кулешов // Проблемы эффективного использования энергоносителей и низкосортных топлив в промышленности: Материалы Междунар. конф., Саратов, 24-25 сент. 1998. - Саратов: СГТУ, 1998. - С. 40-44.
34. Кулешов, О.Ю. Математическая модель и методика расчета газодинамики и теплообмена в хлебопекарных печах / В.М. Седелкин, О.Ю. Кулешов, Т.В. Кузьмичева // Тр. 2-ой Российской нац. конф. по теплообмену, М., 26-30 октября 1998 г. - Т.6. - М., 1998. - С. 204-206.
35. Кулешов, О.Ю. Численное моделирование радиационно-конвективного теплообмена в потоке термически разлагающихся углеводородов / О.Ю. Кулешов, В.М. Седелкин // Тр. 2-ой Российской нац. конф. по теплообмену, М., 26-30 октября 1998 г. - Т.6. - М., 1998. - С. 308-311.
36. Kuleshov, O.Yu. Modern zoning analysis of applied problems of complex heat transfer / O.Yu. Kuleshov, V.M. Sedelkin, T.V. Kuzmicheva // 13-th International Congress of Chemical and Process Engineering CHISAТ98: Summaries 8. Heat transfer processes and equipment, general topics, Czech Republic, Praha, 23-28 August 1998. - Praha, 1998. - P.23.
37. Kuleshov, O.Yu. Mathematical model and calculation technique of conjugate heat-mass transfer in technolgical furnaces / V.M. Sedelkin, O.Yu. Kuleshov, T.V. Kuzmicheva // 13-th International Congress of Chemical and Process Engineering CHISAТ98: Summaries 8. Heat transfer processes and equipment, general topics, Czech Republic, Praha, 23-28 August 1998. - Praha, 1998. - P.21.
38. Kuleshov, O.Yu. Mathematical model of tube-side carbon deposit dynamics in pyrolysis furnace / O.Yu. Kuleshov, V.M. Sedelkin // 13-th International Congress of Chemical and Process Engineering CHISAТ98: Summaries 8. Heat transfer processes and equipment, general topics, Czech Republic, Praha, 23-28 August 1998. - Praha, 1998. - P.22.
39. Кулешов, О.Ю. Современное состояние и перспективы развития методов теплового расчета трубчатых печей газовой и нефтехимической промышленности / О.Ю. Кулешов. - Технол. ин-т Сарат. гос. техн. ун-та. - Энгельс, 1999. - 16 с. - Деп. в ВИНИТИ 29.06.99, № 2085-В 99 // Б.у. Деп. науч. работы. - 1999. - № 8 (331). - б/о 366.
40. Кулешов, О.Ю. Математическое моделирование процессов в трубчатом реакторе пиролиза углеводородов с учетом радиационной составляющей внутреннего теплообмена / О.Ю. Кулешов. - Технол. ин-т Сарат. гос. техн. ун-та. - Энгельс, 1999. - 11 с. - Деп. в ВИНИТИ 29.06.99, №2084-В 99 // Б.у. Деп. науч. работы. - 1999. - № 8 (331). - б/о 358.
41. Кулешов, О.Ю. Разработка физико-математической модели сопряженного теплообмена в промышленных хлебопекарных печах / О.Ю. Кулешов. - Технол. ин-т Сарат. гос. техн. ун-та. - Энгельс, 1999. - 8 с. - Деп. в ВИНИТИ 29.06.99, № 2086-В 99 // Б.у. Деп. науч. работы. - 1999. - № 8 (331). - б/о 384.
42. Кулешов, О.Ю. Методы математического моделирования сложного и сопряженного теплообмена в аппаратах нефтехимической технологии / О.Ю. Кулешов, В.М. Седелкин, Ю.Я. Печенегов // Повышение эффективности тепломассообменных процессов и систем: Матер. 2-ой Междунар. научн. конф., Вологда, 19-22 апреля 2000 г. - Вологда, 2000. - Ч.2. - С. 174-176.
43. Кулешов, О.Ю. Зонально-элементный метод решения прикладных задач сложного теплообмена / О.Ю. Кулешов, В.М. Седелкин // Тепломассообмен ММФ-2000: Тр. 4-го Минского Междунар. форума, Минск, 22-26 мая 2000 г. - Минск, 2000. - Т.2. - С.36-40.
44. Кулешов, О.Ю. Математическая модель и методика расчета сопряженного теплообмена в промышленных хлебопекарных печах / О.Ю. Кулешов, В.М. Седелкин, Т.В. Кузьмичева, В.А. Шныпко // Тепломассообмен ММФ-2000: Тр. 4-го Минского Междунар. форума, Минск, 22-26 мая 2000 г. - Минск, 2000. - Т.11. - С. 125-129.
45. Kuleshov, O.Yu. Combined zone-elements method for solving complex heat transfer problems in highly irregular calculation areas/ O.Yu. Kuleshov, V.M. Sedelkin // 14-th International Congress of Cemical and Process Engineering CHISAТ2000, Czech Republic, Praha, 27-31 August 2000. - Praha, 2000. - P.44.
46. Кулешов, О.Ю. Двухфазная математическая модель теплопереноса в трубчатом реакторе каталитической конверсии углеводородов / О.Ю. Кулешов, В.М. Седелкин // Тепло- и массообмен в химической технологии: Материалы Всеросс. науч. конф., Казань, 19-20 декабря 2000 г. - Казань, 2000. - С.124-125.
47. Кулешов, О.Ю. К расчету теплоотдачи излучающего газового потока / О.Ю. Кулешов, В.М. Седелкин // Тепло- и массообмен в химической технологии: Материалы Всеросс. науч. конф., Казань, 19-20 декабря 2000 г. - Казань, 2000. - С.126-127.
48. Кулешов, О.Ю. Методы расчета сопряженного теплообмена в высокотемпературных трубчатых реакторах / О.Ю. Кулешов, В.М. Седелкин // Тепло- и массообмен в химической технологии: Материалы Всеросс. науч. конф., Казань, 19-20 декабря 2000 г. - Казань, 2000. - С.128-129.
49. Кулешов, О.Ю. Оптимальное управление теплообменом в печах с дифференцированным подводом теплоты / О.Ю. Кулешов, В.М. Седелкин // Прогрессивные процессы и оборудование металлургического производства: Материалы 2-ой Всеросс. науч. техн. конф. - Череповец: ЧГУ, 2001. - С. 33-35.
50. Кулешов, О.Ю. Повышение эффективности факельного сжигания углеводородных топлив в печах / О.Ю. Кулешов, В.М. Седелкин // Энергосбережение в теплоэнергетических системах: Материалы Междунар. науч. техн. конф., Вологда, 24-26 апреля 2001 г. - Вологда, 2001. - С.21-24.
51. Кулешов, О.Ю. Рациональное энергоиспользование в технологических печах / О.Ю. Кулешов, В.М. Седелкин // Энергосбережение в теплоэнергетических системах: Материалы Междунар. науч. техн. конф., Вологда, 24-26 апреля 2001 г. - Вологда, 2001. - С.156-158.
52. Кулешов, О.Ю. Разработка метода контроля температуры экранной поверхности нагрева в трубчатых печах / О.Ю. Кулешов, В.М. Седелкин // Теплофизические измерения в начале XXI века: Материалы 4-ой Междунар. Теплофиз. школы, Тамбов 24-28 сентября 2001 г. - Тамбов 2001. - Ч.1. - С. 58-59.
53. Кулешов, О.Ю. Математическое обеспечение САПР технологических печей / О.Ю. Кулешов, В.М. Седелкин // Теплофизические измерения в начале XXI века: Материалы 4-ой Междунар. Теплофиз. школы, Тамбов 24-28 сентября 2001 г. - Тамбов 2001. - Ч.2. - С. 135.
54. Кулешов, О.Ю. Разработка численного метода расчета тепломассопереноса в хлебопекарных печах / О.Ю. Кулешов, В.М. Седелкин // Сборник докладов Юбилейной Междунар. науч.-практ. конф. Пищевые продукты XXI века. - М.: МГУПП, 2001. - Т.2.- С.61.
55. Кулешов, О.Ю. Математический аппарат анализа режимных характеристик технологических печей / О.Ю. Кулешов, В.М. Седелкин // Моделирование, оптимизация и интенсификация производственных процессов и систем: Матер. Междунар. науч. техн. конф., Вологда, 29-31 октября 2001 г. - Вологда, 2001. - С. 27-28.
56. Кулешов, О.Ю. Математическая модель и методика расчета локального теплообмена в экранированных топках / О.Ю. Кулешов, В.М. Седелкин, Ю.Я. Печенегов // Повышение эффективности теплообменных процессов и систем: Материалы 3-ей Междунар. науч.-техн. конф., Вологда, 20-23 мая 2002 г. - Вологда: ВГТУ, 2002. - С.27-30.
57. Кулешов, О.Ю. Разработка зонального метода расчета характеристик тепломассообмена в промышленных хлебопекарных печах / О.Ю. Кулешов, В.М. Седелкин // Современные научные и информационные технологии: Матер. науч.-метод. конф., посвящ. 25-летию Мех.-маш. ф-та ТИ СГТУ, г. Энгельс, 21 октября 2003. - Саратов, 2003. - С. 33-36.
58. Кулешов, О.Ю. Математическое моделирование сложного и сопряженного теплообмена в промышленных хлебопекарных печах / О.Ю. Кулешов, В.М. Седелкин // Материалы научно-техн. конф., посвященной 50-летию ЭТИ СГТУ: материалы, г. Энгельс, 20-21 ноября 2006 г. - Саратов: СГТУ, 2006. - С. 85-87.
59. Кулешов, О.Ю. Математическое моделирование сложного и сопряженного теплообмена в технологических камерах промышленных хлебопекарных печей / О.Ю. Кулешов, В.М. Седелкин; Энгельс. технол. ин-т Сарат. гос. техн. ун-та. - Энгельс, 2007. - 17 с.: ил. - Библиогр.: 11 назв. - Рус. - Деп. в ВИНИТИ 03.04.07, № 366-В2007.
60. Кулешов, О.Ю. Математическое моделирование сложного и сопряженного теплообмена в трубчатых печах газовой и нефтехимической промышленности / О.Ю. Кулешов, В.М. Седелкин; Энгельс. технол. ин-т Сарат. гос. техн. ун-та. - Энгельс, 2007. - 17 с.: ил. - Библиогр.: 18 назв. - Рус. - Деп. в ВИНИТИ 03.04.07, № 367-В2007.
61. Кулешов, О.Ю. К расчету локальных характеристик радиационного и сложного теплообмена в рамках зонального метода / О.Ю. Кулешов, В.М. Седелкин // Теплофизика в энергосбережении и управлении качеством. Часть 1: материалы Шестой Международной теплофизической школы, г. Тамбов, 1-6 октября 2007 г. - Тамбов: Изд-во ТГТУ, 2007. - C. 133-134.
62. Кулешов, О.Ю. Двухфазная математическая модель тепло-и влагопереноса в выпекаемом тестовом изделии / О.Ю. Кулешов, В.М. Седелкин //Теплофизика в энергосбережении и управлении качеством. Часть 1: материалы Шестой Международной теплофизической школы, г. Тамбов, 1-6 октября 2007 г. - Тамбов: Изд-во ТГТУ, 2007. - C. 135-137.
63. Кулешов, О.Ю. Разработка метода численного расчета локальных характеристик сложного теплообмена в трубчатых печах / О.Ю. Кулешов, В.М. Седелкин // Математические методы в технике и технологиях - ММТТ-21: сб. трудов XXI Международной научной конференции. Том 5. - Саратов: Сарат. гос. техн. ун-т, 2008. - С. 96-98.
64. Кулешов, О.Ю. Математическое моделирование открытого факельного сжигания сбросных газов нефтегазовой промышленности / О.Ю. Кулешов, В.М. Седелкин // Математические методы в технике и технологиях - ММТТ-21: сб. трудов XXI Международной научной конференции. Том 5. - Саратов: Сарат. гос. техн. ун-т, 2008. - С. 185-187.
65. Кулешов, О.Ю. Математическая модель процесса выпечки хлеба в промышленных печах / О.Ю. Кулешов, В.М. Седелкин // Математические методы в технике и технологиях - ММТТ-21: сб. трудов XXI Международной научной конференции. Том 5. - Саратов: Сарат. гос. техн. ун-т, 2008. - С. 192-195.
66. Кулешов, О.Ю. Математические модели и методы расчета сопряженного теплообмена в промышленных огнетехнических установках / О.Ю. Кулешов, В.М. Седелкин // Синтез инноваций: направления и перспективы: материалы научно-практической конференции. - Саратов: Сарат. гос. техн. ун-т, 2009. - С. 52-53.
67. Кулешов, О.Ю. Совершенствование тепловых режимов промышленных печей на базе высокоэффективной зональной методики расчета сложного теплообмена / О.Ю. Кулешов, В.М. Седелкин // Проблемы энерго- и ресурсосбережения: сборник научных трудов. - Саратов: Сарат. гос. техн. ун-т, 2010. - С. 197-205.
68. Кулешов О.Ю. Повышение вычислительной эффективности зонального метода расчета сложного теплообмена на основе методики коррекции базовых оптико-геометрических характеристик излучения / О.Ю. Кулешов, В.М. Седелкин // Труды пятой Российской национальной конференции по теплообмену. - Т.6. - М.: Издат. дом МЭИ, 2010. - С. 227-230.
69. Кулешов, О.Ю. Коррекционный зональный метод расчета радиационного и сложного теплообмена в энергетических установках / О.Ю. Кулешов, В.М. Седелкин // Тепломассообмен ММФ-XIV: Тр. 14-го Минского Междунар. форума. - Минск, 2012. - Т.2. - С. 40-44.
70. Кулешов, О.Ю. Математическая модель и методика численного расчета сопряженного теплообмена в трубчатых печах - реакторах / О.Ю. Кулешов, В.М. Седелкин // Математические методы в технике и технологиях - ММТТ-25: Сборник трудов XXV Международной научной конференции, Волгоград, 2012 г. - Т.2. - Волгоград: ВГТУ, 2012; Харьков: ХПИ, 2012. - С. 68-71.

Авторефераты по всем темам  >>  Авторефераты по техническим специальностям