Hydraulik II Skriptum - Department Wasser-Atmosphäre-Umwelt ...

Universität für Bodenkultur WienDepartment für Wasser-Atmosphäre-UmweltInstitut für Hydraulik und landeskulturelle WasserwirtschaftLV-Nr. 815.300HydrodynamikVortragende: Willibald LOISKANDL und Reinhard NOLZWintersemester 2011/12für den Studiengebrauch an der Universität für Bodenkultur Wien

HYDRODYNAMIK INHALTSVERZEICHNIS Seite IALLGEMEINE ANGABEN ZUR LEHRVERANSTALTUNG ................................................... 1KONTINUUMSMECHANISCHE BEHANDLUNG VON FLÜSSIGKEITSSTRÖMUNGEN 2GRUNDLEGENDE GESETZE DER FLÜSSIGKEITSSTRÖMUNG......................................... 4KONTINUITÄT (GEOMETRISCHES KONZEPT).......................................................................................... 4Transportgleichung .......................................................................................................................... 5BEWEGUNGSGLEICHUNG....................................................................................................................... 5HYDRAULISCHES VERSUCHSWESEN...................................................................................... 8EINLEITUNG........................................................................................................................................... 8ÄHNLICHKEITSBEGRIFFE IM WASSERBAULICHEN VERSUCHSWESEN ..................................................... 9MAßSTABSBEZIEHUNGEN .................................................................................................................... 10DIMENSIONSANALYSE......................................................................................................................... 11BUCKINGHAM'S METHODE (π THEOREM)............................................................................................ 12DIMENSIONSLOSE PRODUKTE ............................................................................................................. 16FEHLERFORTPFLANZUNGSGESETZ....................................................................................................... 17ABFLUSS IN OFFENEN GERINNEN.......................................................................................... 19DIMENSIONSLOSES ENERGIEHÖHENDIAGRAMM FÜR RECHTECKGERINNE ........................................... 20GRAFISCHE ERMITTLUNG DES FLIEßZUSTANDES IN OFFENEN GERINNEN MIT ALLGEMEINERQUERSCHNITTFORM ............................................................................................................................ 21ERMITTLUNG VON PROFILKENNWERTEN............................................................................................. 22Profilkennwerte .............................................................................................................................. 23BESTIMMUNG DER ABFLUSSVERHÄLTNISSE IN NATÜRLICHEN GERINNEN........................................... 24Einteilige Abflussprofile ................................................................................................................. 24Mehrteilige Abflussprofile .............................................................................................................. 24BERECHNUNG DER ABFLUSSVERHÄLTNISSE IN MEHRTEILIGEN PROFILEN........................................... 25PROFILE MIT ZUSAMMENGESETZTEN RAUHIGKEITEN.......................................................................... 26GEGLIEDERTES QUERPROFIL............................................................................................................... 27KORREKTURBEIWERTE........................................................................................................................ 28GÜLTIGKEITSBEREICH DER GAUCKLER/MANNING/STRICKLER-FORMEL NACH NAUDASCHER (1987) 29EINFLUSS DER VEGETATION................................................................................................................ 30GLEICHUNG DER SPIEGELLINIE BEI STATIONÄR UNGLEICHFÖRMIGEM ABFLUSS ................................. 33INSTATIONÄRER ABFLUSS IN OFFENEN GERINNEN....................................................... 35HYDROMECHANISCHE GRUNDLAGEN FÜR DIE MEHRDIMENSIONALE ABFLUSSMODELLIERUNG .......... 35MATHEMATISCHE GRUNDLAGEN ........................................................................................................ 36TURBULENZMODELLIERUNG ............................................................................................................... 39STATISTISCHE TURBULENZMODELLIERUNG ........................................................................................ 39MODELLIERUNGSANNAHMEN.............................................................................................................. 40Flachwasserapproximation ............................................................................................................ 40

HYDRODYNAMIK INHALTSVERZEICHNIS Seite IIIGerinnehydraulik.......................................................................................................................... 110Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 2KONTINUUMSMECHANISCHE BEHANDLUNG VONFLÜSSIGKEITSSTRÖMUNGENDa die Materie (Fluide) aus Molekülen besteht, könnte die Bewegung der Moleküle, ihre gegenseitigeBeeinflussung und die Einwirkung äußerer Einflüsse die Grundlage der Fluidmechanik bilden.⇒ große Komplexität ⇒ großer AufwandLösung: Materie nicht diskret (molekular) sondern kontinuierlich im Raum verteilt angenommen.Eine sinnvolle Definition der Dichte in einem Punkt ist nurdann möglich, wenn um den Punkt ein kleines Volumen dVangenommen wird.Δmdmρ= lim =ΔVdVΔV→0dm....Masse aller Moleküle in ΔVDie Größe eines Elementarvolumens (für ρ = const.) ist nach oben und unten begrenzt.Elementarvolumen = FlüssigkeitsteilchenRichtwert:V = 10 -9 mm 3 (für Gase bei p a = 101330 Pa)= 2,687 10 7 Moleküleρ wird dem Schwerpunkt von dV zugeordnet. Im Raum: ρ = ρ( xyzt , , , )Die Betrachtung als Kontinuum ist dann möglich, wenn das Elementarteilchen wesentlich kleinerals das Strömungsgebiet ist.Diese Betrachtung kann auch für andere Eigenschaften (skalare oder vektorielle) angewandt werden,z. B. Geschwindigkeit.Ideales Gas (Gaskinetik): zusammengesetzt aus Massepunkten mit "willkürlicher" Geschwindigkeitu i (nach Größe und Richtung (Θ i )).1u =NN∑ u ii=1cos Θ 1N .... Anzahl der Moleküle, die durch die Fläche in der Zeiteinheit Δthindurchgehen.u ..... gemittelte Geschwindigkeit des Fluidums in Richtung derFlächennormalenWird die Fläche so gedreht, dass u ein Maximum wird, so erhält man die sogenannte Fluidgeschwindigkeitu als vektorielle Größe u = Max( u), wobei u als kontinuierlich in Raum und Zeitbetrachtet wird.Die so gefundene Geschwindigkeit ist neben der Lage der Fläche noch von der Größe dieser undvom Zeitintervall Δt abhängig.1) Molekülbewegung kleinster makroskopischer MaßstabInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 4GRUNDLEGENDE GESETZE DER FLÜSSIGKEITSSTRÖMUNGKONTINUITÄT (GEOMETRISCHES KONZEPT)Viele hydraulische Prozesse sind mit Transportphänomenen verbunden.Bilanz in der Zeiteinheit dt:Einströmende Masseρ Q dtρQ+∂ ρQ ∂x dx dt[ ( ) ]Ausströmende Masse ( )Masse am Anfang m a = ρ A dx (Anfangsvolumen A dx) und nach einem Zeitschritt dt= [ ρA+ ( ∂( ρA)∂t)dt]dxDa die Masse erhalten bleibt (Erhaltungssatz der Masse) istm eEinströmung − Ausströmung = Speicherung (ΔS)ρ Qdt ∂∂−ρQdt − ( ρQdxdt) = ρAdx+ ( )∂x∂tρ Adtdx−ρ Adx∂∂− ( ρQdxdt) = ( )∂x∂tρ Adtdx∂( ) ( )∂ ρ ∂A +t ∂xρ Q =0Für konstante Dichte:∂A∂t∂Q+ = 0 Kontinuität, Transportgleichung∂xoder mit konstanter Breite:∂h∂t∂q+ = 0∂xStationäre Bedingungen (Einströmung = Ausströmung; ρ =const.)Q= ∫ vdAv.......mittlere GeschwindigkeitQ......DurchflussInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 5Für die Querschnitte (1) und (2) folgt:vA 1 1 = vA 2 2 = QKontinuität für stationäre Bedingungenund inkompressible FlüssigkeitIn der praktischen Anwendung wird häufig die Stromröhre über das ganze Strömungsgebiet angenommen.Die Randbedingungen der Stromröhre sind dann gleich den physikalisch vorgegebenenBerandung des Strömungsgebietes, z. B. Rohrwand, Flussbett usw.Herleitung der Gleichung für den Stofftransport:Bilanz in der Zeiteinheit dt z. B. für Verunreinigung (c)Einströmende Masse:c Q dt[ cQ + ( ∂ cQ ∂x)dx]dtAusströmende Masse: ( )Masse am Anfang ma = c A dxund nach dt m E = c A + ( ∂( cA)∂ t)dt dxcQdt[ ]Ein - Aus = Δ S∂∂− cQdt −∂x∂t∂ ∂( cA) + ( cQ) = 0∂t∂xc∂AA∂cc∂QQ∂c+ + + = 0∂t∂t∂x∂x⎡∂A∂Q⎤ ∂c∂cc⎢+ + A + Q = 0⎣ t x ⎥14243∂ ∂ ⎦ ∂t∂x( cQ) dxdt = cAdx + ( cA) dtdx − cAdxKontinuität = 0∂c∂cA + Q = 0∂t∂x∂c∂c+ u = 0∂t∂xQmit u =ATransportgleichungc = c(x,t)BEWEGUNGSGLEICHUNGDie Grundlage bildet das Newton'sche Kraftgesetz F = m aWird nun diese Gleichung in Richtung der Kraft bzw. der Beschleunigung integriert so folgt,s2s1s21∫Fds= m∫ads = m v2s12 2( − v )die Energiegleichung.Die an einem Körper geleistete Beschleunigungsarbeit ist gleich der Änderung seiner kinetischenEnergie.21Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 7Abb. Bewegung - FluidelementEntsprechend dem Volumen zwischen 1-1' wird der Impuls des betrachteten Flüssigkeitsvolumenverringert2( ρ Δx) v = ρ A v Δt= QρvtA1 1 1 1 11Δ14243und entsprechend dem Volumen Abschnitt 2 2 ′ ( Q ρ v2 Δ t)m− vergrößert.Die Änderung des Impulses muss durch eine Kraft ausgeglichen werden, daher ist( )F= ρQΔt v2 −v1 .Bei dieser Betrachtungsweise ist die Änderung des Impulses nur vom Ort und nicht von der Zeitabhängig.Für instationäre Verhältnisse muss die Gleichung um ein Glied der Form∂( )dV∂t ∫ ρverweitert werden.Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 8HYDRAULISCHES VERSUCHSWESENEINLEITUNGGrundsätzlich kann zwischen mathematischen Modellen und Modellversuch (physikalische Modelle)unterschieden werden. Der Fortschritt im Bereich der EDV, führte zu einem verstärktem Einsatzvon mathematischen Modellen.Mathematische Modelle können dann angewandt werden, wenn die Problemstellung genügenddetailliert mathematisch beschrieben werden kann. Modellversuche sind möglich, falls die grundlegendenphysikalischen Phänomene in verkleinertem Maßstab reproduziert werden können.Entsprechend den Anwendungsbereichen kann folgende Modelleinteilung erfolgen:0-dimensionale Modelle:1-dimensionale Modelle:2-dimensionale Modelle:3-dimensionale Modelle:Weitere Anwendungen:'lumped system'- einfache Speicherprobleme- Bevölkerungswachstum.- Abfluss in offenen Gerinnen- Rohrhydraulik- Bewässerungssysteme- Dammbruchwellensimulation- Hochwasserwellenabfluss in Gerinnen und/oder Speichern.- Details im Abfluss in offenen Gerinnen- Umströmung, z. B. einer Schwelle- vertikale Beschleunigung (short waves)- Hauptanwendungsgebiet der Modellversuche- spezielle hydraulische Probleme, Einbauten in Gerinnen- Turbulenz- Stofftransport und Dispersion- Wärmetransport und Dispersion- Grundwasserströmungen- Morphologische Änderungen.Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 9ÄHNLICHKEITSBEGRIFFE IM WASSERBAULICHEN VERSUCHSWESENa) Geometrische Ähnlichkeit:Festes Verhältnis der messbaren, homologen geometrischen Größen zweier Systemeb) Kinematische Ähnlichkeit :Festes Verhältnis der in zwei Systemen verschiedener Abmessung homologer Zeitenc) Dynamische Ähnlichkeit:Festes Verhältnis der messbaren homologen Kräfte bzw. Kraftwirkungen zweier Systemed) Ähnlichkeit der Stoffübertragungsprozesse (z. B. Sauerstoffaufnahme)e) Reaktionskinetische Ähnlichkeit (z. B. Reaktionsverlauf und Auswirkung)f) Biologische Ähnlichkeit (z. B. Biologische Fermentationsprozesse)Eine vollständige Ähnlichkeit zweier Systeme ist nicht erzielbar, so dass man sich mit partiellerÄhnlichkeit begnügen muss. Es können daher nebeneinander ablaufende physikalische, chemischeund biologische Vorgänge im Gesamten modellmäßig nicht nachgebildet werden, im Rahmen desModellversuchwesens erfolgt daher eine Beschränkung auf die Behandlung von Teilproblemen.Kalibrierung:Verifizierung:Modellanpassung an gemessene Prototypdaten. Modell reproduzierteine spezifische Situation.Vergleich der Modellergebnisse mit einer bekannten Situation desPrototyps, ohne das Modell zu verändern.Kalibrierung und Verifizierung bedingen Prototypdaten.Der Maßstab eines Parameters ist definiert als Verhältniszwischen Prototypwert und Modellwert dieses Parameters.z. B.: Längenmaßstab nL=LLpmDurch diese Festlegung ist der Maßstab meist > 1.Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 10MAßSTABSBEZIEHUNGENMaßstabsgesetz (Scale laws): Der Maßstab eines Produkts von Parametern ist gleich dem Produktder Maßstäbe dieser Parameter.RepvpDp νmnvnLz. B.:nRe= ==Re ν v D nmpmWenn nun in einer Fragestellung die Trägheitskräfte und die zähen Reibungskräfte überwiegen, somuss das Verhältnis dieser Kräfte in Natur und Modell gleich sein.Dies führt zu einer Maßstabsbedingung (scale condition).−1Re v L =z. B.: n = n n nν 1"Reynold'schesmνÄhnlichkeitsgesetz"Im Prinzip sind Abweichungen von den Maßstabsbedingungen möglich oder unter Umständennotwendig, jedoch treten dabei Maßstabsfehler auf.In vielen Fällen von praktischen Anwendungenwerden die Maßstäbe aus der Zusammenfassungvon Parametern bestimmt.Beispiel: stationärer Abfluss mit freier Oberfläche (Reibung vernachlässigt).2Bernoulli H = h + v 2g123hgMaßstabsbestimmungHphp+ hgpnh+ hgphmnh+ nhg( hgmhm)nH= = ==H h + hg 1+hg h 1+hg h1.) Für nh= nhg: → nH= nh= nhgAlle Längen haben gleichen Maßstab.mmm2.) Für nh ≠ nhg nH = f( nh, nhg und hgm hm)mn H ist also nicht nur eine Funktion von anderen Maßstäben sondern zusätzlich von dem Parameterhg m h m abhängig. Daher treten Maßstabsfehler auf.n = n folgt, da ng = 1 istFür den Fall 1 ( H hg)2p2mmv 2gn h = = nv 2goder nv=2vnhmmInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 11oder mit der Froudezahl3 22Trägheit ρLv / L v 2== = Fr3Schwerkraft ρgLgL( )( ) voder Fr = mit c = ghc..........Wellengeschwindigkeit1/ 2vpghm−1/2für n Fr =1 folgt 1 == n1/ 2vnLghpvmAnmerkung:1) Betrachtung der Kräfte (Hydrodynamik)als Bedingung für die Ähnlichkeit liefert als Nebenproduktnoch einen Einblick in den Maßstabsfehler falls n h ≠n hg ist.2) Soll zusätzlich die Reibung berücksichtigt werden, müsste auch das Reynold’sche Modellgesetzerfüllt sein.vpLp−1nRe= 1 = ⇒ nv= nLn ν =1 (gleiche Flüssigkeit im Prototyp und Modell)v LmTriviale Lösung n v = n L = 1mghDIMENSIONSANALYSEDer Zweck der Dimensionsanalyse ist es, die Anzahl der Variablen, die in einem Problem auftreten,durch Bildung von untereinander unabhängigen dimensionslosen Produkten zu reduzieren. Diedimensionslosen Produkte, die anstelle der einzelnen Variablen auftreten, enthalten neben denVariablen auch alle dimensionsbehafteten physikalischen Konstanten, die für das Problem relevantsind. Die Zusammenfassung einzelner Variablen zu Produkten reduziert die Anzahl von Experimentenund erleichtert die Interpretation von theoretischen Überlegungen.Bemerkung: Ist eine Größe eine Funktion von nur einer Variablen, dann kann das Ergebnis einerVersuchsreihe (z. B. 10 Experimente) durch eine einzige Kurve dargestellt werden. Ist die Größevon zwei Variablen abhängig, so wird das Ergebnis durch eine Kurvenschar dargestellt. Für z. B. 10benötigte Kurven sind 100 Experimente notwendig. Drei Variablen führen zu Seiten mit Kurvenscharen.Bei 10 benötigten Seiten steigt die Anzahl der Experimente auf 1000, usw.Grundlage der Dimensionsanalyse ist das Prinzip der Homogenität in den Dimensionen:Eine Gleichung, die einen physikalischen Vorgang richtig darstellt, ist homogen in den Dimensionen.Es wird erwartet, dass die Lösungen dieser Gleichung auch homogen in den Dimensionen sind.Die Dimensionsanalyse liefert partielle Lösungen der Probleme, dabei ist es nicht notwendig dieGleichungen selbst zu kennen. Es müssen aber alle Variablen und physikalischen Konstanten undnur die, die den physikalischen Vorgang beeinflussen, berücksichtigt werden.Eine Gleichung ist homogen, falls sie unabhängig von den verwendeten Grundeinheiten ist.Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 12Beispiel: Überfall über ein WehrPoleni-Formel: Q =μb 322g hQ.......Abfluss [L 3 /T]b........Wehrbreite [L]g........Erdbeschleunigung [L/T 2 ]h........Wasserspiegellage über der Wehrkrone [L]μ .......dimensionslos solange für L,T die gleichenGrundeinheiten verwendet werdenAndererseits sind verschiedene Gleichungen in ihrer allgemeinen Form nicht homogen (z. B.Stricklerformel).Der wichtigste Schritt ist die Auswahl der unabhängigen Variablen, die das betrachtete Problembeeinflussen. Anschließend sind die dimensionslosen Produkte zu bilden.In der Hydraulik werden 3 Arten von Variablen unterschieden, bezogen auf:- die Geometrie,- den Durchfluss,- das Strömungsmedium.BUCKINGHAM'S METHODE (π THEOREM)Die Beziehung des betrachteten Problems mit den Variablen a,b,c....( a,b,c, K ) 0f =wird in eine einfachere Beziehung mit einer kleineren Anzahl von variablen dimensionslosen Argumentenπ 1 ,π 2 ,..., die aus den Variablen, die das Problem beschreiben, hergeleitet werden, transformiert.( π , π , π , K ) 0F 1 2 3 =Die Anzahl der Grundeinheiten in der Hydraulik ist r ≤ 3 (Länge, Masse, Zeit).Die dimensionslosen Argumente π1,π2,....,πn-r sind Produkte von Variablen mit variierenden Potenzen.1 1 1π 1 = a x b y czx y z2 2 2π 2 = a b cJedes Argument π sollte ( r +1 ) Variable enthalten, wobei 2 Bedingungen erfüllt sein müssen:1. Alle Grundeinheiten müssen enthalten sein2. Die Variablen selbst dürfen kein dimensionsloses Argument seinIm allgemeinen Fall mit r = 3 geschieht dies indem 3 von 4 Variablen in jedem Argument wiederholtwerden, z. B. charakteristische Länge, Geschwindigkeit und Dichte. Die 4. Variable ist injedem Argument verschieden (mit Exponent ±1), so dass in der Lösung alle n-Variablen enthaltensind.etc.Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 13Im allgemeinen sind dimensionslose Argumente π einfache Zahlen.[ ] [ A] [ B] [ C] [ N]x y z vπ = = 1Werden nun die Dimensionen A,B,C durch die Grundeinheiten, z. B. L,M,T ausgedrückt und dieExponenten der Grundeinheiten aufsummiert, so ist die Summe gleich 0, da das Produkt der Grundeinheitenπ = 1 ist.Es ergeben sich r Gleichungen für die unbekannten Exponenten x,y,z. Die verbleibenden unbekanntenExponenten sind frei (Erfahrung) wählbar. Als Ergebnis können alle dimensionslosen Ausdrückebestimmt werden.Beispiel 1: Widerstand eines Körpers (z. B. Schiffs), der sich mit gleichbleibender Geschwindigkeitauf einer unendlich ausgedehnten Fläche, einer idealen Flüssigkeit mit unendlicher Tiefe, bewegt.Variablen: Widerstand R (MLT -2 )Geschwindigkeit v (LT -1 )Erdbeschleunigung g (LT -2 )Dichte ρ(ML -3 )Länge des Körpers l (L)→ n=5, r=3 (M,L,T)→ (n-r) dimensionslose Argumente können gebildet werden.x1 y1 z1 1π1= l v ρ gx2 y2 z2 1π 2 = l v ρ RVariable durch Grundeinheiten L,M,T ausdrücken und aufsummieren der Exponenten.x1π 1 = LFür L x1 + y1 − 3z1+ 1 = 0T −y1− 2 = 0Mz1= 0x1 = 1, y1 = − 2,z1= 0−2glπ 1 = 1vg =2vFür L x2 + y2 − 3z2+ 1 = 0T −y2− 2 = 0Mz2+ 1 = 0x =− 2, y =− 2,z =−12 2 2π2x2π 2 = L−2 −2 −1⎛ L⎞⎜ ⎟⎝ T⎠⎛ L⎞⎜ ⎟⎝ T⎠= 1 v ρ R =y1y2⎛⎜⎝⎛⎜⎝M⎞3⎟L ⎠M⎞3⎟L ⎠R2 2l v ρz1z2⎛ L ⎞⎜⎝ 2⎟T ⎠⎛⎜⎝1(Kehrwert von Fr 2 )ML⎞2⎟T ⎠1(Newtonzahl Ne)F R =vghInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 14⎛ R g1Der Zusammenhang zwischen den dimensionslosen Argumenten F2 2 201 v ρ ,⎞⎜ ⎟ =⎝ v ⎠kann experimentell gefunden werden, dazu wird ein Modell (Länge l) mit verschiedenen Geschwindigkeitenv gezogen und der dazugehörige Widerstand R gemessen.Für Wertepaare −v, R → Ne , Fr2 (berechnet) - kann ein Graph aufgestellt werden als FNeFe ( , ).Wird nun ein geometrisch ähnlicher Körper, mit Länge l p mit der Geschwindigkeit v p unter den2 2gleichen Bedingungen wie vorher in einer Flüssigkeit gezogen und die Froudezahl ( Frp= v p g lp)berechnet, so kann mit Hilfe des Graphen die zugehörige Newtonzahl Ne p bestimmt werden.Für den Widerstand folgt R p = Ne p l 2 p v 2 p ρ.Beispiel 2: Bestimmung der Reibung eines Partikels mit Durchmesser dL das sich mit der GeschwindigkeitvLT −1 in einen inkompressiven Medium unendlicher Ausdehnung−3−1−1( [ ML ] und μ [ ML T ])ρ bewegt.Variable: Widerstand R (MLT -2 )Geschwindigkeit v (LT -1 )Dichte ρ (ML -3 )Zähigkeit μ (ML -1 T -1 )Durchmesser d (L)n = 5, r = 3 (M,L,T)(n-r) dimensionslose Argumente können gebildet werden.ππx1 y1 z1 11 = ρ v d μx2 y2 z2 12 = ρ v d RVariable durch Grundeinheiten L,M,T ausdrücken und aufsummieren der Exponenten−x1−1y1z − −( ML ) ( LT ) L ( ML T )3π 1 =1 1 1 1LTM− 3x1+ y1Für:+ z1−1− y1−1x1+ 1x = −1,y= −1,z= −11−1−x−y z −( ML ) ( LT ) L ( MLT )1−1−11= 0= 0= 0( ρvd) 1 Reπ1 = ρ v d μ = μ =3 1 1 1 1 2π 1 2 = = 0Für:FürL − 3x + y + z + 1 = 0TM2 2 2−y− 2 = 02 2 22x2+ 1 = 0x =− 1, y =− 2, z =−2( )π = ρ v d R = R ρv d2−1 −2 −2 2 2⎛Der Zusammenhang zwischen den dimensionslosen Argumenten F vD Rν , ⎞⎜2 2⎟ 0⎝ ρ v d ⎠=Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 15kann experimentell gefunden werden und in Form einer einzigen Kurve als Funktion der Reynolds-Zahl dargestellt werden. Mit Re = vDν(Reynolds-Zahl) und C D (Widerstandskoeffizient)2wird der Querschnitt der Kugel A = 14π D und ρv 2 2 die kinetische Energie der Strömung zurFestlegung von CD berücksichtigt, so ergibt sich für R aus der Newton Zahl2 2 1 2 π 2R = ρvd f ( Re) = ρv( d ) CDwobei CD = 8 f( Re ) ist.2 123 4πAAnmerkung: Bei der Darstellung von experimentell ermittelten Daten besteht die Gefahr einer„spurious correlation“, deshalb ist zu beachten, dass die abhängige Variable nur in einemdimensionslosen Produkt π vorkommt.Widerstandskoeffizient versus Reynoldszahlfür verschiedene KörperWiderstandskoeffizient versus Reynoldszahlfür verschiedene FormbeiwerteInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 16DIMENSIONSLOSE PRODUKTEFroude ZahlFr2v=gLTrägheitGravitationρvL vL TrägheitReynolds Zahl Re = =μ νReibungskräfteEuler ZahlEu= ρ v2 v2,Pa gLTrägheitDruckkraftWeber ZahlWevL= ρ 26TrägheitKapillarkräfteStrouhal ZahlSt= ω L Lv, vTOszillationMittlere GeschwindigkeitKavitations ZahlCa=P − Pρv2vDruckkraftTrägheitCauchy-Mach'sche ZahlMa=vEρelastTrägheit. Formänderungskräftep v : DampfdruckInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 17FEHLERFORTPFLANZUNGSGESETZAuswertung von MessergebnissenMessungen physikalischer Größen sind im Normalfall mit Fehlern behaftet (Messwert weicht vom"wahren" Wert ab). Die Fehler sind entweder systematischer oder stochastischer (zufälliger) Natur.Systematische Fehler sind durch Eichung der Messgeräte behebbar. Von zufälligen Fehler sindweder das Vorzeichen noch der Betrag angebbar.Grundlegende statistische Größen:Das arithmetische Mittel xn∑ x ii=11x =xi ............ Messwertnn ............. Anzahl der MessungenNach Gauß ist x als Schätzwert für den "wahren" Wert x als Minimum der Summe der Quadratedes absoluten Fehlers festgelegt (Methode der kleinsten Quadrate).Die Standardabweichung σ i (mittlerer Fehler der n Einzelmessungen)σ i=∑ f i2 1n −f = x − x Fehler der EinzelmessungiiDer mittlere Fehler des Mittelwertes der Messreihe:∂X∑ fi2= σ =X n n −1∂ ( )oder als relativer Fehlerr = σ 100 %xcvGrößen die nicht direkt gemessen werden, sondern aus Messwerten berechnet werden, sind mitderen Fehlern behaftet. Unter der Voraussetzung einer Normalverteilung und für kleine Standardabweichungen(σi

HYDRODYNAMIK SEITE 18Mit X = F (x1, x2,...)z. B. X = a1 x 1 + a 2 x 2 ......21folgt σ = ( a σ ) + ( a σ ) +xoder für X = a x1 p1 x2 p2 .....1⎛ σ ⎞ ⎛ σ ⎞folgt σ =⎜1⎟ +⎜2X X p1p2⎟ +⎝ x1⎠ ⎝ x2⎠oder mit dem relativen Fehlerrx=p2 21 r x1+ p2 22rx22+K2222Beispiel:Dreieckswehr (Thomsonwehr)für kleinere bis mittlere DurchflüsseQ = Cvtan α h52 Cv ≈ 14 , m s, α meist 45°Wird der Fehler für α vernachlässigt, so folgt für den relativen Fehler für Q( 52)2 2 2 2rQ = rCv + rhr Cv ..........von der Eichungr h ............von MessungFür die Kalibrierung sind Q und h zu messen:QCv =h 52tan α( 52)Cv 2 Q 2 2 h2r = r + − r12Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 19Definition der verwendeten Variablen.ABFLUSS IN OFFENEN GERINNENNach Bernoulli:2v pEnergiehöhe H = + + z =2g ρg123hPmit z hρ g+ = und v Qm = ABestimmung des Fließzustandes:Q2ga22( h)+ hIm allgemeinen ist H = H(Q,h); für konstantes Q folgt H = H(h).EnergiehöhendiagrammInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 20DIMENSIONSLOSES ENERGIEHÖHENDIAGRAMM FÜR RECHTECKGERINNERechteckquerschnitt: A = B h und Q = q BqEnergiehöhe H = 22gh+h2Für den kritischen Abfluss folgt mitVhgrmgr=h mgr = 3g hmgr2 2gr Q2 2gB hmgr2V= =gH min = 3 2qgh grFür den allgemeinen Fall H = H(Q,h) benötigt man zur Auswertung der Energiehöhengleichungeine Kurvenschar. Durch eine dimensionslose Darstellung des Energiehöhendiagramms kann dieKurvenschar durch eine einzige Kurve ersetzt werden.H hgrh= ⎛ 21h ⎝ ⎜ ⎞⎟ +2 h ⎠ hgrgrInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 21GRAFISCHE ERMITTLUNG DES FLIEßZUSTANDES IN OFFENEN GERINNEN MITALLGEMEINER QUERSCHNITTFORMInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 22ERMITTLUNG VON PROFILKENNWERTEN"Streifenwerten" durch Aufsummierung gewonnen:(a) Einzelfläche im Inneren:1A = ( h − z ) + ( h−z + 1)x + 1 −x2[ ]( )i i i x ibzw. an den Rändern (z. B. links):1Ai = h z2( − )2xz− x− z2 2 11 2(b) Benetzter Umfang im Inneren:2( ) ( )2i i+ 1 i i+1 iU = x − x + z −zbzw. an den Rändern (z. B. links):Ui=( h − z )22+⎡⎢⎣( h − z )( x − x ) 221z22− z1⎥ ⎦⎤Gesamtfläche und gesamter benetzter Umfang:nA = ∑ A iU = ∑i=1nU ii=1Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 23ProfilkennwerteProfildaten dargestellt in tabellarischer Form:h U(h) A(h) R(h)SUBROUTINE ITER (F,FI,NN,TABV,FTAB,ITAB)C DATE: 88-12-08C LINEAR - INTERPOLATION OF DATASETCCPARAMETER LIST:C F INDEPENDENT INTERPOLATIONVALUEC FI DEPENDENT INTERPOLATIONVALUEC ITAB PARAMETERC NN NUMBER OF DATAPAIRSC TABV INDEPENDENT TABULATED VALUESC FTAB DEPENDENT TABULATED VALUESCREAL F,FI,TABV(ITAB),FTAB(ITAB)INTEGER NN,JDO 10 J=1,NNIF(F .GE. TABV(J) .AND. F .LT. TABV(J+1))THENFI=FTAB(J)+(FTAB(J+1)-FTAB(J))/(TABV(J+1)-TABV(J))1 *(F-TABV(J))GOTO 99999END IF10 CONTINUE99999 RETURNENDInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 24BESTIMMUNG DER ABFLUSSVERHÄLTNISSE IN NATÜRLICHEN GERINNEN(aus Hydraulik und Hydrodynamik)Einteilige Abflussprofilea) Kompaktes Profil⎛ h 1 ⎞b) Breites Profil ⎜ < Rh→h⎟⎝ B 30 , ⎠Mehrteilige AbflussprofileInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 25BERECHNUNG DER ABFLUSSVERHÄLTNISSE IN MEHRTEILIGEN PROFILENBerechnungsverfahren nach Schmutterer für den Abfluss im Vorland. Das Vorland wird als Ersatzflächein die Berechnungen einbezogen.Für den Abfluss im Vorland folgt: Qv = Avvv = As'vsAs'Ersatzfläche für das Vorland, bei welcher der Abfluss Qv mit der mittleren Abflussgeschwindigkeitvs des Stromschlauches stattfindet.Zufolge hManning/Strickler:woraus sichMit IvA= bzw. h = Rh gilt für breite GerinneBA′= As= I bzw. I Ivssvkksvss⎛ hv⎞⎜h⎟⎝ s ⎠kk⎛⎜⎝hB

HYDRODYNAMIK SEITE 26PROFILE MIT ZUSAMMENGESETZTEN RAUHIGKEITEN(aus Hydraulik und Hydrodynamik)Für einfache Profilformen, aber mit unterschiedlichen Rauhigkeiten entlang des benetzen Umfanges,ohne jedoch eine Unterteilung der Abflussfläche einzuführen, kann durch Bestimmung einesäquivalenten Abflussbeiwertes wie mit einer konstanten Rauhigkeit gerechnet werden.Zur Bestimmung des äquivalenten Rauhigkeitsbeiwertes (z. B. k ST -Wert) wird die durchflosseneQuerschnittsfläche in n-Teile geteilt, von denen der benetzte Umfang U1, U2.....Un und die Rauhigkeitsbeiwertek ST1 , k ST2 .... k STn bekannt sind. Horton und Einstein trafen die Annahme, dass jedeTeilfläche mit der gleichen mittleren Geschwindigkeit durchflossen wird( v v = = v v)= K .1 2 n =Durch die Annahme folgt für den äquivalenten Rauhigkeitsbeiwert:kSTä=⎡⎢⎢⎢⎢⎢⎣n∑i=1UiU( k )3 2STi⎤⎥⎥⎥⎥⎥⎦2 3λä=n∑i=1λ UiUiInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 27GEGLIEDERTES QUERPROFILIst die Berechnung mit einem äquivalenten Abflussbeiwert nicht mehr möglich, wie bei stark gegliedertenQuerprofilen oder bei Querprofilen mit unterschiedlich verteiltem Bewuchs, dann ist eszweckmäßig die Querprofile in einzelne Bereiche zu gliedern, für die eine annähernd gleiche Geschwindigkeitsverteilungund Rauhigkeit angenommen werden kann (z. B. Vorland-Hauptstromschlauch).Die Berechnung des Normalabflusses erfolgt dann mit23 121 ST 1 123 122 ST 2 2v = k R I , v = k R I , ,v = kiSTiR2 3 1 2iIund der Kontinuitätsgleichungn∑Q = A i v i n..... Anzahl der Bereichei=1Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 28KORREKTURBEIWERTE(aus Hydraulik und Hydrodynamik)Im allgemeinen ist die Geschwindigkeit nicht gleichmäßig über den Querschnitt verteilt. Korrekturbeiwertfür die Geschwindigkeitshöhen α (Coriolis Beiwert).31 ∫ v dAAα =vm........ mittlere Geschwindigkeiten,A v3mentsprechend Q/AWird nun die Geschwindigkeit v i der einzelnen Streifen, unter der Annahme eines konstanten EnergieliniengefälleI für den gesamten Querschnitt, herangezogen so gilt für jeden Streifen:v = kiSTiR2 3i1 2Iund der α-Wert folgt näherungsweise mitα =n∑ 1 i=1A3( v A )vi3miAnmerkung: Die Geschwindigkeitsverteilung über den Querschnitt wird analog zu den einzelnenWassertiefen der Streifen angenommen.Boussinesq-Beiwert.β =1An2 viAi∫ v dA ∑A 1 i=1≈2 2v mA v mβ ergibt sich aus der Integration über den Abflussquerschnitt (ähnlich wie α) → 1 je gleichförmigerdie Geschwindigkeitsverteilung ist.Für den ganzen Querschnitt:α =β =1An∑i=1n1A1A1Aα∑i=1n∑i=1nβ∑i=1ivivv3im3imv2imA2imAAiAiiiInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 29GÜLTIGKEITSBEREICH DER GAUCKLER/MANNING/STRICKLER-FORMELNACH NAUDASCHER (1987)a) GMS-Beiwerte kSt sind nicht dimensionslos, deshalb gültig nur bei Gerinnegrößen, für die kStermittelt wurden,b) Zähigkeitseinfluss ist nicht berücksichtigt, deshalb gültig nur für sehr große Reynolds-Zahlen,c) Widerstandsgesetz ist nicht beachtet, deshalb gültig nur für mittlere relative Rauhigkeiten,d) Querschnittsform ist nicht berücksichtigt, deshalb gültig nur für Formen des Gerinnequerschnitts,für die kSt ermittelt wurde,e) Einfluss unterschiedlicher Rauhigkeit und Gliederung des Gerinnequerschnitts sind nicht berücksichtigt,f) Einfluss des Sedimenttransports und der veränderlichen Sohlenform bei beweglicher Sohle sindnicht berücksichtigt,g) Einfluss der zur Wellenbildung führenden Instabilität ist nicht berücksichtigt,h) Einfluss der Luftaufnahme bei extrem hohen Geschwindigkeiten ist nicht berücksichtigt.Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 30EINFLUSS DER VEGETATIONVegetationsbestände bewirken veränderliche Rauhigkeiten.Widerstand einer Pflanze Fw = 0,5 Cw Ap ρ Vm,p 2Vm,p.....................mittlere Anströmgeschwindigkeit [m/s]Cw und Ap ...........f(Form der Pflanze).τ′′ =FwawWand- bzw. Sohlenschubspannungim Sinne eines Formwiderstandesa ...... Elementabstand in Fließrichtungw ..... Elementabstand quer zur Fließrichtungτ = τ'+ τ"τ'........... Wandanteil der Rauheitselemente (Flächenwiderstand)Klassifizierung des Bewuchses bei natürlichen oder naturnah gestalteten Gerinnen nach Schröder(1990),Brettschneider und Schulz(1985).Großbewuchs− Stammbewuchs über Wasserspiegel meist relativ starr− hydraulisch durch Zylinder mit äquivalentem Durchmesser genähert.Mittelbewuchs (Buschbewuchs):− maximale Höhe bis zum Wasserspiegel− flexibel− hydraulisch noch als Einzelelemente zu identifizieren.Kleinbewuchs (Gräser):− Höhe klein gegenüber der Wassertiefe− sehr flexibel− hydraulisch als flächenhafter Rauhigkeitsbelag betrachtetMikrobewuchs (Algen):− Höhen gegenüber Wassertiefe vernachlässigbar− Korngerüst der Sohle verklammernd− hydraulisch als flächenhafter Rauhigkeitsbelag betrachtet.Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 31Großbewuchs (nach PETRYK/BOSMJIAN)kst=1+∑ickw,i2gALst,WAvP,ik2st, wR4 3A..........Abflussquerschnittkst,W ...Strickler-Beiwert der gleichen aber unbewachsenen Wand [m1/3/s]cw,i......Widerstandskoeffizient der Pflanze i [-]Ap,i......angeströmte Fläche der Pflanze i [m]LV .......charakteristische Länge einer Vegetationszone in Fließrichtung [m]Anmerkung: Einfluss der Pflanze und des Sohlenmaterials berücksichtigt.Orientierungshilfe für Cw-Werte für Großbewuchselemente mit äquivalentem Zylinderdurchmesserd.ElementarordnungRelativer Abstand a/dim Gerinne 5 10 20 30Einzelne Reihe längs, mittig 0,44 0,60 0,87 1,06angeordnetVersetzt (gestaffelt)0,58 0,77 1,02 1,17auf halber GerinnebreiteVersetzt (gestaffelt)auf ganzer Gerinnebreite0,79 0,95 1,16 1,23Tabelle: C w -Werte nach Garbrecht, G.,Abflussberechnungen für Flüsse und Kanäle.Die Wasserwirtschaft S1 (1961) Nr. 2, S 40 - 45Nr. 4, S 72 - 771) C w -Werte Mittelwerte mit Streubreiten bis ± 25 %Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 32Mittelbewuchs (nach Schulz)Der Ansatz ist ähnlich wie bei dem Großbewuchs.kst=1+∑icw,i2gALAvkp,ist,wk2st, wR4 3⎛ v⎜⎝ vm,Pm⎞⎟⎠2Zusätzlich Quotient aus Vm,p/Vmvm,p..... mittlere Anströmungsgeschwindigkeit eines Bewuchselementesvm........ mittlere Fließgeschwindigkeit im GerinneDer Kleinbewuchs kann als einziger als Wandrauhigkeit aufgefasst werden.Mikrobewuchs kann zur Bildung von Filmen oder Matten führen die unter Umständen die Oberflächeglätten.Klassifizierung lebender Baustoffe nach Bretschneider und Schulz (1985)Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 34Δxh = h −( z − z ) + β ( h − h ) + ( I + I ) + h2i+ 1 i i i−1 k, j k, j− 1 E, j+1 E, j v, örtl.hk2v= ,2gIE2 2QQ= oder2K A 2 (1/ λ)8gRfolgt für die einzelnen Größen der Energiehöhe für die gegliederten Querschnitte mit den Teilquerschnitteni die Wassertiefe für die Stationierung j+1⎡⎤nn3 3v2 2ji ,Aji ,vj 1, iA−−⎢∑∑ + j+1, i⎥ 2 ⎡ n nR1 1ji ,R⎤i= i=Q Δx⎛ ⎞ ⎛ ⎞j+1, ii+ 1=i−( i−i− 1 ) + β ⎢ − ⎥+ ⎢ji , + j+1, i ⎥ +vörtl , .2gQ 2gQ 16g⎜∑i= 1 λ ⎟ ⎜∑⎢ji , i= 1 λ ⎟j+1, i⎥h h z z A A h⎢ ⎥⎣⎝ ⎠ ⎝ ⎠⎢⎥⎦⎣⎦Nach Naudascher wird für die Geschwindigkeitshöhehn∑h Qki , i ni=13k= = ∑ vi AiQ 2gQi=11gesetzt.Es wird davon ausgegangen, dass das Energiegefälle in allen Teilquerschnitten gleich ist. Der ß-Wert berücksichtigt eventuelle Verluste in Erweiterungen oder Verengungen und ist beiv j ≥ v j+1 mit β= 1 undv j < v j+1 mit β= 2/3 für allmähliche Aufweitungen kleiner als 1:7 anzusetzen.Die Verluste aus plötzlichen Erweiterungen sind über örtliche Verluste zu berücksichtigen undkönnen zum Beispiel über die Bordasche Verlusthöhengleichung quantifiziert werdenMit den vorgestellten Annahmen folgt die Arbeitsgleichung, mit der der Wasserspiegelverlauf fürgegliederte Profile berechnet werden kann.⎡⎤⎢⎥h h z z ⎣ ( h ) ( h ) ⎦ h2 g 16 g ⎢⎛ ⎞ ⎛ ⎞ ⎥⎢⎜∑Aji ,Rji ,/ λji , ⎟ ⎜∑Aj+ 1, iRj+ 1, i/λj+1, i ⎟ ⎥⎢⎣⎝ i= 1 ⎠ ⎝ i=1⎠ ⎥⎦2 2Q Q Δx1 11 ( 1)i += − − i i i −+ β ⎡α− αj j + 1⎤ + ⎢+ ⎥ +2 2 vörtl , .nnwobeiα( h)3/2n⎛R⎞i∑ Ai⎜i 1 λ⎟= i=⎝ ⎠1/23⎡ n⎛R⎞ ⎤i⎢∑Ai⎜ ⎥i=1 λ⎟i⎢⎣⎝⎠⎥⎦Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 36MATHEMATISCHE GRUNDLAGENDen Ausgangspunkt der theoretischen Überlegungen für die dreidimensionale Strömungsbeschreibungbilden schon seit jeher die Kontinuitätsgleichung und die Impulsgleichungen nach Navier-Stokes. Letztere sind sowohl für die laminare als auch turbulente Bewegung, die den Regelfall imAbflussgeschehen in künstlichen und natürlichen Gerinnen darstellt, geeignet. Eine Grundannahmein der Abflussmodellierung ist ein inkompressibles Strömungsverhalten.∂v∂vx y ∂vzrKontinuitätsgleichung + + = div v = 0∂x∂y∂zImpulsgleichungen nach Navier-Stokes2 2 2∂vx∂vx∂vx∂vx 1 ∂p⎛ ∂ vx∂ vx∂ v ⎞x+ vx+ vy+ vz= X − + νt x y z x⎜ + +2 2 2x y z⎟∂ ∂ ∂ ∂ ρ ∂ ⎝ ∂ ∂ ∂ ⎠2 2 2∂vy∂vy∂vy∂vy 1 ∂p⎛ ∂ vy∂ vy∂ v ⎞y+ vx+ vy+ vz= Y − + ν ⎜ + + ⎟2 2 2∂t∂x∂y∂zρ ∂yx y z⎝ ∂ ∂ ∂ ⎠2 2 2∂vz∂vz∂vz∂vz 1 ∂p⎛ ∂ vz∂ vz∂ v ⎞z+ vx+ vy+ vz= Z − + ν2 2 2t x y z14444244443 { {z⎜ + +x y z⎟∂ ∂ ∂ ∂ ρ ∂14 ⎝ ∂ ∂ ∂444 244443⎠Advektion Kraft Druck Zähigkeitx,y,z.....t ...........v x, , v y , v zX,Y,ZpνμρKoordinatenrichtungenZeitGeschwindigkeitskomponente in KoordinatenrichtungRichtungsvektoren der Kraft KDruckdynamische Zähigkeitkinematische ZähigkeitDichteDie Kraft K(X,Y,Z) beinhaltet die konservativen Anteile wie Schwerkraft, Corioliskraft und weitereäußere Einwirkungen wie z. B. den Windeinfluss.Damit ist ein System von vier Gleichungen zur Lösung von vier Unbekannten v x , v y , v z und p gegeben,das theoretisch auch lösbar ist. Im Falle der turbulenten Bewegung jedoch darf nicht die mittlereGeschwindigkeit eingesetzt werden, es müsste die tatsächliche (instationäre) Bewegung betrachtetwerden. Die turbulente Bewegung ist aber unregelmäßig und zufallsbedingten Schwankungenunterworfen, sodass nur einzelne Zustände beschrieben werden können, die zwar theoretisch möglichsind, jedoch in der Praxis nicht auftreten müssen. Einer der ersten der eine Lösung dieses Problemsvorschlug war BOUSSINESQ, der einen Ansatz in Form einer weiteren Zusatzkraft, die alsFolge einer Scheinzähigkeit entsteht, einführte.Da die Beschreibung der Kraftwirkung der viskosen Spannungen eine zentrale Rolle spielt solldiese, ausgehend von den an einem Volumenelement angreifenden Vektoren U,V,W, näher erläutertwerden. Der innere Spannungszustand infolge Zähigkeit ist, unter Zugrundelegung des Newton’schenZähigkeitsgesetztes, durch einen Spannungstensor, gebildet aus den Normalspannungenund den Tangentialspannungen, bestimmt.Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 37∂Ux − Richtung : dxdydz∂xFür die anderen Koordinatenrichtungen giltanalog:∂Vy − Richtung : dxdydz∂y∂Wz − Richtung : dxdydz∂zσ iτ i,jNormalspannungTangentialspannungKraftwirkung der viskosen Spannung in x-Richtungund zugehörige Komponenten.Durch die Zerlegung der Vektoren U, V, W in ihren Komponenten kann der viskosen Spannungstensorgeschrieben werden:U =V =W =( σx, τxy, τxz)( τyx, σy,τyz)( τ , τ , σ )zxzyzDer erste Index gibt die Richtung der Flächennormalen an, der zweite die Richtung der Spannung.Der Spannungstensor enthält nur sechs Unbekannte, da die Tangentialspannungen τ jeweils paarweisegleich sind, z. B. τyx= τxyusw.. Für die laminare Bewegung – nach Newton gilt z. B.σ x = μ dv/dx – liefern die Richtungsableitungen der einzelnen Komponenten die Zähigkeitsterme inder Bewegungsgleichung. Für die x-Richtung folgt:⎛2 2 2∂σ ∂τ∂ ⎞x xy ∂τxz ⎜∂ v vx y ∂ vz+ + = μ + + ⎟μν =∂ ∂ ∂2 2 2x y z⎝ ∂x∂y∂z⎠ρDie anderen Zähigkeitsterme in y- und z-Richtung ergeben sich analog.Für die turbulente Strömung müssen zu den Normalspannungen und Tangentialspannungen nochdie sogenannten scheinbaren turbulente Spannungen hinzugefügt werden.Reynoldsspannungstensor______ _________ _________⎛2⎞⎜ v' v'xv'y v'xv' ⎟xz⎜⎟_________ ______ _________⎜2v'⎟yv'x v' v'y yv'z⎜⎟⎜ _________ _________ ______ ⎟2⎜ v' v'zv'x zv'y v'z ⎟⎝⎠v’ i Fluktuation oder zufällige SchwankungFür die vollständige Beschreibung des turbulenten viskosen Spannungstensors folgt somit:σσσxyz⎛ ∂vx ⎞= 2μ⎜⎟ − ρ⎝ ∂x⎠⎛ ∂vy ⎞= 2μ⎜ − ρx⎟⎝ ∂ ⎠⎛ ∂vz ⎞= 2μ⎜⎟ − ρ⎝ ∂x⎠______2v'x______2v'x______2v'zτττxyyzzx= τ= τ= τyxzyxz.⎛ ∂vy ∂v⎞x= μ⎜ +x y⎟⎝ ∂ ∂ ⎠⎛ ∂vy ∂v⎞z= μ⎜ +z y⎟⎝ ∂ ∂ ⎠⎛ ∂vz∂vx ⎞= μ⎜+ ⎟⎝ ∂x∂z⎠− ρ v'− ρ v'__________________y− ρ v'x_________xv'v'v'zyzInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 38Neben den Grundgleichungen beinhaltet das vollständige Gleichungssystem zusätzlich die Formulierungder Randbedingungen, die an den Ein- und Ausströmquerschnitten (dynamische Randbedingung;z. B. Hydrograph), der freien Oberfläche und der unteren Berandung festgelegt werdenmüssen.Besondere Aufmerksamkeit benötigt die Nachbildung der freien Oberfläche ist, ist diese doch selbstTeil des Lösungsprozesses. Die Lage der freien Oberfläche wird mit der in der kinematischen Randbedingungintegrierten Kontinuitätsgleichung ermittelt De Strömung an der freien Oberfläche hatauch die dynamischen Randbedingungen zu erfüllen. Hierdurch werden die an der Oberflächewirkenden Windschubspannungen einbezogen. Vielfach wird der Anteil auf Grund der Oberflächenspannungvernachlässigt.DynamischeRandbedingungUntere BerandungFreie OberflächeDynamischeRandbedingungKinematische Randbedingung:Dynamische Randbedingung:sξn sσv AC Dρ Avzs∂ξ r= + v∂txs∂ξ r+ v∂x⎛⎞p n n P p n σ ⎜ v v C z v z v v z v∂ x ∂⎟⎝ y ⎠ ρIndex für Geschwindigkeit und Druck an der freien OberflächeWasserspiegellageNormalvektorOberflächenspannungWindgeschwindigkeitWiderstandbeiwertDicht der Luft2 2r r r ∂ξ r ∂ξ ρAr r r rs s−s=A s− + ( ) ( ) ( )2 x+s 2 D A−zsA−zsys∂ξ∂y( )( )Die untere Berandung ist durch das Stokes’sche Wandhaftungsgesetz festgelegt, dadurch ist dieGeschwindigkeit für eine starre Sohle gleich Null (v B = 0). Für eine bewegliche Sohle wird v B durchdie Bewegung der Sohle v sohle ersetzt, diese entspricht der Geschwindigkeit der Wasserschicht überBoden (v sohle = v).Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 39TURBULENZMODELLIERUNGDas räumliche und zeitliche irreguläre Verhalten ist ein wesentliches Merkmal der Turbulenz, diezusätzlich rotationsbehaftet ist. Die Turbulenz stellt keine eigen Bewegungsform dar und ist auchkeine Materialeigenschaft. Im allgemeinen ist die turbulente Scheinzähigkeit größer als die Zähigkeitμ t (x,t) > μ. Durch eine Energiekaskade erfolgt eine Energiedissipation von großen zu kleinenWirbel, bis diese schließlich aufgelöst werden.Entsprechend der Turbulenzauflösung (Diskretisierung) kann die Turbulenzmodellierung eingeteiltwerden in:• direkte Simulation, direkte Lösung der Grundgleichungen;• LES (Large Eddy Simulation),dabei werden kleine Wirbel vernachlässigt, die Effekte werden jedoch in einem subgridModel berücksichtigt (kleine Wirbel sind isotrop).Die Diskretisierung muss zur Darstellung der anisotropen Wirbelstrukturen entsprechendklein genug gewählt werden.• Reynoldsmittelung (Statistische Turbulenzmodellierung),• Vereinfachung zu 2-D und 1-D Modellen.•Die direkte Simulation und die LES dienen vor allem der Nachbildung der Feinstrukturen der Turbulenz,auch als kohärente Strukturen bezeichnet. Unter einer kohärenten Struktur verstehen Nezuund Nakagawa (1993) ein zusammenhängendes Fluidpaket, in dem gewisse Strömungsgrößen einengewissen räumlichen Zusammenhang über eine gewisse Lebenszeit beibehalten. Die wichtigstenkohärenten Strukturen entstehen in der wandnahen Strömung und werden als ’Bursts’ bezeichnet.STATISTISCHE TURBULENZMODELLIERUNGDer Aufwand der direkten Lösung ist auch heute noch für die Ingenieurpraxis zu hoch. Deshalbwurden Wege gesucht, die sich auf die wesentlichen Merkmale der Strömung konzentrieren. Sowurde schon sehr früh die statistische Betrachtung der Turbulenz eingeführt. Die StrömungsgrößenGeschwindigkeit und Druck werden in eine mittlere Geschwindigkeit bzw. einen gemittelten Druckund entsprechende turbulente Schwankungsgrößen zerlegt (Nezu and Nakagawa 1993), dadurchentstehen die sogenannten Reynolds-Gleichungen.Momentane Fliessgeschwindigkeitv(t) = v m +v’v(t) momentane Geschwindigkeitv m mittlere Geschwindigkeitv’ Fluktuation oder zufälligeSchwankunganalog gilt:p = p m + p‘ρ = ρ m + ρ‘ , wenn die Dichte nichtkonstant istWie bei jeder Mittelung geht Information verloren, diese erscheint aber wieder in der Impulsgleichungals zusätzlicher Spannungsterm (Abbott und Basco 1989). Das Problem der Reynoldsgleichungenist daher, dass diese mehr Unbekannte als Gleichungen enthalten.Die Lösung dieses Problems, oft auch als Schließungsproblem bezeichnet, ist die Bestimmung derturbulenten ViskositätInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 40τ = −ρv'v'ij_________ij⎛⎜∂v= ρνt⎝ ∂xij∂v+∂xIn Analogie zur molekularen Viskosität definierte bereits Boussinesq für die turbulente Strömungeine Scheinzähigkeit, wobei diese, der Natur der Turbulenz entsprechend, keine Stoffeigenschaftist, sondern vom Abflussgeschehen selbst abhängt. Den meist verwendeten Ansatz zur Modellierungder turbulenten Strömung stellen die sogenannten k-ε Modelle dar (Rodi 1993), wobei k dieturbulente kinetische Energie und ε die turbulente Dissipation beschreibt.MODELLIERUNGSANNAHMENBei der vollständigen dreidimensionalen Modellierung tritt das Problem der Ermittlung der freienOberfläche auf. Durch die Annahme einer hydrostatischen Druckverteilung (Flachwasserapproximation)kann die Neigung der freien Oberfläche für die Bestimmung des Druckgradienten herangezogenwerden.Eine Integration der Kontinuitäts- und der Impulsgleichung in vertikaler Richtung liefert die tiefengemitteltenFlachwassergleichungen. Zwischen den 2-D und 3-D Modellen stehen die Mehrschichtmodelledie besonders der Geschwindigkeitsänderung nahe der Sohle Rechnung tragenkönnen.Die zusätzliche Integration über die Breite liefert schlussendlich die klassischen instationären, 1-DAbflussmodelle, die sogenannten Flachwasser- oder Saint-Venant-Gleichungen. Bei 1-D Modellekondensieren die Wirkung der Turbulenz, der Dispersion, der Sohlschubspannung und der Einflussder Sekundärströmungen zu einem einzigen Parameter, dem Energieliniengefälle (Malcherek 2001).ji⎞⎟⎠FlachwasserapproximationEs wird angenommen, dass die vertikalen Schubspannungen einen vernachlässigbaren Einfluss aufdie horizontale Geschwindigkeit haben. Es entsteht ein Gleichungssystem aus drei Gleichungen mitdrei Unbekannten. Aus dieser Annahme folgt einen hydrostatische Druckverteilung über die Tiefe .∂v∂ vx y ∂vz+ + = 0∂ x ∂ y ∂z∂ vx ∂vx ∂vx ∂vx 1 ∂ p 1 ⎛∂σ ∂τx xy ∂τ⎞xz+ vx + vy + vz= X − + ⎜ + + ⎟∂t ∂ x ∂ y ∂z ρ ∂ x ρ ⎝ ∂x ∂y ∂z⎠∂ vy ∂vy ∂vy ∂vy 1 ∂ p 1 ⎛∂τ yx∂σ y∂τyz⎞+ vx + vy + vz= Y − + ⎜ + + ⎟∂t ∂ x ∂ y ∂z ρ ∂ y ρ ⎝ ∂x ∂y ∂z⎠p = p + ρg( ξ + z)oDer Spannungstensor und die Randbedingungen vereinfachen sich ebenso.2D - horizontale - hydrodynamische GleichungenDieses Gleichungssystem entsteht durch eine Mittelbildung über die Wassertiefe, die vertikaleGeschwindigkeitsverteilung wird zur Gänze eliminiert. Die drei Gleichungen beinhalten die dreiUnbekannten v x , v y , ξ. Zum Unterschied zu der Flachwasserapproximation ist die Lage der freienOberfläche als zeitlich veränderliche Größe einbezogen.Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 41∂vx ∂vx ∂vx ∂ξ τbx1 ⎡ ∂ ∂ ⎤+ vx + vy + g + − ( hτxx ) ( hτxy) 0∂t ∂ x ∂ y ∂ x ρh ρh ⎢ + =∂ x ∂ y⎥⎣⎦∂vy ∂vy ∂vy ∂ξ τby1 ⎡ ∂ ∂ ⎤+ vx + vy + g + − ( hτxy ) ( hτyy ) 0∂t ∂ x ∂ y ∂ x ρh ρh ⎢ + =∂ x ∂ y⎥⎣⎦( hvx) ∂ ( hvy) 0∂ξ ∂+ + =∂t ∂ x ∂ yv x ......... Geschwindigkeitskomponente in x-Richtungv y ......... Geschwindigkeitskomponente in y-Richtungξ........... Wasserspiegellageτbi ........ Sohlenschubspannungτii......... Schubspannung= über die Tiefe gemittelte WerteDie Lage der freien Oberfläche ist festgelegt durchξ = zb + hzb ............Lage der Sohlenh..............WassertiefeDie effektive Schubspannung τ xy ist nach Oking (1985) definiert als⎛∂u ∂ v⎞τxy= ρvt+⎜∂y ∂ x ⎟⎝ ⎠Die dynamische Austauschgröße ν t beinhaltet Zähigkeits-, Turbulenz- und Dispersions- Komponenten.Wird zur Darstellung des Zusammenhanges der Sohlschubspannung mit der effektiven Schubspannungein gleichförmiger Abfluss und ein kompaktes Profil angenommen, kann geschrieben werdenals∂ξ τbx1 ∂g + − ( hτxy ) = 0∂ x ρh ρh∂yτ bx Sohlschubspannungτ xy Schubspannung in x-Richtung oder effektive Schubspannungg Erdbeschleunigung∂ghI h∂ ySohlgefälleFür h = const. folgt τbx= ρS+ ( τxy)I SGrundsätzlich ist der Abfluss in einem gegliederten Profil (tiefer Hauptkanal mit anschließendemVorland) ein dreidimensionales Problem. Zwischen den schnell fließenden Hauptstrom und densich langsam bewegenden Abfluss im Vorland erfolgt ein Impuls- und Massenaustausch. DieserAustausch zeigt sich besonders deutlich in der Übergangszone Fluss - Vorland, durch einen Geschwindigkeitsgradientenquer zur Fließrichtung. Die dynamische Austauschgröße wird vielfachdurch ein analytisches Modell bestimmt. Beispielhaft sei hier das Konzept von Karausev (zitiertvon Bogardi 1974) nahm Hanxiang (1985) erwähnt.Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 42ANALYTISCHES MODELLAusgehend von Karausev nahm Hanxiang an, dass der Impulsaustausch in Fließrichtung als auchquer dazu von der Geschwindigkeit abhängig ist.ρgh Karausev:ε = vtMCε.........turbulente Dissipationν t ........turbulente Austauschgröße, kinematic eddy viscosityC........Chezy's RauhigkeitsbeiwertM .......Konstante f(C)M = 0,7C + C 10 65h.........Wassertiefeg.........Erdbeschleunigungv.........Geschwindigkeit f(z)2⎛ Z ⎞v = vo1−P⎜1−⎟⎝ h ⎠2Mvm33 ,und P = oder P = 057 , +2CvoCDie Beziehung zw. der Geschwindigkeit an der Oberfläche v ound der v mmittleren GeschwindigkeitvmC−1ist:= 0.9 .voCv oGeschwindigkeit an der Oberflächev mmittlere GeschwindigkeitFür die über die Wassertiefe gemittelte dynamische Austauschgröße folgt durch Integrationh+aρ gh ρghν = ρ ε = vdz vm DLvmh∫ = =MC MCaD L...... Austauschkoeffizient (Längsrichtung)und in Querrichtungghν = DHvmDH = 48 , ρD H = 4,8 D LMCDer größere Austauschkoeffizient ist eine Folge der zwei Wasserkörper, die sich mit unterschiedlicherGeschwindigkeit bewegen.Bemerkung: Hanxiang änderte die Limits für M.M = 0,6C + C C ≤ 60M = 48 C > 60Aus dem vorhergegangenen folgt eine Geschwindigkeitsverteilung quer zur Fließrichtung in exponentiellerForm.2vf= vfo+ C1 exp( α y)( β )2v expm= vmo+ C2 − y2 232 12C = ( vmo − vfo )⎡( Cm Cf ) ( Mm Mf ) +⎣32 12C2 = C1( Cm Cf) ( Mm Mf)α= ( 1 hf)Mf 24 , Cfβ= 1 h M 24 , C11( )m m m⎤⎦Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 43Die Indizes f und m geben das Strömungsgebiet an, f für Vorland (flood plain) und m für Hauptkanal(main channel).Definition sketchDie berechnete Geschwindigkeitsverteilung ermöglicht die Ermittlung von Abflussparametern, wiez. B. die Schubspannung, die Austauschgröße oder die Interaktionsbreite.Velocity distributionInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 44SAINT-VENANT-GLEICHUNGENDifferentialgleichungen, die den instationären Abfluss in offenen Gerinnen beschreiben, wurdenbereits 1871 von de Saint Venant abgeleitet. Die 1-D Modelle betrachten den Zusammenhangzwischen dem Wasserspiegel und dem Durchfluss in der Hauptrichtung (Abbildung 1), sie setzeneinen schwach gekrümmten Flusslauf und kompakte Querschnitte voraus. Die Stärke liegt in derUntersuchung großer Flussabschnitte. Zur Berücksichtigung von Veränderungen quer zur Fliessrichtung,wie bei gegliederten Profilen (Hochwasserabfluss) wurden gekoppelte Modelle entwickelt.Die Teilquerschnitte sind über die Energielinie oder den Wasserspiegel verbunden, dabei istdas Hauptproblem der turbulente Queraustausch zwischen den parallel verlaufenden Hauptstromund dem Vorlandabfluss.Zusätzliche Voraussetzungen:• inkompressibles Strömungsmedium• Die Druckverteilung längs einer beliebigenVertikalen ist hydrostatisch(Stromfadenkrümmung ist vernachlässigbar)• Über einen Querschnitt ist die Geschwindigkeitkonstant.• kleine Wasserspiegelneigung∂ hsin δ / = / tan δ / = / /∂ x∂Adt∂t∂QQ + ∂ xdxAbbildung 1: KontinuitätsbedingungUnter den beschriebenen Voraussetzungen lautet dann das 1-D Gleichungssystem::Kontinuitätsbedingung∂A∂Q+ = 0 (23)∂t∂xBewegungsgleichung2∂Q∂ ⎛ Q ⎞ ⎛ ∂h⎞+ + gA⎜− Is ⎟ + gAIE= 0t x⎜ βA⎟(24)∂ ∂ ⎝ ⎠ ⎝ ∂x⎠Beschleunigungs-, Trägheits-, Schwerkraft-, ReibungsgliedA FließquerschnittQ Durchflussv mittlere Geschwindigkeith Wassertiefeβ Stossbeiwertg ErdbeschleunigungI s SohlgefälleI E Energieliniengefällemuss die Durchfluss-Strömungswiderstands-Beziehung nicht von vornherein festgelegt werden.Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 45Die einzelnen Terme der Bewegungsgleichung stellen der Reihe nach die Kraftwirkungen infolgeder Beschleunigung, Trägheit, Schwerkraft und Reibung auf ein Fluidpaket dar.Die Gleichung 24 kann auch direkt in Form der Beschleunigung geschrieben und interpretiert werden(Gleichung 25).∂v∂v⎛ ∂h⎞+ v + g⎜− Is ⎟ + gIE= 0(25)∂t∂x⎝ ∂x⎠Die einzelnen Terme beschreiben dann die lokale, zeitliche Änderung der Geschwindigkeit, diekonvektive Beschleunigung, den ungleichförmigen Abfluss und wieder den Reibungsanteil alsEnergiehöhenverlust.Durch die Einführung des sogenannten "Hydraulischen Leitwertes" (Conveyance) Q = K I E mussdie Durchfluss - Strömungswiderstands-Beziehung nicht von vornherein festgelegt werden.12E23Z. B. ist für die Formel von Manning/Strickler Q = k R I AK = k St R 2/3 A. (26)k st ....... StricklerbeiwertR ........ hydraulischer RadiusDie vollständige dynamische Welle (Gleichung 24 und 25) kann durch weitere Vereinfachungen, ineine quasi stationäre Welle, eine diffusive Welle, kinematische Welle bis hin zu einem stationären,gleichförmigen Abfluss übergeführt werden. Dies bildet seit jeher die klassischen Ansätze derHydromechanik.StFließgesetze für den stationären gleichförmigen AbflussManning/Strickler:23 12Q = KsTR I AK = k sT R 2/3 AChezy:12Q=A C( R I)K = A C R 1/2Darcy/Weisbach:Q = ⎛ 12g R I A⎝ ⎜ 1 ⎞ 8 ⎟λ ⎠K = ⎛ 12g R A⎝ ⎜ 1 ⎞ 8 ⎟λ ⎠Die Fließgesetze können miteinander über die Widerstandskoeffizienten in Beziehung gesetztwerden.23 12 8g ksTR = C R = ⎛ R⎝ ⎜ ⎞⎟λ ⎠1216 goder C = ksTR = ⎛ ⎝ ⎜ 8 ⎞⎟λ ⎠12Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 46ZUSAMMENSTELLUNG DER MODELLANWENDUNGENAus Tabelle 1 ist mit einem Blick das Anwendungsgebiet der einzelnen Ansätze ersichtlich. Fetthervorgehoben sind die zur Zeit gebräuchlichen Einsatzgebiete in der Fliessgewässermodellierung.Tabelle 1: Modelleinsatzbereiche nach Malcherek (2001) und Nezu & Nakagawa (1993)Raumdimens.Diff.glei-chungenAuflösung AnwendungDNSDirect numerical Simulation3 4 1 mmLES(Large-eddy Simulation)3 4 1 cmRANSReynolds averaged NS 3 4 1 dm3D hydrostatischeSimulation 3 3 1 mTiefenintegriert2-D Simulation 2 3 10 mFliessgewässermodellierung,Morphologie, gegliederte QuerschnitteKlassische instationäre FliessgewässermodellierungDe Saint VenantQuerschnittsintegrierte1-D Simulation 1 2 100 mFeinstrukturen der Turbulenz,kohärente Strukturen, Starkidealisiert, kleine Re-Zahlen,Maßstabspektrum engGroßskalige Turbulenzbewegung,kohärente Strukturenk-ε Modelle, Statistische TurbulenzmodellierungFlachwasserapproximation.KüstengewässerInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 47LÖSUNGSANSÄTZE DER 1-D INSTATIONÄREN GRUNDGLEICHUNGENDer Ansatz des vollständigen Gleichungssystems erfordert, insbesondere bei großen Flusssystemeneinen hohen Rechenaufwand. Für einige Anwendungsfälle ist eine Modellreduktion möglich, wobeidie Näherungslösung dem Ergebnis der Lösung der Bewegungsgleichung bei ungekürztem Ansatzentsprechen soll.Schwere WelleDie Wellengeschwindigkeit beträgt in der einen Richtung w v ( gh)Richtung w v ( gh)212112= + und in der anderen= − . Bei FROUDE-Zahlen Fr < 1 liegt strömender Abfluss vor, die Störungpflanzt sich stromaufwärts fort. Ist Fr > 1 bewegt sich die Störung in Fließrichtung. Werden dieReibungskräfte vernachlässigt, so entsteht eine Wasserwelle, die keinerlei Dämpfung unterliegt.Das Fehlen der Reibungskräfte bewirkt, dass sich Wasserwellen bzw. Störungen ungedämpft fortpflanzen.Kinematische WelleBei diesem Ansatz werden die Trägheitsterme (∂v/∂t und ∂v/∂x) und das Druckglied (∂h/∂x) vernachlässigt.Der Abfluss Q hängt nur von h und nicht von der Zeit oder von der Ableitung von hnach x ab. Der Durchfluss ist eine eindeutige Funktion der Wasserspiegellage.Die Schnelligkeit (celerity) w der Welle ist von der FROUDE-Zahl unabhängig und ergibt sich zuw = 1,5v (v = Fließgeschwindigkeit). Auch bei dieser Wellenform unterliegt eine Störung keinerDämpfung.Die schwere und die kinematische Welle sind Grenzfälle der in natürlichen Flussläufen beobachtetenPhänomene.DiffusionswelleWerden nur die Trägheitsterme (∂v/∂t und ∂v/∂x) in der Bewegungsgleichung unterdrückt, erhältman mit der Kontinuitätsbedingung eine Beziehung, die einer Diffusionsgleichung entspricht.Diffusionswellen pflanzen sich wie kinematische Wellen in der Fließrichtung fort. Ein Unterschiedbesteht jedoch im Verhalten. Diffusionswellen flachen mit der Zeit ab. Die Stärke der Dämpfung isteine Funktion der Wellenlänge.Dynamische Wellestellt die Lösung des vollständigen Gleichungssystems dar. Die Einbeziehung der Trägheitskräftebewirkt, dass sich die Welle entlang zweier Wege fortpflanzen kann. Bei Fr < 1 wandert die Primärwelleflussabwärts, die Sekundärwelle flussaufwärts. Ist Fr > 1 bewegen sich beide Wellen ingleicher Richtung. Das Ausmaß der Abflachung ist bei einer dynamischen Welle von derFROUDE-Zahl und der Länge der Welle abhängig.BERECHNUNGSMETHODENDas Verfahren der Charakteristiken beruht auf einer Umformung der beiden partiellen Differentialgleichungennach DE SAINT VENANT in ein System von vier gewöhnlichen Differentialgleichungen.Das so erhaltene Gleichungssystem bezeichnet man als die sogenannten Gleichungen derCharakteristik. Die Integration erfolgt entlang dieser Charakteristiken (natürliche Methode) oder ineinem festem Gitter (lokale Methode).Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 48METHODE DER FINITEN DIFFERENZEN (FDM)Bei direkten Differenzenverfahren werde die Differentialquotienten der DE-SAINT-VENANT-Gleichungen zu Differenzenquotienten umgeformt und die Gleichungen unter Berücksichtigung derAnfangs- und Randbedingungen gelöst.Einführung in die Methode der finiten Differenzen anhand einer einfachen Transportgleichung.Die Differentialgleichung∂h∂h+ c = 0∂t∂xbeschreibt den Transport eines Mediums mit der Geschwindigkeit c. Diese Gleichung gibt eineAussage über einen bestimmten Punkt der Lösungsebene (eine infinitesimale kleine Region). DurchEinführung von finiten Differenzen erweitert sich die Gültigkeit auf eine finite Region.Darstellung der Lösungsebene (nach Abbott 1989)Projektion und SchnittDie Ableitung ∂h/∂x im Punkt P entspricht der Neigung der Tangente TT.Näherungsmöglichkeiten:Neigung der Sehne AP ("backward difference"):⎛ ∂h⎞⎜ ⎟⎝ ∂x⎠njh≈nj− hΔxnj − 1bzw.⎛ ∂h⎞⎜ ⎟⎝ ∂t⎠njh≈nj− hΔtn−1jSteigung der Sehne PB ("forward difference"):⎛ ∂h⎞⎜ ⎟⎝ ∂x⎠njh≈nj+1− hΔxnjbzw.⎛ ∂h⎞⎜ ⎟⎝ ∂t⎠nj≈hn +1j− hΔtnjSteigung der Sehne AB ("central difference"):⎛ ∂h⎞⎜ ⎟⎝ ∂x⎠njh≈n nj+ 1 − hj−12Δxbzw.⎛ ∂h⎞⎜ ⎟⎝ ∂t⎠nj≈hn + 1 n −1j − hj2ΔtInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 49Weitere Näherungen können durch gewichtete Kombinationen dieser drei Möglichkeiten gebildetwerden, wobei Werte ungleich j und Zeitabschnitte ungleich n Δt einzuführen sind.Die Transportgleichung kann z. B. durch eine "backward difference" in der x-Richtung und durcheine "forward difference" in t-Richtung genähert werden:oderhhn + 1jn + 1j− hΔtEs bedeuten:j..........Index Berechnungspunktn .........Index Zeitschrittnj+ c hnj− h Δxnj−1≈ 0⎛ c Δt⎞x h n c t≈⎜− ⎟⎝ ⎠j + ⎛ ⎝ ⎜ Δ ⎞1⎟Δ Δx⎠hnj−1Der Funktionswert n+1 eines Zeitabschnittes kann aus den zeitlich vorhergehenden Werten berechnetwerden.Explizite DifferenzenverfahrenDie verschiedenen expliziten Verfahren unterscheiden sich nach dem Schema des Lösungsansatzes,das für die Approximation der Orts- und Zeitableitungen gewählt wird. In jedem Fall ist bei explizitenVerfahren die numerische Stabilität zu untersuchen. Es ist zu überprüfen, ob die Lösungen derDifferenzengleichungen eine bestimmte Grenzabweichung von denjenigen der Differentialgleichungennicht überschreiten. Wachsen diese Abweichungen (Truncation Error) zwischen den Differenzen-und den Differentiallösungen unbegrenzt an, so versagt die gewählte Lösungsmethode.Implizite DifferenzenmethodenDie Impliziten Differenzenmethoden unterscheiden sich von den expliziten Verfahren dadurch, dasssowohl in der Approximation der örtlichen und zeitlichen Ableitungen bekannte Zustände an derStelle j zur Zeit t als auch unbekannte Zustände an den Stellen ...,j-1, j, j+1,... zur Zeit t+1 herangezogenwerden. Bei den impliziten Verfahren ist die numerische Stabilität von vornherein gegeben;zu überprüfen ist die Konsistenz. Grundsätzlich darf bei impliziten Verfahren der Zeitschritt beliebiggewählt werden. Haben die Trägheitskräfte einen starken Einfluss auf die Strömung, so führtdie Wahl von großen Zeitschritten bei impliziten Schemata zu erheblichen Fehlern. Sind die Trägheitskräfteklein, jedoch nicht vernachlässigbar, so können die Zeitschritte bedeutend größer gewähltwerden. Dadurch liegen die erforderlichen Rechenzeiten erheblich unter denen der explizitenVerfahren.Ein weiterer Vorteil der impliziten Verfahren besteht in der vom Ortsschritt unabhängigen Wahldes Zeitschrittes. Dies erlaubt die Einführung von verschieden großen Ortsschritten, wodurch eineAnpassung der Dichte der Berechnungspunkte an die unterschiedlichen Verhältnisse ermöglichtwird.Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 50WASSERSPIEGELLAGEBERECHNUNG⎛ ∂gAh I g A Q Qs∂ x+ ⎞ / /⎜ ⎟ + =⎝ ⎠ K20mit K = kST R2/3 A∂h∂xQQ+ Is+ =k R A02 4 3 2STNäherung: breites Gerinne R ≈ hdh Q⇒ = −Idx2 2 10 3k b hST2sKontinuumformExakte Lösungz. B. BresseDruckliniengefälle = WasserspiegelgefälleNäherung durch finite Differenzendhh j+ 1 ≈ h j + Δ xdxdh h j+1− h j≈dx Δxh j+ − h21 j Q≈Δxk b hST2 10 3j−IsQhj+ ≈ hj+ ⎛ 2x⎜Is⎝ k b h−⎞1 Δ ⎟2 2 10 3Euler Methode (Diskrete Form)⎠STjInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 51Eine Verbesserung der Näherung kann erreicht werden, wenn die Steigung in einem Punkt j + 1/2herangezogen wird.⎛ dh ⎞⎜ ⎟⎝123dx ⎠h j+1−hjΔx⎪⎧= ⎨⎪⎩ b2kQ22STh10 3⎪⎫− Is⎬⎪⎭j+1 2⎧ 2⎪ Q ⎫⎪hj+1 = hj+ Δ x⎨2 2 10 3−Is⎬⎩⎪ b kSTh ⎭⎪j+12Vorerst sind nicht alle Werte zur Berechnung in j + 1/2 vorhanden, deshalb wird in zwei Schrittenvorgegangen (Predictor - Corrector - Method)1.) Euler-Methode⎧ 2∗ ⎪ Q ⎫⎪hj+1 = hj+ Δ x⎨2 2 10 3−Is⎬⎩⎪ b kSTh ⎭⎪jhj+1/ 2=hj+ h2*j + 12.) Korrekturhj+1= hj⎪⎧+ Δx⎨⎪⎩ b2kQ22STh10 3⎪⎫− Is⎬⎪⎭j+1/ 2Schritt 1 und 2 werden als verbesserte Euler-Methode bezeichnet.Diskretisierung in h:im Bereich großer Krümmungen der Oberfläche sind kleinere Schritte Δx von Vorteil → InversesSchema.Δhxj+1 = xj+⎛ 2Q ⎞⎜− Is⎟2 2 10 3⎝ k b h ⎠STj+12Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 52Beispiel BehälterausgleichBehälter ABehälterausgleichBehälter BAnnahmen:Keine Speicherung im RohrVerlustbeiwert ξ = 2,0g = 10m/s25 m10.000 m 21 m 2ΔH2 m5.000 m 2ΔHΔHQ=ξ2gAQ=gA2Q ... Durchfluss222A.... RohrquerschnittH.... Energiehöheξ.... Verlustbeiwert= Δ und Q f( t z)1) Q A g H= , z = ΔH2) dV = Qdt = A gz dt und dV = AA dx = AB dy KontinuitätdV⇒ dx = =AAdV⇒ dy = =ABQAAQAB3) z = hA − hB − x − ydtdt4) dz = −(dx + dy)⎛ 1 1 ⎞5) dz =− ⎜ + ⎟⎝ A A ⎠Qdt⎛ A=−⎜⎝ ABA+ AAAOBAB⎞⎟ A g⎠z dt1 dz 1t =−K∫ =− 2z K zΔHOΔH( AAAB)( + )2 2t = ΔH=Δ HK A A A g=B A pipe( )dt=− 1Kdzz210000 . . 5000= ≈10000 + 5000 1 10 3 3650 s 1 h 1 min.Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 531) Q = A g Δ H2) ΔH = h A − h B3) ΔV = Q Δt (dv = Q dt)4) ΔV = A A Δh AnΔhA = hn + 1 n + 1 n−h ⇒ h = h −Δh5)hn+1A= hnAΔV−AA= hnAQΔt−A6) A Δh = A Δh ⇒ Δh = ΔhAAAA A B B B AAB1B B n +1B n7) hB n + = hB n + Δ hBΔh = h −hn + 128) h A ( AnAn + 1= h + h ) 2⎫⎪⎬ ⇒ Δ Hn + 129) h B ( hBn hBn + 1= + ) 2⎭⎪Δt = 200sTimesteph A [m] h B [m]Q n[m3/s]h An+1[m]h Bn+1[m]AAh An+1/2[m]Ah Bn+1/2[m]AQ n+1/2[m3/s]Ah An+1[m]0 5,0 2,00 5,48 4,89 2,22 4,95 2,11 5,33 4,89 2,211 4,89 2,21 4,18 4,79 2,42 4,84 2,32 5,02 4,79 2,412 4,79 2,41 4,88 4,69 2,605 4,74 2,51 4,72 4,69 2,603 4,69 2,60 4,57 4,60 2,78 4,65 2,69 4,43 4,60 2,784 4,60 2,78 4,27 4,51 2,95 4,56 2,87 4,11 4,52 2,945 4,52 2,94 3,97 4,44 3,10 4,48 3,02 3,82 4,44 3,096 4,44 3,09 3,67 4,37 3,24 4,41, 3,17 3,51 4,37 3,237 4,37 3,23 3,38 4,30 3,37 4,34 3,30 3,22 4,31 3,368 4,31 3,36 3,08 4,25 3,48 4,28 3,42 2,93 4,25 3,489 4,25 3,48 2,77 4,19 3,59 4,22 3,54 2,61 4,20 3,5810 4,20 3,58 2,49 4,15 3,68 4,175 3,63 2,33 4,15 3,6811 4,15 3,68 2,18 4,11 3,76 4,13 3,72 2,02 4,11 3,7612 4,11 3,76 1,87 4,07 3,84 4,09 3,80 1,70 4,08 3,8313 4,08 3,83 1,58 4,05 3,89 4,065 3,86 1,43 4,05 3,8914 4,05 3,89 1,26 4,02 3,94 4,035 3,915 1,10 4,03 3,9315 4,03 3,93 1,0 4,01 3,97 4,02 3,95 0,89 4,01 3,9716 4,01 3,97 0,63 4,00 3,99 4,005 3,98 0,50 4,005 3,9817 4,005 3,98 0,50 3,095 4,00 4,00 3,99 0,316 4,00 3,99h Bn+1[m]18 4,00 3,99 0,316 3,995 4,00 ~4,0 3,995 0,22 3,996 3,998_______Σ 3600sInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 54HOCHWASSERWELLE (FLOOD ROUTING)Die grundlegende mathematische Beschreibung der Hochwasserwellenberechnung ist die Erhaltungder Masse.∂Q∂+ bh =0∂x∂tDer Durchfluss in jedem Querschnitt wird als eindeutige Funktion des Wasserstandes angenommen.Q = Q(h)∂QdQ=∂xdh∂h∂xEingesetzt in die Kontinuitätsgleichung folgtdQ ∂h∂+ bh = 0dh ∂x∂t∂hdQ ∂∂tdh b h+ ⎛ ⎝ ⎜ ⎞⎟ = 0⎠ ∂x∂h∂t∂h+ c∗ = 0 c∂x∗=dQdhbAnmerkung:Die „Fortpflanzungsgeschwindigkeit c*“ ist nicht zu verwechseln mit der Fortpflanzungsgeschwindigkeitc = gh (celerity)Beispiel: Flood routing1.) Differentialgleichungh ............. Wassertiefe∂h∂h+ c∗ = 0 c* ........... Fortpflanzungsgeschwindigkeit∂t∂xQ ............ DurchflussdQc dhb = 1, 5 m s b ............. SpeicherbreiteInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 552.) Diskretisierung:Vorwärts in der Zeit, rückwärts in x.hn + 1jΔt− hnj+ c∗hnj− h Δxnj−1= 0hn + 1jn n( hjhj−)n c t− hj+ ∗ Δ− =Δx1 0mitrc t= ∗ΔΔxr... Courent Zahln + 1jnjn n( j j−1)h = h −r h −h3.) Finites Differenzen Schema:n + 1jnj−1 1( )h = r h + −r hnjDifferenzenoperatorh-Werte für folgende Rastereinteilung sind zu berechnen:a) Δx = 3.000 m Δt = 1.000s => r = 1,5 1 . 0003.000= 0,5b) Δx = 3.000 m Δt = 2.000s => r = 1,5 2 . 0003.000= 1,0c) Δx = 3.000 m Δt = 3.000s => r = 1,53.0003.000= 1,5Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 56Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 57Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 58Definition of VariablesINPUT OF DATA:TIMESTEPS(M), STEP Δx(N),RIVERLENGTH(RL)TIMERANGE,TIMESTEPS-OUTPUTTABLE:TIME BOUNDARY - DEPTHTABLE:C-VALUES - DEPTHINITALIZE:DX = RL/NDT = TIME/MRR(0) = 0H(0,0) = HT(1)TSTEP(0) = HT(1)RMAX = 0RMIN = 1BOUNDARY-VALUES RIVERRIDGEDO I=1,NRR(I) = RR(I-1)+DXH(0,I) = HT(1)BOUNDARY-VALUES TIMEDO I = 1,MTSTEP(I) = TSTEP (I-1) +DTCALL SUBROUTINE ITERH(I;0) = HHCALCULATE C-VALUES AND WATERDEPTHDO I = 0,M-1do J = 1,NCALL SUBROUTINE ITERR = CC * DT/DXCALCULATE RMAX AND RMINH(I+1,J) = (1-R)*H(I,J)+R*H(I,J-1)OUTPUT: RMAX, RMINFLOODROUTING:DISTANCE AND TIMETABLEENDInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 59SUBROUTINE ITERDEFINITION OF VARIABLESINTERPOLATE (F) -> (FI)DO J = 1,NNF>TABV(J) AND F

HYDRODYNAMIK SEITE 60PROGRAM FLOODC DATE 88-12-08C COMPUTATION OF A FLOOD WAVE ALONG A CHANNELC PARAMETER LIST:C C C ACCORDING TO HCC HC DEPTH FOR C-VALUESC CC INTERPOLATED C-VALUESC DX LENGTHSTEPSC DT TIMESTEPCH(I,J) DEPTH IN EACH GRIDPOINTC HH INTERPOLATED DEPTH - BOUNDARYC HT DEPTH - BOUNDARYC TH TIME - BOUNDARYC IT OUTPUT FOR SELECTED TIMESTEPSC K NUMBER OF VALUES FOR BOUNDARYC L NUMBER OF C-VALUES FOR INTERPOLATIONC M NUMBER OF TIMESTEPSC N NUMBER OF DXC R R=CC*DT/DXC RL RIVERLENGTHCRMAX MAXIMUM R - VALUECRMIN MINIMUM R - VALUEC RR RIVERRIDGECTIME TIMERANGE CONSIDEREDCTSTEP TIMESTEPSCPARAMETER (IMAX=100,ITAB=50)REALHH,H(0:IMAX,0:ITAB),C(ITAB),HC(ITAB),HT(ITAB),TH(ITAB),RLREAL CC,DX,DT,TIME,R,RR(0:IMAX),TSTEP(0:IMAX),RMAX,RMINCHARACTER LINE*6INTEGER I,J,K,L,N,M,ITOPEN(10,FILE='FLODAT')OPEN(11,FILE='FLORES')C -------------------------- READ DATA ------------------READ(10,*) N,M,RL,TIME,ITREAD(10,*) KREAD(10,*) (TH(I),I=1,K)READ(10,*) (HT(I),I=1,K)READ(10,*) LREAD(10,*) ( C(I),I=1,L)READ(10,*) (HC(I),I=1,L)C -------------------------- INITIALIZE -----------------LINE = '------'DX = RL/NDT = TIME/MRR(0) = 0.H(0,0)= HT(1)TSTEP(0) = 0.RMAX = 0.RMIN = 1.C -------------------------- BOUNDARY VALUES RIVERRIDGE ---DO 20 I=1,NRR(I)=RR(I-1)+DXH(0,I)=HT(1)Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 6120 CONTINUEC -------------------------- BOUNDRY VALUES TIME (LEFT BOUNDARY)--DO 40 I=1,MTSTEP(I)=TSTEP(I-1)+DTCALL ITER (TSTEP(I),HH,K,TH,HT,ITAB)H(I,0)=HH40 CONTINUEC -------------------------- CALCULATE C AND WATERDEPTH ----------DO 60 I=0,M-1DO 60 J=1,NCALL ITER (H(I,J),CC,L,HC,C,ITAB)R=CC*DT/DXIF(R .GT. RMAX) RMAX=RIF(R .LT. RMIN) RMIN=RH(I+1,J)=(1-R)*H(I,J)+R*H(I,J-1)60 CONTINUEC ------------------------- OUTPUT -------------------------------WRITE(11,99993)WRITE(11,99994) RMAX,RMINWRITE(11,99995) (LINE,I=1,N+1)WRITE(11,99996)WRITE(11,99995) (LINE,I=1,N+1)WRITE(11,99997) (RR(I),I=0,N)DO 80 I=0,M,ITWRITE(11,99998) TSTEP(I),(H(I,J),J=0,N)80 CONTINUE99993 FORMAT(' FLOODROUTING - DISTANCE AND TIMETABLE',/)99994 FORMAT(' RMAX =',F4.2,' RMIN =',F4.2,' STABLE FOR 0

HYDRODYNAMIK SEITE 62Beispiel: Abfluss in einem Gerinne, mit kompaktem ProfilDazu sollen die vollständigen Gleichungen nach de Saint Venant herangezogen werden.Kontinuitätsgleichung:∂Q∂+ bh = 0 (1)∂x∂tBewegungsgleichung:∂Q∂ ⎛ 2Q ⎞ ∂+ ⎜βgAh I∂t∂x⎝A ⎠⎟ + ⎛⎜⎝ ∂x+s⎞ gQQ⎟ +⎠ 2= 0 (2)C ARQ........ Durchflussh ........ Wassertiefeb ........ SpeicherbreiteA........ FließquerschnittR........ hydraulischer RadiusI s ........ SohlgefälleC........ Widerstandsbeiwert nach Chezyß ........ Boussenesq-KoeffizientDie Lösung des Problems erfolgt mit Hilfe eines impliziten Differenzschemas. Im wesentlichen gibtes dafür zwei Vorgangsweisen:1. In jedem Gitterpunkt ist sowohl h und Q als unbekannte Größen definiert.2. Die unbekannten Größen werden abwechselnd in den Gitterpunkten definiert (staggered-gridscheme).Falls eine Flussstrecke ( 0 < x≤)L in n Teilstrecken zerlegt wird so hat man 2n+2 Unbekannte zurZeit t + Δt (z. B. n+1 Wasserstände, n+1 Geschwindigkeiten).Anfangsbedingungen:Randbedingungen:Zur Zeit t = 0 in jedem Punkt der Strecke bekanntIm Punkt x = 0 ist z. B. der Wasserstand oder der Durchfluss alsFunktion der Zeit gegeben.Im Punkt x = L ist z. B. eine Schlüsselkurve gegebenInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 63ad 1: Preismann Schema:ad 2: Abbott-Ionesco Schema:staggered gridInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 641) h, Q in jedem Gitterpunkt (Preismann-Verfahren).In aufeinanderfolgenden Gitterpunkten in Raum und Zeit wird sowohl die Kontinuitäts- als auch dieBewegungsgleichung angesetzt.Die allgemeine Form der Gleichungen (Ableitung Anhang 3) lautet:n+1 n+1 n+1 n+1j + B1 jh+ C1 1 Q j + 1 + D1 j h j + 1 = E1jA1 QA2jjn+1 n+1 n+1 n+1j Q j + B2 jh+ C2 1 Q j + 1 + D2 j h j + 1 = E2jj(3)(4)wobeiA1j = −Θ/ΔxB1j = b s (1−Ψ)/ΔtC1j = Θ/ΔxD1j = b s Ψ/Δt1 j =Θ −1ΔxnQ j 1 − Q( ) b( )s+ [( 1 − Ψ)h + h ]n n nE + jj Ψ j + 1A2B2C2D2E2jjjjj=( 1 − Ψ)Δt= −ΘgAΨ= +Δtβ−ΔxAn+1 2j+1 2Qx nβΔ= +ΘgA=( 1 − Ψ)Δtn+1 2j+1 2QΔxnj + 1n + 1 2j+1njΔxnQ jn+1 2jΨ+ QΔtΔtg+2 Cg Q+2 CQnj2 n+1 2( AR) nj + 12 n + 1 2( AR) j+1jn+1 2gAn j+1 2j + 1 −+ +Δxn n n+1 2( 1 − Θ)( h j 1 − h j ) gA IsDie Lösung der Gleichungsmatrix erfolgt durch einen vorwärts und einen rückwärts Eliminations-Schritt (double sweep algorithm)j+1 2Dazu werden folgende Beziehungen eingeführt (Anhang 4):n + 1j j n + 1j jQ = F h + G(5)n + 1jh = H Q + I h Jj n j+ + 1 1 j n j+ +1 1 j(6)Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 65wobei (recurrence relations):−C1jH j =A1 F + B1IJjjj j j−D1j=A1 F + B1j j j= − E1 − A1GA1 F + B1j j jj j j( A2jFj Ij D2j( A2jFj + B2j)Hj + C2j2 − 2 − ( 2 + 2 )( A2jFj + B2j)Hj + C2jFj+ = − + −1Gj+ 1 =E A F A F B Jj j j j j j jist.Weiters sind für den strömenden Abfluss Randbedingungen an beiden Modellenden notwendig.Linke Randbedingung:QB bekannt: FB = 0 und GB = QB B ....... Boundaryoder hB bekannt:FB = ∞ und GB = −FB hBzur Berechnung FB >> angenommenoder mit 1. Gitterpunkt Berechnungbegonnen.3a) A2BA1B QB + A2BB1B hB + A2BC1B Q1 + A2BD1B h1 = E1B A2B4a) A1BA2B QB + A1BB2B hB + A1BC2B Q1 + A1BD2B h1 = E2B A1BE13a) - 4a)BA2B− E2BA1BA2BB1 B − A1BB2B⇒ Q1 =−h BA214BC1444444 B − A1BC2BA244 24BC144444444B − A1BC2B3G1A2BD1B− A1BD2B−h1A214444BC1B− A124444BC2B3F1QB wird dann bei der Rückwärtselimination berechnet.Rechte Randbedingungen:Q jj bekannt: hjjQ jj G= −F Foder h jj bekannt: Q jj = F jj h jj + G jjjjjjjjInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 66Zusammenfassung - LösungsalgorithmusDie Anfangsbedingungen (Durchfluss, Wassertiefe) zum Zeitpunkt n sind bekannt.Weiters sind die Randbedingungen (links und rechts) zum Zeitpunkt n+1 gegeben.1.) Randbedingungen:Q 1 gegeben => F1 = 0 und G1 = QBh 1 gegeben => mit F2 und G2 begonnen (siehe linke Randbedingung)2.) Eliminationsschritt von links nach rechts.3.) Rechte Randbedingung:hjj bekannt => Qjj berechnetQjj bekannt => hjj berechnet4.) Rückwärtseliminationsschritt.Berechnung der fehlenden Werte des linken Randes5.) IterationsschrittBerechnung neuer Koeffizienten.6.) Neuer Eliminationsschritt7.) Zeitschritt beendet; berechnete Werte sind Anfangswerte für nächsten Zeitschritt.Anmerkung: In der Regel sind 2 Iterationsschritte ausreichend.Die Iteration ist notwendig da die Koeffizienten auch vom Zeitpunkt n+1 abhängig sind.Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 67Beispiele für Nichtlinearitäten:TrägheitsgliedA = f(h)n + 12j+1nj+1n + 12j1 Schritt: A = A und A = An + 12 1 n n+1Iteration: j+1 ( )Reibungsglied:A,R,C = f (h)A = A + A j+12A n + 12 1j ( A n A n + 1= + ) j22 n + 12 21 Schritt: ( C AR) = ( C AR)j2 n + 12 2( C AR) = ( C AR)j+12n 1222 Iteration: ( ) = ( )njnj+1C AR C A Rnj+ n + 12 n + 12 n + 12jjj j2n + 122( ) = ( )C AR C A Rj+1n + 12j+1n + 12j+1n + 12j+1Schwerkraftglied: A = f (h)1 Schritt:1A = ( Aj+ Aj+1 )21 nA = Aj+12+ A+ 12 n + 122 Iteration: ( j )nInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 68TE = O(Δt,Δx)Δtr = uΔxn+1Δtnj-1 jC.P.Explizite Schemata:rückwärts in xvorwärts in tr = 1 ....................... genaue Lösungr ≤ 1 .......................stabiln+1nn-1nn-1nj-1n+1j-1nC.P.Δxj j+1C.P.C.P.jj j+1j C.P. j+1vorwärts in xrückwärts in tr = 1 ....................... genaue Lösungr ≥ 1 .......................stabilrückwärts in x und tr = -1......................genaue Lösungr ≥ 0 .......................stabilcelerity +/- r ≥ -1 ...stabilvorwärts in x und tr = 1 ...................... genaue Lösungr ≥ -1 ......................instabilr ≤ -1 ......................Dämpfung der Lösungzentriert in xvorwärts in tr = 1/6 ...................beste Lösungr ≤ 1/2 ...................stabilTE = O(Δt,Δx 2 )Δtr = D⋅ Δx2Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 69Implizite Schematan+1PreismannΨΔxΘΔtRundungsfehler = 0stabilψ = 0,5 ® Θ < 0,5 instabilΘ = 0,5 genaueste LösungΘ > 0,5 dissipativn+1n j j+1→≤Ψ≥Abbott - Ionescostabil für Θ = 0,5nj-1j C.P. j+1ΘΔtImplizitVorteile:• Durch bessere Zentrierung in Raum und Zeitnormalerweise größere Genauigkeit• Stabilität• größere Zeitschritte Δt2 2( , )TE= O Δt ΔxNachteile:• Zwei oder mehrere Unbekannte im folgendenZeitschritt sind in denLösungsalgorithmus einbezogen• Abhängig von der Reihenfolge und der Richtungder Berechnung• Benötigt gut definierte Randbedingungen(obere und untere)ExplizitVorteile:• Einfacher Lösungsalgorithmus, neuer Wertdirekt aus vorhergehenden Zeitschritt berechenbar• schnell• Unabhängig von der Reihenfolge und derRichtung der Berechnung• Für bestimmte Fälle genaue Lösung reproduzierbarNachteile:• bedingt stabil• eingeschränkt durch die Courent-Friedrich-Lewy BedingungTE= OΔt,Δx( )Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 70Fließende Retention (Channel Routing)RETENTIONHierbei ist stets eine zeitliche Verschiebung der Flutmaxima und eine Deformation der Hochwasserwellezwischen zwei Stellen im Gewässerlauf zu beobachten, die abhängig sind von dem imbetrachteten Gewässerabschnitt vorhandenen Retentionsraum und den instationären Beschleunigungsverhältnissendes Abflusses – charakterisiert durch die Fließgeschwindigkeit v und dem Wasserstandh als Funktion von Ort und Zeit.Berechnung des Wellenablaufes (Flood Routing)Ziel: Aus der Ganglinie h 1 (t) bei Pegel P 1 eine Prognose für die Ganglinie bei Pegel P 2 zu erstellen– sind Daten der Gerinnegeometrie, Randbedingungen h 1 (t) und Anfangsbedingungen h (x, t=0)notwendig. Eine Gliederung der Verfahren ist bei BRETSCHNEIDER (1982) zu finden, wobeidiese grundsätzlich ina) Verfahren der Charakteristiken• Natürliche Methoden• Lokale Methodenb) Direkte Differenzenverfahren• Explizite Methoden• Implizite Methodeneingeteilt werden können.Die Integration erfolgt entlang dieser Charakteristiken (natürliche Methode) oder in einem festenGitter (lokale Methode). Dieses Verfahren kann als weitestgehend exakt gelten, weil der instationäreWellenablauf dabei lediglich unter Vernachlässigung des konvektiven Beschleunigungsanteilesin der Vertikalrichtung , ansonsten aber unter Erfassung sämtlicher beteiligter Einflüsse beschriebenwird (näheres dazu bei SCHRÖDER 1972; RADLER 1989).Stehende Retention (Reservoir Routing)Künstliche Retentionsräume werden durch Rückhaltebecken charakterisiert. Zu den natürlichenRetentionsräumen zählen• Seen• natürliche Teichanlagen• ausgedehnte Sumpf- und MoorlandschaftenDer Einfluss von Sümpfen oder Hochmooren auf die Abflussdämpfung ist von der Einzugsgebietsflächeabhängig und zudem vom Grad der anthropogenen Beeinflussung. Dies ist sinnvoll für sehrdetaillierte Modelle und für Einzugsgebiete, in denen ein signifikanter Teil der Gesamtoberflächeaus Seen besteht.Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 71Als Annahme für die Berechnung von Retentionsräumen gilt, dass die Fließgeschwindigkeit instehenden Gewässern vernachlässigbar ist, d. h. der Wasserspiegel ist annähernd horizontal. Durchdiese Annahme reduzieren sich die für die instationäre Strömung maßgebenden hydrodynamischenGleichungen aufdie Kontinuitätsgleichung unddie Ausflusscharakteristik.Das Retentionsverhalten R folgt aus der Differenz Qzu - Qab und ist gleich der zeitlichen Änderungdes Beckeninhaltes V.Q = Q − Q =R zu abdVdtAls geometrische Grundlagen genügen das Becken- bzw. Seeinhaltsdiagramm (Beziehungzwischen Wasserstand und Speichervolumen).Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 72Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 73Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 74Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 75NUMERISCHE LÖSUNG DER RETENTIONSGLEICHUNG( Q Q ) −12 ( Q + Q ) = ( V − V ) Δt1 2 zui zui 1 abi abi 1 i 1 i+ ++ +Wird die Gleichung entsprechend den unbekannten Gliedern zum Zeitpunkt ti+1 = ti+Δt aufgelöst,so folgt( )12 Q + Q −Q − V Δt = 12Q + V Δ tzui zui+ 1 abi i abi+ 1 i+1mit V = V(Qab) lässt sich die Gleichung numerisch auflösen.Beispiel: Beckenentleerungscharakteristik durch Polygonzug genähert.Polygonseite (m) durch Gleichungdargestellt.( m ) ( m) ( m)( )v Q = Vo + k Δ Qaab( m) ( m)Es folgt ( 2)( + − ) + = ( 2) + +Δt Q Q Q V Δt Q Vo k ΔQzui zui+ 1 abi i abi+1ΔQ = Q − Qabi +1Beachte:Nach jedem Schritt muss die für die Berechnung maßgebende Polygonseite bestimmt werden.abiabInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 76Beispiel: Rückhaltebecken mit zusätzlichem ÜberlaufbeckenSchnitt 1-1Q initial = 35 m 3 /sQ OUT,initial = 35 m 3 /sΔt = 1 hourQin (m3/s)400350300250200150100500Inflow hydrograph1 2 3 4 5 6 7 8 9 10 11 12 13t (hours)Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 77Modell:dV = ( Q −Q ) ⋅dtzuAAdhdt( )= Q −Q ⋅ dtzuAUS( )ΔV = A ⋅ Δh = Q −Q −Q ⋅ Δt(1)A A A zu AUS ABΔV = A ⋅ Δh = Q ⋅ Δt(2)B A B ABn 1 nAus (1) ( ) ( )h + −h ⋅ A = Q −Q −Q ⋅ΔtAAA zu AUS AB+ 1 n+12 /A A A zu AUS ABh −h ⋅ A = Q −Q −Q ⋅Δtn noder ( ) ( )n[Euler Methode][Verbesserte Euler Methode]t + 1 Δn= ⋅ Qzu −QAUS − QAB+ hAn⇒ h( )n nAus (2) ( B B)AA+ 1 n ( oder n+1/ 2)h −h ⋅ A = Q ⋅ΔthnBBABt + 1 ΔA Q n= ⋅ h AB+BnBnAAusflussbedingung: Q = 12,5 h32 /AUSAbis 100 m 3 / s (Maximum)⇒ 100 = 12,5 32 /h A⇒ h A = 4,0 m⇒ Q AUS = 12,5 32 /h Ah A < 4,0 mQ AUS = 100 m 3 / sh A > 4,0 mQ AB = 36,5 h − h für h A > 4,5 mABQ AUS INITIAL = 35 m 3 / s ⇒ h A INIT =23 /⎛ 35 ⎞⎜ ⎟ = 199 , m⎝12,5⎠B anfangs trocken ⇒ h B = 0 + 4,5 (Vergleichsniveau) = 4,5 mInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 78Beispiel: Konstruktive und hydraulische Bemessung der Grundablässe und der Hochwasserentlastungfür RHB WartschenbrunnKonstruktionsprinzipDas Auslassbauwerk wird in Form einer Trogsperre mit Tauchwand ausgeführt. Der Grundablasswird in Form einer rechteckigen Öffnung am Fuß der Tauchwand ausgeführt. Zur späteren Nachregulierungdes Grundablasses ist wasserseitig ein verstellbarer Schieber (scharfkantig) vorgesehen.Die Hochwasserentlastung soll über die Krone der Tauchwand (abgerundet, mit lotrechter OW-Seiteund geneigter UW-Seite) ausgeführt werden. Als Tosbecken der Hochwasserentlastung dient gleichzeitigdas Tosbecken des Grundablasses, das innerhalb der Trogsperre angeordnet ist. Wasserseitigist eine Rechenkonstruktion vorgesehen.Auslassbauwerk Rückhaltebecken, Konstruktionsprinzip. Konstruktion und Design: HOLZINGERGerhard; WLSInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 79INSTATIONÄRE ROHRSTRÖMUNGDurch Regelvorgänge (z. B. mit Schieber, Ventile, Pumpen) kommt es zu einer zeitlichen Änderungdes Durchflusses ( ∂Q ∂t ≠ 0 bzw. ∂v ∂t≠ 0). Diese Änderung des Durchflusses bewirkt bei vollfließendenRohrleitungen den sogenannten Druckstoß. Nicht nur Regelvorgänge sondern Erscheinungenwie die Kavitation in Druckrohrleitungen, ungenügend entlüftete Saugrohrleitungen odereine unregelmäßige Förderung von Pumpen bewirken Druckstöße. Der Druckstoß resultiert aus derTrägheit der Flüssigkeitsmasse, die diese eine Änderung ihres Bewegungszustandes entgegensenkt.Generelle Annahmen für die Druckstoßberechnung:• Geschwindigkeit und Druck gleichmäßig über den Querschnitt• Geschwindigkeitshöhe

HYDRODYNAMIK SEITE 80mit A = τ( μ )μ RegelfunktionRt A Ro⇒ vt= voΔH1 +HoAARtRo= v τoΔH1+Hoτ ......... zeitabhängige RegelfunktionFür gleichförmiges Schließen gilt:undAARoRo− A− ARtRT=vttT⇒ τ = 1−AtT⎛ A⎜⎝= τ− AARoRt A RoRoRT⎞⎟ = 1−⎠⎡ t ⎛ vo− vT⎞⎤ΔH= vo⎢1−⎜ ⎟⎥ 1+⎣ T ⎝ v o ⎠⎦HotT⎛ v⎜⎝o− vvoT⎞⎟⎠Wird nun dv/dt gebildet und in die Ausgangsgleichung für ΔH eingesetzt, weiters diese Gleichungnach dΔH/dt aufgelöst so folgt für die maximale Druckhöhe (dΔH/dt=0).ΔH max K= ± K +KH 2 4o( )⎡Lvo− vwobei K = ⎢⎣⎢gH oTT2⎤⎥⎦⎥2+ Schließen− ÖffnenBei diesem Ansatz wird an sich kein Druckstoß berechnet (Elastizität vernachlässigt), daher nur fürkurze Leitungen und sehr langsamen Regelorganänderungen verwendbar.( − )ΔHmaxfür KK Lv o< 005 , : ≅ =vHogHoT1) Gleichförmiges Schließen2) Elastizität vernachlässigt3) langsame Regelorganverstellung[ m]LFaustformel T[ s ] > 300TFür ungleiche Rohrdurchmesser kann mit einer Ersatzlänge L´ gerechnet werden.L2A1LnA1L`=L1 + + ...... +AAN2Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 81THEORIE DER ELASTISCHEN WASSERSÄULEABLEITUNG DER GRUNDGLEICHUNGENVoraussetzungen bzw. vereinfachende Annahmen:- Vernachlässigung der Rohrreibung- Konstanter Leitungsdurchmesser, konstante Wanddicke und gleiches RohrleitungsmaterialA) Bewegungsgleichung (Impulssatz)m dvdt= ∑ Rpr 2 πdxdxT = 2τrπdxvr 2ρgr 2 πdx⎛ ∂p⎞π⎜p + ⎟dx⎝ ∂x⎠2ρrπdx dv 2r π ∂ p=dt ∂x dxdv ∂pρ =dt ∂x( )v = v x,tdv ∂v∂vdx= + v ist klein gegen ∂ vdt ∂t∂xdt∂t(v∼5 m/s bei Wasserkraftanlagen)Daher:∂v∂t∂p= 1ρ ∂xB) Kontinuitätsbedingunga) Geschwindigkeitsänderungd xEinströmendes Wasservolumen in der Zeit dt: V3v V2 r d x∂vv + dx ∂ xdr2π v dtAusströmendes Wasservolumen in der Zeit dt:Volumenbilanzr2 πΔV1⎛ ∂v⎞⎜ v+dx⎟dt⎝ ∂x ⎠2 v= r π ∂ dx dt∂xInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 82Dies ist nur möglich, wenn sich einerseits die Flüssigkeit (Wasser) unter einem erhöhtem Druckkomprimiert und wenn sich andererseits die Rohrleitung ausdehnt (quer zur Längsachse).b) Verkürzung der Wassersäule infolge DruckerhöhungΔdxHooke’sches Gesetz: σ = E wε ε =dxσ Spannung ε Dehnungσ = dp (zusätzliche Druckänderung)Δdxdp = E wdx∂p∂pdp = dt + dxp = p(x,t)∂t{ ∂xklein∂p∂t dt E Δdx= wΔdxdxp= 1 ∂ dx dtE ∂t wΔV = r π Δdx=222r π∂pdx dtE ∂t wc) Dehnung des Rohres infolge Druckerhöhungπ2∫0P = Z P = dp sin ϕ r dA = − dp cos ϕ r / = dp rσ s = dp rdp rσ =sπ20Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 83Hooke’sches Gesetz (Elastische Verformung):dp rs= σ = E ε = E∂p∂pdp = dt +∂t∂x∂prΔ r = dt∂tE sΔV= πr2rrdx {Δrrklein→0bzw.ErsΔrdp =2r2 2 ⎡ 22 2 ⎤[( r + Δr)− r ] dx = πdxr + 2rΔr+ Δ{ r − ⎥⎦3 rklein→0⎢⎣Δ V32πr=E sr3∂pdtdx∂tGesamte Kontinuitätsbedingungr2Δ V1 = ΔV2+ ΔV3v r p 2 rπ ∂ π∂ πdx dt = dx dt +∂x E ∂t Es2 3W∂v∂p⎛=⎜∂x∂t⎝1E+2rsW E r⎞⎟⎠r∂pdx dt∂t Wird⎛ 1 d ⎞ 1ρ ⎜ + =E2w sE⎟ eingesetzt so folgt:⎝ r ⎠ a∂v∂x1=ρa∂pt2 ∂a =⎛ 1ρ⎜⎝ Ew1d ⎞+ ⎟sE ⎠Geschwindigkeit der Druckfortpflanzung in Rohrleitungen( )a = a ρ, E , d, s,Ewrrρ, E w Stoffkennwerte der Flüssigkeitd, s, E r Kennwerte der RohrleitungInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 84Fortpflanzungsgeschwindigkeit der Schall- (Druck-)wellen in Wasser:d1Mit = 0in a =folgtsEr⎛ 1 d ⎞ρ⎜ +⎟⎝ Ew sE r ⎠Geschwindigkeit der Schallfortpflanzung in Wasser:Ewa =ρNach Messungen von Colladon und Sturm im Genfer See mit a = 1435 m/s ergibt sich der Elastizitätsmodulfür Wasser:E a ρ ≈ 210 . N mw =2 9 2Im Weiterem gilt für:Stahlrohre: EST = 2-2.2 1011N/m2Schleudergussrohre: EG = 1.0 1011N/m2Betonrohre: EB = 2-3 1010N/m2Für generelle Vorausberechnungen kann die Druckfortpflanzungsgeschwindigkeit in Stahlrohrleitungena ~ 1000 m/sangenommen werden.ALLGEMEINE LÖSUNGEN DES DIFFERENTIALGLEICHUNGSSYSTEMSDie Bewegungsgleichung∂v∂t∂p= 1ρ∂xund die Kontinuitätsgleichung∂v⎛ 1 2r⎞ ∂p1= ⎜ + ⎟ =∂x⎝ E dE ⎠ ∂tρawr2∂p∂tvermitteln als System linearer, partieller Differentialgleichungen den örtlichen und zeitlichen Zusammenhangzwischen der Strömungsgeschwindigkeit v (x,t) und dem Druck p (x,t).Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 85Die allgemeinen Lösungen dieses Systems lauten:p(x,t) = po + Φ(x-at) + Ψ(x+at)1 1vxt , = vo− Φ x− at + Ψ x+atρaρa( ) ( ) ( )Φ(x-at) und Ψ(x+at) sind zunächst beliebige Funktionen der Variablen x-at bzw. x+at, welche erstaus den jeweiligen Randbedingungen des speziell vorliegenden Problems bestimmt werden.Kontrolle für den Lösungssatz:pv( x, t) = P + Φ( x − at) + ψ( x + at)o1ρa1ρa( x, t) = v − Φ( x − at) + Ψ( x + at)o∂∂x∂∂tliefert:∂p∂x∂v∂t( x − at)( x − at)∂Φ ∂= 0 +∂142431 4243= 0 −1ρaΦ′( x − at)( x − at)( x − at) ∂Ψ( x + at)+∂x∂( x + at)∂Φ ∂∂142431 4243Φ′1( x + at)∂142431 4243 ∂xΨ′( x − at) 1 ∂Ψ( x + at)+∂tρa∂( x + at)−a1( x + at)∂142431 4243 ∂tΨ′+ a∂p∂v−11= Φ′ + Ψ′ und = Φ′ − a + ψ′∂x∂tρaρaund somit: ( ) aHieraus folgt∂v∂t1 ∂p=ρ ∂x(Bewegungsgleichung)Analog kann∂ v 1 ∂p=∂xρa2 ∂t(Kontinuitätsbedingung)hergeleitet werden.Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 86Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 87PHYSIKALISCHE BEDEUTUNG DER LÖSUNGSINTEGRALEFür einen Beobachter, der mit der Geschwindigkeit a in derpositiven x-Richtung (Speicherbecken x Regelorgan) fortschreitet, bleibt der aus Φ(x − at) resultierendeDruckanteil konstant.Beweis:p(x,t) = po + Φ(x - at) + (Ψ = 0, Annahme) = const.d. h.: dp = 0∂( x − at) ∂( x − at)dp = Φ′ dx + Φ′ dt = 0 (Totales Differential)14243∂x14243∂t1−adxdp = Φ′=dt( dx − adt) = 0 → aFür die Ortsveränderung mit dx/dt = a (Druckwellenfortpflanzungsgeschwindigkeit) ist p(x,t) konstant.Analog gilt:Für einen Beobachter, der mit der Geschwindigkeit a in dernegativen x - Richtung (Speicherbecken --> -x |x| (B) Regelorgan), fortschreitet, bleibt deraus Ψ(x+at) resultierende Druckanteil konstant.Beweis:p(x,t) = po + (Φ = 0, Annahme) + Ψ(x + at) = const.d. h.: dp = 0∂( x − at) ∂( x − at)dp = Ψ′ dx + Ψ′ dt = 014243∂x14243∂t1a(Totales Differential)dxdp = Ψ′=dt( dx − adt) = 0 → aFür die Ortsveränderung mit dx/dt = -a (Druckwellenfortpflanzungsgeschwindigkeit in negativer x-Richtung) ist p(x,t) konstant.Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 88Die Funktion Φ(x − at) stellt eine beliebig geformte Druckwelle dar, die sich ohne Deformation mitder Geschwindigkeit a in der positiven x-Richtung fortpflanzt. (A ← B)Die Funktion Ψ(x + at) stellt gleichfalls eine beliebig geformte Druckwelle dar, die sich ohne Deformationmit der Geschwindigkeit a in der negativen x - Richtung fortpflanzt. (A → B)Der gesamte Überdruck bzw. Unterdruck im Leitungsquerschnitt xi zum Zeitpunkt ti ergibt sich ausder Summe der Druckwellen zu:Δp(xi,ti) = p(xi,ti) - po = Φ(xi - ati) + Ψ(xi + ati)Die Geschwindigkeitsänderung Δv = v - v o im Querschnitt x i zum Zeitpunkt t i ist proportional derΨ x + at −Φx − at zufolge:Differenz der Druckwellen ( i i) ( i i)1Δvx ( , t) = vx ( , t) − v = Ψ( x + at ) −Φ( x −at)ρa[ ]i i i i o i i i iFür einen Beobachter im Querschnitt xi ändern sich naturgemäß andauernd sowohl Φ(x - at) undΨ(x + at).Zusammenhang zwischen der zeitlichen Veränderung des Durchflusses (Schließ- bzw. Öffnungsgesetz)und der Form der Druckwelle.Allmähliche Drosselung des Durchflusses:Man kann sich diesen stetigen Abschlussvorgang aus einer unzähligen Anzahl unstetiger, elementarerVorgänge vorstellen.Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 89Die Druckwelle pflanzt sich mit der Druckfortpflanzungsgeschwindigkeit vom Ort der Durchflussänderung(B) zum Speicherbecken (A) fort.2 LIst die Veränderungszeit τ < tr = ,dann tritt an Ort B die volleaDruckänderung Δ p = Φ (x - at) auf. (Direkter Stoß).DIREKTER STOß UND REFLEXIONSSTOßAm Einlaufquerschnitt A wird der statische Druck po aufrecht erhalten. Dies ist aber nur möglich,wenn eine Entlastungswelle (Gegenwelle) auftritt, welche die erste (direkte) Druckwelle egalisiert.Es muss also für x = L im Punkt A gelten:Δp = p − po = 0 = Φ(L−at) + Ψ(L+at) bzw.Φ( L − at) = −Ψ( L + at)Die Gegenwelle Ψ ist das im Abstand L gespiegelte Bild der direkten Welle Φ, jedoch mit negativenOrdinaten.Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 90Vollkommener Schieberabschluss: τ> t =r2LaIst τ > t =auf.r2Laso tritt im Punkt B infolge der Reflexion der abgeminderte Überdruck Δp < ΔpmaxVeränderungszeit (Direkter Stoß):Veränderungszeit (Stoß mit Reflexion):2Lτ ≤ tr = KΔp maxa2Lτ ≥ t = K Δp≤ Δar p maxGRAPHISCHE METHODE ZUR ERMITTLUNG DER DRUCKÄNDERUNG INFOLGE DESWASSERSTOßES NACH SCHNYDER, BERGERONGegeben:Strömungszustand im Leitungsquerschnitt m im Abstand l vonder Störungsstelle (z. B. Schieber) zum Zeitpunkt T.Strömungsgeschwindigkeit v 1,Tbzw. Durchfluss Q l,T = v l,T Ap ,TDruck p l,Tbzw. Druckhöhe y ,T = 11ρgGesucht:P x,t bzw. y x,t und v x,t bzw. Q x,t für Leitungsquerschnitt m 1im Abstand x von der Störungsstelle zum Zeitpunkt t.Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 91Zustandsgleichungen für m1:1y ,t − yo= Φ x − atρg[ ( ) + Ψ( x at)]x +− AQx ,t − Qo=+ρa[ Φ( x − at) − Ψ( x at)]Für einen Beobachter der von Punkt m mit Geschwindigkeit a in Richtung +x wandert, ist Φ konstant,denn für t > T istx = l + a(t-T) und somitΦ( x − at) = Φ[ 1 + a( t − T)− at] = Φ( 1 − aT) = Φim Querschnitt m.Hingegen ändert sich für diesen Beobachter Ψ!Diese veränderliche Funktion kann man eliminieren!Zustandsgleichung für m:y1,TQ1,THieraus folgt:ρgρg− yo− Qo1=ρg− A=ρg[ Φ( 1 − aT) + Ψ( 1 + aT)][ Φ( 1 − aT) − Ψ( 1 + aT)]ρa( y − y ) − ( Q − Q ) = 2Φ( − aT)1, T o 1,T o 1Aρa( y - y ) − ( Q − Q ) = 2Φ( x − at)x,toAx,to⎫ρg⎪⎬ +− ρa⎪A ⎪⎭⎫für m,⎪⎬ −bzw. für m 1⎪⎭Letztlich erhält man:ρa( y − y ) − ( Q − Q ) 0 bzw.ρg,t 1 ,T x,t 1Ax , T =Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 92( Q Q )a− Stoßgerade 1gAyx,ty1,T= x,t − 1,TStoßgerade 1 durch m y(l,T) und Q(l,T) mitagAtan α =Q = 1ZustandsebeneNeigungswinkel:⎛ ⎞α= arctan⎜a ⎟⎝ gA⎠Analog gilt für einen Beobachter der von Punkt m mit a gegen die Strömung nach m2, also in negativerx-Richtung wandert:( + at) = Ψ( + aT) = ΨΨ x 1 im Punkte m,da er zur Zeit t nach m1 gelangt ist:x =1 − ( t − T)aund somit( x + at) = Ψ( 1 − ( 1 − T)a + at) = Ψ( + aT)Ψ 1Daraus folgt analog dem Vorhergehenden:yx,t( Q − Q )− a− y 1 ,T = x,t 1,T Stoßgerade 2gAInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 93−aStoßgerade 2 durch m[y(l,T) undgAQ(l,T)] mit tan α=Q = 1Zur Behandlung konkreter Aufgaben müssen noch Randbedingungen gegeben sein, z. B. Schließeneines Schiebers, mit bekanntem Schließgesetz:ϕ( t )=ARt ()A RoQ(t) = A vo = A Ro v ARo = A R(t) v R(t)v R(t) = () 2gy tQ(t) = A R(t) 2gy()t == ϕ(t) A Ro 2gy()t == ϕ(t) AvvoARo( )( ) ( ) ( )Qt = Qϕ t,yt2gy()tParabelschar mit ϕ(t) als Parameter:Q = ϕAvv14243Roo() t ⋅ 2g y()tkQ =ϕ () t k y()tInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 94Graphische Methode Druckänderung nach Schnyder, BergeronInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 95Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 96Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 97Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 98Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 99ANHANGANHANG 1: LISTE DER VERWENDETEN FORMELZEICHENINTERNATIONALES Einheitssystem SI (Systeme Internationale) Größen, Namen undFormelzeichen nach ÖNORM A6401 und aus dem ON Handbuch 1, "Größen und Einheitenin Physik und Technik"Formelzeichen SI-Einheit GrößeZeichenNameA m 2 Quadratmeter FlächeninhaltV m 3 Kubikmeter Rauminhalt (Volumen)α,β,ν rad Radiant Ebener Winkelu,v,w,c ms -1 Geschwindigkeita ms -2 Beschleunigungω(Omega) rads -1 Winkelgeschwindigkeitρ(Rho) kgm -3 DichteFGN Newton (kgms -2 ) KraftGewichtp Pa Pascal (=Nm -2 = Drucktτ(Tau) =kgm 2 s -3 ) Schubsp.W J Joule (=Nm= ArbeitE =kgm 2 s -2 ) EnergieP W Watt (=J.s -1 =kgm 2 s -2 ) Leistungη(Eta)μ(My)(Pa.s)DynamischeViskositätν(Ny) (m 2 s -1 ) KinematischeViskositätl m Meter Längem kg Kilogramm Masset s Sekunde ZeitT K Kelvin Temperaturt °C GradCelsiusCelsius-TemperaturInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 100ANHANG 2: PHYSIKALISCHE EIGENSCHAFTEN UND KENNGRÖßEN VON REINEMWASSERTabelle 1: Physikalische Kenngrößen von reinem Wasser. Aus Studienunterlagen Hydraulics I, IHE(International Institute for Hydraulic and Enviromental Engineering)Dichte:ρ= dmdVZähigkeit (Viskosität):0,00178dynamische Viskosität (nach Poiseuille) ημ ( ) =1 + 0. 0337t+ 0,00022t2ηkinematische Viskosität ν =ρΔVVolumsänderung infolge Druckänderung: Δp=− EwV1 dV−Volumsänderung infolge Temperaturänderung: β = β = 18 10 5 je ° CV dtInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 101ANHANG 3: STOFFGESETZEStoffgesetze vermitteln den Zusammenhang zwischen den kinematischen Größen und den Spannungenbedingt durch die Deformation des Flüssigkeitselementes im Zuge des Strömungsablaufes(Volumenänderung und Gestaltänderung).Gasgesetz: p = ρ R T T.........absolute Temp.R ........Gaskonstanteρ.........Massendichtep.........DruckGrundgleichungen für die Strömung zäher, inkrompressibler Flüssigkeiten:Volumenänderung:Gestaltänderung:ρ = const τ = μdnInkompressibilitätNewton’scher Reibungsansatzdvtg β=dnn ... Normalabstand zu vGrundgleichungen:ρ dv = ρK− gradp + −μ∇2 v ( Navier/Stokes)dtund div v = O.∇22∂=∂x22∂+∂y22∂+∂z2∇........ Laplace’sche OperatorK........ MassenkraftGrundgleichungen für die Strömung idealer Flüssigkeiten :Volumenänderung:ρ = constInkompressibilitätGrundgleichungen:dvρ = ρK− gradp (Euler)dtund div v = O.Gestaltänderung:τ = OReibungsfreiheitInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 102ANHANG 4: PREISMANN SCHEMAKontinuität ∂Q∂h+ bs = 0∂x∂t2∂Q∂ ⎛ Q ⎞ ⎛ ∂h⎞ g Q QMoment + ⎜ β ⎟ + gA⎜+ Is⎟ +∂t∂xAx2⎝ ⎠ ⎝ ∂ ⎠ C ARDifferenzennäherung:⎛ n n ⎞ ⎛ n+1 n+1∂h⎞( ) ⎜h j + 1 − h j ⎟ ⎜h j + 1 − h j≅ 1 − Θ+ Θ⎟∂x⎜ Δx⎟ ⎜ Δx⎟⎝ ⎠ ⎝⎠⎛ n + 1 n ⎞ ⎛ n+1 n+1∂h⎞( ) ⎜h j − h j ⎟ ⎜h j + 1 − h j≅ 1− Ψ+ Ψ⎟∂t⎜ Δt⎟ ⎜ Δt⎟⎝ ⎠ ⎝⎠⎛ n n ⎞ ⎛ n+1 n+1∂Q⎞( ) ⎜Q j + 1 − Q j ⎟ ⎜Q j + 1 − Q j≅ 1 − Θ+ Θ⎟∂x⎜ Δx⎟ ⎜ Δx⎟⎝ ⎠ ⎝⎠⎛ n+1 n ⎞ ⎛ n+1 n+1∂Q⎞( ) ⎜Q j − Q j ⎟ ⎜Q j + 1 − Q j≅ 1 − Ψ+ Ψ⎟∂t⎜ Δt⎟ ⎜ Δt⎟⎝ ⎠ ⎝⎠⎛ 2 ⎡ n n+1 n n+1∂ ⎞ β⎥ ⎥ ⎤⎜Q Qβ ⎟ ≅ ⎢ j + 1 Q j + 1 Q j Q j−∂x⎝ A ⎠ Δxn+1 2 n+1 2⎢⎣Aj+1Aj ⎦⎡ n n+1 n n+1 ⎤g Q Q 1Q≅⎢ j Q j + 1 Q j + 1 Q j + 1g+2 ⎢ 2 n+1 2 2 n+1 2C AR 2⎢( C AR) ( ) ⎥ ⎥⎥ ⎣ j C AR j+1 ⎦Kontinunität:⎛ n n n 1 n 1n 1 nQ 1 Q j⎞ + +Q j 1 Q j⎛ +h j h j⎞ ⎛ h1 ⎜ + −⎟ + −bs1 ⎜−− Θ+ Θ+ − Ψ⎟ + b ⎜sΨ⎜ Δx⎟ Δx⎜ Δt⎟ ⎜⎝ ⎠⎝ ⎠ ⎝⎛ − Θ ⎞ n+ 1 bs( 1 − ψ)n+1 Θ n+1 ψ n+1⇒ ⎜ ⎟Qj + h j + Q j + 1 + bsh j + 1 =⎝123Δx⎠ 14243 Δt{ Δx123 Δt( ) j ( ) j + j + ⎟ = 0Bewegungsgleichung:⎛ Q1−Ψ⎜⎝( )A1jB1j( −1)( 1− Θ ) b( )sQ − Q + [( 1− Ψ)h + h ]=jj +1444444444Δx24Δ44444444t3n + 1jC1jD1n n n+ 1 jj Ψjn1n+11E1 jj n j 1j n j 1j n j+ 1 1 j n j+ 1 1j− hA1 Q + B1 h + C1 Q + D1 h = E1(3)− QΔtnj⎞ ⎛⎟ +QΨ⎜⎠ ⎝⎞n+ + 1 1 nn nj − Qj+ 1 1 1 1⎟ + β ⎜Qj+ Qj+ +Δt12⎠ Δx⎜n +A j+1⎛⎝−QAΔtn n+1j Qjn+12jn n 1Q Qn n 1⎡ n n n 1 n 1h h h hQ QgA n 1 2⎛( 1 ) j 1 j⎞ ⎛ + +⎡ ++j 1 j⎞ ⎤⎤+j jj 1 j 1⎢ ⎜ + −⎟ ⎜ + −Isg 2⎢+ ++⎟⎥= 0j 1 2− ⎥+− Θ+ Θ++xx⎢ 2 n 1 2 2 n 1 2⎢ ⎜ Δ ⎟ ⎜ Δ ⎟++⎥⎥⎣ ⎝ ⎠ ⎝ ⎠ ⎦ ⎢( C AR) j( C AR)⎣j+1 ⎥⎦⎞⎟⎟⎠n1⎞⎟⎠Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 103⎡⎢( − ψ)⇒⎢⎢⎣14 2n + 1 2( − Θ) gAj+1 2 1n n1 β Q jg Q j n + 1n +− +⎥Qn 1 2n 1 2( )j +h jΔtΔx+A2 + ⎥xj 2 C ARΔj ⎥1 44243 44444444 44444444⎦ 3 B2 jA2 j⎡n n n 1 2Q j 1g Q ⎤+j + 1Θ⋅gA⎢ ψ β +n+1 j+1 2+ +⎥Qn 1 2n 1 2( )j + 1 + h⎢ΔtΔx+A2 + ⎥xj+1 2 C ARΔ⎢⎣j+1 ⎥142431444444424444444⎦ 3 D2 j( − ψ)c2 jn+1 2⎤( 1 − Θ)( h − h )n+1j + 11 n Ψ gAn j+1 2n n n 1 2= Q j + Q j + 1 −j + 1 j + gA+j+1 2Is14Δ44444444444t ΔtΔx24444444444443E2 jj n j + 1j n j + 1j n j+ + 1 1 j n j+ +1 1jA2 Q + B2 h + C2 Q + D2 h = E2(4)Recurrence RelationA1 Q + 1+ B1 h + 1+ C1 Q + + + D1 h + + = E1(3)j n j j n j j n j 1 1 j n j 1 1jj n j + 1j n j + 1j n j+ + 1 1 j n j+ +1 1jn + 1j j n + 1= j + jA2 Q + B2 h + C2 Q + D2 h = E2(4)Substitution Q F h Geingesetzt in (3)( )( )j j n j + 1j j n j + 1j n j+ + 1 1 j n j+ +1 1jA1 Fh + G + B1 h + C1 Q + D1 h = E1j j j n +j 1j n j+ + 1 1 j n j+ +1 1j j jA1 F + B1 h + C1 Q + D1 h = E1 −A1Gn+1− C1j +−+−⇒ h =++++A F + B Q n 1D1jA F + B h n 1E1j A1jGjjj 1j 11j j1j1j j1jA1jFj + B1j1424314243142443H j I j J jn + 1jh = H Q + I h + J(5) und (6) in (4) eingesetztA2 Fh + 1 1+ G + B2 h + + C2 Q + + + D2 h + + = E2j n j+ + 1 1 j n j+ +1 1 j(6)j ( j n j j)j n j j n j 1 1 j n j 1 1j( A jFj B j)h n +2 + 2 j 1 + C2j Q n j+ + 1 1 + D2j h n j+ +1 1 = E2j−A2jGj( )( n + 1 n + 11) n + 1 n +A2j F j + B2 j H j Q j + 1 + I j h j + 1 + J j + C2 j Q j + 1 + D2 j h j + 1 = E2 j − A2 jGj( )+ + 1 1 + +1 1[ 2 + 2 + 2 ] + [( 2 + 2 ) + 2 ] = 2 − 2 − ( 2 + 2 )A F B H C Q A F B I D h E A G A F B Jj j j j j j n j j j j j j n j j j j j j j( )( )− − ( + )( A2jFj + B2j)Hj + C2j⇒ = − A2 F + B I − D++Q+++A F + B H + C h E A G A F B Jn j j2j j2j2n j2j j2j j211j jj 1j 12j j j j2j144442444431444442444443Fj+1nj+ + 1 1 j+ 1 n j+ + 1 1 j+1Q = F h + GG j+1+(5)Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 104LITERATURVERZEICHNISABBOTT M. B.: Computational Hydraulics; Elements of the Theory of Free Surface Flows, Pitman,1979.ABBOTT M. B., BASCO D. R.: Computational Fluid Dynamics an Introduction for Engineers,Longman Scientific & Technical, 1990.ABBOTT M. B.: HydroinformaticBLOß S.: Numerische Modelle von Flüssen, Seen und Küstengewässern, DVWK Schriften 127,1999.BOGARDI J.: Sediment transport in alluvial streams, Akademia Kiado, Budapest, 1974.BOLLRICH G.: Technische Hydromechanik/2, VEB Berlin,BOUVARD M.: Mobile barrages and intakes on sediment transporting rivers, Balkema, 1992.BREINER H.: Neue Gesichtspunkte zur Wasserspiegellageberechnung bei stationären Abflüssen inoffenen Gerinnen, Teil 1, Theoretische Grundlagen, Österr. Wasserwirtschaft, Heft5/6, 1989.BREINER H., HENGL M.: Neue Gesichtspunkte zur Wasserspiegellageberechnung bei stationärenAbflüssen in offenen Gerinnen, Teil 2, Praktische Berechnungsmethodik, Österr. Wasserwirtschaft,Heft 5/6, 1989.BREINER H.: Die durch Integration diskretisierten Grundgleichungen zur eindimensionalen Beschreibungvon Abflussvorgängen in offenen Gerinnen. Österr. Wasserwirtschaft, Heft 1/2, 1990.CHADWICK A., MORFETT J.: Hydraulics in Civil Engineering, Collins, 1986.CUNGE H. Jr., VERWEY: Practical Aspects of Computational River Hydraulics, Pitman, 1980.DE VRIES M.: Scale Models in Hydraulic Engineering, IHE-Delft, 1986FOX J.A.: Hydraulic Analysis of Unsteady Flow in Pipe Networks, Macmillan, 1979.FRANKE P. G.: Instationäre Strömung in Druckleitungen, Abriss der Hydraulik 9, BauverlagWiesbaden und Berlin, 1974.HANXIANG X.: Characteristics of the overbank flow and related hydraulic computations, Proceedings21. Intern. Congress of IAHR, 1985KASTANEK F., LOISKANDL W.: Kriterien der Spiegellagenberechnung regelmäßiger und unregelmäßigerProfile. Wiener Mitteilungen Band 79, Wien 1989.KOBUS H.: Hydraulic Modelling, Parey, 1980.KRISHNAPPAN B. G.: Unsteady, Nonuniform, Mobile Boundary Flow Model-Mobed, EnviromentCanada, 1981.KRISHNAPPAN B. G.: Mobed Users Manual Update I, Enviroment Canada, 1983.KRISHNAPPAN B. G.: Mobed Users Manual Update II, Enviroment Canada, 1986.KUNDU P.: Fluid Mechanics, Academic Press, 1990.Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 105LASARTE A. E.: A mathematical model for flood plain sedimentation Master Thesis,IHE-Delft,1989LOISKANDL W.: A mathematical model for the flow in a single channel with a rectangular crosssection,IHE-Delft,1989.LOISKANDL W.: Diffusive -Sediment transport towards flood plain, IHE-Delft, 1989.LOISKANDL W.: Exercise Free Surface Flow Modelling, IHE-Delft, 1989.LORENZ P.: Adaption und Kritische Beurteilung eines numerischen Modells zur Simulation desAbflussgeschehens in offenen Gerinnen, Diplomarbeit BOKU, 1990.MALCHEREK A.: Physik und Numerik der Oberflächengewässer, Bundesanstalt für Wasserbau,Hamburg, BRD, 2001.NAUDASCHER E.: Hydraulik der Gerinne und Gerinnebauwerke, Springer Verlag, 1987.NAUDASCHER E.: Hydrodynamic Forces, Balkema, 1991.NEZU I. and NAKAGAWA H.: Turbulence In Open-Channel Flows, IAHR-Monograph, Balkema,Rotterdam, Brookfield, 1993.NINOMIYA; O.: Flow Analysis Using a PC, CRC Press, 1991.NOVAK P.: Water hammer and surge tanks, IHE-Delft, 1983.NOVAK P.: Hydraulic Structures, Hyman, 1990.PEYRET R., TAYLOR T. D.: Computational Methods for Fluid Flow, Springer-Verlag, 1985.PREISSLER, BOLLRICH: Technische Hydromechanik Band 1, VEB Berlin,RODI W.: Turbulence Models and Their Application in Hydraulics, A State –of-the-art review,IAHR-Monograph 3 rd ed., Balkema, Rotterdam, Brookfield, 1993.ROUVE G.: Hydraulische Probleme beim naturnahen Gewässerausbau, Forschungsbericht, VCH,1987SCHRÖDER R. C. M.: Hydraulische Methoden zur Erfassung von Rauheiten, DVWK Schriften 92,Parey, 1990.STEPHENSON: Pipeflow Analysis, Developments in Water Science Elsevier, 1984.TOWNSON J. M.: Free-Surface Hydraulics, Hyman, 1991.VENNARD, Street D.: Elementary Fluid Mechanics, 5. Auflage, WileyVERWEY A.: International Courses in Hydraulic and Sanatary Engineering, Delft, 1973.VISCHER, HAGER: Hochwasserrückhaltebecken, VDF, 1992.VOGEL H., 1995: Gerthens Physik, Springer-Verlag.VREUGDENHIL C. B.: Computational Hydraulics, Springer-Verlag, 1989.ZIELKE W.: Elektronische Berechnung von Rohr- und Gerinneströmungen, Erich Schmidt Verlag,1974.Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 106ZUSAMMENFASSUNG DER WICHTIGSTEN FORMELN DER VORLESUNG815.100 HYDRAULIK UND HYDROMECHANIKPhysikalische Kenngrößendm3Dichte : ρ = ρ⎡kg / m ⎤dV ⎣ ⎦ xrelative Dichte : ρrel=ρ2H OViskositätF dvτ = = η⋅A dnWOberflächenspannung:⎡Nm⎤σ =2A ⎢⎣ m ⎥⎦σ = σ − σ = σ ⋅ cosαw GW FW GF wKapillardruckKapillare Steighöhe:hKP K1 1= σ ⎛⎜ +⎞⎟⎝r1 r2⎠4σ= ⋅ cosαρgdkGrundgesetze der HydrostatikFlüssigkeitsdruckDruckkraft dFp = = DFläche dA( x y z )p x = p y = p z = pdp = ρ a dx + a dy + a dz 1. Grundgleichung von Eulerp ü = p - p 0 = ρg hph = ü .... Druckhöhe .ρgEbene Flächen:F = ρ g h s A2x=s+sI I A z ( Satz von Steiner)zD2Is+ AzsIs=zAs= + zz AssInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 107gekrümmte FlächenA X A x∫ ∫F = dF = ρ g z d A = ρ g z Ax x x sx xAyAy∫ ∫F = dFy = ρ g z d A = ρ g z Ay y sy yKraft und einem MomentzAzAz∫ ∫F = dFz = ρ g z d A = ρ gVz2 2 2x y zF = F + F + F cos α = F x / F; cos β = F y / F; cos γ = F z / FTranslationΔ h aa L= ⇒ Δ h= ⋅ RotationL 2gg 22ω = g⋅ΔhroGrundgesetze der Fluidmechanik∂ s∂ Qv = limv = lim∂ t →0∂ t∂ A→0∂ Adv ∂ v= + ( v∇) vdt ∂ tEuler’schen Bewegungsgleichung:∂ v1+ ( v∇ ) v = K − grad p∂tρ,Navier-Stokes;dv 1= K − grad p+Z;dt ρKontinuitätsbedingung∂ρEin − Aus = ΔS Δ S = dx dy dz∂ t∂ρ ∂( ρv) ∂ ( ρvy)x∂( ρvz)∂ρ+ + + = + div ( ρv) = 0∂t ∂ x ∂ y ∂z ∂t∂v∂ vx y ∂vz+ + = div v = 0∂ x ∂ y ∂zQ = v A = v A = constm1 1 m22.Die Bernoulli'sche Druck- und Energiegleichungdv ∂ z 1 ∂ pEulergleichung in Richtung der Stromlinie: =−g−dt ∂ s ρ ∂s2v p+ + gz = E = const ( Energie)2 ρInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 1082v ρ+ p + ρ g z = p = const ( Druck)22v pz H const2 g+ ρg+ = =2v2gpρgzH: Geschwindigkeitshöhe: Druckhöhe: geodätische Höhe, auf einen Vergleichshorizont bezogen (Ortshöhe).: mechanische Energiehöhe.Beschleunigung quer zur Stromlinienrichtung (Zentripetalbeschleunigung)Newton'sche Kraftgesetz in Richtung der Hauptnormalen:2 2v F v dvFn= m ⇒ = =r m r dt⎛ dv ∂z 1 ∂ρ v⎜⎞ ⎟ =− g − =⎝ dt ⎠ ∂n ρ ∂n rn2vr21 ∂ p=− (bei Vernachlässigung der Schwere)ρ ∂ nAusfluss aus kleinen Öffnungenv = 2 g h TORICELLI´sches THEOREM → Q = v A=μ A 2 g h.Schwingung einer Wassersäule21 d ξ= ξ + h2 2rh r ......Reibungsverlusteα dtbei turbulenter Wasserbewegung2 22ξ ⎛dξ⎞− m + μξ = 02ddt⎜⎝dt⎟⎠h r ∼ v⎛ ξ ⎞≈ ⎜ ⎟⎝ dt ⎠2 dgedämpfte Schwingung2Erweiterte Bernoulli-Gleichung2v p+ + z+ hv= konstant = H0.2gρgdHNeigung der Energielinie I =− =dxdh vdxInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 109hvi2v= ξ . und Formel nach DARCY-WEISBACH2gI2v= λ2 gdhv2v l= λ (quadratisches Widerstandsgesetz für turbulente Rohrströmung),2g dLaminar: λ= 64Reλ : Widerstandszahl.kλ(Re, )dImpulssatzdIdt= ∑ R (Impuls pro Zeiteinheit)Impulssatz in der HydromechanikdIdtdIdt− = ρ ( 2 − 1)=∑1 2Qv v R∑ R = FG+ FR+ p1 A1 + p2 A2 F G GewichtskraftF Rp i A iReaktionskraft der WandWasserdruckkraft in den QuerschnittenDer Stützkraftsatz besagt, dass die vom Mantel einer Stromröhre auf die strömende Flüssigkeitausgeübten Kräfte zuzüglich der Massenkraft (Gewicht) mit den in den geschnittenen Fließflächenangreifenden Stützkräften im Gleichgewicht stehen. Die Summe aus Druckkraft und Impulsstromergibt die Stützkraft F = SpA +ρ Qv .DE SAINT VENANT KoeffizientenImpulssatz: β=2∫ v dAA2v A. Werte für den Korrekturbeiwert β:für laminare Bewegung in Röhren: β = 1,33für turbulente Bewegung in Röhren: 1,01 < β < 1,07Freispiegelgerinne: 1,01

HYDRODYNAMIK SEITE 110Laminare und turbulente StrömungvmdREYNOLD´schen Zahl Re =νvt () = ν + v′() t ,v m ......... mittlere Geschwindigkeitd........... Rohrdurchmesserν........... kinematische ZähigkeitTurbinenleistung P = ρ g Q H n η.Förderhöhe HGerinnehydraulikvQρHman= H + + h Pa= .2gηman geo vges(als Strickler-Formel bezeichnet): v = k R I2STQ = A k ST R 2/3 I 1/213 / −126kST[ m s ] =16 /k2/ 3 1/2Schubspannungsgeschwindigkeit v * = g h I⎛ ksCuniverselle Fließformel v =− C1log ⎜ +⎝ C R 4 Rev ⎛ ksC ⎞3oder =− 2,828 C*1log ⎜ + ⎟v ⎝C2R 4Re λ ⎠23⎞⎟λ ⎠2 2 gRIQ = K I Ebhg = 2Strömender und schließender Abflussvorgang2Q B=gA31Froudezahl FRv= = StrömungsgeschwindigkeitghFortpflanzungsgeschwindigkeiteiner kleinen StörungInstitut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7

HYDRODYNAMIK SEITE 111Rechteckquerschnittqhgr = 3 α max2vggrg= h gr αFließzustand Wassertiefe Fließgeschwindigkeit Gefälle FroudezahlStrömen h > h gr v < v gr I E < I EKRIT F r < 1Kritischer Fließzustand h = h gr v = v gr I E = I EKRIT F r = 1Schießen h < h gr v > v gr I E > I EKRIT F r > 1Grundgleichung Wasserspiegellage:dhdxIo− I=2Q B1−3gFIo− I=21−FrPRAKTISCHE V0RGANGSWEISE zur WASSERSPIEGELANALYSE1. Unterteilung des Längenprofils des Gerinnes in Abschnitte gleicher Charakteristika.2. Berechnung der kritischen Tiefe h gr (für Fr 2 =l) und der Normalabflusstiefen h N (mit bekanntenAbflussformeln).3. Zuordnung der entsprechenden Profiltypen (M,S,C,H,A) zu den einzelnen Bereichen;gegebenenfalls sind noch die konjugierten Tiefen zu ermitteln. Überprüfung des Längenprofilsauf Kontrollquerschnitte.3. Voranalyse der möglichen Ausbildungsformen der Wasseroberflächen mit Hilfe der12 theoretisch möglichen Kurvenausbildungsformen. Danach Berechnung der Wasserspiegellagenin den einzelnen Bereichen.In Bereichen strömenden Abflusses hat die Berechnung gegen die Fließrichtung zu erfolgen, inBereichen schießenden Abflusses in Fließrichtung.Formel von Poleni2Q = μ B 2g h33/2Institut für Hydraulik und landeskulturelle Wasserwirtschaft Version 1.7