传热学第四章培训教案.ppt

* 传 热 学 第四章 导热问题的数值解法 数值解法 有限差分： 物理概念明确实施方法简便 边界元法： 将力学中的微分方程的定解问题化为边界积分方程的定解问题，再通过边界的离散化与待定函数的分片插值求解 有限元法： 将连续的求解域离散为一组单元的组合体，用在每个单元内假设的近似函数来分片的表示求解域上待求的未知场函数，近似函数通常由未知场函数及其导数在单元各节点的数值插值函数来表达。从而使一个连续的无限自由度问题变成离散的有限自由度问题。 是在微分方程中用差商代替偏导数，得到相应的差分方程，通过解差分方程得到微分方程解的近似值 §4-1 导热问题数值求解的基本思想 数值求解基本思想： 把原来在时间、空间坐标系中连续的物理量的场，用有限个离散点上的值的集合来代替，通过求解按一定方法建立起来的关于这些值的代数方程，来获得离散点上被求物理量的值。这些离散点上被求物理量值的集合称为该物理量的数值解。 y x n N M §4-1 导热问题数值求解的基本思想 具体作法： 1）建立控制方程及定解条件 2）区域离散(连续域离散成不连续节点) 3）建立节点物理量的代数方程 4）求解代数方程组 5）解的分析 前处理 后处理 x y x y n N m M (m,n) (m,n+1) (m,n-1) (m+1,n) (m-1,n) §4-2 稳态导热的有限差分法 研究问题： 二维、矩形、无内热源、常物性、稳态导热。 1 区域离散 节点：（m,n） 空间步长： △ x，△y 当△ x=△y均匀网格 单元，控制容积←→（m，n） 网格线： △ x，△y任意 2 建立节点物理量的代数方程 两种方法 物理法（能量守恒） 数学法 1）内节点（m，n） 由能量守恒：净导入微元体热量=0 x y n m (m,n) (m,n+1) (m,n-1) (m+1,n) (m-1,n) 由能量守恒：净导入微元体热量=0 时： 2） 边界节点 a）平直边界上的节点（m，n） 净导入（m,n）单元体热量 由能量守恒： 对流流入（m,n）单元热量 + =0 x y n m (m,n) (m,n+1) (m,n-1) (m-1,n) h,tf b）外部角点（m，n） 净导入（m,n）单元体热量 由能量守恒： 对流流入（m,n）单元热量 + =0 x y n m (m,n) (m,n-1) (m-1,n) h,tf 时： x y n m (m+1,n) (m,n-1) (m-1,n) (m,n+1) (m,n) c）内部角点（m，n） 净导入（m,n）单元体热量 由能量守恒： 对流流入（m,n）单元热量 + =0 h,tf 时： 3 代数方程的求解 求解方法 直接：（高斯消元） 迭代：（高斯—赛德尔迭代法） 迭代法：假设初值，迭代初场，重复迭代 判据：当 停止迭代 网格的无关性 1. 判断 2. 在其它条件相同的情况下，判断下列哪种情况先达到热平衡 3. 在其它条件相同的情况下，判断下列哪种情况先达到热平衡 λ1=236w/mK;ρ=2710ke/m3;Cp=902J/kgK。 λ1=39.2w/mK;ρ=7570ke/m3;Cp=470J/kgK。 4. 试述Fo,Bi的物理意义。 λ1 λ2 t3 t2 t1 δ δ x t δ δ x t 0 h,tf h,tf δ δ λ1，λ2谁大？ §4-3 非稳态导热的数值解法 1 区域离散化（空间、时间） x为空间坐标，将区域划分N-1等份，得N个节点。 i为时间坐标，将时间坐标上计算域划分I个节点，I-1等份。 x τ i I n N (n,i ) (n,i+1) (n,i-1) (n+1,i ) (n-1,i ) 一维、非稳态 时间步长 2 建立节点物理量的代数方程 物理法——热平衡法 非稳态的离散有三种格式 向前差分 向后差分 中心差分 以一维问题为例 内部节点： 净导入n的热量=内能增加 △x n-1 n n+1 △y 1 左侧： 右侧： 净导入n的热量： 内能的增加 相应净导入热量 向前差分 向后差分 中心差分 内部节点： 净导入n的热量=内能增加 △x n-1 n n+1 △y 1 若内能增加项选前差格式，净导热量选i层 为“显式差分格式”——一旦i时层上个节点的温度已知，可立即算出（i+1）时层上个内点温度，不必求解联立方程式。 其优点是计算工作量小 缺点是对时间步长及空间步长有一点限制 若内能增加项选前差格式，净导热量选（i+1）层 △x n-1 n n+1 △y 1 3 稳定性判据 步长限制——稳定性判据 对一维、非稳态、显式格式 n=1,2,3…N-1 隐式格式的缺点是计算工作量大 优点是对步