OpenFOAM weirOverflowを自作STLとSnappyHexMeshでカスタマイズしてみた

まくまく

OpenFoamを使って流体解析をしよう！今回はinterFoamの使い方を学ぶため公式チュートリアルのweirOverflowを実行してみようと思います。

コンテンツ

1 weirOverflowをコピー
2 実行
3 自作STLで実行
4 処理の流れ
- 4.1 注意点：alpha.water
- 4.2 注意点：フォルダ削除

weirOverflowをコピー

cp -r $FOAM_TUTORIALS/multiphase/interFoam/RAS/weirOverflow

チュートリアルをコピーします。interFoamのRASフォルダの中になります。

実行

Allrunで実行するとこのような感じになります。左側赤色のinletから流入し、堰を超えて右側青色のoutletへ流れていきます。

アニメーションはこちら。

改訂新版 OpenFOAMの歩き方 (技術の泉シリーズ（NextPublishing）)

created by Rinker

自由表面流れについては、こちらの書籍の「第五章二層流」に詳細な記載があります。damBreakの基礎が学べると思います。

#include        "include/initialConditions"

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    inlet
    {
        type            variableHeightFlowRateInletVelocity;
        flowRate        $inletFlowRate;
        alpha           alpha.water;
        value           uniform (0 0 0);
    }

    outlet
    {
        type            zeroGradient;
    }

    lowerWall
    {
        type            noSlip;
    }

    atmosphere
    {
        type            pressureInletOutletVelocity;
        value           uniform (0 0 0);
    }

    defaultFaces
    {
        type            empty;
    }
}

Uファイルはこのようになっていて、variableHeightFlowRateInletVelocityで流入量を設定しています。実際の数値はinitialConditionsにあり75となっていました。

inletFlowRate        75;
pressure             0;
turbulentKE          4.14e-03;
turbulentEpsilon     4.39e-05;

試しに、流入量をvolumetricFlowRateで75m^3/sにして実行してみると見た目では大きな差はなさそうです。もっと水面が変化する条件だとvariableHeightFlowRateInletVelocityの効果を感じれるのかもしれませんね。このあたりは今後要検証ということにしておきましょう。