哈尔滨工业大学 计算传热学 第四章 扩散方程的数值解法及其应用-2013.ppt

§4.1 一维稳态导热 §4.2 一维非稳态导热 §4.3 多维非稳态导热问题的全隐格式 §4.4 源项及边界条件的处理 §4.5 TDMA和ADI方法（交替方向隐式） §4.6 管内充分发展对流换热的数值求解 §4.6.1 管内充分发展对流换热 §4.6.2 纵向内肋片管内充分发展对流换热 X向动量方程。（r、θ向化掉） 边界条件：固体表面：w=0，对称线： ， 假设： 1. 肋厚不计。 2. 为定值。 3. 不计流向导热。 R θ r w是x向流动速度 注意 能量方程 边界条件： （？号—方向） ， ， ； ， （在肋片上） （无肋处） ， 肋片数 将流动方程转化为导热型方程 1）充分发展区 2）设 则动量方程为 即 （A） 带源项的二维（柱坐标）导热型方程 固面上： 对称线 将能量方程化为导热方程 1）定义无量温度，需一个已知的参考温度， 或是一个具有温度量纲的已知量。 令 由 得 则 2） A─横截面积 对流换热量 壁面换热量 能量平衡 即 （B） （肋片） （无肋） 3）由于边界条件均为第二类，故解不确定，需附加条件 （C） 求解过程： 1）给定N、h，求解（A）。其中Γ=1、S=1，得 2）选定w，求解（B）。其中Γ=1、 3）检查上述 解是否符合（C），若不满足（C），则 计算结果讨论 1）阻力系数f 层流光管（达西公式） 而纵肋管 2） 有一峰值，有最优值 * * 第四章 扩散方程的数值解法及其应用 F(x)：与坐标系和截面形状有关的计算因子 S：内热源。 W P E w e △x W P E （线性分布） 设 为温度的线性函数，沿x轴阶梯分布 W P E 整理： b 得： 一维稳态导热方程的离散形式 （1）由 热阻 看出 是热阻的倒数，即热导，导热能力。 （2） →控制容积中的发热量 （3） 是界面值，实际上只有节点值： a)算术平均线性分布 e W P E 讨论： b)调和平均 当 算术平均 调和平均 不合理 合理 此时热阻应取决于导热系数小的 用控制容积法积分可得： 显式： 隐式： C–N： 一、格式 二、求解 显式：可以逐个节点计算，但受 步长限制。 1. 稳及半稳：需所有节点联立求解。 写成矩阵形式 系数矩阵为三对角矩阵，可以采用“追赶法”即TDMA算法求解。 2. 二维、直角 w e N P W E S 整理 ， ， ， 热惯性 对边界节点补充代数方程，适用于外节点离散法. 在内节点法中，边界节点W的控制容积厚度为零。 由 壁面导热系数 目的是消除边界节点温度 注意 N P E S △x 处理边界条件的方法 附加源项法，适用于内节点离散法 2. 1. 将上式两边— 得 而 一、第二类边界条件 附加源项 qB已知 二、第三类边界条件 则 所以对第三类边界条件不仅有附加常数源项，而且还有附加源项的斜率项 qB未知 三、编通用程序 采用附加源项法，很容易以通用的节点离散方程计算内部和邻近边界节点的温度。而边界节点的温度为： 或 四、源项线性化处理 a. 只要S不是常数，采用S=f（T）加快迭代。 b. 计算收敛。 c. 大，计算稳定，收敛慢。 小，变化快，容易发散。 局部线化处理： 常数项 S在P点随T变化曲线的斜率 2. 广义源项：所有不能包括到非稳态项、对流项和扩散 项中的其它项。 1. 满足线性方程组对角占优条件!!! 随着维数增加，若采用隐式计算则系数矩阵为五对角或七对角阵，不能直接用TDMA方法，但可以在每个方法逐步使用。 一、TDMA法 最后： 消元： 1. 回代： 2. 二、ADI方法 有许多种，介绍P–R方法 记 为 … ADI–PR 每一个方向均可用TDMA法。 稳定条件： 全显 全隐 Peaceman-Rachford ADI 定义：无量纲过余温度不变 充分发展区 简单： →层流→导热型方程 复杂： 紊流 速度进口段：L 温度进口段： 平均 x 只有在一定的热边界条件下才能实现充分发展对流换热。 注意 采用柱坐标，控制方程： x r R 把对流问题化为导热问题 已知： 假设： 1）物性为常数。 2） 3）管壁热阻不计 方程化简： 1）充分发展的简单管流 层流 则 2）管内 则 3）因 则： 边界条件: （对称） 引入无量纲过余温度将能量方程常微分方程 则 又 由 ， ，常物性 得 又 即 令 无因次径向坐标 令 仅与x有关 仅与 有关 →两端必等于同一常数 取 是因为 永远与 异号 叫特征值 右端 这类似于边界层流动的相似变换 把偏方化成常方 边界条件: 记 则 由于方程和边界都是齐次的，必须附加条件才能得到θ的确定唯一解。 由流体平均温度？量加权平均 （A） A式为圆柱坐标的一维导热型方程，其中 为源项 注意