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

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

旋流数对分层旋流流场影响的大涡模拟研究

为了深入理解分层旋流流场特征和燃烧稳定性,采用OpenFOAM对分层旋流燃烧器的冷态和燃烧流场进行了大涡模拟。研究了旋流数对分层旋流流场结构和非稳态特性的影响。采用Q准则显示了流场中的瞬时涡结构;利用功率谱分析了流场中的进动特征。结果表明：在冷态工况下,旋流对回流区的位置和大小影响较小。随着旋流数增大,出口气流受到旋流诱导的离心作用,流动发散,流场扩张角变大,流场下游出现二次回流区。平均流场的三维流线与螺旋涡在空间中均表现成正交关系,表明螺旋涡是由剪切层Kelvin-Helmholtz不稳定性产生。在燃烧工况下,随着旋流数增大,回流区的面积增大,平均温度分布不断沿径向扩张,火焰锋面脉动增强,涡旋发生破碎的位置明显向上游移动。     

分层旋流燃烧器冷态流场的大涡模拟

针对旋流数为0.25,0.45和0.79三种工况下带有中心钝体的分层旋流燃烧器的冷态流场进行了大涡模拟(LES),选取动态Smagorinsky涡黏模型作为亚格子尺度的湍流模型,研究旋流强度对钝体回流区、涡旋破碎和进动特征的影响.模拟结果与实验结果符合得较好.模拟结果表明:3个旋流数下钝体回流区的大小没有明显改变,轴向长度都约为20mm.Q准则用来显示涡旋结构,结果表明螺旋涡产生于旋流剪切层的Kelvin-Helmholtz不稳定性;增大旋流强度,涡旋破碎发生的位置向上游移动.功率谱密度(PSD)表明流场出现进动特征,进动运动沿流向逐渐衰减;旋流数为0.45和0.79时,钝体回流区末端出现进动特征;3个旋流数下,进动频率都约为78Hz.     

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

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

高雷诺数湍流横侧射流的大涡模拟

采用大涡模拟(LES)方法研究横向主流与壁面射流均为高雷诺数的壁面横侧射流(JICF),采用动态Smagorinsky涡粘模型对亚格子尺度进行封闭。数值模拟结果中,近场和远场的平均速度、均方根值均与实验结果符合很好。瞬时压力等值面在射流迎风侧与射流迹线成正交关系,表明JICF近场迎风涡是由于射流出口上游剪切层Kelvin Helmholtz不稳定性引起的。基于长度尺度RD,Broadwell和Breidenthal提出了射流迹线公式,通过拟合数值模拟结果得到公式中常数A=1.35,B=0.3。采用平均流线、平均场涡量与Q准则分析三维涡旋结构。横向主流遇到射流后发生夹带现象,在射流出口下游受涡旋结构的影响发生卷起;平均涡量等值面表明,反向旋转涡对(CVP)的涡旋方向与射流出口下游马蹄涡两个分支的涡旋方向相反;由于横向主流和射流雷诺数较高,Q准则表征的涡旋结构在射流出口下游一倍射流直径位置开始发生破碎。     

燃料分级对旋流火焰的动态特性影响研究---“两机”专刊

为深入研究分级旋流火焰特性,本文以分级旋流模型燃烧室为研究对象,对四个不同燃料分级比(FSR)条件下的分级旋流火焰进行了数值研究,在时均燃烧场特性分析的基础上进一步对燃料分级比为1和3两个工况进行了大涡模拟(WMLES)研究。结果表明：燃料分级比的改变会影响中心回流区(CRZ)的长度和宽度。燃烧室中截面的散点分布图能够显示出不同燃料分级比条件下的燃烧特征。燃料分级比为1时,燃烧室剪切层仅存在零散的涡破碎区;而燃料分级比为3时,伴随涡破碎区还出现了单螺旋分支进动涡核(PVC)。通过FFT变换获得的燃烧室内剪切层速度能谱主频与进动涡核的旋转频率相同,表明内剪切层速度脉动的产生与进动涡核有关。另外进动涡核会使流场内的燃料分布和燃烧模式发生周期性的变化,进而影响燃烧过程。调整燃料分级比在1附近,能够使分级火焰达到稳定燃烧降低排放的目标。     

文氏管出口张角对旋流杯流动不稳定性的影响

为探究旋流杯内流动不稳定现象,针对带不同出口张角文氏管的旋流杯冷态流场进行了大涡模拟,并结合PIV试验结果进行了验证。研究结果表明,文氏管出口张角对时均流场的影响主要在旋流杯内部区域,对下游恢复区域基本无影响;旋流杯瞬态流场中存在旋进涡核（PVC）和小涡旋结构,出口张角对PVC的运动频率和形态基本无影响,对小涡旋结构影响较大,相同Q值下,56°张角旋流杯方案的小涡旋结构更丰富;经本征正交分解方法（POD）分析,旋流杯的强脉动结构为PVC和小涡旋结构。56°张角旋流杯方案的小涡旋结构能量更强,脉动幅度较大,能更好地对燃油进行剪切破碎,有利于强化雾化效果。     

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

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

超声速高度欠膨胀冲击射流的大涡模拟

冲击射流广泛应用于短距起飞垂直降落飞行器(SVTOL)等航空航天领域.本文采用大涡模拟方法对高度欠膨胀的超声速冲击射流的流场进行了数值模拟.本文数值模拟得到了高度欠膨胀冲击射流流场中的激波结构和内外剪切层中不同尺度的涡结构.数值结果观察到了马赫盘的振荡,以及在斜激波、马赫盘及大尺度涡结构的共同作用下,射流内外剪切层之间的环形激波的生成与消失的周期过程.并对流场内剪切层的涡结构的演化进行了研究,数值结果显示内剪切层的大尺度涡结构的形成与马赫盘的振荡相关,在内外剪切层的作用下形成了壁射流区内外交错的涡结构.     

强欠膨胀圆管射流中二次涡环的形成机理

基于大涡模拟方法与高精度混合格式,对出口压力比(n)为1.4与4.0的弱与强欠膨胀圆管射流初始流场进行了数值模拟。结果清晰地描述了两种压力比条件下欠膨胀射流初始流场的结构特征,包括涡与激波结构及它们演变过程,发现强欠膨胀射流中二次涡环是由三波点向下游形成的滑移层失稳而直接卷起,此与弱欠膨胀情况不同。另外,还讨论了二次涡环在绕主涡环外缘向后翻转过程中的演变过程及其对主涡核的影响。当n=4.0时,得到射流域内桶形激波长度与马赫盘直径分别为24mm与0.84mm,与经验公式的计算结果相符。     

吹吸扰动对壁湍流边界层摩擦阻力的影响

对风洞中平板湍流边界层施加局部周期性吹吸扰动,并用IFA-300热线风速仪测量了其近壁区域不同法向位置瞬时速度的时间序列信号.利用湍流边界层对数律平均速度剖面,计算湍流的壁面摩擦阻力.对平板湍流脉动速度信号用子波分析进行多尺度分解,通过条件采样以及相位平均的方法获得不同尺度下的相干结构波形和能量分布,来研究吹吸扰动对平板湍流边界层的摩擦阻力的影响.结果表明,扰动使得窄缝附近流动阻力增加,而下游区域流动阻力减小.      

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

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

伴流速度对平行喷口射流影响的数值研究

采用大涡模拟方法计算研究平行喷口出口马赫数为0.9,伴流速度比分别为0.1、0.3和0.5时的喷射流场特性。计算时采用高精度的数值模拟方法,并结合Smagorinsky亚网格尺度模型。考察了流场统计平均特性、脉动特性以及射流流场中涡结构的发展演变过程,结果表明:伴流速度的增大使得势流核长度变长,减缓了空间剪切层的发展,转捩延迟。喷射流场速度分布具有自相似性,而湍流强度的分布则不具有相似性。通过分析剪切层中轴向速度脉动、径向速度脉动、压力脉动在空间任意两点上的时空相关性,发现随着伴流速度的增大,脉动量在空间上的相关性减弱,而脉动量向下游的传递速率增加。该研究结果为进一步揭示伴流速度对喷流声场的影响提供基础。     

火花型激励合成射流瞬时流场测试

采用粒子图像测速仪(PIV)相位锁定采样方法对一个特定结构的火花型合成射流的非定常流场进行了实验测量，得到了火花型合成射流形成的涡对在外场的发展过程；同时采用热线风速仪对喷孔下游固定位置的瞬时速度进行了测试，对火花型合成射流激励参数(放电器储能、激励频率)的影响进行了初步的对比分析。结果表明:合成射流在孔口喷射初期呈现球面扩散状，在发展过程中不断剪切和卷吸外部气流，形成一系列旋向相反的涡串；合成射流的最大瞬时峰值法向速度出现在放电后大约005T时刻，随后其作用范围迅速扩大，瞬时峰值速度和涡量逐渐降低。在该激励器结构参数和激励参数范围内，合成射流随激励频率和放电器储能的增大而有所增加。      

Precessing motion in stratified radial swirl flow

Vortex/flame interaction is an important mechanism for unsteady combustion in a swirl combustion system. Technology of low emission stirred swirl(Te LESS), which is characterized with stratified swirl flow, has been developed in Beihang University to reduce NOXemission. However,large-scale flow structure would be induced in strong swirl flow. Experiments and computational fluid dynamics(CFD) simulation were carried out to investigate the unsteady flow feature and its mechanism in Te LESS combustor. Hotwire was firstly applied to testing the unsteady flow feature and a distinct mode with 2244 Hz oscillation frequency occurred at the pilot swirl outlet.The flow mode amplitude decayed convectively. Large eddy simulation(LES) was then applied to predicting this flow mode and know about its mechanism. The deviation of mode prediction compared with hotwire test was 0.8%. The spiral isobaric structure in pilot flow passage indicates that precessing vortex core(PVC) existed. The velocity spectrum and phase lag analysis suggest that the periodic movement at the pilot outlet was dominated by precessing movement. Negative tangential momentum gradient reflects that the swirl flow was unstable. Another phenomenon was found out that the PVC movement was intermittently rotated along the symmetric axis.