轴流轩轮机三维粘性非定常流动计算的双时间方法.pdfVIP

中国工程热物理学会 热机气动热力学学术会议 编号 轴流叶轮机三维粘性非定赏流动计算的双时间方法 刘前智 周新海 胡运聪 西北工业大学七系(西安71∞72) (柞者联系电话：02争8495224，Email：liuqz‰霄pu．edu．cn) 摘要 3 本文以双时间方法计算了单级跨音速压气机的三维非定常粘丝流堑，应用LU-sGs 隐式解法快速获得虚拟时间域的定常解，有效地减少了计算时间．对某跨音速单级压气 ● 机同时进行了双时间方法和显式方法的非定常流数值计算，两种方法计算结果是一致 ● 的，但从计算时问上来说．双时{可法与LU．sGs隐式解法相结合的非定常粘流计算模式 ● 具有明显的优点． 一、前言 叶轮机内部由于转动叶捧与静止叶捧的相邻捧列，相互之间产生了强烈的非定常干 扰，探索这种固有现象的物理本质是促进叶轮机设计技术发展的重要途径。随着计算流 体力掌的不断发屁特别是计算机技术的进步，为数值模拟叶轮机内部的三维非定常流 动提供了必要豹条件．已有的研究表明，叶片排之间的非定常干扰对叶轮机气动性能有 较大影嗨，而对叶片所受到的气动力影晌就更太。有效地预测叶轮机中动静时捧相干的 非定常影响是改进目前的定常流设计方法和设计体系的一个重要方向，特别是在高性能 多缀跨音速风扇、压气机的设计中，叶片气动弹性这一与非定常效应密切相关的闯劂导 会更加突出，迫切要求发展实用性强的数值分析手段，来处理叶轮机非定常流问题。 由于叶轮机三维非定常粘性流动的计算量巨大，当前所面临的主要问题是如何快速 获得满足周期性条件的非定常解，以便于工程实际的应用。在以往的定常流计算时，由 于所要求的是非时间精确解，所以在显式方法的时间推进求解过程中，可以采用诸如当 地时问步长、多网格法等措施来有效地加速收敛，而这些措施对于非定常时间精确显式 方法的求解中则不能加以应用。另一种选择，则是采用隐式格式使时间步长限制大大放 ● 宽。但经过对角化、因式分解等步骤，会降低时间精度，特别是对于粘性非定常流，不 适当的简化将使时间精度降至无法接受的地步。然而若不作上述处理，隐式计算的复杂 t 程度及计算工作量都将大大增加。九十年代初，JⅫes∞flj提出的双时闯方法(nIal Tin壕M酏od)在非定常流动计算中具有显著的优越性，得到了广泛的重视口～。该方 法在橱理时间域采用时间=阶精度隐式格式进行时间离教，在每一时间层上引入一个虚 拟时间域，通过在虚拟时间域上求解一个定常解来得到物理时间域的时间二阶精度隐式 解，这样，既避免了隐式直接求解的巨大计算工作量，又可以在虚拟时间域的定常流场 求解过程中采用各种加速收敛措施，减少整个过程的计算时间．应该指出的是，在该方 法中影响整个过程计算时间的因素主要是虚拟时间域定常流场求解的收敛速度，如果定 常流计算不采用十分有效的加速收敛措施，则无法体现该方法的优越性。 本文采用双时间方法计算叶轮机三维非定常粘性漉动。在虚拟时间域采用LIJ．SGs 3硒 1t | 誉—■—_彳—『翻r—————————————————————一 隐