0_1464584218751_upload-39a9dbfb-d412-4c58-9b21-9a957ca37b6e 据说是因为solid和porous两个域没有接起来，求解决

每个网格内的压力不一样，怎样使源项在网格中不同，而不是用整体的源项呢？

@cfd-china 网格应该没太大问题，主要是设置这块怎么检查，0_1464338494250_upload-1e7bd2e0-4b82-4053-adad-bad9dc0e6fdd 有人说是solid和porous这两个区域没有接起来，但是我不知道怎么接

我并没有查看相关的代码，不顾我猜测是在求解pEqn的时候读取fvSolution里面的pimple设置的最后迭代的p求解器来求解.

谢谢各位的回复！

2 sp su 怎样赋值给每个网格，而不是整体的源项？

@buxiangle2012 我看了看simple.loop代码，表示没看懂，不知道怎么和fvsolution设置残差联系起来的。。

我感觉系统不同，时间也会多少有点不同吧 就算同样的配置。 虚拟机win系统，ubuntu直接跑 感觉ubuntu跑得快。虚拟机毕竟你要分配内存给他，而你用ubuntu直接跑直接调用系统内存。所以我觉得直接跑省时间。

@cfd-china 这个东西之前是写给DPMFoam的。。然后我搬到我这边了。。就出现这样了。

@cfd-china 非常感谢您的指导，但依旧希望您有空的时候帮忙设置一下我的那个算例，对我们这些新手的参考作用确实非常大

@李东岳 东岳老师我想说一下自己的理解，不知道对不对：最后速度更新的时候，U = HbyA - rAUfvc::grad(p)中HbyA是通过HbyA = rAUUEqn.H();//公式(15)得出的吧，也就是Ur,预测的速度而不是n+1步的速度对不？我还没看SIMPLE，不过按照我的理解是不是HbyA要随着速度变化不断更新？那内迭代的最终截止条件是什么呢？非常感谢！

@cfd-china 我在文件里找到了#include "compressible/turbulentModel/turbulenceModel.H这行代码，修改成#include "turbulenceModel.H这样后出现了差不多的报错0_1464072322810_upload-7effb60a-4303-444d-8e08-203e98037205

最近找了好多方法来求浓度，都不行，但是controlDic下面的function功能好像可以实现这个功能，有谁有相关经验么？@散漫守望2016

@sejabs 主要是就找到这个例子是比较接近我要做的模型，其他的尚未找到

@cfd-china 我试试，谢谢

@wwzhao 好没关系 谢谢您

怎样在主程序中添加correct( ); 来在更新类中的　rho1()rho2() p()

暂时还没研究过这个类:sad: ....

@gzh3369

如果公式是这样的话：$x=\frac{\sum V_i x_i}{\sum V_i}$。可以这样：

forAll(x, celli)//对所有的体x进行计算 { x[celli] = x[celli] *mesh.V()[celli];//此处即为每个网格的x*网格体积 } scalar sumx = gSum(x);//对所有的x*网格体积进行加和 scalar V = gSum(mesh.V());//计算总网格体积 scalar Final = sumx/V;//计算体均湍流动能耗散率 Info << runTime.timeName() <<" volume average x is "<<x<<endl;//输出结果

楼主好人 w-y10@foxmail.com