摄动差分法对有源对流扩散方程计算精度的提高.pdfVIP

曾方军等：摄动差分法对有源对流扩散方程计算精度的提高  摄动差分法对有源对流扩散方程 计算精度的提高 1,2 1 1 2 曾方军 ，舒适 ，李明军 ，高智 （1. 湘潭大学，数学与计算科学学院，湖南湘潭，411105） （2. LHD，中国科学院力学研究所，北京海淀，100190） 摘要 对流扩散方程中函数源项的存在能使许多格式产生假扩散，许多格式计算有源对流扩散方程所 产生的数值误差比计算齐次对流扩散方程的误差要大很多。本文研讨高智数值摄动差分法能否减轻 上述假扩散的问题；通过对流项和源项同时摄动，摄动后格式的整体精度与未摄动格式相比有了很 大的提高，这是摄动格式较其他格式的一条优良特性；本文最后通过数值实验证实了数值摄动差分 法的这条优良特性。由于具有工程价值的对流扩散型方程（例如 Navier-stokes 方程组）都是有源的， 因此将数值摄动差分法的这条优良特性应用到工程计算中去将会产生重要的实用价值。 关键词 数值摄动差分法, 计算流体力学, 对流扩散方程, 函数源项 引言 于在格式中加入了反扩散，易生假振，也存在 数值误差。目前人们对函数源项的离散处理大 对流扩散型方程是流体力学中一类基本 都只是对源项进行线性化处理[2,4]，但这种处理 的数学物理方程，如流体运动的NS方程组就是 只能提高格式的稳定性和加快解的收敛速度， 对流扩散型方程，在描述实际问题时基本上所 对源项所产生的假扩散起不到减轻作用。因此 有的方程都含有函数源项。已有对流扩散型方 对有源对流扩散方程研究能够提高计算精度的 程计算格式的研究，对函数源项离散方法的研 算法和格式显得十分必要；由于数值摄动差分 究还不很深入且看法不同，有人认为函数源项 法对精度的提高是整体精度的提高[8-10]，高智 的存在对离散方程的计算精度毫无影响；但是 教授建议利用数值摄动差分法应该能提高有源 已有研究[1-3]表明，函数源项的存在能够在计 对流扩散方程的计算精度，本文通过研究和数 算过程中引起假扩散,有时甚至还会引起较大的 值计算证实了数值摄动差分法[5-7]确实能对有 数值误差；文献[1]中Leonard写道：“One would 源对流扩散方程的计算精度起到一定的提高作 have to take great care to generalize a 用。 finite-difference algorithm based on partial 1994年，高智教授提出数值摄动差分法[7], first-order upwinding to take account of the 从起初的三节点格式的全域摄动，到现在的三 unsteady effects,cross-grid transport,and explicit