Der Aufruf der stationären Strömung kann auch direkt ohne die Eingabemaske über die Kommandozeile erfolgen, sofern eine Batch-Datei *.bsi für das Modul SITRA oder eine Batch-Datei *.bge für das Modul GEONEU im Verzeichnis vorhanden ist. Durch eine Eingabe in der Kommandozeile wird in das Verzeichnis gewechselt, in dem die Berechnung ausgeführt werden soll. Die Eingabe von sitra und dem Bestätigen mit der Enter-Taste führt dazu, dass eine Berechnung mit dem Modul SITRA gestartet wird. Der voreingestellte Batch-Datei-Name für das Modul SITRA ist sitra.bsi. Diese Datei wird im Verzeichnis gesucht. Mit dem Aufruf sitra Dateiname kann aber auch eine andere Batch-Datei angegeben werden (die Erweiterung ".bsi" wird bei Bedarf automatisch angehängt).
Die Eingabe von geoneu und dem Bestätigen mit der Enter-Taste führt dazu, dass eine Berechnung mit dem Modul GEONEU gestartet wird. Der voreingestellte Batch-Datei-Name für das Modul GEONEU ist geoneu.bge. Diese Datei wird im Verzeichnis gesucht. Mit dem Aufruf geoneu Dateiname kann aber auch eine andere Batch-Datei angegeben werden (die Erweiterung ".bge" wird bei Bedarf automatisch angehängt).
Die Batch-Dateien können auch von Hand mit einem beliebigen Editor erstellt oder verändert werden.
Beispiel einer Batch-Datei für eine stationäre Strömungsberechnung mit dem Modul GEONEU
out.g # Ausgabedatei
1 1.0e-8 1.0e-6 # Gleichungsloeser(1-iterativ,0-direkt),Residuum,Iter.Differenz
0 0 # Kontrollinien? Geschw. (0-mittel,1-rechts,-1,links)
1 # Berechnung der Geschwindigkeiten? (1-ja,0-nein)
0 # K-Werte aus Eichen uebernehmen? (1-ja,0-nein)
1 5 0 0 1.0 # Iter.Maec.? (1/-1/0)(ja/fort/nein),Anz.,Pit,Mit,Daempf
1 # Bahnlinienberechnung? (1-ja,0-nein)
0.0 # max Wegstrecke (0.0 -> keine Beschr.)
0.0 # max Fliesszeit (0.0 -> keine Beschr.)
0 # Bahnlinie gegen(0) oder in(1) Stroemungsrichtung
0.2 # Porositaet
0.5 # Fangradius
27 0 0 # Anz. Startpunkte in Elementen, um Knoten, einz. Punkten
1350 1668 1861 1960 1885 1664 1437 1346 1345 1547 # Elementnummern
1854 1967 1843 218 1833 1564 1285 1148 1033 1029 # Elementnummern
1135 1200 1030 1143 1445 1674 208 # Elementnummern
Beispiel einer Batch-Datei für eine stationäre Strömungsberechnung mit dem Modul SITRA
out.s # Ausgabedatei
1 # ausfuehrliches Protokoll (1-Ja/0-Nein)
1 0 # Gleichungsloeser P/U (1-iterativ, 0-direkt, -1=Mehrgitter)
1 0 # Kontrolllinienber. (1-Ja, 0-Nein), IMITTEL
0 0 # Pecletzahlen, K-Werte aus Eichen (1-Ja, 0-Nein)
STROEMUNG #
1 # ges./unges. Rechnung (1-ja,0-nein)
0 0 0 # gleiche Konz., dichteabh. Rech., p oder h (1-Ja/0-Nein)
5 0.5 # Anz. Iterationsschritte, Daempfungsfak.
0 0 # Stroemung/Transport stat./0 inst./1 inst.+gesp./unges./2
0 0 # Start. pote/Konz aus ASAT(3) null(2) EICH/AKON(1) =0 (0)
CONDITION NUMBER
In einer Zeile werden alle Eingaben nach dem Zeichen '#' ignoriert, d.h., damit können Kommentare in die Batch-Dateien geschrieben werden.
Die abgebildeten Batch-Dateien für die stationäre Strömungsberechnung stehen auf unserer Homepage zum Download bereit.
Hier findet sich eine vollständige Beschreibung der Batch-Datei für eine Strömungsberechnung mit GEONEU
1. Zeile: Name der Ausgabedatei (beginnend ab 1. Zeichen!)
2. Zeile: Flag fuer Gleichungsloeser: 0=direkt, 1= iterativ,
3. Zeile: 3 Flags: 1. Flag für Kontrolllinienrechnung (1=ja/0=nein), 2. Flag für Geschwindigkeiten für Kontrolllinienrechnung: 0= Mitteln der Geschwindigkeiten, 1= Geschwindigkeiten rechts der Kontrolllinie, -1= Geschwindigkeiten links der Kontrolllinie, 3. Flag für Berechnung vertikaler Durchflussmengen (nur bei 3d-Modell ohne wegfallende Elemente!)
4. Zeile: Flag für Berechnung der Geschwindigkeiten (1=ja/0=nein) (automatisch bei Kontrolllinien oder Durchflussmengenberechnung und automatisch bei Bahnlinienberechnung (s.u.))
5. Zeile: Flag für Übernahme von K-Werten aus Kalibrierung (1=ja/0=nein) (geht nur, wenn vorher Kalibrierung gelaufen ist!)
6. Zeile: Flag und Parameter für Iteration der Mächtigkeit / freien Oberfläche:
Bei 3D-Modellen oder Vertikalmodellen:
Flag für Iteration der freien Oberfläche (1 oder -1 =ja/0=nein), Anzahl Iterationen (bei Flag =1 oder -1), Flag für Startwerte aus EICH (bei Flag=1), Dämpfungsfaktor für Iteration
Bei Horizontalmodellen mit UNTE:
Flag für Iteration der Mächtigkeit (1 oder -1 =ja/0=nein), Anzahl Iterationen, Flag für Protokoll der iterierten Mächtigkeiten (1=ja/0=nein), Flag für Protokoll der iterierten Potentiale (1=ja/0=nein), Dämpfungsfaktor für Iteration
Bei Flag=-1 wird auch iteriert. Ist dann vorher schon eine Iteration gelaufen, wird diese fortgesetzt! (Übernahme von K-Werten aus Kalibrierung ist dann allerdings nicht mehr möglich!)
7. Zeile: Flag für Bahnlinienberechnung (1=ja/0=nein)
ab 8. Zeile: Eingabeparameter für Bahnlinienberechnung
8. Zeile: Max. Wegstrecke einer Bahnlinie (in m)
9. Zeile: Max. Fliesszeit einer Bahnlinie (in ZE)
10. Zeile: Bahnlinien gegen (0) oder in (1) die Fliessrichtung
11. Zeile: Porosität
12. Zeile: Fang'radius' an Brunnen (zw. 0.0 und 1.0)
13. Zeile: 3 Flags für Anzahl Bahnlinien mit Startpunkten in: in Elementmittelpunkten, um Knoten, in beliebigen Punkten:
Ab 14. Zeile:
bei Start in Elementmittelpunkten: Elementnummern (max. 10 Stück pro Zeile)
bei Start um Knoten: eine Zeile mit Anz. der Startpunkte pro Knoten u. Radius (m), weitere Zeilen mit Knotennummern (max. 10 Stück pro Zeile)
bei Start in beliebigen Punkten: je eine Zeile mit (x,y) (2D) bzw. (x,y,z)-Koord. (3D-Modell) für jeden Startpunkt
Hier findet sich eine vollständige Beschreibung der Batch-Datei für eine Strömungsberechnung mit SITRA
1. Zeile: Name der Ausgabedatei (beginnend ab 1. Zeichen!)
2. Zeile: Flag für ausführliches Protokoll (1=ja/0=nein)
3. Zeile: Flag für Gleichungslöser: 0=direkt, 1=iterativ, -1= Mehrgitter-
4. Zeile: Flag für Kontrolllinienrechnung (1=ja/0=nein), wenn Flag=1:
Flag für Wahl der Geschwindigkeiten für die Kontrolllinienrechnung:
0 = Mitteln der Geschwindigkeiten
1 = Geschwindigkeiten rechts der Kontrolllinie
-1 = Geschwindigkeiten links der Kontrolllinie
5. Zeile: Flag für die Berechnung der Pecletzahlen und (bei inst. Rechnung) Courantzahlen (1=ja/0=nein), Flag für Übernahme der K-Werte aus EICHEN (1=ja/0=nein)
6. Zeile: 'STROEMUNG' = Nur Strömungsrechnung oder 'TRANSPORT' = Strömungs- und Transportrechnung
7. Zeile: Für 3D-Modell oder Vertikalmodell: Flag für ges./ungesättigte Rechnung (1=ja/0=nein)
Für Horizontalmodell mit UNTE: Flag für die Iteration der Mächtigkeit (1=ja/0=nein), sonst irrelevante Zeile
8. Zeile: 1. Flag für Transportrechnungen: Berücksichtigung von gleichen Konzentrationen/Temperaturen bei GLEI (1=ja/0=nein), 2. Flag nur wichtig bei 3D-Modell oder Vertikalmodell: dichteabhängige Rechnung (1=ja/0=nein), 3. Flag entscheidet über Berücksichtigung der Druck- oder Potentialgleichung (1=Druck/0=Potential)
9. Zeile: Bei ges./unges. Rechnung oder Iteration der Maechtigkeit und/oder dichteabhängiger Rechnung: Anzahl der Iterationen, Dämpfungsfaktor, sonst irrelevante Zeile.
Bemerkung: bei inst. Rechnung mit Iteration der Mächtigkeit wird die Anzahl der Iterationen pro Zeitschritt immer auf 1 gesetzt (keine Iteration!). Der Dämpfungsfaktor sollte dann 1.0 sein (nur update der Mächtigkeiten)
10. Zeile: Flag für stationäre oder instationäre Strömung und Flag für stationären oder instationären Transport: bei 3D-Modell und Vertikalmodell: 0=stationär, 1=instationär
bei Horizontalmodell mit definierter Grundwassererfüllter Mächtigkeit:
0=stationär, 1=instationär (ohne Berücksichtigung von gespannten/ungespannten Verhältnissen), 2=instationär (mit Berücksichtigung von gespannten/ungespannten Verhältnissen)
11. Zeile: 2 Flags für Startwerte der Potentiale (bei ges/unges. Rechnung oder inst. Strömung) u. Konzentrationen/Temperaturen (bei dichteabh. Rechnung oder inst. Transport): 2=Startwerte aus 'null'-Datei ('null'-Datei nicht vorhanden: Startwert = Eich/Akon), 1=Startwerte aus Eichpotentialen bzw. AKON (Eich bzw. AKON nicht vorhanden: Startwert =0), 0=Startwerte = 0.0
ab 12. Zeile: Eingabeparameter für instationäre Berechnung
12. Zeile: Name der inst. Eingabedatei (beginnend ab 1.Zeichen!)
13. Zeile: Flag für Warmstart mit 'null'-Datei (1=ja/0=nein) ('null'-Datei nicht vorhanden Eingabe wird ignoriert!)
14. Zeile: Flag für die Festlegung der instationären Zeitschritte:
0 = Flag für Zeitschritte aus der inst. Eingabedatei, weitere Eingabe: Verkleinerungsfaktor für Zeitschritte, zus. Verkleinerungsfaktor für die Transportrechnung (optional)
1 = Flag für Rechnung mit fester Zeitschrittweite, weitere Eingaben (getrennt durch Blank oder Tabulator):
Zeitschrittweite, Zeiteinheit für die Zeitschrittweite (0=Sek., 1=Min., 2=Std., 3=Tag, 4=Monat, 5=Jahr), zus. Verkleinerungsfaktor für die Transportrechnung (optional)
15. Zeile: Interpolationsflags für Interpolation inst. Randbedingungen an Zwischenzeitpunkten, bzw. explizte Behandlung inst. Mengenrandbedinungen für POTE, KNOT, FLAE, VORF und EFLA
bei POTE, KNOT, FLAE, KONZ, 1KON, VORF und EFLA: 0= implizite Behandlung ohne Interpolation, 1= implizite Behandlung mit Interpolation
bei FLAE, EFLA und KNOT: -1= explizite Behandlung ohne Interpolation, -2= explizite Behandlung mit Interpoltaion, Flag fuer inst. Bergsenkungen (1=ja/0=nein)
16. Zeile: Anzahl Zeitschritte; oder Eingabe der Randbedingungen für Grubenflutung: GRUB POTE oder MENGE
17. Zeile: Flag für das ausführliche Abspeichern des Ergebnisses (Datei out66) zum Fortsetzen der inst. Rechnung (bei Flag=1)
18. Zeile: Flag für Abspeichern der Zwischenergebnisse in den Hintergrund: 0= Abspeichern nach jedem n.ten Zeitschritt, weitere Eingabe: n, 1= Abspeichern von n festen Zeitschritten, weitere Eingabe: n, 2= Abspeichern von n festen Ganglinien, weitere Eingabe: n
bei Flag = 1 in 18. Zeile:
ab 19. Zeile: Zeitschritte in aufsteigender Reihenfolge, max. 10 pro Zeile, ohne Start- (0) und letzten Zeitpunkt!
bei Flag = 2 in 18. Zeile:
ab 19. Zeile: Knotennummern, max. 10 pro Zeile
Letzte Zeile: CONDITION NUMBER: es wird die Konditionszahl berechnet
Stand: 10.11.2023 | 