##############################################################
２次元円柱周り流れのシミュレーション (Navier-Stokes Eq.)
##############################################################

=========================================================
条件
=========================================================

* 2次元
* 解くべき方程式はNavier-Stokes方程式
* 領域は 縦 0.050 (m) 横 0.2 (m) の２次元矩形領域 ( (-0.1,-0.025)-(+0.1,+0.025) の領域 )
* (x,y)=(-0.4, 0.0)に直径 0.010 (m) の円柱が存在する．

=========================================================
メッシュ
=========================================================

* メッシュ生成コードは以下．

  .. literalinclude:: ../../code/fluid/flow__around_cylinder_NS_XY2D/mesh.py
     		    :caption:  ２次元円柱周り流れのメッシュ生成コード
     		    :language: python

  .. literalinclude:: ../../code/fluid/flow__around_cylinder_NS_XY2D/mesh.conf
     		    :caption:  mesh.conf

  .. literalinclude:: ../../code/fluid/flow__around_cylinder_NS_XY2D/phys.conf
     		    :caption:  phys.conf


* メッシュの様子は以下．

  .. image:: ../../fig/flow__around_cylinder_NS_XY2D/mesh.png
     :width:  600px
     :align:  center
  

=========================================================
シミュレーション .sif ファイル
=========================================================

* elmer シミュレーションコードは以下．

  .. literalinclude:: ../../code/fluid/flow__around_cylinder_NS_XY2D/ns.sif
                      :caption:  ２次元円柱周り流れのシミュレーションコード

                                 
                                 
* 境界条件は、

  + 左側＝上流面 :math:`v_x=0.05, v_y=0, v_z=0`
  + 右側＝下流面：圧力ゼロ ( :math:`p=0` )
  + 側面-固定癖：流速ゼロ ( :math:`v_x,v_y,v_z=0` )


  
=========================================================
シミュレーション結果
=========================================================

* シミュレーション結果の例．
* webサイトではカルマン渦ができているが、自分の環境ではカルマン渦ができていない．

  + おそらくレイノルズ数の関係．レイノルズ数が低くて（粘性を下げてもカルマン渦ができないことから数値粘性ではないかと思われる）
  + メッシュを細かくきる必要があるかも．