Неустойчивость Рэлея-Тейлора: тяжёлое поверх лёгкого

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

Суть явления

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

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

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

Число Атвуда

Ключевой безразмерный параметр задачи - число Атвуда, описывающее контраст плотностей:

$$A = \frac{\rho_2 - \rho_1}{\rho_2 + \rho_1}$$

Здесь $\rho_2$ - плотность верхней (тяжёлой) жидкости, $\rho_1$ - нижней (лёгкой). Число Атвуда меняется от 0 до 1:

• $A \to 0$ - плотности почти равны (например, две смешивающиеся жидкости), неустойчивость развивается медленно и почти симметрично: пузыри и струи похожи.
• $A \to 1$ - тяжёлая жидкость много плотнее лёгкой (вода над воздухом, $A \approx 1$), процесс бурный и асимметричный: тяжёлые струи проваливаются узкими «шипами», лёгкие пузыри всплывают широкими «шапками».

При $A < 0$ тяжёлое находится снизу - конфигурация устойчива, неустойчивости нет.

Рассчитать для своей задачи

Инкремент роста

В линейной стадии (амплитуда мала по сравнению с длиной волны) поверхность раскладывают в гармоники $\propto e^{i k x + n t}$, где $k = 2\pi/\lambda$ - волновое число возмущения, $n$ - инкремент роста. Для идеальной (невязкой) жидкости без поверхностного натяжения классический результат Тейлора:

$$n = \sqrt{A\, g\, k}$$

Возмущение растёт экспоненциально: $a(t) = a_0\, e^{n t}$. Важное следствие - чем короче длина волны (больше $k$), тем быстрее рост: $n \propto \sqrt{k}$. Поэтому в идеальной модели мелкомасштабные возмущения растут стремительнее всего, а спектр стремится к коротким волнам - без стабилизирующего механизма задача была бы плохо поставлена.

Стабилизирует коротковолновую часть поверхностное натяжение $\sigma$. С его учётом инкремент:

$$n = \sqrt{A\, g\, k - \frac{\sigma\, k^3}{\rho_1 + \rho_2}}$$

Подкоренное выражение положительно (неустойчивость есть) только при $k < k_c$, где критическое волновое число

$$k_c = \sqrt{\frac{A\, g\, (\rho_1 + \rho_2)}{\sigma}}$$

Возмущения короче критической длины волны $\lambda_c = 2\pi/k_c$ гасятся натяжением и не растут. Быстрее всего растёт мода с $k_{max} = k_c/\sqrt{3}$ - именно её масштаб задаёт характерный размер первых «грибов».

Стадии развития

Реальный распад границы проходит несколько качественно разных фаз:

1. Линейная стадия. Амплитуда мала, моды растут независимо и экспоненциально по закону $a_0 e^{n t}$. Работает приведённая выше теория.
2. Нелинейная стадия. Когда амплитуда сравнивается с длиной волны, форма перестаёт быть синусоидой: лёгкие пузыри всплывают с почти постоянной скоростью, тяжёлые струи («спайки») заостряются и ускоряются.
3. Грибовидные структуры. На боках струй включается вторичная неустойчивость Кельвина-Гельмгольца - концы заворачиваются в характерные «шляпки грибов» (вихревые валики).
4. Турбулентное перемешивание. Структуры разных масштабов взаимодействуют, граница превращается в развивающийся перемешанный слой, толщина которого растёт как $h \sim \alpha A g t^2$ - квадратично по времени.

Рассчитать для своей задачи

Где встречается

Неустойчивость Рэлея-Тейлора - не лабораторная экзотика, а массовое явление в природе и технике:

• Астрофизика. При взрыве сверхновой ударная волна гонит плотные внешние слои внутрь, и граница распадается на знаменитые «пальцы» - их прекрасно видно в Крабовидной туманности. Тот же механизм формирует структуру остатков сверхновых.
• Инерциальный термоядерный синтез. При сжатии мишени плотная оболочка ускоряется лёгким горячим газом - возмущения растут и могут разрушить симметрию сжатия. Подавление этой неустойчивости - одна из центральных инженерных задач ICF.
• Геофизика. Поднимающиеся плюмы мантии, соляные диапиры (лёгкая соль всплывает сквозь более плотные осадочные породы), грибовидные облака взрывов - всё это проявления того же эффекта.
• Повседневность. Капля густого сиропа, тонущая в воде; «лавовая лампа»; перевёрнутый стакан воды.

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

Связь с неустойчивостью Рихтмайера-Мешкова

Близкий родственник - неустойчивость Рихтмайера-Мешкова, когда границу плотностей пересекает не постоянное ускорение, а импульсное (ударная волна). Это можно считать предельным случаем Рэлея-Тейлора с очень коротким, но сильным ускорением. Отличие: после прохождения удара возмущение растёт не экспоненциально, а линейно по времени, и неустойчивость развивается даже при «неправильном» (стабилизирующем для Рэлея-Тейлора) знаке градиента плотности. Оба эффекта обычно работают вместе в задачах термоядерного синтеза и астрофизики.

Частые ошибки

• «Неустойчивость возникает при любом контакте разноплотных жидкостей». Нет. Нужна именно конфигурация «тяжёлое над лёгким» с ускорением от тяжёлого к лёгкому. Лёгкое сверху (масло на воде) - устойчиво.
• «Длинные волны растут быстрее». Наоборот: в идеальной модели $n \propto \sqrt{k}$, короткие волны растут быстрее. Поверхностное натяжение лишь обрезает совсем короткие, оставляя промежуточный максимум.
• «Поверхностное натяжение всегда мешает неустойчивости». Оно стабилизирует только коротковолновую часть; длинные волны ($k < k_c$) растут независимо от натяжения.
• «Число Атвуда - это отношение плотностей». Нет, это нормированная разность: $A = (\rho_2 - \rho_1)/(\rho_2 + \rho_1)$, а не $\rho_2/\rho_1$.
• «Рэлея-Тейлора и Рихтмайера-Мешкова - одно и то же». Механизм родственный, но первая вызвана постоянным ускорением (рост экспоненциальный), вторая - импульсным ударом (рост линейный).

FAQ

Почему граница распадается именно на грибы? Тяжёлые струи проваливаются вниз заострёнными «шипами», а на их боках, где жидкости движутся друг относительно друга, включается вторичная неустойчивость Кельвина-Гельмгольца. Она заворачивает концы струй в вихревые валики - отсюда форма гриба со шляпкой.

От чего зависит скорость роста неустойчивости? От числа Атвуда $A$ (контраст плотностей), ускорения $g$ и волнового числа $k$ возмущения: $n = \sqrt{Agk}$. Чем больше контраст плотностей и ускорение и чем короче волна - тем быстрее рост (до предела, поставленного поверхностным натяжением и вязкостью).

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

Коротко

Неустойчивость Рэлея-Тейлора развивается на границе, где тяжёлая жидкость лежит над лёгкой в поле ускорения: любое возмущение нарастает, понижая потенциальную энергию системы. Скорость роста задаёт инкремент $n = \sqrt{Agk}$, где $A = (\rho_2 - \rho_1)/(\rho_2 + \rho_1)$ - число Атвуда; короткие волны растут быстрее, пока их не обрежет поверхностное натяжение (критическое $k_c = \sqrt{Ag(\rho_1+\rho_2)/\sigma}$). Граница проходит линейную стадию, нелинейные струи, грибовидные структуры и турбулентное перемешивание со слоем $h \sim Agt^2$. Явление управляет взрывами сверхновых, термоядерным синтезом, мантийными плюмами и тонущей каплей сиропа.