綺麗な円柱メッシュの作成

円柱を立方体メッシュで切るには領域を分割するなど気を配らないといけない。
ただ、公開されているマクロを使えば簡単に切る事が出来る。

確認バージョン
OpenFOAM v1812

１．以下のURLからblockMeshDict.m4をダウンロードする
https://www.ehsanmadadi.com/cylinder-mesh/

２．blockMeshDict.m4をcaseディレクトリ内のsystem内に置く

３．m4がなければインストールする
sudo apt install m4

４．blockMeshDict.m4の設定を変更する

以下のように円柱メッシュは円形断面において曲率を持った四角形で分割されている。
そのため四角形についても形状を定義する

４．１　幾何学的形状

blockMeshDict.m4


// Inner square side half　内部の四角形の1辺の長さ
define(s, 0.35)

// Inner square side curvature　内部の四角形の膨らみ具合を表す。「内部の四角形の1辺の長さ」より若干大きいぐらいが良い。
define(sc, 0.4)

// cylinder radius　半径
define(r, 0.7)

// Height of cylinder　軸方向の長さ
define(z, 20.0)

４．２　メッシュサイズ

blockMeshDict.m4


// Number of cells at inner square 内部の四角形の一片の分割数
define(Ns, 20)

// Number of cells between inner square and circle　内部の四角形と円柱外面の間の分割数
define(Ni, 10)

// Number of cells in the cylinder height　円柱の軸方向の分割数
define(Nz, 100)

単位系をmmにするためconvertToMetersを使用する

convertToMeters 0.001

４．以下を実行
m4 blockMeshDict.m4 > blockMeshDict
cd ../
blockMesh

５．paraFoamかParaViewにてメッシュを確認する

mshメッシュの出力

caseファイル直下で以下を実行
foamMeshToFluent
caseファイル内にfluentInterfaceディレクトリが出来ているので、そのディレクトリの直下にmshファイルがあることを確認する
しかし、このmshファイルは以下のソフトでは読み込めなかった
GMesh
ElmerGUI, ElmerGrid
FreeCAD
Salome-Meca Win版