В рамках курса Массачусетского технологического института (MIT OpenCourseWare) профессор Джан Паоло Беретта представляет глубокий анализ математических методов термодинамики, связывающих фундаментальные уравнения стабильных равновесных состояний с реальными физическими процессами. Лекция охватывает переход от сложных математических преобразований, таких как определители Якоби, к физическим концепциям доступной энергии и знаменитому принципу Ле Шателье — Брауна. Главной практической иллюстрацией теоретических выкладок становится работа современных воздушных пружин, широко применяемых в транспортной индустрии.
🧮 Математический аппарат термодинамики: Якобианы и их свойства 0:12
Изучение свойств стабильных равновесных состояний неизбежно требует развитого математического аппарата для извлечения информации из фундаментальных термодинамических соотношений системы. По словам профессора Беретты, чистые математические манипуляции обретают ценность только тогда, когда им даётся надлежащая физическая интерпретация. Профессор приводит в пример экономику, где аналогичные методы математического анализа принесли Полу Самуэльсону, работавшему в MIT, Нобелевскую премию, поскольку тот осознал универсальность инструментов манипулирования функциями.
Вслед за преобразованиями Лежандра, позволившими получить свободную энергию Гельмгольца, свободную энергию Гиббса, энтальпию и свободную энергию Эйлера, а также соотношения Максвелла, Беретта предлагает использовать детерминанты Якоби (якобианы).
Если имеется система из двух функций двух переменных, матрица их частных производных называется матрицей Якоби, а её определитель — якобианом, который обозначается специальным символом, напоминающим частную производную вектора функций по вектору независимых переменных:
$$\frac{\partial(f, g)}{\partial(x, y)} = \det \begin{pmatrix} \frac{\partial f}{\partial x} & \frac{\partial f}{\partial y} \ \frac{\partial g}{\partial x} & \frac{\partial g}{\partial y} \end{pmatrix}$$
Профессор Беретта выделяет ключевые свойства якобианов, напрямую следующие из линейной алгебры:
- Перестановка строк или столбцов меняет знак определителя на противоположный.
- Если одна из функций тривиальна (например, $g = y$), то вторая строка матрицы содержит производные независимой переменной $y$ по $x$ (что равно $0$) и по $y$ (что равно $1$), превращая якобиан в стандартную частную производную $\left(\frac{\partial f}{\partial x}\right)_y$.
Данное свойство позволяет легко переходить от обычных частных производных к эквивалентным якобианам для их дальнейшего упрощения. Кроме того, для систем функций справедлива теорема об обратной функции: якобиан обратной системы равен обратной величине исходного якобиана. Из этого строго доказывается термодинамическое правило переворота дроби для частных производных при фиксации одной переменной.
Наконец, Беретта напоминает цепное правило для параметрического представления: если независимые переменные $x$ и $y$ сами являются функциями новых параметров $a$ и $b$, то новый якобиан равен произведению исходного якобиана на якобиан перехода:
$$\frac{\partial(f, g)}{\partial(a, b)} = \frac{\partial(f, g)}{\partial(x, y)} \cdot \frac{\partial(x, y)}{\partial(a, b)}$$
Информально это свойство можно представить как символьное «сокращение» промежуточных переменных в числителе и знаменателе.
📉 Связь термодинамических свойств и матрица Гессе 10:24
В качестве примера использования якобианов Беретта рассматривает определитель матрицы Гессе (гессиан) для фундаментального уравнения энергии по энтропии $S$ и объёму $V$. Матрица составляется из вторых частных производных энергии. Поскольку первая производная энергии по энтропии — это температура $T$, а по объёму — давление с минусом ($-p$), гессиан можно переписать через сопряжённые переменные. Вынося знак минус за пределы определителя, мы получаем чистый якобиан температуры и давления по энтропии и объёму:
$$\det H = -\frac{\partial(T, p)}{\partial(S, V)}$$
Используя искусственный приём со вставкой промежуточных переменных в пустое пространство числителя и знаменателя, профессор раскладывает это выражение двумя путяюми. Первый путь — введение пары $(S, p)$. Из-за повторяющихся переменных в числителях и знаменателях выражение распадается на произведение двух частных производных: $\left(\frac{\partial T}{\partial S}\right)_p$ и $\left(\frac{\partial p}{\partial V}\right)_S$.
Опираясь на ранее введённое определение теплоёмкости при постоянном давлении $C_p$ через энтальпию, где при постоянном давлении $dH = TdS$, Беретта показывает, что:
$$\left(\frac{\partial S}{\partial T}\right)_p = \frac{C_p}{T}$$
Следовательно, первая производная представляет собой обратную величину $\frac{T}{C_p}$. Вторая производная связана с изоэнтропийной (адиабатической) сжимаемостью материала $\kappa_S$, определяющей сжатие при постоянной энтропии. Так как определитель матрицы Гессе должен быть неотрицательным по условиям термодинамической устойчивости, это накладывает строгие ограничения на знаки термодинамических коэффициентов.
Второй путь разложения — введение пары переменных $(T, V)$. Аналогичные подстановки выражают гессиан через теплоёмкость при постоянном объёме $C_V$ и изотермическую сжимаемость $\kappa_T$:
$$\det H = \frac{T}{C_V \cdot V \cdot \kappa_T}$$
Поскольку оба подстановки преобразуют одно и то же исходное выражение, профессор приравнивает полученные группы переменных. После сокращения $T$ и $V$ выводится фундаментальное термодинамическое тождество:
$$\frac{C_p}{C_V} = \frac{\kappa_T}{\kappa_S}$$
Для воздуха или любого двухатомного газа отношение теплоёмкостей $C_p/C_V$ общеизвестно и составляет 1.4 (или 7/5). По словам лектора, это автоматически означает, что отношение изотермической сжимаемости к изоэнтропийной для воздуха также строго равно 1.4. Данный факт имеет колоссальное прикладное значение для создания воздушных рессор и систем амортизации.
🔄 Циклические соотношения и уравнение Майера 18:06
Если система описывается уравнением состояния, связывающим давление, объём и температуру ($pVT$-взаимосвязь), её можно представить в виде трёхмерной поверхности. В термодинамике такую поверхность для однофазных состояний можно эквивалентно рассматривать с разных точек зрения, выбирая любые две переменные в качестве независимых. Из свойств якобианов для циклической перестановки трёх переменных выводится так называемое циклическое соотношение:
$$\left(\frac{\partial x}{\partial y}\right)_z \left(\frac{\partial y}{\partial z}\right)_x \left(\frac{\partial z}{\partial x}\right)_y = -1$$
Применяя это правило к параметрам $p, V, T$, можно связать температурные изменения энтропии с легко измеряемыми свойствами. Одно из соотношений Максвелла позволяет оценивать изменение энтропии через изменение объёма и давления, что профессор Беретта метафорически называет основой для гипотетического прибора — «измерителя энтропии». Циклическое соотношение доказывает, что искомая производная равна отношению изобарного коэффициента расширения $\alpha_p$ к изотермической сжимаемости $\kappa_T$.
Для идеального газа, где $V = \frac{nRT}{p}$, эти коэффициенты принимают простой вид:
- Изобарный коэффициент расширения: $\alpha_p = \frac{1}{V}\left(\frac{\partial V}{\partial T}\right)_p = \frac{1}{T}$.
- Изотермическая сжимаемость: $\kappa_T = -\frac{1}{V}\left(\frac{\partial V}{\partial p}\right)_T = \frac{1}{p}$.
В результате их отношение для идеального газа упрощается до величины $p/T$.
Ещё одним важным результатом детерминантного анализа является уравнение Майера, связывающее разность теплоёмкостей $C_p$ и $C_V$ с коэффициентами $\alpha_p$ и $\kappa_T$. Переводя производные в якобианы и подставляя промежуточные параметры $T$ и $V$, Беретта выстраивает матричную структуру, где компоненты симметричны благодаря свойствам гессиана и соотношениям Максвелла. Итоговое уравнение Майера имеет вид:
$$C_p - C_V = \frac{V T \alpha_p^2}{\kappa_T}$$
Переходя к удельным (молярным) теплоёмкостям $c_p = C_p/n$ и $c_V = C_V/n$, для идеального газа уравнение Майера трансформируется в классическую форму $c_p = c_V + R$, где $R$ — универсальная газовая постоянная.
⚡ Функции доступной энергии и термодинамические резервуары 29:56
Профессор Беретта возвращается к концепции доступной энергии (available energy). В системе, взаимодействующей с резервуаром, протекает процесс с поднятием или опусканием груза (weight process). Если процесс обратим, из системы можно извлечь максимальную полезную работу при переходе в состояние взаимного равновесия с окружающей средой. В общем случае рассматриваются два произвольных состояния (начальное и конечное), и рассчитывается максимальная получаемая или минимальная затрачиваемая механическая работа.
Задача всегда решается совместным составлением баланса энергии, баланса энтропии и уравнений свойств участвующих сред. Для классического теплового резервуара фундаментальное уравнение выражает линейную связь между его энергией и энтропией через постоянную температуру $T_R$. В физике такой резервуар часто называют термической баней (heat bath).
Как объясняет профессор Беретта, это название оправдано тем, что резервуар готов обеспечивать тепловое взаимодействие: при извлечении энергии он отдаёт энтропию в строго пропорциональном количестве. Однако лектор подчёркивает, что резервуар не является «хранилищем тепла», поскольку тепло — это не форма энергии, а форма её передачи.
«Свободное обращение с терминами допустимо, если вы понимаете, что делаете. Но если придать бытовой смысл строгому физическому понятию, возникнет путаница. Наша цель — быть максимально аккуратными с языком науки», — напоминает Беретта.
Максимальная обратимая работа равна разности канонических функций доступности $\gamma$ для двух состояний системы:
$$\gamma = E - T_R S$$
Графически функция $\gamma$ имеет абсолютный минимум в точке взаимного равновесия системы с резервуаром. Из факта минимальности этой функции выводится условие выпуклости фундаментального уравнения: при любом виртуальном смещении из равновесия функция $\gamma$ должна возрастать. Это доказывает, что вторая производная энергии по энтропии строго положительна, а вторая производная энтропии — отрицательна, что гарантирует неотрицательность теплоёмкостей. Данный вывод составляет первую часть принципа Ле Шателье — Брауна.
Физический смысл заключается в следующем: если систему, находящуюся в равновесии с резервуаром, возмутить (например, сообщить ей энергию), она отреагирует увеличением температуры. Температура представляет собой «стремление избежать удержания энергии» (escaping tendency). Возросшая температура инициирует спонтанный теплообмен с резервуаром, направленный на компенсацию внешнего воздействия и возвращение к равновесию.
📦 Расширение концепции: Резервуары с переменным объёмом и осмотические системы 41:17
В реальной инженерной практике термодинамические резервуары не ограничиваются только тепловым обменом. Окружающая атмосфера, моря или реки могут изменять свой объём при взаимодействии с технической системой. Если эти изменения относительно малы по сравнению с масштабами окружающей среды, её температура $T_R$ и давление $p_R$ остаются неизменными. Для резервуара с переменным объёмом фундаментальное уравнение описывает плоскость, наклоны которой равны $T_R$ и $-p_R$.
При расчёте балансов для комбинированной системы добавляется уравнение сохранения общего объёма. В результате решения системы уравнений выводится новая функция доступности $\Psi$:
$$\Psi = E - T_R S + p_R V$$
Беретта обращает внимание, что функционально $\Psi$ напоминает свободную энергию Гиббса, но тождественно совпадает с ней только в точке полного взаимного равновесия с окружающей средой, когда температура и давление системы в точности равны $T_R$ и $p_R$. В любом другом состоянии, особенно неравновесном, у системы может вообще не быть единой определённой температуры, однако макроскопическая функция доступности $\Psi$ остаётся строго определённой и задаёт лимиты работы.
Математический анализ устойчивости в трёхмерном пространстве (с учётом изменения объёма) требует, чтобы при любом отклонении от равновесия изменение $\Psi$ было положительным. Это приводит к требованию, чтобы матрица вторых производных энтропии (связанная с квадратичной формой для приращений $dE$ и $dV$) была отрицательно полуопределённой. Из линейной алгебры известно, что для матрицы размера $2 \times 2$ отрицательная полуопределённость означает отрицательность элементов на главной диагонали, но её определитель (произведение двух отрицательных собственных значений) должен быть строго больше нуля.
Профессор иллюстрирует это на классической задаче о создании вакуума в цилиндре под поршнем, находящемся в атмосфере с давлением $p_R$. Интуитивно понятно, что минимальная обратимая работа, которую необходимо затратить на вытягивание поршня для получения пустоты объёмом $\Delta V$, в точности равна $W_{min} = p_R \Delta V$. Термодинамическая формула полностью подтверждает этот механический вывод.
Аналогичный подход Беретта применяет к осмотическим резервуарам, способным обмениваться с системой только определённым типом вещества через полупроницаемые мембраны. Актуальным примером является проектирование систем улавливания углекислого газа ($CO_2$) из атмосферы для расчёта их термодинамической эффективности по второму закону.
Один из гипотетических путей — создание вакуума под поршнем с последующим открытием селективной мембраны для $CO_2$. Поскольку изначально под поршнем концентрация углекислого газа равна нулю, процесс его диффузии через мембрану до выравнивания химических потенциалов будет самопроизвольным и необратимым, что приведёт к генерации энтропии. Для этой системы вводится осмотическая функция доступности, включающая химический потенциал резервуара $\mu_R$:
$$\Omega = E - T_R S - \mu_R n$$
🌐 Обобщённый резервуар Эйлера и устойчивость равновесия 1:11:44
Финальным обобщением теории становится «большой резервуар» (grand reservoir), с которым система может свободно обмениваться энергией, энтропией, объёмом и массами всех $r$ химических компонентов. Изменение фундаментальных свойств такого резервуара описывается расширенным линейным соотношением. Предельная полезная работа в этом случае определяется через разность значений функции доступности Эйлера:
$$\Xi = E - T_R S + p_R V - \sum_{i=1}^r \mu_{Ri} n_i$$
Матрица Гессе для проверки устойчивости расширяется до размеров $(2+r) \times (2+r)$. Находясь в точке взаимного равновесия, мы можем виртуально изменять состояние системы в различных направлениях.
Профессор Беретта показывает, что если зафиксировать температуру и давление на уровне резервуара, но изменять количества компонентов, то требование минимизации функции Эйлера переходит в требование минимизации свободной энергии Гиббса системы (её второй дифференциал должен быть неотрицательным). Если же зафиксировать температуру и объём, условие стабильности трансформируется в требование минимума свободной энергии Гельмгольца.
⚖️ Принцип Ле Шателье — Брауна и его физический смысл 1:19:43
Чтобы строго доказать вторую, более глубокую часть принципа Ле Шателье — Брауна, Беретта обращается к каноническим квадратичным формам. Для функции двух переменных квадратичная форма содержит квадраты приращений независимых переменных и их перекрёстное произведение, что затрудняет мгновенное определение знака выражения. Алгебраический приём позволяет выделить полный квадрат и привести форму к каноническому виду, где перед вторым дифференциалом возникает коэффициент $\lambda$.
Для положительно определённой матрицы коэффициенты $\lambda$ также должны быть неотрицательными, причём коэффициент $\lambda_y$ всегда строго меньше исходной второй производной по $y$ из-за вычитания положительной величины, связанной с перекрёстными производными.
Математическое раскрытие этих коэффициентов показывает, что они представляют собой частные производные сопряжённых потенциалов, но при принципиально разных условиях фиксации параметров. В одном случае фиксируется естественная переменная, в другом — сопряжённый ей потенциал.
Это даёт цепочку термодинамических неравенств:
$$C_p \ge C_V \ge 0$$ $$\kappa_T \ge \kappa_S \ge 0$$
Физический смысл принципа, сформулированного Анри Луи Ле Шателье и Карлом Фердинандом Брауном, заключается в качественной оценке силы ответной реакции физической системы на внешнее воздействие. Если в систему добавляется энтропия, её температура растёт, стимулируя отдачу тепла наружу для ослабления внешнего возмущения. Однако интенсивность этого температурного ответа существенно зависит от наложенных связей: при постоянном объёме реакция системы (рост температуры) будет значительно сильнее, чем при постоянном давлении.
При фиксации давления система сохраняет частичное механическое равновесие с окружающей средой, поэтому её отклик на температурное возмущение оказывается более «мягким». Напротив, при фиксации объёма нарушаются сразу два условия взаимного равновесия (и по температуре, и по давлению), что заставляет систему реагировать гораздо жёстче, стремясь вернуться в исходную точку.
🚗 Практическое применение: Воздушные пружины в машиностроении 1:37:38
В качестве яркого примера действия выведенных законов Беретта подробно разбирает устройство воздушной пружины (рессоры). Конструктивно она представляет собой герметичный цилиндр с поршнем, внутри которого находится воздух. Инженерный расчёт жёсткости такой пружины (коэффициента $k$, связывающего приращение возвращающей силы с величиной деформации) показывает, что жёсткость обратно пропорциональна сжимаемости газа, находящегося в знаменателе формулы.
Здесь термодинамика разделяет работу устройства на два полярных режима:
- Медленное сжатие (изотермический режим): если на пружину давить очень медленно, воздух успевает обмениваться теплом с окружающей средой, оставаясь в термическом равновесии. Жёсткость пружины в этом случае определяется изотермической сжимаемостью $\kappa_T$.
- Быстрый удар (изоэнтропийный режим): если автомобиль или поезд наезжает на дорожное препятствие (кочку), сжатие поршня происходит мгновенно. У системы физически нет времени на теплообмен с резервуаром, и процесс идёт при постоянной энтропии. Жёсткость пружины начинает определяться изоэнтропийной сжимаемостью $\kappa_S$.
Поскольку отношение этих двух видов сжимаемости для воздуха строго равно отношению теплоёмкостей $C_p/C_V = 1.4$, отношение динамической жёсткости к статической также равно 1.4.
«Это означает, что при быстром динамическом ударе воздушная пружина автоматически становится на 40% жёстче, чем при медленном плавном сжатии», — констатирует профессор Беретта.
Производство воздушных рессор сегодня представляет собой огромный мировой бизнес. Благодаря уникальному термодинамическому свойству автоматически менять жёсткость в сочетании с гидравлическими демпферами, эти устройства являются незаменимым элементом подвесок современных высокоскоростных поездов, грузовых автомобилей, автобусов и даже некоторых премиальных моделей велосипедов, обеспечивая эффективное гашение вибраций и высокий комфорт передвижения.
