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

ch06

好像是gfortran要求每个变量都要初始化，intel fortran对没有初始化的变量赋值0

entropies

我不是很清楚 gfortran 的工作方式，看来代码还得写地更严谨一些，以使得可移植性更好。

好像这里有个链接可以检测一下数据：
http://adl.stanford.edu/docs/display/SUSQUARED/Tutorial+6+-+Turbulent+RAE+2822

ch06

Initialize OK ......
  Start ......
*** glibc detected *** ./opencfd: double free or corruption (out): 0x00007f67a0a0e0a0 ***
*** glibc detected *** ./opencfd: double free or corruption (out): 0x00007f679b7fa0a0 ***
*** glibc detected *** *** glibc detected *** ./opencfd: double free or corruption (out): 0x00007f679aff80a0 ***
Aborted (core dumped)
不知道为什么CentOS_OpenCFD在我这儿用intel fortran的openmp总是出问题，gfortran就没问题。
而李老师的程序两个编译器都没问题。
我的是centos6.2

ch06

先测试到这儿，以后有时间在弄。

entropies

你QQ是什么昵称？

ch06

呵呵，你已经拒绝加我为好友了。

[ 本帖最后由 ch06 于 2012-6-15 11:18 编辑 ]

entropies

给我发email吧，seeing_2004@163.com

lixl-imech

原帖由 ch06 于 2012-6-9 19:37 发表

                               
登录/注册后可看大图

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)    ...

OpenCFD-EC采用虚网格(Ghost Cell 或Dumy Cell)来实现边界条件。 在计算域外，扩展一层或两层虚网格（OpenCFD-EC采用两层虚网格）。 对于无滑移固壁边界，虚网格上的速度值与对应内点处速度值符号相反。
这样，计算出的边界处的无粘通量刚好为0. 粘性通量也可得到正确的值。

ch06

原帖由 lixl-imech 于 2012-6-18 11:46 发表

                               
登录/注册后可看大图

OpenCFD-EC采用虚网格(Ghost Cell 或Dumy Cell)来实现边界条件。 在计算域外，扩展一层或两层虚网格（OpenCFD-EC采用两层虚网格）。 对于无滑移固壁边界，虚网格上的速度值与对应内点处速度值符号相反。
这样， ...

多谢李老师指点。

cj1206

李老师：您好！
    在您的程序（OpenCFD-EC2D v1.4.3）的sub_LU_SGS.f90文件中，为计算DF调用如下函数：
   comput_DFn(NVAR1,DF,B%U(:,i-1,j),B%DU(:,i-1,j),B%ni1(i,j),B%ni2(i,j),gamma)
       请问其中参数B%DU(:,i-1,j)代表什么含义？
   另外这里B为 Block_TYPE类型，但是我在这个结构体的定义中没有找到DU这样一个成员，还请您指点。
    谢谢。

lixl-imech

原帖由 cj1206 于 2012-6-22 15:57 发表

                               
登录/注册后可看大图

李老师：您好！
    在您的程序（OpenCFD-EC2D v1.4.3）的sub_LU_SGS.f90文件中，为计算DF调用如下函数：
   comput_DFn(NVAR1,DF,B%U(:,i-1,j),B%DU(:,i-1,j),B%ni1(i,j),B%ni2(i,j),gamma)
       请问其中参数 ...

Block_type 里面包含了du(:,:,: )这个指针成员。 在主程序opencfd_ec2d_v1.4.3.f90第66行有定义。
du（:,:,: ）是LU-SGS中使用的变量， du=U-Un, 即两个时间步物理量的差。

entropies

李老师，你这里有没有k-e turbulence model 或者 k-omega 的？ 我仿照 SST 编了一个 k-omega, 总是有错 。。。

cj1206

李老师：您好！

    感谢您的回复，我最近在学习您的代码，关于边界条件处理上，我有一个问题希望向您请教。
    例如，计算Naca0012亚音速(Ma=0.8)机翼绕流问题的远场边界的设定问题，计算区域大小为机翼弦长20倍左右，Euler方程。

    方法一：只区分入口、出口，但不区分亚音速，超音速，即入口边界的虚拟单元使用远场值，出口边界的虚拟单元使用内部外推。

     方法二：区分入口、出口，也区分亚音速，超音速，完全按照黎曼不变量进行处理，即远场使用特征边界条件。
      
     请问：方法一的处理方式相对于方法二会影响计算定常解的收敛性吗？

lingjie111111

李老师，解N-S方程的三维程序今年和明年这两年的时间里能公布吗？若不能公布，那估计什么时候能公布？

[ 本帖最后由 lingjie111111 于 2012-7-28 13:28 编辑 ]

merryfish

when your 3 dimensional code will be available?