openfoamquest

OpenFOAMや流体解析の使い方を紹介しています

OpenFOAM 多孔質体 設定(4) porous model

 

 前回の続きです。

 

openfoamquest.hatenablog.jp

 

 

例えば円錐でcylindricalの軸をZ軸にして30°の多孔質モデルの場合

半径(板厚)方向のみに流れると考えた場合、dとfの係数は1番目

以外を係数1000倍とする。

1番目が半径、2番目が接線、3番目が軸方向の係数となる。

(X軸をcylindricalの軸にする場合は、dとfの係数は

1番目が軸、2番目が半径、3番目が接線方向となる)

 

porosity1
{
    type            explicitPorositySource;

    explicitPorositySourceCoeffs
    {
        selectionMode   cellZone;

        cellZone        fff-fluid_porous;

        type            DarcyForchheimer;

        // D 100;  // Very little blockage
        // D 200;  // Some blockage but steady flow
        // D 500;  // Slight waviness in the far wake
        D 10E9; // Fully shedding behavior 1000=alpha
        F 1000;   // beta

        d   ($D 10E12 10E12);
        f   ($F 1000000 1000000);


        coordinateSystem
        {
            type cylindrical;
            origin  (0 0 1);
            note "Any note can be written here.";
            rotation
            {
                type axisAngle;
                axis      (0 0 1);
                angle     30;
        }
        }
    }
}

参考

http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2008/HaukurElvarHafsteinsson/haukurReport.pdf