多孔介质对圆柱-翼型干涉噪声的影响

采用LES（large eddy simulation）+FW-H（Ffowcs Williams-Hawkings）方程，研究了圆柱表面使用不同PPI（pore number per inch）和厚度的多孔介质对圆柱尾迹及圆柱-翼型干涉噪声的影响，探索了多孔介质的降噪规律和机理。结果表明：多孔介质能稳定圆柱表面的剪切层，抑制旋涡脱落，从而削弱尾迹对下游翼型的影响，圆柱单音噪声最大可降79 dB，翼型单音峰值降低13.22 dB，宽频噪声降低20 dB；多孔材料PPI的变化对降噪效果影响较小，而厚度是影响流场模态、降噪效果及气动性能的一个关键参数；多孔材料厚度合适时，圆柱-翼型流场形态为“剪切层模态”，可有效降低湍流干涉噪声；多孔材料厚度较小时，发现了一种流场形态，即“剪切层-尾迹模态”，导致翼型噪声增大；合适的多孔介质厚度不仅降噪效果显著，对圆柱-翼型的气动性能也有改善作用。     

圆盘绕流近尾迹实验研究

采用烟线和PIV实验研究Re=22400时直径与厚度比D／H=5圆盘近尾迹（x／D〈5）流场结构。流向平面烟线显示表明圆盘尾迹有稳定的回流区和随机扭曲倾斜的三维涡旋结构。对烟线平面PIV测速数据统计平均,发现该回流区具有较好的对称性,且长度为2．1D。对PIV数据进行POD重构,发现流向雷诺正应力和切应力峰值出现在回流区两侧剪切层,横向雷诺正应力峰值出现在回流区驻点附近,其正是烟线显示涡旋结构脱落区域。圆盘尾迹涡旋结构的产生和脱落源于剪切层的不稳定性及其与回流区的相互作用。圆盘尾迹前两个模态含能仅为10．9％和10．2％,表明湍流结构的随机性;前10个模态含能45．9％可较好描述流场的平均雷诺切应力。     

横流中空心锥形喷雾的扩散特性与影响因素

为了更好地理解和认识空心锥形喷雾射流与横向气流的掺混过程,基于可视化实验测量,针对空心锥形喷雾的初始雾化状态、喷射角度和雾化锥角以及横流速度等因素对空心锥形喷雾射流在横流中的扩散特性进行了系统分析。研究结果表明,随着横流速度的提高,流场中对称反旋涡对（CVP）结构变小,但其稳定性及产生的卷吸气流则是先增强而后减小;提高喷雾初始液滴的动量和数流率,均可以增大射流剪切层的贯穿深度和射流尾迹的伸展高度,同时使流场中CVP的涡量强度增大。基于实验测量结果,建立了空心锥形喷雾垂直入射横流条件下掺混流场液滴群CVP结构特征尺寸和剪切层轨迹的预测关联式。另外针对喷嘴雾化锥角和喷射角度的分析表明,当喷雾初始雾化状态相近时,随着雾化锥角的减小,流场中CVP的水平尺度减小但竖直方向尺度增大且结构更为稳定,同时喷雾液滴的贯穿深度增大;相比喷雾垂直于横流入射,当喷嘴以一定角度逆向横流入射时,CVP结构稳定性减弱,流场剪切层涡结构变大且剪切层区域液滴富集现象减弱,射流尾迹紊乱程度增加,反之,则流场剪切层涡拟序结构变小,射流尾迹现象减弱,CVP结构变小。     

圆柱跨声速绕流中的激波/湍流相互作用大涡模拟研究

采用大涡模拟方法研究了圆柱跨声速绕流中的激波/湍流相互作用问题,来流马赫数M∞取为0.75,基于圆柱直径D的雷诺数为2×105。计算结果表明,圆柱分离点处出现一道斜激波,并且以与涡脱泻Strouhal数一致的特征频率向上游传播。激波运动导致流场中出现反对称的流动模态,剪切层中压力信号的功率谱曲线中存在0.4、-1和-5次方的斜率关系,剪切层中的剪切应力角约为0°,脉动速度以流向脉动速度为主,并且沿剪切层的大尺度结构组织性减小。     

尾流撞击效应对轴流压气机下游叶排涡面的影响

针对某高负荷高压压气机末级,通过研究上游转子尾迹涡面周期性扫掠对下游静子涡面的影响,从涡面非定常作用的角度分析流场时空结构.结果表明,流场性能的改善与涡面时空结构相对应,尾迹涡面同分离涡面的相互作用可显著降低流动损失;与上游尾迹涡面发生耦合的是叶片尾缘旋涡交替脱落频率而非吸力面前缘分离涡不稳定频率.      

基于DDES方法的叶栅分离旋涡的非定常流动数值研究

为了研究叶栅分离旋涡的非定常流动特性,选取亚声速叶栅为研究对象,进行了基于延迟分离涡模拟方法 (DDES)的非定常数值模拟。针对进口Ma0.3的情况,对比分析了进口攻角、叶栅稠度、叶片弯角、叶片中弧线弯度分布和叶片厚度等因素对于叶栅分离旋涡运动的影响。结果表明,当攻角增大时,尾迹附近流场会逐渐由定常模式转化为非定常对涡脱落模式,同时对涡脱落频率不断减小;随着攻角的进一步增大,叶背分离强度不断增加,叶背分离涡开始单独脱落,并主导流场旋涡运动。强烈的叶背分离涡脱落引起了叶片气动力的剧烈脉动,其幅值是小攻角状态下由尾迹对涡脱落引起叶片气动力脉动的30倍以上。叶栅几何参数的改变对于流场旋涡运动同样有很大影响,流场同样呈现出"定常尾迹"、"对涡脱落"、"叶背分离涡脱落"这三种主要的旋涡运动形式之一。叶背分离点的相对位置则是影响旋涡非定常运动形式和旋涡脱落频率的主要内在因素。     

MVG控制斜激波/湍流边界层干涉的大涡模拟

采用浸没边界法(IBM)对带有微型涡发生器(MVG)控制器的激波/湍流边界层干涉流动进行了大涡模拟(LES)。以来流马赫数为2.3的斜激波(由平板上方8°楔产生)为基本流动入射平板湍流边界层,通过在干涉区前布置MVG阵列来控制激波诱导的边界层分离。采用浸没边界法处理MVG的复杂几何,分析了MVG尾迹区平均流速度剖面,雷诺应力,瞬态旋涡结构。结果表明:时均流场显示MVG尾迹区存在一对对转的主流向涡,流向涡加剧了边界内的动量交换从而增加了边界层抗分离能力,而瞬态流场则反映出MVG尾迹区的剪切层由于Kelvin-Helmholtz(K-H)不稳定性会卷起为一列展向旋涡。     

低速轴流压气机转子尾迹三维动态测量

对一台低速轴流压气机的转子尾迹进行了三维动态测量,采用锁相平均、方均根和功率谱密度等分析方法对低速轴流压气机转子在节流过程中尾迹流场的稳态和动态特征进行了研究.结果显示:在绝对坐标系下,转子尾迹有高总压和高方均根特征,且在不同叶高的强度不同;转子尾迹与主流的掺混作用导致平均总压升系数的快速衰减;功率谱密度分析表明,在接近失速点时,转子尾迹90%叶高以上出现0.5倍叶片通过频率的特征频带且其功率随着流量减小而增大.      

非定常来流压力下基于DMD方法的预冷器换热特性

采用大涡模拟方法对非定常来流压力条件下叉排管束预冷器换热特性进行了研究,同时运用动力学模态分解方法对流场主控流动结构进行了识别,探讨了来流压力周期性变化频率对预冷器内部流动、换热性能和熵产的影响。结果表明:来流压力变化频率对预冷器时均和瞬态换热性能影响均不显著,但当来流压力变化频率增大至流场固有频率950 Hz时,流场发生共振,换热性能发生剧烈振荡;管束壁面剪切层运动和绕流脱落涡结构为主控流动结构,其时空演化过程对瞬时换热性能起决定作用;当流场发生共振时,剪切层的生长和演化与来流速度的脉动密切相关,前排管束的绕流涡脱落周期与来流压力/速度变化周期一致,而壁面剪切层的生长周期则为来流压力/速度变化周期的两倍。此外,叉排管束流场的换热熵产决定于主控流动结构,其时空演化特征与主控流动结构演化规律完全一致。     

基于合成射流的二维后台阶分离流主动控制

作为一种新的流动控制激励器,合成射流技术在流动分离控制、降低压力脉动和抑制噪声等方面具有广阔的应用前景。实验利用合成射流主动控制技术对二维后台阶湍流分离再附流动控制进行了研究,应用表面测压、粒子图像测速(PIV)和热线等多种实验测试技术对后台阶表面压力分布、流场结构以及剪切层特性进行了测试。结果表明,在台阶前缘施加合成射流可有效减小回流区范围和降低再附长度,当合成射流的动量系数为0.301×10^-3时,可使再附点长度减小25%。合成射流控制使得沿台阶下游的湍动能和雷诺应力增强,提高了台阶下游流场的混合效率。热线动态结果表明频率是后台阶分离流动控制的关键参数,当频率为260 Hz、激励频率与剪切层涡脱落频率之比为1.32、激励频率等同于旋涡脱落频率时,合成射流控制效果最好,仅需消耗较小的能量即可实现流动控制的目的。     

汽车后视镜非定常流场的大涡模拟

以桑塔纳3000后视镜为例,使用大涡模拟方法,研究汽车后视镜非定常流场。通过分析后视镜表面平均压力因数、中截面平均速度分布发现,平均流场的大涡模拟结果与雷诺平均结果相当一致,两者可以相互验证。比较有无后视镜两种情况,紧邻后视镜的汽车侧窗压力因数标准差值以及后视镜尾部速度脉动标准差值,了解到后视镜不仅自身较大压力脉动,而且影响汽车侧窗和其它区域。后视镜的剪切层以及尾部旋涡脱落导致后视镜尾部速度脉动在个别位置达到来流0.2倍。测点的自功率谱密度以及相关分析有助于从更深角度明确压力、速度脉动的能量分布以及后视镜的非定常流场结构特征,从而更好揭开后视镜导致气动噪声的流动机理。     

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

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

基于大涡模拟的圆转矩喷管尾喷流强化掺混机理研究

为揭示2元圆转矩喷管尾喷流强化掺混的内在机制,应用大涡模拟(LES)方程对2种相同进、出口直径的喷管模型(轴对称、2元圆转矩)在Ma=0.8、高雷诺数(2×10~5)条件下进行了数值模拟计算。结果表明:与轴对称喷管相比,圆转矩喷管射流掺混效应增强,速度衰减快,核心区长度和高温区域面积减小。同时尾喷流拟序结构变化说明:2种喷管主要拟序结构均包含涡环、涡辫、发卡涡、螺旋涡等相似结构;但圆转矩喷管在射流近场诱导出的涡旋更丰富,边角剪切涡发展更快,形成明显的CVP结构,导致其射流柱失稳时刻提前、距离缩短;同时,喷管形式的改变使得射流剪切层内雷诺应力增大,速度脉动增强。拟序结构发展及雷诺剪切应力变化说明在射流流场中涡旋发展耗散速度增大、速度边界层脉动增强、射流柱易失稳是导致射流掺混增强的本质因素,为异形喷管的强化掺混机理提供了依据。     

合成双射流控制水下圆柱绕流流动分离数值模拟研究

为了探究合成双射流激励器及出口射流参数对圆柱绕流流动分离的控制效果,首次对合成双射流控制水下圆柱绕流流动分离进行了数值模拟。数值计算结果显示:保持激励器出口射流振幅不变的条件下,出口射流频率等于尾迹涡脱落特征频率时,射流控制作用与绕流流场耦合效果最好,控制流动分离效果最佳;保持出口射流频率为尾迹涡脱落特征频率时,在数值计算测试范围内,随着射流振幅的增大,射流对于流场的动量掺混能力增强,控制效果也随之增强。机理分析表明:合成射流位于前驻点的控制,主要通过在圆柱前缘形成虚拟气动外型来达到减阻控制的效果;而合成射流位于后驻点的控制,主要是通过增强回流区的动量掺混来提高回流区抑制分离的能力,从而达到减阻控制的效果。     

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

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

趋于临界马赫数的圆柱跨声速绕流特性分析

通过深化认识趋于临界马赫数Ma_cr的圆柱跨声速绕流特性,明确新型飞行器增升减阻设计的空气动力学理论依据。采用大涡模拟方法数值研究了来流马赫数Ma_∞为0.75和0.85、雷诺数Re为2×10⁵的圆柱跨声速绕流。结果表明:当Ma_∞趋于临界马赫数(Ma_cr≈0.9)时,圆柱的阻力下降且升力系数振荡被抑制;通过力的分解,得知圆柱的阻力减小来自旋涡力的影响,而非可压缩性;圆柱的阻力减小与其背压上升有关;剪切层初始阶段的对流马赫数Ma_c随Ma_∞的增加而增大,而增长率相反,这使得剪切层更为稳定、柱体背压更高。此外,由于Ma_∞=0.85时边界层分离点处的激波和尾迹处的激波向下游推移,使得近尾迹处的湍流脉动减弱,进而导致柱体的表面压力振荡和升力系数振荡被抑制。     

连续尺度变换方程在稳态及非稳态流动模拟中的应用

为了提高涡粘性假设的湍流模型对于非稳态流动的求解精度,同时兼顾其对于稳态流动的求解性能,将雷诺应力项与连续变换方程(CSSE)结合而形成新的应力项,使其根据流场尺度、网格尺度及Kolmogorov尺度来自动调节当地的应力雷诺应力模化水平,避免网格因素在流场模拟中产生不利影响,改正了混合RANS/LES方法的速度型偏离对数率问题;同时,该方法并未引入显式亚格子模型(SGS),因此回避了亚格子系数确定对于流场模拟精度产生的影响,改善了湍流模型对于流动不稳定性的辨识精度。在湍流平板算例中,CSSE方法计算的边界层速度型精度与雷诺平均方法(RANS)相当,而对于圆柱尾迹的模拟则证明了CSSE方法具有混合RANS/LES方法的优点,即能够准确模拟流动的不稳定性特征。     

支板喷射超声速湍流燃烧的大涡模拟

为了发展可行的超声速湍流燃烧大涡模拟方法,将设定型PDF(Probability Density Function)模型与LES(Large Eddy Simulation)相结合以封闭亚格子湍流-燃烧相互作用,并将模型用于支板喷射超声速湍流燃烧流场的数值模拟。分别对冷流及燃烧流场进行了模拟,计算结果与实验测量符合较好,表明了所采用方法及模型的可行性。冷流条件下,大尺度湍流涡通过卷吸、拉伸运动主导支板尾迹区的近场混合,并通过破碎过程影响远场混合。燃烧条件下,回流区尺度扩大,剪切层中形成的高温燃烧产物通过大涡卷吸以及回流区末端对流作用进入回流区并与其中的燃料喷流相互作用,使部分燃料预热升温并进入回流区两侧剪切层与主流空气混合、燃烧,从而实现火焰稳定。在薄反应剪切层及大尺度反应涡的边界区域,LES网格不足以直接求解湍流与燃烧的相互作用,PDF模型给出了较强的亚格子脉动。     

低雷诺数NACA0012平面叶栅流场直接数值模拟

采用具有7阶精度的weighted essentially non-oscillatory(WENO)差分格式,直接求解可压缩二维非定常N-S方程组,研究了NACA0012翼型平面叶栅低雷诺数流动的特征.直接模拟及与文献对比的结果表明:叶栅尾缘涡脱落的形成过程与圆柱绕流涡脱落的形成过程非常相似.平面叶栅尾迹区的2阶统计量与孤立翼型尾迹区的2阶统计量具有相同的分布特征,但前者的强度显著大于后者.周期性的涡脱落不仅在上下翼面形成非定常分离,也使得尾迹区某点的总压发生准周期性的变化.随着栅距的减小,翼型上的平均分离位置向前缘移动;尾迹区某点的总压变化频率及其幅值均显著地增加;而且栅距越小,速度脉动2阶统计量反而越大.      

大涡模拟研究钝体有旋流流场的拟序结构

对带有中心钝体的分层旋流燃烧器的冷态流场进行了大涡模拟(LES),选取动态Smagorinsky涡黏模型作为亚格子模型,研究旋流数为0.45时旋流场的大尺度拟序结构。模拟结果表明:瞬时压强等值面显示的内外螺旋涡及进动涡核(PVC)均与平均速度场流线在空间上呈正交关系,表明两种涡结构均由剪切层Kelvin-Helmholtz不稳定性产生。Q准则等值面显示内外螺旋涡在下游20mm左右开始发生破碎。PVC发源于环形旋流与环形射流剪切层附近的下游区域。瞬时周向速度的功率谱密度(PSD)出现明显的特征峰,表明PVC影响附近的流体,使之出现进动特征。采用本征正交分解(POD)重构湍流脉动速度场,不同模态下功率谱密度结果表明前两个模态的大尺度结构具有进动特征,使用前两个模态的周向脉动速度等值面显示了PVC周围流体的大尺度涡旋结构。