Dynamische Energiebilanzen für einen Rührkesselheizer werden erstellt, indem die stationäre Bilanz von der transienten Bilanz subtrahiert wird. Für einen ideal durchmischten Tank mit konstantem Volumen, konstantem Durchfluss und konstanten Stoffeigenschaften lautet das Temperaturmodell

[ \frac{dT'}{dt}+aT'=\frac{1}{\tau}T_i'+KT_{st}' ]

wobei (T'), (T_i') und (T_{st}') die Abweichungen der Tanktemperatur, der Einlasstemperatur und der Heizleistung von ihren nominalen stationären Werten darstellen. Diese Form isoliert nur die dynamischen Änderungen, die für die Reaktionsanalyse und den Regelungsentwurf relevant sind.

Kernpunkt: Definieren Sie jede Variable als ihre Abweichung von einem stationären Arbeitspunkt, setzen Sie diese Definitionen in die dynamische Energiebilanz ein und verwenden Sie die stationäre Gleichung, um die konstanten Terme zu eliminieren. Das resultierende Abweichungsmodell beschreibt das transiente Temperaturverhalten, ohne die nominalen Arbeitspunkt-Offsets mitzuführen.

Beginnen Sie mit der transienten Energiebilanz

Definieren Sie das Kontrollvolumen und die Annahmen

Für einen ideal durchmischten Rührkesselheizer wird angenommen, dass die Flüssigkeitstemperatur im gesamten Tank einheitlich ist. Ein gängiges Modell für Pilotanlagen geht von konstantem Flüssigkeitsvolumen (V), konstantem Volumenstrom (F), konstanter Dichte (\rho) und konstanter Wärmekapazität (c_p) aus.

Unter diesen Annahmen werden kinetische Energie, potenzielle Energie, Wellenarbeit und Volumenänderungseffekte vernachlässigt. Die dominierenden Terme sind die Energieakkumulation, der Enthalpiestrom am Ein- und Ausgang sowie die extern zugeführte Wärme.

Schreiben Sie die dimensionale Energiebilanz

Für einen nicht-reaktiven Flüssigkeitstank mit Wärmezufuhr (Q) gilt:

[ \rho c_p V\frac{dT}{dt}

\rho c_p F(T_i-T)+Q ]

wobei:

(T) die Temperatur der Tankflüssigkeit ist,
(T_i) die Einlasstemperatur ist,
(F) der Volumenstrom ist,
(V) das Flüssigkeitsvolumen ist,
(Q) die Heizleistung ist.

Die Auslasstemperatur ist (T), da der Tank perfekt durchmischt ist.

Bringen Sie die Bilanz in die dynamische Standardform

Die Division durch (\rho c_p V) ergibt

[ \frac{dT}{dt} + \frac{F}{V}T

\frac{F}{V}T_i + \frac{Q}{\rho c_p V} ]

Der Koeffizient (F/V) ist der Kehrwert der Verweilzeit. Mit der Definition

[ \tau=\frac{V}{F} ]

erhält man

[ \frac{dT}{dt} + \frac{1}{\tau}T

\frac{1}{\tau}T_i + \frac{Q}{\rho c_p V} ]

Wenn der Heizeingang durch eine Variable (T_{st}) dargestellt wird, kann der Wärmezufuhr-Term stattdessen als (KT_{st}) geschrieben werden. Das resultierende Modell wird üblicherweise ausgedrückt als

[ \frac{dT}{dt}+aT=\frac{1}{\tau}T_i+KT_{st} ]

wobei (a), (1/\tau) und (K) die relevanten Parameter für Durchfluss, Volumen, Wärmekapazität und Wärmeübertragung zusammenfassen.

Konvertieren Sie das Modell in Abweichungsvariablen

Definieren Sie den nominalen Arbeitspunkt

Seien die stationären Betriebsbedingungen

[ T_s,\qquad T_{i,s},\qquad T_{st,s} ]

Im stationären Zustand ändert sich die Tanktemperatur nicht, daher gilt

[ \frac{dT_s}{dt}=0 ]

Die stationäre Gleichung lautet daher

[ aT_s=\frac{1}{\tau}T_{i,s}+KT_{st,s} ]

Diese Gleichung beschreibt den Arbeitspunkt, um den herum das dynamische Modell entwickelt wird.

Definieren Sie Abweichungen vom stationären Zustand

Die Abweichungsvariablen sind

[ T'=T-T_s ]

[ T_i'=T_i-T_{i,s} ]

[ T_{st}'=T_{st}-T_{st,s} ]

Äquivalent dazu gilt:

[ T=T_s+T' ]

[ T_i=T_{i,s}+T_i' ]

[ T_{st}=T_{st,s}+T_{st}' ]

Da (T_s) konstant ist, gilt

[ \frac{dT}{dt}=\frac{dT'}{dt} ]

Einsetzen in die dynamische Bilanz

Ausgehend von

[ \frac{dT}{dt}+aT=\frac{1}{\tau}T_i+KT_{st} ]

ergibt die Substitution

[ \frac{dT'}{dt} +a(T_s+T')

\frac{1}{\tau}(T_{i,s}+T_i') + K(T_{st,s}+T_{st}') ]

Durch Umstellen erhält man

[ \frac{dT'}{dt}+aT' + aT_s

\frac{1}{\tau}T_{i,s} + \frac{1}{\tau}T_i' + KT_{st,s} + KT_{st}' ]

Die konstanten Terme erfüllen die stationäre Beziehung

[ aT_s=\frac{1}{\tau}T_{i,s}+KT_{st,s} ]

und heben sich daher auf. Die Gleichung in Abweichungsvariablen lautet

[ \boxed{ \frac{dT'}{dt}+aT'

\frac{1}{\tau}T_i' + KT_{st}' } ]

Dies ist die gewünschte dynamische Energiebilanz in Abweichungsform.

Interpretation des Abweichungsmodells

Trennen Sie Stellgrößen und Störgrößen

Der Heizeingang (T_{st}') oder äquivalent (Q') ist typischerweise die Stellgrößenabweichung. Die Einlasstemperaturabweichung (T_i') wird üblicherweise als Störgröße behandelt.

Die Tanktemperaturabweichung (T') ist der Prozessausgang. Die Gleichung beschreibt daher, wie sich sowohl Änderungen am Heizer als auch Einlassstörungen auf die Auslass- oder Tanktemperatur auswirken.

Verwenden Sie die physikalische Form der Wärmeleistung, wenn angemessen

Wenn die Stellgröße die Heizleistung anstelle einer Heizmediumtemperatur ist, definieren Sie

[ Q'=Q-Q_s ]

Das Abweichungsmodell mit konstanten Stoffeigenschaften wird zu

[ \rho c_p V\frac{dT'}{dt}

\rho c_p F(T_i'-T')+Q' ]

oder

[ \boxed{ \frac{dT'}{dt} + \frac{F}{V}T'

\frac{F}{V}T_i' + \frac{Q'}{\rho c_p V} } ]

Diese Form ist für Berechnungen in Pilotanlagen oft transparenter, da (Q') die direkte physikalische Bedeutung einer Änderung der Heizleistung hat.

Erkennen Sie die Zustandsvariable

Für einen Heizer mit konstantem Volumen ist (T') die primäre thermische Zustandsvariable. Wenn sich der Flüssigkeitsstand oder das Volumen signifikant ändert, muss auch die Massenbilanz einbezogen werden:

[ A\frac{dh}{dt}=F_i-F ]

und die Energiebilanz kann geschrieben werden als

[ Ah\frac{dT}{dt}

F_i(T_i-T)+\frac{Q}{\rho c_p} ]

In diesem Fall sind die Temperaturdynamik und die Füllstandsdynamik gekoppelt, und eine einfache Temperaturgleichung erster Ordnung ist möglicherweise nicht mehr ausreichend.

Übertragungsfunktionen für die dynamische Analyse erhalten

Wenden Sie die Laplace-Transformation an

Die Laplace-Transformation von

[ \frac{dT'}{dt}+aT'

\frac{1}{\tau}T_i'+KT_{st}' ]

ergibt

[ s\bar T'(s)-T'(0)+a\bar T'(s)

\frac{1}{\tau}\bar T_i'(s) + K\bar T_{st}'(s) ]

Bei einer anfänglichen Temperaturabweichung von Null,

[ T'(0)=0 ]

gilt also

[ (s+a)\bar T'(s)

\frac{1}{\tau}\bar T_i'(s) + K\bar T_{st}'(s) ]

Leiten Sie die Übertragungsfunktion der Einlasstemperatur ab

Das Setzen der Heizerabweichung auf Null ergibt

[ \boxed{ \frac{\bar T'(s)}{\bar T_i'(s)}

\frac{1/\tau}{s+a} } ]

Diese Übertragungsfunktion beschreibt die Reaktion der Tanktemperatur auf Störungen der Einlasstemperatur.

Leiten Sie die Übertragungsfunktion des Heizeingangs ab

Das Setzen der Einlasstemperaturabweichung auf Null ergibt

[ \boxed{ \frac{\bar T'(s)}{\bar T_{st}'(s)}

\frac{K}{s+a} } ]

Wenn die Heizleistung (Q') direkt verwendet wird, lautet die entsprechende Übertragungsfunktion

[ \boxed{ \frac{\bar T'(s)}{\bar Q'(s)}

\frac{1}{\rho c_p V} \frac{1}{s+F/V} } ]

Diese Beziehungen unterstützen die Analyse von Sprung-, Impuls- und Rampenantworten und bilden eine Grundlage für den Entwurf von Rückkopplungsreglern.

Erweitern Sie die Bilanz, wenn der Heizer ein Mantel oder eine Schlange ist

Beziehen Sie die treibende Kraft der Wärmeübertragung ein

Wenn Wärme über einen Mantel oder eine Schlange zugeführt wird, kann die Leistung modelliert werden als

[ Q=UA(T_{st}-T) ]

wobei (U) der Gesamtwärmedurchgangskoeffizient und (A) die Wärmeübertragungsfläche ist.

Die transiente Energiebilanz wird zu

[ \rho c_p V\frac{dT}{dt}

\rho c_p F(T_i-T) + UA(T_{st}-T) ]

Durch Umstellen erhält man

[ \frac{dT}{dt} + \left( \frac{F}{V}+\frac{UA}{\rho c_p V} \right)T

\frac{F}{V}T_i + \frac{UA}{\rho c_p V}T_{st} ]

In dieser Darstellung gilt also

[ a=\frac{F}{V}+\frac{UA}{\rho c_p V} ]

und

[ K=\frac{UA}{\rho c_p V} ]

Die Abweichungsgleichung behält die gleiche Struktur bei:

[ \frac{dT'}{dt}+aT'

\frac{F}{V}T_i' + \frac{UA}{\rho c_p V}T_{st}' ]

Berücksichtigen Sie bei Bedarf Reaktionseffekte

Für einen reagierenden Rührkessel muss die Energiebilanz auch die durch Reaktionen erzeugte oder verbrauchte Wärme enthalten. Eine allgemeine Enthalpiebilanz kann dargestellt werden als

[ V\frac{d}{dt} \left(\sum_i C_i\bar H_i\right)

F\left[ \sum_i(C_i\bar H_i)_{\text{in}}

\sum_i C_i\bar H_i \right] +\dot Q ]

Reaktionswärme kann die Bilanz nichtlinear machen, da Reaktionsgeschwindigkeiten und Wärmetönungen von Temperatur und Zusammensetzung abhängen können.

Ein lineares Abweichungsmodell wird dann durch Linearisierung der nichtlinearen Bilanz um den gewählten stationären Zustand erhalten. Die Koeffizienten in der Abweichungsgleichung repräsentieren lokale Empfindlichkeiten an diesem Arbeitspunkt.

Verständnis der Kompromisse

Abweichungsmodelle sind lokale Modelle

Die Abweichungsgleichung ist normalerweise in der Nähe des gewählten stationären Zustands genau. Große Temperaturänderungen, sich ändernde physikalische Eigenschaften, variabler Durchfluss, Phasenänderungen oder stark temperaturabhängige Reaktionsgeschwindigkeiten können die Näherung mit konstanten Koeffizienten ungültig machen.

Ein neuer Arbeitspunkt erfordert möglicherweise eine neue Linearisierung und neue Modellparameter.

Der stationäre Zustand muss konsistent sein

Die Aufhebung der konstanten Terme funktioniert nur, wenn (T_s), (T_{i,s}) und (T_{st,s}) die tatsächliche stationäre Energiebilanz erfüllen. Wenn die Basiswerte inkonsistent sind, verbleiben Restkonstanten in der Abweichungsgleichung.

Dies ist eine häufige Quelle für scheinbare Modellierungsfehler.

Eine Annahme von konstantem Volumen kann Füllstandseffekte verbergen

Wenn Einlass- und Auslassdurchflussraten unterschiedlich sind, ändert sich der Tankfüllstand. Da das Flüssigkeitsvolumen die thermische Kapazität bestimmt, verändert die Füllstandsdynamik die Temperaturantwort.

Für füllstandsgeregelte Pilotanlagen ist das Modell mit konstantem Volumen oft angemessen. Während Füllstandstransienten sollten jedoch sowohl Massen- als auch Energiebilanzen beibehalten werden.

Heizertemperatur und Heizleistung sind keine identischen Eingänge

Die Verwendung von (T_{st}') setzt voraus, dass die Temperatur des Heizmediums die gewählte Eingangsgröße ist und die Wärmeübertragungsbeziehung in (K) integriert wurde. Die Verwendung von (Q') behandelt die Wärmeleistung als die direkte Stellgröße.

Die richtige Wahl hängt vom physischen Aktor und den in der Pilotanlage verfügbaren Messungen ab.

Anwendung auf Ihr Projekt

Die geeignete Formulierung hängt davon ab, ob die Pilotanlage für eine grundlegende Reaktionsanalyse, einen Regelungsentwurf oder eine detaillierte Anlagensimulation modelliert wird.

Wenn Ihr Hauptfokus auf der Temperaturreaktion auf Einlassstörungen liegt: Verwenden Sie (\frac{dT'}{dt}+aT'=\frac{1}{\tau}T_i'+KT_{st}'), wobei (T_i') als Störgröße behandelt wird.
Wenn Ihr Hauptfokus auf der Heizerregelung liegt: Verwenden Sie die Heizleistungsabweichung (Q') und schreiben Sie (\frac{dT'}{dt}+\frac{F}{V}T'=\frac{F}{V}T_i'+\frac{Q'}{\rho c_pV}).
Wenn Ihr Hauptfokus auf dem Design von Mantel oder Schlange liegt: Behalten Sie (UA) explizit bei, damit das Modell zeigt, wie Wärmeübertragungsfläche und Koeffizient sowohl Verstärkung als auch Ansprechzeit beeinflussen.
Wenn Ihr Hauptfokus auf der gemeinsamen Regelung von Füllstand und Temperatur liegt: Formulieren Sie gekoppelte Abweichungs-Massen- und Energiebilanzen mit (h') und (T') als Zustandsvariablen.
Wenn Ihr Hauptfokus auf einer reagierenden Pilotanlage liegt: Beginnen Sie mit der vollständigen Enthalpiebilanz und linearisieren Sie Reaktions-, Wärmeübertragungs- und Stoffeigenschaftsterme um den Arbeitspunkt.

Ein korrekt formuliertes Abweichungsmodell verwandelt das Arbeitspunktproblem des Rührkesselheizers in eine klare, handhabbare Beschreibung, wie Störungen und Regelungsmaßnahmen die Temperatur verändern.

Zusammenfassungstabelle:

Hauptaspekt	Zusammenfassung
Hauptannahme	Konstantes Volumen, Durchflussrate, Dichte und Wärmekapazität in einem ideal durchmischten Tank
Dynamisches Modell	`dT/dt + aT = (1/tau)Ti + KTst`
Abweichungsvariablen	`T' = T - Ts`, `Ti' = Ti - Ti,s` und `Tst' = Tst - Tst,s`
Abweichungsgleichung	`dT'/dt + aT' = (1/tau)Ti' + KTst'`
Häufige Störgröße	Einlasstemperaturabweichung, `Ti'`
Häufige Stellgröße	Heizmediumtemperatur, `Tst'`, oder Heizleistungsabweichung, `Q'`
Mantel- oder Schlangenmodell	Beinhaltet `UA` zur Darstellung von Wärmeübertragungseffekten
Modellbeschränkungen	Große Betriebsänderungen, variable Eigenschaften, Füllstandsänderungen und Reaktionseffekte können ein nichtlineares oder gekoppeltes Modell erfordern

Sind Sie bereit, die Ausbildung, Forschung oder Prozessentwicklung im Bereich der Temperaturregelung zu verbessern? Kontaktieren Sie LABPARK, um eine auf Ihre Bedürfnisse zugeschnittene Ausbildungs- oder Berufsbildungs-Pilotanlage für Verfahrenstechnik zu besprechen. LABPARK bietet praktische Pilotanlagensysteme für Chemieingenieurwesen, Bioprozess- und Biotechnologie sowie Anwendungen in der Umwelt- und Wasseraufbereitung. Universitäten, Forschungsinstitute und Unternehmen können unsere Anlagen für praktische Ausbildung, Prozessanalyse, Entwicklung von Regelsystemen, Scale-up-Studien und technische Schulungen nutzen, mit Konfigurationen, die auf Ihre Prozessziele und Laboranforderungen abgestimmt sind.