CFD中文网

@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.》，这本书也有对应中文版《湍流初级教程》。

简单说区别是：

1. 时间谱（频谱）：对时间相关函数做傅里叶变换，任意位置(x,y,z)处的u、v、w时间序列都可以计算对应的时间谱
2. 一维空间谱（波谱）：一维空间相关函数做傅里叶变换。比如u沿着x方向可以计算空间相关函数$R_{uu}^x(r)$，对其傅里叶变换得到u沿着x方向的一维空间谱$S_{uu}^x(k)$。类似的，uvw与xyz组合，一维空间谱总共有9个
3. 三维能谱：
   (1) 三维空间速度u(x,y,z)可以计算三维空间相关函数$R_{uu}(r_1,r_2,r_3)$，对其三维傅里叶变换得到三维空间谱$\Phi_{uu}(k_1,k_2,k_3)$。类似的还可求$\Phi_{vv}(k_1,k_2,k_3)$,$\Phi_{ww}(k_1,k_2,k_3)$
   (2)记$\Phi(k_1,k_2,k_3)=\Phi_{uu}(k_1,k_2,k_3)+\Phi_{vv}(k_1,k_2,k_3)+\Phi_{ww}(k_1,k_2,k_3)$，三维能谱就等于$\Phi(k_1,k_2,k_3)$沿半径为$\kappa$球面积分再乘以1/2系数，记作$E(\kappa)$，具体看下面书中的介绍：
   
   

@东方白杨 不清楚什么原因，你看在程序最后再重新计算一遍保存，看看结果有无区别

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的作用