Неустойчивость Рэлея

Развитие нестабильности Рэлея — Тейлора.

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

Основным параметром, определяющим скорость развития этой нестабильности является число Атвуда.

Аналитическое описание

Задача о неустойчивости Рэлея — Тейлора имеет аналитическое решение в рамках линейной теории устойчивости.

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

$\frac{\partial \vec{v}}{\partial t} + \left( \vec{v} \cdot \nabla \right) \vec{v} = - \frac{1}{\rho} \nabla P + \vec{g},$

$\operatorname{div} \vec{v} = 0.$

В дальнейшем компоненты скорости определяются как . Вполне очевидно, что равновесное решение () удовлетворяет модели, при этом из уравнения Эйлера для давления получается следующее:

$\frac{\partial P}{\partial x} = 0, \quad \frac{\partial P}{\partial y} = 0, \quad \frac{\partial P}{\partial z} = - \rho g$

Откуда определяется равновесное распределение давления (известный результат для давления столба жидкости):

$P_0 = - \rho g z.$

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

$\frac{\partial \vec{v}}{\partial t} = - \frac{1}{\rho} \nabla P,$

$\operatorname{div} \vec{v} = 0.$

Граничные условия задаются исходя из соображений равенства z-компонент скорости жидкостей 1 и 2 на границе раздела и наличия поверхностного натяжения. На верхней и нижней границах, т. к. жидкость идеальная, работают условия непротекания. Удобно принять координату границы раздела в равновесии за 0. На ней выполняется кинематическое условие

$\quad \frac{\partial \zeta}{\partial t} = w,$

и динамическое условие

$\left(P_1 - P_2\right) - \left( \rho_1 - \rho_2 \right) g \zeta = \sigma \Delta \zeta.$

Условие непротекания верхней и нижней границ:

$z=\pm h: \quad w = 0,$

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

Положим, что возмущения имеют вид:

$\vec{v}, P, \zeta \sim e^{\lambda t} e^{i \left( k_x x + k_y y \right)},$

где - скорость роста (инкремент) возмущения, - компоненты волнового вектора возмущения границы.

Из уравнения Эйлера выражается :

$\lambda w = - \frac{1}{\rho} \frac{\partial P}{\partial z},$

а условие несжимаемости даёт уравнение Лапласа для давления. В итоге, скорость течения из задачи удаётся исключить. Остаётся линейное уравнение:

$\frac{\partial^2 P}{\partial z^2} - k^2 P = 0,$

с граничными условиями:

$z=0: \quad \left( P_1 - P_2 \right) - \left( \rho_1 - \rho_2 \right) g \zeta = -\sigma k^2 \zeta,$

$z=0: \quad \frac{1}{\rho_1} \frac{\partial P_1}{\partial z} - \frac{1}{\rho_2} \frac{\partial P_2}{\partial z} = 0,$

$z=\pm h: \quad \frac{\partial P}{\partial z}=0.$

Решение уравнения Лапласа для давления:

$P_1 = C_1 \cosh k \left( h - z \right),$

$P_2 = C_2 \cosh k \left( h + z \right).$

Константы определяются из кинематического условия. Динамическое условие даёт связь между инкрементом и модулем волнового вектора

$\lambda^2 = \frac{\left( \rho_1 - \rho_2 \right)g - \sigma k^2 }{\rho_1 + \rho_2} k \tanh kh,$

откуда непосредственно следует выражение для критического волнового числа возмущений (при ):

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

В предельном случае бесконечно глубоких слоёв () наибольшая скорость роста возмущений достигается при волновом числе

В тонких слоях ():

В природе

Ярким известным проявлением неустойчивости Рэлея — Тейлора являются вымеобразные облака, а также ядерный гриб.

Литература

Лабунцов Д.А., Ягов В.В. Механика двухфазных систем. // М.: Изд-во МЭИ, 2000. - с. 143-146.
Векштейн Г.Е. Физика сплошных сред в задачах. // М.: Институт компьютерных исследований, 2002. - с. 109-111.

Ссылки

http://www.astronet.ru/db/msg/1188634

Lt304888.ru

Туристические услуги

Рекомендуем

Неустойчивость Рэлея — Тейлора

Содержание

Аналитическое описание

В природе

Литература

Ссылки