Научная статья на тему 'Термообработка угольного слоя СВЧ-энергией: аналитическое исследование в условиях теплосброса II и III рода'

Термообработка угольного слоя СВЧ-энергией: аналитическое исследование в условиях теплосброса II и III рода Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

CC BY
243
97
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
УРАВНЕНИЕ ТЕПЛОПРОВОДНОСТИ ФУРЬЕ / УГОЛЬ / МИКРОВОЛНОВАЯ ЭНЕРГИЯ / ЗАКОН БУГЕРА / СМЕШАННЫЕ ГРАНИЧНЫЕ УСЛОВИЯ / ТЕПЛОВОЙ РЕЖИМ / ТЕПЛООТВОД / THE FOURIER HEAT EQUATION / COAL / MICROWAVE ENERGY / THE BOUGUER LAW / MIXED BOUNDARY CONDITIONS / THERMAL REGIME / HEAT REMOVAL

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Саломатов Владимир Васильевич, Карелин Вадим Александрович, Заворин Александр Сергеевич, Саломатов Василий Владимирович, Пащенко Сергей Эдуардович

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по электротехнике, электронной технике, информационным технологиям , автор научной работы — Саломатов Владимир Васильевич, Карелин Вадим Александрович, Заворин Александр Сергеевич, Саломатов Василий Владимирович, Пащенко Сергей Эдуардович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Relevance of the research is related to the necessity of scientific study of microwave exposure technology parameters for coal, especially when drying, at thermal treatment, burning intensification and others. Experiments show that use of microwave radiation help reduce harmful emissions from coal burning, increase energy efficiency, reduce process time, etc. For optimal conditions of microwave treatment it is very popular to search for theoretical approaches, in particular, the engineering practice requires analytical solutions for heat treatment of coal array, which was done in the study. The aim of the research is to develop strict analytic solutions of the problems of coal layer heating by microwave radiation under the heat removal of the II and III order, which allow carrying out the parametric analysis of microwave exposure and searching for the most advantageous modes of heating coal layer. Methods. The tasks of coal bed heating by microwave radiation with boundary conditions of the II and III order were set with a number of simplifications, such as persistence of thermal and electrical properties of the material, they are isotropy, one-dimensionality. This physical process itself is plane electromagnetic wave energy absorption and formation of heat sources in coal array, modeled on the Bouguer law. The energy equation in the form of Fourier in this case is considered independently of the Maxwell equations of electrodynamics. A new dependent variable was determined. It modifies the boundary conditions into homogeneous, and simplifies the search of solutions. As a result, the basic system of equations is divided into two subsystems, and the final solution was found making a superposition of these two sub-tasks. Results. The authors have found analytically strong dependences of temperature distribution over the cross section and time for carbon bed heated by microwave radiation, under heat removal of II and III order of boundary surfaces. These solutions are the basis for optimum parameters of microwave heating technology, a tool to reduce the spread of temperature field on the desired temperature distribution within the layer. They help determine quickly: maximum temperature inside a coal array, coordinates of its location, thermo-destructive voltage, limits of discharged heat fluxes and other.

Текст научной работы на тему «Термообработка угольного слоя СВЧ-энергией: аналитическое исследование в условиях теплосброса II и III рода»

УДК 536.24

ТЕРМООБРАБОТКА УГОЛЬНОГО СЛОЯ СВЧ-ЭНЕРГИЕЙ: АНАЛИТИЧЕСКОЕ ИССЛЕДОВАНИЕ

В УСЛОВИЯХ ТЕПЛОСБРОСА II И III РОДА

Саломатов Владимир Васильевич1,

vvs@itp.nsc.ru

Карелин Вадим Александрович2,

vad2hen@mail.ru

Заворин Александр Сергеевич3,

zavorin@tpu.ru

Саломатов Василий Владимирович1,

vvsalomatov@mail.ru

Пащенко Сергей Эдуардович2,

pashenko-ecolog@yandex.ru

1 Институт теплофизики СО РАН,

Россия, 630090, г. Новосибирск, пр. акад. Лаврентьева, 1.

2 Национальный исследовательский Новосибирский государственный университет, Россия, 630090, г. Новосибирск, ул. Пирогова, 2.

3 Национальный исследовательский Томский политехнический университет, Россия, 634050, г. Томск, пр. Ленина, 30.

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

Цель работы: построение строгих аналитических решений задач нагрева угольного слоя микроволновым излучением в условиях теплосброса II и III рода, позволяющих провести параметрический анализ СВЧ-воздействия, а также найти наивыгоднейшие режимы нагрева угольного слоя.

Методы исследования. Задачи нагрева угольного слоя микроволновым излучением с граничными условиями II и III рода были поставлены с рядом упрощений, таких как постоянство теплофизических и электрофизических свойств материала, их изотропность, одномерность. При воздействии СВЧ происходит объемное поглощение тепла с появлением объемных источников тепла, которые в данной задаче моделировались по закону Бугера. В ходе решения была найдена новая зависимая переменная, трансформирующая граничные условия в однородные и упрощающая поиск самого решения. В результате основная система уравнений была разбита на две подсистемы, а итоговое решение было найдено методом суперпозиции решений этих двух подзадач.

Результаты. Найдены аналитически строгие зависимости распределения температуры по сечению и времени для угольного слоя, нагреваемого микроволновым излучением, в условиях несимметричного теплоотвода II и III рода по граничным поверхностям. Приведенные решения являются основой для получения оптимальных параметров технологии СВЧ-нагрева, инструментом для уменьшения разброса температурного поля от требуемого распределения температуры внутри слоя. С их помощью оперативно определяются: значение и координата максимальной температуры внутри угольного массива, терморазрушающие напряжения, предельные значения сбрасываемых тепловых потоков и др.

Ключевые слова:

Уравнение теплопроводности Фурье, уголь, микроволновая энергия, закон Бугера, смешанные граничные условия, тепловой режим, теплоотвод.

Введение

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

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

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

В настоящее время имеется целый ряд фундаментальных экспериментальных и теоретических работ о нагреве микроволновым излучением различных материалов, включая уголь. Среди работ по изучению воздействия СВЧ-излучения к наиболее значимым можно отнести теоретические исследования [1-7]. Также с использованием различных установок были реализованы экспериментальные работы о СВЧ-нагреве диэлектрических материалов [8-17]. Математическое моделирование электромагнитных и тепловых полей при СВЧ-нагреве диэлектрических сред выполнили В.И. Анфино-гентов, Г.А. Морозов, Т.К. Гараев и др. специалисты Казанского государственного технического университета им. А.Н. Туполева [18, 19]. В Иркутском государственном техническом университете Н.П. Коноваловым, П.Н. Коноваловым, В.В. Фед-чишиным, А.А. Хайдуровой проводились экспериментальные исследования по повышению технических характеристик бурого угля воздействием СВЧ-излучения [20, 21]. Полученные в перечисленных работах результаты вносят весомый вклад в данную проблематику, но полностью не закрывают потребности микроволной обработки материалов. Требуется дальнейшее развитие этого направления исследований, включая, в частности, поиск аналитически строгих зависимостей нагрева материала энергией СВЧ-излучения.

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

Данная статья является продолжением работы [22] тех же авторов, в которой приведены решения для задачи микроволнового нагрева со смешанными условиями теплоотвода I и II рода. Рассмотрены более сложные смешанные несимметричные условия II и III рода, что привело к определенным трудностям при построении итоговых решений.

Математическая постановка задач

и построение их решений

Приведем две задачи из этого класса:

Задача 1. Схема задачи приведена на рис. 1.

Рис. 1. Схема микроволнового нагрева угольного слоя с граничными условиями II и III рода

Fig. 1. Scheme of microwave heating of coal layer with boundary conditions of II and III order

Система уравнений микроволнового нагрева угольного слоя:

dT(х, t) _ л d2T (х, t) x

— a "г e ,

dt dx cp

T (x,0) _ f (x),

— _a.IT(0,t) - TC((t)],

dx

(S, t)

dx

(1.1) (1.2)

(1.3)

(1.4)

где T(x,t) - текущая температура, К; x - поперечная координата, м; t - время, с; a - температуропроводность, м2/с; c - удельная теплоемкость, Дж/кг-К; к - коэффициент поглощения в законе Бугера, 1/м; Я - теплопроводность угля, Вт/м-К; Tc1(t) - переменная по времени температура окружающей среды, К; q2(t) - нестационарная плотность теплового потока на нижней поверхности, Вт/м2; а1 - коэффициент теплоотдачи на верхней поверхности, Вт/м2-К; qv0 - максимальное значение внутреннего источника тепла, Вт/м3.

Данная система введением новой зависимой переменной:

и(х,Г) = Т(х,Г) - ТС1 (Г) - ^х2, (1.5)

1 ¿оЯ

приводится к более простой системе уравнений с однородными граничными условиями:

du( х, t) d2u( х, t)

■_ a—т^- + q(х, t),

dt

dx2

u (х,0) _ F(х),

du (0, t)

dx

_ axu(0, t),

xdTSA _ 0,

dx

(1.6)

(1.7)

(1.8)

(1.9)

где

q( х, t) _^ e - * - T>) - ^ х2 + S

cp 2Sa Sa

- новое распределение источников тепла; - первая производная по времени;

^(х) = / (х) - Г (0) - х2

- новое начальное распределение температуры.

Новую систему с неоднородным уравнением (1.6) и неоднородным начальным условием (1.7) можно разбить на две [1]. Решение будет иметь следующий вид: и(х^)=и!(х^)+и2(х^), для каждой части которого должны выполняться следующие уравнения: для и1

ды1 (х, t) д2щ (х, t) dt ~ дх2

щ (х,0) = 0,

и для u2

ды2 (х, t) du2 (х, t)

Ху

X

2) а%(у5) = —у - из этого уравнения вычисля-

X

ются собственные значения ук. Собственные функции выражаются следующим образом:

а,

(х) = С05(у„х) ^51П(уих),

хУп

(х) = ^-(^Уп )2',

где ап - члены ряда Фурье.

Этого достаточно, чтобы найти частные решения второй задачи (для и2):

= ane-(ГаУп )21 j^cos(y х) +XyLsin(y х)^

'J'

Далее составим ряд:

« = Е ^п е~{4~аГп )2' ( соь(у п х) + -р-вшу х) 1. п=1 ^ хУп )

Коэффициенты ап можно найти из начального условия (1.13):

F (х) = 11, а

а,

cos(yn х) + —1-sin(y, х)

ХУп

+ q( х, t), (1.10) (1.11)

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

|кп|| а

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

д( " дх2 ' и2(0 = F (х), (1.13)

при «старых» граничных условиях (1.8), (1.9).

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

Будем искать функции и1 в следующем виде:

«(х 0 = у (х)Ф,- (/).

Далее находим собственные значения у и собственные функции щ(х) однородной задачи, что приводит к двум уравнениям:

(О + (ТХу)>, (Г) = 0, (1.14)

у" (х) + у2у, (х) = 0. (1.15)

Решения для щ(х) в общем виде выражаются следующим образом:

х) = с1 б1п(у х) + с2 соэ(ух).

Подставляя их в граничные условия, находим:

1\ а1 1) С1 = С2'

(1 12) гДе ШЭ) - собственная функция

и,

со$(у пх) + —^Ш^х);

ХУп

г(£) - весовая функция (в данном случае равная 1). Пределы интегрирования: а=0, Ъ=8.

2 5( (X Л2

|Ы12 = | ^С0«(Упх) + X)J ^ =

1

4Уп

2уп5(1 + а2) + (1 - а2)sin(2y,,5) --4a(cos(2yn5) -1)

где а =

а

Далее выражаем искомую функцию в следующем виде:

™ 1

2(х, t) = e

п=1 Wn\\

S

(•Tay,, )21

cos(УnX) + —Чт(у,х)

ХУп

<i F^)

cos(yпЭ + i~^sin(y

ХУп

d Э

Для поиска решения системы (1.10)-(1.11) относительно щ воспользуемся применимой в данном случае теоремой Стеклова VI [1]. В итоге имеем:

™ 1 u1( хt) = Ег~Ц2

п=1 Ш

t S

<JJe"^)2(t-z) q(£z)

cos(Уnх) +^-^п(Упх) ХУп

а

cos(y,l) +T^sin(yn$

ХУп

d^dz.

Окончательное решение u(x,t) выглядит так:

да

i( x, t ) =

"=1 Л

t s xj je

0 0

cos(r„x) a

X7

~sin(7nx)

<-Tayn )2(t-г)

X q(4,T)

cos(YnO +-T^sin(7n^)

X7n

dÇdr ■

да i

n=1 Л

s

<j F (t)

- (4^7 n )2

cos(^x) ■

a

--^sin(7nx)

X7n

X7n

(1.16)

После возвращения к исходным переменным получаем итоговое решение системы (1.1)-(1.4):

Т(х, 0 = Т (() + ^ х2 + ^ ' 7 ^ 25Я

м 1 а

-=1 Ы ЯУп

xjje

0 0

-(4aïn )2(t-г)

- kt

ср

-Т (г) -

п (г) S2 , q2(T)

2SX

SX

*[cos(7nt) +-T^sin(YnÇ)]dÇdT ■

X7n

да i

n=1 Л

cos(Ynx) -

a

sin(7nx)

X7n

xj\f (t) - rci(0) - qSX t

a

xtœs^n^) +ir^sin(7nt)] dt

X7n

(1.17)

Ниже приведены расчетные графики распределения температуры по глубине слоя и в отдельных областях угольного слоя по времени, построенные по выведенным формулам. Рис. 2 демонстрирует распределение температуры по глубине слоя при значении числа Фурье, равном 0,35. При таком значении времени достаточно ограничится лишь первым членом ряда. Как видно из расчетов, наибольшая температура достигается в плоскости, несколько заглубленной в угольном массиве. Также с помощью полученного решения можно отследить временные изменения в интересующих нас сечениях (рис. 3): зону максимальной температуры, температуру в центре, на поверхности и др. На основе этих данных могут быть оценены термо-разрушающие напряжения, момент зажигания и др. параметры. Из графиков распределения темпе-

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

420

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

400

380

360

340

0.00 0.02 0.04 0.06 0.08 0.10

Рис. 2. Распределение температуры по глубине слоя

Fig. 2. Temperatures distribution on depth of layer

38>

Рис. 3. Временная зависимость температуры в разных сечениях угольного слоя: 1 - максимальная температура по слою; 2 - температура в центре; 3 - температура на поверхности

Fig. 3. Time dependence of temperature in different sections of the coal layer: 1 is the maximum temperature of the layer; 2 is the temperature in the center; 3 is the surface temperature

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

Задача 2. Схема задачи приведена на рис. 4.

Рис. 4. Схема микроволнового нагрева угольного слоя со смешанными граничными условиями II и III рода

Fig. 4. Scheme of microwave heating of coal layer with mixed boundary conditions of II and III order

e

2

Система уравнений микроволнового нагрева угольного слоя:

дТ (х, г) д2Г (х, г) -кх —= а—+ —еы, (2.1)

дг дх2 ср у '

Т (х ,0) = / (х), (2.2)

„ дТ(0, г)

-X дх = Ч1(г), (2.3)

хдТзх^ = (2[Т б, г) - Тс2 (г)]. (2.4)

Данная система упрощается введением новой подстановки:

qx(t)

5^

и(х,г) = Т(х,г)-ТС2(г)х-—--^, (2.5)

и приводится к системе уравнений относительно и(х^) с однородными граничными условиями

ди( х, г) д2ы( х, г)

= а—^^ + ч(х, г), (2.6)

dt дх2

u( х,0) = F (х),

xduiOii)=о,

(2.7)

(2.8)

(2.9)

2 5 Л

Ч(х,г) = -Т '(г)-х+—+-1 + -41-ср С2 X \ 25 2) Х8

- новое распределение источников тепла; ' - знак первой производной по времени;

где

дх

„ дТ(5,t) ,, ч X—^_iZ = a2 u(5, t), дх

„-*х хq1'(t) Г

qx(0)Г

х

5^

(х) = /(х) - ТС2(0)х+-+ ^)

- новое начальное распределение температуры.

Новую систему с неоднородным основным уравнением (2.6) и неоднородным начальным условием (2.7) можно разбить на две [1], если искать решение в виде суммы двух частных решений и(х^)=щ(х^)+и;(х^), для каждой части которого должны выполняться следующие уравнения: для и1

ди, (х, г) д2и, (х, г) —= а—+ Ч(х, г), (2.10)

и для u2

" дх2 u (х,0) = 0,

д^2 (х, t) d2u 2( х, t) dt ~ дх2 ''

u2 (t) = F (х),

(2.11)

(2.12) (2.13)

при «старых» граничных условиях (2.8), (2.9).

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

« С^г) = У(х)Ъ(г).

Далее находим собственные значения у и собственные функции щ(х) однородной задачи, что приводит нас к двум уравнениям:

Ъ (г) + (л/Ху)>, (г) = 0, (2.14)

Щ (х) + у2у. (х) = 0. (2.15)

Решения для щ(х) в общем виде выражаются следующим образом:

х) = С1 БШ(ух) + С2 СОБ(ух).

Подставляя в граничные условия, находим:

1) с:=0;

-X

2) ^(уб) = —у - из этого уравнения находятся

а2

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

Уп(х) = С0я(Упх);

Ъ (х) = ^)2 г,

где ап - постоянная интегрирования.

Этого достаточно, чтобы найти частные решения краевой задачи для и2:

«2п = апе~(ГаУп )2г С0«(Упх).

Далее составляем ряд:

ОТ

«2 =Е апе " (4~ау")2 г С08(у,х).

п=1

Коэффициенты ап можно найти методом ортогональности, используя начальное распределение и2(х,0):

F (х) = Е ап С0$(упх).

п=1

Это есть разложение заданной функции в тригонометрический ряд Фурье по синусам.

Коэффициенты разложения находятся по формуле

ап =7-^2 Ж,

где Ъп(^) - собственная функция ео8(ух); г(^) - весовая функция (в данном случае равна 1).

а =

0, b = 5, \\фп\\2 = J cos(ynх)2dх

= _L (У+sin(2y,5)

Уп ^ 2 4 Далее находим искомую функцию:

и2 (х, t) =

да 1 О

= е<4~ау")21С05(7их)[ F(^СО^)

"=1 ЫЛ

Для поиска решения системы (2.10)-(2.11) относительно щ воспользуемся применимой в данном случае теоремой Стеклова VI, согласно которой это решение может быть представлено в виде ряда Фурье по собственным функциям однородной задачи (4.1)-(4.2):

да

и\( х, t ) = X Ьп ^ (х).

п= 1

Для нахождения коэффициентов Ъ^) разложим заданную функцию qv(x,t) в ряд Фурье по собственным функциям ц/п(х):

члx >')=Е 11 Чп с )w„ (x )>

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

где

1 cä

q"(t)=FT 4v(x'' (x)dx

Далее подставляем выражения для q(x,t) и u1(x,t) в (3.1):

XЬП«П(х) =£аЬп(t)Д^п(х) +ХЯп (х).

п= 1 п=1 п =1

Затем, используя (2.15) и сравнивая коэффициенты перед \ifXx), находим уравнение для Ъ„(^:

Ьп ^)' + (4ау )2 Ьп ^ ) = я ^)

с начальным условием (2.11), из которого следует: Ъ.(0)=0 .

Решение линейного неоднородного дифференциального уравнения имеет вид:

К« ) = / /^)2(t-1) Яп (т)йт,

да 1

и1(x,t) = Хп2С08(^пх) х

п= 1 ||Ып||

хЩУ )2а-1) Я(£т)с08(7п£) йтй^. (2.16)

Промежуточное решение щ^Д):

да 1

м(х,t) = X—¡¡2С08(7пх) х п= 1 Ып||

х/о /О 1о е"^)2( т)cos(7n^)dтd^ +

е-)21 С08(7пх) х

п= 1 |Ып||

J F©cos^)

(2.17)

После возвращения к «старым» переменным получим итоговое решение системы (1.1)-(1.4):

xfj 'e jo jo

q,(') x2 S

T (x'') = T*(')+T(x " S 2) +

+^F^i2cos(r„x) x

1 \\F\

^ e"«- ТСг' (t) -cp

-(га?„ )2('-t) J - 4l (T) f. + Ü + S -

Я l 2S 2 J

4i(t) " Я6

x cos(7 п£)йтй^ +

. 1 -(.4äyn )2' ^Г-jTe ) x п=1 Fn||

7 © - th (0) -

f S

: cos(7nx) Ii q1(0) Г S^

Я l 2S 2J

x cos(jn^)d^.

(2.18)

Обсуждение

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

Заключение

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

x

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

Параметры СВЧ-нагрева, которые можно найти из полученных решений: • Исследовав на экстремум найденные температурные функции, можно найти координату и значение максимальной температуры.

СПИСОК ЛИТЕРАТУРЫ

1. Угольные технологии с использованием СВЧ-излучения / Вл.В. Саломатов, C.O. Сладков, С.Э. Пащенко, Вас.В. Саломатов // Горение и плазмохимия. - 2014. - T. 12. - № 2. -C. 77-96.

2. Диденко А.Н. СВЧ-энергетика: теория и практика. - М.: Наука, 2003. - 446 с.

3. Conversion of coal fly ash to zeolite utilizing microwave and ultrasound energies: a review / Syed Salman Bukharia, Jamshid Behina, Hossein Kazemiana, Sohrab Rohania // FueJ. - 15 January 2015.- V. 140. - P. 250-266.

4. Torrefaction of wheat and barley straw after microwave heating / Sangram Kishor Satpathya, Lope G. Tabilb, Venkatesh Medab, Satya Narayana Naika, Rajendra Prasada // Fuel. - 15 May 2014. - V. 124. - P. 269-278.

5. Faisal Mushtaqa, Ramli Mata, Farid Nasir Ani. A review on microwave assisted pyrolysis of coal and biomass for fuel production // Renewable and Sustainable Energy Reviews. - November 2014. - V. 39. - P. 555-574.

6. A comparative study using direct hydrothermal and indirect fusion methods to produce zeolites from coal fly ash utilizing singlemode microwave energy / Syed Salman Bukhari, Jamshid Behin, Hossein Kazemian, Sohrab Rohani // Journal of Materials Science. - December 2014. - V. 49. - Iss. 24. - P. 8261-8271.

7. Саломатов В.В., Сладков С.О., Пащенко С.Э. СВЧ-технологии в угольной энергетике // Инженерно-физический журнал. -2012. - № 85 (3). - С. 535-549.

8. Иваницкий В.В. О возможности применения сверхвысокочастотных методов нагрева материалов // Известия высших учебных заведений. Проблемы полиграфии и издательского дела. - 2000. - № 3-4. - С. 42-46.

9. Семенов А.С., Байбурин В.Б. СВЧ-энергия и ее применение. Особенности, оборудование, технологические процессы. - Саратов: Изд-во Саратовского ун-та, 1999. - 114 с.

10. Comparison and modeling of microwave tempering and infrared assisted microwave tempering of frozen potato puree / N. Se-yhun, Н. Ramaswany, G. Sumnu, S. Sahin, J. Ahmed // Journal of Food Engineering. - 2009. - № 92 (3). - Р. 339-344.

11. Chemkhi S., Zagrouba F., Bellagi A. Modeling and simulation of drying phenomena with rheological behavior // Brazilian Journal of Chemical Engineering. - 2005. - № 22 (2). - Р. 153-163.

12. A kinetic study of microwave and fluidized-bed drying of a Chinese lignite / Arash Tahmasebia, Jianglong Yua, Yanna Hana, Huan Zhaob, Sankar Bhattacharya // Chemical Engineering Research and Design. - January 2014. - V. 92. - Iss. 1. - P. 54-65.

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

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

• Другие технологические параметры по оптимизации, автоматизации.

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

13. The catalytic methanation of coke oven gas over Ni-Ce/Al2O3 catalysts prepared by microwave heating: Effect of amorphous NiO formation / Zhifeng Qin, Jun Ren, Maoqian Miao, Zhong Li, Ji-anying Lin, Kechang Xie // Applied Catalysis B: Environmental. - March 2015. - V. 164. - P. 18-30.

14. Pickles C.A., Gao F., Kelebek S. Microwave drying of a low-rank sub-bituminous coal // Minerals Engineering. - July 2014. -V. 62. - P. 31-42.

15. Fast deposition of porous iron oxide on activated carbon by microwave heating and arsenic (V) removal from water / Alp Yürüma, Züleyha Özlem Kocaba§-Ataklib, Meltem Sezena, Raphael Semi-atc, Yuda Yürüm // Chemical Engineering Journal. - 15 April 2014. - V. 242. - P. 321-332.

16. Рутинец Р.Г., Брыков С.К., Лохару Э.Х. Тепловые режимы при сверхвысокочастотном нагреве диэлектриков // Инженерно-физический журнал. - 1990. - № 59 (5). - С. 853-858.

17. Application of Carbon Nanocatalysts in Upgrading Heavy Crude Oil Assisted with Microwave Heating / Kewen Li, Binchi Hou, Lei Wang, Yi Cui // Nano Lett. - 2014. - № 14 (6). - P. 3002-3008. DOI: 10.1021/nl500484d.

18. Анфиногентов В.И. Математическое моделирование СВЧ-на-грева диэлектриков. - Казань: Изд-во Казанского гос. техн. унта, 2006. - 140 с.

19. Анфиногентов В.И., Гараев Т.К., Морозов Г.А. Об одной задаче теории СВЧ-нагрева диэлектриков // Вестник КГТУ им. А.Н. Туполева. - 2002. - № 3 (27). - С. 21-22.

20. Хайдурова А.А., Федчишин В.В., Коновалов Н.П. Микроволновая сушка бурых углей и повышение их технологических характеристик // Известия вузов. Проблемы энергетики. -2010. - № 1-2. - С. 31-35.

21. Хайдурова А.А., Коновалов Н.П., Федчишин В.В. Воздействие микроволновой энергии на бурый уголь для улучшения его технологических характеристик // Теплофизические основы энергетических технологий. - Томск: Изд-во ТПУ, 2009. -С. 109-113.

22. Угольный слой при микроволновом нагреве: аналитическое исследование при смешанных граничных условиях I и II рода / В.В. Саломатов, В.А. Карелин, А.С. Заворин, Вас.В. Саломатов, С.Э. Пащенко // Известия Томского политехнического университета. - 2015. - Т. 326. - № 11. - С.135-140. URL: http://www.lib.tpu.ru/fulltext/v/Bulle-tin_TPU/2015/v326/i11/14.pdf (дата обращения: 23.11.2015).

Поступила 10.03.2016 г.

Информация об авторах

Саломатов В.В., доктор технических наук, профессор, главный научный сотрудник Института теплофизики СО РАН.

Карелин В.А., магистрант Национального исследовательского Новосибирского государственного университета.

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

Саломатов В.В., кандидат физико-математических наук, ведущий программист Института теплофизики СО РАН.

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

UDC 536.24

HEAT TREATMENT OF THE COAL LAYER BY MICROWAVE ENERGY: AN ANALYTICAL STUDY IN CONDITIONS OF HEAT REMOVAL OF II AND III ORDER

Vladimir V. Salomatov1,

vvs@itp.nsc.ru

Vadim A. Karelin2,

vad2hen@mail.ru

Alexander S. Zavorin3,

zavorin@tpu.ru

Vasiliy V. Salomatov1

vvsalomatov@mail.ru

Sergey E. Pashchenko2,

pashenko-ecolog@yandex.ru

1 Institute of Thermal Physics SB RAS,

1, Academician Lavrentiev Avenue, Novosibirsk, 630090, Russia.

2 National Research Novosibirsk State University,

2, Pirogov street, Novosibirsk, 630090, Russia.

3 National Research Tomsk Polytechnic University, 30, Lenin Avenue, Tomsk, 634050, Russia.

Relevance of the research is related to the necessity of scientific study of microwave exposure technology parameters for coal, especially when drying, at thermal treatment, burning intensification and others. Experiments show that use of microwave radiation help reduce harmful emissions from coal burning, increase energy efficiency, reduce process time, etc. For optimal conditions of microwave treatment it is very popular to search for theoretical approaches, in particular, the engineering practice requires analytical solutions for heat treatment of coal array, which was done in the study.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

The aim of the research is to develop strict analytic solutions of the problems of coal layer heating by microwave radiation under the heat removal of the II and III order, which allow carrying out the parametric analysis of microwave exposure and searching for the most advantageous modes of heating coal layer.

Methods. The tasks of coal bed heating by microwave radiation with boundary conditions of the II and III order were set with a number of simplifications, such as persistence of thermal and electrical properties of the material, they are isotropy, one-dimensionality. This physical process itself is plane electromagnetic wave energy absorption and formation of heat sources in coal array, modeled on the Bou-guer law. The energy equation in the form of Fourier in this case is considered independently of the Maxwell equations of electrodynamics. A new dependent variable was determined. It modifies the boundary conditions into homogeneous, and simplifies the search of solutions. As a result, the basic system of equations is divided into two subsystems, and the final solution was found making a superposition of these two sub-tasks.

Results. The authors have found analytically strong dependences of temperature distribution over the cross section and time for carbon bed heated by microwave radiation, under heat removal of II and III order of boundary surfaces. These solutions are the basis for optimum parameters of microwave heating technology, a tool to reduce the spread of temperature field on the desired temperature distribution within the layer. They help determine quickly: maximum temperature inside a coal array, coordinates of its location, thermo-destructive voltage, limits of discharged heat fluxes and other.

Key words:

The Fourier heat equation, coal, microwave energy, the Buger law, mixed boundary conditions, thermal regime, heat removal.

REFERENCES

1. Salomatov Vl.V., Sladkov S.O., Pashchenko S.E., Salomatov Vas.V. Microwave technologies in coal power engineering. Gore-nie i plazmokhimiya, 2014, vol.12, no. 2, pp. 77-96. In Rus.

2. Didenko A.N. SVCh-energetika: teoriya i praktika [SHF-energet-ics: theory and practice]. Moscow, Nauka Publ., 2003. 446 p.

3. Syed Salman Bukharia, Jamshid Behina, Hossein Kazemiana, Sohrab Rohania. Conversion of coal fly ash to zeolite utilizing microwave and ultrasound energies: A review. Fuel, 15 January 2015, vol. 140, pp. 250-266.

4. Sangram Kishor Satpathya, Lope G. Tabilb, Venkatesh Medab, Satya Narayana Naika, Rajendra Prasada. Torrefaction of wheat and barley straw after microwave heating. Fuel, 15 May 2014, vol. 124, pp. 269-278.

5. Faisal Mushtaqa, Ramli Mata, Farid Nasir Ani. A review on microwave assisted pyrolysis of coal and biomass for fuel production. Renewable and Sustainable Energy Reviews, November 2014, vol. 39, pp. 555-574.

6. Syed Salman Bukhari, Jamshid Behin, Hossein Kazemian, Sohrab Rohani. A comparative study using direct hydrothermal and indirect fusion methods to produce zeolites from coal fly ash utilizing

single-mode microwave energy. Journal of Materials Science, December 2014, vol. 49, Iss. 24, pp. 8261-8271.

7. Salomatov Vl.V., Sladkov S.O., Pashchenko S.E. Microwave technology in coal energy. Journal of Engineering Physics and Thermophysics, 2012, no. 85 (3), pp. 535-549. In Rus.

8. Ivanitskiy V.V. O vozmozhnosti primeneniya sverkhvysokochas-totnykh metodov nagreva materialov [Possibility of using of microwave heating methods of materials]. Proceedings of the institutions of higher education. Issues of the graphic arts and publishing, 2000, no. 3-4, pp. 42-46.

9. Semenov A.S., Bayburin V.B. SVCh-energiya i eeprimenenie. Oso-bennosti, oborudovanie, tekhnologicheskie protsessy [Microwave energy and its application. Features, equipment, technological processes]. Saratov, Publ. House of Saratov University, 1999. 114 p.

10. Seyhun N., Ramaswany H., Sumnu G., Sahin S., Ahmed J. Comparison and modeling of microwave tempering and infrared assisted microwave tempering of frozen potato puree. Journal of Food Engineering, 2009, no. 92 (3), pp. 339-344.

11. Chemkhi S., Zagrouba F., Bellagi A. Modeling and simulation of drying phenomena with rheological behavior. Brazilian Journal of Chemical Engineering, 2005, no. 22 (2), pp. 153-163.

12. Arash Tahmasebia, Jianglong Yua, Yanna Hana, Huan Zhaob, Sankar Bhattacharya. A kinetic study of microwave and fluidiz-ed-bed drying of a Chinese lignite. Chemical Engineering Research and Design, January 2014, vol. 92, Iss. 1, pp. 54-65.

13. Zhifeng Qin, Jun Ren, Maoqian Miao, Zhong Li, Jianying Lin, Kechang Xie. The catalytic methanation of coke oven gas over Ni-Ce/Al2O3 catalysts prepared by microwave heating: Effect of amorphous NiO formation. Applied Catalysis B: Environmental, March 2015, vol. 164, pp. 18-30.

14. Pickles C.A., Gao F., Kelebek S. Microwave drying of a low-rank sub-bituminous coal. Minerals Engineering, July 2014, vol. 62, pp. 31-42.

15. Alp Yuruma, Zuleyha Ozlem Kocaba§-Ataklib, Meltem Sezena, Raphael Semiatc, Yuda Yurum. Fast deposition of porous iron oxide on activated carbon by microwave heating and arsenic (V) removal from water. Chemical Engineering Journal, 15 April 2014, vol. 242, pp. 321-332.

16. Ruginets R.G., Brykov S.K., Lokharu E.Kh. Teplovye rezhimy pri sverkhvysokochastotnom nagreve dielektrikov [Thermal regimes at microwave heating of dielectrics]. Journal of Engineering Physics and Thermophysics, 1990, no. 59 (5), pp. 853-858.

17. Kewen Li, Binchi Hou, Lei Wang, and Yi Cui. Application of Carbon Nanocatalysts in Upgrading Heavy Crude Oil Assisted with Microwave Heating. Nano Lett., 2014, no. 14 (6), pp. 3002-3008. DOI: 10.1021/nl500484d.

18. Anfinogentov V.I. Matematicheskoe modelirovanie SVCh nagre-va dielektrikov [Mathematical modeling of microwave heating of dielectrics]. Kazan, Kazan State University Publ. House, 2006. 140 p.

19. Anfinogentov V.I., Garaev T.K., Morozov G.A. Ob odnoy zadache teorii SVCh nagreva dielektrikov [A problem in the theory of microwave dielectric heating]. Vestnik KGTU im. AN. Tupoleva, 2002, no. 3 (27), pp. 21-22.

20. Khaydurova A.A., Fedchishin V.V., Konovalov N.P. Microwave drying of brown coal and increasing its technological characteristics. Izvestiya vysshykh uchebnykh zavedeniy. Problemy energe-tiki, 2010, no. 1-2, pp. 31-35. In Rus.

21. Khaydurova A.A., Fedchishin V.V., Konovalov N.P. Vozdeystvie mikrovolnovoy energii na bury ugol dlya uluchsheniya ego tekhnologicheskikh kharakteristik [Effects of microwave energy on brown coal to improve its processing characteristics]. Teplofi-zicheskie osnovy energeticheskikh tekhnologiy [Thermophysical basics of energy technologies]. Tomsk, TPU Publ., 2009. pp. 109-113.

22. Salolomatov Vl.V., Karelin V.A., Zavorin A.S., Salomatov Vas.V., Pashchenko S.E. Coal layer under microwave heating: analytical study under mixed boundary conditions I and II of the genus. Bulletin of the Tomsk Polytechnic University. Geo assets Engineering, 2015, vol. 326, no. 11, pp. 135-140. Available at: http://www.lib.tpu.ru/fulltext/vZBulle-tin_TPU/2015/v326/i11/14.pdf (accessed 23 November 2015).

Received: 10 March 2016.

Information about the authors

Vladimir V. Salomatov, Dr. Sc., professor, chief research worker, Institute of Thermal Physics SB RAS. Vadim A. Karelin, graduate student, National Research Novosibirsk State University.

Alexander S. Zavorin, Dr. Sc., professor, head of the department, National Research Tomsk Polytechnic University.

Vasiliy V. Salomatov, Cand. Sc., programming supervisor, Institute of Thermal Physics SB RAS.

Sergey E. Pashchenko, Cand. Sc., senior research worker, National Research Novosibirsk State University.

i Надоели баннеры? Вы всегда можете отключить рекламу.