中心差分格式

以一维、稳态、无源项的对流-扩散问题作为讨论对象，

 

 

 

可根据图示建立相应的控制方程如下：

 

为获得对流扩散问题的离散方程，必须对

 

因此，F和D在控制体积界面上的值分别为：

 

在此基础上，定义一维单元的Peclet（其物理意义为对流速率与扩散速率之比）数如下：

 

Pe表示对流与扩散的强度之比。因此，当Pe数为0时，对流-扩散问题变为纯扩散问题，物理含义即为：流场中没有流动，只有扩散。当Pe>0，时，流体沿x正方向流动，Pe<0时，情况正好相反，而Pe值很大时，演变为对流问题，扩散作用可以忽略。

此外，再引入两条假定：

（1）在控制体积界面e和w处的界面面积存在如下关系：

 

（2）方程右端的扩散项，用中心差分格式表示。

所以上述对流-扩散方程可以写为下列形式：

 

 

 

将上述方程整理后可得：

 

 

将上式中的各项前的系数分别用一下形式表示为：

 

可以证明，当Pe<2时，中心差分格式的计算结果与精确解基本吻合，但当Pe>2时，中心差分格式所得的解就不具有物理意义，这是由于当Pe>2时，系数aE<0造成的，我们从公式中可以得到系数aE和aw代表相邻两点处的物理量通过对流扩散作用对P点产生影响的大小，负的系数会导致物理上不真实的解（个人形象理解是，当上面系数小于0时，随着时间推进，方程左边的项很可能在某个时间点突然变为负值，明显物理世界，该值没有意义）。

因此为保证Pe<2，根据：

 

对于给定的ρ与Г，必须保证网格数尽可能的小或者速度u很小，因此基于该限制，必须创建其他更合适的离散格式。

参考文献：

【1】计算流体动力学分析：CFD软件原理与应用.