圆柱和直桨叶突然启动瞬态流动的数值研讨.pdfVIP

中固工程热物理学会 流体机械 2006年学术会议论文 编号：067004 圆柱和直桨叶突然启动瞬态流动 的数值研究 李志峰吴大转王乐勤 10027杭州) (浙江大学化工机械研究所3 Tel：0571 Email：hjwlq2@zju．edu．cn 摘要：采用基于动网格方法的有限体积法对圆柱和直桨叶突然启动引起的二维非定常不可压粘性流 进行了数值模拟，给出了计算方法。通过计算得到了Re=5000和9500时圆柱突然启动后的流场分布： 分析了采用不同湍流模式计算高雷诺数下的圆柱启动流动，并对结果进行比较分析：采用大涡模拟 对高雷诺数下的直桨叶突然旋转启动的流场进行了计算，得到了完成启动后叶片内部流场随时问的 演化过程。所采用的方法和结果为进一步求解三维叶轮启动过程内部非定常流动提供了基础和依据。 关键词：瞬态流动：数值模拟；有限体积法：动网格：突然启动 1 引 言 在白然界和工程实际问题中普遍存在的1F定常流动问题，通常情况下它可以分为 两种类型，第一种类型为边界移动引起流场的非定常变化，在工程上例如流体机械的 启动和停机过剧h羽，第二种类型是由于流体自身诱发的周期性非定常流动。其中圆 柱瞬时启动和流体机械叶轮突然加速旋转过程引起的非定常流动属于第一种非定常流 动现象，此类问题的研究对于理论分析、实验测试以及解决工程问题都具有十分重要 的意义，特别是对于许多存在突然启动现象的工程问题。 对于第一类非定常流动问题，在实验研究方面，Bouard等u’对圆柱的突然启动过 程进行了实验，得出40Re10‘范围内分离区内详细的流场分布情况，已经成为验证 数值计算结果的一个重要参考。然而对于叶轮突然加速旋转问题，目前还没有较系统 的实验研究结果。在突然启动问题的数值计算方面，目前采用大多采用差分法和涡方法 14-61，这两类方法在求解几何形状简单的问题具有较好效果，但在处理复杂结构形状方 面存在边界条件处理困难等问题，并且向三维复杂问题的拓展方面存在较大困难。 由于流体机械加速旋转引起的流动问题比较复杂，采用以上方法目前还存在许多 的困难，例如边界条件的处理以及三维的求解问题等。而有限体积法在处理复杂流动 问题方面具有明显的优势，由丁其在流体机械内部流动分析等复杂工程问题上取得很 好的效果，同时，考虑到采用动网格方法H’使得启动过程移动边界的处理更加灵活，本 文采用基于动网格方法的有限体积法对平动和旋转两类突然启动引起的流场进行数值 求解，在验证基于动网格方法的有限体积法求解突然启动问题的有效性基础上，为求 解流体机械启动过程的内部流动打F基础。 幕金项日：国家自然科学基金项目(No l 2控制方程和数值方法 2．1控制方程 ’ 考虑由移动边界洲所包围的任意控制体积内积分形式的非定常不可压黏性流的 体积守恒方程、质量守恒方程与动量守恒方程： 旦dt fd矿一￡，u。积=o’ ‘D 昙工∥矿+L加．u。)删-0， (2) ∞’ 式中，t为时间，y为控制体积的体积，U。为控制体边界的移动速度矢量，n为边界 彳上外法向单位矢量，p为密度，U为速度矢量，∥为动力黏度，墨为源项。 2．2数值求解方法 上述方程中的时间导数项用一阶向后差分格式表示为 竺fudV：坐芝尘堕． (4) dl毋