離散化スキームのエラー

[ogura]

はじめましてoguraと申します。

自身が作成したモデルをcfMeshを用いてメッシュにし、実行したのですが以下のようなエラーが表示されました。
エラー内容が離散化スキーム内のものであることは確認できるのですが、どのように直したら良いのか分かりません。よろしければ、アドバイスお願いいたします。


以下は表示されたエラーです。
--> FOAM FATAL IO ERROR:
keyword div((nuEff*dev2(T(grad(U))))) is undefined in dictionary "/home/custom/Desktop/aaa/system/fvSchemes.divSchemes"

file: /home/custom/Desktop/aaa/system/fvSchemes.divSchemes from line 30 to line 37.

    From function const Foam::entry& Foam::dictionary::lookupEntry(const Foam::word&, bool, bool) const
    in file db/dictionary/dictionary.C at line 441.

FOAM exiting


また、該当箇所かと思われる離散スキームです。

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default         none;
    div(phi,U)      Gauss limitedLinearV 1;
    div(phi,k)      Gauss limitedLinear 1;
    div(phi,epsilon) Gauss limitedLinear 1;
    div(phi,R)      Gauss limitedLinear 1;
    div(R)          Gauss linear;
    div(phi,nuTilda) Gauss limitedLinear 1;
    div((nuEff*dev(T(grad(U))))) Gauss linear;
}


OpenFOAMについて初心者でして、分からないところだらけです。ぜひともお助けください。
よろしくお願いいたします。

[Youhei Takagi]

ogura様

高木と申します。エラーにあるとおり、キーワード"div((nuEff*dev2(T(grad(U)))))"に
対する設定が記述されていません。divSchems{ }の中に一行追加して、例えば、

divSchemes
{
    default         none;
    div(phi,U)      Gauss limitedLinearV 1;
    div(phi,k)      Gauss limitedLinear 1;
    div(phi,epsilon) Gauss limitedLinear 1;
    div(phi,R)      Gauss limitedLinear 1;
    div(R)          Gauss linear;
    div(phi,nuTilda) Gauss limitedLinear 1;
    div((nuEff*dev(T(grad(U))))) Gauss linear;

    div((nuEff*dev2(T(grad(U))))) Gauss linear;
}

などにすれば良いかと思います。ただし、指定したスキームが正しい/計算が
正常に行えるかどうかは別問題なので、ご自身でご検討ください。

2017年5月30日 18:28 ogura <djt...@gmail.com>:

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

[ogura]

高木様

お早い返信ありがとうございます。oguraです。
ご指摘の通りに、divScheme{}を変更しましたところ、うまく計算が行われました。
丁寧な説明ありがとうございました。また、指定したスキームも正しい計算が行われたと思われます。
本当にありがとうございました。

失礼します。