| 研究概要 |
離散要素法(DEM)による流動層内流動解析がかなり行われるようになり,実際の流動をよく再現できることが示されている.しかし,この手法では個々の粒子運動を計算するため,計算機能力の制約から,Geldartの流動様式分類でA粒子やC粒子に分類されるような小さな粒子に対しては適用が困難である.そこで本課題ではDEMによる流動解析が小さな粒子の流動予測にも適用できるのかを調べ,さらに相似則モデルを開発することにより,大きな粒子を用いた計算で小さな粒子の流動を予測することを最終目的としている.
本年度はまず,A粒子からなる非常に小規模な2次元流動層に対して,DEMによる数値シミュレーションを行った.A粒子の流動様式の最大の特徴は,粒子層が一様に膨張する均一流動化状態が存在することであるが,本研究室現有の既存の計算コードでもA粒子に対して均一流動化状態が再現されることを確認した.しかし比較のために実験を行ったところ,圧力損失特性において計算と実験との間で定量的な差異が見られた.この差異の原因として,潤滑効果や付着力の影響を考えた.そこで既存の計算コードで考慮している重力,接触力,流体抵抗力に加えて,潤滑効果や付着力も考慮できるように計算コードを改良した.その結果,付着力の影響を考慮した場合に計算と実験の圧力損失特性は定量的に一致した.ただ,付着力と潤滑効果の両方を考慮した場合には計算が不安定になるというモデルの問題点も明らかになっており,モデルの改良が必要である.
また,付着性の強い微小粒子を流動化させる脈動流動層と呼ばれる装置に対してDEMによる数値シミュレーションを行った.ただし,気流への脈動付加が粒子流動に及ぼす影響のみを調べるため,付着力などは考慮していない.その結果,適当な周波数で気流に脈動を加えた場合,粒子の対流が大きく促進されることが分かった.
|
