应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 自编CFD程序怎样做后处理
52/1返回列表
查看: 3196  |  回复: 22
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 计算强帖 ，点击这里进行查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zhanghuang10

捐助贵宾 (初入文坛)

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

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

» 猜你喜欢
1楼 2013-05-11 18:28:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhanghuang10

捐助贵宾 (初入文坛)
引用回帖:
10楼: Originally posted by onesupeng at 2013-05-13 19:58:35
好吧，我不用跟你辩论了，等你确信你已经找到答案，告诉我一下。

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

你既然推荐陶的书，那么你也看看第一章第 ...

谢谢您的指正。CFD我确实才入门，有很多东西还不知道。但是CFD不是说谁用的久谁就熟，并且谁也不敢保证说在CFD领域他什么都清楚。我是带着求教的态度和您讨论的，您也没必要来挖苦我不懂CFD之类的。
您跟我解释了什么是结构化网格和非结构化化网格，谢谢您的解释。
但是您的解释和教科书上不一样，我提出疑问，这是很正常的。
正如您说，陶先生《计算传热学近代进展》第一章第3节P13页：
“在结构化网格中，任意一个位置可以通过一定的规则予以命名。”
“在非结构化网格中，节点的位置无法用一个固定的法则予以有序地命名”
这样的描述中，没有一个是和网格规则与否的描述相关的。陶老师上课的时也从没讲过结构化网格和非结构化网格是按规不规则区分的，而且用ICEM对一个正方形区域画网格，用结构化和非结构化网格画出来是一样的。
您说您做CFD十年了，请问您在Compter&Fluids, Journal of Computational Physics, JFM，PRL，PRE上发过文章嘛？或者为Ansys公司开发过什么大型的CFD程序嘛？
如果有，说明您这十年确实值得。如果没有，那您这十年也就白费了。
赞一下 回复此楼
11楼2013-05-13 20:51:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 23 个回答

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

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和看手册，很容易找到答案。主要是动手！！！
赞一下(4人) 回复此楼
长期招收博士生，参见http://fsl-unsw.com
4楼2013-05-12 02:20:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 23 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭