CFD for kids - GeoGen2D

shirazbj

GeoGen2D是个前处理程序。

它用一个简单的文本文件描述流场边界及边界条件。然后通过该程序输出网格生成程序需要的输入文件，再用相应程序计算分析流场。
它输出gmsh的geo文件，供用gmsh划分网格。然后可用其他支持gmsh网格的程序计算流场，比如openfvm.
同时它也输出西交nht的2d非结构delauney网格划分教学程序的input.txt文件，供其划分非结构网格。然后用其非结构的教学程序计算流场。

文件下载后，解压倒c:\un123里。在该目录里执行1。2。3三个程序。

1.exe为本程序，GeoGen2D
2.exe为西交nht的2d非结构网格生成教学程序+gmsh网格文件输出
3.exe为西交nht的2d非结构流场计算教学程序+进出口边界+gmsh后处理文件输出
4.exe可以把迭代过程中保存的结果数据转化成gmsh的动画格式，在gmsh下看迭代过程。

流场边界及边界条件文件以绕流三角形为例：给出外边界4个点，内边界（三角形）三个点，每个边的边界条件，就可以开始算了。

case文件名只能为4个字母，比如叫做：tr13.txt

1.exe  要输入case文件名
2.exe  要指定区域内部内点的个数, 1000~2000个，运行完观察生成的网格，酌情加密。
3.exe 不用输入

参考：
gmsh  http://geuz.org/gmsh/
OpenFVM http://openfvm.sourceforge.net/
nht下载 http://nht.xjtu.edu.cn/nht/downpage.asp

[ 本帖最后由 shirazbj 于 2011-7-8 21:07 编辑 ]

shirazbj

算例三角形绕流的结果，对不对另说。呵呵。

shirazbj

un123.txt为计算程序的控制文件，只有3行：
100 - 迭代次数
1   - 对流项一阶迎风， =0 中心格式
0


grid.msh                      网格文件

*.geo                         gmsh的几何描述文件，可用gmsh划分网格

GMSH_ani_p.POS                p的迭代伪彩图动画
GMSH_ani_vel.POS              vel的迭代矢量图动画  
GMSH_bnd.POS                  流场边界
GMSH_FEILD00_VP.POS           流场初值，供检查进口设置
GMSH_FEILD_ST_p.POS           p伪彩图
GMSH_FEILD_VP_point.POS       vel矢量图，位于节点
GMSH_FEILD_VP_v.POS           vel矢量图，位于控制体中心
GMSH_outlet.POS               出口，显示位置，显示的数值没意义
GMSH_RES.POS                  迭代残差


使用方法：
只须在资源浏览器中把要看的gmsh文件托放到桌面上gmsh的快截图标上，gmsh就会启动，并打开相应的图，不用作任何其他操作，非常简便。
动画文件打开后，要点击左下角的箭头形播放图标。

[ 本帖最后由 shirazbj 于 2011-6-12 12:02 编辑 ]

shirazbj

要在三角形后边再加一个矩形，只要增加这几行就可以了：

4
20
4, 2,0   
5,2,0
5,4,0
4,4,0

再执行1，2，3，

结果如图

[ 本帖最后由 shirazbj 于 2011-6-4 18:32 编辑 ]

kaitian2012

谢谢楼主，学习了！

Gnimou

谢谢~~~

shirazbj

瞎算了一个4进口边，4出口边的例子。出口可以在y方向上了。边界也可以是圆了。
下一步希望进出口可以在斜边上。

[ 本帖最后由 shirazbj 于 2011-6-2 21:21 编辑 ]

shirazbj

进口在斜边上的例子，进口给法向和切向速度。出口不用指定系数了。

down.txt

[ 本帖最后由 shirazbj 于 2011-7-8 15:11 编辑 ]

shirazbj

出口在斜边上的例子

room.txt

[ 本帖最后由 shirazbj 于 2011-7-8 15:11 编辑 ]

shirazbj

进口和出口均在圆上的例子。

这个例子只算了很少的布数，算多了就发散了。

cir2.txt

[ 本帖最后由 shirazbj 于 2011-7-8 15:45 编辑 ]

xiaoshuai537

楼主很厉害呀
谢谢楼主分享

shirazbj

边界从文件读入的例子。这里模拟闸阀的四边形是从文本文件valve_gate.dat读入的。

valve_gate.dat内容如下：

1. 4
   
2. 1,0,0
   
3. 2,1,0
   
4. 3,1,0.05
   
5. 4,0,0.05
   

复制代码

shirazbj

ver_d.

shirazbj

这是读入naca_0012翼型的例子。

有翼型后，nht的分网格教学程序经常分不出来。用gmsh一般都分的很好。

需要加个远场条件了。目前不支持。