―1―
向流移動層内の粒子挙動と温度分布に及ぼす炉下部排出構造の影響

向流移動層内の粒子挙動と温度分布に及ぼす
炉下部排出構造の影響
夏 井 貴 臣 *・ 一 田 守 政Inﬂuence of furnace lower exhaust structure on particle behavior
and thermal distribution in counter-current moving bed
Takaomi NATSUI* and Morimasa ICHIDA（平成25年10月25日受理）

The influence of the furnace lower exhaust structure on the particle behavior and the thermal
distribution in the counter-current moving bed where the shaft reduction furnace was assumed for wasexamined with two dimensional reduced scale model，CCD camera，and infrared camera and the followingﬁndings were obtained. A peripheral exhaust and an increase in the outlet number reduce the stagnation
layer，expand the descent area (expansion of the gas circulation area)，and make the thermal distribution
uniform consequentially.
Keywords: counter-current moving bed，exhaust structure，particle behavior，PIV，thermal distribution

1.

緒言は，平均的な還元炉の 1/25 縮尺の二次元模型 ( 長さ325mm ×幅 40mm ×高さ 510mm) である。

向流移動層は高炉などの製鉄還元炉や下水汚泥溶
融炉，排煙処理システム等に広く用いられている1)。
還元炉プロセスは固体粒子と還元ガスの向流移動層であり，粒子降下挙動とそれに関連したガス流れによる温度分布の制御が安定稼動に不可欠とされている。そこで本研究では，還元炉プロセス内の粒子挙動と温度分布に及ぼす炉下部排出構造の影響を，シャフト還元炉の二次元縮尺模型，CCD カメラ，熱画像測定装置を用いて検討した。

2.

実験装置

2.1. 実験装置の構成
実 験 装 置 の 構 成 を 図 1 に 示 す。 実 験 装 置 は 還 元
炉 縮 尺 二 次 元 模 型， 温 風 発 生 装 置， ブ ロ ワ ー， 流量 計， 熱 画像 測 定 装 置， ノ ー トパ ソ コ ン に よ っ て構成される。本実験で使用される還元炉縮尺模型* 秋田高専専攻科修了生

秋田高専研究紀要第 49 号

図 1 実験装置の構成

2.2. 還元炉縮尺二次元模型
図 2 に還元炉縮尺二次元模型を示す。模型全面をアクリル板製にすることで装入試料の可視化を可能とした。装入試料として実際の還元炉で充填されるペレットを模擬したφ 3mm のアルミナ球を用い，高温ガスの代わりに温風を用いた。アルミナ球は装置上部の漏―2―夏井貴臣・一田守政斗から手動で装入する。送風羽口からは約 20℃の温風が吹き込まれる。送風羽口近傍に背面から温度計を差し込むことにより，羽口先端の温度を測定する。
また実験装置前面のアクリル板表面の中心温度を熱画像測定装置で測定する。実験装置に装入されたア
ルミナ球はホッパーを介し，直下に設置されたターンテーブルにより装置外へ排出される。ターンテーブルの回転数によりアルミナ球排出量を制御する。

2.3.

排出構造条件

図2(a) 全面排出(b) 周辺排出(c) 中心排出(d) 中間排出(e) 周辺中心排出(f ) 周辺中間排出還元炉縮尺二次元模型概観

排出構造条件を図 3 に示す。障害物を配置しない(a) 全面排出，左右に排出口を設けた (b) 周辺排出，中央に１本の排出口を設けた (c) 中心排出，中央に 2本の排出口を設けた (d) 中間排出，中心排出と中間排出にそれぞれ周辺排出を付加した (e) 周辺中心排出，(f ) 周辺中間排出の 6 水準で実験を行った。(b)~(f) は排出口の断面積を一定とし，また (e)，(f ) は炉壁近傍の排出口下部を塞がないように障害物の一角をホッパーの角度に合わせて加工した。

3.

基準実験条件の導出 2)還元炉内の装入物の降下挙動は重力とガスの慣性力に支配される。したがって支配的な物理法則はフ

図3

排出構造条件

表1

基準実験条件

ルード数 ( 以下 Fr 数とする ) であるので，実炉と模型で Fr 数が一致するように求めた 3) 基準実験条件を表１に示す。
これを参考にして今回の実験では送風量0.180Nm³/min， ア ル ミ ナ 球 排 出 量 は 1.62kg/minと設定した。アルミナ球排出速度を基準実験条件(0.303kg/min) と比較して約 5 倍にしたのは，ターンテーブル排出装置の排出速度下限値以下の値であり制御できないためであるが，この程度のアルミナ球排出速度の差異は許容範囲と考えた。

平成 26 年２月

―3―
向流移動層内の粒子挙動と温度分布に及ぼす炉下部排出構造の影響

4.

実験方法
本実験では粒子挙動測定と温度分布測定という 2種類の測定実験を行う。以下，4.1 に粒子挙動の測定
方法，4.2 に温度分布の測定方法を示す。

4.1. 粒子挙動の測定方法
2.1 で説明した様に，本実験装置は温風を送風する機能を備えているが，粒子挙動測定実験では無送風の条件下で実験を行った。粒子挙動測定は CCD カメラを用いる。二次元縮尺模型上部までアルミナ球を充填後，アルミナ球排出を開始する。排出開始後 10秒時，10 分時，20 分時，30 分時に CCD カメラでそれぞれ 10 秒間動画を撮影する。撮影した動画を 1 秒間に 30 枚の静止画像としてノート PC に保存する。
適切な時間間隔 ( 今回は 2 秒 ) の 2 枚の静止画像から，PIV 法により粒子速度解析を実施する。実験は前述の
図 4 に示す 6 水準で行った。
4.1.1 PIV 法 4)
PIV 法は，PC に取り込んだ流体の可視化画像データを利用して輝度分布によるパターン追跡からベクトル解析を行う方法である。図 4 のように，ある時間 to とその時間からΔ t たった画像があるとき，toに任意の点 P をとり，P の周辺にマトリクスをとる。
点 P のマトリクス内の輝度分布と類似しているものをΔ t 後の画像からとり点 Q とし，数学的に注目しているマトリクスと相互相関係数が最も高いものを類似しているとする。そしてこのＰＱを移動量として速度を求める。
上記からわかるように PIV 法の特徴である静止画像二枚から解析可能である点，粒子の高分布密度時に適している点，粒子以外（気流等）の画像に対しても有効である点から本実験では PIV 法を解析に用いた。
今回の実験条件 ( アルミナ球排出速度 1.62kg/min ，各種排出構造条件 ) では，PIV 法による粒子の速度ベクトルの解析条件として，解析範囲全域に速度ベクトルが表示されている静止画像のインターバル 60(2 秒 )の条件を採用した。
4.2. 温度分布の測定方法
二次元縮尺模型上部までアルミナ球を充填し，温風発生装置により昇温された空気をブロワーから送風羽口を介して実験装置本体へ送風し，実験装置およびアルミナ球を昇温する。この際の温風温度は 20℃
秋田高専研究紀要第 49 号

(1)t= to(2)t= to ＋Δ t
図4

PIV 法による輝度分布によるパターン追跡前後になる｡ この状態でアクリル板表面の中心温度( 図 3）が 17℃になるまで昇温する。
実験装置アクリル板表面の中心温度が約 17℃に到達したら，熱画像測定装置で温度分布を測定し，これを固定層の温度分布とした。その後アルミナ球を排出開始し，向流移動層温間実験を開始する。
アルミナ球排出後，5 分後ごとに温度分布を測定して，これを 30 分間行う。この温度分布を向流移動層の温度分布とした。実験は図 4 に示す 6 水準の炉下部排出構造条件で行った。

5.

実験結果

5.1. 粒子挙動の時間的変化の比較
全面排出 ( 図 3 の (a)) における 10 秒時，30 秒時，5 分時，10 分時，20 分時，30 分時の粒子速度分布を図 5 に示す。10 秒時，30 秒時には降下速度分布が径方向でばらつきがある。5 分時，10 分時には降下速度分布が高さ方向でばらつきがある。このため排出開始から 10 分時までは非定常状態であると判断できる。20 分時，30 分時には降下速度分布が径方向，高さ方向において共に一様となっているため，定常状態であると判断できる。よって今後の各実験条件における解析結果の比較には 30 分時の定常状態のデータを採用する。本二次元縮尺模型内の粒子は 1回の入れ替わりに約 6 分かかり，アルミナ球が 5~6回入れ替わることに相当している。

―4―夏井貴臣・一田守政(a)10 秒時(b)30 秒時(a) 粒子速度ベクトル(b) 粒子降下速度分布

図 6 全面排出における粒子降下挙動(c)5 分時(d)10 分時(a) 粒子速度ベクトル
図7(e)20 分時
図5(f)30 分時

全面排出時の粒子速度分布の時間変化

5.2. 粒子降下挙動に及ぼす排出構造の影響
5.2.1 全面排出
全面排出時 ( 図 3 の (a)) の粒子速度ベクトル解析結果と粒子降下速度分布を図 6 に示す。粒子の降下速度が径方向にほぼ均一である（図 6 の (b)）。また，速度ベクトルのほとんどが下向きである。炉壁近傍を除く粒子があたかも固定層の塊のように降下している（図 6 の (a)）。
5.2.2 周辺排出周辺排出時 ( 図 3 の (b)) の粒子速度ベクトル解析
結果と粒子降下速度分布を図 7 に示す。粒子は断面積の縮小した周辺の排出口に向かって速度分布を持ちながら降下している（図 7 の (b)）。また障害物直上に非常に大きな停滞層を形成している。

(b) 粒子降下速度分布

周辺排出における粒子降下挙動

5.2.3 中心排出
中心排出時 ( 図 3 の (c)) の粒子速度ベクトル解析結果と粒子降下速度分布を図 8 に示す。粒子は断面積の縮小した中心部の排出口に向かって速度分布を持ちながら降下している（図 8 の (b)）
。また，周辺部 ( 左右の炉壁近傍 ) に大きな停滞層を形成している（図 8の (a)）
。

(a) 粒子速度ベクトル

図8(b) 粒子降下速度分布

中心排出における粒子降下挙動
平成 26 年２月

―5―
向流移動層内の粒子挙動と温度分布に及ぼす炉下部排出構造の影響

5.2.4 中間排出
中間排出時 ( 図 3 の (d)) の粒子速度ベクトル解析結果と粒子降下速度分布を図 9 に示す。粒子は断面積の縮小した中間部の排出口に向かって速度分布を持ちながら降下している（図 9 の ( ｂ )）。また，中心に小さな停滞層，周辺に大きな停滞層を形成している（図 9 の (a)）
。

が降下している為，周辺の停滞層が大幅に縮小している。また，粒子の流れが周辺部に分散されるため，中心部の停滞層も若干縮小している。

(a) 中心排出(b) 周辺中心排出(1) 粒子速度ベクトル(a) 粒子速度ベクトル
図95.3.

(b) 粒子降下速度分布

中間排出における粒子降下挙動

停滞層の縮小対策(c) 中心排出(d) 周辺中心排出5.2 より周辺排出，中心排出，中間排出のいずれの場(2) 粒子降下速度分布

合も中心部または周辺部に大きい停滞層を形成するこ

図 10 中心排出，周辺中心排出時の粒子降下挙動とが示された。周辺部に大きい停滞層が形成される中心排出時・中間排出時には周辺排出を付加し，周辺部停滞層の縮小効果を確認する実験を行った。そして中心部に大きい停滞層が形成される周辺排出時には中心排出・中間排出を付加し，中心部停滞層の縮小効果を確認する実験を行った。

5.3.1

周辺排出付加効果の比較１中心排出（図 3 の (c)）と中心排出に周辺排出を付加した周辺中心排出（図 3 の (e)）の粒子速度ベクトル解
析結果と粒子降下速度分布を図 10 に示す。周辺部に大きい停滞層が形成されている中心排出時（図 10 の (a)）(a) 中間排出(b) 周辺中間排出(1) 粒子速度ベクトルに比べて，周辺中心排出時（図 10 の (b)）は炉壁近傍で粒子が降下しているため，周辺部の停滞層が大幅に縮小している。

5.3.2

周辺排出付加効果の比較 2

中間排出（図 3 の (d)）と中間排出に周辺排出を付加した周辺中間排出（図 3 の (f)）の粒子速度ベクトル解
析結果と粒子降下速度分布を図 11 に示す。周辺部に大(c) 中間排出(d) 周辺中間排出

きい停滞層が形成されている中間排出（図 11 の (a)）(2) 粒子降下速度分布に比べ周辺中間排出（図 11 の (b)）は炉壁近傍で粒子

図 11 中間排出，周辺中間排出時の粒子降下挙動

秋田高専研究紀要第 49 号

―6―
夏井貴臣・一田守政

5.3.3 中心・中間排出付加効果
周辺排出 ( 図 3 の (b)) 周辺排出（図 3 の (b)）に中心排出を付加した周辺中心排出（図 3 の (e)），中間排出を付加した周辺中間排出（図 3 の (f)）の粒子速度ベクトル解析結果と粒子降下速度分布結果を図 12に示す。排出口本数の増加に伴い，停滞層が縮小している。これは個々の障害物の幅が縮小しているためと推定される。

(a) 周辺排出(b) 周辺中心排出(a) 全面排出時(b) 周辺排出時(c) 中心排出時(d) 中間排出時(c) 周辺中間排出(1) 粒子速度ベクトル

図 13(d) 周辺排出(e) 周辺中心排出

全面排出時，周辺排出時，中心排出時，中間排出時の温度分布(f) 周辺中間排出(2) 粒子降下速度分布
図 12 周辺排出，周辺中心排出，周辺中間排出の粒子降下挙動

は還元炉内壁の付着物への進展が予想され望ましくない。そこで温度分布を均一化させることを目的として，周辺部に排出口を付加した周辺中心排出時 ( 図3 の (e))，周辺中間排出時 ( 図 3 の (f)) の温間実験を

5.4. 温度分布に及ぼす排出構造の影響
5.4.1 温度分布に及ぼす排出構造の影響
全面排出時（図 3 の (a)），周辺排出時（図 3 の (b)），中心排出時（図 3 の (c)）
，中間排出時（図 3 の (d)）の温度分布を図 13 に示す。全面排出時には周辺部の温度は高く中心部の温度は低い ( 図 13 の (a))。周辺排出時には大きな停滞層を形成している中心部の温度が高い ( 図 13 の (b))。中心排出時（図 13 の (c)）
，中間排出時（図 13 の (d)）ともに大きな停滞層が形成されている周辺部の温度が高く中心部の温度が低い，不均一な温度分布である ( 図 13 の (c) と (d))。

行い，それぞれ中心排出時 ( 図 3 の (c))，
中間排出時 ( 図3 の (d)) と温度分布を比較した。
中心排出時 ( 図 3 の (c)) と周辺中心排出時 ( 図 3 の(e)) の温度分布測定結果を図 14 に示す。中心排出時に周辺排出を付加した周辺中心排出時には，周辺部( 停滞層 ) の高温領域が縮小した結果，温度分布は均一化される ( 図 14 の (b))。
中間排出時 ( 図 3 の (d)) と周辺中間排出時 ( 図 3 の(f )) の温度分布測定結果を図 15 に示す。図 14 と同様に中間排出時に周辺排出を付加した周辺中間排出時には周辺部の高温領域が縮小し，中心部温度が上昇した結果，温度分布が均一化される ( 図 15 の (b))。

5.5. 温度分布の均一化対策
5.5.1 中心排出時及び中間排出時に周辺排出を付加した場合の温度分布
図 13 の (c)，(d) で見られたような周辺部の高温化
平成 26 年２月

―7―
向流移動層内の粒子挙動と温度分布に及ぼす炉下部排出構造の影響

5.6.

各実験条件による温度分布結果差異のメカニ
ズム
5.6.1 渦度から見た排出構造と粒子挙動の関係
5.2，5.3 で評価した各排出構造条件時の粒子の動きを，粒子流線の渦度解析により検討した。
一般的に渦度ベクトルΩは流速ベクトル V=(Vx，Vy，Vz) により (1) 式のように表される。
(a) 中心排出時
図 14(b) 周辺中心排出時(1)

中心排出時と周辺中心排出時の温度分布

本実験は2次元であるので(2)式になる。
(2)各排出構造条件における渦度分布解析結果を図 17に示す。周辺排出時 ( 図 3 の (b)) には羽口レベルから羽口上 150mm の範囲で周辺部の粒子流線の渦度の絶対値が増加し，停滞層を形成している中心部の(a) 中間排出時
図 15(b) 周辺中間排出時中間排出時と周辺中間排出時の温度分布

渦度は 0 である ( 図 17 の (b))。中心排出時 ( 図 4 の(c)) には羽口レベルから羽口上 150mm の範囲で中心部の粒子流線の渦度の絶対値が増加し，停滞層を形成している周辺部の渦度は 0 である ( 図 17 の (c))。
中間排出時 ( 図 3 の (d)) には羽口レベルから羽口上5.5.2

周辺排出時における排出口本数を増加した
場合の温度分布
周辺排出時 ( 図 3 の (b))，周辺中心排出時 ( 図 3 の(e))，周辺中間排出時 ( 図 3 の (f )) の温度分布測定結果を図 16 に示す。 周辺排出時には大きな停滞層が形成されている中心部が高温領域となっているが ( 図16 の (a))，周辺中心排出時，周辺中間排出時には高温領域は縮小している ( 図 16 の (b) と (c))。排出口本数が増加するにしたがって，高温領域の停滞層が縮小し，温度分布が均一化される。

150mm の範囲で中間部の粒子流線の渦度の絶対値が増加し，停滞層を形成している周辺部の渦度は 0 である ( 図 17 の (d))。渦度の絶対値の増加は，粒子が高さ方向だけでなく，断面積の小さい排出口へ向かって径方向に移動しながら降下している。
排出口を増加させた周辺中心排出時 ( 図 3 の (e)) と周辺中間排出時 ( 図 3 の (f )) には，粒子流線が複雑になり高さ方向と径方向の全体にわたり渦度の絶対値が増加している ( 図 17 の (e) と (f))。
一方，全面排出時には渦度の絶対値が径方向においてきわめて小さい ( 図 17 の (a))。これは，炉壁近傍を除いた粒子があたかも固定層のように降下しており，プラグフローに近いことを示している。

5.6.2(a) 排出口 2 本

図 16(b) 排出口 3 本(c) 排出口 4 本

周辺排出時と周辺中心排出時と周辺中間排出時の温度分布

秋田高専研究紀要第 49 号

周辺への排出口付加 ( 排出口本数増加 ) による温度分布均一化のメカニズム
周辺部に排出口を付加 ( 排出口本数増加 ) したことによる温度分布均一化のメカニズムを，中心排出と周辺中心排出を例として図 18 に示す。中心排出時には 5.6.1 の渦度解析 ( 図 17 の (c)) で説明したように左右の障害物上の周辺部に非常に大きな停滞層が形成される ( 図 18 の (a))。一方中心排出に周辺排出を―8―
夏井貴臣・一田守政(a) 全面排出(b) 周辺排出(a) 中心排出時
図 18(c) 中心排出(b) 周辺中心排出時

周辺に排出口を付加 ( 排出口本数増加 ) による温度分布均一化のメカニズム①周辺排出付加により炉壁近傍で粒子が降下し，(d) 中間排出周辺の停滞層が大幅に縮小する。
②排出口本数増加により個々の障害物の幅が縮小し，停滞層が縮小する。
⑵

温度分布は粒子の降下挙動と密接に関連し，停⑶周辺排出付加と排出口本数増加は，停滞層を縮滞層は高温領域となる。
小，かつ降下領域を拡大 ( ガス流通領域の拡大 )(e) 周辺中心排出(f) 周辺中間排出

図 17 各実験条件における渦度分布

付加した周辺中心排出時には，5.6.1 の渦度解析 ( 図17 の (e)) で説明したように，左右の障害物上の周辺部に形成される停滞層が大幅に縮小される ( 図 18 の(b))。
周辺部の排出口付加は，総排出口断面積が同じであっても，停滞層縮小と粒子降下領域の拡大に効果的である。実用化を想定した三次元のシャフト還元炉 ( 円形 ) では，周辺部の面積比率が大きいため，上記の周辺排出口付加効果はより大きいと考えられる。停滞層の縮小と粒子降下領域の拡大により，高温の停滞層が縮小され，拡大した降下領域のガス流させるため，温度分布を均一化する。
⑷

渦度解析により粒子の動きを定量的に評価できることを確認した。

7．参考文献1)

篠原邦夫，高橋洋志，中村正秋：移動層工学，北海道大学図書刊行会，札幌，(2000)，11.

2)

江守一郎：模型実験の理論と応用，技報堂出版，東京，(1985)，85.

3)

金田啓佑，一田守政：秋田工業高等専門学校，

4)

加賀昭和，井上義雄，山口克人：可視化情報，

47(2012)，1.
14(1994)，53.

通量が増加することから，結果的に装置内の温度分布が均一化されると考えられる。

6.

結言
シャフト還元炉を想定した向流移動層における粒子挙動と温度分布に及ぼす炉下部排出構造の影響を二次元縮尺模型，CCD カメラ，熱画像測定装置を用いて検討し，以下の知見を得た。
⑴

周辺排出付加，排出口本数増加は周辺部，中心部の停滞層縮小対策として有効である。
平成 26 年２月