CFD程序设计.ppt

第四章 热传导 4.1 本章的对象 学习目的： 通过对热传导方程的研究找到一种普适的解决物理问题的方法 如：位势流动，质量扩散，通过多孔介质的流动，以及某些充分发展的通道流。 此外，电磁场理论，热辐射的扩散模型，以及润滑流动是由传导型方程控制的现象中的另一些例子。 4.2 一维稳态热传导 物理模型简化： 为了简化问题突出重点，略去对流项，考虑稳态状况，进而得到简化后的基本方程。 4.2-1 基本方程 作为解释四项基本法则的工具，在3．3节中介绍说明性例子的过程中，我们已经推导了稳态一维导热问题的离散化方程。对这个问题的控制微分方程是： 由该方程推导出离散化方程（控制容积法） 式中 说明： 至于分布假设，梯度dT/dx已经由T对x的分段性的变化算，而对于线性化的源项，假设Tp值代表整个控制容积的值。 当然，应当记住：只要不违背四项基本法则，选择其它形式的分布曲线是可能的，也是允许的。 这里，选择分布曲线的基本原则是在四项基本法则的约束范围内宁可采纳简单一些的分布曲线，而只是在那些必要的场合，才采用高级而复杂形式的分布曲线。 4.2-2网格间距 确定网格间距的思想： 根据因变量随x变化程度来确定网格的划分状况。即，在T~x变化比较陡（或者称梯度较大）的区域内需要进行细分网格。 确定网格间距的方法： 采用仅仅几个网格点进行试探性计算，达到给定精度后停止。 4.2-3 界面导热系数 研究主要内容： 如何通过一定的方法使非线性的导热系数线性化 在方程4.3中，已经把导热系数ke用来代表通过控制容积面e的k值；类似地，kw代表界面w的k值。当导热系数k是x的函数时，我们往往只知道在网格点W 、P、 E上的k值。下列讨论并不是针对均匀导热系数的情况而言的。 不均匀导热系产生的原因： 1 材料的不均匀性（如组合的材料板）引起； 2 温度分布的不均匀，导热系教随温度而变化。 说明：在处理对的通用微分方程时，扩散系数将起着与导热系数相同的作用，因此根据相似理论可以做类似的处理。 这里面介绍两种常用的处理方法： 算术平均值法和调和平均值法。前者直截了当，但不能精确处理组合材料中导热系数突然变化，后者利用了在热传导过程中q的不变性可以很好地解决这一问题。 （一）算术平均值法： 假设k在P点和E点之间呈线性变化，线 性插值得到： 其中插入因子fw是用图4．I所示的距离所定义的一个比值： 如果界面e位于两个网格点之间的中点，那么fe将是0．5，而ke就是kp与kB的算术平均值。 （二）调和平均值法： 我们主要关心的不是导热系数在界面e上的局部值。我们的主要目标是要得到描述的界面热流密度q的良好的表达式： 让我们来讨论这样一种情况：围绕着网格点P的控制容积由具有均匀导热系数的材料所填满，而围绕着E点的控制容积则由导热系数为的材料所填满．对于在P点与E点之间的组合板，根据稳态无内热源一维导热的分析可以得出下列的结果： 把方程（4.6）-（4.8）合并在一起，就得到： 当界面l位于P和E之间的中点时，我们有 于是有： 将方程（4.9）应用于系数的定义[式(4.3)]中，可以得到表达式， 这一公式的效能可以由下列两种极限情况很快看出 这意味着在一个绝热层的表面上热流密度为零，正像它理应存在的情况那样．相反，算术平均公式在这种情况下地均会给出非零的热流密度。 这个结果具有两方面的含义： 1）方程（4.13）表明界面 的导热系数完全与无关。因为围烧着P点的材料导热系数高，其热阻与围绕着E点的材料相比可以忽略不计（算术平均公式却保留Kp对Ke的影响）。 2）ke不等于kE，而是它的1/fe。我们的目的是通过方程（4．7）得到一个正确的qe值。应用方程（4.13）得出： 当 时，温度Tp将一直扩展到界面e， 而温降 将实际上发生在距离 内。于是正确的热流密度将由方程（4．14）给定. 换句话说， 在方程（4.13）中的因子fe可以看成是对方程（4，7）中所用的名义距离 的补偿．