Главная/Уроки/Урок 4

Первое начало термодинамики

Работа, теплота, внутренняя энергия и первое начало.

Внутренняя энергияПервое началоРабота сил давленияТеплопередачаФункции состоянияПолный дифференциалКвазистатические процессыАдиабатический процессЭкстенсивность

В предыдущей лекции были введены термодинамические системы, равновесные состояния и функции состояния. Теперь мы ответим на центральный вопрос: как включить энергию системы в это макроскопическое описание?

В некоторых курсах этот вопрос решается слишком быстро: первое начало представляют как простое сохранение энергии, а затем сразу приводят формулу энергетического баланса без дальнейших обоснований. Мы встретим эти формулы в разделах 3 и 4.3, и спешащий читатель действительно может перейти прямо к ним (см. также краткие итоги в разделе 7).

Однако построение, приводящее к этим формулам, в действительности тонко. Прежде всего, первое начало утверждает больше, чем сохранение энергии, которое (в наши дни) кажется очевидным с микроскопической точки зрения. Мы покажем, что это начало на самом деле заполняет логический разрыв между существованием энергии UmicroU_{\mathrm{micro}}, определённой механически на пространстве примерно с 6N6N координатами, и существованием функции состояния энергии UthU_{\mathrm{th}}, зависящей лишь от нескольких макроскопических переменных.

Затем мы покажем, что аддитивность энергии является дополнительным предположением, необходимым для перехода от сохранения энергии изолированной системы к балансу энергии на границе замкнутой системы. Именно этот последний шаг позволяет наконец уточнить эпистемический и логический статус теплоты: она определяется как всякий перенос энергии, остающийся в этом балансе после учёта макроскопической работы, совершённой над замкнутой системой: Q  =def  ΔUWQ \equiv \Delta U - W.

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

1. Микроскопическая внутренняя энергия

Рассмотрим классическую, неквантовую и нерелятивистскую систему в инерциальной системе отсчёта1. Рассмотрим NN точечных частиц с массами mim_i, положениями ri\vec r_i и скоростями vi\vec v_i. В каждый момент времени её микроскопическое состояние полностью определяется набором

Note 1 : Эти случаи требуют особого подхода, выходящего за рамки данного курса.
(r1,...,rN,p1,...,pN)Γ,гдеpi=mivi,\boxed{ (\vec r_1,...,\vec r_N, \vec p_1,...,\vec p_N) \in\Gamma, \, \text{где} \,\,\, \vec p_i=m_i\vec v_i, }

где Γ\Gamma обозначает классическое фазовое пространство системы, имеющее в этой простой модели размерность dimΓ=6N\dim\Gamma=6N. Предполагается, что все внутренние силы, например силы притяжения или отталкивания между молекулами, консервативны. Им можно сопоставить полную потенциальную энергию Epint(r1,...,rN),гдеFiint=riEpintE_p^{\mathrm{int}}(\vec r_1,...,\vec r_N), \, \text{где} \,\,\, \vec F_i^{\,\mathrm{int}} =-\nabla_{\vec r_i}E_p^{\mathrm{int}}. Система может также находиться под действием внешних консервативных сил, например силы тяжести. Им точно так же сопоставляется внешняя потенциальная энергия Epext(r1,...,rN)E_p^{\mathrm{ext}}(\vec r_1,...,\vec r_N).

Система не обязательно покоится в выбранной системе отсчёта. В каждый момент времени можно определить её центр масс и его скорость V\vec V, вообще говоря зависящую от времени. Вводя полную массу системы и скорости относительно центра масс

M=i=1Nmi,V=1Mi=1Nmivi,ui=viV,M=\sum_{i=1}^{N} m_i, \qquad \vec V=\frac{1}{M}\sum_{i=1}^{N} m_i\vec v_i, \qquad \vec u_i=\vec v_i-\vec V,

полную кинетическую энергию можно записать в виде

Ec=12i=1Nmivi2=12i=1Nmi(V+ui)2.E_c=\frac12 \sum_{i=1}^{N}m_i \vec {v_i}^2 = \frac12 \sum_{i=1}^{N}m_i (\vec V + \vec u_i)^2.

Раскрывая квадрат, замечаем, что перекрёстные члены обращаются в нуль, поскольку imiui=0\sum_i m_i\vec u_i=\vec 0. Этот результат известен в классической механике как теорема Кёнига:

Ec=12MV2+12i=1Nmiui2+Vi=1NmiuiimiviMV=0E_c =\frac12MV^2 +\frac12 \sum_{i=1}^{N}m_i \vec{u_i}^2 +\vec V\mathbin{\cdot} \underbrace{\sum_{i=1}^{N} m_i\vec u_i}_{\sum_i m_i\vec v_i-M\vec V=\vec0}

Отсюда:

Ec=Ecmacro+Ecmicro\boxed{E_c = E_c^{\mathrm{macro}} + E_c^{\mathrm{micro}}}

(1)

Полная кинетическая энергия, таким образом, всегда разлагается2 на макроскопическую кинетическую энергию центра масс и на энергию остаточных относительных, «беспорядочных» движений, которые и составляют то, что называют тепловым движением.

Note 2 : Поступательное движение центра масс не обязательно исчерпывает макроскопическую кинетическую энергию. Система может обладать упорядоченным вращением как целое. В этом случае можно записать разложение vi=V+ω×(riR)+ui,\vec v_i = \vec V +\vec\omega\times(\vec r_i-\vec R) +\vec u_i', где R\vec R — положение центра масс, а ui\vec u_i' обозначает теперь остаточное движение после вычитания поступательного движения и вращения системы как целого. Выбирая ω\vec\omega так, чтобы это остаточное движение обладало нулевым суммарным моментом импульса относительно центра масс, получаем Ec=12MV2+12ωICMω+12i=1Nmiui2,E_c = \frac12 MV^2 + \frac12\vec\omega\cdot \mathbf I_{\mathrm{CM}}\vec\omega + \frac12\sum_{i=1}^{N}m_i {u_i'}^2, где ICM\mathbf I_{\mathrm{CM}} — тензор инерции относительно центра масс. Разделение макроскопической и микроскопической кинетических энергий, следовательно, остаётся возможным.

Собирая различные слагаемые, полную механическую энергию можно окончательно записать так:

Etot=Ecmacro+Epext+[Ecmicro+Epint(r1,...,rN)+Eintautres]внутренняя энергия.E_{\mathrm{tot}} = E_c^{\mathrm{macro}} +E_p^{\mathrm{ext}} +\underbrace{\left[ E_c^{\mathrm{micro}} +E_p^{\mathrm{int}}(\vec r_1,...,\vec r_N) +E_{\mathrm{int}}^{\mathrm{autres}} \right]}_{\text{внутренняя энергия}}.

Замечание: член EintautresE_{\mathrm{int}}^{\mathrm{autres}} объединяет вклады внутренних степеней свободы, которые эта минимальная точечная модель не описывает, например вращательных или колебательных степеней свободы молекул. Если их моделировать явно, фазовое пространство Γ\Gamma должно быть соответственно расширено.

Определение 1 (Микроскопическая внутренняя энергия)
Микроскопическую внутреннюю энергию мы будем обозначать Umicro=Ecmicro+Epint+EintautresU_{\mathrm{micro}}= E_c^{\mathrm{micro}} +E_p^{\mathrm{int}} +E_{\mathrm{int}}^{\mathrm{autres}}. Это функция, определённая на фазовом пространстве, Umicro:ΓR,U_{\mathrm{micro}}:\Gamma\longrightarrow\mathbb R,

которая зависит, вообще говоря, от огромного числа переменных и удовлетворяет соотношению

Etot=Ecmacro+Epext+Umicro.E_{\mathrm{tot}} =E_c^{\mathrm{macro}}+E_p^{\mathrm{ext}}+U_{\mathrm{micro}}.

(2)

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

Etot=Umicro.E_{\mathrm{tot}}=U_{\mathrm{micro}}.

Теперь мы введём функцию состояния — термодинамическую энергию, временно обозначаемую UthU_{\mathrm{th}}, — и поймём, почему её не следует a priori смешивать с UmicroU_{\mathrm{micro}}.

2. Термодинамическая внутренняя энергия

Возьмём в качестве системы кастрюлю с водой, поставленную на нагревательную плиту. Мы остаёмся в простых рамках, описанных выше: кастрюля макроскопически неподвижна, и её внешняя потенциальная энергия не изменяется. Поэтому можно записать Etot=UmicroE_{\mathrm{tot}}=U_{\mathrm{micro}} с точностью до постоянной.

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

ΔEtot=ΔUmicro.\Delta E_{\mathrm{tot}}=\Delta U_{\mathrm{micro}}.

Термодинамика стремится описать в точности то же самое изменение. Поэтому хотелось бы построить величину UthU_{\mathrm{th}} такую, что

ΔUth=ΔUmicro.\Delta U_{\mathrm{th}}=\Delta U_{\mathrm{micro}}.

Решение кажется очевидным. Не достаточно ли отождествить термодинамическую внутреннюю энергию с микроскопической внутренней энергией?

На самом деле такое отождествление было бы математически неверным, поскольку у этих двух функций разные области определения. UmicroU_{\mathrm{micro}} — функция, определённая на фазовом пространстве и зависящая по меньшей мере от 6N6N переменных, тогда как UthU_{\mathrm{th}} зависит лишь от небольшого числа dd независимых переменных, описывающих равновесное состояние X=(x1,...,xd)EX=(x^1,...,x^d)\in\mathcal E (напомним, что xix^i — макроскопические переменные состояния, такие как температура, давление или объём).

Тем не менее эти две функции должны быть связаны так, чтобы Uth(X)U_{\mathrm{th}}(X) представляла на макроскопическом уровне внутреннюю энергию, соответствующую равновесному состоянию XX. Нужно объяснить, как работает это понижение размерности. Вопрос этот глубок и выходит за рамки вводного курса. Микроскопическое описание необходимо, если мы хотим вывести этот переход между двумя масштабами. Подраздел 6.2, выходящий за рамки программы, излагает в общих чертах то, как статистическая физика решает эту задачу.

Но со строго термодинамической точки зрения обосновать это микроскопическое отождествление невозможно, и существование UthU_{\mathrm{th}} как функции лишь нескольких макроскопических переменных приходится постулировать. Это и будет ролью первого начала термодинамики, которое, следовательно, утверждает не только сохранение энергии, но также и прежде всего то, что эта энергия может быть представлена в равновесии функцией состояния, определённой на макроскопическом пространстве E\mathcal E.

Как только выбраны dd независимых переменных состояния xix^i, например (TT, VV, NN) для идеального газа, функция U(T,V,N)U(T,V,N) задаёт скалярное поле на пространстве состояний; в расширенном пространстве координат U,T,V,NU,T,V,N её график представляет собой гиперповерхность, называемую поверхностью равновесия, роль которой в дальнейшем изложении будет решающей, ср. рис. 1.

Гиперповерхность энергии, определённая на пространстве равновесных состояний.
Рисунок 1. Гиперповерхность энергии, определённая на пространстве равновесных состояний.

Отметим в заключение непосредственное следствие, хорошо иллюстрирующее различие между UthU_{\mathrm{th}} и UmicroU_{\mathrm{micro}}: в развиваемой здесь равновесной термодинамике термодинамическая внутренняя энергия, вообще говоря, не определена вне равновесных состояний. Например, при резком процессе в газе возможно и даже обычно, что ему уже нельзя приписать единую температуру или единое давление. В этом случае ему нельзя непосредственно сопоставить и термодинамическую внутреннюю энергию UthU_{\mathrm{th}}, поскольку она определена лишь на пространстве равновесных состояний E\mathcal E. Проблема не в том, что энергии больше не существует, — ведь микроскопическая энергия остаётся прекрасно определённой в каждый момент времени, — а в том, что система просто больше не представляется точкой XEX\in\mathcal E.

После всего сказанного, начиная со следующего раздела, механическая функция UmicroU_{\mathrm{micro}} больше не будет появляться непосредственно: мы будем просто писать UU для термодинамической внутренней энергии UthU_{\mathrm{th}}, подразумевая, что она определена лишь на пространстве равновесных состояний.

Замечание 1 (Расширения вблизи равновесия)
Отметим, что для систем, достаточно близких к равновесию, термодинамическое описание можно расширить, увеличив пространство состояний. Например, в уже встречавшемся примере металлического стержня, помещённого между двумя источниками температуры, можно ввести локальное стационарное поле температуры T(x)T(x) вместо единой температуры и определить локальную плотность внутренней энергии. Полная энергия становится тогда функционалом этих полей. Мы не будем здесь развивать эту так называемую неравновесную термодинамику вблизи равновесия.

3. Первое начало термодинамики

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

Первое начало термодинамики для замкнутой системы
Для всякой замкнутой термодинамической системы существует функция состояния U:ER,U:\mathcal E\longrightarrow\mathbb R,

называемая внутренней энергией, определённая лишь с точностью до аддитивной постоянной и по меньшей мере дважды дифференцируемая, такая, что для всякого равновесного состояния XEX\in\mathcal E

Etot=Ecmacro+Epext+U(X).E_{\mathrm{tot}} = E_c^{\mathrm{macro}} +E_p^{\mathrm{ext}} +U(X).
(3)

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

ΔEtot=0.\Delta E_{\mathrm{tot}}=0.

Гладкость функции UU — тоже нетривиальное предположение. На практике UU можно считать функцией класса C\mathcal C^\infty. Поскольку UU определена лишь с точностью до постоянной из-за потенциальных энергий, входящих в равенство 3, физический смысл имеют только её изменения. Для простоты мы ограничимся далее процессами, при которых макроскопические кинетическая и потенциальная энергии не изменяются. Тогда ΔEtot=ΔU\Delta E_{\mathrm{tot}}=\Delta U.

Рассмотрим теперь замкнутую систему SS, которая может обмениваться со своим окружением энергией, но не веществом. Полную систему

S=S+ExtS'=S+\mathrm{Ext}

всегда можно включить в изолированную систему. Первое начало налагает тогда условие

ΔES=0.\Delta E_{S'}=0.

При дополнительном предположении, что энергия аддитивна, то есть что энергия их взаимодействия пренебрежимо мала, имеем: ES=ES+EExtE_{S'}=E_S+E_{\mathrm{Ext}}, так что

ΔES=ΔEExt.\Delta E_S=-\Delta E_{\mathrm{Ext}}.

Предположение аддитивности позволяет тем самым выразить сохранение энергии в виде баланса между подсистемами: всякая энергия, полученная системой, теряется её окружением, и наоборот. Остаётся теперь уточнить, как эта энергия может передаваться через границу системы. Мы различим два способа передачи: работу, обозначаемую WW, которая обладает независимым определением, унаследованным от механики, и теплоту, обозначаемую QQ.

Энергетический баланс примет тогда вид

ΔES=ΔU=Q+W\boxed{ \Delta E_S = \Delta U = Q+W }
(4)

где WW обозначает работу, полученную системой, а QQ — полученную ею теплоту.

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

Теперь нам остаётся определить эти два способа передачи. Но прежде необходимо отметить весьма стандартное соглашение о знаках, которым мы будем пользоваться на протяжении всего курса.

Запомнить (Банковское соглашение о знаках)
Во всей этой книге перенос энергии считается положительным, когда он получен системой, и отрицательным, когда он отдан ею. Так, Q>0, W>0Q>0,\ W>0

отвечают энергии, полученной системой, тогда как

Q<0, W<0Q<0,\ W<0

отвечают энергии, отданной наружу.

4. Работа, теплота и энергетический баланс

4.1. Работа

Одно из преимуществ работы состоит в том, что она обладает определением, не зависящим от термодинамики и унаследованным от механики. Когда внешняя сила Fext\vec F_{\mathrm{ext}} действует на точку границы системы и эта точка смещается на drd\vec r, элементарная работа, совершённая этой силой и, значит, полученная системой, равна

δW=Fextdr.\delta W = \vec F_{\mathrm{ext}}\cdot d\vec r.

Для протяжённой системы нужно просуммировать вклады всех внешних сил, совершающих работу на её границе.

Случай, который чаще всего будет встречаться нам на практике, — это жидкость или газ в сосуде с подвижной стенкой, что приводит к вычислению работы сил давления. Рассмотрим газ, заключённый в цилиндре плоским поршнем сечения S\mathcal S, ср. рис. 2. Ось xx направлена наружу от газа. Когда поршень смещается на величину dxdx, объём изменяется на

dV=Sdx.dV=\mathcal S\,dx.

Обозначим через PextP_{\mathrm{ext}} внешнее давление, действующее на границу газа. Сила, полученная газом, направлена внутрь и равна

Fext=PextSex.\vec F_{\mathrm{ext}} = -P_{\mathrm{ext}}\mathcal S\,\vec e_x.

Элементарная работа, полученная газом, равна, следовательно,

δWpression=Fext(dxex)=PextSdx,\delta W_{\mathrm{pression}} = \vec F_{\mathrm{ext}}\cdot(dx\,\vec e_x) = -P_{\mathrm{ext}}\mathcal S\,dx,

то есть

δWpression=PextdV\boxed{ \delta W_{\mathrm{pression}}=-P_{\mathrm{ext}} dV }
(5)

Банковское правило позволяет сразу же проверить знак. При сжатии dV<0dV<0, значит δW>0\delta W>0: газ получает работу. При расширении dV>0dV>0, значит δW<0\delta W<0: газ совершает работу над своим окружением.

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

Работа, полученная газом при смещении поршня сечения S. При расширении dV>0, тогда как сила, создаваемая внешним давлением, противоположна смещению: полученная работа, следовательно, отрицательна.
Рисунок 2. Работа, полученная газом при смещении поршня сечения S\mathcal S. При расширении dV>0dV>0, тогда как сила, создаваемая внешним давлением, противоположна смещению: полученная работа, следовательно, отрицательна.

Интегрируя, получаем следующую формулу для конечного процесса ABA \to B, в ходе которого внешнее давление известно:

Wpression[AB]=ABPextdV\boxed{ W_{\mathrm{pression}}[A \to B] = -\int_{A}^B P_{\mathrm{ext}} dV }
(6)

Этот интеграл показывает, что для вычисления полной работы нужно знать значение PextP_{\mathrm{ext}} на всём протяжении процесса. Отсюда важнейшее следствие: работа, вообще говоря, не может быть определена по одним лишь начальному и конечному состояниям AA и BB: её выражение зависит от совершённого процесса.

Здесь важно отметить три частных случая.

  1. Если PextP_{\mathrm{ext}} постоянно, то Wpression[AB]=PextΔVW_{\mathrm{pression}}[A\to B] = -P_{\mathrm{ext}}\Delta V.
  2. Вдоль изохорического процесса dV=0dV=0, и работа сил давления равна нулю: δWpression=0\delta W_{\mathrm{pression}}=0 и Wpression[AB]=0W_{\mathrm{pression}}[A\to B]=0.
  3. При расширении в вакуум Pext=0P_{\mathrm{ext}}=0, поэтому работа сил давления также равна нулю.

Замечание 2 (Другие виды работы)
Работа сил давления — лишь один пример среди прочих. Для растягиваемой нити, если LL обозначает её длину, а FextF_{\mathrm{ext}} — полученную силу натяжения, δW=FextdL.\delta W=F_{\mathrm{ext}}\,dL.

Для поверхности раздела, площадь AA которой увеличивают, поверхностное натяжение γ\gamma приводит в обычных условиях к работе вида

δW=γdA.\delta W=\gamma\,dA.

Энергию можно передавать и без видимого механического перемещения. Например, перемещение заряда dqdq через разность потенциалов может приводить к электрической работе вида

δWeˊlec=Vextdq,\delta W_{\mathrm{élec}}=V_{\mathrm{ext}}\,dq,

где VextV_{\mathrm{ext}} — внешний электрический потенциал, часто обозначаемый в этом контексте Φext\Phi_{\mathrm{ext}}, чтобы не путать его с объёмом. Знак dqdq здесь также выбран в соответствии с банковским соглашением.

Подраздел 6.3 уточнит общий вид внешних работ.

4.2. Теплота

Работа не может быть единственным способом передачи энергии, поскольку нагревание газа, заключённого в жёстком сосуде, повышает его температуру и его энергию без совершения какой-либо внешней работы. Значит, энергия пересекла границу системы, не будучи переданной в форме работы. Этот второй способ передачи энергии называется теплотой, или тепловым переносом.

В принятом здесь построении (см. подраздел 6.1 об эквивалентном альтернативном построении) для физического процесса ABA \to B, связывающего два равновесных состояния AA и BB, после того как выявлены все полученные системой работы, полученная теплота определяется как

Q[AB]  =def  U(B)U(A)W[AB].\boxed{ Q[A \to B] \equiv U(B)-U(A)-W[A \to B]. }
(7)

так что мы вновь получаем объявленное выше соотношение

ΔU=Q+W\boxed{ \Delta U=Q+W }

Поскольку работа зависит от пути процесса ABA \to B, а изменение энергии ΔU=U(B)U(A)\Delta U = U(B)-U(A) от него не зависит, количество обменянной теплоты необходимым образом зависит от него тоже.

На языке термодинамики говорят, что внутренняя энергия UU является функцией состояния (по постулату), тогда как работа WW и теплота QQ ими не являются. Напомним, что это означает: то, что WW и QQ не могут быть функциями состояния системы в данный момент времени, означает, что система не обладает «некоторым количеством работы или теплоты», а скорее что WW и QQ суть переносы энергии на её границе в ходе произвольного процесса ABA \to B, — в точности то, чего мы хотели добиться ввиду различных опытов, подробно разобранных в лекции №2.

Отметим, что при циклическом процессе AAA \to A система возвращается в своё начальное состояние, поэтому

ΔUcycle=0.\Delta U_{\mathrm{cycle}}=0.

Первое начало налагает тогда условие

Qcycle+Wcycle=0.\boxed{Q_{\mathrm{cycle}}+W_{\mathrm{cycle}}=0.}

Циклическая машина, следовательно, не может бесконечно совершать работу, не получая равного количества энергии от своего окружения. Исторически это называют невозможностью вечного двигателя первого рода.

4.3. Дифференциальная форма первого начала

Соотношение

ΔU=Q+W\Delta U=Q+W

связывает два равновесных состояния AA и BB. Оно всегда верно в том смысле, что не предполагает, будто все промежуточные состояния сами могут быть представлены точками пространства равновесных состояний E\mathcal E. Если процесс резкий, система может временно покинуть поверхность равновесия, определённую на E\mathcal E, тогда как U(A)U(A) и U(B)U(B) остаются прекрасно определёнными.

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

γE\gamma\subset\mathcal E

составленным из бесконечно близких равновесных состояний и не покидающим поверхность равновесия. Поскольку UU — дифференцируемая функция состояния на E\mathcal E, её изменение между двумя бесконечно близкими состояниями является полным дифференциалом, обозначаемым dUdU.

Соответствующие элементарные переносы теплоты и работы обозначают тогда δQ\delta Q и δW\delta W. Первое начало принимает в этом случае вид

dU=δQ+δW\boxed{ dU=\delta Q+\delta W }
(8)

Различие в обозначениях между dUd U и δW\delta W или δQ\delta Q не является чисто косметическим. Оно указывает, что первое есть полный дифференциал, а два других — нет, что представляет собой математический перевод того факта, что изменение энергии не зависит от пройденного пути, тогда как QQ и WW от него зависят. Следующий раздел напомнит начинающему читателю математику, необходимую для правильного понимания этого важнейшего пункта.

В частном случае, когда единственная работа — это работа сил давления и когда PextP_{\mathrm{ext}} можно отождествить с давлением PP системы, это соотношение принимает вид

dU=δQPdV.dU=\delta Q-P\,dV.
Запомнить (Две формулировки первого начала)
Интегральная форма ΔU=Q+W\Delta U=Q+W

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

dU=δQ+δWdU=\delta Q+\delta W

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

5. Точные и неточные дифференциальные формы

Этот раздел собирает минимум математики, необходимый для придания точного смысла различию между dUdU, с одной стороны, и δQ\delta Q и δW\delta W — с другой. Читатель, которому это уже известно, может его пропустить.

5.1. Дифференциал функции

Пусть дана функция ff переменных x1,...,xdx^1,...,x^d, предполагаемая дифференцируемой. Её дифференциалом называется выражение

df=i=1dfxidxi,df=\sum_{i=1}^{d}\frac{\partial f}{\partial x^i}\,dx^i,

которое измеряет изменение ff при переходе от точки XX к соседней точке X+dXX+dX. Важно следующее: если следовать вдоль пути γ\gamma, идущего от точки AA к точке BB, то, суммируя все эти элементарные изменения, получаем

γdf=f(B)f(A).\int_\gamma df=f(B)-f(A).

В частности, вдоль замкнутого пути:

df=0.\oint df=0.

5.2. Дифференциальные формы

Рассмотрим теперь выражение того же вида, которое, однако, ни в коем случае не следует смешивать с дифференциалом функции. Пусть даны произвольные функции Ai(x1,...,xd)A_i(x^1,...,x^d). Определим тогда ω\omega как

ω=i=1dAidxi.\omega=\sum_{i=1}^{d}A_i\,dx^i.

Такой объект называют дифференциальной формой. Если задано выражение функций AiA_i, её без труда можно проинтегрировать вдоль пути γ\gamma. Но ничто не гарантирует существования функции ff, дифференциалом которой была бы ω\omega, то есть такой, что Ai=f/xiA_i=\partial f/\partial x^i для всех ii.

  • Если такая функция ff существует, форму называют точной и пишут ω=df\omega = d f. Интеграл γω=γdf=f(B)f(A)\int_\gamma\omega = \int_\gamma df = f(B)-f(A) зависит тогда только от концов пути.
  • В противном случае форму называют неточной и обозначают её тогда δω\delta\omega, а не ω\omega, чтобы напомнить, что она не является дифференциалом никакой функции. (Математики записали бы форму просто как ω\omega.)

5.3. Критерий Шварца

Простой критерий, позволяющий определить, является ли дифференциальная форма точной, — критерий Шварца. Ограничимся для простоты двумя переменными (обобщение непосредственно), положив

ω=A(x,y)dx+B(x,y)dy.\omega=A(x,y)\,dx+B(x,y)\,dy.

Если бы ω\omega была точной, то мы имели бы A=f/xA=\partial f/\partial x и B=f/yB=\partial f/\partial y, а значит

Ay=2fyx=2fxy=Bx,\frac{\partial A}{\partial y} =\frac{\partial^2 f}{\partial y\,\partial x} =\frac{\partial^2 f}{\partial x\,\partial y} =\frac{\partial B}{\partial x},

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

AyBx    ω не является точной.\boxed{ \frac{\partial A}{\partial y}\neq\frac{\partial B}{\partial x} \;\Longrightarrow\; \omega\ \text{не является точной.} }

С точностью до тонкостей, которые в термодинамике не встретятся, верно и обратное.

5.4. Пример: работа сил давления

Когда элементарные переносы могут быть выражены через переменные состояния вдоль квазистатического пути, δQ\delta Q и δW\delta W становятся дифференциальными формами на E\mathcal E, и можно проверить, точны они или нет.

В качестве примера рассмотрим nn молей идеального газа, испытывающего квазистатический процесс, при Pext=PP_{\mathrm{ext}} = P. Тогда δW=PdV\delta W=-P\,dV. Но так как P=nRT/VP=nRT/V, элементарная работа записывается в виде

δW=PdV=nRTVdV=0×dTnRTVdV,то естьA=0,B=nRTV,\delta W=-P dV = -\frac{nRT}{V}\,dV = 0\times dT-\frac{nRT}{V}\,dV, \qquad\text{то есть}\qquad A=0, \quad B=-\frac{nRT}{V},

при отождествлении x=Tx=T и y=Vy=V. Это действительно дифференциальная форма, определённая на пространстве состояний. Критерий Шварца даёт

AV=0,BT=nRV0.\frac{\partial A}{\partial V}=0, \qquad \frac{\partial B}{\partial T}=-\frac{nR}{V}\neq0.

Две смешанные производные не равны: элементарная работа действительно не является полным дифференциалом. Следовательно, в этом случае не существует никакой функции состояния W(T,V)W(T,V), изменением которой была бы δW\delta W, и полученная между двумя состояниями работа и в самом деле зависит от пройденного пути. То же рассуждение справедливо для δQ=dUδW\delta Q=dU-\delta W: так как dUdU точна, а δW\delta W — нет, их разность тоже не может быть точной.

Запомнить (почему δ\delta, а не dd)
dUdU — полный дифференциал: её изменение зависит только от начального и конечного состояний, а её интеграл по циклу равен нулю. δQ\delta Q и δW\delta W — неточные формы: они зависят от пройденного пути, и не существует ни функции QQ, ни функции WW состояния системы таких, что δW=dW\delta W = d W и δQ=dQ\delta Q = dQ.

6. Для дальнейшего изучения

Этот раздел собирает более продвинутые соображения, объявленные ранее. При первом чтении его можно пропустить.

6.1. Два возможных построения первого начала

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

Q=ΔUW.Q=\Delta U-W.

Отсюда вытекает замечательный факт. Если процесс адиабатический, то ΔU=W\Delta U = W, а это требует, чтобы адиабатическая работа WadW_{\mathrm{ad}}, в свою очередь, всегда не зависела от пройденного пути.

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

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

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

Полученная работа, в свою очередь, измеряется чисто механическим или электрическим способом. В этом всё эмпирическое содержание опытов Джоуля, проведённых между 1843 и 1850 годами и повторявшихся с тех пор со всё возрастающей точностью: масса mm, падающая с высоты hh, доставляет mghmgh без какого-либо теплового переноса. Все устройства, которые он придумал и осуществил в адиабатическом сосуде в указанном выше смысле, — лопастная мешалка, нагревательное сопротивление или сжатие газа, — привели его к одному и тому же выводу: одна и та же совершённая работа производит одно и то же изменение состояния, что и составляло для него «механический эквивалент теплоты», как мы объясняли в лекции №2. Джоуль был, таким образом, первым, кто показал (экспериментальное указание на то), что адиабатическая работа не зависит от пройденного пути.

Именно Каратеодори в 1909 году [2], а затем прежде всего Борн в 1921 году предложили возвести этот экспериментальный факт в ранг постулата. Если его принять, становится легко построить функцию состояния UU, полагая

U(B)U(A)=Wad(AB),U(B)-U(A)=W_{\mathrm{ad}}(A\to B),

а затем определить теплоту по разности для общих процессов, как мы это и сделали.

Такое построение воспроизведено во многих сочинениях, в частности у Пиппарда [3] и Каллена [4]. Оно эквивалентно нашему: постулировать UU, а затем вывести независимость адиабатической работы, — то же самое, что постулировать независимость этой работы и построить UU. Подробный рассказ об этой концептуальной эволюции, от Джоуля к Каратеодори и Борну, можно найти у Розенберга [5].

6.2. Микро- и макросостояния

Если не выбирать описанный выше операциональный путь, то раздел 2 выявил фундаментальный вопрос: как энергия, определённая на пространстве с 6N6N координатами, может сводиться к функции лишь нескольких макроскопических переменных? Статистическая физика строит этот переход между двумя масштабами. Вот его основные черты.

В статистической физике точку XX пространства равновесных состояний E\mathcal E называют макросостоянием. Одному и тому же макросостоянию отвечает, вообще говоря, гигантское число микросостояний γΓ\gamma\in\Gamma, совместимых с одними и теми же макроскопическими ограничениями.

Затем доказывается, что энергии этих микросостояний распределены чрезвычайно узко вокруг своего среднего значения при NN\to\infty, — это так называемый термодинамический предел. Иначе говоря, относительные флуктуации энергии стремятся к нулю:

σ(Umicro)Umicro0.\frac{\sigma(U_{\mathrm{micro}})} {\langle U_{\mathrm{micro}}\rangle}\longrightarrow0.

Тогда можно положить Uth(X)=UmicroXU_{\mathrm{th}}(X)= \langle U_{\mathrm{micro}} \rangle_X, где усреднение ведётся по всем микросостояниям, совместимым с данным макросостоянием. Упомянутая выше концентрация и есть то, что придаёт смысл этому определению. Система в каждый момент времени занимает единственное микросостояние, а не среднее: лишь потому, что почти все они несут одну и ту же энергию, можно говорить об энергии макросостояния и о том, что два одинаковых приготовления дают одно и то же измерение.

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

Следующий рисунок подытоживает эту процедуру.

Одно и то же макросостояние X совместимо с необъятным числом микросостояний. Их энергии, в нашем масштабе, чрезвычайно узко сосредоточены вокруг одного макроскопического значения, что и позволяет определить U_ th(X).
Рисунок 3. Одно и то же макросостояние XX совместимо с необъятным числом микросостояний. Их энергии, в нашем масштабе, чрезвычайно узко сосредоточены вокруг одного макроскопического значения, что и позволяет определить Uth(X)U_{\mathrm{th}}(X).

6.3. Обобщённые работы

Различные виды работы, встреченные выше, обладают общей структурой. Мы видели выражения

δW=PextdV\delta W=-P_{\mathrm{ext}}\,dV

для работы давления,

δW=FextdL\delta W=F_{\mathrm{ext}}\,dL

для растяжения нити, а также

δW=γextdA\delta W=\gamma_{\mathrm{ext}}\,dA

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

δW=iνiextdxi,\boxed{ \delta W=\sum_i \nu_i^{\mathrm{ext}}\,dx^i, }

где νiext\nu_i^{\mathrm{ext}} обозначает внешнюю обобщённую силу (обычно интенсивную), сопряжённую координате xix^i (обычно экстенсивной), со знаком, отвечающим нашему соглашению.

В случае, когда внешние обобщённые силы совпадают с соответствующими термодинамическими силами системы (пример: Pext=PP_{\mathrm{ext}}=P), первое начало можно записать в виде

dU=δQ+iνidxi.\boxed{ dU=\delta Q+\sum_i\nu_i\,dx^i.}

6.4. Экстенсивность внутренней энергии

В предыдущей лекции мы представили внутреннюю энергию как экстенсивную величину: умножение размера однородной системы на множитель λ\lambda при фиксированных интенсивных переменных умножает её энергию на λ\lambda,

UλU.U\longrightarrow\lambda U.

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

Уточним это. Рассмотрим NN составляющих, распределённых с постоянной плотностью в пространстве размерности DD, и предположим, что их потенциальная энергия взаимодействия ведёт себя как

Ep(r)1rα.E_p(r)\sim\frac{1}{r^\alpha}.

При фиксированной плотности ρ=N/V\rho=N/V типичный линейный размер системы растёт как

LN1/D.L\sim N^{1/D}.

Оценим тогда полную энергию взаимодействия в три приёма.

Подсчёт соседей.

Зафиксируем одну составляющую и спросим, сколько у неё других на расстоянии между rr и r+drr+dr. Это плотность, умноженная на объём соответствующего слоя:

dn(r)=ρSDrD1dr,dn(r)=\rho\,S_D\,r^{D-1}\,dr,

где SDS_D обозначает площадь единичной сферы в размерности DD, то есть 4π4\pi в трёх измерениях.

Суммирование по расстояниям.

Каждый из этих соседей вносит вклад Ep(r)rαE_p(r)\sim r^{-\alpha}. Энергия взаимодействия одной составляющей со всеми остальными равна, следовательно,

uρaLrD11rαdr=ρaLrD1αdr.u\sim\rho\int_a^L r^{D-1}\,\frac{1}{r^{\alpha}}\,dr =\rho\int_a^L r^{D-1-\alpha}\,dr .

Нижний предел aa — это минимальное расстояние сближения, ниже которого закон 1/rα1/r^\alpha перестаёт быть справедливым и которое не даёт интегралу расходиться при r0r\to0. Верхний предел LL — размер системы: за ним соседей нет.

Суммирование по составляющим.

Умножаем на число составляющих (деля на два, чтобы не считать каждую пару дважды):

UintN2uNρaLrD1αdr,U_{\mathrm{int}}\sim\frac{N}{2}\,u \sim N\rho \int_a^L r^{D-1-\alpha}\,dr,

где геометрический множитель SD/2S_D/2 поглощён оценкой порядка величины. Посмотрим тогда на значение интеграла:

aLrD1αdr=[rDαDα]aL(αD).\int_a^L r^{D-1-\alpha}\,dr =\left[\frac{r^{D-\alpha}}{D-\alpha}\right]_a^L \qquad(\alpha\neq D).

Если α>D\alpha>D, показатель отрицателен: интеграл сходится при LL\to\infty и определяется своим нижним пределом, то есть имеет порядок постоянной aDα/(αD)a^{D-\alpha}/(\alpha-D). Если же α<D\alpha<D, он, напротив, определяется своим верхним пределом и имеет порядок LDα=N1α/DL^{D-\alpha}=N^{1-\alpha/D}. Наконец, в случае α=D\alpha=D первообразная является логарифмом и интеграл равен ln(L/a)=1DlnN\ln(L/a)=\tfrac1D\ln N. Подставляя, получаем тем самым

Uintcste×{N,α>D,NlnN,α=D,N2α/D,α<D.\begin{aligned} U_{\mathrm{int}} \sim \text{cste} \times \begin{cases} N, & \alpha>D,\\[1mm] N\ln N, & \alpha=D,\\[1mm] N^{\,2-\alpha/D}, & \alpha<D. \end{cases} \end{aligned}

Физически истолкование ясно. При α>D\alpha>D мы имеем дело с короткодействующими силами: убывание EpE_p пересиливает рост числа соседей. Каждая составляющая чувствует лишь своё ближайшее окружение, так что её энергия не зависит от размера системы, а полная энергия пропорциональна NN, то есть экстенсивна.

При αD\alpha\leq D всё наоборот: пересиливают далёкие соседи, которых гораздо больше, каждая составляющая чувствует систему целиком, и её собственная энергия растёт с размером системы. Полная энергия растёт тогда быстрее, чем NN, и экстенсивность теряется. Достаточно короткодействующие взаимодействия, таким образом, естественно приводят к экстенсивной энергии, тогда как дальнодействующие взаимодействия могут это свойство разрушить.

Например, в трёхмерной ньютоновской гравитации

Ep(r)1r,D=3,α=1,E_p(r)\sim-\frac1r, \qquad D=3, \qquad \alpha=1,

и предыдущая оценка даёт при фиксированной плотности

UgravN5/3.U_{\mathrm{grav}}\sim-N^{5/3}.

Отрицательный знак, которого сама по себе оценка порядка величины не даёт, отражает притягивающий характер гравитации.

Экстенсивность энергии будет широко использоваться в дальнейшем изложении. Гравитация, конечно, универсальна, но это очень слабая сила: если изучать два приведённых в контакт объёма газа, ею на практике можно пренебречь. С другой стороны, внутренние силы в газе, ван-дер-ваальсова типа, убывают очень быстро (как 1/r61/r^6), так что экстенсивность в этом случае почти точна.

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

7. Итоги

Подведём итог тому, что установила эта лекция и что будет использоваться в дальнейшем.

  • Первое начало постулирует существование функции состояния U:ERU:\mathcal E\longrightarrow\mathbb R, внутренней энергии, и утверждает сохранение энергии изолированной системы.
  • В дальнейшем изложении мы будем считать её также экстенсивной и достаточно непрерывной и дифференцируемой по всем своим переменным.
  • Дополненное аддитивностью энергии, первое начало позволяет сформулировать сохранение как энергетический баланс между замкнутой системой и её окружением. Поскольку работа определена независимо механикой, теплота определяется тогда как оставшийся перенос, что и приводит к полезной формуле ΔU=Q+W\Delta U=Q+W.
  • Для квазистатического процесса этот энергетический баланс принимает дифференциальную форму dU=δQ+δWdU=\delta Q+\delta W, в которой лишь dUdU является полным дифференциалом.

Анализ открытых систем будет рассмотрен в другом месте курса.

8. Литература

О концептуальной эволюции первого начала, от Джоуля к Каратеодори и Борну, см. Розенберг [5]. Классическое изложение см. у Пиппарда [3] и Каллена [4].

  1. J. P. Joule, “On the Mechanical Equivalent of Heat,” Phil. Trans. R. Soc. Lond. 140, 61—82 (1850)
  2. C. Carathéodory, “Untersuchungen über die Grundlagen der Thermodynamik,” Math. Ann. 67, 355—386 (1909)
  3. A. B. Pippard, Elements of Classical Thermodynamics for Advanced Students of Physics, Cambridge University Press (1957)
  4. H. B. Callen, Thermodynamics and an Introduction to Thermostatistics, 2nd ed., Wiley (1985)
  5. R. M. Rosenberg, “From Joule to Caratheodory and Born: A Conceptual Evolution of the First Law of Thermodynamics,” J. Chem. Educ. 87, 691—693 (2010)
  6. A. M. Steane, “First Law, internal energy,” chap. 7 in Thermodynamics: A Complete Undergraduate Course, Oxford University Press (2017)
  7. E. A. Gislason and N. C. Craig, “Cementing the foundations of thermodynamics: Comparison of system-based and surroundings-based definitions of work and heat,” J. Chem. Thermodynamics 37, 954—966 (2005)