第4章导热问题的数值解法.pptVIP

4.3边界节点离散方程的建立及代数方程的求解对于第一类边界条件的热传导问题，处理比较简单，因为已知边界的温度，可将其以数值的形式加入到内节点的离散方程中，组成封闭的代数方程组，直接求解。对于第二类或第三类边界条件的导热问题，所有内节点的离散方程组成的代数方程组是不封闭的，因未知边界温度，因此应对边界上的节点补充相应的代数方程，才能使方程组封闭，以便求解。第30页，共66页，星期日，2025年，2月5日xyqw为了求解方便，将第二类边界条件及第三类边界条件合并起来考虑，用qw表示边界上的热流密度或热流密度表达式。为使结果更具一般性，假设物体具有内热源Φ(不必均匀分布)。第31页，共66页，星期日，2025年，2月5日边界节点(m,n)只代表半个元体，若边界上有向该元体传递的热流密度为qw，据能量守恒定律：4.3.1边界节点离散方程的建立(1)平直边界上的节点qwxyqw第32页，共66页，星期日，2025年，2月5日(2)外部角点如图所示，外角点仅代表1/4个以为边长的元体。假设边界上有向该元体传递的热流密度为，则据能量守恒定律得其热平衡式为：xyqw第33页，共66页，星期日，2025年，2月5日(3)内部角点内部角点代表了3/4个元体，在同样的假设条件下xyqw第34页，共66页，星期日，2025年，2月5日讨论关于边界热流密度的三种情况：（1）绝热边界即令上式即可。（2）值不为零（3）对流边界此时，将此表达式代入上述方程，并将此项中的与等号前的合并。对于的情形有：流入元体，取正，流出元体，取负第35页，共66页，星期日，2025年，2月5日（a）平直边界（b）外部角点（c）内部角点第36页，共66页，星期日，2025年，2月5日4.3.2处理不规则区域的阶梯型逼进法当计算区域出现曲线边界或倾斜边界时，常常采用阶梯形的折线来模拟真实边界，然后用上述方法建立边界节点的离散方程。第37页，共66页，星期日，2025年，2月5日一个大的工业炉由一根很长的耐火黏土砖柱支撑，其截面尺寸为1m×1m。在稳定运行时，柱子的三个表面保持在500K，另一个表面暴露在300K，h=10W/㎡·K的空气流中。利用Δx=Δy=0.25m的网格，确定柱中的二维温度分布及单位长度柱子对空气的散热速率。第38页，共66页，星期日，2025年，2月5日第39页，共66页，星期日，2025年，2月5日第40页，共66页，星期日，2025年，2月5日4.3.3代数方程的求解方法2）迭代法：先对要计算的场作出假设（设定初场），在迭代计算中不断予以改进，直到计算前的假定值与计算结果相差小于允许值为止的方法，称迭代计算收敛。1）直接解法：通过有限次运算获得精确解的方法，如：矩阵求解，高斯消元法。第41页，共66页，星期日，2025年，2月5日2迭代法目前应用较多的是：1）雅可比迭代法（简单迭代）：每次迭代计算，均用上一次迭代计算出的值。2）高斯——赛德尔迭代法：每次迭代计算，均是使用节点温度的最新值。第42页，共66页，星期日，2025年，2月5日在计算后面的节点温度时应按下式（采用最新值）例如：根据第k次迭代的数值可以求得节点温度：第43页，共66页，星期日，2025年，2月5日设有一三元方程组：其中（i=1,2,3；j=1,2,3）及是已知的系数（均不为零）及常数。第44页，共66页，星期日，2025年，2月5日采用高斯——赛德尔迭代法的步骤：（1）将三元方程变形为迭式方程：第45页，共66页，星期日，2025年，2月5日（2）假设一组解（迭代初场），记为：并代入迭代方程求得第一次解每次计算均用最新值代入。（3）以新的初场重复计算，直到相邻两次迭代值之差小于允许值，则称迭代收敛，计算终止。第46页，共66页，星期日，2025年，2月5日判断迭代是否收敛的准则：k及k+1表示迭代次数；—第k次迭代得到的最大值当有接近于零的t时，第三个较好第47页，共66页，星期日，2025年，2月5日迭代能否收敛的判据1）对于一个代数方程组