CC BY
54
УДК 532.546+536.24: 666.3.047
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ТЕПЛО- И МАССОИЕРЕНОСА ПРИ АДСОРБЦИОННО-КОНТАКТНОЙ СУШКЕ КЕРАМИЧЕСКИХ МАТЕРИАЛОВ (МАТЕМАТИЧЕСКАЯ МОДЕЛЬ)
Дорняк Ольга Роальдовна д.т.н., доцент
ФГБОУВПО «Воронежская государственная лесотехническая академия», Воронеж, Россия
Усачев Александр Михайлович к.т.н., доцент
Шмитько Евгений Иванович д.т.н., профессор
ФГБОУВПО «Воронежский государственный архитектурно-строительный университет»
Аль Сарраджи Салах Хусейн Мела аспирант
ФГБОУ ВПО «Воронежская государственная лесотехническая академия», Воронеж, Россия
В рамках механики многофазных систем построена математическая модель процесса тепло- и массопереноса при адсорбционно-контактной сушке керамического кирпича
Ключевые слова: МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ, ТЕПЛО- И МАССОПЕРЕНОС, АДСОРБЦИОННО-КОНТАКТНАЯ СУШКА, КАПИЛЛЯРНО-ПОРИСТЫЕ МАТЕРИАЛЫ
Введение
UDC 532.546+536.24: 666.3.047
MATHEMATICAL MODELLING OF HEAT AND MASS TRANSFER DURING
ADSORPTION-CONTACT DRYING OF CERAMIC MATERIALS (MATHEMATICAL MODEL)
Domyak Olga Roaldovna Dr.Sci.Tech., associate professor
Voronezh State Forestry Academy, Voronezh, Russia
Usachev Aleksandr Mikhailovich Cand.Tech.Sci., associate professor
Shmitko Evgeniy Ivanovich Dr.Sci.Tech., professor
Voronezh State University of Architecture and Civil Engineering, Voronezh, Russia
Al' Sarradzhi Salakh Hussein Mela postgraduate student
Voronezh State Forestry Academy, Voronezh, Russia
The mathematical model of heat and mass transfer during adsorption-contact drying of ceramic bricks within the framework of the mechanics of multiphase systems is constructed
Keywords: MATHEMATICAL MODELLING, HEAT AND MASS TRANSFER, ADSORPTION-CONTACT DRYING, CAPILLARY AND POROUS MATERIALS
Математическое моделирование процессов сушки капиллярно-пористых материалов развивается в течение многих десятилетий по нескольким направлениям. Основные из них связаны как с изучением фундаментальных закономерностей теплофизических процессов, так и с решением практических проблем прогнозирования результатов обезвоживания материалов [1,2]. Развитие математического моделирования процессов сушки сдерживается недостатком
знаний о механизмах взаимодействия фаз, а также о структуре материала, изменяющейся во взаимосвязи с процессами тепло- и массопереноса. Как отмечено в работе [2], моделирование процесса сушки капиллярно-пористых тел все еще остается близким к искусству.
Математические модели развиваются применительно к традиционным промышленным методам сушки (конвективному, радиационно-кондуктивному, кондуктивному, сублимационному, токами высокой частоты и другим), а также к новым способам. Оригинальный метод адсорбционно-контактной сушки (АКС) был разработан в Институте катализа имени Г.К. Борескова СО РАН. Метод был разработан применительно к сушке гранул различных материалов, в том числе, включая гранулы взрывчатых веществ. Основная идея метода заключается в том, что при контакте пористых гранул, насыщенных жидкостью (донора), с сухими пористыми гранулами (акцептором), возможен массообмен с переносом флюида от донора к акцептору. Интенсивность переноса существенно зависит от градиента химических потенциалов между донором и акцептором [1]. Массоперенос возможен в широком температурном диапазоне, включая обычную комнатную температуру.
Авторами [3] разработан метод АКС глиняного сырца (в работах [3-4] он называй диффузионно-контактным). Следует отметить, что метод АКС особенно интересен, как метод сушки термолабильных материалов, свойства которых ухудшаются в условиях обычной термической сушки [5-6]. Условия АКС могут обеспечить эффективную сушку при невысоких температурах. Изготовление керамического кирпича из высокочувствительных глин (к ним относятся, в том числе, глины семилукского месторождения), требует значительных материальных затрат для создания мягкого режима термообработки, позволяющего снизить число образцов, которые растрескиваются на этапе сушки сырца.
влагонепроницаемый колпак высушиваемый сырец
рамки
Рисунок 1 - Схема
адсорбционно-контактной сушки
глиняного сырца.
Схема сушки контактно-диффузионным способом представлена на рисунке 1. Отформованный сырец укладывается на посредник-подложку, выполненную из влагоемкого материала, накрывается влагонепроницаемым колпаком и помещается в сушильную камеру [4]. Посредник с тыльной стороны подвергается конвективному высушиванию движущимся теплоносителем с заданными температурно-влажностными параметрами. Таким образом, влага из сырца диффундирует в объем материала-посредника, влажность которого поддерживается на уровне, безопасном для объекта высушивания, а затем происходит удаление влаги уже из посредника [3-4].
Исследования [6-9] позволили определить основные механизмы массообмена при адсорбционно-контактной сушке. Если донор насыщен раствором компонента, нелетучего при температуре массообмена, то по газовой фазе происходит перенос паров только растворителя до момента установления равновесного состояния. Через зону контакта донора и акцептора возможен перенос преимущественно растворителя до установления равновесия. Жидкофазный перенос происходит вследствие перепада давления, вызванного различием кривизны менисков жидкости в пределах жидкофазного домена (при условии образования "жидкостного" мостика в зоне контакта), а также по
смачивающим пленкам за счет градиента расклинивающего давления.
В работах, представленных в [1], проведено детальное экспериментальное исследование кинетики адсорбционно-контактной сушки пористых материалов, обоснованы уравнения кинетики влагопереноса применительно к АКС гранул с деформируемой и недеформируемой структурой (А.Д.Симонов, Н.А.Языков, В.Н.Пармон). На основе современного подхода к описанию текстурно-чувствительных массообменных процессов, протекающих при сушке капиллярно-пористых материалов, определены основные механизмы массообмена при АКС (В.Б.Фенелонов).
При изучении процессов тепло- и массопереноса при АКС объектов с достаточно большими размерами, для которых характерно неравномерное и нестационарное развитие полей теплофизических переменных, использование кинетических уравнений ограничено. В данной работе сформулирована математическая модель адсорбционно-контактной сушки керамических материалов в рамках механики многофазных систем [10], следуя [11]. Данный подход не предполагает модельного представления пористой структуры материала, но требует описания динамики полей целого ряда текстурных характеристик пористой системы материалов. Методология механики гетерофазных систем позволяет совершенствовать модели тепло- и массопреноса в капиллярно-пористых средах, учитывая различные механизмы физико-химического взаимодействия фаз.
Математическая модель тепло- и массопереноса при адсорбционно контактной сушке материала
Пусть выполнены следующие предположения. Материал образцов донора и акцептора изотропен. Начальная влажность по сечению образцов распределена в общем случае неравномерно. Жидкая фаза практически несжимаема, ее плотность постоянна. Среднее значение начальной влажности сырца может быть как выше, так и ниже предела гигроскопичности материала. Процессы деформирования материала вследствие тепловых воздействий в данной работе не
Научный журнал КубГАУ, №89(05), 2013 года рассматриваются.
Рисунок 2 - Расчетная схема Учитывая наличие плоскостей симметрии для полей теплофизических переменных, рассмотрим, как объект исследования - тело, состоящее из половины кирпича, разрезанного вдоль направления движения сушильного агента, и фрагмента платформы, боковые вертикальные грани которого являются плоскостями симметрии для полей влагосодержания и температуры в посреднике (рисунок 2). Это составное тело имеет поверхность разрыва для всех текстурных параметров, а также теплофизических параметров твердой фазы. Поверхность разрыва совпадает с плоскостью контакта донора (О) и акцептора
Математическая формулировка задачи, которая ставится для описания различных процессов в многофазной системе, включает уравнения, записанные для теплофизических переменных, усредненных по объемам всех имеющихся в системе фаз - 1 (газообразной), 2 (жидкой ), 3 (твердой). При моделировании процессов сушки капиллярно-пористой системы с широким распределением
радиусов пор можно учитывать различные механизмы переноса воды, которые определяются формой связи влаги с материалом твердой фазы. В данной работе задача, гидродинамики поставлена для нескольких типов воды.
В уравнениях состояния и сохранения, записанных ниже, знак усреднения по фазе опущен. В уравнениях символ/ обозначает переменную/ отнесенную к объему фазы с номером / , и усредненную по этому объему, т.е. </,>„ следуя обозначениям [10].
Газовая фаза состоит из двух компонент - неконденсирующегося газа и водяного пара. Примем предположение об идеальности паровой и газовой составляющих. Параметры, относящиеся к первой компоненте, обозначим нижним индексом ко второй - IV. Плотность, концентрация, давление и удельная внутренняя энергия компонент и парогазовой смеси определяются следующим образом [10]:
Л п Л п р\ ~Р\у+Р\ё. р! . р! . 9 9 9 (1)
Р1*= рГ/1%. Р1„=Р1°уВД„. 9 9 (2)
= С^Т\. и1у = . 9 9 (3)
Р\ =р"ВД. -в, =хвъ.+(\-г)в^ (4)
Здесь В - индивидуальная газовая постоянная, Дж/(кг-К); су - теплоемкость при постоянном объеме, Дж/(кгК)\ р - плотность, кг/м3, Т - температура, К; р -давление, Па; % - концентрация пара в парогазовой смеси; и - удельная внутренняя энергия, Дж/кг. Верхний индекс "о" у величины плотности указывает на ее истинное значение в отличие от приведенного значения рг=рг°аг (аг - объемная концентрация /-фазы).
Значения скорости паровой и газовой компонент у1у и у1ё могут быть
различны. Для их описания введены среднемассовая скорость смещений элементарных макрообъемов первой фазы Vj и диффузионные скорости пара и газа - wlv и wlg. Относительное движение компонент определяется законом бинарной диффузии Фика [10]:
Параметр О - коэффициент бинарной диффузии, зависящий в общем случае от температуры газа, (м21 с). Уравнения сохранения массы для парогазовой смеси и неконденсирующегося газа имеют вид [10]:
где j -поток массы пара, обусловленный фазовыми переходами, отнесенный к единице времени и единице площади, кг (м2 с)\ - удельная поверхность
раздела 1 и 2 фазы, лг1; УА, У° - области пространства, занятые акцептором и донором соответственно; - объемное содержание /-ой фазы в материале. Случай />0, соответствует испарению, у<0 - конденсации.
Уравнение движения и теплопроводности парогазовой фазы при сделанных предположениях можно записать следующим образом [10]:
V1 = Xviv +(l-5C)vig .
(5)
Wig - vlg vl Wlv = Vlv _ Vl
(6)
(7)
d(pfai)
dt
_ s{2,xgVA + V(p1°a1v1 ) = sl2j' Sl2~[sD7XeVD
(8)
+ V[p1°a1(l-x)(v1 +wg)] = 0
dt
(9)
аг +a2 +a3 =1
(10)
1 _ а 1^1 у1
^13^(01)
*13 =
К1А3,хеУА у = К^,хеУ° . 1
х¥А,хе¥А Ц)0 хеУВ
1 ’ • (П)
ср1р°ах^-У(а1Х1УТ1)+а1В1Т1-^- + су\$^(Т1\112-Т1) + (212+(2\т,
ср\~%ср\+^ р\ё (2^ = ~Т}) ; у = 152,3.
Здесь К13 - коэффициент проницаемости 1-ой фазы в 3-ей при полном насыщении, м2; Т, - относительная фазовая проницаемость для /-ой фазы; 0І -насыщенность объема порового пространства /-ой фазой; ср - теплоемкость при постоянном давлении, Дж/(кгК)\ X -коэффициент теплопроводности, Вт/(м-К); |1
- вязкость, Па-с; л,, - удельная площадь поверхности раздела /-ой и /-ой фаз, лг1;
аь-
]1 - коэффициент теплоотдачи /-ой фазы через межфазную поверхность 11р
Он~
Вт/(м2-К); р - интенсивность теплообмена через межфазную поверхность X,, к /-ой фазе, Дж/(м3-с). Верхние индексы А и I) - обозначают принадлежность к акцептору и донору соответственно. Переменные с индексом Еу относятся к границам раздела фаз / и у.
Задача для парогазовой фазы является сопряженной и должна решаться для объемов поровой системы донора, акцептора и пространства под колпаком. В работе [12] построена балансовая математическая модель термодинамического состояния парогазовой смеси под влагоизолирующим колпаком применительно к исследуемому процессу АКС.
Отметим, что большинство параметров уравнений (1-13) одинаковы для
донора и акцептора, за исключением тех, которые непосредственно связаны со структурными особенностями материала. Среди них - удельная поверхность капиллярно-пористого материала и коэффициент проницаемости.
Жидкая фаза - вода, заполняющая поры, свободная и связанная. Для более точного описания гидродинамических процессов в рамках того же метода усреднения по отдельным фазам, введем дополнительное усреднение по объемам свободной и связанной воды. Причем в объеме связанной воды выделим объем воды в смачивающих граничных слоях и тонких прослойках, где вода не имеет контакта с газовой фазой [11].
Для некоторого мгновенного микропараметра жидкой фазы х'2, усредняемого по жидкому макрообъему (IV2, имеем:
\х2/2 а2свх2св а2сшх2сш "^а2т.слЛ:2т.сл
Здесь х2св, х2см, х2тсл - рассматриваемый параметр, усредненный по объемам свободной воды, воды в смачивающих пленках и тонких слоях соответственно:
О^св У2а/(1V2^ С^т.сл (''зт.с:/(Хсм (1 V2cм/(1 К2.
Введем дополнительно параметры объемного содержания в образце свободной воды /7?еи=с1 К2еи/(1V, воды в тонких капиллярах ттл=6У2тл^^У и воды в смачивающих пленках тсм=дУ2смМУ. Тогда
{х2)? = (т2свх2св т2смх2см ~>гт2 т.сл л:2т.сл )
х ао
^2св ^Х2св>2св, ^2см ^Х2см>2см? ^2т.сл <'Х2іхл>2т.сл-
где а2 - объемное содержание воды в материале. Очевидно,
т2св = <^2св^2- т2см =<^2см^2- т2т.сл =^2т.сл^2
Насыщенность пористой структуры материала свободной водой, водой, содержащейся в тонких слоях и смачивающихся пленках, представляется в виде следующих выражений:
^2св "*"^2см "*"^2т.сл =®2
а _ т2св а _ т2см а _ т2т.с
^2св — ^2см — ^2т.с —
а2 +OL\ . OL2 +а1 • ОС2 +OL\
9 9
Уравнения сохранения массы. Рассмотрим уравнения сохранения массы для каждого из видов воды, содержащейся в капиллярно-пористом теле. Учтем, что вода в тонких слоях практически неподвижна, вследствие малости здесь градиента давления жидкости по сравнению с силами сопротивления, во-первых, и не испытывает фазовых переходов, во-вторых [11]:
Vo = 0
2т. сл ^ (14)
а2т.сл — 1 — а2св ~ а2см (15)
Для воды свободной и в смачивающих пленках получим, полагая
Р2 = const .
рО[Э(м2св) + V(m2cBv )] = -JSl2cB
8t . (16)
Р2[Э<”2СМ) + V (ш2см v2 )] = -j-ц 2см
91 , (17)
^12 = ^12св "*"^12см (18)
Если сложить почленно уравнения (16-17), учитывая (14-15), получим уравнение сохранения жидкой фазы, усредненное по всему занимаемому ею объему [10]:
'~2 = (*2}, =-7('»2св^2св +т1см^2 ш +«>2т слетел)
где 2 - вектор скорости
жидкой фазы, усредненные по ее объему.
Уравнение изменения количества движения. Уравнение изменения количества движения запишем без учета конвективных членов, динамических эффектов фазовых переходов и массовых сил. Для тонких смачивающих пленок, перетекающих по пористой системе образцов под действием градиента расклинивающего давления, после усреднения уравнений движения по объему Усм имеем:
Р°2>П2шд^ = ->п^Р2ш -^^см—У3> + у*( ^ )
31 к23^2ф2ш) .
К23
К*,хеУА
К^,хєУВ 2
х¥А,хеУА
(21)
*7 о к) )
То =2ЦоЄ0 2см О Vу 7 ^2см ^ у ^^2см/
2см 2 2см ; 2 ; ^=1,2,3. (22)
Второй член справа от знака равенства в уравнении (20) соответствует силе
хк
сопротивления фильтрационному переносу жидкой фазы в пористой системе;
- вектор напряжений, Па; К2з - коэффициент проницаемости 2-ой фазы при полном насыщении, м2.
Если образец находится в поле внешних сил, создающих градиент давления, то уравнения для определения кинематических характеристик течения свободной воды могут быть записаны в виде аналогичном (20):
к] , А7 е\' = —СУ -V? + У,у=! )
^=2^2» ; 2“ 2 7 * 2ов'’;ко=1АЗ. (24)
Если влиянием внешних СИЛ МОЖНО пренебречь, ТО /?2св=СОП81 и У2еи=0.
Отметим, что после почленного сложения уравнений движений воды трех типов
не удается получить, как в случае уравнения неразрывности, уравнения движения
жидкой фазы, усредненного по всему объему, занимаемому этой фазой.
Уравнение теплопроводности жидкой фазы используется в форме [10]
ВТ
с2Р2а2^~ = ^(а2^2^^2) ~ С28\2]^2\у'[-> ~ ^2) + б12 + 0.
(25)
Давление воды. Усредненное давление в жидкой фазе представляется в
виде:
(р 2)2 = (т2свР2св т2смР2см т2т.спР2т.сп )
а2
Давление свободной воды в недеформируемой системе, когда влиянием массовых сил можно пренебречь, при изотермических условиях постоянно и равно давлению окружающей среды рс. При деформировании
капиллярно-пористого материала давление свободной воды в ненасыщенных областях по-прежнему будет оставаться равным рс. В зонах полного насыщения жидкой фазой давление свободной воды должно определяться в процессе решения задачи из условия совместного деформирования твердого скелета и воды (равенство нормальных напряжений на границе жидкости и деформируемых стенок поровых каналов):
Рз =-^<®'з'>3
3 , (26)
где о3 - тензор напряжений твердой фазы.
Расклинивающее давление в смачивающих пленках можно определить,
следуя методике работы [13], позволяющей это сделать в рамках равновесной
термодинамики двухфазных многокомпонентных систем. Выражения для р2см
вычисляется из равенства химических потенциалов жидкости и пара по формуле Кельвина с помощью уравнения изотерм сорбции
Р2см = РваЛТ1,12) + Ь Р1У
РяМ^ІЛі) (27)
а = /'4(Г,ф) и = /°(Т,ф)
где Я - универсальная газовая постоянная Длс (моль-К) \ у - молярный объем
т
жидкости, мъ!молъ; и - влагосодержание, кг/кг; - температура границы фаз 1 и 2, К; - давление насыщенного пара у плоской поверхности жидкости, Па; Р\у - давление паровой компоненты газообразной фазы, Па; ср - относительная влажность воздуха.
Принято, что сопротивление течению в тонких слоях существенно выше имеющегося здесь градиента давления. В этом случае величина Ръ&А не используется при решении задачи переноса воды.
Для твердой фазы, деформации которой в данной работе не а- = сот1
рассматриваются ^ , а уравнение теплопроводности можно записать в
виде:
дТ
сЗРзаЗ~- = біЗ +023 +^(а3^3^^з)
(28)
сз
с^,хе¥А
РзА,хє¥а Рз° °,хеУ°
(29)
Математическая модель включает уравнения сохранения на межфазных поверхностях. На границе раздела жидкость-пар (в поверхностной фазе Е12) в общем случае следует учитывать неравновесность фазовых переходов, связанную с тем, что количество пара, испаряющегося с поверхности раздела фаз зависит от кинетических возможностей паровой фазы. Кинетика неравновесных фазовых переходов описывается уравнением Герца-Кнудсена-Ленгмюра без учета скачка
температур в граничном кнудсеновском слое пара [10]:
_ к[р$ш(^21^12^ _
Здесь к - коэффициент аккомодации.
Уравнение сохранения энергии в поверхностной Е12 фазе без учета тепловой инерции Е12 фазы и влияния искривленности межфазных поверхностей на ее температуру можно представить в виде
512./^ = ОС 12^12 (^1 _ ^11x12 ^ + а12‘у12 (^2 (31)
Условия теплообмена на остальных поверхностях раздела фаз могут быть записаны в виде балансовых соотношений [10]:
513аз 1(^213 ~тз) = _513а1з(^2пз ~ТХ). 52заз2(Г22з -Т3) = -523а23(Г223 -Т2)
9
(32)
Краевые условия на внешних границах. Система уравнений в частных производных, уравнений состояния и соотношений на границе фаз (1-32) для
переменных Тъ Т2 ,Т3, У| , р\, х, У2см, а2 должна быть дополнена граничными и
начальными условиями. Краевые условия для составного объекта АКС поставим для четырех групп внешних границ (рисунок 2) - ^ - и и ^4 и ^5 ?
ГГ> = £*2 и ^3 и = ^3 и ^2 “ А и ^1 и Л6 и ГЖ
р
На множестве ^ справедливы условия симметрии для всех усредненных теплофизических переменных:
, /-1,2,3; ; ; ; . (33)
Множество объединяет внешние грани донора, под
влагоизолирующим колпаком. Здесь практически мгновенно устанавливается
равновесная влажность соответствующая температуре Тк и относительной влажности паровоздушной смеси под колпаком ср^ (нижний индекс к относится к параметрам парогазовой смеси в пространстве под влагоизолирующим колпаком). Примем, что
я, дт<
2см
дп
= «,Гс№
г,
I)
г
Тк)
и
/=1,2,3:
I
г =1к Р1\г =рк
^Р2см
дп
V
2см
= 0
«2|г = /°(Тк’(Рк)аЗР°3/Р2
(34)
Здесь и ниже —коэффициенты теплоотдачи /-ой фазы с окружающей средой
на границе Г, (/=А, Э, Ак), Вт (мгК). Верхние индексы т и п обозначают касательную и нормаль к границе области.
На границах акцептора, соответствующих объединению областей ^, практически мгновенно устанавливается равновесная влажность мА, соответствующая текущей температуре сушильного агента Тс и его относительной влажности фс. Для внешних поверхностей тела посредника имеем
я,- дт<
дп
= ^А(Т,
Гл
Гл
Тс)
/=1,2,3:
_ПдЛ_
дп
-1с)
дрх
дп
= Р1Г(А
п
2см
др2
СМ
г'А
дп
Гл ^2
V
2см
= 0
(35)
где (3,г1 (лг1) и ур (м/с) - коэффициенты массоотдачи /-ой фазы с окружающей средой на границе Г| (/-А, Э, Ак).
Часть поверхности акцептора
оказывается под влагоизолирующим
колпаком. Здесь краевые условия могут быть представлены следующим образом
-А*
дЪ
дп
Г
Ак
г
-тк)
Ак
/=1,2,3:
%
Г -Кк рх\т =рк
Ак ■ '1Ак
2см
Ф2
см
ГА
дп
К23^2^в2см)
1 Ак
V
2см
о
- Ак
а
1Л(Тк,(рк)а3р;/р°2
21г
Ак
(36)
Начальные условия.
В начальный момент времени тела донора и акцептора имеют заданную температуру и влажность.
Т\(0,х1,х2,х3) = Т2(0,х1,х2,х3) = Т3 (0,х1,х2,х3) = Т0(л^,х2,х3)
По распределению влажности в пористой системе конкретного материала,
а
зная плотность его твердой фазы р3 и пористость 1°, можно вычислить объемное содержание каждой из трех фаз, в том числе, определить поле
а,
концентраций 20.
а
20
а^п,хе УА
"20' а^0,хеГВ
(37)
а2(°,*1,*2,*з) =а20 (*Ъ*2>*з) •
В начальный момент пар и жидкость неподвижны,
?1(0,*ь*2>*з) = 0 ^ ^2(03х1,х2’Хз) = 0 (38)
Объемная концентрация пара в образцах в начальный момент времени /0 определяется с использованием известных термодинамических соотношений по равновесному значению относительной влажности воздуха и заданной температуре Т0
*0 =1
Х(®>х1>х2’хз) = %()(х1’х2’хз) .
Хо-^1'А уг>хе Vе*
Х° ’ . (39)
Давление парогазовой смеси при /=0 вычисляется с использованием уравнений состояния по найденному распределению /о-
Начальные распределения по сечению объемных содержаний свободной воды т2св, воды в смачивающих пленках т2см и тонких слоях т2тс зависят от температуры и влажности среды. Для их вычислений пока нет точных методик, поэтому придется принять дополнительные гипотезы и расчетные схемы, позволяющие определить концентрацию различных типов вода в материале. В [11] предложена расчетная схема, которая связывает данные структурные параметры со значением влажности материала на примере материала древесины.
Предварительные оценки. Направление капиллярных потоков зависит от соотношения значений давления жидкой фазы на поверхности контакта двух
О А
материалов. Из (27) следует, что если на границе контакта Э) (Р >(Р , то материал донора (/)) обладает более высоким, чем материал акцептора (А) уровнем расклинивающего давления. Если граница контакта капиллярно-пористых тел проходит, в том числе по жидкой фазе, то имеет место процесс обезвоживания донора и увлажнение акцептора за счет течения смачивающих пленок.
Рисунок 3 - Изотермы сорбции различных материалов при 71=293 К.
Как видно из рисунка 3, в равновесном состоянии при влажности материала, меньшей значения максимальной гигроскопической влажности, для обезвоживания глиняного сырца при температуре 71=293 К за счет жидкофазного переноса в качестве акцептора (посредника) может выступать цементный камень (В/Ц=0,4). Данный механизм не работает, если в качестве посредника выступает, например, теплопор, корундовый или пеношамотный легковесы и другие (для указанных условий).
Интенсивность влагопереноса зависит от градиента расклинивающего давления в граничных слоях донора и акцептора. Оценку этой величины можно дать, дифференцируя уравнение (27), в направлении, перпендикулярном плоскости контакта образцов. Для определенности обозначим это направление х. Для однородного температурного поля в приближении однотемпературной модели
ф2см = КТ 1 дф ди дх V (р ди дх ^40^
Из (40) видно, что интенсивность жидкофазного переноса зависит от градиента влагосодержания в материале, а также его сорбционных свойств -
. 1 д ф
_ д
комплекса ср ои и температуры. Скорость стационарного
фильтрационного течения жидкости можно получить из уравнения движения (17) и (40):
(41)
Формула (41) позволяет оценить кинетику влагопереноса в предельных ситуациях. С использованием опытных изотерм сорбции и известной эмпирической формулы Аверьянова [14] получено, что в изотермических
условиях для глиняного сырца функция является убывающей, а функция
относительной фазовой проницаемости ^22см) _ возрастающей. Таким образом, с ростом насыщенности капиллярно-пористой системы, интенсивность переноса связанной влаги будет зависеть от вклада двух конкурирующих факторов, определяемых, прежде всего особенностями структуры материала.
Влияние температуры на кинетику влагопереноса, как видно из (41), определяется явной линейной зависимостью от Т, а также вкладом величины вязкости воды ц2(1) - убывающей функции температуры [15] и значения упругости водяного пара ср(Т) - возрастающей функции Т при данном влагосодержании [16]. Расчеты по формуле (41) показывают, что с увеличением температуры происходит интенсификация капиллярных перетоков.
Перенос влаги в виде пара происходит наряду с жидкофазным переносом, а также без последнего при потере жидкофазного контакта между образцом и посредником. Заметим, что величина давления парогазовой фазы зависит от теплофизических и структурных факторов, в том числе, от удельной поверхности пор, объемного содержания фаз, коэффициента проницаемости. Направление переноса пара, а также интенсивность переноса могут быть определены путем исследования математической модели (1-39).
Заключение. Процесс адсорбционно-контактной сушки является многофакторным и нестационарным. Для изучения развития полей теплофизических переменных в пористой системе донора и акцептора, работы отдельных механизмов влагопереноса, а также выбора рациональных режимов АКС на основе методологии механики многофазных систем синтезирована математическая модель (1-39). Основным методом анализа данной нелинейной модели могут быть численные методы. Анализ сопряженной математической модели тепломассопереноса в составном теле донор-акцептор позволяет дать оценку физических и структурных факторов, определяющих характер процессов тепло- и массопереноса в рассматриваемом объекте при различных
температурно-влажностных воздействиях.
Литература
1. Современные подходы к исследованию и описанию процессов сушки пористых тел / Под ред. В.Н Пармона. - Новосибирск : Изд-во СО РАН, 2001. - 300 с.
2. Wei Wang, Guohua Chen, Arun S. Mujumdar Physical Interpretation of Solids Drying: An Overview on Mathematical Modeling Research // Drying Technology. - 2007. - V.25, №4. P. 659-668.
3. Шмитько Е.И., Суслов А.А, Усачев A.M. Новый способ повышения эффективности процессов сушки керамических изделий // Строительные материалы. - 2006. - №5. - С. 20-22.
4. Усачев А.М. Повышение эффективности процесса сушки сырца керамического кирпича пластического формования : Автореф. дис... канд. техн. наук. - Воронеж, 2006. - 24 с.
5. Лыков А.В. Теория сушки. - М. : Энергия, 1968. - 471 с.
6. Фенелонов В.Б. Введение в основы адсорбции и текстурологии.
URL: http://www.catalysis.ru/block/index.php9ID=5&SECTION_ID=207 . Дата обращения:
01.04.2013
7. Боресков Г.К. Гетерогенный катализ. М. : Наука, 1986. - 303 с.
8. Прозоров Е.И., Ставров В.М., Чураев Н.В. Кинетика удаления жидкости из капиллярно-пористых тел, псевдоожиженных пористыми частицами // Инженерно-физический журнал. - 1978. - Том 34. - С. 423-429.
9. Языков Н.А., Симонов А.Д., Фенелонов В.Б. Механизм массопереноса в процессе адсорбционно-контактной сушки материалов // Теоретические основы химической технологии. -1997. - Том 31, №4. - С. 409 - 415.
10. Нигматулин Р.И. Основы механики гетерогенных сред. - М. : Наука, 1978. - 336 с.
11. Дорняк О.Р. Тепломассоперенос в ненасыщенных коллоидных капиллярно-пористых анизотропных материалах: Автореф. дис... д-ра техн. наук. - Воронеж, 2007. - 32 с.
12. Дорняк О.Р. Аль Сарраджи С. X. М. Математическая модель температурно-влажностного состояния парогазовой смеси в негерметичной камере с границами различной степени проницаемости // Лесотехнический журнал. - 2012. - №4. - С. 7 - 14.
13. Гринчик Н.Н. Процессы переноса в пористых средах, электролитах и мембранах. -Минск: АНК «Институт тепло-и массообмена» АН Беларуси, 1991. - 251 с.
14. Полубаринова-Кочина П Я. Теория движения грунтовых вод [Текст] / П. Я. Полубаринова-Кочина. - М. : Наука, 1977. - 664 с.
15. Краткий физико-технический справочник: Справочник / Под редакцией К.П. Яковлева. Том 1. Математика. Физика. - М. : Изд-во физ-матем. лит., 1960. - 446с.
16. Никитина Л.М, Термодинамические параметры и коэффициенты массопереноса во влажных материалах. - М. : Энергия, 1968 - 500 с.
References
1. Sovremennye podhody k issledovaniju i opisaniju processov sushki poristyh tel / Pod red. V.N. Parmona. - Novosibirsk : Izd-vo SO RAN, 2001. - 300 s.
2. Wei Wang, Guohua Chen, Arun S. Mujumdar Physical Interpretation of Solids Drying: An Overview on Mathematical Modeling Research // Drying Technology. - 2007. - V.25, №4. P. 659-668.
3. Shmit'ko E.I., Suslov A.A., Usachev A.M. Novyj sposob povyshenija jeffektivnosti processov sushki keramicheskih izdelij // Stroitel'nye materialy. - 2006. - №5. - S. 20-22.
4. Usachev A.M. Povyshenie jeffektivnosti processa sushki syrca keramicheskogo kirpicha plasticheskogo formovanija : Avtoref. dis... kand. tehn. nauk. - Voronezh, 2006. - 24 s.
5. Lykov A.V. Teorija sushki. - M. : Jenergija, 1968. - 471 s.
6. Fenelonov V.B. Vvedenie v osnovy adsorbcii i teksturologii. URL:
http://www.catalysis.ru/block/index.php?ID=5&SECTION_ID=207 . Data obrashhenija: 01.04.2013
7. Boreskov G.K. Geterogennyj kataliz. M. : Nauka, 1986. - 303 s.
8. Prozorov E.I., Stavrov V.M., Churaev N.V. Kinetika udalenija zhidkosti iz kapilljarno-poristyh tel, psevdoozhizhennyh poristymi chasticami // Inzhenemo-fizicheskij zhurnal. -1978. - Tom 34. - S. 423-429.
9. Jazykov N.A., Simonov A.D., Fenelonov V.B. Mehanizm massoperenosa v processe adsorbcionno-kontaktnoj sushki materialov // Teoreticheskie osnovy himicheskoj tehnologii. - 1997. -Tom 31, №4. - S. 409-415.
10. Nigmatulin R.I. Osnovy mehaniki geterogennyh sred. - M. : Nauka, 1978. - 336 s.
11. Dornjak O.R Teplomassoperenos v nenasyshhennyh kolloidnyh kapilljarno-poristyh anizotropnyh materialah: Avtoref. dis... d-ra tehn. nauk. - Voronezh, 2007. - 32 s.
12. Dornjak O.R. АГ Sarradzhi S. H. M. Matematicheskaja model' temperaturno-vlazhnostnogo sostojanija parogazovoj smesi v negermetichnoj kamere s granicami razlichnoj stepeni pronicaemosti // Lesotehnicheskij zhurnal. - 2012. - №4. - S. 7 - 14.
13. Grinchik N.N. Processy perenosa v poristyh sredah, jelektrolitah i membranah. - Minsk: ANK «Institut teplo-i massoobmena» AN Belarusi, 1991. - 251 s.
14. Polubarinova-Kochina P.Ja. Teorija dvizhenija gruntovyh vod [Tekst] / P. Ja. Polubarinova-Kochina. -M. : Nauka, 1977. - 664 s.
15. Kratkij fiziko-tehnicheskij spravochnik: Spravochnik / Pod redakciej K.P. Jakovleva. Tom 1. Matematika. Fizika. - M. : Izd-vo fiz-matem. lit., 1960. - 446s.
16. Nikitina L.M., Termodinamicheskie parametry i kojefficienty massoperenosa vo vlazhnyh materialah. - M. : Jenergija, 1968 - 500 s.
