冲压发动机助推段压强振荡现象数值分析

为了揭示串连式冲压发动机飞行助推段进气道和补然室内压强振荡产生的机理,针对冲压发动机的助推段开展了进气道和补燃室内流动的大涡模拟,获得了进气道和补然室内压强振荡的频率和振幅,并和飞行实验数据进行了对比。结果表明,在助推段串连式冲压发动机进气道和补燃室内发生的压强振荡现象属于典型的整体模式不稳定现象;助推飞行马赫数、飞行高度的变化对整体模式不稳定现象压强振荡的频率和振幅影响不大。     

轴流风扇自噪声强度特性的数值研究

利用流场/声场混合模型计算分析轴流风扇自噪声的强度特性。使用大涡模拟(LES)数值计算方法获取湍流流场信息,以此为输入,使用声类比方法噪声预测模型获得了不同来流条件下的自噪声频谱。0°攻角条件下,叶片噪声源主要位于叶片尾缘区域,并且叶背上的噪声强度大于叶盆;攻角增大和速度增大都会引起叶片自噪声强度的增大;攻角增大引起的噪声增大主要集中在低频区域。     

涡破裂诱导的垂尾抖振数值模拟

 采用数值模拟方法研究了大后掠三角翼前缘涡破裂诱导的垂尾抖振问题,分析了大迎角条件下的垂尾抖振特性。采用Navier-Stokes方程求解非定常气动力、耦合结构动力学方程,建立了气动弹性方程,在时域内采用松耦合方式推进以得到垂尾结构响应。研究结果表明:涡破裂流的脉动频带覆盖了垂尾扭转模态的固有频率,诱发了垂尾抖振现象;与传统的颤振频域响应特性不同,垂尾抖振响应的各阶位移与加速度响应主频均位于各阶结构模态固有频率附近。此外,弯曲与扭转响应存在耦合效应,且耦合作用的频率与提取的垂尾表面气动载荷脉动频率一致。垂尾的位移响应由一阶弯曲模态主导,振幅不大;加速度响应主要由扭转模态产生,量级较大,使结构持续遭受严重的附加惯性载荷作用。     

高超声速进气道强制转捩流动的大涡模拟

吸气式高超声速飞行器常在进气道边界层内布置粗糙颗粒或涡流发生器强制流动转捩为湍流以确保发动机正常启动。为了清晰认识强制转捩过程,采用隐式大涡模拟方法,对强制转捩问题开展了数值模拟研究。针对钻石形和斜坡形涡流发生器,计算得到涡流发生器诱导的流动结构,显示出强制转捩流动由涡流发生器产生的反向旋转流向涡对主导。扰动沿流向增长和发展,导致流向涡对以偶模式或奇模式失稳,偶模式失稳产生对称形式的涡对破碎,而奇模式失稳则导致非对称（弯曲）形式的涡对破碎。流向涡对破碎后产生一系列发卡涡并最终促使边界层转捩为湍流。最后就计算网格和数值耗散对隐式大涡模拟结果的影响以及计算的收敛性进行了讨论。     

侧风下的汽车风振噪声研究与控制

采用大涡模拟的计算方法,对某轿车在侧风工况下的风振噪声特性进行了研究。首先,通过实车道路试验验证仿真方法的准确性;其次,采用上述计算方法分析不同侧风速度、角度对风振噪声的影响;最后,提出在 B 柱内壁上使用 V 型沟槽抑制风振噪声的方案。结果表明,数值仿真与试验结果比较吻合;在侧风速大于5 m/s、侧风角度小于140°时的风振噪声比没有侧风高,根据侧风角度不同最高可达26 dB,不适合开窗;所采用的控制方案降噪幅度可达7 dB,效果明显。     

基于混合方法的膨胀腔非稳态流动及气动噪声分析

采用混合方法对气流作用下的膨胀腔的非稳态流动和气动噪声进了数值仿真分析。建立基于有限体积法的流体仿真模型,采用大涡模拟方法对膨胀腔内部非稳态流动进行计算。数值仿真结果捕捉到了腔体内部非稳态流动以及涡量的发展变化过程,并与实验观测结果基本吻合。采用积分插值的方法将流场结果插入声学网格上,基于气动声学理论计算腔体内部的流噪声源。建立声学仿真模型,将流噪声源导入声场网格中计算远场响应点的声压,并与实验结果进行对比,声压级以及共振频率都吻合良好。结果表明:混合仿真方法准确直观揭示出管内流噪声产生的机理和过程,尾管辐射噪声的仿真结果与实验结果及理论计算结果的偏差都在5%以内。      

超声速喷流混合流场大涡模拟

以光学窗口外冷喷流为研究背景,采用大涡模拟方法对后台阶外形切向喷流混合流场进行了研究。数值方法基于隐式亚格子模型,采用高精度WENO格式进行空间离散,并通过超声速平面混合层流动对数值方法进行了考核验证。喷流混合流场计算模型与试验一致,来流和喷流马赫数分别为3.4和2.5。数值模拟清晰地捕捉到了流场波系以及混合剪切层、壁面边界层等典型流场结构,并精细预测了混合层发生失稳、转捩及发展为充分发展湍流的时空发展过程。数值模拟得到的湍流大尺度结构的位置和形态与实验图像一致。通过对瞬时流场、统计平均流场和脉动参数的分析,揭示了流场结构特征及其时空演化规律,并获得了流场密度脉动特性。      

不同喷压比下欠膨胀射流对特征频率激励的响应

为了揭示欠膨胀激励射流的流动机理,以及考察不同喷压比下射流对相同激励的流动响应,采用大涡模拟方法,对喷压比NPR=5.60和9.34的欠膨胀定常射流和激励射流进行了三维数值计算。激励频率为定常射流中固有的轴对称频率=14.569kHz,激励形式为在射流喷管入口处施以正弦压力扰动。结果表明,特征频率激励影响射流的声场特征,缩小射流核心区的范围,减少射流近场的激波胞格数目,并影响射流气体与环境空气的混合。同时,激励射流的特征频率转变为激励频率及其高阶倍频,激励射流的主不稳定模态均为轴对称模态。其中,NPR=9.34的欠膨胀射流的主不稳定模态和外加压力扰动的形式相一致,射流与外加激励发生了更加剧烈的流动耦合和响应。这使得在NPR=9.34时,射流核心区长度减小得更多,压力脉动的振幅更大,激励对射流混合的增强作用更加明显。     

Numerical study of transonic separation flow using an anisotropic turbulence model

A transonic turbulent separation flow in a converging-diverging transonic diffuser was studied, when there existed a separation bubble on the top wall of the diffuser triggered by strong shock-wave-boundary-layer-interaction (SWBLI). To capture the essential behavior of this complex flow, the current study utilized an anisotropic turbulence model developed on the basis of a statistical partial average scheme. The first order moment of turbulent fluctuations, retained by a novel average scheme, and the turbulent length scale, can be determined from the momentum equations and mechanical energy equation of the fluctuation flow, respectively. The two physical quantities were readily used to construct the nonlinear anisotropic eddy viscosity tensor and to significantly improve the computational results. Comparisons between the computational results and experimental data were carried out for velocity profiles, pressure distribution, skin friction coefficient, Reynolds stress as well as streamline vectors distribution. Without using any empirical coefficients and wall functions, the numerical results were in good agreement with the available experimental data, further confirming that the nonlinear anisotropic eddy viscosity tensor is the decisive factor for the success of the computational results.     

低旋流预混燃烧稳燃机理的大涡模拟

采用大涡模拟(LES)方法对比了高低旋流的预混燃烧特性,并分析了低旋流的稳燃机理。结果表明:LES耦合部分搅拌反应燃烧模型(PaSR)可以准确的捕捉旋流预混燃烧的细小火焰结构和流场分布。不同旋流强度引发不同的流场特性,低旋流诱发较弱的回流区,而高旋流会诱发很强的回流区,这种强烈的流场特征往往会引起振动频率极低的进动涡核拟序结构(PVC)。通过涡脱落和燃烧温度分布的瞬时耦合关系,发现火焰剪切层引起的内外涡脱落频率差是低旋流火焰稳定的最主要因素,低旋流燃烧受涡旋结构的影响交替出现W和V型火焰。