Modifizierte sphärische Besselfunktionen liefern das Konzentrationsprofil, während die Damköhler-Zahl angibt, wie stark die Reaktion mit der Diffusion konkurriert. Für ein sphärisches Katalysatorkügelchen mit dem Radius (R), einer konstanten Oberflächenkonzentration und einer irreversiblen Reaktion erster Ordnung lautet die stationäre dimensionslose Konzentration

[ \Theta(r)=\frac{C(r)}{C_s} =\frac{\sinh(\phi r/R)}{(r/R)\sinh(\phi)}, \qquad \phi=\sqrt{\mathrm{Da}}, ]

wobei

[ \mathrm{Da}=\frac{kR^2}{D_{\mathrm{eff}}}. ]

Hierbei ist (k) die Geschwindigkeitskonstante erster Ordnung und (D_{\mathrm{eff}}) die effektive intrapartikuläre Diffusivität.

Die Damköhler-Zahl ist der wichtigste Screening-Parameter: Eine kleine (\mathrm{Da}) führt zu einer nahezu gleichmäßigen Konzentration im Kügelchen, während eine große (\mathrm{Da}) starke interne Konzentrationsgradienten und eine geringere Katalysatorausnutzung erzeugt. Die modifizierte sphärische Besselfunktion stellt die mathematisch korrekte reguläre Lösung für die stationäre Diffusion gekoppelt mit einer Reaktion in einer Kugel dar.

Wie das Modell ein Katalysatorkügelchen darstellt

Grundlegende Transportbilanz

Für eine stationäre, sphärisch symmetrische Diffusion mit einer Reaktion erster Ordnung gilt

[ D_{\mathrm{eff}}\frac{1}{r^2} \frac{d}{dr}\left(r^2\frac{dC}{dr}\right)-kC=0. ]

Mit dem dimensionslosen Radius (\rho=r/R) und der Konzentration (\Theta=C/C_s) ergibt sich

[ \frac{1}{\rho^2}\frac{d}{d\rho} \left(\rho^2\frac{d\Theta}{d\rho}\right) -\mathrm{Da},\Theta=0. ]

Das Modell setzt ein sphärisches Kügelchen, radiale Symmetrie, eine konstante effektive Diffusivität, isothermen Betrieb und eine Reaktionskinetik erster Ordnung im gesamten Volumen des aktiven Partikels voraus.

Randbedingungen

Im Zentrum des Kügelchens erfordert die Symmetrie

[ \left.\frac{d\Theta}{d\rho}\right|_{\rho=0}=0. ]

An der äußeren Oberfläche führt eine vorgegebene Oberflächenkonzentration zu

[ \Theta(1)=1. ]

Die Bedingung im Zentrum ist wesentlich, da sich die scheinbaren (1/\rho)-Terme in der allgemeinen Lösung so kombinieren müssen, dass eine endliche Konzentration entsteht, wenn (\rho) gegen null geht.

Warum modifizierte sphärische Besselfunktionen auftreten

Die allgemeine Lösung kann geschrieben werden als

[ \Theta(\rho) =A\frac{\sinh(\sqrt{\mathrm{Da}},\rho)}{\rho} +B\frac{\cosh(\sqrt{\mathrm{Da}},\rho)}{\rho}. ]

Der Term mit (B) divergiert im Zentrum, daher erfordert die physikalische Endlichkeit (B=0). Durch Anwendung der Oberflächenbedingung erhält man

[ \boxed{ \Theta(\rho)= \frac{\sinh(\sqrt{\mathrm{Da}},\rho)} {\rho\sinh(\sqrt{\mathrm{Da}})} }. ]

Äquivalent dazu ist dies das Verhältnis modifizierter sphärischer Besselfunktionen nullter Ordnung:

[ \Theta(\rho)= \frac{i_0(\sqrt{\mathrm{Da}},\rho)} {i_0(\sqrt{\mathrm{Da}})}, \qquad i_0(x)=\frac{\sinh x}{x}. ]

Die elementare hyperbolische Form ist für Berechnungen in Pilotanlagen meist praktischer, während die Schreibweise (i_0) die Lösung mit der umfassenderen Familie sphärischer Transportprobleme verknüpft.

Wie die Damköhler-Zahl das Profil steuert

Kleine Damköhler-Zahl

Wenn

[ \mathrm{Da}\ll 1, ]

ist die Diffusion im Vergleich zur Reaktion schnell. Die Konzentration bleibt im gesamten Kügelchen nahezu gleichmäßig, sodass der größte Teil des Katalysatorvolumens nahe der Oberflächenkonzentration arbeitet.

In diesem Bereich kann eine Erhöhung der intrinsischen Katalysatoraktivität die Gesamtrate direkt beeinflussen, da die interne Diffusion nur einen geringen Widerstand verursacht.

Mittlere Damköhler-Zahl

Wenn (\mathrm{Da}) von der Größenordnung eins ist, haben Reaktion und Diffusion einen vergleichbaren Einfluss. Zwischen der Oberfläche des Kügelchens und seinem Zentrum bildet sich ein messbarer Konzentrationsgradient.

Dies ist häufig der aussagekräftigste Bereich für Pilotanlagenexperimente, da Änderungen der Kügelchengröße, der Temperatur, der Lösungsmitteleigenschaften oder der Katalysatorbeladung sowohl bei den Konzentrationsprofilen als auch beim Umsatz beobachtbare Verschiebungen erzeugen können.

Große Damköhler-Zahl

Wenn

[ \mathrm{Da}\gg 1, ]

verbraucht die Reaktion das Edukt schnell in der äußeren Region des Kügelchens. Das Edukt dringt weniger tief ein, sodass der innere Bereich unzureichend genutzt wird.

Die charakteristische Eindringtiefe skaliert näherungsweise mit

[ \ell \sim \sqrt{\frac{D_{\mathrm{eff}}}{k}}, ]

sodass eine Erhöhung des Kügelchenradius bei konstant gehaltenem (k) und (D_{\mathrm{eff}}) im Allgemeinen die Einschränkungen durch die interne Diffusion verstärkt.

Konzentration im Zentrum

Die Konzentration im Zentrum folgt direkt aus dem Profil:

[ \Theta(0)= \lim_{\rho\to 0}\Theta(\rho)

\frac{\sqrt{\mathrm{Da}}} {\sinh(\sqrt{\mathrm{Da}})}. ]

Für kleine (\mathrm{Da}) nähert sich dieser Wert eins. Für große (\mathrm{Da}) wird er klein und liefert damit ein einfaches Maß dafür, wie tief das Edukt in das Katalysatorkügelchen eindringt.

Berechnung der Katalysatorwirksamkeit

Definition des Wirksamkeitsfaktors

Der interne Wirksamkeitsfaktor vergleicht die tatsächliche Reaktionsrate im Kügelchen mit der Rate, die eintreten würde, wenn das gesamte Kügelchen der Oberflächenkonzentration ausgesetzt wäre:

[ \eta= \frac{\text{tatsächliche gesamte Reaktionsrate}} {\text{Rate bei }C=C_s\text{ im gesamten Kügelchen}}. ]

Für eine Reaktion erster Ordnung in einer Kugel mit konstanter Oberflächenkonzentration gilt

[ \boxed{ \eta= \frac{3}{\phi^2} \left(\phi\coth\phi-1\right) }, \qquad \phi=\sqrt{\mathrm{Da}}. ]

Dieser Faktor liegt zwischen null und eins. Ein Wert nahe eins weist auf eine wirksame Nutzung des inneren Katalysatorvolumens hin, während ein kleiner Wert auf einen starken intrapartikulären Diffusionswiderstand hindeutet.

Grenzverhalten

Für (\phi\ll 1) gilt

[ \eta\approx 1-\frac{\phi^2}{15}. ]

Das Kügelchen verhält sich somit nahezu wie ein kinetisch kontrollierter Katalysator.

Für (\phi\gg 1) gilt

[ \eta\approx \frac{3}{\phi} =\frac{3}{\sqrt{\mathrm{Da}}}. ]

Die Gesamtrate wächst dann weniger stark mit dem Katalysatorvolumen, da nur ein Bereich nahe der äußeren Oberfläche wesentlich beiträgt.

Abschätzung der Reaktionsrate

Wenn (V_p) das Volumen des Kügelchens ist, kann die gesamte Reaktionsrate geschrieben werden als

[ \dot N_A

\eta,k,C_s,V_p. ]

Für ein sphärisches Kügelchen gilt

[ V_p=\frac{4}{3}\pi R^3. ]

Dieser Ausdruck ermöglicht es Ingenieuren in Pilotanlagen, eine experimentell gemessene Oberflächenkonzentration oder ein externes Stoffübergangsmodell mit der vorhergesagten Verbrauchsrate des Kügelchens zu verknüpfen.

Erweiterung des Modells um den externen Stoffübergang

Durch einen Film kontrollierte Oberflächenkonzentration

Die Bedingung einer festen Oberflächenkonzentration ist nützlich, wenn der externe Stoffübergang schnell ist. Im Betrieb einer Pilotanlage kann sich die Bulk-Konzentration (C_b) jedoch von der Oberflächenkonzentration des Kügelchens (C_s) unterscheiden, da das Kügelchen von einem Flüssigkeitsfilm umgeben ist.

Der externe molare Fluss wird üblicherweise ausgedrückt als

[ N_A=k_f(C_b-C_s), ]

wobei (k_f) der externe Stoffübergangskoeffizient des Films ist.

Der aus der sphärischen Lösung erhaltene interne Fluss muss an der Oberfläche diesem Filmfluss entsprechen.

Biot-Zahl für den Film

Die relative Bedeutung der externen und internen Diffusion kann durch eine Stoffübergangs-Biot-Zahl dargestellt werden:

[ \mathrm{Bi}m=\frac{k_fR}{D{\mathrm{eff}}}. ]

Eine große (\mathrm{Bi}_m) weist auf einen kleinen externen Filmwiderstand im Vergleich zum internen Diffusionswiderstand hin. Eine kleine (\mathrm{Bi}_m) zeigt an, dass der Film die Oberflächenkonzentration deutlich verringern kann, bevor das Edukt in das Kügelchen eindringt.

Die interne Damköhler-Zahl und die externe Biot-Zahl beantworten daher unterschiedliche Fragen: (\mathrm{Da}) misst die Reaktion relativ zur intrapartikulären Diffusion, während (\mathrm{Bi}_m) den externen Filmtransport mit der internen Diffusion vergleicht.

Interpretation der Gesamtrate

Wenn beide Widerstände relevant sind, ist die Rate vom Bulk zum Partikel geringer als die allein auf Grundlage von (C_b) vorhergesagte Rate. Eine nützliche ingenieurtechnische Darstellung lautet

[ \frac{1}{k_{\mathrm{overall}}} \approx \frac{1}{k_f} + \frac{1}{k_{\mathrm{internal}}}, ]

wobei der interne Beitrag aus dem Wirksamkeitsfaktormodell bestimmt wird.

Dieser Rahmen hilft zu unterscheiden, ob eine Pilotanlage kleinere Kügelchen verwenden, die Durchmischung verbessern, die Temperatur erhöhen oder die Katalysatorformulierung verändern sollte.

Anwendung des Modells in einer Pilotanlage

Auswahl der dimensionslosen Eingangsgrößen

Schätzen Sie für jede Betriebsbedingung

[ \mathrm{Da}=\frac{kR^2}{D_{\mathrm{eff}}}. ]

Die Geschwindigkeitskonstante (k) kann stark von der Temperatur und der Zusammensetzung des Reaktionsgemischs abhängen. Die effektive Diffusivität (D_{\mathrm{eff}}) spiegelt Porenstruktur, Porosität, Tortuosität, Lösungsmitteleigenschaften und gegebenenfalls Adsorptionseffekte wider.

Erzeugung von Konzentrationsprofilen

Wählen Sie Werte von (\rho) zwischen null und eins und berechnen Sie

[ \Theta(\rho)= \frac{\sinh(\sqrt{\mathrm{Da}},\rho)} {\rho\sinh(\sqrt{\mathrm{Da}})}. ]

Der Wert bei (\rho=0) sollte mithilfe des Grenzausdrucks und nicht durch direkte numerische Substitution berechnet werden, da die Formel im Zentrum die hebbare Form (0/0) besitzt.

Vergleich von Kügelchenauslegungen

Die Wiederholung der Berechnung für unterschiedliche Kügelchenradien zeigt den Einfluss der Partikelgröße unmittelbar. Da (\mathrm{Da}) mit (R^2) skaliert, erhöht eine Verdopplung des Kügelchenradius die interne Damköhler-Zahl um den Faktor vier, sofern (k) und (D_{\mathrm{eff}}) unverändert bleiben.

Dies macht den Kügelchendurchmesser zu einer wichtigen Auslegungsvariable, obwohl kleinere Kügelchen den Druckverlust erhöhen oder die Feststoffrückgewinnung in Festbett- oder Wirbelschichtsystemen erschweren können.

Verknüpfung der Profile mit dem gemessenen Umsatz

Das Modell sagt die lokalen Reaktionsraten voraus durch

[ r_A(r)=kC(r). ]

Die Integration dieser Rate über das Kügelchenvolumen ergibt die gesamte Partikelverbrauchsrate, die anschließend in Bilanzen für Festbett-, Suspensions- oder immobilisierte-Enzym-Reaktoren aufgenommen werden kann.

Verständnis der Zielkonflikte

Die Oberflächenkonzentration ist möglicherweise nicht bekannt

Das normierte Profil setzt voraus, dass (C_s) vorgegeben ist. Wenn Messungen nur die Bulk-Konzentration liefern, ist ein externer Filmansatz erforderlich, um (C_s) zu bestimmen.

Die direkte Verwendung von (C_b) in der Lösung mit fester Oberflächenkonzentration kann die Reaktionsrate überschätzen, wenn der Filmwiderstand erheblich ist.

Die effektive Diffusivität ist häufig unsicher

Die Zuverlässigkeit des Modells hängt wesentlich vom Wert von (D_{\mathrm{eff}}) ab. Die Annahme, dass die molekulare Diffusivität der intrapartikulären Diffusivität entspricht, kann bei porösen, gequollenen, gewundenen oder adsorbierenden Katalysatormaterialien ungenau sein.

Daten aus der Pilotanlage sollten daher nach Möglichkeit zur Validierung oder Schätzung effektiver Transportparameter verwendet werden.

Kinetik erster Ordnung kann eine Näherung sein

Die analytische Lösung gilt für eine irreversible Reaktion erster Ordnung. Sie beschreibt nicht direkt eine starke Substratinhibierung, Sättigungskinetik, reversible Reaktionen, Wärmefreisetzung oder eine konzentrationsabhängige Diffusivität.

Diese Effekte können ein nichtlineares numerisches Randwertmodell erfordern.

Der stationäre Zustand beschreibt kein Anfahrverhalten

Die modifizierte sphärische Bessellösung beschreibt das stationäre Konzentrationsfeld. Das transiente Diffusions-Reaktions-Verhalten erfordert eine zeitabhängige Gleichung und geeignete Anfangsbedingungen.

Bei transienter Diffusion ohne den Reaktionsterm können gewöhnliche sphärische Besselfunktionen oder trigonometrische Eigenfunktionen auftreten. Diese sollten nicht mit den modifizierten sphärischen Besselfunktionen verwechselt werden, die im stationären Diffusions-Reaktions-Problem verwendet werden.

Hohe Reaktionsraten können thermische Einschränkungen verursachen

Eine hohe Damköhler-Zahl kann auf einen starken Eduktverbrauch hindeuten, stellt jedoch für sich genommen keine thermische Kontrolle fest. Bei exothermen oder endothermen Reaktionen müssen Wärmeleitung und Wärmeübertragung separat modelliert werden, gegebenenfalls mit temperaturabhängigem (k) und (D_{\mathrm{eff}}).

Die richtige Wahl für Ihr Ziel treffen

Verwenden Sie zunächst das dimensionslose Profil und den Wirksamkeitsfaktor und ergänzen Sie externe Film- oder transiente Modelle, wenn die Bedingungen der Pilotanlage dies erfordern.

Wenn Ihr Hauptaugenmerk auf der Abschätzung des internen Diffusionswiderstands liegt: Berechnen Sie (\mathrm{Da}=kR^2/D_{\mathrm{eff}}), werten Sie das sphärische Profil aus und geben Sie den Wirksamkeitsfaktor (\eta) an.
Wenn Ihr Hauptaugenmerk auf der Auswahl des Kügelchendurchmessers liegt: Vergleichen Sie die (\mathrm{Da})-Werte für verschiedene mögliche Radien, da der interne Diffusionswiderstand mit dem Quadrat des Kügelchenradius zunimmt.
Wenn Ihr Hauptaugenmerk auf der Diagnose einer Begrenzung durch den externen Stoffübergang liegt: Kombinieren Sie das interne Modell mit einem Filmkoeffizienten und der Stoffübergangs-Biot-Zahl, anstatt Bulk- und Oberflächenkonzentration als gleich anzunehmen.
Wenn Ihr Hauptaugenmerk auf der Interpretation transienter Pilotanlagendaten liegt: Verwenden Sie die stationäre modifizierte Bessellösung nur als Grenzwert für lange Zeiten und formulieren Sie für Anfahr- oder Sprungänderungsexperimente ein zeitabhängiges Diffusions-Reaktions-Modell.
Wenn Ihr Hauptaugenmerk auf der Verbesserung der Katalysatorproduktivität liegt: Erhöhen Sie die intrinsische Aktivität nur dann, wenn das Kügelchen kinetisch limitiert bleibt; bei großer (\mathrm{Da}) können eine Verringerung der Partikelgröße oder eine Verbesserung der Diffusivität einen größeren Nutzen bringen.

Ein sphärisches Modell mit modifizierten Besselfunktionen verwandelt Diffusion und Reaktion auf Kügelchenebene in messbare Auslegungskriterien für Katalysatorausnutzung, Stoffübergangsdiagnose und die Maßstabsübertragung auf Pilotanlagen.

Zusammenfassungstabelle:

Parameter oder Konzept	Ingenieurtechnische Erkenntnis	Anwendung in der Pilotanlage
Damköhler-Zahl, `Da = kR²/D_eff`	Vergleicht die Reaktionsstärke mit der intrapartikulären Diffusion	Ermittelt, ob kinetische oder diffusionsbedingte Kontrolle vorliegt
Modifiziertes sphärisches Bessel-Profil	Sagt die Konzentration im gesamten Kügelchen voraus	Schätzt lokale Konzentrations- und Reaktionsratengradienten
Kleine `Da`	Nahezu gleichmäßige Konzentration und hohe Katalysatorausnutzung	Verbesserungen der intrinsischen Aktivität können die Produktivität steigern
Große `Da`	Starker interner Stoffabbau und unzureichend genutzte Kügelchenbereiche	Kleinere Kügelchen oder eine höhere effektive Diffusivität in Betracht ziehen
Wirksamkeitsfaktor, `η`	Quantifiziert den Anteil der genutzten Katalysatorkapazität	Korrigiert ideale Berechnungen der Reaktionsrate
Stoffübergangs-Biot-Zahl, `Bi_m`	Vergleicht den externen Filmwiderstand mit der internen Diffusion	Diagnostiziert, ob Durchmischung oder Partikelauslegung die Leistung begrenzt
Kügelchenradius	`Da` steigt mit dem Quadrat des Radius	Gleicht den Diffusionswiderstand gegen Druckverlust und Feststoffrückgewinnung ab
Modellgrenzen	Setzt stationären Zustand, Kinetik erster Ordnung und konstante Eigenschaften voraus	Bei Bedarf numerische transiente, thermische oder nichtlineare Modelle verwenden

Optimieren Sie die Leistung von Katalysatorkügelchen mit LABPARK

Möchten Sie Diffusionsbegrenzungen, Katalysatorausnutzung oder die Leistung bei der Maßstabsübertragung bewerten? Kontaktieren Sie LABPARK, um eine auf Ihre Anforderungen zugeschnittene Ausbildungs- oder berufsbildende Pilotanlage für verfahrenstechnische Grundoperationen zu besprechen. LABPARK unterstützt Universitäten, Forschungseinrichtungen und Unternehmen mit Pilotanlagen für Chemieingenieurwesen, Bioprozesstechnik und Biotechnologie sowie für Umwelt- und Wasseraufbereitungsanwendungen. Unsere Systeme helfen Ihnen, die Auswirkungen der Damköhler-Zahl, Konzentrationsprofile, Wirksamkeitsfaktoren, externen Stoffübergang, die Auswahl der Kügelchengröße und den Reaktorumsatz unter praxisnahen Betriebsbedingungen zu untersuchen. Nutzen Sie zuverlässige experimentelle Daten, um D_eff und kinetische Parameter zu validieren, Stofftransportengpässe zu identifizieren, Ingenieure auszubilden und fundiertere Entscheidungen bei der Prozessauslegung zu treffen.