华北电力大学导热问题数值解法.ppt

第四章 导热问题的数值解法 内容简介： 本章介绍了数值求解物理问题的基本思想，以二维稳态导热问题和一维非稳态导热为例，介绍了采用有限差分法数值求解导热问题的基本过程。 4-1 导热问题数值求解的基本思想 一、 数值解法的本质 数值解法是用物理问题所涉及的空间和时间区域内有限个离散点(称为节点)的物理量近似值来代替物体内实际连续的物理量分布，将连续物理量分布函数的求解问题转化为各节点物理量值的求解问题。 二、数值解法的应用 数值解法是目前各行各业解决工程问题和进行科学研究的一种重要方法。如：航空航天、冶金、电力、电子、化工、建筑、能源动力、机械等。 三、大型的数值计算商用软件 目前市场上有许多针对某一领域的商用软件，在我们专业应用的有： FLUENT——是目前国际上比较流行的商用CFD软件包，在美国的市场占有率为60%。凡跟流体，热传递及化学反应等有关的工业均可使用。 CFX——是由英国AEA公司开发，是一种实用流体工程分析工具，用于模拟流体流动、传热、多相流、化学反应、燃烧问题。 此外还有：Phoenics和 STAR-CD 等。 四、数值解法的基本步骤 建立控制方程及定解条件 确定节点（区域离散化） 建立节点物理量的代数方程 设立物理量的初值 求解代数方程 是否收敛 解的分析 改进初场 是 否 五、二维稳态导热问题的数值求解 1、物理问题 接下来以一个二维稳态导热问题为例，介绍数值求解导热问题的的具体过程，重点是节点离散方程的建立和代数方程组的迭代求解。 二维矩形区域内的稳态、含有内热源、常物性的导热问题。 y x 0 H W h,tf q=0 h,tf q=0 2、建立控制方程及定解条件 3、区域的离散 将二维矩形区域按照下图的格式进行离散化. x y m n 节点——网格线的交点。x方向，M各节点；y方向，N个节点 区域内部的节点称为内节点；边界上的节点称为外节点。 步长——相邻两个节点间的距离，记为Δx、Δy 。 M N (m,n) (M,n) (M,N) (1,1) 4、建立节点的离散方程 ① 内节点 （1）泰勒级数展开法 根据泰勒级数展开式，在x方向，用节点(m,n)的温度t (m,n)来表示节点(m+1,n)的温度t（m+1,n） 截断误差 两式相加： 同理，y方向 由数学描述 离散方程 若Δx=Δy ——均分网格 （2）热平衡法 假定节点（m,n）所代表的控制区域四个边界分别为n、s、e、w ，节点（m,n）所代表的控制体的能量守恒方程为： (m,n) (m+1,n) (m,n+1) (m-1,n) (m,n-1) n e s w 由于温度分布未知，假设热量由四个邻点传入 从四个边界导入控制体的能量可应用傅立叶定律得到: 代入能量守恒方程式，并整理得: ② 边界节点离散方程的建立 边界节点离散方程只能采用热平衡的方法建立。 对于平直边界（如右边）上的节点 热平衡法的优点 （1）物理意义明确； （2）适用于内节点和边界节点。 对于其它外节点，如外部角点和内部角点等可按同样的方法推导，区别是其所代表的控制体体积不同。 写出所有内节点和边界节点的离散方程，n个未知节点温度，n个代数方程式： 代数方程组的求解方法：直接解法、迭代解法 5、节点方程组的求解 直接解法：通过有限次运算获得代数方程精确解; 如：矩阵求逆、高斯消元法等。 缺点：所需内存较大、方程数目多时不便、不适用于非线性问题。 迭代解法：先对要计算的场作出假设、在迭代计算过程中不断予以改进、直到计算结果与假定值的结果相差小于允许值。称迭代计算已经收敛。 迭代解法有多种：高斯-赛德尔迭代是最常用的。 取n=3的三元方程组，并作变换 假设一组初始温度， 开始迭代计算：第一轮 得到一组新的温度 以此组温度开始新一轮迭代 重复上述计算，直到相邻两次迭代值之差小于允许值 判断迭代是否收敛的准则： 6 3 10 ~ 10 - - 取 值一般 相对偏差 允许的偏差； — e e k及k+1表示迭代次数； —第k次迭代得到的最大值 对解的结果进行进一步的分析，获得定性或定量上的一些新的结论。 5.进行解的分析