Gesetz von Hagen-Poiseuille

Mit dem Gesetz von Hagen-Poiseuille [poaː'zœj] (nach Gotthilf Heinrich Ludwig Hagen, 1797–1884 und Jean Léonard Marie Poiseuille, 1797–1869) wird der Volumenstrom ${\dot {V}}$ – d. h. das geflossene Volumen V pro Zeitspanne – bei einer laminaren stationären Strömung eines homogenen Newton’schen Fluids durch ein Rohr (Kapillare) mit dem Radius $r$ und der Länge $l$ beschrieben.

Formulierung

Das Gesetz lautet

${\dot {V}}={\frac {dV}{dt}}={\frac {\pi \cdot r^{4}}{8\cdot \eta }}{\frac {\Delta p}{l}}=-{\frac {\pi \cdot r^{4}}{8\cdot \eta }}{\frac {\partial p}{\partial x}}$

oder als Vektor:

${\dot {\mathbf {V} }}=-{\frac {\pi \cdot r^{4}}{8\cdot \eta }}{\frac {\partial p}{\partial \mathbf {x} }}$

mit

Variable	Bedeutung	SI-Einheit
${\dot {V}}$	Volumenstrom durch das Rohr	m3·s−1
${\dot {\mathbf {V} }}$	Volumenstromvektor	m3·s−1
$r$	Innenradius des Rohres	m
$l$	Länge des Rohres	m
$\eta$	dynamische Viskosität der strömenden Flüssigkeit	Pa·s
$p$	Druck	Pa
$\Delta$	Mathematischer Operator für die Differenz einer Größe über die Länge $l$	-
$x$	Lagekoordinate	m
$\mathbf {x}$	Lagevektor	m

Laminares Strömungsprofil

Dieses Gesetz folgt direkt aus dem stationären, parabolischen Strömungsprofil durch ein Rohr, das aus den Navier-Stokes-Gleichungen hergeleitet werden kann – oder direkt aus der Definition der Viskosität, siehe unten. Bemerkenswert ist die Abhängigkeit des Volumendurchflusses von der vierten Potenz des Radius des Rohres. Dadurch hängt der Strömungswiderstand sehr stark vom Radius des Rohres ab, so würde beispielsweise eine Verringerung des Rohrdurchmessers auf die Hälfte den Strömungswiderstand auf das 16fache erhöhen.

Das Gesetz gilt nur für laminare Strömungen. Bei größerem Durchfluss einer Rohrleitung, verbunden mit höheren Strömungsgeschwindigkeiten bzw. größeren Abmessungen, kommt es zu turbulenten Strömungen mit wesentlich höherem Strömungswiderstand als nach Hagen-Poiseuille zu erwarten wäre. Die konkreten Verhältnisse turbulenter Strömungen werden u. a. mit den Formeln von Blasius, Nikuradse bzw. Prandtl-Colebrook beschrieben.

Das Gesetz von Hagen-Poiseuille gilt grundsätzlich nur bei vollständig ausgebildetem hydrodynamischen Strömungsprofil (parabolisches Geschwindigkeitsprofil). Strömt bspw. Flüssigkeit aus einem Tank über ein Rohr aus, dann muss das Rohr hinreichend lang sein, damit das Gesetz von Hagen-Poiseuille gültig ist. Im Anlauf der Strömung muss das parabolische Strömungsprofil unter zusätzlichem Druckverlust ("Energieaufwand") nämlich erst ausgebildet werden. Die Druckdifferenz in der oberen Formel für den Volumenstrom bezieht sich deshalb auf die Druckdifferenz einer vollständig ausgebildeten Strömung.

In sehr dünnen Röhren, in denen die Grenzschicht maßgeblich das Strömungsprofil beeinflusst und nicht sehr klein gegenüber dem Radius ist, lässt sich dieses stark vereinfachte mathematische Modell der Strömung ebenfalls nicht anwenden.

Für kompressible Fluide (wie z. B. Gase) gilt ein modifiziertes Gesetz.

Herleitung

Hier ist die Überlegung, aus der das Hagen-Poiseuille-Gesetz und das ihr zugrundeliegende Strömungsprofil folgt: Bezeichne $v(s)$ die Strömungsgeschwindigkeit an der Stelle $s$ eines kreisförmigen Rohres mit Radius $R$. Betrachten wir einen Hohlzylinder der Länge $l$ und der Wanddicke $ds$ zwischen den Radien $s$ und $s^{+}=s+ds$. Der Zylinder solle sich im Gleichgewichtszustand befinden, also keine Beschleunigung erfahren, daher ist die Summe aller auf die Flächen wirkenden Kräfte gleich null. Aus der Reibung auf die Außen- bzw. Innenfläche $A^{+}$ bzw. $A$ mit der dynamischen Viskosität $\eta$ sowie der Druckdifferenz $\Delta p$ auf die Hohlzylinder-Grundfläche $dA=2\pi \,s\,ds$ ergibt sich die Kraftgleichung:

$\eta (A^{+}dv^{+}-A\,dv)/ds+dA\Delta p=0$.

Dabei ist $\eta A^{+}dv^{+}/ds$ die Reibung mit dem nach außen benachbarten Strömungszylinder, der den Radius $s+ds$ hat. Die Geschwindigkeitsdifferenz ${\frac {dv^{+}}{ds}}={\frac {dv}{ds}}\vert _{s+ds}$ verteilt sich auf die Schichtdicke $ds$ und wirkt entlang der Außenfläche $A^{+}=2\pi s^{+}l$. Analog gilt dies für die Reibung an der Innenfläche mit dem nach innen benachbarten Strömungszylinder.

Im Grenzübergang $ds\to 0$ ergibt sich eine Differentialgleichung zweiter Ordnung für $v(s)$:

${\begin{aligned}2\pi \,l\,\left(s+ds\right){\frac {dv}{ds}}\vert _{s+ds}-2\pi \,l\,s\,{\frac {dv}{ds}}\vert _{s}+{\frac {2\pi \,s\,ds\,\Delta p}{\eta }}&=0\\\left(s+ds\right){\frac {dv}{ds}}\vert _{s+ds}-s\,{\frac {dv}{ds}}+{\frac {s\,ds\,\Delta p}{\eta \,l}}&=0\\s\,{\frac {dv}{ds}}+s\,{\frac {d^{2}v}{ds^{2}}}\,ds+{\frac {ds\,dv}{ds}}+{\mathcal {O}}\left(ds^{2}\right)-s\,{\frac {dv}{ds}}+{\frac {s\,ds\,\Delta p}{\eta \,l}}&=0\\s\,{\frac {d^{2}v}{ds^{2}}}\,ds+{\frac {ds\,dv}{ds}}+{\frac {s\,ds\,\Delta p}{\eta \,l}}&=0\\{\frac {d^{2}v}{ds^{2}}}+{\frac {dv}{s\,ds}}+{\frac {\Delta p}{\eta \,l}}&=0\\\end{aligned}}$

Die Lösung muss die Randbedingung $v(R)=0$ (Haftbedingung) erfüllen und ist dadurch eindeutig bestimmt:

$v(s)={\frac {\Delta p}{4\eta \,l}}\,(R^{2}-s^{2})$.

Dies ist genau das genannte quadratische Strömungsprofil. Durch Integration folgt dann das Gesetz von Hagen-Poiseuille:

${\dot {V}}=\int _{0}^{2\pi }\int _{0}^{R}v(s)s\,ds\,d\varphi ={\frac {\pi R^{4}}{8\eta }}{\frac {\Delta p}{l}}$.

Alternative Herleitung über das Spatprodukt:

I) Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \overrightarrow{s} = \begin{pmatrix} R\,x\cos(\varphi) \\ R\,x\sin(\varphi) \\ 0 \end{pmatrix} }

II) Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \overrightarrow{v} = \frac{\Delta p}{4 \eta\, l} \,(R^2 - |\overrightarrow{s}|^2) \begin{pmatrix} 0 \\ 0 \\ 1 \end{pmatrix}}

III) Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \dot{V} = \int_0^{2\pi} \int_0^1 {\partial\over\partial{x}}\overrightarrow{s}\times{\partial\over\partial{\varphi}}\overrightarrow{s}\bigcirc\overrightarrow{v} \,dx\,d\varphi}

Durch Einsetzen von I in II folgt:

IV) Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \overrightarrow{v} = \frac{\Delta p}{4 \eta\, l} \,(R^2 - R^2x^2) \begin{pmatrix} 0 \\ 0 \\ 1 \end{pmatrix}}

Durch Einsetzen von I, II und IV in III folgt:

V) Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \dot{V} = \int_0^{2\pi} \int_0^1 {\partial\over\partial{x}}\begin{pmatrix} R\,x\cos(\varphi) \\ R\,x\sin(\varphi) \\ 0 \end{pmatrix}\times{\partial\over\partial{\varphi}}\begin{pmatrix} R\,x\cos(\varphi) \\ R\,x\sin(\varphi) \\ 0 \end{pmatrix}\bigcirc\left[\frac{\Delta p}{4 \eta\, l} \,(R^2 - R^2x^2) \begin{pmatrix} 0 \\ 0 \\ 1 \end{pmatrix}\right] \,dx\,d\varphi = }

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle = \int_0^{2\pi} \int_0^1 \frac{\Delta p}{4 \eta\, l}(R^2 - R^2x^2)\left[\begin{pmatrix} R\cos(\varphi) \\ R\sin(\varphi) \\ 0 \end{pmatrix}\times\begin{pmatrix} -R\,x\sin(\varphi) \\ R\,x\cos(\varphi) \\ 0 \end{pmatrix}\bigcirc \begin{pmatrix} 0 \\ 0 \\ 1 \end{pmatrix}\right] \,dx\,d\varphi = }

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle = \int_0^{2\pi} \int_0^1 \frac{\Delta p}{4 \eta\, l} R^2x(R^2 - R^2x^2) \,dx\,d\varphi = \int_0^{2\pi} \frac{\Delta p R^4}{16 \eta\, l}\,d\varphi = \frac{\pi \Delta p R^4}{8 \eta\, l} }

Nicht kreisförmige Kanalquerschnitte

Rechteck-Kanal

Für einen Rechteck-Kanal mit den Abmessungen Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle b} und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle h} lässt sich dieses Gesetz in der folgenden Form angeben:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \dot V = \frac{K\cdot\min(b,h)^3\cdot\max(b,h)}{12\eta l}\cdot \Delta p}

Hierbei ist

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle K=1-\sum_{n=1}^{\infty}\frac{1}{(2n-1)^5}\cdot\frac{192}{\pi^5}\cdot\frac{\min(b,h)}{\max(b,h)}\tanh\left((2n-1)\frac{\pi}{2}\frac{\max(b,h)}{\min(b,h)}\right)}

Die Abweichung vom exakten Wert bei Berechnung von K in erster Näherung (n=1) beträgt maximal 0,67 %, in zweiter Näherung 0,06 %, in dritter Näherung 0,01 %.

Einige Beispielwerte, berechnet in dritter Näherung:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \frac{\min(b,h)}{\max(b,h)}}	0	1/10	1/5	1/4	1/3	1/2	2/3	3/4	1
Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle K}	1	0,9370	0,8740	0,8425	0,7900	0,6861	0,5873	0,5414	0,4218

Formeln für weitere Querschnittsformen werden in vielen Lehrbüchern hergeleitet.

Elliptischer Querschnitt

Für elliptische Querschnitte ergibt sich

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \dot V = \frac{\pi}{4\eta} \cdot \frac{a^3 \cdot b^3}{a^2 + b^2} \cdot \frac{\Delta p}{l}\,,}

wobei Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle a} und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle b} die beiden Halbachsen der Ellipse repräsentieren.

Man beachte den Spezialfall Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle a=b=r} ,

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \dot V = \frac{\pi \cdot r^4}{8 \cdot \eta}\frac{\Delta p}{l}\,,}

bei dem sich die Gleichung auf die Gleichung für zylindrische Röhren reduziert.

Anwendungen

Im Gültigkeitsbereich des Gesetzes bewirkt etwa die Verengung eines runden Leitungsradius um 10 % einen Durchsatzrückgang um 1 − 0,9⁴ = 34 %. Um den ursprünglichen Durchfluss bei verkleinertem Radius wieder zu erreichen, muss somit die Druckdifferenz um über 52 % steigen.

Außerdem bildet das Gesetz von Hagen-Poiseuille die Grundlage einer Vielzahl von Modellgleichungen bei der Durchströmung von Schüttgütern.

Eingeschränkte Gültigkeit im Blut

Das Gesetz von Hagen-Poiseuille bezieht sich auf Newtonsche Flüssigkeiten. Bei Newtonschen Flüssigkeiten ist die Viskosität keine Funktion der Scherrate. Ein Beispiel für eine solche Flüssigkeit ist Wasser. Das Blutplasma ist auch eine Newtonsche Flüssigkeit, nicht aber das Blut: Es ist eine inhomogene Suspension aus verschiedenen Zellen in Plasma. Hier ist die Viskosität von der Höhe der Scherrate (also der Strömungsgeschwindigkeit) abhängig. Weiterhin spielt auch die Deformierbarkeit der Erythrozyten eine Rolle. Diese können sich beispielsweise „geldrollenartig“ in dünnen Gefäßen aggregieren. Im Übrigen handelt es sich hier eher nicht um laminare, sondern turbulente Strömungszustände.

Dieses spezielle Fachgebiet der Rheologie des Blutes wird als Hämorheologie (englisch hemorheology) bezeichnet.

Literatur

• Wolfgang Beitz; Karl-Heinrich Grote (Hrsg.): Dubbel. Taschenbuch für den Maschinenbau. 20. Auflage. Springer-Verlag, Berlin / Heidelberg / New York 2001, ISBN 3-540-67777-1
• James P. Hartnett; Milivoje Kostic: Heat Transfer to Newtonian and Non-Newtonian Fluids in Rectangular Ducts. In: Advances in Heat Transfer, Volume 19, 1989
• Rainer Klinke (Hrsg.): Physiologie. Zahlreiche Tabellen. 5. Auflage. Georg Thieme Verlag, Stuttgart / New York 2005, ISBN 3-13-796005-3