传热学第4章报告.pptVIP

4.1 导热问题数值求解的基本思想 4.2 内节点离散方程的建立 将上两式相加  对于二维稳态导热问题，在直角坐标中，其导热微分方程为： 4.3 边界节点离散方程的建立 及代数方程的求解 判断迭代是否收敛的准则： k及k+1表示迭代次数； —第k次迭代得到的最大值 当有接近于零的t 时，第三个较好 4.3 边界节点离散方程的建立 及代数方程的求解 主对角线占优：对于常物性导热问题所组成的差分方程组,迭代公式的选择应使每一个迭代变量的系数总是大于或等于该式中其他变量系数绝对值之和,此时用迭代法求解代数方程一定收敛. 4.4 非稳态导热问题的数值解法 1.一维非稳态导热时间—空间区域的离散化 非稳态项的离散 向前差分： 向后差分： 中心差分： 4.4 非稳态导热问题的数值解法 2.一维非稳态导热微分方程的离散方法 （1）泰勒级数展开法 中心差分 向前差分 显式差分格式 优点：计算工作量小； 缺点：受时间及空间步长的限制。 4.4 非稳态导热问题的数值解法 中心差分 向后差分 优点：不受时间及空间的步长影响； 缺点：计算工作量大。 隐式差分格式 4.4 非稳态导热问题的数值解法 优点：①不受网格是否均匀限制； ?????????? ②不受物性是否为常数限制。 （2）热平衡法 4.4 非稳态导热问题的数值解法 网格傅立叶数 网格毕渥数 4.4 非稳态导热问题的数值解法 * * 第4章 热传导问题的数值解法 重点内容： ① 掌握导热问题数值解法的基本思路； ?????② 利用热平衡法和泰勒级数展开法建立节点 的离散方程。 掌握内容：数值解法的实质。 热平衡法建立节点的离散方程。 了解内容：了解非稳态导热问题的两种差分格式 及其稳定性。 本章重点内容 作业 4-9 引言 求解导热问题的三种基本方法： 理论分析法； (2)数值计算法； (3)实验法 一.数值解法的基本概念 1.实质： 把原来在时间、空间坐标系中连续的物理量的场，如导热物体的温度场等，用有限个离散点上的值的集合来代替，通过求解按一定方法建立起来的关于这些值的代数方程，来获得离散点上被求物理量的值。 ???? 数值解：离散点上被求物理量值的集合。 2.基本思路：（1）物理模型简化成数学模型是基础； （2）建立节点离散方程是关键； （3）一般情况微分方程中，某一变量在某一坐标方向所需边界条件的个数等于该变量在该坐标方向最高阶导数的阶数。 二.数值求解的步骤 建立控制方程及定解条件 确定节点（区域离散化） 建立节点物理量的代数方程 设立温度场的迭代初值 求解代数方程 是否收敛 解的分析 改进初场 是 否 二.数值求解的步骤 （1）建立控制方程及定解条件 控制方程： 定解条件： 二维矩形域内稳态无内热源，常物性的导热问题 二.数值求解的步骤 （2）区域离散化（确立节点） 步长：相邻两节点间的距离 元体：每个节点都可以看成是以它为中心的一个小区域的代表，把节点代表的小区域称为元体，又叫控制容积。 节点：网格线的交点，用（ ， ）表示 二.数值求解的步骤 （3）建立节点物理量的代数方程（离散方程） ?对节点（ , ）的代数方程，当 时，有： 二.数值求解的步骤 （4）设立迭代初场 （5）求解代数方程组 求解时遇到的问题： ① 线性； ② 非线性； ③ 收敛性等。 线性方程组：代数方程中各项系数在整个求解过程中不再变化； ②非线性方程组：代数方程中各项系数在整个求解过程中不断更新。 ③是否收敛判断：是指用迭代法求解代数方程是否收敛，即本次迭代计算所得之解与上一次迭代计算所得之解的偏差是否小于允许值。 （6）解的分析 热流量 热应力 热变形 建立离散方程常用的方法 （1）Taylor（泰勒）级数展开法 (2)控制容积平衡法(也称为热平衡法) x y n m (m，n) M N 二维矩形域内稳态无内热源，常物性的导热问题 一.泰勒级数展开法 根据泰勒级数展开式，用节点(m,n)的温度tm,n 来表示节点(m+1,n)的温度tm+1,n 用节点(m,n)的温度tm,n来表示节点(m-1,n) 的温度tm-1,n x y n m (m，n) M N 截断误差 移项整理即得二阶导数的中心差分： 同样可得： 一. 泰勒级数展开法 = 其节点方程为： 若 ，则 一. 泰勒级数展开法 二. 热平衡法 基本思想：对每个有限大小的控制容积应用能量守恒，从而获得温度场的代数方程组，它从基本物理现象