気体物性値の出力方法に関して

Kazuki Sakaki

unread,

Dec 23, 2015, 1:33:39 AM12/23/15

to OpenFOAM

この度はお世話になります。

現在、OpenFoamを用いた燃焼解析(reactingFoamやsprayFoam)を行っておりますが、
気体の定圧比熱等の出力がうまくできずに困っています。

計算自体はCHEMKIN formatのファイルを用いて問題なく行えているのですが、
圧力、温度、化学種濃度のように時刻毎の定圧比熱、定積比熱、比熱比の3つを出力しようとしたのですがうまく行きません。

createFields.Hにおいて密度にならって

volScalarField Cp
    (
        IOobject
        (
            "Cp",
            runTime.timeName(),
            mesh,
            IOobject::NO_READ,
            IOobject::AUTO_WRITE
        ),
        thermo.Cp()
    );

という記述を追加してみたのですが、意図したとおりには動きませんでした。
（コンパイルはとおり、計算自体は問題なく実行できました）
同様の問題をCFD onlineにおいても発見したのですが、いまいち解決策が分からず途方にくれています。

本件に関しまして知見をお持ちの方がいらっしゃいましたら、お時間のある際にご教示いただければ幸いです。

以上、よろしくお願いいたします。

Hiroaki Fujio

unread,

Dec 23, 2015, 3:46:38 AM12/23/15

to open...@googlegroups.com

先に"Cp"の名前でIOobjectにNO_WRITEで登録されている変数があるため出力されません。

下記heThermo.Cの333行目をご覧ください。

https://github.com/OpenFOAM/OpenFOAM-2.3.x/blob/master/src/thermophysicalModels/basic/heThermo/heThermo.C

そこで、例えば下記のようにするとCp_ファイルとして出力されます。

volScalarField Cp_

(

IOobject

(

"Cp_",

runTime.timeName(),

mesh,

IOobject::NO_READ,

IOobject::AUTO_WRITE

),

thermo.Cp()

);

また、このままではCp_の中身が更新されませんので、毎タイムステップのエネルギーと圧力の式を解いた後に

Cp_ = thermo.Cp();

を行う必要があります。

それから特に大した問題ではありませんが、NO_READよりREAD_IF_PRESENTの方が良いかと思います。

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

Kazuki Sakaki

unread,

Dec 23, 2015, 8:56:30 PM12/23/15

to OpenFOAM

早速のご対応ありがとうございました。

また丁寧に補足のご説明をいただきありがとうございます。

初歩的なトラブルで申し訳ありませんでしたが、お陰様で無事解決いたしました。