simple算法的疑惑

Aurora827

我看书上说源项b=0可以作为收敛条件，我编程的时候用while(err<abs(b(j,i))结果一循环，所得结果是NAN，不知道怎么回事，希望高手指点。在编程的时候，这个收敛判断如何编写？在求压力修正的时候是不是还要再假设一下P‘的压力分布，然后迭代算出修正压力？

周华

是不是数值格式设置问题。SIMPLE方法可以参考versteeg和malalasekera合著的那本书，或者看Anderson的书也可以。