超声速来流边界层厚度对浅腔声学特性的影响

通过分析不同来流边界层厚度与空腔深度比(δ/D)下腔内中心线上的脉动声压级的分布和不同测点的声压频谱特性,讨论了超声速来流边界层厚度对浅腔(长深比分别为12和15)声学特性的影响.试验来流马赫数为1.5,基于每米的雷诺数为2.26×107.结果表明,δ/D减小导致浅腔内的噪声更加强烈,腔前后部的声压级分布更不均匀;除了个别离散频率外,腔内不同测点其余离散频率对应的声压级都有不同程度的增大.δ/D减小引起空腔前部和后部区域几乎整个离散频率范围内的噪声声压级有明显升高;因超声速浅腔流动,腔中部产生的激波的干扰因素的影响,边界层流动特性对浅腔中部区域的声学特性影响较小.      

高亚声速空腔绕流气动噪声特性研究

通过分析空腔底面中心线上声压级分布与不同测点声压频谱特性,着重研究了高亚声速空腔绕流的气动噪声特性。空腔模型长深比分别为6、10和15,自由来流马赫数为0.8,基于每米的雷诺数为1.55×107,测量的空腔前缘的边界层厚度为0.034m。结果表明:空腔后缘处于噪声产生区,声压级较高;闭式和过渡式空腔因深度较小,来流剪切层触及了空腔底面,干扰了从腔后壁向腔前壁的噪声反馈回路,限制了腔内流动自激振荡的形成;开式空腔深度较大,剪切层直接跨过空腔中部、撞击腔后壁,并产生强烈噪声,噪声从腔后壁通过空腔向前壁的反馈回路未受到干扰,故腔内流动出现自激振荡和多个声压峰值频率。     

基于DDES的有／无舱门腔体气动噪声仿真分析

采用延迟分离涡仿真（ DDES）的方法研究了空腔和带有舱门的腔体的气动噪声特性。仿真在自由来流M＝0．85下,长宽深比值为5：1：1的矩形开放腔体展开。空腔仿真的声压级幅值结果与试验数据相比误差在4％以内,频率位置的预测较好。 DDES方法对于气动噪声问题具有很好的仿真准确性。空腔仿真结果和带有舱门的结果对比后表明,舱门迫使剪切层气流更多的撞击后壁进入腔体内部,使腔体内部尤其是中部压强波动增大,声压级在低频部分增量显著,达到了5 dB左右。     

湍流边界层厚度对三维空腔流动的影响

采用脱体涡模拟(DES)方法开展了不同湍流边界层厚度(TTBL)下的三维空腔非定常流动数值计算。空腔长、宽、深比例为5:1:1,来流马赫数为0.85,雷诺数为13.47×10⁶ m^-1,各工况湍流边界层厚度比值为1:2:4:8。研究结果表明,湍流边界层厚度对自由剪切层的发展、空腔底部静态压力分布、脉动压力及空腔流动类型均有重要影响,且随着边界层厚度的增大,下游剪切层覆盖的范围会增大,但是剪切层增长率降低;空腔前后静态压力压差减小、压力梯度下降;腔内局部测点的脉动压力声压级下降,各阶声压峰值频率向低频方向偏移;空腔流动类型往开式流动方向转换。     

混合动理学通量WENO方法拓展

在马赫数2.5来流条件下，开展了高频微秒脉冲放电控制压缩折角激波/边界层干扰非定常性的风洞实验，放电位于压缩折角上游沿流向布局的6对电极之间，所选取的放电频率为14 kHz,接近于来流边界层的特征频率。采用高速纹影成像技术记录流场的动态变化，并基于纹影图像灰度值的时间序列采用平均、均方根、本征正交分解、动态模态分解、傅里叶变换等方法进行处理，对比研究有/无控制情形下激波/边界层干扰的非定常特性。研究发现，对于无控情形的基准流场，流动的低频特性表现为分离激波的振荡及边界层大尺度涡经过激波的脱落行为，中、高频特性表现为边界层小尺度涡与激波的相互作用；对于受控情形，来流边界层内的大涡尺度在放电作用下增大，大尺度涡与分离激波相互作用使得激波的振荡转变为稀疏压缩波的脉动，流经激波的边界层脉动更强，分离激波的低频振荡(10～300 Hz)有所改善。此时，流动的低频特性主要表现为边界层大尺度涡经过激波的脱落行为，而中、高频特性与基准流场相似。     

基于 DDES 算法的有扰流片腔体气动噪声分析北大核心CSCD

为了进一步了解腔体流场的物理特性以及锯齿形扰流片这种被动控制方法对腔体噪声振荡的抑制效果,对马赫数0.85下、长深比为5∶1的矩形腔体进行了气动声学分析。分析过程采用了基于Spalart-Allmaras单方程模型的延迟分离涡模拟(DDES)的CFD方法。计算结果以声压级(SPL)形式与英国QinetiQ基地进行的实验研究结果进行比较。计算准确地预测出了第二阶和第三阶主导模态的声调噪声幅值,其误差均不超过2dB。为了验证腔体前缘安装扰流片这种广泛应用的腔体噪声被动控制方法,计算重点分析了带有锯齿形扰流片的腔体结构噪声等级。计算分析结果与空腔结果对比分析后认为前缘锯齿形扰流片可以降低第二阶主导模态的峰值噪声等级10dB以上,其它频率下的腔体噪声等级可降低接近5dB,降噪效果较好。     

基于DDES算法的有扰流片腔体气动噪声分析

为了进一步了解腔体流场的物理特性以及锯齿形扰流片这种被动控制方法对腔体噪声振荡的抑制效果,对马赫数0.85下、长深比为5∶1的矩形腔体进行了气动声学分析。分析过程采用了基于Spalart-Allmaras单方程模型的延迟分离涡模拟(DDES)的CFD方法。计算结果以声压级(SPL)形式与英国QinetiQ基地进行的实验研究结果进行比较。计算准确地预测出了第二阶和第三阶主导模态的声调噪声幅值,其误差均不超过2dB。为了验证腔体前缘安装扰流片这种广泛应用的腔体噪声被动控制方法,计算重点分析了带有锯齿形扰流片的腔体结构噪声等级。计算分析结果与空腔结果对比分析后认为前缘锯齿形扰流片可以降低第二阶主导模态的峰值噪声等级10dB以上,其它频率下的腔体噪声等级可降低接近5dB,降噪效果较好。     

机翼尾迹脱落涡噪声的理论研究

Ｍｉｃｈｅｌ等人１９９８年应用平面传声器阵列对飞机过顶噪声进行的测量研究首次发现，机翼尾迹脱落涡噪声是某些类型飞机重要的噪声源。为发展一种预测这种噪声源的理论预测模型，应用ｖｏｎ Ｋａｒｍａｎ涡街模型模拟二维机翼下游尾迹脱落涡，尾迹涡的强度和脱落频率应用这个模型进行计算。基于Ｈｏｗｅ后缘噪声理论，并结合尾迹模型，本文发展了一种预测脱落涡噪声声压级和指向特征的气动声学模型。对６架现代商用飞机的机翼尾迹脱落涡噪声的计算表明，本文理论模型预测的涡脱落频率、声压级以及噪声的指向性等与实验测量结果有较好的一致性。     

螺旋桨滑流对带后缘襟翼机翼气动特性影响的数值分析

采用动量理论结合有限体积法求解N-S方程的方法,进行了螺旋桨滑流对带后缘襟翼机翼气动特性干扰的数值模拟分析.螺旋桨模型简化为一个激励盘,即一个无厚度可穿透的圆盘,以给定的圆盘拉力载荷,即压强差来模拟螺旋桨滑流的效果.采用嵌套网格技术处理桨盘与机翼间复杂的几何关系.算例计算表明,该计算方法能够较好的模拟滑流对带后缘襟翼机翼气动特性的影响.     

平流层飞艇尾部形状对气动阻力的影响

为了研究平流层飞艇尾部动量边界层厚度与尾涡结构,应用LES(大涡模拟)方法计算了零攻角工况下飞艇绕流场,并对LOTTE和M-LOTTE两种飞艇进行了对比分析.采用Q分布和涡量描述回转体尾涡结构,根据Q分布可以确定M-LOTTE飞艇较LOTTE飞艇尾部分离区显著减小;并分析了回转体的轴对称曲面动量边界层厚度对飞艇气动阻力的影响,随着飞艇尾部厚度逐渐减小,动量边界层厚度逐渐增大,M-LOTTE飞艇尾部动量边界层厚度明显小于LOTTE飞艇.飞艇尾部动量边界层厚度分布说明了M-LOTTE飞艇的总阻力系数较LOTTE飞艇降低17.2%的原因,同时也表明飞艇尾部形状对飞艇气动阻力影响较大.      

超声速空腔流激振荡与声学特性研究

基于高速风洞试验研究了超声速时空腔流激振荡与声学特性.试验马赫数为1.5,基于每米的雷诺数为2.26×107,来流边界层厚度为0.024 m,试验空腔长深比分别为15,12和6.结果表明:空腔内形成的剪切层与腔后壁相撞诱发腔内较强烈噪声,噪声从腔后缘向腔前缘传播时受到腔内流动的干扰,故同频率下腔后缘处的声压均高于腔前中部区域的声压.闭式和过渡式空腔长深比较大,剪切层与腔底面相撞在腔内形成的压缩波或激波,干扰了从腔内声波反馈回路、限制了流激振荡的形成,故腔内未出现明显的声压峰值激振频率;开式空腔长深比较小,剪切层直接跨过空腔中部与腔后壁相撞,产生的噪声向腔前缘传播,腔内形成流激振荡,并出现多个声压峰值激振频率.      

基于PaSR模型的低旋流燃烧大涡模拟

为了深入理解低旋流流场特征和燃烧稳定性,基于OpenFOAM平台,采用动态k方程模型和有限速率PaSR模型对甲烷/空气预混低旋流燃烧进行了大涡模拟,研究了气流入口速度、当量比和压力等流场参数对流场结构和燃烧非稳态特性的影响,分析了流场大尺度结构与火焰相互作用。结果表明,流场结构和火焰抬升高度受入口速度影响较小,流场和火焰形态能够保持自相似性;随着当量比和压力提高,流场扩张性增强并在燃烧区下游产生回流区,火焰稳定不依赖回流区,根部火焰锋面形状由U形转变为W形,火焰抬升高度降低。火焰锋面稳定在剪切层,剪切层产生的周期性有序涡结构引起当地流场速度脉动和火焰表面褶皱,反映了流场非稳态特性;通过剪切层监测点瞬时轴向速度分析,涡结构特征频率随速度增大而提高,由250Hz提高至300Hz,随当量比和压力提高而降低,由250Hz降低至125Hz。     

扰动波相位差对剪切层中相干结构的影响

当自由剪切层中存在多个扰动时,扰动波的相位差对大尺度相干结构的发展具有一定的作用。本文用数值模拟的方法对此进行了探讨,发现相位差的变化可使大尺度相干结构的演化表现为旋涡的合并或撕裂,并在一定程度上影响剪切层动量厚度的发展。     

有/无尾迹作用下低压涡轮叶栅分离边界层转捩的大涡模拟

采用经过大量算例验证的可压缩大涡模拟求解器对雷诺数为60154、马赫数为0.402的低压涡轮叶栅T106D-EIZ进行了细致模拟,计算了定常来流和周期性尾迹来流两种工况.对计算结果的分析表明:定常来流工况下,叶片吸力面后部出现大尺寸的层流分离泡,分离剪切层的转捩过程受Kelvin-Helmholtz （K-H）不稳定性控制;尾迹来流工况下,由于来流尾迹的周期性扫掠,时均分离泡尺寸变小,叶栅总压损失降低.对相位平均和瞬态流场的分析表明,尾迹引起的逆射流使分离点后移,形成卷起涡结构,逆射流掠过卷起涡的过程中与其发生强烈的相互作用,产生大量气动损失,而后卷起涡破碎,流动转捩为湍流.      

涡发生器控制平板边界层分离的大涡模拟

为了研究涡发生器(VGs)间距λ对控制边界层分离效果的影响,选取了4种涡发生器间距,λ/H(H为涡发生器高度)分别为5,7,9,11.采用大涡模拟(LES)方法对带逆压梯度的平板边界层分离流动及VGs控制分离流动进行了数值模拟.分析了有无VGs控制时,湍流场中大尺度相干结构及其演化规律,分别从旋涡间距、边界层内流体动能、压差损失等方面考察了VGs间距对控制流动分离效果的影响.研究结果表明当λ/H为5时,VGs间距过小抑制了旋涡的展向发展,λ/H为9,11时,VGs间距过大边界层内流体动能偏低,当间距λ/H为7时流动控制效果更优,此时计算域压差损失最小,相比较无VGs控制时,压差损失降低了30.95%.      

反旋喷嘴进气旋转盘腔压力损失特性

针对旋转盘腔的压力损失问题,利用反旋喷嘴改变进入旋转盘腔的气流的切向速度,减少旋转腔内涡的生成,从而减小气体流经旋转盘腔的压力损失.通过实验的方法测得反旋进气条件下的流阻特性,应用k-ω SST模型模拟了不同工况下腔内的流场.研究表明:旋转盘腔内的压降由湍流参数和掺混后旋流系数共同控制.喷嘴自身压降所占比例为极大值时,腔内无大涡出现,而在非极值点的工况下,腔内均为对涡的流动结构.      

基于简化射线模型的剪切层相位修正方法

给出了基于简化射线模型的三维有限厚度风洞剪切层相位修正方法,利用数值模拟结果对该方法的物理依据进行了解释,并将其与传统的Amiet方法和平均马赫数法进行了对比研究,通过风洞试验对简化射线法的有效性进行了验证,最后给出了简化射线法对空间弯曲的剪切层的相位修正方法。结果表明:在来流马赫数小于等于0.3、测量角在40°~140°的范围内,当声源和传声器在同一个平面内时,对于风洞剪切层引起声波相位畸变修正,简化射线法、平均马赫数法和Amiet方法之间没有明显的差别,三种方法之间的相对误差小于1%。      

用于可压缩自由剪切流动的湍流混合长度

抓住可压缩流动变密度特性,构造出基于有效涡量的三维von Karman混合长度。湍流模型采用仅依赖湍动能k的单方程KDO(Kinetic Dependent Only)模型,引入新构造的混合长度替换旧尺度得到CKDO模型。为了验证其描述可压缩自由剪切湍流的能力,选择无壁面束缚、密度梯度大和可压缩效应强的自由剪切混合层为算例,其对流马赫数Ma_c=0.8。计算结果表明,KDO模型对混合层的速度分布有着良好的控制和模拟,而经可压缩修正后的CKDO模型与原模型及其他可压缩修正模型相比,所计算的速度分布、主雷诺剪切力和混合层厚度与试验结果更加接近,说明了该混合长度对可压缩混合层这种自由剪切湍流有着良好的刻画能力。