Материал: Тишин ВБ Новоселов АГ Процессы переноса в технолог аппаратах

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

ента k1 равно 6,87 10 6 . Данная величина получена при условии, что обработка опытных данных производилась в системе СИ. Следует помнить, что коэффициент k1 имеет размерность, так как он, в отли-

чие от безразмерного коэффициента в уравнении (2.130), получен чисто эмпирически.

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

Уравнение (5.66) может быть использовано при вычислении коэффициентов массоотдачи в жидкой фазе; при растворении газов в жидкостях, близких по свойствам к воде; в аппаратах колонного типа в пределах изменения приведенной скорости газа от 0,005

до 0,1 м/с.

5.4. Гидродинамика и теплообмен в условиях естественной конвекции

5.4.1. Гидродинамика

Движение жидкости при естественной конвекции обусловлено изменением еѐ плотности с изменением температуры. Такие течения возникают в аппаратах пищевых и микробиологических производств, в которых протекают процессы, сопровождающиеся выделением теплоты тем или иным способом: механическим, химическим и биологическим, например при анаэробном культивировании микроорганизмов, в процессах брожения и т. п.

Рассмотрим для примера вертикальный цилиндроконический бродильный аппарат (ЦКТ), представленный на рис. 5.4. Аппарат предназначен для сбраживания пивного сусла. В процессе брожения дрожжи потребляют субстрат, размножаются, выделяя продукты метаболизма, в частности диоксид углерода и теплоту, которую приходится постоянно отводить. Для этой цели аппарат снабжѐн распределѐнными по высоте несколькими теплообменными секциями (рубашками), в которые подаѐтся хладоноситель. Высота аппарата может

171

готово-

достигать 10 м. При таком уровне жидкости абсолютное давление

внижней части аппарата будет примерно в два раза больше, чем

вверхней части.

x

1

R

2

r

rm

ux

r0

Рис. 5.4. Схема цилиндроконического бродильного аппарата: 1 – бродильный аппарат; 2 – теплообменные секции

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

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

Рассмотрим вариант, когда температура у стенки аппарата ниже, чем в центре. В этом случае у стенки жидкость будет двигаться вниз, а в центре – вверх (см. рис. 5.4). В таком случае на некотором расстоянии r = r0 от оси x скорость ux = 0.

172

Регулирование скорости движения по высоте можно осуществлять подачей хладоносителя в теплообменные секции.

Для решения задачи теплообмена между теплопередающей поверхностью и жидкостью необходимо знать распределение скорости движения жидкости по сечению аппарата.

Решение различных задач, связанных с движением жидкости и газов вдоль плоских вертикальных нагретых поверхностей в условиях естественной конвекции, можно найти в работах [16, с. 305; 2, с. 303]. Различие будет заключаться в том, что в нашем случае решение будет производиться в цилиндрических координатах.

В своѐм решении поставленной задачи будем следовать работе [16, с. 305], приняв движение жидкости в аппарате ламинарным, а распределение температуры – линейным, т. е.

 

 

 

 

 

 

 

 

 

 

1

 

 

r

,

(5.67)

 

 

 

 

 

 

 

T T

T

 

 

 

 

 

 

 

2

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T1 T2

 

 

 

 

 

 

 

 

 

где

T и

T

– текущие и средние значения температуры;

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

T1

T2 , здесь T1

и T2 – температура в центре и на стенке аппара-

та соответственно; r

– текущий радиус;

R – радиус аппарата.

 

 

Сразу же отметим, что принятие условий ламинарного течения

и линейности профиля температуры по живому сечению аппарата вызвано лишь стремлением упростить решение задачи. Следует иметь в виду, что температура в центре аппарата будет определяться мощностью теплового источника, находящегося внутри потоков жидкости. Этим источником являются дрожжевые клетки, более или менее равномерно распределѐнные в объѐме жидкости [см. уравнение (3.27)]. Поэтому температурный профиль в аппарате следовало бы искать из уравнения (3.7). Задача сильно усложнится. Поскольку нас, прежде всего, интересует путь решения поставленной задачи, а именно, влияние физических свойств среды на скорость движения жидкости в аппарате, ограничимся принятыми условиями.

Рассмотрим установившееся движение жидкости вдоль оси x под действием массовой силы, сил трения и давления. Считая движение жидкости осесимметричным, можно использовать для решения поставленной задачи несколько видоизменѐнное уравнение (2.22а):

173

1 d

r

dux

 

dp

g.

(5.68)

 

 

 

 

 

 

r dr

dr

 

dx

В уравнении (5.68) плотность жидкости является функцией температуры и вычисляется по уравнению

 

 

 

 

(1 (T T )) ,

(5.69)

где – среднее значение плотности жидкости в аппарате; – коэф-

фициент объѐмного расширения жидкости.

Изменение давления по оси x происходит в результате изменения гидростатического давления, т. е

 

 

 

 

 

 

 

dp

 

g.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С учѐтом последнего выражения и равенства (5.69) уравне-

ние (5.68) примет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 d

 

dux

 

 

 

 

 

 

 

 

 

 

(5.70)

 

r

 

 

 

g

 

 

 

g(1 (T T )).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r dr

dr

 

 

Определив из уравнения (5.67) разность температур и подставив еѐ в уравнение (5.70), после несложных преобразований получим

d r

dux

 

 

 

 

g

 

T

 

r

2

dr.

(5.71)

 

 

 

dr

 

 

2R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Представим уравнение (5.71) в следующем виде:

 

d r

 

dux

 

 

 

Ar

2

dr,

 

(5.72)

 

 

dr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

g

T

.

 

 

 

(5.73)

 

 

 

 

 

 

 

 

 

 

 

2R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

174

Для нахождения профиля распределения локальных скоростей движения жидкости по радиусу аппарата необходимо проинтегрировать равенство (5.72) дважды.

После первого интегрирования получим

r

dux

 

A

r3 C .

(5.74)

 

 

 

dr

3

1

 

 

 

 

Знак минус перед производной означает, что положительному приращению радиуса соответствует отрицательное приращение ско-

рости. Постоянную интегрирования C1

находим из следующего ус-

ловия: при r 0,

ux um

const производная dux / dr

0 и C1 0 .

Второе интегрирование производим при следующем гранич-

ном условии: r

r0 , ux

0 :

 

 

 

 

 

 

 

 

 

 

u

 

 

A

r

3

C .

(5.75)

 

 

 

x

 

 

 

 

 

 

 

9

 

 

2

 

 

 

 

 

 

 

 

 

 

Постоянную интегрирования C2 находим из условий: r = r0,

ux = 0, С2

Ar3 / 9.

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

Подставив значение C2

 

в уравнение (5.75), получим уравне-

ние, описывающее профиль скорости по сечению аппарата в восходящей части потока:

u

 

A

r3

r3 .

(5.76)

xy

 

 

9

0

 

 

 

 

 

 

 

Для нахождения уравнения, описывающего профиль скорости в нисходящей части потока жидкости, произведѐм двойное интегрирование (5.72). После первого интегрирования получим

 

r

d ux

 

A

r3

C .

(5.77)

 

 

 

 

 

d r

3

 

3

 

 

 

 

 

 

Постоянную интегрирования

C3

находим из такого условия:

при r rm ( rm

расстояние от оси аппарата до середины нисходяще-

го потока), ux

um производная dux / dr

0:

 

 

175