Tworzenie siatki w BlueKenue

Wpis będzie jeszcze zmieniany. Dopiero się uczę Telemaca i BlueKenue, więc jest duża szansa, że są błędy w tym co napisałem niżej. Postaram się wrzucić później jakieś obrazki, poprawić tekst i lepiej to sformatować, żeby było czytelniejsze. Na końcu znajdziecie link do filmu na YT, na którym bazowałem oraz adres strony, którą też wykorzystywałem do nauki.

Przygotowanie siatki w BlueKenue do modelowania 2D i 3D w oprogramowaniu Telemac:

  1. Przygotowanie plików shapefile (poligon) z obszarem modelu. Ważne, żeby we wczytywanym pliku była tylko i wyłącznie JEDNA kolumna, która będzie wypełniona wartością liczbową, mówiącą o wielkości boku oczka siatki.
  2. Przygotowanie plików shapefile z obszarami o większej rozdzielczości siatki (większa rozdzielczość = gęstsza siatka). Może to być jeden plik, można też rozdzielić na kilka. Nakładanie się poligonów nie jest problemem, ponieważ zawsze będzie tworzona siatka o większej gęstości.
  3. Przygotowanie plików shapefile z obszarami wyłączonymi z modelowania, np. wyspami. BlueKenue ma problem z poprawną interpretacją poligonów, które mają dziury, więc nie można tego zrobić na wcześniejszym etapie.
  4. W Telemacu pliku wczytujemy przez File -> Import -> ArcView Shape File. Żeby zobaczyć co wczytaliśmy, musimy z panelu po lewej chwycić o przesunąć nasz plik do zakładki View ->2D view.
  5. Shp zapisujemy do plików i2s (2D Line Set). Zaznaczamy warstwę, którą chcemy wyeksportować, następnie File -> Save Copy As… i zapisujemy.
  6. Shp możemy usunąć, po czym wczytujemy właśnie wyeksportowane warstwy przez File -> Open, zmieniamy tyl pliku na Line Sets i dodajemy.
  7. Przygotowujemy generator siatki File -> New -> T3 Mesh Generator.
  8. W automatycznie otwartym nowym oknie z ustawieniami generatora wpisujemy nową wartość “Default Edge Length“. Co tam wpisać? Zależy od obszaru, ale modelując obszar wielkości Małego i Wielkiego Bełtu, śmiało możemy zacząć od 1000 (szczerze powiedziaszy to nawet nie wiem czy jest to potrzebne skoro i tak mamy te informacje z pliku shp).
  9. Obszar modelu chwytami i przesuwamy do “Outline“.
  10. Obszary o zmienionej rozdzielczości przesuwamy do “Density“.
  11. Wyspy przesuwamy do “SoftLines“.
  12. Wracamy do generatora i klikamy “Run“. Zapoznajemy się z podstawowymi informacjami i jeżeli nie ma niczego dziwnego, to potwierdzamy, że jest ok.
  13. Wygenerowaną siatkę przesuwamy do widoku jak wcześniej. Klikamy PPM -> Properties -> Data i tam widzimy ile węzłów zostało stworzonych.
  14. Ponownie zaznaczamy w panelu naszą siatkę i zapisujemy ją do formatu t3s File -> Save Copy As….
  15. Docelowym formatem jest Selafin. Konwertujemy naszą siatkę File -> New -> Selfafin Object.
  16. Najpierw tworzymy model dna, więc klikamy PPM na nowy obiekt i wybieramy Add Variable…, a New Variable Properties ustawiamy Name jako BOTTOM, Units powinno się automatycznie wypełnić jako M.
  17. Musimy jeszcze dodać brakujący plik z batymetrią File -> Open… i polecam plik xyz, bo z GeoTiff mam jakiś problem.
  18. Dane należy interpolować do siatki, więc File -> New -> 2D Interpolator.
  19. Od razu przesuwamy do nowego obiektu wcześniej dodany plik xyz.
  20. Z obiektu Selafin zaznaczamy BOTTOM i idziemy do Tools -> Map Object i wybieramy interpolator (domyślnie newInterpolator2D).
  21. Wynik można sprawdzić przesuwając BOTTOM do karty widoku.
  22. Z listy ponownie wybieramy obiekt Selafin, a następnie File -> Save Copy As….
  23. Teraz należy dodać lokalizację warunków brzegowych File -> New -> Boundary conditions (Conlim)…. Do listy naszych obiektów została dodana nowa rzecz i przesuwamy ją do widoku, wyłączając (lub usuwając z niego) wszystkie inne.
  24. Przesuwamy widok w miejsce, które nas interesuje i dwukrotnie klikamy LPM na węzeł, wciskamy i przytrzymujemy klawisz Shift i ponownie dwukrotnie klikamy LMP na kolejny węzeł, gdzie chcemy zakończyć.
  25. W nowym oknie wpisujemy nazwę oraz definiujemy typ warunku (w moim przypadku jest to Open Boudary with prescribed H o kodzie 5 4 4). To samo robimy w pozostałych miejscach, a na końcu zapisujemy do pliku bc2.
  26. Z listy obiektów klikamy na nasze warunki brzegowe (ikona tabeli) i zapisujemy do pliku cli.

 

Link do filmu, który pokazuje to krok po kroku: Generating mesh and geometry (selafin, *.slf) file for TELEMAC-2D using Bluekenue
Adres strony z materiałami do nauki: TELEMAC — Hydro-Informatics

De Bever Piotr

Autor bloga, z wykształcenia geodeta, hydrolog, modelarz hydrauliczny i GISowiec

Leave a Reply