книги из ГПНТБ / Дьярмати, И. Неравновесная термодинамика. Теория поля и вариационные принципы
.pdfИнтегральный принцип термодинамики |
217 |
уравнения баланса (6.12) формула (6.47) приводит к ва риационной задаче
Грс |
дТ |
- ~ (ѴГ)2 dV = О, |
|
ѵ ~дТ |
|||
|
|
уравнением Эйлера—Лагранжа для которой является дифференциальное уравнение
pCt, - ^ + V-(LrVr) = 0. |
(6.49) |
Оно описывает теплопроводность в твердых телах в уни версальном виде. Поэтому при Г = \/Т уравнение (6.49) включает в себя энтропийное представление (6.36), при Г = —- ln Т — энергетическое представление (6.32) и, на конец, при Г = — Т — уравнение Фурье в традиционной форме (6.15). Еще раз подчеркиваем, что локальные ве личины а г и Ч;г, определяемые соотношениями (6.44) и (6.45), имеют размерность производства энтропии толь ко при Г = 1/Т и Lp 2= Lqq, т. е. в энтропийном представ лении. Это очевидно из выражений (6.5) и (6.6), лежа щих в основе подробно рассмотренных нами трех пред ставлений.
Метод описания, применяемый в универсальном «Г»- представлении, очень важен. Действительно, он не только содержит в сжатой форме основные результаты различ ных представлений, что открывает возможность единого теоретического описания, но и представляет интерес в практическом отношении. Например, при некоторых до пущениях, зависящих от самой задачи, он позволяет про стейшим способом исследовать проблемы теплопровод ности, которые могут описываться только нелинейным уравнением
Р ^ ^ = Ѵ-[Л(7’)ѴП
так как коэффициент теплопроводности к’ = К(Т) зави сит от температуры. Если нам удастся установить зави симость Х(Т) (например, экспериментально в конкрет ных случаях), то мы можем построить такое «Г»-пред- ставление, в котором коэффициент Lr , согласно (6.43), постоянен. В этом случае нелинейная вариационная
2 1 8 |
Глава VI |
задача в универсальном «Г»-представлении сводится к линейной вариационной проблеме.
Естественно, что нелинейные вариационные задачи, справедливые в том случае, когда X зависит от темпера туры, и приводящие к нелинейной форме уравнений пе реноса, подобной последнему уравнению, представляют собой тип задач, значительно отличающихся от тех, воз можная теория которых упоминалась в гл. V, § 5. Это достаточно очевидно, поскольку нелинейные конститутив ные уравнения (5.82) определяют нелинейные соотноше ния между потоками и силами. Если же назвать это не линейностью в точном смысле слова, то необходимо ска зать, что при исследовании проблемы в универсальном «["»-представлении можно исключить только нелинейно сти типа X = Х(Т) (или другие подобные), т. е. нелиней ности более слабые. Конечно, введение универсального «["»-представления возможно не только для теплопровод ности, но (при выполнении соответствующих условий) и для других проблем переноса. Вот почему «["»-представ ление, разработанное Фархашем [85], очень полезно в практическом отношении.
§2 . Формулировка интегрального принципа
Впредыдущем параграфе, исходя из интегральной формы принципа наименьшего рассеяния энергии (6.1), заданного в представлении через силы для частного слу чая теплопроводности, мы сформулировали новый инте гральный принцип. Этот принцип выражен в различных
представлениях вариационными условиями (6.22), (6.30), (6.38) и (6.47). Рассмотрим подробную запись (6.37) плотности лагранжиана, относящуюся к интегральному принципу, сформулированному в энтропийном представ лении (6.38). Используя последние выражения в (6.6) и (6.10), ее можно записать в компактной форме
& U T = p s - y q. |
(6.50) |
Сравнивая это выражение с исходным (6.33), мы видим, что плотность лагранжиана совпадает с подынтеграль ным выражением в объемном интеграле в (6.33). В (6.33)
Интегральный принцип термодинамики |
219 |
поверхностный интеграл влияет только на граничные условия (его обращение в нуль означает, что Js или Jq и Т фиксированы вдоль граничной поверхности во время б-процесса). Следовательно, граничные условия в выра жении (6.50) для плотности лагранжиана уже учтены. То же самое, очевидно, справедливо и для интеграль ного принципа, записанного в форме (6.38).
Теперь обобщим наши результаты, учитывая усло вия, заданные для случая теплопроводности, и таким образом сформулируем интегральный принцип термоди намики. Величину s всегда можно задать соотношением (3.25), которое следует из соотношения Гиббса для си стем, находящихся в состоянии целлулярного равнове сия; потенциал рассеяния Ч7 также можно считать из вестным. Отсюда с очевидностью следует, что в энтро
пийном представлении функцию |
|
2 = ps — W, |
(6.51а) |
или в более общем случае функцию |
|
2 ’ = а — Чг, |
(6.516) |
можно рассматривать как термодинамическую плотность лагранжиана. Другими словами, это означает, что тер модинамическая плотность лагранжиана равна разности между «кинетической» частью ps, содержащей произ водную по времени, и «потенциальной» частью, которая определяется потенциалом рассеяния ЧГ
Обозначим через Гг интенсивные параметры, которые входят в энтропийное представление (3.25) и градиента ми которых являются соответствующие силы Хі. Тогда в случае f независимых скалярных параметров Г, для интегральной функции Лагранжа L, зависящей от этих параметров, справедлив интегральный принцип
L = |
J |
9? dV = шах. |
(6.52) |
|
к |
|
|
Ему соответствует вариационное условие |
|
||
6Ц Г„ Г2, |
... , |
Tf] = 6 j' S d V = О, |
(6,53) |
V
2 2 0 |
Глава Vf |
|
||
которому отвечают уравнения Эйлера — Лагранжа |
||||
д2 |
3 |
|
|
|
- 1 - дха д(дГі/дха) |
(»•= 1, 2, .... f). |
(6.54) |
||
<ЗГ; |
||||
|
а=1 |
|
|
|
В |
вариационной задаче |
(6.53) варьирование |
должно |
|
проводиться исключительно по параметрам Гг, градиен ты которых определяют силы. Такой способ варьирова ния становится понятным, если принять во внимание, что интегральный принцип, записанный в форме (6.52) или (6.53), тесно связан с силовым представлением принципа наименьшего рассеяния энергии, так как он относится к тем параметрам, градиентами которых яв ляются силы. При формулировке принципа наименьшего рассеяния энергии в представлении через силы [см. (4.21) —(4.25)] варьирование проводилось только по си лам при постоянных потоках. Уравнения баланса пока зывают, что входящие в них (а также в s) производные параметров по времени определяются потоками [см., на пример, (6.12)]. Таким образом, условия, согласно кото рым производные по времени не изменяются, по суще ству являются следствием постоянства потоков. Послед нее вариационное условие, однако, не является новым и независимым, потому что оно следует из принципа Онсагера, сформулированного в представлении через силы.
Прежде чем анализировать дальше теоретические во просы, докажем на примерах, что принцип, сформули рованный в основном с помощью соотношений (6.51) — (6.54), действительно является подлинным и точным ин тегральным принципом термодинамики. Он эквивален тен полной системе уравнений переноса, описывающих течение необратимых процессов во времени и простран стве.
§ 3. Вывод уравнений Фика для изотермической диффузии
Рассмотрим континуум из К компонентов в отсут ствие внешних сил. Пусть температура и давление одно родны, а градиенты химических потенциалов не обра щаются в нуль. Если на какое-то время пренебречь хи-
Интегральный принцип термодинамики |
221 |
мическими реакциями, то в системе, в которой выпол няются указанные выше условия, происходят только процессы диффузии, приводящие к локальному рассея
нию энергии
к- 1
Tad — 2 Jk • XI, |
(6.55) |
k=l |
|
где
х ; ^ - ѵ ( ^ - ц к) (/1=1, 2, ... , к - 1) (6.56)
— силы диффузии. Приведенные выше выражения, вы текающие из общего выражения (3.87), даны в энерге тическом представлении, которое удобно для решения изотермических проблем. Целесообразно ввести незави симые параметры
= |
(6 = 1 , 2, ... , К - 1), (6.57) |
при использовании которых потенциал рассеяния в энер гетическом представлении с учетом (4.15а) принимает вид
к - 1 |
к - I |
S |
L * i k X \ - X l = \ V L h V n - V n . (6.58) |
i, k = l |
i, k =l |
Теперь, следуя Верхашу [65, 79], запишем правую часть выражения (6.1) для принципа наименьшего рас сеяния энергии в энергетическом представлении, т. е. будем исходить из следующего вариационного условия:
6 J [ T o d - 4 r * 0,d ] d V |
= (6.59) |
V |
|
которое в силу (6.55) и (6.58) связано с принципом экстремума
к- 1 |
к-і |
1 |
dV — max. (6.60) |
Lft= l |
i, k=\ |
(Заметим, что, изменив знак, можно записать также и принцип минимума [65, 79], однако в интересах единства изложения сохраняется формулировка с максимумом,
2 22 Глава VI
которая использовалась до сих пор.) |
Применяя |
извест |
||||
ное из векторного анализа тождество |
|
|
|
|||
Ѵ -(/*П) = ПѴ-/* + /*-ѴП |
( k = \ , |
2, |
K - D |
(6.61) |
||
и преобразуя |
дивергенцию j |
V ■(/*Г*) dV |
при помощи |
|||
|
|
V |
|
|
получаем |
|
теоремы Гаусса в поверхностный интеграл, |
||||||
|
'К- 1 |
к- 1 |
|
|
|
|
- j |
V ПѴ./* + і £ L b v n - v n dV + |
|
||||
V с ft=i |
i, k=i |
|
|
|
|
|
|
|
+ |
J TUI ■dQ e= max. |
(6.62) |
||
Воспользуемся теперь уравнениями баланса компонен
тов |
V - /fe = 0 |
( £ = 1, 2, |
.. .. К), |
(6.63) |
pc* + |
||||
справедливыми |
для нашего |
случая, |
и потребуем |
(как |
мы это делали в случае теплопроводности), чтобы по токи /й, а также величины pèft [в силу (6.63)] были фик сированы при варьировании по параметрам Г£, т. е. чтобы выполнялись условия б/ft = 0 и Ö(pCft)= 0; кроме того, примем, что вдоль граничной поверхности системы варьирование не производится. При этих условиях экс тремальная задача (6.62) уже полностью определяется объемным интегралом
к-1 |
к- 1 |
ѴП d V : : шах. (6.64) |
[ 2 |
S |
|
L f t=l |
І, k=\ |
|
Здесь подынтегральное выражение есть не что иное, как плотность лагранжиана (в энергетическом представле нии) для диффузионной задачи
к-і |
к- |
1 |
|
і?Ь = р £ r kck - ~ |
J ] |
LlftVnѴП, |
(6.65) |
k=\ |
І, k=\ |
|
|
где параметр, по которому производится варьирование, указан в индексе. Используя плотность лагранжиана
Интегральный принцип термодинамики |
2 2 3 |
(6.65), получаем из уравнений Эйлера — Лагранжа (6.54) дифференциальные уравнения
Pck - |
S V |
• [L\kV (ц{ - |
цк)] = 0 |
(6.66) |
|
( |
k |
= 2,l , |
К |
-1), |
|
где мы в соответствии с (6.57) |
заменили параметры Г*% |
||||
химическими потенциалами. Уравнения (6.65) фактиче ски идентичны уравнениям Фика для многокомпонент ных изотермических систем.
Это легко показать. Действительно, функции
Ий= М Р і» Р2>•••> Рк-\) (6 = 1, 2, . .., К — 1) (6.67)
в изотермическом и изобарическом случаях удовлетво ряют тождествам
к- i д( |
_ |
. |
У ( ц , . - ^ ) = Ѵ |
^ |
- І х^ Ѵру ( / = 1 , к 2, ...1)- ., (6.68) |
/=1 |
1 |
|
.С помощью этих тождеств и используя следующие вы
ражения для коэффициентов диффузии: |
|
|||
D h = % |
~ |
( * . / = ! » |
2............К - 1 ) , |
(6.69) |
<=і |
1 |
|
|
|
получаем из (6.66) дифференциальные уравнения |
|
|||
p c ^ s V |
^ V p , ) |
( k = U |
2, ... , К - 1). |
(6.70) |
Если наша диффузионная система удовлетворяет усло вию, согласно которому плотность р постоянна и не за висит от времени, то из уравнения баланса масс (2.17) следует, что дивергенция скорости движения центра масс обращается в нуль, \-ѵ = 0. С другой стороны, если граничные условия для диффузионной системы та
ковы, что |
нормальная компонента скорости |
ѵ равна |
нулю на границе системы, то скорость ѵ равна |
нулю по |
|
всей системе. Этот так называемый «бесконвекционный» случай кажется весьма специальным с теоретической точки зрения, однако он очень часто встречается на
2 2 4 |
Глава VI |
практике. Вместо (6.70) (принимая во внимание пре образование коэффициентов диффузии) можно записать уравнения
- i f = Е Ѵ ' (d ;,Vc/) |
2........./С-1), |
(6.71) |
/=I |
|
|
которые представляют собой известную полную систему диффузионных уравнений Фика [3].
Для дальнейших рассуждений несколько преобра зуем плотность лагранжиана (6.65). Записывая в нашем изотермическом и изобарическом случаях уравнение
Гиббса — Дюгема (3.28а) |
в виде |
|
|
|
|
|
||||
|
к - I |
|
/ |
|
к - 1 |
\ |
|
|
|
|
|
2 |
ck diik + |
1 - |
2 |
ck)diiK = 0 |
|
(6.72) |
|||
|
ft=i |
|
\ |
|
ft=i |
/ |
|
|
|
|
и затем с помощью (6.57) |
переходя |
к параметрам П, |
||||||||
получаем уравнение |
|
|
|
|
|
|
|
|
||
|
|
к- |
1 |
|
|
|
|
|
|
|
|
|
2 |
ck dY\ = |
|
|
|
|
(6.73) |
||
где обе |
части |
являются |
полными |
дифференциалами. |
||||||
Воспользуемся соотношением |
|
|
|
|
|
|||||
|
К - 1 |
(Зс |
|
( Ä = l , 2, |
.. .. |
К - |
1) (6.74) |
|||
|
/=і |
(ЭГ |
|
|||||||
|
|
|
|
|
|
|
|
|
||
которое |
вытекает |
из |
выражения |
для |
функции |
|||||
ck ~ ck(y\> П> •••> rj_,), |
а также соотношениями сим |
|||||||||
метрии для смешанных вторых производных |
|
|
||||||||
|
д*к |
дсі |
|
dck |
|
д*»к |
|
|
(6.75) |
|
§ik — <эг; дгі |
дТІ |
|
дг: |
|
|
|
Ski |
|||
|
|
|
|
|
||||||
|
|
(i, k = l , |
2, |
... , |
ТС- |
1). |
|
|
|
|
С помощью (6.74) и (6.75) запишем плотность лагран жиана (6.65) в общем виде
К-1 |
К-1 |
S4 .-P £ |
£ Ч,ѵг; •ѵг;. <в.7в) |
і, k=\ |
i, k=l |
Интегральный принцип термодинамики |
22 5 |
Из определения параметров П (6.57) и из правила варьирования, которое мы использовали, следует, что эта плотность лагранжиана зависит только от химиче ских потенциалов и их градиентов, т. е. не содержит производных по времени от Г* в качестве независимых переменных. С помощью плотности лагранжиана (6.76) получаются уравнения Эйлера — Лагранжа
p 2 , gifcr; + |
s ' v . ( L ; Ävr;) = o ( / - 1, 2 ,..., * - і ) , (6.77) |
k = I |
fe= 1 |
которые являются общими уравнениями для многоком понентной диффузии. Эти уравнения очень просто за писать с помощью плотностей или произвольных кон центраций (например, мольных долей).
Все сказанное выше можно обобщить, если еще учесть существование химических реакций между ком понентами. Для этого необходимо лишь определить в энергетическом представлении соответствующий потен циал рассеяния. Он аналогичен последним двум членам в подынтегральном выражении (5.66). Таким образом, потенциал рассеяния имеет вид
/ С -1 |
R |
/С -1 |
|
|
ѵ г : |
ѵ ѵ г ; г ; . |
(6.78) |
і , k = \ |
I. r—I |
1, fc=! |
|
С его помощью, |
используя (6.51) и (6.76), получаем для |
||
плотности лагранжиана |
|
|
|
г-г— р 2 |
2 Ѵ Щ ѵ г ; - |
|
|
і , ft=l |
i, k = l |
|
|
|
-4 2 L i r 2 |
(6-79) |
|
|
j, r= l |
i, k = l |
|
Производя дифференцирование, указанное в (6.54),при ходим к уравнениям Эйлера — Лагранжа
/С-1 |
2 |
V •(/.;,№ •) = |
|
||
р S г „ г ; - |
|
||||
k=1 |
f t = l |
|
' |
|
|
|
R |
|
/ / C |
- l |
(6.80) |
|
= 2 |
|
L)r[ 2 |
vkr |
|
|
j , r—1 |
, r \ k = |
l Rr |
|
|
8 Зак. 787
226 |
Г лава VI |
Эти уравнения описывают эволюцию в пространстве и времени многокомпонентных изотермических диффу зионных и реагирующих систем, т. е. характеризуют пе ренос компонентов, обусловленный диффузией, а также происходящие в системе химические реакции.
Заметим, что из уравнения (6.80) можно получить стационарные уравнения, представляющие собой част ный случай уравнения (5.696). Хотя в гл. V использова лось энтропийное представление, а здесь— энергети ческое, различие между ними (если и в гл. V ограничиться изотермическим случаем) несущественно. Действительно, параметры Г/4 и П, определяемые соот ветственно соотношениями (5.65) и (6.57), отличаются только постоянным множителем Т. Поэтому плотности лагранжиана в соответствующем энтропийном представ лении отличаются от выведенного раньше множителем Г-1, т. е. справедливо соотношение
2 г = Т~'2т*. |
(6.81) |
§4. Вывод общего уравнения движения гидродинамики
Вкачестве подтверждения интегрального принципа выведем теперь общее уравнение движения вязкого по тока. Особенность этого примера состоит в том, что при вязком течении термодинамическая сила определяется градиентом скорости \ѵ , следовательно, варьирование должно проводиться по величине ѵа (а = 1, 2, 3), кото рая, согласно Казимиру, является типичным «ß»-napa- метром.
Общеизвестное уравнение Навье — Стокса для вяз кого течения было впервые выведено из интегрального принципа Верхашем [65, 79]. Выводом обобщенной фор мы уравнения Навье — Стокса мы обязаны Бэрэцзу [80], который, принимая во внимание асимметрическую часть тензора гидродинамического давления, получил более
общую форму уравнения Навье — Стокса, включающую и член, описывающий вращательную вязкость. Ниже, исходя из интегрального принципа термодинамики, вы водится наиболее общая форма гидродинамического уравнения движения. Сначала, однако, мы приведем
