extrudeMeshを用いた曲がり円管メッシュ作成
915 views
Skip to first unread message
H.Kusaka
Nov 23, 2014, 4:20:16 AM11/23/14
to open...@googlegroups.com
お世話になっております．日下と申します．
先日は「buoyantPimpleFoamを用いた非圧縮、層流の円管内振動流による熱輸送解析」
で大変お世話になりました．おかげさまで直円管内振動流については理論解とよく一致する結果が得られました．
（これについては，近く上のスレッドで報告させていただくつもりです）
現在，解析対象を曲がり管に移すため，bolckMeshで曲がり管メッシュ作成を試みておりましたが，
CFD-ONLINEでも報告されているように，blockMeshで曲がりに沿った断面形状（円形）を維持できないというトラブルが発生してしまいました．
extrudeMeshを用いて押出でメッシュを作成すればよいとの投稿がされていましたので，
blockMeshで作成した曲がり管の入口面などを曲がった軸に沿って押し出せば良いかと思っていたのですが，
下記アドレスを参照すると押し出しに用いる軸は直線のみであり，曲がりに沿った軸を定義することは難しいのではないかと思います．
extrudeMeshを用いて曲がりに沿った軸を定義することで，曲がり管メッシュを作成することは可能なのでしょうか．
併せてケースファイルを添付させていただきます．
なお，現在のメッシュは管半径a=5 mm,管長L=314 mm, 曲率半径b=100 mmで，
曲率中心が原点，管の片端の中心が(b,0,0)になるように作成しております．
また簡単のため半断面についてのみ扱っています．
どんなに些細な事でも構いません．
皆様ご教授のほどどうかよろしくお願いいたします．
ohbuchi
Nov 25, 2014, 11:34:05 PM11/25/14
to open...@googlegroups.com
OpenFOAM-2.0.xでextrudeMeshが機能拡張され、曲がり管の様なメッシュも作成できる様になっています。
extrudeModeでwedgeを使います。
例
<blockMeshDict>
convertToMeters 1;
vertices
(
(4 -0.6 -0.5) // 0
(5 -0.6 -0.5) // 1
(5 0 -0.5) // 2
(4 0 -0.5) // 3
(4 -0.6 0.5) // 4
(5 -0.6 0.5) // 5
(5 0 0.5) // 6
(4 0 0.5) // 7
);
blocks
(
hex (0 1 2 3 4 5 6 7) (5 3 5) simpleGrading (1 1 1)
);
edges
(
);
boundary
(
inlet
{
type patch;
faces
(
(0 1 5 4)
);
}
walls
{
type wall;
faces
(
(0 3 2 1)
(4 5 6 7)
(1 2 6 5)
(0 4 7 3)
);
}
outlet
{
type patch;
faces
(
(2 3 7 6)
);
}
);
<extrudeMeshDict>
constructFrom mesh;
sourceCase "$FOAM_CASE";
sourcePatches (outlet);
flipNormals false;
extrudeModel wedge;
nLayers 50;
expansionRatio 1.0;
wedgeCoeffs
{
axisPt (0 0 0);
axis (0 0 1);
angle 90;
}
mergeFaces false;
2014年11月23日日曜日 18時20分16秒 UTC+9 H.Kusaka:
extrudeModeでwedgeを使います。
例
<blockMeshDict>
convertToMeters 1;
vertices
(
(4 -0.6 -0.5) // 0
(5 -0.6 -0.5) // 1
(5 0 -0.5) // 2
(4 0 -0.5) // 3
(4 -0.6 0.5) // 4
(5 -0.6 0.5) // 5
(5 0 0.5) // 6
(4 0 0.5) // 7
);
blocks
(
hex (0 1 2 3 4 5 6 7) (5 3 5) simpleGrading (1 1 1)
);
edges
(
);
boundary
(
inlet
{
type patch;
faces
(
(0 1 5 4)
);
}
walls
{
type wall;
faces
(
(0 3 2 1)
(4 5 6 7)
(1 2 6 5)
(0 4 7 3)
);
}
outlet
{
type patch;
faces
(
(2 3 7 6)
);
}
);
<extrudeMeshDict>
constructFrom mesh;
sourceCase "$FOAM_CASE";
sourcePatches (outlet);
flipNormals false;
extrudeModel wedge;
nLayers 50;
expansionRatio 1.0;
wedgeCoeffs
{
axisPt (0 0 0);
axis (0 0 1);
angle 90;
}
mergeFaces false;
2014年11月23日日曜日 18時20分16秒 UTC+9 H.Kusaka:
0 new messages