Quelltext der Seite Idealreaktor

Chemische Reaktionen finden in der [[Chemische Reaktionstechnik|Chemischen Reaktionstechnik]] in [[Chemischer Reaktor|Reaktoren]] statt. Da die [[Hydrodynamik]] solcher Reaktoren äußerst kompliziert ist, werden diese Reaktoren [[Ideales Objekt|idealisiert]] betrachtet. Dies nennt man '''Idealreaktor''' (abgekürzt&nbsp;IR).

Das bedeutet, dass [[Mathematik|mathematisch]] schwierig erfassbare reale Bedingungen, die in der Praxis jedoch immer in einem gewissen Rahmen vorhanden sind, z.&nbsp;B. unvollständige Rückvermischung oder [[turbulente Strömung]] mit [[Laminare Strömung|laminaren]] Bereichen an der Reaktorwand, durch einfache Bedingungen ersetzt werden, z.&nbsp;B. durch vollständige Rückvermischung oder ideale [[Pfropfenströmung]]. Dies ermöglicht es, Konzentrations- oder Temperaturverläufe als [[Analytische Funktion|analytische Lösung]] der [[Bilanz]]-[[Differentialgleichung]]en zu gewinnen.

== Die Grundlagen der Berechnung von Konzentrations- und Temperaturverläufen ==
Als Basis der Reaktorberechnung dienen die Erhaltungsgleichungen für Masse und Energie, welche zuerst in allgemeiner Schreibweise für eine Komponente dargestellt sind:<ref name=":0">{{Literatur |Autor=G. Herbert Vogel |Titel=Verfahrensentwicklung: Von der ersten Idee zur chemischen Produktionsanlage |Auflage=1 |Verlag=Wiley |Datum=2002 |ISBN=978-3-527-28721-5 |DOI=10.1002/9783527611232}}</ref>

Masse: <math>\frac{\partial c}{\partial t} = -\nabla\cdot (c\cdot\vec u) + \nabla\cdot (D\cdot\nabla c) + \nu\cdot r</math>

Dabei steht <math>c</math> für die molare Konzentration, <math>\vec u</math> für die Strömungsgeschwindigkeit, <math>D</math> für den Diffusionskoeffizienten und <math>\nu \cdot r</math> für eine Reaktion <chem>\nu A ->[k_A] P</chem>, wobei <math>\nu</math> der stöchiometrische Koeffizient der Komponente A und <math>r </math> die Reaktionsgeschwindigkeit ist. Im einfachsten Fall einer Reaktionskinetik erster Ordnung entspricht der Term <math>\nu \cdot r</math> = <math>-k_A\cdot c_A</math>. Wird obige Gleichung auf eine Dimension reduziert ergibt sich folgender Ausdruck:

:<math>\frac{\partial c}{\partial t} = -u\cdot \frac{\partial c}{\partial x} + D\cdot\frac{\partial^2 c}{\partial x^2} + \nu\cdot r </math>

Dabei wird angenommen, dass sowohl die Strömungsgeschwindigkeit, als auch der Diffusionskoeffizient nicht vom Ort abhängen. Mithilfe dieser Gleichung können die Konzentrationsverläufe von dem kontinuierlichen/diskontinuierlichen Rührkessel und dem idealen Strömungsrohr hergeleitet werden.

Energie: <math>\frac{\partial T}{\partial t} = -\nabla\cdot (T\cdot\vec u) + \nabla\cdot \left(\frac{\lambda}{c_\mathrm{p}\cdot\rho }\cdot\nabla T\right) + \frac{(-\Delta_RH)\cdot r}{c_\mathrm{p}\cdot\rho } + \frac{k\cdot A}{c_\mathrm{p}\cdot\rho\cdot V_R} \cdot (T_K - T)</math>

Dabei steht <math>T</math> für die Temperatur des Reaktionsmediums, <math>\lambda</math> für die Wärmeleitfähigkeit des Reaktionsmediums, <math>c_\mathrm{p}</math> für die Wärmekapazität des Reaktionsmediums, <math>\rho</math> für die Dichte des Reaktionsmediums, <math>\Delta_RH</math> für die Reaktionsenthalpie der ablaufenden Reaktion, <math>k</math> für den Wärmedurchgangskoeffizienten zwischen Reaktionsmedium und Kühlmittel, <math>A</math> für die Wärmeaustauschfläche, <math>V_R</math> für das Volumen des Reaktionsmediums und <math>T_K</math> für die Temperatur des Kühlwassers. Auch diese Gleichung kann auf eine Dimension reduziert werden, wodurch sich folgender Ausdruck ergibt.

:<math>\frac{\partial T}{\partial t} = -u\cdot \frac{\partial T}{\partial x} + \frac{\lambda}{c_\mathrm{p}\cdot\rho}\cdot\frac{\partial^2 T}{\partial x^2} + \frac{(-\Delta_RH)\cdot r}{c_\mathrm{p}\cdot\rho } + \frac{k\cdot A}{c_\mathrm{p}\cdot\rho\cdot V_R} \cdot (T_K - T) </math>

Dabei wird angenommen, dass sowohl die Strömungsgeschwindigkeit, als auch die Wärmeleitfähigkeit nicht vom Ort abhängen.

Im Allgemeinen wird zwischen folgenden idealen Reaktoren unterschieden:

== Diskontinuierlicher Betrieb ==
Beim [[Chargenprozess|diskontinuierlichen Betrieb]] wird der leere Reaktor mit den [[Reaktant|Ausgangsstoffen]] und den gegebenenfalls notwendigen [[Lösemittel]]n befüllt, und die chemische Reaktion läuft an. Nach Abschluss der Reaktion wird der Inhalt mit den [[Produkt (Chemie)|Produkten]] und den übriggebliebenen [[Edukt]]en entnommen. Der Reaktor wird gesäubert, bevor er wieder eingesetzt wird. Wird ein Reaktor so verwendet, so bezeichnet man ihn auch als ''Satzreaktor'' (engl.: ''batch reactor'').

=== Der diskontinuierliche, ideale Rührkessel ===
In einem ''idealen'' diskontinuierlichen [[Rührkessel]] (engl.: ''Stirred tank reactor''&nbsp;(STR) oder auch ''Batch-Reaktor'') wird die [[Reaktionsmasse (Chemie)|Reaktionsmasse]] vorgelegt. Die Zusammensetzung im Reaktor bzw. die [[Konzentration (Chemie)|Konzentrationen]] verändern sich dann über die Zeit:

:<math>\frac{\partial c}{\partial t} \neq 0</math>.

Es erfolgt eine ''ideale'', also beliebig schnelle Durchmischung des Reaktorinhalts. Dadurch sind im Reaktionsraum zu jedem diskreten Zeitpunkt die Konzentrationen und die Temperatur an jedem Ort gleich:

:<math>\frac{\partial c}{\partial x} = 0</math> und <math>\frac{\partial T}{\partial x} = 0</math>,

der Zustand ist also räumlich, nicht aber zeitlich, „[[Gradient (Mathematik)|gradienten]]<nowiki />frei“.

== Kontinuierlicher Betrieb ==
Beim [[Kontinuierlicher Prozess|kontinuierlichen Betrieb]] wird der Reaktor an einer Stelle kontinuierlich mit neuen Edukten befüllt und diese an anderer Stelle zusammen mit den gebildeten Produkten kontinuierlich entnommen.

=== Kontinuierlicher, idealer Rührkessel (KIK) ===
Beim kontinuierlichen, idealen Rührkessel (engl.: ''Continuously stirred Tank reactor''&nbsp;(CSTR)) findet eine vollständige Rückvermischung statt. Unter „ideal“ sei hier gleichzeitig auch verstanden, dass die zufließenden Edukte mit unendlich hoher Geschwindigkeit mit der Masse im Reaktionsraum „Rührkessel“ vermischt werden.

Die [[Massenkonzentration]]en aller Stoffe (also auch der Produkte) sind sowohl im Kessel selbst als auch im Auslauf örtlich ''und zeitlich'' konstant (örtlich und zeitlich [[Stationärer Vorgang|stationär]] = „raumzeitstationär“):

:<math>\frac{\partial c}{\partial x} = 0</math> und <math>\frac{\partial c}{\partial t} = 0</math>

Dies ist für reaktionstechnische und [[Kinetik (Chemie)|kinetische]] Untersuchungen positiv, weil stationär anfallende Messwerte besser oder genauer erfasst werden können.

Weil die jeweiligen Edukte im Reaktor miteinander reagieren und damit Edukt verbraucht wird, wird die Menge pro Volumeneinheit (d.&nbsp;h. die Konzentration) der Edukte im Reaktionsraum niedriger sein als im Zulaufstrom. Daher liegt an der Eintrittsstelle ein „Konzentrationssprung“ vor, d.&nbsp;h. sowie die Edukte in den Kessel kommen, ist ihre Konzentration sofort auf die im Kessel herrschende niedrigere Konzentration abgesunken:

:<math>c_\text{Kessel} < c_\text{Zulauf}</math>

==== Anwendung des Massenerhalts auf eine Reaktion erster Ordnung ====
Wird die, auf eine Dimension reduzierte, Gleichung des Massenerhalts auf einen CSTR, für eine Reaktion erster Ordnung <chem>A ->[k_A] P</chem>, angewendet, können die oben beschriebenen Vereinfachungen vorgenommen werden. Zusätzlich kann der Diffusionsterm, aufgrund der vollständigen Rückvermischung, vernachlässigt werden. Aufgrund der Annahme der zeitlichen Stationarität hängt die Differentialgleichung nur noch von einer Variablen ab, wodurch die partiellen Ableitungen durch Ableitungen für Funktionen einer Variablen ersetzt werden.<ref name=":0" />

:<math>D\cdot \frac{\mathrm{d}^2 c}{\mathrm{d} x^2} = 0</math>

Wird wie oben das Innere des Reaktors betrachtet, so ist die örtliche Änderung der Konzentration gleich null. Wird der Reaktor jedoch als Bilanzraum betrachtet, in den lediglich ein Eduktstrom ein und ein Produktstrom ausfließt, kann der Teil der örtlichen Ableitung als Differenzquotient geschrieben werden.

:<math>-u \cdot \frac{\mathrm{d}c}{\mathrm{d}x} = -u \cdot \frac{\Delta c}{\Delta x} = u \cdot \frac{c_0 - c}{L} = \frac{1}{\tau}\cdot(c_0 - c)</math>

Wobei <math>c_0</math> die Eingangskonzentration des Edukts A und <math>\tau</math> die hydrodynamische Verweilzeit des Reaktors ist, welche sich wie folgt berechnet:

:<math>\tau = \frac{L}{u} = \frac{V_R}{\dot{V}}</math>

Die hydrodynamische Verweilzeit gibt an, wie lange das Fluid eines gegebenen Volumenstroms <math>\dot V</math> braucht, um durch das Volumen des Reaktors <math>V_R</math> zu strömen. Mit den gegebenen Relationen lässt sich die Gleichung des Massenerhalts wie folgt schreiben:

:<math>0 = \frac{1}{\tau}\cdot(c_0 - c) - k_A\cdot c </math>

Wird die nun algebraische Gleichung nach der Konzentration <math>c</math> des Edukts A aufgelöst, kann folgender Ausdruck für die Konzentration im Reaktor bzw. am Ausgang des Reaktors erhalten werden:

:<math>c = \frac{c_0}{1 + k_A\cdot \tau} = \frac{c_0}{1 + \mathrm{Da_I}}</math>

Hierbei wurde die [[Damköhler-Zahl]] mit der Definition <math>\mathrm{Da_I} = k_A \cdot \tau \cdot c_0^{n-1}</math> und der Reaktionsordnung <math>n = 1</math> verwendet. Dieses Vorgehen kann für Reaktionen beliebiger Reaktionsordnung angewendet werden, indem lediglich der Reaktionsterm <math>\nu \cdot r</math> auf die jeweilige Reaktion angepasst wird. Reagiert die betrachtete Komponente in mehr als einer Reaktion, so setzt sich der Gesamt-Reaktionsterm additiv aus den einzelnen Reaktionen zusammen.

Bedingt durch den [[Volumenstrom]] durch den Kessel, hat die Reaktionsmasse eine begrenzte [[Verweilzeit (technischer Prozess)|Verweildauer]] im Reaktionsraum – und dies „zu allem Überfluss auch noch“ mit einer Verlassens-Wahrscheinlichkeit geprägt, also ''nicht'' diskret, sondern in einer „Bandbreite“ = Verweilzeit[[Verteilung einer Zufallsvariablen|verteilung]]. Für Reaktionen [[Reaktionsordnung|positiver Ordnung]] ist die spezifische Produktleistung daher immer kleiner als in einem Idealreaktor.

==== Rührkesselkaskade ====
Durch die Hintereinanderschaltung mehrerer&nbsp;KIK zu einer Kaskade nähert sich ein derartiges System mit steigender Anzahl der Rührkessel dem Verhalten des Idealreaktors an: ab zehn Rührkesseln ist im Realfall kein Unterschied mehr messbar.

{{Anker|Rohrreaktor}}
=== Ideales Strömungsrohr (IR) ===
In einem idealisierten, [[Rohr (Technik)|rohr]]<nowiki />förmigen Reaktor [engl.: ''Plug-flow reactor''&nbsp;(PFR) bzw. ''Plug-flow tubular reactor''&nbsp;(PFTR)] liegt eine [[Pfropfenströmung]] vor. Man kann sich diese Strömung als eine Wanderung einer sehr langen Reihe sehr dünner Scheibchen der Reaktionsmasse durch das Rohr vorstellen, die untereinander ''keinen'' Stoff- oder Wärmeaustausch haben.

Hier finden die Stoffumwandlungen entlang des Fließweges&nbsp;x statt, die Konzentrationen der Stoffe ändern sich entlang des Rohres:

:<math>\frac{\partial c}{\partial x} \neq 0</math>

An einem bestimmten Ort&nbsp;x im Rohr hingegen liegen zu jedem Zeitpunkt die gleichen Konzentrationen vor (da die Reaktion innerhalb aller Scheibchen gleich schnell abläuft):

:<math>\left.\frac{\partial c}{\partial t}\right|_x = 0</math>

Betrachtet man den Reaktionsverlauf ''innerhalb eines bestimmten Scheibchens während seiner Wanderung durch das Rohr'', so lässt er sich genauso betrachten wie in einem diskontinuierlichen, idealen Rührkessel:

:<math>\frac{\partial c}{\partial t} \neq 0</math>

Innerhalb eines einzelnen Scheibchens (Radius&nbsp;R, Laufvariable&nbsp;r) sind zu einem bestimmten Zeitpunkt alle Konzentrationen und die Temperatur räumlich konstant:

:<math>\frac{\partial c}{\partial r} = 0</math> und <math>\frac{\partial T}{\partial r} = 0</math>

==== Anwendung des Massenerhalts auf eine Reaktion erster Ordnung ====
In diesem Beispiel wird eine Reaktion erster Ordnung <chem>A ->[k_A] P</chem> in einem PFR betrachtet und ein Ausdruck hergeleitet, der den Konzentrationsverlauf entlang des Rohrreaktors beschreibt. Werden die obigen Annahmen der zeitlichen Stationarität und der ausbleibenden axialen Vermischung auf die Gleichung des Massenerhalts angewendet, kann folgender Ausdruck erhalten werden:<ref name=":0" />

:<math>0 = -u \cdot \frac{\mathrm d c}{\mathrm d x} - k_A \cdot c</math>

Wobei der Diffusionsterm, durch die Forderung ausbleibender axialer Vermischung, gleich null gesetzt wird. Diese homogene Differentialgleichung erster Ordnung kann über die Trennung der Variablen integriert werden. Dabei kann als Randbedingungen der Ein- und Ausgang des Rohrreaktors verwendet werden um die Konzentration <math>c_e</math> am Ende des Reaktors zu erhalten.

:<math>u\cdot \int_{c_0}^{c_e} \frac{1}{c} \, \mathrm dc = -k_A \int_0^L \, \mathrm dx</math>

:<math>\ln\left(\frac{c_e}{c_0}\right) = -k_A \cdot \frac{L}{u} = -k_A \cdot\tau</math>

:<math>c_e = c_0 \cdot e^{-k_A\cdot\tau}</math>

Wird als Randbedingung der Reaktoreingang und ein beliebiger Ort <math>x</math> im Reaktor verwendet, ergibt sich dadurch der Konzentrationsverlauf entlang des Rohrreaktors wie folgt:

:<math>c(x) = c_0 \cdot e^{-k_A \cdot \frac{x}{u}}</math>

== Einzelnachweise ==
<references />

[[Kategorie:Chemischer Reaktor]]