基于Peclet数判别法的一维对流扩散方程分类研究.docxVIP

基于Peclet数判别法的一维对流扩散方程分类研究

摘要：采用Peclet数的绝对值大小来判别一维对流扩散方程为对流占优型或是扩散占优型方程，运用三种隐式差分格式—中心隐式格式、对流C-N型格式和扩散C-N格式，对不同Peclet数的算例进行离散和求解。然后，将计算区域中所有节点的解析解与数值解表示成矩阵形式，并求解出它们的矩阵2范数之后作比较，两者越接近则代表差分格式精度越高。通过比较得出了当方程Peclet数的绝对值小于等于0.5时，方程为扩散占优型方程。在离散方法选取方面，针对扩散项的离散可以采用更高精度的差分格式，如扩散C-N格式；当Peclet数的绝对值大于等于20时，方程为对流占优型方程。此时，针对对流项可以采用更高精度的差分格式，如对流C-N格式；当Peclet数的绝对值介于0.5与20之间时，无法用Peclet数判断方程类型,不过可以选择折衷的离散格式减小误差，如中心隐式格式。

关键字：一维对流扩散方程Peclet数判别法有限差分方法数值模拟MR(2010)主题分类号：39A14；65M06 中图分类号：O242.2

文献标识码：A

1.引言

一维对流扩散方程是描述流体流动和传热问题的一类线性化模型方程。 土壤、大气等多

孔介质中水、盐分、温度以及污染物质的对流扩散问题都会遇到此类方程。 在一维对流扩散

方程的求解过程中，反映流体对流和扩散两种物理作用的分别是对流项和扩散项。所以，根据方程中对流项还是扩散项占主导作用，通常可将方程分为对流占优型和扩散占优型两类方程。然而，要想得到精确度较高的数值结果，这两种类型方程的离散方法不能采用相同的离散格式。因此，需要有一种判别方法来判断方程的类型，关于对流占优型和扩散占优型方程的判别方法一直是近年来研究的热点问题。这对研究不同类型的方程使用合适的差分格式进行离散具有实际的意义。由于Peclet数的绝对值表示了对流作用相对扩散作用的大小，即绝

大，扩散所起的作用就可以忽略。反之，当Peclet数为零时，方程就为纯扩散方程。本文选用一维定解非稳态对流扩散方程为例，通过考察Peclet数的绝对值大小来对方程进行分类，

方程一般形式如下：

uau2uf(x,t),x1xx2,t0

tx x212

u(x,0)g(x)(1)u(x1,t)(t),u(x2,t)(t)

uu(x,t)

其中a和分别代表对流项系数和扩散项系数。假定求解区间长度为s，Peclet数的绝对值

计算公式为：

计算公式为：

(2)Peas

(2)

从公式（2）中可以看出，当计算区间长度给定， Peclet数是由对流和扩散系数确定的[1]。

面介绍方程（1）的离散方法。

2.离散方法

2.1显式格式离散

对于上述方程（1），需要离散非稳态项（简称U项）、对流项（简称C项）和扩散项（简称D项）。常见的离散方法有显式格式和隐式格式两种。显式格式有：中心显式格式、修正中心显式格式、迎风差分格式等。比如，以中心显式格式为例，即使用向前差分格式、一阶中心差分格式与二阶中心差分格式组合分别离散U项、C项和D项。其离散形式如下：

(3)n1nnnnnn

ujuj uj1uj1 uj12ujuj1

(3)

2ha2

h2

2h

其截断误差为o（h2）。然而，由vonNeumann判别条件判断此种格式将受到稳定性条件的限制，即：2,1。相应其它显式格式同样有稳定性条件限制。所以，显示格

a2,h22

式时间步长和空间步长h取值将受到限制。因此，若采用显式格式求解一维非稳态对流扩散方程问题，得到的数值解精度将受到限制，甚至误差很大。所以，显式格式的离散效果欠

佳，为了弥补它的缺陷，尝试采用无条件稳定的隐式格式离散（ 1）式。

2.2隐式格式离散

常见的隐式格式有三种：向后差分格式、一阶中心差分与二阶中心差分组合（简称中心隐式格式）；向后差分格式、C-N型一阶中心差分与二阶中心差分组合（简称对流C-N型格

式）；向后差分格式、一阶中心差分与C-N型二阶中心差分组合（简称扩散

式）；向后差分格式、一阶中心差分与

C-N型二阶中心差分组合（简称扩散C-N

型格式）。

三种隐式格式的离散形式如下：

1）中心隐式格式：

n1nun

n1nunj1ujn

n1n1

uj1 uj1

a(j1j1)

2h

n1n1

unj112unj1