|
u |
ж |
, |
(5.27) |
m |
|
|||
|
2 |
|||
|
|
|
||
|
|
|
|
где – зазор между пластинами.
Трудность решения поставленной задачи обусловлена сложностью определения касательного напряжения на стенке и, следовательно, динамической скорости. Эта сложность связана с тем, что в гофрированных каналах два источника турбулентности – стенка и гофры. Определить раздельно вклад каждого из них в потери энергии, которые возникают при движении жидкостей в их каналах, пока не представляется возможным. Поэтому при расчѐте гидравлических потерь вычисляют так называемые условные касательные напряжения по уравнению (5.19) с заменой в нѐм значения на у .
В подобных случаях динамическую скорость целесообразнее выражать не через касательное напряжение, а через среднее значение диссипации энергии (удельную мощность).
Согласно уравнению энергии (1.41), диссипация энергии происходит в основном в пристенном слое, где градиенты скоростей максимальны, и описывается уравнением
E |
2 , |
(5.28) |
0 |
|
|
где – скорость сдвига.
С другой стороны, касательные напряжения на стенке связаны со скоростью сдвига зависимостью
|
|
0 |
. |
|
(5.29) |
|
|
|
|
|
|
|
|
Из совместного решения уравнений (5.28) и (5.29) следует |
||||||
|
0 |
( E )0,5. |
|
(5.30) |
||
|
|
0 |
|
|
|
|
Приняв допущение Прандтля о постоянстве касательных на- |
||||||
пряжений, после подстановки |
0 |
в уравнение (2.82) получим извест- |
||||
|
|
|
|
|
|
|
ную зависимость |
|
|
|
|
|
|
|
|
|
0,25 |
|
|
|
|
u |
|
E0 |
. |
(5.31) |
|
|
|
2 |
|
|||
|
|
|
|
|
|
|
ж
151
Допустим, с определенной степенью достоверности, что турбулентность в гофрированных каналах изотропна. В таком случае, согласно работе А.Н. Колмогорова, величину E0 можно выразить через среднее значение диссипации энергии E в виде равенства
E |
4 |
E |
4 |
|
N |
, |
(5.32) |
|||
|
|
|
|
|
||||||
0 |
|
|
|
|
|
|
Vж |
|
|
|
|
|
|
|
|
|
|
|
|
||
где N – мощность источника турбулентности; Vж – объѐм жидкости |
||||||||||
в аппарате. |
|
|
|
|
|
|
|
|
|
|
Подставив выражение (5.32) в уравнение (5.31), запишем |
|
|||||||||
|
|
|
E |
0,25 |
|
(5.33) |
||||
u |
|
|
|
|
|
|
, |
|
||
|
2 |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
||
|
|
|
ж |
|
|
|
|
|
|
|
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
N |
|
|
0,25 |
|
|
||
u |
|
|
|
|
|
|
. |
(5.34) |
||
|
|
|
|
|
|
|
|
|||
|
|
2 V |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
||
|
|
|
ж ж |
|
|
|
|
|
|
|
Мощность источника турбулентности N (мощность потока) при движении жидкости в межпластинном зазоре найдем по уравнению
|
|
|
l |
|
|
ж w3 |
|
|||
N |
pQ |
у f |
|
|
|
|
|
, |
(5.35) |
|
dэ |
2 |
|||||||||
|
|
|
|
|
||||||
где l и f – длина канала и площадь его поперечного сечения. |
|
|||||||||
Так как объѐм жидкости |
в канале Vж f l , то из |
уравне- |
||||||||
ний (5.34) и (5.35) следует |
|
|
|
|
|
|
|
|
|
|
|
|
y w3v |
0,25 |
|
(5.36) |
|||||
|
u |
|
|
|
. |
|
||||
|
2dэ |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
||
В окончательном виде уравнение для расчѐта коэффициентов теплоотдачи в пластинчатых аппаратах, с пластинами сетчатопоточного типа с пересекающимися гофрами «в ѐлку», примет вид
152
0,95 |
|
|
|
|
||||||
|
|
0,36 |
жu*dэ |
|
|
|
Pr0,33 . |
(5.37) |
||
|
|
|||||||||
|
|
|
||||||||
Таким образом, чтобы рассчитать |
|
|
по уравнению (5.37), не- |
|||||||
|
||||||||||
обходимо знать значения коэффициентов |
и |
у |
в равенстве (5.36). |
|||||||
|
|
|
|
|
|
|
|
|
|
|
Их можно определить, имея опытные данные по теплообмену и гидравлическим сопротивлениям или эмпирические уравнения для их расчета для того или иного типа пластин. Значения у , взятые из ра-
боты [27] для пластин сетчато-поточного типа с пересекающимися гофрами «в ѐлку»: ПР 03, ПР 05, ПР 05М и «Розенблат» – приведены
вприл. 5. Пластины подобных типов наиболее часто применяются
всовременных пластинчатых теплообменниках.
Опытные значения критерия Нуссельта рассчитывались по
уравнению типа |
|
|
|
|
|
|
|
|
Nu |
0 |
A Rem2 |
Pr0,43 |
; |
(5.38) |
|||
|
|
2 |
|
|
|
|
|
|
коэффициенты гидравлического сопротивления – по уравнению |
|
|||||||
|
|
|
|
B |
|
, |
|
(5.39) |
|
|
у |
|
|
|
|
||
|
|
|
Re0,25 |
|
||||
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
||
где значения A2 , B и m2 зависят от типа пластин и приведены в [27]. Поиск значений коэффициента осуществлялся в пределах
изменения Pr от 2 до 50 и |
m |
от 80 до 1000. В результате обработки |
|
|
опытных данных было установлено, что в уравнении (5.36) коэффициент = 1,8.
Отклонение опытных значений от вычисленных по уравнению (5.38) для указанных типов пластин не превышает 14 %. Можно предположить, такая же точность будет и при расчете теплообменников, собранных из других пластин сетчато-поточного типа с пересекающимися гофрами. Однако следует иметь в виду, чтобы воспользоваться уравнением (5.37), необходимо предварительно провести гидравлические испытания теплообменника в целях определе-
ния коэффициента |
у |
в равенстве (5.36). |
|
|
153
Выполненные решения позволяют сделать заключение о том, что эмпирические уравнения являются частными случаями уравнений, полученных полуэмпирическим способом на основе аналогии между переносом импульса и теплоты.
Решив сравнительно простые задачи, прейдѐм к решению более сложных задач теплообмена при аэробном культивировании микроорганизмов в аппаратах различных конструкций.
5.2.Теплообмен между газожидкостной смесью
итеплопередающей поверхностью
Вданном подразделе рассматриваются решения нескольких задач, связанных с определением коэффициентов теплообмена, в аппаратах, наиболее часто применяемых для нагрева или охлаждения газожидкостных смесей в микробиологической и других отраслях пищевой промышленности. Отвод теплоты от культуральной среды
вних осуществляется с помощью тепловой рубашки, встроенных или выносных теплообменных устройств.
Наиболее простой вариант представляет решение задачи теплообмена между стенкой аппарата и газожидкостной смесью в барботажных культиваторах колонного типа.
5.2.1.Культиваторы колонного типа
степловой рубашкой
Пустотелые аппараты подобного типа широко используются в различных отраслях микробиологической промышленности, в частности для культивирования дрожжей, производства лимонной кислоты и т. п. Особенностью культиваторов колонного типа является отсутствие принудительной циркуляции жидкости через аппарат. Однако даже если такое и имеет место, например в установках непрерывного культивирования микроорганизмов, приведенная скорость жидкой фазы в них бывает настолько мала, что она не оказывает заметного влияния на скорость переноса импульса. Иными словами, источником турбулентности будет являться только движение газа относительно культуральной среды, т. е. приведѐнная скорость газа. Это, конечно, не означает полного отсутствия циркуляции жидкости внутри аппарата, но скорость еѐ будет определяться всѐ той же при-
154
ведѐнной скоростью газа и может быть рассчитана по уравнению, предложенному автором работы [33]:
wж 4(vг gdп )1/ 3. |
(5.40) |
Вклад касательных напряжений, обусловленных движением жидкости вдоль теплопередающей поверхности, в развитие турбулентности, как показывают расчѐты, настолько мал по сравнению с долей, приходящейся на движение газа относительно жидкой фазы, что им можно пренебречь.
Влияние на скорость циркуляции вязкости культуральной среды учитывается еѐ влиянием на размер пузырей. Согласно работе [7], диаметр пузырей и динамическая вязкость связаны пропорциональностью dп ~ 0,25 . В таком случае wж ~ 0,083 , т. е. влияние вязкости
на скорость циркуляции невелика.
При решении рассматриваемой задачи будем считать, что воздух барботируется через культуральную жидкость, представляющую собой суспензию, в которой твѐрдой фазой являются клетки микроорганизмов, а жидкой (в большинстве случаев) – водные растворы различных питательных веществ и продуктов метаболизма. Физические свойства таких суспензий берутся либо из таблиц справочной литературы, либо рассчитываются по имеющимся в литературе формулам [34].
Для решения интеграла (5.15) необходимо найти значение динамической скорости. Использование для этой цели уравнения (2.68) представляет определенную сложность, так как источниками турбулентности служат относительная скорость фаз и стенка. Влиянием стенки в данной задаче можно пренебречь, потому что касательные напряжения, обусловленные скоростью циркуляции суспензии, невелики. Вследствие чего воспользуемся тем же приѐмом, что и в предыдущей задаче.
Входящую в уравнение (5.34) величину N найдем, исходя из следующих соображений. Мощность источника турбулентности бу-
дет |
определяться сопротивлением |
газожидкостного |
слоя |
высо- |
той H гж при движении газовой фазы относительно жидкой фазы. Так |
||||
как |
в данном случае приведѐнная |
скорость жидкости |
wж |
0 , то |
мощность источника турбулентности будет определяться только приведѐнной скоростью газа и может быть выражена зависимостью
155