Неравновесная термодинамика

Неравновесная термодинамика — раздел термодинамики, изучающий системы вне состояния термодинамического равновесия и необратимые процессы. Возникновение этой области знания связано главным образом с тем, что подавляющее большинство встречающихся в природе систем находятся вдали от термодинамического равновесия.

Содержание

История

Необходимость в создании новой теории возникла в первой половине двадцатого века. Пионером в этом направлении стал Ларс Онзагер, в 1931 году опубликовавший две работы, посвященные неравновесной термодинамике.[1][2] В дальнейшем значительный вклад в развитие неравновесной термодинамики внесли Эккарт[3], Майкснер и Райк[4], Д. Н. Зубарев[5], Пригожин[6], Де Гроот и Мазур[7], Гуров К. П. и другие. Следует отметить, что теория неравновесных систем активно развивается и в настоящее время.

Классическая формулировка неравновесной термодинамики

Основные положения

Классическая неравновесная термодинамика основана на фундаментальном предположении о локальном равновесии. Концепция локального равновесия заключается в том, что равновесные термодинамические соотношения справедливы для термодинамических переменных, определенных в элементарном объеме, то есть рассматриваемая система может быть мысленно разделена в пространстве на множество элементарных ячеек, достаточно больших, чтобы рассматривать их как макроскопические системы, но в то же время достаточно малых для того, чтобы состояние каждой из них было близко к состоянию равновесия. Данное предположение справедливо для очень широкого класса физических систем, что и определяет успех классической формулировки неравновесной термодинамики.

Концепция локального равновесия подразумевает, что все экстенсивные переменные (энтропия, внутренняя энергия, массовая доля компонента k{displaystyle k}

 ) заменяются своими плотностями:

s=s(x,t),u=u(x,t),ck=ck(x,t).{displaystyle s=s(mathbf {x} ,t),;;;;u=u(mathbf {x} ,t),;;;;c_{k}=c_{k}(mathbf {x} ,t).} 

В то же время все интенсивные переменные, такие как температура, давление и химический потенциал должны быть заменены соответствующими функциями координат и времени:

T=T(x,t),p=p(x,t),μk=μk(x,t),{displaystyle T=T(mathbf {x} ,t),;;;;p=p(mathbf {x} ,t),;;;;mu _{k}=mu _{k}(mathbf {x} ,t),} 

при этом они определяются так же, как и в равновесном случае, т.е. T−1=∂s∂u,T−1p=∂s∂v,T−1μk=∂s∂ck{displaystyle T^{-1}={frac {partial s}{partial u}},;T^{-1}p={frac {partial s}{partial v}},;T^{-1}mu _{k}={frac {partial s}{partial c_{k}}}}

 .

Далее, посредством введенных выше функций переписываются законы и соотношения из равновесной термодинамики в локальной форме. Первое начало (закон сохранения энергии):

∂e∂t+∇⋅Je=0{displaystyle {frac {partial e}{partial t}}+nabla cdot {boldsymbol {J}}^{e}=0} , e{displaystyle e}  — сумма плотностей кинетической и внутренней энергий, Je{displaystyle {boldsymbol {J}}^{e}}  — поток энергии.

Второе начало:

производство энтропии в каждой части системы, вызванное необратимыми процессами неотрицательно, то есть ∂s∂t+∇⋅Js=σ,σ⩾0{displaystyle {frac {partial s}{partial t}}+nabla cdot {boldsymbol {J}}^{s}=sigma ,;sigma geqslant 0} .

Важную роль в классической неравновесной термодинамике играет локальная форма уравнения Гиббса—Дюгема:

Tds=du+pdv−∑kμkdck{displaystyle Tds=du+pdv-sum _{k}mu _{k}dc_{k}} 

Переписав на последнем соотношении с учетом локальной формы закона сохранения энергии, массы, и сравнив с локальной формой второго начала, нетрудно получить следующий вид для производства энтропии:

σ=q⋅∇(1T)−∑kJk⋅∇(μkT)−1Tpν∇⋅ν−1TPν0:V0+ρT∑lAlξ˙l+1Tε⋅i.{displaystyle sigma ={boldsymbol {q}}cdot nabla left({frac {1}{T}}right)-sum _{k}{boldsymbol {J}}_{k}cdot nabla left({frac {mu _{k}}{T}}right)-{frac {1}{T}}p^{nu }nabla cdot {boldsymbol {nu }}-{frac {1}{T}}{overset {0}{mathbf {P} ^{nu }}}:{overset {0}{mathbf {V} }}+{frac {rho }{T}}sum _{l}A_{l}{dot {xi }}_{l}+{frac {1}{T}}{boldsymbol {varepsilon }}cdot {boldsymbol {i}}.} 

Здесь:

  • q{displaystyle {boldsymbol {q}}}  — поток теплоты,
  • ν=1ρ∑kρkνk{displaystyle {boldsymbol {nu }}={frac {1}{rho }}sum _{k}rho _{k}{boldsymbol {nu }}_{k}}  — скорость центра масс,
  • Jk=ρk(νk−ν){displaystyle {boldsymbol {J}}_{k}=rho _{k}({boldsymbol {nu }}_{k}-{boldsymbol {nu }})}  — поток диффузии,
  • тензор вязких напряжений разложен следующим образом: P=pU+pνU+Pν0{displaystyle mathbf {P} =pmathbf {U} +p^{nu }mathbf {U} +{overset {0}{mathbf {P} ^{nu }}}} , где тензор вязкого давления Pν{displaystyle mathbf {P} ^{nu }}  разложен на объемное вязкое давление pν{displaystyle p^{nu }}  и девиатор с нулевым следом Pν0{displaystyle {overset {0}{mathbf {P} ^{nu }}}} ,
  • аналогично, тензор скоростей деформации может быть разложен следующим образом: V=13(∇⋅ν)U+Vν0{displaystyle mathbf {V} ={frac {1}{3}}(nabla cdot {boldsymbol {nu }})mathbf {U} +{overset {0}{mathbf {V} ^{nu }}}} ,
  • двоеточие :{displaystyle :}  — двойное скалярное произведение тензоров,
  • Al{displaystyle A_{l}}  — химическое сродство реакции l{displaystyle l} , ξl{displaystyle xi _{l}}  — соответствующая степень полноты реакции,
  • ε=E+ν×B{displaystyle {boldsymbol {varepsilon }}={boldsymbol {E}}+{boldsymbol {nu }}times {boldsymbol {B}}}  — электрическое поле в системе координат, движущейся со скоростью ν{displaystyle {boldsymbol {nu }}} , i=∑kzkJk{displaystyle {boldsymbol {i}}=sum _{k}z_{k}{boldsymbol {J}}_{k}}  — ток проводимости.

Потоки и силы

В рамках классической неравновесной термодинамики описание необратимых процессов происходит при помощи термодинамических сил и термодинамических потоков. Основанием для введения данных величин является то, что через них производство энтропии выражается в простой форме. Дадим явные выражения для различных сил и потоков. Из приведенного выше выражения для производства энтропии видно, что σ{displaystyle sigma }

  представляет собой билинейную форму:

σ=∑αJαXα{displaystyle sigma =sum _{alpha }J_{alpha }X_{alpha }} ,

где Jα{displaystyle J_{alpha }}

  — термодинамический поток, Xα{displaystyle X_{alpha }}  — термодинамическая сила. Следует особо подчеркнуть произвольность разделения на термодинамические потоки и силы. Например, множитель 1T{displaystyle {frac {1}{T}}}  можно отнести не к силе, а к потоку. Силы и потоки можно даже поменять местами, однако всё же естественно считать, что термодинамические силы порождают термодинамические потоки, как градиент температуры порождает поток теплоты. Пример разделения сил и потоков показан в таблице:

Сила Xα{displaystyle X_{alpha }}  ∇1T{displaystyle nabla {frac {1}{T}}}  −∇μkT{displaystyle -nabla {frac {mu _{k}}{T}}}  −1T∇⋅ν{displaystyle -{frac {1}{T}}nabla cdot {boldsymbol {nu }}}  −1TV0{displaystyle -{frac {1}{T}}{overset {0}{mathbf {V} }}}  −1TAl{displaystyle -{frac {1}{T}}A_{l}}  −1Tε{displaystyle -{frac {1}{T}}{boldsymbol {varepsilon }}} 
Поток Jα{displaystyle J_{alpha }}  q{displaystyle {boldsymbol {q}}}  Jk{displaystyle {boldsymbol {J}}_{k}}  pν{displaystyle p^{nu }}  Pν0{displaystyle {overset {0}{mathbf {P} ^{nu }}}}  ρξ˙{displaystyle rho {dot {xi }}}  i{displaystyle {boldsymbol {i}}} 

Как видно, потоки и силы могут быть не только скалярами, но также векторами и тензорами.

Линейные материальные уравнения

Основная статья: Линейная неравновесная термодинамика

Потоки являются неизвестными величинами, в отличие от сил, которые представляют собой функции от переменных состояния и/или их градиентов. Экспериментально установлено, что потоки и силы связаны друг с другом, причем заданный поток зависит не только от своей силы, но может зависеть также от других термодинамических сил и от переменных состояния:

Jα=Jα(X1,…Xα;T,p,ck).{displaystyle J_{alpha }=J_{alpha }(X_{1},dots X_{alpha };,T,p,c_{k}).} 

Соотношения такого вида между потоками и силами называются феноменологическими соотношениями или материальными уравнениями. Они в совокупности с уравнениями баланса массы, импульса и энергии представляют замкнутую систему уравнений, которая может быть решена при заданных начальных и граничных условиях. Так как в положении термодинамического равновесия силы и потоки обращаются в нуль, то разложение материального уравнения вблизи положения равновесия принимает следующий вид:

Jα=∑βLαβXβ,Lαβ=∂Jα∂Xβ.{displaystyle J_{alpha }=sum _{beta }L_{alpha beta }X_{beta },;;;;L_{alpha beta }={frac {partial J_{alpha }}{partial X_{beta }}}.} 

Величины Lαβ{displaystyle L_{alpha beta }}

  называются феноменологическими коэффициентами и в общем случае зависят от переменных состояния T{displaystyle T} , p{displaystyle p}  и ck{displaystyle c_{k}} . Важно отдавать себе отчет в том, что, например, такая сила, как ∇1T{displaystyle nabla {frac {1}{T}}}  способна вызывать не только поток теплоты q{displaystyle {boldsymbol {q}}} , но электрический ток i{displaystyle {boldsymbol {i}}} . На феноменологические коэффициенты накладывается ряд ограничений, подробнее о них изложено в соответствующей статье.

Другим важным результатом, полученном в рамках линейной неравновесной термодинамики является теорема о минимуме производства энтропии:

В линейном режиме полное производство энтропии в системе, подверженной потоку энергии и вещества, в неравновесном стационарном состоянии достигает минимального значения.

Так же в это случае (линейный режим, стационарное состояние) показано, что потоки с собственными нулевыми силами равны нулю. Таким образом, например, при наличии постоянного градиента температуры, но при отсутствии поддерживаемого градиента концентрации система приходит к состоянию с постоянным потоком тепла, но с отсутствием потока вещества.

Системы вне локального равновесия

Несмотря на успехи классического подхода, у него есть существенный недостаток — он основывается на предположении о локальном равновесии, что может оказаться слишком грубым допущением для достаточно обширного класса систем и процессов. Примеры включают в себя такие явления, как распространение ультразвука в газах, суспензии, растворы полимеров, гидродинамика фононов, ударные волны, разреженные газы и т.д. С целью преодоления этого недостатка был создан целый ряд различных подходов.

Рациональная термодинамика

Расширенная необратимая термодинамика

Гамильтоновы формулировки неравновесной термодинамики

См. также

Неравновесный статистический оператор

Примечания

  1. L. Onsager, Phys. Rev. 37 (1931) 405
  2. L. Onsager, Phys. Rev. 38 (1931) 2265
  3. C. Eckart, Phys. Rev. 58 (1940) 267, 269, 919
  4. J. Meixner and H. Reik, Thermodynamik der Irreversiblen Prozesse (Handbuch der Physik III/2), (S. Flugge, ed.), Springer,Berlin, 1959.
  5. D. N. Zubarev, Double-time Green functions in statistical physics, Sov. Phys. Uspekhi, 1960, 3(3), 320—345.
  6. I. Prigogine, Introduction to Thermodynamics of Irreversible Processes, Interscience, New York, 1961.
  7. S.R. de Groot and P. Mazur, Non-equlibrium Thermodynamics, North-Holland, Amsterdam, 1962.

Литература