CC BY
7УДК: 537.86:536.33
МЕТОД ОЦЕНКИ ПЕРЕПАДОВ ТЕМПЕРАТУРЫ И ВЛАГОСОДЕРЖАНИЯ ПРИ ЭЛЕКТРОМАГНИТНОЙ СУШКЕ ОБРАЗЦА ПРОИЗВОЛЬНОЙ ФОРМЫ
А.М. Афанасьев, Г.А. Попов, Б.Н. Сипливый
Анатолий Михайлович Афанасьев *, Глеб Александрович Попов
Кафедра Информационной безопасности, Институт приоритетных технологий, Волгоградский государственный университет, пр. Университетский, 100, Волгоград, Российская Федерация, 400062 E-mail: [email protected]*, [email protected]
Борис Николаевич Сипливый
Кафедра теоретической физики и волновых процессов, Физико-технический институт, Волгоградский государственный университет, пр. Университетский, 100, Волгоград, Российская Федерация, 400062 E-mail: [email protected]
Аналитическими методами исследован процесс удаления влаги из образца произвольной формы, находящегося в рабочей камере установки для сушки электромагнитными волнами. Образец обдувается потоком воздуха и подвергается воздействию двух видов излучений, имеющих большую и, соответственно, малую по сравнению с размерами образца глубину проникновения; для повышения равномерностей обдува и тепловой обработки осуществляется вращательное движение образца. В качестве математической модели сушки использованы уравнения теории тепломассопереноса А.В. Лыкова в среде с постоянными теплофизическими коэффициентами. С помощью преобразований трех видов - осреднения по времени оборота, осреднения по координатам и перехода к асимптотике по времени -исходные уравнения приведены к виду, содержащему лишь числовые характеристики функций. На основе этих соотношений предложены формулы для величин, названных характерными перепадами температуры и влагосодержания, которые характеризуют степени неоднородности соответствующих полей. Введенные перепады зависят от теплофизических параметров материала, характерного размера образца (для шара - диаметр), а также от мощностей генераторов электромагнитных волн с большой и с малой глубиной проникновения. В случае шара характерные перепады дают точные значения разности температур и разности влагосодержаний между поверхностью шара и его центром. Результаты работы позволяют организовать электромагнитную сушку, в которой при высокой интенсивности испарения влаги с поверхности образца не происходит перегрева материала (возникает за счет большого перепада температуры) или его разрушения от механических деформаций (возникает за счет большого перепада влагосодержания).
Ключевые слова: уравнения Лыкова; электромагнитная сушка; образец произвольной формы; характерные перепады температуры и влагосодержания
METHOD OF EVALUATING CHANGES IN TEMPERATURE AND MOISTURE CONTENT USING ELECTROMAGNETIC DRYING OF SAMPLE OF ARBITRARY SHAPE
A.M. Afanasiev, G.A. Popov, B.N. Siplivy
Anatoly M. Afanasiev*, Gleb. A. Popov
Department of Information Security, Institute of Priority Technologies, Volgograd State University, University
ave., 100, Volgograd, 400062, Russia
E-mail: [email protected]*, [email protected]
Boris N. Siplivy
Department of Theoretical Physics and Wave Processes, Physico-Technical Institute, Volgograd State University, University ave., 100, Volgograd, 400062, Russia E-mail: [email protected]
Analytical methods were used to study the process of removing moisture from the sample of arbitrary shape, located in the working chamber of the installation for drying by electromagnetic waves. The sample is blown by the air flow and exposed to two types of radiation, having a large and, accordingly, small compared to the size of the sample penetration depth. To increase the uniformity of blowing and heat treatment, the sample is rotated. The equations of heat and mass transfer theory are used as a mathematical model of drying Lykov in a medium with constant thermal coefficients. Using transformations of three types - averaging over turnaround time, averaging over coordinates, and transition to time asymptotics - the original equations are reduced to a form that contains only numerical characteristics offunctions. On the basis of these relations, formulas for the values called characteristic temperature and moisture content differences, which characterize the degree of heterogeneity of the corresponding fields, are proposed. The introduced differences depend on the thermophysical parameters of the material, the typical size of the sample (for the ball - diameter), as well as on the power of generators of electromagnetic waves with a large and small penetration depth. In the case of a ball, the characteristic differences give the exact values of the temperature difference and the moisture content difference between the surface of the ball and its center. The results of the work allow to organize electromagnetic drying, in which at a high intensity of moisture evaporation from the surface of the sample there is no overheating of the material (due to a large temperature drop) or its destruction from mechanical deformations (due to a large difference in moisture content).
Key words: Lykov equations, electromagnetic drying, arbitrary shape sample, characteristic temperature and moisture content differences
Для цитирования:
Афанасьев А.М., Попов Г.А., Сипливый Б.Н. Метод оценки перепадов температуры и влагосодержания при электромагнитной сушке образца произвольной формы. Изв. вузов. Химия и хим. технология. 2019. Т. 62. Вып. 10. С. 131-139
For citation:
Afanasiev A.M., Popov G.A., Siplivy B.N. Method of evaluating changes in temperature and moisture content using electromagnetic drying of sample of arbitrary shape. Izv. Vyssh. Uchebn. Zaved. Khim. Khim. Tekhnol. 2019. V. 62. N 10. P. 131-139
ВВЕДЕНИЕ
риеи - пластины, цилиндра с прямоугольным сечением, прямоугольного параллелепипеда, цилиндра Вопросу моделирования процессов тепло- с круговым сечением и шара - алгоритмы расчета и массопереноса при сушке влажных материалов перепадов уже разработаны авторами, соответ-электромагнитным излучением уделяется большое ственно в статьях [7-11]. В настоящей статье точ-внимание в литературе по химическим техноло- ные формулы для перепадов из этих работ впервые гиям [1 -6]. Одна из проблем организации режимов обобщаются в виде приближенных формул для обсушки с наперед заданными свойствами состоит в разцов с произвольной геометрией. следующем. Для повышения скорости сушки при- Схема установки
ходится повышать интенсивность излучения. Это Одна из применяемых на практике схем
сопровождается увеличением перепадов темпера- сушки изображена на рисунке (например, такую
туры и влагосодержания по глубине материала. Но конструкцию имеют микроволновые бытовые печи
большие перепады температуры могут привести к с грилем). Рабочей камерой устройства является
перегреву материала, а большие перепады влагосо- объемный резонатор, в котором расположены два
держания - к его разрушению от механических де- генератора электромагнитных волн, сверхвысоко-
формаций. Ввиду этого ставится задача об органи- частотного (СВЧ) и инфракрасного (ИК) диапазо-
зации таких режимов, в которых при максимально нов. Обычно можно считать, что СВЧ излучение,
возможной интенсивности сушки перепады темпе- обладая большой проникающей способностью, ин-
ратуры и влагосодержания снижены до безопас- дуцирует источники тепла, распределенные по
ного уровня. Чтобы такой режим осуществить, сле- всему объему образца; в противоположность
дует располагать алгоритмами, аналитическими этому, источники тепла, вызванные поглощением
или численными, для расчета указанных величин. ИК излучения, распределяются в тонком поверх-
В настоящей статье будут рассматриваться анали- ностном слое материала. Находящийся под воздей-
тические алгоритмы. Для тел с простейшей геомет- ствием этих двух видов излучения влажный обра-
зец с помощью специального устройства приводится во вращение, чем достигается равномерность термообработки. Для отвода окружающих образец водяных паров его обдувают потоком воздуха.
Рис. Установка для сушки электромагнитным излучением: 1 - влажный образец; 2 - вращающийся столик; 3 - объемный резонатор; 4 - СВЧ излучатель; 5 - ИК излучатель; 6 - воздушный поток
Fig. Set up for drying by electromagnetic radiation: 1 - wet sample; 2 - rotating table; 3 - volume resonator; 4 - microwave emitter; 5 - IR emitter; 6 - air flow
Математическая модель электромагнитной сушки
Обратимся к процессам, происходящим в рабочей камере установки для сушки излучением. Влажный образец - тело V с границей Е - обдувается воздушным потоком и подвергается воздействию электромагнитных волн. Переменные точки в области V и на поверхности Е будем обозначать соответственно, как Ми Р. За счет поглощения электромагнитной энергии внутри образца возбуждаются процессы нагревания материала и испарения жидкости. Вследствие этого возникают тепловые потоки и потоки влаги как внутри образца, так и с его поверхности. Внутренний перенос тепла происходит по закону теплопроводности Фурье, а внутренний перенос влаги - по закону влагопро-водности А.В. Лыкова: часть влаги перемещается в виде жидкости, часть - в виде пара, и вызывается это движение градиентами влагосодержания и температуры (диффузия и термодиффузия). Движение влаги может происходить также и под действием градиента давления влажного воздуха внутри пор (массоперенос по закону фильтрации Дарси). Однако эта составляющая потока влаги вносит заметный вклад в общий поток лишь при температуре материала, близкой к температуре кипения воды, и учитывать ее здесь мы не будем.
Обозначим
T = T(M, т), U = U(M, т) -температура тела [°С] и его влагосодержание
SU/ St = am5-V2T + am -V2U; M eV; X-ST/dn = r(1 -y)J + Q - S; P eE;
[кг/кг] как функции точки М и времени т [с]. Положение точки будем задавать ее координатами в системе, жестко скрепленной с вращающимся телом. При условии, что теплофизические характеристики материала могут считаться постоянными, начально-краевая задача для расчета полей Т и и будет формулироваться следующим образом [12, 13]: ср0(дТ/дт) = Я-У2Т + гур0(<5и/дт)+Ж; МеV; (1)
(2) (3)
8-дТ/ дп + ди/ дп = J| (атр0 ); Р еЕ; (4)
Т(М, 0) = Л(М); и(М, 0) = /2 (М); т = 0. (5)
Здесь (1) и (2) - уравнения распространения тепла и влаги в области V; (3) и (4) - краевые условия на поверхности Е; (5) - начальные условия при т = 0, где/1(М) и/2(М) - некоторые заданные функции координат. В этих уравнениях: Ж [Вт/м3] и [Вт/м2] - плотности объемных и поверхностных источников тепла, возникающих за счет работы генераторов электромагнитных волн СВЧ и ИК диапазонов соответственно (представление о поверхностных источниках тепла введено в статье авторов [13]); с, ро, у, Я, От, 5 - теплофизические характеристики материала (удельная теплоемкость [Дж/(кг°С)], плотность в сухом состоянии [кг/м3], критерий испарения, коэффициент теплопроводности [Вт/(м°С)], коэффициент диффузии влаги [м2/с], относительный коэффициент термодиффузии влаги [1/°С]); г - удельная теплота парообразования воды [Дж/кг]; ддп - символ производной по направлению внутренней нормали к границе Е; Q и J - интенсивность теплообмена [Вт/м2] и интенсивность массообмена [кг/(м2с)] поверхности Е с воздушной средой, которые в каждой точке Р рассчитываются по формулам
(6) (7)
В этих формулах Тв и ф - температура и влажность воздуха вблизи точки Р за пределами пограничного слоя; о - постоянная Стефана-Боль-цмана [Вт/(м2°С4)]; А - коэффициент теплового излучения; О(Т) - функция, моделирующая зависимость относительного парциального давления насыщенного водяного пара от его температуры Т; Т1=273 °С - постоянная; Оу(Р) и ат(^) - коэффициенты тепло- и массообмена поверхности образца с воздушной средой (соответственно [Вт/(м2°С)] и [кг/(м2с)]), которые сложным образом зависят от скорости воздуха V вблизи точки Р. Вид функции G(Т) можно найти в [14]; там же приводятся и фор-
Q(P, t) = «w [T (P, T)-Tв ]+o2[(T (P, T)+T )4 -(TB + Ti )4 J (P, x) = «m [G(T (P, х))-ф. G(TB )].
мулы для коэффициентов aw(V) и am(V) для некоторых частных случаев.
Поставим задачу дать приближенный анализ электромагнитных процессов и процессов теп-ломассопереноса в изображенной на рисунке ситуации. С этой целью подвергнем исходные уравнения преобразованиям, с помощью которых из уравнений последовательно будут исключены координаты и время, и конечным итогом будут новые уравнения, содержащие лишь числовые характеристики всех функций. Рассмотрим подробнее эти процедуры.
Осреднение по времени оборота, по координатам и переход к асимптотике по времени
Пусть Ат - время одного оборота, AV -объем тела, АЕ - площадь его поверхности, и F(M, т) - некоторая функция координат и времени. Операции осреднения этой функции по времени оборота, по объему тела и по его поверхности введем следующим образом:
1 т+Лт/2
F (M, т)= — JF(M, Cd,
ЛТ т-Лт/2
<F}(т) = ^ fF(M, T)dVM '
< F )(т) = -^ ff F (P, T)dS p •
ЛЕ E
В результате выполнения этих операций все непрерывные функции координат и времени, которые мы рассматриваем, будут заменены на новые непрерывные функции.
Будем считать время Ат малым по сравнению с характерными временами изменения полей температуры и влагосодержания. Тогда в результате выполнения операции осреднения по времени оборота получим, что функции Т, U, их производные по времени, их лапласианы в области V и нормальные производные на поверхности Е существенных изменений не приобретут, т.е. они сохранят свой прежний вид.
Этот результат мы отразили в таблице; так же будем поступать и с последующими результатами.
Переходя к плотности СВЧ потерь W, примем, что время оборота, в отличие от ситуации с полями температуры и влагосодержания, является большим по сравнению со временем релаксации электромагнитных процессов. Кроме того, будем считать, что способность материала поглощать энергию электромагнитных волн в каждой точке пространства не зависит от температуры и влагосо-держания в этих точках. Тогда плотность потерь W(M, т) в каждый момент времени будет опреде-
ляться только текущим расположением образца относительно СВЧ излучателя [15]. Ясно, что зависимость от времени будет периодической, с периодом, равным времени оборота. Выполнив осреднение, получим новую функцию Ж(М). Аналогичным образом из функции S(P, т) в результате осреднения получится функция 5 (Р).
Обсудим теперь условия тепло- и массооб-мена на поверхности образца. Интенсивность теплообмена Q, согласно формуле (6), определяется теплообменом конвекцией по закону Ньютона и теплообменом излучением по закону Стефана-Больцмана, а интенсивность массообмена J, согласно закону испарения Дальтона (7), зависит от перепада парциального давления водяного пара по толщине пограничного слоя. Осредняя функцию J(Р,т) по времени оборота, мы учтем, что функция G(Т(Р,т)) за это время не изменяется, переменными оказываются только характеристики воздушного потока вблизи точки Р: его влажность ф, температура Тв и скорость V, от которой зависит коэффициент массообмена ат. Выполнив осреднение, получим функцию J (г (Р, т)). Аналогичным образом из функции Q(Р,т) после осреднения получится функция Q(г (р, т)).
Таким образом, за время оборота функции Т и и не изменяются; функции Ж и S изменяются за счет изменения ориентации тела относительно излучателей; функции Q и J изменяются за счет изменения ориентации тела относительно воздуховодов.
Осреднение по времени оборота позволило исключить из всех функций осцилляции, связанные с вращением образца. Следующее за ним осреднение по координатам, т.е. по поверхности или по объему, в зависимости от области задания функции, исключит из функций переменные точки М и Р, и оставит в них только время т (разумеется, если только оно осталось после осреднения по времени оборота). Результаты представлены в таблице и дополнительных пояснений не требуют.
Рассмотрим теперь операцию, позволяющую исключить из осредненных по координатам функций время. Как известно, сушка влажного материала начинается с периода прогрева, за которым следует период постоянной скорости сушки [12]. Исследование периодов постоянной скорости имеет большое значение для практики, поскольку во многих случаях именно в эти промежутки времени из материала удаляется основная часть влаги. Характерными признаками периода постоянной скорости являются стационарность поля темпера-
туры и квазистационарность поля влагосодержа-ния (величина V уменьшается во времени по линейному закону со скоростью, одинаковой во всех точках и во все моменты). В математической физике решения такого рода называют квазистационарными, или установившимися [16]. Формально поля Т и V становятся такими лишь асимптотически, т.е. при т^да, но фактическое время их установления всегда является конечным. Таким образом, для этого режима мы имеем:
Т(М,т) = ГШ(М), и(М,т)=иш(м)+0. (8)
Здесь V» = ди/дт = const < 0 - скорость сушки, а Uco(M) есть распределение влагосодержа-ния в момент т = 0, который в приведенной формуле мы принимаем за начало установившегося режима. Благодаря стационарности поля Т, отвечающие асимптотическому режиму интенсивности тепло- и массообмена перестают изменяться во времени. Например, функция (Q}(t), как это отражено в таблице, преобразуется в постоянную (Q)
. Аналогичным образом решается вопрос и об асимптотике функции (J }(т).
Таблица
Результаты осреднения функций и асимптотические значения этих величин
Исходные функ- Осреднение по времени обо- Осреднение по координа- Асимптотика по вре-
ции рота там мени
ST/ St(M, т) cT/ 5т (M, т) (ST/ ST}(T) 0
dU/ Sr(M, т) SU/ 5т (M, т) (SU/ Sx}(x) V„
V2T (M, т) V2T (M, т) (V2T }(т) (V2T }м
V2U (M, т) V2U (M, т) (V2U }(т) (V2U }m
ST/ Sn (P, т) ST/ Sn (P, т) (дг/ а«}(т) (ST/SnL
SU/ Sn (P, т) dU/ Sn (P, т) (SU/ Sn}(x) (SU/ Sn}M
W (M, т) W (m ) (W} (W}
S (P, т) S (p) (S} (S}
Q(T (P, т), т) Q(T (p, т)) (Q}(T) (Q}„
J (T (P, т), т) J (T (P, т)) ( J }(т) (J }„
Выполнив над уравнениями системы (1)-(4) преобразования, перечисленные в таблице, получим новую систему следующего вида:
о = Х-<У2т >„ + гур0^+<ж >; (9)
ат5-<у2г>ю + ат-<у2и^; _ (Ю) X - <8Т/8п>ю = г(1 - у) - < 1 >ю + <0Ю - <5>; (И) З-^Т/дП^ -Кди/дп^ =<1> ^(атРо ). (12) Подчеркнем, что в отличие от исходных уравнений, где присутствуют функции координат и времени, здесь мы имеем соотношения между числовыми (интегральными) характеристиками этих функций.
Характерные перепады АТ* и Д V* Согласно плану работы, прежде всего нас будут интересовать перепады температуры АТ и влагосодержания Ди между внутренними областями материала и его поверхностью. Уточним эти понятия. Прежде всего, будем говорить только об асимптотических режимах, когда, согласно фор-
муле (8), указанные перепады от времени не зависят. Примем, далее, что
ДТ = Т(М)-Т(Р), где меК, Ре! -- какие-нибудь две точки, соответственно внутри тела и на его поверхности. Аналогичным образом определим и перепад Ди. Ясно, что введенные перепады будут зависеть от выбора точек М и Р, что затрудняет использование этих величин как некоторых параметров, с помощью которых можно судить о степени неоднородности полей Т и V. Переопределим тогда понятия перепадов так, чтобы при новой формулировке возникшее затруднение уже не имело места. С этой целью обратимся к решенной авторами задаче об электромагнитной сушке однородного шара [11]. Если плотности объемных и поверхностных источников тепла являются постоянными, а условия тепло- и массооб-мена на поверхности шара не зависят ни от координат, ни от времени, то в режиме постоянной скорости сушки стационарное распределение темпе-
ратуры и квазистационарное распределение влаго-содержания будут параболическими. Например, для поля влагосодержания мы будем иметь и (г, т)=ахг2 +а2 + Утт .
Здесь г - расстояние до центра шара, а1 и а2 - постоянные, V» - введенная выше скорость сушки, а время, как и в формуле (8), отсчитывается от начала установившегося режима. Перепадом влагосодержания между центром шара и его поверхностью г = Я назовем величину
Д и = и (0, т) - и (Я, т) = -ахЯ 2.
Чтобы дать обобщение перепада А и для тела произвольной формы, придадим этой формуле другой вид. Вычислим нормальную производную поля и на поверхности шара (нормаль внутренняя):
ди/ д«|2 = -ди/ дг\г=к = -2ахЯ.
Сравнивая это с формулой для А и, получим, что
Ди = (£/ 4)-ди/ д^.
Здесь В = 2Я - диаметр шара. Эта формула связывает существование перепада влагосодержа-ния между центром и поверхностью шара с наличием вполне определенного градиента влагосодер-жания на поверхности. Имея в виду такой смысл полученного результата, и обращаясь теперь к телу произвольной формы, мы определим величину Аи* - характерный перепад влагосодержания -формулой
Ди* = {Б */4><ди/ дп)^. (13)
Здесь В* - характерный размер тела (для шара - диаметр), а <ди/ди)|Е - среднее значение
нормальной производной ди/дп на поверхности тела. Для тела в форме шара это соотношение имеет точный смысл, как об этом говорилось выше.
Аналогичной формулой введем и характерный перепад температуры:
ДТ* = (Б */4><дТ/дя)|Е. (14)
Отметим, что по смыслу всех величин, будем иметь А и* > 0 и АТ* > 0, если влагосодержание и, соответственно, температура в среднем возрастают при переходе от поверхности к внутренним областям.
Формулы для характерных перепадов
Чтобы найти перепады по формулам (13) и (14), мы должны определить входящие в них средние нормальные производные полей Т и и на поверхности тела. Обратимся для этого к уравнениям (9)-(12). Решающим обстоятельством для дальнейшего является прием, позволяющий исключить
входящие в эти уравнения лапласианы полей. Пусть ф(М) - произвольное достаточно гладкое скалярное поле. Для участвующей в наших рассуждениях области V с границей Е (нормаль внутренняя) запишем интегральную теорему для лапласиана [16]:
Ш^2Ф ^^ = -#|фйЕ.
V Е дП
Здесь слева мы имеем умноженный на объем ЛV средний лапласиан <У2ф), а справа -умноженную на площадь ЛЕ среднюю нормальную производную <дф/ди). Применяя эту теорему к полям Т и и, в асимптотическом режиме будем иметь:
ДЕ де
<у2Т)» = -Т^Т/ дп)да; <У2и )х = -—<ди/ дп)да.
ДV ДV
Исключим из уравнения (9) величину V» с помощью (10), и в получившееся уравнение подставим на место лапласианов полей полученные выше их выражения через нормальные производные. Рассматривая новое уравнение совместно с (12), мы будем иметь два уравнения с двумя неизвестными величинами -нормальными производными полей Т и и. Разрешая эту систему, получим:
<дТ/ дп)ю =1Г^ <Ж )-уг< J )0 X V ДЕ
-I
ХГ 4 "" ДЕ
Здесь мы ввели новую постоянную материала Х = у + Х1(атр0-).
Вместо величин ) и < S), которые не могут быть измерены непосредственно, введем Рсвч и Рик - потребляемые установкой мощности соответствующих генераторов и общую потребляемую мощность Р = Рсвч + Рик. Очевидно, что
<Ж) = Рсвч/ДV, <5) = Р^/ДЕ.
Выражая в (15) величину ) через Рсвч, и подставляя получившиеся выражения для нормальных производных в (13) и (14), получим итоговые формулы для характерных перепадов:
<ди1дп)ж=-{хг^)ж-Д^ Ж)]. (15)
ДТ * = — (с, -уг{ J >„-Де)
X
(16)
ди
* = J >„-ДЕ-Рсвч)
Здесь обозначено С*=В */(4 ДЕ) - геометрическая характеристика образца.
РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ
Результат проведенного исследования заключается в формулах (16). Они связывают характерные перепады температуры и влагосодержания
АТ* и А и* с параметрами В*, ЛЕ (геометрические характеристики образца); от, ро, 5, X, у (теплофизи-ческие характеристики материала); Рсвч (выходная мощность СВЧ генератора); < J) (интенсивность
сушки).
Чтобы воспользоваться этими формулами, необходимо знать интенсивность сушки <J) , для
которой у нас нет алгоритма расчета, и, следовательно, она должна быть получена в результате измерений. Это обстоятельство можно обойти следующим образом. Подставляя в (11) значение <дТ/дп)х из (15), получим после преобразований
уравнение, имеющее смысл баланса энергетических потоков на поверхности образца:
г< J )„= Р/ ДЕ. (17)
Заметим теперь, что в режимах с высокой интенсивностью сушки потери энергии за счет теплообмена с воздушным потоком оказываются обычно небольшими по сравнению с затратами энергии на испарение воды [17, 18]. Это означает, что <0Ж << г< J)ж, и тогда из предыдущей формулы получим
г<J)Ю•ДЕ« Р. (18)
При выполнении этого условия формулы (16) принимают следующий приближенный вид:
лг* = C* (^свч -yP); ли* = C*5 (xp - рвч) • (19)
В этих формулах переменная (J у исключена, и роль заменяющей ее переменной выполняет общая потребляемая мощность Р.
Заметим, что перепады AT* и AU* получаются независящими от характеристик воздушного потока. Это согласуется с точным решением задачи о сушке шара [11]: варьирование характеристик воздушного потока приводит лишь к повышению или понижению температуры поверхности шара, а перепады температуры и влагосодержания между поверхностью шара и его центром остаются при этом неизменными.
Важно отметить следующее обстоятельство. Рассмотрим множество режимов, когда Рсвч и Рик одновременно изменяются, но так, что Рсвч + Рик = Р = const.
Согласно (18), интенсивность сушки (J у
будет в любом таком режиме одной и той же, в то время как перепады AT* и AU*, подсчитанные по формулам (19), различными. Пользуясь этим обстоятельством, можно варьированием мощности генераторов подобрать режим, в котором при заданной интенсивности сушки характерные пере-
C * 5,
пады температуры и влагосодержания будут соотноситься нужным для нас образом. Другой способ установления необходимого соотношения между перепадами изложен авторами в работе [19]. Вместо двух генераторов, с большой и с малой глубиной проникновения, там используется только один генератор электромагнитных волн, а выбор оптимального режима сушки осуществляется согласованным варьированием мощности этого генератора и его частоты.
Как уже отмечалось, в случае шара рассмотренная здесь задача о расчете асимптотических характеристик сушки будет иметь точное аналитическое решение, и вследствие этого полученные выше приближенные формулы приобретут точный смысл. В частности, формулы (19) будут определять перепады температуры и влагосодер-жания между поверхностью шара и его центром. В свою очередь, аналитические решения задач сушки для тел с простой геометрией имеют многочисленные прямые и косвенные опытные подтверждения [7, 10-12, 14, 17]. Значит, можно надеяться, что для тел, форма которых не слишком отличается от шарообразной, формулы (19) окажутся вполне пригодными для инженерных оценок.
ВЫВОДЫ
Сформулирована начально-краевая задача для расчета полей температуры и влагосодержания в образце произвольной формы, находящемся в рабочей камере установки для сушки электромагнитным излучением. Образец расположен на вращающемся столике, он обдувается воздушным потоком и находится под воздействием излучений генераторов электромагнитных волн ИК и СВЧ диапазонов. Исследование начально-краевой задачи проведено при следующих упрощающих предположениях: теплофизические и электромагнитные характеристики влажного материала считаются постоянными; фильтрационное движение влаги, происходящее по закону Дарси, не учитывается; время оборота образца является небольшим по сравнению с характерными временами изменения полей температуры и влагосодержания и большим по сравнению со временем релаксации электромагнитных процессов; сушка рассматривается лишь в установившемся режиме, когда ее скорость почти не изменяется (режим постоянной скорости сушки). Для приближенного анализа процесса сушки выполнено осреднение исходных уравнений по времени оборота и по координатам, а также осуществлен переход к асимптотике по времени; в результате этих операций уравнения для искомых функций
были заменены соотношениями для их числовых (интегральных) характеристик. На основе этих соотношений, и с привлечением точного решения задачи об электромагнитной сушке шара, были введены параметры, названные характерными перепадами температуры и влагосодержания, которые позволяют судить о степени неоднородности данных полей. Контроль за характерными перепадами позволяет при интенсивных способах сушки избежать перегрева материала или его разрушения от
ЛИТЕРАТУРА
1. Афанасьев А.М., Никишова А.В., Сипливый Б.Н. Переходные процессы при сушке конвекцией и сушке инфракрасным излучением. Изв. вузов. Химия и хим. технология. 2017. Т. 60. Вып. 10. С. 94-101. Ш1: 10.6060Лсс!20176010.5568.
2. Рудобашта С.П., Зуева Г.А., Зуев Н.А. Влияние термодиффузии на кинетику осциллирующей инфракрасной сушки. Изв. вузов. Химия и хим. технология. 2016. Т. 59. Вып. 4. С. 83-87.
3. Натареев О.С., Кокина Н.Р., Натареев С.В. Теплопе-ренос в процессе конвективной сушки влажного материала. Изв. вузов. Химия и хим. технология. 2015. Т. 58. Вып. 2. С. 67-72.
4. Натареев С.В., Венкин Е.Н., Натареев О.С. Массопе-ренос в процессе сушки влажного материала при конвективном подводе теплоты. Изв. вузов. Химия и хим. технология. 2012. Т. 55. Вып. 3. С. 104-108.
5. Рудобашта С.П., Зуева Г.А., Карташов Э.М. Тепломассо-перенос при сушке сферической частицы в осциллирующем электромагнитном поле. Теор. основы хим. технологии. 2016. Т. 50. № 5. С. 539-550. Б01: 10.7868/80040357116050109.
6. Рудобашта С.П., Карташов Э.М., Зуев Н.А. Тепломассопе-ренос при сушке в осциллирующем электромагнитном поле. Теор. основы хим. технологии. 2011. Т. 45. № 6. С. 641-647.
7. Афанасьев А.М., Сипливый Б.Н. Обобщение понятия температуры мокрого термометра на случай электромагнитной сушки. Инж.-физ. журн. 2016. Т. 89. № 3. С. 610-616.
8. Афанасьев А.М., Сипливый Б.Н. Теория электромагнитной сушки: асимптотическое решение начально-краевой задачи для прямоугольной области. Физика волн. проц. ирадиотех. системы. 2012. Т. 15. № 1. С. 77-83.
9. Афанасьев А.М., Сипливый Б.Н. Асимптотические распределения температу-ры и влагосодержания при электромагнитной сушке образца, имеющего форму прямо-угольного параллелепипеда. Изв. вузов. Электромеханика. 2012. № 3. С. 3-8.
10. Афанасьев А.М., Сипливый Б.Н. Теория электромагнитной сушки: асимпто-тическое решение начально-краевой задачи для цилиндра. Теор. основы хим. технологии. 2014. Т. 48. № 2. С. 222-227. Б01:10.7868/80040357114020018.
11. Афанасьев А.М., Сипливый Б.Н. Задача о сушке шара электромагнитным излучением. Инж.-физ. журн. 2013. Т. 86. № 2. С. 322-330.
12. Лыков А.В. Теория сушки. М.-Л.: Энергия. 1968. 471 с.
13. Афанасьев А.М., Сипливый Б.Н. Концепция поверхностных источников тепла в теории сушки электромагнитным излучением. Изв. вузов. Электромеханика. 2017. Т. 60. № 2. С. 13-20.
14. Афанасьев А.М., Сипливый Б.Н. О краевых условиях массообмена в виде законов Ньютона и Дальтона. Инж.-физ. журн. 2007. Т. 80. № 1. С. 27-34.
механических деформаций. Результаты работы могут быть использованы как средство инженерных оценок при разработке технологических линий и устройств для нагревания и сушки электромагнитными волнами.
Исследование выполнено при финансовой поддержке РФФИ и Администрации Волгоградской области в рамках научного проекта № 18-48340001 р_а.
REFERENCES
1. Afanasiev A.M., Nikishova A.V., Siplivy B.N. Transients during drying convection drying and infrared radiation. Izv. Vyssh. Uchebn. Zaved. Khim. Khim. Tekhnol. 2017. V. 60. N 10. P. 94-101 (in Russian). DOI: 10.6060/tcct.20176010.5568.
2. Rudobashta S.P., Zueva G.A., Zuev N.A The influence of diffusion on the kinetics of the oscillating infrared drying. Izv. Vyssh. Uchebn. Zaved. Khim. Khim. Tekhnol. 2016. V. 59. N 4. P. 83-87 (in Russian).
3. Natareev S.V., Kokina N.R., Natareev O.S. Heat transfer in the convective drying process of moist material. Izv. Vyssh. Uchebn. Zaved. Khim. Khim. Tekhnol. 2015. V. 58. N 2. P. 67-72 (in Russian).
4. Natareev S.V., Venken E.N., Natareev O.S. Mass transfer in drying process of moist material at convective supply of heat. Izv. Vyssh. Uchebn. Zaved. Khim. Khim. Tekhnol. 2012. V. 55. N 3. P. 104-108 (in Russian).
5. Rudobashta S.P., Zueva G.A., Kartashov E.M. Heat and mass transfer during drying of a spherical particle in an oscillating electromagnetic field. Teor. Osnovy Khim. Tehnol. 2016. V. 50. N 5. P. 539-550 (in Russian). DOI: 10.7868/S0040357116050109.
6. Rudobashta S.P., Kartashov E.M., Zuev N.A. Heat and mass transfer during drying in an oscillating electromagnetic field. Teor. Osnovy Khim. Tehnol. 2011. V. 45. N 6. P. 641-647 (in Russian).
7. Afanasiev A.M., Siplivy B.N. Generalization of the Wet-Bulb Temperature Notion to the Case of Electromagnetic Drying. J. Eng. Phys. Thermophys. 2016. V. 89. N 3. P. 620-626. DOI: 10.1007/s10891-016-1419-3.
8. Afanasiev A.M., Siplivy B.N. Electromagnetic drying theory: asymptotic solution of the initial boundary value problem for a rectangular region. Fizika Voln. Prots. Radiotekh. Sistemy. 2012. V. 15. N 1. P. 77-83 (in Russian).
9. Afanasiev A.M., Siplivy B.N. The asymptotic distribution of temperature and moisture content using electromagnetic drying of the sample, having the shape of a rectangular parallelepiped. Izv. Vyssh. Uchebn. Zaved. Elektromekhanika. 2012. N 3. P. 3-8 (in Russian).
10. Afanasiev A.M., Siplivy B.N. Theory of electromagnetic drying: an asymptotic solution of an initial-boundary value problem for a cylinder. Theor. Found. Chem. Eng. 2014. V. 48. N 2. P. 206-211.
11. Afanasiev A.M., Siplivy B.N. Problem of electromagnetic drying of a sphere. J. Eng. Phys. Thermophys. 2013. V. 86. N 2. P. 340-348.
12. Lykov A.V. Theory of drying. M.-L.: Energiya. 1968. 471 p. (in Russian).
13. Afanasiev A.M., Siplivy B.N. The concept of surface heat sources in the theory of drying by electromagnetic radiation. Izv. Vyssh. Uchebn. Zaved. Elektromekhanika. 2017. V. 60. N 2. P. 13-20 (in Russian).
15. Шимони К. Теоретическая электротехника. М.: Мир. 1964. 775 с.
16. Тихонов А.Н., Самарский А.А. Уравнения математической физики. М.: Наука. 1966. 724 с.
17. Афанасьев А.М., Сипливый Б.Н. Оптимизация процесса электромагнитной сушки капиллярно-пористых материалов. Изв. вузов. Электромеханика. 2006. N° 5. С. 3-10.
18. Календерьян В.А., Бошкова И.Л., Волгушева Н.В. Кинетика микроволновой сушки сыпучего органического материала. Инж.-физ. журн. 2006. Т. 79. №№ 3. С. 123-127.
19. Афанасьев А.М., Сипливый Б.Н. Оптимизация сушки электромагнитным излучением: аналитическое исследование проблемы. Изв. вузов. Химия и хим. технология. 2019. Т. 62. Вып. 6. С. 69-76. БО1: 10.6060ЛуШ201962ф.5870.
20. Падохин В.А., Зуева Г.А., Кокурина Г.Н., Кочкина Н.Е., Федосов С.В. Комплексное математическое описание тепло-и массопереноса в процессе сушки неограни-ченного тела цилиндрической формы аналитическими методами теории теплопроводно-сти. Теор. основы хим. технологии. 2015. Т. 49. № 1. С. 54-64. БО1: 10.7868/80040357115010108.
14. Afanasiev A.M., Siplivy B.N. On the boundary conditions of mass transfer in the form of Newton and Dalton laws. Inzh.-Fiz. Zhurn. 2007. V. 80. N 1. P. 27-34 (in Russian).
15. Shimoni K Theoretical electrical engineering. M.: Mir. 1964. 775 p. (in Russian).
16. Tikhonov A.N., Samarskiy A.A. Equations of mathematical physics. M.: Nauka. 1966. 724 p. (in Russian).
17. Afanasiev A.M., Siplivy B.N. Optimization of electromagnetic drying of capillary-porous materials. Izv. Vyssh. Uchebn. Zaved. Elektromekhanika. 2006. N 5. P. 3-10 (in Russian).
18. Kalender'jan V.A., Boshkova I.L., Volgusheva N.V. Kinetics of microwave drying of bulk organic material. Inzh.-Fiz. Zhurn. 2006. V. 79. N 3. P. 123-127 (in Russian).
19. Afanasiev A.M., Siplivy B.N. Optimization of drying by electromagnetic radiation: analytical study of the problem. Izv. Vyssh. Uchebn. Zaved. Khim. Khim. Tekhnol. 2019. V. 62. N 6. P. 69-76 (in Russian). DOI: 10.6060/ivkkt201962fp.5870.
20. Padokhin V.A., Zueva G.A., Kokurina G.N., Kochkina N.E., Fedosov S.V. Complex mathematical description of heat and mass transfer in the drying process of an unlimited cylindrical body by analytical methods of the theory of thermal conductivity. Teor. Osnovy Khim. Tekhnol. 2015. V. 49. N 1. P. 54-64 (in Russian). DOI: 10.7868/S0040357115010108.
Поступила в редакцию 14.03.2019 Принята к опубликованию 07.08.2019
Received 04.03.2019 Accepted 07.08.2019
