加载中…
个人资料
浪迹天大
浪迹天大
  • 博客等级:
  • 博客积分:0
  • 博客访问:5,657
  • 关注人气:8
  • 获赠金笔:0支
  • 赠出金笔:0支
  • 荣誉徽章:
相关博文
推荐博文
谁看过这篇博文
加载中…
正文 字体大小:

OpenFOAM中圆柱体的二维画法

(2014-12-25 10:22:36)
标签:

openfoam

wedge

makeaxialmesh

分类: OpenFOAM
首先声明,使用的方法是makeAxialMesh,makeAxialMesh工具,请自行去WIKI下载,编译openfoamwiki.net/index.php/Contrib/MakeAxialMesh
来张UserGuide上的插图:
OpenFOAM中圆柱体的二维画法

在OpenFOAM(其他平台不知道)中,圆柱体的二维表现为一个夹角不超过5度的楔形。那么怎样把这个楔形画出来呢?
在system文件夹下,需要两个文件,rotationDict and collapseDict

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
// ************************rotationDict************************************** //
// Use the same parameters as from the command line
// makeAxialOldMode 1;

// used in old and in new mode
makeAxialAxisPatch center;
makeAxialWedgePatch frontAndBack;

// used only in old mode
makeAxialOffset 0.01;

// used only in new mode
rotationVector (0 1 0);
originVector (0 0 0); //offset
// originVector (0 0); // origin

wedgeAngle 5;

// revolve option
// old and default mode, points are projected on wedges
// points are revolved
revolve 0;

// ************************************************************************* //

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
// *collapseDict* //

// If on, after collapsing check the quality of the mesh. If bad faces are
// generated then redo the collapsing with stricter filtering.


collapseEdgesCoeffs
{
    // Edges shorter than this absolute value will be merged
    minimumEdgeLength   1e-6;

    // The maximum angle between two edges that share point attached to
    // no other edges
    maximumMergeAngle   179;
}


// ************************************************************************* //



在blockMeshDict文件中,网格照常定义。只是需要注意一下几点:
1.把需要坍塌的那个面(此面合并成一条线,即圆柱体的轴线),定义为center , type symmetryPlane(定义边界条件是也用此type)。
2.需要frontAndBack。
3.生成楔形网格后,frontAndBack会变成两个    frontAndBack_pos, frontAndBack_neg。定义边界条件时,用wedge类型。

然后,按如下指令操作:
1
2
3
4
blockMesh
makeAxialMesh -overwrite
collapseEdges -overwrite
checkMesh
具体原理还没研究,以后补充……
也欢迎各位大神指点~

0

阅读 评论 收藏 转载 喜欢 打印举报/Report
  • 评论加载中,请稍候...
发评论

    发评论

    以上网友发言只代表其个人观点,不代表新浪网的观点或立场。

      

    新浪BLOG意见反馈留言板 电话:4000520066 提示音后按1键(按当地市话标准计费) 欢迎批评指正

    新浪简介 | About Sina | 广告服务 | 联系我们 | 招聘信息 | 网站律师 | SINA English | 会员注册 | 产品答疑

    新浪公司 版权所有