Тепломассообмен. Нестационарная теплопроводность

1. Тепломассообмен 6

Нестационарная теплопроводность:
● Охлаждение (нагревание) пластины
● Охлаждение (нагревание) цилиндра
● Регулярный режим

2. Внутренняя задача

● Частный случай (А):
Bi
(практически Bi >100): Bi – число (критерий) Био:
соотношение конвективной
/
1
(Bi
) ; 0; , теплоотдачи снаружи и теплопроводности внутри тела.
1/
В данном случае очень интенсивное наружное охлаждение,
поэтому температура поверхности пластины, погруженной
в жидкость, сразу становится равной температуре жидкости.
Распределение температур в пластине зависит от ее
теплопроводности λ и геометрических размеров
от условий внутри пластины (внутренняя задача).
, то есть

3. А) Внутренняя задача В) Внешняя задача

А) Внутренняя задача Bi
В) Внешняя задача Bi 0
Fo1 Fo2 Fo3
Fo 0
1
Fo 0
1
Fo1
Fo1
А)
Fo2
Fo2
Fo3
Fo3
0
В)
x
0
x

4. Внешняя задача

● Частный случай (В):
Bi 0
(практически Bi < 0,1),
Bi ( / ) 0; 0;
1/
теплопроводность (λ)
значительная.
Из-за высокого коэффициента теплопроводности пластины
температуры в ней быстро выравниваются. Охлаждение
слабое и все зависит от внешнего коэффициента
конвективной теплоотдачи
(внешняя задача).
Обозначения: - половина толщины пластины, м;
- теплопроводность пластины, Вт/(мК);
- коэффициент конвективной теплоотдачи, Вт/(м²К).

5. Средний случай

● Частный случай (С): 0,1 Bi 100.
Интенсивность охлаждения зависит и от внутреннего
термического сопротивления
( / )
и внешнего -
(1/ )
.
Распределение температур в пластине для этого случая
показано на следующем слайде.
Из уравнения (13 ТМО 4) следует, что для любого момента
времени
распределение температур имеет вид
симметричной кривой с максимумом на оси пластины (Х=0).
Касательные к кривым в точках
+А и –А на расстоянии
этих касательных
X 1 проходят через точки
X 1/ Bi
tg X 1
X0
и тангенс угла наклона
(см. следующий слайд).

6. Температурное поле в пластине

Fo 0
1
Fo1
Fo2
X 0
X 1
A
0
X0
2
X0
A
x

7. Теплота, отданная от пластины к жидкости

Теплота, отданная с обеих сторон пластины к окружающей
, равна изменению
внутренней энергии пластины, Дж: Q 2 f c(t t ), (1)
п
0
ж
где 2 f V объем пластины, м³; V m ее масса, кг;
- половина толщины пластины, м; f - ее поперечное
сечение, м²; - плотность материала пластины, кг/м³;
ее жидкости за время от
0
до
с – теплоемкость материала пластины, Дж/(кгК).
Тогда за любой промежуток времени от
0
в безразмерной форме: (от Fo = 0 до Fo1).
до
1
или
За это время внутренняя энергия пластины изменится на, Дж:
t1ср tж
Q Qп Q1 2 f c(t0 tж )(1
),
t0 t ж
(2)

8. Средняя по толщине пластины безразмерная избыточная температура

Q Qп (1 1ср ),
или
где
ср
t
tж
ср
1
1
t0 t ж
(3)
средняя по толщине пластины
безразмерная избыточная температура
в момент времени
1 .
В соответствии с теоремой о среднем, средняя безразмерная
избыточная температура
пластины найдется как:
где
1X
dx,
X0
ср
по уравнению (13 ТМО 4), тогда при
(4)
Fo 0,3 можно
ограничиться только первым членом ряда, то есть:
1ср
2sin 2 1
2
exp( 12 Fo).
1 1 sin 1 cos 1
(5)

9. Охлаждение (нагревание) бесконечного цилиндра

Бесконечный цилиндр
( 2r0 ) радиусом r0 отдает теплоту
окружающей его жидкости при:
Const;tж Const; 0 t t0 Const.
Дифференциальное уравнение теплопроводности в полярных
(цилиндрических) координатах для бесконечного цилиндра:
2 1
a( 2
).
r r
r
Начальные условия: при
Граничные условия: при
(6)
0;0 r r0 0 t0 tж Const
r 0; 0 ( )r 0 0;
r
r r0; 0 ( )r r r r .
0
r 0
.

10. Числа подобия для охлаждения (нагревания) цилиндра

Решение в общем виде: f ( R, Bi, Fo),
где R r безразмерный радиус цилиндра;
(8)
r0
r0
Bi
Fo
a
2
r0
число (критерий) Био, который представляет
собой соотношение конвективной теплоотдачи
снаружи и теплопроводности внутри цилиндра.
число (критерий) Фурье – безразмерное время.
По аналогии с пластиной, при Bi 0 (практически
при Bi <0,1) степенные ряды становятся настолько быстро
сходящимися, что можно ограничиться только первым
членом ряда .
1

11. Безразмерные избыточные температуры

Безразмерные избыточные температуры:
R 0 N0 (Bi)exp( 12 Fo);
● на оси цилиндра
● на поверхности цилиндра P (Bi)exp( 2 Fo).
R 1
0
(9)
(10)
1
Функции N ( Bi); P ( Bi) табулированы и приведены в
0
0
справочниках.
По аналогии с бесконечной пластиной зависимости (9) и (10)
линейные в логарифмических координатах и по ним можно
найти:
t
t
r 0
r r
ж
tr 0 t ж
; r r 0
,
0
t0 t ж
t0 t ж
или в безразмерном виде, по графикам:
R 0 Ф1(Bi, Fo); R 1 Ф2 (Bi, Fo).

12. Теплота, отданная от цилиндра к окружающей его жидкости

Теплота, отданная цилиндром за время от 0 до ,
равна изменению внутренней энергии цилиндра, Дж:
Qц r02 c(t0 tж ),
ср
а за время от 0 до 1 : Q Qц (1 1 ),
где
1ср
При
(11)
(12)
t1ср tж
средняя по цилиндру безразмерная
t0 t ж
Fo 0,25:
избыточная температура в момент времени .
1
1ср
4Bi 2
2 Fo).
2 2
exp(
1
1 ( 1 Bi 2 )
(13)
Аналогично есть решение и для охлаждения (нагревания) шара.

13. Регулярный режим охлаждения

n( )
Iст.
IIстадия
X 0
n( 1)
X 1
n( 2 )
0
1
2

14. Регулярный режим охлаждения (нагревания) тел

Анализ решений для охлаждения (нагревания) тел разной
формы показывает, что все они представляют сумму
бесконечного ряда, члены которого соответствуют быстро
убывающим экспоненциальным функциям.
Например, для бесконечной пластины
при
Const;tж Const
было получено:
где
x
An cos( n )e
n 1
a
n2
2
,
An константа для каждого члена ряда, которая
находится из начальных условий.
Множитель cos(
n
x
) U n зависит только от координаты Х.
(1)

15. I – неупорядоченная стадия охлаждения

n2 a
mn является постоянным, положительным,
Комплекс
2
вещественным числом:
Тогда уравнение (1)
запишется в виде:
m1 m2 m3 ... mn
AnU n e
, где
mn
n 1,2,3,...
.
.
(2)
n 1
Уравнение (2) справедливо для тел разной геометрии,
которая учитывается видом сомножителей A ,U .
n
n
При малых значениях времени от 0 до
1
изменение температур зависит от начального распределения
температур в теле.
В этом случае поле температур будет определяться не только
первым, но и последующими членами ряда (2)
«I – неупорядоченная стадия охлаждения».

16. II стадия охлаждения – регулярный режим

Но начиная с некоторого момента времени 1 начальные
условия играют второстепенную роль, процесс определяется
интенсивностью охлаждения и физическими свойствами тела.
Тогда температурное поле достаточно точно описывается
только первым членом ряда «II стадия охлаждения –
регулярный режим»,
m1
AU
e
1 1
для которого:
(3)
Логарифмируя (3) и опуская индексы, получим:
n m n( AU ), или
n m c( x, y, z),
то есть в полулогарифмических координатах эта
зависимость – прямолинейная.
(4)

17. III стадия охлаждения – стационарный режим

При длительном охлаждении (
или
Fo
)
все точки тела принимают одинаковую температуру, равную
температуре окружающей жидкости
tж .
Это III стадия охлаждения – стационарный режим.
Для регулярного режима после дифференцирования
уравнения (2) имеем:
1
m Const,
то есть относительная скорость изменения температуры
равняется константе «m», не зависящей от координат
и времени.
«m», 1/с – темп охлаждения.
(5)

18. Темп охлаждения

Если есть экспериментальный график изменения
избыточной температуры тела во времени (см. слайд 13), то
темп охлаждения в стадии
регулярного режима, 1/с:
m
n 1 n 2
tg .
2 1
(6)
Зависимость темпа охлаждения от физических свойств
тела, его геометрии, размеров и условий теплообмена на
поверхности можно найти из теплового баланса.
Изменение внутренней
энергии тела, Дж:
где
vср
vср
dQ c V
d ,
средняя по объему избыточная температура, К.
Теплота (7) отдается от поверхности тела к окружающей
его жидкости.
(7)

19. Первая теорема Кондратьева

По уравнению конвективной
dQ ср Fср Fd .
теплоотдачи, Дж:
(8)
ср
Здесь средняя по поверхности избыточная
F
температура; ср средний коэффициент теплоотдачи.
Приравнивая (7) и (8) с учетом того, что c V C - полная
теплоемкость тела, Дж/кг; F / v коэффициент
неравномерности распределения температуры в теле, имеем:
1 vср
ср F (9) то есть при ср Const темп
ср
m
,
темп охлаждения однородного
C
v
изотропного тела (относительная скорость охлаждения)
пропорционален коэффициенту теплоотдачи, поверхности
тела и обратно пропорционален его полной теплоемкости
(первая теорема Кондратьева).

20. Коэффициент неравномерности распределения температуры

Итак коэффициент неравномерности
распределения температуры
m
в теле из (9):
Как же он зависит от числа Био?
А)
ср
C
.
Bi 0
(практически Bi < 0,1) – внешняя задача:
распределение температур не зависит
от геометрических
размеров тела и его
физических свойств
Fср
vср
Fср
ср 1.
v
(10)

21. Диапазон изменения коэффициента

В)
Bi
(практически Bi > 100) – внутренняя задача:
распределение температур зависит только от
геометрических размеров тела и его физических свойств.
Из-за высокого внешнего коэффициента теплоотдачи
;t t ; 0 .
c
ж
Следовательно, в общем случае, коэффициент
изменяться от (1 при Bi = 0) до (0 при Bi = ∞).
См. следующий слайд.
будет

22. Зависимость

f (Bi)
1
0
Bi

23. Вторая теорема Кондратьева

При Bi ;( ) темп охлаждения тела «m»
становится пропорциональным его коэффициенту
температуропроводности "a"
a km .
(вторая теорема Кондратьева)
Коэффициент пропорциональности зависит только от
геометрии и размеров тела.
2
m
a
Для бесконечной пластины: m
,
ср
2
где половина толщины пластины,
тогда с учетом того, что:
получим: Bi ;ctg 0; /2;
,
Bi 0;ctg ; 0
то есть в диапазоне Bi = 0 - ∞:
a
ctg
0 /2
.
Bi
,
(11)
(12)

24. Регулярные режимы I, II, III родов

При Bi (практически при Bi>100) из (12) для
/2: m ( /2 )2 a
откуда: k ( /2 ) 2
, то есть
a ( /2 ) 2 m km
(13)
- коэффициент пропорциональности для пластины.
Есть также свои выражения для цилиндра и шара.
На основе теории регулярного режима разработаны
экспериментальные методы определения теплопроводности
и коэффициентов температуропроводности тел.
Const - регулярный режим I рода;
tж tж0 b - регулярный режим II рода;
tж tж0 tm cos( ) - регулярный режим III рода (
При: t
ж
tm
- амплитуда колебаний температуры жидкости).
- частота и