乙醇-空气稀液雾值班火焰的大涡模拟研究

为加深对气液两相湍流燃烧现象的认知,检测火焰面模型在液雾燃烧中的适用性,本文在欧拉-拉格朗日架构下使用火焰面/反应进度变量模型(FPV)数值模拟了乙醇-空气稀液雾值班火焰。欧拉坐标系下的气相湍流场使用大涡模拟方法模拟,离散液相则使用拉格朗日颗粒轨道模型进行描述,考虑了相间质量、动量、能量交换。燃烧模型中采用碳元素定义混合物分数,在混合物分数方程源项中体现液相对燃烧模型的影响。模拟得到的气相温度分布和液相统计值均和实验数据较好吻合,验证了该燃烧模型对稀液雾扩散类型火焰的适用性。分析瞬时图发现,该稀液雾火焰的最高燃烧温度往往位于当量混合物分数附近,在出口下游20倍直径处火焰完全点着,此处上游FPV模型能给出局部点火熄火现象。蒸发作用在剪切层和点火区域较强,而液滴和火焰作用较弱,单一液滴很少被火焰包围。     

稳定器流动、燃烧实验及冷态大涡模拟研究

蒸发式V型稳定器是一种适用于加力式涡扇发动机加力燃烧室的外涵稳定器,本文在冷态实验中测量其流阻,同时在热态实验中测量其工作范围,并用大涡模拟对其冷态流场进行了三维全尺寸非稳态的数值模拟。分析表明,多处的回流区的组织有利于增加稳定器的稳定工作范围。      

脉冲横向射流的大涡模拟

为了研究脉冲气膜冷却的流动特性,采用基于高阶对称加权本质无振荡(WENO)-Roe格式和隐式牛顿迭代时间推进、动态亚格子模型和预处理技术下的大涡模拟(LES)计算方法,对吹气比为1.0下的斜管横向射流进行了稳态和脉冲两种不同工况下的数值模拟.计算结果显示了射流流场中存在发卡涡结构,研究表明射流孔进口位置的脉冲改变了射流出口拟序结构,即改变了射流与主流的掺混过程.      

基于大涡模拟的离心式喷嘴雾化过程研究

为实现离心式喷嘴雾化过程的精确数值仿真,探究喷嘴内部流动特性与外部液膜破碎形式,采用基于大涡模拟的仿真方法,对一种典型的四进口离心式喷嘴进行研究,仿真结果揭示了喷嘴内部相界面的振荡现象与外部液膜的破碎细节,并通过耦合流体体积法(VOF)与离散相模型(DPM),获得液滴粒径的空间分布特征.研究结果表明:在液体填充过程中,...     

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

为了深入理解分层旋流流场特征和燃烧稳定性,采用OpenFOAM对分层旋流燃烧器的冷态和燃烧流场进行了大涡模拟.研究了旋流数对分层旋流流场结构和非稳态特性的影响.采用Q准则显示了流场中的瞬时涡结构;利用功率谱分析了流场中的进动特征.结果表明:在冷态工况下,旋流对回流区的位置和大小影响较小.随着旋流数增大,出口气流受到旋流...     

贫油预混预蒸发燃烧室回火现象的大涡模拟

为研究贫油预混预蒸发(LPP)燃烧室的回火现象及形成机理,基于大涡模拟(LES)及部分预混燃烧模型,对Schneider试验进行数值计算,分别将冷态及燃烧温度场与试验结果进行对比,验证了模型的准确性。将该模型应用于LPP燃烧室,并利用详细的化学机理对航空煤油火焰传播速度进行修正,成功重现了回火现象的发生并向主燃级上游传递的动态过程,通过流场分析发现局部形成的不规则涡结构在回火现象中扮演着重要角色。进一步从涡量方程的角度对涡结构进行分析,发现燃烧引起的斜压及体积膨胀是诱发回火现象的主要原因。     

方孔横向射流的大涡模拟研究

采用大涡模拟方法,对射流雷诺数ReD=4700、吹风比R=1.5的方孔横向射流进行了计算,为处理完全三维的复杂流动,亚格子模型使用了局部平均的动力模型.通过对时间平均流场的分析,了解了射流的穿透、扩散性质以及肾形涡的形成原理;而瞬态流场的时间演化序列则揭示了大尺度相干结构在流场演化中的主导作用.计算结果表明:使用局部平均动力模型的大涡模拟方法健壮性好,具有较高的数值精度,能够准确给出方孔横向射流的复杂流动图画,适合进行这类复杂流动问题的机理研究.     

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

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

超声速混合层液滴两相流动的大涡模拟

为研究超声速气流中液滴与气流的混合及液滴蒸发对混合的影响,采用大涡模拟(LES)方法数值仿真了超声速混合层内液滴两相流流场结构,气相流场采用亚格子(SGS)模型和切应力输运(k-ωSST)湍流模型,液相模拟采用轨道模型和单液滴蒸发模型。在混合层前缘入口处均匀持续地投放液滴,并在液滴入口处下方添加非周期小扰动,并观察液滴蒸发过程对该小扰动造成的影响。分析了入口小扰动在流场中不同的发展情况,发现液滴的蒸发过程使混合层厚度增加并加速混合层的发展,对气相流场扰动较强,可能导致流动失稳,对混合过程有很大的促进作用。     

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

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

基于湍动能谱的自适应大涡模拟方法

为了在计算量有限的前提下实现对不同尺度上湍流燃烧结构的直接求解,提出了新的基于湍动能谱的自适应大涡模拟的概念。在这一概念下,利用湍动能谱和自适应网格技术建立了跨尺度算法,从而实现DNS和LES两种方法的耦合。在湍流非预混氢气射流火焰的模拟中,该方法既能捕捉大尺度湍流结构,又能分辨小尺度燃烧结构,计算出的物理量的瞬态和时均分布和DNS结果吻合得较好。计算偏差主要体现在DNS和LES区域的过渡处,这与亚格子湍动能在过渡处的变化不够平滑有关,表明跨尺度算法需要优化,以完善对不同尺度间湍动能传输的模拟。      

合成气旋流预混燃烧的大涡模拟

如何高效洁净地燃烧合成气是整体煤气化联合循环(IGcc)系统中面临的主要问题.利用旋流贫燃预混燃烧可以达到强化传热传质,提高燃烧反应效率,降低NOx等目的,但由于合成气中所含有的H2成分,同样面临燃烧不稳定问题.为了对合成气旋流预混燃烧室的燃烧不稳定性及污染物排放提供理论和数值上的指导,采用大涡模拟方法,对典型的旋流预混燃烧室中的合成气燃烧进行了模拟.结果显示,随着H2含量的增加,火焰缩短向上游发展,燃烧室下游涡破碎更加集中.较高的入口雷诺数,加剧了涡团破碎的程度和速度,同时影响中心回流区和几何回流区的大小及分布.随着当量比的减小,火焰褶皱程度加大,燃烧趋向不稳定.增大压力火焰长度会明显缩短,NOx排放量也会随之增大,但NOx的反应速率只有在较低压力下才会对压力变化比较敏感.     

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

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

基于跨声速高压涡轮凹槽叶顶间隙模型的大涡模拟研究

为研究跨声速高压涡轮叶顶间隙非定常流动特性及流场结构,基于跨声速高压涡轮凹槽叶顶间隙的几何特征,建立了可行的简化数值模型,并通过大涡模拟对叶顶凹槽间隙内部的非定常流动进行了数值计算,探讨了凹槽叶顶几何参数对流动稳定性和气动性能的影响.研究结果表明,当跨声速泄漏流动流经凹槽叶顶时,在凹槽入口处形成激波,同时凹槽前分离泡发...     

浅槽孔气膜冷却的大涡模拟研究

为了探究浅槽孔气膜冷却的强化冷却机理,针对平板浅槽型气膜孔冷却射流与主流相干作用进行大涡模拟研究,并与圆孔射流进行对比分析。研究结果表明:浅槽孔内部会形成一对循环涡旋,将冷却气向展向卷吸;浅槽孔下游会同时出现肾形涡对和反肾形涡对,反肾形涡对抑制了气膜冷却射流的抬升,改善了冷气贴壁性;与圆孔相比,浅槽孔气膜冷却射流的流场拟序结构分布更为无序,浅槽前缘未出现马蹄涡旋,射流下游也未形成串列的发卡涡;当吹风比增长到1.5时,浅槽孔在低频区(400Hz处)出现主频,在高频区的振幅则远低于圆孔,表明高吹风比下浅槽孔气膜冷却射流流场中低频小尺度涡旋占主导作用。整体而言,浅槽孔提升了气膜的展向扩展能力,冷却性能远优于传统圆孔。     

湍流大涡数值模拟的现状与展望

近年来,随着计算机技术和计算流体力学(CFD)技术的发展,大涡数值模拟(LES)越来越多地应用于解决湍流流动问题.本文介绍了当前LES研究与应用的主要领域以及取得的成果,分析了LES的局限及实施的难点,并对LES未来的研究方向和重点进行了总结.     

超音速底部流动的大涡模拟研究

超音速底部流动是CFD的难点问题之一。采用基于可压缩流动的大涡模拟方法求解N-S方程,对超音速底部流动进行了模拟研究,分析了其流动结构。计算结果表明：得到的速度型、底面压力分布与实验值吻合,优于雷诺平均方法的结果。