$FOAM_TUTORIALS/lagrangian/MPPICFoam/injectionChannel
2つの流入口 lowerInlet、upperInlet からそれぞれ流速 (18.7939, 6.8404, 0) m/s、(18.7939, -6.8404, 0) m/s で空気が流入し、流路内で衝突します。衝突後の流れは流出口 outlet から計算領域外へ流出します。全体で0.1秒間分の計算を行ないます。
離散粒子の物性はディレクトリ constant 内のファイル kinematicCloudProperties で定義されています。粒子の流入位置、流入開始時刻、流入持続時間、速度、サイズ分布などは subModels の injectionModels で以下のようにそれぞれ model1、model2 として定義されています。
injectionModels
{
model1
{
type patchInjection;
parcelBasisType fixed;
patchName lowerInlet;
U0 (18.7939 6.8404 0);
nParticle 1;
parcelsPerSecond 1390885;
sizeDistribution
{
type normal;
normalDistribution
{
expectation 650e-6;
variance 25e-6;
minValue 500e-6;
maxValue 800e-6;
}
}
flowRateProfile constant 1;
massTotal 0;
SOI 0;
duration 60;
}
model2
{
type patchInjection;
parcelBasisType fixed;
patchName upperInlet;
U0 (18.7939 -6.8404 0);
nParticle 1;
parcelsPerSecond 1390885;
sizeDistribution
{
type normal;
normalDistribution
{
expectation 650e-6;
variance 25e-6;
minValue 500e-6;
maxValue 800e-6;
}
}
flowRateProfile constant 1;
massTotal 0;
SOI 0;
duration 60;
}
}
メッシュは以下の通りです。メッシュ数は9776です。
メッシュ
計算結果は以下の通りです。
時刻0.001秒での離散粒子
時刻0.005秒での離散粒子
時刻0.01秒での離散粒子
時刻0.02秒での離散粒子
時刻0.1秒での離散粒子
また最終時刻での流れ場の様子は以下の通りです。
時刻0.1秒での流速(U)
時刻0.1秒での圧力(p)
13分27.9秒 ※シングル、Inter(R) Core(TM) i7-2600 CPU @ 3.40GHz 3.40GHz