Смешивание в вязкоупругих жидкостях с использованием упругой турбулентности

Основная информация

• ID статьи: 2409.06391
• Название: Mixing in viscoelastic fluids using elastic turbulence
• Авторы: Reinier van Buel, Holger Stark (Technische Universität Berlin)
• Классификация: physics.flu-dyn
• Дата публикации: сентябрь 2024
• Ссылка на статью: https://arxiv.org/abs/2409.06391

Аннотация

В данной работе исследуется влияние упругой турбулентности на смешивание поля скалярной концентрации в вязкоупругих жидкостях в двумерной геометрии Тейлора-Куэтта с использованием численного решения модели Олдройда-B. Состояние потока определяется параметром порядка вторичного потока, указывающим на область упругой турбулентности. При снижении числа Вейссенберга ниже $Wi_c$ с начального турбулентного состояния наблюдается слабохаотический поток. Конвекция как в турбулентном, так и в слабохаотическом состояниях индуцирует смешивание, что иллюстрируется временной эволюцией стандартного отклонения концентрации растворённого вещества относительно однородного распределения. В частности, смешивание в состоянии упругой турбулентности является интенсивным и количественно оценивается через скорость смешивания, время смешивания и эффективность смешивания. Все три величины подчиняются степенным законам. Важно отметить, что параметр порядка сильно коррелирует со скоростью смешивания и, следовательно, является хорошим индикатором смешивания в жидкости.

Научный контекст и мотивация

Проблемная постановка

1. Проблема смешивания при низких числах Рейнольдса: В ньютоновских жидкостях при низких числах Рейнольдса течение является ламинарным, а процесс смешивания протекает медленно и контролируется в основном молекулярной диффузией. Это представляет значительное ограничение в микрофлюидных устройствах.
2. Открытие упругой турбулентности: Добавление высокомолекулярных полимеров в жидкость делает её вязкоупругой, позволяя ей проявлять явления упругой турбулентности даже при пренебрежимо малых инерционных эффектах.
3. Значимость смешивания в микромасштабе: В устройствах типа "лаборатория на кристалле" (lab-on-a-chip) эффективное смешивание жидкостей критически важно, особенно в микронном масштабе.

Научная мотивация

• Исследование потенциала упругой турбулентности для повышения эффективности смешивания
• Понимание физической природы механизмов смешивания в вязкоупругих жидкостях
• Предоставление теоретического руководства для микрофлюидных приложений

Основные вклады

1. Определение природы перехода: Идентификация перехода от ламинарного течения к упругой турбулентности как субкритического перехода и обнаружение существования состояния слабохаотического потока
2. Установление методов характеризации смешивания: Введение трёх количественных показателей — скорости смешивания, времени смешивания и эффективности смешивания, все подчиняющихся степенным законам
3. Обнаружение сильной корреляции между параметром порядка и смешиванием: Доказательство того, что параметр порядка вторичного потока сильно коррелирует со скоростью смешивания и может служить хорошим индикатором степени смешивания жидкости
4. Верификация механизма смешивания Бэтчелора: Наблюдение экспоненциального затухания в состоянии упругой турбулентности, соответствующего области смешивания Бэтчелора, согласующегося с экспериментальными результатами

Методология

Геометрическая установка и управляющие уравнения

Исследование использует двумерную геометрию Тейлора-Куэтта с двумя концентрическими цилиндрами:

• Радиус внешнего цилиндра: $r_o = 10\mu m$
• Радиус внутреннего цилиндра: $r_i = 2.5\mu m$
• Внешний цилиндр вращается против часовой стрелки с угловой скоростью $\Omega = 2\pi s^{-1}$

Движение жидкости описывается обобщённым уравнением Навье-Стокса: $\rho\left(\frac{\partial \mathbf{u}}{\partial t} + \mathbf{u} \cdot \nabla \mathbf{u}\right) = -\nabla p + \eta_s \nabla^2 \mathbf{u} + \nabla \cdot \boldsymbol{\tau}$

Определяющая модель Олдройда-B

Тензор полимерного напряжения $\boldsymbol{\tau}$ описывается моделью Олдройда-B: $\boldsymbol{\tau} + \lambda \stackrel{\nabla}{\boldsymbol{\tau}} = \eta_p[\nabla \otimes \mathbf{u} + (\nabla \otimes \mathbf{u})^T]$

где $\stackrel{\nabla}{\boldsymbol{\tau}}$ — верхняя конвективная производная.

Моделирование смешивания

Эволюция поля скалярной концентрации $c(\mathbf{r}, t)$ управляется уравнением конвекции-диффузии: $\frac{\partial c}{\partial t} + (\mathbf{u} \cdot \nabla)c = D\nabla^2 c$

В данном исследовании принято $D = 0$, что пренебрегает молекулярной диффузией.

Ключевые безразмерные параметры

• Число Рейнольдса: $Re = \rho\Omega r_o^2/\eta_s = 0.628 \times 10^{-3}$
• Число Вейссенберга: $Wi = \lambda\Omega$
• Отношение вязкостей: $\beta = \eta_p/\eta_s = 1.5$

Экспериментальная установка

Численный метод

• Использование открытого решателя конечных объёмов OpenFOAM®
• Специализированный решатель RheoTool для обработки вязкоупругих течений
• Сетка: радиальная сетка со 100 элементами и 120 элементами в окружном направлении
• Временной шаг: $\delta t = 10^{-5}s$

Определение параметра порядка

Интенсивность вторичного потока определяется как параметр порядка: $\Phi = \sigma = \sqrt{\langle[\mathbf{u}(\mathbf{r},t) - \mathbf{u}_0(\mathbf{r})]^2\rangle}/u_0^{\max}$

Показатели характеризации смешивания

Нормализованное стандартное отклонение: $\Delta c = \sqrt{\frac{\langle[c(\mathbf{r},t) - c_\infty]^2\rangle}{\langle[c(\mathbf{r},0) - c_\infty]^2\rangle}}$

где $c_\infty$ — концентрация в состоянии полного смешивания.

Экспериментальные результаты

Переходы состояния потока

1. Субкритический переход: Обнаружено, что переход от ламинарного течения к упругой турбулентности имеет субкритический характер с явлением гистерезиса
2. Критические числа Вейссенберга:
   • Восходящая ветвь: $Wi_c = 9.95$
   • Нисходящая ветвь: $Wi_c^I = 9.26$
3. Слабохаотическое состояние: Наблюдается слабохаотический поток при $Wi_c^I$ ниже определённого значения

Анализ характеристик смешивания

Состояние слабохаотического потока ($5.03 < Wi \leq 8.17$)

• Степенной закон затухания стандартного отклонения с показателем, близким к $-1/6$
• Энергетический спектр демонстрирует степенной закон с показателем $-5$ или выше

Состояние упругой турбулентности ($Wi \geq 9.42$)

• Начальное степенное затухание (показатель около $1/5$), переходящее в экспоненциальное затухание
• Форма экспоненциального затухания: $\Delta c \sim \exp(-\alpha_{mix}t)$
• Энергетический спектр демонстрирует характеристики спектра Бэтчелора (масштабирование $k^{-1}$)

Степенные законы

1. Время смешивания: $t_{mix} \sim (Wi - Wi_c^I)^{-0.52}$
2. Концентрация в области II: $c_{II}/c_\infty \sim (Wi - Wi_c^I)^{0.28}$
3. Масштабное соотношение: $c_{II} \sim t_{mix}^{-1/2}$

Корреляция между параметром порядка и скоростью смешивания

Экспериментальные результаты показывают сильную корреляцию между параметром порядка $\Phi_d$ и скоростью смешивания $\alpha_{mix}$, что указывает на то, что усиление колебаний скорости непосредственно повышает эффективность смешивания в жидкости.

Связанные работы

Исследования упругой турбулентности

• Первоначальное открытие упругой турбулентности в растворах полимеров Гройсманом и Штайнбергом
• Анализ линейной и нелинейной устойчивости в геометрии Тейлора-Куэтта
• Наблюдение спиральных и полосчатых структур потока в трёхмерном численном моделировании

Исследования усиления смешивания

• Смешивание при растяжении в четырёхвалковой системе
• Упругая неустойчивость в геометрии крестообразного канала
• Усиление конвективного теплообмена в змеевидных каналах

Выводы и обсуждение

Основные выводы

1. Природа перехода: Подтверждение субкритического характера перехода упругой турбулентности в течении Тейлора-Куэтта
2. Механизм смешивания: Доказательство того, что упругая турбулентность обеспечивает эффективное смешивание при экстремально низких числах Рейнольдса
3. Область Бэтчелора: Наблюдение характеристик смешивания Бэтчелора в состоянии упругой турбулентности
4. Предсказательный индикатор: Параметр порядка может служить надёжным предсказательным индикатором эффективности смешивания

Ограничения

1. Двумерное ограничение: Исследование ограничено двумерной геометрией; трёхмерные эффекты могут дать различные результаты
2. Зависимость от сетки: Явление гистерезиса исчезает при более мелкой сетке, что указывает на влияние численной точности
3. Чувствительность к начальным условиям: Характеристики смешивания зависят от выбора начального распределения концентрации

Направления будущих исследований

1. Влияние упругих волн: Исследование влияния упругих волн на эффективность смешивания, скорость смешивания и время смешивания
2. Трёхмерное расширение: Верификация обнаруженных степенных законов в трёхмерной геометрии
3. Практические приложения: Исследование потенциала практического применения в микрофлюидных устройствах

Глубокая оценка

Преимущества

1. Систематическое исследование: Комплексный анализ характеристик смешивания в процессе перехода от ламинарного течения к упругой турбулентности
2. Количественная характеризация: Введение нескольких количественных показателей для характеризации степени и эффективности смешивания
3. Теоретическая верификация: Численные результаты согласуются с теорией смешивания Бэтчелора и экспериментальными наблюдениями
4. Практическая ценность: Предоставление важного теоретического руководства для микрофлюидных приложений

Недостатки

1. Проблемы численной точности: Зависимость явления гистерезиса от сетки указывает на возможные ограничения численного метода
2. Понимание механизма: Недостаточное объяснение физического механизма степенного затухания в состоянии слабохаотического потока
3. Диапазон параметров: Относительно ограниченный диапазон исследуемых параметров, отсутствие более широкого исследования пространства параметров

Влияние

1. Теоретический вклад: Предоставление важных сведений для понимания механизмов смешивания в вязкоупругих жидкостях
2. Перспективы применения: Значительное руководство для проектирования микрофлюидных устройств
3. Методологическая ценность: Установленные методы характеризации могут быть распространены на другие геометрии потока

Применимые сценарии

• Оптимизация смешивания в микрофлюидных устройствах
• Контроль смешивания в процессах переработки полимеров
• Смешивание в микромасштабе для биомедицинских приложений

Библиография

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