Die instationäre Diffusion bietet eine praktische Möglichkeit, sowohl Konzentrationsprofile als auch zeitabhängige Stoffübergangsraten in Pilotanlagen abzuschätzen. Für eine halbunendliche Flüssigkeitsschicht, die plötzlich einer konstanten Oberflächenkonzentration ausgesetzt wird, verwendet man die Ähnlichkeitsvariable (\eta=x/[2(Dt)^{1/2}]), um ein Konzentrationsprofil in Form einer Fehlerfunktion zu erhalten. Die daraus resultierende Eindringtiefe beträgt für einen Fortschritt von 95 % ungefähr (\delta(t)=2.8(Dt)^{1/2}), während der momentane Oberflächenfluss mit (t^{-1/2}) abnimmt. Diese Ergebnisse helfen bei der Abschätzung von Kontaktzeiten, der Entwicklung von Grenzschichten, instationären Flüssen und dem Maßstabsübertragungsverhalten bei Absorption, Membranverfahren, Extraktion und verwandten Grundoperationen.

Die zentrale Erkenntnis ist, dass die instationäre Diffusion eine Stoffübergangszone erzeugt, deren Dicke mit (\sqrt{Dt}) wächst, während der momentane Fluss mit (1/\sqrt{t}) abnimmt. Pilotanlagenmodelle können diese Beziehungen nutzen, um das Eindringverhalten zu frühen Zeitpunkten vom stationären Langzeitbetrieb zu unterscheiden und aussagekräftige Verweilzeiten sowie effektive Stoffübergangskoeffizienten zu berechnen.

Wie Ähnlichkeitslösungen den instationären Stoffübergang beschreiben

Ausgangspunkt ist das zweite Fick’sche Gesetz

Für die eindimensionale Diffusion in eine Flüssigkeitsschicht wird das Konzentrationsfeld beschrieben durch

[ \frac{\partial \Theta}{\partial t}

D\frac{\partial^2 \Theta}{\partial x^2} ]

wobei (D) der Diffusionskoeffizient und (\Theta) eine dimensionslose Konzentrationsvariable ist.

Eine häufig verwendete Definition lautet

[ \Theta=\frac{C-C_s}{C_0-C_s} ]

wobei (C_0) die anfängliche Konzentration im Volumen und (C_s) die vorgegebene Oberflächenkonzentration ist.

Diese Formulierung ist geeignet, wenn die Flüssigkeitsschicht zunächst homogen ist und sich die Grenzflächenkonzentration im Vergleich zur Diffusionszeit durch die Schicht schnell ändert.

Anwendung der Ähnlichkeitsvariable

Definiert wird

[ \eta=\frac{x}{2(Dt)^{1/2}} ]

Diese Variable kombiniert Weg und Zeit zu einer einzigen Koordinate. Sie berücksichtigt, dass die instationäre Diffusion keine Grenzschicht mit fester Dicke erzeugt; die beeinflusste Region breitet sich näherungsweise mit ((Dt)^{1/2}) aus.

Durch Einsetzen in das zweite Fick’sche Gesetz wird die partielle Differentialgleichung reduziert zu

[ \frac{d^2\Theta}{d\eta^2} + 2\eta\frac{d\Theta}{d\eta}=0 ]

Das instationäre Problem wird dadurch in eine gewöhnliche Differentialgleichung überführt.

Verwendung des Fehlerfunktionsprofils

Für ein halbunendliches Medium mit einer plötzlich aufgeprägten Oberflächenkonzentration lautet die Lösung

[ \Theta(\eta)=\operatorname{erfc}(\eta) ]

Äquivalent dazu kann die normierte Konzentrationsänderung geschrieben werden als

[ \frac{C-C_0}{C_s-C_0}

\operatorname{erfc} \left( \frac{x}{2\sqrt{Dt}} \right) ]

Dieses Profil liefert die Konzentration an jedem Ort und zu jedem Zeitpunkt, ohne die instationäre Gleichung für jede Betriebsbedingung separat lösen zu müssen.

Umwandlung des Profils in eine Eindringtiefe

Definition der technischen Eindringtiefe

Eine Pilotanlage benötigt ein praktisches Kriterium, um zu bestimmen, wann ein Flüssigkeitsbereich durch die Diffusion wesentlich beeinflusst wurde. Wenn 95 % der Konzentrationsänderung abgeschlossen sind, beträgt die verbleibende dimensionslose Abweichung ungefähr 0,05.

Für die Lösung mit der komplementären Fehlerfunktion gilt

[ \operatorname{erfc}(\eta)\approx 0.05 ]

bei ungefähr

[ \eta=1.4 ]

Daher gilt

[ \delta(t)=2\eta\sqrt{Dt} ]

und

[ \boxed{\delta(t)\approx 2.8\sqrt{Dt}} ]

Diese Eindringtiefe ist ein definierter technischer Kennwert und keine scharf abgegrenzte physikalische Grenzfläche. Konzentrationsänderungen setzen sich jenseits von (\delta) fort, sind jedoch im Verhältnis zur Differenz zwischen Oberflächen- und Volumenkonzentration gering.

Abschätzung der erforderlichen Kontaktzeit

Wenn ein Prozess erfordert, dass die Diffusionszone eine Entfernung (L) erreicht, setzt man (\delta=L):

[ L\approx 2.8\sqrt{Dt} ]

Dies ergibt

[ \boxed{t\approx \frac{L^2}{7.84D}} ]

Diese Abschätzung kann als erste Näherung für die Zeit verwendet werden, die erforderlich ist, damit die Diffusion einen gesamten Film, Membranbereich, eine poröse Schicht oder eine stagnierende Flüssigkeitszone beeinflusst.

Die Beziehung zeigt außerdem eine wichtige Auslegungsbeschränkung: Die Diffusionszeit steigt mit dem Quadrat der Entfernung. Eine Verdopplung der charakteristischen Diffusionsstrecke erfordert unter gleichen Bedingungen ungefähr die vierfache Zeit.

Verwendung der Fourier-Zahl für endliche Systeme

Für eine endliche Schicht mit der Dicke (L) definiert man die Fourier-Zahl

[ Fo=\frac{Dt}{L^2} ]

Die Näherung für ein halbunendliches Medium ist bei kurzen Zeiten am nützlichsten, wenn die Eindringtiefe deutlich kleiner als die Schichtdicke ist. Sobald die Diffusion die gegenüberliegende Begrenzung erreicht, muss die endliche Geometrie mithilfe einer Lösung für eine endliche Platte, einen Zylinder oder einer anderen geeigneten Lösung modelliert werden.

Für zylindrische Feststoffe oder Pellets ist die relevante dimensionslose Zeit analog

[ \tau=\frac{Dt}{R_0^2} ]

wobei (R_0) der charakteristische Radius ist. Lösungen für Zylinder enthalten im Allgemeinen Reihen mit Besselfunktionen anstelle des einfachen Fehlerfunktionsprofils.

Berechnung der instationären Stoffübergangsrate

Bestimmung des Oberflächenflusses mit dem ersten Fick’schen Gesetz

Der Diffusionsfluss ist

[ N_x=-D\frac{\partial C}{\partial x} ]

Durch Differenzieren des Konzentrationsprofils in Form der Fehlerfunktion und Auswertung an der Oberfläche erhält man einen Fluss proportional zu

[ \boxed{ N_x\big|_{x=0}

\left(\frac{D}{\pi t}\right)^{1/2} K C_0 } ]

für die durch (K C_0) dargestellte Konzentrations- und Gleichgewichtsbeziehung.

Die genaue Konzentrationstriebkraft hängt von den gewählten Randbedingungen ab. Bei einer Gasabsorption kann (K C_0) beispielsweise eine Konzentrationsdifferenz an der Grenzfläche oder eine gleichgewichtsangepasste Triebkraft darstellen und nicht eine universelle Konstante.

Interpretation des Verhaltens (t^{-1/2})

Der momentane Fluss ist unmittelbar nach der Beaufschlagung am größten, da der Konzentrationsgradient an der Grenzfläche steil ist. Mit dem Wachstum der Eindringzone wird der Gradient flacher, und der Fluss nimmt gemäß

[ N_x\propto t^{-1/2} ]

ab. Dieses Verhalten ist von zentraler Bedeutung für die Penetrationstheorie und den instationären Kontakt zwischen kontinuierlich erneuerten Fluidelementen.

Berechnung des mittleren Flusses über eine Kontaktperiode

Der momentane Fluss ist nicht immer die in einer Pilotanlage gemessene Größe. Wenn ein Fluidelement für eine Zeit (t_c) exponiert wird, beträgt der mittlere Fluss

[ \overline{N}_x

\frac{1}{t_c}\int_0^{t_c}N_x(t),dt ]

Für einen Fluss der Form

[ N_x(t)=A t^{-1/2} ]

ergibt sich für den Mittelwert

[ \overline{N}_x

\frac{2A}{\sqrt{t_c}} ]

Somit ist der mittlere Fluss über ein Kontaktintervall bei diesem idealisierten Diffusionsproblem in einem halbunendlichen Medium doppelt so groß wie der momentane Fluss am Ende des Intervalls.

Diese Unterscheidung ist wichtig, wenn im Labor Messungen auf Grundlage der insgesamt übertragenen Masse mit Korrelationen verglichen werden, die auf momentanen Stoffübergangskoeffizienten beruhen.

Anwendung der Methode auf Grundoperationen in Pilotanlagen

Gasabsorption und -desorption

In einem Gas-Flüssigkeits-Kontaktor kann die Grenzfläche als Oberfläche betrachtet werden, an der sich die Flüssigkeitskonzentration schnell ändert. Die Ähnlichkeitslösung schätzt ab, wie weit der gelöste Stoff während einer durch Turbulenz, Wirbelneuerung oder Grenzflächenkontakt bestimmten Expositionszeit in die Flüssigkeit eindringt.

Eine Abschätzung nach der Penetrationstheorie kann durch einen effektiven Koeffizienten ausgedrückt werden:

[ k_L\sim \left(\frac{D}{t_c}\right)^{1/2} ]

mit einem Proportionalitätsfaktor, der durch die gewählten Randbedingungen und die Definition des Flusses bestimmt wird.

Messungen der Gas- und Flüssigkeitszusammensetzungen in Pilotanlagen können anschließend mit Vorhersagen auf Grundlage von

[ N_A=k_L(C_{A,i}-C_{A,b}) ]

verglichen werden, wobei (C_{A,i}) und (C_{A,b}) die Konzentrationen des Stoffes A an der Grenzfläche beziehungsweise im Volumen der Flüssigkeit sind.

Das Modell ist besonders nützlich, um zu untersuchen, wie Rühren, Gasgeschwindigkeit, Flüssigkeitsgeschwindigkeit und Grenzflächenerneuerung die effektive Kontaktzeit und Stoffübergangsrate verändern.

Membrandiffusion

Zu kurzen Zeiten nach einem Konzentrationssprung verhält sich eine Membran näherungsweise wie eine halbunendliche Platte. Das instationäre Konzentrationsprofil lautet

[ \frac{c-c_{0}}{c_{\infty}-c_{0}}

\operatorname{erf} \left( \frac{z}{\sqrt{4Dt}} \right) ]

Diese Form ist äquivalent zur Darstellung mit der komplementären Fehlerfunktion, wobei die Konzentrationsnormierung anders gewählt ist.

Die Eindringtiefe hilft bei der Bestimmung, ob der gelöste Stoff die gegenüberliegende Membranoberfläche erreicht hat. Bei längeren Zeiten nimmt die gesamte Membrandicke am Transport teil, und das System nähert sich dem stationären Verhalten eines Dünnfilms.

Das geeignete Diffusionsmodell hängt auch von der Membranmorphologie ab. Dichte Membranen werden im Allgemeinen mit molekularer Diffusion modelliert, während poröse Strukturen möglicherweise Beschreibungen durch konvektive Diffusion oder Knudsen-Diffusion erfordern.

Flüssig-Flüssig-Extraktion

Bei der Flüssig-Flüssig-Extraktion kann die instationäre Diffusion das Eindringen eines gelösten Stoffes in dispergierte Tropfen oder die Ausbildung von Konzentrationsgrenzschichten in der Nähe der Grenzfläche beschreiben. Die charakteristische Kontaktzeit bestimmt, ob sich der Tropfen wie ein halbunendliches Medium verhält oder ob die Diffusion bereits sein Zentrum erreicht hat.

Für einen Tropfen mit dem Radius (R) vergleicht man die Kontaktzeit mit

[ t_D\sim\frac{R^2}{D} ]

Kurze Kontaktzeiten begünstigen ein Verhalten nach dem Penetrationstyp. Längere Kontaktzeiten erfordern Modelle für endliche Tropfen und können innere Zirkulation, Grenzflächenwiderstand und sich ändernde Gleichgewichtsbedingungen einschließen.

Feststoff-Flüssigkeits-Extraktion und Auslaugung

Bei einer Feststoffmatrix, die mit einer gut durchmischten Flüssigkeit kontaktiert wird, ändert sich die Flüssigkeitskonzentration, während der gelöste Stoff den Feststoff verlässt. Die Diffusion im Feststoff muss mit einer Stoffbilanz der Flüssigphase gekoppelt werden.

Die halbunendliche Ähnlichkeitslösung ist zu frühen Zeiten nützlich, doch sobald sich die Eindringfront dem Zentrum oder der gegenüberliegenden Begrenzung nähert, ist ein Modell für eine endliche Platte oder einen Zylinder erforderlich. Das resultierende Modell kann Extraktionszeit, Badkonzentration und den Einfluss der Partikelgröße abschätzen.

Adsorption und Imprägnierung

Die instationäre Diffusion beschreibt auch die Bewegung gelöster Stoffe in poröse Partikel, Katalysatorpellets und andere Feststoffmatrizen. Die relevante Geometrie kann planar, zylindrisch oder sphärisch sein.

Wenn der Feststoff im Verhältnis zur Eindringtiefe groß ist, bietet die Lösung für ein halbunendliches Medium eine zweckmäßige Näherung für frühe Zeiten. Wenn das Partikel klein ist oder die Exposition lange dauert, müssen innere Diffusion und Gleichgewichtsaufnahme gemeinsam modelliert werden.

Verknüpfung der Diffusion mit Messungen in Pilotanlagen

Bestimmung eines effektiven Stoffübergangskoeffizienten

Pilotanlagen geben den Stoffübergang typischerweise durch Koeffizienten und nicht allein durch Konzentrationsprofile an. Aus dem gemessenen Fluss kann ein effektiver Koeffizient abgeleitet werden:

[ k_{\text{eff}}

\frac{N_A}{\Delta C_{\text{driving force}}} ]

Für eine instationäre Penetration gilt

[ k_{\text{eff}}\propto \sqrt{\frac{D}{t_c}} ]

Dies verknüpft die molekulare Diffusivität und die Expositionszeit mit der experimentell beobachteten Stoffübergangsrate.

Das Ergebnis ist kein universeller Anlagenkoeffizient. Es handelt sich um einen effektiven Wert, der auch Hydrodynamik, Grenzflächenfläche, Phaseneigenschaften und den gewählten Mittelungszeitraum widerspiegelt.

Kombination des lokalen Stoffübergangs mit der Grenzflächenfläche

Die gesamte Stoffübergangsrate wird üblicherweise geschrieben als

[ \dot n_A

k_{\text{eff}},a,V,\Delta C ]

wobei (a) die Grenzflächenfläche pro Volumeneinheit und (V) das kontaktierte Volumen ist.

Eine hohe Diffusivität garantiert keine hohe Anlagenleistung, wenn die Pilotanlage eine unzureichende Grenzflächenfläche oder eine schlechte Phasendispersion aufweist. Umgekehrt kann eine verstärkte Durchmischung den Stoffübergang verbessern, indem sie die Grenzfläche erneuert und die effektive Expositionszeit einzelner Fluidelemente verkürzt.

Validierung des Modells mithilfe von Zusammensetzungsdaten

Bei der Gasabsorption oder Extraktion werden die Einlass- und Auslasszusammensetzungen gemessen und die übertragene Menge aus einer Phasenbilanz berechnet. Vergleichen Sie diese gemessene Rate mit der integrierten instationären Vorhersage.

Abweichungen können darauf hindeuten, dass die Annahme eines halbunendlichen Mediums ungültig ist, die Grenzflächenkonzentration nicht konstant ist, sich die Diffusivität ändert oder hydrodynamische Effekte das ideale Diffusionsmodell dominieren.

Verständnis der Zielkonflikte

Die Annahme eines halbunendlichen Mediums gilt nur in einem begrenzten Bereich

Die Lösung mit der Fehlerfunktion setzt voraus, dass die diffundierende Region noch keine zweite Begrenzung erreicht hat. Sie wird unzuverlässig, wenn (\delta) sich der Filmdicke, Membrandicke, dem Partikelradius oder dem Tropfenradius nähert.

In diesem Stadium sollten Lösungen für endliche Gebiete oder numerische Lösungen mit den tatsächlichen Randbedingungen verwendet werden.

Eine konstante Diffusivität kann unrealistisch sein

Das Grundmodell nimmt an, dass (D) konstant ist. Temperaturgradienten, konzentrationsabhängige Viskosität, Quellung, Wechselwirkungen in Mehrkomponentensystemen und Änderungen der Membranstruktur können diese Annahme ungültig machen.

Für mehrkomponentige Gasgemische können Stefan-Maxwell-Gleichungen erforderlich sein. Für verdünnte Flüssigkeitslösungen kann eine pseudobinäre Näherung ausreichen, da das Lösungsmittel das Transportverhalten dominiert.

Ideale Randbedingungen entsprechen möglicherweise nicht der Anlage

Eine plötzlich aufgeprägte, konstante Oberflächenkonzentration ist eine nützliche mathematische Bedingung, doch reale Pilotanlagen können einen endlichen Grenzflächenwiderstand und sich ändernde Konzentrationen im Volumen aufweisen. Die Grenzfläche kann außerdem kontinuierlich erneuert werden, anstatt stationär zu sein.

Die Filmtheorie ist für relativ stationäre Grenzschichten nützlich, während sich die Penetrationstheorie besser für turbulente oder instationäre Grenzflächen eignet. Die Wahl der falschen Theorie kann zu einem plausiblen, aber irreführenden Stoffübergangskoeffizienten führen.

Geometrie und Anisotropie beeinflussen die Maßstabsübertragung

Planare Diffusion ist nicht mit radialer Diffusion in einem Pellet, Tropfen oder einer zylindrischen Membran gleichzusetzen. Die charakteristische Länge muss der dominierenden Transportrichtung entsprechen.

Anisotrope Materialien erfordern richtungsspezifische Diffusivitäten und Randbedingungen. Eine Skalierung, die nur auf der geometrischen Größe basiert, kann daher die tatsächliche Transportzeit falsch darstellen.

Höhere Stoffübergangsraten können die Kontaktzeit verkürzen

Da der instationäre Fluss umgekehrt proportional zur Quadratwurzel der Expositionszeit skaliert, führen kürzere Expositionszeiten im idealen Penetrationsmodell zu höheren momentanen Flüssen. Ein kurzer Kontakt kann jedoch auch die insgesamt übertragene Menge begrenzen.

Bei der Optimierung einer Pilotanlage müssen daher Fluss, Verweilzeit, Grenzflächenfläche, Phasendurchsatz und die erforderliche Auslassspezifikation gegeneinander abgewogen werden.

Die richtige Wahl für Ihr Ziel treffen

Verwenden Sie das instationäre Modell als Berechnung aus ersten Prinzipien und verfeinern Sie es anschließend anhand der Geometrie, Hydrodynamik und des Phasenverhaltens der tatsächlichen Pilotanlage.

Wenn Ihr Hauptaugenmerk auf der Abschätzung des Stoffübergangs zu frühen Zeiten liegt: Verwenden Sie (\eta=x/[2\sqrt{Dt}]), das Konzentrationsprofil mit der Fehlerfunktion und (\delta\approx2.8\sqrt{Dt}), um den beeinflussten Bereich und den momentanen Fluss zu berechnen.
Wenn Ihr Hauptaugenmerk auf der Wahl der Verweilzeit liegt: Vergleichen Sie die erforderliche Diffusionsstrecke mit (\delta(t)) oder schätzen Sie (t\approx L^2/(7.84D)) für eine Penetration von 95 % ab.
Wenn Ihr Hauptaugenmerk auf Gasabsorption oder Flüssig-Flüssig-Extraktion liegt: Kombinieren Sie die auf der Penetration basierende Schätzung (k_{\text{eff}}\sim(D/t_c)^{1/2}) mit Grenzflächenfläche, Gleichgewichtsdaten und gemessenen Phasenbilanzen.
Wenn Ihr Hauptaugenmerk auf der Membranleistung liegt: Verwenden Sie zu kurzen Zeiten die Lösung für ein halbunendliches Medium und wechseln Sie zu einem stationären Modell für eine endliche Platte, sobald die Diffusion die gegenüberliegende Membranoberfläche erreicht.
Wenn Ihr Hauptaugenmerk auf der Maßstabsübertragung liegt: Bewahren Sie die relevanten dimensionslosen Gruppen, insbesondere (Fo=Dt/L^2), und bewerten Sie gleichzeitig Strömungsregime, Temperatur, Druck, Grenzflächenfläche und Mehrkomponenteneffekte.
Wenn Ihr Hauptaugenmerk auf der Modellgenauigkeit liegt: Ersetzen Sie die Ähnlichkeitslösung durch Modelle für endliche Geometrien, Mehrkomponentensysteme oder numerische Modelle, sobald Begrenzungen, Konzentrationsänderungen, Anisotropie oder Grenzflächenwiderstand relevant werden.

Prinzipien der instationären Diffusion verwandeln Konzentrationsdaten aus Pilotanlagen in verständliche Vorhersagen zu Penetration, Verweilzeit und Stoffübergangsrate.

Übersichtstabelle:

Schlüsselkonzept	Technische Anwendung
Ähnlichkeitsvariable	Modelliert instationäre Konzentrationsprofile mit (\eta=x/[2\sqrt{Dt}])
Lösung mit der Fehlerfunktion	Prognostiziert Konzentrationsänderungen in halbunendlichen Schichten, Filmen, Membranen und Partikeln
Eindringtiefe	Schätzt die beeinflusste Entfernung mit (\delta\approx2.8\sqrt{Dt}) für einen Fortschritt von 95 % ab
Diffusionszeit	Nähert die erforderliche Kontaktzeit mit (t\approx L^2/(7.84D)) an
Instationärer Fluss	Zeigt, dass der momentane Stoffübergangsfluss mit (t^{-1/2}) abnimmt
Effektiver Koeffizient	Verknüpft die Stoffübergangsleistung mit Diffusivität und Kontaktzeit: (k_{eff}\propto\sqrt{D/t_c})
Prüfung der Maßstabsübertragung	Verwendet Fourier-Zahl, Geometrie, Hydrodynamik, Grenzflächenfläche und Phasenverhalten
Modellgrenzen	Erfordert Modelle für endliche Gebiete, Mehrkomponentensysteme oder numerische Modelle, wenn Begrenzungen und Nichtidealitäten relevant sind

Bessere Stoffübergangssysteme für Pilotanlagen entwickeln

Sie möchten die Diffusionstheorie in eine zuverlässige Anlagenleistung überführen? Kontaktieren Sie LABPARK bezüglich Bildungs- und Berufsbildungs-Pilotanlagen für Grundoperationen, die für Chemieingenieurwesen, Bioprozesstechnik und Biotechnologie sowie Umwelt- und Wasseraufbereitungsanwendungen entwickelt wurden.

LABPARK unterstützt Universitäten, Forschungsinstitute und Unternehmen mit Systemen im Pilotmaßstab zur Untersuchung von Absorption, Desorption, Membrantransport, Flüssig-Flüssig-Extraktion, Adsorption, Auslaugung und verwandten Grundoperationen. Unsere Anlagen helfen Ihnen, Konzentrationsprofile, Einflüsse der Verweilzeit, instationäre Flüsse, effektive Stoffübergangskoeffizienten und das Maßstabsübertragungsverhalten unter realistischen Betriebsbedingungen zu messen.

Arbeiten Sie mit LABPARK zusammen, um Ihren Prozess zu bewerten, Diffusions- und hydrodynamische Modelle zu validieren, die Bedienerschulung zu verbessern und eine Pilotanlagenkonfiguration auszuwählen, die zuverlässige Forschung, Lehre und kommerzielle Prozessentwicklung unterstützt. Besprechen Sie Ihre Anwendung mit LABPARK.