
Литература от Абакумова ч2 / Осмотический массоперенос
.pdf
Рис. 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μ |
|
dр |
Sdτ, |
(7.1.2) |
|
|
||||
|
RT0 dx |
|
где Т0 – начальная температура грунта и внешней парогазовой среды, принимается постоянной в течение всего процесса сублимации.
С учетом соотношений (7.1.1) и (7.1.2) на границе при x =ξ(τ) получим условие
− |
Dμ 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τ)] |
|
|
ср |
cр |
|
|
, |
(7.1.7) |
||||
|
|
s |
|
erf z |
|||||
|
|
|
|
|
|
|
|
|
где безразмерный коэффициент z = β(2
D ). Если воспользоваться дополнительным граничным условием (7.1.3), то можно определить
β = |
2 |
(ps − pcр ) 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 − pcр ) Dμ |
|
|
|
|
|
|
|
|
|
|
ic = − |
Dμ d p |
|
γ0u0β |
|
|
γ0u0β |
|
γ0u0β2 |
|||||||
|
|
|
= |
|
= |
|
|
|
) |
≈ |
|
= |
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, τ) =pcр + (ps −pcр)(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б, кривые 1′– 3′) из-за плохого подвода тепла через обезвоженные пористые слои торфа средняя интенсивность сушки 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