第四节 导热问题的数值解法 Numerical Methods of Heat Conduction.docVIP

第四章 导热问题的数值解法 Numerical Methods of Heat Conduction § 4-1 导热问题数值求解的基本思想 及内节点离散方程的建立 导热问题一般为: 上述问题的解法有以下两种: 1. 理论解(analytical method): 通过对上述方程积分求得(有限情况)。 2. 数值解(numerical method): 用某种方式把微分方程化为关于各个离散点(节点)的代数方程,通过解代数方程获得问题近似解的方法。 连续——离散(任意情况) 数值求解的基本步骤 1.数学描述 2.区域离散化 3.建立节点物理量的代数方程 4.设立迭代初场 5.求解代数方程组 6.解的分析 1.数学描述 导热问题一般为: 无限长棱柱（如图）导热、沿高度各截面的温度分布相同，可简化为二维问题。 2. 区域离散化 有限差分法原理 finite difference 有限元法 finite element 边界元法 boundary element 有限分析法 finite analysis 网格划分 grid 节点(node): 网格线交点. 控制容积(control volume): 节点代表的区域 ,其边界位于两点之间. 界面(interface): 控制容积的边界. 网格划分方法: practice A 先确定节点，后定界面 practice B 先确定界面，后定节点 均分网格: 节点编号: 从小往大排 3. 代数方程的建立 （1）Taylor 级数展开法 对点(m,n)作Taylor 展开: 两式相加得: 同理 代入微分方程得: 对于正方形网格 则有: （2）热平衡法（热力学第一定律） 说明: 1 用此方法所得的边界方程是有O((x2)精度 2解析解是温度（物理量）的连续函数 3数值解得出离散点上的数值 § 4-2 边界节点离散方程的建立及代数方程的求解 1. 边界上离散方程的建立 对于边界节点要根据边界条件来确定。 1) 第一类边界条件, y=b 处将边界温度直接代入即可,方程封闭。 2)对于第三类边界条件, y=0对控制体直接应用热力学第一定律 当时,上式为: 3) 对于第二类边界条件 4）不规则边界的处理 ? 折线法 ? 坐标变换 2. 代数方程的求解 ? 直接求解 (内存大) ? 矩阵求逆 ? 消元法 ? 迭代法 (使用较多) ? Gauss—Seidel迭代 ? 点迭代 ? 线迭代 ? 块迭代 直接求解法 ? 内存大 迭代法 ? 使用较多 Gauss---Seidel 迭代 线性方程组: 可以写成: 上角标为计算序号,计算时先给出ti 的初值,然后用上式进行迭代。 终止计算的方法： 3. 例题 针肋如右图所示，碳钢 (=43.2W/(m.K)，求其温度分 布及换热量。 解： 以上是精确解，现在我们用数值方法求解： 网格划分如右图： 该问题的数学描述为 节点2： 同理得节点3: 节点4 用热力学第一定律，导入的热量应等于对流散出的热量，固有： 得 三种情况的计算结果如下温度分布 X 0 10 15 20 30 θ 175 139.5 127.9 119.7 113.4 θ 175 139.8 120.13 113.8 θ 175 128.13 114.29 热量计算： 误差 精确解 (=15.06 W 四节点 ( =11.94 W 21% 三节点 ( =10.52 W 30% 如取5 节点, 则( 的误差为 19% § 4-3 非稳态导热问题的数值解法 ? 多非稳态项 ? 扩散项的处理方法与前一样 空间坐标 x 1(N ?x 空间步长 时间坐标 ( 1(I ?( 时间步长 （n,i）代表了时间空间区域中的一个接点位置 t(i)n 将温度函数 t 在节点（n,i+1）和（n,i-1）对点（n,i）作泰勒级数展开 从第一式得出 向前差分 forward difference 从第二式得出 向后差分