大気解析の境界条件

2019年2月13日

はじめに

大気解析の境界条件について。

バージョン

OpenFOAM v1812

大気解析の境界条件

U

    minX
    {
        type            atmBoundaryLayerInletVelocity;
        flowDir         constant (1 0 0);
        zDir            constant (0 0 1);
        Uref            constant 2;
        Zref            constant 10;
        z0              0.1;
        zGround         0;
    }

大気境界層の流束分布を与える。flowDir は流れの方向、zDir は高さ方向のベクトル、Uref は Zref 高さでの流束、z0 は粗度、zGround は地面の高さである。粗度は長さの単位を持ち、地面の状態によって値が変わる。障害物の高さが高いほど値は大きい傾向にはあるが、高さそのものを表しているわけではない。

k

    minX
    {
        type            atmBoundaryLayerInletK;
        flowDir         constant (1 0 0);
        zDir            constant (0 0 1);
        Uref            constant 2;
        Zref            constant 10;
        z0              0.1;
        zGround         0;
    }

epsilon

    minX
    {
        type            atmBoundaryLayerInletEpsilon;
        flowDir         constant (1 0 0);
        zDir            constant (0 0 1);
        Uref            constant 2;
        Zref            constant 10;
        z0              0.1;
        zGround         0;
    }

速度に応じた k, epsilon を設定する。

nut

    ground
    {
        type            nutkAtmRoughWallFunction;
        z0              uniform 0.1;
        value           $internalField;
    }

z0 に応じた粗さの壁関数を用いる。