fvMatrix 的系数正负求解稳定性
-
最近修改方程发现
-fvc::div(phiHbyA)+fvm::laplacian(rAUf, p_rgh)
fvc::div(phiHbyA)-fvm::laplacian(rAUf, p_rgh)
求解的稳定性完全不一样前者稳定，后者发散了，觉得很奇怪，大佬们有何高见@东岳 -
这俩个不是一个加一个减么
-
@东岳 对的，按照理论来讲只是方程的系数变了一个符号，相当于
solve(pEqn==0) 和solve（-pEqn==0) -
后面UEqn也跟着处理了么
-
@东岳 Ueqn没有处理，是通量的符号出了问题吗？
-
@yfclark pEqn.flux()跟里面的写法有关系，如果你写反了，那组建U的时候，也要反一下
-
@东岳 ok，明白了，谢谢东岳老师
