发布OpenCFD-EC2D; 更新 Ver 1.5.4 ； 三维软件使用协议

ch06

原帖由 lixl-imech 于 2012-6-9 11:13 发表

                               
登录/注册后可看大图

如果采用等距网格，则不需读取网格文件。
但平板边界层流动应当采用非均匀的边界层网格为好，在近壁区网格需要非常密集。

我用李老师提供的mesh.f90生成了一个矩形计算域，但是初始化时壁面的速度与主流大小相等方向相反


f_no, face, istart,iend, jstart, jend,  neighb, subface orient
     1     1     1      1     1     17    -20     0     0
     2     2     1     49     1     1     -10     0     0
     3     3     49   49     1     17    -20     0     0
     4     4     1    49     17    17    -10     0     0
这个是我改的边界条件

[ 本帖最后由 ch06 于 2012-6-9 13:12 编辑 ]

ch06

这个是流线图

ch06

if(Twall .le. 0) then   ! 绝热壁
B%U(1,i2,j2)= B%U(1,i1,j1)       ! d(0)=d(1)   对称
           B%U(2,i2,j2)=-B%U(2,i1,j1)       ! u(0)=-u(1)  -> d*u 反对称
           B%U(3,i2,j2)=-B%U(3,i1,j1)       ! v(0)=-v(1)  -> d*v 反对称
           B%U(4,i2,j2)= B%U(4,i1,j1)       ! E(0)=E(1)   对称

李老师壁面条件，为什么要这样处理，没看明白这句：B%U(2,i2,j2)=-B%U(2,i1,j1)
是要抵消壁面的速度吗？

ch06

原帖由 lixl-imech 于 2012-6-9 11:13 发表

                               
登录/注册后可看大图

如果采用等距网格，则不需读取网格文件。
但平板边界层流动应当采用非均匀的边界层网格为好，在近壁区网格需要非常密集。

谢谢，李老师指点，目前是先试着使用openec，网格还不知如何根据需要调整。

ch06

原帖由 entropies 于 2012-5-26 03:41 发表

                               
登录/注册后可看大图

Try this and see what happens.

CentOS_OpenCFD算的结果居然出现分离了（用openmp算的结果）

[ 本帖最后由 ch06 于 2012-6-13 12:59 编辑 ]

ch06

openec就没有这种现象

entropies

Thank you. Let me look into that. By the way, how many iterations you've run for this case on centos_opencfd?

entropies

By the way, how did you plot the streamlines? By using the streamtraces in Tecplot and click on the plot?

ch06

原帖由 entropies 于 2012-6-13 17:51 发表

                               
登录/注册后可看大图

Thank you. Let me look into that. By the way, how many iterations you've run for this case on centos_opencfd?

54000步

ch06

原帖由 entropies 于 2012-6-13 17:57 发表

                               
登录/注册后可看大图

By the way, how did you plot the streamlines? By using the streamtraces in Tecplot and click on the plot?

是的，我用tecplot做的流线图

ch06

原帖由 ch06 于 2012-6-9 11:36 发表

                               
登录/注册后可看大图

我用李老师提供的mesh.f90生成了一个矩形计算域，但是初始化时壁面的速度与主流大小相等方向相反


f_no, face, istart,iend, jstart, jend,  neighb, subface orient
     1     1     1      1     1      ...

在壁面出现这个问题，应该是我输出ghost点问题，去掉ghost的值就好了。

entropies

你有没有QQ或者EMAIL？ 还有你用的是gfortran?

ch06

已经给发qq邀请了，我这儿的短消息打开时，很多时候无法显式！

[ 本帖最后由 ch06 于 2012-6-14 07:07 编辑 ]

ch06

原帖由 entropies 于 2012-6-14 06:59 发表

                               
登录/注册后可看大图

你有没有QQ或者EMAIL？ 还有你用的是gfortran?

是的，我用的是gfortran

entropies

try to use intel fortran, it's free to download from the website if you're running linux.

gfortran is a little bit fuzzy, I don't know why