大迎角下俯仰振荡对细长体气动特性的影响

通过对头部无扰动和有扰动的两个细长体模型进行数值计算,分析了各个流场中涡系结构和气动力变化历程,并对比两者之间的差别.发现在高频率小振幅的俯仰振荡作用下,头部无扰动和有扰动的细长体的流场随时间的发展慢慢由非对称转变为对称.相对于头部有扰动的细长体,在第一个振荡周期内非定常俯仰运动对头部无扰动的细长体作用更明显.     

气泡稳定性问题

在一维球对称解的基础上,研究了气泡表面二维扰动发展的稳定性问题。用WKB方法给出了扰动所满足的微分方程的一致有效渐近解,描绘了不同区域中扰动发展的行为。缝合相应的渐近展开,得到了在坍缩和膨胀阶段中扰动发展的整体解,讨论了扰动发展的完整的图象。对有关结果进行推广,讨论了流体的粘性和可压缩性对气泡稳定性的影响。     

关于跨声速流动问题的多重解

本文讨论了跨声速流小扰动方程Riemann问题的多重解及物理上可实现的解的判别。     

中心分级燃烧室雾化特性试验研究

为了研究中心分级燃烧室雾化特性,设计中心分级燃烧室头部进行雾化特性试验研究。头部预燃级采用贫油直接喷射,主燃级采用预混预蒸发。试验采用了相位多普勒粒子分析仪测量液滴粒径及速度,用10μm以下的小液滴速度近似流场速度。实验结果表明:流场具有中心回流区(PRZ)、唇口回流区、角落回流区(CRZ)、预燃级高速射流区及主燃级高速射流区等结构;流场结构对称且受工况改变影响较小;大液滴集中在中心回流区及角落回流区。中心分级燃烧室不同于其他分级燃烧室,其具有特殊的流场结构及燃油雾化分布规律;在近场流场区域65μm～75μm液滴集中在PRZ区域和CRZ区域;近场区域内Case 4雾化D32是Case 3的60%。     

流向强迫作用下的液体初始雾化机制及动力学特征

使用Volume of fluid (VOF)方法和基于树形数据结构的自适应算法来研究射流雾化的破碎过程以及扰动对射流破碎机理产生的影响。在无扰动情况下,液体射流的头部、液丝和液滴随着射流时间的发展不断发生演变。射流头部先呈现蘑菇状外形,随后液丝生成,并慢慢转变成网兜状,直至断裂形成小液滴。在周期性流向强迫的作用下,射流液柱的表面会形成周期波,其液丝破裂形成液滴的时机与稳定射流情形相比会有所提前,射流形成的头部更趋于扁平,最终生成的液滴数量更多。低中频阶段,随着扰动频率的增大,射流未扰动液柱长度(L)逐渐缩短,液滴直径的概率密度分布(PDF)趋于尖锐,液滴平均直径(SMD)增大。在高频阶段,随着扰动频率的增大,L会随之增大,液滴直径的PDF分布变得平缓,SMD会减小。     

跨声速压气机转子旋转失速现象的非定常模拟

为研究跨声速压气机转子失速机理,全周非定常数值模拟了某跨声速压气机单转子的失稳过程。结果表明:该转子由叶尖Spike扰动诱发旋转失速。在小流量稳定工作状态,压气机转子叶尖区域存在"旋转不稳定"(Rotating Instability,RI)流动现象。压气机节流过程中,转子进出口的流量降低,叶尖区流场非定常波动幅值增大。近失速状态时,RI扰动团的典型流场结构"径向涡"在叶尖区域形成堵塞,导致相邻叶片前缘间歇性地出现溢流现象。随着压气机进一步节流,转子叶尖的负荷达到极值,叶片通道尾缘逆压力梯度过大,出现倒流。尾缘倒流的出现又进一步增加通道内的堵塞,最终形成Spike扰动。失速先兆对应的流场结构是沿叶片前缘额线向相邻叶片压力面周向运动的"径向涡"结构。     

跨、超音速流动的区域分解方法与并行算法

研究二维跨、超音速无粘性流动的 Euler方程区域分解方法、并行算法及其应用。通过内边界耦合条件实现相邻子区域解的光滑过渡,以得到总体流场的数值解。发展了一种多块区域之间守恒型的有效内边界耦合方法,对二维翼型跨音速流动和钝头体超音速流动等进行了分区数值求解,分区计算结果与其他单区计算结果作了比较,并讨论了多种区域分解数目的分区计算效率。并行计算采用纯结点并行编程方式和“先进先出”的同步控制等待机制,利用 PVM并行环境对二维绕翼型跨音速流动做了二区和四区分区并行计算。     

头部和后体对钝头体侧向力的影响

钝头体大迎角飞行时会出现随机的非对称流动现象,引起不确定的较大侧向力,进而使其偏离运行轨道。通过在钝头体头部施加人工扰动块可以固定其大迎角下的非对称流场结构,得到确定的侧向力,以利于改善钝头体的大迎角飞行特性及机动性。本文讨论了在头部人工扰动块主控流场结构的基础上,模型后体对侧向力影响的存在性问题,在迎角为50°、雷诺数为1.54×105的条件下,利用实验对周向角为90°和270°、子午角为10°的扰动位置的球形扰动主控下的侧向力影响因素进行了研究。发现钝头体大迎角下的非对称流动结构在头部主控的基础上,后体对非对称流动的影响不会消失,且其为影响头部扰动主控作用的重要因素。尽管模型后体的影响不会改变钝头体头部对于流场结构的主控地位,但会影响头部扰动控制的精准程度。所以在通过钝头体头部施加扰动进而得到确定的侧向力的同时,还需要减小模型后体对流场的影响,对其结构和加工质量进行优化,以更好地通过人工扰动主控流场结构。     

叶根失速先兆触发跨声速压气机失速的机制研究

实验中发现了一种起始于叶根的新型压气机失速先兆,为了探索这种先兆的初始扰动特征及失速的动态发展机制,对该压气机开展全周非定常数值研究。数值结果表明,在压气机进入失速的过程中,静子叶根区域首先形成堵塞区,验证了失速初始扰动起源于叶根区域。该扰动由六个轴对称分布、以45%转子转速旋转的堵塞团组成。随着失速程度的加深,静子叶根堵塞团在低叶高范围内合并为整圈的分离区,造成整个叶根区域流动的堵塞。叶根的堵塞导致气流从堵塞区上方绕行,加快了叶尖区域的流动速度,激波的增强导致叶顶泄漏涡破碎程度加剧,形成转子叶尖堵塞区并发展为旋转失速团,最终导致压气机失速。     

头部几何微扰对流动不对称性影响的数值研究

使用低耗散的Roe格式,数值模拟了头部几何微扰对细长体不对称绕流的影响.计算发现,周向网格的加密本身也是一种微扰动,会对流场产生很大的影响.如果再加上头部几何微扰,随其高度的变化,侧向力分布显示出类似于几何扰动周向角变化时产生的双稳态现象.在对应的攻角下,头部扰动的方式决定了头部分离的方向,进而决定了后体的流场发展情况.     

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

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

高压压气机叶栅的高速模型和低速模型的相似变换准则

本文以相似理论为基础,从线性小扰动方程、跨音小扰动方程和全扰动势方程出发,推导出一组在进口马赫数不同时,保持两组叶栅气动力相似的变换准则。变换涉及跨音与亚音,跨音和不可压,亚音与不可压及其全扰动的情形。应用二维和准三维无粘数值计算方法进行数值验证,证明了变换关系在无粘下的正确性。文中使用了局部线化概念和一些假设,其合理性一并由数值试验得到证实。本文在理论和数值试验上为压气机叶栅的高速模型和低速模型之间的相似变换提供了必要的基础。      

计算三维粘性／无粘性干扰流动的准联立分区算法

 应用摄动理论导出了粘性／无粘性流动的匹配条件，与边界层方程联立形成能计及粘性／无粘性干扰效应的模型方程。同时引用等效吮吸理论对该模型方程与全位势方程进行交替迭代求解，建立了计算高雷诺数流动的一种快速收敛的粘性／无粘性分区算法。计算结果表明该算法适用于包括小尺度分离的流动；与传统的半逆解法相比，计算效率明显提高     

跨音速透平叶栅速度面设计方法的研究

 低能量损失系数跨音速透平叶栅的设计和实现是提高跨音速透平性能的关键问题之一。本文给出了基于有限面积法的跨音速透平叶栅速度面反问题设计方法。首先推导了基本方程的对称形式,用于上游奇点处理,并对速度面上不连续边界条件作了分析和数值验证说明。然后详细叙述了流函数场的有限面积解法包括主方程的积分变换,基本有限面积格式和边角有限面积格式处理以及解精确度的粗略分析。所编制的程序简明快速。计算结果表明解相当稳定并且不需要在音速线附近等处加密网格或作其他处理。算例也指明音速线上临界奇点位置的选取,对流场及叶型型线形状有一定影响,应在设计时注意。     

旋翼三维非定常跨音速绕流的数值模拟

 利用小扰动位势方程，用有限差分方法求解了旋翼跨音速流场。导出了严格的三维非定常位势方程远场无反射边界条件，并成功地用于数值计算中，缩小了物理求解域的空间尺度。采用了贴体坐标和网格加密技术，以得到合理的非均匀网格分布。应用较少的网格点数成功地计算了无升力情况下旋翼的三维非定常亚、跨音速绕流流场，得到了与实验结果符合良好的数值解     

全速位积分方程跨音速绕流的数值解

从全速位方程出发,利用Green公式将其化为激波捕获积分方程和激波装配积分方程,然后离散进行数值解。流场出现激波时,对激波捕获积分方程应用上风技术捕捉到激波,然后应用激波装配技术计算,得到了满意的结果。经算例考核,该方法具有计算区域小,收敛快和CPU时间较少等优点。     

跨音速三维机翼及翼身组合体非守恒型全位势流AF3算法

 本文对跨音速非守恒型全位势流,用Baker的快速有限差分法,改变其网络拓扑,研制出分析机翼及翼身组合体跨音速绕流的实用程序。其中包括:对“C-H”型贴体坐标网格提出AF3算法的具体实施方案取得了快速收敛效果;采用张量形任意曲线坐标上的非守恒全位势方程,使得程序中跨音速求解部分既能用于分析单独机翼又能用于翼身组合体跨音速绕流计算。     

一种求解核函数方程的新方法—辅助核法

本文提出了求解核函数方程的辅助核法。此法克服了偶极子栅格法对网格分布的限制,具有更大的适用范围。数值计算表明,此法对求解绕翼型的跨声速激波流动是很有利的。     

翼-身组合体跨声速流场计算

本文把Holst的AF2因式分解格式推广到求解机翼-机身跨声速绕流的全速势方程。利用求解椭圆型偏微分方程组来数值生成贴体的计算网格。算例表明本方法运算可靠,费时较少,可以向工程实用部门推广使用。     

跨音速自适应风洞的任意初始壁方案

提出一种新的二维柔壁自适应风洞进行跨音速三维模型实验时,基于任意初始壁状态的洞壁调整方案,讨论了在较高跨音速实验Ma数的实验条件下,采用洞壁预调效应处理的洞壁自适应方案的实验结果。