[工学]No08 1013 4 导热问题的数值解法.ppt

数值求解的基本步骤： 物理问题数值求解过程的框图： (1) 泰勒级数展开法 根据泰勒级数展开式，用节点(m,n)的温度tm,n 来表示节点(m+1,n)的温度tm+1,n 若取上面式右边的前三项，并将两式相加 移项整理即得二阶导数的中心差分： 对于二维稳态导热问题，在直角坐标中，其导热微分方程为： 、在图中所示的有内热源的二维导热区域中，一个界面绝热，一个界面等温温度为t0（包括节点4），其余两个界面与温度为tf的流体对流换热，λ，h均匀，内热源强度为Φ。试列出节点1，2，5，6，9，10的离散方程式。 、直径为d=0.5mm的球，其材料的密度ρ＝8920kg/m3, 比定压热容cp＝400 J/(kg.K)，导热系数为25 W/(m.K)。其初温为20℃，突然将其置入温度t =125℃的气流中，表面传热系数为420 W/(m2.K)，问热电偶指示温度达到100℃时所需的时间；在这期间热电偶吸收的热量是多少？ 、质量为M、比热为c、表面面积为A的小型导电体，初温等于环境温度Ta ，通电加热后温度逐渐升高，最终达到平衡温度Tb , 试证：  并说明 的意义（ ，其中h是导电体与环境的对流换热系数）。 从第二式得出 向后差分 back difference 二级数相减得到: 中心差分 central finite difference 常物性一维非稳态问题: 非稳态项-向前差分 扩散项-中心差分 ——显示格式 explicit finite difference scheme可直接求解. 如果扩散项也用(i+1)时层上的值来表示,则 不能写成与显示格式类似的形式， 称为隐示格式 implicit finite difference scheme 不能直接求解，必需联立求解. 与显示格式比,计算工作量大,但它是时间步长没有限制. 边界节点的处理 式中: — 网格傅里叶数 一维无限大平板 边界: 内节点: 计算机实习： 1、补充材料中计算机实习指导书: 练习题一、二、三 5 4 13 405机房 42 传热学（能动学院） 5 4 13 404机房 42 传热学（能动学院） 5 4 12 405机房 34 传热学（能动学院） 5 4 12 406机房 56 传热学（能动学院） 5 4 11 405机房 30 传热学（能动学院） 5 4 11 409机房 60 传热学（能动学院） 单 元 日 次 周 次 机 房 上机人数 班级名称 2、作业上交： 请于第14周发到：zhhjia@yeah.net 邮件名称以姓名学号命名，附件采用压缩文件； 、等截面直肋如右图所示,碳钢?=43.2W/(m.K), 写出1，2，3，4节点的离散方程。 、 图中示出了常物性、有均匀内热源 、二维稳态导热问题局部边界区域的网格配置，试用热平衡法建立节点0的离散方程式（设 ）。 图 、 直径为d、单位长度电阻为R、发射率为ε的金属棒，初始时与温度为T∞的环境处于热平衡状态，后通过电流I，已知棒与环境的表面传热系数为 h。试导出通电流期间金属棒温度随时间变化的规律，并写出处于新的热平衡状态的条件。（不用求解） 、 图中为一扇形区域中的二维、稳态、无内热源的导热问题，其一侧平直表面维持均匀恒定温度 的表面绝热，其余表面与温度为 tf的环境换热，表面传热系数为h，为用数值法确定区域中的温度分布，在径向和周向分别用定步长 和 划分计算区域，节点布置如图中所示。假定区域的物性参数已知且为常数，试用热平衡法建立节点 1和 7的离散方程（有限差分方程式）。（不要求整理） 该方法基本思路是：从导热微分方程出发，采用泰勒级数展开方法，获得相邻节点温度和空间步长表示的一阶和二阶导数 * * 第四章　热传导问题的数值解法 导热问题一般为: 1. 分析解(analytical method): 通过对上述方程积分求得(有限情况)，也称理论解。 2. 数值解(numerical method): 用某种方式把微分方程化为关于各个离散点(节点)的代数方程,通过解代数方程获得问题近似解的方法。 ???? 把原来在时间、空间坐标系中连续的物理量的场，如导热物体的温度场等，用有限个离散点上的值的集合来代替，通过求解按一定方法建立起来的关于这些值的代数方程，来获得离散