残差のプロットについて

[k]

openfoam初心者のkです。

現在、foam-extend-3.1をインストールし
チャネルでの解析をしていますが、解析自体はうまくいくのですが
function objects を使った残差のプロットがうまく行きません。

現在行っている方法は、
http://www.geocities.jp/penguinitis2002/study/OpenFOAM/proc_results.html
を参照し

system/controlDict に下に示されているものを追加し、
計算後にexecFlowFunctionObjectsを実行していますが
うまく行きません。
(バージョン3.0x以前　ではない方を使用。)

functions
{
    residuals
    {
        type residuals;
        libs ("libutilityFunctionObjects.so");
//        functionObjectLibs ("libutilityFunctionObjects.so");
        fields (U p);
    }
}

現状、エラー内容も理解できていません。

--> FOAM FATAL ERROR:
Unknown function type residuals

Table of functionObjects is empty


    From function functionObject::New(const word& name, const Time&, const dictionary&)
    in file db/functionObjects/functionObject/functionObject.C at line 73.

FOAM exiting


また、別の方法として、同URL内にある logファイル保存し、
plotというスクリプトを用意して行う方法も試しましたがうまく行きません。

このスクリプトというのは新しいドキュメントという
解釈であっているでしょうか。

この際、下のように出てきます。
bash: ./plot: 許可がありません


初歩的な質問で申し訳ございません。
エラー内容及び解決策のご教授お願いいたします。

[すもも]

foam-extend は使っていないのでわかりませんが、
単純に residuals という functionObjects がない可能性があります。

plot のほう

>bash: ./plot: 許可がありません

は、実行権限がないだけと思われるので、

$ chmod +x plot

として実行許可を設定してから実行すればよいと思います。


2017年8月31日木曜日 15時57分29秒 UTC+9 k:

[k]

すもも様

openfoamにはあって、foam-extendにはないことも
あるのですね。


plot

言葉足らずで申し訳ありません。
URLにある

実行権限を設定する。

$ chmod +x plot

つぎのように実行する。

$ ./plot

という操作を行っていますが、うまく行きません。
スクリプトを用意したあと、この２つのみ行いましたが
間違っているでしょうか。 

よろしくお願いいたします。

[k]

すもも様

residualsというfunctionObjectsがないのではないか
という事に関してですが、

src/postProcessing内のディレクトリにresidualsがない場合は、
functionObjectsではできないということでしょうか。

今確認したところ、postProcessing内にresidualsというディレクトリはありませんでした。

[すもも]

>openfoamにはあって、foam-extendにはないことも

>あるのですね。

foam-extend が residuals が入る前の OpenFOAM をベースにしていれば、

ないかもしれません (わざわざ追加していなければ)。

plot については、gnuplot が使われていますが、

gnuplot 自体は実行できる状態でしょうか？

[k]

すもも様

なるほど．ありがとうございます．

gnuplotは問題なく実効することが出来るので，

実効権限を与える方法ありますでしょうか．

また，スクリプトというのは，

gedit plot

と打ち込み作成したのですが，間違えでしょうか．

初心者で初歩的な質問になってしまい申し訳ございません．

よろしくお願いいたします．

[すもも]

plot の中身が gnuplot のコマンドになっているなら、代わりに

$ gnuplot -p plot

でも実行できるはずです。

あるいは、logs の中身を gnuplot で手動で表示するとか、
LibreOffice Calc でグラフ表示してもよいと思います。

[oda]

すもも様

gnuplot -p plot を実行することで残差のプロットを行うことが出来ました．

初歩的な質問にも関わらずご教授頂きありがとうございます．

また，別件で質問なのですが，

現状では，k-εモデルを使って計算を回していましたが，

LRRモデルに切り替えて乱流モデル以外変えずに

計算を行ったところ上手くいきません．

blockMeshは出来ますが，simpleFoamを実行したところ

Segmentation fault (コアダンプ)

と出てしまいました．

過去の投稿を見て緩和係数を小さくしたりしましたが

上手くいきません．

何か注意点などありますでしょうか．

oda

[すもも]

>また，別件で質問なのですが，

タイトルを変えて別件として再投稿したほうがよいと思います。
どなたかにコメントもらえるかもしれません。

Segmentation falut 直前のログデータがあると、解決のヒントになるかもしれません。

[nakagawa]

oda様

すでに解決済みのようですが，foam-extendでのresidualグラフ化について見かけた情報を共有させてください。

foam-extendを使っていないため動作未確認ですが，logファイルからpylab利用してグラフを作成するツールが同梱されているようです。

ご参考まで。

foam-extend-foam-extend-3.2/bin/plotResidual.py

https://github.com/Unofficial-Extend-Project-Mirror/foam-extend-foam-extend-3.2/blob/master/bin/plotResidual.py

なかがわ

2017年9月5日火曜日 14時31分57秒 UTC+9 oda:

[oda]

すもも様

初歩的な質問にも関わらず最後まで丁寧にありがとうございます．

もう少し考えた上で投稿したいと思います．

oda

[oda]

nakagawa様

pylabを初めてお聞きしたので，よく調べ動作確認してみたいと思います．

ありがとうございます．

oda