Gesetz von Hagen-Poiseuille

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

Formulierung

Das Gesetz lautet

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/“:): \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:

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/“:): \dot \mathbf{V} = -\frac{\pi \cdot r^4}{8 \cdot \eta}\frac{\partial p}{\partial \mathbf{x}}

mit

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.^[2]

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 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/“:): v(s) die Strömungsgeschwindigkeit an der Stelle 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/“:): s eines kreisförmigen Rohres mit Radius $ R $. Betrachten wir einen Hohlzylinder der Länge 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/“:): l und der Wanddicke 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/“:): ds zwischen den Radien 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/“:): s 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/“:): 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 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/“:): A^+ bzw. 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/“:): A mit dem dynamische Viskosität 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/“:): \eta sowie der Druckdifferenz 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/“:): \Delta p auf die Hohlzylinder-Grundfläche 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/“:): dA = 2\pi \,s\,ds ergibt sich die Kraftgleichung:

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/“:): \eta (A^+ dv^+ -A\,dv)/ds + dA \Delta p = 0 .

Dabei 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/“:): \eta A^+ dv^+/ds die Reibung mit dem nach außen benachbarten Strömungszylinder, der den Radius 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/“:): s+ds hat. Die Geschwindigkeitsdifferenz 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/“:): \frac{dv^+}{ds} = \frac{dv}{ds}\vert_{s+ds} verteilt sich auf die Schichtdicke $ ds $ und wirkt entlang der Außenfläche 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/“:): 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 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/“:): ds \to 0 ergibt sich eine Differentialgleichung zweiter Ordnung fü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/“:): v(s) :

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/“:): \begin{align} 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^2v}{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^2v}{ds^2}\,ds+\frac{ds\,dv}{ds} + \frac{s\,ds\,\Delta p}{\eta\, l} &= 0 \\ \frac{d^2v}{ds^2}+\frac{dv}{s\, ds} + \frac{\Delta p}{\eta\, l} &= 0 \\ \end{align}

Die Lösung muss die Randbedingung 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/“:): v(R) = 0 (Haftbedingung) erfüllen und ist dadurch eindeutig bestimmt:

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/“:): 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:

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/“:): \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/“:): \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/“:): \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/“:): \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/“:): \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/“:): \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/“:): = \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/“:): = \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/“:): 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/“:): 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/“:): \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/“:): 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:

Formeln für weitere Querschnittsformen werden in vielen Lehrbüchern^[3] 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/“:): \dot V = \frac{\pi}{4\eta} \cdot \frac{a^3 \cdot b^3}{a^2 + b^2} \cdot \frac{\Delta p}{l}\,,

wobei $ 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/“:): 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/“:): 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/“:): \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

Einzelnachweise