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

24小时热门版块排行榜    

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

zhanghuang10

捐助贵宾 (初入文坛)

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

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

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

仿真模拟

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-05-11 18:28:44
已阅   回复此楼   关注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的回帖

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

onesupeng

金虫 (职业作家)

【答案】应助回帖

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

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

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

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

onesupeng

金虫 (职业作家)
引用回帖:
11楼: Originally posted by zhanghuang10 at 2013-05-13 20:51:15
谢谢您的指正。CFD我确实才入门，有很多东西还不知道。但是CFD不是说谁用的久谁就熟，并且谁也不敢保证说在CFD领域他什么都清楚。我是带着求教的态度和您讨论的，您也没必要来挖苦我不懂CFD之类的。
您跟我解释了 ...

我在第8楼中的定义是这样的，是为了直观的表述我的理解：
"直观讲，结构网格就是横七竖八规则排列的网格，就像军训列队一样。
       非结构网格就是除结构网格以外的所有网格。"

你觉得有什么问题？再直接一点，可以用i，j，k表述的三个维度的坐标体系，他就是结构网格。。。如我给你的tecplot的前两种，否则他就是非结构，如我给你的tecplot标记的后述情形。请别继续说那1维只能用i不能用ijk这样的话。

结构和非结构当然不以几何图形为标记，但确实以编号为标记，一个结构网格完全可用非结构网格的方式标记，这没什么好较真的。反之却不一定。

至于我的水平如何，那我只能说，你导师一个研究组一年的发文数量和质量，不见得比我高。谢谢。

8楼给出的文章，是我其中一个，欢迎批评指正，谢谢。
赞一下 回复此楼
长期招收博士生，参见http://fsl-unsw.com
12楼2013-05-13 21:46:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

onesupeng

金虫 (职业作家)
http://www.cfd-online.com/Forums ... tructured-mesh.html

我附上一个帖子，第三楼讲的比较详细。和上面第二段、第三段是一个意思。
赞一下 回复此楼
长期招收博士生，参见http://fsl-unsw.com
13楼2013-05-13 22:06:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

onesupeng

金虫 (职业作家)
引用回帖:
14楼: Originally posted by pswpswpsw at 2013-12-15 02:21:36
你这样和别人交流的态度。可以看得出你在平时生活中是一个怎样的人了。

我只能说，你要先学如何做人，再学好做学问。

不过，学生不教，过错不在于你，而在于你的导师。...

前辈教训的太有道理了，按照您的意思，您觉得我应该怎样帮助别人呢？您在生活中又如何为我等不会做人不会做学问的人做表率的呢？让我等有个标杆嘛，要不我等生活没有目标，谢谢。
赞一下 回复此楼
长期招收博士生，参见http://fsl-unsw.com
15楼2013-12-15 05:07:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

onesupeng

金虫 (职业作家)
引用回帖:
14楼: Originally posted by pswpswpsw at 2013-12-15 02:21:36
你这样和别人交流的态度。可以看得出你在平时生活中是一个怎样的人了。

我只能说，你要先学如何做人，再学好做学问。

不过，学生不教，过错不在于你，而在于你的导师。...

另外，您老可以看看最开始的回帖，您觉得哪里可以改进？怎么改进？谢谢。
赞一下 回复此楼
长期招收博士生，参见http://fsl-unsw.com
16楼2013-12-15 05:10:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

onesupeng

金虫 (职业作家)
引用回帖:
17楼: Originally posted by lan5902108 at 2013-12-15 21:57:26
顶你。
...

雪中送炭啊~
回复此楼
长期招收博士生，参见http://fsl-unsw.com
18楼2013-12-16 10:23:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

onesupeng

金虫 (职业作家)
引用回帖:
19楼: Originally posted by xjman at 2015-03-03 12:07:57
支持一记，不过我觉得楼上的可能回复错了人。
不要怪我挖坟，我是搜索非结构网格的连接关系挖出来的。我的tecplot里面居然没有数据的例子。
另外，e=Equd(ee)和e=Etri(ee)是您自己写的子程序吧？...

就是三角形或者四边形单元的单元编号，比如有n个三角形单元和m各四边形单元，第ee个单元在总单元中的位置。如果只有三角形或者四边形，XX（ee）=ee，1、2、3、4、5这样子。。。
赞一下 回复此楼
长期招收博士生，参见http://fsl-unsw.com
20楼2015-03-03 13:57:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

onesupeng

金虫 (职业作家)
引用回帖:
21楼: Originally posted by xjman at 2015-03-04 04:05:33
非常感谢！您这个我还是没看明白，但我试着把单元的每个顶点编号写出来，可以画出云图了。不过不知道用程序怎么写，网格数如果多了的话手动写就不可能了。
你这写的是与该单元相邻的单元编号吗？三角形单元和四边 ...

我的做法是分开的。但是应该可以，比如三角形单元，你当成四边形，第四个点和第三个点一样试一试。

至于程序怎么去实现，属于比较初级层面的问题，可能我没空具体去教你怎么搞
赞一下 回复此楼
长期招收博士生，参见http://fsl-unsw.com
22楼2015-03-21 18:10:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhanghuang10 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭