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

книги из ГПНТБ / Шубин Г.С. Физические основы и расчет процессов сушки древесины

.pdf
Скачиваний:
11
Добавлен:
24.10.2023
Размер:
11.89 Mб
Скачать

Из изложенного вытекает, что при данной температуре t суще­ ствует не единственное значение cpi, а диапазон Д<р, в котором влагосодержание древесины изменяться не будет, что облегчает, на­ пример, кондиционирование влагосодержания древесины и, следо­ вательно, применение системы допусков и посадок, а также взаимо­ заменяемых деталей.

 

 

 

1.3.3. РАВНОВЕСНОЕ

ВЛАГОСОДЕРЖАНИЕ ДРЕВЕСИНЫ

 

 

 

 

 

ПРИ ПОВЫШЕННОМ

ДАВЛЕНИИ

Для расчета

ряда

технологических процессов, например сушки

со сбросом

давления,

сушки с пропиткой, сушки

в гидрофобных

жидкостях,

тепловой

обработки,

требуется

знание

равновесного

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

при давлениях, превышающих

атмосферное.

На рис. 1.3.4

крайняя

правая линия — это линия равновесного вла­

госодержания

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

( р „ = 1 бар).

Она

соответствует

линии фмакс при атмосферном давлении из рис. 1.3.1. На эту диа-

Рис. 1.3.6. Принципиальная схема экспери­ ментальной установки для исследования равновесного влагосодержания древесины

при повышенном

давлении:

 

/ — образец; 2 — пружинные

весы; 3 — стеклянный

цилиндр; 4 — испытательная

камера;

5—паропе­

регреватель; 6 — контактный

манометр;

7 — галь­

ванометр; 8 регулятор напряжения

грамму можно, пользуясь зависимостью (1.3.4), нанести аналогич­

ные линии

при

других

давлениях, которые

будут

лежать

правее

изобары р п

=

1

бар.

 

 

 

 

 

 

 

 

Возникает

вопрос,

можно

ли экстраполировать

линии Wv

в об­

ласть повышенного

давления

до пересечения

с соответствующими

линиями ф =

фмако-

В работе

[8] линии Wv

при повышенных

давле­

ниях были

построены

по эмпирическим

уравнениям,

полученным

для р п > 1 бар.

 

 

 

i

 

 

 

 

 

Надежный ответ на вопрос о величинах Wv

при р п > 1

бар может

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

Такое

исследование на

специальной установке (рис. 1.3.6) было проведено

в МЛТИ * [25].

Образец

1

(сосна,

береза)

в виде тонкого среза подвешивался

к кварцевым пружинным весам 2, которые крепились в стеклянном цилиндре S, герметично соединенном с испытательной камерой 4У

* Работа под руководством проф. П. С. Серговского и автора выполнялась инж. Э. Б. Щедриной.

2 L

выполненной из трубы.

Сюда из автоклава АГ-2 по змеевиковому

пароперегревателю 5 подавался пар, параметры которого

регулиро­

 

 

вались

при помощи

контактного

is

• 1X

манометра 6 и гальванометра 7.

 

Образец

выдерживался

до устой­

 

*|* 1 <

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

в процес­

 

 

се сорбции. Изменение массы об­

 

 

разца

(длины пружинки)

фикси­

 

VZ

ровалось

оптически

при

помощи

 

д

горизонтального

микроскопа

 

МГ-2, чтр позволило

обеспечить

 

высокую точность замеров. Опыты

h

X

С

*" "Hi.

Рис. 1.3.7. Равновесное влагосодержание древесины в перегретом паре:

ТОО' 110° 120° 130° П0° 150° 160° 170° 180'

/ — Р п = 1 , 5

бар; 2 — рп=2 бар

Температура t°, с

 

 

 

проводились при давлении р'а =

1,2;

1,4; 1,5;

2 бар в среде насы-

щенного и перегретого пара.

 

 

 

A l l

Для суждения о возможности экстраполяции линий равновес­ ного влагосодержания (при положительном ответе количество весьма сложных опытов может быть сокращено) результаты экс­ периментов сопоставлялись с линиями равновесного влагосодержа­ ния в координатах Wv — t при р = const, построенными на основа­ нии экстраполяции линий Wv. Результаты такого сопоставления приведены на рис. 1.3.7. Опытные точки вполне удовлетворительно укладываются на линии экстраполированных значений равновес­ ного влагосодержания древесины при давлении выше 1 бар. Это позволяет построить обобщенную диаграмму равновесного влагосо­ держания древесины (рис. 1.3.8).

Глава 2

ИСХОДНЫЕ УРАВНЕНИЯ И ХАРАКТЕРИСТИКИ ПЕРЕНОСА ВЛАГИ И ТЕПЛА ПРИ СУШКЕ ДРЕВЕСИНЫ

2.1. ОСНОВНЫЕ ЗАКОНЫ ПЕРЕНОСА ВЛАГИ И ТЕПЛА ВНУТРИ ВЛАЖНОГО ТЕЛА

При изучении переноса влаги и тепла внутри тела могут быть два подхода.

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

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

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

личиной

изменения энтропии можно

оценивать результаты пере­

носа.

 

 

В термодинамике необратимых процессов, к которым относится

и сушка,

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

характеризуется скоростью а

изменения энтропии 5 во времени т [26]:

23-

Скорость изменения энтропии определяется суммой произведе­ ний плотности потока переносимой субстанции на термодинами­ ческую силу X* для каждого из видов переноса:

п

1 = 1

Выражение (2.1.2), являющееся основным соотношением Онзагера, может служить основанием для определения потоков и дви­ жущих сил переноса (энергии, массы), в частности в процессах сушки.

Сопоставление уравнения (2.1.2) с преобразованным термодина­ мическим уравнением Гиббса [27] показывает, что для полной ана­ логии между классической термодинамикой и термодинамикой не­ обратимых процессов уравнение 2.1.2 целесообразно представить в виде

п

^ =

2

^ ,

(2.1.3)

где Хг = ТХ*.

i =

i

 

 

 

 

По Онзагеру [26] одному

из авторов теории

термодинамики

необратимых процессов — основой

этой теории являются принципы

линейности и взаимности.

 

 

 

Первый фиксирует многолетней практикой установленный факт,

что плотность переноса q в каком-либо

процессе пропорциональна

термодинамической силе этого процесса

X:

q=LX,

(2.1.4)

где L — коэффициент пропорциональности.

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

диент температуры может

вызвать,

помимо

основного

эффекта —•

потока тепла, также побочный эффект — поток массы).

 

 

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

ние на поток всех взаимосвязанных явлений. Отсюда

феноменоло­

гические законы принимают вид системы линейных уравнений:

п

 

 

 

 

 

 

qi=^LikXk

( / = 1 , 2,

 

п).

 

(2.1.5)

k=i

 

 

 

 

 

 

Здесь Lih — так называемые кинетические

коэффициенты. Коэф­

фициенты Ьц связаны с действием прямого

эффекта

для

данного

потока, а Ьщ (при i^k)

с налагающими

эффектами

(коэффици­

енты увлечения).

 

 

 

 

 

 

Согласно принципу взаимности

влияние

различных

потоков

друг на друга равнозначно, в связи

с чем перекрестные коэффици­

енты (не связанные с прямым эффектом)

должны быть равны

 

Lik=Lu.

 

 

 

 

(2.1.6)

24

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

Предположим первоначально, следуя Ю. А. Михайлову [28]г что перенос массы осуществляется в жидкой фазе и химические превращения в теле отсутствуют. Тогда на основании уравнения Гиббса [27] и закона сохранения скорость возрастания энтропии

 

 

" = -Я{-^--Ч^-т,

 

(2Л.7)

где qt и

qu-—плотность

потоков энергии и вещества

(влаги) во

влажном

материале;

Т — температура;

— химический

потенциал

переноса

вещества.

 

 

 

Очевидность данной записи можно

показать, анализируя раз­

мерности на примере скорости возрастания энтропии за счет пере­ носа только энергии — первый член правой части уравнения (2.1.7).

Размерность энтропии — единица энергии

на

°К

(например,

дж/°К),

а скорость изменения энтропии в

единице

объема—•

втГК-мг.

Плотность потока энергии qt характеризуется

размер­

ностью вт/м2. Тогда движущая сила переноса

тепла

в виде

X* (ко­

торая включает градиент температуры °К/м) должна иметь размер­ ность °К/м(°К)2, что приводит к записи

J2

Из сопоставления уравнений (2.1.3) и (2.1.7) видно, что термо­

динамической силой переноса энергии является

 

 

Xt=

L y y .

 

 

 

(2-1.8)

а переноса вещества

 

 

 

 

 

^ = - 7

V - f .

 

 

 

(2-1.9)

На основании (2.1.5) с учетом (2.1.8) и (2.1.9) система феноме­

нологических уравнений потоков энергии qt

и

вещества qu

будет

иметь вид:

 

 

 

 

 

qt=-Ltt^—LtUTX7^r-

 

 

 

(2.1.10)

q u = L u t ^ f

- L u u T

V J

^

,

(2.1.11)

где

 

 

 

(2.1.12)

Ltu=Lut.

 

 

 

Отсюда видно взаимное влияние переноса тепла и массы. В связи

с тем что химический потенциал [х в общем

случае является

функ­

цией температуры, влагосодержания и давления, т. е. fx = f (t,

0, р),

от градиента температуры V T

зависит

не

 

только прямая

сила

25

переноса тепла, но и сила переноса массы. В связи с этим более удобно изменить уравнения (2.1.10) и (2.1.11) таким образом, чтобы каждое из них имело независимые силы, одна из которых была бы

связана

только с неравномерностью

температурного

поля,

а дру­

г а я — с

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

потенциала.

Преобразования, выполненные, например, в работе

[28], приво­

дят к системе уравнений

 

 

 

 

 

 

 

 

 

qq

= - L q

g

^ — L q U

^ ;

 

 

(2

 

 

q u

= - L U q

^ — L u v V * ,

 

 

 

(2.1.14)

где в = (лт — изотермический химический потенциал, а

 

 

 

 

 

v H - a H , v < / .

 

 

 

< 2 - U 5 >

В уравнении (2.1.13) величина

qq

характеризует

 

поток

тепла,

связанный с потоком энергии qt,

соотношением

 

 

 

 

 

 

 

qq=qt-iqu>

 

 

 

 

 

.16)

где / — удельная энтальпия системы.

 

 

 

 

между

Уравнения

(2.1.13)

и (2.1.14)

позволяют установить связь

 

 

 

 

 

 

 

 

 

 

(2.1

кинетическими

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

и

истинными

(не

осложненными

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

Так, истинный

коэффициент теплопроводности

получается

Яд =

> массо- и влагопроводности

Ки = Ьци, а

коэффициент

б ' ,

характеризующий

в неизотермических условиях

поток

вещества,

вызванный градиентом температуры [28] \

 

 

 

 

 

8 ' = —

1

и ч

 

 

 

 

 

Т L

ии

 

 

 

 

 

 

 

 

 

 

Тогда выражения (2.1.13)

и (2.1.14)

примут вид:

 

(2.1.17)

 

qq=-\S?T-\uVTS7Q;

 

 

 

 

fc,=

- W e - X t

 

, 8 ' V 7 \

'

(2.1.18)

Если массо- и теплоперенос'осуществляются не в жидкой фазе (как это было принято нами вначале), а в газовой или двухфазном потоке, выражения (2.1.17) и (2.1.18) сохраняют свой вид. При

1 По новой терминологии А. В. Лыкова [5], член —jr- , являющийся коэф­

фициентом перед V Т в законе переноса (2.1.14), называется коэффициентом термовлагопроводности, а отношение этого коэффициента к коэффициенту массопроводности Lw — относительным коэффициентом термовлагопроводности 6' (размерность 1/единица потенциала переноса влаги).

26

этом коэффициенты тепло- и массопроводности становятся эквива­ лентными. Так,

где индексы 1 и 2 относятся соответственно к газовой и жидкой фа­ зам. Уравнение переноса тепла записывается в этом случае так:

qg = -hVT+Iq„,

(2.1.19)

где kq и / — эквивалентные коэффициенты теплопроводности и эн­ тальпия.

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

 

Яр=-^РУР,

 

(2Л.20)

где кр

эквивалентный коэффициент молярного переноса.

будет

Тогда

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

влаги

иметь вид*:

 

 

 

qu=-baVe-WVT-KVp.

'

(2.1.21)

2.2.ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ ПЕРЕНОСА ВЛАГИ

ИТЕПЛА

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

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

Учет источников влаги и тепла может быть выполнен по-раз­ ному [29]. Одним из методов их учета является введение критерия

1 В последующем изложении плотность потоков qt и qv и коэффициенты пе­ реноса Хч и Хи обозначаются соответственно q, q', X, X'.

27

фазового превращения г, под которым понимают отношение изме­ нения влагосодержания за счет фазового превращения dll$ (испа­ рение— конденсация) к общему изменению влагосодержания dU:

 

. = - 2 g f .

(2.2.1)

При е = 0 фазовые

превращения

отсутствуют, влагосодержание

изменяется только за

счет переноса

жидкости. При е = 1

влагосо­

держание изменяется

только за счет

испарения жидкости,

перенос

которой отсутствует.

 

 

 

Количество влаги в единице объема при фазовых превращениях

 

 

dU

 

где 6U

скорость изменения влагосодержания (т-— время); ро —

плотность абсолютно сухого тела.

 

Источник тепла

dU

 

 

 

 

где го — теплота фазового

перехода.

 

Если учесть, что влагосодержание связано с потенциалом пере­

носа влаги

соотношением

 

 

 

 

39 = Щ-,

(2.2.2)

в котором

С" —удельная

влагоемкость, являющаяся

аналогом

удельной теплоемкости и показывающая количество влаги, которое необходимо сообщить телу, чтобы поднять его потенциал на одну единицу (подробнее см. разделы 2.3 и 2.4), система уравнений пе­ реноса тепла и влаги будет иметь вид [31]:

-^=а\/Ч+гйе

(2.2.3)

•=а' V 2 € > + « ' 3 ' V 2 ^

(2.2.4)

При наличии внутреннего избыточного давления р система ус­ ложняется:

£=aV*t+r^-g+-%£-VpVt.,

 

(2.2.5)

64 =a'\72e-^a'b'\74-{-aJb'pv2p;

 

(2.2.6)

dp

' „ 2

дв

 

_ = a

P V Р~г

ST.

(2.2.7)

28

Здесь а — коэффициент

температуропроводности;

С

и С р удель­

ная теплоемкость

тела

и парогазовой

смеси;

К — коэффициент

воздухопроницаемости тела; а'—коэффициент

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

сти молекулярного

и капиллярного

переноса

влаги;

6^ — отноше­

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

влагопереноса А / / А / ;

а' — коэффициент потен-

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

 

 

 

 

 

а'р = - ^

,

 

 

 

(2.2.8)

 

 

ЧРо

 

 

 

 

 

где С'в —удельная

влагоемкость тела по отношению к влажному

воздуху в процессе

молярного движения

парогазовой

смеси.

Использование дифференциальных уравнений переноса воз­ можно лишь при наличии сформулированных краевых условий: на­

чальных, характеризующих распределение

потенциалов

переноса

(например, температуры)

при т = 0, и

граничных,

отображающих

закон

взаимодействия

поверхности

тела

с

окружающей средой.

Для

конвективной сушки в общем случае

применительно

к си­

стеме

уравнений

(2.2.3—2.2.7)

граничные

условия

можно

запи­

сать:

 

 

 

 

 

 

 

 

 

 

 

« Сс -

Q - Ч

V On -

(1 - е) г0ав (9„ - в р ) = 0 ;

 

(2.2.9)

 

« e ( e „ - e p ) + x ' ( v e ) n + X ' 8 ' ( V O n + X p ( v p ) n = 0 ;

(2.2.10)

 

 

 

/ > n = / » = const.

 

 

 

(2.2.11)

В этих выражениях индекс п обозначает поверхность. Первый член граничного условия для переноса тепла (2.2.9) определяет ко­ личество тепла, подведенного к поверхности тела (tc — температура среды; а — коэффициент теплообмена), второй член — количество тепла, отведенного от поверхности внутрь тела, третий член — за­ траты тепла на испарение влаги на поверхности (а' — коэффициент

массообмена между средой и поверхностью, отнесенной к разности

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

В п

и равновесного

по­

тенциала массопереноса в среде 6 Р ) . При е =

1

(испарение

происхо­

дит только

внутри тела) этот член уравнения

(2.2.9)

равен 0,

при

8 = 0 он максимален, так как вся влага в жидком виде

подводится

к поверхности.

 

 

 

 

 

 

 

 

 

 

Первый

член граничного

условия

для

переноса

влаги

(2.2.10)

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

в окружаю­

щую среду,

второй,

третий и четвертый члены — количество влаги,

подводимое

к поверхности соответственно под действием

градиен­

тов потенциала переноса влаги V 9 , температуры Vt

и избыточного

давления

V p .

 

 

 

 

 

 

 

 

 

 

Когда

градиент

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

может

рассматриваться

как

потенциал

 

переноса

влаги — к древесине

это

относится

для

зон,

в которых

 

U<Un.r

(W<Wn.v)

—подробнее

см. раздел

2.4,диф­

ференциальные уравнения переноса и условия

однозначности

(при

V p ' = 0) примут вид [31]:

 

 

 

 

 

 

 

 

 

29

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

sr0

дЦ

 

 

(2.2.12)

 

С

дх

 

 

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

 

 

 

dU

a'42U+a'b

уЧ-

(2.2.13)

дх

 

 

 

граничные условия

 

 

 

а ( * с - tB) -1(

v t)a - (1 - е) г 0 а и Р о (£/„ - Up) = 0;

(2.2.14)

*'VPo{Ua-С/р)

+ а'ро ( V f / ) n

+ a ' p 0 8 ( V 0 n = 0 .

(2.2.15)

Здесь а' —коэффициент влагообмена, отнесенный к разности

влагосодержания на поверхности Uu и равновесного влагосодержания в среде £ / р ; б — относительный коэффициент термодиффузии (размерность IIград), равный отношению коэффициента термодиф­ фузии a't к коэффициенту потенциалопроводности а' при записи

закона переноса влаги в виде

q' = -a'p0 v U - a ' t p 0 v * = — a ' p 0 ( V £ / + 8 Vt).

(2.2.16)

Уравнения (2.2.9—2.2.10) или (2.2.14—2.2.15) являются гранич­ ными условиями третьего рода. При больших значениях коэффици­ ентов внешнего тепло- и массообмена ( а - > « , а^—>-оо) граничные

условия превращаются в их частный случай — граничные условия первого рода.

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

ный в какой-либо зоне тела (хотя

бы в минимальной

дозе), мгно­

венно отражается во всем теле (говорят, что

период

релаксации

теплового %г или влажностного %' напряжения равен нулю).

 

Однако в некоторых явлениях

переноса

период

релаксации

имеет конечную, сравнительно большую величину — например,

при

капиллярном переносе влаги т ' = (0,4-М ,2) Ю - 4

сек [5]. Тогда в

об­

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

ное уравнение при переносе тепла q=XVt

преобразуется:

< 7 = — X Vt — trq,

(2.2.17)

 

где величина хг связана со скоростью переноса щ соотношением

30

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