微缝共振器的非线性吸声机理

为模拟高声压级时微缝处涡声相互作用现象并计算微缝声阻抗,发展了一种离散涡模型.联合微缝离散涡模型和一维声传播模型,计算了微缝共振器的声阻抗和吸声系数.数值结果表明,微缝声学性质主要取决于涡脱落和脱落涡层的演化;共振器的吸声系数随着入射声压级的增加而增大;入射声波频率对共振器吸声系数有显著影响.研究结果和先前直接数值模拟...     

一种预测平板尾迹噪声的时域无网格方法

基于离散涡方法和涡声理论建立了一种预测二维平板尾迹发声的时域无网格方法。该方法应用解耦方式完成声场计算,首先使用离散涡方法计算了均匀来流作用下的平板尾迹流场,得到了流场中点涡的涡量、位置和速度等关键参数,然后基于涡声理论建立了自由空间中点涡发声模型,并引入了时域边界元方法来模拟平板表面对声场的散射作用,计算得到了平板尾迹涡发声的偶极子声场分布和指向性等关键特征。通过对上下表面涡以及平板散射对声场贡献的深入分析表明,从平板尾缘上下角点脱落并卷起的涡团均为偶极子源,平板的散射作用使得声场在一定程度上得到加强,并且使声场具有极大值方向垂直于平板表面的偶极子指向性特征。所建立的无网格方法计算快速,能同时获得流场和声场分布的关键特征,可提升对气动噪声产生机理的基本认识,同时还为尾迹噪声的理论研究提供了一种具有工程应用价值的可靠计算方法。     

钝体绕流有旋流中回流区与进动涡核的大涡模拟

针对旋流数为0.57、0.68、0.91和1.59四种工况下的悉尼旋流燃烧器的冷态流场进行了大涡模拟,选取动态Smagorinsky涡黏模型作为亚格子尺度的湍流模型,研究不同旋流数下的流场结构、进动频率和进动涡核。模拟结果表明,旋流数为0.68时,钝体回流区长度最短。随着旋流数的增加,中心射流出口的旋流剪切层不断衰减,而下游的旋流剪切层不断增强。功率谱分析表明,进动现象的出现和消失对应于与旋流剪切层的增强和衰减;中心射流与下游区域具有不同的进动频率,表明流场中存在着两个独立的大尺度涡旋结构。不同取值位置的周向速度相关性分析进一步佐证了两个涡旋结构的存在。轴向位置70 mm处的横截面上瞬时流线和压强分布证实了下游流场存在着进动涡核。瞬时压强等值面显示了中心射流出口和下游流场进动涡核的三维螺旋形结构。下游流场的进动涡核均与平均速度场流线在空间上成正交关系,表明进动涡核是由剪切层Kelvin-Helmholtz不稳定性产生。     

旋翼非定常平行桨-涡干扰流场的数值模拟

建立了一个基于Euler方程的旋翼平行桨-涡干扰流场的数值模拟方法。该方法在时间离散上使用隐式的双时间法,以避免因物理时间步长过小而引起的数值不稳定问题;为减小由于空间离散格式精度和网格密度引起的涡数值耗散,采用了预定涡方法。以Kitaplioglu的平行桨-涡干扰试验模型为算例,计算了表面瞬时压强分布,并与可得到的试验结果相比较,验证了方法的有效性,且描述了干扰流场的气动力变化特性以及桨叶表面扰动压力波的产生和传播过程,然后分析了马赫数、涡强度、干扰距离等参数对桨-涡干扰流场特性的影响。在此基础上,得出了一些有意义的结论。     

大涡模拟模型环形燃烧室污染特性

在任意曲线坐标系下,采用大涡模拟模型环形燃烧室两相燃烧流场中污染物的生成.用亚网格涡破碎(EBU)燃烧模型估算化学反应速率;分别用NO亚网格动力学模型与CO亚网格模型来模拟NO和CO的生成;采用随机离散模型模拟气液两相湍流流动,两相之间耦合采用PSIC算法.计算结果与瞬态速度场、出口温度和污染物质量分数分布的实验数据比较吻合,表明大涡模拟方法和亚网格模型可以用来预估实际燃烧室两相喷雾燃烧流场和污染物生成过程.     

后视镜罩边缘结构对流场和气动噪声的影响

采用大涡模拟(LES)的计算方法,对车外后视镜不同边缘结构引起的外部流场和镜后车身表面监测点处的气动噪声进行了数值仿真.研究表明:不同的镜罩边缘结构在较大程度上影响了流经后视镜罩的气流速度和流线方向,对后视镜后部流场和监测点处声压级产生较大影响.相较于原模型的光滑边缘结构,模型1后部流场涡团更远离车身表面,有利于降低气动噪声,某监测点1/3倍频程中心频率处声压级最大降幅接近10dB;模型2后部流场涡团产生分离,并更靠近车身表面,反而使气动噪声增大.模型1的镜罩边缘结构改进方案对车外后视镜-A柱区域的流场和气动噪声都有较好的改善.     

双缝激波诱导推力矢量喷管的非定常流动特性研究

采用分离涡模拟方法对双缝激波诱导推力矢量喷管的三维流场进行数值模拟,分析原始流场的流线、压力系数分布。运用傅里叶模态分解技术对喷管对称面压力系数进行模态分解和时间演化。结果表明:双缝激波诱导推力矢量喷管的流场中包含两个模态,一阶模态主要表现为下游二次流口激波处和尾流上剪切层处的压力振荡,二阶模态主要表现为剪切层处涡的脱落。本文对研究矢量喷管非定常流场特定频率下对应的潜在动力学信息具有重要意义。     

基于火焰面模型的超声速燃烧混合LES/RANS模拟

为了明晰超燃冲压发动机燃烧室内部燃烧过程的细节,建立了超声速湍流燃烧稳态火焰面亚格子模型,并采用混合LES/RANS方法对氢燃料超燃冲压发动机进行算例验证.控制方程对流项用五阶精度WENO格式离散,时间方向采用二阶Runge-Kutta方法.研究表明:(1)冷流流场中燃料分布与大尺度结构分布相似,说明混合过程受大涡控制;(2)燃烧流场中涡的尺寸明显变大,且仅存在于火焰面上,另外温度分布和主要生成物分布与涡量云图基本相同,说明燃烧过程也由大涡控制;(3)时均计算结果与实验阴影基本符合,速度剖面和温度剖面与实验测量值定性一致,说明本文的数值模拟方法和燃烧模型可以较好地描述和预测超声速流动燃烧过程.     

超声速弱欠膨胀冲击射流流场结构

为了对超声速弱欠膨胀冲击射流的流场结构细节进行研究,使用大涡模拟方法对其进行了数值模拟。利用三阶迎风和四阶对称紧致格式对无量纲化轴对称可压缩滤波N-S方程进行空间离散,时间上推进采用的是三阶精度的TVD型Rugge-kutta法。亚格子尺度模型采用的是修正Sm agorinsky涡粘性模型。通过与经典的冲击射流实验比较,证明了程序的可靠性。数值模拟得到了剪切层以及壁面射流中的涡结构和主射流中的激波结构,并且在此基础上对涡合并和板前激波和涡干扰现象进行了深入研究。发现涡合并现象主要出现在流场的上游,越往下游出现的几率越小;涡和板前激波的相互作用会引起激波位置和强度以及冲击平板上冲击区的压强的显著变化,同时也会导致涡的变形。     

分区LES/DES混合方法及缝翼三维流场模拟

提出一种分区大涡/脱体涡模拟混合方法,采用高阶空间有限差分方法和二阶隐式LU-SGS时间推进方法,对有限展长缝翼的三维流场进行了数值计算。该混合方法吸收了脱体涡模拟方法和大涡模拟方法各自的优点,在非核心区域使用脱体涡方法进行计算,相对于大涡模拟方法具有较高的计算效率;相对于分区脱体涡模拟中采用的雷诺平均方法,在计算量没有明显增加的条件下提高了对流场各尺度流动的模拟能力。计算结果与风洞实验数据吻合良好,同时缝翼流场主要观测点的压力脉动数据为后续缝翼噪声分析和低噪声优化提供了基础。     

高声强下多狭缝共振腔的吸声性能

为深入理解高入射声强下多狭缝共振腔吸声机理,采用低频散低耗散的计算气动声学方法对二维多狭缝共振腔开展直接数值模拟研究.首先对标准单狭缝共振腔计算结果进行验证,随后相同的数值模拟方法被应用于相同穿孔率的多狭缝共振腔的数值模拟中.结果显示:高声强下涡脱落对吸声系数的贡献占据了主导地位,各入射频率下均超过68%.多狭缝低频时会导致脱落涡总能量的下降,而高频时升高;而黏性耗散作用随着共振腔狭缝数目的增加而增强.因此综合作用下低频时多狭缝共振腔对吸声效果影响不大,但在高频时多狭缝共振腔有更好的吸声效果.      

大迎角翼—身组合体涡流特性计算

论述了用非线性离散涡法来模拟大迎角下翼－身组合体涡流绕流的计算模型及计算方法；给出了一个典型翼－身组合体的涡流流态计算结果及非线性气动特性和截面压强分布。     

卡门涡街与翼型的相互作用

本文利用保角变换,用离散涡方法求解了卡门涡街与翼型的相互作用。较详细地分析了它作为偶极子声源的特性,并模拟了在涡翼相互作用过程中尾迹的发生和演化。计算结果表明:在相互作用时,尾迹所吸收的声能可以是正值也可以是负值;升力和声偶极子强度出现了周期性脉动,这种脉动与翼型的振动和噪声有直接的关系。     

圆柱尾流中的二次分离及环量衰减问题

在圆柱绕流的非定常流动中,相应于旋涡的周期性脱落柱后的二次分离是周期性间断发生。二次涡的影响不能忽略。但对长时间流动的模拟还必须考虑尾涡的环量衰减。模拟Bluff Body流动有效的离散涡模型应包括二次分离及尾涡衰减两种机因。本文数值计算结果证明了上述结论。     

涡和翼型的强相互作用

基于涡云模型研究涡和翼型的强相互作用，重点模拟主涡在贴近翼型表面经过时涡变形对声偶极子强度的影响，并与单点涡模型的研究结果进行了比较。     

环形中心钝体驻涡燃烧室驻涡腔有无喷射的对比

采用三维雷诺平均Navier-Stokes(N-S)方程、renormalization group (RNG) k-ε湍流模型和标准壁面函数对驻涡腔有无喷射的环形中心钝体驻涡燃烧室的冷态流场进行了数值仿真,分析了驻涡腔有无喷射对环形中心钝体驻涡燃烧室涡系结构、驻涡腔流动参数和燃烧室总体性能的影响.结果表明:相比于无喷射时,驻涡腔添加喷射可以使驻涡腔内形成相对稳定的双驻涡结构;驻涡腔喷射的存在使得环形中心钝体驻涡燃烧室出口截面总压损失系数降低约9.2%,并能提高驻涡腔内的平均气流参数;驻涡腔喷射对环形中心钝体驻涡燃烧室出口截面流动参数沿流道高度方向的变化趋势影响不大.      

穿孔板切向流效应的理论和实验研究

 建立了一种有关切向流对小孔声阻抗影响的小扰动势流模型,其要点是使用法向质点速度连续而非广为采用的质点位移连续边界条件来匹配小孔剪切层两侧的声场,并在小孔前缘施加Kutta条件以反映切向流效应的机理———锐缘处的声涡转化。模型中还包含了实际应用中经常遇到的小孔形状和穿孔板厚度两个重要因素。理论预测的声阻和声抗与实验结果良好吻合。     

有限长圆管内刚体旋转流的直接数值模拟

通过对三维不可压缩粘性流体的直接数值模拟,研究了有限长圆管内刚体旋转流在初始扰动下,扰动的增长和动力学演化过程.首先,通过改变流动的Reynolds数和旋转角速度ω,分析了轴对称流动中Reynolds数和ω对流动失稳的影响.对于高Reynolds数轴对称流动,利用数值模拟得到的不稳定结果与无粘线性稳定性理论预测的结果一致.最后,为进一步研究流动不稳定发展的三维效应,又进行了完全三维流动的模拟.数值结果表明,三维流动中不稳定的主导模态为螺旋模态,模态在线性段指数增长,导致流动产生螺旋型的漩涡破裂,之后经过非线性段的增长后达到饱和,流动最终发展为轴对称泡型旋涡破裂的稳定状态.