Модель решеточного метода Больцмана для динамики сжимаемых неидеальных жидкостей

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

• ID статьи: 2510.14712
• Название: Lattice Boltzmann model for non-ideal compressible fluid dynamics
• Авторы: S. A. Hosseini, M. Feinberg, I. V. Karlin (ETH Zurich)
• Классификация: physics.flu-dyn
• Журнал публикации: Journal of Fluid Mechanics
• Ссылка на статью: https://arxiv.org/abs/2510.14712

Аннотация

В данной работе предложен метод решеточного Больцмана для моделирования течений сжимаемых неидеальных жидкостей. Метод использует решетку первого соседства и вводит согласованный набор поправочных членов через квазиравновесный аттрактор, обеспечивая положительную определённость и галилеевскую инвариантность диссипативной функции Навье-Стокса. Такая конструкция исключает необходимость в расширенных шаблонах или временной регуляризации, сохраняя при этом численную устойчивость и термодинамическую согласованность в широком диапазоне состояний течения. Полученная модель точно воспроизводит гидродинамику на уровнях Эйлера и Навье-Стокса. В качестве строгой верификации авторы впервые в рамках метода решеточного Больцмана продемонстрировали количественно точное моделирование взаимодействия капли с ударной волной при числе Маха до 1,47.

Исследовательский контекст и мотивация

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

Динамика сжимаемых неидеальных жидкостей — это развивающееся и быстро растущее направление гидромеханики, обусловленное появлением методов и технологий, работающих при близкритических, трансритических и сверхкритических состояниях. К таким жидкостям относятся:

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

Значимость

1. Потребности промышленных приложений: органические циклы Ренкина, сверхкритические CO₂-турбины и другие энергетические технологии
2. Дефицит экспериментальных данных: сложность экспериментальных установок затрудняет получение данных
3. Потребность в численных инструментах: необходимы согласованные и эффективные численные инструменты для понимания сложных физических явлений

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

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

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

1. Предложена новая двухраспределительная модель LBM: использует классическую решетку первого соседства (D3Q27) для моделирования неидеальных сжимаемых течений
2. Введен механизм квазиравновесного аттрактора: поправочные члены fᵢ и gᵢ обеспечивают положительную определённость и галилеевскую инвариантность диссипативной функции
3. Реализовано независимое управление объёмной вязкостью: решена проблема возможной отрицательной объёмной вязкости в структуре BGK
4. Впервые реализовано моделирование взаимодействия капли с ударной волной при высоком числе Маха: достигнута количественно точная модель при Ma=1,47 в рамках метода LBM
5. Сохранена термодинамическая согласованность: поддерживается численная устойчивость в широком диапазоне состояний течения

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

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

Моделирование сжимаемых неидеальных жидкостей, удовлетворяющих системе уравнений Навье-Стокса-Кортевега (NSK):

• Сохранение массы: ∂ₜρ + ∇·(ρu) = 0
• Сохранение импульса: ∂ₜ(ρu) + ∇·(ρu⊗u) + ∇P + ∇·T^NS + ∇·T^K = 0
• Сохранение энергии: ∂ₜ(ρE) + ∇·((ρE+P)u) + u·(∇·T^K) + ∇·(u·T^NS) + ∇·q = 0

где T^K — тензор поверхностного натяжения Кортевега, T^NS — тензор вязких напряжений.

Архитектура модели

1. Двухраспределительная конструкция

Используются два набора функций распределения {fᵢ, gᵢ}, i=1,...,27:

• fᵢ: обработка сохранения массы и импульса
• gᵢ: обработка сохранения энергии

2. Уравнение решеточного Больцмана

fᵢ(x+cᵢδt, t+δt) = fᵢ + 2β(fᵢᵉᵠ - fᵢ) + (1-β)(fᵢ* - fᵢᵉᵠ)
gᵢ(x+cᵢδt, t+δt) = gᵢ + 2β(gᵢᵉᵠ - gᵢ) + (1-β)(gᵢ* - gᵢᵉᵠ)

где β = Pδt/(2μ + Pδt) — параметр релаксации.

3. Функции равновесного распределения

Используется произведение формы:

fᵢᵉᵠ = ρ ∏_α Ψᵢα(uα, P/ρ + uα²)

где Ψᵢα — базовые функции:

Ψᵢα(ξα, ζαα) = 1 - cᵢα² + ½[(3cᵢα² - 2)ζαα + cᵢαξα]

4. Поправочные члены квазиравновесия

Поправки для распределения f:

• Силовой член: Fα = -κρ∂α∇²ρ (сила Кортевега)
• Поправка напряжений: Φαα включает поправку диагонального тензора напряжений и независимое управление объёмной вязкостью

Поправки для распределения g:

• Смещённая скорость: u* = u + δt F/ρ
• Смещённая температура: T* = T - δt²F·F/(2ρ²cᵥ)
• Поправка теплового потока: учитывает эффекты градиента энтальпии для неидеальных газов

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

1. Механизм квазиравновесного аттрактора

Введением fᵢ и gᵢ система релаксирует не просто к локальному равновесию, а к модифицированному квазиравновесному состоянию, что обеспечивает:

• Корректную макроскопическую динамику
• Численную устойчивость
• Галилеевскую инвариантность

2. Независимое управление объёмной вязкостью

В традиционной модели BGK объёмная вязкость η' = ((D+2)/D - ρcs²/P)μ может быть отрицательной. Независимое управление реализуется через поправочный член Φ':

Φ' = P((D+2)/D - ρcs²/P - η/μ)(∇·u)

3. Термодинамическая согласованность

Тщательно разработанные смещённая температура T* и поправки энергии обеспечивают термодинамическую согласованность при нелинейных уравнениях состояния.

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

Верификационные случаи

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

Параметры вещества

Основная верификация проводилась на азоте (N₂):

• Критическое давление: Pc = 3,4×10⁶ Па
• Критическая плотность: ρc = 241,96 кг/м³
• Критическая температура: Tc = 126,2 K
• R/cv = 0,4

Уравнение состояния

Используется уравнение состояния Ван-дер-Ваальса:

P(ρ,T) = ρRT/(1-bρ) - aρ²

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

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

1. Верификация базовых физических величин

• Скорость звука: на ветвях насыщенной жидкости и пара идеально совпадает с аналитическим решением в диапазоне Tr∈0,7,1
• Сдвиговая вязкость: сохраняет галилеевскую инвариантность при различных числах Маха с ошибкой <1%
• Тепловой поток Куэтта: согласуется с аналитическим решением при Pr∈{0,6,1,2,4,9} и Ma∈{0,8,1,2,1,6}

2. Свойства многофазного течения

• Плотности коэкзистенции: полностью совпадают с конструкцией равных площадей Максвелла в диапазоне Tr∈0,3,0,99
• Сходимость интерфейса: демонстрирует хорошие сходимостные свойства при различных разрешениях сетки

3. Неидеальная ударная труба

Успешное моделирование трёх различных конфигураций неидеальной ударной трубы:

• Случай I: разреженная ударная волна в области Γ<0
• Случай II: классическое поведение волны сжатия
• Случай III: смешанная разреженная волна, динамически пересекающая линию Γ=0

Все случаи показывают отличное совпадение с эталонными данными.

4. Взаимодействие ударной волны с каплей

• Впервые реализовано моделирование при высоком числе Маха Ma=1,47 в рамках метода LBM
• Точно захвачены сложные волновые структуры: падающая ударная волна, прошедшая волна, отражённая ударная волна, стержень Маха
• Эволюция деформации капли высоко согласуется с экспериментальными данными

Абляционные исследования

Верификация сходимости модели при различных разрешениях сетки на интерфейсе:

• Четыре уровня разрешения: δx от 5 мкм до 0,1 мкм
• Демонстрирует хорошие свойства сходимости по сетке

Ключевые находки

1. Численная устойчивость: сохраняется при высоких числах Маха и сильных неидеальных эффектах
2. Вычислительная эффективность: использование стандартной решетки первого соседства обеспечивает контролируемые вычислительные затраты
3. Физическая точность: точно воспроизводит сложные физические явления неидеальных жидкостей

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

Традиционные многофазные модели LBM

1. Модель цветовых градиентов (Gunstensen et al. 1991)
2. Псевдопотенциальная модель (Shan & Chen 1993)
3. Модель свободной энергии (Swift et al. 1996)

Развитие сжимаемого LBM

• He & Doolen (2002): первые попытки сжимаемых неидеальных течений
• Существующие модели в основном ориентированы на приложения кипения
• Отсутствие систематических моделей для сжимаемых неидеальных течений

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

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

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

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

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

Ограничения

1. Модель с единственным временем релаксации: в настоящее время используется простая модель столкновений BGK, которая может быть недостаточно робастна для моделирования турбулентности
2. Ограничения уравнения состояния: основная верификация проведена для уравнения Ван-дер-Ваальса, требуется дальнейшая верификация для других сложных уравнений состояния
3. Ограниченная трёхмерная верификация: большинство верификаций сосредоточено на одномерных и двумерных конфигурациях

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

1. Расширение на более робастные модели столкновений: для моделирования турбулентных конфигураций
2. Более сложные уравнения состояния: такие как Пенга-Робинсона и другие
3. Крупномасштабные трёхмерные приложения: применение к практическим инженерным задачам

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

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

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

Недостатки

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

Влияние

1. Академический вклад: открывает новое направление для метода LBM в области сжимаемых неидеальных течений
2. Инженерная ценность: предоставляет важный численный инструмент для технологий сверхкритических жидкостей
3. Методологическое значение: концепция квазиравновесного аттрактора может повлиять на другие численные методы

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

1. Технологии сверхкритических жидкостей: CO₂-турбины, органические циклы Ренкина
2. Многофазные течения: разрушение капель, процессы фазовых переходов
3. Газодинамика ударных волн: сложные волновые взаимодействия в высокоскоростных течениях
4. Прямое численное моделирование: фундаментальные физические исследования, требующие высокой точности

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

1. Guardone, A. et al. (2024). Nonideal compressible fluid dynamics of dense vapors and supercritical fluids. Annual Review of Fluid Mechanics, 56, 241-269.
2. Hosseini, S.A. & Karlin, I.V. (2023). Lattice Boltzmann for non-ideal fluids: Fundamentals and practice. Physics Reports, 1030, 1-137.
3. Shan, X. & Chen, H. (1993). Lattice Boltzmann model for simulating flows with multiple phases and components. Physical Review E, 47(3), 1815.