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

iSswV45XbM

.pdf
Скачиваний:
1
Добавлен:
15.04.2023
Размер:
3.67 Mб
Скачать

Рисунок 9. Изменение температурного напора между теплоносителями по длине поверхности теплообмена

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

 

 

 

 

 

 

t =

t'·e­m·k·F

(3.14)

где m = (

 

 

 

 

 

 

 

для прямоточной схемы движения и

 

 

 

 

 

 

m = (

 

 

 

 

для противоточной.

 

 

 

 

Следовательно, в теплообменниках с противоточной схемой движения t

уменьшается в случае, когда W1 < W2 , т.е. m > 0, а при W1 > W2 , т.е. когда m < 0, температурный напор t увеличивается (см. Рис. 9). Выполняя

усреднение температурного напора по всей поверхности теплообмена получают расчетную формулу среднелогарифмического температурного напора:

= ; (3.15)

или же выражение (1.15) записывают в таком виде

= , (3.16)

20

где – большая разность температур между теплоносителями;

– меньшая разность температур между теплоносителями.

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

поверхности теплообмена (когда отношение

/

≤ 2) средний

температурный напор

между

теплоносителями

определяют, как

среднеарифметическое от

и

:

 

 

=

 

(

 

.

 

(3.17)

 

 

 

Для сложных схем движения теплоносителей определение температурного напора по формулам (3.15 - 3.17) будет не точным, поэтому для учета особенностей теплообмена в этом случае используют экспериментально полученные зависимости между величинами P и R для

определения

поправочного коэффициента

к

рассчитанному

среднелогарифмическому напору

. При этом, под

вспомогательной

величиной P

понимается тепловая

эффективность,

а

R - представляет

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

P =

 

 

,

(3.18)

 

 

R =

 

 

.

(3.19)

 

 

Величина P всегда меньше единицы, а R может быть и меньше, и больше единицы в зависимости от соотношения водяных эквивалентов W1 и W2 . На рисунках Приложений 27 – 32 даны такие зависимости для перекрестного, одновременно прямоточного и противоточного, и многократно перекрестного схем движения.

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

21

Q =

 

,

(3.20)

 

 

 

 

 

 

 

 

 

 

 

 

где l -длина трубки;

λ – теплопроводность материала трубки;

dвн и dнар – внутренний и наружный диаметры трубки;

α1 и α2 – соответственно коэффициенты теплоотдачи внутри и снаружи трубки;

t1 и t2 – соответственно температуры горячего и холодного теплоносителей.

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

k = . (3.21)

В практических расчетах по формуле (3.21) руководствуются следующим правилом:

при α1 > α2 за средний диаметр принимают наружный диаметр;

при α1 < α2 за средний диаметр принимают внутренний диаметр;

при α1 ≈ α2 за средний диаметр принимают dср = (dвн + dнар).

Имеющаяся практика расчета тонкостенных труб, в которых dнар/dвн <2, а определение коэффициента теплопередачи ведется по формуле справедливой для плоской стенки

k =

 

,

(3.22)

 

 

 

 

 

 

 

 

 

 

приводит к дополнительной ошибке порядка 4% и поэтому при современном уровне использования вычислительной техники не является оправданной.

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

22

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

Использование критериальных уравнений предполагает следующий стандартный алгоритм расчета коэффициентов теплоотдачи:

1.В зависимости от вида конвективного теплообмена выбирают конкретное уравнение из справочной литературы.

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

3.Определяют режим течения теплоносителя по критерию Рейнольдса Re для вынужденного движения или Релея Ra при свободном движении и уточняют вид критериального уравнения.

4.Рассчитывают по критериальному уравнению число Нуссельта Nu, являющееся безразмерным коэффициентом теплоотдачи.

5.Вычисляют коэффициент теплоотдачи, используя

функциональную связь α = Nu.

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

- со стороны горячего теплоносителя температура стенки определяется, как разность между температурой горячего теплоносителя и ½ сренелогарифмического температурного напора

:

t1ст = t1 -

 

,

(3.23)

 

-со стороны холодного теплоносителя, учитывая высокую теплопроводность материала стенки, выбирают температуру стенки на 1-3 градуса меньше, чем со стороны горячего теплоносителя:

23

t2ст = t1ст – (1···3)0C.

(3.24)

Теплопроводность материала стенки определяют из справочных данных при средней температуре стенки, равной ½(t1ст+t2ст). Дальнейшее уточнение t1ст и t2ст производят, используя метод последовательного приближения. Для этого по критериальным уравнениям рассчитывают коэффициенты теплоотдачи α1 и α2. Затем вычисляют коэффициент теплопередачи по формуле (3.21) или (3.22) и определяют плотность теплового потока q через поверхность теплообмена:

q = k (t1 – t2).

(3.25)

Зная плотность теплового потока, определяют температуры стенки на горячей и холодной стороне:

t1ст = t1 - , (3.26)

t2ст = t2 +

 

.

(3.27)

 

Процентное расхождение между полученными значениями и заданными в первом приближении оценивают допускаемой величиной, обычно в 1 – 5 %. Если расхождение больше установленной нормы, то расчет повторяют, задаваясь уже полученными из расчетов значениями t1ст и t2ст.

Возможен и другой способ определения температур t1ст и t2ст. Для этого задаются в первом приближении ожидаемыми значениями коэффициентов теплоотдачи α1 и α2, выбирая из справочной литературы диапазоны, реально существующие для случаев нагрева и охлаждения жидкостей, газов или кипения и пленочной конденсации жидкостей. В этом случае теплопроводность материала стенки находят из справочников при температуре, равной среднему значению температур горячего и холодного теплоносителей ½(t1+t2). Учитывая, что разность температур между теплоносителями складывается из суммы разностей температур между горячим теплоносителем и стенкой, разности температур между горячей и холодной стороной стенки и разности между температурой стенки с холодной стороны и холодным теплоносителем, можно составить два уравнения, которые определяют связь между этими температурными

24

перепадами

 

и

соответствующими

им

термическими

сопротивлениями:

 

 

 

 

 

 

 

 

t1 – t2 = (t1 -

t1ст) + (t1ст t2ст) + (t2ст - t2),

(3.28)

 

 

 

 

 

=

 

 

=

 

,

 

(3.29)

 

 

 

 

 

 

 

 

 

где

=

 

 

– термическое сопротивление конвекции с горячей

 

 

стороны;

=– термическое сопротивление теплопроводности стенки;

=– термическое сопротивление конвекции с холодной

стороны;

толщина стенки трубы.

Решая совместно эти два уравнения, находят выражения для температур стенки с горячей стороны t1ст и с холодной t2ст:

t1ст = t1 -

 

 

;

(3.30)

 

 

 

t2ст = t1ст – (t1 - t1ст )

 

.

(3.31)

 

 

 

Далее определяют

коэффициенты теплоотдачи

и

по

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

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

1) при ламинарном вязкостном движении, когда произведение критериев Грасгофа Gr и Прандтля Pr меньше 8·105, среднее значение критерия Нуссельта Nuср определится уравнением [2, С. 212],

Nuср = 1,55· (Pe

 

)0,33· (

 

-0,14·

(3.32)

 

 

в котором l – длина трубки, – поправка на гидродинамический начальный участок, d - определяющий размер, которым может быть либо внутренний диаметр трубы , либо эквивалентный диаметр канала сложной формы dэкв:

25

dэкв=

 

,

(3.33)

 

где f – площадь поперечного сечения канала,

 

– полный периметр канала,

 

либо (для кольцевого канала)

 

dэкв = D – dнар,

(3.34)

где D – внутренний диаметр наружной трубы,

 

либо (для внешнего прохода теплоносителя в межтрубном пространстве кожухотрубного теплообменника)

dэкв = , (3.35)

где – число трубок;

физические свойства жидкости, входящие в критерии уравнения,

выбираются при температуре t = 0,5(tст

+

 

,

значение динамической

 

вязкости жидкости

находится при

средней

температуре стенки, а

коэффициент εl, учитывающий поправку на гидродинамический начальный участок, выбирают из следующей Табл. 1 [2, С. 213].

Таблица 1

Значения εl при ламинарном режиме

l/d

1

2

5

10

15

20

30

40

50

εl

1,90

1,70

1,44

1,28

1,18

1,13

1,05

1,02

1

2) при ламинарном вязкостно-гравитационном режиме движения, когда произведение критериев Грасгофа и Прандтля Gr· Pr>8·105 [6, С. 374]

Nuср = 0,15·Re0,33·Pr0,43· Gr0,1· (Prж/Prст)0,25εl ,

(3.36)

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

3) при турбулентном режиме движении (Re >104) [5, С. 89]

 

Nuср = 0,021·Re0,8·Pr0,43· (Prж/Prст)0,25εl ,

(3.37)

За определяющую температуру в этом уравнении принимают среднюю температуру соответствующего теплоносителя на входе и выходе из

26

теплообменника tср = ( + ), а за определяющий размер – эквивалентный

диаметр dэкв равный учетверенной площади канала, деленной его полный периметр (для круглых труб dэкв равен геометрическому диаметру трубы).

4) при переходном режиме движения, когда 2300 <Re <104 [6, С. 378]

Nuср = K0·Pr0,43· (

 

 

)0,25,

 

 

 

 

 

 

(3.38)

 

 

 

 

 

 

 

 

 

 

где комплекс K0

находят из Табл. 2 [6, С. 377].

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Re·10-3

2,2

2,3

 

2,5

 

3,0

3,5

4,0

5

6

7

8

 

9

10

K0

2,2

3,6

 

4,9

 

7,5

10

12,2

16,5

20

24

27

 

30

33

При расчете теплообмена в изогнутых трубах в вышеприведенных уравнениях полученное число Nuср умножают на поправочный коэффициент εl , который определяется выражением [2, С. 219]

εl = 1+ 1,8

 

,

(3.39)

 

где – радиус изгиба трубы.

5) при движении теплоносителя в пучке труб под углом атаки ψ = 90о

предложены следующие зависимости для расчета среднего значения числа

Nu [1,5, С. 104]:

при коридорной компоновке и Re < 1000

 

Nuср = 0,56Re0,5·Pr0,36·(Prж/Prст)0,25;

(3.40)

при коридорной компоновке и Re > 1000

 

Nuср = 0,22Re0,65·Pr0,36·(Prж/Prст)0,25;

(3.41)

при шахматной компоновке и Re < 1000

 

Nuср = 0,56Re0,5·Pr0,36·(Prж/Prст)0,25;

(3.42)

при шахматной компоновке и Re > 1000

 

Nuср = 0,40Re0,60·Pr0,36·(Prж/Prст)0,25;

(3.43)

6) при пленочной конденсации водяного пара и ламинарном стекании пленки под действием силы тяжести, которое происходит, если dнар < 20·( )0,5, используется формула для расчета среднего значения

коэффициента теплоотдачи, предложенная Нуссельтом [2, С. 272]:

27

αср = a·

 

,

(3.44)

 

где a = 0,943, b = h – т.е. высоте вертикальной стенки или тубы, a = 0,728, b = dнар для горизонтальной трубы,

поверхностное натяжение пленки,

плотность пленки,

ускорение свободного падения,

температура насыщения водяного пара.

6)теплоотдача при развитом пузырьковом кипении воды в условиях свободного движения преимущественно зависит от тепловой нагрузки q и

давления насыщения pн, поэтому в практических расчетах ограничиваются применением эмпирических формул, справедливых в диапазоне давлений от 0,1 до 4 Мпа, и полученных при анализе и обработке экспериментальных данных [2, С. 311]:

αср = 3,0 q0,7 pн0,15,

(3.45)

αср = 38,7 t2,33 pн0,5,

(3.46)

где t соответствует разности температур стенки tст и насыщения водяного пара .

7) теплоотдача при пузырьковом кипении и вынужденной конвекции в трубах должна учитывать гидродинамическое влияние; такое влияние учитывается на практике отношением коэффициента теплоотдачи по формулам развитого кипения, когда отсутствует воздействие скорости потока на теплообмен, т. е. по формулам (3.45, 3.46), и коэффициента теплоотдачи по формулам вынужденного конвективного теплообмена однофазной жидкости, например, по формуле (3.37). Интерполяционная зависимость при этом имеет вид [2, С. 317]:

= . (3.47)

Если это отношение меньше 0,5, то среднее значение коэффициента теплоотдачи αср рассчитывают по формулам (3.45, 3.46), при величине отношения больше 2 среднее значение αср рассчитывают по формуле (3.37), а в случае, когда отношение лежит в диапазоне 0,5 – 2 коэффициент теплоотдачи αср определяют по формуле:

28

αср = αст·

 

.

(3.48)

 

8) теплоотдача при использовании в качестве теплоносителя излучающего газа представляет собой сложный теплообмен и должна учитывать суммарный коэффициент теплоотдачи, что предполагает метод раздельного учета теплообмена конвекцией и излучением; в этом случае суммарное значение среднего коэффициента теплоотдачи будет включать конвективную составляющую αк и лучистую αл :

αср = αк + αл .

(3.49)

Коэффициент αл определяют после расчета плотности лучистого теплового потока qл, который согласно закону Стефана-Больцмана может быть определен по приближенной формуле [5, С. 176]:

qл =

[(

 

)4 - (

 

)4],

(3.50)

 

 

αл =

 

,

 

 

 

 

(3.51)

 

 

 

 

 

где – коэффициент излучения абсолютно черного тела = 5,67 Вт/(м2·K4),

– приведенная степень черноты, =

 

,

 

степень черноты газа,

степень черноты поверхности теплообмена,

среднее значение абсолютной температуры газа, К,

среднее значение абсолютной температуры стенки, К.

Вид уравнения (3.50) предполагает передачу тепла от объема газа к более холодной стенке, поэтому степень черноты нужно выбирать при температуре . Если же поток тепла направлен от стенки к газу, то величину нужно выбирать при температуре . Это позволяет учесть то обстоятельство, что поглощательная способность газа по отношению к излучению от стенки не равна его степени черноты.

В практических технических расчетах наибольший интерес представляет степень черноты газового объема. Для углекислого газа и водяного пара выполнены многочисленные исследования и их результаты для удобства пользования сведены в номограммы (в частности, номограмма для СО2 представлена в Приложении 24). В ней степень черноты СО2 приведена в зависимости от произведения парциального давления газа p и длины пути луча l. Номограмма построена с учетом того, что излучение проходит через единичную площадку из окружающей её газовой полусферы с

29

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]