基于SIMPLE的紧致方法求解对流换热问题研究.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：083447 基于SIMPLE 的紧致方法求解对流换热 问题 张昆，杨茉，王治云，郑建城 （上海理工大学动力学院，上海， 200093 ） (Tel Email:zhangkun52015@) 摘要：通过延迟修正的方法建立了基于 SIMPLE 的有限容积紧致方法求解对流换热问题，对底部加热 的方腔内自然对流换热问题进行数值模拟，证明了方法的可行性和高精度。同时也说明了不同网格系 统对数值解精度的影响。 关键词：有限容积紧致格式；SIMPLE 算法；对流换热； 非等距网格。 0 前言 在自然界中，许多真实的物理流场都具有多种的空间尺度和时间尺度，有的还具有 流动的非定常性，非稳定性和耗散性。对于这些复杂流场很难用低精度的数值方法获得 满意的数值解[1] 。紧致格式是一类高精度的格式，能在更宽的谐波范围内逼近真实解的 谐波成份，更好地反映流动的小尺度结构。紧致格式的特点在于结点的导数采用隐式求 解，使用尽可能少的网格点达到更高的精度。虽然网格基架点少，精度高，但是任意结 点的导数值不仅与相邻结点有关，而且与相邻结点的导数值有关，需要求解有关结点导 数的代数方程组才能获得。所以求解任意结点的导数都涉及整个流场或一个方向上所有 结点的信息，需要偏微分方程的离散采用显式格式。若采用隐式格式，直接将紧致格式 代入偏微分方程中进行离散，则每一个离散方程中的未知数均涉及所有结点，离散方程 组相当繁琐，当网格结点数目较多时，甚至是不可能实现的。 近年来，如何构造基于各种算法的求解对流换热问题的紧致方法的研究取得了很大 的进展。Ming Li, Tao Tang 建立了基于涡量流函数的紧致方法求解对流换热问题 [2](1995)，这种方法涡量边界条件的确定是比较困难的，不像在原始变量法中确定速度 与压力的边界条件那样可以由物理概念上的考虑而得出明确的结果[3] 。K.K.Q.Zhang, B.Shotorban等建立了基于人工可压缩性法的紧致方法求解对流换热问题[4](2006)，人工 可压缩性系数的选择对这种方法的影响很大，不合适的可压缩性系数得不到收敛的解。 SIMPLE算法是一种发展比较成熟的求解对流换热问题的方法，但由于SIMPLE算法 在时间格式上可以采用隐式或半隐式方法，所以紧致格式不能像经典格式那样直接用于 SIMPLE算法中控制方程的离散。1999 年，杨茉等通过延迟修正方法将QUICK格式引入 SIMPLE算法中获得成功[5]。受到这种方法的启发，本文通过延迟修正方法建立了一种基 于SIMPLE的有限容积紧致方法求解对流换热问题。用该方法计算底部加热的方腔 自然对流换热问题，数值模拟结果证明了该方法的可行性和高精度。 资助项目：国家自然科学基金项目(No.；上海市重点学科建设项目资助(No. J50501) 1 数值方法 1.1 有限容积紧致格式 以模型方程为例说明基于非等距网格的有限容积紧致格式，建立如图 1 所示的网格 系统。