Файл: Быстрай Г.П. Термодинамика открытых систем Часть 2 (2006).pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 29.06.2024
Просмотров: 356
Скачиваний: 1
Термодинамика открытых систем
Глава 4. ТЕРМОДИНАМИКА ЛОКАЛЬНО−НЕРАВНОВЕСНЫХ ПРОЦЕССОВ
ПЕРЕНОСА 4.1. Основные положения
Одной из наиболее последовательных и детально разработанных термодинамических теорий, не опирающихся на принцип локального равновесия является так называемая “расширенная необратимая термодинамика” (РНТ) (см. обзор [25]). В рамках РНТ рассматриваются следующие дифференциальные уравнения для диссипативных потоков релаксационного типа [25]:
q + τ |
T |
∂q |
= −λ T , |
(4.1) |
||||
∂t |
||||||||
j + τ |
D |
∂j = −D C , |
(4.2) |
|||||
∂t |
||||||||
p + τϑ |
∂p |
= −ζ ϑ |
(4.3) |
|||||
∂t |
||||||||
P + τ2 p |
∂P |
ϑ |
• |
|||||
= −2ηϑ , |
(4.4) |
|||||||
∂t |
||||||||
где λ − коэффициент теплопроводности, D − коэффициент диффузии, τT , τD , τϑ , τ2 p − времена релаксации соответствую-
щих диссипативных потоков; P=pδ+Pυ, δ − единичный тензор, p − вязкое давление (1/3 следа тензора P), Pυ−часть тензора P со следом, равным нулю, ζ − объемная вязкость, η − сдвиговая вяз-
•
кость, υ − симметрическая часть градиента скорости. При этом потоки уже не определяются градиентом соответствующего термодинамического потенциала переноса, а являются решениями эволюционных уравнений (4.1) – (4.4). Эти уравнения описывают процессы релаксации диссипативных потоков к своим локальноравновесным значениям. Например, в системе с нулевым градиентом концентрации начальное значение массового потока j0 релаксирует к равновесному значению j=0 по экспоненциальному закону:
65
Термодинамика открытых систем
t |
||||||
j(t) = |
j0 |
− |
||||
exp |
τT |
. |
||||
Уравнение Максвелла-Катанео (4.1) может быть представлено как приближение первого порядка при разложении в ряд Фурье
по τT более общего соотношения: q(t + τT ) = −λ T . Послед-
нее означает, что между тепловым потоком и градиентом температуры существует временной сдвиг, равный времени релаксации. Уравнения (4.1)-(4.4) описывают простейшие случаи одноступенчатой (или одностадийной) релаксации и не учитывают как перекрестных, так и пространственно-нелокальных эффектов.
Учет перекрестных эффектов в приближении РНТ позволяет представить уравнения возмущенного движения в виде
Ji |
+ τi |
dJi |
= aie |
∂S |
+ aii |
∂S |
; |
|||||||||||
dt |
∂ξe |
∂ξi |
||||||||||||||||
J e + τe |
dJ e |
= aee |
∂S |
+ aei |
∂S |
; |
(4.5) |
|||||||||||
dt |
∂ξi |
|||||||||||||||||
∂ξe |
||||||||||||||||||
здесь τi, τe – время релаксации внутренних и внешних термодинамических потоков. Характерные пространственно-временные масштабы L, h, τe=t0 и τi определяют две характерные скорости
[25]
ϑe = τLe , ϑi = τhi .
Скоростьϑe , представляющая собой отношение макромасштабов
рассматриваемого процесса, характеризует линейную скорость изменения параметров системы, вызванную внешними причинами. Например, это может быть скорость перемещения изотерм при движении источника тепловыделения в теплопроводящей
среде. Отношение микропараметровϑi является внутренней характеристикой самой системы и не зависит от внешних условий. Величина ϑi − скорость распространения возмущений потенциала переноса для внутреннего потока. Например, в газах ха-
66
Термодинамика открытых систем
рактерными микропараметрами среды как для процессов теплопереноса, так и процессов массопереноса, являются средняя длина свободного пробега h и время между двумя последовательны-
ми столкновениями молекул τ. Поэтому ϑi – средняя скорость молекул газа, причем ϑi =3D/h=3a/h (D – коэффициент диффу-
зии), поскольку в газах a=D. В расплавах металлов коэффициент диффузии примеси D~10–9−10–8 м2с–1 значительно меньше коэффициента температуропроводности a~10–5−10–4 м2с–1. В результате скорость распространения концентрационных возмущений ϑD ~1−20 мс–1 много меньше скорости распространения тепло-
вых возмущений ϑT ~10–3– 104 мс–1. В такой системе сначала
устанавливаются локально-равновесные значения потока, обладающего минимальным временем релаксации, а только в последующем – локально-равновесные значения другого потока. При
этом характерное время τD ~h/ ϑD много больше, чем время тепловой релаксации τT ~h/ ϑT . Это означает, что в такой системе, сначала, через время порядкаτT устанавливаются локально– равновесные значения температуры и только через время порядка τD – локально-равновесные значения концентрации. После-
довательная релаксация к тепловому, а лишь затем к диффузионному равновесию может возникнуть в системах со сложной структурой, например, в полимерах и капиллярно-пористых средах.
В результате для локально-неравновесных систем скорость изменения энтропии, объединяющая все внешние и внутренние потоки, также будет зависеть от времени релаксации τr, которое связано с одним из наибольших времен релаксации потоков (с самым длительным лимитирующим процессом):
G + τr |
dG |
= −Je X e + Ji X i +σ. |
(4.6) |
|
dt |
||||
Прежде чем использовать полученное уравнение (4.6), обратимся к уравнению теплопроводности для локально-неравновесных систем.
67
Термодинамика открытых систем
Следует отметить, что для локально-неравновесных систем, описываемых уравнениями возмущенного движения - нестационарными уравнениями Онзагера (4.5) можно также сформулировать теоремы 1-3.
4.2.Гиперболическое уравнение теплопроводности
систочником тепла
Если приближения локального равновесия не выполняются τ<<t0, то этим эффектом, учитывающим время релаксации, пренебречь уже нельзя. В этом случае процесс переноса тепла описывается уравнением гиперболического типа
∂T |
+ τT |
∂2T |
= a 2T + |
W |
, |
(4.7) |
|
∂t |
∂t 2 |
CV ρ |
|||||
отличающимся от параболического уравнения теплопроводности наличием второй производной температуры по времени и со-
держащим время релаксации теплового потока τT , здесь
[W ] = Вт/ м3 , [W / Cvρ] = K / c , [ a ] = м2 / c , [λ] = Вт/ мK .
Справедливость данного уравнения теплопроводности также можно доказать в рамках термодинамики неравновесных процессов.
T, K
ϑT
х
Рис.4.1. Фронт повышения температуры при отсутствии локального равновесия.
68
Термодинамика открытых систем
Однако, прежде чем это делать, отметим, что уравнение теплопроводности гиперболического вида описывает два семейства характеристик x ± ϑT t = const . Это означает, что тепловой
сигнал (или высокочастотные тепловые возмущения) распространяются в локально-неравновесных условиях с конечной ско-
ростью ϑT =( a / τT )1/ 2 (рис.4.1)[26].
Иными словами, уравнение (4.7) говорит о том, что изменения температуры на поверхности полубесконечного тела будут
распространяться в его объем с конечной скоростью ϑT , в отли-
чие от параболического уравнения теплопроводности. Для уравнения (4.7) можно ввести понятие теплового пограничного слоя
δT ~ ϑT τT = (aτT )1/ 2 − именно на такое расстояние распространится тепловое возмущение за характерное время τT . На фронте
распространяющегося температурного возмущения терпят сильный разрыв не только температура и энтропия, но и свободная
• ••
энергия F, а также ее производные F, F . Это и свидетельству-
ет, что для таких процессов принцип локального равновесия не выполняется.
1. Термодинамическое обоснование гиперболического уравнения теплопроводности. При выводе гиперболического уравнения будем также исходить из закона сохранения энергии для неравновесных систем с источником (1.11). После дифференцирования по времени этого уравнения, имея в виду что
•
J q = Lqq X q − τT J q ,
получаем с учетом (3.3) дифференциальное уравнение сохранения энергии [18] в случае фиксированных потоков:
∂S ∂T |
∂ |
2 |
T |
∂X q |
• |
• |
∂σ |
e |
|||||||||||
+ S |
=T J |
−T |
X |
q |
J |
q |
τ |
T |
+T |
,(4.8) |
|||||||||
∂t ∂t |
2 |
q ∂t |
∂t |
||||||||||||||||
∂t |
0 |
0 |
0 |
||||||||||||||||
или −
69
Термодинамика открытых систем
∂T |
S + δS ∂2T |
2 |
||||||||||||||
+ |
=T0 |
λT |
− |
1 dT |
− |
|||||||||||
∂t |
• |
∂t 2 |
T 2 |
|||||||||||||
0 |
dx |
|||||||||||||||
S |
0 |
|||||||||||||||
+T |
∂σe |
∂t |
. |
|||||||||||||
∂t |
||||||||||||||||
0 |
∂S |
|||||||||||||||
12 d 2T ∂∂t +
T0 dtdx S
Здесь второе слагаемое в правой части (4.8) перенесено в левую часть и оно вошло в структуру второго члена в левой части. Ло- кально-неравновесная энтропия равна S(t)=S+δS, где приращение δS находится из выражения
• |
2 |
|||||||||||||
• |
J q |
λ ∂ |
T |
|||||||||||
δS = T X |
τ |
= |
τ |
. |
||||||||||
2T0 ∂x2 |
||||||||||||||
0 |
q •• |
T |
T |
|||||||||||
T |
||||||||||||||
Таким образом, локально-неравновесная |
энтропия |
является |
||||||||||||
• • |
||||||||||||||
функцией параметров неравновесия X q ,J q , или |
||||||||||||||
S( t ) = Seq |
+ |
λτT |
2T |
, |
(4.9) |
|||||||||
2T0 |
||||||||||||||
здесь [ S ] = Дж / м3K . В результате после деления правой и
•
левой частей уравнения (4.8) на S получаем гиперболическое уравнение теплопроводности (4.7), в котором время релаксации определяется как характеристика изменения энтропии
•
τT = δS / S . В результате источник тепла как и в параболическом уравнении теплопроводности определен в явном виде:
W |
∂σe |
T 2 |
∂σe |
||||||||||
= T |
= |
0 |
, |
||||||||||
C |
ρ |
∂S |
C |
∂T |
|||||||||
0 |
V |
ρ |
|||||||||||
V |
V |
V |
|||||||||||
или
W= T02 ∂∂σTe V .
Вслучае независимости источника от температуры из последнего выражения после интегрирования получаем равенство, связы-
70
Термодинамика открытых систем
вающее введенную нами величину σe с интенсивностью внутренних источников тепла
σe = W .
T0
4.3. Термодинамика процессов переноса тепла при отсутствии локального равновесия
Знание пространственно распределенной температуры, которая находится из гиперболического уравнения теплопроводности (4.7), позволяет для локально-неравновесных процессов переноса массы вычислить скорость изменения энтропии, свободной энергии, производство энтропии, которые также зависят от координат и времени. Материал, изложенный в этом параграфе, следует сравнить с изложенным в параграфе 3.2.
1. Скорость изменения энтропии. Скорость изменения энтропии для гиперболического уравнения теплопроводности
(4.7) равна
dS |
= σe + J q X q = |
W |
λ |
( T )2 − |
λ |
• |
||
+ |
T (T )τT .(4.10) |
|||||||
2 |
2 |
|||||||
dt |
T0 |
|||||||
T0 |
T0 |
|||||||
При этом |
производство энтропии является знакоположительной |
|||||||
функцией. Функция источников в такой задаче также знакоопределена – она знакоположительна. В уравнении (4.10)
[ σe ] = [ σi ] = Дж / м3 Kc .
Второй закон термодинамики применительно к уравнению теплопроводности (4.7) выражается неравенством
σi = J q X q = |
λ |
( T )2 − |
λ |
• |
|
T (T )τT ≥ 0 . |
|||||
T 2 |
T 2 |
||||
0 |
0 |
||||
2. Скорость изменения свободной энергии. Скорость изменения свободной энергии для локального объема при суще-
71
Термодинамика открытых систем
ствующем градиенте температуры для уравнения теплопроводности равна
dF |
= −T ( σe |
+ J |
q |
X |
q |
) , |
||||||||||
dt |
0 |
|||||||||||||||
или |
||||||||||||||||
• |
||||||||||||||||
dF |
W |
λ |
( T )2 − |
λ |
||||||||||||
= −T |
+ |
T (T )τ |
, |
|||||||||||||
T |
||||||||||||||||
dt |
0 |
2 |
2 |
|||||||||||||
T0 |
||||||||||||||||
T0 |
T0 |
|||||||||||||||
здесь температура в трехмерной задаче является функцией времени и координат T =T( x, y,z,t ) , от координат также зависит
• |
• |
||
скорость изменение свободной энергииF |
= F( x, y,z,t ) . |
||
3. Теорема Пригожина. |
Производство энтропии |
||
( σi |
• |
||
= J q X q , J q = Lqq X q − τT J q ) |
при постоянной мощности |
||
теплового источника (W = const ) стремится убывать и принимает минимальное положительное значение в стационарном состоянии в соответствии с термодинамическим уравнением
dσ |
i |
1 |
∂W |
1 ∂ |
2 |
F |
• |
|||||||||||||
= − |
− |
, |
[ F ] = Вт / м3 . |
|||||||||||||||||
dt |
T0 ∂t 2 |
|||||||||||||||||||
T0 ∂t |
||||||||||||||||||||
4. Функция Релея (мера рассеивания полной энергии). |
||||||||||||||||||||
При теплопроводности для уравнения (1.21) |
||||||||||||||||||||
dΛ |
= −2Φ , Λ = ΛF + Λe , |
(3.8) |
||||||||||||||||||
dt |
||||||||||||||||||||
функция Релея (1.22) примет вид |
• |
|||||||||||||||||||
T |
λ |
λ |
||||||||||||||||||
Φ = |
0 |
J q X q = |
( T )2 |
− |
J q T τT . |
|||||||||||||||
2T0 |
||||||||||||||||||||
2 |
2T0 |
|||||||||||||||||||
Здесь Λe термодинамический потенциал внутренних источ-
• •
ников Tσe ≡ Λe , Λe =W а ΛF = F − F0 - термодинамический потенциал (свободная энергия) неравновесного состояния.
72
Термодинамика открытых систем
4. Вариационный принцип. Несложно показать,
что уравнение (3.8), на котором по сути и основан вывод уравнения теплопроводности, тождественно некоторому термодинамическому вариационному принципу
δ Λ• + 2Φ Jq = 0 .
5. Энтропия, свободная энергия для локально-
неравновесных состояний. Как ведет себя свободная энергия для описываемых локалъно-неравновесных состояний с потоками тепла? Она равна в соответствии с (4.9)
F(t) = U0 −T0 S( t ) = F0 − τT a 2T . 2T0
Для процессов, связанных с понижением температуры
•
(T < 0 ), из уравнения теплопроводности (4.7) в самой простой задаче, когда a 2T >>W/Cvρ, τT / ∆t <<1 следует неравенство
a 2T < 0 . Поэтому для таких локально-неравновесных процессов энтропия меньше, а свободная энергия F(t) =U0 −Т0 S(t) больше, чем для локально-равновесных:
S( t ) = S |
eq |
− λτT |
2T |
, |
F( t ) = F + |
τT a |
2T |
; |
|||||
2T0 |
0 |
2T0 |
|||||||||||
F0 = U0 −T 0S0 . |
|||||||||||||
Это и означает, что при низких температурах все |
равновесные |
||||||||||||
системы должны быть так или иначе упорядочены, а свободная энергия для них принимает минимальное значение.
Для процессов, связанных с повышением температуры
•
(T > 0 ) из уравнения теплопроводности для выше указанных
условий следует неравенство a 2T > 0 . Поэтому для таких ло- кально-неравновесных процессов энтропия больше, а свободная энергия F(t) =U0 −Т0 S(t) меньше, чем для локально-
равновесных:
73
Термодинамика открытых систем |
|||||||||||||
S( t ) = S |
eq |
+ λτT |
2T |
, |
F( t ) = F − |
τT a |
2T |
. |
|||||
2T0 |
0 |
2T0 |
|||||||||||
Это и означает, что при высоких температурах все равновесные системы должны быть упорядочены меньше, чем при низких, свободная энергия принимает при этом максимальное значение.
6. Температура для локально-неравновесных состоя-
ний. Принимая во внимание, что локально-равновесная темпера- |
|||||||
тура T определяется из соотношения T −1 |
= (∂S |
eq |
/ ∂U ) , ло- |
||||
кально-неравновесная температура θ |
0 |
V |
|||||
может быть найдена из из |
|||||||
(4.9) дифференцированием S(t) по U: |
|||||||
1 = |
1 |
+ |
τT a |
2T . |
(4.11) |
||
T0 |
2T 2 |
||||||
θ |
|||||||
0 |
|||||||
Для процессов, связанных с высокоинтенсивным понижением температуры из уравнения теплопроводности следует
неравенство a 2T < 0 , в результате получаем, что локальнонеравновесное значение температуры θ<T, т.е. меньше чем ло- кально-равновесное:
T0 |
• |
|||||
θ = |
, T |
< 0 . |
||||
τT a |
2T |
|||||
1 + |
||||||
2T |
||||||
0 |
||||||
Для процессов, связанных с высокоинтенсивным повышением температуры из уравнения теплопроводности следует неравенство a 2T > 0 , в результате получаем, что локальнонеравновесное значение температуры θ>T, т.е. больше чем ло- кально-равновесное:
T0 |
• |
|||||
θ = |
, T > 0 . |
|||||
1 − |
τT a |
2T |
||||
2T |
||||||
0 |
Следует подчеркнуть [25], что, строго говоря, в уравнении для потока тепла (4.1) должен стоять градиент локальнонеравновесной температуры θ, а не T . Обычно в практиче-
74