OpenFOAMの1.5の解析ソルバー「rhoTurbFoam」を使いRANSの解析を行いたいと思っているのですが、壁の条件で困っています。透
過壁にしたいのですが圧力が反射してかえってきてしまいます。以下は、テスト用のblockMeshDictとPです。おかしなところがあったら指摘し
ていただけないでしょうか。
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
* * * //
convertToMeters 0.1;
vertices
(
(0 0 0) //0
(1 0 0) //1
(1 1 0) //2
(0 1 0) //3
(0 0 0.1) //4
(1 0 0.1) //5
(1 1 0.1) //6
(0 1 0.1) //7
);
blocks
(
hex (0 1 2 3 4 5 6 7) (300 50 1) simpleGrading (1 1 1)
);
edges
(
);
patches
(
wall walls
(
(3 7 6 2)
(1 5 4 0)
)
patch inlet
(
(0 4 7 3)
)
patch outlet
(
(1 2 6 5)
)
empty frontAndBack
(
(0 3 2 1)
(4 5 6 7)
)
);
mergePatchPairs
(
);
//
************************************************************************* //
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
* * * //
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 100000;
boundaryField
{
walls
{
type zeroGradient;
}
inlet
{
type fixedValue;
// Field Value
value 100100;
}
outlet
{
type waveTransmissive;
field p;
phi phi;
rho rho;
psi psi;
gamma 1.3;
fieldInf 1e5;
lInf 0.3;
value uniform 1e5;
}
frontAndBack
{
type empty;
}
}
//
************************************************************************* //
controlDictやfvSchemeの設定など、他の条件がわからなかったので、
rhoTurbFoamのチュートリアルケース(cavity)の設定をそのまま流用し、
境界の名前を整合させ、圧力だけを貴条件に変更して計算してみました。
但し、メッシュがこのままでは計算時間がかかりすぎるので、
block分割数は(30x30x1)と減らしました。
⇒計算初期で、わずかに圧力反射は見られましたが、すぐに収まりました。
ただ気になった点として、
lInf(0.3)が系全体の長さ(0.1)に対して大きすぎると思うのですが、
何か意図があったんでしょうか?
入口圧が10100で、pInf(100000)なので、出口圧が100075にしかなってくれず、
これはこれで、waveTransmissiveの仕様通りの状態になってくれてはいたんですが。。。
ちなみに、このケースでは、waveTransmissiveを使わないで、
fixedValueでも、すんなり計算できることも確認しました。
以上、ご参考まで。