基于全局线性稳定性分析的翼尖双涡不稳定特征演化机理

相比于机翼产生的孤立翼尖涡，加装小翼之后的翼尖涡表现出双涡甚至多涡结构，并且呈现出更加复杂的不稳定特征。为揭示翼尖双涡结构不稳定特征及其演化机理，采用体视粒子图像测速（SPIV）技术和全局线性稳定性分析（LAS）方法对不同雷诺数和攻角下带双叉弯刀小翼的M6机翼产生的翼尖涡结构在尾迹区的不稳定特征进行研究。试验结果表明，对称布置的双叉弯刀小翼产生的翼尖涡包含上/下小翼产生的主涡（上/下主涡）结构，两者构成近似等强度的同转涡对，在相互靠近的同时以20 rad/s的角速度相互缠绕。对上/下主涡瞬时涡核位置的统计分析表明，翼尖涡摇摆幅值随流向位置逐渐增大，随雷诺数的增加而增大，随攻角的增加先增大后减小。对16倍弦长的尾迹截面处的翼尖双涡结构进行全局时间稳定性分析，不同工况下，上/下主涡最不稳定模态（模态P/模态S）的稳定性曲线变化规律与摇摆幅值的变化规律相一致，表明翼尖涡的摇摆源自于其内在的不稳定性特征。增加流向扰动波数，发现模态P切向波数逐渐增加；而模态S则是径向波数逐渐增加。不同工况下，模态P的切向波数为5~6，扰动波数分布在[2.75，5]的区间内，所对应的不稳定放大率均大于模态S，而不稳定放大率最大的模态扰动范围作用在上主涡的整个涡核区域，表明这种大切向波数的扰动模态在翼尖涡流控中的潜在价值，也意味着加装小翼会增加涡结构的个数，增强不稳定性的发展，有助于翼尖涡的快速失稳衰减。     

采用DMD方法研究叶栅不同攻角的拟序结构

为分析平面叶栅分离流非定常拟序流动特征,对三个不同攻角下的叶栅进行了单通道的大涡模拟仿真,并采用动力模态分解（DMD）三个工况的流场结构进行了分析。DMD方法对包含复杂时空信息的叶栅分离流流场进行了解耦,剥离出了反映流场主要动力信息的模态,获得了其频率和与之对应的空间结构。并且通过DMD方法,将原本需要研究大量不同时刻的流场,转移到仅需要对少量模态的研究即可,实现了保留主要动力特征的低维近似。通过DMD分析表明:气流经过叶片前缘产生流动分离,形成不稳定的剪切涡结构,它和尾迹区脱落涡相互耦合,并形成新的拟序结构。随着攻角的增大,前缘剪切涡及其与尾迹涡的耦合也同时增强,流场变得更加复杂。      

基于DMD方法的翼型大迎角失速流动稳定性研究

基于剪切应力输运湍流模型的SST-DDES混合方法对NACA0012翼型大迎角分离流动进行非定常数值模拟,采用动力学模态分解(Dynamic Mode Decomposition,DMD)数学工具对失速初始状态、浅失速状态以及深失速状态的流场进行稳定性分析。结果表明:DMD方法准确地提取了翼型大迎角流动中的主频和高阶倍频及对应的流场模态结构;与FFT分析结果相比,频率最大差异小于0.16%;且发现两者提取的频率在流动中的主导作用顺序也一致。通过特征值对相应的模态进行稳定性分析,所有模态的放大率均非常小,所有模态处于弱发散、弱收敛或稳定极限环状态。DMD提取的一阶模态主要表现为分离涡演化过程中最主要的静止分离涡结构,前三阶低频对应的模态涡结构与流动中以此频率进行演化的涡结构比较一致,更高阶的倍频主要表现为尾涡和尾迹区的涡结构。且发现不同模态系数之间存在相位差,说明分离涡流动中不同频率对应的涡结构运动不同步。     

展向不同间距壁面粗糙元对近壁湍流拟序结构的影响

采用氢气泡流动显示技术,研究了壁面不同展向间距粗糙元对近壁湍流拟序结构的影响。实验中粗糙元布置在氢泡丝前,其间距按如下两种排列设置：粗糙元直径的0、1和3倍,等间距1、2和4cm;每种情况粗糙元直径均为2.3、4和6mm。基于平均流速和水力直径的流动雷诺数从14000到48000变化。实验观测了光滑壁面和粗糙元壁面的近壁湍流拟序结构条带特征。发现：对于两种不同排列,相同雷诺数和相同粗糙元间距条件下条带有量纲间距均随着粗糙元直径增大而减小,条带有量纲高度均随粗糙元直径的增大而增大;相同雷诺数和相同直径粗糙元条件下条带有量纲间距均随着展向间距增大而增大,条带有量纲高度均随展向间距的增大而减小。粗糙元壁面湍流条带有量纲间距均比光滑壁面小,条带有量纲高度均比光滑壁面大。研究结果对近壁湍流的流动控制具有重要指导意义。     

吹吸扰动控制猝发过程中雷诺应力的实验

采用热线测速技术和双丝边界层热线探针,在风洞中以高于对应最小湍流时间尺度的分辨率,精细测量了平板湍流边界层施加不同频率周期性质量引射前后,不同流向和法向位置的瞬时流向、法向和展向速度分量的时间序列信号.采用流向脉动速度局部平均结构函数模极大值法检测了壁湍流猝发过程中拟序结构不同速度分量和雷诺应力分量的条件相位平均波形.发现在扰动源下游一定范围内,摩擦因数和雷诺应力幅值有所减小,扰动改变了法向或展向脉动速度与流向脉动速度的条件相位平均波形的相位差,使其保持在π/2的奇数倍左右,此时雷诺应力几乎为0,改变了湍流边界层原有雷诺应力的产生和扩散输运机制,抑制了湍流雷诺应力的产生和输运.      

淹没射流湍流场的 TR-PIV 测量及流场结构演变的 POB 分析

利用平面激光诱导荧光流动显示技术（LIF）和二维高速粒子图像测速技术（TR-PIV）对淹没射流和自由面相互作用湍流场进行了细致的测量,并利用本征正交分解方法（POD）对测得的流场进行分解,提取流场中含能大尺度结构并分析其与自由面相互作用特点。射流出口雷诺数 Re ＝5600,淹没深度 H/D ＝4。LIF 实验结果表明：自由液面的存在使得射流在向下游扩散发展中存在向上运动的趋势。时均流场也表明射流上半部分扩散速率更快,在相互作用区域,最大速度位置逐渐向自由液面靠近,类似自由射流的流向湍流度双峰值分布形式随流向距离增大而逐渐消失。POD 分解得到的空间模态表明：射流中有序的相干结构在上游逐渐发展,随后迅速向上发展并与自由液面发生相互作用,在下游,由于自由液面的限制,大尺度结构又开始向下发展。     

来流边界层拟序结构对超声速混合层流场结构的影响研究

超声速混合层涉及可压缩湍流的根本问题,具有重要的应用背景。通过设计超声速混合层实验装置、应用新近提出的高分辨率NPLS测试技术,拍摄了来流边界层分别为层流和湍流流态下混合层的流向和展向流动图像。根据流动图像的特征,分析了混合层的流向与展向流场切面中拟序结构的成因;深入讨论了来流边界层中拟序涡结构与混合层涡结构的相互作用问题;比较了层流和湍流来流条件下混合层拟序结构的异同及其对混合效率的影响。结果表明：当来流边界层为湍流时,对应的混合层具有较高的混合效率。     

压气机动静叶干扰的非定常流动特性

采用大涡模拟(large eddy simulation,LES)方法计算动叶尾迹对静叶干扰的流场信息。利用涡量分布揭示动叶尾迹在静叶通道内的演化过程,利用压力梯度识别激波结构及波振源,运用动力学模态分解(dynamic mode decomposition,DMD)方法对静叶通道内流场的时空结构进行模态分解。结果表明:流场中存在3处波振源,分别位于动叶尾缘、静叶前缘和静叶尾缘处;发现静叶通道内流场的频谱具有多峰值现象,模态分解的第1阶流动代表动叶尾迹在通道内随时间迁移,对应频率为动叶通过频率(blade passing frequency,BPF)是通道内旋涡非定常波动的主导频率;第2阶流动是动叶通过频率的2倍频流动,旋涡的空间尺度为1阶模态的1/2,为更小尺度的扰动。      

超声速/高超声速边界层转捩后期流场的模态分析

将动力模态分解(DMD)方法应用到超声速和高超声速边界层转捩后期的流场分析中,通过获得流场主要的相干结构和对模态的重构,研究了相干结构与壁面阻力和热流的关系。结果表明超声速和高超声速边界层转捩的流场结构存在明显差别。超声速转捩流场由低频流向涡的模态主导,这些模态对转捩后期的壁面阻力和热流有重要贡献;高超声速转捩流场中存在多个不同量级频率的模态,在DMD频谱上表现为多个不同的分支,通过对不同分支能量最高的模态进行考察,我们发现低频模态的结构为流向条带,高频模态的结构为二维扰动波,这些模态对壁面阻力和热流的影响与模态的结构形式类似。     

混合层中由次谐波导致的失稳过程的数值研究

用拟谱方法对随时间发展的二维混合层流场进行了直接数值模拟,给出了高雷诺数下在Ｋｅｌｖｉｎ－Ｈｅｌｍｈｏｌｔｚ波卷起后存在的限振幅次振波地大尺度拟序结构进一步的失稳演变过程,包括大涡的配对,撕裂以及二次配对,同时探讨了初始的次谐波强度,雷诺数相位差等对大涡配对和合并过程的影响。     

Detached-eddy simulation of supersonic flow past a spike-tipped blunt nose

Space vehicle in atmosphere travels mostly at supersonic speed and generates a very strong bow shockwave around its blunt nose. Oblique shock and conical separated flow zone generated by a forward disk-tip spike significantly reduce the drag by reducing the high pressure area on the blunt nose. This study employs improved delayed detached eddy simulation to investigate the characteristic flow structures around a spike-tipped blunt nose at Mach number of 3 and Reynolds number(based on the blunt-body diameter) of 2.72x10~8. The calculated time-averaged quantities agree well with experimental data. Characteristic frequencies in different flow regions are extracted using fast Fourier transform. It is found that two distinct instability modes exist: oscillation mode and pulsation mode. The former is related to the foreshock/turbulence interaction with nondimensional frequency at around 0.004. The latter corresponds to the interaction between turbulence and shock structures around the blunt nose, with a typical coherent structure shedding frequency at 0.092.     

带斜孔肋大宽高比矩形通道的强化传热特性

为了获得带斜孔肋大宽高比矩形通道的强化传热特性,并寻求最佳的孔排倾斜角度,调节孔排倾斜角度和通道雷诺数,使其分别在0°~30°和3×104~9×104范围内变化,通过数值计算系统分析了通道摩擦因数和带肋壁努塞尔数等参数的变化规律.研究发现:相比于常规肋片,新型斜孔肋有效改善了肋片后方紧邻肋片的局部区域的壁面换热,并降低通道的摩擦因数,但传热增强因子有所减小;随着孔排倾斜角度的增大,通道的相对摩擦因子单调升高,传热增强因子则呈现出先升高后降低的变化过程,因此存在着最优孔排倾斜角度为15°,此时斜孔肋的强化传热综合指标达到最大值;随着通道雷诺数的增大,斜孔肋通道的摩擦因数小幅减小,换热则逐渐增强.      

出口侧向扩张孔流量系数试验

采用试验研究的方法,探讨了不同结构参数出口侧向扩张孔的流量系数在不同工况下的变化规律。通过试验研究发现:出口侧向扩张孔圆柱段孔径为0.6 mm时的流量系数比孔径为0.8 mm时大3%;当孔倾角由30°减小至25°时,流量系数随之增大,当倾角继续减小至20°时,流量系数不变;当扩张角由20°增大至40°时,流量系数先增大后减小,且在扩张角为25°时最大;流量系数与次流通道雷诺数成正比,雷诺数越大,流量系数达到极限值所需的压降系数越小,且当压降系数大于0.9后,不同次流雷诺数下的流量系数差别在2%左右。      

HyTRV流场特征与边界层稳定性特征分析

高超声速转捩研究飞行器（HyTRV）是为研究三维复杂外形的边界层转捩问题而设计的一款具备真实飞行器典型特征的升力体标模。为支撑更加全面系统的理论分析、数值模拟、风洞试验和飞行试验研究,采用高精度数值模拟方法、线性稳定性理论（LST）和e^N方法对HyTRV标模的典型流动特征和边界层失稳特征进行了分析。研究表明,HyTRV展现出多个相对独立的横流区域和多个流向涡结构;HyTRV的边界层存在横流失稳模态、第二模态、附着线失稳模态等常见模态。横流失稳模态出现在周向高低压区之间的横流区域,能够主导转捩发生;横流区域同时也存在第二模态,其N值普遍比横流失稳模态小;附着线失稳模态呈现出第二模态特性,且频率非常高。还研究了攻角和单位雷诺数的影响。结果表明,随着攻角增加,标模下表面中心线的流向涡结构逐渐消失,横流雷诺数逐渐减小;上表面流向涡结构逐渐从腰部移向顶端,并出现新的流向涡结构。增加攻角,所有失稳模态的N值总体上逐渐减小;增加单位雷诺数,N值显著增加。基于研究结果,针对流向涡失稳、横流失稳、第二模态和附着线失稳等给出了研究建议。     

斜拉索的阻力系数研究

对几种斜拉索的阻力系数进行了研究,得到阻力系数对于倾角及雷诺数变化的曲面拟合函数,曲面图直观地显示了阻力系数变化.表面花纹索和表面螺旋线索的阻力系数变化较为相似,都比光洁索的阻力系数变化要平缓,但表面花纹索的阻力系数要略小于螺旋线索.     

引入DMD方法研究有/无控气流分离的动态结构

为分析非定常流动控制技术抑制分离流的机理,对弯曲扩压通道的试验模型进行了数值模拟,针对扩压通道在无控和采用最佳射流频率状态下的计算结果引入了动力模态分解(DMD)技术进行分析。通过DMD技术能够将包含时空信息的扩压通道复杂流场进行分解,捕获流场包含的动力信息和对应的拟序流动结构。将无控和有控流场分解的结果进行对比分析后表明:采用有效激励措施时,和脱落涡频率一致的涡系对流场的影响更加突显,流场整体上表现得更加有序;非定常控制抑制了一部分涡的增长,使得各模态整体上更加稳定;而有控流场占主导地位的涡系结构相比无控流场较为有序,且对主流区未形成明显的直接影响。     

平行四边形截面的旁路式双喉道气动矢量喷管数值研究

将异型出口设计与气动推力矢量喷管结合,提出了平行四边形截面的旁路式双喉道气动矢量喷管(bypass dual throat nozzle,BDTN),并与基准矩形截面的BDTN进行流场结构及气动性能的对比研究。三维数值计算表明:平行四边形构型与矩形构型具有相同的气动性能变化规律;相较于基准矩形BDTN构型,由于壁面倾斜导致流场结构变化,平行四边形构型在矢量状态下的俯仰矢量角及推力系数有所降低,但对非矢量状态下的推力系数和流量系数影响不大;壁面倾角是性能变化的重要因素,相同落压比时壁面倾角越小,矢量角越小,壁面倾角不小于60°时喷管的稳定矢量角均可达10°以上,最大矢量角均可达15°以上;平行四边形出口增强了尾喷流与环境气流的掺混,出口射流中心线速度衰减大大加快,有利于提高红外隐身特性。