壁の設定について

612 views
Skip to first unread message

unread,
Apr 5, 2010, 5:58:34 AM4/5/10
to OpenFOAM
Openfoam初心者です。

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;
}
}

//
************************************************************************* //

E.Mogura

unread,
Apr 15, 2010, 4:27:38 PM4/15/10
to OpenFOAM
E.Mogura です。

controlDictやfvSchemeの設定など、他の条件がわからなかったので、
rhoTurbFoamのチュートリアルケース(cavity)の設定をそのまま流用し、
境界の名前を整合させ、圧力だけを貴条件に変更して計算してみました。
但し、メッシュがこのままでは計算時間がかかりすぎるので、
block分割数は(30x30x1)と減らしました。

⇒計算初期で、わずかに圧力反射は見られましたが、すぐに収まりました。

ただ気になった点として、
lInf(0.3)が系全体の長さ(0.1)に対して大きすぎると思うのですが、
何か意図があったんでしょうか?

入口圧が10100で、pInf(100000)なので、出口圧が100075にしかなってくれず、
これはこれで、waveTransmissiveの仕様通りの状態になってくれてはいたんですが。。。

ちなみに、このケースでは、waveTransmissiveを使わないで、
fixedValueでも、すんなり計算できることも確認しました。

以上、ご参考まで。

unread,
Apr 23, 2010, 12:03:55 AM4/23/10
to OpenFOAM
雑な質問に対して返信していただきありがとうございます。
また、返信が遅れてしまい申し訳ありませんでした。

>入口圧が10100で、pInf(100000)なので、出口圧が100075にしかなってくれず、
>これはこれで、waveTransmissiveの仕様通りの状態になってくれてはいたんですが。。。

これは、圧力が透過しなかったため
わかりやすいように入力したもので深い意味はありませんでした。

lInfについてよく調べ、もう一度解析を行いたいと思います。
また、fixedValueでも、試したいと思います。

ご助言ありがとうございました。

--
このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
このグループに投稿するには、open...@googlegroups.com にメールを送信してください。
このグループから退会するには、openfoam+u...@googlegroups.com にメールを送信してください。
詳細については、http://groups.google.com/group/openfoam?hl=ja からこのグループにアクセスしてください。

Reply all
Reply to author
Forward
0 new messages