Корреляционные функции и спектральный анализ

Понятие корреляции в турбулентных потоках

В турбулентной гидродинамике ключевую роль играет статистическое описание поля скоростей. Поле u(x, t) характеризуется сложной пространственно-временной структурой, и для его количественного анализа используются корреляционные функции.

Двухточечная корреляционная функция скорости определяется как среднее произведение отклонений скорости в двух точках пространства:

Rij(r) = ⟨ui(x) uj(x + r)⟩,

где r — вектор разделения, i, j — компоненты вектора скорости, а угловые скобки ⟨⋅⟩ обозначают статистическое усреднение по ансамблю или времени.

Основные свойства корреляционных функций:

  1. Симметрия: Rij(−r) = Rji(r).
  2. Нормировка: при r = 0 корреляция достигает максимума Rij(0) = ⟨uiuj, что связано с локальной кинетической энергией турбулентности.
  3. Зависимость от масштаба: корреляции быстро убывают с увеличением |r|, что отражает диссипативные процессы на малых масштабах.

Для изотропной турбулентности корреляционная функция упрощается и может быть разложена на продольную и поперечную составляющие:

$$ R_{ij}(\mathbf{r}) = \left[ R_L(r) - R_N(r) \right] \frac{r_i r_j}{r^2} + R_N(r) \delta_{ij}, $$

где RL(r) и RN(r) — продольная и поперечная корреляции, соответственно.


Функция структуры

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

S2(r) = ⟨|u(x + r) − u(x)|2⟩.

Связь с корреляционной функцией выражается через:

S2(r) = 2[⟨u2⟩ − R(r)].

Функция структуры позволяет выделять характерные масштабы турбулентности, в том числе инерциальный диапазон, где выполняется известная зависимость Колмогорова:

S2(r) ∼ (εr)2/3,

где ε — скорость диссипации турбулентной энергии на единицу массы.


Спектральное представление турбулентности

Для более детального анализа используется энергетический спектр E(k), который описывает распределение кинетической энергии по волновым числам k = 2π/λ:

$$ \frac{1}{2} \langle u_i u_i \rangle = \int_0^\infty E(k) \, dk. $$

Связь между корреляционной функцией и спектром задается преобразованием Фурье:

$$ E(k) = \frac{1}{2\pi^2} \int_0^\infty R(r) \frac{\sin(kr)}{kr} \, r^2 dr, $$

что позволяет переходить от пространственного описания к волновому.

Основные особенности спектра:

  1. Интервал больших масштабов: энергия концентрируется на малых k, где формируется крупномасштабная структура потока.
  2. Инерциальный диапазон: диапазон k, в котором энергия переносится от больших к малым масштабам без диссипации. Согласно гипотезе Колмогорова, для изотропной турбулентности:

E(k) = CKε2/3k−5/3,

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


Временные корреляции и спектр

Помимо пространственной структуры, важны временные корреляции:

Rij(τ) = ⟨ui(t)uj(t + τ)⟩,

где τ — временной лаг. Преобразование Фурье по времени дает спектр частот E(ω), который позволяет анализировать динамику турбулентного потока.

Для стационарной турбулентности часто используется гипотеза Тейлора, связывающая временные и пространственные корреляции через “замороженный” поток:

ω = kU,

где U — характерная скорость потока.


Кросс-корреляции и перенос энергии

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

Ruv(r) = ⟨u(x)v(x + r)⟩.

На основе спектрального анализа кросс-корреляций определяется коэффициент переноса энергии между компонентами и выявляются анизотропные эффекты.


Практическое применение корреляций и спектров

  1. Определение характерных масштабов турбулентности: по корреляционной длине L = ∫0R(r)/R(0) dr.
  2. Расчет коэффициентов диффузии для моделирования турбулентного перемешивания.
  3. Сравнение с экспериментальными данными: спектры энергии, измеренные методом ЛДА или PIV, позволяют оценивать достоверность моделей турбулентности.
  4. Анализ анизотропии и структурных особенностей: использование продольных и поперечных корреляций.