第5讲 离散方程求解.pptx

第五讲 离散方程的求解 屠基元 教授 清华大学 墨尔本皇家理工大学 2 控制方程向代数方程的转化I 均匀发热无限大平板的稳定导热问题。 x y z A B TA=100 oC TB=400 oC L q =500 W/m3 大平板区域的有限体积离散 3 控制方程向代数方程的转化 II 对于恒定的传热系数和热流量，方程变为： 引入线性近似梯度： 将上述方程进行重新整理，得到代数方程： 4 控制方程向代数方程的转化III 对于节点2和节点3，系数为 ： 对于节点1，系数为： 对于节点4，系数为： 5 控制方程向代数方程的转化 IV 控制体每一节点的系数值 代数方程的矩阵形式： 代入 6 二维离散化I 一维网格 二维网格 7 二维离散化 II 因为 代数方程式 以上方程适用于内部点: 对于节点或点 8 矩阵求解 I 对于线性代数： 反演方法效率不高 高斯消元法并非最好的方法 对于方程个数较多的方程组，可以运用迭代法求解 9 矩阵求解 II 每个节点未知变量 的一般形式可以写成： 重新排列上述方程，用雅可比方法迭代： k+1 为新迭代级 对于高斯-赛德尔迭代方法 k 为上一步迭代级 收敛性判据 10 迭代法 – 雅可比方法 利用雅可比迭代方法求解方程组的演示 11 迭代方法– 高斯-赛德尔方法 利用高斯-赛德尔迭代方法求解方程组示例 12 二维热传导– 高斯-赛德尔方法 利用高斯-赛德尔迭代方法求解方程组 差分方程的求解 将前面得到的差分方程改写为， 或者简单的写成矩阵的形式， 其中， 差分方程的求解 差分方程的求解 与方程（1）对比，知， 由方程（1）系数阵[A]的特殊性，通常称之为三对角方程（Tri-diagonal equation) 三对角方程可以采用非常高效的追赶法（TDMA法）求解 基于矩阵分解 属于必须掌握的内容 差分方程的求解 TDMA法Fortran源程序 计算机实现：算例 求解下面的一维稳态导热问题： 计算机实现：算例 求解区域的离散化： 内节点法：先划分控制容积，在确定节点 均匀网格：x=x 将整个求解区域划分为（N-2）个控制容积，N个节点（包括2个边界节点） 内部节点的差分方程 计算机实现：算例 其中， 计算机实现：算例 注意：采用内节点法划分网格时，近边界节点与其它内部节点不尽相同，所以必须单独考虑。 当i＝2时 (x)w=½ x w= W= 1 所以，当i＝2时 计算机实现：算例 所以，当i＝2时， 计算机实现：算例 当i＝3，4，…，N-2时， 计算机实现：算例 同样，当i＝N-1时 (x)e=½ x e= E= N 所以，当i＝ N-1时 计算机实现：算例 所以，当i＝ N-1时， 计算机实现：算例 最后得到由（N－2）个方程构成的方程组为 求解上面的方程，即可得到（N－2）个未知数，即，T2, T3, T4,……., TN-1。 计算机实现：算例 注意：上面的方程组是非线性的，必须用迭代法求解 求解方法： 假定一个温度分布：Ti，i＝1，2，3，。。。，N 计算i ，i＝1，2，3，。。。，N 计算a, b, c, d 用TDMA法求解方程组，得到新的温度分布： Ti’ 计算：Max{abs(Ti -Ti’), i＝1,2,3,……,N} 判断： abs(Ti -Ti’)是否小于（精度要求） 如果不能满足精度要求，令Ti ＝Ti’，重复上面的计算 满足精度要求：计算结束 计算机实现：算例 希望大家用计算机完成上面的计算，并与下面的分析解结果比较： 特别提示 计算机实现的基础地位 关键：掌握循环变量的使用 基础：对算法清晰透彻的把握 保障：细心细心再细心 29 迎风格式I 包含对流项和扩散项的控制方程： 网格雷诺数 30 迎风格式 II 收敛条件 为了避免数值振荡解 如果 31 带有中心差分的数值振荡 32 迎风格式III 例子如下： 然而 如果 使用向后差分 使用向前差分 迎风格式 一阶精度