Традиционный вывод уравнения Ван-дер-Ваальса
Состояние определенной массы любого вещества можно описать с помощью трех параметров:
давления \(p\),
объема \(V\) и
температуры \(T\). Эти параметры связаны между собой. Их взаимосвязь описывается
уравнением состояния, которое в общем случае имеет вид: \[F\left( {p,V,T} \right) = 0.\] Конкретный вид уравнения зависит от свойств вещества. Например, разреженный газ при достаточно высокой температуре хорошо описывается моделью
идеального газа. Уравнением состояния для него является известное
уравнение Клапейрона (\(1799-1864\)), предложенное в \(1834\) году: \[pV = \frac{m}{M}RT.\] Здесь \(m\) − масса газа, \(M\) − молярная масса (т.е. масса одного моля данного газа), \(R\) − универсальная газовая постоянная. Для одного моля газа это уравнение принимает следующий вид: \[pV = RT.\] Проведенные позднее эксперименты выявили отклонение в поведении реальных газов от законов идеального газа. Эти результаты были обобщены голландским физиком
Яном Дидериком Ван-дер-Ваальсом (\(1837-1923\)), который в \(1873\) году предложил более точное уравнение состояния реального газа. Оно называется
уравнением Ван-дер-Ваальса и в расчете на один моль записывается в виде \[\left( {p + \frac{a}{{{V^2}}}} \right)\left( {V - b} \right) = RT.\] Данное уравнение учитывает силы притяжения и отталкивания, действующие между молекулами. Силы притяжения учитываются благодаря
пристеночному эффекту. Действительно, для частиц, находящихся во внутренней области, силы притяжения со стороны других молекул в среднем скомпенсированы. Однако для частиц вблизи стенок сосуда возникает нескомпенсированная сила притяжения \(f,\) направленная внутрь сосуда. Эта сила, с одной стороны, пропорциональна концентрации частиц \(n\) в сосуде, а с другой стороны − пропорциональна концентрации частиц в пристеночном слое. В результате получаем: \[f \sim {n^2} \sim \frac{1}{{{V^2}}},\] где \(n\) − концентрация молекул в сосуде, \(V\) − объем \(1\) моля газа.
Рассмотренный эффект притяжения молекул пристеночного слоя приводит к уменьшению давления на стенки сосуда. При формальном переходе от уравнения Клапейрона к уравнению Ван-дер-Ваальса это соответствует замене \[p \to p + \frac{a}{{{V^2}}},\] где \(a\) − коэффициент, зависящий от конкретного газа и размеров сосуда.
Силы отталкивания между молекулами в модели Ван-дер-Ваальса учитываются очень просто: предполагается, что молекулы имеют форму шара радиуса \(r\) и не могут приблизиться друг к другу на расстояние между центрами, меньшее чем \(2r.\) Можно считать, что вокруг одной из двух молекул существует "запрещенный" (исключенный) объем (рисунок \(1\)), равный \[\frac{4}{3}\pi {\left( {2r} \right)^3} = 8 \cdot \frac{4}{3}\pi {r^3}.\] Следовательно, в расчете на одну молекулу исключенный объем равен \[{b_0} = 4 \cdot \frac{4}{3}\pi {r^3} = 4{V_0},\] где \({V_0}\) − объем одной молекулы.
В результате , если в уравнении Клапейрона объем пространства, доступного для движения молекул, был равен \(V,\) то теперь он становится равным \[V - {N_A}{b_0} = V - b,\] где \({N_A}\) − число Авогадро (равное числу молекул в одном моле газа), \(b\) − исключенный объем, обусловленный отталкиванием молекул.
|
|
|
Рис.1
|
|
Рис.2 Ян Дидерик Ван-дер-Ваальс
|
Предполагая, что силы притяжения и отталкивания молекул действуют совместно и заменяя \(p\) и \(V\) новыми выражениями, получаем из уравнения Клапейрона уравнение Ван-дер-Ваальса: \[\left( {p + \frac{a}{{{V^2}}}} \right)\left( {V - b} \right) = RT.\]
Изотермы уравнения Ван-дер-Ваальса
При фиксированной температуре уравнение Ван-дер-Ваальса описывает зависимость \(p\left( V \right).\) На плоскости \(pV\) такая зависимость изображается в виде семейства
изотерм, каждая из которых соответствует определенной температуре. Для исследования характера этой зависимости преобразуем уравнение Ван-дер-Ваальса к следующему виду: \[ {\left. {\left( {p + \frac{a}{{{V^2}}}} \right)\left( {V - b} \right) = RT\;} \right| \cdot {V^2},}\;\; {\Rightarrow \left( {p{V^2} + a} \right)\left( {V - b} \right) = RT{V^2},}\;\; {\Rightarrow p{V^3} + aV - pb{V^2} - ab - RT{V^2} = 0,}\;\; {\Rightarrow \left. {p{V^3} - \left( {pb + RT} \right){V^2} + aV - ab = 0\;} \right|:p,}\;\; {\Rightarrow {V^3} - \left( {b + \frac{{RT}}{p}} \right){V^2} + \frac{a}{p}V - \frac{{ab}}{p} = 0.} \] При фиксированном значении \(p\) полученное уравнение является уравнением третьей степени по
переменной \(V\). Известно, что кубическое уравнение может иметь \(1\) или \(3\) действительных корня. Первый случай встречается при высоких температурах \(T\) (зеленая изотерма \(AB\) на рисунке \(3\)). При понижении температуры на изотерме появляется волнообразный участок. В этом случае существует три корня (синяя изотерма \(EFGHIJL\)). Переход между двумя типами изотерм происходит при определенной температуре \({T_K},\) которая называется
критической температурой.
Критическая точка
В критической точке газ характеризуется критическими параметрами \({T_K}\), \({p_K}\) и \({V_K},\) которые определяются лишь свойствами газа. С алгебраической точки зрения изотерма в критической точке имеет три действительных корня, которые равны между собой. Это обстоятельство позволяет легко вычислить значения \({T_K}\), \({p_K}\) и \({V_K}.\) Действительно, в этом случае уравнение состояния газа \[{V^3} - \left( {b + \frac{{RT}}{p}} \right){V^2} + \frac{a}{p}V - \frac{{ab}}{p} = 0\] должно записываться как \[{\left( {V - {V_K}} \right)^3} = 0.\] Раскрывая куб разности и приравнивая коэффициенты при соответствующих степенях \(V,\) находим выражения для критических параметров: \[ {{\left( {V - {V_K}} \right)^3} = 0,}\;\; {\Rightarrow {V^3} - 3{V^2}{V_K} + 3V{V_K} - V_K^3 = 0.} \] Учтем, что в критической точке \(p = {p_K}\), \(T = {T_K}.\) Следовательно, получаем: \[ {\left\{ \begin{array}{l} - \left( {b + \frac{{RT}}{p}} \right) = - 3{V_K}\\ \frac{a}{{{p_K}}} = 3V_K^2\\ - \frac{{ab}}{{{p_K}}} = - V_K^3 \end{array} \right.,}\;\; {\Rightarrow \left\{ \begin{array}{l} b + \frac{{RT}}{p} = 3{V_K}\\ \frac{a}{{{p_K}}} = 3V_K^2\\ \frac{{ab}}{{{p_K}}} = V_K^3 \end{array} \right..} \] Разделим третье уравнение на второе: \[\require{cancel} {\frac{{\cancel{a}b \cdot \cancel{p_K}}}{{\cancel{p_K} \cdot \cancel{a}}} = \frac{{V_K^{\cancel{3}}}}{{3\cancel{V_K^2}}},}\;\; {\Rightarrow {V_K} = 3b.} \] Из второго уравнения определяем \({p_K}:\) \[ {{p_K} = \frac{a}{{3V_K^2}} } = {\frac{a}{{3 \cdot {{\left( {3b} \right)}^2}}} } = {\frac{a}{{27{b^2}}}.} \] Наконец, из первого уравнения находим критическую температуру \({T_K}:\) \[ {{T_K} = \frac{{\left( {3{V_K} - b} \right){p_K}}}{R} } = {\frac{{\left( {3 \cdot 3b - b} \right) \cdot \frac{a}{{27{b^2}}}}}{R} } = {\frac{{8a}}{{27bR}}.} \] Итак, критические параметры газа в модели Ван-дер-Ваальса зависят лишь от величин \(a, b\) и определяются формулами \[ {{V_K} = 3b,}\;\;\; {{p_K} = \frac{a}{{27{b^2}}},}\;\;\; {{T_K} = \frac{{8a}}{{27bR}}.} \] Заметим, что мы нашли критические величины без использования производной. Можно показать, что в критической точке с параметрами \({V_K}\), \({p_K}\), \({T_K}\) первая и вторая производная функции \(p\left( V \right)\) равны нулю (смотрите Пример \(1\)).
Приведенное уравнение состояния
Используя критические параметры \({V_K}\), \({p_K}\) и \({T_K}\) можно перейти к безразмерным переменным: \[ {\varphi = \frac{V}{{{V_K}}},}\;\;\; {\pi = \frac{p}{{{p_K}}},}\;\;\; {\tau = \frac{T}{{{T_K}}}} \] и переписать уравнение Ван-дер-Ваальса
в приведенной форме: \[ {V = \varphi {V_K},\;p = \pi {p_K},\;T = \tau {T_K},}\;\; {\Rightarrow \left( {\pi {p_K} + \frac{a}{{{{\left( {\varphi {V_K}} \right)}^2}}}} \right)\left( {\varphi {V_K} - b} \right) = R\tau {T_K},}\;\; {\Rightarrow \left( {\frac{{\pi a}}{{27{b^2}}} + \frac{a}{{{{\left( {3b} \right)}^2}{\varphi ^2}}}} \right)\left( {3b\varphi - b} \right) = R\tau \cdot \frac{{8a}}{{27bR}},}\;\; {\Rightarrow \left( {\pi + \frac{3}{{{\varphi ^2}}}} \right)\left( {3\varphi - 1} \right) = 8\tau .} \] Данное уравнение является более универсальным, чем исходное уравнение Ван-дер-Ваальса. Изотермы для различных веществ, построенные в соответствии с приведенным уравнением, будут совпадать при одинаковых значениях \(\tau.\) Состояния разных тел, которые описываются одинаковыми приведенными параметрами \(\pi,\) \(\varphi,\) \(\tau,\) называются
соответственными состояниями. Если два каких-либо тела или вещества обладают любыми двумя одинаковыми приведенными параметрами, то они имеют и третью одинаковую величину. Это универсальное свойство называется
законом соответственных состояний.
Рассмотрим еще один интересный инвариант. Поскольку в модели Ван-дер-Ваальса критические величины \({V_K}\), \({p_K}\), \({T_K}\) выражаются лишь через два параметра \(a\) и \(b,\) мы можем составить комбинацию этих величин, не зависящую от \(a\) и \(b.\) Можно записать следующее соотношение: \[ {\frac{{{p_K}{V_K}}}{{{T_K}}} } = {\frac{{3b \cdot \frac{a}{{27{b^2}}}}}{{\frac{{8a}}{{27bR}}}} } = {\frac{{3\cancel{a}\cancel{b} \cdot \cancel{27}\cancel{b}R}}{{\cancel{27}\cancel{b^2} \cdot 8\cancel{a}}} } = {\frac{{3R}}{8}}\;\;\; {\text{или}\;\;\;\frac{{{p_K}{V_K}}}{{{T_K}R}} } = {\frac{3}{8} = 0,375.} \] Это отношение называется
критическим коэффициентом и не зависит ни от каких величин, характеризующих то или иное вещество, т.е. является инвариантной величиной.
Фазовый переход "газ-жидкость" и метастабильные состояния
Рассмотрим снова уравнение Ван-дер-Ваальса при температурах не выше критической: \(T < {T_K}.\) В этом случае изотермы имеют волнообразный характер. Участок такой изотермы между минимумом и максимумом характеризуется положительной производной \(\large\frac{{dp}}{{dV}}\normalsize > 0,\) что соответствует
нестабильному состоянию вещества. При малой флуктуации объема \(dV > 0\) давление газа будет также возрастать (так как \(dp > 0\)), что приведет к взрывному расширению газа. При отрицательной флуктуации \(dV < 0,\) наоборот, будет происходить быстрый коллапс вещества.
На самом деле в этой области на \(pV\)-диаграмме происходит фазовый переход "газ-жидкость", сопровождающийся значительным изменением объема \(V\) при постоянном давлении \(p\) и постоянной температуре \(T.\) Такой переход на \(pV\)-диаграмме изображается прямолинейным горизонтальным участком (отрезок \(FJ\) на рисунке \(4\)). В сосуде при этом возникает двухфазное состояние
газ + жидкость. При перемещении фазовой точки влево вдоль отрезка \(JF\) объем газа уменьшается, а объем жидкости, соответственно, возрастает.
Расположение горизонтального участка \(FJ\) определяется из термодинамических соображений − на основании т.н.
правила Максвелла, согласно которому площади криволинейных фигур \(FGH\) и \(HIJ\) должны быть равными (смотрите выше рисунок \(4\)).
Множество начальных и конечных точек горизонтальных участков при разных температурах \(T\) (где \(T \le {T_K}\)) образует кривую, которая называется
бинодалью. Кривая, проходящая через точки минимума и точки максимума всех изотерм при \(T \le {T_K},\) называется
спинодалью.
Нестабильные состояния (для которых \(\large\frac{{dp}}{{dV}}\normalsize > 0\)) находятся внутри области, ограниченной спинодалью. Область диаграммы между спинодалью и бинодалью в принципе удовлетворяет критерию устойчивости \(\large\frac{{dp}}{{dV}}\normalsize < 0.\) Поэтому состояния вещества, которые реализуются в этой области, называются
метастабильными состояниями. В левой половине указанной области образуется
перегретая жидкость, а в правой половине −
переохлажденный пар.
Удивительно, что такая простая модель реального газа, в которой лишь грубо учитываются межмолекулярные силы, способна отражать реальный набор различных состояний газа и жидкости. Кроме регулярного газа и жидкости, фазового перехода "газ-жидкость" и рассмотренных метастабильных состояний, модель Ван-дер-Ваальса фиксирует также такое экзотическое явление, как
растянутая жидкость (участок \(MNS\) на рисунке \(4\)). Это состояние жидкости является неустойчивым.
В приведенных ниже задачах мы акцентируем внимание на исследовании уравнения Ван-дер-Ваальса с помощью производной.