时隔九个月,今天重新看这个帖子,终于搞清楚了,记录一下。 确实是像李老师说的那样,Usource()就是源项,这里是手动加的源项,之前一直以为是系数矩阵的非对角元。 Udiag()代表U方程的对角线元素,来自IduMatrix,而Usource()是源项矩阵,来自fvMatrix。上面的代码之所以会出现Udiag(),其实就是像李老师说的那样,是为了增加对角占优。而多孔隙模型的源项就是从Usource加入的: Usource[celli] -= V[celli]*((Cd - I*isoCd) & U[celli]); 这里Cd就是公式(2)中括号中的部分。