Рис.10.1. К выводу балансовых уравнений переноса
В газовой или жидкой среде, находящейся в движении, выделим произвольный объѐм V, ограниченный поверхностью А (рис. 10.1). На поверхности А выделим элемент поверхности dА и представим его в векторной форме, умножив на единичный вектор ,
нормальный к этому элементу и направленный из объѐма, |
. |
Составим балансовое уравнение по типу: |
|
Накопление внутри объёма = Вход – Выход + Образование
Примем, что в произвольном объеме нет источников субстанции или стоков, т.е.
образование равно нулю.
Плотность потока субстанции через элементарную площадку будет Знак “–“ в этом произведении делает входящие потоки положительными, а
выходящие – отрицательными.
Результирующий поток массы или энергии (Вход минус Выход) будет получен суммированием всех потоков через замкнутую поверхность A:
|
|
q dA |
(10.12) |
A
Таким образом, физически этот интеграл представляет разницу между входящими и выходящими потоками субстанции через всю поверхность А.
4
Если в объѐме V происходит накопление субстанции, то это вызовет изменение потенциала переноса во времени ddt , которое для элементарного объѐма dV можно
представить как |
d |
dV , а для всего объема |
как интеграл: |
|||||
|
||||||||
|
dt |
|
|
|
|
|
|
|
|
|
M |
d |
dV |
|
(10.13) |
||
|
|
dt |
|
|||||
|
|
V |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Приравняв выражения (10.12) и (10.13), получим: |
||||||||
|
|
|
|
|
d |
|
|
|
|
|
q dA |
|
|
dV |
(10.14) |
||
|
|
dt |
||||||
|
|
A |
|
|
V |
|
|
|
|
|
|
|
|
|
|
||
Согласно теореме Остроградского-Гаусса, дающей преобразование интеграла,
взятого по объѐму V, ограниченному поверхностью A, в интеграл, взятый по этой поверхности, будем иметь:
|
|
|
q dA div q dV |
(10.15) |
|
AV
Сучѐтом (10.15) соотношение (10.14) примет вид:
|
|
|
(10.16) |
|
|
t |
div q dV |
||
V |
|
|
|
|
|
|
|
|
|
Интеграл, взятый по произвольному объѐму, может быть равен нулю только в случае равенства нулю подынтегральной функции:
|
|
|
|
t |
div q |
0 |
(10.17) |
Полученное выражение (10.17) и есть основное дифференциальное уравнение переноса субстанции – массы или энергии, как будет показано далее. В случае изотропных сплошных сред с его помощью можно получать поля температур или концентраций в однофазной среде. Искомой величиной является плотность потока субстанции , которая определяет удельный поток массы или энергии.
5
Дифференциальное уравнение конвективного теплообмена
(Уравнение Фурье-Кирхгофа)
Дифференциальное уравнение конвективного теплообмена является частным случаем уравнения переноса энергии в форме теплоты в однофазной сплошной изотропной среде. При этом теплоѐмкость сP, теплопроводность T и плотность среды считаются постоянными; отсутствует также перенос энергии в форме теплового излучения и объѐмное источники (стоки) теплоты.
Как было отмечено выше, потенциалом переноса теплоты является энтальпия единицы объѐма среды (ур. 10.4):
|
|
|
|
|
cP T |
|
|
(10.4) |
|||
Тогда с учѐтом выражений (10.7) и (10.9), будем иметь: |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
qT qMT qKT |
a grad cP T cP T v |
(10.18) |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
Где qT |
– плотность потока теплоты, представляющая векторную сумму молекулярной |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
qMT и конвективной qKT компонент. |
|
|
|
||||||||
Основное уравнение переноса субстанции (10.17) в этом случае примет следующий |
|||||||||||
вид: |
|
|
|
|
|
|
|
|
|
|
|
|
cP T |
|
|
|
|
|
|
|
|
(10.19) |
|
|
t |
|
div a grad cP T cP T v 0 |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
T |
|
|
|
|
|
|
|
|
|
|
|
t |
div |
a grad T |
T v 0 |
|
|
(10.20) |
||||
При a = const получим: |
|
|
|
|
|||||||
|
div a grad T a 2T |
|
|
(10.21) |
|||||||
2 - оператор Лапласа. |
|
|
|
|
|||||||
Дифференциальная операция |
|
сопоставляет скалярную функцию и |
|||||||||
скалярную функцию |
|
|
|
|
|
|
|
|
|||
|
|
|
|
2 |
|
2 |
2 |
|
|
|
|
|
div grad |
|
|
|
|
z 2 2 |
|
|
(10.22) |
||
|
x 2 |
y2 |
|
|
|||||||
Примем также, что гидродинамически |
среда является стационарной, |
тогда с |
|||||||||
учѐтом |
|
уравнение неразрывности имеет вид: |
|
||||||||
|
|
|
0 |
|
|
|
|
|
|
|
(10.23) |
|
div v |
|
|
|
|
|
|
|
|||
6
Дивергенцию от T v , как произведения векторной и скалярной величины, можно
представить в виде:
|
|
|
(10.24) |
div T v T divv |
v gradT |
||
С учѐтом (10.21), (10.23) и (10.24) выражение (10.20) примет вид: |
|
||
T |
|
|
(10.25) |
t |
v grad T a 2T |
||
|
|
|
|
Полученное выражение (10.25) называется дифференциальным уравнением конвективного теплообмена или уравнением Фурье-Кирхгофа. Оно является частным случаем дифференциального баланса энергии в форме теплоты в движущейся среде, где имеет место перенос энергии теплопроводностью.
Полная форма уравнения конвективного теплообмена в скалярном виде будет:
T |
v |
|
T |
v |
|
T |
v |
|
T |
|
2T |
|
2T |
|
2T |
|
(10.26) |
|||
|
x |
|
y |
|
z |
|
a |
|
2 |
|
2 |
|
2 |
|
||||||
t |
|
x |
|
y |
|
z |
|
x |
|
y |
|
z |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Левая часть этого соотношения представляет собой субстанциональную
производную: |
DT |
T vx |
T |
vy T |
vz T |
|
Dt |
t |
x |
y |
z |
Решением дифференциального уравнения конвективного теплообмена в общем |
|||||
виде является функция |
|
, которая представляет собой нестационарное поле |
|||
температур в движущейся среде. |
|
|
|||
В неподвижной среде |
|
и выражение (10.25) принимает вид: |
|||
T |
a 2T |
|
|
(10.27) |
|
t |
|
|
|
|
|
Уравнение (10.27) описывает нестационарное температурное поле в неподвижной среде, может применяться также для твѐрдых тел и называется уравнением
нестационарной теплопроводности.
Уравнение для установившегося процесса в неподвижной среде или в твердом теле имеет вид:
2T 0 |
(10.28) |
Отметим, что согласно уравнению (10.27), локальное изменение температуры T
t
пропорционально коэффициенту температуропроводности a, который, таким образом,
характеризует теплоинерционные свойства среды. При прочих равных условиях быстрее нагреется или охладится то тело, которое имеет бόльший коэффициент температуропроводности.
7
Граничное условие
Уравнения Фурье-Кирхгофа на практике используется совместно с граничным условием, т.е. условием на границе среды у неподвижной твердой стенки. Вблизи твердой стенки теплота передается только теплопроводностью внутри пограничного слоя.
Следовательно, по закону Фурье:
Q |
|
T |
dA gradT dA |
(10.29) |
|
||||
|
T n |
T |
|
|
В тоже время, количество теплоты, передаваемой из ядра потока к твердой стенке,
можно выразить законом Ньютона (уравнение теплоотдачи):
Q T Tст dA |
(10.30) |
|
|
Если перенос тепла стационарный, это один и тот же тепловой поток: |
|
|
|
Q T gradT dA T Tст dA |
|
T gradT T Tст |
(10.31) |
Это и будет граничным условием, дополняющим уравнение Фурье-Кирхгофа.
Элементы теории подобия в теплообмене
Рассмотрим гидродинамически одномерный поток жидкости. Запишем уравнение Фурье-Кирхгофа:
T |
vz |
T |
a 2T |
(10.32) |
t |
|
z |
|
|
Получим |
приближенное решение этого |
уравнения методами теории |
||
подобия.
Для этого зададим константы подобия, выражающие отношения величин,
входящих в уравнение Фурье-Кирхгофа: al , at , aT, aa , av .
Умножим каждый из элементов дифференциального уравнения (10.32) на
соответствующую константу подобия, причем последняя как постоянная величина,
выносится за знак дифференциала.
|
aT |
|
T |
avaT |
v |
|
T |
|
aaaT |
a 2T |
(10.33) |
|
|
z |
|
|
|||||||
|
|
|
|
|
|
|
al2 |
|
|||
|
at |
t |
al |
|
z |
|
|
||||
Для сохранения тождественности полученного и исходного уравнений необходимо
выполнение следующего условия:
8