Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Лебедев Д.П. Тепло- и массообмен в процессах сублимации в вакууме

.pdf
Скачиваний:
29
Добавлен:
25.10.2023
Размер:
14.77 Mб
Скачать

заполненных газом. S-образный характер кривых зави­ симостей теплопроводности от давления можно объяс­ нить различным вкладом теплопроводности газа в об­ щую теплопроводность капиллярно-пористого тела в про­ цессе сублимации при изменении давления в камере и ориентировочно установить следующие области: 1 —

Рис. 2-18. Зависимость теплопроводности пищевых про­ дуктов от давления.

1 — говядина; 2 — персик.

область континуума (760—100 мм рт. ст.); 2 — область газовой проводимости (100—5- 10~2 мм рт. ст.)\ 3 — область кондуктивной проводимости (от 5-10~2 мм рт. ст.).

Ширина представленных областей определяется структурой каркаса капиллярно-пористого тела и может изменяться в зависимости от диаметра пор (числа Кнудсена).

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

50

со скольжением. Резкое уменьшение теплопроводности

при давлениях в камере ниже 1 мм рт.

ст., очевидно,

является следствием уменьшающегося

вклада пара

в перенос тепла,

о котором упоминалось выше.

Как видно из

рис. 2-18, смена режимов течения пара

внутри капиллярно-пористого тела в значительной сте­ пени изменяет теплопроводность пористого каркаса.

Для всех газов, использованных в [Л. 2-39], тепло­ проводность достигает постоянного значения для дав-

Рис. 2-19. Расчетная структура пористого материала.

а — расчетная модель [Л. 2-22); б — круглая структура (куриный белок); в — ячеечная структура (картофельный крахмал); г — расчетная модель [Л. 2-58].

лений ниже 0,02 мм рт. ст. При этих давлениях газ не дает значительного вклада в теплопроводность, и, сле­ довательно, практически все тепло проводится твердой структурой образца. Рассматриваемые результаты рабо­ ты [Л. 2-39] были получены при отсутствии течения газа через поры, в то время как при реальной сублимацион­ ной сушке имеет место поток водяного пара через вы­ сушенный слой продукта в вакуум.

Для теоретического расчета в [Л. 2-22, 2-44] была предложена модель (рис. 2-19,а) пористого материала, положенная в основу вывода уравнения, описывающего зависимость теплопроводности от давления.

4*

51

Теплопроводность рассматриваемой системы описы­ вается уравнением

1 — d

1— П

П

(2-38)

 

А5 (1 - Я ) + А р/ 7 + 8^ Г + 8 Ар

 

 

где Xs — теплопроводность

твердой

части

материала;

ХР — теплопроводность

пустой и

пористой

части;

П

пористость материала;

б — часть

слоя

материала,

пред­

ставляющая гидродинамическое

сопротивление потоку

пара.

ХР = ХЯ + ХТ,

где

XR— эквивалент­

Теплопроводность

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

Для большинства разностей температур эквивалент­

ный

коэффициент теплопроводности (для засыпки из

стеклянных шариков)

 

 

 

i R=

c ^ d . o m ( M ,

(2-39)

где

Тт =0,5(Гі + 7'2) ;

d — среднее расстояние

между

стенками

пор; ф « 1 — угловой коэффициент;

С12 = 3,5—

4,5 ккал/(м2ч- град) — коэффициент

излучения; 7\ и

Т2— температуры противолежащих

стенок

единичного

капилляра

(поры).

 

 

 

Теплопроводность ХГопределяется как

 

 

Я.,. —■-

А,

2l/d 2 — а '

(2-40)

 

1 +2 .Кп ■

 

 

 

где Аго — теплопроводность газа, не зависящая от давле­

ния;

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

Кп — число

Кнуд­

сена;

I —длина свободного пробега

молекул газа

в по­

рах материала.

(область низкого вакуума и большие

Если

Кп<1

диаметры пор материалов), то X не зависит от давления;

если

же

Кп>1,

что имеет место при высоком вакууме

и малых диаметрах пор материала, то X зависит от дав­ ления газа, заполняющего поры.

Расчетные данные по уравнению (2-38) для коэффи­

циента теплопроводности

трех исследуемых

веществ

в зависимости от давления

приведены на рис.

2-20,а, и

5 2

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

Систематические отклонения вычисленных по форму­ ле (2-38) кривых, у которых падающая ветвь все время проходит несколько более круто, чем у эксперименталь­ ных, можно объяснить тем, что действительные поры материала имеют различную величину, для расчетов бралась средняя величина диаметра пор.

Рис. 2-20. Влияние вакуума на коэффициенты теллоіпіртшстости Я и диффузии D'.

а — зависимость коэффициента теплопроводности от

давления: сплошные кри­

вые — эксперимент;

пунктирные кривые — расчет;

^ —стеклянные

шарики

диаметром 0,305 мм;

О — стеклянные

шарики

диаметром 1,88 мм; X — мате­

риал

Ютсжг;

б — влияние

неконденеирусмых

газов

на коэффициенты

тепло­

проводности и

диффузии;

пунктирные

кривые — коэффициент теплопроводно­

сти;

сплошные

кривые — коэффициент

диффузии; 1 — гелий; 2 — азот.

 

Исходя непосредственно из микроскопического ана­ лиза структуры высушенных сублимацией пищевых про­ дуктов, в [Л. 2-58] предлагаются их наиболее типичные модельные ячейки (рис. 2-19,6, в), для которых в [Л. 2-59] в соответствии с рис. 2-19,г выведена формула для рас­ чета коэффициента теплопроводности:

X= сЯГ bXs -f-

( ^ 1)

где Яг и Xs — теплопроводности газа и твердой части материала.

Как показано в [Л. 2-39, 2-58], состав газа, заполняю­ щего поры материала, может существенно изменять его теплопроводность. Например, для мяса, начиная с ва­ куума р> 1 мм рт. ст., наблюдается увеличение коэффи-

53

циента теплопроводности при заполнении пор продукта азотом и гелием (рис. 2-20,6).

Теплопроводность тел, имеющих капиллярно-пори­ стую волокнистую структуру, может быть различной в зависимости от направления волокон; например, в ра­ ботах [Л. 2-55 и 2-56] были проделаны эксперименты по определению теплопроводности мяса при тепловом по­ токе, как параллельном, так и перпендикулярном его волокнам. Содержание влаги в образцах при этом было почти одинаковым.

Вобразцах, где тепловой поток был направлен па­ раллельно волокнам, теплопроводность была на 10— 16%’ выше, чем в образце, где тепловой поток перпен­ дикулярен.

Втабл. 2-8 представлены значения теплопроводности некоторых пищевых продуктов, высушенных сублима­ цией.

Таблица 2-8

Коэффициент теплопроводности пищевых продуктов, высушенных сублимацией

р, мм

Материал

 

бтІ(М'град)

Источник

pm. cm.

 

0,02

Говядина,

высушенная

0,0375

[Л.

2-46]

0,17

сублимацией

 

0,057

[Л.

2-39]

Необработанная говядина,

 

высушенная

сублима­

 

 

 

0,17

цией

 

 

0,065

[Л.

2-39]

Приготовленная говядина,

 

высушенная

сублима­

 

 

 

0,145

цией

 

 

 

 

 

Лосось

 

 

0,0145

I

2-37]

0,08

Треска

 

 

0,019

[Л.

0,08

Окунь

фаршированный,

0,0225

J

 

0,4

Судак

0,045—0,06

1

2-16]

0,29

*==(—15-Н+5) °С

0,045—0,052

[Л.

Судак

бланшированный

J

 

 

(кусочки)

 

 

 

 

t = (—20н— 5)

0,0256—0,034

[Л.

2-37]

Соки (влажность

5%)

Эффективный коэффициент диффузии. На рис. 2-20,6 показано изменение эффективного коэффициента диффу­ зии D' в зависимости от вакуума для пищевого продукта (мяса), заполненного гелием и азотом (Л. 2-39].

Анализ рис. 2-18, 2-20,а, б показывает, что скорость сублимационной сушки, очевидно, должна быть с одной

54

стороны ограничена при давлениях ниже 1 мм рт. cf. теплообменом (уменьшение коэффициента теплопровод­ ности материала), а с другой стороны — при давлениях выше 1 мм рт. ст. массопереносом (уменьшение эффек­ тивного коэффициента диффузии). Эти физические усло­ вия механизма процесса необходимо учитывать при установлении режимных параметров сублимационной сушки (особенно для третьего периода сушки).

2-4. СТРУКТУРНЫЕ И ТЕПЛО ФИЗИЧЕС КИЕ СВОЙСТВА ПОРИСТЫХ МЕТАЛЛОВ

Теоретическое исследование процесса двойного фазового перехода вода — лед, лед — пар в проницаемых метал­ лах в процессах сублимации (см. гл. 5) требует всесто­ роннего знания зависимости их теплофизических свойств (теплопроводности, теплоемкости) от целого ряда фак­ торов: структурных характеристик материала, насыще­ ния его влагой, льдом, характера фазового перехода вблизи 0°С и т. п.

Структурно-физические свойства. При выборе метода исследования структурно-физических свойств проницае­ мых пластин принималось во внимание, что процесс сублимации льда в вакууме в пористой металлокерамике происходил при давлении питательной воды в 1 кгс/см2. Это давление определяло минимальный (рабочий) раз­ мер пор, «раскрывающихся» в пластине и участвующих в процессе сублимации.

Воздухопроницаемость. Эксперимент по определению воздухопроницаемости пористых пластин проводился на

установке,

описанной

в (Л.

2-5].

Прокаленный до тем­

пературы

200 °С в

электропечи

образец

продувался

сухим воздухом, при этом

определялась

зависимость

расхода воздуха через проницаемую пластину от дав­ ления.

Водопроницаемость. Многие авторы [Л. 2-1, 2-5, 2-6, 2-17, 2-18] отмечают нестабильность процесса течения жидкости через пористые тела, но, однако, до сих пор нет единого мнения о причинах, которые вызывают это явление. Водопроницаемость металлокерамических пла­ стин практически не изучалась.

Нами была исследована проницаемость пластин из порошков металлов с различной пористостью. Образцы пористых пластин насыщались бидистиллированной во­

дой под вакуумом и исследовались на водопроницае­ мость на установке, описанной в [Л. 2-2]. Пример полу­ ченных данных представлен на рис. 2-21. Анализ экс­ периментальных результатов позволил сделать вывод, что режим течения жидкости в титановых пластинах

близок к ламинарному. В

пластинах

никеля режим те­

 

 

чения

является

 

более

 

 

сложным. Для

определе­

 

 

ния коэффициентов

водо­

 

 

проницаемости

для

тита­

 

 

новых пластин можно ис­

 

 

пользовать

 

уравнение

 

 

Дарси, для никелевых пла­

 

 

стин — двучленные

 

урав­

 

 

нения [Л. 2-2].

 

пор по

 

 

Распределение

 

 

размерам. Для установле­

 

 

ния распределения пор по

 

 

размерам в проницаемых

 

 

пластинах

традиционно

Рис. 2-21.-Фильтрация воздуха и

используются метод ртут­

воды через проницаемую пласти­

ной парометрии

[Л. 2-2,

ну из «икеля 26.

 

2-20],

метод

определения

1 — воздух; 2 — вода.

 

распределения пор по раз­

образца [Л. 2-20] и метод

 

мерам

по

микрошлифам

максимального

давления пу­

зырьков [Л. 2-6].

 

 

 

 

 

 

 

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

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

С целью дальнейшего использования проницаемых пластин после исследования их структурных характери­ стик в теплофизических экспериментах нами был выбран для установления распределения пор по размерам метод максимального давления пузырьков [Л. 2-5]. Этот метод основан на измерении давления, которое необходимо, чтобы вытеснить воздух через поры проницаемой пла­ стины, заполненной жидкостью с известным поверхност­ ным натяжением.

56

В работах |Л. 2-5, 2-7, 2-17, 2-18], применяя метод максимального давления пузырьков и закон Пуазейля, получены диаметры и распределения пор по размерам для целого ряда капиллярно-пористых материалов — огнеупоров, керамики и т. п. Эксперименты проводились на установках, описанных в і[Л. 2-5]. Нами методом максимального давления пузырьков получены для целого ряда проницаемых металлических пластин распределе­ ния нор по размерам на установке, близкой к описанной

см3

см3

смг-сек

смг-сек

Рис. 2-22. Явление гистерезиса у пластин различ­ ной пористости.

а — титан-1А; б — бронза-13; в — алюминий-4.

в [Л. 2-2]. В качестве примера на рис. 2-22 представлен

график

зависимости расхода воздуха от давления W=

= іf(p)

для трех проницаемых пластин

металлокерами­

ки: титана-ІА, бронзы-13 и алюминия-4.

методу, предло­

Обработка полученных кривых по

женному в [Л. 2-6] и модифицированному автором, по­ зволила построить для рассмотренных проницаемых пла­ стин распределение пор по диаметрам n = f(d) и срав­ нить его с максвелловским распределением n — F{r), где г — радиус пор (рис. 2-23).

57

Для титана

ІА было получено уравнение

F (г) =

0,307 [ехр(—0,73* 108 г2)]г2 • ІО15.

Для бронзы 13

F (г) =0,159 [ехр (—0,033 • 108 г^г2• 1012.

Явление гистерезиса. При течении газа через насы­ щенную водой проницаемую пластину наблюдается явле­ ние гистерезиса, заключающееся в отставании изменения

 

 

 

расхода

газа

при увели­

 

 

 

чении и

уменьшении его

 

 

 

движения (рис. 2-22,аи в).

 

 

 

Как видно из рис. 2-22,а

 

 

 

и в,

для

алюминиевых

 

 

 

пластин

площадь

петли

 

 

 

гистерезиса

 

А В—С—

 

 

 

E—DG—F при много­

 

 

 

кратном

изменении

дав­

 

 

 

ления газа растет и в ко­

 

 

 

нечном

счете

стабилизи­

а)

 

руется. Это, по-видимому,

 

объясняется

тем,

что ма­

Рис. 2-23. Распределение пор по

териал имеет большое ко­

личество

пор

со сложной

диаметрам

в проницаемых пла­

конфигурацией.

 

 

стинах и сравнение его с максвел­

из

рис.

ловским

распределением

(пунк­

Как

видно

тирная линия).

 

2-22,6 гистерезис у прони­

а — бронза-13; б — титан-1А.

 

цаемых пластин

из

брон­

зы не наблюдается. По-видимому,

это

в

определенной

степени

связано с тем,

что размер

пор

у бронзы лежал

в пределах 30—90 мкм, в то время как у титана состав­ лял 3—7 мкм.

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

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

Для определения теплофизических свойств пористых металлов нами [Л. 2-61] использовался как стационар-

58

пый метод продольного теплового потока, так и неста­ ционарный метод монотонного разогрева.

Схема экспериментальной установки для определе­ ния теплопроводности пористых металлов стационарным методом показана на рис. 2-24. Исследуемые образцы 9 располагались симметрично между основным нагревате­ лем 4 и теплоприемниками 2. Толщина основного на­ гревателя составляла 1,5 мм.

Рис. 2-24. Прибор для определения X (симметричная схема).

а — конструкция нагревателя.

Использование двух образцов сокращало время экс­ перимента практически в 2 раза и обеспечивало эффек­ тивную тепловую защиту нагревателя.

Контакт образцов с нагревателем производился с по­ мощью винтов 6 через стальной шарик 7. С помощью нагревателя 3 и холодильников 5 на охранном цилиндре 8 поддерживалось такое же распределение температур, что и на испытуемых образцах.

Коэффициент теплопроводности определялся из со­

отношения

 

 

 

 

Я -

V

 

 

(ДГ, + ДT2)F

где

Q — мощность

основного нагревателя; А74—пере­

пад

температур в

первом

образце; АТ2— перепад тем­

ператур во втором

образце

(разность между АTt и АТг

59

Соседние файлы в папке книги из ГПНТБ