２次元円柱周りの流れ計算で[[snappyHexを使って・・・という記事>]]があった。これまでsnappyで2次元メッシュということを考えたことがなかったが、確かに細分割指定しなければ、2次元メッシュも作成できるわなぁ。。。と新たな発見（？）をした思い。

⇒ DEXCS-OpenFOAMでチャチャッとやってみた。~
simpleFoamは定常計算なので、仮にカルマン渦が計算できたとして、物理的にどういう意味があるのか？という問題もありますが、計算パラメタを色々変えると結果がそれなりに変化して、数値計算の色んな事象を理解する助けにもなりそうです。

形状作成

–
-↑とりあえずデフォルトのモデルをオープン
–
-↑中味を全部削除し、外枠だけ残す
–
-↑スペースキーを押して、円柱を配置
–
–
-↑数値パッドを使って、円柱の位置を決める
–
-↑サイズを1/2（直径＝1）にしました⇒dexcs.stlとしてNamedAscii形式のエキスポートしておく
–
-↑外枠（基礎メッシュ）のサイズ（厚み）も円柱と同じになるようにした
-↓外枠（基礎メッシュ）の位置をずらした（下右）

基礎メッシュ作成

–
-blockMeshExporter GUI を使って、blockMeshDictを作成（分割数は、ｘ方向１００、、ｙ方向１、ｚ方向５０）
–
-blockMesh作成⇒autoPatch実行して、とりあえずパッチの名前を表示確認
–
-↑createPatchDictを編集。上下面をwall2、サイド面をwallに変更。wall2のtypeはemptyに。
–
-↑createPatchを実行して、パッチの名前を表示確認

snappyHexメッシュ作成

–
-↑ポイントは、円柱（dexcs.stl）をにチェックマークを付けるが、細分割指定をしない（Smax,Smin=0）ということ。レイヤーメッシュは付加すること可能。
–
-↑基礎メッシュの一を動かしたので、メッシュ制御点も変更

–
-↑メッシュ確認。メッシュ数：4,988、節点数：10,300作成に要した計算時間はほんの数秒。

流体特性

–
-↑動粘性係数は1.0e-4
–
-↑乱流モデルは使用しない

境界条件

-DEXCSのデフォルトから境界名が変わっているので、patch名取得ボタンを押してから境界条件を設定
–

-乱流モデルは使用しないが、ｋ、epsillonの初期値は必要
-流入風速を0.1としたので、レイノルズ数1000の流れ計算ということになる~
（ Re=UD/ν ： U=0.1,D=1,ν=1.0e-4）