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

entropies

这个file里 init_flow_read 感觉只考虑了 LAP=1 的情况？

CentOS_OpenCFD 的存储方式是 node based 和 cell based，没有考虑 init_flow_read 读取 buffer cells 和 ghost cells 的情况， sorry for that.  但我还是认为我的存储更好。

entropies

有些情况下，比如 internal flow，知道某部分是层流，某部分是湍流，是不是可以考虑利用 turbulence 的 block属性，强制某部分是层流，而只在某些部分解湍流？

shirazbj

FFLAGS  = -g -O0 -check bounds
-check bounds:     bounds overflow
这个check bounds gfortran支持么？

图中的k怎么声明成real的了。

[ 本帖最后由 shirazbj 于 2012-5-27 10:35 编辑 ]

shirazbj

回复 184# entropies 的帖子

还有这个

entropies

gfortran 有自己的溢出检查flag, 是 -fbounds-check，其他的flags可以查阅 manpage;

另如果你用的是 gfortran, 那应该用的是 linux, intel 的 fortran compiler 也可以免费装啊。

internal boundary 是从 1.4.1 里的 farfield1 改过来的， k 应该声明成 integer， good catch!

我想 step case 还有很多地方要改，现在还没有测完，它需要用到 pressure inlet 和 outlet boundary，需要写相应的 bc 。

186＃的问题肯定是我IO没有写好，我再改改看，让gfortran编译后出来正确的显示。


PS： 这周末我估上去会把精力放在改nondimensional为 dimensional, 为我自己的研究用。

shirazbj

我是在win xp下用gcc 4.5.1里的gfortran.

I1.2改成I4.2就能算了。

输出我用gmsh, 得自己写。

outlet boundary我也想知道怎么写，想把流量算平衡了。

entropies

J. Blazek 的书上有 其实 ver 1.4.1 里的 farfield1 也写了 outflow的 处理 但是我不喜欢 把 inlet 和 outlet 写在一个 farfield1里

cj1206

李老师：您好。
    我在单独编译sub_Finite_Difference.f90编译器会报错，我查看了一下错误的位置和原因，大体是由于B被定义为Block_TYPE的结构体，但Block_TYPE内并没有诸如B%Akx(1,1),B%Aky(1,1),B%Aix(1,1),B%Aiy(1,1)的成员变量，不知是何原因，还望指教？

    疑似错误的位置：sub_Finite_Difference.f90 34行，357行。

lixl-imech

原帖由 cj1206 于 2012-5-27 16:23 发表

                               
登录/注册后可看大图

李老师：您好。
    我在单独编译sub_Finite_Difference.f90编译器会报错，我查看了一下错误的位置和原因，大体是由于B被定义为Block_TYPE的结构体，但Block_TYPE内并没有诸如B%Akx(1,1),B%Aky(1,1),B%Aix(1,1),B%A ...

      这段代码是以前版本的，还没有调试通过。 这是有限体积程序中使用有限差分模块作为“补丁”的一段代码。即将发布的1.5.0版就可以使用了。

lixl-imech

发布OpenCFD-EC 2D ver 1.5.0,  见主贴。
该版本支持有限差分-有限体积混合方法。 该方法是作者提出来的一种新方法，可能会存在一些问题。可作为方法探索使用。

shookware

你可以用标准长度（米）做无量纲单位，就可以计算有量纲问题了。

shookware

计算温度越界了，通常出现这种报错都是由于计算发散导致的（个人经验）

shookware

RK算法竟然和LUSGS算法求定常流的结果产生一定的差异，原因是什么？看着好奇怪。按照道理讲定常问题的解是唯一的啊（对空间格式算法相同的前提下）

entropies

3rd RK 是不是更精确些？

ch06

感谢李老师！