離散化スキームの変更

[OKUYAMA]

chtMultiRegionSimpleFoamを使って円管内を流体が流れるモデルの定常解析をしています。

調べたいのは円管の出入口での温度と圧力の変化です。

これまでは対流項の離散化スキームを、以下の通り1次精度風上差分(upwind)にしていたのですが、レイノルズ数が約25万と大きめのため、高次精度のスキームに変えようと考えています。

divSchemes
{
    default          none;

    div(phi,U)       bounded Gauss upwind;

    div(phi,K)       bounded Gauss upwind;
    div(phi,h)       bounded Gauss upwind;
    div(phi,k)       bounded Gauss upwind;
    div(phi,omega)   bounded Gauss upwind;
    div(phi,epsilon) bounded Gauss upwind;
    div(phi,R)       bounded Gauss upwind;
    div(R)           bounded Gauss linear;
    div((muEff*dev2(T(grad(U))))) Gauss linear;
}

この場合、どの変数のスキームを変えるのがいいのでしょうか。

パターンとしては、

・とにかく全部(div(phi,U)～div((muEff*dev2(T(grad(U))))) 変える。

・対流項(div(phi,*))だけ変える。

・物理量に直接影響しそうな div(phi,U) とdiv(phi,h) だけ変える。

くらいかなと考えています。

初歩的な質問で申し訳ございませんが、よろしくお願いいたします。