Generowanie losowych przestrzennie skorelowanych profili wiatrowych do symulacji wiatru CFD

1. Wprowadzenie

Dokładna ocena działania wiatru na konstrukcje zależy nie tylko od zastosowanego modelu obciążenia, ale fundamentalnie od sposobu, w jaki reprezentowana jest nadchodząca warstwa graniczna atmosfery (ABL). W praktycznym projektowaniu inżynierskim wiatr jest zjawiskiem wysoce niestabilnym, turbulentnym i przestrzennie skorelowanym. Jednak większość podejść opartych na normach, w tym standardowy model wiatru Eurokodu, redukuje tę złożoność do równoważnych wartości statycznych, aby zapewnić prostotę, solidność i konserwatywne marginesy bezpieczeństwa.

W symulacjach wiatru opartych na obliczeniowej dynamice płynów (CFD) definicja warunku brzegowego na wlocie jest decydującym czynnikiem wpływającym na dokładność przewidywanych nacisków powierzchniowych, obciążeń szczytowych i efektów dynamicznych.

Metoda udokumentowana tutaj wdraża losowy, przestrzennie skorelowany generator turbulentnego wlotu na pionowej płaszczyźnie wlotowej (płaszczyzna YZ) i jest przeznaczona do:

• Przejściowych symulacji wiatru CFD

• Oceny ciśnienia i sił istotnych dla obciążenia

• Walidacji w porównaniu z warunkami wlotowymi podobnymi do tunelu aerodynamicznego

• Zaawansowanych zastosowań wykraczających poza konserwatywny zakres Eurokodu

To podejście znacznie wykracza poza standardowy model wiatru Eurokodu (EC), który jest przede wszystkim zaprojektowany do wyprowadzenia statycznego obciążenia projektowego, nie zaś do przepływów z rozdzielczością czasową.

2. Przegląd zastosowanej metody

Generacja wlotu składa się z pięciu sprzężonych komponentów:

2.1 Płaszczyzna wlotowa z dyskretyzacją przestrzenną

• Wlot zdefiniowany na płaszczyźnie YZ przy stałym 𝑥=𝑥₀

• Regularna siatka: 𝑛_𝑦 × 𝑛_𝑧 punktów (np. 20 × 100 → 2,000 punktów pomiarowych)

• Każdy punkt siatki reprezentuje sygnał prędkości zależny od czasu

• Nie stosuje się uśredniania bocznego (ważne dla zachowania fluktuacji)

Pozwala to na fizyczne przedstawienie efektów koherencji pionowej i bocznej.

2.2 Średni profil wiatru (prawo potęgowe, wymuszone fizycznie)

Średni profil wiatru definiowany jest za pomocą formuły prawa potęgowego, empirycznego przybliżenia inżynierskiego, które reprezentuje wzrost prędkości wiatru z wysokością w warstwie granicznej atmosfery. Model ten wiąże prędkość z warunkami wiatru odniesienia i eksponentem zależnym od terenu. Choć prostszy niż logarytmiczna formuła warstwy granicznej wywodząca się z teorii podobieństwa, podejście prawa potęgowego zapewnia praktyczną i szeroko stosowaną definicję wlotu dla zastosowań inżynierskich. Formuła zapewnia zgodny rozkład prędkości pionowej z odniesieniem do określonej wysokości.

2.3 Model turbulencji: Spektrum Kaimal

Turbulencja w każdym punkcie wlotowym jest generowana za pomocą spektrum Kaimal, które reprezentuje zachowanie warstwy granicznej atmosfery i definiuje, jak energia turbulentna jest dystrybuowana w różnych częstotliwościach. To podejście uwzględnia realistyczne charakterystyki częstotliwościowe fluktuacji wiatru, zamiast zakładać stałe zachowanie. Model odzwierciedla intensywność turbulencji i charakterystyczne skale długości, aby opisać strukturę wariacji prędkości. W rezultacie generowana turbulencja zawiera fizycznie spójny wkład energetyczny zmieniający się w zależności od częstotliwości. W każdym punkcie wlotowym turbulencja jest generowana za pomocą autospektrum Kaimal, zgodnie z teorią warstwy granicznej atmosfery [1]:

To prowadzi do energetycznych treści zależnych od częstotliwości, które są nieobecne w statycznym podejściu Eurokodu.

2.4 Koherencja przestrzenna i macierz przekroju widmowego

Model koherencji przestrzennej stosowany jest, aby zapewnić fizycznie realistyczną korelację między punktami wlotowymi poprzez opisanie, jak podobieństwo turbulencji maleje z odległością i częstotliwością. Funkcja koherencji jest używana do definiowania relacji między sygnałami w różnych lokalizacjach, odzwierciedlając przestrzennie skorelowane fluktuacje wiatru. Na podstawie tej koherencji konstruowana jest pełna macierz gęstości przekroju widmowego, łącząca poszczególne spektra z ich korrelacją przestrzenną. Uzyskana macierz reprezentuje zarówno rozkład energii, jak i sprzężenie przestrzenne fluktuacji prędkości na wlocie. Aby zapewnić fizycznie realistyczną korelację między punktami wlotowymi, stosowany jest model koherencji zależny od odległości:

To prowadzi do pełnej macierzy gęstości przekroju widmowego:

Kluczowe cechy:

• W pełni skorelowana przestrzennie turbulencja

• Zanik korelacji zależny od częstotliwości

• Uwzględnia realistyczne struktury podmuchów na wlocie

2.5 Solidna dekompozycja Cholesky'ego (stabilność numeryczna)

Zastosowana jest solidna dekompozycja Cholesky'ego, aby utrzymać stabilność numeryczną, gdy macierz przekroju widmowego staje się źle uwarunkowana, zwłaszcza blisko ziemi. Macierz jest najpierw symetryzowana, aby zapewnić zgodność, a następnie stosowane jest adaptacyjne ładowanie diagonali dla poprawy stabilności. Stosowana jest automatyczna strategia ponownej próby z rosnącą regularizacją, jeśli to konieczne, zapewniając dodatnią określoność. To podejście umożliwia stabilne i niezawodne generowanie skorelowanych danych szeregów czasowych. Ponieważ macierz przekroju widmowego może stać się źle uwarunkowana (szczególnie blisko gruntu), stosowana jest solidna faktoryzacja Cholesky'ego:

• Wymuszenie hermitowskie:

• Adaptacyjne ładowanie diagonali

• Automatyczna strategia ponownej próby z rosnącą regularizacją

To gwarantuje dodatnią określoność i stabilne generowanie szeregów czasowych.

2.6 Konstrukcja szeregów czasowych (odwrotna FFT)

Wprowadzane są losowe fazy, aby reprezentować turbulencję w sygnale, podczas gdy spektrum jest skalowane i strukturyzowane, aby zapewnić fizycznie realistyczny wynik rzeczywisty. Następnie stosowana jest odwrotna FFT, aby przekształcić informacje z domeny częstotliwości w sygnał prędkości zależny od czasu. Ostateczny wynik to pole prędkości złożone z przepływu średniego połączonego z fluktuującymi komponentami.

• Losowe fazy stosowane dla każdej częstotliwości

• Jednostronne spektrum skalowane przez sqrt(2Δ𝑓)

• Wymuszona symetria hermitowska

• Odwrotna FFT → sygnał prędkości z rozdzielczością czasową w każdym punkcie wlotowym

Rezultat:

3. Wizualizacja i diagnostyka

Metoda umożliwia:

• Chwilowe profile pionowe 𝑈(𝑧,𝑡_o)

• Porównania wielu linii pomiarowych (bez artefaktów uśredniania)

• Animacje konturowe pełnej płaszczyzny wlotowej 𝑈(𝑦,𝑧,𝑡)

To sprawia, że fizyka przepływu jest bezpośrednio obserwowalna, w przeciwieństwie do statycznych profili EC.

1 Rozkład punktów początkowych na płaszczyźnie YZ ukazujący siatkę przestrzenną używany do generowania losowych wzburzonych przepływów

2 Porównanie średniego profilu wiatru i chwilowych profili fluktuacji w różnych miejscach rozpiętości

3 Pionowe chwilowe pole wiatru przy wlocie pokazujące przestrzennie skorelowany rozkład prędkości turbulencji

4. Porównanie ze standardowym modelem wiatru Eurokodu

Tabela 1 podsumowuje fundamentalne różnice koncepcyjne i metodologiczne między standardowym modelem wiatru Eurokodu (EC) a optymalną metodą generacji wlotu CFD opartą na losowym modelu. Choć oba podejścia dążą do reprezentowania efektów wiatru na konstrukcje, bazują na różnych filozofiach modelowania i służą różnym celom inżynierskim.

Tabela 1: Porównanie ze standardowym modelem wiatru Eurokodu

5. Kluczowe różnice w interpretacji

Podejście Eurokodu jest głównie przeznaczone do zapewnienia bezpiecznych obciążeń projektowych poprzez uwzględnienie efektów turbulencji poprzez częściowe współczynniki bezpieczeństwa, ale nie dostarcza szczegółowych informacji na temat czasu lub mechanizmów rozwoju obciążenia. W przeciwieństwie do tego, obecna metoda umożliwia głębsze zrozumienie fizyki poprzez rozwiązywanie, kiedy, gdzie i dlaczego występują szczyty ciśnienia, w tym penetrację podmuchów, szczyty ciśnienia indukowane separacją oraz przestrzennie skorelowane efekty obciążenia. To czyni ją szczególnie wartościową dla lekkich konstrukcji, elementów elewacyjnych i okładzinowych, systemów membranowych, wyposażenia dachowego i badań obejmujących dynamiczną reakcję konstrukcyjną:

📘Eurokod

• Zaprojektowany w celu zapewnienia bezpiecznych obciążeń projektowych

• Efekty turbulencji uwzględnione w częściowych współczynnikach bezpieczeństwa

• Brak informacji o tym, kiedy i jak występują obciążenia

🌬️Obecna metoda

• Rozwiązuje, kiedy, gdzie i dlaczego formują się szczyty ciśnienia
• Penetracja podmuchów
• Szczyty ciśnienia indukowane separacją
• Skorelowane obciążenie na dużych obszarach

⭐Szczególnie ważne dla:

• Lekkich konstrukcji

• Elementów elewacyjnych i okładzinowych

• Systemów membranowych i wyposażenia dachowego

• Badań reakcji dynamicznej

6. Implikacje inżynierskie

Proponowana metoda nie ma na celu zastąpienia Eurokodu w projektowaniu zgodnym z przepisami; powinna być rozumiana jako podejście uzupełniające, które zapewnia dodatkowe spojrzenie fizyczne wykraczające poza założenia zwykle zawarte w procedurach Eurokodu. Jest szczególnie odpowiednia do badań wrażliwości, celów badawczych i walidacyjnych, procesów optymalizacyjnych oraz wyjaśniania rozbieżności, które mogą wystąpić między wynikami opartymi na Eurokodzie a symulacjami CFD. W praktycznych zastosowaniach inżynierskich podejście to często pomaga wyjaśnić, dlaczego wartości ciśnienia wyprowadzone z CFD różnią się od przewidywań Eurokodu, zwłaszcza w odniesieniu do lokalnych szczytów ciśnienia i przejściowych efektów przepływu, bez konieczności wskazywania błędu w którejkolwiek z metodologii.

7. Podsumowanie

Przedstawiona metoda generacji wlotu wprowadza fizycznie spójne modelowanie turbulencji połączone z koherencją przestrzenną i reprezentacją pola wiatru z rozdzielczością czasową, wspierane przez solidną implementację numeryczną. W porównaniu do standardowego podejścia Eurokodu, metoda ta przesuwa perspektywę analityczną z konserwatywnych kopert obciążeniowych w kierunku zrozumienia zachowania przepływu bazującego na fizyce, umożliwiając tym samym głębszy wgląd i bardziej zaawansowane analizy inżynierii wiatrowej.