24小时热门版块排行榜    

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

zhanghuang10

捐助贵宾 (初入文坛)


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

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

» 收录本帖的淘帖专辑推荐

仿真模拟

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

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的回帖

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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 08工学调剂 +4 用户573181 2026-03-20 4/200 2026-03-20 10:32 by 朗月清风1
[考研] 材料专硕英一数二306 +6 z1z2z3879 2026-03-18 6/300 2026-03-20 08:49 by xingguangj
[考研] 一志愿苏州大学材料求调剂,总分315(英一) +3 sbdksD 2026-03-19 3/150 2026-03-19 23:21 by fmesaito
[考研] 324分 085600材料化工求调剂 +3 llllkkkhh 2026-03-18 3/150 2026-03-19 14:22 by houyaoxu
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +3 @taotao 2026-03-19 3/150 2026-03-19 14:07 by peike
[考研] 0703化学调剂,求各位老师收留 +10 秋有木北 2026-03-14 10/500 2026-03-19 05:52 by anny19840123
[考研] 330求调剂 +3 小材化本科 2026-03-18 3/150 2026-03-18 21:55 by 无懈可击111
[考研] 354求调剂 +4 Tyoumou 2026-03-18 7/350 2026-03-18 21:45 by Tyoumou
[考研] 材料专业求调剂 +5 hanamiko 2026-03-18 5/250 2026-03-18 20:19 by 星空星月
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 331求调剂(0703有机化学 +7 ZY-05 2026-03-13 8/400 2026-03-18 14:13 by 007_lilei
[考研] 312求调剂 +8 陌宸希 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
[考研] 268求调剂 +6 简单点0 2026-03-17 6/300 2026-03-18 09:04 by 无际的草原
[考研] 考研化学学硕调剂,一志愿985 +4 张vvvv 2026-03-15 6/300 2026-03-17 17:15 by ruiyingmiao
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
信息提示
请填处理意见