データの出力に関して

[Daiki Kono]

大学での研究において、OpenFOAMを使用しております。
２次元軸対称衝突噴流のよどみ流れにおいて壁面に働くせん断応力を算出するため、非常に簡単なモデルを作成し計算を行っています。
得られた速度場データから壁面せん断応力を求めるために、OpenFOAMの計算結果を出力して、MATLABで解析を行っているのですが、その出力方法についてご教授願いたいと思い、投稿いたしました。
現在、データの出力について知っていることは以下です。
・sampleDictを用いた出力。
・controlDict内のfunctionｓを用いた出力。
上記の方法でrawファイルを作成し、スペース区切りをタブ区切りに直し、それをexcelで並び替え、MATLABで読み込んでいます。現在のモデルは非常に簡単なので、excelで並び替えて、MATLABで読み込むという作業が可能なのですが、モデルが複雑になった時におそらく難しいと思います。また、現在の出力法では壁面のデータが得られていないので、その点も改善したいと思っています。

参考にcontrolDict内の記述を以下に示します。

application    simpleFoam;
startFrom        startTime;
startTime        0;
stopAt            endTime;
endTime        30000;
deltaT            1;
writeControl    timeStep;
writeInterval    1000
purgeWrite        2;
writeFormat    ascii;
writePrecision    6;
writeCompression    off;
timeFormat        general;
timePrecision    6;

runTimeModifiable    true;

libs
(   
    "libOpenFOAM.so""libgroovyBC.so"
    "libIOFunctionObjects.so"
    "libcloudFunctionObjects.so"
    "libfieldFunctionObjects.so"
    "libforces.so"
    "libFVFunctionObjects.so"
    "libjobControl.so"
    "libsystemCall.so"
    "libutilityFunctionObjects.so"
);

functions
{
    func_dist
    {
        type surfaces;
        fields (U p);
        interpolationScheme cellPoint;
        surfaceFormat raw;
        surfaces
        (
            dist
            {
            type plane;
            basePoint (0 0 0);
            normalVector (0 0 1);
            triangulate no;
            interpolate no;
            }
        );
        outputControl outputTime;
        outputInterval 1;
    }
}

[Masashi Imano]

今野です．

壁面せん断応力のようにCFD解析において重要な値は他のユーザも算出するはずなので，
"wall shear stress OpenFOAM"や"壁面せん断応力 OpenFOAM"で検索すれば，
OpenFOAMで壁面せん断応力を効率的に算出するのに必要な情報が得られると思います．

なお，壁面せん断応力の値は壁関数や乱流モデルにも依存することにご注意ください．

2015年1月25日 9:30 Daiki Kono <k.dai...@gmail.com>:

> --
> このメールは Google グループのグループ「OpenFOAM」に登録しているユーザーに送られています。
> このグループから退会し、グループからのメールの配信を停止するには openfoam+u...@googlegroups.com
> にメールを送信してください。
> このグループに投稿するには open...@googlegroups.com にメールを送信してください。
> http://groups.google.com/group/openfoam からこのグループにアクセスしてください。
> その他のオプションについては https://groups.google.com/d/optout にアクセスしてください。



--
IMANO Masashi, Ph.D.