中科大计算热物理讲义05对流扩散方程的数值方法.docxVIP

97

第五章对流扩散方程的数值方法

前章所讲的导热问题，主要针对介质没有宏观迁移的导热体，所解的只是通用微分方程组中简化了的能量方程——扩散方程。虽然也涉及有介质宏观流动的充分发展的管流问题，但最终要解的方程还是归结到扩散方程，它只不过是扩散方程解法的一种应用。解扩散方程，空间导数项均为二阶，离散方法相对简单，用中心差分就可以了。本章要讲的对流扩散方程，则是针对热物理中介质的宏观流动和扩散相互耦合的一类更为普遍的问题，必须联立求解连续、动量和能量方程。由于对流项的引入，要比单纯求解扩散方程复杂得多。非线性的对流项和作为动量方程源项的压力梯度项都是一阶导数项，它们所包含的物理意义不一样，数值求解这类问题的一些主要困难都出自对它们的离散。为了循序渐进，条理清晰，本章先不追究流场速度如何得到，而假定速度已知，讨论通用方程中求解函数的对流项如何离散，解决在引入对流项后，一般的对流扩散方程的求解途径问题。下一章进而讲述如何求解流场，并解决压力和速度的耦合问题。

5.1合适选择对流项离散格式的重要性[1,2,3,4,5]

对流扩散方程中的对流项，代表的是流动介质的宏观定向迁移，具有明确的物理意义。合适的对流项离散格式理应能充分体现该项的物理本质。要做到这点，对于一阶导数形式的对流项，迎风型离散格式应是一种好的选择。但一阶迎风格式，与扩散项通常采用的二阶精度的中心差分离散格式不匹配，从而使整个对流扩散方程的离散精度降低。如果对流项采用中心差分格式离散，在精度上能达到与扩散项匹配的要求，但中心差分格式不能体现对流项的物理本质，常会引起数值解的振荡，导致因为对流项离散格式不当的不稳定性。为了要解决前两者间的矛盾，近年来发展了一些高阶格式，精度不低于二阶而又不产生数值振荡，但由于格式构造复杂，所形成的离散代数方程的求解又相对困难，计算机时消耗较多，这又涉及到数值解的经济性问题。

准确性、稳定性和经济性成了选择对流项离散格式形式需要综合考虑的因素。下面先以最简单的一维稳态无源项对流扩散方程为例，介绍几种常用离散格式构成及其数值解对精确解的偏离程度，分析这些格式的特性和使用条件，从而为讨论多维一般对流扩散方程的离散方法准备条件，并为改进对流项离散格式的设计打下基础。

5.2一维稳态对流扩散问题[6]

为能与精确解作比较，讨论一种最简单的一维稳态无源项对流扩散问题，其模型方程的守恒形式为

(5-1a)

5.2.1模型方程的精确解

设ρ,u,Γ为常数，方程(5-1a)在下列边界条件下

x=0,φ=φ0；x=L,φ=φL（5-1b）

图5-1模型方程(5-1)的精确解

98

具有如下形式的精确解析解

其中Peclet数Pe=ρuLΓ。图5－1示出了在不同Pe数下，φ随x变化的曲线。从图易见，Pe=0时，φ随x成线性变化，相应于无流动的纯扩散问题。当Pe的绝对值从零增大时，φ随x逐渐偏离线性变化，在Pe1时，偏离线性不大；直到Pe=5的中等数值，φ随x在整个解域仍然变化比较平稳。但随着Pe绝对值进一步增大到Pe10的数值时，在0≤x≤L的大部分范围内，上游值φ0(对Pe0)或φL(对Pe0)占了优势，仅在靠近外（或内）边界L（或0）的薄层内，φ值才由φ0(或φL)迅速上升（或下降）至φL(或φ0)，这是一种类似于边界层类型问题的特性。精确解的变化特性为设计某些离散格式和分析格式性能提供了比较依据。

5.2.2中心差分格式

1.控制容积积分离散

按图3-2所示的一维网格系统，将方程(5-1)对控制体P自w→e积分，有

设求解函数φ在对流项和扩散项中都为分段线性空间分布，则有

将界面流量ρu记为F,界面单位面积上扩散系数Γ/δx记为D，移项整理上式，得aPφP=aEφE+aWφW(5-5)

其中

(5-6a)对流扩散方程，连续方程是必须要满足的。对于本例，连续方程满足时，必有Fe=Fw，因此离散格式系数aP等于邻点系数之和。

为简化格式写法并便于分析格式性质，定义一个以网格步长为特征尺度的