Файл: Быстрай Г.П. Термодинамика открытых систем Часть 2 (2006).pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 29.06.2024
Просмотров: 82
Скачиваний: 0
Термодинамика открытых систем
Глава 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