| 查看: 3196 | 回复: 22 | |||
| 本帖产生 1 个 计算强帖 ，点击这里进行查看 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
自编CFD程序怎样做后处理
|
|||
|
大家好！ 我现在在做CFD计算。网格是用ICEM划分的，将ICEM生成的文件读取后，导入自编的CFD程序进行流场的计算得出结果。计算结果的值（如速度、压力、温度）储存在每个网格点(node)上，请问我该怎样对得到的结果进行图像显示，得到流线、等高线图等信息？ 谢谢 |
回复此楼» 猜你喜欢
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
物理学I论文润色/翻译怎么收费?
已经有164人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有23人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
onesupeng
金虫 (职业作家)
- 计算强帖: 13
- 应助: 256 (大学生)
- 贵宾: 1.36
- 金币: 2336.2
- 散金: 9224
- 红花: 92
- 帖子: 4583
- 在线: 1303.8小时
- 虫号: 394701
- 注册: 2007-06-07
- 专业: 流体力学
|
我在第8楼中的定义是这样的，是为了直观的表述我的理解： "直观讲，结构网格就是横七竖八规则排列的网格，就像军训列队一样。 非结构网格就是除结构网格以外的所有网格。" 你觉得有什么问题？再直接一点，可以用i，j，k表述的三个维度的坐标体系，他就是结构网格。。。如我给你的tecplot的前两种，否则他就是非结构，如我给你的tecplot标记的后述情形。请别继续说那1维只能用i不能用ijk这样的话。 结构和非结构当然不以几何图形为标记，但确实以编号为标记，一个结构网格完全可用非结构网格的方式标记，这没什么好较真的。反之却不一定。 至于我的水平如何，那我只能说，你导师一个研究组一年的发文数量和质量，不见得比我高。谢谢。 8楼给出的文章，是我其中一个，欢迎批评指正，谢谢。 |
12楼2013-05-13 21:46:08
windzhu
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 6688.5
- 散金: 204
- 红花: 3
- 帖子: 173
- 在线: 137.8小时
- 虫号: 2204544
- 注册: 2012-12-25
- 性别: GG
- 专业: 流体力学
2楼2013-05-11 19:50:37
3楼2013-05-11 20:54:52
onesupeng
金虫 (职业作家)
- 计算强帖: 13
- 应助: 256 (大学生)
- 贵宾: 1.36
- 金币: 2336.2
- 散金: 9224
- 红花: 92
- 帖子: 4583
- 在线: 1303.8小时
- 虫号: 394701
- 注册: 2007-06-07
- 专业: 流体力学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
csgt0: 金币+2, 计算强帖+1, 多谢指导 2013-05-13 11:39:30
zhanghuang10: 金币+5, ★★★很有帮助 2013-05-13 21:09:14
感谢参与，应助指数 +1
csgt0: 金币+2, 计算强帖+1, 多谢指导 2013-05-13 11:39:30
zhanghuang10: 金币+5, ★★★很有帮助 2013-05-13 21:09:14
|
我列一个详细的tecplot的规矩给你看看吧。 打开open(unit=70,file='field.dat') 首先是结构网格的输出： 二维时： write(70,*)'VARIABLES="X","Y","U","V","P"' write(70,*)'ZONE F=POINT, I=',nx,', J=',ny do j=1,ny do i=1,nx write(70,*)xc(i),yc(j),u(i,j),v(i,j),p(i,j) enddo enddo 三维类似： write(70,*)'VARIABLES="X","Y","Z","U","V","W","P"' write(70,*)'ZONE F=POINT, I=',nx,', J=',ny,' K=',nz do k=1,nz do j=1,ny do i=1,nx write(70,*)xc(i),yc(j),zc(k),u(i,j,k),v(i,j,k),w(i,j,k),p(i,j,k) enddo enddo enddo 非结构网格： 二维4边形单元 write(70,*)'variables= "x" "y" "u" "v" "p" "w" ' write(70,*)'zone f=fepoint et=quadrilateral n=',np,' e=',Nequd do i=1,np write(70,*)x(i),y(i),u(i),v(i),p(i),w(i) enddo do ee=1,Nequd e=Equd(ee) write(70,*)(ep(e,j),j=1,4) enddo 三角形单元： write(70,*)'zone f=fepoint et=triangle n=',np,' e=',Netri do i=1,np write(70,*)x(i),y(i),u(i),v(i),p(i),w(i) enddo do ee=1,Netri e=Etri(ee) write(70,*)(ep(e,j),j=1,3) enddo 其他情况类似。多google、baidu和看手册，很容易找到答案。主要是动手！！！ |
4楼2013-05-12 02:20:50