www.wikidata.de-de.nina.az
Die Stromung eines Fluids Flussigkeit oder Gas ist eine Potentialstromung wenn das Vektorfeld der Geschwindigkeiten mathematisch so geartet ist dass es ein Potential besitzt Das Potential kann man sich anschaulich als die Hohe in einer Reliefkarte vorstellen wo dann die Richtung der grossten Steigung in einem Punkt der dortigen Geschwindigkeit entspricht Ein solches Potential ist in einem homogenen Fluid vorhanden wenn die Stromung rotationsfrei wirbel bzw vortizitatsfrei ist und keine Zahigkeitskrafte Reibungskrafte auftreten oder diese vernachlassigbar klein sind Jede aus der Ruhe heraus beginnende Stromung eines homogenen viskositatsfreien Fluids besitzt ein solches Potential 1 Eine Potentialstromung ist der rotationsfreie Spezialfall der Stromung eines homogenen viskositatsfreien Fluids das durch die Euler schen Gleichungen beschrieben wird diese gelten auch fur Stromungen mit Rotation Wirbelstromung Wenn jedoch bei Scherbewegungen die Zahigkeit berucksichtigt werden muss wie z B in Grenzschichten oder im Zentrum eines Wirbels so ist mit den Navier Stokes Gleichungen zu rechnen Potentialstromungen konnen als sehr gute Naherung von laminaren Stromungen bei niedrigen Reynolds Zahlen verwendet werden wenn die fluiddynamische Grenzschicht an den Randern der Stromung keine wesentliche Rolle spielt In der stationaren Potentialstromung inkompressibler Fluide gilt die bernoullische Druckgleichung global die technische Rohrstromungen gut beschreibt Wegen ihrer einfachen Berechenbarkeit werden Potentialstromungen auch als Anfangsnaherung bei der iterativen Berechnung der Navier Stokes Gleichungen in der numerischen Stromungsmechanik verwendet Inhaltsverzeichnis 1 Definition 2 Anwendungsbereiche und Limitierungen 3 Bestimmungsgleichungen fur die Stromung 3 1 Bilanzgleichungen 3 2 Barotropes Fluid in konservativem Schwerefeld 3 3 Stationare Stromung 3 4 Schallwellen 3 5 Randbedingungen 4 Zirkulation 5 Potentialstromung eines inkompressiblen Fluids 5 1 Bestimmungsgleichungen 5 2 Analogon der Warmeleitung 5 3 Ebene Potentialstromung 5 3 1 Komplexes Geschwindigkeitspotential 5 3 2 Naherungsweise rechnerische Losung 5 3 3 Zeichnerische Ermittlung von Potentialstromungen 5 3 4 Konstruktion von Geschwindigkeitspotentialen 5 3 5 Kraftwirkungen auf umstromte Korper 5 3 6 Potentialwirbel 6 Beispiel 7 Weblinks 8 Einzelnachweise 9 LiteraturDefinition BearbeitenEine Potentialstromung ist eine Stromung die ein Geschwindigkeitsfeld v x t displaystyle vec v vec x t nbsp besitzt das vom Ort x displaystyle vec x nbsp und von der Zeit t abhangt und sich aus dem Gradient grad eines Geschwindigkeitspotentials ϕ x t displaystyle phi vec x t nbsp berechnet v grad ϕ v i ϕ x i i 1 2 3 displaystyle vec v operatorname grad phi quad Leftrightarrow quad v i frac partial phi partial x i i 1 2 3 nbsp Deshalb werden solche Stromungen Potentialstromungen genannt Die linke Vektorgleichung ist die koordinatenfreie Version wahrend die rechten Gleichungen in kartesischen Koordinaten x 1 2 3 displaystyle x 1 2 3 nbsp gelten Jedes Gradientenfeld ist rotationsfrei weshalb Potentialstromungen immer rotationsfrei sind Umgekehrt existiert nach dem Poincare Lemma immer dann ein Geschwindigkeitspotential wenn die Stromung rotationsfrei ist Anwendungsbereiche und Limitierungen BearbeitenPotentialstromungen enthalten nicht alle Charakteristika von realen Stromungen Die inkompressible Potentialstromung macht eine Reihe von falschen Voraussagen wie beispielsweise das d Alembert sche Paradoxon demgemass von der Stromung auf einen Korper keine Kraft in Richtung der Stromung ausgeubt wird Alle Phanomene die eine hydrodynamische Grenzschicht oder turbulente Stromung nebst Dissipation von Energie beinhalten wie z B Stromungsabrisse konnen mit Potentialstromungen nicht abgebildet werden Nichtsdestotrotz ist das Verstandnis von Potentialstromungen in vielen Bereichen der Stromungsmechanik hilfreich So konnen fur nicht zu komplizierte Geometrien analytische Losungen ihrer Umstromung berechnet werden und auch die von der Stromung geleistete Auftriebskraft angegeben werden So wird ein tieferes Verstandnis der Stromung erreicht Potentialstromungen finden viele Anwendungen in der Auslegung von Flugzeugen Wie eingangs erwahnt konnen Potentialstromungen als sehr gute Naherung von laminaren Stromungen bei niedrigen Reynolds Zahlen verwendet werden wenn die hydrodynamische Grenzschicht an den Randern der Stromung keine wesentliche Rolle spielt Eine Technik in der numerischen Stromungsmechanik koppelt eine viskose Grenzschichtstromung an eine Potentialstromung ausserhalb der Grenzschicht In einer Potentialstromung kann jede Stromlinie durch eine Wand ersetzt werden ohne die Stromung dadurch zu verandern eine Technik die im Flugzeugdesign angewendet wird Bestimmungsgleichungen fur die Stromung BearbeitenNicht jedes Geschwindigkeitspotential reprasentiert eine physikalisch plausible Stromung Damit das Geschwindigkeitspotential unter den getroffenen Annahmen im Einklang mit den physikalischen Gesetzen ist muss es der Bilanzgleichung fur den Impuls in Form der Euler schen Gleichungen gehorchen und die Massenbilanz erfullen Bilanzgleichungen Bearbeiten Die Impulsbilanz in Form der Euler schen Gleichungen lautet v t grad v v 1 r grad p k displaystyle frac partial vec v partial t operatorname grad vec v cdot vec v frac 1 rho operatorname grad p vec k nbsp und die Massenbilanz r r div v 0 displaystyle dot rho rho operatorname div vec v 0 nbsp Darin ist r die Dichte und der aufgesetzte Punkt die substantielle Zeitableitung div die Divergenz eines Vektorfeldes p der Druck und k displaystyle vec k nbsp ein Beschleunigungsfeld z B Schwerebeschleunigung Dieses System aus vier Gleichungen mit funf Unbekannten drei Geschwindigkeiten Druck und Dichte wird durch eine Zustandsgleichung geschlossen die die Dichte als Funktion des Drucks darstellt Einsetzen des Potentials in die Euler schen Gleichungen liefert mit der Grassmann Entwicklung 2 v t grad v v 1 r grad p v t 1 2 grad v v v rot v 1 r grad p grad ϕ t 1 2 grad ϕ 2 1 r grad p k displaystyle begin aligned frac partial vec v partial t operatorname grad vec v cdot vec v frac 1 rho operatorname grad p amp frac partial vec v partial t frac 1 2 operatorname grad vec v cdot vec v vec v times operatorname rot vec v frac 1 rho operatorname grad p amp operatorname grad left frac partial phi partial t frac 1 2 operatorname grad phi 2 right frac 1 rho operatorname grad p amp vec k end aligned nbsp weil die Rotation eines Gradientenfeldes verschwindet Die Massenbilanz formt sich mit dem Geschwindigkeitspotential und dem Laplace Operator D d i v g r a d displaystyle Delta operatorname div circ grad nbsp um zu r r d i v g r a d ϕ r r D ϕ 0 displaystyle dot rho rho operatorname div grad phi dot rho rho Delta phi 0 nbsp Nun handelt es sich um ein System aus vier Gleichungen mit drei Unbekannten Potential Druck und Dichte das nur losbar ist wenn die Differenz f k 1 r grad p displaystyle vec f vec k tfrac 1 rho operatorname grad p nbsp die Integrabilitatsbedingung rot f 0 displaystyle operatorname rot vec f vec 0 nbsp erfullt Dies ist in barotropen Fluiden in einem konservativen Schwerefeld der Fall Barotropes Fluid in konservativem Schwerefeld Bearbeiten In einem barotropen Fluid ist die Dichte nur eine Funktion des Drucks Dann kann die Integration der Euler Gleichungen vorab erfolgen was die Berechnung wesentlich vereinfacht Barotrope Stromungen liegen dann vor wenn das Temperaturfeld homogen ist oder wenn im technisch besonders wichtigen Fall die Stromung homentrop ist 3 In einem barotropen Fluid gibt es eine Druckfunktion P mit der Eigenschaft P d p r p d P d p r grad P 1 r grad p displaystyle P int frac mathrm d p rho p quad Leftrightarrow quad mathrm d P frac mathrm d p rho quad Leftrightarrow quad operatorname grad P frac 1 rho operatorname grad p nbsp In einem inkompressiblen Fluid ist die Dichte konstant und P p r Einsetzen der Druckfunktion in die Euler Gleichung erlaubt es den Gradientenoperator heraus zu ziehen grad ϕ t 1 2 grad ϕ 2 1 r grad p grad ϕ t 1 2 grad ϕ 2 P k displaystyle operatorname grad left frac partial phi partial t frac 1 2 operatorname grad phi 2 right frac 1 rho operatorname grad p operatorname grad left frac partial phi partial t frac 1 2 operatorname grad phi 2 P right vec k nbsp Damit diese Gleichung erfullt werden kann muss das Beschleunigungsfeld k displaystyle vec k nbsp der Integrabilitatsbedingung rot k 0 displaystyle operatorname rot vec k vec 0 nbsp entsprechen Dann gibt es nach dem Poincare Lemma ein Potential V mit der Eigenschaft k grad V displaystyle vec k operatorname grad V nbsp und die obige Gleichung fuhrt auf ein Gradientenfeld grad ϕ t 1 2 grad ϕ 2 P V 0 displaystyle operatorname grad left frac partial phi partial t frac 1 2 operatorname grad phi 2 P V right vec 0 nbsp das uberall verschwindet weshalb der Term in den grossen Klammern ortsunabhangig ist ϕ t 1 2 grad ϕ 2 P V C displaystyle frac partial phi partial t frac 1 2 operatorname grad phi 2 P V C nbsp Die Integrationskonstante C konnte noch von der Zeit abhangen aber diese Zeitabhangigkeit kann dem Potential ϕ displaystyle phi nbsp zugeschlagen werden ohne dass sich dessen physikalische Bedeutung grad ϕ v displaystyle operatorname grad phi vec v nbsp andern wurde 4 Nun handelt es sich um ein geschlossenes System aus drei Gleichungen inklusive r r p und der Massenbilanz mit drei Unbekannten Potential Dichte und Druck Stationare Stromung Bearbeiten In einer stationaren Potentialstromung ohne Beschleunigungsfeld gilt 5 1 M x 2 2 ϕ x 2 1 M y 2 2 ϕ y 2 1 M z 2 2 ϕ z 2 2 M x M y 2 ϕ x y 2 M y M z 2 ϕ y z 2 M z M x 2 ϕ z x 0 displaystyle begin aligned amp left 1 M x 2 right frac partial 2 phi partial x 2 left 1 M y 2 right frac partial 2 phi partial y 2 left 1 M z 2 right frac partial 2 phi partial z 2 amp 2M x M y frac partial 2 phi partial x partial y 2M y M z frac partial 2 phi partial y partial z 2M z M x frac partial 2 phi partial z partial x 0 end aligned nbsp mit den Machzahlen M x 1 c ϕ x M y 1 c ϕ y und M z 1 c ϕ z displaystyle M x frac 1 c frac partial phi partial x qquad M y frac 1 c frac partial phi partial y qquad text und qquad M z frac 1 c frac partial phi partial z nbsp und der lokalen Schallgeschwindigkeit c Die Stromungsgeschwindigkeit v grad ϕ displaystyle vec v operatorname grad phi nbsp ergibt sich aus dem Gradient des Geschwindigkeitspotentials Diese Gleichungen sind bei sub trans und supersonischer Stromung anwendbar so lange ihre Rotationsfreiheit gegeben ist 5 BeweisIn einer stationaren Stromung ohne Beschleunigungsfeld lautet die Eulergleichung und die Massenbilanz grad v v 1 r grad p c 2 r grad r div r v grad r v r div v 0 displaystyle begin aligned operatorname grad vec v cdot vec v amp frac 1 rho operatorname grad p frac c 2 rho operatorname grad rho operatorname div rho vec v amp operatorname grad rho cdot vec v rho operatorname div vec v 0 end aligned nbsp Hier tritt die Schallgeschwindigkeit c p r S displaystyle c sqrt left frac partial p partial rho right S nbsp auf die die Wurzel aus der Ableitung des Drucks nach der Dichte bei konstanter Entropie S ist Skalare Multiplikation der Euler Gleichung mit der Geschwindigkeit und Einsetzen der Massenbilanz liefert v grad v v c 2 r grad r v c 2 r r div v c 2 div v displaystyle vec v cdot operatorname grad vec v cdot vec v frac c 2 rho operatorname grad rho cdot vec v frac c 2 rho rho operatorname div vec v c 2 operatorname div vec v nbsp Division durch c 2 displaystyle c 2 nbsp ergibt mit dem Geschwindigkeitspotential und der Machzahl M v c grad ϕ c displaystyle vec M vec v c operatorname grad phi c nbsp M g r a d g r a d ϕ M d i v g r a d ϕ D ϕ i 1 3 2 ϕ x i 2 i j 1 3 M i M j 2 ϕ x i x j 0 displaystyle begin aligned vec M cdot operatorname grad grad phi cdot vec M operatorname div grad phi amp Delta phi Leftrightarrow quad sum i 1 3 frac partial 2 phi partial x i 2 sum i j 1 3 M i M j frac partial 2 phi partial x i partial x j amp 0 end aligned nbsp Im sub oder supersonischen Bereich aber nicht im transsonischen kann bei kleinen Anstromwinkeln von schlanken Korpern eine weitere Annahme getroffen werden Das Geschwindigkeitspotential wird in eine Parallelstromung v x displaystyle v infty x nbsp und ein kleines Storungsfeld f displaystyle varphi nbsp aufgeteilt ϕ v x f displaystyle phi v infty x varphi nbsp In diesem Fall kann die linearisierte Potentialgleichung eine Naherung der kompletten Gleichung benutzt werden 5 1 M x 2 2 f x 2 2 f y 2 2 f z 2 0 displaystyle 1 M x 2 frac partial 2 varphi partial x 2 frac partial 2 varphi partial y 2 frac partial 2 varphi partial z 2 0 nbsp Schallwellen Bearbeiten Schallwellen mit geringer Amplitude konnen mit dem folgenden Potentialstromungsmodell approximiert werden 6 2 ϕ t 2 c 2 D ϕ displaystyle frac partial 2 phi partial t 2 bar c 2 Delta phi nbsp Darin ist c displaystyle bar c nbsp die mittlere Schallgeschwindigkeit im homogenen Fluid Diese Gleichung ist eine Wellengleichung die fur das Geschwindigkeitspotential F und gleichfalls fur den Druck sowie die Dichte gilt Randbedingungen Bearbeiten Die Dirichlet Randbedingungen geben den Wert des Geschwindigkeitspotentials auf einer Flache vor woraus sich die Tangentialgeschwindigkeit des Fluids auf der Flache ergibt Mit Neumann Randbedingungen werden die Ableitungen des Potentials senkrecht zu Flachen festgelegt was hier der Vorgabe von Geschwindigkeitskomponenten senkrecht zur Flache entspricht Ansonsten gilt das bei den Euler schen Gleichungen zu den Randbedingungen gesagte Zirkulation BearbeitenIn einem stationaren Stromungsgebiet in dem die Bernoulli Gleichung auf allen Stromlinien denselben Wert hat ist die Zirkulation im inneren des Gebiets gleich Null und die Stromung eine Potentialstromung Trotzdem konnen bei einer aus der Ruhe heraus entstehenden Bewegung einer homogenen reibungslosen Flussigkeit Wirbel entstehen Denn alle im Ruhezustand im Innern der Flussigkeit gezogenen materiellen Linien weichen bei einsetzender Stromung einer sich bildenden Trennflache aus und schneiden diese niemals Der Kelvin sche Wirbelsatz bezieht sich auf materielle Linien und sagt demnach uber die Beziehungen der Gebiete auf beiden Seiten der Trennflache nichts aus Deshalb konnen aus der Ruhe heraus an Kanten Trennflachen und Wirbel ohne Widerspruch zum Wirbelstatz entstehen 7 Da auch der Potentialwirbel eine Potentialstromung ist eine von null verschiedene Zirkulation aufweist und inkompressible Potentialstromungen superponierbar sind kann mit einem Potentialwirbel in oder ausserhalb des Stromungsgebietes eine Zirkulation eingebracht werden Allgemeiner kann ein divergenzfreier Geschwindigkeitsanteil superponiert werden der sich aus der Rotation eines Vektorfeldes ergibt dessen Berechnung die Kenntnis der Rotationsverteilung erfordert und die Losung zusatzlich erschwert 8 Potentialstromung eines inkompressiblen Fluids BearbeitenDie Annahme der Inkompressibilitat ist fur Flussigkeiten bei moderaten Drucken und fur Gasstromungen weit unterhalb der Schallgeschwindigkeit eine haufig sinnvolle Naherung In diesem Fall entkoppeln sich die Bestimmungsgleichungen fur die Geschwindigkeit und den Druck Aus der Laplace Gleichung bestimmt sich mit den Randbedingungen das Potential und aus den Euler Gleichungen die sich auf die Bernoulli sche Druckgleichung reduzieren der Druck Bestimmungsgleichungen Bearbeiten Bei Inkompressibilitat ist das Geschwindigkeitsfeld divergenzfrei div v d i v g r a d ϕ D ϕ 0 displaystyle operatorname div vec v operatorname div grad phi Delta phi 0 nbsp weswegen das Potential die Laplace Gleichung erfullt Auf Grund der Massenbilanz r r D ϕ r 0 displaystyle dot rho rho Delta phi dot rho 0 nbsp ist die Dichte zumindest zeitlich konstant bei Inkompressibilitat ist sie auch raumlich konstant Bei konstanter Dichte vereinfacht sich die bei den barotropen Fluiden eingefuhrte Druckfunktion P zu P p r Multiplikation der Euler Gleichungen mit der Dichte liefert dann r ϕ t r 2 grad ϕ 2 p r V E displaystyle rho frac partial phi partial t frac rho 2 operatorname grad phi 2 p rho V E nbsp In einer stationaren Stromung entfallt der erste Term und es verbleibt die Bernoulli sche Druckgleichung die hier im gesamten Stromungsgebiet gilt Die Konstante Energie E ist die Summe aus kinetischer innerer Druck p und ausserer Energie Schwerefeld rV eines Fluidpartikels und diese Gesamtenergie ist fur alle Partikel in der stationaren Stromung gleich Analogon der Warmeleitung Bearbeiten Hauptartikel Warmeleitung nbsp Vergleich der uber Warmeleitung berechneten Stromung farbig und schwarz mit den analytischen Stromlinien weiss nbsp Das Betragsquadrat des Warmestromes entspricht der DruckverteilungEine Potentialstromung kann mit Finite Elemente Programmen berechnet werden in denen die Warmeleitungsgleichung implementiert ist was in vielen gangigen Programmen der Fall ist Die stationare Warmeleitung in einem homogenen Festkorper ohne innere Warmequellen gehorcht der Laplace Gleichung D T 0 displaystyle Delta T 0 nbsp weshalb die Temperatur T dort derselben Bestimmungsgleichung unterliegt wie das Geschwindigkeitspotential F in einem inkompressiblen Fluid Wird die Temperatur in einem Festkorper als Geschwindigkeitspotential aufgefasst dann entspricht der Temperaturgradient zu dem der Warmestrom proportional ist der Stromungsgeschwindigkeit Die Dirichlet Randbedingungen legen die Temperaturen auf Oberflachen fest wodurch sich der Temperaturgradient und mithin die tangentiale Stromungsgeschwindigkeit in der Flache vorgeben lasst Neumann Randbedingungen bestimmen die Warmestromdichte die in der Potentialstromung die Normalgeschwindigkeit zur Wand einstellt Das obere Bild zeigt einen Ausschnitt der Temperaturverteilung aus einer Warmeleitungsrechnung mit der Finite Elemente Methode in einer quadratischen Scheibe mit einem mittig ausgesparten Profil das im Beispiel unten behandelt wird Die Temperatur ist farbkodiert von rot nach blau rot hoch blau gering gezeichnet und die von oben nach unten verlaufenden schwarzen Linien sind ihre Hohenlinien Die von links nach rechts orientierten kurzen schwarzen Striche stellen Warmestrome dar Um Randeffekte gering zu halten ist die Scheibe zehnmal so gross wie das Profil denn die im Beispiel unten analytisch berechneten weiss gezeichneten Stromlinien gelten fur eine Stromung in der gesamten wandfreien Ebene An den Grenzen der Scheibe ausserhalb des Bildes oben und unten sowie am Profil wurden keine Randbedingungen definiert Am linken Rand wurde die Warmestromdichte von eins vorgegeben und am rechten Rand die Warmestromdichte von minus eins sowie mittig die Temperatur null festgelegt Mit einer konstanten Temperatur an der rechten Begrenzung ware ein senkrechter Ausfluss bestimmt worden Die Hohenlinien der Temperatur entsprechen den Hohenlinien des Geschwindigkeitspotentials und sind im Einklang mit der Theorie augenscheinlich senkrecht zu den weissen Stromlinien und die Warmestrome sind mit den Stromlinien vertraglich Das untere Bild zeigt das Betragsquadrat des Warmestromes rot gross blau klein von dem nach der Bernoulli Gleichung ohne Schwerebeschleunigung p E r 2 grad ϕ 2 displaystyle p E frac rho 2 operatorname grad phi 2 nbsp der Druck eine lineare Funktion ist Wegen des negativen Vorzeichens herrscht also in den blauen Bereichen ein hoher Druck und in den roten ein niedriger Ebene Potentialstromung Bearbeiten Wenn die Stromung in der Ebene stattfindet dann konnen die Eigenschaften komplexer Funktionen ausgenutzt werden Es wird ein komplexes Geschwindigkeitspotential definiert dessen Realteil das reelle Geschwindigkeitspotential und dessen Imaginarteil die Stromfunktion ist Beide Funktionen beschreiben dieselbe Stromung Deshalb gelten die aus der Stromfunktion abgeleiteten Aussagen Die Hohenlinien der Stromfunktion sind Stromlinien zwischen denen der Volumenstrom uberall gleich ist Komplexe Geschwindigkeitspotentiale konnen aus einfacheren zusammengesetzt und transformiert werden so dass mit einfachen Mitteln komplizierte Stromungen untersucht werden konnen Insbesondere kann auch die von der Stromung auf den Korper ausgeubte Kraft exakt berechnet werden Komplexes Geschwindigkeitspotential Bearbeiten Die Ebene wird als komplexe Zahlenebene aufgefasst in der der Wert des Potentials als Realteil einer holomorphen Funktion f dargestellt wird f z ϕ x y i ps x y mit z x i y C displaystyle f z phi x y mathrm i psi x y quad text mit quad z x mathrm i y in mathbb C nbsp Die Menge C displaystyle mathbb C nbsp enthalt alle komplexen Zahlen und i ist die imaginare Einheit Die Funktion f ist das komplexe Geschwindigkeitspotential aus dem sich die Geschwindigkeit uber die Ableitungen v x ϕ x v y ϕ y oder w v x i v y d f d z displaystyle v x frac partial phi partial x quad v y frac partial phi partial y quad text oder quad w v x mathrm i v y frac mathrm d f mathrm d z nbsp ergibt Hier ist w die komplexe Geschwindigkeit Der Imaginarteil des Potentials ist die Stromfunktion deren Hohenlinien Stromlinien sind Weil die Funktion f holomorph ist gelten die Cauchy Riemann schen Differentialgleichungen ϕ x ps y und ϕ y ps x displaystyle frac partial phi partial x frac partial psi partial y quad text und quad frac partial phi partial y frac partial psi partial x nbsp weswegen die Gradienten des reellen Geschwindigkeitspotentials und der Stromfunktion zueinander senkrecht sind grad ϕ grad ps 0 displaystyle operatorname grad phi cdot operatorname grad psi 0 nbsp Die Gradienten sind wiederum orthogonal zu den Hohenlinien so dass auch diese sich im rechten Winkel schneiden Weitere Differentiation offenbart 2 ϕ x 2 2 ps y x 2 ϕ y 2 D ϕ 0 2 ps y 2 2 ϕ x y 2 ps x 2 D ps 0 displaystyle begin aligned frac partial 2 phi partial x 2 amp frac partial 2 psi partial y partial x frac partial 2 phi partial y 2 quad Rightarrow quad Delta phi 0 frac partial 2 psi partial y 2 amp frac partial 2 phi partial x partial y frac partial 2 psi partial x 2 quad Rightarrow quad Delta psi 0 end aligned nbsp weswegen sowohl der real als auch der imaginare Teil des komplexen Geschwindigkeitspotentials die Laplace Gleichung in der Ebene erfullt Die obere Gleichung legt die Divergenzfreiheit und die untere die Rotationsfreiheit der Stromung fest Aus den Randbedingungen berechnet sich das Potential aus dem sich uber die Ableitungen das Geschwindigkeitsfeld berechnet und der Druck wird aus der Bernoulli schen Druckgleichung ermittelt Naherungsweise rechnerische Losung Bearbeiten Der Laplace Operator hat in der Ebene in einem regelmassigen Netz mit Maschenweite x y h displaystyle triangle x triangle y h nbsp den Differenzenquotient D ps 2 ps x 2 2 ps y 2 ps x h y ps x h y 2 ps x y h 2 ps x y h ps x y h 2 ps x y h 2 0 ps x y ps x h y ps x h y ps x y h ps x y h 4 displaystyle begin aligned Delta psi amp frac triangle 2 psi triangle x 2 frac triangle 2 psi triangle y 2 amp frac psi x h y psi x h y 2 psi x y h 2 amp frac psi x y h psi x y h 2 psi x y h 2 stackrel displaystyle 0 rightarrow psi x y amp frac psi x h y psi x h y psi x y h psi x y h 4 end aligned nbsp Die Differenzen wurden mit displaystyle triangle nbsp geschrieben um eine Verwechselung mit dem Laplace Operator zu vermeiden Die Formel gestattet eine einfache Losung der Laplace Gleichung 9 In das Stromungsgebiet wird ein Netz mit quadratischen Maschen mit Kantenlange h gelegt Der Wert der Stromfunktion am Knoten mit Koordinaten x y ist dann der Mittelwert der in x und y Richtung benachbarten Knotenwerte Bei einem durchflossenen Kanal wird auf einem Rand die Stromfunktion zu null gesetzt und auf dem gegenuberliegenden Rand derjenige Wert festgelegt der dem verlangten Volumenstrom entspricht Fur die Knotenwerte im Stromungsgebiet werden zunachst geschatzte Werte eingesetzt und dann iterativ mittels obiger Gleichung angepasst bis eine zufriedenstellende Genauigkeit erreicht ist Zeichnerische Ermittlung von Potentialstromungen Bearbeiten Hauptartikel Stromungsnetz nbsp In sechs Schritten manuell erstellte Zeichnung einer Potentialstromung um eine offene Hubschutze mittels XfigPotentialstromungen konnen auch von Hand skizziert werden was in der Geotechnik und Hydromechanik angewendet wird 9 Die oben aufgefuhrten Eigenschaften von ebenen Potentialstromungen geben dazu Anhaltspunkte Undurchlassige Rander oder freie Oberflachen sind Stromlinien Stromlinien durfen sich nicht schneiden Potentiallinien durfen sich nicht schneiden Potentiallinien und Stromlinien schneiden sich rechtwinklig und das Netz aus Potential und Stromlinien kann so gestaltet werden dass es aus annahernd quadratischen Maschen besteht Zu Anfang werden die Rander des Stromungsgebiets gezeichnet Schritt 1 im Bild und die Richtung der Stromung am Ein und Auslass festgelegt 2 Dann werden 3 zwischen den Randern des Stromungsgebiets einige wenige parallele Stromlinien mit gleichen Abstanden gelegt unter Beachtung der Tatsache dass die Rander ebenfalls Stromlinien sind Die Potentiallinien werden so eingezeichnet 4 dass sie diese Stromlinien im rechten Winkel kreuzen und quadratische Maschen entstehen Durch Eintragung weiterer Strom und Potentiallinien 5 6 wird das Netz bis zum gewunschten Mass verdichtet Konstruktion von Geschwindigkeitspotentialen Bearbeiten Weil die Laplace Gleichung linear ist ist die Stromung die sich aus der Summe zweier Geschwindigkeitspotentiale ergibt wieder eine Potentialstromung So lassen sich komplexe Stromungen durch Superposition aus einfachen Stromungen zusammensetzen von denen einige siehe die unten stehenden Bilder angegeben seien Eine Parallelstromung mit konstanter komplexer Geschwindigkeit w v x i v y displaystyle w v x mathrm i v y nbsp ergibt sich aus f z w z displaystyle f z wz nbsp Eine Staupunktstromung mit Staupunkt im Ursprung ergibt sich mit f z z n n 2 3 4 displaystyle f z z n n 2 3 4 dots nbsp Eine Multipolstromung im Ursprung hat das Potential f z z n n 1 2 3 4 displaystyle f z z n n 1 2 3 4 dots nbsp Quellen haben ein Potential der Form f z Q 2 p ln z z 0 displaystyle f z frac Q 2 pi ln z z 0 nbsp Der Ort der Quelle ist z 0 displaystyle z 0 nbsp ihre Starke ist Q und ln ist der Naturliche Logarithmus An der Stelle z z 0 displaystyle z z 0 nbsp befindet sich eine Singularitat in der die Laplace Gleichung verletzt wird Senken sind Quellen mit negativer Starke Der Potentialwirbel siehe unten ergibt sich durch Vertauschung der Real und Imaginarteile bei einer Quellstromung was durch Multiplikation des Potentials mit i geschieht f z i G 2 p ln z z 0 displaystyle f z mathrm i frac Gamma 2 pi ln z z 0 nbsp Strudel ergeben sich aus der Uberlagerung von Quellen Senken und Potentialwirbel f z Q i G 2 p ln z z 0 displaystyle f z frac Q mathrm i Gamma 2 pi ln z z 0 nbsp nbsp Parallelstromung nbsp Staupunktstromung nbsp Multipolstromung nbsp Quelle Senke nbsp Wirbel nbsp StrudelIn den obigen Bildern sind die Hohenlinien des reellen Potentials rot und die der Stromfunktion blau gezeichnet Der Abstand der roten Linien gibt einen Eindruck von der Stromungsgeschwindigkeit wobei in Bereichen mit geringen Abstanden die Geschwindigkeit hoch ist Die blauen Linien sind Stromlinien Die Funktionen deren Hohenlinien bei den Quellen Senken Wirbeln oder Strudeln radial verlaufen machen irgendwo einen Sprung was eine Folge der Nicht Differenzierbarkeit der Logarithmusfunktion ist Mit der Methode der Bildladungen kann durch geschicktes Einbringen von gedachten Quellen und Senken ausserhalb des durchstromten Gebiets die Stromung so umgelenkt werden dass sie vorgegebene Randbedingungen erfullt Mit konformen Abbildungen konnen Stromungsfelder um einfache Grundkorper auf komplizierte Geometrien ubertragen werden Die Ubertragung geschieht mit einer zweiten komplexen Funktion w gemass P z w f z displaystyle Pi z w f z nbsp die nach dem Riemannschen Abbildungssatz konform ist Die Kutta Schukowski Transformation ubertragt die Umstromung eines Kreiszylinders auf ein Flugelprofil siehe das Beispiel unten Mit der Schwarz Christoffel Transformation kann die Parallelstromung in der oberen Halbebene auf beliebige mit geraden Linien berandete einfach zusammenhangende Gebiete und auch das innere von Polygonen ubertragen werden Kraftwirkungen auf umstromte Korper Bearbeiten Das d Alembert sche Paradoxon besagt dass auf einen beliebig geformten Korper keine Kraft in Richtung der Stromung wirkt In einem komplex differenzierbaren Geschwindigkeitsfeld ubt die Stromung auch keine Kraft senkrecht zur Stromungsrichtung auf von ihr umstromte Korper aus was eine Folgerung aus dem Integralsatz von Cauchy ist Wenn das Geschwindigkeitspotential irgendwo innerhalb der Kontur eines Korpers nicht komplex differenzierbar ist dann kann die Zirkulation der Geschwindigkeit langs der Kontur ungleich null sein und der Korper erfahrt nach dem Satz von Kutta Joukowski eine zu dieser Zirkulation proportionale Auftriebskraft Mit der 1902 entdeckten Formel konnten erste auftriebserzeugende Flugelprofile entwickelt werden Eine bei der Berechnung der Kraftwirkung auf umstromte Korper wichtige Grosse ist die Zirkulation G der Geschwindigkeit langs eines Weges W die mit einem Kurvenintegral berechnet wird G W v x d x v y d y displaystyle Gamma oint W v x mathrm d x v y mathrm d y nbsp Wenn die Kurve eine umstromte Kontur ist dann kann die Zirkulation mit der komplexen Geschwindigkeit berechnet werden G W w z d z R displaystyle Gamma oint W w z mathrm d z quad in mathbb R nbsp Der Integralsatz von Cauchy besagt dass das Kurvenintegral einer komplexen Funktion zwischen zwei Punkten wegunabhangig ist wenn die Funktion holomorph also komplex differenzierbar ist Das Kurvenintegral verschwindet demnach entlang einer geschlossenen Kontur wenn die Funktion in dem von der Kontur umschlossenen Gebiet holomorph ist Die Zirkulation kann nur dann von null verschieden sein wenn das Geschwindigkeitsfeld irgendwo innerhalb der Kontur nicht komplex differenzierbar ist 10 Die komplexe Kraft F F x i F y displaystyle F F x mathrm i F y nbsp wobei F y displaystyle F y nbsp in positiver y Richtung wirkt die auf einen von einer ebenen Potentialstromung umstromten Korper wirkt dessen konstanter Querschnitt die ebene Kontur W ist und dessen Ausdehnung senkrecht zum Querschnitt gleich L ist berechnet sich mit der ersten Blasius schen Formel F i L r 2 W w 2 z d z displaystyle F mathrm i L frac rho 2 oint W w 2 z mathrm d z nbsp mit dem Kurvenintegral des Geschwindigkeitsquadrats entlang der Kontur Der Satz von Kutta Joukowski besagt dass die Kraft die auf den umstromten Korper wirkt proportional zu seiner Zirkulation ist F i L r G w L r G e i p 2 w displaystyle F mathrm i L rho Gamma w infty L rho Gamma e mathrm i frac pi 2 w infty nbsp Weil die Kraft immer senkrecht zur Stromungsgeschwindigkeit w displaystyle w infty nbsp im Unendlichen ist wird die Kraft auch Auftriebskraft genannt In wirbelfreien Stromungen verschwindet die Zirkulation so dass eine wirbelfreie Potentialstromung keine Kraft auf von ihr umstromte Korper ausubt Potentialwirbel Bearbeiten Der Potentialwirbel oder freie Wirbel ist eine echte rotationsfreie Potentialstromung die dennoch kreist d h in einem topologisch zweifach zusammenhangenden Gebiet wie beispielsweise dem Luftraum in einem Saal mit Mittelsaule eine Zirkulation aufweist Ein besonderer Potentialwirbel ist an einer freien Wasseroberflache zu beobachten wenn der Druck im Zentrum so gering wird dass sich die Oberflache merklich einsenkt und einen Wirbeltrichter Strudel bildet Reicht der Trichter unbegrenzt weit in die Tiefe so herrscht im ganzen Flussigkeitsgebiet Potentialstromung nicht jedoch im luftgefullten Kern Beim freien Wirbel bewegen sich alle Fluidpartikel auf konzentrischen Kreisbahnen mit Geschwindigkeits betragen v r displaystyle v r nbsp die ausser im Kerngebiet dem Abstandsgesetz mit Konstante c displaystyle c nbsp und Achsenabstand r displaystyle r nbsp entsprechen v r c r displaystyle v r frac c r nbsp Es handelt sich um eine reziproke Proportionalitat d h je weiter entfernt von der Achse die Partikel sind desto langsamer bewegen sie sich und die innersten Partikel bewegen sich am schnellsten Dadurch ist eine vollig andere Geschwindigkeits und Druckverteilung gegeben als beim quasi starren Wirbel der kein Potential besitzt und dessen Geschwindigkeit proportional zum Achsenabstand ist v w r displaystyle v omega cdot r nbsp In realen Fluiden bilden sich freie Wirbel nur naherungsweise als Potentialstromungen aus da in ihrer Mitte Zahigkeitskrafte zu einer quasi starren Rotation fuhren und das Geschwindigkeitsfeld hier Vortizitat besitzt Nach Bestehorn ist der Potentialwirbel eigentlich ein Punktwirbel 11 Das komplexe Potential des Wirbels lautet mit z z 0 r e i f displaystyle z z 0 re mathrm i varphi nbsp wie oben angedeutet f z i G 2 p ln z z 0 i G 2 p ln r e i f G 2 p f i G 2 p ln r displaystyle f z mathrm i frac Gamma 2 pi ln z z 0 mathrm i frac Gamma 2 pi ln left re mathrm i varphi right frac Gamma 2 pi varphi mathrm i frac Gamma 2 pi ln r nbsp Auf der negativen reellen Achse bei f n p n 1 3 5 displaystyle varphi pm n pi n 1 3 5 ldots nbsp ist der komplexe Naturliche Logarithmus nicht differenzierbar weswegen das Potential gar keines ist Wohl benutzbar ist aber die Stromfunktion ps G 2 p ln r displaystyle psi frac Gamma 2 pi ln r nbsp die den Wirbel vollstandig beschreibt Die Wirbelstarke ist w 1 2 D ps 1 2 G d z z 0 displaystyle omega frac 1 2 Delta psi frac 1 2 Gamma delta z z 0 nbsp mit dem Dirac Delta d weshalb es sich beim Potentialwirbel um einen unendlich starken Punktwirbel an der Stelle z z 0 displaystyle z z 0 nbsp handelt siehe Hamel Oseenscher Wirbel Beispiel Bearbeiten nbsp Ebene Potentialstromung um einen Zylinder und ein FlugelprofilDas komplexe Geschwindigkeitspotential f z z 1 z x i y 1 x i y x x x 2 y 2 ϕ i y y x 2 y 2 ps displaystyle begin aligned f z amp z frac 1 z x mathrm i y frac 1 x mathrm i y amp underbrace x frac x x 2 y 2 phi mathrm i underbrace left y frac y x 2 y 2 right psi end aligned nbsp beschreibt die Umstromung eines Zylinders denn dessen Umfang bei x 2 y 2 1 displaystyle x 2 y 2 1 nbsp ist eine Stromlinie mit ps 0 siehe den oberen Bildteil Die Stromlinien sind Hohenlinien der Stromfunktion ps und haben die Gleichung y y x 2 y 2 ps x y y ps y 2 displaystyle y frac y x 2 y 2 psi quad Rightarrow quad x pm sqrt frac y y psi y 2 nbsp und sind im Bild im Abstand von Dps 0 3 weiss gezeichnet Die Stromlinie mit ps 0 ist rot gezeichnet Die Geschwindigkeit ergibt sich durch Ableitung des Geschwindigkeitspotentials v x i v y d f d z 1 1 z 2 v 2 1 2 y 2 x 2 1 x 2 y 2 2 displaystyle begin aligned v x mathrm i v y amp frac mathrm d f mathrm d z 1 frac 1 z 2 Rightarrow quad v 2 amp 1 frac 2 y 2 x 2 1 x 2 y 2 2 end aligned nbsp In Abwesenheit eines Schwerefeldes liefert die Bernoulli Gleichung den Druck r 2 v 2 p E p E r 2 1 2 y 2 x 2 1 x 2 y 2 2 displaystyle frac rho 2 v 2 p E quad Rightarrow quad p E frac rho 2 left 1 frac 2 y 2 x 2 1 x 2 y 2 2 right nbsp Mit der konformen Abbildung w z z z 0 R 2 z z 0 mit R ℜ z 0 1 ℑ z 0 2 displaystyle w z z z 0 frac R 2 z z 0 quad text mit quad R Re z 0 sqrt 1 Im z 0 2 nbsp wird der Zylinder auf eine Kontur transformiert die einem Flugelprofil ahnelt 12 siehe den unteren Bildteil wo z 0 0 15 e i 5 6 p displaystyle z 0 0 15 e mathrm i frac 5 6 pi nbsp und R 0 867 2 displaystyle R 0 8672 nbsp ist Die Funktionen ℜ displaystyle Re nbsp und ℑ displaystyle Im nbsp extrahieren den Real bzw Imaginarteil ihres komplexen Arguments Der Betrag von z 0 displaystyle z 0 nbsp kontrolliert die Bauchigkeit wohingegen die Phase die Wolbung des Profils beeinflusst Die Stromung wird wie die Kontur transformiert so dass sich die Stromlinien in der w Ebene mit w z z z 0 R 2 z z 0 und z y y ps y 2 i y displaystyle w z z z 0 frac R 2 z z 0 quad text und quad z pm sqrt frac y y psi y 2 mathrm i y nbsp ergeben Sie sind im unteren Bildteil weiss gezeichnet Die Stromlinien dieses Profils wurden oben im Abschnitt Analogon der Warmeleitung mit den Ergebnissen aus einer Warmeleitungsanalyse verglichen Weblinks BearbeitenK Luders R O Pohl G Beuermann K Samwer Stromlinien Modellversuche IWF Wissen und Medien 2004 abgerufen am 3 Januar 2021 Einzelnachweise Bearbeiten Oertel 2012 S 75 grad f f 1 2 grad f f f rot f displaystyle operatorname grad vec f cdot vec f frac 1 2 operatorname grad vec f cdot vec f vec f times operatorname rot vec f nbsp Spurk 2010 S 118 Greve 2003 S 147 a b c John D Anderson Modern compressible flow With historical perspective McGraw Hill New York NY 2002 ISBN 0 07 242443 5 S 358 359 Horace Lamb Hydrodynamics 6th edition Cambridge University Press Cambridge 1993 ISBN 0 521 05515 6 S 492 495 Oertel 2012 S 75 f Spurk 2010 S 353 a b A Malcherek Hydromechanik fur Bauingenieure pdf Universitat der Bundeswehr Munchen S 48ff abgerufen am 9 Oktober 2016 deutsch Bestehorn 2006 S 95 Bestehorn 2006 S 87 Bestehorn 2006 S 91 ff Literatur BearbeitenJ H Spurk Stromungslehre Einfuhrung in die Theorie der Stromungen 8 uberarbeitete Auflage Springer Verlag Heidelberg Dordrecht London New York 2010 ISBN 978 3 642 13142 4 S 343 ff doi 10 1007 978 3 642 13143 1 eingeschrankte Vorschau in der Google Buchsuche abgerufen am 2 Januar 2021 Michael Bestehorn Hydrodynamik und Strukturbildung Springer Berlin u a 2006 ISBN 3 540 33796 2 Ralf Greve Kontinuumsmechanik Ein Grundkurs fur Ingenieure und Physiker Springer Berlin u a 2003 ISBN 3 540 00760 1 Hans J Lugt Wirbelstromung in Natur und Technik G Braun Karlsruhe 1979 ISBN 3 7650 2028 1 Herbert Oertel Hrsg Prandtl Fuhrer durch die Stromungslehre Grundlagen und Phanomene 11 uberarbeitete und erweiterte Auflage Vieweg Braunschweig u a 2002 ISBN 3 528 48209 5 Heinz Schade Ewald Kunz Stromungslehre 2 durchgesehene und verbesserte Auflage de Gruyter Berlin u a 1989 ISBN 3 11 011873 4 Jurgen Zierep Grundzuge der Stromungslehre 4 uberarbeitete Auflage G Braun Karlsruhe 1990 ISBN 3 7650 2039 7 Abgerufen von https de wikipedia org w index php title Potentialstromung amp oldid 229067764