simpleFoamの計算結果を初期値にした場合のrhoPimpleFoamでの実行エラー

[MURAI Kazuhiko]

村井です。　いつもお世話になります。

昨日、rhoPimpleFoamでのメッシュ不具合での実行エラーについて、

解決して頂いたばかりでの、追加質問で恐縮ですが、お願いします。

rhoPimpleFoamで流れ場無しで計算をしようとしましたが、

流れの発達までに時間がかかるため、simpleFoamで計算した結果（U,P,phi）を

初期条件にしようと思い、

simpleFoamで計算した結果（U,P,phi）をrhoPimpleFoamにコピーして計算開始したところ、

今度は下記のエラーが出てきました。

--> FOAM FATAL ERROR: 

Different dimensions for =

     dimensions : [1 -1 -1 0 0 0 0] = [0 2 -1 0 0 0 0]

    From function dimensionSet::operator=(const dimensionSet&) const

    in file dimensionSet/dimensionSet.C at line 165.

FOAM aborting

次元が違うというのは分かりますが、dimensionsの左辺はどうやらalphaSgsとmuSgsのようです。

速度と圧力を持ってきただけでエラーが出るのは、よくわかりません。

もちろん、境界条件はrhoPimpleFoam用に書き換えております。

問題のケースファイルは、下記のURLにアップロードしました。

https://copy.com/xwE2E1LNzu0A

どなたか、ご存知の方がいらっしゃいましたら、アドバイスをお願いします。

よろしくお願いします。

[nozaki]

非圧縮性用のソルバー(simpleFoam) と圧縮性用のソルバー(rhoPimpleFoam)

とでは、圧力p 及び流束phi の次元がそれぞれ密度の分だけ異なるので、このような

次元が合わないというエラーが出ているのだと思います。

p, phi の dimensions をそれぞれ

[1 -1 -2 0 0 0 0]、[1 0 -1 0 0 0 0]

に変更すればこのエラーは回避できると思います

（本来は密度の分だけこれらの変数の internalField も修正するべきですが）。

以上です。

野崎

P.S. simpleFoam の代わりに rhoSimpleFoam を使用すればこのような次元の

書き換えが不要になると思います。

2013年12月29日日曜日 17時19分33秒 UTC+9 MURAI Kazuhiko:

[MURAI Kazuhiko]

nozakiさん

村井です。　ご回答ありがとうございます。

頂いたアドバイスの通り、次元を修正したら、無事に計算が流れました。

しかし…圧縮性の計算は、なかなか収束性が厳しいですね。

商用ソルバーで経験はしていますが、なかなか大変です。

ありがとうございます。

2013年12月29日 23:08 nozaki <fumiya...@gmail.com>:

--
このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
このグループから退会し、メールの受信を停止するには、openfoam+u...@googlegroups.com にメールを送信します。
このグループに投稿するには、open...@googlegroups.com にメールを送信してください。
http://groups.google.com/group/openfoam からこのグループにアクセスしてください。
その他のオプションについては、https://groups.google.com/groups/opt_out にアクセスしてください。

--
☆。.:*:・'°★。.:*:・'°☆。.:*:・'°★。.:*:・'°☆。.:*:・'°
Name          村井  和彦  （MURAI Kazuhiko）
Phone         0563-75-0320
Mobile         080-5103-9878
Zip-code      445-0073
Address       愛知県西尾市寄住町灯籠下36番地1
E-mail          murai.k...@gmail.com
★。.:*:・'°☆。.:*:・'°★。.:*:・'°☆。.:*:・'°★。.:*:・'°