【2017年整理】数值传热学四-数值计算.ppt

第四章 热传导;§4.1 本章的对象;二、以导热问题的数值解作为学习起点的原因;④. 把热传导用作流体流动计算方案的基本组成部分 的做法有助于理解动量传递与热量传递之间的类 似性（用某种方法把速度与温度相比拟）。;§4.2 一维稳态热传导;分布假设： 由T 对 x 的分段线性的变化算得； 源项的线性化?TP代表整个控制容积内的值，即采用阶梯性分布进行计算的。;§4.2-2 网格间距;①.对所要得到的解进行某些定性的预计，使设计得到某些指导；;§4.2-3 界面导热系数; k 值不均匀性产生的原因 由材料的不均匀性引起（如组合材料板）； 材料均匀，T分布的不均匀性也会导致k 的不均匀。;显然，这相当于线性插值。当界面e位于两个节点之间的中点时，fe=0.5, 此时;②.调和平均法;③.两种方法的比较; P控制容积是良导热体 ,按算术平均法计算, ;调和平均法可以合理求解上述两种极限情况;调和平均公式的推导是对于稳态、无内热源、k在相邻的两个控制容积之间发生阶梯式变化导出的。从定性上，串联热阻叠加的适应性不受上述条件的限制。;§4.2-4 非 线 性;收敛：计算达到最终不变的状态称之为迭代的收敛。;§4.2-5 源项的 线 性化;几点说明：;离散化方程的一般式：;源项线性化方法举例：;eg2. 已知：S=3+7 T, 可能的线性化形式有：;(3) 推荐的方法;结论：在所有负斜率的直线中，与已知曲线相切的直线通常为最佳(图中2线)；较陡的直线是可以接受的(图中3线) ；欠陡的直线是不希望采纳的(图中1线) ，它不能体现已知的S 随T 的下降速度。;§4.2-6 边 界 条 件;先定控制容积，后定节点位置 B类网格* ;3. 构建边界节点的附加方程 A类网格 ;将qi 代入上式得：;3. 构建边界节点的附加方程 B类网格* ;③. 第三类边界条件 已知 h、 tf;§4.2-7 线性代数方程的求解;正过程：消元的目的是将每个包含三个待求变量的方程，消去一个变量，使之成为两个待求变量的方程，继而把所有这些两个变量的方程新系数值计算出来，直到边界上原来只有两个待求变量的方程，经消元后变成了单变量方程，并能直接求出边界变量值为止。;TDMA算法的数学表述;3. TDMA的数学推演 ; 对比③与②即可得到Pi、Qi的递推公式：;4. TDMA的计算机程序说明;§4.3 非稳态一维热传导;2. 对整个控制容积积分上述微分方程; 一些假设可归纳为一般化的式子：;经整理得： （将上标1去掉）;§4.3-2 显式、克兰克-尼科尔森模式以及全隐式模式;?;显式格式： 前的系数有可能变为负值( 可看成是 在时间方向上的一个相邻点)，从而违背正系数法则。实际上，要使系数为正，时间步长必须足够小，以使;2. 全隐式格式 f =1;3. 克兰克-尼克尔森模式 (C-N) f =0.5 ;两点说明：;§4.3-3 全隐式离散化方程;3. 全隐格式的主要原则:;§4.4 二维与三维问题;P、E两点间e控制面上的热流密度;式中各系数为：;§4.4-2 三维问题的离散化方程;§4.4-3 代数方程的求解;迭代法：对于非线性问题，方程必须用最新的系数来重复求解，采用直接解法是很不经济的。;4. 迭代求解的两种方法: 高斯-赛德尔逐点计算法和逐行法;?计算机存贮器中所存在相邻点的温度值。本次;这个准则是充分条件，并非必要条件。即有时虽然违背了这个条件，但迭代仍可能收敛。 对于满足正系数法则的方程，如非稳态全隐格式 , 由于各系数均大于零，且;④. 对于稳态无内热源的问题，对内节点只有 然而可以从边界邻点的离散方程中找出使得不等号成 立的条件。;高斯—赛德尔点迭代计算的缺点 收敛速度慢，特别对于网格节点数多的情况尤为明显。 因为它是以一次迭代传递一个网格间距的速率来传递边界所给予的信息的(仅将扫描起始边的影响传到整个区域)。;迭代方法简介：;实际上，采用逐线迭代时，二维导热的离散化方程可改写为：;逐线迭代方向选择的原则（扫描方向）;§4.5 超松弛与欠松弛;2. 方程式的写法;3. ?值的选取;对TP产生的变化加以修正;5. i 值的选取;?? 越小， 越大，欠松弛作用愈强； 当 ??很大时， 很小，不再有欠松弛作用，而直接转化为稳态问题的求解。;线迭代：自左向