研究実績の概要 |
流れ場の最小長さスケールと同程度の径をもつ粒子が希薄に混入した場合の乱れの増減は,基本的かつ重要な乱流変調問題として知られており,数値シミュレーションによる解明が期待されている. よく用いられている質点モデルは,粒子からの直接的な影響を受ける流れ場の空間スケールや粒子表面の応力分布が考慮されないため,乱流変調の基礎研究には不十分である. 粒子境界層を解像できる細かい計算格子を用いると,過大な計算負荷のため,小さな領域しか解析できなくなる. したがって広範なスケール幅をもつ乱流に対する粒子混入の効果を解析するためには,粒子径と同程度の幅の計算格子を用いて粒子と流体の相互作用を扱うモデルが求められる. そこで本研究では,粒子サイズと同程度の領域内での体積平均に基づいて, 粒子表面応力を反映させた流体相への作用力モデルの提案を行った.まず流れの基礎方程式に体積平均を施した方程式を導出し, 流体相への作用力の表式を示した.境界適合格子を用いた数値解析結果に基づいて流体相への作用力のモデリングを行った. まず, 一様流および相対速度の無い剪断流の中における粒子表面応力分布について, 遅い流れを仮定した解析解との比較を行い, 定性的な一致を確認した. この結果より作用力の関数形を決定し, 定量的な不一致を係数調整で補うことでモデリングを行った. 従来の質点モデルと異なり, 本モデルでは流体場へ与える作用力が粒子との距離だけでなく, 粒子から見た向きにも依存するという非等方性を導入できる. このことはせん断の効果を作用力モデルに取り入れられるという意味でも本質的な重要性をもつ. いくつかの剪断流に対しては, 一様流および相対速度のないせん断流に対するそれぞれのモデルの重ね合わせが, 作用力の再現に有効であることを数値解析の結果から確認した.
|