fvMatrix 的系数正负求解稳定性

[yfclark]

最近修改方程发现
-fvc::div(phiHbyA)+fvm::laplacian(rAUf, p_rgh)
fvc::div(phiHbyA)-fvm::laplacian(rAUf, p_rgh)
求解的稳定性完全不一样前者稳定，后者发散了，觉得很奇怪，大佬们有何高见@东岳

[东岳]

这俩个不是一个加一个减么

[yfclark]

@东岳 对的，按照理论来讲只是方程的系数变了一个符号，相当于
solve(pEqn==0) 和solve（-pEqn==0)

[东岳]

后面UEqn也跟着处理了么

[yfclark]

@东岳 Ueqn没有处理，是通量的符号出了问题吗？

[东岳]

@yfclark pEqn.flux()跟里面的写法有关系，如果你写反了，那组建U的时候，也要反一下

[yfclark]

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