而
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 )