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

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

 

 

 

 

 

m

 

 

 

 

 

(

)d

 

 

 

 

 

0

 

 

;

 

(5.12)

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

(

)d

 

 

 

 

0

 

 

 

.

(5.13)

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

Подынтегральные функции (

)

и m (

) находим из решения

интегралов (4.37) или (4.38), в зависимости от условий теплообмена или массообмена, совместно с системой уравнений (4.42) (или любыми другими, взятыми из прил. 6). Для трѐхслойной модели турбулентного потока решение интегралов (5.12) и (5.13) выразится суммой трѐх интегралов.

В случае теплообмена решение будет иметь вид

6

0

 

 

d

30

 

 

d

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

.

(5.14)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

(0,124 )4

 

 

1

 

 

 

 

 

1

 

 

 

 

6

 

 

 

1

30

 

 

 

1

 

 

Pr

Pr 5

Pr

2,5

 

 

 

 

 

 

 

 

 

 

 

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

из уравнения (5.14) в уравнение (5.12),

найдѐм среднюю величину безразмерной разности температур:

 

 

 

Pr

 

m

6

d

 

 

 

30

 

 

 

 

d

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

. (5.15)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m 0

0 1 Pr (0,124

)4

 

6 1

 

Pr (0,2

1)

301 Pr (0,4

1)

 

 

 

 

Аналогичным образом можно найти среднее значение безраз-

мерной концентрации

 

 

m из уравнения (5.13), заменив Pr на PrD .

 

 

 

 

 

 

Чтобы определить коэффициенты теплоотдачи из уравне-

ния (5.3),

необходимо

 

предварительно

рассчитать

m

по

уравне-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

нию (5.4) и найти значение

 

 

по уравнению (5.15). Для упрощения

 

решения задачи функцию

Pr/

 

f (

m ,

Pr) в уравнении (5.3) ап-

 

проксимируют зависимостью вида

 

 

 

 

 

 

146

Pr

C mn Prm .

(5.16)

 

Подставив зависимость (5.16) в уравнение (5.3), получим

Nu A 1m n Prm .

(5.17)

Уравнения (5.16) и (5.17) пригодны для расчѐта коэффициентов теплообмена при движении как однофазных, так и многофазных потоков в любых аппаратах, так как они получены на основе общих положений полуэмпирической теории турбулентного переноса.

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

5.1. Теплообмен между потоком однофазной жидкости и гладкой теплопередающей поверхностью

5.1.1. Теплообмен в каналах с различной формой поперечного сечения

Теплообмен в трубах с круглым поперечным сечением. Ре-

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

Обработка экспериментальных данных по теплообмену при движении в трубах однофазных [9, 10] и многофазных потоков позволила установить численные значения коэффициентов пропорциональности и показателей степеней в уравнениях (5.16) и (5.17).

В пределах изменения

m

от 100 до 2000 и

Pr

от 4 до 300

 

 

 

 

можно принять C 0,09; n

0,05 ; m 0,33.

 

 

После подстановки значений C, n и m в уравнение (5.17) по-

лучим

 

 

 

 

Nu

0,18 m0,95 Pr0,33 .

 

(5.18)

 

 

147

 

 

Основным источником турбулентности в трубах является стенка, а величина турбулентных пульсаций будет определяться касательными напряжениями на стенке, которые рассчитываются по уравнению

w8 2 . (5.19)

Значение динамической скорости найдѐм из уравнения (2.68)

путѐм подстановки в него уравнения (5.19):

 

 

 

 

 

 

(5.20)

u w

 

.

8

Коэффициент гидравлического трения

рассчитывается по

известным в гидравлике формулам. Считая трубы гидравлически гладкими, запишем

 

0,316

.

 

(5.21)

 

 

 

 

Re0,25

 

Подставляя значение λ из формулы (5.21) в равенство (5.20),

получим

 

 

 

 

 

 

0,2

 

 

u w

 

.

(5.22)

Re0,125

Из уравнений (5.4) и (5.22) следует

 

m

0,1 Re0,875 .

(5.23)

 

 

 

 

 

Таким образом, мы установили связь между максимальным масштабом турбулентности и критерием Re. Если принять за границу перехода от турбулентного режима течения к ламинарному значе-

ние Reкp

2300, то из уравнения (5.23) получим граничное значение

 

 

mкр

87.

(5.24)

 

 

 

 

 

Критическое значение максимального масштаба турбулентно-

сти

mкр

может служить оценкой отношения сил инерции, обуслов-

 

 

 

 

148

ленных турбулентными пульсациями, и сил вязкости не только в однофазных потоках, но и многофазных. Только при mкр < 87 говорить

следует не о ламинарном течении, а о вязкостно-турбулентном, в котором силы вязкости играют первостепенную роль в формировании профиля скорости по сечению потока.

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

После подстановки уравнения (5.23) в уравнение (5.18) полу-

чим

Nu 0,0208Re0,83 Pr0,33 .

(5.25)

Уравнение (5.25) по виду схоже с уравнением Михеева (4.3). Примерно такого же вида уравнения можно найти в работе [8]. Строго говоря, показатель степени при Pr величина переменная и зависит

от того, в каких пределах изменения Pr и

m

аппроксимируется

 

 

 

 

функция Pr/

 

f (Pr, m ) уравнением вида (5.16).

 

 

 

Следует обратить внимание на отсутствие в уравнении (5.25) параметра, учитывающего направление теплового потока, как это сделано в уравнении (4.3) и др. Это связано с тем, что выбор или расчѐт физических свойств жидкостей производится при температуре стенки. Кроме того, при турбулентном режиме движения среды, с чем мы и будем в дальнейшем иметь дело, отношение средней вязкости к вязкости у стенки близко к единице и ею можно пренебречь

5.1.2.Теплообмен при движении жидкостей

вгофрированных каналах пластинчатых аппаратов

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

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

149

Наиболее полные сведения по таким зависимостям можно найти в работе [27], изданной более тридцати лет тому назад. Однако к настоящему времени различными фирмами выпущены теплообменники, скомпонованные из пластин c иными формами гофр, но ка- ких-либо сведений об уравнениях, необходимых для их расчетов, в литературе нет.

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

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

В своем решении будем опираться на работу [30], в которой обобщены результаты исследований по теплообмену в аппаратах, различающихся по конструкции и принципам действия. Решение задачи основано на уравнении (5.3).

Уравнение (5.3) приводится, применительно к пластинчатым аппаратам, к следующему виду:

 

 

Nu

4 m

 

 

Pr

.

 

 

 

 

 

 

 

 

 

 

(

)

 

 

Приведя данное уравнение к виду (5.19), запишем

 

 

Nu

0,36

0,95m

Pr0,33,

(5.26)

где Nu

 

dэ / – критерий Нуссельта, здесь dэ

2 – эквивалент-

 

ный диаметр канала.

Безразмерное расстояние от стенки до оси канала (максимальный масштаб турбулентности)

150