Historia użytkownika
Poniższy przykład opisuje eksperymenty tunelu aerodynamicznego przeprowadzone przez Laboratorium Tunelu Wiatrowego Środowiska (EWTL) na Uniwersytecie w Hamburgu [1] jako przykład walidacyjny w części 9.4 WTG-Merkblatt M3. Zamierzamy użyć zmierzonych pól prędkości i danych o chropowatości modelu miasta Michel (Case BL3-3) do walidacji numerycznych symulacji CFD w złożonych strukturach miejskich. Przykład ten może należeć do Grupy 2, zgodnie z Rysunek 2.2 w WTG-Merkblatt-M3, opierając się na badaniu średniej wartości prędkości wiatru:
- G2: Wartości absolutne z wymaganiami średniej dokładności: Obszar zastosowania może obejmować parametry lub wstępne badania, gdy planowane są późniejsze badania z wyższą dokładnością (np. badania tunelowe klasy G3).
- R2: Samotność: wszystkie istotne kierunki wiatru z wystarczająco drobną rozdzielczością kierunkową.
- Z2: Statystyczne wartości średnie i odchylenia standardowe: pod warunkiem, że dotyczą stacjonarnych procesów przepływowych, dla których wystarczająca jest statystyczna weryfikacja fluktuacji za pomocą współczynnika szczytowego.
- S1: Efekty statyczne: Są wystarczające do przedstawienia modelu konstrukcyjnego z niezbędnymi szczegółami mechanicznymi, ale bez właściwości masowych i tłumiących.
Opis
Badanie koncentruje się na zidealizowanym, ale geometrycznie szczegółowym modelu miasta umieszczonym w przepływie warstwy granicznej atmosfery. Pomiary w tunelu aerodynamicznym zostały przeprowadzone w obiekcie WOTAN, który ma sekcję testową o długości 18 m, szerokości 4 m i wysokości 2,75-3,25 m. Odpowiednie pole chropowatości zostało scharakteryzowane długością chropowatości z0=1,53 m i wykładnikiem profilu α=0,27, reprezentującym warunki terenu "bardzo chropowatego". Zarejestrowano łącznie 1838 punktów pomiarowych dla kilku konfiguracji dachów. Czasowe komponenty poziome prędkości u i v, w tym wartości średnie, wariancje, korelacje i spektra, uzyskano za pomocą anemometru laserowego Dopplera (LDA) w 500-600 Hz. Punkty pomiarowe rozmieszczono w profilach pionowych i poziomych, w kanionach ulicznych i w określonych lokalizacjach powtarzalnych. Zbiór danych miasta Michel służy jako przypadek walidacyjny (C5) według wytycznych VDI 3783 część 9 [2]. Do walidacji, oprócz wskaźnika trafień, stosuje się odchylenie względne D=0,25 i odchylenie bezwzględne W=0,08, aby uwzględnić powtarzalność i niepewność pomiarową. Ten zbiór danych został zweryfikowany i przyjęty przez kilka instytucji (np. KalWin [3]) do celów walidacji CFD i porównania modeli.
Wymagania dotyczące dokładności WTG-Merkblatt M3
WTG-Merkblatt M3 dostarcza dwóch kluczowych metod weryfikacji wyników symulacji. Metoda Wskaźnika Trafień ocenia, ile z symulowanych wartości Pi poprawnie odpowiada wartościom referencyjnym Oi w obrębie zdefiniowanej tolerancji, przy użyciu podejścia binarnej klasyfikacji (trafienie albo pudło). To podejście ocenia wiarygodność symulacji przez obliczenie wskaźnika trafień q, podobnie jak funkcje ufności stosowane w teorii niezawodności. Z kolei metoda Znormalizowanego Średniego Błędu Kwadratowego (e2) oferuje bardziej szczegółową ocenę dokładności, kwantyfikując średnie kwadratowe odchylenie między wartościami symulowanymi a referencyjnymi, znormalizowane w celu uwzględnienia różnic skali. Razem te metody dostarczają jakościowych i ilościowych miar dla walidacji symulacji.
Wyniki i dyskusja
Porównanie między znormalizowanymi wartościami prędkości (U/Uref) uzyskanymi z symulacji RWIND a pomiarami eksperymentalnymi pokazuje umiarkowany poziom zgodności w analizowanym zbiorze danych. Analizowano łącznie 43 punkty walidacyjne, w których wartości odchyleń wynosiły od około 2% do niemal 50%, co wskazuje, że chociaż symulacja wychwytuje ogólną wielkość i trend pola prędkości, lokalne różnice pozostają znaczące w niektórych regionach. Analiza wskaźnika trafień dodatkowo podkreśla to zachowanie: tylko 18,60% punktów danych mieści się w ścisłej tolerancji ±10%, zwiększając się do 37,21%, gdy tolerancja jest rozluźniona do ±20%. Wartość znormalizowanego błędu średniego 𝑒2=0,2498 potwierdza umiarkowane ogólne odchylenie między przewidywanymi a zmierzonymi wartościami.
Przestrzennie, niższe odchylenia obserwuje się w obszarach prawdopodobnie charakteryzujących się bardziej stabilnymi lub przytwierdzonymi warunkami przepływu, podczas gdy większe różnice występują tam, gdzie oczekuje się silniejszych gradientów lub złożonych zjawisk przepływowych, takich jak strefy separacji czy regiony za falą. Różnice te mogą być przypisane ograniczeniom związanym z modelowaniem turbulencji RANS w stanie stacjonarnym, założeniom dotyczącym funkcji ściennej, efektom rozdzielczości siatki lub czułości na warunki wlotowe. Pomimo tych ograniczeń, symulacja wykazuje odpowiednią wydajność dla analizy inżynieryjnej i przewidywania trendów. Jednak w przypadku zastosowań wymagających większej miejscowej dokładności, może być konieczne dalsze udoskonalenie modelowania turbulencji, rozdzielczości siatki lub definicji warunków brzegowych.
Tabela 1: Porównanie znormalizowanej prędkości (U/Uref) między RWIND a danymi eksperymentalnymi
| X | Y | Z | U/Uref – RWIND | U/Uref – Eksperymentalne | Odchylenie (%) |
|---|---|---|---|---|---|
| -67.25 | -41.76 | 2 | 0.113 | 0.0942 | 16.64 |
| -70.41 | -27.84 | 2 | 0.120 | 0.0787 | 34.42 |
| -73.58 | -13.92 | 2 | 0.096 | 0.0693 | 27.81 |
| -76.75 | 0 | 2 | 0.098 | 0.0959 | 2.14 |
| -63.23 | 0 | 2 | 0.144 | 0.1305 | 9.38 |
| -61.71 | -12 | 2 | 0.134 | 0.1285 | 4.10 |
| -59.44 | -21.99 | 2 | 0.113 | 0.1399 | 23.81 |
| -57.17 | -31.98 | 2 | 0.124 | 0.1385 | 11.69 |
| -56.89 | -44.45 | 2 | 0.158 | 0.1186 | 24.94 |
| -46.54 | -47.15 | 2 | 0.149 | 0.1393 | 6.51 |
| -48.19 | -23.16 | 2 | 0.117 | 0.156 | 33.33 |
| -48.19 | -12 | 2 | 0.105 | 0.154 | 46.67 |
| -49.71 | 0 | 2 | 0.177 | 0.1347 | 23.90 |
| -36.19 | 0 | 2 | 0.117 | 0.1547 | 32.22 |
| -36.19 | -16.61 | 2 | 0.210 | 0.1633 | 22.24 |
| 31.81 | 127.75 | 2 | 0.091 | 0.086 | 5.49 |
| 58.81 | 119.75 | 2 | 0.096 | 0.093 | 3.13 |
| 72.31 | 127.75 | 2 | 0.065 | 0.0558 | 14.15 |
| 72.31 | 115.54 | 2 | 0.082 | 0.0423 | 48.41 |
| 85.81 | 118.91 | 2 | 0.149 | 0.0748 | 49.80 |
| 85.81 | 127.75 | 2 | 0.192 | 0.0978 | 49.06 |
| -149.89 | -124.58 | 2 | 0.103 | 0.1115 | 8.25 |
| -156.52 | -106.63 | 2 | 0.632 | 0.4036 | 36.15 |
| -142.43 | -99.38 | 2 | 0.658 | 0.334 | 49.24 |
| -141.58 | -112.44 | 2 | 0.476 | 0.4192 | 11.93 |
| -130.02 | -143.76 | 2 | 0.117 | 0.1723 | 47.26 |
| -119.53 | -120.2 | 2 | 0.313 | 0.384 | 22.68 |
| -127.31 | -105.09 | 2 | 0.634 | 0.3833 | 39.54 |
| -114.27 | -84.87 | 2 | 0.119 | 0.0894 | 24.87 |
| -105.26 | -112.85 | 2 | 0.466 | 0.3084 | 33.82 |
| -100.19 | -77.62 | 2 | 0.174 | 0.1187 | 31.78 |
| -36.27 | -94.25 | 2 | 0.483 | 0.2455 | 49.17 |
| -35.24 | -109.91 | 2 | 0.249 | 0.2782 | 11.73 |
| -48.97 | -103.35 | 2 | 0.366 | 0.2337 | 36.15 |
| -55.52 | -120.02 | 2 | 0.220 | 0.284 | 29.09 |
| -69.23 | -113.45 | 2 | 0.265 | 0.2553 | 3.66 |
| -103.49 | -58.25 | 2 | 0.270 | 0.1459 | 45.96 |
| -118.54 | 21.4 | 2 | 0.024 | 0.0224 | 6.67 |
| -121.46 | 20.73 | 2 | 0.062 | 0.0361 | 41.77 |
| -120.54 | 30.17 | 2 | 0.117 | 0.0952 | 18.63 |
| -117.61 | 30.84 | 2 | 0.139 | 0.0896 | 35.53 |
| -122.53 | 38.95 | 2 | 0.170 | 0.0907 | 46.65 |
| -125.45 | 38.28 | 2 | 0.158 | 0.1128 | 28.61 |
Tabela 2: Metryki walidacyjne porównania znormalizowanej prędkości (U/Uref)
| Metryka | Wartość |
|---|---|
| Łączna liczba punktów danych (N) | 43 |
| Wskaźnik trafień q (tolerancja 10%) | 18.60% |
| Wskaźnik trafień q (tolerancja 20%) | 37.21% |
| Forma odchylenia średniego błędu (e²) | 0.2498 |