轴对称问题贴体坐标系中控制方程的有限容积法地研究.pdfVIP

中国工程热物理学会 流体机械 2004年学术会议 编号；047024 轴对称问题贴体坐标系中控制方程 的有限容积法的研究 曾敏陶文铨李增耀 西安交通大学能源与动力工程学院动力工程多相流国家重点实验室西安710049 Tel：029 E-mailwqtao@,mail．xjtu．edu．cn 摘要：圆柱轴对称问题在很多工程技术中经常会遇到的．本文采用有限容积法，通过张量分析导出了 圆柱曲线坐标系下控制方程的离教形式，并对所开发的程序进行了考核和验证． 关键词： 轴对称，贴体坐标，有限容积法。张量分析 1引言 在大多数实际的工程问题中，流动和换热过程都是在复杂的几何形状内发生的，对 于复杂区域内的流动，贴体坐标的方法得到广泛的应用。但是对于轴对称问题采用贴体 坐标方法的计算在公开的文献中不多见，本文采用CFD领域中广为采用的有限容积法， 开发出了适合计算各种轴对称问题的贴体坐标程序。 2坐标变换及速度的选取 z，{ 图l基向量 图2圆柱曲线坐标系 在曲线坐标系中(图1)，存在如下两个基向量【1l： 对于图2所示圆柱曲线系，我们可以获得它的协变基向量为： · OZ艺．防： 出t 防t 叶 ： e(05琵“琵，，气z)2瓦¨面7铴一丌 61 有了这个协变基向量后，其它的几何变量值如：；‘”、岛、矿、彰、j(J’和‘，等可以 用张量的形式很简洁的表达出来【l】． 如图3所示，，点的速度矢量可以按如下三 种方式表示： 孑=，r=叫圳=扩爵。 (3) 在计算平面中，作为求解变量的速度分量有三种 选择。根据文献【2】的分析，采用逆变速度作为计 算平面上连续性方程的速度分量，而以笛卡儿速 度作为计算平面上动量方程的速度分量是一种较 图3速度分量 为合适的选择，本文采用的也是这种组合。 l、笛卡儿速度2、逆变速度3、协变速度 3方程的离散 对于轴对称问题，取一个弧度作为圆周方向的计算区域。为简化起见，以其一个截 面来表示，如图4所示。采用有限容积方法，在图4所示的控制容积中对如下的控制方 程进行离散，当对流项采用不同的格式来处理时，会得出不同的离散方程。采用规正变 V·(p旷妒一I’．V妒)=墨 (4) 卜一(A0——叫 I。!盟．一Ij组l 卜—(碍卜一 图4曲线坐标系中控制容积 图5计算平面中的几何因子 量的方法“1可以得到最终统一形式的离散结果如下： 4办=∑4。丸+品∥+＆+瓦。 其中： 4=4+4+4+鸣一·写以=咒一《+嚣一鼍