传热学导热问题的数值解法课件.pptx

《传热学导热问题的数值解法课件.pptx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第四章　导热问题的数值解法主讲人：孙晴本章知识结构第四章导热问题的数值解法1、导热问题数值求解的基本思想2、稳态导热离散方程的建立和求解方法3、非稳态导热的数值解法24.1导热问题数值求解的基本思想3数值解法的基本思想:用导热问题所涉及的空间和时间区域内有限个离散点的温度近似值来代替物体内实际连续的温度分布，将连续温度分布函数的求解问题转化为各节点温度值的求解问题。通过对各离散节点建立代数离散方程，将导热微分方程的求解问题转化为节点温度代数方程的求解问题。4数值解法的基本内容与步骤:（1）对实际导热

2、问题的几何、物理性质进行分析，建立控制方程及定解条件。（2）求解域离散化：用与坐标轴平行的网络线将所涉及的空间和时间区域划分成有限个子区域，将网络线的交点作为节点,每个节点就代表以它为中心的子区域（元体或称为控制容积），节点温度就代表子区域的温度。（3）建立节点温度代数方程。（4）设立迭代初场。4.1导热问题数值求解的基本思想5（6）对计算结果进行分析，若不符合实际情况，则修正上述步骤，重复进行计算，直到结果满意为止。（5）求解节点温度代数方程组，得到所有节点的温度值；4.1导热问题数值求解的基本思

3、想6是否求解域离散化改进初场建立控制方程及定解条件设立迭代初场建立节点温度代数方程求解代数方程组解的分析是否收敛？导热问题数值解法的流程图4.1导热问题数值求解的基本思想4.2稳态导热离散方程的建立和求解7用有限差分近似微分，用有限差商近似微商（导数）建立节点离散方程的方法有两种：泰勒级数展开法热平衡法1.泰勒级数展开法根据泰勒级数展开式，用节点(i,j)的温度ti,j来表示节点(i+1,j)而温度ti+1,j用节点(i,j)的温度ti,j来表示节点(i-1,j)的温度ti-1,j84.2稳态导热离

4、散方程的建立和求解一、内节点若取上面两式右边的前三项，并将式两相加，得二阶导数的中心差分：同样可得：9截断误差未明确写出的级数余项中的ΔX的最低阶数为24.2稳态导热离散方程的建立和求解对于二维稳态导热问题，在直角坐标中，其导热微分方程为：其节点方程为：10向前差分格式向后差分格式中心差分格式4.2稳态导热离散方程的建立和求解2.热平衡法114.2稳态导热离散方程的建立和求解根据傅里叶定律，L,R,T,B各节点向P节点的导热量：12Z方向取单位长度4.2稳态导热离散方程的建立和求解13则内热源发热量

5、在稳态导热下：则4.2稳态导热离散方程的建立和求解二、边界节点离散方程的建立1.位于平直边界上的节点144.2稳态导热离散方程的建立和求解红色框内为平直边界上P点的网格单元。边界条件：第三类边界热平衡式为：15若有内热源，则可在上式左边加上内热源项：绝热边界，只要令h=0，即可。4.2稳态导热离散方程的建立和求解对于第二类边界条件，设为qw则上式中的对流换热项,用热流密度qw来替代164.2稳态导热离散方程的建立和求解2.外部角点174.2稳态导热离散方程的建立和求解3.内部角点184.2稳态导热离

6、散方程的建立和求解19求解代数方程的迭代法代数方程组的求解直接解法——矩阵求逆，高斯消元法等迭代法——简单迭代法，高斯-赛德尔迭代法，4.2稳态导热离散方程的建立和求解对于常物性导热问题所组成的差分方程组，迭代公式的选择应使每一个迭代变量的系数总是大于或等于该式中其它变量系数绝对值之和，即对角占优，此时用迭代法求解代数方程一定收敛。20迭代过程是否已经收敛的判据：迭代过程能否收敛的判据：K代表迭代次数4.2稳态导热离散方程的建立和求解214.3非稳态导热问题的数值解法对非稳态导热，从能量关系来看网格

7、单元不仅与相邻的网格单元之间有热量传递，本身的内能也将随时间发生变化。一维非稳态导热问题用有限差分法求解空间，时间的坐标的划分22以无内热源，物性参数为常数的一维非稳态导热为例，导热微分方程为：温度对坐标的二阶导数：非稳态项，温度对时间的一阶导数有三种不同的格式。根据泰勒级数展开，若用节点(i,k)的温度参数来表示节点（i.k+1）的温度4.3非稳态导热问题的数值解法23舍去上式中左边第三项及以后各尾项，移项整理此式是一阶导数向前差分的表达式。类似可得向后差分的表达式：同理可得中心差分的表达式：4.

8、3非稳态导热问题的数值解法24若温度对时间的一阶导数采用向前差分，则导热微分方程可改写为：上式移项整理4.3非稳态导热问题的数值解法25稳定性条件对于点i上k+1时刻的温度是由该点在第k时刻的基础上综合相邻点温度后得出。所以如k时刻i点温度较高，则其下一时刻的温度也较高，反之亦然。其表现在差分方程的稳定性条件是：各项的系数必须大于或等于零。即：若温度对时间的一阶导数采用向后差分，则4.3非稳态导热问题的数值解法26整理后可得隐式差分格式：优缺点：隐式格式计算工作量大，