Оценка статистической согласованности синтетической турбулентности в условиях ограничений волновых чисел

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

• ID статьи: 2510.13089
• Название: Evaluation of Statistical Consistency in Synthetic Turbulence under Wavenumber Bounds
• Авторы: Hongyuan Lin, Yi Liu, Shizhao Wang, Chun-Hian Lee
• Классификация: physics.flu-dyn
• Дата публикации: 15 октября 2025 г.
• Ссылка на статью: https://arxiv.org/abs/2510.13089

Аннотация

Метод случайных рядов Фурье (RFM) широко используется для синтеза турбулентности благодаря его математической ясности и простоте. Однако между предустановленными входными параметрами и синтезированными результатами сохраняются расхождения, источники которых еще не получили должной оценки. Данное исследование направлено на систематическую оценку влияния калибровки спектральных коэффициентов, сеточных ограничений и границ волновых чисел на точность турбулентности, генерируемой методом RFM. Результаты показывают, что переколибровка спектральных коэффициентов имеет решающее значение для обеспечения согласованности кинетической энергии турбулентности. Верхняя граница волнового числа, определяемая разрешением сетки, контролирует общий уровень кинетической энергии турбулентности, тогда как нижняя граница, установленная размером расчетной области, контролирует точность энергетического спектра в диапазоне низких волновых чисел. Кроме того, расширение диапазона волновых чисел за пределы сеточных ограничений может повысить точность кинетической энергии турбулентности, но одновременно усиливает спектральные искажения.

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

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

Синтетическая турбулентность имеет решающее значение в вычислительной гидродинамике (CFD) и вычислительной аэроакустике (CAA), с приложениями, включающими исследование широкополосного шума, генерацию турбулентного входного потока для высокоточного моделирования и исследование дисперсии частиц. Недостаточная точность синтетической турбулентности может привести к введению искусственного шума или искажению изучаемых физических механизмов.

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

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

1. Ошибки на уровне энергии: проявляются как отклонения полной кинетической энергии турбулентности (TKE)
2. Ошибки на спектральном уровне: проявляются как искажения энергетического спектра турбулентности E(κ) в дискретизированном диапазоне волновых чисел

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

Хотя предыдущие исследования уже рассматривали определение коэффициентов или предлагали эмпирические критерии для волновых чисел, мало внимания уделялось источникам ошибок, что может привести к неожиданным отклонениям в практических приложениях. Данное исследование систематически оценивает влияние калибровки спектральных коэффициентов, сеточных ограничений и границ волновых чисел на статистическую согласованность турбулентности, генерируемой методом RFM.

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

1. Систематический анализ ошибок: впервые систематически оценено влияние калибровки спектральных коэффициентов, сеточных ограничений и границ волновых чисел на точность синтетической турбулентности
2. Идентификация источников ошибок: четко определено, что ошибки TKE в основном происходят из усечения высоких волновых чисел, а ошибки E(κ) — из недостаточного покрытия низких волновых чисел
3. Совершенствование теоретической базы: установлена теоретическая база, согласно которой выбор границ волновых чисел должен учитывать сеточные ограничения
4. Практические руководящие принципы: предоставлены практические рекомендации для выбора параметров в приложениях CFD и CAA

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

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

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

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

Основы метода случайных рядов Фурье

Компонента пульсирующей скорости u' в произвольном пространственном положении x может быть представлена как:

u'(x) = 2∑(m=1 to M) u_m cos(k_m·x + ψ_m)σ_m

где:

• M: общее количество мод
• u_m: амплитуда m-й моды
• k_m: волновой вектор m-й моды
• ψ_m: случайный угол фазы
• σ_m: единичный вектор направления

Спектр фон Кармана-Пао

Используется широко применяемая форма спектра фон Кармана-Пао:

E(κ) = αu'²/κ_e · (κ/κ_e)⁴/[1+(κ/κ_e)²]^(17/6) · exp[-2(κ/κ_η)²]

где α — спектральный коэффициент, κ_e — волновое число, содержащее энергию, κ_η — волновое число Колмогорова.

Калибровка спектральных коэффициентов

Для учета эффектов конечного числа Рейнольдса вводится локальное число Рейнольдса:

Re_L = TKE²/(εν)

Откалиброванные спектральные коэффициенты α и волновое число, содержащее энергию κ_e, определяются путем решения уравнения, содержащего конфлюэнтную гипергеометрическую функцию.

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

1. Систематическая аналитическая база: установлена комплексная аналитическая база, учитывающая как калибровку спектральных коэффициентов, так и границы волновых чисел
2. Теория сеточных ограничений: четко определены границы волновых чисел, обусловленные сеточными ограничениями: κ_min = 2π/L, κ_max = π/Δx
3. Механизм разделения ошибок: разделены ошибки TKE и спектральные ошибки, проведен раздельный анализ их факторов влияния

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

Базовые параметры

• Пульсация скорости: u' = 0,25 м/с
• Скорость диссипации: ε = 0,5402 м²/с³
• Предустановленная кинетическая энергия турбулентности: TKE_pre = 0,09375 м²/с²

Конфигурация сетки

• Базовая сетка: 128³ точек
• Длина расчетной области: 0,18π м
• Количество мод: 10³ (линейное распределение)

Показатели оценки

Относительная ошибка TKE

Err_TKE = |TKE_syn - TKE_pre|/TKE_pre

Относительная спектральная ошибка

Err_E(κ_j) = |E_syn(κ_j) - E_pre(κ_j)|/E_pre(κ_j)

Результаты экспериментов

Эффективность калибровки спектральных коэффициентов

Откалиброванные спектральные коэффициенты снизили ошибку TKE с 17,28% до 10,14%, подтверждая важность калибровки коэффициентов.

Эффект расширения границ волновых чисел

При расширении верхней границы до 10κ_max:

• Ошибка TKE значительно снижается до 1,239%
• Но спектральная ошибка резко возрастает до 28,50% (среднее значение) и 41,15% (RMS)
• Это указывает на то, что выход за пределы сеточных ограничений приводит к эффектам наложения спектра

Влияние верхней границы волнового числа

С увеличением разрешения сетки с 32³ до 512³:

• Err_TKE снижается с 32,8% до 0,64%
• Спектральная ошибка остается в диапазоне 5-8%
• Подтверждается, что верхняя граница в основном контролирует согласованность TKE

Влияние нижней границы волнового числа

Эффект увеличения расчетной области (снижение κ_min):

• Err_TKE снижается с 36,12% до примерно 8%
• Спектральная ошибка снижается с 17% до менее 5%
• Подтверждается, что нижняя граница в основном контролирует точность спектра низких волновых чисел

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

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

• Kraichnan (1970): первоначально предложил метод случайных рядов Фурье
• Karweit et al. (1991): выбор диапазона волновых чисел на основе экспериментальных данных
• Béchara et al. (1994): введение интегрального масштаба и волнового числа Колмогорова
• Bailly & Juvé (1999): связь между расчетной областью и разрешением сетки
• Saad et al. (2017): связь диапазона волновых чисел с теорией дискретизации

Сравнение вклада данной работы

По сравнению с существующими исследованиями, данная работа впервые систематически анализирует источники ошибок, а не просто предоставляет эмпирические критерии, обеспечивая теоретическую основу для выбора параметров.

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

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

1. Калибровка спектральных коэффициентов является предпосылкой согласованности TKE: при неограниченном диапазоне волновых чисел переколиброванные коэффициенты могут устранить отклонения между расчетными и предустановленными значениями TKE
2. Верхняя граница волнового числа в основном контролирует согласованность TKE: повышение разрешения сетки, тем самым повышая теоретическую верхнюю границу, может значительно снизить ошибку TKE
3. Нижняя граница волнового числа контролирует точность спектра низких волновых чисел: расширение расчетной области, тем самым снижая нижнюю границу, может улучшить согласованность E(κ) в диапазоне низких волновых чисел
4. Выбор диапазона волновых чисел должен учитывать сеточные ограничения: когда границы выходят за пределы сеточных ограничений, происходит эффект наложения спектра, приводящий к значительному увеличению ошибки E(κ)

Ограничения

1. Единственная спектральная форма: исследование в основном основано на спектре фон Кармана-Пао; применимость других спектральных форм требует дальнейшей проверки
2. Предположение об изотропности: исследование предполагает изотропную турбулентность; случаи анизотропной турбулентности в практических приложениях не рассматриваются
3. Фиксированное распределение мод: используется линейное распределение мод; влияние других распределений требует дальнейшего изучения

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

1. Оптимизация количества и распределения мод: изучение влияния количества и распределения мод в пределах фиксированных границ волновых чисел
2. Расширение на анизотропную турбулентность: распространение метода на сценарии анизотропного турбулентного потока
3. Верификация в практических приложениях: проверка руководящих принципов выбора параметров в конкретных приложениях CFD и CAA

Углубленная оценка

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

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

Недостатки

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

Влияние

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

Применимые сценарии

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

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

Данная работа цитирует 21 важную ссылку, охватывающую историческое развитие и последние достижения в методах синтеза турбулентности, обеспечивая прочную теоретическую основу для исследования. Ключевые ссылки включают новаторскую работу Kraichnan, анализ теории дискретизации Saad и др., а также последние методы калибровки коэффициентов Guglielmi и др.