@zju_klh
尝试在controlDict中增加一行看看,添加lib,其中XXX替换为编译模型库的名字
libs ("libXXX.so");
@zju_klh
尝试在controlDict中增加一行看看,添加lib,其中XXX替换为编译模型库的名字
libs ("libXXX.so");
@东方白杨 这是输入 uvw三维湍流场 画出的 球面积分能谱 吗?
@zju_klh options文件最后回车加个空行试试
建议看下书了解几种谱的区别,《Tennekes, H. and J. Lumley, A first course in turbulence. 1972: MIT Press.》,这本书也有对应中文版《湍流初级教程》。
简单说区别是:
@东方白杨 不清楚什么原因,你看在程序最后再重新计算一遍保存,看看结果有无区别
GradCenPot = - fvc::grad(CenPot);
GradCenPot.write();
@Voynich 目前看到的处理都是把k忽略掉,这个帖子OpenFOAM 不可压缩湍流模型的 divDevReff 函数 | Giskard's CFD Learning Tricks中的2个观点,其实本质都是k被忽略掉,只是不知清楚k忽略对结果影响具体是多大,我也没具体搜过有无人做过这方面的比较
@coolhhh 如果还有问题,就用最笨的方法,用probes
,人工生成监测点位置,这样结果就没问题。但probes没法监测壁面边界的数据,比如wallShearStress
,见帖子probes和sample提取wallShearStress问题
@lizhisongsjtu
我之前是用of2206, axis
用distance
试下。另外网格划分有1000个吗,1000个点是否会很密了?
x_0_1
{
type uniform;
axis distance;
start (0.00010000 0.00000100 -1.72500000);
end (0.00010000 1.99999900 -1.72500000);
nPoints 100;
}
@lb0331 参考这个帖子如何控制使用pisoFoam计算时的库朗数?,计算时候严格让CFL完全小于1是很难的,可以适当大一点,只要最后结果计算是合理的。
@郑学习 好的,感谢回复
@郑学习 现在我也遇到了buff/cache过大导致计算变慢问题。想问下现在把swap拉高到256,后面还会出现buff/cache过大问题吗?
@18835721626
试试把PIMPLE
中参数改为:momentumPredictor yes;
@东方白杨 感谢回复,学习了
@东方白杨 如果存在window=20,那么21s之后按照前20秒求平均。这句话是否是指21s输出的fieldAverage,是1~21s的平均结果?还是指21s输出的还依然是0~20s的平均结果吗?
@Afr1yne 如果模型什么都没变,只是网格加密, u k e正常也是一样的。把fvSchemes换成最稳定的格式计算看看,再逐个改二阶格式试算,看哪个设置导致的发散
@Afr1yne 如果普通的均匀网格都发散,有点难搞,得再仔细查一查什么原因。fvSchemes和fvSolution改为更稳定的格式试试
基本不用修改,发散大概率是网格问题。看看fvSchemes
和fvSolution
怎么设置的,有的格式在好点的网格能算,复杂的网格就发散
对于不可压缩湍流,of输出的wallShearStress还要乘以空气密度。另外是否可能实验的数据无量纲化处理了?找下类似的研究,别人的文献数量级是多少,判断哪个计算错了
我没具体操作过,正常画网格时候就要像你那个图下面挖空一样,额外挖空一个小长方形
@Afr1yne 设置边界要掏空一个区域,源项参考这个帖子:关于OpenFOAM中的fvOptions的作用