rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
-
我把我的求解文件放到这里,大家有兴趣的可以帮忙瞧一下。
链接:https://pan.baidu.com/s/15zFl76UsfD_XyWwUYEdyGQ
提取码:o3ae -
最近将离散格式全部改为中心差分格式(原先全部是二阶迎风),计算的结果还比较让人满意,但是这个差分格式还是不太让人满意,并不稳定,对于这个算例最终可以收敛,收敛的结果如下图所示。但是只要换一个叶片马上又发散了,现在真的是很疑惑,前缘的这类现象究竟要怎么调,或者需要动求解器内部的结构了吗?
divSchemes { default none; div(phi,U) Gauss limitedLinearV 1;//Gauss linearUpwind limited; div(phi,k) Gauss limitedLinear 1; div(phi,omega) Gauss limitedLinear 1; div(phi,e) Gauss limitedLinear 1; div(phi,h) Gauss limitedLinear 1; div(phi,K) Gauss limitedLinear 1; div(phi,Ekp) Gauss limitedLinear 1; div(phiv,p) Gauss limitedLinear 1;//Gauss linearUpwind limited; div(phid,p) Gauss limitedLinear 1;//Gauss linearUpwind limited; div((phi|interpolate(rho)),p) Gauss limitedLinear 1;//Gauss linearUpwind limited; div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; }
-
@东岳 首先非常感谢东岳老师
就我目前能找到的试验数据,Naca 65-810,fluent 计算和试验值相比较准一些(以下是该叶型计算值与试验值的对比),但是OpenFOAM 的计算结果在趋势上就差很多,所以个人倾向于fluent的计算结果。下面给出NACA65-810的计算网格以及fluent的计算结果与试验值的对比。
NACA65-810网格(y+~30-60) 链接:https://pan.baidu.com/s/1tZPl2NfspkOqjmoV4rSpvQ
提取码:jhg4 -
最后一次自己顶一下这个帖子
最近尝试了很多,从差分格式到求解和算法控制,一直都没有得到合理的结果,应该是求解器本身的问题了。但是自己时间限制,马上研三了,再做下去心里实在没底,准备去用自己ANSYS 那一系列的软件做事情去了。
但是还是有些不甘心,毕竟了尝试了那么久。我的求解文件和一些网格都在上面,如果有人知道我的错误在哪里,请一定要告诉我,不论时间过去多久,我都愿意拿自己现在一个月的助研费作为酬劳(400,不多哈)。
最后还是感谢东岳老师一直以来的帮助。
-
-
你需要把OpenFOAM的数据、Fluent数据、实验数据放在一起看区别
-
对不上的时候需要找原因。如何找对不上的原因需要从数据上分析,Cp对不上可能别的能对上。还需要看别的能不能对上,比如壁面剪切力,轴线湍流动能、进口速度等重要的变量。
我近期要备课OKS,只能在7月后期抽出时间给你debug。你这个要写SCI么?还是就放在硕士论文中?我觉得你的研究很系统,如果你做出来,我愿意赞助 http://www.cfd-china.com/topic/1716/cfd验证算例支持计划
-
15/53