差分思想说课.doc

10.2导热数值解法 10.2.1数值求解基础 什么是数值解？ 就是把难以作出解析解的导热问题的微分方程描述作某种区域离散化后变成在离散点上可以求解的代数方程组，再通过对这些离散区域的代数方程组求解得到各个离散点上的温度值。这些温度点是分布在研究区域上的温度，如果这些点足够密集，它们可以描述物体的空间温度分布，或说温度场。我们把这样求得的离散点上的温度值叫温度场数值解。这种求解方法叫数值化方法。 目前常用的数值化方法有：有限元法、有限差分法、有限体积法。这里我们介绍有限差分法。 空间离散化 把一个研究的连续区域划成多个不重叠的小区域，即用网格划分区域，网格的交点叫节点。网格划分有两种方法：1.先选定节点，然后确定控制体界面；2.先定控制体界面，然后以控制体中心定为节点。下面介绍第一种方法如图：代表边界节点，代表内部节点 节点标记 如右图： 差分方程的获得 直观法 如图令 在处的估计值 1 同样 2 如果足够小，上面两式可以作为处的近似值。我们把1式叫做点处的向后差分，2式叫做点处的向前差分，另外定义 3 叫做点处的中心差分。 同样函数在处的二阶导数值： 点处的中心差分 同样向前和向后差分： 点处的向前差分 点处的向后差分 这种方法可以很直观地写出差分方程，但缺点是不能估计截断误差，从而对如何选择差分时没有依据。 泰勒级数法 同样考察前面的函数。 在点展开成泰勒级数 4 再对在点展开成泰勒级数 5 由4得到： 由5得到： 4-5得到： 或者用和展开代入 这样通过泰勒级数展开得到： 向前差分 向后差分 中心差分 这样不仅得到一阶导数差分格式，同时得到截断误差估计。 同样我们用4+5可以得到： ……中心差分 我们看到用泰勒级数展开法可以同时估计截断误差。 10.2.2稳态导热有限差分方程建立 稳态导热微分方程： (1)内节点差分方程 为了推导简单，我们这里取等距离网格则代入微分方程中得到： 这就是稳态导热内节点的差分方程（分别是三维、二维、一维问题的差分格式），如果写成编制程序的通用形式： 三维： 这样根据节点编号；；，根据这样的节点编号后内节点共有个方程，而未知数共有个，为了使这样一个节点上的差分代数方程组封闭，必须要同时建立边界节点差分方程[]个,才能保证方程组封闭。 (2)边界节点差分方程的导出 A.对第一类边界条件，边界节点温度已知（当然这就是边界节点差分方程） B.对第二、三类边界条件，边界节点方程处理方式有两种：补充边界节点方程（代数方程）；另一种方法是设法从内节点方程中消去未知的边界点温度，这样只有个方程，也只有个未知数，方程组封闭。结果在临近边界节点方程中增加了源项，也称附加源项法。 补充边界节点方程法： 一般是对边界节点控制体作能量平衡。如右图一个二维问题的第二类边界条件，对边界节点P控制体建立能量平衡： 同样考虑 第二类 推广到三维： 第二类 对第三类边界条件，这时上面的，这样代入上式得到： 第三类边界条件 顶角节点控制体补充方程： 二维 三维 附加源项法： 是把所有边界都按第一类边界条件处理，对非第一类边界条件的情况，把边界条件上所规定的界面热流以附加源项的形式送入与边界相邻的控制体P节点的差分方程： 节点P的差分方程就是内节点方程： 其中是边界节点，对非第一类边界条件是未知的；现在我建立W点能量平衡方程： 实际上每个边界节点可以写出一个这样的关系式， 这些关系式联立就可以求出用由内节点温度和表示的边界节点表达式，代入相邻边界节点的内节点方程中，从而消除了边界节点温度未知数，同时使内节点方程多了源项。故称附加源项法。 (3)差分方程求解 对稳态导热，得到的差分方程组是一个线性代数方程组。求解有高斯消元法：对变量不是很多的方程组，把系数矩阵消成一个三角形矩阵，而后由三角形矩阵依次求得方程的解；迭代求解法：首先任意给定初值，然后由初值求得第一次近似解，判断误差，不合格，再由第一次解代入求第二次近似解，….，以此类推，直到： 其中是根据工程实际人为规定的误差。达到这一值后能够满足工程精度要求。 把计算得到的空间节点上的温度值作出等值线（二维问题）或等值面（三维问题）图就是一个很好的数值解。 10.2.3非稳态导热有限差分方程建立 （1）以一维不稳态导热为例： 在和范围积分（用有限体积积分法得到差分方程） 右边： 令 右边= 左边= 上式整理得： 其中：，