HYDRODYNAMIK SEITE 36MATHEMATISCHE GRUNDLAGENDen Ausgangspunkt der theoretischen Überlegungen für die dreidimensionale Strömungsbeschreibungbilden schon seit jeher die Kontinuitätsgleichung <strong>und</strong> die Impulsgleichungen nach Navier-Stokes. Letztere sind sowohl für die laminare als auch turbulente Bewegung, die den Regelfall imAbflussgeschehen in künstlichen <strong>und</strong> natürlichen Gerinnen darstellt, geeignet. Eine Gr<strong>und</strong>annahmein 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ähigkeit............. x,y,z Koordinatenrichtungent ........... Zeitv x, , v y , v z Geschwindigkeitskomponente in KoordinatenrichtungX,Y,Z Richtungsvektoren der Kraft Kp Druckν dynamische Zähigkeitµ kinematische Zähigkeitρ DichteDie Kraft K(X,Y,Z) beinhaltet die konservativen Anteile wie Schwerkraft, Corioliskraft <strong>und</strong> 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 <strong>und</strong> 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 <strong>und</strong> 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 soll diese,ausgehend von den an einem Volumselement angreifenden Vektoren U,V,W, näher erläutertwerden. Der innere Spannungszustand infolge Zähigkeit ist, unter Zugr<strong>und</strong>elegung des Newton’schenZähigkeitsgesetztes, durch einen Spannungstensor, gebildet aus den Normalspannungen<strong>und</strong> den Tangentialspannungen, bestimmt.<strong>Institut</strong> für <strong>Hydraulik</strong> <strong>und</strong> <strong>landeskulturelle</strong> <strong>Wasserwirtschaft</strong> Version 1.0
HYDRODYNAMIK SEITE 37∂Ux − Richtung : dxdydz∂xFür die anderen Koordinatenrichtungen giltanalog:∂Vy − Richtung : dxdydz∂y∂Wz − Richtung : dxdydz∂zσ iτ i,jNormalspannungKraftwirkung der viskosen Spannung in x-RichtungTangentialspannung<strong>und</strong> 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- <strong>und</strong> z-Richtung ergeben sich analog.Für die turbulente Strömung müssen zu den Normalspannungen <strong>und</strong> Tangentialspannungen nochdie sogenannten scheinbaren turbulente Spannungen hinzugefügt werden.______ _________ _________⎛2⎞⎜ v' v'xv'y v'xv' ⎟xz⎜⎟_________ ______ _________Reynoldsspannungstensor ⎜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'zyz<strong>Institut</strong> für <strong>Hydraulik</strong> <strong>und</strong> <strong>landeskulturelle</strong> <strong>Wasserwirtschaft</strong> Version 1.0