第六章 对流扩散方程的差分格式.pptVIP

* 第六章 对流扩散方程的差分格式 导热型方程：（原始或经过变换的） 二阶导数项（扩散），源项 对流扩散方程：（动量或能量） 二阶导数项（扩散），源项 一阶导数项（扩散），压力梯度。 一维稳态无内热源的对流扩散方程: 密度， 扩散系数。 精确解 对流热能量方程 贝克立数 对流传热量/导热量 纯导热 上游信息对流到下游，下游信息无法通过扩散传到上游。 纯对流 以上精度解为例讨论各种差分格式的性能。 x TL 1 -5 Pe=0 1 5 L T0 0 §6-1 中心差分与迎风差分。 连续方程 常数 而 对控制容积采用分断线性进行积分。 记 W w e E 则 记 即 记 以网格间距 为特征尺寸的贝克列数。 在均匀网格和常物性条件下，上式为: 即 （B） （A） (5)中心差分格式只重于低雷诺问题 (2)式（A）为中心差分格式 (4)导致不合理的原因 而使 在 大时， 分布偏离线性很远，积分假设不合理。且 大时，在x=L/2处， 上游值而非平均值。 (3)存在的问题，用（B）式计算。 (1)当满足连续性条件时， 讨论 取 （a） 得 （b） 得 显然不合理，真实值 在 和 之间。 为了避免上式造成物理不真实，而构造迎风差分格式：（介绍第二类，有守恒性）在控制容积分时，界面上的未知量恒取上游节点的值，而不象中心差分那样取两边节点的平均值。 取最大值 保证 二阶导数的扩散项仍采用中心差分， 整理后得 但 当u0 当u0 （4）当 很大时，在x=L/2处, 而近风格式仍取中心差分，放大了扩散作用。 讨论 （3）由于近风格式在界面上取 等于取 储水箱→搅拌均温→控制容积→壁面导热→扩散→管道→对流 进入管道的水温等于上游水箱的温度，但它不知道将流入的箱中的任何情况。 放大了对流作用。 故该格式没有考虑大小的影响。 储水箱 管道 （1）方程（c）的系数没有负数，解在物理上总是真实的。 （2）格式的物理含义 ?一、指数格式 记总通量 对流和扩散， 则 控制容量积分 把精度解代入J中 § 6-2 指数格式混合格式与乘方格式 故 而 (D) 指数 区别就在函数 二．混合格式 虽然指数格式是精确解,但计算过繁,通过对 随 变化及其三条切线 斯波尔丁提出 (e) 讨论 (1)在 之间， ，就是中心差分格式。 (2)在 的区域里 即扩散项取零的逆风格式。 由于其综和了中心和迎风两种格式的优点，故称混合格式。但最好还是看成精确解的三条直线的近似，包络线。 三、乘方格式 由于混合格式在 附近，偏离真值较远，帕坦卡提出乘方格式。 讨论 (1)该格式计算量比指数小，且与指数格式的解差别很小。 (f) A和B的性质的讨论 (1)当 时，扩散量=0， 完全由对流造成，即 § 6-3 通用表达式 为了在讨论中引入 记 界面i+ 上的值可以有界面两侧节点值表示 （y） 得B-A= 此为和差特性。 i i+1 i+1/2 x （2）对称特性 对坐标I，D在界面前。 坐标II相对I反向，流向不变，而此时D在界面后。（界面的前后相对坐标方向而言，并非对流向） 由于坐标描述的是同一量，故 I C C D D 是 的函数，而 是与之对称的。 上述若对任何成立，必得 即 的值以 =0 的轴对称的。 根据对称特性可以说明为什么前面讨论格式特性只研究函数 因为 A A B B A，B系数性质的含义 根据对称及和差特性，我们仅须知道红线的值，即在 则所有 和 均确定了。 在 通用格式 对B 各种格式的通用表达式 A B 1 P W w e E 而 根据通量守恒 即 显然不论那种格式，仅仅是 表达式的区别。 中心近风混合指