非线性导热方程高精度离散格式.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：093011 非线性导热方程高精度离散格式 王辉，徐明海 （中国石油大学（华东）储运与建筑工程学院，山东青岛经济技术开发区 266555） (Tel: Email: minghai@, 316328693@) 摘要：对于常系数导热方程采用传统的控制容积积分法离散就可以获得较高精度的数值解，而对于变 系数和非线性导热方程则会产生较大的误差，为获得满意的结果需要加密网格，因此会大量消耗存储 空间和运算时间从而增加计算成本。针对上述问题，本文基于微元体平衡法并结合控制容积积分法， 由能量守恒定律重新推导了关于节点温度的差分方程，给出了离散方程系数计算公式。通过几个代表 性的算例对本文方法进行了考核并与文献离散方法的计算结果进行了对比。数值试验结果表明采用本 文方法在较少的网格数下就可以获得较高的精度。 关键词：导热方程；非线性；变系数；元体平衡法；控制容积积分法 0 前言 导热现象在实际工程中随处可见，例如埋地管道原油停输温降、金属材料铸造加工 过程中的温度分布、食品的冷冻过程以及注蒸汽井筒周围地层温度分布等，以上导热现 象通过建立物理模型和数学推导可以统一地由导热微分方程来描述。对于简单的导热方 程，可以直接通过数学求解得到显式解析解。而具体工程问题需要考虑的因素很多，其 抽象建立的模型就较为复杂，此时建立的导热方程就难以直接求解，而在实验条件不允 许的情况下，需要用数值方法通过计算机编程进行求解。目前流动与传热数值计算领域 内主流方法首推有限容积法，该方法将计算区域分成一系列控制容积，每个控制容积以 节点作为代表，通过将守恒型的控制方程对控制容积进行积分从而导出离散方程，利用 该方法导出的离散方程可以保证守恒性且离散方程的系数物理意义明确和易于理解。对 于几何结构简单、常系数导热方程，采用传统的控制容积积分离散就可以获得较高精度 的数值解。而实际上导热材料的热物理性质（如导热系数、密度、比热）是随温度或者 空间几何坐标变化的，且导热截面积也可能并非常数而是依赖坐标而变化。诸如此类非 线性及变系数导热方程若采用传统的控制容积积分法离散方程可能会产生较大的误差， 因此为了获得令人满意的结果不得不加密网格，但是这样做会消耗大量的内存空间和 CPU 运算时间从而导致计算成本增加。为此人们一直致力于开发一种新的算法力求用最 少的节点获得最高精度的解。对于导热系数依赖温度变化的情形，Patankar 教授在其著 [1] 作 中推荐采用调和平均法来处理变导热系数，该方法由热阻串联原理而得出，尤其对 于分层介质适应性较好，但对于强非线性导热问题则精度较差。Voller[2,3]和Swaminathan 基于Kirchhoff 变换发展出了界面导热系数的积分平均插值法方法，并运用该法成功求解 了相变传热问题，从而大大提高了変物性导热问题的数值计算精度，但该方法只涉及了 资助项目：国家自然科学基金 (No 物性依赖温度变化的情形，而对于物性依赖坐标变化和变截面导热问题并未充分考虑。 笔者从傅里叶导热定律出发同样推导得出了上述界面导热系数的积分平均插值公 式，并在此基础上进行了扩展，采用微元体平衡法，依据能量守恒定律，对傅里叶定律 表达式分离变量并积分，综合考虑了変物性和变截面等非线性因素，由此构造了关于节 点温度的新的差分方程，给出了离散方程的详细计算公式，通过几例有精确解的一维和 二维空间非线性及变系数导热问题对本文离散方法进行了考核并与传统的文献离散方法 的计算结果进行了对比，结果表明采用本文方法只需少数几个节点就可以获得较高精度 的解，从而为导热问题的高精度数值计算开辟了一条新的途径。 1 离散方程的建立 元体平衡法的基本思想是：把物理上的守恒定律直接应用于微元控制体，并把节点 看做控制容积的代表，从而导出节点上待求变量之间的代数关系式，可以将它看做控制 容积积分法的一种变形与补充。下面主要讨论一维、二维空间下的导热问题离散方程的 建立，对于三维空间导热同样可以采用本文方法建立离散方程，在此不做讨论。 1.1一