传热学-第四章-热传导问题的数值解法资料.pptVIP

* * * * § 4-4 非稳态导热问题的数值解法 非稳态导热数值解法的特点： （1）非稳态导热微分方程多了非稳态项，因此单值性条件中增加了初始条件； （2）除了对空间域进行离散外，还需要对时间域进行离散； （3）利用热平衡法导出节点温度方程时需要考虑控制容积的热力学能随时间的变化； （4）由于时间和空间同时离散，在有些情况下空间步长和时间步长不能任意选择，否则会带来节点温度方程求解的稳定性问题。 * * 第三类边界条件下，常物性、无内热源无限大平壁的一维非稳态导热问题为例。 1) 求解域的离散 2) 节点温度差分方程的建立 运用热平衡法可以建立非稳态导热物体内部节点和边界节点温度差分方程。 1. 一维非稳态导热的数值求解： * * 内部节点温度差分方程 内部节点n所代表的控制容积在i时刻的热平衡： 时间项一阶导数采用向前差分，扩散项采用当前（i）时层上的值来表示，热平衡方程式可写成 左侧导入热量 右侧导入热量 计算第i+1时层温度值时，用第i时层的已知值 ——内部节点温度方程的显式差分格式 * * 两点结论： (a) 任意一个内部节点n在(i+1)时刻的温度都可以由该节点及其相邻节点(n-1) 、(n+1)在i 时刻的温度由上式直接求出，不必联立求解方程组，这是显式差分格式的优点。这样就可以从初始温度出发依次求出各时刻的节点温度； (b) 必须满足显式差分格式的稳定性条件，即 稳定性条件说明，一旦空间步长?x或时间步长??的数值确定之后，另一个步长的数值的就不能任意选择，必须满足稳定性条件。 物理意义: * * 隐式差分格式： 如果节点温度取下一时层(i+1)的温度值，内部节点n所代表的控制容积的热平衡方程式可写成: ——内部节点温度方程的隐式差分格式 * * 隐式格式与显式格式的区别： * * (2) 边界节点温度差分方程 边界节点0所代表的控制容积在k 时刻的热平衡： 时间项一阶导采用向前差分，扩散项采用前一时层（i）温度（显示格式），热平衡方程式可写成： 上式写成显函数的形式 * * 同内部节点温度方程的显式差分格式的道理相同，需将上一时层的信息传递到下一时层去，因此必须满足： * * 第四章小结 （2）掌握有限差分法的原理; 重点掌握以下内容： （3）能够根据导热问题的特点，合理地进行求解域的离散； (1) 理解数值解法的基本思想，熟悉数值解法的基本步骤； （4）重点掌握热平衡法建立节点温度差分方程； （5）会利用计算机求解节点温度差分方程组，并掌握求解结果表达方法。 * * 作业： 第四版： 4-3，4-10，4-15 * * x y n=1 m m=M n=N m=1 w s e n (m-1,n) (m,n-1) (m+1,n) (m,n+1) (m,n) 上海交通大学 * 上海交通大学 * 上海交通大学 * 第四章　导热问题的数值解法 第四章　导热问题的数值解法 第四章　导热问题的数值解法 第四章　导热问题的数值解法 第四章　导热问题的数值解法 第四章　导热问题的数值解法 第四章　导热问题的数值解法 作业： 第四版：3-2，3-31，3-48 * * 第四章导热问题的数值解法 Numerical Method for Heat Conduction * * 主要内容（重点掌握）： 导热问题数值求解的基本思想 内外节点离散方程的建立 非稳态导热问题的数值解法 * * 特点： (1) 分析法 a 能获得所研究问题的精确解，可以为实验和数值计算提供比较依据； b 分析解具有普遍性，各参数的物理意义及影响清晰 c 局限性很大，对复杂的问题无法求解； 前述2、3章对导热问题的求解思路： 导热微分方程+边界条件+初始条件 稳态问题：直接积分法 非稳态问题：分离变量法 解析解（analytical solution） 局限性： 简单几何形状及边界条件 工程实际中面临的大部分问题几何形状和边界条件要复杂的多，由于数学上的困难还不能给出解析解，导致目前解析解只能作为某些简单问题的参照依据，不能解决实际问题。 研究传热学问题的三种基本方法：(1)理论分析法； (2)实验法；(3)数值计算法 §4-0 引言 * * (2) 实验法: 是传热学的基本研究方法： a 偏向于机理研究； b.受场地，燃料动力源等因素的影响，无法完全复现研究对象，具有时间、空间上的局限性 c.费用昂贵 (3) 数值方法 数值方法：把原来在时间和空间连续的物理量的场，用有限个离散点上的值的集合来代替，通过求解按一定方法建立起来的关