Задача 5.4.8

Решение

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

Режим молекулярных пуль (разреженный газ, $\lambda \gg d$)

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

Поскольку длина свободного пробега $\lambda$ колоссально превосходит расстояние между дисками $d$, эта молекула летит сквозь зазор как микроскопическая пуля, не встречая препятствий на своем пути. Она бьет точно в верхний подвешенный диск, отдавая ему приобретенный тангенциальный импульс целиком и полностью. Крутящий момент, передаваемый в таком баллистическом режиме, определяется исключительно плотностью «огня» — количеством молекул, пересекающих зазор в единицу времени. Чем выше давление газа $P$, тем пропорционально больше концентрация частиц в объеме, тем интенсивнее бомбардировка мишени. Удвоение давления означает строгое удвоение числа независимых переносчиков импульса, поэтому передаваемый момент сил и, как следствие, угол закручивания $\varphi$ растут строго линейно с ростом давления.

Режим вязкого континуума (высокое давление, $\lambda \ll d$)

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

Здесь на сцену выходит один из самых поразительных и контринтуитивных результатов кинетической теории — великий парадокс Максвелла. В 1866 году Джеймс Клерк Максвелл в своей знаменитой Бейкеровской лекции математически предсказал и экспериментально доказал, что динамическая вязкость идеального газа $\eta$ абсолютно не зависит от его давления. Физическая природа этого макроскопического парадокса элегантна в своей микроскопической симметрии: при повышении давления мы увеличиваем плотность газа, то есть количество переносчиков импульса в единице объема линейно возрастает. Однако строго пропорционально этому росту сокращается длина свободного пробега $\lambda$ — дистанция, на которую каждый переносчик способен этот импульс транспортировать без рассеяния. Рост плотности идеально и без остатка компенсируется падением дальнодействия молекул.

Газ становится гуще, но передача импульса в нем становится ровно настолько же более локализованной. В результате динамическая вязкость среды фиксируется на жесткой константе, зависящей исключительно от температуры и квантово-механического сечения рассеяния самих молекул. Передаваемый крутящий момент перестает чувствовать нагнетаемое давление, и угол закручивания нити выходит на горизонтальную асимптоту — неизменную константу $\varphi_2$.

Зависимость угла от давления при $\varphi \ll \varphi_2$
Условие $\varphi \ll \varphi_2$ геометрически означает, что система находится на самом начальном участке своего эволюционного пути. Малые углы отклонения реализуются исключительно в режиме сверхнизких давлений, когда газ работает в баллистическом режиме независимых молекулярных пуль. Поскольку в этой асимптотической области момент сил растет строго пропорционально давлению, зависимость угла имеет вид прямой пропорциональности. Зная опорную точку ($P_1$, $\varphi_1$), мы получаем строгую линейную зависимость.

$$\varphi = \varphi_1 \frac{P}{P_1}$$

Для тех, кому интересны математические детали модели

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

Ультраразреженный предел и аккомодация импульса

Обратимся к кинетической теории. Пусть давление $P$ исчезающе мало, так что межмолекулярными столкновениями в зазоре можно пренебречь. Выделим на нижнем вращающемся диске элементарное кольцо радиуса $r$ и ширины $dr$, площадь которого составляет $dS = 2\pi r dr$.

Кинетический поток частиц, то есть число молекул, ударяющихся о единицу площади поверхности в единицу времени, определяется фундаментальной формулой:
$$\nu = \frac{1}{4} n \langle v \rangle$$
где $n$ — объемная концентрация молекул, а $\langle v \rangle = \sqrt{\frac{8kT}{\pi m}}$ — средняя арифметическая скорость теплового движения.

При столкновении молекулы с поверхностью происходит обмен импульсом, эффективность которого характеризуется коэффициентом аккомодации тангенциального импульса Максвелла (TMAC), обозначаемым $\sigma$. Физический смысл этого коэффициента заключается в том, какая доля молекул полностью термализуется поверхностью (диффузное рассеяние), а какая отскакивает абсолютно упруго, сохраняя свой тангенциальный импульс (зеркальное отражение).

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

Для построения идеализированной аналитической модели мы примем гипотезу полного диффузного рассеяния, положив $\sigma = 1$. Это означает, что каждая молекула, покидающая нижний диск, уносит макроскопический тангенциальный импульс $p_\tau = m v_{disk} = m \omega r$. Достигая верхнего неподвижного диска, она передает ему этот импульс целиком.

Элементарная сила вязкого трения, действующая на выделенное кольцевое сечение верхнего диска, равна потоку импульса через это сечение:
$$dF = p_\tau \cdot \nu \cdot dS = (m \omega r) \left( \frac{1}{4} n \langle v \rangle \right) 2\pi r dr$$

Элементарный крутящий момент вычисляется как произведение силы на плечо $r$ ($dM = r dF$):
$$dM = \frac{\pi}{2} m n \langle v \rangle \omega r^3 dr$$

Интегрирование по всей площади диска от центра до внешнего радиуса $R$ дает точное значение момента:
$$M_1 = \int_0^R \frac{\pi}{2} m n \langle v \rangle \omega r^3 dr = \frac{\pi}{8} m n \langle v \rangle \omega R^4$$

Используя уравнение состояния идеального газа $P = n k T$, мы можем выразить концентрацию $n = P/kT$ и получить искомую функциональную зависимость момента от макроскопического давления:
$$M_1(P) = \left( \frac{\pi m \langle v \rangle \omega R^4}{8 k T} \right) P$$
Математика неумолимо подтверждает физическую интуицию: в ультраразреженном пределе крутящий момент строго линеен по давлению.

Газодинамический предел и течение Куэтта
В противоположной крайности высоких давлений газ представляет собой непрерывную сплошную среду. Геометрия задачи сводится к классическому сдвиговому течению, известному в гидродинамике как течение Тейлора-Куэтта в плоском пределе. При малых числах Рейнольдса профиль скорости между дисками строго линеен. Вектор скорости направлен по касательной к концентрическим окружностям, а его модуль меняется от $\omega r$ на вращающемся нижнем диске до нуля на подвешенном верхнем по закону:
$$v(z) = \omega r \frac{z}{d}$$
где координата $z$ отсчитывается от неподвижного верхнего диска вниз.

Закон вязкого трения Ньютона связывает сдвиговое напряжение $\tau$ с градиентом скорости через динамическую вязкость $\eta$:
$$\tau = \eta \frac{\partial v}{\partial z} = \eta \frac{\omega r}{d}$$

Вычислим суммарный крутящий момент, интегрируя это сдвиговое напряжение по всей площади диска:
$$M_2 = \int_0^R r \cdot \tau \cdot 2\pi r dr = \int_0^R r \left( \eta \frac{\omega r}{d} \right) 2\pi r dr = \frac{2\pi \eta \omega}{d} \int_0^R r^3 dr = \frac{\pi \eta \omega R^4}{2 d}$$

Полученное выражение представляет собой предельную константу. Приравнивая этот момент к моменту упругости нити $\gamma \varphi_2$, мы получаем аналитическое выражение для максимального угла закручивания:
$$\varphi_2 = \frac{\pi \eta \omega R^4}{2 d \gamma}$$

Максвелловский парадокс в формулах
Чтобы строго доказать независимость асимптотического момента $M_2$ от давления, необходимо вскрыть внутреннюю структуру динамической вязкости идеального газа $\eta$. Согласно кинетической теории, транспорт импульса в газе описывается формулой:
$$\eta \approx \frac{1}{3} \rho \langle v \rangle \lambda$$

Плотность газа тривиально выражается через концентрацию и массу одной молекулы: $\rho = m n$. Длина свободного пробега $\lambda$, в свою очередь, обратно пропорциональна концентрации частиц и их эффективному сечению рассеяния $\sigma_{scat} = \pi d_{mol}^2$ (где $d_{mol}$ — газокинетический диаметр молекулы):
$$\lambda = \frac{1}{\sqrt{2} \sigma_{scat} n}$$

Подставим эти фундаментальные микроскопические величины в макроскопическое выражение для вязкости:
$$\eta = \frac{1}{3} (m n) \langle v \rangle \left( \frac{1}{\sqrt{2} \sigma_{scat} n} \right) = \frac{m \langle v \rangle}{3 \sqrt{2} \sigma_{scat}}$$

Концентрация $n$, а вместе с ней и термодинамическое давление $P$, математически аннигилируют, сокращаясь в числителе и знаменателе. Динамическая вязкость $\eta$, а следовательно и предельный угол закручивания $\varphi_2$, оказываются зависящими исключительно от массы молекул, их поперечного сечения и температуры системы.

Сшивка через число Кнудсена (Kn)
Истинная красота теоретической физики проявляется не в анализе крайностей, а в способности описать непрерывный переход между асимптотиками. Для построения единой теории нам потребуется безразмерный критерий подобия — число Кнудсена, определяемое как отношение длины свободного пробега к характерному макроскопическому размеру системы: $Kn = \lambda / d$.

В переходном режиме, когда длина свободного пробега $\lambda$ становится сопоставимой с зазором $d$, вступает в игру феномен вязкого проскальзывания на границах (boundary slip). Впервые этот эффект был экспериментально обнаружен Кундтом и Варбургом в 1875 году при наблюдении за затуханием осцилляций диска в разреженном газе, а в 1879 году получил теоретическое обоснование от самого Максвелла.

Суть феномена заключается в том, что газ, находящийся непосредственно у поверхности твердого тела, не покоится относительно нее (классическое условие прилипания нарушается). Формируется так называемый кнудсеновский слой, в котором макроскопическая скорость газа у стенки $v_{slip}$ отлична от скорости самой стенки. Величина этой скорости проскальзывания прямо пропорциональна градиенту скорости в объеме и длине свободного пробега:
$$v_{slip} = \zeta \lambda \left\vert \frac{\partial v}{\partial z} \right\vert$$
где $\zeta$ — безразмерный коэффициент порядка единицы (часто называемый коэффициентом скольжения), который зависит от того самого коэффициента аккомодации импульса $\sigma$. Для абсолютно диффузного рассеяния $\zeta \approx 1$, но для реальных поверхностей он может быть больше.

Наличие проскальзывания на обоих дисках приводит к тому, что профиль скорости становится более пологим. Градиент скорости уменьшается, словно физический зазор между дисками расширился за свои геометрические пределы. Строгое решение дифференциального уравнения со скользящими граничными условиями дает переопределенный градиент:
$$\frac{\partial v}{\partial z} = \frac{\omega r}{d + 2\zeta \lambda}$$

Подставив этот градиент в интеграл крутящего момента, мы получаем обобщенную формулу, описывающую передачу импульса при любых давлениях:
$$M(P) = \frac{\pi \eta \omega R^4}{2(d + 2\zeta \lambda)}$$

Эта формула — подлинный триумф аналитической мысли. Проверим ее на прочность в наших двух пределах.

При высоких давлениях $P \to \infty$ длина свободного пробега стремится к нулю ($\lambda \to 0$). Слагаемое $2\zeta\lambda$ в знаменателе исчезает, и формула вырождается в чисто гидродинамическое решение $M \to M_2$, не зависящее от давления.

При сверхнизких давлениях $P \to 0$, напротив, $\lambda \to \infty$. Геометрическим зазором $d$ в знаменателе можно пренебречь по сравнению с гигантской длиной свободного пробега. Тогда уравнение принимает вид:
$$M(P) \approx \frac{\pi \eta \omega R^4}{4\zeta \lambda}$$

Здесь необходимо вспомнить, что макроскопическая вязкость $\eta$ содержит в себе длину свободного пробега неявно ($\eta \propto \rho \lambda \propto n \lambda$). Отношение $\eta / \lambda$ оказывается величиной, строго зависящей только от концентрации частиц $n$, а значит, линейно зависящей от давления $P$. Подстановка кинетических выражений для $\eta$ и $\lambda$ мгновенно и без потерь трансформирует эту дробь в формулу молекулярных пуль $M \propto P$.

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