Ar3
С3 3m .
Подставляя значение C3 в равенство (5.77), запишем
|
dux |
|
|
|
|
A |
r2 |
|
|
|
A |
|
rm |
. |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
dr |
|
|
|
|
3 |
|
|
|
|
|
|
3 |
|
r |
|
||||
После интегрирования данного выражения получим |
|
|||||||||||||||||||
u |
|
A |
r |
3 |
|
A |
r |
3 |
ln r C . |
(5.78) |
||||||||||
|
9 |
|
|
|
|
3 |
|
|
||||||||||||
|
x |
|
|
|
|
|
|
|
m |
|
|
4 |
|
|||||||
Постоянную интегрирования |
|
|
C4 |
|
находим из условия |
r R, |
||||||||||||||
ux 0 : |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Ar 3lnR |
|
|
|
AR3 |
|
||||||||||
С4 |
|
|
|
|
|
m |
|
|
|
|
|
|
|
|
|
|
. |
|
||
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
9 |
|
|
|
||||
Подставив значение C4 в равенство (5.78), запишем уравнение профиля скоростей в нисходящем потоке:
u |
|
A |
r |
3 |
R |
3 |
A |
r |
3 |
ln |
R |
. |
(5.79) |
||
|
|
|
|
|
|
|
|
|
|||||||
|
xn |
9 |
|
|
|
|
|
3 m |
|
r |
|
|
|||
Уравнения (5.76) и (5.79) описывают профили локальных скоростей от оси аппарата до стенки. В них требуется найти значения радиусов r0 и rm. Нетрудно доказать, что
r R r0 . (5.80)
m |
2 |
|
Радиус r0 можно найти из условия равенства объѐмных расхо-
дов жидкости в восходящем и нисходящем потоках.
Расход жидкости через любое сечение определяется уравнением (13). После подстановки в него уравнения (5.76) или (5.79) получим уравнения для расчѐта объѐмного расхода жидкости в восходящем и нисходящем потоках:
176
|
|
|
r0 |
|
|
|
|
A r0 |
|
|
|
|
|
Q |
2 u |
xv |
rdr |
2 |
|
r3 |
r3 rdr; |
|
|
(5.81) |
|||
|
|
|
|||||||||||
v |
|
|
|
|
|
|
9 |
0 |
|
|
|
|
|
|
|
|
0 |
|
|
|
|
0 |
|
|
|
|
|
Q 2 |
R u |
|
rdr |
|
2 |
A r0 |
r3 R3 |
r3 ln |
R |
rdr . |
(5.82) |
||
xn |
|
|
|
|
|||||||||
n |
|
|
|
|
9 0 |
|
|
m |
r |
|
|
||
|
r |
|
|
|
|
|
|
|
|
|
|||
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
Для определения радиуса r0 воспользуемся системой MathCad. Для наглядности рассмотрим конкретный пример движения жидко-
сти |
в |
условиях |
естественной |
конвекции в аппарате радиу- |
||
сом R |
1,5 м при разности температур T 2,5 К. Жидкость облада- |
|||||
ет |
следующими |
физическими |
свойствами: |
0,1 10 3 град 1, |
||
|
1000 кг/м3, |
1,3 10 3 Па∙с. В расчѐтах используем следующую |
||||
программу: |
|
|
|
|
||
|
|
|
|
r0 |
0,0001 |
|
|
|
|
|
Given |
|
|
|
|
|
0 |
Qv r0 Qn r0 |
|
|
|
|
|
r0 |
MinErr r0 |
|
|
|
|
|
|
r0 |
0,901, |
|
где r0 = 0,0001 – начальная произвольно заданная величина. Результаты расчѐтов, представленные на рис. 5.5 в виде графи-
ческого изображения профиля скорости по радиусу аппарата, позволяют сделать определѐнные выводы.
Во-первых, средние значения скоростей движения жидкости в восходящем и нисходящем потоках различны ( wв 3,33 м/с; wн 2,076 м/с), и они настолько высоки, что течение не может быть ламинарным. Ламинарный режим возможен при очень малых значениях T .
Во-вторых, при турбулентном режиме движения жидкости, как уже было сказано в подразд. 2.5, профиль скорости по сечению потока отличается от профиля при ламинарном течении. Турбулентные потоки, согласно двухслойной модели Прандтля, делятся на две зоны, в которых профили скоростей описываются системой двух уравнений (2.76). Поэтому температура по сечению потока не может изменяться по линейному закону. Линейно она меняется только в при-
177
стенном слое. В турбулентном ядре она изменяется мало, так же как и скорость.
В-третьих, при турбулентном режиме движения не будет чѐткой границы между восходящим и нисходящим потоками, как это предполагается при ламинарном. Здесь мы имеем дело с так называемой свободной турбулентностью, в условиях которой будет постоянно происходить обмен турбулентными вихрями между противоположно направленными потоками. Желающие более подробно ознакомиться с решением задач процессов переноса при свободной турбулентности могут обратиться к работе [2, с. 673–700].
ux, м/с
R
rm
r0
|
|
r |
Рис. 5.5. Изменение локальных скоростей по радиусу аппарата |
||
Линии соответствуют: |
уравнению (5.76); |
уравнению (5.79) |
Указанные обстоятельства приводят к тому, что радиус r0 , устанавливающий границу между потоками, будет постоянно меняться. Найти можно лишь его некое осреднѐнное значение r0 . Решение за-
дачи по определению радиуса сильно усложняется. Мы не будем пытаться еѐ сейчас решать, а наметим лишь путь решения.
Так же как и при ламинарном режиме движения, r0 находится из условия равенства расходов Qв и Qн . Однако входящие в уравнения (5.81) и (5.82) значения локальных скоростей uxv и uxn следует
178
находить из уравнений (2.76). Но здесь возникает другая сложность, связанная с определением динамической скорости u . В данном случае можно воспользоваться уравнением (5.33).
Динамическая скорость является мерой интенсивности турбулентных пульсаций [2, стр. 543] и зависит от источника их возникновения. В рассматриваемом нами аппарате источников турбулентности два: касательные напряжения на стенке и относительная скорость движения восходящего и нисходящего потоков. Но интенсивность этих источников (при заданной температуре теплообменной поверхности) определяется разностью температур T . Последняя, в свою очередь, зависит от тепловой энергии Ek , выделяемой дрожжевыми
клетками в процессе сбраживания сусла и рассчитываемой либо по уравнению (3.26), либо по (3.27). Для пивных дрожжей численный коэффициент в уравнении (3.27) должен быть иной. Таким образом, аналогично уравнению (5.47), суммарную диссипацию энергии можно выразить равенством
|
Ek |
Eo E , |
(5.83) |
где E |
– диссипация энергии, обусловленная касательным напряже- |
||
нием |
на стенке аппарата, |
которое рассчитывается по |
уравне- |
нию (1.41а); Eo – диссипация энергии, обусловленная относительным движением потоков.
Значение Ek находят экспериментально. Зная Ek и E , можно
определить и E .
Исходя из указанных соображений, динамическую скорость вычислим по уравнению (5.33), заменив в нѐм E на Ek . В итоге получим
|
|
|
0.,25 |
|
|
u |
χ |
μEk |
. |
(5.84) |
|
ρж2 |
|||||
|
|
|
|
5.4.2. Теплообмен
Аналитическое решение задачи по теплообмену между стенкой аппарата и жидкостью при естественной конвекции значительно затруднено, поэтому в литературе приводятся уравнения для расчѐта
179
коэффициентов |
теплообмена |
в |
условиях естественной конвекции |
||||
в критериальном виде [16, с. 308]: |
|
|
|
|
|||
|
|
Nu |
C(Gr Pr)n, |
(5.85) |
|||
где C и n – |
|
коэффициенты, |
|
определяемые |
экспериментально |
||
(табл. 5.1). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 5.1 |
|
|
|
|
|
|
|
|
Режим движения |
|
Gr Pr |
|
|
C |
|
n |
Ламинарный |
|
500 |
|
|
1,18 |
|
0,125 |
Переходный |
|
500–2∙107 |
|
|
0,54 |
|
0,25 |
Турбулентный |
|
Более 2 107 |
|
|
0,135 |
|
0,33 |
В уравнении (5.85) Gr g R2
T / 2 – критерий Грасгофа. Он связан с критерием Рейнольдса зависимостью
Gr B Re,
где коэффициент B зависит от формы поперечного сечения канала. При турбулентном режиме движения коэффициент теплооб-
мена между стенкой аппарата и жидкостью можно определить и по уравнениям (5.4), (5.18) и (5.84), из которых следует, что в условиях естественной конвекции уравнение для расчѐта коэффициента теплообмена принимает вид
|
|
|
|
|
0,25 |
0,95 |
|
|
|
|
|
|
Ek |
|
|
|
|
|
|
|
ж |
|
da |
|
||
|
|
|
|
|
|
|||
0,18 |
т |
|
|
ж |
|
|
Pr0,33. |
(5.86) |
dа |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
Из уравнений (5.84) и (5.86) наглядно видна связь динамической скорости u , определяющей гидродинамические условия в аппарате, коэффициента теплообмена с удельной скоростью прироста биомассы
и со скоростью выделения дрожжевыми клетками теплоты Ek в условиях естественной конвекции.
180