扩散系数界面插值方法的比较研究.pdf

中国工程热物理学会 传热传质学 学术会议论文 编号：093013 扩散系数界面插值方法的比较 王辉，徐明海 （中国石油大学（华东）储运与建筑工程学院，山东青岛经济技术开发区 266555） (Tel: Email: minghai@, 316328693@) 摘要：在有限容积法的基础上，针对计算区域内扩散系数变化的情形，分别采用内节点法和外节点 法 2种区域离散化方法，通过数值求解几类典型的変物性导热问题，考察了算术平均、调和平均、 待求变量插值以及积分平均等4种扩散系数界面插值方法的计算特性。结果表明，对于不同的変物 性问题，采用不同的空间离散方法时，各插值方法计算特性不尽相同，并没有综合性能最优的界面 插值方法，应该视具体情况分别采用相应的处理方法。 关键词：界面插值；扩散系数；区域离散；有限容积法 0 前言 在实际的工程计算中，经常会遇到扩散系数变化的情形，如变物性导热、湍流流 动、多相渗流以及原油流变性等问题。使用有限容积法对描述以上工程问题的控制方 程进行离散后，方程系数中出现了界面扩散系数项。由于扩散系数在界面位置事先没 有定义，因此必须通过节点处的值进行插值处理。NHT学术界经常采用的界面插值方 法主要包括：算术平均、调和平均、待求变量界面插值、积分平均等。算术平均法[1,2] 相当于线性插值，由于处理简单，易于编程，该方法在早期的导热问题数值计算中曾 [3] [4,5] 广为采用 。Patankar教授在 1978年基于热阻串联原理提出了调和平均法 ，因其 物理意义清晰，尤其适用于复合介质的传热计算，因而成为导热系数界面插值的主流 [6] 算法。Chung和Payne通过进一步研究，指出把物性阶跃面设置在节点位置处较之设 在界面位置处时其热流密度计算结果更准确。笔者通过理论分析指出了文献[3-5]中关 于边界界面与内部单元界面导热系数插值格式不一致的问题，并采用调和平均法做出 了修正，算例表明修正方法有效地提高了调和平均法计算変物性导热问题的数值精度 [7] [8,9] 。Voller和Swaminathan基于Kirchhoff变换提出了积分平均插值法，该方法是 目前所有插值方法中精度最高的一种，由于它需要进行积分运算，工作量大，尤其当 扩散系数的积分式无法通过解析式表达时，则需要进行数值积分，由此所增加的运算 [10] 工作量更为明显，因此不适用于大型网格计算。刘中良教授 提出了一种新的待求变 量界面插值方法，即先通过线性插值求出界面温度，再依据导热系数与温度的函数关 系式求出界面导热系数，因其与积分平均法精度相当且实施方便也倍受学界关注，但 是该法同样要求扩散系数具有解析表达式，因此当扩散系数不依赖于待求量或空间坐 标变化时，该方法就不再适用。 以上4种界面插值方法各有优劣，对于不同的实际问题，它们各自具有不同的适 用范围，尤其对于不同的空间离散格式其计算精度也不尽相同。本文首先对4种界面 插值方法的原理及计算公式做了简要分析并进行了拓展。然后分别采用内节点法和外 资助项目：国家自然科学基金 (No 节点法2种空间离散格式，通过几个具有精确解的変物性扩散问题对各种插值方法进 行了考核，分析比较了各种插值方法的计算结果的精度。最后提出了适于工程数值计 算有效建议。 1界面插值方法 在计算离散方程的系数时，需要用到界面扩散系数，通过节点处的扩散系数进行 界面插值处理是最为有效的办法。参照图 1所示的等截面一维空间区域离散示意图， 用 和 分别表示节点W和P处的扩散系数，单元界面w处的界面扩散系数用 表 W P