|
|
|
|
|
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