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

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

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

где Fo

a 2 – число урье, представл ющее собой безразмерное время;

Bi

– число Био, представляющее собой отно ение термического сопро-

тивления стенки к термическому сопротивлени ю переда и тепла с поверх-

ности стенки к окружающ ей среде; ni – собственные числа задачи, о ределяе-

мые как n Bi ct g n.

 

 

 

 

 

ni 2i 1 2.

Для чисел

Bi n1

 

2,

n2 3

2, n3 5 2,...,

Для чисел

Bi 0n1

0,

n2

, n3

2 ,...,

n i i 1 .

В большинстве практических задач требуется определен ие температуры в характерных точ ах тела (либо на поверхности x , либо в средне плоско-

сти x ). В этих случая безразмерная координата

x x принимает значе-

ние 1 или 0. При значениях числа Фурье o 0,3

ряд становится бы стросхо-

дящимся и с достаточной точностью может быть заменен первым членом ряда:

 

 

2sinn1

cos 1x exp n12Fo .

(3.8)

n1

sinn1cosn1

 

 

 

Таким образом, при заданных координатах безразмерн

й перепад темпе-

ратур является функцией чисел Фурье и Био. Решение уравнения представлено в виде графиков в логарифмических координатах (рис. 3.2, 3.3).

Рис. 3.2. Зависимость безразмерного перепада темпера ур от чисел Фурье и Био для поверхности пластины [8]

40

Рис. 3.3. Зависимость безразмерного перепада темпера ур от чисел Фурье и Био для середины пластины [8]

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

- при заданной продолжительности нагрева и интенсивности теплоотдачи определять перепад температур на поверхностях и в ередине пластин ;

-

при заданных пер падах температу и интенсивност теплоотдачи оп-

ределять продолж ительность нагрева;

-

ри задан ных времени нагрева и температурах опр делять интенсив-

ность теплоотдачи с поверхности.

Уравнение соответствует регулярному теплов му режиму, а поле температур остается подобным себе во все последующие моменты времени. Такие

процесс ы называют автомодельны и во времени.

 

0

Количество теплоты, поступающее в

ластину со всех сторон за время от

до , равно изменению энтальпии пластины за этот период

темпера-

тура во всех точках пластины достигает температуры жидкости)

 

 

Q 2 f c tж

t0 ,

(3.9)

 

а для произвольного промежутка времени от 0 до 1

 

 

Q 2 f c t1 t0 2 f c tж t0 1

 

Q 1

 

,

(3.10)

 

 

 

где t0– температура пластины в начальный момент времени; f – площадь боковой пов рхности стенки; t1 – средняя температура по толщине пластины в мо-

41

мент времени 1; t1 tж t0 tж – средний безразмерный перепад темпе-

ратур в момент времени 1.

Характер изменения температурного поля в пластине при различных числах Био представлен на рис. 3.4.

а

б

Рис. 3.4. Изменение температурного поля в пластине в зависимости от числа Био: а –Bi 100; б –Bi 0,1

3.2. Нестационарное температурное поле бесконечно длинного цилиндра

Рассмотрим задачу по определению температурного поля в неограниченном цилиндре радиусом r0 с начальной температурой t0. Введем следующие допущения:

-рассматривается нестационарный процесс;

-внутренние источники тепловыделения отсутствуют;

-теплофизические свойства стенки изотропны по всем направлениям, а значения параметров равны средним значениям;

-длина (или высота) цилиндра существенно превышает ее толщину.

Для удобства вычислений расчет будем вести от температуры окружающей среды tж t . Дифференциальное уравнение теплопроводности для ци-

линдрической системы координат примет вид

а

 

2

1

 

,

0,

0

r r ,

(3.11)

 

 

r

2

r

r

 

 

 

 

0

 

 

 

 

 

 

 

 

 

а граничные условия запишутся в виде

42

 

 

 

 

 

 

 

0 tж t0

при

0;

 

(3.12)

 

 

 

 

 

 

 

 

 

 

0

при r r0;

 

(3.13)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0;

 

при r 0.

 

(3.14)

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Точное решение получено методом разделения переменных и может быть

записано в виде

 

 

 

 

 

 

 

2J1 ni

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

J0 nir

r0 exp ni

Fo ,

(3.15)

 

n

J 2

n

J 2

n

0

 

i

0

 

 

 

i

 

1

i

 

 

 

 

 

 

 

 

 

 

где J0, J1 – функции Бесселя;

 

J0 ni

J1 ni ni

Bi ;Bi r0

.

 

 

Количество теплоты, поступающее в цилиндр за время от 0

до ,

равно изменению энтальпии цилиндра за этот период (температура во всех точ-

ках цилиндра достигает температуры жидкости)

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

 

r 2l c t

ж

t

0

,

 

 

(3.16)

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

а для произвольного промежутка времени от 0 до 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q Q 1

 

,

 

 

 

 

(3.17)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

r

 

2 rdr 21

 

 

 

 

 

 

где

 

tж t tж t0

 

 

0

rdr – средний безразмерный пе-

 

 

r

2

 

 

 

 

 

 

 

0

 

0

 

 

 

 

 

 

0

 

 

 

 

 

 

репад температур в момент времени 1,

t

– средняя температура цилиндра;

r rr0.

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

 

 

 

 

 

 

 

 

4J12 ni

 

 

 

 

2

 

 

 

 

4Bi2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

exp

 

 

 

 

 

 

exp

 

2

 

2

 

 

2

 

 

i

 

2

 

2

 

2

 

i

 

n

 

J

 

n

J

 

 

n

 

 

n Fo

 

i 1 n

n

 

Bi

 

 

 

n Fo . (3.18)

 

 

 

i 1

 

0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

i

 

i

 

 

 

 

 

 

 

 

i

i

 

 

 

 

 

 

 

При значениях числа Фурье F o 0,3

ряд становится быстросходящимся

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4Bi2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

exp

n1 Fo .

 

 

 

 

 

 

(3.19)

 

 

 

 

 

 

 

 

 

 

n12 n12 Bi2

 

 

 

 

 

 

Решение уравнения представлено в виде графиков в логарифмических координатах (рис. 3.5, 3.6).

43

Рис. 3.5. Зависимость безразмерного перепада темпера ур от чисел Фурье и Био на оси цилиндра [8]

Рис. 3.6. Зависимость безразмерного перепада темпера ур от чисел Фурье и Био на поверхности цилиндра [8]

44

3.3. Нестационарная теплопроводность тел конечных размеров

Рассмотрим температурное поле прямоугольного бруса, состоящего из однородного изотропного материала. Необходимо найти распределение температуры в брусе для любого момента времени, а также среднюю температуру для определения подведенной (отведенной) теплоты (рис. 3.7).

Рис. 3.7. Физическая область для задачи нестационарной теплопроводности прямоугольного бруса

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

 

 

2

 

2

 

2

,

 

a

x

2

y

2

z

2

 

 

 

 

 

 

 

 

 

аначальные и граничные условия следующим образом:

0 tж t0 при 0; 0при ; для поверхности при 0

 

 

при x x;

 

 

 

x

 

 

 

при y y ;

 

 

 

y

 

z при z z ,

(3.20)

(3.21)

(3.22)

(3.23)

(3.24)

45

а для центра бруса при

0

 

 

 

0

при x 0;0 y y ;0 z z ;

(3.25)

 

x

 

 

 

 

 

0

при y 0;0 x x;0 z z ;

(3.26)

 

y

 

 

 

 

 

0

при z 0;0 x x ;0 y y .

(3.27)

 

z

 

 

 

Как показано в [8], безразмерный перепад бруса равен произведению безразмерных перепадов температур на пластинах. Рассмотренный метод известен под названием теоремы о перемножении решений:

 

 

 

x y z ,

(3.28)

где х

t x, tж

; y

t y, tж

; z

 

t z, tж

.

 

t0 tж

 

 

 

 

 

 

 

 

 

t0 tж

 

 

 

 

t0 tж

 

Средняя безразмерная температура бруса может быть определена как

 

 

 

 

 

 

 

x

 

y

 

z ,

(3.29)

 

 

 

 

 

 

 

где х t t0 x tжtж ; y t t0 y tжtж ; z t t0 z tжtж .

Величины, входящие в уравнения, могут быть определены по соотношениям, изложенным выше, а также с помощью графиков на рис. 3.2, 3.3, 3.5, 3.6.

При расчета цилиндра конечной длины (рис. 3.8) рассматривают температурное поле бесконечного цилиндра диаметром 2r0 и плоской пластины тол-

щиной 2 z (длина цилиндра принимается равной l 2 z ).

Рис. 3.8. Физическая область для задачи нестационарной теплопроводности цилиндра конечных размеров

46

Безразмерную температуру для такого цилиндра можно записать как

 

 

 

z r ,

(3.30)

где z

t z, tж

; r

t r, tж

,

 

 

 

 

 

t0 tж

 

 

 

 

 

t0 tж

 

а среднюю температуру для любого момента времени как

 

 

 

 

 

 

 

 

z

 

r ,

(3.31)

 

 

 

 

 

 

где z t t0 z tжtж ; r t t0 r tжtж .

Величины, входящие в уравнения, могут быть определены по соотношениям изложенным выше, а также с помощью графиков на рис. 3.2, 3.3, 3.5, 3.6.

При этом в качестве характерных размеров выбирается радиус цилиндра r0 и половина высоты цилиндра z .

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

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

1.Что является целью решения нестационарной задачи теплопроводности?

2.Каков физический смысл критериев Био и Фурье?

3.Как влиет значение числа Био на характер распределения температуры в плоской пластине?

4.Какие теплофизические расчеты можно проводить с помощью графиков на рис. 3.2, 3.3, 3.5, 3.6?

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

6.Каков алгоритм решения задачи нестационарной теплопроводности для цилиндра конечных размеров?

7.Как изменится методика расчета нестационарной теплопроводности, если телу свойственна тепловая анизотропия?

47

4. ФУНДАМЕНТАЛЬНЫЕ ОСНОВЫ КОНВЕКТИВНОГО ТЕПЛООБМЕНА

4.1. Основные понятия и определения

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

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

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

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

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

Re w0l . При этом ламинарный режим течения в трубах наблюдается при

Re 2300, а устойчивый турбулентный режим при Re 10000. Промежуточные значения числа Рейнольдса относятся к так называемому переходному режиму течения. Для внешних задач обтекания плоской пластины, переход от ла-

минарного к турбулентному течению происходит при Re 3 105.

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

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

48

нимать в виде непрерывного поля. Исключение составляют разреженные газы. Подробно теория теплообмена в таких средах рассмотрена в работах [5, 8].

Условия теплообмена определяются законом Ньютона-Рихмана согласно которому

q

t

c

t

ж

 

t

 

 

 

,

(4.1)

 

 

 

 

c

 

 

 

n

 

c

 

 

 

 

 

 

 

 

 

 

 

где – коэффициент теплоотдачи, Вт м2 К . В свою очередь,

 

 

t

 

 

 

 

 

(4.2)

 

 

 

t

n n 0.

Если задача решается в одномерной постановке (например, задана координата y), имеется естественный линейный масштаб l0, то выражение для определения коэффициента теплоотдачи можно записать в безразмерном виде:

 

Nu l0 Y Y 0 ,

(4.3.)

где Nu – безразмерный коэффициент теплоотдачи; Y y l0

– безразмерная ко-

ордината; t tж

tc tж – безразмерная температура.

 

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

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

Таблица 4.1

 

для различных видов конвективного теплообмена

Вид конвективного теплообмена

 

, Вт м2 К

Естественная конвекция в газах

 

2…30

Естественная конвекция жидкостей

 

50…1100

Вынужденная конвекция газов

 

5…500

Вынужденная конвекция жидкостей

 

200…2·104

Кипение воды

 

2·103…5·105

Теплоотдача от жидких металлов

 

102…3·104

Пленочная конденсация водяных паров

 

4·103…2·104

Капельная конденсация водяных паров

 

4·104…2·105

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

49