攪拌槽式ヒーターの動的エネルギー収支は、過渡的な収支式から定常状態の収支式を差し引くことで記述されます。 定容積で完全混合されたタンクにおいて、流量と物性が一定である場合、温度モデルは以下のようになります。

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

ここで、(T')、(T_i')、(T_{st}') はそれぞれタンク温度、入口温度、加熱入力の公称定常値からの偏差です。この形式は、応答解析や制御設計において重要な動的変化のみを分離します。

要点： すべての変数を定常運転点からの偏差として定義し、それらの定義を動的エネルギー収支式に代入し、定常状態の式を用いて定数項を相殺します。得られた偏差モデルは、公称運転点のオフセットを含まずに過渡的な温度挙動を記述します。

過渡的なエネルギー収支から始める

制御体積と仮定の定義

完全混合攪拌槽式ヒーターでは、タンク内の液温は均一であると仮定されます。一般的なパイロットプラントモデルでは、液量（V）、体積流量（F）、密度（\rho）、比熱（c_p）が一定であると仮定します。

これらの仮定の下では、運動エネルギー、位置エネルギー、軸仕事、体積変化の影響は無視されます。支配的な項は、エネルギーの蓄積、入口および出口のエンタルピー流、および外部から供給される熱です。

次元付きエネルギー収支式の記述

熱入力（Q）がある非反応性液体タンクの場合、

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

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

ここで：

(T) はタンク内の液温、
(T_i) は入口温度、
(F) は体積流量、
(V) は液量、
(Q) はヒーターの熱負荷です。

タンクは完全に混合されているため、出口温度は (T) となります。

収支式を標準的な動的形式にする

(\rho c_p V) で割ると、

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

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

係数 (F/V) は滞留時間の逆数です。ここで

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

と定義すると、

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

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

ヒーター入力が変数 (T_{st}) で表される場合、熱入力項は (KT_{st}) と書くことができます。得られるモデルは一般的に以下のように表現されます。

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

ここで (a)、(1/\tau)、(K) は、流量、体積、比熱、熱伝達に関するパラメータをまとめたものです。

モデルを偏差変数に変換する

公称運転点の定義

定常状態の運転条件を

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

とします。定常状態ではタンク温度は変化しないため、

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

したがって、定常状態の式は

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

となります。この式は、動的モデルが展開される運転点を示しています。

定常状態からの偏差の定義

偏差変数は以下の通りです。

[ T'=T-T_s ]

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

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

同様に、

[ T=T_s+T' ]

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

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

(T_s) は定数であるため、

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

動的収支式への代入

以下の式から始めて、

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

代入すると、

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

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

整理すると、

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

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

定数項は定常状態の関係式

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

を満たしているため、相殺されます。偏差変数による式は

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

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

となります。これが偏差形式による目的の動的エネルギー収支式です。

偏差モデルの解釈

操作量と外乱の分離

加熱入力 (T_{st}')、あるいは同等に (Q') は、通常操作量の偏差です。入口温度の偏差 (T_i') は一般的に外乱として扱われます。

タンク温度の偏差 (T') はプロセス出力です。したがって、この式はヒーターの変化と入口の外乱の両方が、出口またはタンク温度にどのように影響するかを記述します。

適切な場合に物理的な熱負荷形式を使用する

操作入力が加熱媒体の温度ではなくヒーターの熱負荷である場合、

[ Q'=Q-Q_s ]

と定義します。定物性の偏差モデルは以下のようになります。

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

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

または

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

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

(Q') はヒーター負荷の変化という直接的な物理的意味を持つため、この形式はパイロットプラントの計算においてより明確であることが多いです。

状態変数の認識

定容積ヒーターの場合、(T') が主要な熱状態変数です。液位や体積が大きく変化する場合は、質量収支も考慮する必要があります。

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

また、エネルギー収支は以下のように書くことができます。

[ Ah\frac{dT}{dt}

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

その場合、温度のダイナミクスは液位のダイナミクスと結合されるため、単一の一次温度式では不十分な場合があります。

動的解析のための伝達関数の取得

ラプラス変換の適用

以下の式をラプラス変換すると、

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

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

以下が得られます。

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

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

初期温度偏差がゼロの場合、

[ T'(0)=0 ]

したがって、

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

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

入口温度伝達関数の導出

ヒーター偏差をゼロに設定すると、

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

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

この伝達関数は、入口温度の外乱に対するタンク温度の応答を記述します。

ヒーター入力伝達関数の導出

入口温度偏差をゼロに設定すると、

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

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

ヒーター負荷 (Q') を直接使用する場合、対応する伝達関数は

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

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

これらの関係は、ステップ応答、パルス応答、ランプ応答の解析をサポートし、フィードバックコントローラー設計の基礎を提供します。

ヒーターがジャケットやコイルの場合の収支式の拡張

熱伝達の駆動力を含める

ジャケットやコイルを通じて熱が供給される場合、熱負荷は以下のようにモデル化できます。

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

ここで (U) は総括伝熱係数、(A) は伝熱面積です。

過渡的なエネルギー収支式は

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

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

整理すると、

[ \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} ]

したがって、この表現では

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

および

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

となります。偏差方程式は同じ構造を維持します。

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

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

必要に応じた反応効果の考慮

反応を伴う攪拌槽の場合、エネルギー収支には反応によって生成または消費される熱を含める必要があります。一般的なエンタルピー収支は以下のように表せます。

[ 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 ]

反応熱は、反応速度や反応熱が温度や組成に依存する場合があるため、収支式を非線形にする可能性があります。

その場合、選択した定常状態の周囲で非線形収支式を線形化することで、線形偏差モデルが得られます。偏差方程式の係数は、その運転点における局所的な感度を表します。

トレードオフの理解

偏差モデルは局所的なモデルである

偏差方程式は、通常、選択した定常状態の近くで正確です。大きな温度変化、物理的性質の変化、流量の変動、相変化、または温度に強く依存する反応速度は、定数係数の近似を無効にする可能性があります。

新しい運転点では、新しい線形化と新しいモデルパラメータが必要になる場合があります。

定常状態は一貫していなければならない

定数項の相殺は、(T_s)、(T_{i,s})、(T_{st,s}) が実際の定常状態のエネルギー収支を満たしている場合にのみ機能します。基準値に一貫性がない場合、偏差方程式に残留定数項が残ります。

これは、明らかなモデリングエラーの一般的な原因です。

定容積の仮定は液位の影響を隠す可能性がある

入口流量と出口流量が異なる場合、タンクの液位が変化します。液体体積は熱容量を決定するため、液位のダイナミクスは温度応答を変化させます。

液位制御されたパイロットプラントでは、定容積モデルが適切な場合が多いです。ただし、液位の過渡状態では、質量収支とエネルギー収支の両方を保持する必要があります。

ヒーター温度とヒーター負荷は同一の入力ではない

(T_{st}') を使用することは、加熱媒体の温度が選択された入力であり、熱伝達関係が (K) に組み込まれていることを前提としています。(Q') を使用することは、熱負荷を直接の操作量として扱います。

正しい選択は、物理的なアクチュエータとパイロットプラントで利用可能な測定値によって異なります。

プロジェクトへの適用方法

適切な定式化は、パイロットプラントが基本的な応答解析、制御設計、または詳細な機器シミュレーションのいずれのためにモデル化されているかによって異なります。

主な焦点が入口外乱に対する温度応答である場合： (\frac{dT'}{dt}+aT'=\frac{1}{\tau}T_i'+KT_{st}') を使用し、(T_i') を外乱入力として扱います。
主な焦点がヒーター制御である場合： ヒーター負荷の偏差 (Q') を使用し、(\frac{dT'}{dt}+\frac{F}{V}T'=\frac{F}{V}T_i'+\frac{Q'}{\rho c_pV}) と記述します。
主な焦点がジャケットやコイルの設計である場合： (UA) を明示的に保持し、伝熱面積と係数がゲインと応答時間の両方にどのように影響するかをモデルに示します。
主な焦点が液位と温度の同時制御である場合： (h') と (T') を状態変数として、結合された質量およびエネルギー収支の偏差モデルを定式化します。
主な焦点が反応を伴うパイロットプラントである場合： 完全なエンタルピー収支から始めて、運転点の周囲で反応、熱伝達、物性の各項を線形化します。

正しく定式化された偏差モデルは、攪拌槽式ヒーターの運転点の問題を、外乱や制御アクションがどのように温度を変化させるかという明確で扱いやすい記述へと変換します。

要約表：

主要な側面	要約
主な仮定	完全混合槽における定容積、流量、密度、比熱
動的モデル	`dT/dt + aT = (1/tau)Ti + KTst`
偏差変数	`T' = T - Ts`, `Ti' = Ti - Ti,s`, `Tst' = Tst - Tst,s`
偏差方程式	`dT'/dt + aT' = (1/tau)Ti' + KTst'`
一般的な外乱	入口温度の偏差、`Ti'`
一般的な操作量	加熱媒体温度 `Tst'` またはヒーター負荷の偏差 `Q'`
ジャケット/コイルモデル	熱伝達効果を表すために `UA` を含める
モデルの制限	大きな運転変化、物性変化、液位変化、反応効果には非線形または結合モデルが必要な場合がある

温度制御のトレーニング、研究、またはプロセス開発を改善する準備はできていますか？お客様のニーズに合わせてカスタマイズされた教育用または職業訓練用の単位操作パイロットプラントについて、LABPARKまでお問い合わせください。LABPARKは、化学工学、バイオプロセスおよびバイオテクノロジー、環境および水処理アプリケーション向けの、実践的なパイロットプラントシステムを提供しています。大学、研究機関、企業は、当社の機器を使用して、実践的な教育、プロセス解析、制御システム開発、スケールアップ研究、技術トレーニングを行うことができ、構成はお客様のプロセス目標と実験室の要件に合わせて調整されます。