第五章 对流-扩散方程的离散格式
第五章——对流-扩散问题的有限体积法
混合格式兼具中心差分格式和迎风差分格式的优 点，具有守恒性、有界性和迁移性，其缺点是按 Taylor级数展开后截断误差为一阶，精度不高
第五章 对流-扩散问题的有限体积法
边界条件处理：
第五章 对流-扩散问题的有限体积法
上机课：高速流计算
第五章 对流-扩散问题的有限体积法
一维稳态对流-扩散问题的有限体积法 举例：考虑一维无源项的稳态对流-扩散问题： （核心区的稳态能量方程） d d d d ( u ) ( ) ( u ) 0 dx dx dx dx
第五章 对流-扩散问题的有限体积法
W
uw
w
P
x
xwp
ue
e
E
d d d ( u ) ( ) dx dx dx
P点中心差分
d ( u ) 0 dx
x pe
xWP
xPE
d d ( uA ) e ( uA ) w (A ) e (A ) w dx dx
设： F u
D x
w e De 有： Fw ( u) w Fe ( u ) e Dw xWP xPE
aPP aWW aEE
第五章 对流-扩散问题的有限体积法
当速度较大时，采用中心差分格式处理边界值， 下游边界条件对数值计算法
离散格式的性质： （1）在数学上，一个离散格式必须要引起很小 的误差才能收敛于精确解，即要求离散格式必 须稳定或网格必须满足稳定性条件。 （2）在物理上，离散格式所计算出的解必须要 具有物理意义，对于得到物理上不真实的解的 离散方程，其数学上精度再高也没有价值
FeE FwP De (E P ) Dw (P W )
[(Dw ) ( De Fw ) ( Fe Fw )]P DwW ( De Fw )E
热流问题数值计算Chapter 5(1)
主讲陶文铨西安交通大学能源与动力工程学院热流中心CFD-NHT-EHT CENTER 2007年11月29日, 西安热流问题数值计算第五章有回流的流动与换热第流场数值计算概述5.1.1两类主要流动与两类数值解法5.1.4两种构造对流项离散格式的方法1.两类主要流动2.两类数值求解方法5.1 流场数值计算概述5.1.2强制对流的涡量方程5.1.3一维模型方程5.1.1 两类主要流动与两类数值解法回流型，其基本区别在于是否存在漩涡(vortex)vorticity) 的区别漩涡是一种宏观的流动形态，特点是流体速度发生反转；涡量是粘性流体的基本特性，只要是粘性流体流动中必有涡量。
动力工程中大多为回流型（椭圆型）流动。
本章仅介绍回流型流动的数值解法。
2. 两类数值求解方法数值求解回流型的流动可以大别为原始变量法与涡量流函数法。
原始变量法u,v,p为求解变量，由于不可压缩流体没有关于压力的独立的方程，数值求解时需要做特殊处理；5.1.3一维模型方程为研究离散格式基本特点又不使过程复杂化，5.1.4两种构造对流项离散格式的方法1. Taylor控制容积积分法－给出界面上被求函数的插值方式对同一种格式，如控制容积积分法得可以认为是控制容积内导数积分中5.2.1 中心差分5.2.2 迎风差分5.2.3 混合格式5.2.4 指数格式5.2.5 乘方格式5.2对流扩散方程的离散格式本节中通过将一维模型方程在取分段线性型线，经整理可得：()eexδΓ+−EaWa做如下变化：()e e x δΓ++为保证代数方程迭代求解的收敛性，我们要求计算中质量守恒一定要满足，于是下列两点边值问题：Pe 随当当当得出结果如右。
，4P =100,W φ=5.5.2 一维对流-扩散方程的迎风控制容积法的定义－界面上未知函数永远取上游Patankar教授提出一种专门符号表示FORTRAN 的Max：,X Y，于是有：(),0,0e P e E eu F Fρφφφ=−−类似地有：(),0,0w W w P wu F Fρφφφ=−−3.对流项一阶迎风、扩散项中心差分的离散方程P P E E W Wa a aφφφ=+,0E e ea D F=+−()P E W e wa a a F F=++−,0W w wa D F=+由于0,0E W a a ≥≥因此FUD 总可以得出物理上合理的解（physically plausible solution )，自五十年代提出以来，半个世纪中得到广泛地采用。
对流扩散方程数值方法
4
情况
E100,W50
E50,W100
P 25 P 125
适用条件有限
➢ 当P 2 时，系数变成负数，会产生振荡
➢ 小适，用或条只件能P处理2 流速这比意较味小着的：问网题格步长必须很
5.2.3 一阶迎风格式
➢ 1.积分离散和格式推演
d(u)
dx
ddx
d
dx
(u)e(u)w d d x e d d x w
5.2.4 指数格式
➢ 1．指数格式推演，根据精确解表达式
L 0 0 e e x x p p ( ( u u L x// ) ) 1 1 e x p e ( x P p e (P x e ) / L ) 1 1
对流-扩散总通量密度
J u d
dx
L 0 0 e e x x p p ( ( u u L x// ) ) 1 1 e x p e ( x P p e (P x e ) / L ) 1 1
1Pe 2
,
0
当 P e2 aED e0
同理可得
aWDwPw,
1Pw, 2
0
2．混合格式评述
➢ 当Pe 2 时，与中心差分格式相同 ➢ 当Pe 2 时，为扩散项=零的迎风格式 ➢ 形式简单，兼有中心差分和迎风格式的优点，
较为适用
5.2.6 乘方律格式
➢ 1．乘方律格式提出和构成[Patankar] 四条包络线
a P a E a W (F e F w )
2. 一阶迎风格式评述
➢ 符合对流项物理意义，离散系数恒>零，无数 值振荡。得到广泛的应用。
➢ 迎风思想：二阶迎风、三阶迎风和QUICK格式 都很好地吸取了迎风思想：让信息主要来自流 动上游方向
对流扩散问题有限体积法
流体仿真与应用第八讲二、对流-扩散问题的有限体积法◆中心差分格式（例子）节点增加到20个结果◆离散格式的性质在数学上，一个离散格式必须要引起很小的误差（包括离散误差和舍入误差）才能收敛于精确解，即要求离散格式必须要稳定或网格必须满足稳定性条件。
在物理上，离散格式所计算出的解必须要有物理意义，对于得到物理上不真实的解的离散方程，其数学上精度再高也没有价值。
通常，离散方程的误差都是因离散而引起，当网格步长无限小时，各种误差都会消失。
然而，在实际计算中，考虑到经济性（计算时间和所占的内存）都只能用有限个控制容积进行离散。
因此，格式需要满足一定的物理性质，计算结果才能令人满意。
主要的物理性质包括：守恒性、有界性和迁移性。
◆离散格式的性质——守恒性满足守恒性的离散方程不仅使计算结果与原问题在物理上保持一致，而且还可以使对任意体积（由许多个控制容积构成的计算区域）的计算结果具有对计算区域取单个控制容积上的格式所估计的误差。
◆离散格式的性质——迁移性③当Pe 为有限大小时，对流和扩散同时影响一个节点的上、下游相邻节点。
随着Pe 的增加，下游受的影响逐渐增大，而上游受的影响逐渐变小。
①，即纯扩散，无对流。
②，即纯对流，无扩散。
0=Pe ∞=Pe◆迎风格式迎风格式（Upwind Differencing Scheme ）在确定控制容积界面上的值时就考虑了流动的方向性，其思想为：在控制容积界面上对流项的取上游节点处的值，称之为第二类迎风格式。
中心差分格式的缺点是，它不能识别流动的方向，控制容积界面上的值取相邻上、下游节点的平均值。
当对流作用较强时，这样的处理就与其物理特征（某点的值受上游的影响，而不受下游的影响）不一致了。
φφφ◆迎风格式◆迎风格式在控制容积界面上对流项的取其上游节点处的值EW →φWw φφ=Pe φφ=()()W P w P E e W w P e D D F F φφφφφφ−−−=−()()[]()Ee W w w P w e e w w D F D F F D F D φφφ++=−+++WE →Pw φφ=Ee φφ=()()[]()Ee e W w Pw e e e w F D D F F F D D φφφ−+=−+−+◆迎风格式通用形式WW E E P P a a a φφφ+=()w e E W P F F a a a −++=EW →ww W F D a +=eE D a =W E →w W D a =ee E F D a −=◆迎风格式的特点迎风格式满足守恒性。
第五章对流扩散问题(假扩散 高阶格式差分方程的求解)
第五章 对流扩散问题———假扩散
在P点的控制容积上对上边的 NW 微分方程进行积分，则：
1 6 N 5 w P E NE
u(Te Tw ) v(Tn Ts ) 0
W 2 SW
A 迎风格式
u(Te Tw ) v(Tn Ts ) 0 ue / w 0
n bb a EE n a EE WW WW b
第五章 对流扩散问题———假扩散
A5.3 PDMA算法 对五对角阵，有没有类似三对角阵TDMA那样的直接 求解方法呢？实际上对五对角阵，人们也可以找到相 应的直接求解方法，这个方法就是PDMA算法。下边 以一维为例来介绍这个算法。
第五章 对流扩散问题———假扩散
由此可见，对一维而言，所得到的差分方程不再是可以 用TDMA直接求解的三对角矩阵方程，而变成一个五对 角阵方程。对二维或三维而言，逐线联立求解的方程组 也不再是可以用TDMA直接求解的三对角矩阵方程，也
变成一个五对角阵方程。
那么，针对这样一种五对角方程，通常我们如何来求解 它们呢？
假设经过代入后得到的上三角方程为
i Ai i 2 Bi i 1 Ci
*#
问题的关键就变成为：找出系数Ai, Bi与Ci和系数a, b, c, d
, e及f之间的关系。为此，写出i-2点的上三角方程如下：
i 2 Ai 2i Bi 2i 1 Ci 2
****
(****)*(di+ eiBi-2)，有：
(di ei Bi 2 )i 1 Ai 1 (di ei Bi 2 )i 1 Bi 1 (di ei Bi 2 )i Ci 1 (di ei Bi 2 )
哈尔滨工业大学 计算传热学 第五章 对流-扩散方程的离散格式-2013
aPP aEE aWW
Fe Fw exp( Pw ) aE , aW exp( Pe ) 1 exp( Pw ) 1
(D)
aP aE aW (Fe Fw )
区别就在函数 aE和aW
aE De
Pe aE De exp( Pe ) 1
aE Pe De
该格式计算量比指数小，且指数格式的解差别很小。
§ 5-3
为了在讨论中引入 PE 记
通用表达式
x
i
J*
i+1 i+1/2
x
1 界面i+ 上的值可以用界面两侧节点值表示 2
J * Bi Ai 1 （y）
系数A和B的性质的讨论 (1)当 i i 1 时，扩散量=0， J *完全由对流造成，即
即
aPP aEE aW W
显然不论那种格式，仅仅是 A(| P |) 表达式的区别。
A( P )
A(|P |)
中心 1 0.5 | P | 迎风 1 混合 [| 0,1 0.5 | P | |] 指数 | P | [exp(| P |) 1]
1.0
迎风
指数 乘方
乘方 | 0, (1 0.1| P |)5 |
中心
混合
P
§ 5-4
原始的假扩散概念
关于假扩散的讨论
一维非稳态对流方程（纯对流,没有扩散）
u t x
显示迎风差分格式
in1 in
t
u
in in 1
x
, o(x, t )
将上式在（i，n）点做Taylar级数展开，保留二阶。
上述若对任何成立，必得
B( P ) A( P ) A( P ) B( P )
对流扩散方程解析解
对流扩散方程解析解对流扩散方程（Convection-DiffusionEquation，CDE）是描述物理系统中物质扩散和热对流运动的方程。
它源于20世纪30年代真空磁体理论中发现的电子运动方程，在50年代被普及应用于各种工程、物理学和化学领域，如电子、热传输、水力学等，具有不可缺少的重要意义。
一般来说，对流扩散方程可以被描述为：$$frac{partial y}{partial t}=afrac{partial^2 y}{partial x^2}+bfrac{partial y}{partial x}+cfrac{partial y}{partial y}+d$$其中，a、b、c和d是常数，t和x分别代表时间和物理位置。
若把空间坐标投射到它们的平面上，则可以用更具体的形式表述为： $$frac{partial y}{partial t}=afrac{partial^2 y}{partial x^2}+bfrac{partial y}{partial x}+cfrac{partial y}{partial y}+d+frac{partial y}{partial z}$$其中，z是投射后的空间坐标，a、b、c和d也可以改变以适合不同的实际应用场景。
对于对流扩散方程的解析解，有两种基本方法：一种是用不定积分法；另一种是用微分平面法，也称作渐进分析方法。
从一般的原理上来看，不定积分法是把对流扩散方程拆解成多个简单的可求解的微分方程，然后分别求解它们，最后再综合求得总解。
此外，它还可以运用标准积分法来近似求解，特别有利于解复杂的多变量方程。
而渐进分析（Perturbation Analysis）是把复杂的问题划分成几个渐进步骤，每一步把问题简化为可以近似解决的状态，依此不断迭代，最终求得近似解。
这种技术通常用来求解非线性方程，对于对流扩散方程求解也非常有效，能有效地提高准确度和计算速度。
此外，还有其他一些求解方法，比如拉格朗日法（Lagrange Method）、拉普拉斯正则化（Laplace Regularization）以及偏微分方程的泛函理论方法（Functional Theory of Partial Differential Equations）等。
2.4常用的离散格式
低阶格式的假扩散特性
迎风格式,指数格式,混合格式及乘方格式等 一阶格式应用于实际问题时都可能引起较严 重的假扩散,这在HVAC领域的高大空间流体 流动及传热计算中尤为明显. 因此,为了有效地克服或减轻假扩散所带来的 计算误差,空间导数应当采用二阶或更高阶的 格式(如QUICK格式,二阶迎风差分格式等).
离散格式
假设速度场已知，则为求解离散方程，需计算广义未 知量在边界e和w处的值。
为完成这一任务，必须决定界面物理量如何通过节点 物理量的插值表示。
各种不同的插值方法就构成了不同的离散格式。
中心差分格式
一阶迎风格式
混合格式
指数格式
乘方格式
1
2.4.1术语的约定
对离散格式的讨论以一维稳态对流扩散方程为例，不 涉及瞬态项。
3
Central differencing scheme 中心差分格式
(x) P
P
interpolated value
e E
eE
We determine the value of at the face by linear
interpolation between the cell centered values.就是界 面上的物理量采用线性插值公式来计算。
基于此限制，中心差分格式不能作为对于一般 流动问题的离散格式，需创建其它更合适的格 式（对纯扩散稳态，如热传导是适用的）。
5
对流扩散方程的精确解
6
精确解随Pe数的变化
（Pe＝0纯扩散，Pe增大对流增强）
7
具体算例
（不同计算工况意味着不同Pe数）
8
第一种工况Pe=0.2
尽管网格粗糙，但数值解与精确解非常接近。
对流方程及其解法
对流方程及其解法对流方程是描述流体运动的最基本方程之一，涉及热、动量、物质等的传递现象，对于各种物理问题的研究都具有重要意义。
本文将从对流方程的基本形式和意义出发，探讨其常见解法及相关应用。
一、对流方程的基本形式与意义对流方程是描述流体中质量、热量和动量传递的方程，其基本形式可以写作：$$ \frac{\partial\phi}{\partial t} + (\mathbf{v}\cdot\nabla)\phi =\nabla\cdot(\Gamma\nabla\phi) $$其中，$\phi$为描述流体量的变量，如温度、密度、浓度等；$\mathbf{v}$为流体的流速，$\Gamma$为扩散系数。
对该方程的解析求解较为困难，故通常采用数值方法进行求解。
下面介绍几种常见的数值解法。
二、有限差分法有限差分法是在连续方程的基础上，利用有限差分代替导数，将微分方程变为代数方程组，从而利用计算机求解的方法。
其基本思想是将求解区域划分为有限个网格，对每个网格内的量用差分代替导数，从而得到有限差分方程。
以简单的二维对流扩散为例，其对流方程为：$$ \frac{\partial\phi}{\partial t} + u\frac{\partial\phi}{\partial x} + v\frac{\partial\phi}{\partial y} = \Gamma\frac{\partial^2\phi}{\partial x^2} + \Gamma\frac{\partial^2\phi}{\partial y^2} $$其中，$u$和$v$分别代表$x$和$y$方向的流速。
对该方程进行离散，假设$\phi_{i,j}$为$x=i\Delta x$，$y=j\Delta y$处的$\phi$值，则可以得到：$$ \frac{\phi^{k+1}_{i,j} - \phi^k_{i,j}}{\Delta t} +u\frac{\phi^k_{i+1,j} - \phi^k_{i-1,j}}{2\Delta x} +v\frac{\phi^k_{i,j+1} - \phi^k_{i,j-1}}{2\Delta y} $$$$ = \frac{\Gamma\Delta t}{(\Delta x)^2}(\phi^k_{i+1,j} -2\phi^k_{i,j} + \phi^k_{i-1,j}) + \frac{\Gamma\Delta t}{(\Deltay)^2}(\phi^k_{i,j+1} - 2\phi^k_{i,j} + \phi^k_{i,j-1}) $$其中，$k$为时刻，$\Delta x$和$\Delta y$分别为$x$和$y$方向的网格间距。
第5章-对流-扩散方程的离散格式
uL
0
Pe表示对流与扩散作用 的相对大小。
0
4/59
传热与流体流动的数值计算
二、对流项的中心差分
d d d u 采用控制容积积分法 对方程 dx dx dx e u e w u w P 2 2 x w x e
aE De Fe ,0 , aW Dw Fw ,0
对流项一阶迎风：
aW i 1 aE i 1 P ,0 1 P ,0 P D D
12/59
传热与流体流动的数值计算
A P P
B P A P P A P P ,0 P B P A P P ,0
24/59
传热与流体流动的数值计算
四、aE、aW的通用表达式
* Je B Pe P A Pe E
J d J P D d x x
*
18/59
传热与流体流动的数值计算
一、通量密度及其离散表达式（续）
J*的离散表达式：
J * Bi Ai 1
Behind Ahead 界面后的项 界面前的项 以坐标轴正方向为依据的“前”、“后”。
19/59
传热与流体流动的数值计算
负系数会导致物理上不真实的解。
7/59
传热与流体流动的数值计算
三、对流项的迎风格式
Taylor展开法
d i i 1 , ui 0 dx i x
i 1 i , ui 0 x
控制容积积分法 e界面 ue 0 , P ; ue 0 , E w界面 uw 0 , W ; uw 0 , P
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
d i i1 dx i x
i1 i x
， ui 0 ui 0
对多维问题，用此方法构造的对流 项的离散格式，只有在求解区域内 流速不发生逆向时，所形成的离散 方程才具有守恒性。
2、控制容积积分法定义
规定界面上的未知量恒取上游节点的值
e界面上： ue 0 , p ； ue 0 , E
J * Pi Pi1
B A P
2、对称特性
对于坐标系I，C位于界面之后，而D位 于界面之前，于是：
J * B(P )C A(P )D
对于坐标系II，D位于界面之后，而C
位于界面之前，于是：
由于
J * B(P )D A(P )C
J * J *'
对Je： 0 P ， L E , L xe
对Jw：
Je
Fe [ P
P E ]
exp(Pe ) 1
0 W ， L P , L xw
Jw
Fw[W
W P ]
exp( Pw ) 1
对于控制容积P，代入对Je、 Jw的表 达式整理得：
P [Fe
exp( Pe ) exp( Pe ) 1
Fw
exp(
1 Pw
)
] 1
E
exp(
Fe Pe
)
1
W
Fw exp( Pw ) exp( Pw ) 1
令： 则
； aE
Fe exp(Pe ) 1
aW
Fw exp(Pw ) exp(Pw ) 1
aP aE aW (Fe Fw)
对于控制体P可写出：
J
* e
B(Pe )P
A(Pe )E
J
* w
B(Pw )W
A(Pw )W
将上两式代入通量密度守恒方程，得
Je
Jw
De
J
* e
Dw
J
* w
0
整理得 P{De B(Pe ) Dw A(Pw )} De A(Pe )E DwB(Pw )W (2)
0
0,(1 0.1Pe )5
0,Pe
Pe , Pe 10
见下页表格：
5.3.5 5种三点格式系数计算式的汇总
不同格式离散方程的形式相同，但 系数不同。具体见下表5－1：
5.4 对流-扩散方程5种3点格式系 数特性的分析
5.4.1 通量密度及其离散表达式
由于
J* J
1、在对流项中心差分的数值解不出现振荡的 参数范围内，在相同的网格节点数下，采 用中心差分的计算结果要比采用迎风差分 的结果误差更小；
2、一阶迎风格式离散方程的系数aE及aW永远 大于零，因而无论在任何条件计算下都不 会引起解的振荡，永远可以得出看似合理 的解；
3、由于一阶迎风的截差阶数低，除非相当密 的网格，其计算结果的误差较大；
w界面上： uw 0 , W ; uw 0 , P
与中心差分格式的区别：迎风差分界面上的 未知量恒取上游节点的值，而中心差分取的
是上、下游节点的算术平均值。
(u)e Fee P max( Fe,0) E max( Fe,0) P[ Fe,0]E Fe,0 uw Fww W max( Fw,0) P max( Fw,0) W[ Fw,0]W Fw,0
程：
d (u) d ( d )
dx
dx dx
若边界条件为：x = 0， 0 ；x＝L， L。
则方程的解为：
0 exp( ux / ) 1 exp( Pex / L) 1 L exp( uL) 1 exp( Pe) 1
其中Peclet数Pe uL。
足，aP仍为相邻各系数的和。aE, aW包括了 扩散与对流作用的影响。
2、特性分析
网格Pe数：
P
ux
常物性下(1)式可写为：
P
(1
1 2
P )E
(1 2
1 2
P )W
5.2.3 对流项的迎风格式 1、两种离散方式下的迎风格式
⑴ Taylor展开法 （如下图） 以流动方向而言，P点的一阶导数永远 是该方向上的向后差分，永远从上游 获取构成一阶导数所必须的信息
C[B(P ) A(P )] D[A(P ) B(P )]
要使此式对任何，的组合都成立，只有 ： B(P ) A(P ) 0 ，即： B(P ) A(P ) A(P ) B(P ) 0，即： A(P ) B(P ) 如下图：
混合格式综合了中心差分和考虑迎 风作用两方面的因素，定义式为：
0 Pe 2
aE De
1
1 2
Pe
1
2 Pe 2 Pe ,
1 2 Pe ,
0
Pe Pe 2
5.3.3 指数格式（exponential scheme）
w
(x)w
]
E
[
e
(x)
e1Βιβλιοθήκη 2(u)e
]
w[
w
(x)
w
1 2
(
u)
w
]
令 F u ，D （扩导）则上式可变为：
x
aPP aEE aWW 式⑴
a a a aE
De
1 2
Fe
aW
Dw
1 2
Fw
p
EW
在数值计算中，若连续性方程始终得到满
动量方程的压力梯度项处理涉及到 压力与速度的耦合问题。
5.1.1 对流项离散格式的重要性
对流项离散格式是否合适将会影响：
⑴ 数值解的准确性（假扩散误差） ； ⑵ 数值解的稳定性 ； ⑶ 数值解的经济性 。
5.1.2 构造对流项离散格式的两种方式
1、Taylor展开方式 对于节点上的一阶导数给出其相应的离散
2、节点值表示的界面总通量密度计算式 将分析解
0 exp( ux / ) 1 exp( Pex / L) 1 L exp( uL) 1 exp( Pe) 1
代入通量密度定义式得：
J
F[0
0 L ]
exp( Pe) 1
把式（2）用于计算界面总通量密度Je, Jw:
5.4.3 系数特性的推论
对5种3点格式的任何一种，若在P >0时， A(P)的计算式为已知，则在 P P P 的 范围内，A(P)，B(P)的计算式均可得出。 对于A(P)：当P <0 ，按和差特性和对称 性有：
A(P ) B(P ) P A(P ) P A( P ) P
中心差分(CD):
aE
De
1 2
Fe
De (1
1 2
Pe )
aW
Dw
1 2
Fw
Dw (1
1 2
Pw )
对同一界面 Pe Pw P ，De Dw D 于是有：
aW
(i 1) D
aE (i) D
(1
1 2
P
)
(1
1 2
P
)
P
迎风差分(FUD)：
aE De Fe,0 De 1 Pe,0
aW Dw Fw,0 Dw 1 Pw,0
对同一界面，有：
aW (i 1) aE (i) 1
D
D
P ,0
1
P ,0
P
只要知道
利用系数A，B得特性，把B用A得表达式写出，
得:
A(Pe ) A( Pe ) Pe ,0
(3)
B(Pw) A(Pw) Pw) A( Pw ) Pw,0 (4)
Pe数表示对流与扩散作用的相对大小.当Pe的 绝对值很大时，导热或扩散作用就忽略.
5.2.2 对流项的中心差分
1、定义及系数的构成
用控制容积积分法时，中心差分相当于界 面上取分段线性的型线。将控制方程对控 制容积P作积分，对均分网格，离散方程为：
P
[
1 2
(u)e
e
(x)e
1 2
(u)w
( /x)
J
u
d dx
x
[P
d ] d(x /x)
所以
J*
J
( / x)
P
d d(x / x)
如图 界面上的表达式为： J * Bi Ai1
5.4.2 系数A、B间关系的分析 1、和差特性
当 i i1 时，界面上的扩散通量零，
热流问题的数值计算
Numerical Simulations of Thermal & Fluid Problems
第五章 对流—扩散方程的离散格式
主讲 李炎锋
2008年7月 北京
5.1 对流项离散格式的重要性及两 种离散格式
非线性对流项的处理涉及到对流项 的离散格式（物理过程观点：对流作 用带有强烈的方向性）；
方式，如表5-1。
例如：对一维均分网格，节点P一阶导数的 中心差分为：
E W i1 i1
x P 2x
2x