第五章对流扩散方程的离散格式
第五章——对流-扩散问题的有限体积法
混合格式兼具中心差分格式和迎风差分格式的优 点，具有守恒性、有界性和迁移性，其缺点是按 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
对流扩散问题有限体积法
流体仿真与应用第八讲二、对流-扩散问题的有限体积法◆中心差分格式（例子）节点增加到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 )
输运方程对流扩散方程
输运方程对流扩散方程输运方程是描述物质传输过程的数学模型，常见的有对流扩散方程。
对流扩散方程是由对流和扩散两种机制共同产生的输运过程来描述的，它的一般形式为：∂c/∂t+∇·(v*c)=∇·(D*∇c)其中，c表示物质的浓度或者响应变量，t表示时间，v表示流体的速度场，D表示物质的扩散系数，∇表示梯度运算符。
对流项描述了物质的对流运动，即物质随着流体的移动而移动。
对于三维坐标系来说，对流项可以表示为∇·(v*c)。
具体来说，对流项的每一项分别表示了物质在x、y和z方向上的携带速度与浓度梯度的乘积。
扩散项描述了物质由浓度高处至浓度低处的扩散现象，即物质自发性地从高浓度区域向低浓度区域传播。
扩散项可以表示为∇·(D*∇c)，其中D是扩散系数，表示物质扩散的速率与浓度梯度的乘积。
对流扩散方程的物理意义是描述了物质在流体中传输的速率与物质浓度梯度之间的关系。
通过对流项，方程能够描述物质随着流体的运动快速传输的现象；而通过扩散项，方程能够描述物质由浓度高处向浓度低处传输的现象。
综合考虑对流和扩散的作用，对流扩散方程能够比较准确地描述物质在流体中的传输过程。
对流扩散方程在科学和工程领域有广泛的应用。
例如，在污染物传输和扩散模拟中，对流扩散方程可用于描述污染物由源区到周围空气或水体的传输过程。
在热传导模拟中，对流扩散方程可用于描述热量由高温区域到低温区域的传导过程。
在物质传递过程中，对流扩散方程也被广泛应用于描绘物质的传输行为。
总结起来，对流扩散方程是一种常见的输运方程，它能够描述物质由流体传输并扩散的过程。
通过对流项和扩散项的综合作用，对流扩散方程能够比较准确地描述物质在流体中的传输行为，所以在科学和工程领域有着广泛的应用。
对流扩散方程clank
对流扩散方程clank标题：对流扩散方程的概述引言概述：对流扩散方程是数学中常见的描述物质传输过程的方程。
它在众多领域中都有广泛的应用，如流体力学、热传导、质量传输等。
本文将从五个大点出发，详细阐述对流扩散方程的相关内容。
正文内容：1. 对流扩散方程的基本概念1.1 对流扩散方程的定义1.2 对流扩散方程的一般形式1.3 对流扩散方程的物理意义2. 对流项与扩散项的影响2.1 对流项的作用2.2 扩散项的作用2.3 对流项与扩散项的相互作用3. 对流扩散方程的解析解与数值解3.1 解析解的求解方法3.2 数值解的求解方法3.3 解析解与数值解的比较4. 对流扩散方程的边界条件和初值条件4.1 边界条件的选择与影响4.2 初值条件的确定与影响4.3 边界条件和初值条件的耦合效应5. 对流扩散方程的应用领域5.1 流体力学中的应用5.2 热传导中的应用5.3 质量传输中的应用总结：对流扩散方程是描述物质传输过程的重要方程，其基本概念包括方程的定义、形式和物理意义。
对流项和扩散项是方程中的两个关键因素，它们分别对物质传输起到对流和扩散的作用，并且相互作用影响着传输过程。
对流扩散方程的求解可以采用解析解和数值解两种方法，它们各有优劣，需要根据具体情况选择。
边界条件和初值条件是方程求解中必要的条件，它们的选择与确定对结果有重要影响。
对流扩散方程在流体力学、热传导和质量传输等领域都有广泛应用，它为我们理解和解决实际问题提供了重要的数学工具。
总之，对流扩散方程是一个复杂而重要的数学方程，它在物质传输过程中起着关键作用。
深入理解和研究对流扩散方程，对于解决实际问题具有重要意义。
第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
一类二维稳态对流——扩散方程的有限差分法
一类二维稳态对流——扩散方程的有限差分法一维稳态扩散方程描述了物质在一维空间中的扩散行为。
然而，在某些情况下，我们需要研究物质在二维平面中的扩散行为，例如热传导、流体传输等。
本文将介绍一类二维稳态对流-扩散方程的有限差分法。
二维稳态对流-扩散方程可以写作：∇·(D∇u) + ∇·(cu) + fu = 0 —— (1)其中，D是扩散系数，c是速度场，u是待求解的物理量，f是源项。
在这个方程中，第一项表示物质的扩散项，第二项表示对流项，第三项表示源项。
我们需要求解方程(1)，找到u的分布。
为了应用有限差分法来求解二维稳态对流-扩散方程，需要将二维空间离散化为一个网格。
假设我们将x方向离散为Nx个等距的节点，y方向离散为Ny个等距的节点，那么我们可以得到一个(Nx+1)×(Ny+1)的网格。
我们在网格节点上定义未知量u，然后将方程(1)对节点处的u进行离散化。
首先，我们对方程(1)的扩散项进行离散化。
我们使用五点差分格式来近似二维Laplace算符∇·(D∇u)。
对于网格节点(x,y)，我们可以得到以下差分格式：(Dij(xi+1,yj)ui+1,j + Dij(xi-1,yj)ui-1,j +Dij(xi,yj+1)ui,j+1 + Dij(xi,yj-1)ui,j-1 -4Dij(xi,yj)ui,j) / ∆x^2 + (Dij(xi,yj)ui,j) / ∆y^2其中，∆x和∆y是网格步长，Dij是扩散系数。
接下来，我们对方程(1)的对流项进行离散化。
我们使用中心差分格式来近似二维梯度算符∇·(cu)。
对于网格节点(x,y)，我们可以得到以下差分格式：(cxi+1/2,yj(ui+1,j - ui,j)) / ∆x + (cxi-1/2,yj(ui,j - ui-1,j)) / ∆x + (cyi,j+1/2(ui,j+1 - ui,j)) / ∆y + (cyi,j-1/2(ui,j - ui,j-1)) / ∆y其中，cxi+1/2,yj、cxi-1/2,yj、cyi,j+1/2和cyi,j-1/2是速度场在节点(x,y)处的中心点处的x和y分量。
第五章对流扩散问题(假扩散)
该问题的数值解如下：
1 n i
u t n u t n (1 ) i i 1 x x
MUD : du d ux d (( ) ) dx dx 2 dx
第五章 对流扩散问题———假扩散
由图可以看出，在区间 P 2 ，中心差分格式预报的 P 值优于迎风格式的预报值。对比这两种格式，其扩散项 的处理是完全相同的，所不同的仅仅是对流项的处理上 ，在中心差分格式中对流项的差分格式具有二阶精度， 而在迎风格式中对流项的差分格式只具有一阶精度。在 区间 P 2 ，两种格式预报 P值所表现出的差异性恰恰是 这两种格式精度不同的体现。观察上图，迎风格式所预 报的 P值具有该高不高和该低不低的特点，这一特点正 是由一阶精度迎风格式所引起的扩散系数为 ux / 2 的 假扩散项造成的。也反映了假扩散项的影响。
n n ux ut 2 n ( ) i u( ) i (1 )( 2 ) i O( x 2 , t 2 ) t x 2 x x
由此可以看出，我们前边得到的差分方程所逼近的是 一个非稳态对流扩散问题，而非原型问题所要求的非 稳态对流问题。
第五章 对流扩散问题———假扩散
1 n (1 P
ut n ut n ) P W x x
用编号法表示
1 n i
1 n , n 在点 (i, n) i i 1
u t n u t n (1 ) i i 1 x x
做Taylor展开
n n u 2 n 1 2 n ( ) i u( ) i ( 2 ) i x ( 2 ) i t O( x 2 , t 2 ) t x 2 x 2! t
第五章 对流扩散问题———假扩散
将 i 1 和 i 1 台劳 级数展开代入
第五章对流扩散方程
• 比混合格式复杂，计算量增加，但准确性 提高
• 稳定性：若用中心差分格式不能体现对流 项的物理本质，常会引起数值解的振荡
• 经济性：若用高阶格式，无数值振荡，但 格式复杂，求解相对困难，机时消耗较多
5.2 一维稳态对流扩散问题
d (u)
dx
d dx
d
dx
5.2.1 模型方程的精确解
d (u)
dx
d dx
d
dx
边界条件：
x 0, 0；x L, L
采用迎风思想：从来流上游方向找依赖区
在界面e上,若 ue 0,则e P ;若 ue 0,则e E ; 在界面w上,若 uw 0,则w W ;若 uw 0,则w P
界面流量
• 引入符号 • 对流：
a1, a2 max(a1, a2 )
(u)e Fee P Fe,0 E Fe,0 , (u)w Fww W Fw,0 P Fw,0
• 控制方程变为： dJ 0; 或 J const dx
J
F
0
0 L
exp(Pe) 1
界面上通量
Jw
Fw
W
W P
exp(Pw )
1Hale Waihona Puke JeFeP
P
exp(
E
Pe )
1
Fe exp( exp(Pe
Pe )
) 1
Fw exp(Pw )
1 P
Fe exp(Pe )
1E
Fw exp(Pw exp(Pw )
) 1
W
合并整理结果
aPP aEE aWW
• 系数
aE
第五章 对流-扩散方程的离散格式
见下页表格：
5.3.5 5种三点格式系数计算式的汇总 不同格式离散方程的形式相同，但 系数不同。具体见下表5－1：
5.4 对流-扩散方程5种3点格式系 数特性的分析
5.4.1 通量密度及其离散表达式
J J ( / x )
*
由于 所以
d d J u [ P ] dx x d ( x / x)
i i 1 d ， ui 0 dx i x i 1 i x
ui 0
对多维问题，用此方法构造的对流 项的离散格式，只有在求解区域内 流速不发生逆向时，所形成的离散 方程才具有守恒性。
2、控制容积积分法定义
规定界面上的未知量恒取上游节点的值 e界面上： ue 0 , p ； ue 0 , E
把式（2）用于计算界面总通量密度Je, Jw: 对Je： ， , L x
0 P L E e
P E J e Fe [ P ] exp( Pe ) 1
对Jw：
0 W ， L P , L xw
W P J w Fw [W ] exp( Pw ) 1
对于坐标系I，C位于界面之后，而D位 于界面之前，于是： J * B( P )C A( P ) D 对于坐标系II，D位于界面之后，而C 位于界面之前，于是：
J B( P ) D A( P )C
*
由于
J J
*
*'
C [ B( P ) A( P )] D [ A( P ) B( P )]
exp( Pe ) 1
Fe ；
Fw exp( Pw ) aW exp( Pw ) 1
对流扩散方程.
A对流扩散方程的求解对流扩散问题的有效数值解法一直是计算数学中重要的研究内容，求解对流扩散方程的数值方法主要是有限差分法(FDM)、有限元法(FEM)、有限体积法(FVM)、有限解析法(FAM)、边界元法(BEM)、谱方法(SM) 等多种方法。
但是对于对流占优问题，用通常的差分法或有限元法进行求解将出现数值震荡。
为了克服数值震荡，80年代，J.Douglas,Jr.和T.F.Russell 等提出特征修正技术求解对流扩散占优的对流扩散问题，与其它方法相结合，提出了特征有限元方法、特征有限差分方法、特征混合元方法；T.J.Hughes和A.Brooks提出过一种沿流线方向附加人工黏性的间断有限元法，称为流线扩散方法(SDM)。
有限差分法、有限元法、有限体积法是工程应用中的主要方法。
对流扩散方程的特点对流扩散方程右端第一项为扩散项，左端第二项则是对流项。
由于其方程本身的特点，给建立准确有效的数值求解方法带来一定的困难。
对流和扩散给流体中由流体携带的某种物理量的变化过程，可以通过一个无量纲的特征参数(Peclet数)来描述，Peclet数Pe的定义为：Pe=|ν|L/D。
这里v是来流速度，L是特征长度，D是物质的扩散系数。
如果Pe数较小，即对流效应相对较弱，这类问题中，扩散占主导地位，方程是椭圆型或抛物线型；如果Pe数较大，即溶质分子的扩散相对于流体速度而言是缓慢的，这类问题中，对流占优，方程具有双曲型方程的特点。
对于对流占优问题的求解，采用常规的Galerkin有限元方法，为了避免求解结果产生数值振荡，获得稳定解，则应使每个单元的局部Peclet数，Peh=|ν|h/D≤2，这里h为单元的最大尺寸，|v|为单元中的最大速度分量值。
因此，用本文方法求解对流占优对流扩散问题，要得到稳定解，则要通过加密有限元网格来实现。
第五章对流扩散问题(假扩散 高阶格式第一个内节点的差分方程)
第五章 对流扩散问题———假扩散
1.0 0.8 0.6 0.4 0.2 0.0 0.0 t=0 AS t=4 UD c=1.0 t=8
0.6 0.8 1.0 x c值进一步增大，我们看到：c=0.010.050.0750.40.81.0 ，而 FD/FDc=0.01=10.960.9340.6060.2020，可见此时UD格式的整体假扩散 已经衰减了100%，使其计算结果和精确解完全一样。
第五章 对流扩散问题———假扩散
以上结果再次证明了我们关于时间项和对流项采用一阶
格式时，对其假扩散的分析。 也证明了我们提出了二阶
迎风格式和QUICK格式确实能够消除由对流项离散而引
起的假扩散。
第五章 对流扩散问题———假扩散
例 设有下列一维非稳态对流问题
u , t x 0 x 1, u 0.1
n i 1 / 2
n i 1 / 2
)
n i
ut c x
第五章 对流扩散问题———假扩散
其中，
n 1 i
c(
n i 1 / 2
n i 1 / 2
)
n i
ut c x
对SUD
n n i 1 / 2 i 1 / 2
n n 3 n 4 i 1 i2 i 2
第五章 对流扩散问题———假扩散
UD格式
x 1 / 100
xi xi 1 x
t cx / u
t=0
i= 1
2
3
4
98
99
100
101
0 i 由初始条件确定
t=t
i= 1
2
3
4
对流扩散方程
所 以 迎 风 格 式 的 稳 定 性条 件 是
h2 2 ah
当a 0时 ， 情 况 类 似 ， 稳 定 性条 件 是
h2 2 a h
也可以利用中心显格式来讨论稳定性,于是将上面格式改为:
u n1 j
u
n j
a
un j 1
只需验证 G 1,由于 1 cos wh 0，条件转化为：
4 4 2 (1 cos wh) r 2 (1 cos wh) 0
即 4 - 2r 2 (4 2 r 2 )(1 cos wh) 0
由于 1 cos wh 0,1,上述不等式转化为
2
4 2r 2 0, 4 2r 2 2(r 2 4 2 ) 0
为了简单方便,设a>0,先对方程作扰动,得到另外一对流
扩散方程
u t
a u x
1
1 R
2u x 2
其中R 1 ha
2
对上面的方程构造迎风格式
u n1 j
u
n j
a
u
n j
u
n j 1
h
1
1 R
un j 1
2u
n j
h2
u
n j 1
称为逼近对流扩散方程的Samarskii格式.
n j
u n1 j
2 得到如下差分格式：
u n1 j
u
n j
a
u
n j 1
u
n j 1
2h
（
2
a
2）u
n j 1
2u
n j
h2
u
n j 1
稳 定 性 分 析 完 全 类 似 于中 心 差 分 格 式 ， 显 然 有
数值传热第五章课件2陶文铨
主讲陶文铨西安交通大学能源与动力工程学院热流中心CFD-NHT-EHT CENTER2010年10月18日, 西安数值传热学第五章对流扩散方程的离散格式(2)对流项离散格式的重要性及两种离散方式5.5.1假扩散的含义与成因5.5.2一阶截差格式引起严重假扩散举例1.本来的含义2.扩充的含义3.Taylor 展开法的分析5.5关于假扩散的讨论5.5.3网格倾斜交叉引起的计算误差5.5.4 非常数源项引起的假扩散5.5.5 两个名例以一维非稳态纯对流过程为例俩分析，其中有两n nφφ2(,O x φΔΔ其中关于时间的二阶导数项可做如下变化：时才没有这部分的计算误差。
2. 扩充的含义现有文献中常常将较大的计算误差都称为假扩散，大致有以下几项原因：(1) 一阶导数的一阶截差格式；(2) 流动方向与网格线呈倾斜交叉；(3) 离散格式未计及非常数源项的影响。
5.5.2一阶截差格式引起严重假扩散举例1.一维稳态对流扩散问题对流项用FUD，扩散项用CD，当Pe较大时，数值计算结果严重偏离精确解。
Physically plausible solution纯对流传递纯对流传递由离散方程：1n−1此时只有对流，没有扩散！时则有严重假扩散！0.8C =0.8C =当时，产生了严重的扩散作此种误差称为流向假扩散Γ≠Γ气流01. 设UE对P 控制容积，有2. 设控制容积，此时：计算误差纯对流传递三个对流问题的归纳这就是假扩散纯对流传递3)网格倾斜交叉引起的计算误差E冷热流体之间产生了温度均匀化的过程，即交叉5.5.5 已知流场计算温度场232(1),2(1)u y x v x y =−=−−参考解xT严重假扩散2) Leonard细高方腔中的自然对流换热5.6.1采用高阶格式克服流向假扩散5.6可以克服或减轻假扩散的格式与方法5.2.2 克服、减轻交叉假扩散的方法1. 采用二阶迎风2.采用三阶迎风3. 采用QUICK 格式1. 采用有效扩散系数2.采用自适应网格4. 采用SGSD 格式可以克服或减轻假扩散的格式与方法相当于界面上的中心差分)W WWxφ+Δ如型线上凹，则(2) FVM向上游取两点定义界面插值2.采用三阶迎风展开定义－一阶导数的三阶偏差分格式3. 采用定义－界面的插值在中心差分基础上考虑曲中心差分插值率修正？需要满足两个条件：插值的正确修正：相邻(2)0W PE φφφ−+<型线下凹8Cur −对e-界面u e 小于零时，取,,W P φφφu e 大于零时，取怎样相邻的三点?QUICK(2)e φφ=1/2w i φφ−=有：4. 采用CD条件稳定，但没有二阶假扩散；二阶迎风绝对稳定，组合起来，但是：如何确定值，特别是如何由计算结果来5. 高阶格式实施中的问题f u f计算边界：固o2) 代数方程的求解：等时，5.6.2用减小扩散系采用自适应网格（以减轻流5.7 对流-扩散方程离散形式稳定性分析5.7.1 数值计算中常见的三种不稳定性5.7.2 分析对流项格式不稳定性的“符号不变原则”5.7.3 稳定性分析结果讨论5.7.4 对流项格式问题讨论小结2.“符号不变”原则的基本思想3. “符号不变”原则的实施步骤4. “符号不变”原则的实施例子1. 研究背景扩散方程离散形式稳定性分析也会产生振荡的解，称为对流项离散格式的不稳态定性，研究目的是，找出产生振荡的临界Peclet 数。
哈尔滨工业大学 计算传热学 第五章 对流-扩散方程的离散格式-2013
乘方 | 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级数展开，保留二阶。
Pe x T T0 exp( cux / k ) 1 exp( L ) 1 TL T0 exp( cuL / k ) 1 exp( Pe ) 1
贝克立数 Pe
cuL
k
cuAT
kA T L
对流传热量/导热量
TL
1
纯导热 P e 0 上游信息对流到下游， 纯对流 下游信息无法通过扩 P e
Pe 10
aE Pe DE aE (1 0.1Pe )5 Pe DE aE (1 0.1Pe )5 DE aE 0 DE
10 Pe 0
0 Pe 10
Pe 10
(f)
aE [| 0, (1 0.1| Pe |)5 |] [| 0, Pe |] De 讨论
而
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
上述若对任何成立必得根据对称特性可以说明为什么前面讨论格式特性只研究函数ab系数性质的含义根据对称及和差特性我们仅须知道红线的值即在0105中心迎风混合指数乘方10迎风指数乘方中心混合54关于假扩散的讨论原始的假扩散概念一维非稳态对流方程纯对流没有扩散差分格式没有考虑非常数源项的影响
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
aP = aE + aW
aE = De – Fe / 2 aW = Dw + Fw / 2
在流场的实际求解过程中， 每一个迭代层次上，即使速度 场尚未收敛，也要保证连续方 程是满足的。
3. 对流项的中心差分与迎风差分
3.2 对流项的中心差分
三点说明：
系数 aE , aW 包含了对流 F 与扩散 D的作用的影响；
对均分网格：
2. 对流项离散格式的重要性及两种离散方式
2.2 构造对流项离散格式的两种方式 （2）控制容积积分
给出界面上被求函数的插值方式
2. 对流项离散格式的重要性及两种离散方式
2.2 构造对流项离散格式的两种方式 （3）两种定义之间的关系
对某种对流项离散格式，都可以用两种方法给出其相应 的定义；
3.1 一维对流-扩散问题模型方程的精确解
上游优势
3. 对流项的中心差分与迎风差分
3.1 一维对流-扩散问题模型方程的精确解
希望所构建的离散方程形式也具有这样的物理特性
3. 对流项的中心差分与迎风差分
3.2 对流项的中心差分 Central Scheme （CS）
分段线性
均分网格
令
对流项
----界面上的流量
1. 简 介
对流与扩散作用在物理本质上的区别
从物理过程来看，扩散作用与对流作用在传递信息或扰动方面 的特点有很大区别：
扩散是由分子的不规则热运动所致，分子不规则热运动对空间不同方向
的几率都是一样的，因而扩散过程可以把发生在某一地点上的扰动的影响 向各个方向传递。
对流是流体微团宏观的定向运动，带有强烈的方向性。在对流的作用下，
两种定义方式的截断误差阶数是一致的，均为二阶截差 （中心差分，分段线性）；
Taylor 展开： 逼近 P点导数值 控制容积积分：逼近 控制容积内导数的积分平均值
3. 对流项的中心差分与迎风差分
3.1 一维对流-扩散问题模型方程的精确解
该数学描写的精确解（解析解）：
3. 对流项的中心差分与迎风差分
可以证明，扩散项具有二阶截差的中心差分确实具有 把扰动向四周均匀传递的特性。中心差分离散格式能够很 好地反映扩散过程的的特点；
对大多数有实际意义的问题，这一离散格式能很好地反 映扩散过程的特点，且已完全能满足需要。数值计算结果 误差的主要来源，在于对流项的离散格式。
1. 简 介
证明：扩散项的中心差分可以将扰动均匀地向四周传递
3.2 对流项的中心差分
特性分析： 将x=L/2带入，得P点精确解：
中心差分解：
3. 对流项的中心差分与迎风差分
3.2 对流项的中心差分
特性分析：
P∆ 精确解 P 中心差分 P
3. 对流项的中心差分与迎风差分
3.2 对流项的中心差分
令
扩散项
----界面上的扩导
3. 对流项的中心差分与迎风差分
3.2 对流项的中心差分
aP
aE
aW
aP = aE + aW + ( Fe - Fw )
3. 对流项的中心差分与迎风差分
3.2 对流项的中心差分
求解过程中连续方程始终满足： Fe = Fw
对节点 i：
− =
∆
0
0
−+
∆
=
−
对节点 i+1：
0
−
∆
=
=
=−
00
−
+
∆=ຫໍສະໝຸດ 1. 简 介0对节点 i-1：
− =
∆
=
00
−
+
∆
=
若取
=0.25 <0.5（稳定性条件），则n时刻的扰动到n+1时刻如图：
n 时刻
n +1 时刻
/2
i-2 i-1 i i+1 i+2 x
/4
i-2 i-1 i i+1 i+2 x
3. 对流项的中心差分与迎风差分
3.2 对流项的中心差分
特性分析：
3. 对流项的中心差分与迎风差分
3.2 对流项的中心差分
特性分析：
在一维模型方程离散求解的均匀网格中，已知 W 100 ， E 200 。试对P∆=0, 1, 2, 4 四种情形，按中心差分计算 的 P 值，并解释结果。
3. 对流项的中心差分与迎风差分
第五章 对流-扩散方程的离散格式
王娴
西安交通大学航天学院
内容
简介 对流项离散格式的重要性及两种离散方式 对流项的中心差分与迎风差分 对流离散格式假扩散特性 克服或减轻假扩散的格式 对流项离散格式讨论小结
1. 简 介
描写对流换热的控制方程包括质量守恒，动量守恒，能量 守恒。这些微分方程中最高阶导数项为二阶（扩散项）： 离散 中心差分 （无方向性，容易）
扩散作用下扰动的传递
1. 简 介
两个一阶导数项：非线性对流项、压力梯度项。数值处理 方法远比二阶扩散项复杂的多。
不可压缩流场数值解的关键问题都是由这两个一阶导数项 的离散所引起的。
1. 非线性对流项处理 对流项的离散格式 （本章） 2. 动量方程中压力梯度项的处理 压力与速度之间的耦
合关系 （下章）
对数值计算及其结果而言，对流项的离散方式构造得是否 合适，影响到： 1. 数值解的准确性 （一阶截差严重的数值结算误差，假扩散问题） 2. 数值解的稳定性（一定条件下解的震荡，中心差分的条件稳定） 3. 数值解的经济性（复杂格式，求解耗内存与时间，加速）
2. 对流项离散格式的重要性及两种离散方式
2.2 构造对流项离散格式的两种方式 （1）Taylor展开方式
De , Dw 部分是由扩散项的中心差分所形成， 代表了扩散作用的影响； 与流量有关的部分则是界面上分段线性型线 在均匀网格下的表现，体现了对流的作用。
3. 对流项的中心差分与迎风差分
3.2 对流项的中心差分
三点说明：
所谓不同格式，具体表现在这两部分表达形式 的不同上；
本章如无特别说明，扩散项均取中心差分，其 表达式即为
发生在某一地点上的扰动只能向其下游方向传递而不会逆向传播。
为某一物理量的扰动，t0 为初始 时刻， t1 ，t2 表示相继时刻；虚 线所示图形表示在扩散或对流作用 下扰动的传递情形。
2. 对流项离散格式的重要性及两种离散方式
2.1 重要性
数学观点：对流项是一阶导数项，离散不存在困难。 物理过程观点：对流作用带有强烈的方向性 难 ！
离散扰动分析法：采用非稳态显式的某种格式来研究该格式传递扰动的特性。
研究一维非稳态扩散方程
=
的显式格式:
− =
∆
−+ ∆
开始时刻物理量场均匀， 处处相等，且假定其为零。在某一时刻 n， 某一节点 i 上突然有一个扰动 ，其余各点上扰动为零，随着时间的推
移，这一扰动传递的情形可按上述差分方程来确定：
1. 简 介