多级轴流压气机湍流和二次流径向掺混的简单统一模型及其应用

本文构造出一套新的统一表示湍流和二次流径向掺混作用的多级轴流压气机子午通流流场的模型方程组。除了掺混系数外,新的方程组没有其它的未知关联项,使得含径向掺混的多级轴流压气机子午通流的计算简单得多,相应计算机程序的计算时间只是无粘计算机时的约1.5倍。新方程组的数值结果与试验结果以及与则Gallimore掺混理论结果比较,表明新方程组的合理性和优越性。新方程组有严格的理论基础,不仅能用于亚音压气机,还可能用于跨音和超音压气机。      

轴流压气机考虑展向掺混的数值计算

为了更加准确地反映压气机内部流场的展向分布特性并提高计算精度,基于流线曲率法发展了1种考虑压气机内部展向掺混作用的计算方法。应用Gallimore-Cumpst y(G-C)掺混模型模拟紊流扩散的作用,在2维Cart er公式基础上引入改善的Rober t s修正模型对落后角进行预估,通过落后角反映二次流的作用,并对3级轴流压气机3S1进行了计算分析和比较。结果表明:不考虑掺混因素时所得计算数据的展向分布情况与实际流场相差较大;只有同时考虑紊流掺混和二次流2种因素,才能得到最优结果。该方法能够较好地模拟轴流压气机的展向掺混作用,得到更加贴近实际的气动参数分布,可为压气机设计和优化提供参考。     

多级轴流压气机的径向掺混

 随着压气机向着更高的级负荷和更低的展弦比方向发展,径向掺混问题在气动分析和气动设计中越来越重要,当前对决定径向掺混的主要因素众说不一。本文用严格的理论,得到了用统一形式表示含二次流和湍流在内的各个径向掺混因素的多级轴流压气机通流基本方程组,并讨论了二维通流流场湍流掺混与实际流场三维剪切结构的关系。基本方程组的数值结果与试验结果以及其它掺混理论结果比较表呀,该方程组能够满意地模拟多级轴流压气机中的径向掺混特性。     

考虑径向掺混的轴流压气机通流彻体力模型计算及验证

为研究考虑径向掺混对通流计算结果的影响,在通流彻体力模型中引入径向掺混项,并以某型低速轴流压气机为算例开展了有无径向掺混对通流计算结果的影响分析。结果表明:在落后角和损失系数相同时,径向掺混模型的加入使叶排出口气流参数沿径向的分布更加均匀;模型对于压气机后面级端壁区流动参数展向分布的预估更加准确。通过与该台压气机试验测量数据的对比,验证了本模型的可靠性。     

跨声压气机内确定应力耦合问题研究

基于三维定常Denton程序发展了二维定常、二维非定常、三维非定常计算程序,以二维跨声压气机叶栅、三维跨声单级压气机近失速点为算例,进行确定应力耦合问题的研究.结果表明:在确定应力应用到时均方程的耦合处理中存在一些因素导致时均方程的解和非定常解的时均值有一定偏差,它们分别是空间差分离散引入的误差、黏性相关项引入的误差、掺混面处理引入的误差.在多排耦合问题中掺混面的处理方式至关重要,对耦合计算结果有较大的影响,简单的掺混面处理并不能获得理想的效果.      

跨声轴流压气机转子尾迹的总压平均方法

对一台跨声轴流压气机转子的尾迹流场进行了数值模拟研究和总压实验测量,结合理论分析方法,研究了面积平均和质量平均这两种常用的总压平均方法在跨声轴流压气机转子尾迹测量中的联系与差异.结果表明:在靠近转子尾缘的测量截面,由于转子尾迹流场中存在总压高而轴向速度低的区域,导致面积平均总压大于质量平均总压;而随着尾迹与主流在下游逐渐掺混均匀,面积平均总压又小于质量平均总压.两种总压平均方法的这一规律性关联可以在研究压气机转子尾迹特征尤其是进行数值计算和实验测量结果对比分析时提供一些有价值的参考.      

跨声速轴流压气机三维黏性流场的数值模拟

改进了动静叶级间的掺混面处理方法,提出了用参数分布修正方法来保证参数在掺混面上的周向不均匀性,并使用亚松弛的参数分布系数来解决程序在动静叶轴向间距较小时的计算稳定性和收敛性问题。采用经改进的NAPA软件对NASA37级压气机的设计点和在不同转速下的流场进行了数值模拟,计算结果与试验结果吻合较好。     

考虑端区黏性和径向掺混的流线曲率通流算法研究

高负荷多级压气机的技术难点之一在于如何快速、完整、准确地获得全工况下叶片前尾缘气动参数分布和整机特性结果,而通流计算凭借速度快、指向性明确的优点,一直是多级压气机设计的核心技术。为了提高通流代码的适用范围和预测精度,推导了统一适用于轴流、斜流、离心压气机的流线曲率方程,引入端区黏性和径向掺混模型,应用Fortran语言针对轴流压气机全新开发了其通流特性计算程序,对NASA Rotor37、Rotor67转叶和GY1-2J压气机进行了计算,计算结果表明:新模型的加入提高了程序的计算精度,初步验证了程序的有效性。并提出提高程序的通用性是今后通流计算的发展方向。     

进气畸变对压气机性能影响的三维彻体力模型

基于三维彻体力模型理论初步发展了一个能有效预测进气畸变下压气机性能的压气机三维稳定性分析模型(CSAC).该模型通过求解带有源项的可压三维Euler方程组,对压气机内部全环三维流场进行模拟.叶片对气流的作用力和做功源项根据压气机三维稳态Navier-Stokes(N-S)解的相关气动参数计算得到,并给出了源项与叶片进口气动参数之间的数据库关联方法.利用该模型计算和分析了NASA Rotor 37在均匀进气下的气动性能,并与三维稳态N-S计算结果对比,验证了该模型的准确性.最后利用该模型模拟了Rotor 37在稳态周向总压畸变下的性能.结果表明:总压畸变不但降低了压气机的气动性能和稳定裕度,而且在压气机转子出口诱发总温畸变,研究表明该模型在消耗较少的计算资源下正确地反映了进气畸变对压气机性能的影响,是目前分析进气畸变对压气机性能影响的有力工具.      

多级轴流压气机非设计性能的数值预估

采用LU-SGS隐式解法以及叶排间掺混面模型数值模拟压气机三维粘性流动,研究了多级轴流压气机的设计与非设计状态性能的预估问题。对一台五级高压压气机进行了详细的流场数值计算,包括了设计转速的100%,90%及80%等多个转速及各种工作状态。结果表明,本文所建立的数值方法与计算机程序对多级轴流压气机各种非设计状态的计算都具有收敛容易和快速的优点,计算得到的压气机总性能和各分级性能与设计参数进行了对比,证明该方法和程序的预估结果是比较准确的。根据设计转速下的计算结果,预估了压气机设计点的喘振裕度。此外,还对不同状态下叶片进口气流角的变化进行了分析。      

超声速气流中轴向涡强化混合及燃烧的数值研究

通过求解N－S方程，数值计算了扩张圆管内的非平衡流动，并研究了超声速混合和燃烧的强化技术。结果表明，在超声速领域，预旋燃料方式是一种比较有效的强化技术。     

燃烧室流动的数值计算

提出了多连通域流场统一计算的处理方法,采用三维贴体曲线网格,非交错网格法和k-ε湍流模型。对包括扩压器,内外环通道和火焰筒的整个燃烧室流场进行了数值计算。可以计算流过主燃孔、掺混孔和气膜孔的流量。计算结果定性合理。      

大速差射流燃烧室内三维湍流气固两相流的实验研究

本文用二维激光多普勒测速仪和探针取样法对大速差射流燃烧室中冷态三维湍流气固两相流动进行了实验研究,测量了气相速度,颗粒速度、颗粒质量流分布及两相各自的湍流特性,探讨了高速射流作用下两相流场的主要物理特征及其对火焰的稳定和强化燃烧的作用。     

应用TVD格式分离求解不可压N-S流动

提出一种新的将连续方程与动量方程分离求解的TVD引用方法。对动量方程中的对流项采用三阶精度、迎风TVD格式离散, 粘性项和压力项采用中心差分格式离散。压力采用压力替代法求解。计算所用的网格是具有交错网格主要特点的修正非交错网格, 并采用一种压力修正法和G-S中心对称迭代法加速收敛。给出了二维管流和二维空穴流动在不同状态下的计算结果, 并与其它解法的计算结果或精确解进行比较。结果表明, TVD格式不仅能用于不可压流数值计算, 而且还能得到比非TVD差分法精度更高的结果。      

积分方程法计算翼型的跨音速绕流

从跨音速小扰动方程出发推导积分方程的过程中,本文用任意形状的封闭曲线C_Q（其极限趋于零）挖去奇点Q,最终得到无奇性（指无穷奇性,不包括Cauchy奇性）的积分方程。 对于跨音速流中的圆头翼型的前缘问题,提出了一种解决办法。 证明了Nixon给出的反演公式对于超临界有激波的小扰动流动也成立。 关于积分方程法中的人工粘性方法,对Sachdev和Lobo提出的方法做了改进。 最后给出了NACA0012翼型在有无升力和有无激波各种状态下的计算结果。比较表明,本方法的计算结果与其它方法的计算结果符合得较好,且计算量很小。     

内外函混合加力燃烧室扩压器内液雾湍流两相流场和浓度场数值模拟

 本文用κ-ε湍流模型和液雾轨道模型,对二维内外函混合加力燃烧室扩压器内有蒸发的液雾两相流动和混合进行了数值预报。其中对气相速度场和温度场的计算结果与实验符合令人满意,液雾两相流的计算结果在趋势上也是合理的,表明本文的模拟方法可用于加力燃烧室性能估算及优化设计。     

时间推进积分方程法求解非定常跨声速流动

本文提出了一种求解非定常跨声速流动的新方法——时间推进积分方程法,此法克服了时间线化积分方程法的限制,能较好地模拟激波的运动。本文首先用一维波(?)问题——模型问题阐明此法的基本思想,然后将它应用于二维非定常跨声速流动中。本文还首次引入拟速度位的概念,使时间推进积分方程式得到简化,尾涡条件和Kutta条件更易处理。数值计算表明时间推进积分方程法是合理可靠的。     

混合流场控制的大涡模拟

采用大涡模拟方法研究了三维小扰动下,混合流场大尺度拟序结构的产生和演化过程,捕捉了展向涡的卷起、配对、合并,以及二次流向涡的出现等大尺度的三维拟序结构,分析了拟序结构与入口扰动方式之间的内在联系,再现了涡卷自身撕裂而引发转捩的现象。