3.3 速度u差分方程系数的产生
按照(传热5.57与5.58),系数的先要计算控制体积四个面上的通量Fe,Fw,Fn,Fs与扩散强度
De,Dw,Dn,Ds。程序中以flow与diff来存储这些量。
在得到通量flow与扩散强度diff之后,各个系数都可以通过一个统一的计算过程计算得到。
通量的计算涉及到控制体四个控制面上的速度与密度。而由于交错网格的存在,它们的信息不总是存储在同一几何面的中心位置上,因而需要分段或者插职计算。这是阅读setup2( )子程序一个主要的困难来源。
下面我们以u的离散方程介绍一下每一个u(i,j)控制体的e,w,s,n四个面上的flow与diff的计算。
表3 u(i,j)控制体的e,w,s,n四个面上的flow与diff的计算
在具体的介绍之前,我们先了解u方程系数设定的流程(图8)。首先,记住
aim(i,j)=aW(i,j);
aip(i,j)=aE(i,j); ajm(i,j)=aS(i,j);
ajp(i,j)=aN(i,j)