fvSolution内の収束判定値の設定について

[Daiki Kono]

大学の研究にて層流の定常解析をsimpleFoamを用いて行っています。

素人質問で申し訳ないのですが、fvSolution内の収束判定条件について質問させていただきます。

収束判定値の設定はsolversとSIMPLEにおいてそれぞれ設定しなければならないようですが、

どういった違いがあるのでしょうか。

またrelaxationFactorsはどういったものなのでしょうか。

ざっくりとした質問で申し訳ありませんが、ご教授いただけたら幸いです。

参考にfvSolutionの記述を以下に掲載します。

*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    p
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-06; //絶対残差許容値
        relTol          0.01; //相対残差許容値
    }

    U
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-05;
        relTol          0.01;
    }
}

SIMPLE
{
    nNonOrthogonalCorrectors 0;

    residualControl
    {
 p               1e-5;
 U               1e-5;
    }
}

relaxationFactors
{
    fields
    {
        p               0.3;
    }
    equations
    {
        U               0.7;
    }
}

// ************************************************************************* //

[nozaki]

はじめまして。

simpleFoam では、

運動方程式と連続の式を離散化してできる２つの連立方程式 (UEqn と pEqn) を

交互に解き、修正を繰り返しながら、両方を同時に満たす流束場と圧力場を求めます。

離散化した方程式の左辺と右辺の差を残差 (residual) と言いますが、

この残差の大きさで収束判定を行います。

各変数について残差の大きさが、"residualControl" で指定した値より

小さくなると、計算は終了します。

運動方程式と連続の式を離散化すると、UEqn と pEqn の２つの連立方程式ができますが、

OpenFOAM では、この連立方程式の解を反復法で求めます。

この反復計算の繰り返しをコントロールするのが、"solvers" の設定です。

"relaxationFactors" は、緩和係数と訳されて、計算の安定性や収束性をコントロールします。

私の資料（http://www.slideshare.net/fumiyanozaki96/openfoam-32087641）の

6章にも説明があるので、読んでみて下さい。

以上ご参考になれば。

野崎

2015年1月29日木曜日 18時15分52秒 UTC+9 Daiki Kono:

[Daiki Kono]

野崎様

返信が遅くなり、申し訳ありません。

アドバイスありがとうございます。

もう一度じっくり基礎から学ぼうと思います。

2015年1月30日金曜日 22時58分31秒 UTC+9 nozaki:

[nozaki]

もし理解しにくい点などがあれば、ぜひ教えてください。

資料をより良くするために、皆様からのご指摘・アドバイスをお待ちしてます。

2015年2月12日木曜日 16時49分05秒 UTC+9 Daiki Kono:

[Cicko LoL]

野崎様の詳細な資料によって
fvSolution内SIMPLEのU,PのtoleranceがSIMPLE定式のUの運動方程式とPのポアソン方程式の反復計算に使われるということはわかりました
しかしfvSolution内solversのU,Pのtoleranceはどこに使われるのでしょうか？