CFD моделирование пристеночной области

Физическая структура пристеночной области

При течении вязкого теплоносителя вдоль твердой поверхности формируется турбулентный пограничный слой. Его внутренняя область, непосредственно примыкающая к стенке, характеризуется доминирующим влиянием молекулярной вязкости и резким изменением динамических и тепловых параметров в направлении нормали к стенке. Толщина этой внутренней области составляет примерно 10–20% от общей толщины пограничного слоя, но именно здесь сосредоточена основная доля теплового сопротивления (до 80–90%). Это объясняется тем, что в непосредственной близости от стенки турбулентный перенос импульса и тепла подавлен, и перенос осуществляется неэффективными молекулярными механизмами (теплопроводность и вязкость).

Принято описывать структуру этой области в безразмерных координатах, известных как «единицы стенки» (wall units). Ключевым масштабом является динамическая скорость (скорость трения), определяемая через касательное напряжение на стенке:

u_τ = √(τ_w / ρ)

Безразмерное расстояние от стенки и безразмерная продольная скорость определяются как:

y⁺ = (y · ρ · u_τ) / μ
U⁺ = U / u_τ

где y — нормальное расстояние от стенки, μ — динамическая вязкость, ρ — плотность. В этих координатах выделяют три характерных подслоя:

• Вязкий подслой (Viscous sublayer) (безразмерное расстояние от стенки менее 5): Турбулентное трение пренебрежимо мало, справедливо соотношение:

  μ_turbulent << μ

  Течение ламинарное. Профиль скорости линеен:

  U⁺ = y⁺

  Перенос тепла осуществляется только молекулярной теплопроводностью, что создает основное термическое сопротивление.

• Буферный слой (Buffer layer) (безразмерное расстояние от стенки от 5 до 30): Область сосуществования и сопоставимого влияния молекулярных и турбулентных механизмов переноса:

  μ_turbulent ≈ μ

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

• Логарифмический слой (Log-law region) (безразмерное расстояние от стенки более 30 при условии, что расстояние от стенки много меньше толщины пограничного слоя): Полностью турбулентный режим, где доминирует турбулентная вязкость:

  μ_turbulent >> μ

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

  U⁺ = (1/κ) · ln(y⁺) + B

  где κ ≈ 0,41 (константа Кармана), а B ≈ 5,0–5,5 (аддитивная постоянная для гладкой стенки). Аналогичный логарифмический профиль существует и для безразмерной температуры, что лежит в основе тепловых пристеночных функций.

Метод пристеночных функций (Wall Functions)

Метод пристеночных функций является стандартным инженерным подходом для моделирования пристеночной турбулентности, позволяющим избежать чрезмерного сгущения сетки к стенке. Его суть заключается в том, что вязкий подслой и буферная область не разрешаются численно. Вместо этого первый пристеночный узел сетки располагается в логарифмическом слое (безразмерное расстояние от стенки более 30), а профили скорости и температуры между этим узлом и стенкой восстанавливаются по полуэмпирическим формулам — «функциям стенки» (Launder & Spalding, 1974).

Ключевые положения метода:

• Предполагается локальное равновесие турбулентности, при котором генерация турбулентной кинетической энергии равна ее диссипации. Это позволяет связать касательное напряжение на стенке с кинетической энергией турбулентности в пристеночном узле.
• Граничные условия для моделей турбулентности задаются не на стенке, а из предположений о структуре течения в логарифмическом слое.
• Тепловой поток рассчитывается по тепловой пристеночной функции (аналог логарифмического закона для температуры) с учетом турбулентного числа Прандтля Pr_t ≈ 0,85–0,9.

Стандартный и масштабируемый подходы:

• Стандартные пристеночные функции: Наиболее экономичны, но сильно зависят от расположения первого узла. При безразмерном расстоянии от стенки менее 11 они могут давать значительные ошибки, так как узел попадает в буферный слой, где логарифмический закон не выполняется.
• Масштабируемые пристеночные функции (Scalable Wall Functions): Разработаны для преодоления проблемы неоптимальной сетки. Они ограничивают минимальное значение безразмерного расстояния от стенки снизу значением 11,06. Это делает решение менее чувствительным к локальным переизмельчениям сетки.
• Неравновесные пристеночные функции: Позволяют частично учесть влияние градиента давления на профиль скорости, что делает их предпочтительными для течений с отрывом, рециркуляцией или сильными перепадами давления.

Полное разрешение вязкого подслоя (Low-Reynolds Number Models)

Альтернативой пристеночным функциям является прямое численное моделирование всей пристеночной области, включая вязкий подслой. Этот подход требует использования «низкорейнольдсовых» версий моделей турбулентности и чрезвычайно мелкой сетки у стенки.

• Требования к моделям турбулентности:
  • Стандартная высокорейнольдсовая k-ε модель и большинство RSM моделей невалидны в вязком подслое. Для них обязательно использование пристеночных функций.
  • Модели, пригодные для расчета вплоть до стенки, содержат демпфирующие функции и дополнительные члены, учитывающие эффекты молекулярной вязкости. К ним относятся: низкорейнольдсовые k-ε модели (Launder-Sharma, Chien, Abe-Kondoh-Nagano), k-ω модель (в особенности SST-версия) и модель переноса касательных напряжений Transition SST.
• Двухслойный подход (Two-Layer Model / Enhanced Wall Treatment):
  Это гибридный метод, реализованный во многих CFD-кодах. Область течения разделяется по турбулентному числу Рейнольдса, определяемому через расстояние от стенки и кинетическую энергию турбулентности. В пристеночной зоне (при значениях этого числа менее 200) решается однопараметрическая модель, а во внешней — стандартная k-ε модель. Это позволяет использовать одну модель турбулентности на сетках с различным разрешением у стенки, плавно переключаясь от формализма пристеночных функций (при безразмерном расстоянии от стенки более 30) к разрешению вязкого подслоя (при безразмерном расстоянии порядка 1).

Практические рекомендации по выбору стратегии и построению сетки

Выбор между пристеночными функциями и полным разрешением — это всегда компромисс между точностью и вычислительными затратами, диктуемый спецификой задачи.

1. Рекомендации при использовании пристеночных функций:

• Целевой диапазон безразмерного расстояния от стенки: Центр первой пристеночной ячейки должен находиться строго в логарифмической области. Рекомендуемый диапазон от 30 до 150. Значения до 500 допустимы, но могут снизить точность расчета теплоотдачи. Важно избегать размещения первого узла в буферном слое (диапазон от 5 до 30).
• Минимальное количество ячеек: В логарифмическом слое должно быть не менее 8–10 ячеек для адекватного описания профиля.
• Шероховатость: При наличии заданной эквивалентной песочной шероховатости необходимо убедиться, что пристеночная функция включает модифицированный закон стенки, сдвигающий аддитивную постоянную B. Этот эффект важен для расчета гидросопротивления в промышленных теплообменниках.

2. Рекомендации при полном разрешении (Low-Re подход):

• Требование к безразмерному расстоянию от стенки: Расстояние до центра первой ячейки должно быть порядка 1. Считается хорошей практикой, когда безразмерное расстояние первых 3–5 ячеек не превышает 5.
• Разрешение вязкого подслоя: Внутри вязкого подслоя (безразмерное расстояние менее 5) должно находиться минимум 3–5 ячеек. Суммарно в пристеночной области рекомендуется иметь 15–20 ячеек для гладкого перехода к турбулентному ядру потока.
• Коэффициент роста: Рекомендуемый коэффициент роста (отношение высот соседних ячеек по нормали к стенке) не должен превышать 1,15–1,2. Резкие переходы размера ячеек в этой зоне приводят к ошибкам численного дифференцирования и снижению точности.

3. Общие замечания:

• Численная схема и безразмерное расстояние от стенки: Важно понимать, как конкретный решатель интерпретирует высоту первой ячейки. Для схем с интегрированием в центрах ячеек (cell-centered) безразмерное расстояние обычно вычисляется по расстоянию от стенки до центра первой ячейки (половина ее высоты).
• Верификация: При возможности следует проводить верификацию на задачах с известным аналитическим или экспериментальным эталоном (например, течение в трубе или на пластине), чтобы убедиться в корректности выбранной стратегии до начала расчета полного теплообменника.