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,明白了,谢谢东岳老师