Strömungen

Beschreibung von Strömungen

Vektorfeld der Strömung

An jedem Punkt hat die Geschwindigkeit $\vec{v}\left( \vec{r}\right)$ einen Betrag und eine Richtung.

Das Vektorfeld $\vec{v}\left( \vec{r}\right)$ ist durch Stromlinien charakterisiert. Wenn $\vec{v}\left( \vec{r}\right)$ nicht von der Zeit abhängt, heisst die Strömung stationär.

Bahnlinien: Bahn eines Teilchens

Bei stationären Strömungen sind Bahnlinien und Stromlinien identisch. Inkompressible Strömungen sind Strömungen mit konstanter Dichte $\rho$.

Fluss

Fluss

Der Fluss ist definiert als

$$d\phi=\rho v\cos\alpha dA$$ (3.220)

oder

$$d\phi=\rho\vec{v\cdot}d\vec{A}$$ (3.221)

Integralform

$$\phi=\int\limits_{A}\rho\vec{v}d\vec{A=}\int\limits_{A}\vec{j}d\vec{A}$$ (3.222)

wobei $A$ beliebige Fläche (auch gekrümmt) ist. $\phi$ ist der Fluss. Er hat die Einheit

$$\left[\phi\right] = \frac{kg}{s}$$

$\vec{j}=\rho\vec{v}$ ist die Stromdichte (analog zum elektrischen Strom). Ihre Einheit ist

$$\left[\vec{j}\right] = \frac{kg}{m^2 s}$$

Bei einer geschlossenen Fläche fliesst netto Masse aus dem Volumen heraus, wenn eine Quelle im Volumen ist. Bei geschlossenen Flächen wird der Fluss nach aussen immer positiv gezählt, nach innen negativ.

Berechnung der Divergenz

Im allgemeinen Fall ändert die Dichte mit dem Ort. Wir haben links die Dichte $\rho(x)$ und die Geschwindigkeit $v_x(x)$ und rechts $\rho(x+dx) = \rho(x)+\frac{\partial \rho}{\partial x} dx$ und $v_x(x+dx) = v_x(x)+\frac{\partial v_x}{\partial x} dx$. Zusammen erhalten wir

$$d\phi_{1}\left( x\text{,} y\text{,} z\right) =-\rho\left( x\text{,} y\text{,} z\right) v_{x}\left( x\text{,} y\text{,} z\right)$$

$$d\phi_{2}\left( x\text{,} y\text{,} z\right) =\rho\left( x+dx\text{,} y\text{,} z\right) v_{x}\left( x+dx\text{,} y \text{,} z\right) dydz$$

$$= \left(\rho\left( x\right) v_{x}\left( x\right) + \left[\rho\lef... ...l\rho}{\partial x}v_x\left( x\text{,} y\text{,} z\right)\right]dx\right) dydz$$

wobei die Zweitordnungsterme mit $dx^2$ vernachlässigt wurden. Der Nettofluss ist

$$d\phi_{1}\left( x\text{,} y\text{,} z\right)+d\phi_{2}\left( x\text{,} y\text{,} z\right) =\left[\rho\left( x\text{,} y\text{,} z\right)\frac{\partial v_... ...{\partial\rho}{\partial x}v_x\left( x\text{,} y\text{,} z\right)\right]dxdydz$$

$$=\frac{\partial}{\partial x}\left[\rho\left( x\text{,} y\text{,} z\right)v_x\left( x\text{,} y\text{,} z\right)\right]dxdydz$$

$$=\frac{\partial}{\partial x}\left[\rho\left( x\text{,} y\text{,} z\right)v_x\left( x\text{,} y\text{,} z\right)\right]dV = d\phi_{x}\left( x\text{,} y\text{,} z\right)$$

Analog bekommt man für die $y$- und die $z$-Richtung:

$$d\phi_{y}\left( x\text{,} y\text{,} z\right) =\frac{\partial}{\partial x}\left[\rho\left( x\text{,} y\text{,} z\right)v_x\left( x\text{,} y\text{,} z\right)\right]dV$$

$$d\phi z\left( x\text{,} y\text{,} z\right) =\frac{\partial}{\partial x}\left[\rho\left( x\text{,} y\text{,} z\right)v_x\left( x\text{,} y\text{,} z\right)\right]dV$$ (3.223)

Der Nettofluss aus dem Volumen ergibt sich zu

$$d\phi\left( x\text{,} y\text{,} z\right) =d\phi_{x}\left( x\text{,} y\text{,} z\right)+d\phi_{y}\left( x\text{,} y\text{,} z\right)+d\phi_{z}\left( x\text{,} y\text{,} z\right)$$

$$=\left( \frac{\partial \left(\rho v_{x}\right)}{\partial x}+\frac... ...ght)}{\partial y}+\frac{\partial \left(\rho v_{z}\right)}{\partial z}\right) dV$$ (3.224)

Ohne Quelle ist $d\phi=0$. Die Grösse

$$\textrm{div} {}\left(\rho \vec{v}\right) =\frac{\partial \left(\... ...o v_{y}\right)}{\partial y}+\frac{\partial \left(\rho v_{z}\right)}{\partial z}$$ (3.225)

(Divergenz) beschreibt die Quellen und Senken in einem Fluss. Sie ist in dem Falle auch gleich null.

Wenn $\textrm{div} {}\left( \vec{v}(\vec{r})\right) \neq 0$ ist, so muss sich die Dichte an der Stelle $\vec{r}$ ändern, oder es muss eine Quelle vorhanden sein.

Wir erinnern uns, der Fluss $\phi$ beschreibt die Massenänderung pro Zeit im Volumen $dV$. Die Masse im Volumen $dV$ nimmt um $d\phi dt$ ab (bei positivem $d\phi$ und nimmt um $\dot{\rho}dV dt$ zu. Zusammen ergibt sich

$$dm = -d\phi dt = -\textrm{div} {}{\left(\rho \vec{v}\right)} dV dt = \dot{\rho} dV dt$$

und damit

$$d\phi= \textrm{div} {}\left(\rho\vec{v}\right)dV =- \dot{\rho}dV$$ (3.226)

Damit bekommt man

$$\textrm{div} {}\left(\rho \vec{v}\right)=-\dot{\rho}$$ (3.227)

Dies ist die Kontinuitätsgleichung.

Wenn $\rho$ nicht vom Ort abhängt, hat man auch.

$$\rho \textrm{div} {}\vec{v}=-\dot{\rho}$$

Dann heisst $\rho \textrm{div} {}\vec{v}=\textrm{div} {} \vec{j}$ die Quelldichte.

Eine quellenfreie inkompressible Strömung hat überall $\textrm{div} {}\vec{v}=0.$

Es gilt:

$$\phi=\underset{\text{Fluss durch A (\glqq Materialmenge\grqq )}}{... ...\text{{\uml A}nderung der Dichte}}{\underbrace{\iiint\limits_{V}\dot{\rho}dV} }$$ (3.228)

da der Satz von Gauss für ein beliebiges Vektorfeld $\vec{x}$ besagt

$$\iint\limits_{A}{}\vec{x}\cdot d\vec{A}=\iiint\limits_{V}(\textrm{div} {}\vec{x})dV$$ (3.229)

Stromlinien in einer inkompressiblen Flüssigkeit

Die Stromlinien durch $A$ definieren einen Schlauch, die Stromröhre, die keinen Austausch mit der Umgebung hat. Also ist in einer inkompressiblen Flüssigkeit

$$A_{1}v_{1}=A_{2}v_{2}$$ (3.230)

Dies ist die makroskopische Kontinuitätsgleichung.

Geschwindigkeitsgradient und Rotation

Wenn die Strömung inhomogen ist, werden mitgeführte Teilchen gedreht.

Wasser strömt mit $\pm\ell\frac{dv_{y}}{dx}$ am Würfel vorbei. Wir haben ein Gleichgewicht, wenn

$$\omega_{z}=\frac{\partial v_{y}}{\partial x}-\frac{\partial v_{x}}{\partial y}$$ (3.231)

$\omega$ zeigt in die $z$ Richtung. Im Allgemeinen ist

$$\vec{\omega}=\left( \frac{\partial v_{z}}{\partial y}-\frac{\part... ...l y},\frac{\partial v_{y}}{\partial x}-\frac{\partial v_{x}}{\partial y}\right)$$ (3.232)

mit $\vec{\nabla=}\left( \frac{\partial}{\partial x},\frac{\partial }{\partial y},\frac{\partial}{\partial z}\right)$ wird

$\textrm{rot} {} \vec{v=\nabla\times v}$ die Rotation des Strömungsfeldes.

Es gilt dann

$$\underset{\text{Bahnkurve s}}{\underbrace{\oint\vec{v}d\vec{s}} }... ... berandete Fl{\uml a}che}}{\underbrace{\int \textrm{rot} {} \vec{v}d\vec{A}}}$$ (3.233)

Falls $\textrm{rot} {} \vec{v}=0$ ist kann $\vec{v}$ aus dem Geschwindigkeitspotential $U$ abgeleitet werden.

$$\vec{v}=-\textrm{grad} {}U$$ (3.234)

Dann gilt $\textrm{rot} {} \vec{v}=0$

Für inkompressible Flüssigkeiten gilt

$$\textrm{div} {}\vec{v}=-\textrm{div} {}\textrm{grad} {}U=-\Delta U=0$$ (3.235)

Wir haben also drei unterschiedliche physikalische Phänomene, die durch die gleiche Mathematik beschrieben werden:

Strömung $\longleftrightarrow$ Graviation $\longleftrightarrow$ Elektrostatik

Lokale und totale Ableitungen

Mitbewegtes System

Sei $S^{\ast}$ das Laborsystem, $S$ das mitbewegte System,das $\Delta m$ folgt. Seien $s^{\ast}$ lokale zeitliche Ableitungen und $s$ totale zeitliche Ableitungen

Das 2. Newtonsches Gesetz beschreibt die Bewegung von $\Delta m,$ aber nur in $S\left( t\right)$

(in $S^{\ast}$ betrachtet man Volumina, nicht Massen)

Also ist

$$\vec{F}_{m}=\frac{\vec{F}}{m}=\vec{a}=\frac{d\vec{v}} {dt} in S\left( t\right)$$ (3.236)

Lokale Ableitung:

$\vec{r}$ ist in $S^{\ast}$ fest

$$\frac{\partial\rho}{\partial t} =\frac{\partial}{\partial t}\rho\left( \vec{r},t\right) =\frac{\partial}{\partial t}\rho\left( x,y,z,t\right)$$

$$\frac{\partial\vec{v}}{\partial t} =\frac{\partial}{\partial t}\vec{v}\left( \vec{r},t\right) =\frac{\partial}{\partial t}\vec{v}\left( x,y,z,t\right)$$ (3.237)

Totale zeitliche Ableitungen

In $S\left( t\right)$ beschreiben die physikalischen Grössen das gleiche Teilchen.

$$\vec{r}=\vec{r}\left( t\right) =\left( x\left( t\right) ,y\left( t\right) ,z\left( t\right) \right)$$ (3.238)

wobei $x\left( t\right)$ die Koordinaten in $S^{\ast}$ sind

$$\rho =\rho\left( \vec{r}\left( t\right) ,t\right) =\rho\left( x\left( t\right) ,y\left( t\right) ,z\left( t\right) t\right)$$

$$\vec{v} =\vec{v}\left( \vec{r}\left( t\right) ,t\right) =\vec{v}\left( x\left( t\right) ,y\left( t\right) ,z\left( t\right) t\right)$$ (3.239)

Ableitung: totale Ableitung auf der Bahn $\vec{r}\left( t\right)$

$$\frac{d\rho}{dt} =\frac{d}{dt}\rho\left( x\left( t\right) ,y\left( t\right) ,z\left( t\right) ,t\right)$$

$$\vec{a} =\frac{d\vec{v}}{dt}$$ (3.240)

Zusammenhang:

Dichte

$$\frac{d}{dt}\rho=\frac{\partial}{\partial t}\rho+\left( \textrm{grad} {} \rho\right) \cdot\vec{v}$$ (3.241)

Beweis:

$$\frac{d}{dt}\rho=\frac{\partial\rho}{\partial x}\cdot\frac{\parti... ...{\partial z}\cdot\frac{\partial z}{\partial t} +\frac{\partial\rho}{\partial t}$$ (3.242)

qed.

Beschleunigung:

$$\vec{a}=\frac{d\vec{v}}{dt}=\frac{\partial\vec{v}}{\partial t}+\l... ... {}\left( \frac{1}{2}\vec{v}^{2}\right) -\vec{v}\times \textrm{rot} {}\vec{v}$$ (3.243)

Kontinuitätsgleichung

Grund Massenerhaltung

$$\frac{\partial\rho}{\partial t}+\textrm{div} {}\left( \rho\vec{v}\right) =\frac{\partial\rho}{\partial t}+\rho \textrm{div} {} \vec{v}=0$$ (3.244)

Beweis: Ortsfests Volumen

$$\Delta V =\Delta x\cdot\Delta y\cdot\Delta z$$

$$\delta\left( \Delta m\right) =\Delta x\cdot\Delta y\cdot\Delta z\cdot\frac{\partial}{\partial t}\rho\left( z+\frac{1}{2}\Delta z\right) \delta t$$

$$=-\rho\left( z+\Delta z\right) \cdot v_{z}\left( z+\Delta z\right... ...rho\left( z\right) v_{z}\left( z\right) \cdot\delta t\cdot\Delta x\cdot\Delta y$$

$$\frac{\delta\left( \Delta m\right) }{\partial t} =\Delta x\Delta y\Delta z\frac{\partial}{\partial t}\rho\left( z\... ...{z}\left( z\right) \right) }{\partial z}\cdot\Delta z\cdot\Delta x\cdot\Delta y$$

$$\frac{\partial\rho\left( z\right) }{\partial t} =-\frac{\partial }{\partial t}\left( \rho\left( z\right) v_{z}\left( z\right) \right) usw.$$ (3.245)

Stationäre Strömung

Sei $\frac{\partial}{\partial t}f=0$ für beliebige $f$. Dann ist die Dichte $\frac{dp} {dt}=\left( \textrm{grad} {} \rho\right) \cdot\vec{v}$ und die Kontinuitätsgleichung $\textrm{div} {}\left( \rho\vec{v}\right) =0$.

$$\vec{a}=\frac{d\vec{v}}{dt}=\textrm{grad} {}\left( \frac{1}{2}\vec{v}^{2}\right) -\vec{v}\times \textrm{rot} {}\vec{v}$$ (3.246)

Im Stromfaden gilt

$$A_{1}\rho_{1}v_{1}=A_{2}\rho_{2}v_{2}=const$$ (3.247)

Inkompressible Flüssigkeiten

$$\frac{d\rho}{dt} =\frac{\partial\rho}{\partial t}=0$$

$$\textrm{div} {} \vec{v} =0$$ (3.248)

dann gilt: $A_{1}v_{1}=A_{2}v_{2}=const.$

Innere Reibung

Innere Reibung in einer Flüssigkeit

Moleküle haben am Rand im Mittel die Geschwindigkeit der Wand.

$\vec{v}\left( z\right)$ ist parallel zur Wand.

Für die Kraft gilt

$$\vec{F}=\eta A\frac{\vec{v}}{z}$$ (3.249)

$\eta:\left[ \frac{Ns}{m^{2}}\right]$ heisst Viskosität (Scherviskosität)

Beispiel

Wasser: $1.8\cdot10^{-3}\frac{Ns}{m^{2}}$

Glyzerin $1\frac{Ns}{m^{2}}$

Allgemein:

$$\vec{F}=\eta A\frac{d\vec{v}}{dz}$$ (3.250)

wobei $A$ klein sein soll.

Temperaturabhängigkeit:

Moleküle müssen ihren Platz wechseln ( $\rightarrow$ Bolzmannstatistik)

$$\eta=\eta_{\infty}e^{\frac{b}{T}}$$ (3.251)

Laminare Strömung

Volumenkräfte

Laminare Stromlinien sind dadurch charakterisiert, dass benachbarte Stromlinien benachbart bleiben (Beispiel: Blut)

$$d\vec{F}_{1} =\left. -\eta\frac{\partial\vec{v}}{\partial x}\right\vert _{1}dydz$$

$$d\vec{F}_{2} =\eta\left. \frac{\partial\vec{v}}{\partial x}\right\vert _{2}dyd... ...ial x}\right\vert _{1}+\frac{\partial^{2}\vec{v}}{\partial x^{2}}dx\right) dydz$$

$$d\vec{F}_{R} =d\vec{F}_{1}+d\vec{F}_{2}=\eta\frac{\partial ^{2}\vec{v}}{\partial x^{2}}dxdydz=\eta\frac{\partial^{2}\vec{v}}{\partial x^{2}}dV$$ (3.252)

Allgemein:

$$d\vec{F}_{R}=\eta\left( \frac{\partial^{2}\vec{v}}{\partial x^{2}... ...v}}{\partial^{2}y^{2}}+\frac{\partial ^{2}\vec{v}}{\partial^{2}z^{2}}\right) dV$$ (3.253)

oder

$$\vec{F}_{V_{R}}=\frac{d\vec{F}}{dV}=\eta\Delta\vec{v}$$ (3.254)

die Volumenkraft der Reibung

Die Druckkraft ist

$$d\vec{F}_{p}=pdydz-\left( p+\frac{\partial p}{\partial x}dx\right) dydz=-\frac{\partial p}{\partial x}dV$$ (3.255)

also

$$\vec{F}_{V_{p}}=-\textrm{grad} {} p$$ (3.256)

$\vec{F}_{V_{R}}$ und $\vec{F}_{V_{p}}$ beschreiben die Dynamik

Strömung durch einen Spalt

Strömung durch einen Spalt

$v=0$ an der Wand

$v=v_{0}$ in der Mitte

$\frac{dv}{dx}\Rightarrow$ Reibungskraft $F_{R}=2\ell b\eta\frac{dv}{dx}$

Druck: $F_{p}=2xb\ell\frac{dp}{dz}$

$\Rightarrow\frac{dv}{dx}=\frac{1}{\eta}\frac{dp}{dz}x$

$\Rightarrow v\left( x\right)$ ist eine Parabel

$$v=v_{0}-\frac{1}{2\eta}\frac{dp}{dx}x^{2}=v_{0}-\frac{p_{1}-p_{2}}{2\eta\ell }x^{2}$$

Am Rand ist $v=0$

$$\Rightarrow\hspace{1cm} v_{0}=\frac{p_{1}-p_{2}}{2\eta\ell }d^{2}$$

Rohrströmung

Rohrströmung

Rand: $v=0$

$$F_{R}=2\pi r\ell\eta\frac{dv}{dr}$$

$$F_{p}=\pi r^{2}\left( p_{1}-p_{2}\right)$$

also

$$\frac{dv}{dr}=\frac{p_{1}-p_{2}}{2\eta \ell}$$

$$\Rightarrow\hspace{1cm} v=v_{0}-\frac{p_{1}-p_{2}}{2\eta l}r^{2}$$

und

$$v_{0}=\frac{p_{1}-p_{2}}{4\eta l}R^{2}$$

infinitesimal

$$v_{z}\left( r\right) =-\frac{1}{4\eta}\frac{dp}{dz}\left( R^{2}-r^{2}\right)$$

Volumenstrom: $d\dot{V}=2\pi rdr\cdot v\left( r\right)$

also:

$$\dot{V}=\int\limits_{0}^{R}2\pi rv\left( r\right) dr=\frac{\pi\left( p_{1}-p_{2}\right) }{8\eta l}R^{4}=\frac{\pi}{8\eta}\frac{dp}{dz}R^{4}$$ (3.257)

Das ist das Gesetz von Hagen-Poiseuille, wobei $\frac{dp}{dz}=-8\eta\frac {\bar{v}}{R^{2}}$ und $\bar{v}=\frac{\dot{V}}{\pi R^{2}}$

Der Strömungswiderstand ist: $\frac{8\eta\ell}{\pi R^{4}}$

Druck und Volumenstrom

Druckkraft $F_{p}=\pi R^{2}\left( p_{1}-p_{2}\right) =\frac{8\eta l}{R^{2} }\dot{V}$

Strömung um Kugeln

Strömung um eine Kugel

Die Kugel hat im Abstand $r$ keinen Einfluss mehr auf die Strömung

$$\Rightarrow\hspace{1cm}-\frac{dv}{dz}\sim\frac{v}{r}$$

Oberfläche: $4\pi r^{2}$

$$F\approx\eta\frac{dv}{dz}\cdot A=-\eta\frac{v}{r}4\pi r^{2}\sim-4\pi\eta vr$$

Genauer erhält man

$$F=-6\pi\eta vr$$ (3.258)

das Stokes-Gesetz.

Das Navier-Stokes-Gesetz, Bewegungsgleichung einer Flüssigkeit

Das Navier-Stokes-Gesetz, die Bewegungsgleichung für viskose Medien lautet

$$\rho\frac{d\vec{v}}{dt} = \rho\left(\frac{\partial \vec{v}}{\partial t}+\textrm{grad} {}\frac{\vec{v}^2}{2}-\vec{v}\times\textrm{rot} {}\vec{v}\right)$$ (3.259)

$$= \vec{F}_V -\textrm{grad} {}p +\eta \Delta \vec{v}+\left(\zeta+\frac{\eta}{3}\right) \textrm{grad} {}\textrm{div} {} \vec{v}$$

$$= \vec{F}_V-\textrm{grad} {}p -\eta \textrm{rot} {}\left(\textr... ...ght)+\left(\zeta+\frac{4\eta}{3}\right)\textrm{grad} {}\textrm{div} {}\vec{v}$$

wobei $\eta$ die Scherviskosität, $\zeta$ die Volumenviskosität, $\vec{F}_V$ die Volumenkraft (siehe Gleichung (3.57) ), $\vec{v}$ das Geschwindigkeitsfeld und $\rho$ die Dichte der Flüssigkeit ist.

Die Volumenviskosität kommt von der Dissipation beim Komprimieren oder Expandieren von Flüssigkeiten oder Gasen. Wenn wir eine Volumenelement $dxdydz$ betrachten, so messen wir neben der statischen Druckkraft

$$p\; \vec{n}= \frac{F}{A}$$

noch eine dynamische Spannung

$$\sigma^*\;\vec{n}= \frac{F^*}{A}=-\zeta\frac{1}{\rho}\frac{d\rho}{dt}\;\vec{n}= -\zeta\left(\frac{d}{dt}\ln\rho\right)\vec{n}$$

Die Einheit der Volumenviskosität oder zweiten Viskosität ist $[\zeta]=\frac{kg}{m s}= [\eta]$.

In Spezialfällen vereinfacht sich die Navier-Stokes-Gleichung.

Reibungslose Medien

Hier ist $\eta = 0$ und $\zeta = 0$. Die Navier-Stokes-Gleichung vereinfacht sich zu

$$\rho\frac{d\vec{v}}{dt} = \rho\left(\frac{\partial \vec{v}}{\partial t}+\textrm{grad} {}\frac{\vec{v}^2}{2}-\vec{v}\times\textrm{rot} {}\vec{v}\right)$$ (3.260)

$$= \vec{F}_V -\textrm{grad} {}p$$

Inkompressible Medien

Hier ist $\textrm{div} {} \vec{v}=0$. Die Navier-Stokes-Gleichung wird also

$$\rho\frac{d\vec{v}}{dt} = \rho\left(\frac{\partial \vec{v}}{\partial t}+\textrm{grad} {}\frac{\vec{v}^2}{2}-\vec{v}\times\textrm{rot} {}\vec{v}\right)$$ (3.261)

$$= \vec{F}_V -\textrm{grad} {}p +\eta \Delta \vec{v}$$

$$= \vec{F}_V-\textrm{grad} {}p -\eta \textrm{rot} {}\left(\textrm{rot} {} \vec{v}\right)$$

Da das Volumen sich nicht ändert, hat die Volumenviskosität keinen Einfluss.

Potentialströmungen inkompressibler Flüssigkeiten

Hier gilt $\textrm{div} {} \vec{v}=0$ und $\textrm{rot} {} \vec{v}=0$. Ebenso verschwindet die partielle Ableitung nach der Zeit. Die Navier-Stokes-Gleichung wird zu

$$\rho\frac{d\vec{v}}{dt} = \rho \textrm{grad} {}\frac{\vec{v}^2}{2}$$ (3.262)

$$= \vec{F}_V -\textrm{grad} {}p$$

Dies ist auch die Bernoulli-Gleichung.

Potentialströmung

Hier gilt $\textrm{rot} {} \vec{v}=0$. Ebenso verschwindet die partielle Ableitung nach der Zeit. Die Navier-Stokes-Gleichung wird zu

$$\rho\frac{d\vec{v}}{dt} = \rho\left(\textrm{grad} {}\frac{\vec{v}^2}{2}\right)$$ (3.263)

$$= \vec{F}_V-\textrm{grad} {}p +\left(\zeta+\frac{4\eta}{3}\right)\textrm{grad} {}\textrm{div} {}\vec{v}$$

Prandtl-Grenzschicht

Prandtl-Grenzschicht

Reibungskraft: $F_{R}=\eta A\frac{V}{D}$

Verschieben um $\ell$:

$$W=F_{R}\cdot\ell=\eta A\frac{V}{D}\ell$$ (3.264)

Kinetische Energie in der Grenzschicht:

$$E_{kin}=\frac{1}{2}\int\limits_{0}^{D}A\rho dz\left( v\cdot\frac{z}{D} ^{2}\right) =\frac{1}{6}A\rho v^{2}D$$ (3.265)

mit $W=E_{kin}$wird$$

$$D=\sqrt{\frac{6\eta l}{\rho v}}$$

Reynolds-Kriterium

$$D \ll\ell$$

$$\Rightarrow\frac{\ell}{D}\gg1$$

$$\sqrt{\frac{\rho v\ell^{2}}{6\eta l}} =\sqrt{\frac{\rho v\ell}{6\eta}}\gg1$$ (3.266)

$$\mathop{\mathcode\lq \'''27\mathcode\lq \*''2A\mathcode\lq \.''613A \mathc... .../''2F\mathcode\lq \:''603A \kern0pt\fam0 Re}\nolimits =\frac{\rho v\ell}{\eta}\gg1$$ (3.267)

dabei ist $\ell$ eine typische Dimension und $v$ die mittlere Geschwindigkeit.

wenn

• $Re\gg1$: turbulente Strömung mit laminarer Grenzschicht

• $Re\ll1$: laminare Strömung (Grenzschicht macht keinen Sinn)

Allgemein gilt: Es gibt für jede Geometrie eine kinetische Reynoldszahl. $Re_{krit}$ mit

$$\mathop{\mathcode\lq \'''27\mathcode\lq \*''2A\mathcode\lq \.''613A\mathcode\lq \-''2D\mathcode\lq \/''2F\mathcode\lq \:''603A \kern0pt\fam0 Re}\nolimits >\mathop{\mathcode\lq \'''27\mathcode\lq \*''2A\mathcode\lq \.''613A\mathc... ...mathcode\lq \/''2F\mathcode\lq \:''603A \kern0pt\fam0 Re}\nolimits _{krit}\Rightarrow$$

$$\mathop{\mathcode\lq \'''27\mathcode\lq \*''2A\mathcode\lq \.''613A\mathcode\lq \-''2D\mathcode\lq \/''2F\mathcode\lq \:''603A \kern0pt\fam0 Re}\nolimits <\mathop{\mathcode\lq \'''27\mathcode\lq \*''2A\mathcode\lq \.''613A\mathc... ...mathcode\lq \/''2F\mathcode\lq \:''603A \kern0pt\fam0 Re}\nolimits _{krit}\Rightarrow$$ (3.268)

Bem. Strömungen mit der gleichen Reynoldszahl sind ähnlich $\Rightarrow$ Windkanal

Bei $D\ll l$ ist $F_{R}=\eta A\frac{V}{D}=A\sqrt{\frac{V^{3}\eta\rho}{l}}$

$F_{R}$ ist etwa der Mittelwert aus

$$\left. \begin{array}[c]{c} \text{Stokes}\sim\eta v\ell\text{... ...ent, zu gross} \end{array}\right\} \text{f{\uml u}r Flugzeuge}$$

Strömung idealer Flüssigkeiten

Ideale Strömung

ideal: keine Reibung, laminar

Volumenerhaltung: $A_{1}\Delta x_{1}=A_{2}\Delta x_{2}=\Delta V$

Arbeit: $\Delta W_{1}=p_{1}A_{1}\Delta x_{1}$

$$\Delta W_{2}=p_{2}A_{2}\Delta x_{2}$$

Energiebilanz: $\Delta W=\Delta W_{1}-\Delta W_{2}=\Delta E_{kin}$

$$\left( p_{1}-p_{2}\right) \Delta V=\frac{1} {2}\rho\Delta V\left( v_{2}^{2}-v_{1}^{2}\right)$$

also

$$p+\frac{1}{2}\rho v^{2}=p_{0}=const$$ (3.269)

Dies ist die Bernoulli-Gleichung.

Ist die Flüssigkeit im Gravitationsfeld, muss noch $\rho gh$ berücksichtigt werden (allg. $W_{pot})$

Anwendungen

Manometer

$_{}$

Manometer

Ein Manometer misst nur den statistischen Druck

Prandtlsches Staurohr

$_{}$

Prandtlsches Staurohr

$$p_{0}-p=\frac{1}{2}\rho v^{2}$$ (3.270)

Strömung aus einem Loch

$_{}$

Ausströmen aus einem Loch. Der äussere Druck $p_0$, der ja überall gleich ist, wurde vernachlässigt.

$$\frac{1}{2}\rho v^{2} =p$$

$$v =\sqrt{\frac{2p}{\rho}}$$ (3.271)

Wenn der Druck $p$ ein Schweredruck ist, gilt

$$p=\rho gh$$

Es folgt

$$v=\sqrt{2gh}$$

Wenn $v>\sqrt{\frac{2p_{0}}{\rho}}=v_{k}$ wird der statische Druck $<0$. Negative Drucke können nicht existieren. Das System reagiert mit einer Dampfbildung (Kavitation).

Strömungswiderstand

Stromlinie

Nach Bernoulli ist der Druck vorne und hinten gleich. Also gäbe es keinen Widerstand (Paradoxon von d'Alembert).

Reales Bild einer Wirbelstrasse

Def. Wirbel: Wenn ein Boot auf einem geschlossenen Weg angetrieben wird

Def. Zirkulation :

$$\Gamma=\oint\vec{v}d\vec{s}=\int \textrm{rot} {} \vec{v }d\vec{a}\neq0$$ (3.272)

Potentialwirbel

$$v_{r} =0$$ (3.273)

$$v_{\varphi} =\frac{\Gamma}{2\pi r}$$ (3.274)

Beim Potentialwirbel gilt:

$$\begin{array}[c]{c} \textrm{rot} {} \vec{v=0 }\text{f{\um... ...textrm{rot} {} \vec{v\neq0}\text{ f{\uml u}r }r=0 \end{array}$$

Für $r\neq 0$ existiert also ein Geschwindigkeitspotential $\phi=\frac{\Gamma}{2\pi}\varphi$

Druck und Druckgradient

Nach Bernoulli gibt es einen Radius $r_0$, bei dem der dynamische Druck negativ würde. $r_0$ ergibt sich zu

$$p =p_{0}-\frac{1}{2}\rho v^{2}=p_{0}-\rho\frac{\Gamma^2}{8\pi^{2} }\cdot\frac{1}{r^{2}}$$

$$p =0 r_{0}=\frac{\Gamma}{2\pi}\left( \frac{\rho}{2p_{0} }\right) ^{\frac{1}{2}}$$ (3.275)

d.h. für $r<r_{0}$ ist das Konzept des Potentialwirbels nicht sinnvoll.

Volumenkraft

$$\vec{F}_{V}=-\textrm{grad} {}p=-\rho\frac{\Gamma^{2}}{4\pi^{2}}\frac{\vec{r}}{r^{4}}$$ (3.276)

(nach innen gerichtet)

Helmholtzsche Wirbelsätze

Helmholtzsche Wirbelsätze

In Wirbeln gelten die folgenden Gesetze:

1. Wirbelsatz

Im Inneren eines Gases oder einer Flüssigkeit können keine Wirbel beginnen.

2. Wirbelsatz

Wirbel enthalten zu jeder Zeit die gleichen Teilchen

3. Wirbelsatz

Die Zirkulation

$$\Gamma = \oint \vec{v}\cdot d\vec{s}$$ (3.277)

ist für jeden Wirbelquerschnitt $Q$ senkrecht zum Wirbelfaden konstant.