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

Литература от Абакумова ч2 / Осмотический массоперенос

.pdf
Скачиваний:
47
Добавлен:
25.03.2016
Размер:
2.3 Mб
Скачать

Рис. 6.13. Кинетика сорбции асбестоцементом радиоактивной метки NaJ* во времени τ (ч) без носителя (1); с носителем (2) (нерадиоактивной солью NaJ); N, N0 – текущие и начальные значения радиоактивности

J131 с носителем NaJ снижает сорбцию самого радиоактивного индикатора (кривая

2).

С целью выявления фазового состава переносимой влаги был применен метод гидрофобной прослойки из гидрофобизированного расплавленным окисленным петролатумом кварцевого песка. Образцы содержали 85% портландцемента и 15% асбеста марки БП-5-65 (табл. 6.2).

Таблица 6.2. Характеристики асбестоцемента и результаты экспериментов *

 

P

u0

γ

В

α

Q 103

δ 103

D 107

 

 

10

0,23

1520

0,280

82

1,78

1.3

4,31

 

 

7

0,26

1520

0,287

89

0,93

1,1

4,61

 

*Р – давление

3,5

0,29

1420

0,312

93

0,27

0,7

5,28

 

прессования, МПа; u0 – начальное влагосодержание; γ – плотность, кг/м3; В

предельная влагоемкость, кг воды/кг смеси; α – степень водопоглощения; Q – интенсивность

потока влаги, кг/м2; δ – термоградиентный коэффициент, К1; D – коэффициент диффузии, м2

 

 

 

 

 

 

На рис. 6.14 приведены графи-

 

 

 

 

 

ки распределения температур. Тер-

 

 

 

 

 

моперенос влаги в направлении,

 

 

 

 

 

перпендикулярном

плоскости асбе-

 

 

 

 

 

стоцементных пленок, незначителен.

 

 

 

 

 

Вызвано это тем, что вдоль пленок

 

 

 

 

 

имеется система сплошных водопро-

 

 

 

 

 

водящих каналов. В поперечном на-

 

 

 

 

 

правлении она отсутствует.

 

 

 

 

 

 

В начальный период нагрева,

 

 

 

 

 

как и в торфе (см. рис. 6.12), наблю-

 

 

 

 

 

дается волнообразное распределение

 

 

 

 

 

влаги (рис. 6.15) в левой половине

 

 

 

 

 

образца, примыкающей к нагревае-

 

 

 

 

 

мому торцу. С увеличением времени

Рис. 6.14. Распределение температуры t

«волны» избыточной влажности пе-

(°С): а – по длине образца в начале опыта

ремещаются в глубь образца. За вре-

(1), через 10 мин (2), 30 мин (3), 1 ч (4), 1

мя нагрева 24 ч по длине образца не

ч 45 мин (5), 3 ч (6); б – во времени τ (ч)

устанавливается линейное распреде-

на горячем (1) и холодном (6) торцах и на

ление влагосодержания.

расстоянии 1 (2), 3 (3), 5 (4), 7 (5) см от

 

 

 

 

 

нагревателя

 

 

 

 

 

 

 

 

185

Рис. 6.15. Кинетика термопереноса влаги в образцах свежесформованного асбестоцемента, имеющих начальное влагосодержание u0 0,29 (а), 0,26 (б), 0,23 (в). Изменение влагосодержания u по длине образца l (см) через 2 (1), 4 (2), 24 (3) часа от начала нагрева и во времени τ (ч) в сечениях на расстоянии 1 (1), 3 (2), 5 (3), 7 (4) см от нагреваемого торца

Температура по мере удаления от нагреваемого торца образца быст-

ро снижается (см. рис.

6.14а). Различие расположения максимумов для

 

 

графиков

радиоактивности и

влажности

 

 

(рис. 6.16, кривая 1 и 2) свидетельствует о

 

 

наличии

смешанного механизма

переноса

 

 

влаги. Метка от горячего торца переносится

 

 

влагой в жидкой фазе. Кроме того, имеет

 

 

место пароперенос. Пар конденсируется в

 

 

крупных порах. Конденсат по макропорам

 

 

перемещается в зону испарения (к горячему

Рис. 6.16. Распределение изме-

торцу) и переносит с собой метку в обрат-

ном направлении. При наличии встречных

нений радиоактивности

N

потоков влаги метка распределяется вблизи

(имп./мин) (1) и влагосодержа-

зоны испарения в образце.

 

ния u (2) по длине l (см) образ-

Были проведены опыты поопределению

ца асбестоцемента после 18 ч

зависимости показаний микротензиометров Рк

нагрева (u0 =0,26)

 

 

 

отвлагосодержания асбестоцемента.

 

 

186

 

 

При влагосодержании u >0,7 значения Рк при комнатной температуре мало изменяются (рис. 6.17). При u <0,7 наблюдается увеличение ка- пиллярно-осмотического давления, которое достигает при u =0,16 величины 0,25 МПа.

Через гидрофобную прослойку (рис. 6.18) влага в жидкой фазе не переносится. В гидрофобных порах влага не адсорбируется. Поры заполнены паром при ϕ→1, что существенно осложняет пароперенос через прослойку. В первой половине (до гидрофобной прослойки) распределение влажности подобно начальному участку кривой 2 (см. рис. 6.15), а во второй – всей кривой 2.

Рис. 6.18. Распределение изменений влагосодержания u по длине l (см) асбестоцементного образца началь-

 

ного влагосодержания 0,26 до (1) и

 

после (1) гидрофобной прослойки

Рис. 6.17. Зависимость потенциала влаги

Согласно табл. 6.2, интенсив-

ности потока влаги Q через сечение

Р (МПа) от влагосодержания асбестоце-

мента u при измерении микротензиомет-

образца на расстоянии 2 см от на-

рами с внутренним размером 0,3 (1), 0,4

греваемого торца и термоградиент-

(2), 1,9 (3) мкм

ные коэффициенты δ уменьшаются

 

с увеличением начального влагосо-

держания u0 асбестоцемента. При больших значениях u0 в нем возрастает степень водонасыщения α=u0/В и соответственно объемы защемленного воздуха. При этом снижается объем открытых, связанных между собой поровых каналов (табл. 6.3). В связи с этим уменьшается пароперенос к холодным торцам образцов и массоперенос конденсированной влаги в обратном направлении.

Таблица 6.3. Объемы открытых поровых каналов и защемленного воздуха

Влагосодержание свежесфор-

Объем защемленного

Объем открытых

мованного асбестоцемента

воздуха, %

поровых каналов, %

0,23

4,4

9,9

0,26

7,3

7,7

0,29

10,2

5,5

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

187

7. ТЕПЛОМАССОПЕРЕНОС ПРИ ПРОМЕРЗАНИИ МАТЕРИАЛОВ

7.1. Сублимационное обезвоживание материалов

При значениях температуры Т<273 К и относительной влажности внешней парогазовой среды φ<1 происходит сублимация льда. Перенос пара в стефановском режиме в порах материала подобен массопереносу в одиночных капиллярах (раздел 2.2). Толстые пленки влаги в порах не образуются. Зависимость перемещения зоны сублимации имеет вид ξ=βτ1/2, где β – коэффициент, м/с1/2. Структура материала при сублимационном обезвоживании не претерпевает существенных изменений. Кристаллы льда образуются в макропорах. В микропорах из-за малого количества сорбированных молекул воды образование зародышей кристаллов и их рост маловероятны [51, 52].

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

лекул воды снижается, а льда возрастает.

В первую очередь переносятся молекулы воды, имеющие более низкую энергию связи с центрами сорбции материала. Сорбированная вода, как показали диэлектрические исследования (рис. 7.1), не переходит в лед даже при температуре 160 К [35]. Наблюдается только снижение их под-

Рис. 7.1. Зависимость диэлектрической

вижности в микропорах.

проницаемости ε от температуры Т (К)

Рассмотрим

математическую

при частоте 2 105 Гц торфа с влагосо-

модель процесса сублимации.

держанием 0,245 (1), 0,195 (2), 0,05 (3)

188

Пусть ось координат x направлена от поверхности в глубь полубесконечной среды (грунта). За отрезок времени dτ граница фазового перехода перемещается на величину dξ. Количество влаги, сублимируемой за это время из элементарного объема dV =Sdξ,

dm 0u0Sdξ.

(7.1.1)

Поток пара направлен в противоположном направлении оси x. Из уравнений Фика и Клапейрона-Менделеева следует

dm = −

Dμ

 

Sdτ,

(7.1.2)

 

 

 

RT0 dx

 

где Т0 – начальная температура грунта и внешней парогазовой среды, принимается постоянной в течение всего процесса сублимации.

С учетом соотношений (7.1.1) и (7.1.2) на границе при x (τ) получим условие

Dμ

= γ0u0

dξ

.

(7.1.3)

 

 

 

 

RT dx

 

 

 

dτ

 

 

0

 

 

 

 

 

 

Вследствие переноса пара внутри пор изменяется его парциальное давление р. Уравнение изменения р во времени и толщине слоя сублима-

ции запишется следующим образом:

 

 

 

p(x, τ)

= D

2 p(x, τ)

, τ>0, 0

<x (τ)

(7.1.4)

∂τ

x2

и решим его при начальных и граничных условиях:

 

 

 

p(0, τ) =pср,

 

(7.1.5)

 

р(х, 0) =p(ξ(τ), τ) =ps,

(7.1.6)

где pср, ps – соответственно парциальное давление пара в парогазовой среде и насыщенного пара.

В мерзлой зоне при x давление пара р=рs. Для зоны сублимации краевая задача (7.1.4) – (7.1.6) имеет решение

р(x, τ)= p

 

+(p

p

 

)

erf [x (2

Dτ)]

 

ср

 

 

,

(7.1.7)

 

 

s

 

erf z

 

 

 

 

 

 

 

 

 

где безразмерный коэффициент z = β(2 D ). Если воспользоваться дополнительным граничным условием (7.1.3), то можно определить

β =

2

(ps p) Dμexp(z2 )

.

(7.1.8)

π

γ0u0 RT0erf z

 

 

 

 

 

 

189

 

 

 

Сублимация – продолжительный во времени процесс. Экспериментальные значения z <10–3. Полагая exp(–z2) 1 и erf z (2/π1/2)z, получим β≈[2Dμ(ps pср)/(RT0γ0u0)]1/2. Из этого соотношения при φ<1 коэффициент диффузии

 

 

 

 

 

D =β2 RT0γ0u0

[2μps (1−ϕ)].

 

 

 

(7.1.9)

Интенсивность сублимации на поверхности слоя (x = 0) с учетом

(7.1.8)

 

 

 

 

(ps p) Dμ

 

 

 

 

 

 

 

 

 

 

ic = −

Dμ d p

 

γ0u0β

 

 

γ0u0β

 

γ0u0β2

 

 

 

=

 

=

 

 

 

)

 

=

. (7.1.10)

RT

d x

πτRT0 erf z

z

2

2 τ

0

 

 

 

2exp(

 

τ

 

 

 

Значения β, рассчитанные по уравнениям (7.1.8) и (7.1.10), с по-

грешностью менее 1% совпадают.

 

 

 

 

 

 

πD ),

 

 

Полагая, что

erf (x 2 Dτ)= erf (z x ξ)βx (ξ

решение (7.1.7)

запишем в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

р(x, τ) =p+ (ps p)(x/ξ),

 

 

 

(7.1.11)

Согласно изотермам десорбции влаги из торфа, при положительных по Цельсию температурах в диапазоне 0,1 <φ<0,7 наблюдается линейная зависимость влагосодержания от φ. При Т<273 К вода переходит в лед, но количество незамерзшей влаги пропорционально φ=p/ps. Приближенно можно записать распределение влаги в зоне сублимации:

u(x, τ) =uп + (uξ uп)(x/ξ),

(7.1.12)

где uп и uξ – соответственно влагосодержания при x =0 и x (τ). Объектами исследования служили цилиндрические образцы диамет-

ром 3 см и высотой 8,5 см механически диспергированного торфа. Они имели начальные влагосодержания в диапазоне u0 =2,9 – 7 и плотности твердых органоминеральных компонентов в единице объема мерзлого торфа γ0 =126–253 кг/м3 (табл. 7.1). Экспериментально находили перемещение границы сублимации ξ(τ), распределения полей влагосодержаний в идентичных образцах по высоте и во времени.

Таблица 7.1. Параметры процесса сублимации льда из мерзлого торфа

Т , К

u

0

γ

ϕ

, %

p

, Па

u ,

u ,

β

10

6

,

 

5

,

0

 

0,

 

s

 

п

н

 

 

D 10

 

 

кг/кг

кг/м3

 

 

 

 

кг/кг

кг/кг

м/с1/2

м2

 

272

6,97

132

78,5

561

0,19

1,4

5,5

 

1,45

 

272

4,76

175

78,5

561

0,19

1,4

5,3

 

1,23

 

270

7

126

83

474

0,22

1,4

 

5

 

 

1,66

 

270

6,21

148

83

474

0,22

0,8

3,3

 

1,23

 

270

4,63

184

83

474

0,22

0,8

5,5

 

1,93

 

270

2,94

253

83

474

0,22

0,8

3,3

 

0,99

 

267

4,8

174

82,3

367

0,16

0,64

2,8

 

1,16

 

 

 

 

 

 

 

190

 

 

 

 

 

 

 

 

 

На рис. 7.2 приведены результаты опытов по сублимационному обезвоживанию, полученные при атмосферном давлении, скорости воз-

 

душного потока над открытой торцевой

 

поверхностью образца v = 0,35 м/с и от-

 

носительной влажности воздуха парога-

 

зовой среды ϕ= 0,65.

 

Процесс сублимации льда из

 

мерзлого торфа в атмосферу (рис. 7.2а)

 

сопровождается частичным испарением

 

незамерзшей воды uн. Согласно кало-

 

риметрическим опытам, при темпера-

 

турах 272, 270, 267, 253 К ее значения в

Рис. 7.2. Кривые сублимационной

исследуемом торфе соответственно бы-

сушки (1 – 3), время τ (ч) (а); зависи-

ли равны 1,4; 0,8; 0,64; 0,4. При влаго-

мости средней интенсивности субли-

содержании 1,4 в мерзлом торфе имеет-

мационной сушки iс 106 (кг/(м2 с)

ся осмотическая (~ 1) и сорбированная

торфа от влагосодержания u (1– 3)

влага (~ 0,4). Для меньших влагосодер-

при температурах 272 (1, 1); 270 (2,

жаний uн доля осмотической влаги со-

2); 267 (3, 3) (б)

ответственно снижалась.

 

В течение всего времени сушки

(рис. 7.2б, кривые 13) из-за плохого подвода тепла через обезвоженные пористые слои торфа средняя интенсивность сушки iс постепенно уменьшается. С увеличением температуры и скорости потока воздуха до v 2 м/с значения iс возрастают. При v >2 м/с скорость сублимационной сушки изменяется незначительно.

При снижении барометрического давления до P =0,33 Па средняя интенсивность iс увеличивается, а время вакуумной сушки существенно уменьшается. Как и для отдельных капилляров (раздел 2.2), при атмосферном давлении воздух «мешает» диффузии пара из внутренних слоев образца в окружающую парогазовую среду. При снижении внешнего давления диффузия пара возрастает. Вакуум-сублимационное обезвоживание мерзлых образцов торфа позволяет оценить количество незамерзающей воды. Для этого необходимо определить влагосодержание слоя материала, контактирующего с фронтом сублимации.

Согласно рис. 7.3, в исследованном диапазоне значений ϕ экспериментально подтверждается зависимость ξ(τ) =βτ1/2, в которой средний коэффициент β=4,87 10–6 м/с1/2. Значения β для каждого режима сублимационного обезвоживания и исходного влагосодержания материала изменяются в пределах (2,8 – 5,5) 10–6 м/с1/2 (см. табл. 7.1). В таблице приведены коэффициенты диффузии D, рассчитанные по соотношению (7.1.9).

Расчетные зависимости интенсивности обезвоживания ic от времени τ (рис. 7.4) согласуются с результатами опытов (точки).

191

Рис. 7.3. Зависимость глубины сублимационного иссушения торфа ξ (мм) от величины τ1/2 1/2): 1) Т=272 К, u0 =4,8 – 7,

γ0 =130 – 175 кг/м3, ϕ=0,78; 2) Т=270 К, u0 =2,9 – 7, γ0=125–250 кг/м3, ϕ=0,83; 3) Т=267 К, u0 =4,8, γ0 =175 кг/м3, ϕ=0,82

Рис. 7.4. Сопоставление экспериментальных (точки) и расчетных (сплошные линии) значений интенсивности сублимации

iс 106 кг/(м2 с) льда из мерзлого торфа: а – Т=272 К, ϕ=0,79, v =1,25 м/с, u0 =7, γ0

=130

кг/м3

(1), u0 =4,8, γ0 =175 кг/м3

(2); б – Т=270 К, ϕ=0,83, v =1,25 м/с, u0 =7, γ0

=126

кг/м3

(1),

u0 =6,2,

γ0 =148

кг/м3

(2),

u0

=4,6,

γ0 =184 кг/м3 (3); время τ (сут)

Согласно экспериментам (рис. 7.5), в зоне сублимации устанавливается линейное распределение влагосодержания, а на границе раздела торф и парогазовая смесь – равновесное влагосодержание материала uп. Вторая

192

 

 

 

 

характерная точка uξ находится на пе-

 

 

 

 

ресечении кривых влагосодержания в

 

 

 

 

 

 

 

 

зоне обезвоживания и границы субли-

 

 

 

 

 

 

 

 

мации.

 

 

 

 

 

 

 

 

Интенсивность

обезвоживания

 

 

 

 

зависит от скорости диффузии пара в

 

 

 

 

порах торфа через зону сублимации.

 

 

 

 

При температуре Т< 273 К и влагосо-

 

 

 

 

держании uξ имеет место динамиче-

 

 

 

 

ское равновесие между насыщенным

 

 

 

 

паром в порах торфа, льдом и неза-

 

 

 

 

мерзшей водой uн. В зоне обезвожива-

 

 

 

 

ния из-за низкой средней скорости

 

 

 

 

сублимации (~10–9 м/с) в каждой про-

 

 

 

 

слойке материала устанавливается со-

 

 

 

 

ответствующее динамическое равно-

Рис. 7.5. Распределения влагосодер-

весие между льдом, связанной влагой

жаний u

по высоте x (см) мерзлого

и паром.

 

 

 

 

торфа в

процессе

сублимации

льда

В поверхностных слоях образца

Т=272 К,

v =1,25

м/с, ϕ=0,78:

а –

удаляется часть льда и сорбированной

u0 =6,8, γ=130 кг/м3, τ=1107 (1), 2160

(2), 2860 (3) ч; б – u0 =5, γ0 =175 кг/м3,

влаги.

Влагосодержание

достигает

τ=1250 (1), 2300 (2), 3000 (3) ч

 

значения uп ≈0,2

(см.

табл. 7.1). В

 

 

 

 

мерзлой

зоне в

процессе

иссушения

начальное влагосодержание u0 остается постоянным.

7.2.Физико-математическаямодельпереносатепла

ивлаги при промерзании материалов

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

Для прогнозных расчетов тепломассопереноса во влажных материалах необходима простая модель с минимальным числом входящих в нее

193

коэффициентов, которая основывается на экспериментальных результатах

идопущениях [53].

1.Процессы промерзания влажных материалов (в частности, грунтов) продолжительны по времени. Перемещение границы фазового перехода ап-

проксимируется зависимостью ξ=βτ1/2. Скорость ее перемещения мала (dξ/dτ~ 10–7 м/с). В связи этим возмущенные поля температуры Т и влагосодержания u в грунте релаксируют до установившихся распределений Т и u быстрее, чем перемещаются границы интервала. Только в этом случае могут выполняться заданные граничные условия I – IV рода. Процесс промерзания (оттаивания) можно рассматривать как последовательную смену установившихся состояний. При этом в слое материала толщиной y(τ) в момент времени τ>>0 распределения температуры и влагосодержания оказываются асимптотически близкими к этим величинам в слое постоянной толщины y =y(τ). Такой подход позволяет получить простые решения задач с подвижными границами, полагая в них первые производные y(τ) и производные более высоких порядков, равными нулю.

2.При низких температурах фазового перехода воды в лед термопереносом влаги можно пренебречь. Поле температур на 1 – 3 порядка устанавливается быстрее, чем поле влагосодержания, поэтому во влагонасыщенных материалах можно отдельно рассматривать переносы тепла и влаги.

3.Коэффициенты математической модели являются переменными

Рис. 7.6. Схема распределения температуры Т (К) (1) и влагосодержания u (2) в промерзающем грунте

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

С учетом приведенных допущений для расчета переноса тепла и влаги при фазовых превращениях следует воспользоваться готовыми решениями с постоянными размерами тел.

Найдем распределение температур в промерзающем слое грунта конечной толщины l. Принимаем температуры на поверхности грунта Т1 (0, τ) =Тс (рис. 7.6, кривая 1), границе мерзлой и талой зон Т1(ξ, τ) =

=Т2(ξ, τ) =Тф и нижней границе Т2(l, τ), равнойначальнойпостояннойтемпературеТ0.

194