OpenFoamチュートリアルcavityについて
cavity流れについては以下の記事で書いていますので、まずはこちらを参照ください。＞＞＞OpenFoamではじめての流体解析 まずはキャビティ流れを確認する
前回の記事と今回のゴール
前回までの記事で、cavityをカスタマイズして左から右へ流体が流れるようにしました。こちらの記事を参照ください。本記事ではこの記事に記載しているコードをベースにしています＞＞＞OpenFOAMチュートリアルCavityカスタマイズ：流体を左から右へ流す方法
こんな感じで左が流入側、右が流出側です。上下は壁です。
今回はこの正方形の解析領域の形状を変更します。
尚、本シリーズはOpenFoamの使い方に慣れることを目的としているので、解析の正確性などを保証するものではなく、まずは設定を変更して流体の流れを変えること。この実現をゴールとします。（試行錯誤して行っていますので余計な記載や、変更が十分でないところがあることご了承願います。また環境やバージョンによって結果が異なる可能性もあります。）
cavityのファイル構成を確認
まずはcavityのファイル構成を確認しましょう。
フォルダは「0」「constant」「system」の3つです。
「0」フォルダ：初期条件のファイルを含むフォルダです。初期速度、圧力など。
「constant」フォルダ：定数パラメータを含むフォルダです。物理的な特性など。
「system」フォルダ：計算の制御や設定を含むフォルダです。時間ステップの設定、ソルバーの選択、解析制御ファイルの作成方法など。
system – blockMeshDict
systemフォルダにあるblockMeshDicファイルを編集していきましょう。変更するところはedgesに必要な記載を追加するだけです。
vertices
(
(0 0 0)
(1 0 0)
(1 1 0)
(0 1 0)
(0 0 0.1)
(1 0 0.1)
(1 1 0.1)
(0 1 0.1)
);
blocks
(
hex (0 1 2 3 4 5 6 7) (20 20 1) simpleGrading (1 1 1)
);
edges
(
arc 0 1 (1 -1 0.1)
arc 4 5 (1 -1 0.1)
polyLine 2 3 ((0.5 0.7 0))
polyLine 6 7 ((0.5 0.7 0.1))
);
verticesで指定した頂点の0と1、4と5は円弧で繋いで、2と3、6と7はポリラインで繋いでみました。円弧は内挿点を設定、ポリラインはリストで書きます。言葉では分かりにくいかもしれませんが、つぎの画像のような形状になります。
先日の面の結合とこのedgesによる形状変更を組み合わせると色々な形状が作れそうです。ただ、解析はうまくいかなかったので、その部分については工夫が必要かもしれません。
