流出境界条件について

[Katsunori Ikuma]

今年度から研究でOpenFOAMを使用することになりました大学院のIkumaと申します．
初心者なので初歩的な質問をするかもしれませんが，宜しくお願い致します．

現在，角柱を挿入したLES解析を行なっています．

本日は，流出境界条件について質問があります．
流出境界条件にadvectiveを用いようとしているのですが，
下記のURL（advectiveの概要）を見てもあまり理解できず困っています．
どなたか御教授頂けないでしょうか？

宜しくお願い致します．

http://foam.sourceforge.net/docs/cpp/a00029.html

[ohbuchi]

src/finiteVolume/fields/fvPatchFields/derived/advective/advectiveFvPatchField.H
の冒頭部分にあるDescriptionに詳しい説明があります。いわゆる対流流出境界条件で、
DDt(psi, U)＝０となる様にpsiの境界条件を設定します。
非定常解析でしか使えません（初期U=0だから）。
記述方法は下記の通りです。

outlet
{
    type     advective;
    phi       phi;
}

以上、ご参考まで。



2013年8月7日水曜日 18時28分30秒 UTC+9 Katsunori Ikuma:

[Katsunori Ikuma]

ohbuchi様

Ikumaと申します。
返信が大変遅くなり申し訳ございませんでした。

詳しいご説明ありがとうございました。
advectiveFvPatchField.Hを参照させていただいたところ、
psiとUにはどのような値が入るのか理解できませんでした。
初歩的な質問で大変申し訳ございませんが、ご教授頂けないでしょうか？
宜しくお願い致します。

2013年8月8日木曜日 12時58分28秒 UTC+9 ohbuchi:

[ohbuchi]

DDT(psi,U)=∂ψ/∂t+∂Uψ/∂x=０、実質微分ですね。psiは解くべきスカラー変数で、Uは速度です。

コード内では面積フラックスphiをフェース面積ベクトルmagSfで割って速度ベクトルUを求めています。

　Foam::advectiveFvPatchField<Type>::advectionSpeed() 

計算は陰的に行うので、マトリクスの係数行列を修正することで境界条件が課されます。

　void Foam::advectiveFvPatchField<Type>::updateCoeffs()

で係数を修正しています。これはすべての境界条件にあてはまります。

この様に定常解析では、境界条件のUが更新されないのでadvectiveが機能しません。

尚、fvSchemesでddtSchemeをsteadyStateにするとadvectionが指定しない様にWarningが表示されます。

以上、ご参考まで。

2013年8月23日金曜日 12時45分01秒 UTC+9 Katsunori Ikuma:

[Katsunori Ikuma]

ohbuchi様

Ikumaと申します。

詳しいご説明ありがとうございました。
理解を深めることができました。

流出境界条件にadvectiveを用いようと思います。

2013年8月23日金曜日 19時38分29秒 UTC+9 ohbuchi: