自然对流换热高精度非结构网格有限体积法.pdfVIP

第 30卷 第 3期 计 算 物 理 V0l_30．No．3 2013年 5月 CHINESEJOURNALOFCOMPUTATIONALPHYSICS May，2013 文章编号：1001-246X(2013)03-0337-09 自然对流换热的高精度非结构 网格有限体积法 解 岩， 欧阳洁 ， 周 文， 任朝倩 (西北工业大学应用数学系 ，西安 710129) 摘 要：用非结构网格有限体积法求解 自然对流换热时，传统的对流项离散格式难 以兼顾数值精度与计算效率 ，我 们发展了一种耦合高精度格式的延迟修正方法 ，用于对流项的离散．高 e数下方腔驱动流数值计算验证了该方法 具有较高的计算精度和较好的稳定性．Boussinesq流体的自然对流换热数值模拟，表明该方法能有效克服高 o数 时数值计算发散 ，可准确捕捉 自然对流换热问题中不同偏心率下的等温线和流线分布特征． 关键词：非结构；自然对流换热；Boussinesq流体 中图分类号 ：024 文献标识码 ：A 0 引言 自然对流换热是工程技术领域中广泛存在 的一类物理现象，如动力电站封闭母线散热、绕核反应器空气 腔的自然冷却、电子元件的散热等无不涉及到 自然对流换热．这类 问题的特点是封闭腔壁面对流换热与其 内 部孤立物体的对流换热耦合 ，相互影响，因而该问题 的研究有一定难度．1978年，Kuehn与 Goldstein13采用 涡量一流函数法数值计算研究了水平环形域 中的 自然对流换热 ，给 出了环腔 内的温度与流线分布 ，并与 Mach．Zehnder干涉仪的实验测定结果做 了对 比，二者吻合较好．国内陶文铨等人 采用涡量．流函数法对该 问题进行探讨，在三维的流体流动 中并不存在流函数，当遇到三维流体力学问题时，涡量．流函数法难 以实 施 ． 20世纪后半叶，有限体积法 (finitevolumemethod)因其能保证物质的局部质量守恒原理、存贮空间小、 稳定性强、便于推广到三维情形 以及能适应复杂的求解区域等 优点，被广泛应用于流体流动和传热的研 究．Touazi等 基于结构网格采用 CCFV (collocatedclusteredfinitevolume)格式对方腔 自然对流换热进行了 研究．2011年，Behrooz等人 同样基于结构网格采用人工压缩法数值模拟 了自然对流换热．但是，许多 自然 对流换热问题 中，求解区域不规则、甚至十分复杂，传统的结构网格难 以准确地对其进行剖分，因此发展非结 构网格下的数值求解器，增强对复杂区域的适应性势在必行．然而由于非结构网格本身拓扑关系 比较复杂， 使得结构网格上的高精度格式推广到非结构网格上时非常繁琐，有时甚至难 以实现．基于上述原因，采用高 精度非结构网格有 限体积法求解 自然对流换热问题非常必要． 采用非结构 网格有限体积法求解对流扩散问题时，传统的对流项离散格式难 以兼顾数值精度与计算效 率．我们在文献 [7]研究的基础上 ，发展一种耦合 Q-QUICK (quasi—quadraticupwindinterpolationofeonvective kinetics)格式的高精度延迟修正方法离散对流项 ，该格式简称为延迟修正 Q—QUICK格式 ，即DCQ—QUICK (deferredcorrectionQ—QUICK)格式．在 DCQ-QUICK中，上风侧远邻点处的速度采用二阶迎风插值计算，避免 了传统加权平均方法导致的精度损失，而且，由于DCQ—QUICK中高阶Q．QUICK格式对离散方程所做贡献体 现在源项中，因而在保证数值结果具有较高计算精度的同时，大大降低 了方程求解 的复杂度．为避免运用交 错网格繁琐、复杂的缺点，本文对流场控制方程的整体求解采用了基于 SIMPLEC算法 的非结构 同位网格 收稿日期 ：2012—09—06；修回 日期：2012—12—19 基金项 目：国家重点基础研究发展计划(973)(2012CB025903)资助项 目 作者简介 ：解岩 (1984一)，女 ，安徽 ，硕士生 ，从事有限体积法研究及复杂流体数值模拟 十通讯