Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Методическое пособие 769

.pdf
Скачиваний:
23
Добавлен:
30.04.2022
Размер:
8.47 Mб
Скачать

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

виде

Nu 0,26Re0,6 Pr0,37,

(6.32)

которое справедливо при Re 103 2 105 ;Pr 0,7 300.

Таким образом, использование теории подобия позволяет обобщить массив экспериментальных данных и существенно снизить их требуемый объем.

Контрольные вопросы

1.Каковы требования к подобию физических явлений?

2.Что утверждает первая теорема подобия?

3.Какие числа подобия можно получить из уравнения сохранения импульса и энергии при естественных масштабах?

4.Назовитенеобходимыеидостаточныеусловияподобияфизическихявлений?

5.Каков физический смысл гидродинамических чисел подобия?

6.Что выражают числа подобия в процессах теплообмена?

7.Какие числа подобия принимают в качестве определяющих в процессах теплообмена?

8.Каков физический смысл критериев Нуссельта и Рэлея? Для анализа каких процессов они используются?

9.Почему теория подобия необходима при моделировании?

70

7. КОНВЕКТИВНЫЙ ОДНОФАЗНЫЙ ТЕПЛОПЕРЕНОС

7.1. Теплообмен при наружном обтекании плоской пластины

Рассмотрим конвективный однофазный теплоперенос при движении теплоносителя около плоской пластины. Самой простейшей является внешняя стационарная задача обтекания пластины потоком несжимаемой жидкости. Поток

жидкости направлен вдоль пластины и имеет постоянную скорость w0, температура t0 постоянна во времени, внутренние источники тепловыделения отсут-

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

Рис. 7.1. Развитие гидродинамических и температурных слоев на плоской пластине

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

Теплообмен при ламинарном обтекании плоской пластины является классической модельной задачей. Математически данная задача может быть описана системой уравнений

 

w

 

 

wy

0;

 

 

 

 

x

 

 

 

 

 

y

 

 

 

 

x

 

 

 

 

 

 

w

w

w

 

 

w

 

2w

;

(7.2)

x

 

 

x

x

x

x

 

y

y

 

y2

 

 

w

w

 

 

 

w

 

2w

 

(7.3)

z w

 

z a

z

 

x

x

 

 

y

y

 

y2

 

 

c соответствующими граничными условиями:

71

 

 

 

 

 

 

 

 

 

y 0;

wx

wy 0;

wz 0;

(7.4)

 

 

 

 

 

 

 

 

 

y ;

wx

w ;

wz wz ,

(7.5)

где

 

t tc

 

;

 

 

t tc

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Точное аналитическое решение такой задачи было получено Блазиусом, а впоследствии и Польгаузеном [10]. В результате получено выражение для определения профиля скорости в пограничном слое, коэффициента гидравлического сопротивления и толщины пограничного слоя. Изменение скорости на границе пограничного слоя относительно основного потока не превышает 1%.

C f 2 0,332Re1x

2;

 

(7.6)

5xRe 1 2

5 x w

1 2 ,

(7.7)

x

 

 

 

 

где Rex w x .

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

ентов а. Последнее отношение получило названия безразмерного критерия Прандтля Pr a cp .

В общем случае толщина температурного пограничного слоя может быть выражена как f Pr .

Для обтекания плоской пластины ламинарным потоком жидкости толщина температурного пограничного слоя определится как

5x Rex1 2 Pr 1 3 .

(7.8)

В свою очередь,

 

 

 

 

1 2

 

 

 

 

 

x Re1x

2 Pr1 3

2 3 w

x 1 2 1 6c1p

3.

(7.9)

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

Nux x 0,332Re1x

2 Pr1 3,

(7.10)

а при граничных условиях второго рода qc const

Nux 0,47Re1x

2 Pr1 3 .

(7.11)

72

Данные зависимости справедливы для ламинарного течения потока жид-

кости Rex Rex.кр 3 105,

Pr 0,6 15.

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

Осредненный коэффициент теплоотдачи на участке длиной l может быть определен как

l 0

 

 

1 l

xdx 2 x l ,

(7.12)

а безразмерный коэффициент теплоотдачи при обтекании пластины при граничных условиях первого рода

Nu

l 0,664Re1x

2 Pr1 3 .

(7.13)

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

Для турбулентного режима течения коэффициент трения может быть определен по формуле Шлихтинга:

C f 2 0,0295Rex

0,2,

(7.14)

а локальный безразмерный коэффициент теплоотдачи

Nux 0,296Re0,8x Pr0,43 .

Выражение справедливо для Rex 3 105. Установлено, что коэффициент теплоотдачи сильнее зависит от массовой скорости при турбулентном режиме

73

течения w 0,8, в то время как ослабевает при ламинарном режиме течения w 0,5 . В то же время влияние молекулярного переноса слабее –

0,57(при ламинарном течении ), еще более слабое падение коэффициента теплоотдачи наблюдается по длине пластины – х 0,2 (при ламинарном

течении х 0,5 ).

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

 

Таблица 7.1

Диапазоны значений числа Прандтля для различных сред

Теплоноситель

Диапазон числа Прандтля

Щелочные металлы

0,8·10-2…10-2

Сплавы металлов, ртуть

10-2… 10-1

Смеси газов

0,2…0,8

Газы

0,7…2,0

Спирты, легкие углеводороды, хладоны, вода

2,0…10

Высокомолекулярные спирты

2,0…90

Масла

90…103

Так, например, для идеального газа a

, Pr что1, связано с тем, что у

упругих твердых сфер возможно лишь поступательное движение молекул. Для реальных газов Pr 0,7 1, однако при больших давлениях и вблизи термоди-

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

ностью, поэтому для них характерны малые числа Прандтля.

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

При теплообмене в условиях обтекания пластины ламинарным потоком возможны два предельных случая.

Случай 1. Pr 1,что означает что коэффициент кинематической вязкости

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

Nux 0,339Re1x

2 Pr1 3 .

(7.16)

74

Случай 2. Pr 1,что означает, что теплообмен происходит в условиях

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

Nux 0,564Pe1x

2.

(7.17)

Схема формирования профилей температуры и скорости при критических значениях числа Прандтля представлены на рис. 7.2.

а

б

Рис. 7.2. Принципиальная схема формирования гидродинамического и температурного слоев: а –Pr 1;б –Pr 1

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

7.2.Теплообмен при движении потока внутри каналов различной формы

иобтекании их снаружи

При исследовании теплообмена внутри каналов или обтекании поверхностей рассматривают два типа задач. Первая из них – это внешние задачи. Это классическое обтекание одиночной трубы или пучка труб потоком жидкости или газа. К внутренним задачам относятся движение внутри труб и каналов различного поперечного сечения. На характер течения оказывает влияние шероховатость поверхности.

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

Если обтекание плоской пластины носит более именее прогнозируемый характер, то в случае обтекания цилиндра оно является крайне неустойчивым

(рис. 7.3).

75

а

б

 

в г

Рис. 7.3. Качественное влияние числа Рейнольдса на гидродинамическую картину течения внешнего обтекания цилиндра [10]:

а –Re 1;б – Re 60 300;

в – Re 103 2 105; г

Re 2 105

При значении числа Рейнол дса Re w d 1 набл дается безотрыв-

ное течение и полная симметрия линий тока. При R 1 60 наблюдается фор-

мирование стаци

нарных вихрей. У величение скорости потока, а следовательно

и числа Рейноль

са, приводит к потере устойчивости. При этом наблюдаются

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

ихревые

образования увеличиваются по ме е удаления от поверхности цилинд а.

При значении Рейнольдса

Re 300 происходит возникновение турбу-

лентнос и в следе. При Re 1000

отрыв происходит при отр 2, а при зна-

чении Re 5 103 устанавливается устой ивое состояние с углом

отрыва

отр 82 . На положение точки отрыва оказывает влияние степень турбулентности набегающего потока. Гидра лическо сопротивление при обтекании ци-

линдра связано с наличие

зоны отрыва, так называемое сопротивление формы.

При это м величи на отр в

рьирует я в пред лах 80 88.

76

Закономерности теплообмен в зависимости от числа ейнольдса и угла отрыва отр представлены на рис. 7 4.

Рис. 7. 4. Зависимость коэффициента теплоотдач при обте ании одиночного цилиндра при изменении числа Рейнольдса [10]: 1 – Re 60; 2 – 300 Re 1000; 3 –103 Re 2 105 ; 4 –Re Reкр 2 105; 5 –Re 2 1 5

Эмпирические расчетные формулы носят общий вид:

 

Nu C Rem Prn k ,

(7.18)

T

 

где Nu D ; коэффициенты C, m, n определяются по табл. 7.2, а значение kT - по табл. 7.3. При проведении расчетов свойств жидкости выбираются по средней температуре. Для учета переменно ти свойств А. А. Жукаускасом [11]

был предложен поправочный коэф фициент kT

(табл. 7.3).

Таблица 7.2

 

 

 

 

 

 

Значения коэффициентов для критериального у авнения теплооб

ена

Re

C

 

 

m

 

n

 

64

1,34

 

 

0,32

 

0,

1

64 – 103

0,52

 

 

0,5

 

0,

7

103 – 2·105

0,26

 

 

0,6

 

0,

7

>2·105

0 ,023

 

 

0,8

 

0,4

Таблица 7.3 Значения коэффициента kT для критериального уравнения теплообмена

Ка пельные жидкости

Газ ы

kT Pr

Prc

0,25 при tc t

 

Pr

Prc

 

kT

0,2 при tc t

kT t

tc 0,

 

25 при tc t

kT 1 при tc t

77

Эмпирическая формула (7.18) относится к обтеканию цилиндра под прямым углом. При обтекании под углом, отличным от 2, вводится поправка:

 

 

 

 

при

2.

(7.19)

 

 

1 0,54cos2

 

Наиболее распространенным в инженерной практике является не обтекание одиночного цилиндра, а движение в пучках труб. Поперечное обтекание пучков труб – широко распространенный случай конвективного теплообмена. В случае одновременного обтекания цилиндров в потоке по характеру течения различают шахматные (staggered tube bundle)и коридорные (in-line tube bundle)

пучки труб (рис. 7.5).

а

б

Рис. 7.5. Шахматный (а) и коридорный (б) пучки труб

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

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

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

 

m

 

n

S

n1

 

Nu C1 Re

 

Pr

 

1

 

kT ,

(7.20)

 

S2

 

 

 

 

 

 

 

а эмпирические зависимости представлены в табл. 7.4 [10].

78

Критериальные зависимости для расчета конвективного

Таблица 7.4

 

 

 

 

 

 

теплообмена в пучках труб

 

 

 

 

 

 

 

Re

Коридорный пучок

Шахматный пучок

 

 

 

 

Nu 0,9Re0,4 Pr0,36 kT

Nu 1,04Re0,4 Pr0,36 k

 

 

 

Re 102

 

 

Re 40

 

 

T

 

103

 

 

 

 

 

 

 

Nu 0,52Re0,5 Pr0,36 k

Nu 0,71Re0,5 Pr0,36 k

 

 

 

 

 

 

 

 

T

 

 

 

 

 

T

 

 

Re 103

 

40 Re 103

 

 

 

 

 

102

 

 

 

 

 

 

 

 

 

 

 

Nu 0,35Re0,6 Pr0,36 S

S

 

0,2 k

 

 

 

 

 

 

S1

S2 2

1

2

T

 

103 Re 2 105

Nu 0,27Re0,63 Pr0,36 kT

 

 

 

 

 

 

Nu 0,4Re0,6 Pr0,36 k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S1

S2 2

 

T

 

 

 

 

 

 

 

 

 

 

 

 

Re 2 105

Nu 0,033Re0,8 Pr0,4 kT

Nu 0,031Re0,8 Pr0,4 S

S

2

0,2 k

 

 

 

 

 

 

 

 

1

 

T

 

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

Теплофизические физические свойства жидкостей выбираются по средней температуре в пограничном слое и среднемассовой температуре по длине пучка от ряда к ряду. Это может быть сделано путем осреднения температуры (при незначительных значениях удельного теплового потока) или разбивке труб на участке по длине с вычислением среднего коэффициента теплоотдачи на каждом из участков.

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

Как и в случае одиночного цилиндра, в эмпирическую зависимость вводится поправка на переменность свойств жидкости (табл. 7.5) [10].

Таблица 7.5

Поправочный коэффициент на свойства жидкости

 

 

 

 

 

Капельные жидкости

 

kT

Pr

Prc 0,25

Газы

kT

t

tc 0,25

при tc tж

 

kT

1 при tc t

 

 

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

79