冷却を制御する幾何学的相似条件

「スケールアップ理論を考えてみよう ー 冷却編【冷却モデル】」のページでは，冷却速度に関して2つの考え方を紹介しました。

最初に，スケールアップ前後で冷却速度を一定にすることを目指してみたいと思います。

単位体積あたりの槽壁伝熱量

”単位体積あたりの冷却速度”を”単位体積あたりの槽壁伝熱量q_V”と呼びたいと思います。

単位体積あたりの槽壁伝熱量q_Vは，「単位時間あたりの伝熱量q/製品仕込量V」で表すことができます。

まずは，式の形から考えてみたいと思います。

「スケールアップ理論を考えてみよう ー 冷却編【冷却モデル】」のページで紹介したように，製品中の熱は”乳化物 → 伝熱境膜 → 乳化槽 → 冷却水”の順に移動すると考えることにしました。

ここでは一纏めにした伝熱係数を考えることにして，これを総括伝熱係数Uと定義することにします。

📝[memo]　乳化撹拌装置における乳化槽を使用するときは，総括伝熱係数Uは一定であるとして考えます。

単位時間あたりの槽壁伝熱量q_Vは，総括伝熱係数U，製品と冷却水の温度差ΔT，乳化槽の伝熱面積A/製品仕込量Vに比例します。

冷却条件を等しくしてスケールアップするので，製品と冷却水の温度差ΔTは一定になります。

すると，単位時間あたりの伝熱量qは乳化槽の「伝熱面積A/製品仕込量V」のみに比例することがわかります。

そして，「伝熱面積A/製品仕込量V」は乳化槽の形状に起因する条件であるため，単位時間あたりの伝熱量は”力学・運動学的相似条件”を適用すべきかもしれませんが，結果的に”幾何学的相似条件”として考えたいと思います。

例えば，10倍にスケールアップ（製品仕込量Vが10倍）した事例を考えます。

このとき，乳化槽の伝熱面積Aは10倍にはならず，約4.64倍（= 10^2/3倍）にしかなりません。

すなわち，「伝熱面積A/製品仕込量V」を一定にすることはできません。

その結果，スケールアップ前後で単位体積あたりの槽壁伝熱量q_Vに0.464倍（= 10^-1/3倍）の差異が生じていることがわかります。

📝[memo]　q_V2 = q₂/V₂ = 10^2/3q₁/10V₁ = 10^-1/3q₁/V₁ = 10^-1/3q_V1

このように，10倍にスケールアップしたときは単位体積あたりの槽壁伝熱量q_Vに0.464倍の差異が生じていました。

では，このような差異が生じないためにはどうしたら良いかというと，q_V1 = q_V2を満たせば良いことがわかります。

そこで，x倍にスケールアップしたとして計算をしてみると，「q_V2 = x^-1/3q_V1」の関係式が得られます。

これを満たすためには，x = 1（1倍にスケールアップ＝同じ乳化槽）のときになります。

したがって，幾何学的相似を満たすことで，冷却速度（単位体積あたりの槽壁伝熱量q_V）を一定にすることは困難であるという結論が得られます。

ジャケット冷却とコイル冷却

ジャケット冷却では，乳化槽の伝熱面積Aや製品仕込量Vのバランスが崩れるのでスケールアップが上手くできませんでした。

そこで，乳化槽の伝熱面積Aに相当する部分を大きくすることで，スケールアップが上手くできるように工夫することができます。

その手段の一つとして，コイル冷却があります。

一方で，「スケールアップ理論を考えてみよう ー 冷却編【冷却モデル】」のページで紹介したように，エマルション製品は冷却するにつれて増粘していく傾向が見られます。

すると，乳化槽壁面やコイル表面に伝熱境膜が生成します。

乳化槽の場合は掻取ミキサーを使用してこのような伝熱境膜を掻き取ることができますが，コイル冷却の場合は伝熱境膜は残存したままです。

その結果，コイル冷却では冷却効率が悪化していくことが懸念されます。

したがって，エマルション製品の冷却工程ではコイル冷却を採用することはほとんどありません。

製品と冷却水の温度差を変数とするとき

ところで，製品中の熱は”乳化物 → 伝熱境膜 → 乳化槽 → 冷却水”の順に移動するので，総括伝熱係数Uは定数とみなされます。（総括伝熱係数Uは変化しません。）

一方，製品と冷却水の温度差ΔTは冷却条件によって変化します。

したがって，製品と冷却水の温度差ΔTを変化させることで，約0.464倍（= 10^-1/3倍）に変化する「伝熱面積A/製品仕込量V」を補正できるかもしれません。

これは，冷却速度を一定にするためには，あくまでもq_Vが等しくなれば良いという考え方です。

結論を言うと，製品と冷却水の温度差ΔTが約2.15倍（= 10^1/3倍）にになれば良いことがわかります。

それでは，実務においてこのような対応ができるでしょうか？

例えば，80℃の製品を冷却することを考えます。

冷却水の温度が10℃のときは，製品と冷却水の温度差ΔTは70℃になります。

10倍にスケールアップしたときは製品と冷却水の温度差ΔTが150℃であることが求められるので，-70℃の冷却水を使用しなければなりません。

このような対応は実際には無理な話なので、製品と冷却水の温度差ΔTを変化させるという考え方は採用できません。

冷却を制御するための計算（推算）式

そこで，力学・運動学的相似条件を考えることになります。

このとき登場するのが”冷却を制御するための計算（推算）式”です。

これは，撹拌によって製品を流動させ，“槽壁伝熱量”を変化させるものです。

回転数Nが決まれば，単位体積あたりの槽壁伝熱量q_Vが算出できそうです。

問題点

しかしながら，この式を使用するにあたって問題点があります。

冷却時において，温度によってエマルション製品の密度ρや粘度ηが連続的に変化します。

そのため，密度ρや粘度ηが定数として決まらないため，単位体積あたりの槽壁伝熱量q_Vが一定とは限りません。

したがって，単位体積あたりの槽壁伝熱量q_Vを一定とする（＝冷却速度を一定とする）条件を算出することは困難であるという結論が得られます。

単位体積あたりの槽壁伝熱量の式の導出

ご参考までに，単位体積あたりの槽壁伝熱量の式の導出したいと思います。

ひたすら式変形が続くので，読み飛ばしても問題ありません。

🚩 [引用：社団法人化学工学協会編『化学工学便覧』丸善，1988]

永田の式

hD/λ = α(C_pη/λ)^1/3(η/η_w)^0.14Re^β₁×(d/D)^β₂(ΣC_i/iz)^β₃(Σb_i/D)^β₄×n_p^β₅(sinθ)^β₆(Z/D)^β₇

⇔ h ∝ Re^β₁/D …(1)

インデックスの詳細は明記していませんが，スケールアップ前後で変化しない定数を無視すると上記の式が得られます。

伝熱面積

A = kD² …(2)

面積は長さの2乗に比例します。（kは比例定数）

掻取ミキサーの直径

d = D …(3)

掻取ミキサーの羽根径＝乳化槽径として，ホームページの表記に合わせます。

全伝熱量

q = h・ΔT・A …(4)

伝熱係数hを使用した単位時間あたりの伝熱量を示しています。

式(1), (2), (3)を式(4)に代入すると，

q = h・ΔT・A = (Re^β₁/D)・ΔT・kD² ∝ Re^β₁・ΔT・D …(5)

製品仕込量

V = sD³ …(6)

体積は長さの3乗に比例します。（sは比例定数）

単位体積あたりの槽壁伝熱量

q_V = q/V = Re^β₁・ΔT・D/sD³ ∝ Re^β₁・ΔT/D² = (D²Nρ/60η)^β₁・ΔT/D² = D^{(2β₁ – 2)}(Nρ/60η)^β₁・ΔT

式(5)÷式(6)より算出します。