Historia użytkownika
Poniższy przykład opisuje eksperymenty przeprowadzone w tunelu aerodynamicznym przez Laboratorium Tunelu Wiatrowego (EWTL) na Uniwersytecie w Hamburgu [1] jako przypadek walidacyjny w Części 9.4 WTG-Merkblatt M3. Zamierzamy wykorzystać zmierzone pola prędkości modelu miasta Michel (Przypadek BL3-3) do walidacji numerycznych symulacji CFD w złożonych strukturach miejskich. Przykład może należeć do Grupy 1, zgodnie z Rysunek 2.2 w WTG-Merkblatt-M3, na podstawie badania ś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 o wyższej dokładności (np. badanie w tunelu aerodynamicznym klasy G3).
- R2: Samotne: wszystkie istotne kierunki wiatru z wystarczająco drobną rozdzielczością kierunkową.
- Z2: Średnie statystyczne i odchylenia standardowe: pod warunkiem, że dotyczą one stacjonarnych procesów przepływu, dla których wystarczająca jest statystyczna weryfikacja fluktuacji za pomocą współczynnika szczytowego.
- S1: Efekty statyczne: Są one wystarczające do reprezentacji modelu strukturalnego z niezbędnym szczegółem mechanicznym, ale bez właściwości masy i tłumienia.
Opis
Badanie koncentruje się na zidealizowanym, ale geograficznie szczegółowym modelu miasta umieszczonym w przepływie w warstwie granicznej atmosfery. Pomiary w tunelu wiatrowym zostały przeprowadzone w urządzeniu WOTAN, o sekcji testowej długości 18 m, szerokości 4 m i wysokości 2,75-3,25 m. Odpowiednie pole szorstkości charakteryzowało się długością szorstkości z0=1,53 m i wykładnikiem profilu α=0,27 reprezentującymi "bardzo szorstkie" warunki terenowe. Łącznie zarejestrowano 1,838 punktów pomiarowych dla kilku konfiguracji dachowych. Zmienno-czasowe składowe poziome prędkości u i v, w tym wartości średnie, wariancje, korelacje i widma, uzyskano za pomocą 2D Anemometru Laserowo-Dopplerowskiego (LDA) przy częstotliwości 500-600 Hz. Punkty pomiarowe były rozmieszczone w profilach pionowych i poziomych, w kanionach ulicznych i w określonych miejscach powtarzalności. Zestaw danych miasta Michel służy jako przypadek referencyjny do walidacji (C5) zgodnie z wytycznymi VDI 3783 Część 9 [2]. Do walidacji, oprócz współczynnika trafienia, stosuje się relatywne odchylenie D=0,25 i absolutne odchylenie W=0,08, aby uwzględnić powtarzalność i niepewność pomiarów. Ten zestaw danych został zweryfikowany i przyjęty przez kilka instytucji (np. KalWin [3]) do celów walidacji CFD i porównania modeli.
Wymagania dokładności według WTG-Merkblatt M3
WTG-Merkblatt M3 zapewnia dwie kluczowe metody walidacji wyników symulacji. Metoda Wskaźnika trafienia ocenia, ile z wartości symulowanych Pi poprawnie odpowiada wartościom referencyjnym Oi w określonej tolerancji, wykorzystując podejście klasyfikacji binarnej (trafienie lub nietrafienie). To podejście ocenia niezawodność symulacji, obliczając wskaźnik trafienia q, podobny do funkcji zaufania używanych w teorii niezawodności. W przeciwieństwie do tego, 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 i referencyjnymi, znormalizowane w celu uwzględnienia różnic skali. Razem, te metody dostarczają jakościowych i ilościowych miar do walidacji symulacji.
Wyniki i dyskusja
Porównanie znormalizowanych wartości prędkości (U/Uref) uzyskanych z symulacji RWIND i pomiarów eksperymentalnych wykazuje umiarkowaną zgodność w badanym zbiorze danych. Przeanalizowano łącznie 43 punkty walidacyjne, przy których wartości odchylenia wahają się od około 2% do prawie 50%, co wskazuje, że chociaż symulacja uchwyciła ogólną wielkość i trend pola prędkości, lokalne rozbieżności pozostają znaczne w niektórych obszarach. Analiza wskaźnika trafienia dodatkowo podkreśla to zachowanie: tylko 18,60% punktów danych mieści się w ścisłej tolerancji ±10%, rosnąc do 37,21% przy bardziej liberalnej tolerancji ±20%. Znormalizowana wartość średniego błędu 𝑒2=0,2498 potwierdza umiarkowane ogólne odchylenie między wartościami przewidywanymi i mierzonymi.
Przestrzennie, niższe odchylenia obserwuje się w regionach prawdopodobnie charakteryzujących się bardziej stabilnymi lub przyczepnymi warunkami przepływu, podczas gdy większe rozbieżności występują tam, gdzie można oczekiwać silniejszych gradientów lub złożonych zjawisk przepływu, takich jak strefy separacji lub regiony turbulencji. Te różnice mogą być przypisane ograniczeniom inherentnym w modelowaniu turbulencji RANS, założeniom funkcji ściennych, efektom rozdzielczości siatki lub wrażliwości na warunki wlotowe. Pomimo tych ograniczeń, symulacja wykazuje odpowiednią wydajność do analiz inżynierskich i przewidywania trendów. Jednak w przypadku zastosowań wymagających większej lokalnej dokładności może być konieczna dalsza poprawa 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 – Experimental | 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 walidacji dla porównania znormalizowanej prędkości (U/Uref)
| Metryka | Wartość |
|---|---|
| Całkowita liczba punktów danych (N) | 43 |
| Wskaźnik trafienia q (tolerancja 10%) | 18.60% |
| Wskaźnik trafienia q (tolerancja 20%) | 37.21% |
| Odchylenie formy średniego błędu (e²) | 0.2498 |