円筒内の流れの解析にて

[q10...@shibaura-it.ac.jp]

皆様こんにちは

酒井です．

　この前の円筒内流れについてはシミュレーションは正しい結果を得ることができましたが，ソルバをsimpleFoamに変更して解析を行うと添付ファイルみたいな解析結果になり得たい結果になりません．円筒の長さは2500mm，直径100mm，動粘性係数1e-6，速度をinletから0.001m/s，レイノルズ数は100です．また，初期条件としてnutとnuTidaのフォルダを以下のように設定しました．

　また，コマンドはblockMesh→snappyHexMesh(overwrite)→decomposePar→mpirun -np 4 simpleFoam -parallel→reconstructPar→paraFoamです．

//nutフォルダ--------------------------------------------

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

internalField   uniform (0 0 0);

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

internalField   uniform 0;

boundaryField

{

    inlet_1_solid       

    {

        type            freestream;

freestreamValue uniform 0;

    }

    outlet_1_solid    

    {

        type            freestream;

        freestreamValue uniform 0;

    }

    side_1_solid      

    {

        type            nutUSpaldingWallFunction;

value uniform 0;

    }

    defaultFaces

    {

        type            zeroGradient;

    }

//-------------------------------------------------------------

//nuTidaフォルダ--------------------------------------

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

internalField   uniform 0;

boundaryField

{

    inlet_1_solid       

    {

        type            freestream;

freestreamValue uniform 0;

    }

    outlet_1_solid    

    {

        type            freestream;

        freestreamValue uniform 0;

    }

    side_1_solid      

    {

        type            fixedValue;

value uniform 0;

    }

    defaultFaces

    {

        type            zeroGradient;

    }

//-------------------------------------------------------------

　以前にもこのような投稿があったと思うのですが，私が検索すると検索結果なしになっていました．投稿の回答としては縮流を考慮しなければならないことから計算領域を伸ばした方が良いと書いてあったと思いますが，具体的な方法がわかりませんでした．

　大変恐縮ですが，どのような方法をとればいいのかご教示をお願いします．また長文になってしまい申し訳ありません．

[ohbuchi]

こんにちは。

まず、Re=100だと層流ですので、nut, nuTildaは不要です。

また、freeStreamの使い方も不適切だと思います。今回関係ありませんが、nutの出入口

境界は type calculatedが一般的です。

速度分布が十分に発達しないのは、入り口速度分布が壁直上まで有限値になっているのが

悪影響していると思います。

壁で０になる適当な速度分布を与えて計算してみて下さい。

以上、ご参考まで。

2013年8月18日日曜日 16時43分37秒 UTC+9 q10...@shibaura-it.ac.jp:

[q10...@shibaura-it.ac.jp]

返答ありがとうございます．

最初に記述し忘れましたが，今回は層流を想定して解析を行いました．

nut,nuTidaの修正点についてもありがとうございます．この初期条件は乱流時に使用するのですね．

次回からnutとnuTidaを抜いて考えていきます．

これから，色々と勉強していきますが，また不明点がありましたらよろしくお願い致します．

2013年8月18日日曜日 20時00分22秒 UTC+9 ohbuchi:

2013年8月18日日曜日 20時00分22秒 UTC+9 ohbuchi:

2013年8月18日日曜日 20時00分22秒 UTC+9 ohbuchi:

2013年8月18日日曜日 20時00分22秒 UTC+9 ohbuchi: