两种湍流模式对绕椭球体三维分离流计算的比较

将新近的Shih Lumley(S L)模式应用于雷诺数为7 2×106的6∶1椭球绕流的全N S方程计算,作为参考也计算了Jonse Launder(J L)模式。在没有发生分离时,两种湍流模式计算的平均流场与实验结果符合较好。在分离发生后,J L模式计算的分离结构与实验差别很大,没有计算出二次分离线;S L模式的计算结果较好,计算得到的表面摩擦力线图谱中主要分离结构与实验结果较接近,分离区的表面压力系数分布也与实验值定性相符;在S L模式计算结果中,壁面附近出现小尺度的分离再附结构.     

绕椭球三维流动的分离结构随雷诺数的变化

通过数值计算求解N S方程研究了绕4:1椭球流动结构随雷诺数的变化。计算的雷诺数范围为102-5×103。计算结果系统地显示了分离结构随雷诺数的变化情况及分离涡的空间结构。研究表明,在所计算的雷诺数范围内,10°和30°迎角的分离涡系有本质区别。分离涡系在10°迎角下以闭式分离为主导,在30°迎角下则以开式分离为主导。     

不同雷诺数下翼型气动特性及层流分离现象演化

低雷诺数下空气黏性效应突出,翼型表面普遍存在层流分离现象,相比常规雷诺数情况气动特性显著恶化。采用带预处理的Roe方法求解非定常可压缩Navier-Stokes方程的数值模拟技术和低雷诺数低湍流度风洞油流显示试验技术,对FX63-137翼型不同雷诺数下气动特性和流动结构展开深入研究。通过风洞油流显示试验可以清晰获得低雷诺数层流分离流动的两道油流汇集线。数值模拟结果表明其分别为时均化主分离线和二次分离线,两种结果定性定量均吻合较好,证明了本文的研究方法有效可靠;雷诺数从500 000降至20 000,翼型气动特性和层流分离流动结构均发生显著的变化,伴随阻力系数剧增和升力系数剧降,时均化流动结构从附体至出现经典的长层流分离泡,并最终演化为后缘层流分离泡,相应的两种分离泡的非定常流动结构也存在显著差异;对于阻力系数和升力系数而言,存在不同的临界雷诺数,因为导致阻力系数剧增的机理在于经典长层流分离泡的产生使翼型压差阻力大增,而造成升力系数剧降的主要原因在于后缘层流分离泡使得等效翼型后部弯度减小;非定常结果显示正是由于翼型表面漩涡周期性的生成与脱落,才造成了低雷诺数下升力系数的周期性波动。翼型上表面主分离涡即将脱落时,流线在后缘附近再附,升力系数达到峰值;而当流体从下表面向上卷起二次分离涡时,尾部流线大尺度分离,升力系数降至谷值。     

高超声速下表面凸起干扰气动热实验研究

 对高超声速飞行器表面凸起附近的气流流动和气动加热开展了实验研究和分析。实验在高超声速炮风洞中进行,来流马赫数为8.2、单位雷诺数为9.35×10⁶ m^-1。利用薄膜传热测量方法进行了凸起几何形状和边界层状态对干扰流动加热的影响评估。利用流油图谱和纹影摄像法得到了凸起周围的流动特征:若凸起上游边界层未分离,最大峰值热流发生在凸起侧方附近处;若凸起上游边界层完全分离,最大峰值热流通常发生在凸起的上游表面。实验发现最大峰值热流和来流边界层状态关系不大,原因是流动干扰区表现出较强的三维扰动特性,使得来流层流边界层在干扰区内会转变成过渡甚至完全湍流状态。     

压气机静叶栅层流分离泡转捩与角区分离数值模拟与实验

采用γ-Reθ转捩模型对某可控扩散叶型(CDA)平面叶栅全攻角范围进行了三维数值计算,通过对比数值计算结果与叶栅实验吻合较好。在此基础上,分析了进口来流湍流度和雷诺数变化对叶栅表面层流分离、转捩以及角区分离的影响。结果表明:进口湍流度低于5%时,吸力面存在层流分离,当进口湍流度大于5%后,层流分离移除,但转捩会一直存在;随着进口湍流度或雷诺数增加,吸力面和压力面转捩位置均会前移;随着进口湍流度增加,吸力面角区分离会有所减小,雷诺数增加对角区分离的影响不大。      

小展弦比机翼加装格尼襟翼的低雷诺数试验

通过风洞试验研究了在低雷诺数下加装格尼襟翼的小展弦比机翼气动特性,机翼展弦比为1.67,格尼襟翼为1%~4%弦长高度,试验雷诺数分别为2.0×105和5.0×105.天平测力和表面测压的试验结果表明:低雷诺数下小展弦比机翼加装一定高度的格尼襟翼后,升力系数明显提高,加装1%弦长高度的格尼襟翼还能够提高机翼的升阻比.这是因为在试验雷诺数下,合适高度的襟翼在提高了机翼升力的同时并未显著增大机翼阻力.对比不同试验雷诺数下格尼襟翼的作用效果,表明格尼襟翼能够减少低雷诺数气流分离的不利影响,并且在较小的雷诺数下这种作用更加显著.关于格尼襟翼对低雷诺数层流分离现象的影响,还需要通过细致的流场显示技术进行研究.      

低雷诺数下翼型层流分离泡及吹吸气控制数值研究

在低雷诺数下Eppler387翼型表面会出现层流分离泡现象,因此本文使用Fluent求解器开展吹/吸气控制翼型表面层流分离泡的数值研究,主要探究了射流位置、射流角度、射流速度比三个控制参数对层流分离泡控制的影响规律.研究结果表明,采用与SSTκ-ω湍流模型耦合的γ～(Reθt)转捩模型可以准确预测层流分离泡的位置;吹/吸气可以有效抑制低雷诺数下层流分离泡的发展,明显提高低雷诺数下翼型升阻比.固定射流位置,较大吸气速度比和较小吹气速度比可分别获得较好的流动控制效果,且吸气控制比吹气控制对层流分离泡的抑制作用更加有效.     

安装角对压气机叶栅气动噪声特性的影响

 基于大涡模拟(LES)和边界元方法对轴流压气机叶栅湍流流场以及流场诱导的噪声进行计算,在不同叶栅安装角下研究来流攻角和来流雷诺数对叶栅气动噪声产生、辐射的影响。研究表明:来流雷诺数不变时,同一安装角下,随着来流攻角从-5°~20°变化,叶栅监测曲线上的声压级先减小后增大,在0°来流攻角下声压级达到最小。安装角为45°时,外场总声压级随来流攻角的分布与30°安装角变化趋势相近。但安装角为60°时,总声压级的变化则明显变缓。在0°来流攻角下,总声压级比安装角为30°和45°时增加了近6 dB,但在其他正来流攻角下,变化并不明显。叶栅的最小声压值出现在弦线方向附近,安装角改变时,最小声压级出现的位置也不同。安装角不变,随着来流雷诺数的增大,叶栅表面的分离减小,损失降低。但叶栅表面的压力脉动随着来流雷诺数的增大而增大,使外场辐射噪声增加。     

来流湍流度对超高负荷低压涡轮附面层流动控制手段的影响

实验研究了表面粗糙度耦合上游尾迹的流动控制技术,分析了来流湍流度（FSTI）在流动控制过程中对叶片吸力面附面层分离、转捩特性的影响.实验发现:在速度峰值点至分离点之间布置粗糙高度与弦长之比为1.05×10-4的粗糙条带可以在来流湍流度为0.4%与2.2%的低雷诺数范围内降低叶型损失.在雷诺数为85000的状态下,FSTI影响了尾迹通过区、尾迹诱导转捩区及自然转捩区的附面层动量厚度,造成了叶型损失的差异,但FSTI对抑制区的影响较小.      

低温烧蚀材料的烧蚀实验研究

本文给出了低温烧蚀材料的高超声速风洞烧蚀实验研究的结果。模型外形是半锥角为9°的球锥体,头部半径为15毫米和25毫米,材料为蜂蜡和樟脑,实验条件是M_∞=6,α=0°、2°、6°和8°,Re_∞为(1.90～4.42)×10~7米~(-1)。实验结果表明:随着Re_∞的增加,可以获得层流、转捩和湍流烧蚀表面形态和端头烧蚀外形;证实了在实验过程的后期出现平衡外形,在平衡外形中存在角点;融熔型材料和升华型材料的烧蚀表面形态差别较大;有攻角和零攻角的烧蚀外形有较大差异。     

突然启动流动问题:从不可压到高超声速流动

突然启动流动问题在气动弹性、扑翼飞行和机动飞行中有重要应用价值。鉴于此,对突然启动流动问题涉及的流动与升力演化机制和理论分析方法进行统一介绍并指出目前的理论空白。比较性地介绍了升力随时间变化的原因、预测升力演化的理论方法和主要变化规律。从中发现,不可压缩和可压缩突然启动问题存在各自独立的研究方法与结论,因此进行统一介绍与分析对建立二者之间的关联有指导意义。同时指出,大迎角可压缩突然启动问题尚无理论分析方法和研究结果。作为一项补充研究,采用数值计算发现一个之前尚未报道的现象,即升力系数首先从较低值快速增加至一个峰值,接着快速下降,趋于定常值。该文综述介绍的方法可用于分析突然启动问题升力演化规律。     

涡轮平面叶栅端壁附近的流动和损失

用五孔针对涡轮导向叶栅流场作了详细测量,着重研究了端壁附近的三维流动,分析了叶栅中的二次流现象和流动损失机理。在实验中清楚地探测到了通道涡和马蹄涡的两个分支。实验发现通道涡主宰了整个二次流流场,马蹄涡吸力面分支在与通道涡的相互作用过程中强度逐渐减弱直至消失。在０．６２Ｃｘ截面及其下游截面的吸力面附近,还发现了一个显著的低速高损失核,该核在吸力面上向上爬升并一直延续到叶栅下游。     

不同转速下跨声速轴流压气机内部流动失稳的机理

以跨声速轴流压气机转子NASA Rotor 67为研究对象,采用数值模拟方法,开展100%、80%及60%转速下跨声速轴流压气机内部流动失稳触发机制的机理研究。数值结果与实验数据的对比分析表明:在3个转速下,数值总性能曲线的变化趋势与实验数据符合一致。通过压气机内部流场的详细分析,得出其基本流动机理。在3个转速下,随着压气机节流,叶顶泄漏涡(TLV)的起始位置逐渐向叶片前缘移动,叶顶泄漏涡也逐渐向相邻叶片压力面偏转,相比近峰值效率点,近失速点时在100%、80%以及60%转速下叶顶泄漏涡的偏转角度分别为3°、6°和9°。在100%和80%转速下,叶顶泄漏涡与激波相互作用所导致的堵塞是触发压气机内部流动失稳的机制,并且在80%转速下,叶顶泄漏涡发生破碎;而在60%转速下,泄漏涡在相邻叶片出现的叶顶前缘溢流(LESF)是触发压气机内部流动失稳的主要机制,叶片吸力面尾缘出现的小尺度附面层气流分离(BLFS)不是主要机制。      

雷诺数对轴流压气机转子非定常流动的影响

针对某轴流压气机转子,结合地面试验采用5通道的数值模拟方法在两个海拔高度下对该转子进行数值模拟,分析雷诺数(Re)对压气机中非定常流动的影响。在地面条件下数值模拟绝对坐标系下的频率均落于试验旋转不稳定(RI)的驼峰频率带宽范围内。相对坐标系下相邻通道的数值监测信号互功率谱证明了压气机中周向扰动的存在。该周向扰动的产生是泄漏流与叶顶载荷分布周期性相互作用的结果。高空20 km下压气机中存在特征频率不同的多个周向扰动。此时在叶片近吸力面的尾缘处存在较大范围的流动分离区,由此引起的径向低能流体与叶顶间隙泄漏流相互作用,并共同影响叶顶载荷分布,引起低雷诺数下压气机流动非定常性周向传播。      

Precessing motion in stratified radial swirl flow

Vortex/flame interaction is an important mechanism for unsteady combustion in a swirl combustion system. Technology of low emission stirred swirl(Te LESS), which is characterized with stratified swirl flow, has been developed in Beihang University to reduce NOXemission. However,large-scale flow structure would be induced in strong swirl flow. Experiments and computational fluid dynamics(CFD) simulation were carried out to investigate the unsteady flow feature and its mechanism in Te LESS combustor. Hotwire was firstly applied to testing the unsteady flow feature and a distinct mode with 2244 Hz oscillation frequency occurred at the pilot swirl outlet.The flow mode amplitude decayed convectively. Large eddy simulation(LES) was then applied to predicting this flow mode and know about its mechanism. The deviation of mode prediction compared with hotwire test was 0.8%. The spiral isobaric structure in pilot flow passage indicates that precessing vortex core(PVC) existed. The velocity spectrum and phase lag analysis suggest that the periodic movement at the pilot outlet was dominated by precessing movement. Negative tangential momentum gradient reflects that the swirl flow was unstable. Another phenomenon was found out that the PVC movement was intermittently rotated along the symmetric axis.     

坡面浅层明流流态界定方法之商榷

为了探求坡面薄层水流流态界定方法,从流体力学与泥沙运动学观点出发,对4种不同粗糙尺度床面（光滑玻璃、人工加糙粒径分别为0．075、0．245和0．380mm）、5种不同坡度下坡面薄层水流放水进行试验研究,并对薄层水流阻力规律与明渠水流阻力规律进行了比较分析。在此基础上提出了粘性底层厚度与水深之比作为新的无量纲临界判数,结果表明该判数能更好地反映坡面薄层水流流动型态。     

燃气非平衡流再生冷却流动传热数值模拟

为准确预测液体火箭发动机推力室身部再生冷却换热状况,采用数值模拟方法,对燃气、推力室壁和超临界气氢进行三维流动和换热耦合计算。采用6组分9步反应动力学模型计算燃气的非平衡化学反应,采用DO模型计算燃气辐射换热,考虑超临界气氢物性随温度和压力的变化。获得了室壁温度场、燃气及冷却剂流场。结果表明,Redlich-Kwong方程、Peng-Robinson方程、Lucas法、TRAPP法能分别准确计算超临界氢的密度、定压比热容、粘度、导热系数,采用燃气非平衡流计算所得流场值更符合实际情况。     

低速轴流压气机叶尖非定常流动试验

利用动态压力传感器对一低速轴流压气机转子的叶顶间隙流场进行详细的试验测量,通过对信号特征的分析,对压气机节流过程中叶顶间隙的非定常流动发展演变规律进行了研究。结果表明：压气机完全失速时,叶尖存在一以46.5%转子转速周向传播的失速团;节流过程中,叶尖前缘处的动态压力信号中存在非定常波动的特征频率带,其变化规律与叶顶流场压力非定常波动的能量迁移有关;随着压气机流量减小,叶顶泄漏流影响区域向前缘移动,失速团在叶顶前缘附近产生,并向尾缘方向扩展,最终覆盖叶片全部弦长;近失速工况时,叶顶间隙相邻通道内泄漏流相互作用,造成通道中的低压区“一前一后”交替分布从而形成一个空间上周期约2个叶片通道的扰动波。     

径流及斜流压气机任意曲面叶型长短叶片的造型设计方法

发展了一套径流及斜流式叶轮机的任意曲面叶型长短叶片的造型方法,并在离心式压气机和斜流式压气机的设计中做了具体的应用。新发展了令叶片的厚度按三次样条函数沿流向分布的一种设计思想,得到了比较理想的叶片。同时在径流及斜流式压气机中使用分流叶片,引入了长短叶片的造型设计方法。通过CFD流场检验验证了这种叶片造型方法的可操作性和适用性。这为径流式以及斜流式叶轮机转子、静子的进一步气动优化设计奠定了基础。      

脉动作用下波壁管内非牛顿流体流动特性的研究

在热量和质量传递过程中,流体的流动状态对整个传递过程有着重要的影响,为此,通过实验方法对以聚丙烯酰胺（PAM）溶液为代表的非牛顿流体在波壁管内的流动特性进行压力差和可视化研究。研究发现,PAM溶液存在一个最佳减阻浓度,在该浓度下与牛顿流体相比PAM溶液的转捩点明显提前。此外在脉动流场下,在一个脉动周期内,波壁管内流体存在明显的稳定流动和不稳定流动两种结构;相同净雷诺数条件下,脉动流场下波壁管内流体的流动混合情况比定常流动状态下强烈很多,意味着脉动流场具有更加优越的质量传递特性。