postProcess 求解的div不为0

[vbcwl]

请教一下各位前辈，我做的槽道流，体积平均速度为15.43，然后使用postProcess -func 'div(U)'求出来散度，值都在0.03 0.04这附近，甚至能到0.13，求解格式是div(U) Gauss linear，用的不可压缩pimpleFoam，这个散度这么大是什么情况？

[vbcwl]

@vbcwl 我后面用grad求解出du/dx dv/dy dw/dz这些，相加确实等于这个div求解结果，但是du/dx 与(u2-u1)/(x2-x1)相差巨大，有好几倍，就算这个手算精度差也不至于这么大差距吧，请问有前辈了解么

[李东岳]

你要是验证$\nabla\cdot\bfU=0$，最好验证phi是否守恒，而不是U。OpenFOAM这个保证phi是严格守恒的。你试试？

[vbcwl]

@李东岳 明白了，谢谢东岳老师