24小时热门版块排行榜    

查看: 3418  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 070300化学学硕求调剂 +5 太想进步了0608 2026-03-16 5/250 2026-03-16 14:11 by xwxstudy
[教师之家] 焦虑 +7 水冰月月野兔 2026-03-13 9/450 2026-03-16 10:00 by Quakerbird
[硕博家园] 深圳大学硕士招生(2026秋,传感器方向,仅录取第一志愿) +4 xujiaoszu 2026-03-11 8/400 2026-03-16 09:45 by xujiaoszu
[考研] 中科院材料273求调剂 +3 yzydy 2026-03-15 3/150 2026-03-15 21:15 by ms629
[考研] 268求调剂 +5 一定有学上- 2026-03-14 6/300 2026-03-14 22:20 by 运气yunqi
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 290求调剂 +4 @将就将就看 2026-03-10 8/400 2026-03-14 14:23 by 千千运气
[考研] 环境调剂 +6 晓看天暮看云 2026-03-09 6/300 2026-03-14 01:16 by JourneyLucky
[考研] 311求调剂 +8 zchqwer 2026-03-10 8/400 2026-03-14 00:01 by JourneyLucky
[考研] 279求调剂 +3 Dizzy123@ 2026-03-10 3/150 2026-03-13 23:02 by JourneyLucky
[考研] (081700)化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 281求调剂 +9 Koxui 2026-03-12 11/550 2026-03-13 20:50 by Koxui
[考研] 【0856】化学工程(085602)313 分,本科学科评估A类院校化学工程与工艺,诚求调剂 +7 小刘快快上岸 2026-03-11 7/350 2026-03-13 16:06 by ruiyingmiao
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[考研] 304求调剂(085602一志愿985) +12 化工人999 2026-03-09 12/600 2026-03-13 12:02 by JourneyLucky
[考研] 277求调剂 +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by 白夜悠长
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考研] 321求调剂(食品/专硕) +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
信息提示
请填处理意见