Научная статья на тему 'ИССЛЕДОВАНИЕ КИНЕТИКИ ТЕПЛОВЛАГООБМЕНА ПРИ ТЕРМИЧЕСКОЙ ОБРАБОТКЕ И СУШКЕ ТОНКИХ ВЛАЖНЫХ ТЕПЛОИЗОЛЯЦИОННЫХ МАТЕРИАЛОВ'

ИССЛЕДОВАНИЕ КИНЕТИКИ ТЕПЛОВЛАГООБМЕНА ПРИ ТЕРМИЧЕСКОЙ ОБРАБОТКЕ И СУШКЕ ТОНКИХ ВЛАЖНЫХ ТЕПЛОИЗОЛЯЦИОННЫХ МАТЕРИАЛОВ Текст научной статьи по специальности «Физика»

CC BY
23
9
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ВЛАГОСОДЕРЖАНИЕ / КОЭФФИЦИЕНТ СУШКИ / СКОРОСТЬ СУШКИ / ИНТЕНСИВНОСТЬ СУШКИ / СУШКА ТЕПЛОИЗОЛЯЦИОННЫХ МАТЕРИАЛОВ

Аннотация научной статьи по физике, автор научной работы — Ольшанский А. И., Голубев А. Н.

Проведено исследование кинетики сушки тонких теплоизоляционных плоских материалов. Представлена аппроксимация кривой скорости сушки различными методами. При определении длительности сушки зональными методами использовалось уравнение скорости сушки с коэффициентом сушки. Изложены способы обработки опытных данных зональными методами А. В. Лыкова, В. В. Красникова и методом Б. С. Сажина. Уравнение, полученное обработкой опытных данных по методу Б. С. Сажина для определения длительности процесса сушки, содержит только величину влагосодержания прогрева материала, знать критическое влагосодержание нет необходимости. Приведены зависимости для вычисления коэффициента сушки зональными методами. По результатам эксперимента представлены формулы для установления длительности тепловой обработки теплоизоляционных материалов. Дана зависимость относительной скорости сушки от относительного влагосодержания. Обработка экспериментальных данных обобщенными комплексными переменными создает условия инвариантности, позволяющие переходить от одних переменных к другим, от одной системы координат к другой, что сокращает число экспериментов. Приведены формулы с комплексными переменными для определения длительности сушки материалов. На основании уравнения кинетики сушки А. В. Лыкова и уравнения Г. К. Филоненко для относительной скорости сушки представлены формулы для нахождения плотности тепловых потоков, интенсивности испарения влаги и температуры для периода падающей скорости сушки. Проанализированы решения экспериментальных уравнений. Установлена погрешность, вызываемая обработкой опытных данных. Представлена проверка достоверности полученных эмпирических уравнений и проведено сопоставление расчетных значений основных параметров кинетики сушки с экспериментом.

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

Похожие темы научных работ по физике , автор научной работы — Ольшанский А. И., Голубев А. Н.

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

INVESTIGATION OF THE KINETICS OF HEAT AND MOISTURE EXCHANGE DURING HEAT TREATMENT AND DRYING OF THIN WET THERMAL INSULATION MATERIALS

. The kinetics of drying of thin heat-insulating flat materials has been studied. The approximation of the drying rate curve by various methods is presented. When determining the duration of drying by zonal methods, the equation of the drying rate with the drying coefficient was used. The techniques of processing experimental data by zonal methods of A.V. Lykov, V. V. Krasnikov and by the method of B. S. Sazhin are described. The equation obtained by processing experimental data using the B. S. Sazhin method to determine the duration of the drying process contains only the amount of moisture content of the heating of the material, there is no need to know the critical moisture content. Dependences for calculating the drying coefficient by zonal methods are given. According to the results of the experiment, formulas for determining the duration of heat treatment of thermal insulation materials are presented. The dependence of the relative drying rate on the relative moisture content is presented as well. The processing of experimental data by generalized complex variables creates invariance conditions that allow moving from one variable to another, from one coordinate system to another, which reduces the number of experiments. Formulas with complex variables for determining the duration of drying of materials are given. Based on the drying kinetics equation by A. V. Lykov and the equation for the relative drying rate by G. K. Filonenko, equations are given for determining the density of heat fluxes, the intensity of moisture evaporation and temperature for the period of decreasing drying rate. Solutions of experimental equations are analyzed. The error caused by the processing of experimental data is established. Verification of the reliability of the empirical equations obtained is presented and comparison of the calculated values of the main parameters of the drying kinetics with the experiment is carried out.

Текст научной работы на тему «ИССЛЕДОВАНИЕ КИНЕТИКИ ТЕПЛОВЛАГООБМЕНА ПРИ ТЕРМИЧЕСКОЙ ОБРАБОТКЕ И СУШКЕ ТОНКИХ ВЛАЖНЫХ ТЕПЛОИЗОЛЯЦИОННЫХ МАТЕРИАЛОВ»

https://doi.org/10.21122/1029-7448-2023-66-1-66-79 УДК 66.047.7

Исследование кинетики тепловлагообмена

при термической обработке и сушке

тонких влажных теплоизоляционных материалов

А. И. Ольшанский1*, А. Н. Голубев1*

'-Витебский государственный технологический университет (Витебск, Республика Беларусь)

© Белорусский национальный технический университет, 2023 Belarusian National Technical University, 2023

Реферат. Проведено исследование кинетики сушки тонких теплоизоляционных плоских материалов. Представлена аппроксимация кривой скорости сушки различными методами. При определении длительности сушки зональными методами использовалось уравнение скорости сушки с коэффициентом сушки. Изложены способы обработки опытных данных зональными методами А. В. Лыкова, В. В. Красникова и методом Б. С. Сажина. Уравнение, полученное обработкой опытных данных по методу Б. С. Сажина для определения длительности процесса сушки, содержит только величину влагосодержания прогрева материала, знать критическое влагосодержание нет необходимости. Приведены зависимости для вычисления коэффициента сушки зональными методами. По результатам эксперимента представлены формулы для установления длительности тепловой обработки теплоизоляционных материалов. Дана зависимость относительной скорости сушки от относительного влагосодержания. Обработка экспериментальных данных обобщенными комплексными переменными создает условия инвариантности, позволяющие переходить от одних переменных к другим, от одной системы координат к другой, что сокращает число экспериментов. Приведены формулы с комплексными переменными для определения длительности сушки материалов. На основании уравнения кинетики сушки А. В. Лыкова и уравнения Г. К. Филоненко для относительной скорости сушки представлены формулы для нахождения плотности тепловых потоков, интенсивности испарения влаги и температуры для периода падающей скорости сушки. Проанализированы решения экспериментальных уравнений. Установлена погрешность, вызываемая обработкой опытных данных. Представлена проверка достоверности полученных эмпирических уравнений и проведено сопоставление расчетных значений основных параметров кинетики сушки с экспериментом.

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

Для цитирования: Ольшанский, А. И. Исследование кинетики тепловлагообмена при термической обработке и сушке тонких влажных теплоизоляционных материалов / А. И. Ольшанский, А. Н. Голубев // Энергетика. Изв. высш. учеб. заведений и энерг. объединений СНГ. 2023. Т. 66, № 1. С. 66-79. https://doi.org/10.21122/1029-7448-2023-66-1-66-79

Адрес для переписки Address for correspondence

Ольшанский Анатолий Иосифович Ol'shanskii Anatolii I.

Витебский государственный Vitebsk State

технологический университет Technological University

просп. Московский, 72, 72, Moskovsky Ave.,

220035, г. Витебск, Республика Беларусь 220035, Vitebsk, Republic of Belarus

Тел.: +375 212 47-50-26 Tel.: +375 212 47-50-26

[email protected] [email protected]

Investigation of the Kinetics of Heat and Moisture Exchange during Heat Treatment and Drying of Thin Wet Thermal Insulation Materials

A. I. Ol'shanskii4, A. N. Golubev4

'Vitebsk State Technological University (Vitebsk, Republic of Belarus)

Abstract. The kinetics of drying of thin heat-insulating flat materials has been studied. The approximation of the drying rate curve by various methods is presented. When determining the duration of drying by zonal methods, the equation of the drying rate with the drying coefficient was used. The techniques of processing experimental data by zonal methods of A.V. Lykov, V. V. Krasnikov and by the method of B. S. Sazhin are described. The equation obtained by processing experimental data using the B. S. Sazhin method to determine the duration of the drying process contains only the amount of moisture content of the heating of the material, there is no need to know the critical moisture content. Dependences for calculating the drying coefficient by zonal methods are given. According to the results of the experiment, formulas for determining the duration of heat treatment of thermal insulation materials are presented. The dependence of the relative drying rate on the relative moisture content is presented as well. The processing of experimental data by generalized complex variables creates invariance conditions that allow moving from one variable to another, from one coordinate system to another, which reduces the number of experiments. Formulas with complex variables for determining the duration of drying of materials are given. Based on the drying kinetics equation by A. V. Lykov and the equation for the relative drying rate by G. K. Filonenko, equations are given for determining the density of heat fluxes, the intensity of moisture evaporation and temperature for the period of decreasing drying rate. Solutions of experimental equations are analyzed. The error caused by the processing of experimental data is established. Verification of the reliability of the empirical equations obtained is presented and comparison of the calculated values of the main parameters of the drying kinetics with the experiment is carried out.

Keywords: moisture content, drying coefficient, drying speed, drying intensity, drying of thermal insulation materials

For citation: Ol'shanskii A. I., Golubev A. N. (2023) Investigation of the Kinetics of Heat and Moisture Exchange during Heat Treatment and Drying of Thin Wet Thermal Insulation Materials. Energetika. Proc. CIS Higher Educ. Inst. and Power Eng. Assoc. 66 (1), 66-79. https://doi.org/10. 21122/1029-7448-2023-66-1-66-79 (in Russian)

Введение

Сушка - процесс удаления влаги из твердых тел, связанный с высокими энергозатратами из-за большой теплоты парообразования. На испарение влаги из материала при температуре 90-150 °С расходуется в семь раз больше теплоты, чем на испарение воды. Возможность теоретического расчета длительности сушки влажных материалов в реальных условиях протекания процесса крайне ограничена [1-5]. В этом случае используют экспериментальные исследования с разработкой приближенных эмпирических уравнений для расчета кинетики сушки с минимальным числом постоянных, определяемых из опыта [1-6]. Наиболее эффективны методы обработки опытных данных, основанные на обобщенных переменных, характерных для сушки различных материалов.

Основная часть

Кинетика сушки. Исследование сушки конкретного материала начинается с построения кривых: сушки u = f (т), ее скорости du ¡à т =

= f (u) и температурных Т = f (т), Т = f (u), где u и Т - влагосодержа-

ние и среднеинтегральная температура материала. Вид этих графических зависимостей определяется свойствами материала [1, 3, 5].

На рис. 1 показаны типичные для тонких плоских материалов кривая

сушки u = f (т) и температурные кривые tп = f (т), t = f (т), построенные для войлочной пластины, где tn - температура на поверхности материала.

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

до скорости сушки в первом периоде N, где N = du/dт = const и влагосодержание линейно уменьшается во времени т. При сушке тонких материалов наблюдается период постоянной температуры на уровне /м.т, совпадающий с периодом постоянной скорости сушки до значения критического влагосодержания мкр (рис. 1). Время сушки в первом периоде т1.

1,2

1,0

0,8

0,6

0,4

0,2

Рис. 1. Кривая сушки u = f (т) и температурные кривые tTl = ft), t = f (т) в процессе сушки войлочной пластины (режим сушки: температура среды tc = 90 °C; скорость воздуха v = 5 м/с; относительная влажность воздуха ф = 5 %)

Fig. 1. Drying curve u = f (т) and temperature curves tTl = ft), and t = f (т) during felt plate drying process in drying mode: medium temperature tc = 90 °C, air speed v = 5 m/s, relative humidity of the air ф = 5 %

При достижении влагосодержания икр начинается второй период - падающей скорости сушки. Температура увеличивается от 1шт до температуры среды при длительности сушки Т ^ да. Кривая сушки асимптотически к концу процесса приближается к равновесному влагосодержанию ир,

сушку необходимо заканчивать при кондиционном значении икон > ир во избежание пересушки материала и перерасхода энергии [1, 2]. Исследованиями установлено, что кривая сушки во втором периоде представляет собой

две, а иногда и три сложные сопряженные кривые в точках с влаго-содержаниями йкр2 и йкр3 [1, 3].

К основным кинетическим параметрам относятся время и скорость сушки, температура материала, плотность тепловых потоков, интенсивность испарения влаги. Основная задача кинетики - определение длительности процесса [1-6], при этом аналитический путь очень сложен и не всегда возможен [1, 5].

Зависимость скорости сушки от влагосодержания dü/d т = f (й) во втором периоде имеет сложный вид. По классификации А. В. Лыкова [1-3, 5], для самых разнообразных материалов наблюдаются всего шесть типов кривых скорости сушки. В простейшем случае для расчета длительности сушки применяют линейную аппроксимацию кривой скорости сушки при постоянном коэффициенте K = const по однозональному методу А. В. Лыкова [1-7].

На рис. 2 показаны кривая скорости сушки для войлочной пластины и ее аппроксимация однозональным методом А. В. Лыкова (I), двухзональным методом В. В. Красникова (II) и зональным методом (III) при переменном значении коэффициента K с заменой сложной кривой скорости прямыми. При Ай « 0,1 будет шесть зон, где 1-6 - зоны, в каждой из которых K = const.

0.05

0.04

0,03

0,02

001

0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 u

Рис. 2. Кривая скорости сушки du/dТ = f (U) войлочной пластины

(режим указан на рис. 1): I - аппроксимация методом А. В. Лыкова; II - методом В. В. Красникова; III - зональным методом; 1-6 - аппроксимация по зонам

Fig. 2. Curve of the drying rate du/dт = f (U) in the process of drying the felt plate

(the drying mode is shown in Fig. 1): I - approximation by the method of A. V. Lykov; II - approximation by the method of V. V. Krasnikov; III - approximation by the zonal method;

1-6 - approximation by zones

Уравнение скорости сушки для однозонального метода (рис. 2, прямая I)

du

d т

K (й - uv ).

(1)

При интегрировании (1) при K = const получим длительность сушки в периоде падающей скорости [1-3, 51

!--1П

К

С и — ир ^

и — и

V кР Р У

(2)

где и - текущее влагосодержание материала; тп отсчитывается от нуля.

Длительность сушки с учетом времени в первом периоде по однозо-нальному методу А. В. Лыкова

1 N

(и0 — икр ) —(икр — ир )1П

^ и — ир ^

и — и

V кр р У

Коэффициент сушки К по однозональному методу [1-3, 5]

К

N

икр — ир

■ = СОП81;.

(3)

(4)

Длительность сушки по двухзональному методу В. В. Красникова (рис. 2, ломаная II) определяется уравнением (ир = 0)

Т = Тт

и0 — икр

N

_±1п + -!]п ^

К

и

К

и

(5)

где и1, и11 - текущее влагосодержание соответственно в первой и второй части второго периода

икр2 < «I < икр1; ир < «П < икр2.

Коэффициенты сушки К1 , К2 в первой и второй частях второго периода численно равны тангенсам углов наклона отрезков прямых к оси т. Величина этих коэффициентов зависит от вида влажного тела и скорости сушки в первом периоде N.

В отличие от зональных методов расчета кривой скорости сушки, метод расчета кинетики сушки на основе единого кинетического уравнения, предложенный Б. С. Сажиным [4], позволяет описывать весь процесс, включая первый и второй периоды сушки [3, 4, 8, 9]. В его основу положено обобщенное уравнение массопередачи, которое применительно к конвективной сушке имеет вид [3, 4, 8, 9]

ёи

ё Т

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

К (и0 — и )(ипр — ир ) ,

(6)

где К - эффективный коэффициент сушки.

Интегрирование (6) от начального влагосодержания и0 с учетом влаго-содержания прогрева материала и дает длительность сушки до заданного конечного влагосодержания [4]

К (и0 — ир )

г.

Параметр X по методу Б. С. Сажина определяется выражением [4, 7]

т

и

0

Z = ln (U -U- ). (8)

(U0 - ипр )(U - up )

В системе координат Z = f (т) кривые кинетики представляют собой прямые, угол наклона которых к оси абсцисс определяет эффективный коэффициент сушки K, а отрезки, отсекаемые по оси Z, равны Z0 = U.

Коэффициент сушки Kравен тангенсу угла наклона прямых Z = f (т) [4, 7]

1 Z Zo

tg у = ln-0 « K = const.

т

В реальных условиях K изменяется, поскольку кривая сушки представляет собой две или три сложные экспоненты, сопряженные в точках с влагосодержанием U 2, U 3. Закономерности изменения K вытекают из обобщенной кривой сушки, впервые построенной Г. К. Филоненко в системе координат (U0 - ир ) = f (Nт) [6]. При такой системе обработки

опытных данных все кривые сушки при всех режимах сводятся в одну кривую. Комплекс Nt назван обобщенным временем сушки [1, 2]. В. В. Красниковым и Г. К. Филоненко [2, 6] на основе многих экспериментальных данных показано, что при начальном влагосодержании u0 при любом режиме сушки сохраняется неизменной величина Nt = const. Из обобщенной кривой сушки следует уравнение для кривой скорости сушки [2, 6]

dU

d т

= N * N. (9)

Следовательно, скорость сушки в любой момент второго периода можно определить по величине N для любого режима сушки и величине относительной скорости N [2]. Влияние режима сушки на К отражается величиной скорости в первом периоде N, которая является обобщенной переменной, устанавливающей влияние всех параметров на процесс. Относительная скорость сушки N не зависит от режима и для конкретного материала при данном методе является лишь функцией влагосодержа-ния [1-3], что вытекает из метода обобщения кривых сушки Г. К. Филоненко [2-6]:

1 (Ы ~ ,_,

N =^~Т = / (ы )•

N (т

С целью установления действительной зависимости N = f (^Ыкр)

проведено графическое дифференцирование более 20 кривых сушки для исследуемых материалов с анализом кривых скорости в широком диапазоне режимов = 90-120 °С; V = 3, 5 и 10 м/с). В результате получена зависимость, которая приближенно аппроксимируется уравнением

N *

f - Л и

и

V кР У

(10)

Все кривые скорости сушки для исследуемых материалов относятся к третьему типу кривых по классификации А. В. Лыкова [1-5] и обращены выпуклостью к оси т (рис. 2). Это характерно для капиллярно-пористых материалов [1-4].

На основании (1), (9) и (10)

K !

NN

и - и

N (и/иКр У'

и - и

(11)

р р

Следовательно, в процессах сушки K является переменной величиной.

Кривую скорости сушки для войлочной пластины (рис. 2, линия III) разделим на несколько зон, для каждой из которых K является величиной постоянной, а влагосодержание изменяется по линейному закону. Таким образом, кривая скорости сушки заменяется одной ломаной линией со ступенчатыми прямыми на выбранных промежутках Дм, с большей точностью аппроксимирующей действительную кривую скорости.

Подставим в (2) значения K из (11), тогда с учетом времени сушки в первом периоде длительность процесса при переменных K

1

т = — N

(и0 -икр )-"

и - и„

N *

ln

^ и - и

и - и

V кР р у

(12)

где N найдем по (10).

Ключевым звеном при определении длительности сушки при интегрировании уравнения скорости (1) является коэффициент К. Для однозонального метода относительная скорость сушки

N * =-

и - и„

икр - ир

(13)

В табл. 1 приведены результаты расчета длительности сушки зональными методами и методом Б. С. Сажина.

Результаты расчета длительности сушки войлочной пластины разными методами аппроксимации кривой скорости сушки (рис. 2)

The results of calculating the duration of drying a felt plate using different methods for approximating the drying rate curve (Fig. 2)

Таблица 1

Пластина шерстяного войлока: 135 х 100 х 8 мм; р = 200 кг/м3; и. = 1,14; и = 0,75 ; и = 0, N = 0,054 мин-1. Режим сушки: ¿с = 90 °С; V = 5 м/с; ф = 6 %

Влагосодержание и 0,60 0,50 0,40 0,30 0,25 0,20 0,10

и / и кр 0,80 0,67 0,53 0,40 0,33 0,27 0,14

Коэффициент сушки

Однозональный метод А. В. Лыкова K (4) N , K «-= const; K = 0,072 мин-1; и = 0,75 — ' ' ' кр ' и - и Ч> р

Продолжение табл. 1

Двухзональный метод В. В. Красникова Кь К2 (5) К = 0,065 мин-1; К2 = 0,115 мин-1; u , = 0,75 ; u = 0,23 15 5 L 5 5 кр1 5 5 кр2 5

Метод Б. С. Сажина К (6) К = const; К = 0,115 мин-1; й = 0,75 ' 5 ' кр >

Зональный метод К = var (11) 0,063 0,054 0,047 0,038 0,024 0,023 0,022

Относительная скорость сушки

, u - u N =-^ u - u кр р (13) 0,80 0,67 0,53 0,40 0,33 0,27 0,14

(- v-6 , u N = — экс u V кр / (10) 0,700 0,526 0,362 0,231 0,169 0,124 0,045

Длительность сушки т, мин

Однозональный метод А. В. Лыкова (3) 11,9 14,2 16,8 21,2 23,5 27,6 37,8

Двухзональный метод В. В. Красникова (5) 11,4 13,3 17,4 21,6 24,6 28,4 38,7

Метод Б. С. Сажина (7) 11,2 13,3 16,8 21,6 23,5 26,8 38,2

Зональный метод (12) 12,8 13,3 16,5 20,5 24,5 28,8 38,2

Эксперимент - 12,5 13,5 16,5 20,5 25,0 28,5 38,0

Анализ расчетов времени сушки, представленных в табл. 1, показал следующее: по однозональному методу А. В. Лыкова (рис. 2, прямая I) погрешность в расчетах 7-10 %; по двухзональному методу В. В. Красникова (рис. 2, ломаная II) - 6-8 %; по зональному методу (рис. 2, линия III) -4-5 %; по методу Б. С. Сажина - 5-7 %. Таким образом, более точным является метод при переменных значениях коэффициента сушки К.

Обобщенные комплексные переменные. К обобщенным переменным сушки относятся: N - обобщенное время сушки, тп/тI - относительное

время сушки, й/мкр - относительное влагосодержание, N - относительная

скорость сушки [1-2]. Обработка опытных данных комплексными переменными создает условия инвариантности, позволяющие переходить от одних переменных к другим, от одной системы координат к другой, что значительно сокращает число экспериментов [2, 10-12]. Использование комплекса N придает исследованию сушки обобщенный характер, когда анализируется не единичный частный случай, а множество подобных совокупностей [2, 3].

На рис. 3 даны зависимости Ктп = /(й/йкр) для периода падающей

скорости в процессах сушки керамики, войлока и асбеста при различных режимах.

В результате обработки экспериментальных данных методом наименьших квадратов получено уравнение [10, 11]

N ijj =Р exp

f - \

~ ~ u -2,2—

V ukp J

Длительность сушки с учетом времени в первом периоде

Тт + ттт « —

1 тт N

(u0 -икр) + Pexp

-2,2

_ > и

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

кр

(15)

Рис. 3. Зависимость Nxn = f (ujuKp )

для процессов сушки: 1 - керамики; 2 - войлока; 3 - асбеста; режим сушки: tc = 90, 120 °C; v = 3, 5, 10 м/с

Fig. 3. Dependence of the generalized drying

time Nijj = f (й/йкр ) for the processes

of drying: 1 - ceramics, 2 - felt; 3 - asbestos; drying modes are: tc = 90, 120 °C; v = 3, 5, 10 m/s

Коэффициент в для исследуемых материалов дан в табл. 2.

Значения коэффициентов р, A, A1 в уравнениях (14), (16), (17) в процессах сушки исследуемых материалов

The magnitudes of the coefficient p, A1 and A in equations (14), (16), (17) in the processes of drying of the materials under study

Таблица 2

T

Наименование материала в (14) A (16) A1 (17)

Керамическая плитка 0,36 1 0,55U + 0,19 икр 0,70 - 0,10-^ икр

Листовой асбест 1,76 1 0,39U + 0,18 икр и 0,50 - 0,10-°- икр

Шерстяной войлок 0,65 1 0,32U + 0,14 икр и 0,75 - 0,10-0- икр

На рис. 4 представлена обработка опытных данных для тп/т: = = / (ц/икр) в процессах сушки керамики, асбеста и войлока при различных режимах в виде экспоненциальной и степенной зависимостей [13-16]:

— « A exp тт

-2u

, ^р J

^ ж a,

TI

/ _ \-U5

u

% J

(16)

(17)

где коэффициенты A, A, - нелинейные функции отношений соответственно

Tii/TI

2,5 2,0 1,5 1,0 0.5

0,1

0,2

0,4

0,5

0,6

0,7

0,8 й/«кр

Рис. 4. Зависимость тп/т1 = f (ujuKp) в процессах сушки: 1 - керамики; 2 - асбеста; 3 - войлока (режимы указаны на рис. 3)

Fig. 4. Dependence тп/т1 = f (w/w^) in the processes of drying: 1 - ceramics; 2 - asbestos; 3 - felt (drying modes are given in Fig. 3)

0,2

0,4

0,6

0,8 uu

1,0

Рис. 5. Зависимость коэффициентов A = f (u/uKp) (a), A, = f (ы0/мкр) (b) для сушки:

1 - керамики; 2 - асбеста; 3 - войлока (режимы указаны на рис. 3) Fig. 5. Dependence of coefficients A = f (uju^) (a), A, = f (ы0/ыкр) (b) for drying of: 1 - ceramics; 2 - asbestos; 3 - felt (drying modes are given in Fig. 3)

Коэффициент A, в первом приближении вычисляется по зависимости u

A, « 1,2 - 0,38з^. Более точные уравнения для определения коэффициен-

"кр

тов A, Ai даны в табл. 2.

b

a

0

Тогда уравнение (15) с учетом (16), (17) примет вид

1 + А ехр

т * т: + тп

и0 - икр

N

и0 - икр

(

-2-Ц-

V икр у

N

1 + А

( _ Л-1'25

и

V икр у

Основное уравнение кинетики сушки А. В. Лыкова [1, 3, 8]

^ = N * (1 + ЯЪ), 4:

(18)

(19)

(20)

где 4:, 4П - плотность потока теплоты соответственно в первом и втором периодах.

При сушке тонких материалов расход теплоты на нагревание влажного тела значительно меньше, чем на испарение из него влаги, и число Ребин-дера ЯЪ << 1 [1, 2]. Пренебрегая числом Ребиндера при малых его значениях [1, 8], для плотности потока теплоты во втором периоде запишем

На основании (10) получим 4п * 4:

( _ у,6

и

и

V Ч> у

4:: = 4:N . = Ф ^

( _ у,6

и

и

V Ч> у

(21)

где 4: = грЯ^; г - теплота парообразования; р - плотность; Яу - отношение объема тела к поверхности.

Интенсивность испарения влаги в периоде постоянной скорости сушки [1-4]

7: = РЯД, (22)

а в периоде падающей скорости

4: / \1,6

7:: =

4::

( - ^ 1,6

и

V кР у

г

Г - \ и

V икр у

(23)

Г. К. Филоненко на основании обобщенной кривой сушки (9) получил зависимость, устанавливающую связь влагообмена с теплообменом [6]:

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

N*

'о -*и

' - 7

о м.т

Температура на поверхности материала вычисляется по формуле

'п = 'с -('с - 'м.т ) N * = 'с -('с - 'м.т )

г - V,6

и

и

V кр у

(24)

(25)

х ~ х: + х::

В табл. 3 приведены результаты расчетов основных параметров сушки по полученным уравнениям и сравнение расчетных значений с экспериментом.

Таблица 3

Результаты расчета основных параметров кинематики сушки по полученным уравнениям и их сравнение с экспериментальными значениями

The results of calculating the main parameters of the kinematics of drying according to the obtained equations and comparing the calculated values with the experimental ones

Влаго-содержание u Относительное влаго-содержание u икр Длительность сушки т, мин Температура поверхности Гп, °C Плотность потока теплоты во втором периоде qn, Вт/м2 Интенсивность испарения влаги в периоде падающей скорости j'jj, 103 кг/м2с

(экс) (7) (12) (15) (18) (19) (экс) (25) (21) (23)

Керамическая плитка: 120 х 78 х 5 мм; р = 1860 кг/м3; и0 = 0,2; uкр = 0,1; режим сушки: ^ = 120 °С; V = 5 м/с; ф = 5 %; N = 0,025 мин-1; ^ = 48,5 °С; т, = 4 мин; д1 = 7200 Вт/м2;= 3,5-10-3 кг/м2с

0,08 0,8 6,5 7,2 6,4 6,1 7,1 6,8 65,0 64,0 6380 2,66

0,06 0,6 8,5 8,6 8,3 8,9 8,4 8,8 75,0 73,5 4770 1,98

0,04 0,4 11,5 11,3 11,4 12,2 11,2 11,8 83,0 85,2 9260 1,37

0,02 0,2 16,5 17,2 16,5 15,9 18,6 17,2 92,0 90,5 1690 0,71

Листовой асбест: 120 х 80 х 6 мм; р = 770 кг/м3; м0 = 0,46; мкр = 0,21; режим сушки: ^ = 120 °С; V = 5 м/с; ф = 5 %; N = 0,028 мин-1; ^ = 37 °С; т, = 7,5 мин; д1 = 4840 Вт/м2; ]\= 2,15-10-3 кг/м2с

0,16 0,8 11,5 11,2 11,4 11,3 11,4 11,2 55,0 57,0 3860 1,61

0,14 0,7 12,5 12,4 12,5 12,2 12,8 12,3 59,0 61,2 3420 1,42

0,12 0,6 14,5 15,3 14,5 14,2 13,9 14,2 65,0 67,2 2830 1,18

0,08 0,4 18,5 20,5 18,4 18,3 19,2 19,2 76,0 79,4 1890 0,79

0,04 0,2 24,5 25,6 23,8 26,9 25,7 25,4 85,0 80,7 1050 0,40

0,02 0,1 29,0 29,8 28,7 31,4 32,3 30,8 98,0 101 680 0,26

Шерстяной войлок: 135 х100 х 8 мм; р = 220 кг/м3; м0 = 1,14 ; икр = 0,75 ; режим сушки: ^ = 90 °С; V = 5 м/с; ф = 5 %; N = 0,054 мин-1; = 34 °С; т, = 7,5 мин; = 3350 Вт/м2; ]1 = 1,38-10-3 кг/м2с

0,60 0,80 12,5 11,2 12,8 12,4 12,1 11,8 41,0 40,8 2650 1,07

0,50 0,66 13,5 13,2 13,3 14,2 14,2 13,2 43,0 42,2 2150 0,89

0,40 0,53 16,5 16,8 16,5 16,2 16,2 16,7 46,0 45,7 1810 0,75

0,30 0,40 20,5 21,6 20,5 19,5 21,3 21,4 51,0 50,0 1440 0,61

0,20 0,26 28,5 26,8 28,8 27,6 30,2 28,5 58,0 57,6 840 0,34

ВЫВОДЫ

1. Проведено исследование кинетики сушки тонких теплоизоляционных плоских материалов. Рассмотрены способы обработки опытных данных обобщенными комплексными переменными и методы расчета длительно-

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

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

ЛИТЕРАТУРА

1. Лыков, А. В. Теория сушки / А. В. Лыков. М.: Энергия, 1968. 472 с.

2. Красников, В. В. Кондуктивная сушка / В. В. Красников. М.: Энергия, 1973. 288 с.

3. Акулич, П. В. Расчеты сушильных и теплообменных установок / П. В. Акулич. Минск: Беларус. наука, 2010. 443 с.

4. Сажин, Б. С. Научные основы техники сушки / Б. С. Сажин. М.: Наука, 1997. 447 с.

5. Рудобашта, С. П. Массоперенос в системах с твердой фазой / С. П. Рудобашта. М.: Химия, 1980. 248 с.

6. Филоненко, Г. К. Сушильные установки / Г. К. Филоненко, П. Д. Лебедев. М. - Л.: Гос-энергоиздат, 1952. 263 с.

7. Лыков, А. В. Кинетика теплообмена в процессе сушки влажных материалов / А. В. Лыков, П. С. Куц, А. И. Ольшанский // Инженерно-физический журнал. 1972. Т. 23, № 3. С. 401-406.

8. Ольшанский, А. И. Кинетика тепловлагообмена и метод расчета длительности конвективной сушки натуральной кожи / А. И. Ольшанский, А. Л. Климентьев // Вес. Нац. акад. навук Беларуси Сер. фiз.-тэхн. навук. 2020. Т. 65, № 4. С. 464-475. https://doi.org/10. 29235/1561 -8358-2020-65-4-464-475.

9. Куц, П. С. Некоторые закономерности тепловлагообмена и приближенные методы расчета кинематики процесса сушки влажных материалов / П. С. Куц, А. И. Ольшанский // Инженерно-физический журнал. 1977. Т. 32, № 5. С. 1007-1014.

10. Ольшанский, А. И. Экспериментальное исследование кинетики сушки тонких плоских влажных материалов методом регулярного режима с использованием обобщенных комплексных переменных / А. И. Ольшанский, А. М. Гусаров // Инженерно-физический журнал. 2017. Т. 90, № 3. С. 700-713.

11. Ольшанский, А. И. Регулярный тепловой режим нагревания влажных плоских капиллярно-пористых материалов в процессе сушки / А. И. Ольшанский, В. И. Ольшанский // Инженерно-физический журнал. 2014. Т. 87, № 6. С. 1308-1318.

12. Ольшанский, А. И. Исследование теплообмена в процессе сушки влажных материалов по экспериментальным данным влагообмена с использованием обобщенных комплексных переменных / А. И. Ольшанский // Инженерно-физический журнал. 2014. Т. 87, № 4. С. 887-897.

13. Ольшанский, А. И. Экспериментальные исследования тепловлагообмена в процессе конвективной сушки тонких влажных материалов / А. И. Ольшанский, С. В. Жерносек, А. М. Гусаров // Энергетика. Изв. высш. учеб. заведений и энерг. объединений СНГ. 2018. Т. 61, № 6. С. 564-578. https://doi.org/10.21122/1029-7448-2018-61-6-564-578.

14. Ольшанский, А. И. Исследование теплообмена в процессе сушки методом регулярного режима / А. И. Ольшанский, В. И. Ольшанский // Инженерно-физический журнал. 2012. Т. 85, № 2. С. 385-391.

15. Ольшанский, А. И. Исследование тепломассообмена в процессах тепловой обработки и сушки теплоизоляционных материалов / А. И. Ольшанский, С. В. Жерносек, А. М. Гусаров // Энергетика. Изв. высш. учеб. заведений и энерг. объединений СНГ. 2022. Т. 65, № 2. С. 156-168. https://doi.org/10.21122/1029-7448-2022-65-2-156-168.

16. Лыков, А. В. Теория тепло- и массопереноса / А. В. Лыков, Ю. А. Михайлов. М. - Л.: Госэнергоиздат, 1963. 535 с.

Поступила 09.04.2022 Подписана в печать 14.06.2022 Опубликована онлайн 31.01.2023 REFERENCES

1. Lykov A. V. (1968) Theory of Drying. Moscow, Energiya Publ. 472 (in Russian).

2. Krasnikov V. V. (1973) Conductive Drying. Moscow, Energiya Publ. 288 (in Russian).

3. Akulich P. V. (2010) Calculations of Drying and Heat Exchange Installations. Minsk, Belaruskaya Navuka Publ. 443 (in Russian).

4. Sazhin B. S. (1997) Scientific Fundamentals of Drying Techniques. Moscow, Nauka Publ. 447 (in Russian).

5. Rudobashta S. P. (1980) Mass Transfer in Systems with a Solid Phase. Moscow, Khimiya Publ. 248 (in Russian).

6. Filonenko G. K., Lebedev P. D. (1952) Dryers. Moscow - Leningrad, Gosenergoizdat Publ. 263 (in Russian).

7. Lykov A. V., Kuts P. S., Olshanskii A. I. (1972) Kinetics of Heat Transfer in the Process of Drying Wet Materials. Journal of Engineering Physics, 23 (3), 1082-1086. https://doi. org/10.1007/bf00832214.

8. Alshansky A. I., Klimentyev A. L. (2020) Kinetics of Heat and Moisture Exchange and Method for Calculating the Duration of the Convective Drying Process of Natural Leather. Vestsi Natsyyanal 'nai Akademii Navuk Belarusi. Seryya Fizika-Technichnych Navuk = Proceedings of the National Academy of Sciences of Belarus. Physical-Technical Series, 65 (4), 464-475. https://doi.org/10.29235/1561-8358-2020-65-4-464-475 (in Russian).

9. Kutz P. S., Ol'shanskii A. I. (1977) Some Features of Heat and Moisture Transfer and Approximate Methods of Calculating the Drying Kinetics of Moist Materials. Journal of Engineering Physics, 32 (6), 650-656. https://doi.org/10.1007/bf00862568.

10. Ol'shanskii A. I., Gusarov A. M. (2017) Experimental Study of the Kinetics of Drying Thin Flat Moist Materials by the Regular-Regime Method Using Generalized Complex Variables. Journal of Engineering Physics and Thermophysics, 90 (3), 665-678. https://doi.org/10.1007/ s10891-017-1614-x.

11. Ol'shanskii A. I. (2014) Regular Heat Regime of Heating of Moist Capillary-Porous Materials in the Process of Their Drying. Journal of Engineering Physics and Thermophysics, 87 (6), 1362-1373. https://doi.org/10.1007/s10891-014-1139-5.

12. Ol'shanskii A. I. (2014) Study of the Heat Transfer in the Process of Drying of Moist Materials from Experimental Data on Moisture Transfer. Journal of Engineering Physics and Thermophysics, 87 (4), 917-928. https://doi.org/10.1007/s10891-014-1089-y.

13. Ol'shanskii A. I., Zhernosek S. V., Gusarov A. M. (2018) Experimental Studies of Heat and Moisture Exchange in the Process of Convective Drying of Thin Wet Materials. Em^tika. Izvestiya Vysshikh Uchebnykh Zavedenii i Energeticheskikh Ob 'edinenii SNG = Energeti-ka. Proceedings of CIS Higher Education Institutions and Power Engineering Associations, 61 (6), 564-578. https://doi.org/10.21122/1029-7448-2018-61-6-564-578 (in Russian).

14. Ol'shanskii A. I., Ol'shanskii V. I. (2012) Investigation of Heat Transfer in the Drying Process by the Method of Regular Regime. Journal of Engineering Physics and Thermophysics, 85 (2), 415-422. https://doi.org/10.1007/s10891-012-0667-0.

15. Ol'shanskii A. I., Zhernosek S. V., Gusarov A. M. (2022) Investigation of Heat and Mass Transfer in the Processes of Heat Treatment and Drying of Thermal Insulation Materials. Ern^ti^. Izvestiya Vysshikh Uchebnykh Zavedenii i Energeticheskikh Ob 'edinenii SNG = Energetika. Proceedings of CIS Higher Education Institutions and Power Engineering Associations, 65 (2), 156-168. https://doi.org/10.21122/1029-7448-2022-65-2-156-168 (in Russian).

16. Lykov A. V., Mikhailov Yu. A. (1963) Theory of Heat and Mass Transfer. Moscow - Leningrad, Gosenergoizdat Publ. 535 (in Russian).

Received: 9 April 2022

Accepted: 14 June 2022

Published online: 31 January 2023

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