トップ - 技術ドキュメント - OpenFOAM 付属チュートリアル一覧

タンク内の水のスロッシング(3次元、6自由度)

更新:2017/7/28
OpenFOAM 4.x

ケース

$FOAM_TUTORIALS/multiphase/interDyMFoam/ras/sloshingTank3D6DoF

概要

解析モデル 解析モデル

メッシュ移動機能を使用して揺れるタンクの中に入った水の動きの解析を40秒分、行います。水の体積はタンクの2割程度とし、残りは空気であるとします。

水、空気の各相の物性はファイル constant/transportProperties で指定しています。

メッシュの移動はファイル constant/dynamicMeshDict で以下の様に指定します。

dynamicFvMesh   solidBodyMotionFvMesh;
solidBodyMotionFvMeshCoeffs
{
    solidBodyMotionFunction tabulated6DoFMotion;
    tabulated6DoFMotionCoeffs
    {
        CofG            (0 0 0);
        timeDataFileName "$FOAM_CASE/constant/6DoF.dat";
    }
}

指定されたファイル constant/6DoF.dat には以下の様に40秒までの並進3方向、回転3方向の値が記載され、これによって運動が定義されます。

100
(
(0 ((0 0 0) (0 0 0)))
(0.40404 ((0.401298 0.952899 0.321827) (4.82741 2.79073 2.00649)))
(0.808081 ((0.786273 1.8071 0.635266) (9.52899 5.3597 3.93137)))

……(省略)……

(39.596 ((1.625 0.77368 -0.260102) (-3.90153 5.2878 8.12498)))
(40 ((1.82589 1.65428 -0.575807) (-8.6371 2.70906 9.12945)))
)

またファイル system/controlDict で以下の様に設定し、指定点、指定面での圧力 p を出力しています。出力データはフォルダ postProcessing 以下に保存されます。

functions
{
    probes
    {
        type            probes;
        libs ("libsampling.so");
        writeControl   timeStep;
        writeInterval  1;
        probeLocations
        (
            (0 9.95 19.77)
            (0 -9.95 19.77)
        );
        fixedLocations  false;
        fields
        (
            p
        );
    }
}

メッシュは以下の通りで、メッシュ数は25840です。

メッシュ メッシュ

計算結果は以下の通りです。

液相 water の界面
初期時刻での液相 water の界面 初期時刻での液相 water の界面
最終時刻での液相 water の界面 最終時刻での液相 water の界面
指定座標での圧力 p の時系列 指定座標での圧力 p の時系列。実線:座標 (0, 9.95, 19.77)、破線:座標 (0, -9.95, 19.77)。

実行コマンド

cp -r $FOAM_TUTORIALS/multiphase/interDyMFoam/ras/sloshingTank3D6DoF sloshingTank3D6DoF
cd sloshingTank3D6DoF

m4 system/blockMeshDict.m4 > system/blockMeshDict
blockMesh
cp 0/alpha.water.orig 0/alpha.water
setFields
interDyMFoam

# 3次元描画
paraFoam

# 指定座標での圧力の時系列をグラフ描画
gnuplot
gnuplot>plot "postProcessing/probes/0/p" using 1:2 title 'point 1' with lines lc rgb "black", "postProcessing/probes/0/p" using 1:3 title 'point 2' with lines dt (5,5) lc rgb "black"

計算時間

40分55.38秒 ※シングル、Inter(R) Core(TM) i7-2600 CPU @ 3.40GHz 3.40GHz

参照