CC BY
80
Научный журнал КубГАУ, №103(09), 2014 года
1
УДК 664.727
ИССЛЕДОВАНИЕ ЭФФЕКТА ДАВЛЕНИЯ В ПРОЦЕССАХ ТЕПЛОМАССОПЕРЕНОСА С ТВЕРДОЙ ФАЗОЙ ПРИ СУШКЕ
Подгорный Сергей Александрович к.т.н
Косачев Вячеслав Степанович д.т.н., профессор
Кошевой Евгений Пантелеевич д.т.н., профессор
Кубанский государственный технологический университет, Краснодар, Россия
В статье рассмотрены вопросы математического описания эффекта давления в процессах тепломассопереноса с твердой фазой при сушке
Ключевые слова: ПОТЕНЦИАЛЫ ПЕРЕНОСА, МЕТОД КОНЕЧНЫХ ЭЛЕМЕНТОВ, МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ, ТЕПЛОМАССОПЕРЕНОС
UDC 664.727
STUDYING THE EFFECT OF PRESSURE IN THE PROCESSES OF HEAT MASS TRANSFER AT SOLID PHASE DURING DRYING
Podgornyy Sergey Alexandrovich Cand.Tech.Sci.
Kosachev Vyacheslav Stepanovich Dr.Sci.Tech., professor
Koshevoy Evgeniy Panteleevich Dr.Sci.Tech., professor
Kuban State Technological University, Krasnodar, Russia
The article studies the issues of mathematical description of pressure effect in the processes of heat mass transfer with a solid phase during drying
Keywords: TRANSFER POTENTIALS, FINITE ELEMENTS METHOD, MATHEMATICAL MODELING, HEAT MASS TRANSFER
Теория тепло-массопереноса [1] создана на базе термодинамики необратимых процессов и позволяет изучать переноса энергии (тепла), массы вещества и давления. Интенсивный и быстрый процесс сушки приводит к росту давления в капиллярно-пористых телах, при этом градиент давления в пределах пористой структуры твердого тела становится важным и должен быть включен в формулировку системы уравнений тепло-массопереноса Лыкова [2—5]. В технологии сушки важным является развитие поля давления, что в ряде случаев оказывает определяющее влияние на качество высушиваемого материала. Рассмотрим систему уравнений Лыкова-Михайлова [1] переноса температуры, влаги и давления к капиллярно - пористой среде:
d
—t(x, X) = aq 2
dx dx
J 2 cm J
t(x, X ) + є* r-----0(x, X )
cq dx
d d2 d2 d2
—0(x, x) = am-2 0(x, X) + am-50-2t(x, x ) + am*5p-2P(x x )
dX dx2 dx2 dx2
d
—p(x, x ;>=ap 2 dX dx2
d2 , , cm d 0( ,
p(x, x) - є-------0(x, x)
cq dX
Где t, 0, p - соответственно потенциалы переноса тепла, массы и давления; ap - коэффициент потенциалопроводности фильтрационного движения парогазовой смеси ap=k/cbg; k- коэффициент воздухопроницаемости; cb -коэффициент характеризует емкость капиллярно-пористого тела по
http://ej.kubagro.ru/2014/09/pdf/58.pdf
Научный журнал КубГАУ, №103(09), 2014 года
2
тношению к влажному воздуху в процессе молярного движения парогазовой смеси; g - плотность тела; e - коэффициент фазового превращения; г - удельная теплота испарения; 5 - термоградиентный коэффициент.
Тепловые характеристики а и c отмечены индексом q (aq°a, cq°c). Влажностные характеристики а и c отмечены индексом m (am°a, cm°c). Удачное решение системы уравнений предложено [6] путем переформулировки уравнений Лыкова-Михайлова, с введением дополнительных параметров, регулирующих градиенты давления, влажности и температуры. В данной работе изложена проверка соответствия этих модифицированных уравнений фундаментальным соотношениям Ларса Онзагера [1], основанным на многокомпонентной задаче переноса потенциалов неравновесной термодинамики с одной стороны и определении связи параметров уравнений Лыкова-Михайлова с параметрами, предложенными в статье [6]. Полученная уточненная система уравнений использовалась для сравнительной оценки процесса сушки капиллярно-пористого тела в обычных условиях и в вакууме. Приведем эту систему уравнений к каноническому виду относительно производных по времени, заменяя их в правой части уравнений.
Сначала заменим производную в первом уравнении. Затем объединим производные температуры по координате в результате первое уравнение системы уравнений Лыкова-Михайлова примет вид:
(aq'cq +
5е) d2
iL«* t) = v“q c^am cm ^ -dlt(x, t) + !m:fml! лі ^ t) + am'°m're'8p .4**, x)
dt cq dx cq dx cq dx
q dx q dx q
Выразим производную по времени в третьем уравнении системы уравнений Лыкова-Михайлова.
Объединим производные давления по координате в результате третье уравнение системы уравнений Лыкова-Михайлова примет вид:
—p(x,t) = dt
am'cm'e'8e d2 ^ amcme d2 (ap'cq amcme'8p) d2
. -d^t(x,t) . Ae(x,t) +
dx2 cq dx2
---2p(x, t)
dx2
c
c
q
q
Таким образом, система уравнений Лыкова-Михайлова в матричном виде принимает вид:
http://ej.kubagro.ru/2014/09/pdf/58.pdf
Научный журнал КубГАУ, №103(09), 2014 года
3
(
f 10 0 0 1 0 10 0 1)
f d ^ — t(x,t) dt
— 0(x,t) dt =
d . . —p(x,t) V dt )
[+ am' cm' r'e ' 80) am' cm' r' e am' cm' r' e ' 8p
cq cq cq
am 80 am 3.
am' cm' e ' 80 am' cm' e (ap ' cq — am' cm' e'
cq cq cq
d2
-%t(x,T) dx2
d2
0(x,t)
dx2
d
2
—2P(x,t)
V dx
Согласно принципу симметрии кинетических коэффициентов, послужившему основой феноменологической термодинамики неравновесных процессов Ларса Онзагера, матрица в правой части уравнения должна быть симметричной. Для этого вводим произвольные множители, преобразующие матрицу коэффициентов переноса системы уравнений Лыкова-Михайлова к Онзагеровскому виду, упрощая и приводя подобные:
f 10 0 ^
0 Cq 0
0 0 C
P )
(aq ' cq + am' cm r' e' 80) am' cm' r' e am' cm' r' e'8p am' cm' r' e ' 80 am' cm r' e am' cm' r ' e'8p
c c c aq1 cq c c
q q q q q q
am'80 am am'8p simplify ® am' C0 '80 am C0 am'8p ' C0
am' cm e'80 am' cm' e (ap ' cq — am' cm' e' 8p) Cp' am' cm' e'80 Cp' am' cm' e Cp ' (ap ' cq — am' cm' e ' 8p)
c c c c c c
L q q q J Lq q qJ
Определяем первый множитель (Ce) из условия симметрии:
= am’Cq' dq solve ,Cq ®
am cm'r'e
cm'r'e
q
cq ' 80
Подставляем это множитель в матрицу коэффициентов переноса:
-/ s \ s -| am cm re ' d0 am cmr'e
(aq' cq + amcm r' e ' 8e) am cmr e amcmr e' 8p
1 0 0
cm r e
0 ----- 0
cq 80
0 0 Cp
V p)
cq
am80
am 8p
am cm e ' 80 am cm e (ap ' cq ~ am cm e ' 8p)
simplify ®
cq 80
am cmre'8p cq' 80
Cp'amcme' 80 Cp'am cm e Cp' (ap'cq am cm e'8p)
Определяем второй множитель (Cp) из условия симметрии:
am'cm' r ' e'8p
cq ' 80
Cp'am'cm'e r' 8p
------------solve , Cp ® —
c
q
80
Подставляем этот множитель в матрицу коэффициентов переноса:
am' cm ' r' e'80 am' cm' r' e ar
f 1 0 0 t
cq' 80
r ' 8.
00
(aq' cq + am' cm' r e' 8b) am 'cm' r e am'cm'r' e' 8p
am' cm 'e' 80
m
^ ' cm' e
am 8p
(a
'8p)
simplify ®
q
a c r e
q
a c r e
am' cm ' r ' e'8p
cq cq ' 80 cq ' 80
am ' cm ' r' e' 8p am ' cm ' r' e' 8p r' 8p ' (ap ' cq — am' cm' e'8p)
q
cq 80
cq 80
c
q
q
q
c
c
q
q
am cm r e am cm r e
a
m
c
q
c
c
c
q
q
q
c
c
c
q
q
q
aq +
c
q
c„. ' r e
c
c
c
m
q
q
q
0
0
am 80
8
0
c
c
c
q
q
q
http://ej.kubagro.ru/2014/09/pdf/58.pdf
Научный журнал КубГАУ, №103(09), 2014 года
4
Проверяем полученную матрицу на симметричность:
am' cm ' r'E' 50 am ' cm ' r'E
q q
am' ^,'Г' E a ■ С 'Г E
cq d0
am ' cm ' r' E' Sp
am ' cm ' r' E' Sp cq ' S0
am' cm'r' E' Sp am' cm ' r ' E' Sp r8p'(ap ' cq ~ am' cm' E'Sp)
cq S0
cq S0
am' cm' r'E' S0 am'cm'r' E
qq a c r E a c r E
cq S0
am ' cm ' r ' E' Sp
am ' cm ' r ' E' Sp cq ' S0
am' cm' r' E' Sp am ' cm ' r ' E' Sp r'Sp'(ap ' cq ~ am' cm' E'
Sp)
cq S0
cq S0
simplify ®
( 0 0 0 1
0 0 0 v 0 0 0 j
Таким образом, при умножении левой и правой части системы уравнений Лыкова-Михайлова в матричном виде на Онзагеровские множители:
(1 0 0 1
cq ' S0
r SL
00
(10 0 1
0 1 0 V 0 0 1 j
(L л
dt d
t(x,t) 0(x,t) p(x,t)
simplify ®
-t(x,t)
dt
cm' r E' A 0(x,t)
dt
cq S0
r' Sp '—p(x't)
dt
S0
(1 0 0 1
0 cm' r' E 0
cq ' S0
0 0 r'Sp
—
V S0 J
(aq'cq + am' cm' r' e'8b) am 'cm' r' E am' cm' r' E' Sp
q
am S0
q
am Sp
am' cm' E (ap ' cq ~ am' cm' E' Sp)
simplify ®
am cm r E S0 am cm r E
am' cm'r' E am' cm'r' E am ' cm' r' E'Sp
cq cq S0 cq S0
am' cm ' r ' E' Sp am' cm ' r ' E' Sp r' Sp ' (ap ' cq - am' cm' E'
cq cq ' S0 cq' S0
Получаем систему уравнений удобную для решения задачи переноса трех компонентов в рамках термодинамики неравновесных процессов
dt
-t(x,t) =
aq +
am cm ' r E' 80 1 d2
a
q J
dx
t(x,t) +
a^■ r- E
V q J
dx
-0(x,t) +
( am cm 'r'E' 8p 1
c
qJ
dx
,p(x,t)
( cm r eA
cq 'dq
v q 0 j
( r 8„ 1
dt
-0(x,t)
( am cm ' r E 1 d2
V q J
t(x,t) +
dx
cq S0
0(x,t) +
dx
( am' cm' r 'e Sp 1 cq ' 80
dx
p(x,t)
p
S0
V ° /
■—p(x,t) = dt
(a c r E S 1 2
m m p
c
qJ
dx
t(x,t) +
( am cmre Sp 1 cq ' 80
0(x,t) -
dx
r Sp (ap cq - am cm E Sp)
cq S0
p(x,t)
dx
Аналогичная система уравнений для численного решения этой задачи была представлена в виде:
Po cq' 8 1 d
V cm I dt
---t(x,t)
(kq + E' 1 km)'
dx
t(x,t) +
V
8 1 d 2
e ' 1 km---I 2 0 (x,t) +
cm J dx
8 1 d 2
E'l' kp---I---2p(x,t )
cm J dx
, ( 8 1 ,2 ,2 ,2
(e' 1po cm) 0(x,t) = e'1 km--------------I----2t(x,t) + (e'1 km)------20(x,t) + (e l'kp)-------2p(x>t)
dt I cm J dx2 dx2 P dx2
^ 1rocpkp I d
dt
p(x,t) =
S 1 d2 d2
E ' 1 kp--I---2t(x,t) + (E' 1 kp )--2 0 (x,t) -
cm J dx dx
1(1 - e)'
2p(x>t)
dx
T r
aq +
aq +
c
c
q
q
c
c
q
q
c
q
d
c~, 'r'E
m
0
0
S
0
a„ '^. 'r'E
mm
aq +
c
c
c
q
q
q
c
q
a
m
am' cm ' E 'S0
c
c
c
q
q
q
2
2
d
d
d
2
2
a c r E
d
d
d
mm
2
2
d
d
8
2
d
c
m
2
k
2
d
p
k
k
m
m
http://ej.kubagro.ru/2014/09/pdf/58.pdf
Научный журнал КубГАУ, №103(09), 2014 года
5
Нивелируя эту систему уравнений относительно производной температуры по времени:
ро•cq 51 d -----— 1 — t(x,- )
d-
(kq + e-1- km )•
р о •c q 5
(£-1ро • cm)-0(x, -)
d-
ро•cq 5
e • 1 • k.
р о •c q 5
5 1 ,2 d
d2 t( ) I , k 5 I d2
TTt(x’X ) I £- 1 km — I- —1
dx I cm J dx
--------- + —---
ро•cq 5
m2 cm J dx
ро•cq 5
(x,t) (e-
------ + ---
m2
dx
ро•cq 5
x,t) (e
----- + —
1ро•cp•kp )d , ,
----------- I p(x, - )
dt
ро•cq 5
5 1 d 2 2
Є'!' kp----I---21( x, -) (e-1- kp)---j 0 (x, -)
cm J dx dx
) Ie-1-kp• 5 I d2 ( I 2 p( cm J dx2
' f р о - cq 5 I
cm J
d2 1-kp) 2p(x, -) dx2
( ро•cq 51
1 cm J
i-(1 -e^;^ km d2 2p(x, - ) dx2
ро•cq 5
ро•cq 5
ро•cq 5
Получаем аналог системы уравнений переноса Лыкова - Михайлова:
d-
-t(x,t) =
(kq + Xkm-e)' cq- ро _
dx
t(x, -) +
11kp£ 1 v Чр° j
dx
0(x, - ) +
11 km£ 1
v cqро j
dx
p(x, -)
( 2 1 1cm2e
cq-5
v q J
—0(x, -) d-
( 1-k™-e 1 ,2 d
v Vp° j
( 1-cm -k„ -e 1 ,2
—2t(x, -) + dx2
mp
v Ч5ро J
0(x, - ) +
dx2
' 1cmkm£ 1
v cq •5 ро j
---2p(x,- )
dx2
( 1 cm cp kp 1 v cq • km- 5 J
—p(x,-) d-
( 1-k^ •e 1 ,2 d
p
v Чро J
---2t(x,- ) +
dx2
( 1c™ k„ -e 1 d2
mp
са5ро
vq
---2 0(x,- ) +
dx2
1 cm kp2 (e - !)
VV5^ J
—2p(x, -) dx2
Сравнивая её с полученной нами системой уравнений:
d
•t(x,-)
( amcmre501
d-
( c ■ r- e A
v cq•50 J
aq +
—0(x,-) d-
,2
q J
t(x,-) +
dx
v q J
( Є г Е'І 1 ,2
0(x,-) +
dx
mm p
v cq J
d
dx
p(x,-)
( a • c •r^p'^
v q J
—2t(x,-) + dx2
c-50
v q 0 J
—20(x, -) + dx2
( am' cm'r e ' 5p 1 c_50
v q 0 J
—2 p(x,-) dx2
( 1 v 50 J
—p(x,-) d-
( amcmre-5p 1
c
( e г Е'^ 1 ,2
----2 t(x, - ) +
dx2
m m p
c50
v q 0 J
1 z.
0(x,-) -
dx2
r-5p-(ap cq - am cm e-5p)
cq-50
v q J
Получаем возможность идентификации параметров этих систем относительно друг друга путем приравнивания соответствующих множителей при производных этих систем:
---2 p(x,-)
dx2
5
c
c
m
m
c
c
c
m
m
m
2
d
c
c
c
m
m
m
+
c
c
c
c
m
m
m
m
2
2
2
d
d
d
d
2
2
2
d
mm
2
2
2
a., -c.. -r-е
mm
mm
2
2
http://ej.kubagro.ru/2014/09/pdf/58.pdf
Научный журнал КубГАУ, №103(09), 2014 года
6
Л 2
1cm e cmre cq8 cq89 1 cm cp' kp r 8p cq • km'8 §9
И соответствующих слагаемых правых частей этих систем уравнений с учетом их симметричной структуры относительно главной диагонали (соотношения Онзагера):
(kq + 1 • km' e) cq • Po
- aq +
am' cm' r~ e89
cq
1 • kp 'e am' cm'г e cq • Po cq
1 • cm' kp 'e _ am' cm'r e cq •8 po cq • 89
1 • km' e am' cm'r e ' 8p
cq Po
c
q
1 • cm- km’ e _ f am- cm- r e 8p^
cq •8 po ^ cq • 89 у
1 • cm- kp • (e - !) _ r 8p • (am- cm- e 8p - ap • cq) cq • km' 8 po cq • 89
Полученная система уравнений может быть использована для определения функциональной зависимости параметров этих систем относительно друг друга. Для этого преобразуем полученную систему уравнений к виду, удобному для численного решения, избавившись от множителей при производных по времени в левой части системы уравнений переноса Лыкова - Михайлова, разделив их на множители при соответствующих производных. После приведения подобных получаем соответствующую матрицу кинетических коэффициентов в правой части
http://ej.kubagro.ru/2014/09/pdf/58.pdf
Научный журнал КубГАУ, №103(09), 2014 года
7
5 :q + e1km) cm 5 e-1km cm 5 e-1-V— cm
P ■ c -5 Ko q P ■ c -5 oq P ■ c -5 oq
c c c
m m m
5 e-1-km- — cm e-1km e-1kp
e-1-Po-cm e-1-Pocm e-1-Po-cm
5 e-1kp- — cm e-1-kp kp2 -1(1 -Є)і“ km
1 Po- cp-kp 1pocp-kp 1 - Po- cp-kp
km km km
simplify ®
kq + km-e 1-kme 1-kpe
cqpo cqpo cqpo
km5 km kp
2 cm po cmPo cmPo
km-5-e km-e kp (e - 1)
cm- cppo cppo cppo
Учитывая положения неравновесной термодинамики, о превалировании кинетических коэффициентов главной диагонали матрицы переноса, получаем возможность уточнить их числовые значения при совместном влиянии переносимых компонентов.
Для этого потребуем, чтобы коэффициенты главной диагонали были равны между собой, вводя целевую функцию:
f kq + 1km e km ^ 2 " km kp (e- 1)"
О "O о 1 о 3. "О о V 1 о 3. "О о 1 1 о "О 0 1
Z(d,cm’ Po,kq,km,kp>cp) •
Остальные коэффициенты должны быть меньше них по абсолютной величине, что соответствует следующей системе неравенств:
Po > 1000
2
kp(e-1) f km' e + km'5 e Л
1 о "О о V 1 о "О о 1 о 3. о "О о
> 0
km f km' 5 + kp Л
cm' po V О 3 -"О о cmpo J
> 0
kq +1 ■ km' e f 1 ■ km' e 1 - kp - e Л
О "O о V О "O о О "О о
>0
Используя в качестве начальных приближений наиболее характерные, для зерновых продуктов, значения параметров [3], уточняем их с учетом выполнения данных соотношений и получаем следующий набор данных
http://ej.kubagro.ru/2014/09/pdf/58.pdf
Научный журнал КубГАУ, №103(09), 2014 года
8
( 5 Л
cm
О Q.
kq
km
kp
l
e
cq
v cp J
( - 7 ^
3.053 х 10
1.161 х 10 1.229 х 10 1.375 х 102.851 х 10"
1.157 х 10" 35.313
3
2 х 10 0.056
3
V 4.706 х 10
Используя полученные соотношения, определяем значения параметров, входящих в систему уравнений переноса температуры, влаги и давления к капиллярно - пористой среде системы уравнений переноса Лыкова -Михайлова:
Х'cm e cmre
V5
cq '5Є
X'cm • cp • kp cq• km 5
r 5M
5e
(kq + *•'km e)
c q P o
' - aq +
1kp e am cm r •e
II о Q. cr о c q
X' cm'kp'e am'cm'r'e
c q' 5'P o cq' 5Є
1 km'e am' cm'r 'e '5p
c q' P o cq
X' cm'km 'e ( am'cm r 'e
cq '5'Po V cq 5Є
X' cmkp2 ' (e - k) r' 5 p-(am'cm'£'5p - apcq) cq km' 5'P o c q' 5Є
Таким образом, набору параметров:
5
3
7
7
am' cm 'r 'e 5
e
c
q
p
http://ej.kubagro.ru/2014/09/pdf/58.pdf
Научный журнал КубГАУ, №103(09), 2014 года
9
( 8 > ( - 7 ^ 3.053 x 10
cm 1.161X 10- 5
ro 1.229X 103
kq Цп 1.375X 10- 3 - 7
:= 2.851X 10
kp l -7
1.157X 10 35.313
e cq 2 X 10- 3 0.056
ч cp У ч4.706X 10- 6 ,
соответствует набор параметров системы уравнений переноса температуры, влаги и давления к капиллярно - пористой среде системы Лыкова - Михайлова:
(a > aq ( - 5> 1.99X 10
am 8.10X 10- 6
ap 1.99X 10- 5
8q 0.026
3.4X 10- 5
чr -j ч 35.316 j
Нулевая невязка множителей при производных этих систем показывает их однотипный характер:
http://ej.kubagro.ru/2014/09/pdf/58.pdf
Научный журнал КубГАУ, №103(09), 2014 года
10
(kq + 1 kme) f am cm reS9 'j
cq Pq
aq +
1 cm e cm'r e
cq 5 cq d9
= 0
1 cm'cp ■ kp _ cq ■ km5 5Є
1 kp e am^ cm^r e
cq Pq
0
c
q
0
0
c
q
, ■c ■k„ e a~_ ■ c~ ■ re
mp mm
cq 5 Pq cq 5Є
'■ cm^ kp2(e _ 1) rV
= _0
cq ■ km5Po
cq-69
1 kme am^ cm^ re5p
cq Pq
= 0
1 cm^ kme am^ cm^ re5p
cq 5 Pq
cq 5Є
= 0
a
c
q
= 0
С учетом полученных соотношений было проведено моделирование процессов переноса тепла, влаги и давления из сферического тела для различных начальных и граничных условий методом конечных элементов
[7]. В результате были получены типичные диаграммы изменения температуры, влажности и давления по радиусу сферы во времени, когда влажная (0 = 1) и холодная (t = 0) сфера в которой внутреннее давление равно давлению окружающей среды (p = 1) нагревается этой средой от температуры равной нулю до температуры равной единице (на поверхности сферы) и одновременно сушится за счет удаления влаги с поверхности сферы:
1.0000,9000,8000,700-2_ 0,000S- 0,500-“ 0,40( 0,30002000,100 0,000
1,000 0,800 0,600
11§ Время
0,125 0,300 0,475
0,650 РадИу
0,825 ,000
0,000 0,200 0,400
pip
http://ej.kubagrQ.ru/2014/09/pdf/58.pdf
Научный журнал КубГАУ, №103(09), 2014 года
11
Рисунок 1 - Зависимости изменения температуры, влажности и давления при окружающем давлении.
В дальнейшем было проведено моделирование аналогичного процесса в условиях вакуумного воздействия на сферическую частицу. Типичные диаграммы изменения температуры, влажности и давления по радиусу сферы во времени, когда влажная (0 = 1) холодная (t = 0) сфера в которой внутреннее давление равно атмосферному давлению (р = 1) помещается в вакуум (р = 0) и нагревается от температуры равной нулю до температуры равной единице (на поверхности сферы) и одновременно сушится за счет удаления влаги с поверхности сферы, были также получены методом конечных элементов:
0,096
0,080
0,06
-1,000 0,900 0,800 -0,700 0,600 0,500 ф 0,400 '0,300 0,200 0,100 0,000
0,125 0,300 0,475
0,650 РаЛИу'
0,825 ,000
1,0000,9000,8000,7000,600' 0,500 0,400 0,3000,200 0,100 0,000
0,125
0,300
0,475
Радиус
1,000
Время
Рисунок 2 - Зависимости изменения температуры, влажности и давления под вакуумом.
Сравнительный анализ полученных диаграмм показал, что вакуумирование позволяет интенсифицировать удаление влаги с одновременным предотвращением роста давления во внутренних областях сферической частицы, что в реальных условиях сушки приводит к нежелательным последствиям (растрескиванию семян).
Вывод
При наличии, сравнимых по интенсивности переноса, градиентов давления, влажности и температуры в обычных условиях сушки наблюдается значительный рост давления внутри капиллярно-пористого тела, приводящий к растрескиванию. В тоже время в условиях вакуума такого явления не наблюдается при прочих равных условиях проведения процесса сушки.
Литература
1. Лыков А.В., Михайлов Ю.А. Теория тепло- и массопереноса. Госэнергоиздат, 1963.
http://ej.kubagro.ru/2014/09/pdf/58.pdf
Научный журнал КубГАУ, №103(09), 2014 года
12
2. Lewis R.W., Ferguson W.J. The effect of temperature and total gas pressure on the moisture content in a capillary porous body, International Journal for Numerical Methods in Engineering 29 (1990) 357—369.
3. Irudayaraj J., Haghighi K., Stroshine R.L. Finite element analysis of drying with application to cereal-grains, Journal of Agricultural Engineering Research 53 (1992) 209-229.
4. Irudayaraj J., Wu Y. Analysis and application of Luikov's heat, mass, and pressure transfer model to a capillary porous media, Drying Technology 14 (1996) 803-824.
5. Datta A.K. Porous media approaches to studying simultaneous heat and mass transfer in food processes. I: problem formulations, Journal of Food Engineering 80 (2007) 80-95.
6. Conceicaoa R.S.G., Macedob E.N., Pereira L.B.D., Quaresma J.N.N. Hybrid integral transform solution for the analysis of drying in spherical capillary-porous solids based on Luikov equations with pressure gradient. International Journal of Thermal Sciences 71 (2013) 216—236
7. Kosachev V.S., Koshevoy E.P., Podgorny S.A. Using rounding function in problems of finite-element analysis. Studies in mathematical science. 2012.-V.4.-№2.-pp.17-24.
References
1. Lykov A.V., Mihajlov Ju.A. Teorija teplo- i massoperenosa. Gosjenergoizdat, 1963.
2. Lewis R.W., Ferguson W.J. The effect of temperature and total gas pressure on the moisture content in a capillary porous body, International Journal for Numerical Methods in Engineering 29 (1990) 357—369.
3. Irudayaraj J., Haghighi K., Stroshine R.L. Finite element analysis of drying with application to cereal-grains, Journal of Agricultural Engineering Research 53 (1992) 209-229.
4. Irudayaraj J., Wu Y. Analysis and application of Luikov's heat, mass, and pressure transfer model to a capillary porous media, Drying Technology 14 (1996) 803-824.
5. Datta A.K. Porous media approaches to studying simultaneous heat and mass transfer in food processes. I: problem formulations, Journal of Food Engineering 80 (2007) 80-95.
6. Conceicaoa R.S.G., Macedob E.N., Pereira L.B.D., Quaresma J.N.N. Hybrid integral transform solution for the analysis of drying in spherical capillary-porous solids based on Luikov equations with pressure gradient. International Journal of Thermal Sciences 71 (2013) 216—236
7. Kosachev V.S., Koshevoy E.P., Podgorny S.A. Using rounding function in problems of finite-element analysis. Studies in mathematical science. 2012.-V.4.-№2.-pp.17-24.
http://ej.kubagro.ru/2014/09/pdf/58.pdf
