基于混合笛卡尔网格方法的扑动翼型数值模拟

针对单体扑动翼型与前扑动/后静止串联布置的双翼型发展了适用于运动边界问题的非定常混合笛卡尔网格方法,并进行了数值模拟研究.在物面附近使用贴体结构网格,外部使用自适应笛卡尔网格来填充,使用最近“贡献单元”方法来实现两套网格之间的信息传递.发展了一套非定常可压缩求解器,使用格心格式2阶精度的有限体积方法实现空间离散,使用隐式LU-SGS双时间步方法实现时间离散.物面边界运动过程中,在贴体网格外边界进行笛卡尔网格的动态几何自适应,使用逆距离插值方法进行新鲜网格单元参数的确定.对扑动翼型的研究重点关注了推力系数与推进效率:除却很小的扑动幅值与减缩频率的工况,在单体扑动翼型后部一倍弦长处放置一个静止翼型能够增大推力系数;但推进效率的改变较为复杂,与扑动的幅值以及减缩频率相关.     

直升机主减机匣结构振动噪声分析与优化

针对直升机主减速器机匣的振动噪声问题,对机匣进行基于频率响应和模态贡献量的结构动力学特性分析,给出机匣的振动特性,并确定对结构振动特性起主要影响的模态.利用间接边界元与有限元相结合的方法,应用基于结构面板声功率贡献量的分析方法,进行结构噪声功率分析和结构噪声功率面板贡献量分析,给出机匣的声学辐射特性,找出相应激励频率下对结构噪声功率贡献量最大的面板.以该面板为设计域进行结构拓扑优化,并根据优化结果合理布置加强筋,以提高结构刚度,达减振降噪的目的.结果表明:结构速度频率响应峰值下降了36%,减振效果良好,结构声功率级有了明显的降低,其中声功率级峰值下降了5dB,降噪效果良好,为直升机主减机匣提供了一种可行的减振降噪方法.      

Sajben扩压器复杂流动中湍流模型的性能评估

用8个常用的湍流模型对Sajben扩压器中跨声速流动进行了数值模拟,评估了Spalart-Allmaras, 标准k-ε, RNG (re-normalization group) k-ε,realizable k-ε,标准k-ω,SST(shear stress transport) k-ω,v2-f,Reynolds stress共8个湍流模型对激波/湍流边界层相互作用的模拟预测能力.通过与实验数据比较发现:SST k-ω模型和v2-f模型比其他模型模拟的更准确,其中SST k-ω模型比v2-f更能准确地预测壁面压力,然而对于分离点、再附点以及分离区长度v2-f比SST k-ω预测得更准确.      

后掠机翼人工转捩最佳粗糙带高度数值预测

介绍了基于当地变量的γ-Re_θ转捩模型,并将该模型应用到后掠机翼的转捩预测和人工转捩最佳粗糙带高度以及人工转捩技术能够模拟的大气飞行雷诺数的确定中。为检验γ-Re_θ转捩模型对后掠机翼转捩的预测能力,对ONERA M6机翼和DLR-F4标模机翼进行了边界层转捩预测,采用结构化网格和有限体积法求解雷诺平均Navier-Stokes(RANS)方程,得到了机翼表面的摩擦阻力系数分布,从而可以得到相应的转捩位置,预测得到的转捩位置与试验结果比较吻合,说明该模型对后掠机翼转捩预测是可信的。最后在DLR-F4标模机翼上表面固定了粗糙带,通过相同的方法得到了转捩位置,从而确定了马赫数为0.785、雷诺数为3.0×10⁶时最佳粗糙带高度为0.11 mm;通过不断增大雷诺数使自由转捩位置不断向前缘移动,验证了人工转捩对大气飞行雷诺数的模拟能力。结果表明,在最佳粗糙带高度为0.11 mm下,可以实现对大气飞行高雷诺数的模拟。     

非定常预处理方法在倾转旋翼飞行器悬停状态气动干扰模拟中的应用

针对倾转旋翼飞行器的多种速度尺度并存的复杂流动问题,发展了能够统一求解同时包含不可压缩流动和可压缩流动的非定常预处理方法和一种考虑了流场物理信息传播方向的预处理远场边界条件。结合适用于倾转旋翼飞行器的复杂外形及大位移问题的非结构混合网格技术和重叠插值技术,建立了倾转旋翼飞行器的数值模拟方法。在此基础上,首先通过悬停算例验证了数值方法的有效性,确定了选择预处理截断参数的原则。然后探究了倾转旋翼飞行器直升机模式悬停状态的尾迹涡结构演化过程,研究了“喷泉效应”和“地面效应”对倾转旋翼气动性能的影响。计算结果表明,旋翼尾迹涡的耗散在不同方位角存在较大差异。“地面效应”抑制了旋翼尾迹的发展,提高了旋翼整个方位角区域的升力系数,而不仅限于机翼上方方位角区域。 “地面效应”作用强于“喷泉效应”作用,使旋翼相对于孤立旋翼拉力系数略微增大了5.26%,但同时整个倾转旋翼飞行器的总的向上的力只有孤立旋翼拉力的82.46%。结论表明预处理方法和预处理远场边界条件能够较好地模拟倾转旋翼飞行器悬停状态的干扰流场,兼顾了数值计算的耗散性和稳定性。气动干扰对倾转旋翼飞行器整体气动性能有重大影响。      

不同波系配置的鼓包压缩面引起的激波/边界层干扰特性

为了对不同波系配置的鼓包压缩面引起的锥形激波和来流边界层干扰的流场特性进行对比分析,分别以总偏转角相同的单锥、双锥和等熵锥等3种外锥流场为基准流场,在来流马赫数为2.0的条件下,采用流线追踪法设计了3种鼓包压缩面。在此基础上,采用数值仿真的方法对这3种鼓包压缩面在不同马赫数下的三维激波/边界层干扰流场进行了对比分析,结果表明:相比于单锥鼓包,双锥和等熵鼓包流场中的涡量集中区明显减小,且涡核处的涡量值也有所降低;对于双锥和等熵鼓包,由于受到后续波系的干扰,其表面流动的拓扑结构不再呈现出锥形相似性;双锥与等熵鼓包具有和单锥鼓包相同的边界层排移能力,且这一能力不随来流马赫数的变化而变化。此外,采用多波系配置能提高鼓包外压缩系统的总压恢复系数,波系数目越多,提升的程度就越大,且这一优势在高马赫数下更加明显。      

颤振边界预测的系统稳定性分析方法

对颤振边界预测方法中的颤振裕度法与自回归滑动平均模型（ARMA）稳定性分析方法进行研究,采用数值仿真算例研究这两种方法受参数识别误差的影响,采用风洞颤振试验算例对比分析其预测颤振边界的有效性与准确性,结果表明:①颤振裕度法的判据受阻尼比识别误差的影响较小,相对于阻尼比识别误差更容易受频率识别误差的影响,ARMA稳定性分析方法的判据相对于颤振裕度法更易受阻尼比识别误差的影响,颤振裕度法比ARMA稳定性分析方法鲁棒性更高;②对于机翼弯扭耦合颤振,两方法的判据变化趋势相近,相差采样频率4次方的数量级,两者随着风速的增加具有平缓的下降趋势,有助于较早地预测颤振边界;③对于面内弯曲为主的颤振,ARMA稳定性分析方法比颤振裕度法适用性更好。      

受限混合层的线性稳定性

针对超声速边界层/混合层组合流动,利用可压缩线性稳定性理论研究了流动的线性失稳特性。基本流场选取了具有不同速度特征的两股来流,采用双曲正切的混合层剖面叠加可压缩边界层自相似性解剖面构造。重点考察了混合层中心与壁面距离、对流马赫数等参数对组合流动稳定性特征的影响,其中壁面采用绝热壁面。混合层中心与壁面的距离为5~15倍的边界层厚度,混合层的对流马赫数为0.6~1.2。结果表明:该组合流动中存在独特的多重不稳定模态,并相互影响;且其不稳定模态随着壁面距离及对流马赫数的变化呈现出不同的主导行为。      

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.     

跨声速风洞槽壁干扰评估与修正技术的应用

目前先进跨声速风洞试验段多采用槽壁形式,而国内对槽壁干扰的认识较少,尚未进行过槽壁干扰的评估和修正工作。本文基于理想槽壁均匀边界条件和经典方法,利用先进构型民机标模对槽壁干扰特性进行了评估和修正,在验证方法准度的基础上,对比分析了国内2.4 m跨声速风洞和欧洲跨声速风洞（ETW）槽壁干扰的差异和规律。结果表明,修正后的试验数据与ETW参考数据吻合较好,修正量、干扰因子与ETW评估结果一致。与传统孔壁试验修正方法相比,本文方法原理清晰,计算简便、快速,可方便应用于其他类似槽壁风洞中,可作为壁压信息法的辅助手段以提高国内试验数据质量,并为国内先进跨声速风洞的设计和调试提供借鉴。