Переходный процесс подъёма жидкости по стержню в жидкости Олдройда-Б

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

• ID статьи: 2501.01016
• Название: Transient rod-climbing in an Oldroyd-B fluid
• Авторы: Tachin Ruangkriengsin, Rodolfo Brandão, Katie Wu, Jonghyun Hwang, Evgeniy Boyko, Howard A. Stone
• Классификация: physics.flu-dyn (гидродинамика)
• Дата публикации: 2 января 2025 г. (отправлено на arXiv)
• Журнал: Journal of Fluid Mechanics (в процессе рецензирования)
• Ссылка на статью: https://arxiv.org/abs/2501.01016

Аннотация

В данной работе исследуется эффект Вайссенберга (явление подъёма жидкости по стержню) в неньютоновских жидкостях, то есть подъём границы раздела жидкости вдоль вращающегося стержня. Авторы используют модель Олдройда-Б и применяют асимптотический анализ с безразмерными переменными для изучения зависящей от времени высоты границы раздела. Исследование показывает, что при малых деформациях управляющие уравнения содержат временной пограничный слой, в котором граница раздела быстро поднимается на коротких временных масштабах, а затем постепенно приближается к стационарному состоянию. Путём введения растянутой временной переменной авторы выводят переходные поля скорости и соответствующие профили границы раздела, а также переинтегрируют малые, но конечные инерционные эффекты для исследования их взаимодействия с вязкоупругостью.

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

Определение проблемы

Явление подъёма жидкости по стержню (эффект Вайссенберга) является одним из наиболее характерных поведений неньютоновских жидкостей. Когда стержень вращается в вязкоупругой жидкости, поверхность жидкости поднимается вдоль стержня. Это явление всё ещё недостаточно изучено с теоретической точки зрения.

Значимость

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

Ограничения существующих методов

1. Теоретические ограничения: Предыдущие исследования сосредоточены главным образом на стационарном поведении с использованием моделей жидкостей второго порядка
2. Проблемы модели: Модели жидкостей второго порядка действительны только при малых числах Дебора или Вайссенберга и численно неустойчивы
3. Временная зависимость: Отсутствует теоретический анализ переходной динамики

Мотивация исследования

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

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

1. Первый теоретический анализ: Предоставляет первый полный переходный теоретический анализ явления подъёма жидкости по стержню в жидкости Олдройда-Б
2. Теория временного пограничного слоя: Обнаруживает и анализирует структуру временного пограничного слоя, раскрывая двухэтапный процесс быстрого подъёма и медленного приближения к стационарному состоянию
3. Безразмерная аналитическая база: Устанавливает систематический метод безразмерного анализа на основе числа Вайссенберга
4. Исследование инерционных эффектов: Впервые исследует взаимодействие малых, но конечных инерционных эффектов с вязкоупругостью
5. Критерий подъёма: Предлагает теоретический критерий для определения условий возникновения подъёма жидкости: $Re < 4\beta_p Wi$

Подробное описание методов

Определение задачи

Исследуется временная эволюция высоты свободной поверхности жидкости $h(r,t)$ при вращении бесконечно длинного цилиндрического стержня (радиус $a$) с угловой скоростью $\Omega$ в вязкоупругом полимерном растворе.

Математическая модель

Управляющие уравнения

Используется конститутивная модель Олдройда-Б:

• Уравнение эволюции тензора конфигурации: $\frac{\partial\mathbf{A}}{\partial t} + \mathbf{u} \cdot \nabla\mathbf{A} - \mathbf{A} \cdot \nabla\mathbf{u} - (\nabla\mathbf{u})^T \cdot \mathbf{A} = -\frac{1}{\lambda}(\mathbf{A} - \mathbf{I})$
• Тензор напряжений: $\boldsymbol{\sigma} = -p\mathbf{I} + 2\mu_s\mathbf{E} + G(\mathbf{A} - \mathbf{I})$
• Уравнение количества движения: $\nabla \cdot \boldsymbol{\sigma} = \rho g\mathbf{e}_z$

Обезразмеривание

Введены ключевые безразмерные параметры:

• Число Вайссенберга: $Wi = \lambda\Omega$
• Параметр гравитации: $G = \frac{\rho ga}{(\mu_s + \mu_p)\Omega}$
• Число Рейнольдса: $Re = \frac{\rho a^2\Omega}{\mu_s + \mu_p}$

Характеристическое масштабирование

Авторы выявляют ключевое характеристическое масштабирование: $G = Wi \cdot \mathcal{G}$, где $\mathcal{G} = O(1)$, что упрощает двухпараметрическое асимптотическое разложение до однопараметрического.

Технические инновации

1. Разделение временных масштабов

Обнаруживается наличие двух временных масштабов:

• Внешний временной масштаб: $T = O(1)$, соответствует стационарному поведению
• Внутренний временной масштаб: $\tau = T/Wi$, соответствует переходному пограничному слою

2. Метод возмущения области

Для малых деформаций границы раздела граничные условия разлагаются в ряд Тейлора от $Z = H(R,T)$ к $Z = 0$, что упрощает обработку граничных условий свободной поверхности.

3. Преобразование Вебера

При рассмотрении инерционных эффектов применяется преобразование Вебера для обработки оператора Бесселя дифференцирования на полубесконечной области, преобразуя уравнение в частных производных в обыкновенное дифференциальное уравнение.

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

Теоретическая аналитическая база

Данная работа является чисто теоретическим исследованием, использующим метод асимптотического анализа:

1. Малый параметр: $Wi \ll 1$ (слабая вязкоупругость)
2. Диапазон гравитации: $G \gg Wi^{-1}$ (сильные гравитационные эффекты)
3. Предположение о малых деформациях: Деформация границы раздела относительно небольшая

Стратегия анализа

1. Раздел 4: Внешний анализ – стационарный профиль границы раздела
2. Раздел 5: Внутренний анализ – переходный профиль границы раздела
3. Раздел 6: Анализ с включением малых инерционных эффектов

Результаты исследования

Основные результаты

1. Стационарный профиль границы раздела

$H_{steady}(R) = \frac{2\beta_p Wi^2}{\mathcal{G}R^4}$

Это согласуется с результатами, полученными Joseph & Fosdick (1973) с использованием модели жидкости второго порядка.

2. Переходный профиль границы раздела

$H_{transient}(R,\tau) = \frac{2\beta_p(1-(\tau+1)e^{-\tau})Wi^2}{\mathcal{G}R^4}$

где $\tau = T/Wi$ – растянутая временная переменная.

3. Критерий подъёма

При включении инерционных эффектов условие возникновения подъёма жидкости по стержню имеет вид: $4El\beta_p > 1$ где $El = Wi/Re$ – число упругости.

Характеристики временной эволюции

1. Этап быстрого подъёма: Граница раздела быстро поднимается в течение времени $\tau \approx 6$
2. Приближение к стационарному состоянию: Переходное решение сходится к стационарному решению в пределе больших времён
3. Немонотонное поведение: При включении инерции высота границы раздела может сначала подняться, а затем опуститься

Ключевые выводы

1. Окружное напряжение: Окружное напряжение, возникающее в результате растяжения полимера, является основным механизмом, управляющим подъёмом
2. Конкуренция инерции: Центробежная сила, вызванная инерцией, конкурирует с окружным напряжением вязкоупругости
3. Структура пограничного слоя: Наличие временного пограничного слоя объясняет двухэтапное поведение быстрого переходного процесса и медленного стационарного состояния

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

Историческое развитие

1. Serrin (1959): Первая теоретическая попытка с использованием модели жидкости второго порядка
2. Giesekus (1961): Более общие конститутивные уравнения с пренебрежением инерционными эффектами
3. Joseph & Fosdick (1973): Систематическая база метода возмущения области
4. Oldroyd (1950, 1951): Ранние приложения модели Олдройда-Б

Преимущества данной работы

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

Заключение и обсуждение

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

1. Структура пограничного слоя: Явление подъёма жидкости по стержню демонстрирует явный временной пограничный слой, в котором граница раздела быстро поднимается, а затем медленно приближается к стационарному состоянию
2. Разделение масштабов: Переходный временной масштаб составляет $O(Wi)$, стационарный временной масштаб – $O(1)$
3. Механизм подъёма: Окружное напряжение, возникающее в результате растяжения полимера, является основным механизмом, управляющим подъёмом
4. Влияние инерции: Малые инерционные эффекты могут значительно влиять на поведение подъёма, существует критическое число упругости

Ограничения

1. Предположение о малых деформациях: Применимо только при малых деформациях границы раздела
2. Слабая вязкоупругость: Ограничено диапазоном $Wi \ll 1$
3. Сильная гравитация: Требует условие $G \gg Wi^{-1}$
4. Упрощённая модель: Пренебрегает поверхностным натяжением и эффектами контактного угла

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

1. Высокие числа Вайссенберга: Исследование поведения при $Wi \gg 1$
2. Сложные конститутивные модели: Расширение на модели FENE-P, Giesekus и другие
3. Конечные деформации: Рассмотрение случаев больших деформаций границы раздела
4. Экспериментальная верификация: Сравнение с экспериментальными данными

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

Достоинства

1. Теоретическая строгость: Математические выводы точны, методы асимптотического анализа применены надлежащим образом
2. Физические инсайты: Раскрывает структуру временного пограничного слоя и механизм подъёма
3. Методологические инновации: Искусное сочетание преобразования Вебера и метода возмущения области
4. Практическая ценность: Обеспечивает теоретическую основу для реологических измерений

Недостатки

1. Ограниченный диапазон параметров: Применимость ограничена относительно узким диапазоном параметров
2. Отсутствие экспериментального сравнения: Нет прямого сравнения с экспериментальными данными
3. Численная верификация: Некоторые результаты представлены только численными расчётами без аналитических выражений

Влияние

1. Научный вклад: Заполняет пробел в теории переходного подъёма жидкости по стержню
2. Методологическая ценность: Предложенные методы анализа могут быть распространены на другие задачи с комплексными жидкостями
3. Инженерные приложения: Обеспечивает теоретическую поддержку для проектирования реологических приборов и характеризации материалов

Области применения

1. Разбавленные полимерные растворы: Жидкости Бога и другие вязкоупругие жидкости с постоянной вязкостью при сдвиге
2. Низкие скорости вращения: Условия экспериментов с малыми числами Вайссенберга
3. Реологические измерения: Определение коэффициента нормального напряжения и времени релаксации

Список литературы

1. Joseph, D.D. & Fosdick, R.L. (1973). The free surface on a liquid between cylinders rotating at different speeds. Arch. Rat. Mech. Anal., 49, 321-380.
2. Oldroyd, J.G. (1950). On the formulation of rheological equations of state. Proc. R. Soc. A, 200, 523-541.
3. Bird, R.B., Armstrong, R.C. & Hassager, O. (1987). Dynamics of Polymeric Liquids, Vol. 1, 2nd edn. John Wiley and Sons.

Общая оценка: Это высококачественная теоретическая работа по гидродинамике, впервые систематически анализирующая переходное поведение явления подъёма жидкости по стержню в жидкости Олдройда-Б. Математические выводы строги, физические инсайты глубоки, методы инновационны, и работа вносит важный вклад в теорию сложных жидкостей и приложения реологии. Несмотря на ограничения диапазона параметров, теоретическая база и методы анализа имеют хороший потенциал для обобщения.