紧致差分的matlab程序,对流占优扩散方程的一种新C—N 紧致差分格式

冯立伟+席伟

摘要：给出了对流扩散方程一种新的C-N紧致差分格式，其截断误差为时间二阶空间四阶，且为无条件稳定的。编制了MATLAB程序，数值试验表明了格式的有效性。

关键词：对流扩散方程；紧致格式；C-N格式

中图分类号：TP311 文献标识码：A 文章编号：1009-3044(2015)30-0173-04

A New C-N Compact Difference Scheme for Solving Convection-Diffusion Equation

FENG Li-wei1，XI Wei2

(1. ShenYang University of Chemical Technology ， Shenyang 110142， China； 2. ShenYang University of Chemical Technology ， Shenyang 110142， China)

Abstract： It gives a new C-N compact finite difference scheme for solving Convection-diffusion equation. Its error is two order in time and four order in space， and is unconditionally stable， and writes MATLAB programs. Finally， a numerical experiment shows the effectiveness of the scheme.

Key words：convection-diffusion equation； compact scheme；C-N scheme；

对流扩散方程描述了在自然界中大量出现的对流和扩散现象，在流体力学、环境科学以及能源开发等诸多领域有着广泛的应用。因此研究对流扩散问题的数值计算方法就尤为重要[1-3]。田振夫、葛永斌等[4-6]使用Hennite 插值思路给出了求解对流扩散方程的空间四阶差分格式。杨志峰和陈国谦等[7-9]使用综合变换建立了求解对流扩散方程的一种两层隐式四阶紧致差分格式。肖建英等[10]通过引入指数变换造了一种高精度的紧致隐式差分格式，本文采用指数变换将含源项的对流扩散方程变为纯扩散方程，扩散项使用pade逼近的紧致差分离散，时间层上采用C-N格式，得到了