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

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

基于超大涡模拟的燃烧室气动性能仿真研究进展

航空发动机燃烧室涉及旋流、雾化蒸发、掺混、化学反应、湍流与火焰相互作用等多尺度强耦合物理化学过程,相关的高 精度建模和数值模拟面临极大的挑战。超大涡模拟是近些年发展的兼顾计算精度、计算效率和强鲁棒性的数值模拟新方法,具备 试验室尺度和复杂工程应用场景下湍流流动与燃烧仿真能力。针对航空发动机燃烧室相关流动与燃烧基本特征,阐述了超大涡 模拟的理论方法及特点,从旋流流动、湍流燃烧、液雾雾化、碳烟生成、燃烧不稳定等典型多物理过程,以及双旋流模型燃烧室和高 温升燃烧室气动性能集成仿真等方面介绍了超大涡模拟的研究进展,对涉及的物理机制进行了分析,为超大涡模拟在航空发动机 燃烧室中规模化工程应用提供了坚实支撑。超大涡模拟在较低的计算资源消耗下具备与传统大涡模拟相当的计算精度,是一种 经济可承受的燃烧室高精度气动性能仿真新方法。     

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

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

湍流贫燃预混射流火焰的大涡模拟

为了研究高Karlovitz(Ka)数下湍流预混火焰的驻定机制和火焰结构,采用化学热力学建表方法结合假定概率密度函数(PDF)模型,对值班预混射流火焰(PPJB)进行大涡模拟(LES)研究。选取自点火模型作为建表模型,使用假定PDF模型考虑湍流和化学反应之间的相互作用,其中假定双混合物分数的概率密度分布为Dirichlet分布,反应进度变量的概率密度分布为β分布。模拟结果表明,能较好预测PPJB火焰的温度和主要组分的分布。使用中间组分(CO和OH)的质量分数来定性描述PPJB火焰结构,结果表明高Ka数下火焰结构受到湍流涡旋的显著影响。     

基于蚁群算法的TSP的仿真与研究

对于旅行商问题的研究始于19世纪,从20世纪中叶开始随着计算机技术的发展而不断发展。蚁群算法是一种新型的优化算法,于20世纪90年代提出,最早成功应用于解决旅行商问题。研究表明,蚁群算法有着极强的鲁棒性发现较好解的能力。通过编程实现了用蚁群算法解决旅行商问题,通过仿真实验研究了各参数对算法的影响。     

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

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

高温升燃烧室综合燃烧性能超大涡模拟

为了对高温升燃烧室性能进行计算分析,运用超大涡模拟方法开展数值模拟,并同步开展了雷诺平均数值模拟作为对 比。计算结果表明：旋流器设计与火焰筒开孔设计的匹配合理,Rothstein提出的射流迹线公式能够合理预测主燃孔的射流穿透。超 大涡模拟计算得到的出口径向温度分布系数的剖面曲线趋势和最大值位置均与试验结果符合较好。通过与试验结果进行定量比较 发现,在相同计算网格条件下,超大涡模拟方法预测高温升燃烧室综合燃烧性能的精度明显高于雷诺平均方法的。     

燃烧室与涡轮气动性能全3维联合数值仿真

针对燃烧室与涡轮2大部件相互耦合关系下多物理场的精准辨析问题,开展了基于CFX统一软件平台的部件间气动性 能全3维联合数值仿真研究。选取地面起飞状态点作为数值仿真工况,采用化学热力学建表方法计算化学反应过程,通过源项法 模拟叶片冷却喷射。根据计算结果分析了2大部件相互耦合关系下多物理场的分布规律,结果表明：各流场参数分布趋势合理, 验证了部件间全3维联合数值仿真的可行性;与部件单独仿真结果进行了对比,结果表明：部件联合仿真与单独仿真获取的压力 和马赫数的演化规律基本一致,温度场分布表现出差异,联合仿真克服了部件单独仿真无法准确捕捉热斑的不足。部件联合仿真 为多部件耦合关系下气动性能匹配和优化设计提供了工程可用的3维分析手段。     

拉伸层流预混火焰化学热力学建表方法研究

通过求解物理空间上简化的层流甲烷/空气对撞预混火焰面方程,获得不同当量比下考虑拉伸效应的火焰面高维数据。为更准确地描述流动与燃烧的相互作用,采用基于矩阵分析的主成分分析法(PCA)对火焰面高维数据进行分析,并选用分析出的前两个主成分作为新的建表特征标量,进而构建考虑拉伸效应的层流预混火焰化学热力学表。相比于已有的拉伸层流预混火焰面建表方法,使用PCA方法得到的两个建表特征标量之间不相关,无需条件概率密度信息。同时PCA方法具有一定的通用性,可以推广到不同燃料不同当量比下的拉伸层流预混火焰面建表中。     

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

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