ブジネさん : Boussinesq approx

-----------------
今日辺りから，また，寒くなるとか．

-----------------
BBC ニュース Olympic に関して．
 judo の話とか．
-----------------
-----------------
こないだ，このページに書いた，でんでらりゅうば ．なるほど．
 たしかに，北部九州人は，「とりよーと？」やなぁ．
でも，長崎の じいちゃん，ばーちゃん  の会話とかを思い出すと，
「とっとっと」でもオッケイかなぁと，同意．
-----------------

buoyantBoussinesqPimpleFoam

について考えてみる．

0/ ディレクトリに必要な初期変数は

U, p, k, epsilon と，
T, kappat, p_rgh の計７つ ．

U, p, k, epsilon は simpleFoam と同じ．


計算を始めると，nut, phi  の２変数が勝手に生成される．
中間変数なんだろう．

/home/hoge/OpenFOAM/OpenFOAM-2.1.x/applications/solvers/heatTransfer/buoyantBoussinesqPimpleFoam/buoyantBoussinesqPimpleFoam.C

の但し書きを読んでみる．

 Description
    Transient solver for buoyant, turbulent flow of incompressible fluids

    Uses the Boussinesq approximation:
    \f[
        rho_{k} = 1 - beta(T - T_{ref})
    \f]

    where:
        \f$ rho_{k} \f$ = the effective (driving) kinematic density
        beta = thermal expansion coefficient [1/K]
        T = temperature [K]
        \f$ T_{ref} \f$ = reference temperature [K]

    Valid when:
    \f[
        rho_{k} << 1
    \f]

-----------------

色々精査しながら，気がつく．

CASE02/constant/g

のなかで， 重力加速度が gy=-9.8 になってた．gz=-9.8 とすべきところを．

チュートリアルサンプルでは y 方向が上になってるから，
それを流用したため．

結果表示見て，なんかおかしいなぁとおもってた．


ということで，やり直した結果動画．

なんかよくわからないけど，単純なベクトル表示を行った ogv ファイルは
さくっと，YouTube アップロードができる．
-----------------

で，流線の動画も．

-----------------

-----------------

-----------------

-----------------