Математическое моделирование тепло- и массопереноса при адсорбционно-контактной сушке керамических материалов (вычислительный эксперимент)

Размещено на http://www.allbest.ru

1

Научный журнал КубГАУ, №89(05), 2013 года

УДК 532.546+536.24: 666.3.047

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ТЕПЛО- И МАССОПЕРЕНОСА ПРИ АДСОРБЦИОННО-КОНТАКТНОЙ СУШКЕ КЕРАМИЧЕСКИХ МАТЕРИАЛОВ (ВЫЧИСЛИТЕЛЬНЫЙ ЭКСПЕРИМЕНТ)

Дорняк Ольга Роальдовна

д.т.н., доцент

Аль Сарраджи Салах Хусейн Мела

аспирант

Усачев Александр Михайлович

к.т.н., доцент

Шмитько Евгений Иванович

д.т.н., профессор

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

One-temperature limited version of the general model of adsorption-contact drying for the ceramic brick is considered. Numerical modeling of heat and mass transfer in the donor and the acceptor samples for one-dimensional approach is carried out. It is shown, in particular, that the moisture diffusion coefficient of material in the contact zone can vary significantly due to changes in the conditions for liquid-phase transfer

Ключевые слова: МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ, КЕРАМИЧЕСКИЕ МАТЕРИАЛЫ, АДСОРБЦИОННО-КОНТАКТНАЯ СУШКА

Keywords: MATHEMATICAL MODELING, CERAMIC MATERIALS, ADSORPTION-CONTACT DRYING

Введение

Для сушки кирпичей, изготовленных из высокочувствительного к термообработке сырья, разработана оригинальная технология применения адсорбционно-контактной сушки (АКС), уменьшающая вероятность растрескивания образцов [1-2]. Сущность такого метода заключается в следующем. Отформованный сырец помещается на подложку из пористого материала с высокой влагоемкостью и накрывается сверху влагонепроницаемым колпаком. Подложка с тыльной стороны подвергается высушиванию теплоносителем. Тем самым осуществляется диффузия влаги из сырца в подложку, а уже за счет подачи и отбора сушильного агента происходит интенсивная конвективная сушка самой подложки [1]. Адсорбционно-контактный метод позволяет снизить нежелательную тепловую нагрузку на материал и одновременно усилить эффективность сушки [3-4]. Поиск эффективных режимов АКС требует решения задач тепло- и массопереноса в капиллярно-пористых телах с учетом особенностей технологического процесса и структурных свойств материалов. Задачи данного типа исследуются методами математического моделирования, при этом известны различные подходы созданию моделей. Метод математического моделирования, наиболее известный среди инженеров, связан с именем А.В. Лыкова [5]. Широко используется теория многофазной фильтрации [3,7]. Работы по механике многофазных систем Р.И. Нигматулина [6] позволили создать математические модели сушки дисперсных систем с учетом специфики механизмов переноса и взаимодействия отдельных фаз. Ряд современных направлений в теории сушки основан на идеях иерархического моделирования [4].

В первой части публикуемой работы синтезирована математическая модель тепломассопереноса при АКС керамического кирпича на основе методологии многофазных систем [6]. В данной части исследован предельный, более простой, вариант этой модели, который близок по форме к диффузионно-фильтрационной модели А.В.Лыкова [5] и получен при следующих предположениях [8-9]. Температура всех фаз дисперсной системы одинакова (однотемпературное приближение). Фазовые переходы происходят по равновесной схеме. Вкладом работы внутренних сил по расширению парогазовой фазы и влиянием конвективного переноса тепла в этой фазе можно пренебречь.

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

Рассмотрим часть плиты посредника из тонкопористого материала, на которой уложена заготовка глиняного кирпича. Сверху кирпич прикрыт влагозащитным колпаком. Система находится в конвективной сушильной камере при атмосферном давлении с заданными характеристиками сушильного агента - температурой T_c и относительной влажностью ц_c. Ставится задача изучения особенностей развития полей влажности, температуры, давления в составном образце, содержащем кирпич (донор) и посредник (акцептор). Предположим, что поля влагосодержания, температуры, общего давления и концентрации пара однородны в направлении осей у и z, так что u=u(t,x), T=T(t,x), p=p(t,x), ч=ч(t,x) (рисунок 1). Тогда для изучения процессов переноса достаточно рассмотреть задачу в одномерном приближении.

Рисунок 1 - Расчетная схема адсорбционно-контактной сушки глиняного сырца.

Математическая модель

Уравнение влагопроводности в составном образце запишем в виде

,

; ; . (1)

где коэффициент диффузии a_m зависит от изменяющихся значений влагосодержания, температуры и проницаемости [8-9]:

; . (2)

В данной работе влияние на фильтрационный перенос насыщенности соответствующей фазой не рассматривается, т.е. функции относительной проницаемости и . Как видно из (2), для описания кинетики влагопереноса необходимы уравнения политерм сорбции и для донора и для акцептора -

; . (3)

Уравнение теплопроводности представляется в следующей форме

; (4)

, ;

;;; ; i=1,2,3.

Объемная концентрация отдельных фаз вычисляется по формулам

; . (5)

Коэффициент фазовых переходов зависит от влажности, температуры и давления пара [8-9]:

. (6)

Уравнение давления парогазовой смеси имеет вид

; (7)

.

Для вычисления полей концентраций паровой и газовой компоненты в паровоздушной фазе -

и (8)

служит следующее уравнение:

; (9)

; .

Уравнения состояния идеального газа для смеси и ее компонент имеют вид

; ; .

; . (10)

Интенсивность фазовых переходов, определяемая в многофазной модели с помощью соотношения Герца-Кнудсена-Ленгмюра, после предельного перехода к однотемпературной модели записывается в виде

. (11)

Формула (11) уже не имеет первоначального физического смысла, связанного с оценкой разности молекулярных потоков от поверхности жидкости в пар и обратно [10], и к₁ - не является коэффициентом аккомодации. Для вычисления j, таким образом, имеем аналог формулы Дальтона (уравнения удельной поверхностной влагоотдачи) [4,11], согласно которой массовая скорость испарения жидкости с поверхности осушаемого образца в установившемся процессе

, ,

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

Для определения давления насыщенного пара при данной температуре воспользуемся известной формулой Антуана, которая аппроксимирует кривую фазовых переходов 1 рода в диапазоне изменения температуры 284 - 441 К следующим образом [12]:

. (12)

Здесь величина p_sat имеет размерность мм. рт.ст. Константы в (12) имеют следующие значения; ; [12].

Начальные условия задают произвольное распределение зависимых переменных:

; . (13)

; ; ; .

Отметим, что функция не является независимой функцией. Концентрация пара в образце однозначно определяется величиной относительной влажности воздуха, которая в равновесных условиях связана с влажностью материала уравнениями политермы сорбции (3).

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

. (14)

Для температуры на внешних границах составного образца ставятся условия Ньютона-Рихмана

; . (15)

Граничные значение давления парогазовой смеси и концентрации каждой из компонент в капиллярно-пористом образце полагаем равными соответствующим значениям в соседнем объеме - в сушильной камере для акцептора (подложки) и в пространстве под колпаком для донора (высушиваемого образца) -

; ; ; . (16)

Параметры сушильного агента под влагоизолирующим колпаком T_k, p_k, ч_k изменяются в процессе сушки. Они определяются из решения сопряженной задачи тепломассопереноса в образце и в среде под колпаком. Поведение парогазовой фазы, заполняющей негерметичный объем сложной формы, ограниченный поверхностями материалов с различной степенью проницаемости - весьма сложная и существенно неодномерная задача. Методика ее приближенного решения построена в [13], следуя [14].

Уравнения сохранения массы записаны, в балансовой форме для воды, паровой и газовой компонент:

; (17)

; (18)

; ; . (19)

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

, ;

, .

Для вычисления расхода парогазовой смеси через возможные отверстия между плитой акцептора и колпаком используется формула, основанная на теореме Бернулли:

, (20)

где б_m - коэффициентом расхода сужающего устройства, е_m - поправочный множитель, учитывающий уменьшение плотности вещества при его прохождении через сужающее устройство, S_c - площадь отверстия [15]. Будем считать, что расход паровой и газовой компонент в отдельности из объема под в колпаком в сушильную камеру пропорционален их концентрации:

, ;,. (21)

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

Балансовое уравнение для температуры в паровоздушном слое между образцом и защитной поверхностью имеет вид:

(22)

, , ,

, , , ,

, .

Начальные условия для уравнений (17-19, 22) приняты в предположении, что температурно-влажностное состояние газообразной фазы под влагоизолирующим колпаком при t=0 совпадает с характеристиками сушильного агента в камере:

, ; . (23)

По заданному значению ц_с определяем начальную концентрацию паровой и газовой компонент (ч_с и (1-ч_с) соответственно) с помощью известных соотношений [16]:

, ,

, . (24)

Величина начальной массы паровоздушной смеси под влагоизолирующим колпаком -

; ; . (25)

Масса конденсированной фазы в начальный момент времени

. (26)

Система нелинейных дифференциальных уравнений (1, 4, 7, 9, 17-19, 22), дополненная уравнениями состояния (10-12), начальными (13, 23-26) и граничными условиями (14-16) описывает нестационарный тепло- и массоперенос при адсорбционно-контактной сушке образцов из двух различных материалов в одномерном приближении.

Безразмерная форма начально-краевой задачи

Приведем сформулированную задачу тепломассообмена к безразмерному виду. В качестве единиц величин с независимой размерностью используем характерный размер заготовки l_xap, характерное время процесса t_xap, давление p_xap, температуру T_xap. Безразмерные переменные определим следующим образом:

; ; ;

;

; ;

; ;;

;;;;

; ;;;

; .

Безразмерные уравнения переноса в материале донора и акцептора имеют вид:

; ;

; ; ;

; ; .

;

, ;

; ; ;

; i=1,2,3. ;

.

; ; .

; ,

; ; ; .

;

; ; .

; ; .

;

; ;

; ; ;. (27)

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

; ;

; ;

, ; ,

; ; .

; ;

;

; ;

, , .

, , , . (28)

Вычислительный эксперимент

Система уравнений (27-28) исследована численно. Дискретный аналог уравнений в частных производных построен методом контрольного объема [17]. Обыкновенные дифференциальные уравнения решаются с использованием метода Рунге-Кутта. Вычислительный алгоритм является итерационным в связи с нелинейностью поставленной задачи. Тестирование программного комплекса проведено на основе ряда точных решений уравнений с дополнительными источниковыми членами. Сравнивая расчетные зависимости среднего влагосодержания сырца и данные опытов (рисунок 4), можно сделать вывод об их удовлетворительном совпадении.

На рисунках 2-7 представлены результаты численного анализа, иллюстрирующие особенности процессов тепло- и массопереноса при адсорбционно-контактной сушке образцов, полученных из высокочувствительной глины способом пластического формования. В качестве посредника рассматривается цементный камень (В/Ц=0,4).

Сорбционные свойства материалов описаны с использованием уравнения Л.Б. Цимерманиса, полученного на основе термодинамического анализа [18]. Для типичных капиллярно-пористых тел это уравнение имеет вид

,

где u - равновесное влагосодержание, ц - относительная влажность воздуха, а₀ и k - безразмерные постоянные, называемые соответственно структурной активностью капиллярно-пористого материала и коэффициентом приращения активности связи. Как показано в [18], равновесное влагосодержание цементного камня снижается с увеличением температуры, и для максимального гигроскопического влагосодержания справедливо соотношение

.

Здесь u_S0 - максимальное гигроскопическое содержание при T=273 K, б_ST - постоянный множитель.

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

Водяной пар и воздух: R=8.3144 Дж/(моль^.К); B_1v=461.9 Дж/(кг^.K); B_1g=284 Дж/(кг^.K); ₁=10^-6 Па^.с; ₁=0.0248 Вт/(м^.K); с_{p 1v}=2.034^.10³ Дж/(кг^.K); с_{p 1g}=9.05^.10² Дж/(кг^.K); с_{v 1v}=1.58^.10³ Дж/(кг^.K); с_{v 1g}=7.07^.10² Дж/(кг^.K); p_atm=1.013^.10⁵; ₁=10^_6 Па^.с; D=2.16^.10^_5 м²/c.

Вода: ₂=10³ кг/м³; ₂=10^-3 Па^.с; V'₂=18^.10^_6 м³/моль; ₂=0.648 Вт/(м^.K), с₂=4.2^.10³ Дж/(кг^.K); L=2.25^.10⁶ м³/моль; ₁=1.875^.10^_11; =0.04; ₂₃=10^_17-10^_18м²; ₁₃=10^_11-10^_12 м².

Донор: ₃=2370 кг/м³; =1400 кг/м³; ₃=0.73 Вт/(м^.град), с₃=0.796^.10³ Дж/(кг^.K); u_S0=0.23 кг/кг; б_ST=1.1^.10^_3 K^-1; a₀=0.8862; k=3.12; ₁₃=2.5^.10^_14 м²; ₂₃=10^_18 м²; s₁₂=7.98^.10⁶ м^-1; Н_d=0.015 м; a_а=0.05 м, b_а=0.02 м.

Акцептор: ₃=4800 кг/м³; =1680 кг/м³; ₃=0.88 Вт/(м^.град), с₃=0.838^.10³ Дж/(кг^.K); u_S0=0.0967 кг/кг; б_ST=0.418^.10^_3 K^-1; a₀=0.6640; k=14.8; ₁₃=5^.10^_14 м²; ₂₃=10^_20 м²; s₁₂=0.15^.10⁶ м^-1; Н_а=0.01 м.

Колпак: h_k=0.065 м. Окружающая среда: T_c=293ч393 К; p_с=p_atm; _с=45ч100%.

Начальные значения: u₀^A=0.0169 кг/кг; u₀^D=0.28 кг/кг; T₀=293 K; p₁₀=p_atm; M₂₀=0; Характерные значения: t_xap=1281600 с?356 час; l_xap= Н_d.

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

Результаты и обсуждение и вычислительного эксперимента

а)	б)
Рисунок - 2. Изменение коэффициента влагопроводности в различных точках образцa акцептора - а) и донора - б) со временем при цс= 0,45, Tc=323 K, vc=0.37 м/с: x/lxap=0.0066 - 1a, 0.3317 - 2a, 0.6633 - 3a; x/lxap=1.6567- 1d, 1.1692 - 2d, 0.6716 - 3d; пунктирная линия соответствует коэффициенту влагопроводности, усредненному по объему образца. Величина am0=1.4920.10-7 м2/c.

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

Расчеты коэффициента влагопроводности в характерном режиме адсорбционно-контактной сушки показывают существенную разницу значений a_m в различных точках образца и посредника (рисунок 2). Кривые 3a и 3d относятся к точкам в области, принадлежащей акцептору и донору вблизи зоны контакта, 2a и 2d - в центре, и 1a и 1d - рядом с внешними границами этих тел. Особенности изменения a_m связаны с динамикой развития профиля влагосодержания (рисунок 3).

Пористая структура донора практически весь расчетный период заполнена свободной водой, которая в отсутствие массовых сил неподвижна. Обезвоживание донора происходит за счет испарения в объеме и жидкофазного переноса на границе с акцептором за счет капиллярных эффектов 1 и 2 рода. При этом жидкофазный переход в области соприкосновения плиты посредника и кирпича может происходить при наличии «мостика» жидкости - жидкого кластера [4]. В материале акцептора работают те же механизмы обезвоживания что и в доноре, а при уменьшении влажности ниже предельного гигроскопического значения добавляется механизм влагопереноса за счет градиента расклинивающего давления в смачивающих пленках. Кинетика сушки в доноре и акцепторе различна по интенсивности, так что можно ожидать, что жидкофазный контакт между ними имеет динамический характер. Изменение влагосодержания в акцепторе в зоне контакта может привести к полному отсутствию контакта с донором по жидкому кластеру или резкому снижению его интенсивности. Именно такое положение иллюстрируют немонотонные зависимости коэффициента влагопроводности от времени - кривые 2a, 3a на рисунке 2а. Замедлившийся жидкофазный перенос через границу тел может активизироваться вновь, когда концентрация жидкой фазы вырастет за счет капиллярных течений на границе, при этом коэффициент влагопроводности вновь возрастает.

а)	б)
в)	г)
д)	е)
Рисунок 3. Распределение влагосодержания по сечению комплекса донор-акцептор при цс= 0,45, Tc=323 K, vc=0.37 м/с для t/txap=0.005 - а); 0.025 - б); 0.05 - в); 0.15 - г); 0.20 - д) 0.25 - е).

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



Рисунок - 4. Изменение среднего влагосодержания uav (кг/кг) в доноре (d) и акцепторе (a) со временем при Tc=323 K, vc=0,37 м/с для цс= 0,45 - 1; 0,60 - 2; 0,80 - 3. Маркером + отмечены данные эксперимента [2].	Рисунок - 5. Изменение среднего влагосодержания uav (кг/кг) в доноре (d) и акцепторе (a) со временем при цс= 0,45, vc=0,37 м/с для Tc=353 K - 1; 323 - 2; 293 - 3.

На рисунке 4 приведены кривые кинетики адсорбционно-контактной сушки при различных значениях относительной влажности сушильного агента. Как видно из графиков, влияние данного фактора на скорость обезвоживания, как и при других видах сушки, весьма существенно. Если акцептор высушивается недостаточно интенсивно, то влажность донора также снижается медленно. При ц_с=0.45 акцептор приобретает равновесную влажность, соответствующую температуре окружающей среды, существенно быстрее (почти в 2 раза), чем при ц_с=0.6. Для упругости пара в паровоздушной смеси с ц_с=0.8 средняя влажность акцептора в течение расчетного периода практически не снижается. Влага, заполняющая пористую структуру донора, наиболее интенсивно покидает его при ц_с=0.45, а при ц_с=0.8 длительность процесса сушки образца существенно возрастает.

Графики на рисунке 5 иллюстрируют влияние температуры сушильного агента на кинетику АКС. Повышение температуры на 30 градусов приводит к значительной интенсификации процесса обезвоживания донора. Например, в жестком режиме (T_c=353 K, ц_с= 0,45) заданная конечная средняя влажность высушиваемого материала сырца достигается в 5 раз быстрее, чем при T_c=323 K, ц_с=0,45.

Рисунок - 6. Изменение со временем средней по объему температуры в образце донора (1) и акцептора (2) со временем при цс= 0,45, Tc=323 K, vc=0.37 м/с.	Рисунок - 7. Изменение со временем критерия Кирпичева в образце донора (1) и акцептора (2) со временем при цс= 0,45, Tc=323 K, vc=0.37 м/с.

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



а)	б)
Рисунок - 8. Изменение со временем среднего по объему давления парогазовой смеси (а) и концентрации пара в парогазовой смеси в образце (б) для донора (1) и акцептора (2) со временем при цс= 0,45, Tc=323 K, vc=0.37 м/с.

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

Концентрация паровой компоненты снижается в течение процесса сушки, как в высушиваемом образце, так и в посреднике. Скорость снижения определяется диффузионным числом Пекле (Pe_D=8,1^.10^_6) и числом Дарси. На внешних границах донора и акцептора концентрация пара не одинакова, т. к. в среде сушильного агента, содержание пара существенно ниже, чем под колпаком, где воздух постепенно приближается к состоянию насыщения (при данной температуре). Выравнивание концентрации пара в объеме V^A происходит быстрее, чем в V^D, поскольку при заданных параметрах условия сопротивление переносу газообразной фазы в доноре выше, чем в акцепторе. Для донора число Дарси Da₁₃ имеет в несколько раз меньшее значение, чем для акцептора.

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

Рисунок 9 иллюстрирует развитые профили теплофизических переменных, соответствующих окончанию процесса АКС. Граница между донором и акцептором имеет безразмерную координату x/l_xap=0.66. Здесь имеют место скачки влагосодержания u, концентрации пара ч, избыточного давления парогазовой смеси, интенсивности парообразования j₁₂ и коэффициента диффузии влаги a_m. Величина разрыва связана с текстурными особенностями контактирующих материалов, в частности, объемными концентрациями газообразной, жидкой и твердой фаз в каждом материале - б₁, б₂, б₃. (рисунок 9), а также удельными поверхностями пористой системы этих материалов и .



Рисунок - 9. Распределение по сечению комплекса донор-акцептор безразмерных теплофизических параметров при ц_с= 0,45, T_c=323 K, v_c=0.37 м/с для t/t_xap=0.25: влагосодержания u; температуры T; концентрации пара ч; избыточного давления газовой фазы и пара - p₁-p_atm и p_1v-p_atm; плотности парогазовой среды с₁; интенсивности парообразования j₁₂; относительной влажности ц; коэффициента влагопроводности a_m; объемного содержания газообразной, жидкой и твердой фазы - б₁, б₂, б₃.

Основной проблемой, которую решает использование АКС при высушивании сырца из высокочувствительной глины, является проблема снижения трещинообразования. Экспериментальное исследование [2] позволило авторам обосновать критерий прочности образца из глиняного сырца при сушке. Для прогноза трещиностойкости выбрано число Кирпичева Ki [20]:

.

где q_m - плотность массового потока влаги на поверхности плиты кг/м²ч; h - половина толщины пластины, м; u₀ - начальное влагосодержание, кг/кг; a_m - коэффициент влагопроводности, м²/ч; с_m - средняя плотность сухого материала, кг/м³. Оценивая градиент влагосодержания как разность влагосодержаний между центральным и поверхностным слоем, получаем выражение вида

.

Экспериментально установлено, что при Ki<0.58 растрескивания высушиваемого глиняного сырца не происходит [2]. На рисунке 7 показано изменение числа Кирпичева, вычисленного на основе рассматриваемой математической модели для режима АКС, с параметрами сушильного агента ц_с= 0,45, T_c=323 K, v_c=0.37 м/с. Данный режим, как установлено опытами [2], обеспечивает допустимый интервал изменения критерия Кирпичева для донора. Что касается материала подложки, то он находится в условиях жесткого режима сушки и его усадочные напряжения могут превысить предельный уровень, что допустимо в исследуемом технологическом процессе. Результаты вычислительного эксперимента полностью согласуются с выводами [2].

Заключение

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

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

Список обозначений

a_m - коэффициент влагопроводности, м²/с;

a - коэффициент температуропроводности, м²/с; ;

B - индивидуальная газовая постоянная, Дж/(кг^.К);

с - теплоемкость, Дж/(кг^.К);

D - коэффициент бинарной диффузии, м²/с;

i - энтальпия, Дж/кг;

j -поток массы пара, обусловленный фазовыми переходами, отнесенный к единице времени и единице площади, кг/(м^2.с);

H= H_а+H_d - координата верхней границы донора, м;

H_а - высота акцептора, м;

H_d - высота донора, м;

_j3 - коэффициент проницаемости j-ой фазы при полном насыщении пористой системы, м²;

L - удельная теплота фазовых переходов Дж/кг;

l - длина, м;

M - масса, кг;

m - массовый расход, кг/с;

n - внешняя нормаль к поверхности;

p - давление, Па;

R - универсальная газовая постоянная, Дж/(моль^.К);

s_ij - удельная поверхность раздела фаз i и j, м^-1;

S - площадь поверхности, м²;

T - температура, К;

t - время, с;

q -тепловой поток, Дж/с;

V - объем, м³;

v - скорость, м/с;

- объем моля воды;

u - влагосодержание, кг/кг;

w - диффузионная скорость, м/с;

x, y, z - декартовы координаты, м; - объемное содержание фазы;

- коэффициент теплоотдачи от акцептора к сушильному агенту, Вт/(м²K);

- коэффициент теплоотдачи от донора к газообразной фазе под влагонепроницаемым колпаком, Вт/(м²K);

- коэффициент теплоотдачи от паровоздушной среды под колпаком к сушильному агенту, акцептору и донору соответственно, Вт/(м²K);

µ - вязкость, Па^.с;

- критерий фазовых переходов;

- коэффициент аккомодации;

₁ , ₂ - коэффициенты;

- концентрация пара;

-коэффициент теплопроводности, Вт/(м^.К);

- насыщенность; - плотность, кг/м³;

- относительная фазовая проницаемость;

ц - относительная влажность воздуха.

Нижние индексы относятся: 1,2,3 - к газообразной, жидкой и твердой фазам; с - сушильный агент; g -газовая компонента; k - паровоздушная среда под влагоизолирующм колпаком; v -пар; sat - cостояние насыщения; Г - принадлежность к внешним границам образца; 0 -начальное значение; хар - характерное значение параметра; atm - атмосферное.

Верхние индексы: - истинное значение физической величины; * - безразмерное значение; A - акцептор; D - донор.

Безразмерные комплексы: - тепловое число Пекле; - массообменное число Пекле; - диффузионное число Пекле; - число Кутателадзе; - число Коссовича; ; - число Рейнольдса; - число Дарси; - число Кельвина; , i= D, A, k - число Нуссельта; ; - числа Био.

Литература

1. Шмитько Е.И., Суслов А.А., Усачев А.М. Новый способ повышения эффективности процессов сушки керамических изделий // Строительные материалы - 2006. - №5. - С. 20-22.

2. Усачев А. М. Повышение эффективности процесса сушки сырца керамического кирпича пластического формования. - автореф. Дис. ... Канд. Техн. Наук : 05.23.05 / Воронеж, 2006. - 24 с.

3. Хейфец Л.И., Неймарк А.В. Многофазные процессы в пористых средах. - М. : Химия, 1982. - 320 с.

4. Современные подходы к исследованию и описанию процессов сушки пористых тел / Под ред. В. Н. Пармона. - Новосибирск : Изд-во СО РАН, 2001. - 300 с.

5. Лыков А. В. Тепломассообмен. - М. : Энергия, 1971. - 560 с.

6. Нигматулин Р.И. Динамика многофазных сред. - М.: Наука, 1987. Ч. 1. - 464 с.; Ч. 2. - 360 с.

7. Гринчик H.H. Процессы переноса в пористых средах, электролитах и мембранах. - АНК "Институт тепло- и массообмена" АН Беларуси, г. Минск, 1991г. - 251 с.

8. Дорняк О. Р. Тепломассоперенос в ненасыщенных коллоидных капиллярно-пористых анизотропных материалах. - автореф. Дис. ... Д-ра техн. Наук : 01.04.14 / Воронеж, 2007. - 32 с.

9. Дорняк О. Р. Уравнение общего давления в математической модели диффузионно-фильтрацинного влаготеплопереноса в коллоидных капиллярно-пористых телах // Вестник Воронежского государственного технического университета. 2011. Т. 7. № 12-1. С. 68-71.

10. Исаченко В. П. Теплообмен при конденсации. - М. : Энергия, 1977. - 240 с.

11. Светлов Ю. В. Термовлажностные процессы в материалах и изделиях легкой промышленности. - М. : Издательский центр академия, 2006. - 272 с.

12. Рид Р., Праусниц Дж., Шервуд. Т. Свойства газов и жидкостей. - Л. : Химия, 1982, 463 с.

13. Дорняк О. Р., Аль Сарраджи С. Х. М. Математическая модель температурно-влажностного состояния парогазовой смеси в негерметичной камере с границами различной степени проницаемости // Лесотехнический журнал. - 2012. № 4. С. 7-14.

14. Дейч М. Е. Гидродинамика двухфазных сред. - М. : Энергоатомиздат, 1981. - 472 с.

15. Кремлевский П. П. Расходомеры и счетчики количества: Справочник. -- Л. : Машиностроение. Ленингр. отд-ние, 1989. -- 701 с.

16. Крутов В. Н. Исаев С. И., Кожинов И. А. Техническая термодинамика. - М.: Высш. Шк., 1991. - 384 с.

17. Патанкар С. Численные методы решения задач теплообмена и динамики жидкости. - М. : Энергоатомиздат, 1984. - 149 с.

18. Цимерманис Л.-Х. Б. Сорбция, структурообразование, массоперенос (термодинамика влажного тела) - М. : «Алекс», 2006. - 232 с.

19. Языков, Н.А. Механизм массопереноса в процессе адсорбционно-контактной сушки материалов / Н.А. Языков, А.Д. Симонов, В.Б. Фенелонов // Теоретические основы химической технологии. 1997. - Т.31. - №4. - с. 409 - 415.

20. Колодежнов В. Н. Безразмерные комплексы и критерии подобия в гидроаэромеханике: Справочник. Воронеж: Воронежский госпедуниверситет, 2011. - 580 с.

Размещено на Allbest.ru

...