Вихри в полимерной реке: как упругость меняет порог устойчивости ✛ ✎

Оригинал: "Competing instabilities in wide-gap viscoelastic Taylor-Couette flow: Taylor vs. helical vortices"
· Alexander Proskurin
arXiv:2608.29552v1 · 2026-08-30 · CC BY-NC-ND 4.0 · собственный разбор · 3 мин
Численное моделирование показало, что в широком зазоре между цилиндрами упругость полимерного раствора делает винтовые вихри опаснее осесимметричных.

Контекст

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

Методы

Исследование выполнено методом прямого численного моделирования трёхмерного течения вязкоупругой жидкости в кольцевом зазоре. Геометрия: безразмерный радиус внутреннего цилиндра ξ = 0.049, длина расчётной области 21.3. Уравнения движения включают дополнительный член с запаздыванием, описывающий упругость полимерных молекул. Для дискретизации использован спектральный метод Фурье вдоль оси цилиндров и сетка из 545 элементов в поперечном сечении; полиномиальный порядок p = 5, 50 мод Фурье, шаг по времени 0.025. На стационарное течение накладывался белый шум конечной длительности, и отслеживалась эволюция энергии возмущения.

Результаты

Для ньютоновской жидкости (K = 0) критическое число Рейнольдса, полученное из нелинейных расчётов, составило около 265, что совпадает с линейной теорией. Для вязкоупругой жидкости с K = 0.1 критическое значение снизилось до 202. При Re = 190 возмущение затухало, при Re = 210 энергия выходила на постоянный уровень, а при Re = 250 росла. Длительная эволюция при сверхкритических числах показала формирование либо вихрей Тейлора (для K = 0), либо винтовых вихрей (для K = 0.1). Таким образом, упругость не только снижает порог устойчивости, но и меняет тип вторичного течения.

Значение

Работа стоит на отшибе среди соседей по архиву: рядом исследования вихрей в сверхтекучем гелии, электронных вихрей на нанострунах и течений в пористых средах. Отличие данной работы в том, что она впервые подтверждает прямым нелинейным моделированием предсказание линейной теории о конкуренции мод в вязкоупругом течении Куэтта–Тейлора с широким зазором. Соседние работы либо используют другие модели жидкости, либо рассматривают иные геометрии, поэтому данное исследование заполняет нишу нелинейной верификации для слабоупругих растворов.

Развитие

Естественным продолжением станет расширение диапазона параметров: увеличение упругости K, изменение радиуса внутреннего цилиндра и учёт конечной длины цилиндров. Интересно также сравнить с экспериментами на реальных полимерных растворах, где могут проявляться эффекты, не учтённые в модели Осколкова. Возможно исследование перехода к развитой турбулентности и роли упругих напряжений в формировании вторичных течений.

Влияние

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

Следующие шаги

Провести расчёты при других значениях ξ и K, чтобы построить карту режимов. Сравнить с экспериментальными данными по течению полимерных растворов в зазоре с вращающимся внутренним цилиндром.

Ключевые проблемы

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

Взгляд машины знаний

Этот раздел адресован автору работы и его коллегам. Соседние работы нашёл смысловой поиск по нашему архиву — не по совпадению слов, а по смыслу. Мы не рецензируем: мы показываем поле вокруг.

Что сделано в работе

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

В чём сила

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

Куда можно двигаться

  • Может быть интересно сопоставить полученные пороги устойчивости с подходами, где неустойчивость рассматривается как локальная потеря устойчивости среды, и проверить, насколько универсальны критические параметры для разных геометрий.опора: Когда среда теряет устойчивость arXiv:2607.10584
  • Рядом лежит работа о влиянии тонкой пластины на ранний переход к турбулентности; возможно, перенести идею о пассивных или активных включениях в вязкоупругое течение Куэтта–Тейлора и посмотреть, как они модифицируют порог и структуру вторичных течений.опора: Тонкая пластина будит турбулентность раньше срока arXiv:2601.17341
  • Интересно сравнить модель Осколкова с другими реологическими моделями, например, с моделями, учитывающими конечную растяжимость полимерных цепей, чтобы оценить, насколько вывод о доминировании винтовых мод зависит от выбора модели.опора: Три переводчика и танец вихря arXiv:2501.00037

Что можно попробовать

  • Проверить, сохраняется ли смена доминирующей моды при других значениях радиуса внутреннего цилиндра.
  • Апробировать введение в поток пассивных частиц или нитей и проследить их влияние на порог устойчивости.
  • Сравнить критические числа Рейнольдса для модели Осколкова и модели FENE-P при одинаковых параметрах упругости.

Почему это может быть важно

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

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

\sigma = 2(\nu D + \kappa \frac{\partial D}{\partial t})
Напряжение складывается из вязкой части (пропорциональной скорости деформации) и упругой части (пропорциональной скорости изменения деформации), где κ — время запаздывания.

Ключевые числа

  • критическое число Рейнольдса для K=0: 265
  • критическое число Рейнольдса для K=0.1: 202
  • безразмерный радиус внутреннего цилиндра: 0.049
  • длина расчётной области: 21.3
  • количество элементов сетки: 545
Учёные
Daniel BernoulliLeonhard EulerClaude-Louis NavierGeorge Gabriel StokesLudwig Boltzmann
Понятия
число Рейнольдса уравнения Навье — Стокса неустойчивость Кельвина–Гельмгольца уравнение Пуазёйля вязкость неньютоновская жидкость турбулентная коагуляция вихри в сверхтекучей жидкости неустойчивость уравнение Фоккера — Планка турбулентность вычислительная гидродинамика скорость препятствие метастабильное состояние теорема Колмогорова — Арнольда — Мозера турбулентное перемешивание гармонический осциллятор вихревая решётка завихрённость