книги из ГПНТБ / Шубин Г.С. Физические основы и расчет процессов сушки древесины
.pdfИз изложенного вытекает, что при данной температуре 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], приво |
|||||||||
дят к системе уравнений |
|
|
|
|
|
|
|
|||
|
|
= - 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) величина |
характеризует |
|
поток |
тепла, |
||||||
связанный с потоком энергии 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