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

无忧尊者

我感觉也是。这需要判断边条是I或者J，在另外一个方向循环次数减一。不过无所谓了，因为这样做是处理了第一层虚网格，但是虚网格有可能会在BUFFER里面重新更新，不带影响

北羽

1.5.4版本的SST k-w模型计算时，残差显示NAN，不知道是我编译的问题，还是程序里面有bug。希望李老师解答。
用的网格文件盒配置文件就是压缩包里面的RAE2822算例的文件。
测试发现，是OpenMP的问题。不采用OpenMP时可以正常计算。但使用OpenMP时，如何修改还未找到决解方案。
已经解决了，供大家参考。
将原来的语句
!$OMP PARALLEL DEFAULT(PRIVATE) SHARED(nx,ny,B,Kt,Wt,f1,Qk,Qw,Fluxk,Fluxw,d,uu,v)
更改为
!$OMP PARALLEL DEFAULT(PRIVATE) SHARED(nx,ny,B,Kt,Wt,f1,Qk,Qw,Fluxk,Fluxw,d,uu,v,Re)

最后，差分模块Residual_FDM子程序里面的的OpenMP也有问题，希望李老师解答，也希望已经成功解决该问题的朋友不吝赐教，先谢谢各位！

[ 本帖最后由 北羽 于 2013-10-31 09:15 编辑 ]

lixl-imech

多谢指出程序的Bug.
使用OpenMP编程时，在SHARED列表中容易忽略一些公共变量 （如Re等）。
另外一种方法是，将PRIVATE改为FIRSTPRIVATA。 这样即使Re被当成private属性的变量，由于具有firstprivate属性，也可以获得正确的值。

原帖由 北羽 于 2013-10-31 08:41 发表

                               
登录/注册后可看大图

1.5.4版本的SST k-w模型计算时，残差显示NAN，不知道是我编译的问题，还是程序里面有bug。希望李老师解答。
用的网格文件盒配置文件就是压缩包里面的RAE2822算例的文件。
测试发现，是OpenMP的问题。不采用OpenMP时 ...

tianhao007

请教李老师，gridgen生成的plot3d网格为grd格式，而您在getbcin中要求为plot3d 为.dat格式，请问这一问题应该怎么处理。谢谢！

apango

请问一下，壁面上的剪切应力应该如何输出？

liuluoqin

李老师，您好！
我在尝试使用您的opencfd_ec2d_v1.5.4程序时，出现如下错误，不知道该如何解决，乞盼您能给予帮助：
ifort -openmp -o ompexec opencfd_ec2d_v1.5.4.f90
sub_NS_singlegrid.f90(42): error #6404: This name does not have a type, and must have an explicit type.   [M]
!$OMP PARALLEL DO DEFAULT(SHARED) PRIVATE(i,j,m)
----------------------------------------------^
sub_NS_singlegrid.f90(42): error #7656: Subobjects are not allowed in this OpenMP clause; a named variable must be specified.   [M]
!$OMP PARALLEL DO DEFAULT(SHARED) PRIVATE(i,j,m)
----------------------------------------------^
sub_NS_singlegrid.f90(50): error #7656: Subobjects are not allowed in this OpenMP clause; a named variable must be specified.   [M]
!$OMP PARALLEL DO DEFAULT(SHARED) PRIVATE(i,j,m)
----------------------------------------------^
compilation aborted for opencfd_ec2d_v1.5.4.f90 (code 1)

祝好

liuluoqin

我直接将m删除，问题就解决了！

Entropy

For Opencdf-SC2D, is it possible to get a copy of source code including compact schemes ?  coz  it is "not support yet " in 1.52 version. sorry for chinese type method disable.

chfenger

向大神致敬！

hust_illusion

将主文件中的include语句注释掉，在项目中手动添加进所有需要用到的文件。

wilson732

请问3D版的软件采用OpenMP编程吗？谢谢

Aeronpu

作为一个后来人，刚开始上手，谢谢李老师！

Aerodynamics

顶顶顶！李老师功在当代，利在千秋！！！

tigermeng11

ch06 发表于 2011-1-19 09:56
我现在是在二维程序的基础上，算轴对称问题，因此我程序中的y应该对应成您图中的r
由于是二维，seta只有0 ...

请问，把EC二维程序编程轴对称程序后，除了边界条件的变化，还应该添加了坐标变换产生的源项吧？那么请问你是如何处理的？我加上源项后程序总是出现负温度负压力，不能运行

tigermeng11

lixl-imech 发表于 2011-1-18 09:32
普通的对称边界条件容易处理，无粘流的对称边界条件与滑移壁面条件是相同的。
轴对称边界条件要复杂 ...

李老师，如果EC程序用有限体积法的话，物理量本来就储存在网格中心，那么轴对称问题的话是不是就不存在奇点的问题了？轴对称边界可以直接用对称边界了？