国家自然科学基金(51176178)
- 作品数:11 被引量:48H指数:5
- 相关作者:叶桃红韩超张宏达张济民张培更多>>
- 相关机构:中国科学技术大学更多>>
- 发文基金:国家自然科学基金更多>>
- 相关领域:动力工程及工程热物理航空宇航科学技术理学更多>>
- 稀释对射流扩散火焰抬举及吹熄特性影响的实验研究被引量:3
- 2018年
- 以低热值气体燃料的燃烧利用为背景,以N_2或CO_2稀释的甲烷或丙烷及内径2mm或3mm的喷管作为多重变量,实验研究了射流扩散火焰的抬举速度及吹熄速度随燃料稀释比的变化规律.结果发现,随着稀释比的增大,火焰抬举速度近似呈线性缓慢减小,吹熄速度近似呈指数快速衰减.分析表明,这主要归因于稀释导致了当量混合层流火焰速度、初始燃料质量分数及燃烧放热的减小,并对低热值气体燃料的燃烧器设计提出了建议.此外,大于临界稀释比时,附着火焰随射流速度的增大不能转变为稳定的抬举火焰而会发生直接吹熄,依据甲烷和丙烷临界稀释比对应的射流雷诺数推测出二者的临界稀释比可能受不同机理的控制.
- 严野蒋年鹏张海涛
- 关键词:低热值气体燃料
- 钝体绕流有旋流中回流区与进动涡核的大涡模拟被引量:12
- 2014年
- 针对旋流数为0.57、0.68、0.91和1.59四种工况下的悉尼旋流燃烧器的冷态流场进行了大涡模拟,选取动态Smagorinsky涡黏模型作为亚格子尺度的湍流模型,研究不同旋流数下的流场结构、进动频率和进动涡核。模拟结果表明,旋流数为0.68时,钝体回流区长度最短。随着旋流数的增加,中心射流出口的旋流剪切层不断衰减,而下游的旋流剪切层不断增强。功率谱分析表明,进动现象的出现和消失对应于与旋流剪切层的增强和衰减;中心射流与下游区域具有不同的进动频率,表明流场中存在着两个独立的大尺度涡旋结构。不同取值位置的周向速度相关性分析进一步佐证了两个涡旋结构的存在。轴向位置70 mm处的横截面上瞬时流线和压强分布证实了下游流场存在着进动涡核。瞬时压强等值面显示了中心射流出口和下游流场进动涡核的三维螺旋形结构。下游流场的进动涡核均与平均速度场流线在空间上成正交关系,表明进动涡核是由剪切层Kelvin-Helmholtz不稳定性产生。
- 张济民韩超张宏达叶桃红
- 关键词:大涡模拟钝体
- 湍流贫燃预混射流火焰的大涡模拟被引量:7
- 2015年
- 为了研究高Karlovitz(Ka)数下湍流预混火焰的驻定机制和火焰结构,采用化学热力学建表方法结合假定概率密度函数(PDF)模型,对值班预混射流火焰(PPJB)进行大涡模拟(LES)研究。选取自点火模型作为建表模型,使用假定PDF模型考虑湍流和化学反应之间的相互作用,其中假定双混合物分数的概率密度分布为Dirichlet分布,反应进度变量的概率密度分布为β分布。模拟结果表明,能较好预测PPJB火焰的温度和主要组分的分布。使用中间组分(CO和OH)的质量分数来定性描述PPJB火焰结构,结果表明高Ka数下火焰结构受到湍流涡旋的显著影响。
- 张宏达叶桃红陈靖赵马杰
- 关键词:大涡模拟
- 高雷诺数湍流横侧射流的大涡模拟被引量:3
- 2016年
- 采用大涡模拟(LES)方法研究横向主流与壁面射流均为高雷诺数的壁面横侧射流(JICF),采用动态Smagorinsky涡粘模型对亚格子尺度进行封闭。数值模拟结果中,近场和远场的平均速度、均方根值均与实验结果符合很好。瞬时压力等值面在射流迎风侧与射流迹线成正交关系,表明JICF近场迎风涡是由于射流出口上游剪切层Kelvin Helmholtz不稳定性引起的。基于长度尺度RD,Broadwell和Breidenthal提出了射流迹线公式,通过拟合数值模拟结果得到公式中常数A=1.35,B=0.3。采用平均流线、平均场涡量与Q准则分析三维涡旋结构。横向主流遇到射流后发生夹带现象,在射流出口下游受涡旋结构的影响发生卷起;平均涡量等值面表明,反向旋转涡对(CVP)的涡旋方向与射流出口下游马蹄涡两个分支的涡旋方向相反;由于横向主流和射流雷诺数较高,Q准则表征的涡旋结构在射流出口下游一倍射流直径位置开始发生破碎。
- 赵马杰曹长敏张宏达叶桃红
- 关键词:湍流大涡模拟马蹄涡
- 甲烷/空气射流抬举火焰的大涡模拟计算被引量:14
- 2014年
- 为了研究中高温伴流中的自动着火特性,采用大涡模拟(LES)和均匀搅拌器(PSR)结合假定概率密度函数的建表湍流燃烧模型对甲烷/空气射流在高温伴流中的抬举火焰数值研究,并比较反应进度的不同概率密度分布对结果的影响。结果表明,计算结果与实验值符合得较好。PSR建表方法能够捕捉到高温伴流抬举火焰的点火过程,能够正确预测火焰的抬举高度,但在火焰稳定燃烧区域该建库燃烧模型还存在不足。同时,反应进度的亚格子脉动对准确模拟抬举高度十分重要。
- 韩超张培叶桃红陈义良
- 关键词:大涡模拟
- 超声速H_2/Air湍流扩散燃烧RANS数值模拟被引量:2
- 2015年
- 为了研究修正的火焰面反应进度变量燃烧模型在超声速湍流扩散燃烧问题中的适用性,对德国宇航中心(DLR)超声速燃烧室开展RANS数值模拟。基于Open Foam软件平台中密度求解器分别对三维冷态场和燃烧场进行模拟分析。将网格自适应加密技术用于流场的计算;燃烧场计算中,通过分析不同压力下层流火焰面数据库,引入了反应进度变量源项的压力修正系数,压力修正系数α等于2.2。计算结果表明,冷态场中压力分布、波系分布、速度分布以及燃烧场中波系分布、速度分布、温度分布结果均与实验值符合较好。压力修正方法能够较好地解决超声速湍流扩散燃烧问题。湍流Schmidt数敏感性分析表明,湍流Schmidt数Sct对湍流火焰结构有较大影响,文中Sct等于0.7时能得到与实验值较为一致的分布。
- 曹长敏叶桃红
- 关键词:数值模拟
- 大涡模拟研究钝体有旋流流场的拟序结构被引量:17
- 2014年
- 对带有中心钝体的分层旋流燃烧器的冷态流场进行了大涡模拟(LES),选取动态Smagorinsky涡黏模型作为亚格子模型,研究旋流数为0.45时旋流场的大尺度拟序结构。模拟结果表明:瞬时压强等值面显示的内外螺旋涡及进动涡核(PVC)均与平均速度场流线在空间上呈正交关系,表明两种涡结构均由剪切层Kelvin-Helmholtz不稳定性产生。Q准则等值面显示内外螺旋涡在下游20mm左右开始发生破碎。PVC发源于环形旋流与环形射流剪切层附近的下游区域。瞬时周向速度的功率谱密度(PSD)出现明显的特征峰,表明PVC影响附近的流体,使之出现进动特征。采用本征正交分解(POD)重构湍流脉动速度场,不同模态下功率谱密度结果表明前两个模态的大尺度结构具有进动特征,使用前两个模态的周向脉动速度等值面显示了PVC周围流体的大尺度涡旋结构。
- 张宏达张济民韩超叶桃红
- 关键词:本征正交分解
- 超声速横侧射流混合特性的大涡模拟
- 2015年
- 采用大涡模拟方法研究了Gamba超声速燃烧室内的横侧射流流场中大尺度涡旋结构以及当地混合特性.超声速来流受到音速射流流体的阻碍,形成了复杂的激波和涡旋结构.由计算结果中的平均马赫数分布图可以清楚地看到激波结构,包括弧形激波、λ激波、桶状激波以及马赫盘;采用Q准则表征三维涡旋结构,可以看到稳态的反向旋转涡对(CVP)、尾迹反向旋转涡对(TCVP)、马蹄涡以及非稳态的射流剪切层涡等结构;此外,由平均流线图可以看到,TCVP结构与CVP结构的旋转方向相反,不对称的CVP结构导致燃料质量分数展向分布不均匀.引入概率密度函数方法分析当地混合特性,结果表明射流近场混合主要发生在射流出口上游回流区以及桶状激波下侧和射流剪切层下侧的射流尾迹区内;射流远场混合分数的概率密度分布从β分布逐渐过渡为高斯分布.研究射流浓度衰减特征,结果表明氢气质量分数沿射流浓度最大迹线呈指数-0.7衰减.
- 曹长敏赵马杰周涛涛叶桃红
- 关键词:大涡模拟概率密度函数
- 自由旋流冷态流场中涡旋破碎泡的大涡模拟被引量:2
- 2013年
- 采用大涡模拟方法(LES)模拟了悉尼旋流燃烧器的中等旋流数算例(N29S054和N29S045)的冷态流场,研究了涡旋破碎泡(下游二次回流区)的不稳定模式.LES结果得到的统计矩总体上与实验值符合得较好.研究发现:涡旋破碎泡没有螺旋形结构,其周围也没有出现螺旋形状的进动涡核,表明涡旋破碎泡没有典型的进动特征;瞬时速度分布显示了涡旋破碎泡存在着周期性的收缩/崩塌与膨胀的现象;功率谱的特征峰证实了涡旋破碎泡存在着周期性运动.
- 张济民韩超叶桃红
- 关键词:大涡模拟
- 不同湍流预混燃烧模型在本生灯火焰中的比较被引量:4
- 2014年
- 湍流预混燃烧的大涡模拟数值计算近来受到很大的关注,为了进一步发展湍流预混燃烧模型,利用正在发展的湍流预混燃烧模型等值面集合的G函数方法、层流预混火焰生成流型的化学热力学建表FGM方法、G函数/FGM耦合方法三种燃烧模型,对甲烷/空气湍流本生灯火焰进行了大涡模拟计算,并对上述燃烧模型进行了比较。大涡模拟结果表明,G函数方法能够较为准确捕捉火焰的传播特性,FGM方法能够捕捉到火焰和湍流的相互作用,但无法准确预测火焰的传播特性。而耦合的方法既可以捕捉火焰的传播特性,又能够捕捉到火焰和湍流的相互作用。
- 韩超张培叶桃红陈义良
- 关键词:湍流预混火焰G函数