Турбулентные потоки, сопровождаемые химическими реакциями, представляют собой одно из наиболее сложных направлений современной физики турбулентности. Здесь сталкиваются два нелинейных процесса — неустойчивое вихревое движение и нелинейная динамика химической кинетики. В результате возникают пространственно-временные структуры, которые невозможно описать, исходя лишь из законов гидродинамики или лишь из законов химии: требуется комплексный подход.
Математическое описание базируется на уравнениях Навье–Стокса для несжимаемой или сжимаемой среды, дополненных уравнениями переноса концентраций реагентов и теплового баланса:
$$ \rho \left( \frac{\partial u_i}{\partial t} + u_j \frac{\partial u_i}{\partial x_j} \right) = -\frac{\partial p}{\partial x_i} + \mu \nabla^2 u_i + f_i, $$
где ui — компоненты скорости, ρ — плотность, μ — динамическая вязкость, fi — массовые силы.
$$ \rho c_p \left( \frac{\partial T}{\partial t} + u_j \frac{\partial T}{\partial x_j} \right) = \lambda \nabla^2 T + Q_\text{реакц.}, $$
где T — температура, λ — теплопроводность, Qреакц. — тепловыделение от реакции.
$$ \frac{\partial Y_\alpha}{\partial t} + u_j \frac{\partial Y_\alpha}{\partial x_j} = D_\alpha \nabla^2 Y_\alpha + \dot{\omega}_\alpha, $$
где Yα — массовая доля компонента α, Dα — коэффициент диффузии, ω̇α — скорость химического превращения.
Совокупность этих уравнений образует систему, описывающую реактивные турбулентные течения.
1. Медленные реакции Если характерное время реакции τреакц. значительно превышает турбулентное время смешения τтурб., то турбулентность успевает полностью перемешать реагенты до их превращения. Такой режим называют смешение-контролируемым.
2. Быстрые реакции Когда τреакц. ≪ τтурб., реакция происходит мгновенно на границах смешивания потоков. Это приводит к формированию тонких реакционных зон, где скорость определяется интенсивностью турбулентной диффузии.
3. Сравнимые масштабы При τреакц. ∼ τтурб. возникает сложная взаимосвязь: реакция изменяет локальную температуру, а следовательно и плотность, что влияет на поле скоростей. В таких случаях необходимо учитывать двустороннюю связь между гидродинамикой и химией.
Для анализа реактивной турбулентности вводятся специальные безразмерные параметры:
$$ Da = \frac{\tau_\text{турб.}}{\tau_\text{реакц.}} $$
определяет относительную роль турбулентного перемешивания и химической кинетики.
$$ Ka = \left( \frac{\tau_\eta}{\tau_\text{реакц.}} \right), $$
где τη — время на колмогоровском масштабе. Оно характеризует способность мелкомасштабных вихрей разрушать реакционную зону.
$$ Le = \frac{\alpha}{D}, $$
где α — тепловая диффузивность, D — молекулярная диффузия реагентов. При Le ≠ 1 возникает различие в переносе тепла и вещества, что усиливает неустойчивости.
Особое значение имеют пламенные фронты, формируемые в условиях турбулентности. Их структура зависит от числа Рейнольдса и числа Дамкёлера.
1. Модели усреднения (RANS) Используют уравнения для средних величин, включая дополнительные модели турбулентного горения: eddy dissipation model, flamelet models, PDF-модели.
2. LES-моделирование Метод крупных вихрей позволяет описать крупномасштабную структуру турбулентности, добавляя подрешёточные модели для химических реакций. Особенно эффективен в условиях сложных геометрий и переменных плотностей.
3. Директное численное моделирование (DNS) Позволяет без модели решать всю систему уравнений на сетке, разрешающей все масштабы турбулентности и реакции. Однако чрезвычайно дорого с вычислительной точки зрения.
Турбулентность усиливает локальные неустойчивости пламени:
Эти процессы приводят к появлению сложных фрактальных структур пламени, которые играют ключевую роль в реальных реактивных потоках.
Исследование взаимодействия химических реакций и турбулентности имеет фундаментальное и прикладное значение: