24小时热门版块排行榜    

查看: 3417  |  回复: 22
本帖产生 1 个 计算强帖 ,点击这里进行查看

zhanghuang10

捐助贵宾 (初入文坛)


[求助] 自编CFD程序怎样做后处理

大家好!
我现在在做CFD计算。网格是用ICEM划分的,将ICEM生成的文件读取后,导入自编的CFD程序进行流场的计算得出结果。计算结果的值(如速度、压力、温度)储存在每个网格点(node)上,请问我该怎样对得到的结果进行图像显示,得到流线、等高线图等信息?
谢谢
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

onesupeng

金虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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和看手册,很容易找到答案。主要是动手!!!
长期招收博士生,参见http://fsl-unsw.com
4楼2013-05-12 02:20:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

onesupeng

金虫 (职业作家)

【答案】应助回帖


csgt0: 金币+1, 多谢指导,欢迎交流 2013-05-14 16:13:27
引用回帖:
9楼: Originally posted by zhanghuang10 at 2013-05-13 10:13:19
结构网格和非结构网格的定义并不是根据它们形状规则与否定义的,而是根据它们的生成方式来定义的。一般来讲,结构化网格通过插值映射关系、偏微分方程等代数方法生成;而非结构化网格通过计算几何的方式生成。你如 ...

好吧,我不用跟你辩论了,等你确信你已经找到答案,告诉我一下。

另外,我图片是给你展示一个例子,并用最简单的列队示例进行说明,你都那样理解,那就没有办法了。

你既然推荐陶的书,那么你也看看第一章第三节,那儿有给你定义什么叫结构网格,什么叫非结构网格。我玩了结构网格和非结构网格快十年,今天第一次听说我把结构网格和非结构网格搞错了(是被一个cfd为何物都不清楚的人指出的。。。),我不知道该笑呢还是该哭亦或者是该倍感悲凉。。。
长期招收博士生,参见http://fsl-unsw.com
10楼2013-05-13 19:58:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

windzhu

木虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
可以采用很多后处理软件对计算出来的结果进行处理,如Tecplot、Paraview、Orign等。我最近在学IECM,倒想向您请教怎么学习IECM的,呵呵
锐气藏于胸,和气浮于脸,才气见于事,义气施于人
2楼2013-05-11 19:50:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhanghuang10

捐助贵宾 (初入文坛)


引用回帖:
2楼: Originally posted by windzhu at 2013-05-11 19:50:37
可以采用很多后处理软件对计算出来的结果进行处理,如Tecplot、Paraview、Orign等。我最近在学IECM,倒想向您请教怎么学习IECM的,呵呵

谢谢你的回复。
ICEM我也不是很熟悉,只会一些简单的操作,看看你需要掌握icem到怎样的程度,我们可以随时讨论。
现在我是自编的程序得到的计算结果,而且我需要求解的几何区域不规则。我现在的主要疑问是怎样安排输出数据的结构才能使后处理软件,比如tecplot能够读取,显示出速度、温度云图来。
3楼2013-05-11 20:54:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

supersunyan

金虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
如果数据比较大的话,可以用二进制文件格式,不然文件很大,tecplot中的帮助文件对数据格式有详细的说明,另外tecplot可以输入很多中间格式,你也可以参考的。
默默奋斗,实现梦想
5楼2013-05-12 09:43:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhanghuang10

捐助贵宾 (初入文坛)


引用回帖:
4楼: Originally posted by onesupeng at 2013-05-12 02:20:50
我列一个详细的tecplot的规矩给你看看吧。
打开open(unit=70,file='field.dat')
首先是结构网格的输出:
二维时:
write(70,*)'VARIABLES="X","Y","U","V","P&q ...

谢谢您!请问您所列的tecplot的保存规则是从哪里找的呢,有什么手册里写的有吗?
6楼2013-05-12 12:24:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhanghuang10

捐助贵宾 (初入文坛)


引用回帖:
4楼: Originally posted by onesupeng at 2013-05-12 02:20:50
我列一个详细的tecplot的规矩给你看看吧。
打开open(unit=70,file='field.dat')
首先是结构网格的输出:
二维时:
write(70,*)'VARIABLES="X","Y","U","V","P&q ...

我觉得你所指的结构网格和非结构网格和教科书上的定义不一样,教科书上所指的结构网格和非结构网格只是他们生成的方式不同,而你所指的结构网格和非结构网格我个人认为应该叫做均分网格和非均分网格。
7楼2013-05-12 12:58:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

onesupeng

金虫 (职业作家)

【答案】应助回帖

引用回帖:
7楼: Originally posted by zhanghuang10 at 2013-05-12 12:58:04
我觉得你所指的结构网格和非结构网格和教科书上的定义不一样,教科书上所指的结构网格和非结构网格只是他们生成的方式不同,而你所指的结构网格和非结构网格我个人认为应该叫做均分网格和非均分网格。...

直观讲,结构网格就是横七竖八规则排列的网格,就像军训列队一样。
       非结构网格就是除结构网格以外的所有网格。

下面这个文章图7A,外围像结构网格的样子,内部想非结构结构网格的样子。
http://challenge.vuse.vanderbilt ... ech_Review_2013.pdf
长期招收博士生,参见http://fsl-unsw.com
8楼2013-05-12 20:09:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhanghuang10

捐助贵宾 (初入文坛)


引用回帖:
8楼: Originally posted by onesupeng at 2013-05-12 20:09:05
直观讲,结构网格就是横七竖八规则排列的网格,就像军训列队一样。
       非结构网格就是除结构网格以外的所有网格。

下面这个文章图7A,外围像结构网格的样子,内部想非结构结构网格的样子。
http://challe ...

结构网格和非结构网格的定义并不是根据它们形状规则与否定义的,而是根据它们的生成方式来定义的。一般来讲,结构化网格通过插值映射关系、偏微分方程等代数方法生成;而非结构化网格通过计算几何的方式生成。你如果对一个规则的区域生成四边形的非结构化网格时,生成的网格也有可能是规则的。所以你根据形状的规则来与否来判断结构或者非结构化网格是不对的。建议你看看陶文铨写的《计算传热学近代进展》的网格生成那一章。
9楼2013-05-13 10:13:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhanghuang10 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 311求调剂 +4 26研0 2026-03-15 4/200 2026-03-16 14:46 by 可淡不可忘
[考研] 321求调剂 +4 大米饭! 2026-03-15 4/200 2026-03-16 08:41 by Linda Hu
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+5 NSFC2026我来了 2026-03-10 14/700 2026-03-15 23:53 by 不负韶华的虎
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 308求调剂 +4 是Lupa啊 2026-03-09 4/200 2026-03-14 02:06 by tranquil_ya
[考研] 一志愿北京化工大学材料与化工296分求调剂 +16 稻妻小编 2026-03-09 18/900 2026-03-14 02:00 by JourneyLucky
[考研] 一志愿华中农业大学071010,总分三百二,求调剂 +3 困困困困坤坤 2026-03-10 3/150 2026-03-14 00:35 by JourneyLucky
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 求材料调剂 +5 隔壁陈先生 2026-03-12 5/250 2026-03-13 22:03 by 星空星月
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 277求调剂 +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by 白夜悠长
[考研] 08食品或轻工求调剂,本科发表3篇sci一区top论文,一志愿南师大食品科学与工程 +3 我是一个兵, 2026-03-10 3/150 2026-03-13 10:21 by Yuyi.
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考研] 工科0856专硕化学工程269能调剂吗 +10 我想读研11 2026-03-10 10/500 2026-03-13 10:14 by Yuyi.
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生,2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
[考博] 读博申请 +5 感dd 2026-03-10 7/350 2026-03-11 17:02 by QGZDSYS
[考研] 调剂 +5 呵唔哦豁 2026-03-10 5/250 2026-03-10 22:00 by 28375m
信息提示
请填处理意见