Bewegungsgleichungen für Fluide in Zylinderkoordinaten werden angewendet, indem der Impulstransport in radiale, azimutale und axiale Komponenten aufgelöst und anschließend mit Kontinuitäts-, Energie- und Stoffgleichungen gekoppelt wird. In Rohrreaktoren und Pilotanlagen für Rohrströmungen ermöglicht diese Formulierung die Vorhersage von Geschwindigkeitsprofilen, Druckverlust, Wandschubspannung, Verweilzeitverhalten sowie Wärme- und Stofftransport unter realistischen Betriebsbedingungen. Sie ist besonders wirkungsvoll, weil sich die Geometrie der Anlagen und die Randbedingungen natürlich mithilfe des Radius (r), des Winkels (\theta) und der axialen Position (z) ausdrücken lassen.

Impulsbilanzen in Zylinderkoordinaten überführen die Geometrie eines Rohrs oder Behälters in lösbare Transportgleichungen. Die resultierenden Geschwindigkeits- und Druckfelder bilden die Grundlage für die Auslegung von Rohrleitungen und Pumpen, die Interpretation von Pilotdaten, die Diagnose nichtidealer Strömungen und die Maßstabsübertragung von Rohrsystemen.

Warum Zylinderkoordinaten verwendet werden

Die Koordinaten entsprechen der Anlagengeometrie

Ein Zylindersystem beschreibt die Position mithilfe von (r), (\theta) und (z), mit den Geschwindigkeitskomponenten (v_r), (v_\theta) und (v_z).

Bei einem geraden Rohrreaktor oder einer Rohrleitung folgt (z) der Anlagenlänge, (r) misst den Abstand von der Mittellinie und (\theta) beschreibt die Umfangsposition. Dadurch wird die Rohrwand zu einer einfachen Randbedingung bei (r=R).

Die Differentialoperatoren erfassen die radiale Geometrie

Für eine skalare Größe (f) lautet der Gradient

[ \nabla f = \frac{\partial f}{\partial r}\mathbf{e}r+ \frac{1}{r}\frac{\partial f}{\partial \theta}\mathbf{e}\theta+ \frac{\partial f}{\partial z}\mathbf{e}_z ]

Die Divergenz des Geschwindigkeitsfeldes ist

[ \nabla\cdot\mathbf{v}

\frac{1}{r}\frac{\partial}{\partial r}(rv_r) + \frac{1}{r}\frac{\partial v_\theta}{\partial \theta} + \frac{\partial v_z}{\partial z} ]

Für ein inkompressibles Fluid ist dieser Ausdruck gleich null.

Die radialen Faktoren mit (1/r) sind wesentlich. Ihr Weglassen würde zu falschen Erhaltungsgleichungen und falschen Geschwindigkeits- oder Druckfeldern führen.

Wie die Impulsgleichungen angewendet werden

Die materielle Ableitung beschreibt die Fluidbeschleunigung

Die Bewegung eines Fluidelements wird dargestellt durch

[ \frac{D}{Dt}

\frac{\partial}{\partial t} + \mathbf{v}\cdot\nabla ]

In Zylinderkoordinaten berücksichtigen die konvektiven Terme Geschwindigkeitsänderungen, die durch die Bewegung in radialer, azimutaler und axialer Richtung verursacht werden.

Die allgemeine Impulsbilanz wird geschrieben als

[ \rho\frac{D\mathbf{v}}{Dt}

-\nabla p+\nabla\cdot\boldsymbol{\tau} +\rho\mathbf{g} ]

wobei (p) der Druck, (\boldsymbol{\tau}) der viskose Spannungstensor und (\rho\mathbf{g}) Körperkräfte wie beispielsweise die Schwerkraft darstellen.

Die axiale Gleichung sagt Druckverlust und Strömungsentwicklung voraus

Bei gewöhnlicher Rohrströmung ist die axiale Komponente normalerweise am wichtigsten. Sie bilanziert axiale Konvektion, Druckgradientenkräfte, viskose Spannungen und Schwerkraft.

Bei einer vollständig entwickelten, stationären und achsensymmetrischen Rohrströmung ist die Geschwindigkeit vorwiegend axial:

[ v_r=0,\qquad v_\theta=0,\qquad v_z=v_z(r) ]

Die axiale Impulsbilanz reduziert sich dann auf einen Zusammenhang zwischen dem axialen Druckgradienten und der radialen Änderung der Schubspannung. Für eine laminare newtonsche Strömung ergibt sich daraus das bekannte parabolische Geschwindigkeitsprofil.

Die axiale Gleichung wird verwendet zur Bestimmung von:

Druckverlust entlang des Rohrs
Volumenstrom
Wandschubspannung
Reibungsbedingtem Energieverlust
Erforderlicher Förderhöhe und Pumpenleistung
Der erforderlichen Strecke, bis eine sich entwickelnde Strömung vollständig entwickelt ist

Die radiale Gleichung bestimmt die Druckvariation über den Rohrquerschnitt

Bei nicht rotierender, vollständig entwickelter Rohrströmung ist die radiale Druckvariation häufig vernachlässigbar. Die radiale Impulsbilanz wird jedoch wichtig, wenn die Strömung Drall, Krümmung, Auftrieb, Mehrphaseneffekte oder starke radiale Beschleunigung aufweist.

Die radiale Bilanz umfasst Terme im Zusammenhang mit radialen Geschwindigkeitsgradienten und gegebenenfalls der Zentrifugalbeschleunigung. Diese Terme helfen zu bestimmen, ob der Druck zwischen der Rohrmittellinie und der Wand deutlich variiert.

Die azimutale Gleichung beschreibt Drall und Rotation

Die azimutale Komponente wird benötigt, wenn (v_\theta) ungleich null ist. Dies tritt in Drallreaktoren, statischen Mischern, Zyklonabscheidern, rotierenden Anlagen und einigen gerührten zylindrischen Behältern auf.

Die radiale Impulsbilanz enthält dann den Zentrifugalterm

[ \frac{v_\theta^2}{r} ]

Dieser Term kann einen radialen Druckgradienten erzeugen und die Verweilzeitverteilung, Durchmischung und Wandbelastung verändern. Er sollte nicht in ein einfaches Modell einer nicht rotierenden Rohrströmung aufgenommen werden, sofern keine azimutale Bewegung physikalisch vorhanden ist.

Anwendung der Gleichungen auf Pilotanlagen für Rohrströmungen

Geometrie und Betriebsbedingungen festlegen

Das Modell beginnt mit Rohrradius, Länge, Einlasskonfiguration, Rauheit, Fluideigenschaften und Betriebsvolumenstrom. Dichte und Viskosität können für eine isotherme inkompressible Flüssigkeit als konstant behandelt werden; in Wärmeübertragungs- oder Reaktionssystemen können jedoch temperaturabhängige Stoffeigenschaften erforderlich sein.

Das Strömungsregime wird mithilfe der Reynolds-Zahl bewertet:

[ \mathrm{Re}=\frac{\rho U D}{\mu} ]

wobei (U) die mittlere Geschwindigkeit und (D) der Rohrdurchmesser ist.

Physikalisch korrekte Randbedingungen vorgeben

An einer feststehenden undurchlässigen Wand erfordert die Haftbedingung

[ v_z(R)=0 ]

und, sofern zutreffend,

[ v_r(R)=0,\qquad v_\theta(R)=0 ]

Am Einlass kann das Modell eine gleichförmige Geschwindigkeit, ein gemessenes Geschwindigkeitsprofil oder ein entwickeltes Profil vorgeben. Am Auslass wird üblicherweise der Druck oder eine geeignete konvektive Bedingung festgelegt.

Diese Entscheidungen sind wichtig, da sich ein vorgegebenes gleichförmiges Einlassprofil stromabwärts durch den viskosen Impulstransport entwickelt. Wird am Einlass ein vollständig entwickeltes Profil verwendet, obwohl die Pilotanlage tatsächlich eine Einlaufstrecke besitzt, können der vorhergesagte Druckverlust und die Verweilzeit verfälscht werden.

Geschwindigkeitsprofile und Wandschubspannung berechnen

Die radiale Ableitung der axialen Geschwindigkeit bestimmt die axiale Wandschubspannung:

[ \tau_{rz}

\mu\frac{\partial v_z}{\partial r} ]

für ein newtonsches Fluid in einer einfachen achsensymmetrischen Situation.

Der Geschwindigkeitsgradient ist bei einem symmetrischen Profil an der Mittellinie null und in Wandnähe am größten. Daher bestimmt der wandnahe Bereich einen großen Teil des Reibungsdruckverlusts und kann den Wärme- und Stofftransport stark beeinflussen.

Lokale Felder mit messbaren Größen verknüpfen

Messungen in Pilotanlagen umfassen gewöhnlich Volumenstrom, Einlass- und Auslassdruck, Temperatur und manchmal lokale Konzentration oder Geschwindigkeit. Die Lösung in Zylinderkoordinaten verknüpft diese Messungen mit lokalen Größen wie (v_z(r,z)), Schubspannung und Druckverteilung.

Der Volumenstrom wird beispielsweise aus

[ Q

\int_A v_z,dA

2\pi\int_0^R v_z(r),r,dr ]

ermittelt. Der Faktor (r) ergibt sich aus dem zylindrischen Flächenelement und ist für eine korrekte Integration erforderlich.

Anwendung der Gleichungen auf Rohrreaktoren

Die Hydrodynamik bestimmt die Verweilzeit

Das axiale Geschwindigkeitsfeld bestimmt, wie lange sich das Material im Reaktor aufhält. Eine Pfropfenströmungsnäherung nimmt eine im Wesentlichen gleichförmige axiale Geschwindigkeit an, während eine laminare Strömung aufgrund der langsameren Bewegung des Fluids in Wandnähe gegenüber dem Fluid an der Mittellinie zu einer breiten Verteilung lokaler Verweilzeiten führt.

Dieser Unterschied beeinflusst Umsatz, Selektivität, Temperaturverlauf und die Interpretation von Tracerexperimenten.

Der radiale Transport steuert die nichtideale Strömung

Radiale Geschwindigkeitsgradienten erzeugen Impulsdiffusion, während die molekulare radiale Diffusion und bei turbulenter Strömung die turbulente Durchmischung Wärme und Stoffe über das Rohr transportieren.

Ein Reaktormodell koppelt daher die Impulsgleichungen mit Energie- und Stoffbilanzen. Eine allgemeine Stoffbilanz hat die Form

[ \frac{Dc_i}{Dt}

\nabla\cdot(D_i\nabla c_i)+R_i ]

wobei (c_i) die Konzentration, (D_i) der Diffusionskoeffizient und (R_i) die Reaktionsgeschwindigkeit ist.

Die Impulslösung liefert das Geschwindigkeitsfeld, das im konvektiven Term verwendet wird. Die Stoff- und Energielösungen bestimmen anschließend die Konzentrations- und Temperaturverteilungen.

Druckverlust kann die Reaktorleistung beeinflussen

Druckverlust ist nicht nur eine mechanische Auslegungsgröße. In Gasphasenreaktoren beeinflusst der Druck Dichte, Geschwindigkeit, Verweilzeit und Reaktionsgeschwindigkeiten.

In Festbettrohrreaktoren erhöht das Festbett den Strömungswiderstand. Die Impulsbilanz muss dann den mit dem porösen Medium verbundenen Druckverlust berücksichtigen, der häufig mithilfe einer experimentell validierten Festbettbeziehung dargestellt wird.

Der Wärmeübergang hängt von der wandnahen Lösung ab

In einem beheizten oder gekühlten Rohrreaktor kann die Wandrandbedingung eine Wandtemperatur oder einen Wandwärmestrom vorgeben. Der radiale Temperaturgradient bestimmt den Wandwärmestrom:

[ q_r=-k\frac{\partial T}{\partial r} ]

Das Geschwindigkeitsprofil beeinflusst die thermische Grenzschicht und damit den Gesamtwärmeübergangskoeffizienten. Änderungen des Volumenstroms, der Viskosität oder der Wandbedingung können sowohl das Temperaturfeld als auch das Reaktionsergebnis verändern.

Erweiterung des Modells auf turbulente Strömung

Turbulenz verändert die erforderliche Schließung

Die momentanen Gleichungen in Zylinderkoordinaten bleiben auch für turbulente Strömungen gültig, doch die Auflösung jeder turbulenten Schwankung ist für die Auslegung von Pilotanlagen meist unpraktikabel.

Ingenieurmodelle verwenden daher gemittelte Gleichungen mit zusätzlichen turbulenten Spannungen. Diese erfordern ein Turbulenzmodell oder eine experimentell angepasste Reibungskorrelation.

Reibungskorrelationen ergänzen die Impulsbilanz

In Labor- und Pilotstudien zu Rohrströmungen wird die Druckverlustvorhersage häufig mithilfe einer Reibungsfaktorkorrelation formuliert. Das Buckingham-Pi-Theorem ordnet die relevanten Variablen wie Druckverlust, Rohrlänge und -durchmesser, Geschwindigkeit, Dichte und Viskosität in dimensionslosen Gruppen.

Die resultierenden Korrelationen stellen den Zusammenhang zwischen Reibungsfaktor und Reynolds-Zahl sowie bei rauen Rohren der relativen Rauheit her. Sie sind für schnelle Anlagenberechnungen nützlich, während die differentiellen Impulsgleichungen eine detailliertere lokale Beschreibung liefern.

Das Wandverhalten ist in turbulenten Systemen wichtig

Turbulente Rohrströmungen enthalten einen wandnahen viskosen Bereich und einen äußeren turbulenten Bereich. Der Druckgradient an der Wand und die Wandschubspannung bestimmen die Reibungsgeschwindigkeit und beeinflussen das dimensionslose Geschwindigkeitsprofil.

Pilotanlagen können verwendet werden, um gemessene Druckgradienten und Geschwindigkeitsprofile mit etablierten Reibungsgesetzen zu vergleichen. Dies ist wertvoll bei der Auswahl von Korrelationen für ein bestimmtes Fluid, Rohrmaterial oder einen bestimmten Betriebsbereich.

Die Abwägungen verstehen

Detaillierte Modelle erfordern mehr Daten

Ein multidimensionales zylindrisches Modell kann Einlaufeffekte, sich entwickelnde Profile, Drall, Wärmeübertragung und Stofftransport abbilden. Es erfordert jedoch auch zuverlässige Stoffeigenschaften, Einlassbedingungen, Wandbedingungen, Geometrie und numerische Auflösung.

Eine eindimensionale Druckverlustkorrelation ist einfacher anzuwenden, kann jedoch keine lokalen Rezirkulationen, radialen Durchmischungsbegrenzungen oder nichtgleichförmigen Temperaturfelder aufzeigen.

Achsensymmetrie ist eine Annahme, keine Garantie

Eine achsensymmetrische Modellierung ist geeignet, wenn die Anlage und die Betriebsbedingungen ungefähr gleichförmig über den Umfang sind. Sie wird unzuverlässig, wenn der Einlass asymmetrisch ist, Armaturen starke dreidimensionale Strukturen erzeugen oder der Reaktor interne Elemente enthält, die die Symmetrie aufheben.

Statische Mischer, Bögen, Einspritzstellen und teilweise blockierte Rohre können ein vollständig dreidimensionales Modell erfordern.

Laminare Annahmen können den Pilotbetrieb falsch darstellen

Eine laminare Lösung ist für viskose Fluide und niedrige Volumenströme nützlich. Ihre Anwendung auf einen Übergangs- oder turbulenten Betrieb kann jedoch die Durchmischung deutlich unterschätzen und den Druckverlust falsch berechnen.

Umgekehrt kann eine turbulente Korrelation für hochviskose oder stark nichtnewtonsche Fluide ungeeignet sein, sofern sie nicht für diese Bedingungen validiert wurde.

Fehler bei den Randbedingungen können das Ergebnis dominieren

Falsche Einlassprofile, unrealistische Wandbedingungen oder ein zu nahe am Entwicklungsbereich platzierter Auslass können numerische Ergebnisse erzeugen, die zwar stabil erscheinen, aber die physikalische Pilotanlage nicht repräsentieren.

Die Validierung anhand gemessener Druckverluste, Volumenströme, Austrittstemperaturen und Verweilzeitdaten aus Tracerversuchen ist daher Bestandteil des Modellierungsprozesses.

Die richtige Wahl für Ihr Ziel treffen

Das geeignete Modell hängt von der Entscheidung ab, die die Pilotanlage unterstützen soll.

Wenn Ihr Hauptaugenmerk auf Druckverlust und Pumpenauslegung liegt: Verwenden Sie die axiale Impulsbilanz mit validierten Reibungsfaktor- oder Festbettkorrelationen und berücksichtigen Sie gegebenenfalls Höhendifferenzen und örtliche Verluste.
Wenn Ihr Hauptaugenmerk auf der Verweilzeitverteilung liegt: Lösen Sie das axiale Geschwindigkeitsprofil auf und koppeln Sie es mit einem Tracertransportmodell, insbesondere wenn die Geschwindigkeitsunterschiede zwischen Wand und Mittellinie bei laminarer Strömung erheblich sind.
Wenn Ihr Hauptaugenmerk auf Reaktorumsatz oder Selektivität liegt: Koppeln Sie zylindrische Impuls-, Energie- und Stoffbilanzen, damit Geschwindigkeits-, Temperatur-, Konzentrations- und Reaktionsgeschwindigkeitsfelder konsistent berechnet werden.
Wenn Ihr Hauptaugenmerk auf der Wärmeübertragungsleistung liegt: Lösen Sie die radialen Geschwindigkeits- und Temperaturgradienten in Wandnähe auf und validieren Sie den vorhergesagten Wärmeübergangskoeffizienten anhand von Pilotmessungen.
Wenn Ihr Hauptaugenmerk auf Durchmischung oder Drall liegt: Behalten Sie die azimutale Impulsgleichung und den Zentrifugalterm bei und verwenden Sie ein dreidimensionales Modell, wenn die Anlage nichtachsensymmetrische Einbauten enthält.
Wenn Ihr Hauptaugenmerk auf der Maßstabsübertragung liegt: Vergleichen Sie dimensionslose Gruppen, Randbedingungen, Druckgradienten, Verweilzeitverhalten sowie Wärme- oder Stoffübergangskoeffizienten zwischen Pilot- und Industriegeometrien.

Gleichungen in Zylinderkoordinaten überführen die Geometrie und das Transportverhalten von Rohrsystemen in einen Rahmen zur Vorhersage, Messung und Verbesserung der Leistung von Pilotanlagen.

Zusammenfassungstabelle:

Anwendung	Wichtige Gleichungen oder Ergebnisse	Nutzen für die Pilotanlage
Rohrströmung	Axiale Impulsbilanz, Geschwindigkeitsprofil, Druckverlust, Wandschubspannung	Auslegung von Pumpen und Vorhersage von Reibungsverlusten
Verweilzeit	Axiales Geschwindigkeitsfeld und Tracertransport	Analyse von Umsatz und Selektivität
Wärmeübertragung	Radialer Temperaturgradient und Wandwärmestrom	Temperaturregelung des Reaktors
Stofftransport	Konvektions-, Diffusions- und Reaktionsbilanzen	Vorhersage von Konzentration und Umsatz
Drallströmung	Azimutale Impulsbilanz und zentrifugale Druckeffekte	Analyse von Durchmischung und Wandbelastung
Maßstabsübertragung	Reynolds-Zahl, Reibungsfaktor und Übergangskoeffizienten	Validierung der Auslegung vom Pilot- zum Industriemaßstab

Verbessern Sie Ihre Pilotanlagenauslegung

LABPARK unterstützt Universitäten, Forschungsinstitute und Unternehmen dabei, die Analyse von Transportgleichungen in zuverlässige Versuchsanlagen zu überführen. Unsere Pilotanlagen für Bildungs- und berufliche Ausbildungseinheiten in der chemischen Verfahrenstechnik, Bioprozess- und Biotechnologie sowie Umwelt- und Wasseraufbereitung ermöglichen präzise Untersuchungen von Strömungsverhalten, Druckverlust, Verweilzeit, Wärmeübertragung, Stoffübertragung und Reaktorleistung.

Arbeiten Sie mit LABPARK zusammen, um Anlagen auszuwählen, die zu Ihren Fluideigenschaften, Betriebsbedingungen, Forschungszielen und Anforderungen an die Maßstabsübertragung passen. Unsere Pilotanlagen helfen Ihnen, Modelle zu validieren, Ingenieure auszubilden, Betriebsbedingungen zu vergleichen und verlässliche Daten vor der Umsetzung im Großmaßstab zu gewinnen. Kontaktieren Sie LABPARK, um Ihre Anforderungen an eine Pilotanlage für Rohrreaktoren, Rohrströmungen oder verfahrenstechnische Grundoperationen zu besprechen.