チュートリアルのrhoPimpleFoamおよびsonicFoamにおけるrhoについて

[Choko]

はじめまして

OpenFOAM初心者のChokoと申します。
表題の件についてお伺いしたく、投稿させていただきました。

本当に始めのはじめでつまずいていてお恥ずかしい限りなのですが、
どうかよろしくお願いいたします。

使用バージョンはOpenFOAM2.1.1です。

今後乱流圧縮性遷音速（非定常）気体を扱って研究をしていきたいことから、
チュートリアルのsonicFoamを走らせたところ、
endTimeまで回ったものの発散し、
logを確認するとrhoのresidualおよびiterationがずっと0のままとなっていました。

設定はデフォルトのまま何も変えておらず、ソルバのrhoEpnの確認もしたのですが、
原因がわかりません。
チュートリアルを走らせるに当たり何か行わなくてはならない設定があるのでしょうか。

sonicFoamがもともとバグを含んでいる可能性を考え（そのようなことはないのでしょうが）、
同じく圧縮性乱流rhoPimpleFoamでも試したところ、
可視化した結果としてはそれらしくなりましたが
logを確認するとsonicFoamと同じくrhoのresidualおよびiterationがずっと0のままでした。
（本質問とは直接関係のないことですが、非定常での研究をしていきたいと考えております。）

境界条件やメッシュ等を自分の研究対象のものに変更し、
境界条件として指定したTotalPressureおよびTotalTemparature内のrho等を変更してみるなどしたのですが、
rhoの件は解決しないままです。

どのようにすれば解決できるか、お力を貸していただけますでしょうか。

どうぞよろしくお願いいたします。

[ohbuchi]

こんにちは。

チュートリアルをそのまま動かして発散するとすれば、OpenFOAMがちゃんとセットアップ

できていないか、バグがあるかのどちらかです。

尚、私の環境ではsonicFoamチュートリアルは問題なく計算できます。

具体的にどのチュートリアルが発散したのでしょうか？

また、sonicFoamは圧縮性ソルバで、rhoは状態方程式から計算します。従って反復計算は

必要ありません。偏微分方程式で計算するのは、p,U,Tです。

2013年2月10日日曜日 20時13分35秒 UTC+9 Choko:

[Choko]

ohbuchi様

返信ありがとうございます。

発散したのはsonicFoam/ras/prismです。
試した限りでは他のチュートリアルは正常に動いているようです。

rhoの計算の件、
以前自分でも圧縮性のコードを書いて計算していたはずなのに、
久々の流体力学とバグへの戸惑いとで失念しておりました。

丁寧なご指摘ありがとうございます。

[ohbuchi]

sonicFoam/ras/prismに関するバグ修正履歴やバグレポートは出ておらず

最新版の当該ファイルは2.1.xリリース時から修正されていません。

何が原因か判りませんが、OpenFOAMが正しくセットアップされていない

可能性があります。

バイナリ版もリリースされていますので、そちらをお試しになるのも良いでしょう。

2013年2月11日月曜日 11時09分29秒 UTC+9 Choko:

[Choko]

ohbuchi様

prismの件、情報ありがとうございます。
バイナリ版を試してみようと思います。

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