二次元計算の厚み方向の寸法について

[otisin]

はじめまして。

最近、OpenFOAMを使い始めたものです。

現在矩形領域内に円柱を設置し一方向から一様流量の流体を流すという初歩的なことをしております。

その後に、円柱を角速度ωで回転させ「pisoFoamによる２次元円柱周りの流れの解析」

http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?plugin=attach&refer=%C2%E8%A3%B2%A3%B4%B2%F3%CA%D9%B6%AF%B2%F1%A1%A7H250810&openfile=2Dflow_around_cylinder_with_pisoFoam_rev.pdf

を参考に寸法を設定したのですが、p.25に「２次元問題とするためz方向厚を1mmと小さくする」といった表記があったため、そのまま適用すると抗力係数cdの値に約1000倍の狂いが発生しました。

比較として「The turbulent mixing layer: geometry of large vortices」の論文を使用させていただいたのですが、流れ場は完全に一致したこと、各ディレクトリやファイルもすべて確認したのですが原因が分かりませんでした。

初歩的な質問で申し訳ございませんが、ご教授いただけると幸いです。

宜しくお願い致します。

参考になるかわかりませんが、各ファイルも載せておきます。

他に必要なファイルがあればお申し付けください。

0/p

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0;

boundaryField

{

    ROTOR

    {

type zeroGradient;

    }

    INLET

    {

        type            fixedValue;

value uniform 0;

    }

    OUTLET

    {

        type fixedValue;

        value uniform 0;

    }

    UPANDDOWN

    {

        type fixedValue;

        value uniform 0;

    }

    INTERFACE1

    {

        type cyclicAMI;

        value $internalField;

    }

    INTERFACE2

    {

        type cyclicAMI;

        value $internalField;

    }

}

0/U

    ROTOR

    {

        type            movingWallVelocity;

        value           uniform (0 0 0);

    }

    INLET

    {

        type fixedValue;

        value uniform (1 0 0);

    }

    OUTLET

    {

type inletOutlet;

inletValue uniform (1 0 0);

value $internalField;

    }

    UPANDDOWN

    {

type inletOutlet;

inletValue uniform (1 0 0);

        value $internalField;

    }

    INTERFACE1

    {

type cyclicAMI;

        value $internalField;

    }

    INTERFACE2

    {

type cyclicAMI;

        value $internalField;

    }

}

constant/dynamicMeshDict

dynamicFvMesh   dynamicMotionSolverFvMesh;

motionSolver    solidBody;

cellZone        FLUID1;

solidBodyMotionFunction  rotatingMotion;

origin        (0 0 0);

axis          (0 0 1);

omega         1.0; // rad/s

system/controlDict/functions

    forces

    {

    type            forceCoeffs;

    libs            ("libforces.so");

    writeControl    timeStep;

    timeInterval    0.01;

    log             true;

    patches         ("ROTOR");

    pName           p;

    UName           U;

    rho             rhoInf;      // Indicates incompressible

    rhoInf          1.0;           // Redundant for incompressible

    liftDir         (0 1 0);

    dragDir         (1 0 0);

    CofR            (0 0 0);  // Axle midpoint on ground

    pitchAxis       (0 0 1);

    magUInf         1.0;

    lRef            2.0;        // Wheelbase length

    Aref            2.0;        // Estimated

    }