Адаптация дифференциального уравнения энергии к условиям топочных процессов в котельных агрегатах Текст научной статьи по специальности «Физика»

Теплоэнергетика

УДК 621.1.016 й01: 10.14529/power150101

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

Е.В. Торопов, К.В. Осинцев

Южно-Уральский государственный университет, г. Челябинск

Предложена форма дифференциального уравнения энергии для потока среды в топке котельного агрегата с учетом значимых для данного случая механизмов переноса теплоты: конвекции, излучения, теплопроводности и объемных тепловыделений. Коэффициенты переноса теплоты в уравнении энергии адаптированы к реальным условиям процесса; раскрыт их физический смысл и пределы варьирования, определены числа подобия и интервалы их влияния на пространственное температурное поле топки, режимы течения и теплообмен топочной среды. Показано, что при наблюдаемых на практике значениях оптической плотности среды в топках 2,8-3,6 в математическом описании вектора плотности потока излучения возможен переход от сферической зависимости к градиентной форме. Также градиентная форма описания рекомендована для конвективно-кондуктивного механизма образования температурных молей в топке. Определены участки перехода к турбулентному режиму в пристенном пограничном слое, масштаб турбулентных полей при этом достигает 0,3 м. Уточнено влияние температуры потока среды на ее теплофизические свойства.

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

Введение

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

Уравнение переноса энергии в топке КА

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

X (Ч, ) = о, (1)

г

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

Рассмотрим плоское двумерное течение = 0), где координата х направлена вдоль потока среды, а координата у - поперек течения, причем у = 0 на стенке КА; тогда уравнение энергии (1) примет вид

РСр

= 1

дТ

wx--Ь w

дх

( д 2т

д 2Т

дТ_

у ду \ (

дх2 ду2

дЧ

и, х + дЧг

и, у

дх

ду

+ Чу

(2)

Введем безразмерные переменные

Х = х/10, У = у/10, Яв = w0l0/v, №х = wJw0, Жу = Wy/wо, 0 = Т/Т0, <Ч)и,х = Чи,х/00Тэ4, <ч)и,у = Чи,у /О0Т04;

здесь 10 - характерный размер рассматриваемой системы, м, который рационально принять равным продольному размеру топки КА, тогда все остальные размеры, как и текущие координаты точки, будут составлять доли от 10. В качестве характерной скорости w0, м/с, можно принять скорость в основном течении среды, вдали от пограничных слоев, остальные переменные с нижним индексом «0» характеризуются подобным же образом. При этом подразумевается, что уравнения сплошности, сохранения импульса и теплоотдачи на границах системы записываются как обычно. В этом уравнении также введена объемная плотность тепловыделения чУ, Вт/м3, которая принята постоянной, не зависящей от х и у, но это ограничение не принципиально, можно ввести в уравнение заданные пространственные или температурные функции объемной плотности тепловыделения, это не усложнит задачу анализа.

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

59 w 59 1

Wx-+ W — = —

x 5X y 5Y Pe

( Я2,

5 29 5 +

2

K5X2 5Y

_1_

Bo

5{q)

5X

u,x +5< q)

U, y

5Y

Po

Pe'

(3)

Адаптация коэффициентов переноса

к реальным условиям в топке

Для условий теплообмена в топке КА числа подобия, входящие в уравнение (3), имеют следующее теплофизическое содержание. Число Пекле Ре = ЯеРг характеризует относительный вклад в формирование температурного поля потока среды в топке КА конвективного механизма переноса теплоты по сравнению с теплопроводностью. Число Померанцева Ро = д„/02/ХГ0 характеризует безразмерную плотность тепловыделений.

Число Больцмана Во = (рм>)0ср /сГ03 описывает соотношение потоков энергии конвективного и радиационного характера, а величина 1/Во - участие радиационного механизма в формировании температурного поля потока среды в топке. Роль числа Во в радиационно-конвективном теплообмене аналогична роли числа Ре в кондуктивно-конвективном теплообмене: когда число Ре мало, можно пренебречь переносом тепла конвекцией по сравнению с переносом тепла теплопроводностью. Когда число Во имеет малую величину, можно пренебречь вкладом конвекции в формирование температурного поля потока среды по сравнению с излучением. Число Во легко преобразуется в комплекс чисел подобия Во = (Ке/К)(Рг/Вы), где, кроме отмеченных Яе и Рг, приведены числа N = с0ТЛк и Бугера Ви = к10, где с0 = 5,67-10-8 Вт м2/К4 - коэффициент излучения абсолютно черного тела, к - коэффициент ослабления топочной среды, м-1. Число N является радиационно-кондуктивным параметром, число Ви характеризует излучательную способность топочной среды КА. Приведенная зависимость для Во является определением, получить из нее расчетные зависимости можно только после установления количественных зависимостей между величинами N, Яe, Pг, Bu, Ро - аналитически, численно или экспериментально.

Выбор температур Т и Т0 имеет принципиальное значение для анализа влияния числа Больцма-на на процессы в топке КА и он связан с задачами анализа. Так, при определении результата теплообмена между координатами х1 и х2 можно записать уравнение для результирующего числа Больц-мана АВо = Ворез.

Ворез =pwcp (71 -72 )/ст0 (4 -?24 ).

(4)

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

число Во, основанное на средней по участку температуре, при этом температура для конвективного потока определяется как средняя арифметическая, Тср = (Т1 + Т2)/2, а для потоков излучения необходимо другое усреднение Тср4 = (Т14 + Т24)/2. Расхождение между средними температурами, подсчитываемыми по этим формулам, зависит от разности Т1 - Т2 = ДТ: при ДТ = 100 К разница 0,13 %, при ДТ = 300 К разница достигает 1,2 %, что также можно признать допустимым.

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

Число Во в виде сомножителя 1/Во со знаком «минус» входит параметром во второе слагаемое правой части уравнения (3), это слагаемое описывает влияние потоков излучения на формирование температурного поля топки КА. Входящее в Во число Ви = к10 является параметром для суммы градиентов плотности потоков излучения {д)и х и (д)и у по координатным направлениям:

Ви ~^га^ {д)иху). С увеличением Ви снижается

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

Сравнив характерный размер топочного пространства 10, м, с величиной, обратной коэффициенту ослабления 1/к, получим Ви = /0/(1/к) = к/0; если Ви « 1,0, среда считается оптически тонкой, при Ви» 1,0 среду считают оптически толстой (массивной). В топках КА при Ви = 2,8-3,6 оптическая массивность среды позволяет перейти в математическом описании вектора плотности потока излучения от интеграла интенсивности излучения по сферическому углу к градиентной форме

Торопов Е.В., Осинцев К.В.

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

описания Чи, что и применено при выводе соотношений (2), (3). Приведенные выше соображения подтверждают достоверность примененной формы описания дифференциального уравнения энергии.

Преобразованное нелинейное уравнение (3) не поддается аналитическому решению, как и исходное (2), но позволяет выделить несколько аспектов анализа, связанных с теплофизикой и теплообменом. Левая часть уравнения (3) описывает изменение температуры потока безвихревого течения со скоростью, имеющей компоненты №х, Жу в пространственном поле температур с градиентами д0/дХ, д0/дУ, и это изменение связано с действием конвективного, кондуктивного и радиационного механизмов, а также с тепловыделениями в потоке среды. Ограничение относительно безвихревого течения можно ослабить, если вместо молекулярного коэффициента теплопроводности применить эффективный коэффициент, учитывающий турбулентные пульсации, Хэф = + Хт.

Первое слагаемое правой части (3) в общем соответствует решению уравнения Лапласа, относящегося к теории потенциалов [2], которое в частном случае плоской поверхности дает описание плоских изотермических поверхностей, определяющая температура которых растет при удалении от поверхности экранов КА. Если раскрыть множитель перед скобкой 1/Рв = Xэф/pwcp/0, то становится ясно, что первое слагаемое (3) описывает соотношение лапласовских температурных полей, образованных механизмами теплопроводности и конвекции.

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

ка, отсчитывается от начала течения х = 0, то есть от оси горелочных устройств до оси горизонтального газохода х = 10 , при этом безразмерная координата точки изменяется в пределах 0 < Ьт < 1,0.

Количественные оценки коэффициентов

переноса

Режим движения потока среды вдоль вертикальных стен топки КА определяется двумя критическими значениями Явх,1 и Явх,2, первое соответствует турбулизации ламинарного пограничного слоя, второе - полностью развитому турбулентному пограничному слою. В научной литературе [3] обычно приводятся зависимости числа Рейнольдса от степени турбулентного возмущения потока для идеализированных условий в точке встречи потока с поверхностью х = 0. При этом учитывается влияние шероховатости стенки, условия формирования слоя при х = 0, угол встречи потока с поверхностью и другие факторы. В теплоэнергетических приложениях можно считать, что полностью развитая турбулентность наступает при Явх,кр = 5 105, откуда можно получить зависимость для хкр = 5 105^ А^), результаты расчетов по которой приведены в таблице.

Определение хкр имеет значение для адаптации процесса развития турбулентности и в этом отношении для расчета толщины турбулентного пограничного слоя 5 можно использовать рекомендуемую в [3] зависимость 5 = 0,37х-Явх-0'2, которую для повышения точности графического изображения можно представить как 5усл = 5/х0,8 = = 0,37(^ ^/,2. На основе приведенных данных можно сделать вывод о том, что при х > хкр толщина турбулентного пограничного слоя растет и к точке х = 10 м составляет 0,189 м, а к точке х = 20 м достигает 0,329 м.

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

Зависимость координаты полной турбулизации пограничного слоя хкр, м, (числитель) и условной толщины турбулентного пограничного слоя, 8усл, м0'2, (знаменатель) от плотности потока массы топочной среды, рп>, кг/м2-с, и ее температуры Т, К

pw, кг/м2-с Температура Т, К

1473 1573 1673 1773

15 1,75/0,0300 1,82/0,0302 1,89/0,0305 1,96/0,0307

20 1,32/0,0284 1,37/0,0286 1,42/0,0288 1,47/0,0298

25 1,05/0,0272 1,09/0,0274 1,14/0,0276 1,17/0,0278

30 0,87/0,0262 0,91/0,0264 0,95/0,0266 0,98/0,0268

тем передают ее вихрям меньшего масштаба, и далее передача энергии продолжается с уменьшением масштаба. В самых мелких вихрях энергия диссипирует, то есть превращается в теплоту, вклад этого источника энергии незначителен по сравнению с другими источниками в уравнениях (2) и (3). Считается, что таким образом происходит диссипация кинетической энергии потока газов в топке КА вследствие работы сил вязкости; этот каскадный процесс передачи энергии является характерной чертой турбулентного движения. Другой его особенностью является трехмерность пульсационных движений, когда мгновенные значения скорости в произвольной точке турбулентного потока зависят от трех координат.

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

Режимные оценки, приведенные выше и в таблице, имеют предельный (ш1п1-шах) характер. Фактически, в условиях высокой турбулентности набегающего потока при х < хкр образование ламинарного пограничного слоя маловероятно, поэтому оценки турбулентности в пограничном слое имеют значение минимальных. С другой стороны, плотность потока массы рм> = 10...20 кг/м2-с для показателей среднего течения среды в топке КА недостижимы и могут наблюдаться только в применении к локальным показателям как максимальные.

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

при выводе основных характеристик движения. Действительно, уравнение сплошности div(pw) = 0 позволяет не вводить температурные поправки в характеристики течения, так как pt = p0T0/T, wt = w0T/T0, то (pw)t = (pw)0. Влияние температуры на характеристики потока среды проявляется через температурную зависимость абсолютной вязкости, которая в исследованном диапазоне Т = 1273... 1773 К может быть представлена в виде fit = [47,9 + 0,022(T- 1273)]10-6, кг/мс = Нс/м2 = = Па с, для продуктов сгорания среднего состава (р = 0,101 МПа, СО2 = 13 %, Н2О = 11 %, N2 = 76 %).

Объемная плотность тепловыделений в потоке топочной среды qv, Вт/м3, в уравнении (2) принята постоянной, это допустимо в качестве первого приближения при анализе и адаптации процессов в топке КА. После преобразований уравнения (2) в уравнение (3) этот член преобразуется в отношение Ро/Ре; интегральная оценка этого источника тепло -ты может быть получена делением тепловой нагрузки КА Вренр, Вт, на общий объем топки VT, м3: qv = BpQn/VT, где Вр - расчетный расход топлива, кг/с. Но эта оценка не очень точно отражает функцию тепловыделения, так как сгорание топлива в объеме топки происходит неравномерно, для уточнения этого показателя рекомендуется применить те же соображения, что и для числа Во.

Заключение

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

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

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

Литература

1. Равич, М.Б. Упрощенная методика теплотехнических расчетов / М.Б. Равич. - М.: Изд-во АН СССР, 1966. - 407 с.

2. Корн, Г. Справочник по математике для научных работников и инженеров: определения, теоремы, формулы: пер. с англ. /Г. Корн, Т. Корн. -4-е изд., пер. со 2-го америк. перераб. изд. - М.: Наука, 1977. - 832 с.

3. Кутателадзе, С. С. Теплопередача и гидродинамическое сопротивление: справ. пособие / С.С. Кутателадзе. - М.: Энергоатомиздат, 1990. -367 с.

4. Осинцев, К.В. Расчет теплообмена в топке котельного агрегата: учеб. пособие / К.В. Осинцев, Е.В.Торопов. - Челябинск: Издат. центр ЮУрГУ, 2010. - 179 с.

Торопов Е.В., Осинцев К.В.

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

5. Особенности процессов горения и теплообмена в камерных топках котлов с фронтальной компоновкой горелок / В.В. Осинцев, Е.В. Торопов, М.П. Сухарев, К.В. Осинцев; под ред. проф. Е.В. То-ропова. - Челябинск: ООО «Изд-во Рекпол», 2006. -204 с.

6. Организация факельного сжигания низкосортного твердого топлива и природного газа в топках котлов с фронтальной компоновкой горелок / В.В. Осинцев, Е.В. Торопов, М.П. Сухарев, К.В. Осинцев; под ред. Е.В. Торопова. - Челябинск: Изд-во ЮУрГУ, 2007. - 150 с.

Торопов Евгений Васильевич, д-р техн. наук, заслуженный деятель науки и техники РФ, профессор, кафедра «Промышленная теплоэнергетика», Южно-Уральский государственный университет, г. Челябинск; evtor@mail.ru.

Осинцев Константин Владимирович, канд. техн. наук, доцент, кафедра «Промышленная теплоэнергетика», Южно-Уральский государственный университет, г. Челябинск; osintsev2008@yandex.ru.

Поступила в редакцию 8 декабря 2014 г.

DOI: 10.14529/power150101

ADAPTATION OF DIFFERENTIAL EQUATION OF THE ENERGY TO CONDITIONS OF FURNACE PROCESSES IN THE BOILER UNITS

E.V. Toropov, South Ural State University, Chelyabinsk, Russian Federation, evtor@mail.ru,

K.V. Osintsev, South Ural State University, Chelyabinsk, Russian Federation, osintsev2008@yandex.ru

The paper presents the form of a differential equation of energy for the flow in the furnace of the boiler unit. This form takes into account the mechanisms of heat transfer significant for this case. Examples of such mechanisms are convection, radiation, heat conductivity and volumetric generation of heat. Heat transfer coefficients adapted to the real conditions of the process. Their physical meanings and limits of variation are revealed in the paper. In addition, the authors define the similarity numbers and intervals of their impact on spatial temperature field of the boiler furnace, regimes of the flow, heat transfer in the contents of furnace volume. It is shown that if the observed practical values of the optical density of the medium in the boiler furnaces are 2.8-3.6 then for the mathematical description of the radiation flux density vector the transition is possible from spherical dependence to the gradient dependence. Also, the gradient form is recommended for convection-conduction mechanism of the formation of temperature moles in the boiler furnaces. The areas of transition to the turbulent regime in the boundary layer, the scale of the turbulent fields in this case rises up to 0.3 m. The paper specifies the effect of temperature of the flow on its thermal and physical properties.

Keywords: energy equation, heat transfer, aeromechanics, furnace boiler, boiler unit.

References

1. Ravich M.B. Uproshchennaya metodika teplotekhnicheskikh raschetov [Simplified Method Thermal Calculations]. Moscow, 1962. 303 p.

2. Korn G., Korn T. Spravochnik po matematike dlya nauchnykh rabotnikov i ruinzhenerov [Mathematical Handbook for Scientists and Engineers]. Moscow,1968. 720 p.

3. Kutateladzhe S.S. Teploperedacha i gidrodinamicheskoe soprotivlenie [Heat Transfer and Flow Resistance]. Moscow, Energoatomizhdat Publ., 1990. 367 p.

4. Osintsev K.V., Toropov E.V. Raschet teploobmena v topke kotelnogo agregata [Calculation of Heat Transfer in the Furnace of the Boiler Unit]. Chelyabinsk, South Ural St. Univ. Publ., 2010. 179 p.

5. Osintsev V.V., Toropov E.V., Sykharev M.P., Osintsev K.V. Osobennosti protsessov goreniya i teploobmena v kamernikh topkakh kotlov s frontalnoy komponovkoy gorelok [Features of Combustion and Heat Transfer in the Chamber Boiler Furnaces from the Front Burner Layout]. Chelyabinsk, Rekpol Publ., 2006. 204 p.

6. Osintsev V.V., Toropov E.V., Sykharev M.P., Osintsev K.V. Organizatsiya fakelnogo szhiganiya nizko-sortnogo tverdogo topliva i prirodnogo gaza v topkakh kotlov s frontal noi komponovkoy gorelok [Organization Flaring of Low-Grade Solid Fuels and Natural Gas in Boilers with Front Burner Layout]. Chelyabinsk, South Ural State University Publ., 2007. 204 p.

Received 8 December 2014

БИБЛИОГРАФИЧЕСКОЕ ОПИСАНИЕ СТАТЬИ

Торопов, Е.В. Адаптация дифференциального уравнения энергии к условиям топочных процессов в котельных агрегатах / Е.В. Торопов, К.В. Осинцев // Вестник ЮУрГУ. Серия «Энергетика». - 2015. -Т. 15, № 1. - С. 5-10. Б01: 10.14529/power150101

REFERENCE TO ARTICLE

Toropov E.V., Osintsev K.V. Adaptation of Differential Equation of the Energy to Conditions of Furnace Processes in the Boiler Units. Bulletin of the South Ural State University. Ser. Power Engineering, 2015, vol. 15, no. 1, pp. 5-10. (in Russ.) DOI: 10.14529/power150101