4-1_导热问题的数值解法讲义.pptVIP

4-3 非稳态导热问题的数值解法 (2) 热平衡法 三、建立节点物理量的代数方程（离散方程） 右边第三类边界 显示格式存在稳定性问题：如果节点 tn(i) 前面的系数小于零，则数值解出现不稳定的震荡结果。 显示格式 显示格式：格式右边全部为第 i 时间层的温度值，只要 i 时间层温度已知，即可计算得到 i+1 时间层的温度。 即：空间步长?x和时间步长?τ的选取有限制 显示格式的稳定性条件： 4-3 非稳态导热问题的数值解法 四、一维导热问题显式差分格式的不稳定性讨论 隐式格式 隐式格式：空间离散采用（i+1）时层的值。 隐式格式不存在稳定性问题，对时间步长和空间步长没有限制，但是计算量较大。 非稳态导热节点离散方程的两种格式： 4-3 非稳态导热问题的数值解法 2．讨论一维导热问题显式差分格式稳定性限制的物理意义 由节点显式差分方程可以看出， 前的系数可能是负值。当系数为负时， 时刻节点温度越高 ， 时刻的温度却越低，违背了热力学第二定律。 4-3 非稳态导热问题的数值解法 表面温度在不同 Δτ 时的计算值 4-3 非稳态导热问题的数值解法 4-3 非稳态导热问题的数值解法 1．首先选择空间坐标间隔 Δx ，即距离步长。对二维问题一般使 Δy=Δx　。 2．对显式格式差分方程，根据方程的稳定性条件选择允许的最大时间步长 Δτ ； 在稳定性条件允许范围内，Δτ 越大，计算工作量越小，但精度较差；对一维问题，一般取　　　　 ，即可满足工程精度要求；　　　 对于隐式差分方程，Δx、Δτ 可任意选取，不必进行稳定性条件校核； 4-3 非稳态导热问题的数值解法 五、节点方程组的求解 3. 按题意给定的初始温度分布，确定各节点上的温度初值； 4. 根椐建立的差分方程组，求 Δτ 时刻各节点的温度 tn(1) ； 5. 再由 tn(1) 为初值，换用τ=2Δτ 　　　　（即 k=2 ），重复 4 计算出 tn(2) ，如此反复，最后得到 j 时刻的 tn(j) 　。 4-3 非稳态导热问题的数值解法 * * 区域离散化：用一系列与坐标轴平行的网格线把求解区域划分成若干个子区域，以网格线的交点作为需要确定温度值的空间位置，称为节点 ( 结点 ) ，节点的位置用该节点在两个方向上的标号m，n表示。 以节点为中心，在两个节点的中心处划分界限，定出节点所代表的小区域，叫控制体。二维控制体是一个个的矩形面积。网格的步长在每一个方向上可以均匀划分，也可以不均匀的划分。同样，也可对时间坐标以一定的步长划分出时间网格。 节点(node): 网格线交点。 控制容积(control volume):节点代表的区域,其边界位于两点之间。 界面(interface): 控制容积的边界。 网格划分方法: A 先确定节点，后定界面 B 先确定界面，后定节点 * * * * * * * * 从所有方向流入控制体的总热量＋ 控制体内热源生成热＝ 0 边界内角点 三、建立节点物理量的代数方程 §4-2 稳态导热问题的数值解法 （3）边界节点的有限差分方程 边界节点离散方程中的两个问题： 边界热流密度的具体处理方法 绝热边界 第二类边界 第三类边界 不规则边界的处理方法 多段折线模拟不规则边界，网格越密越接近实际 坐标变换：保角变换 §4-2 稳态导热问题的数值解法 三、建立节点物理量的代数方程 建立节点离散方程的泰勒级数法与热平衡法的比较： 泰勒级数法属于纯数学方法，而热平衡法基于能量守恒原理，物理概念明确，且推导过程简捷； 泰勒级数法对于建立边界节点的离散方程较困难； 当导热物体物性或内热源不均匀时，泰勒级数法不适用，而热平衡法能够方便处理。 §4-1 导热问题数值求解的基本思想 三、建立节点物理量的代数方程 例1： 对如图所示的圆截面直肋的一维稳态、无内热源、常物性导热问题，试分别列出内节点m和端部节点M的离散方程式。已知圆截面直径为d。 四、设立迭代初场 §4-2 稳态导热问题的数值解法 五、求解代数方程组 写出所有内节点和边界节点的温度差分方程式： 直接解法 迭代解法 直接解法：矩阵求逆、高斯消元法等 缺点：所需内存较大、方程数目多时不便、不适用于非线性问题（若物性为温度的函数，节点温度差分方程中的导热系数不再是常数，而是温度的函数。这些系数在计算过程中要相应地不断更新） §4-2 稳态导热问题的数值解法 五、求解代数方程组 迭代解法：简单迭代（Jacobi迭代）、高斯-赛德尔（Gauss-Seidel）迭代、块迭代、交替方向迭代等 先对要计算的场作出假设（给定初始值）、在迭代计算过程中不断予以改进、直到计算结果与假定值的结果相差小于允许值。称迭代计算已经收敛。 Gau