Файл: Быстрай Г.П. Термодинамика открытых систем Часть 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–910–8 м2с–1 значительно меньше коэффициента температуропроводности a~10–510–4 м2с–1. В результате скорость распространения концентрационных возмущений ϑD ~120 мс–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