Создание случайного пространственно скоррелированного профиля ветра для CFD-симуляции ветра

1. Введение

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

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

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

• Переходных CFD симуляций ветра

• Оценки давлений и сил на нагрузки

• Валидации по сравнению с условиями входа, подобными течению в аэродинамической трубе

• Продвинутые приложения за пределами консервативного охвата Еврокода

Этот подход значительно превосходит стандартную модель ветра Еврокода (EC), которая в первую очередь предназначена для статического проектирования нагрузки, а не для временно разрешенной физики потока.

2. Обзор реализованного метода

Генерация входного потока состоит из пяти сопряженных компонентов:

2.1 Пространственно дискретизированная входная плоскость

• Вход определяется на плоскости YZ при постоянном 𝑥=𝑥₀

• Регулярная сетка: 𝑛_𝑦 × 𝑛_𝑧 точек (например, 20 × 100 → 2,000 датчиков)

• Каждая точка сетки представляет временно зависимый сигнал скорости

• Среднее по ширине не применяется (важно для сохранения флуктуаций)

Это позволяет физически представить вертикальные и боковые эффекты когерентности.

2.2 Средний профиль ветра (закон степени, физически обеспеченный)

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

2.3 Модель турбулентности: спектр Каймала

Турбулентность в каждой точке входа генерируется с использованием спектра Каймала, который представляет поведение атмосферного пограничного слоя и определяет, как турбулентная энергия распределяется по частотам. Этот подход захватывает реалистичные частотно-зависимые характеристики флуктуаций ветра, а не предполагает постоянное поведение. Модель отражает интенсивность турбулентности и характерные масштабы длины для описания структуры изменений скорости. В результате сгенерированная турбулентность включает физически последовательное энергетическое содержание, которое меняется с частотой. В каждой точке входа турбулентность генерируется с использованием автоспектра Каймала, соответствующего теории атмосферного пограничного слоя [1]:

Это создает зависящее от частоты энергетическое содержание, которого нет в статическом подходе Еврокода.

2.4 Пространственная когерентность и кросс-спектральная матрица

Модель пространственной когерентности применяется для обеспечения физически реалистичной корреляции между точками входа, описывая, как подобие турбулентности уменьшается с расстоянием и частотой. Эта функция когерентности используется для определения взаимосвязи между сигналами в разных местах, отражая пространственно скоррелированные флуктуации ветра. На основе этой когерентности строится полная матрица плотности кросс-спектра, комбинирующая отдельные спектры с их пространственной корреляцией. Полученная матрица представляет как распределение энергии, так и пространственное соединение флуктуаций скорости по входу. Для обеспечения физически реалистичной корреляции между точками входа применяется модель когерентности, зависящая от расстояния:

Это приводит к полной матрице плотности кросс-спектра:

Ключевые характеристики:

• Полностью пространственно скоррелированная турбулентность

• Зависимое от частоты уменьшение корреляции

• Захватывает реалистичные структуры порывов по входу

2.5 Робустное разложение Холецкого (численная устойчивость)

Применяется робустное разложение Холецкого для поддержки численной устойчивости, когда матрица кросс-спектра становится плохо обусловленной, особенно около земли. Сначала матрица симметризуется для обеспечения согласованности, затем применяется адаптивная диагональная подгрузка для улучшения устойчивости. Используется автоматическая стратегия повторных попыток с увеличивающейся регуляризацией при необходимости, обеспечивая положительное определение. Этот подход обеспечивает стабильную и надежную генерацию скоррелированных временных рядов. Поскольку матрица кросс-спектра может становиться плохо обусловленной (особенно у земли), используется робустная факторизация Холецкого:

• Принудительная гермитовская симметрия:

• Адаптивная диагональная подгрузка

• Автоматическая стратегия повторных попыток с увеличивающейся регуляризацией

Это гарантирует положительное определение и стабильную генерацию временных рядов.

2.6 Построение временных рядов (обратное преобразование Фурье)

Случайные фазы вводятся для представления турбулентности в сигнале, в то время как спектр масштабируется и структурируется для обеспечения физически реалистичного результирующего действительного значения. Затем применяется обратное преобразование Фурье (FFT) для преобразования частотной информации в временно зависимый сигнал скорости. Окончательный результат - это поле скорости, составленное из среднего потока в сочетании с флуктуирующими компонентами.

• Случайные фазы применяются к каждой частоте

• Односторонний спектр масштабируется на sqrt(2Δ𝑓)

• Принудительная гермитовская симметрия

• Обратный FFT → временно разрешенный сигнал скорости в каждой точке входа

Результат:

3. Визуализация и диагностика

Метод позволяет:

• Мгновенные вертикальные профили 𝑈(𝑧,𝑡_o)

• Сравнение нескольких пробных линий (без артефактов усреднения)

• Полные анимации контуров входной плоскости 𝑈(𝑦,𝑧,𝑡)

Это делает физику потока непосредственно наблюдаемой, в отличие от статических профилей на основе Еврокода.

1 Распределение входных точек на плоскости YZ, показывающее пространственную сеть для генерации случайного турбулентного потока

2 Сравнение средней характеристики ветра и мгновенной флуктуирующей характеристик в нескольких местах по пролету

3 Поле вертикального ветра с моментальным впуском, показывающее пространственно коррелированное турбулентное распределение скорости ветра

4. Сравнение со стандартной моделью ветра Еврокода

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

Таблица 1: Сравнение со стандартной моделью ветра Еврокода

5. Ключевые различия в интерпретации

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

📘Еврокод

• Разработан для обеспечения безопасных проектных нагрузок

• Эффекты турбулентности встроены в частичные коэффициенты безопасности

• Нет информации о том, когда или как возникают нагрузки

🌬️Текущий метод

• Разрешает, когда, где и почему формируются пиковые давления
• Проникновение порывов
• Пики давления, вызванные разделением
• Коррелированная нагрузка на большие площади

⭐Особенно важно для:

• Легковесных конструкций

• Элементов облицовки и фасадов

• Мембран и оборудования на крыше

• Исследований динамического отклика

6. Инженерные последствия

Предлагаемый метод не стремится заменить Еврокод для регулирующего проектирования; вместо этого он должен пониматься как дополнительный подход, предоставляющий дополнительное физическое понимание, выходящее за рамки предположений, типично встроенных в процедуры Еврокода. Он особенно подходит для исследований чувствительности, исследовательских и валидационных целей, оптимизационных процессов и объяснения расхождений, которые могут возникнуть между результатами, основанными на Еврокоде и CFD-симуляциях. В практических инженерных приложениях этот подход часто помогает прояснить, почему значения давления, выведенные из CFD, отличаются от прогнозов Еврокода, особенно в отношении локальных пиковых давлений и переходных эффектов потока, причем это не обязательно указывает на ошибку в любой из методологий.

7. Резюме

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