图1最大的作用也许是帮助你熟悉SIMPLE_C的框架。方程求解的细节都被包装在了setup2( )子程序中了。setup2( )子程序的流程见图2。
图2 setup2( )子程序的流程
有一个细节需要注意。图2中u、v、压力校正以及其它标量方程可能求解多次,但是无需求解至收敛。这是因为对于一个非线性的问题,在一组固定的系数下把代数方程的求解一直进行到收敛是不必要也是不明智的。在计算出一组系数之后,我们必须进行充分的迭代以求解方程,以便从这些系数上得到足够的好处。但是把时间过多地花费在这些只是以暂时性的系数为基础的方程是不明智的。求解次数在SIMPLE_C在ntimes[]中设定。这一部分更详细的讨论可以参阅《传热与流体流动的数值计算》4.4-3节。