Verwenden Sie eine Wärmebilanz mit beweglicher Grenzfläche: Die Pseudostationär-Näherung behandelt die Wärmeleitung durch die bestehende Feststoffschicht als instantan im Vergleich zur langsameren Bewegung der Fest-Flüssig-Grenzfläche. Die Feststoffwachstumsrate wird dann ermittelt, indem die an der Grenzfläche freigesetzte latente Wärme mit der durch die Feststoffschicht, die Rohrwand und den internen Kühlmittelfilm übertragenen Wärme gleichgesetzt wird. Mit zunehmender Schichtdicke steigt deren thermischer Widerstand, wodurch die Wachstumsrate sinkt.

Die praktische Berechnung erfolgt über eine widerstandsbasierte Stefan-Bilanz: Bestimmen Sie die Wärmestromrate von der Phasenwechseltemperatur zur Kühlmittel-Mischtemperatur und teilen Sie diese Rate durch die latente Wärme, die zur Bildung des zusätzlichen Feststoffs erforderlich ist.

Annahmen der Pseudostationär-Näherung

Thermische Profile passen sich schneller an als die Grenzfläche

Das Temperaturfeld im Feststoff wird zu jedem Zeitpunkt als stationär angenähert, auch wenn sich die Fest-Flüssig-Grenzfläche bewegt. Dies ist sinnvoll, wenn die thermische Diffusion durch den Feststoff wesentlich schneller erfolgt als die Zeit, die für eine nennenswerte Änderung der Schichtdicke benötigt wird.

Die Näherung bedeutet nicht, dass der Prozess tatsächlich stationär ist. Die Feststoffdicke und der Wärmeübertragungswiderstand variieren weiterhin mit der Zeit.

Latente Wärme dominiert die Bilanz der beweglichen Grenzfläche

Bei Gefrier- oder Kristallisationsprozessen ist die latente Schmelz- oder Kristallisationswärme der wesentliche Energiebedarf im Zusammenhang mit der Grenzflächenbewegung. Die fühlbare Wärmeaufnahme im Feststoff und in der Flüssigkeit wird vernachlässigt oder als sekundäre Korrektur behandelt.

Die Grenzflächenbilanz basiert daher auf

[ \text{abgeführte Wärme}=\text{durch neuen Feststoff freigesetzte latente Wärme}. ]

Die Grenzfläche bleibt nahe der Phasenwechseltemperatur

Der Fest-Flüssig-Grenzfläche wird üblicherweise die Gleichgewichts-Phasenwechseltemperatur (T_m) oder eine effektive Kristallisationstemperatur zugewiesen, wenn Konzentrations- oder Druckeffekte wichtig sind.

Dies liefert die thermische Triebkraft für den Wärmeabzug durch das Rohr und die Rohranordnung.

Aufbau des thermischen Widerstandsmodells

Definition der Rohr- und Feststoffgeometrie

Sei

(R_i) der Innenradius des Rohrs,
(R_t) der Außenradius des Rohrs,
(R_s(t)) der Radius der beweglichen Fest-Flüssig-Grenzfläche,
(k_w) die Wärmeleitfähigkeit der Rohrwand,
(k_s) die Wärmeleitfähigkeit der Feststoffschicht,
(h_i) der interne Wärmeübergangskoeffizient auf der Kühlmittelseite,
(T_c) die Mischtemperatur des Kühlmittels,
(\rho_s) die Feststoffdichte,
(L) die latente Wärme der Erstarrung oder Kristallisation.

Für eine Ablagerung, die von der Rohroberfläche nach außen wächst, nimmt der Feststoff den Bereich zwischen (R_t) und (R_s) ein, und seine Dicke beträgt

[ \delta(t)=R_s(t)-R_t. ]

Ausdruck des Widerstands pro Rohrlängeneinheit

Für ein zylindrisches Rohr sind die Widerstände des Kühlmittelfilms, der Rohrwand und der Feststoffschicht pro Längeneinheit:

[ \mathcal{R}_i'=\frac{1}{h_i,2\pi R_i}, ]

[ \mathcal{R}_w'=\frac{\ln(R_t/R_i)}{2\pi k_w}, ]

und

[ \mathcal{R}_s'=\frac{\ln(R_s/R_t)}{2\pi k_s}. ]

Zusätzliche Widerstände können bei Bedarf einbezogen werden, wie z. B. ein Kontaktwiderstand an der Grenzfläche zwischen Rohr und Feststoff oder ein flüssigkeitsseitiger Widerstand außerhalb der wachsenden Schicht.

Der Gesamtwiderstand pro Längeneinheit ist daher

[ \mathcal{R}_{\text{tot}}'

\frac{1}{h_i,2\pi R_i} + \frac{\ln(R_t/R_i)}{2\pi k_w} + \frac{\ln(R_s/R_t)}{2\pi k_s}. ]

Berechnung der instantanen Wärmeabfuhrrate

Unter der Pseudostationär-Näherung ist die abgeführte Wärme pro Rohrlängeneinheit

[ q'(R_s)

\frac{T_m-T_c}{\mathcal{R}_{\text{tot}}'}. ]

Diese Wärmerate wird unter Verwendung des aktuellen Grenzflächenradius berechnet. Folglich ist die Berechnung quasistationär, erfasst aber dennoch die Zeitabhängigkeit, die durch das Schichtwachstum verursacht wird.

Ableitung der Feststoffwachstumsrate

Aufstellen der Bilanz der latenten Wärme

Wenn sich die Grenzfläche um (dR_s) nach außen bewegt, beträgt das neu gebildete Feststoffvolumen pro Rohrlängeneinheit

[ dV'=2\pi R_s,dR_s. ]

Die entsprechende freigesetzte latente Wärme ist

[ dQ_{\text{latent}}'

\rho_s L,2\pi R_s,dR_s. ]

Die Division durch die Zeit und die Gleichsetzung dieser Rate mit der abgeführten Wärme ergibt

[ q'(R_s)

2\pi R_s\rho_s L\frac{dR_s}{dt}. ]

Das Einsetzen des Widerstandsausdrucks führt zur maßgeblichen Wachstumsgleichung:

[ \boxed{ \frac{dR_s}{dt}

\frac{T_m-T_c} { 2\pi R_s\rho_s L \left[ \frac{1}{h_i,2\pi R_i} + \frac{\ln(R_t/R_i)}{2\pi k_w} + \frac{\ln(R_s/R_t)}{2\pi k_s} \right] } } ]

Diese Gleichung liefert die instantane radiale Wachstumsrate der Feststoffschicht.

Interpretation des Rückgangs der Wachstumsrate

Wenn (R_s) zunimmt, steigt auch der Feststoffwiderstand

[ \frac{\ln(R_s/R_t)}{2\pi k_s} ]

an. Der Wärmestrom, der die Phasengrenze erreicht, sinkt daher, und die Grenzfläche schreitet langsamer voran.

Der Faktor (2\pi R_s) im Term der latenten Wärme ist ebenfalls von Bedeutung: Ein radialer Zuwachs erzeugt mehr Feststoffvolumen pro Rohrlängeneinheit, je größer der Umfang der Grenzfläche wird.

Nützliche Vereinfachungen für technische Schätzungen

Verwendung einer planaren Näherung für eine dünne Schicht

Wenn die Feststoffdicke im Vergleich zum Rohrradius gering ist, (\delta \ll R_t), kann der zylindrische Widerstand durch einen planaren Widerstand angenähert werden. Die Wachstumsgleichung wird zu

[ \frac{d\delta}{dt}

\frac{T_m-T_c} {\rho_s L \left[ \frac{1}{h_i} + \frac{t_w}{k_w} + \frac{\delta}{k_s} \right]}, ]

wobei (t_w) die Rohrwanddicke ist.

Wenn die Widerstände des Kühlmittelfilms und der Wand zu einem konstanten Widerstand (R_0) zusammengefasst werden, dann gilt

[ \frac{d\delta}{dt}

\frac{T_m-T_c} {\rho_s L\left(R_0+\delta/k_s\right)}. ]

Integration zur Abschätzung der Dicke über die Zeit

Bei konstanten Stoffwerten und konstanter Temperaturdifferenz ergibt die Integration

[ t

\frac{\rho_s L}{T_m-T_c} \left[ R_0\delta+\frac{\delta^2}{2k_s} \right]. ]

Dieses Ergebnis zeigt zwei Grenzbereiche.

Bei geringer Dicke dominieren die Widerstände des Kühlmittelfilms und der Wand, sodass das Wachstum zeitlich annähernd linear verläuft. Bei größerer Dicke dominiert die Wärmeleitung im Feststoff, und das bekannte diffusionsähnliche Verhalten tritt auf:

[ \delta \propto \sqrt{t}. ]

Verwendung der zylindrischen Form für erhebliche Ablagerungen

Wenn die Schichtdicke im Verhältnis zum Rohrradius nicht klein ist, kann die planare Näherung die sich ändernde Geometrie und den Widerstand unterschätzen. Die zylindrische Gleichung sollte dann numerisch oder nach Entwicklung der Widerstandsterme analytisch integriert werden.

Für eine rohrbasierte Pilotanlage ist diese Unterscheidung wichtig, da Ablagerungen so dick werden können, dass sowohl die lokale Wärmeübertragungsfläche als auch der logarithmische Leitungswiderstand das Ergebnis beeinflussen.

Anwendung des Modells auf eine rohrbasierte Prozesseinheit

Abschätzung des effektiven Wärmeübergangskoeffizienten

Das Widerstandsmodell identifiziert, welcher Teil des Systems das Wachstum begrenzt. Der interne Kühlmittelfilm wird durch (1/h_i) dargestellt, während die wachsende Feststoffschicht im planaren Grenzfall (\delta/k_s) beiträgt.

Ein effektiver Gesamtwert kann aus dem Gesamtwiderstand definiert werden, bezogen auf eine konsistente Referenzfläche:

[ U=\frac{1}{R_{\text{total}}}. ]

Dies ermöglicht den Vergleich gemessener Wachstumsraten mit berechneten Werten und kann aufzeigen, ob die dominierende Begrenzung die kühlmittelseitige Konvektion, die Rohrwandleitung oder die Feststoffschichtleitung ist.

Bewertung der Kühlmittelbetriebsbedingungen

Die Kühlmitteltemperatur und die Durchflussrate beeinflussen die thermische Triebkraft und den internen Koeffizienten (h_i). Eine Senkung von (T_c) erhöht im Allgemeinen die Triebkraft (T_m-T_c), während eine Erhöhung der Kühlmittelgeschwindigkeit den internen Filmwiderstand verringern kann.

Das Modell kann daher verwendet werden, um vorherzusagen, wie sich Änderungen der Kühlmittelzufuhrbedingungen auf die anfängliche Wachstumsrate und die Zeit bis zum Erreichen einer Zielschichtdicke auswirken.

Berücksichtigung von Temperatur- und Zusammensetzungseffekten

In Kristallisationssystemen kann (T_m) mit der Konzentration, dem Druck oder der Ausscheidung von gelösten Stoffen an der Grenzfläche variieren. Wenn die Flüssigkeit zunehmend konzentriert wird, können sich die Phasenwechseltemperatur und die verfügbare Triebkraft während des Betriebs ändern.

In diesem Fall sollte die Bilanz als Funktion von Ort, Zeit oder lokaler Zusammensetzung aktualisiert werden, anstatt ein einzelnes konstantes (T_m) zu verwenden.

Kopplung des Modells mit dem axialen Rohrverhalten

Eine reale Wärmeübertragungseinheit kann eine Kühlmitteltemperatur aufweisen, die sich entlang des Rohrs ändert. Die lokale Triebkraft ist dann

[ \Delta T(z,t)=T_m(z,t)-T_c(z,t), ]

und die Wachstumsgleichung muss sowohl entlang der Rohrlänge als auch über die Zeit gelöst werden.

Dies ergibt eine lokale Dicke (R_s(z,t)), mit der Bereiche identifiziert werden können, in denen übermäßige Ablagerungen, Blockaden oder ungleichmäßige Kristallisation wahrscheinlich sind.

Verständnis der Kompromisse

Die Näherung vernachlässigt transiente fühlbare Erwärmung

Das Pseudostationär-Modell vernachlässigt die Zeit, die zum Abkühlen von neu gebildetem Feststoff und nahegelegener Flüssigkeit benötigt wird. Dies ist am ehesten vertretbar, wenn die Stefan-Zahl klein ist, was bedeutet, dass die latente Wärme im Vergleich zur fühlbaren Wärme, die mit der Temperaturänderung verbunden ist, groß ist.

Wenn fühlbare Effekte signifikant sind, ist ein vollständiges transientes Modell mit beweglicher Grenzfläche angemessener.

Konstante Stoffwerte können unrealistisch sein

Wärmeleitfähigkeit, Dichte, latente Wärme, Viskosität und Wärmeübergangskoeffizienten können mit Temperatur und Zusammensetzung variieren. Die Annahme, dass sie konstant sind, vereinfacht die Analyse, kann aber die vorhergesagte Dicke und Wachstumsraten verzerren.

Die Variation der Stoffwerte sollte einbezogen werden, wenn der Betriebsbereich groß ist oder wenn experimentelle Daten systematische Abweichungen vom einfachen Modell zeigen.

Oberflächenbedingungen können das Grenzflächenverhalten verändern

Das Modell geht im Allgemeinen von einer wohldefinierten Grenzfläche und einem eindimensionalen radialen Wärmestrom aus. Rauheit, Fouling, Keimbildungsverhalten, natürliche Konvektion, Scherung und Kristallmorphologie können diese Annahmen verletzen.

Bei der Kristallisation kann die Ablagerung auch porös oder dendritisch statt einer dichten Feststoffschicht sein. Ihre effektive Leitfähigkeit und Grenzflächenfläche können dann erheblich von idealen Werten abweichen.

Wachstum kann die nutzbare Prozesskapazität verringern

Eine dickere Feststoffschicht erhöht den thermischen Widerstand und kann die effektive Wärmeübertragungsleistung der Einheit verringern. Sie kann auch Strömungskanäle verengen, den Druckverlust erhöhen oder durch übermäßige Verweilzeit oder ungleichmäßige Kristallgröße zu Änderungen der Produktqualität führen.

Daher ist die Maximierung der instantanen Wachstumsrate nicht immer gleichbedeutend mit der Optimierung des gesamten Einheitenbetriebs.

Die richtige Wahl für Ihr Ziel

Verwenden Sie das Pseudostationär-Modell als erstes technisches Werkzeug und erhöhen Sie die Modellkomplexität nur dann, wenn die Annahmen das Ergebnis wesentlich beeinflussen.

Wenn Ihr Hauptaugenmerk auf der Schätzung der Feststoffdicke über die Zeit liegt: Verwenden Sie die zylindrische Widerstandsbilanz bei erheblichen Ablagerungen und die planare integrierte Gleichung für Schätzungen bei dünnen Schichten.
Wenn Ihr Hauptaugenmerk auf der Identifizierung der Wärmeübertragungsbegrenzung liegt: Trennen Sie die Widerstände von Kühlmittelfilm, Rohrwand und Feststoffschicht, um zu bestimmen, welche Betriebsvariable den größten Effekt hat.
Wenn Ihr Hauptaugenmerk auf der Optimierung des Kühlmittelbetriebs liegt: Bewerten Sie Kühlmitteltemperatur und Durchfluss gemeinsam, da Triebkraft und interne Konvektion beide die Abfuhrrate der latenten Wärme steuern.
Wenn Ihr Hauptaugenmerk auf der Vorhersage von ungleichmäßigem Wachstum in einer Pilotanlage liegt: Koppeln Sie die lokale Wachstumsgleichung mit den Bilanzen für axiale Kühlmitteltemperatur und Zusammensetzung.
Wenn Ihr Hauptaugenmerk auf einer hochpräzisen Vorhersage liegt: Beziehen Sie transiente fühlbare Wärme, variable Stoffwerte, externen Flüssigkeitswiderstand und Ablagerungsmorphologie ein, wenn das experimentelle Verhalten nicht durch das Pseudostationär-Modell erklärt werden kann.

Eine pseudostationäre Widerstandsbilanz wandelt eine bewegliche Erstarrungsgrenzfläche in eine praktische Wärmeübertragungsberechnung um und macht gleichzeitig deutlich, warum sich das Wachstum verlangsamt, wenn die Feststoffschicht dicker wird.

Zusammenfassungstabelle:

Schlüsselelement	Technische Erkenntnis
Kernannahme	Thermische Profile passen sich schneller an als die Fest-Flüssig-Grenzfläche.
Triebkraft	Phasenwechseltemperatur minus Kühlmittel-Mischtemperatur, ΔT = Tm − Tc.
Wärmeübertragungspfad	Interner Kühlmittelfilm, Rohrwand und wachsende Feststoffschicht.
Wachstumsbilanz	Abgeführte Wärme entspricht der durch neu gebildeten Feststoff freigesetzten latenten Wärme.
Zylindrisches Wachstum	Verwenden Sie den aktuellen Grenzflächenradius und den logarithmischen Feststoffwiderstand für dicke Ablagerungen.
Schätzung dünner Schichten	Das Wachstum ist anfangs nahezu linear und nähert sich dann δ ∝ √t an, wenn der Feststoffwiderstand steigt.
Hauptsteuergrößen	Kühlmitteltemperatur, Kühlmittelgeschwindigkeit, thermische Eigenschaften, latente Wärme und Ablagerungsdicke.
Modellbeschränkungen	Transiente fühlbare Erwärmung, variable Stoffwerte, Zusammensetzungseffekte, Fouling und poröse Ablagerungen erfordern möglicherweise eine höhere Modellgenauigkeit.

Vorhersage und Verbesserung Ihres rohrbasierten Prozesses

Müssen Sie Gefrieren, Kristallisation oder Feststoffablagerungen in einer Pilotanlage steuern? Kontaktieren Sie LABPARK, um Ihre Anwendung mit unserem Team zu besprechen.

LABPARK bietet Pilotanlagen für die verfahrenstechnische Ausbildung und Berufsbildung für Universitäten, Forschungsinstitute und Unternehmen, einschließlich Lösungen für Verfahrenstechnik, Bioprozess- und Biotechnologie sowie Umwelt- und Wasseraufbereitung. Wir können Ihnen helfen, Kühlmittelbedingungen, Wärmeübertragungswiderstand, Grenzflächenwachstum, axiale Temperaturvariation und Scale-up-Anforderungen zu bewerten, damit Sie:

Die Dicke der Feststoffschicht und die Wachstumszeit zuverlässiger abschätzen können.
Identifizieren können, ob Kühlmitteldurchfluss, Rohrwandleitung oder Ablagerungswiderstand die Leistung begrenzen.
Risiken durch Blockaden, Druckverlust und ungleichmäßiges Wachstum vor dem Pilotbetrieb bewerten können.
Praktische Instrumentierung und Betriebsbedingungen für Lehre, Forschung und Produktentwicklung auswählen können.

Kontaktieren Sie LABPARK, um Ihr Zielmaterial, die Rohrgeometrie, die Kühlmittelbedingungen und das Prozessziel zu definieren und einen auf Ihre Bedürfnisse zugeschnittenen Pilotanlagenansatz zu entwickeln.