径向导热界面参数插值方法研讨.pdf

中国工程热物理学会 传热传质学 学术会议论文 编号：093012 径向导热界面参数插值方法研究 王辉，徐明海 （中国石油大学（华东）储运与建筑工程学院，山东青岛经济技术开发区 266555） (Tel: Email: minghai@, 316328693@) 摘要：由于文献中的导热离散方程中热导计算没有考虑导热截面积的连续性变化，因此对强非线性 变截面热传导问题进行数值求解时会产生较大的误差。本文以变系数径向导热问题为研究对象，由 傅里叶导热定律提出了较传统控制容积积分法更为合理的热导计算格式，根据流的相容性原理分别 推导得出极坐标、柱坐标、球坐标下的界面半径和界面导热系数插值公式。通过几个变系数和非线 性导热问题对本文插值方法进行了考核，并与文献中的热导计算格式进行了对比。数值实验结果表 明，由于本文方法充分考虑了变截面和変物性的综合影响，因此对于非线性导热问题能得到更精确 的结果。 关键词：热导；界面半径；界面导热系数；插值方法 0 前言 对于导热方程采用控制容积积分法离散之后，对于每个节点而言其离散方程中邻 点系数具有明显的物理意义，它们的大小反映了周围节点对主节点的温度的影响程度 的大小，而每个邻点系数则表示节点之间导热阻力的倒数（热导）,因此如何有效地计 算热导成为导热问题数值计算是否成功的关键所在。由热导计算公式可以看出其主要 由界面导热系数、界面导热面积因子、节点间距等要素组成，其中节点间距可由已知 的节点坐标直接计算得出。界面导热系数的插值方法没有定论，因此一直成为 CFD/NHT 学术界研究的热点之一。早期的导热问题数值计算中曾广泛采用算术平均 法，该方法等同于线性插值，其优点是简单明了，当导热系数是温度的线性或缓变函 数时，该法是有效可行的，但是当导热系数是温度的强函数或导热系数不是温度的连 续函数时（如复合材料的交界面），该方法就会造成很大误差，目前该法只在文献[1] 中推荐使用。Patankar 教授在其著作[2] 中推荐使用调和平均法，因其物理意义清晰， 尤其适用于复合介质的传热计算，因而成为导热系数界面插值的主流算法。对于某些 可以求得精确解的情形，采用上述2 种方法的数值计算结果表明，即使对有热源或导 热系数呈连续性变化的场合，调和平均也要比算术平均好些[3,4] 。Voller[5] 和 Swaminathan 基于Kirchhoff 变换提出了积分平均插值法，该方法在目前所有插值方法 中精度最高，但由于它需要进行积分运算，工作量大，尤其当导热系数的积分式无法 通过解析式表达时，由数值积分所增加的运算工作量更为明显，因此不适用于大型网 格计算。刘中良[6]教授提出了一种新的待求变量界面插值方法，即先通过线性插值求 出界面温度，再依据导热系数与温度的函数关系式求出界面导热系数，并用此法成功 求解了边界零扩导问题。但是当导热系数不依赖于待求量或空间坐标变化时，该方法 就不再适用。Date[7]也在其著作中指出,采用算术平均和待求变量插值法来处理非线性 资助项目：国家自然科学基金 (No 导热系数虽然可以接受，但两者均无法通过物理中的传导理论得到有效验证，并极力 推荐调和平均法以确保热流通量在界面处的连续性。笔者撰文指出了文献[3-5]中关于 边界界面与内部单元界面导热系数插值格式不一致的问题，并采用调和平均法做出了 修正，修正方法有效地提高了调和平均法计算変物性导热问题的数值精度，并成功求 解了边界零扩导问题[8] 。关于界面导热面积因子的插值文献中广泛使用的方法是直接 采用物理界面处的导热面积因子。实际上，对于等截面导热问题，采用上述插值方法 可以获得令人满意的结果，而对于非线性变截面导热问题上述方法会产生较大的误差， 主要是因为上述方法仅仅对导热系数和导热面积因子分别进行界面插值处理而没有综 合考虑。 本文以变系数径向导热问题为研究对象，综合了导热面积因子和导热系数等界面 参数的耦合作用，由傅里叶导热定律提出了较传统控制容积积分法更为合理的热导计 算格式，根据流的相容性原理分别推导得出极坐标、柱坐标、球坐标下的界面半径和 界面导热系数插值公式。通过几个变系数和非线性导热