Исследование излучательной способности конструкционных материалов и пылегазовых сред применительно к расчету теплообмена в котлах-утилизаторах Таймаров, Михаил Александрович

Введение к работе

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

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

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

Опытные данные ' по оптическим константам' вецеств частиц и Фракционному составу служат основой для теоретических расчетов коэффициентов ослабления излучения пылевыми потоками в котлах-утили-

-.? т

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

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

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

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

1.Создание стендовых установок и программ для обработки данных на ЭВМ по исследованию интегральной и спектральной излучательной способности конструкционны:: материалов.

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

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

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

5.Создание экспериментальной методики и установки для изуче
ния поглощательной и излучательной способности пылевых потоков в
стендовых условиях. Проведение экспериментов, обработка результа
тов с применением ЭВМ. ₀

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

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

8. Применение зонального метода к расчету радиационного теплообмена в энерготехнологических котлоагрегатах. Расчеты на ЭВМ, обобщение, анализ результатов. Выработка рекомендаций по оптимизации конструкций и интенсификации теплообмена котлов-утилизаторов.

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

Выполнены исследования нормальной интегральной и спектральной излучательной способности огнеупоров при нагреве их в вакууме и в воздушной среда в диапазоне 600-21 СО К. Получены формулы для расчета интегральной излучательной способности огнеупоров в зависимости от химического состава, температуры при нагреве в Еакууме и в воздушной среде. Для диапазона длин волн 1-6 мкм предложена фор-

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

Путем экспериментов в стендовых условиях установлен характер влияния на интегральную поглощательную способность пылевых потоков среднего диаметра, температуры и концентрации частиц твердых дисперсных фаз из энерготехнологическнх агрегатов: УККС-6/40, КУ-60, КС-450 ВТКУ. РКК-20/40, КУ-125. печь КС, БКЗ-50-39У.

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

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

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

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

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

Результаты работы использованы при тепловых расчетах новых конструкций котлов-утилизаторов на Белгородском заводе "Энерго-маш". в разработках высокотемпературных энергетических установок во Всесоюзном институте огнеупоров (г.Санкт-Петербург), для проектирования энерготехнологических агрегатов в Нефтекамском ЦБПО"Баш-нефть", в расчетах радиационного теплообмена при реконструкции печи обжига керамики на Казанском комбинате строительных материалов, при режимно- и пусконаладочных работах на энергетических котлах Казанской ТЭЦ-2, при проектировании и эксплуатации энррготехноло-гических агрегатов на Бугульминском механическом заводе, в режимно-и пусконаладочных работах в Нефтекамском П0"Иско;к" на котлах типа ДКВР. Использование результатов работы подтверждено документально актами и протоколами внедрения.

Апробация работы. Основное содержание работы опубликовано в статьях автора в журналах "Теплоэнергетика", "Измерительная техника". Инженерно-физический журнал, "Промышленная теплотехника", "Энергомашиностроение", Известия вузов "Энергетика", в отраслевых и межвузовских сборниках. Результаты работы докладывались на: отчетных научно-технических конференциях КХТИ им. СМ. Кироза (1984-1990 г.г.). Седьмой Всесоюзной конференции по тепломассообмену (г.Минск, 1984 г.), Е'ыездном заседании Совета ГКНТ "Надежность котельных поверхностей нагрева и актуальные вопросы теплообмена и гидравлики (г. Подольск, 1984 г.), Научно-техническом совещании "Состояние и перспективы сокращения вредных выбросов предприятиями промышленности в атмосферу" (г. Казань, 1985 г.). Четвертой Всесоюзной конференции "Влияние минеральной части энергетических топлив на условия работы паровых котлов" (г. Таллин. 1986 г.), на заседании секции "Теплообмен излучением" ГКНТ СССР по проблеме "массо- и теплоперенос в технологических процессах" (г. Уфа, 1986 г.)., Пятой Всесоюзной научной конференции "Математическое моделирование сложных химико-технологических систем" (г. Казань, 1988 г.), Республиканской'научно-практической конференции "Пути поЕыае-ния надежности и эффективности теплоснабжения промышленных предприятий" f г. Казань, 1989 г.).

Автором получено 2 авторских свидетельства на изобретения.

Объем работа.Диссертация состоит из введения, четырех глав, заключения (объемом 294 стр.),приложения(53 стр.),списка использованной' литературы из 333 наименований и содержит:65 рис.,53 табл. СОДЕРЖАНИЕ РАБОТЫ

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

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

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

Е(9) - E(AS₃ - AS,)/(AS₂ - AS,), (1)

где Є - угол между нормалью к поверхности образца и направлением излучения:*:⁰- эффективная степень черноты АЧ^ДЗз.ДБ^ДЗг -разности сигналов приемника при регистрациях плотностей излучений от образца,- окружающего фона и АЧТ при открытом и закрытом экране.

Спектральная степень черноты материалов исследовалась на установке, разработанной на базе спектрофотометра ИКС-14 с размещением АЧТ в канале сравнения. Выражение для определения спектральной степени черноты материала в окончательной записи имеет вид

_Х(Т) D_xE⁰x(T⁰)(exp(c₂/n) - 1)/{ехр(с₂/ХТ) - 1). (2) где D_x- пропускание фотометрического клина; Е⁰х,?-спектраль-ная эффективная степень черноты и температура АЧТ; с₂Д-2-я конс-*-В исследовании излучательной способности материалов участвовали: д.т.н. Гарифуллин Ф.А., к.т.н. Русев К.А., к.т.н. Давлетбаева Д.3.

танта излучения Планка и длина волны, мкм; Т-температура образца.

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

t = e'⁰((S₆-S4HS2-Si)-(S5-S4)(S4-S,))/(S5-S4)(Sa-S₄). (3±

где Бб^-сигналы приемника при регистрациях плотностей излучений от образца, нагреваемого пламенем горелки, и самого пламени, mB;S₂,Si-сигналы приемника при регистрациях плотностей излучений АЧТ и окружающего фона, мВ;3₅-сигнал приемника, измеряющего излучение АЧТ на фоне пламени горелки по оптической оси, мВ.

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

Графический способ обработки результатов сводился к применению градуированных функций S,MB=f(Е),Вт/м^г,где Е-излучение АЧТ.

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

F = 1/(1 + 0.000062 (Ra/Sm)²). (4)

где Ra.Sm - среднеарифметическое отклонение (мкм) и средний шаг (мм) неровностей профиля, измеряемые в опытах на двойном микроскопе нис-11 и на профилографе-профилометре П-203.

В диапазоне Т=400-900 К влияние окисной пленки на степень черноты по результатам опытов на 35 марках металлов выражается как

і = е„ + A lg(t/T„), (5)

где е„ -степень черноты в точке, прини»аемой за.начало отсчета; А -режимный коэффициент: т, т„ -текущее и начальное время время окисления (для сталей 2 участка: до к после 650 К)(табл.1).

Для окисленных на воздухе при Т=800-1070 К сталей длина волны не влияет на степень черноты, за исключением легированных сталей типа 1Х18НЮТ, имеющих в области спектра 2-4 мкм повышенные значения степени черноты.' по сравнению с диапазоном 6-12 мкм.

Излучение углеродистых сталей уже в конце первичного нагрева в области Т=900 К является изотропным. Излучение цветных металлов

- .9 -

Таблица 1.

. Вторичный

Тн.МИН А

и сплавов отклоняется от закона Ламберта.

В диапазоне длин волн 1-13 ыкм при нагреве на воздухе из числа исследованных 38 марок огнеупоров, наименьшей селективностью обладают хромитопериклазовые огнеупоры. Совместное влияние длины волны Х(мкм) в области 1-6 ыкм, температуры Т(К) и концентрации основных компонентов К (окислоз Сг, Fe, Si, Hg, Al, Ca в %) на спектральную степень черноты огнеупоров ех записывается формулой _х=1+(ехр(-1.4Х ехр(0,09\^г-0,6Х)/Х_и))(0,008(Ксг2оз⁺ +Кгвгоз+Кгео)+0.003 Ksio₂⁺0.002(Kn_g0+K_Aieo3+Kcao)-l). (6) где Х„=2896/Т-длина волны (мкм) максимума излучения АЧТ при температуре образца Т (К).

Интегральную степень черноты є в диапазоне температур
600-2200 К можно описать формулой: для алюмосиликатных огнеупоров:
е = 1 - 0,0025 Каі2оз - kTⁿ . (^?)

где к=0,0002, п=1-нагрев на воздухе: к=0,008,п=о,5-нагрев в вакууме 0.133 Па. и для нагнезитохромитовых огнеупоров:

Е - 1 - 0,004(K_Mgo - Ксггоз)' - КТ . <8)

где к-0.0001-нагрев на воздухе и к=о. 00005-в вакууме 0,133 Па. Исследования 47 различных химических составов загрязняющих

отложений с поверхностей нагрева котлов-утилизаторов показали, что для них наблюдается аналогичный, как и для огнеупоров, вид эмпирической зависимости степени черноты t_x от длины волны X, температуры и концентрации окислов. Отклонения от закона Ламберта для загрязняющих отложений при 8=70-85 незначительны. Излучательная способность чистых поверхностей нагрева выые. чем у отложений на них.

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

Основные выводы по первой главе:

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

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

а)излучательная способность чистой (без окисной пленки) поверхности практически слабо изменяется в зависимости от продолжительности нагрева т. Верхняя граница области г^и42 мин. Т ^й640 К.

б)после образования тонкой окисной пленки излучательная способность резко возрастает в зависимости от продолжительности нагрева. Верхняя граница области т*69 мин. Т^к840 К.

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

3. Предварительная выдержка в средах углеводородных газов (ра
бочие среды печей пиролиза) при температурах около 1100 К в тече
нии 20 часов образцов хромоникелевых сталей и сплавов типа Х23Н32Т
вызывает, по сравнению с водяным парой и воздухом, возрастание из
лучательной способности в среднем на 10 % для диапазона температур
Т-450-900 К.

4.Наиболее сильно в зависимости от угла в изменяется направ-

- и .-

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

5.-Для спектрального диапазона 1-13 мкм степень черноты окисленных на воздухе при 920-1070 К в течении 2-5 часов сталей и сплавов:Ст.3 сп, ЭИ712. 09Г2С. сталь 20. 40Х2Н2МА, 38ХНЗМА. 12Х1МФ. ЭП182. Д16АТ. 19ГН. 35ХГСА для температур 700-900 К незначительно изменяется в зависимости от длины волны X и их излучение может приниматься серым. Степень черноты сталей 08Х18НЮТ. ІХ18Н10Т. ЭИ925. СП28. ЭП53 для Т=700-900 К наиболее резко изменяется в диапазоне Х=3-6 мкм. У этих сталей при Х=1-3 мкм спектральная степЕ..лЬ черноты на 25-75 % выше, че: степень черноты при Х=6-13 мкм. Общей тенденцией для всех исследованных- сталей является увеличение их спектральной степени черноты с ростом температуры, составляющее в диапазоне 700-900 К в среднем 0,04.

6.Пак при нагреве в воздушной среде, так и при нагреве в вакууме, для интервала температур 600-2000 К излучательная способность огнеупоров понижается с ростом температуры, с увеличением содержания в составе огнеупоров оксидов алюминия, магния, кальция. Оксиды кремния, хрома и железа вызывают повышение степени черноты и меньшую зависимость степени черноты от температуры. В области длин волн 1-7 мкм спектральная степень черноты огнеупоров на 30-70 % ниже, чем в области 8-13 мкм, в которой практически для всех марок огнеупоров спектральная степень черноты может быть принята близкой к 0,9.

7.Излучательная способность динасових и шамотных огнеупоров при нагреве в вакууме увеличивается, по сравнению с нагревом в воздушной среде, соответственно в среднем на 0,2 и на 0,1. Для корундовых (высокоглиноземистых), магнезитовых и хромомагнезитові»: огнеупоров степень черноты увеличивается при нагреве в вакууме, по сравнению с нагревом на воздухе . в среднем на 0,04-0,08. Полученные змпиричесшіе зависимости позволяют, с погрешностью 10 % (при 1600 К) для условий нагрева на воздухе и с погрешностью 13 % при нагреве в вакууме, описать 'Елпяние температуры и химического состаза на интегральную степень черноты огнеупоров.

8. В интервале температур 500-1000 1С очищенные от загрязнений

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

9.Уменьшение интегральной степени черноты поверхностей нагрева вследствие образования на них загрязнений составляет но производствам (и котлам-утилизаторам) : конверторное получение стали (0КГ-100-ЗБ) 31. мартеновское производство (КУ-125) 25%, обжиг цинковых концентратов (УККС-6/40) 22%, медеплавильное поиззодство (ТОП-35/40) 33%, переработка никелевого штейна (РКК-20/40) 12%. шлаковозгоночный процесс получения свинца (РК-75/40) 21%, сажевое производство (ПКК-75/24-150-5) 34%. сжигание сероводорода (ПКС-10/40) 6%. целлюлозно-бумажное производство (СРК-700) 21%.

10.Интегральная степень черноты наружных загрязняющих отложений поверхностей нагрева не зависит от температуры в диапазоне 500-1000 К для агрегатов ПКК-75/24-150-5. КУ-125. ТОП-35/40 (за исключением воздухоподогревателя), РК-75/40, ПКС-10/40. С ростом температуры в диапазоне 50О-1000 К наблюдается увеличение интегральной степени черноты - загрязняющих отложений агрегатов: ОКГ-ЮО-ЗБ на 0.11. УККС-6/40 на 0,06, СРК-700 на 0,10 (диапазон 450-650 К). Уменьшение интегральной степени черноты с ростом температуры имеет место для загрязняющих отложений из воздухоподогревателя котла ТОП-35/40 (на 0,06 в интервале 550-800 К) и для СРК-700 (на 0,18 в интервале 650-900 К).

Во второй главе исследуются*:излучательная и поглощательная способности потоков пылевых частиц (взятых из газоходов промышленных агрегатов) на созданной стендовой установке, эмиссионные свойства двухфазных рабочих сред по результам экспериментов на действующих энерготехнологических агрегатах с помощью оптических приборов и оборудования, разработанных при непосредственном учас-тии автора, первичные радиационные характеристики потоков твердых * -В экспериментах на котлах участвовали:д.т.н.Сагадеев В.И., к.т. н. Кочеров М.И., к.т.н. Зайцев В. А., инж.Горбатенко И.В.. ина.Романов B.C.; к.т.н.Степанов И.Е. (исследование оптических констант).

- 13 .-

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

' Установка для исследования радиационных свойств пылевых потоков представляла собой вертикальную аэродинамическую трубу замкнутого типа, в верхней части которой находился электрообогреваемый дозатор пыли, а в нижней - циклон-пылесборник. Воздух в пылевую камеру средней части подавался тягодутьевым вентилятором через нормальную диафрагму двумя потоками: основная часть поступала по осевому вводу непосредственно в камеру, а оставшаяся часть - в дозатор для обеспечения транспортировки пыли. С одной стороны пылевой камеры размещался приемник теплового излучения ППТ-131. а с другой по оптической оси - АЧТ. Из решения уравнений теплового баланса на приемнике получено, что с приемлемыми методическими погрешностями, окончательные формулы для определения поглощательноП А и излучательной е способности пылевых потоков имеют вид:

А(Т,Т) - (S, - S₂)/S,. е = t⁰S₃/S,. (9)

где Б^Бг-сигналы приемника при регистрации излучения от АЧТ без пылевого потока и с ним. мВ; Т,Т-температуры пылевого потока и АЧТ ;с-эффективная степень черноты АЧТ: Б₃-сигнал приемника при регистрации излучения от пылевого потока, мВ.

В экспериментах применялась графическая обработка результатов с построением графиков градуирования приемника по АЧТ. Эксперименты проводились при концентрациях пыли и в аэропотоке 20-2500 г/нм³ для двух геометрических толщин слоя L: 0,08 и 0,2 м.

Из предположения справедливости закона Бугера-Бера'по опытным данным для поглощательной способности A=l-exp(-kFjiL). определялся безразмерный коэффициент ослабления лучей k=-ln(i-A)/FnL. Его значения оказались зависимыми от /iL. Здесь F- удельная поверхность пыли (м^г/г), рассчитываемая по формуле: п

F * (0,138/р) I (ADj/AX^lgd + ДХі/Хі), (10)

где р-плотность вещества частиц пыли (г/см³);п - число прямолинейных участков, на которые разделена кривая дисперсного состава D(X); ДХі=ДХц.і-Хі - разность меаду максимальным и минимальным размерами частиц на 1-ом прямолинейном участке; uDi-Di»j-Di - количество частиц до размера Х_1рЖ.

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

k =(1- D(ML/260)^p)(B(X_cp.T)ⁿ * С(Т - 290)^т). (И)

где D.р, В.п.С.т - эмпирические параметры;'X_cp.-6/pF - средний по удельной поверхности F диаметр частиц: р - плотность вещества частиц пыли; Т,Т - температуры АЧТ и пылевых частиц, К.

Значения среднего диаметра частиц Х_ср.(мкм) и эмпирических параметров D.р.В,п.С.т для пылевых потоков различных энерготехно-логнческих агрегатов приведены в табл. 2.

Таблица 2.

Представленные в табл.2 энерготехнологические агрегаты, из которых были взяты для стендовых испытаний пробы пылей. установлены в технологических линиях цветной металлургии:БКЗ-50-39У (Сред-неуральский медеплавильный комбинат.г. Ревда). РКК-20/40 (Мончегорский комбинат "Североннкель"). УККС-6/40 (Алмалыкский горнсме-таллургический комбинат), печь КС (Алавердский горнометаллургический комбинат), сернокислотных производств: КС-450ВТКУ (Череповецкий комбинат "Аммофос"), огнеупорных производств: КУ-60 (Бсро-Еичевский огнеупорный комбинат), КУ-125 (Челябинский металлургический комбинат, обжиг доломита), КУ-125 (Саткинский огнеупорной комбинат "Магнезит").

Непосредственные эксперименты на котлоагрегатах ь премиален ных условиях при различных режимах подтвердили установлений и стендовых условиях Факт, что при значениях произведений юнпогпрі-

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

Послойные измерения падающих лучистых потоков проводились на котлах-утилизаторах через лючки в обмуровке, в которых размещались спектрорадиометрические приборы с водоохлаждаемыми моделями АЧТ. Геометрическая толщина L изменялась от 0,05 до 1.5 м с шагом 0.05 м путем перемещения модели АЧТ по оптической оси спектрорадиометра вглубь газохода. Температура рабочей среды измерялась методом двух термопар и отсосными пирометрами. Перед экспериментами спектрора-диометры градуировались по плотности излучения АЧТ.

Степень черноты газопылевого потока е_г„ определялась из зависимости е_гп=Е_гп/Е. где Е_ГП.Е - потоки излучения от исследуемой газопылевой среды и АЧТ при температуре среды Т. По полученным в экспериментах значениям степени черноты е_гп находились интегральные коэффициенты ослабления газопылевого потока К_гп из формулы: E_rn=l-exp(-K_rilL). В общем виде, экспериментально полученное значение интегральной степени черноты газопылевого потока также может быть представлено известным выражением: е_гп - е_г + е_п - гп . где е_г, е_п - степени черноты газовой и твердой фаз. Для случая газой фазы, образующейся при обжиге серосодержащего сырья, имеем г^в so2 + Єн2о - sos Сиге где so2. н2о - интегральные степени черноты двуокиси серы и водяного пара.

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

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

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

Даьные по степени черноты двуокиси серы c_S02 показывают, что увеличение L для фиксированных р и Т сопровождается уменьшением влияния L на t_S02- Характер зависимости e_S02 от ^т объясняется положением основных полос в спектре излучения S0₂ и наличием таких факторов как увеличение числа полос и их расширением с ростом Т.

Степень черноты отдельной полосы t_t для S0₂ при увеличении Т растет до тех пор, пока максимум монохроматической плотности излучения АЧТ не перемещается в область длин волн, в которой расположена данная полоса. При дальнейшем повышении Т максимум переходит в область более коротких длин волн и степень черноты.полосы будет убывать. Пока длина волны Х„ ( максимум спектральной плотности излучения АЧТ) превышает значение Xj_Cp.. соответствующее средней длине волны в полосе, расширение полосы и увеличение числа спектральных ЛИНИЙ будет УСКОРЯТЬ РОСТ E_t С РОСТОМ Т. При \_и<^1ер. эти

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

Для сильных полос S0₂ v_t, v₃. 2Vj, V!+v₃ с центрами на дл*'нах волн 8,68, 7,35, 4.34, 4 мкм эти температуры составляют 334, 394. 667 и 724 К. Поэтому для интервала температур, имеющих место в газоходах котлов-утилизаторов, с повышением Т происходит только убывание Ei ВСЄХ ПОЛОС.

Интегральная степень черноты газовой.фазы при обжм;- серосо^ держащего сырья в зависимости от температуры Т описывается Фоому-лой: _г = 1.29 L^0,5(r р L)⁰-¹⁸ ехр(-0,0013 t). где L - геометрическая толщина излучающего слоя.м; г.р - суммарная объемная доля трехатомных газов и общее давление газовой смеси, равное для продуктов обжига серосодержащего сырья 0.1 МПа. Интегральная погло'".а-тельная способность газовой фазы может быть представлена как

0 0,5

А_г = (Т/Т)^0,5 (1- (1- tp)^(T/T ' ). где Г- температура АЧТ; _г - степень черноты газовой фазы при ее температуре Т.

Эксперименты в стендовых условиях и спектральные измерения на натурных котлоагрегатах подтвердили наличие критерия подобия для описания эмиссионных свойств твердой дисперсной Фазы в виде отношения среднего по удельной поверхности F диаметра частиц Х_ср. и длине волны Х_м. соответствующей максимуму излучения АЧТ при температуре частиц Т. При этом интегральный безразмерный коэффициент ослабления лучей взвешенных частиц выражается в виде k_n = О (Х_ср./Х_и)^п, где С,n-эмпирическне коэффициент и показатель.

Для определения оптических констант веществ твердой дисперсной фазы из энерготехнологических агрегатов использован метод Ши-нона, основанный на измерениии коэффициентов зеркального отраженш поляризованного излучения при углах падения луча на поверхності исследуемого образца 20 и 70. Применено состояние поляризации, і котором электрический вектор колеблется параллельно плоскости падения луча. Для реализации выбранного метода измерений принята оптическая схема спектрофотометра "Спекорд М-80". Подготовка образцов с размерами 15x15x2 мм и оптически гладкой поверхностью заклю чалась в прессовании пылевых частиц в прессформе с приложение] сдвигающего усилия в определенном режиме. Качество оптической по верхности контролировалось на профилографе^профилометре П203 пр: вертикальном увеличении 10000 раз для скоростей движения алмазно головки 1 мм'мин и диаграммной бумаги 80 мм/мин. По профилограмма определялись параметры шероховатости, которые сравнивались с мик рошероховатостью эталонного зеркала.

При проведении экспериментов в канал измерения спектрофото метра сначала устанавливалось эталонное зеркало с алюминиевым пок рыгаем, а затем исследуемый образец с оптически гладкой поверх ностью. Коэффициент отражения от исследуемого образца определяло по формуле R-R,(I/I_S). где R,-коэффициент отражения алюминиевог покрытия, равный 97%; I,Ij-интенсивности отраженного излучения с образца и эталона. Значения коэффициентов отражения выводились ї цифровое печатагдее устройство с шагом v»4 см"¹.В дальнейшей обрг ботке использовались результаты расчетов по Формулам Френеля. Зті

результаты протабулированы с интервалом R=0.005%.

Химический состав веществ пылевых частиц при исследовании оптических констант был следующий ( по массе): ГРК-700: rJa₂S0₄---99: 0КГ-100-ЗБ: Fe₂0₃»73,2; Са0=8.5; Fe0-=7,2;S10₂-2,3;Zn0-1.5;Mn0=O, 9; УККС-6/40: Zn0=30,4; CuO=ll. 9; Fe₂0₃=4.2; Fe0=2.32; ⁿb0=2.28; S10₂-l, 9; КУ-125: Fej,0₃=75,3; Mg0=6.3: Ca0=6; Zn0=3,8; FeO-1,8; S10_z-1.3; Р₈0₅-0.9; КС-450ВТКУ: Fe=49, 4; S10₂=29,G;S=6,1; Al₂0₃=3.49; Zn=0.79; Cu=0,45; БКЗ-210-140ф:5102»61.5;А1г0з=19;Ге_г0з=8.7:Са0-3,7;Ре0-2.67;НЕ0-1.9: РКЖ-25/.'О; Fe=30; Cu=2^; S₀6i*.=15: S10₂=10; Pb=5; Zn=5; As-4; C-2.

При X=2-15 mkm вещество частиц из котлоагрегатов относится к промежуточному классу: между металлами, для которых показатель поглощения эе>10 и полупроводниками, для которых 0,01<эе<0,1. Дисперсия оптических констант зависит от соотношения содержания основных химических компонентов, входящих в состав пылей, и их микроструктуры. Показатель преломления п вещества пыли из агрегатг {С-450 ВТКУ имеет наименьшую дисперсию, описываемую полиномом :едьмой степени относительно длины волны.

Основные выводы по второй главе:

1.С ростом температура частиц интегральная степень черноты галевых потоков из этих частиц имеет тенденцию к некоторому повышению. При постоянной температуре частиц, с ростом температуры источника излучения - АЧТ, поглощательная способность пылевых пото-сов повышается для агрегатов: БКЗ-50-39У, РКК-20/40, КС-450 ВТКУ, 'ККС-6/40 и понижается для: печи ICC, КУ-125, КУ-60. Характер влия-іия температуры во многом определяется зависимостью спектрал' ных юэффициентов ослабления от X и дисперсией оптических констант.

2.Влияние температуры частиц на поглощательную способность шлевых потоков более сильное, чем влияние температуры источника ізлучения- АЧТ. Например, для пылевого потока из агрегата 'ККС-6/40 эффект от повышения температуры-частиц от 490 К до 630 К называется на увеличении поглощательноЯ способности в . 2 оаза ;ильнее. чем от повышения температуры источника излучения ОТ 900 К ;о 1300 К.

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

Ввиду селективности спектра излучения потока пылевых частиц их козффищіент ослабления по отношению к собственному излучению оказывается на 15 % выше, чем по отношению к излучению АЧТ.

1. Запыленные газовые потоки при обжиге серосодержащего сырья являются селективными излучателями. Наибольшая селективность наблюдается при малых толщинах слоя L и обусловлена составляющей излучения двуокиси серы. Спектр излучения твердой дисперсной фазы относительно слабо селективен и его степень черноты Ехп Для диапазона 1-15 мкм убывает с ростом X. С ростом L значения е*_п возрастают при дЬ<160- г/м². При ць>160 г/м² зависимость е_Хп от L. при неизменных концентрациях пыли д. ослабевает.

2. Относительный вклад твердой дисперсной фазы в суммарное излучение пылегазового потока при обжиге серосодержащего сырья изменяется в зависимости от X и pL, увеличиваясь с ростом uL. При цЬ>200 г/м² по всему интервалу 1-Ю мкм эта величина составляет 90%. При рШбО г/м² интегральный безразмерный коэффициент ослабления лучей к_п определяется по формуле кп=0.06(Х_ср./Хн)^0,г-

6.Для диапазона длин волн 2-15 мкм и Т*293 К дисперсия оптических констант веществ пылевых частиц из котлов-утилизаторов сильно зависит от соотношения содержания основных химических компонентов и характеризуется значениями: 1<п<1.6 и 0.1<к<1.1.

В третьей главе отражены результаты расчетного исследования основных параметров радиационного теплоообмена по известным математическим моделям. С помощью теории Ни по полученным в настоящей работе оптическим константам вещг -лз твердых дисперсных фаз энерготехнологических агрегатов для случаев монодисперсных частиц рассчитаны безразмерные спектральные коэффициенты ослабления лучей поглощением и рассеянием, которые оказались очень сильно зависимыми от параметра дифракции р в области его значений от 0 до 4 (р =ЯХ/Х. Х-диаметр частицы, мкм; л-длина волны, мкм). Для больших значений р»б0-100 суммарный коэффициент ослабления приближается к двум, а коэффициент рассеяния в среднем на 30-40? выше, чем коэффициент поглощения. С уьеличением показателя преломления п и длины волны X. коэффициент ослабления возрастает при значениях р около 1. При этих р. уменьшение X и рост показателя поглощения ае приводят к повышению коэффициента ослабления лучей поглощением.

Рассчитанные по экспериментальным значениям оптических констант размерные суммарные коэффициенты ослабления лучей к=к_п+к_р.ко- эФФициенты поглощения к_п. рассеяния к_р и критерий Шустера Sc=k_p/k полидисперсных систем выражаются функциями.в зависимости от длины волны Х(мкм).следующего вида: k,k_n.k_p,Sc = аХ^ь (табл. 3).

Таблица 3 . Значения параметров функций для к. к„. к_р, Sc

к, см"' кп, см"¹ кр. см"' Sc

Агрегат ab aba b ab УККС--6/'.0 0.19 0.0^5 0,086 0.119 0,116-0,109 0.566-0,108 КУ-125(м.П. J 0.С28 0.095 0.011 0.2 0,0173 0.0067 0,624-0.088 КС-450ВТКУ 0,096-0.014 0.038 0,106 О,OS -0.13 0.624-0.416 ОКГ-ЮО-ЗБ 0.16 -0.066 0.0592 0.155 0,11 -0.327 0,696-0.261 РКН-25/40 0,035 0.041 0,0135 0,147 0.0217 -0.053 0,625 -0.095 БКЗ-210-140Ф 0,024 -0,22 0.0098 -0.026 0,0163 -0.506 0,687 -0,289 Примечание:(м.п.-мартеновское производство).

Для описания дисперсного состава применены функция логарифмически нормального распределения и модальный диаметр Х_н. Использованные в расчетах значения Х_м(мкм) и концентрации пыли fi (г/нм³) по агрегатам составляли:

УККС-5/40 КУ-125 КС-450ВТКУ ОКГ-ЮО-ЗБ РКН-25/10 БКЗ-210-140Ф
Х„ 0,564 0.125 14.1 2.45 15.5 6,1
Ц 300 50 300 150 150 20.

При значениях оптических констант в области п = 1-1.5, h= 0.2-0,6 и концентрациях частиц 300 г/нм³ для диапазона спектра 1-15 мкм уменьшение модального диаметра с 14,1 до 0,564 мкм (^.е. ^я в 25 раз) приводит к росту суммарного коэффициента ослабления в 2-2,5 раза. При этом коэффициенты поглощения также возрастают в среднем в 2-3 раза. Различный характер поведения к_п и к_р с увеличением длины волны X не дает более сильного возрастания суммарного коэффициента ослабления лучей ( к_п растет с увеличением X, к_р - падает). Снижение концентраций частиц с 300 до 150 г/нм³ уменьшает суммарный коэффициент ослабления к в 2-2,5 раза.

Увеличение коэффициента ослабления к происходит наиболее интенсивно при снижении показателя преломления п. чем при уменьшении, модального диаметр частиц Х_н при прочих равных услозиях. Обчей ' теденцией является уменьшение критерия Шустера Sc ростом X.

При известных F.u. L осуществляется переход от рассчитанных размерных спектральных коэффициентов ослабления лучеЛ к безразмерным, осредненные значения которых во всем диапазоне длин волн X:

со со со

K_s=;k_xE⁰xdX//E⁰xdX=(aCi/6T⁴)/(l/X^{b + 5})(exp(c_?/XT)-l)-¹c!X=

⁰ со

=((aCjT ^b)/(6⁰c_z^b*⁴))/z^b*³:e^z-l)"¹riz. (12)

где kx=aX^b -функция дисперсии коэффициента ослабления; с_1гс_г-1-я и 2-я константы излучения Планка;6-постоянная СтеФана-Больцмана. Последний интеграл сводится к гамма-функции ст параметра b

; z^b+3(e²-l)~¹dz=1.04r(b+4). о

01'эзначив (асі1.04Г(Ь+4))/(бс₂^ь+4) =п, получаем зависимость, рассчитанного по оптическим константам, среднепланковсксго коэффициента ослабления лучей от температуры Т: k_s= nT ^ь.

Полученные в настоящей работе расчетные данные по критерию Шустера и экспериментальные - по эмиссионным свойствам двухфазных сред в газоходах энерготехнологических агрегатов позволили, по известным математическим моделям, произвести анализ зависимости коэффициента -тепловой эффективности экранов (КТЭ) \|) от полей температур, радиационных характеристик рабочих сред и граничных поверхностей. ВеЛИЧИНа ір = а_Пал-аоОр/Чпад^сЯрез/Чпа_А. ГДЄ Япад. Qoup. q_Pe3~

падающий. обратный и результируюгчй потоки.

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

ется применение понятия эффективной приведенной температуры 9эсгИ_ЭС1. Изменением в_э$ охватываются все встречающиеся на практике сочетания диапазонов изменения X и Т. Конкретное распределение температур з газоходах котлов-утилизаторов по результатам экспериментальных обследований выражается в виде аналитических зависимостей для симметричных профилей с изотермическим ядром.

Для случая двухфазной рассеивающей среды применима формула,
полученная Адзерихо К.С.: д^^ХТ^=с_г{1п{1+1/к))'¹. (13)

где A=(K/(l-exp(-Kt₀>>/ (ехр(-К(х₀-т')))/(ехр(с_г/ХТ(т'))-1).

K=2(i-Sc)^0,s; т₀.т'- общее и текущее значения оптических толщин.

Расчеты показали, что-при %₀>5 с увеличением критерия Шустера Sc от 0 до 0,6 0_э0 возрастает на 11 %, а при т₀<1 - на 1-1,5 %. С ростом полки изотермического ядра т_ц в 2 раза значение 8_эо уменьшается в среднем на 10 % при т₀<5 , температуре центра 8_Ц=3-10~³ м-К и 8ц/8_Ст=2. где О_ст-прпведекная температура стенки. Сравнение корректно вводимой эффективной температуры T_D0 с широко используе-чыми для описания неизотермических объемов среднеарифметической Г_са=0.5(Тц+Т_ст) и среднегеометрической Т_сг=(Т_аТ_Ст)⁰'⁵ температур доказывает, что при оптических толщинах т₀<2 и 0_Ц=3-Ю"³ м-К величины Т_са и Т_сг занижены в среднем на 25%. Для т_о=15-20. Sc=0, т_ц=0 і 8ц/8ст>2 значения Т_са и Т_сг являются завышенными по отношению' к :_э$ на 15-35%. При увеличении критерия Шустера Sc от 0 до 0.6 для ;_о=0,5-2, т_ц=0,5т₀, 8_ц=3-10"³ м-К, 8_Ц/9_СТ=2 эффективная степень іернотьі е_эо уменьшается в среднем на 20 55. Влияние степени черноты >граничиващих стенок е_ст заметно сказывается на величине эффек-чівной степени' черноты _Эо лишь при относительно небольших опти-іескіїх толщинах т₀<1.5. При г₀=0,5, т_ц=0, Sc-О возрастание с_от от 1.40 до 0,75 понгскает с_э0 на 15 й.

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

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

і.При прочих одинаковых характеристиках полидисперсного пылевого потока изменение показателя преломления п оказыгает большее влияние на коэффициенты ослабления лучей, чем изменение модального диаметра. С ростом показателя преломления в 1,3-1,4 раза: коэффициент ослабления лучей за счет поглощения уменьшается в 7-8 раз, а коэффициент ослабления рассеянием уменьшается в 6-7 раз. С ростом .модального диаметра в 6-7 раз коэффициент поглощения уменьшается в 4-5 раз.

2. В области длин волн Х=2-1б мкм дисперсия коэффициента ослабления лучей к* пылевыми потоками энерготехнологических агрегатов ' выражается в виде:кх=аХ^ь, где а, Ь- параметры, определяемые первичными радиационными характеристиками. Температурная зависимость среднепланковского коэффициента ослабления лучей, рассчитанного по оптическим константам, описывается формулой к₃=пТ^ь.

3. Значения употребляемых при расчетах радиационного теплообмена в котлах-утилизаторах среднеарифметической Т_са и среднегеометрической Т_сг температур, при малых оптических толщинах т <5 и приведенных-температурах центрального ядра 6_ц<3-10"³ м-К, занижены, а при т >15. завышены относительно корректно вводимой эффективной температуры Т_Эф неизотермического объема рабочей среды. С увеличением критерия Шустера Sc г~ 0.3 до 0,5 для т_о=0,5-2 эффективная температура Т_Э1> возрастает не более,чем на 1.5%. а эффективная степень черноты Сэо убывает не более, чем на 4%.

4.Для условий, имеющих место в газоходах котлов-утилизаторов, при т₀<4, Тц/Тст=2. 6ц(1.5-5)-10^-3 м-К. коэффициент тепловой эффективности экранов i|) очень мало зависит от т₀ и степени черноты стенки с_ст. Величина гр для исследованных агрегатов составляет 0.4-0.75, сильно уменьшается с ростом длины волны X и возрастает с увеличением Т_э&.

Четвертая глава посвящена применению полученных данных по из-лучательноЯ способности материалов и пылегазсвых сред в зональном методе расчета теплообмена котлов. Как показали экспериментьдгл.2)

для сильной запыленности и высокой концентрации S0₂ коэффициент ослабления лучей понижается с ростом L. Поэтому при укрупнении объемных зон дол::;но соблюдаться условие: макснгэльное значение эффективной длины пути луча L₃₄2: для обжигового газа (при парциальном давлении S0₂ р=и.008-0,014 МПа) pL=15 МПа-мм.

Определение оптимального числа единичных исгатаний Р проводилось на геометрической модели з форме параллелепипеда, заполненного диатермической средой. Сравнение точных аналитических значении угловых коэффициентов с вычисленными методом Монте-Карло в диапазоне р=500-3000. показывает, что средневзвешенная погрешность. о₀р связана с числом испытаний Р зависимостью: б_ср = 73,6/Р^0,426.

Исследование влияния селективности спектра излучения сернистого газа на точность расчета радиационного теплообмена зональным методом проводилось с использованием модели Хоттеля, по которой полосатый спектр трехатомного газа представляется в виде трех серых газов и одного прозрачного газа. Первоначально определялись весовые коэффициенты а„, b_n для каждой полосы спектра п в зависимости от температур газа Т и АЧТ Т. затем - коэффициенты поглощения к„ и коэффициенты полиномов: а_п0 - для излучательной способности, Ьпо - для поглощательной способности сернистого газа при парциальном давлении двуокиси серн 0,012 Мпа, имеющие значения: .

п Излучательная способность __1т

, к„,1/(м-МПа) а„ о а_п ц10³ а_{я 2}^10⁶ а_{п З}іі0⁹ ^

1. 2,06 0,4344 0,2271 -0.2823 0,0554

2. 26.78 0,1471 0,2819 -0,3101 0,6670

3. 810 0,1911 -0,1309 0,0253 -0,0066 .

п. Поглощательная способность при температуре газа 1^1200 К ._

, b_n о b„, JL0 ³ Ь„ _г -JO ^G b_n зЦО⁹

4. 1,1621 -1,9712 1.6071 --,5138

5. 0,3621 -0,4056 0,2733 -0.06S4

3 0,2086 -0, 3091 0,2381 -0, 0694 .

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

При обжиге серного колчедана для двухфазной среды расхождения в результатах расчетов по селективной и серой моделям имеют меньшую величину, наиболее сильно проявляясь при небольших значениях эффективных толщин излучающенго слоя до 0,1 м. Для котла-утилизатора при содержании огарковой пыли в обжиговом газе 200 г/нм³ коэффициент передачи тепла конвекцией составляет по х^ау движения рабочей среды: на входе в ширмезый пучок 11%, ь середине высоты ширм - 15%, на выходе из ширм - 18 от значения коэффициента теплопередачи излучением.

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

Основные выводы по четвертой главе:

1.В расчете теплообмена зональным методом при геометрическом моделировании укрупнение объемных зон и зон кладки для постоянного деления на поверхностные зоны (экраны) не приводит к существенным погрешностям. Расхождение результатов по суммарному тепловосприя-тию для 50.- и 30-зоиной геометрических моделей составляет около 4%, по зональному тепловосприятию максимальное расхождение 9%. Среднее по всей поверхности нагрева суммарное тепловосприятие для 60-зонной модели на 5% превышает тепловосприятие, рассчитанное для 45-зокной модели.

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

3.Сравнение результатов расчета угловых коэффициентов методом статистических испытаний с точньсш аналитическими значениями показывает, что погрешность убывает с увеличением длины выборки случайных чисел (затрачиваемого машинного времени). Оптимальная величина длины выборки 1200-1300 испытаний при погрешности 3%.

4.Суммарный тепловой поток, рассчитанный по селективной модели, на 5-7 % меньше, чем поток, определенный по серой модели. Радиационные потоки по отдельным участкам ширм п;лі расчете по селективной модели также ниже на 5-7 % по сравнению с серой моделью. Пренебрежение селективными свойствами газовой фг^ы приводит к завышению теплоотдачи от пылегазового потока и к занижению теплоотдачи от поверхностей, имеющих сплоыной спектр излучения.

5.Средние температуры объемных и поверхностных зон. полученные в р:зультате расче-а по селективной и серой моделям, различаются между собой в среднем на 20-50. При этом температуры объемных зон по селективной модели выше на 20-50, чем температуры объемных зон. рассчитанные по серой модели. Для поверхностных зон температуры, рассчитанные по селективной модели, інше на 20-30". чем температуры, рассчитанные по серой модели. Машинное время, затрачиваемое для расчета по селективной модели среднем на 20" превышает время, затрачиваемое на расчет теплообмена зональным методом по сеоой модели.

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

7.Для исследованных котлоагрегатов по результатам расчета зо-[альным методом доля тепла, передаваемого конвекцией, составляет в іреднем 20% от общего тепловосприятия (при концентрации пыли в поте д более 100-150 г/нм³). Величина погрешности при раочете обобщенных угловых коэффициентов методом статистических испытаний юзрастает с увеличением количества зон. или иначе, с уменьшением относительных размеров зон. Для кадцого котла-утилизатора сущетву-т' оптимальное количество зон. при котором достигается наименьшая :огрешность расчета для постоянного числа испытаний.

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

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

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

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

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

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

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

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

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

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