黄渊
- 作品数：12 被引量：36H指数：4
- 供职机构：中国空气动力研究与发展中心更多>>
- 发文基金：国家自然科学基金更多>>
- 相关领域：航空宇航科学技术更多>>
- 轴向三级旋流燃烧室流场结构大涡模拟被引量：2
- 2017年
- 为深入了解真实航空发动机燃烧室内部复杂流场结构,在自有CFD平台上采用动态亚网格模型对一种轴向3级旋流燃烧室的单个头部矩形试验模型0.5MPa下冷态流场结构进行了LES(大涡模拟).为避免试验模型简化误差,对包括火焰筒上约2 000个气膜孔在内的燃烧室所有精细结构进行了完全仿真.计算模拟了燃烧室内复杂流场从静止启动到统计定常状态的完整非定常发展过程,成功捕捉到主旋流与横向对冲射流相互影响作用及涡旋破碎等细观结构,获得的测点湍动能谱与湍流经典理论中惯性子区-5/3规律一致,LES时间平均流场结构与已有PIV(particle image velocimetry)试验结果吻合,表明所建立的高仿真网格与LES方法可进一步用于真实航空发动机环形燃烧室流场数值模拟.
- 周瑜乐嘉陵陈柳君黄渊
- 关键词：航空发动机大涡模拟
- 采用CARS试验技术与UFPV数值方法研究航空发动机燃烧室（英文）
- 2017年
- 在自主开发的软件平台上,采用基于URANS的方法计算航空发动机燃烧室的三维两相燃烧流动,考虑了液态燃油从液膜-液滴-燃气-燃烧的完整物理化学过程。其中,颗粒相采用LISA一次破碎模型,KH-RT二次破碎模型和标准的蒸发模型,湍流燃烧模型采用可以考虑非稳态燃烧特性的非稳态火焰面/反应进度变量方法,得到了航空发动机燃烧室中温度、组分浓度和燃油液滴的颗粒直径分布规律。同时,采用CARS光学手段测量燃烧室主燃区的温度分布,并将数值计算结果与光学试验测量值进行比较,数值计算结果和试验值吻合较好,数值计算误差小于7.3%。说明了本文的数值计算方法和UFPV方法在计算航空发动机燃烧室的两相燃烧流动时具有较高的精度。
- 熊模友乐嘉陵黄渊宋文艳杨顺华郑忠华
- 关键词：航空发动机燃烧室两相燃烧
- 采用火焰面/反应进度变量方法模拟湍流燃烧被引量：3
- 2016年
- 在自主开发的软件平台上添加了火焰面/反应进度变量方法,选取了两种反应进度变量分别对甲烷/空气的同轴射流燃烧室进行了数值模拟.将计算结果与稳态火焰面模型、文献参考值以及实验值进行了对比研究,结果表明:采用火焰面/反应进度变量方法能捕捉到稳态火焰面模型所不能预测到的火焰抬举等非稳态现象,得到的结果与文献参考值以及实验值出入很小,能准确反应真实的燃烧过程;另外不同反应进度变量的定义对计算结果有较大的影响,采用火焰面/反应进度变量方法对主要中间产物一氧化碳的预测比文献参考值更好,更接近于实验值.
- 熊模友乐嘉陵黄渊宋文艳杨顺华郑忠华
- 关键词：湍流燃烧
- 直流喷射首次破碎的形变过程研究被引量：4
- 2016年
- 为深入了解直流喷射过程中射流柱的细部结构以及其脱落过程,利用LES结合VOF的方法,对静止大气中垂直射流现象进行模拟。模拟得到的液相喷雾结构和试验结果能够很好地吻合。通过计算观察得到射流柱进入到大气中后由于Rayleigh-Taylor(RT)不稳定性迅速形成伞状的头部,头部的边缘在不断变薄失稳的过程中脱落形成液带,液带脱落产生的液滴具有等距性,间距为0.106mm,并在伞状边缘的下方形成一个气涡,气涡与上游的射流柱相互作用,促进射流柱表面的脱落。射流柱表面由于Kelvin-Helmholtz(K-H)表面波的作用呈现鱼鳞状的结构,并导致射流柱整体断裂、破碎,其表面波波长由初始的0.26mm迅速增长到0.78mm。
- 刘日超乐嘉陵杨顺华郑忠华宋文艳黄渊
- 关键词：VOF雾化
- KH-RT模型在横向来流作用下射流雾化过程的应用被引量：8
- 2017年
- 针对射流柱在横向来流下的破碎过程,采用欧拉-拉格朗日模型对其进行了数值模拟。其中射流柱的一次雾化采用WAVE模型,颗粒的二次雾化采用KH-RT混合模型。计算得到的液滴颗粒平均直径在70μm左右,还得到了雾化场的喷雾结构、速度分布、索特尔直径分布,以及射流柱在不同来流速度下破碎过程等。同时通过采用不同二次破碎模型,分别将计算结果与试验结果进行对比,结果表明采用的KH-RT混合雾化模型在一定程度上能够更加真实地反映雾化过程,并较好地捕捉射流雾化特点。
- 刘日超乐嘉陵杨顺华郑忠华宋文艳黄渊
- 关键词：混合模型射流雾化
- 亚声速横向气流中液体射流破碎过程的直接模拟被引量：9
- 2016年
- 为研究亚声速横向气流中液体射流柱的变形、弯曲以及其破碎过程,利用LES结合VOF的方法,对射流破碎过程进行了直接模拟。通过计算观察得到射流柱进入到横流气体中后由于RayleighTaylor(RT)不稳定性和Kelvin-Helmholtz(K-H)不稳定性的共同作用迅速发生变形,并形成表面波,同时伴随着细小液滴的脱落,气流在射流柱后方形成涡状结构。在射流柱的迎风面上同时存在两种表面波,即RT表面波和KH表面波。模拟所得液柱纵向初始表面波波长为0.22mm,与K-H表面波波长理论公式计算所得相接近。通过对比射流柱的背风面和迎风面两侧,发现在横向来流中背风面上的大量细小液滴主要是由横向来流对液柱的不稳定性扰动造成迎风面上的液体向背风面方向挤压并剪切,最终脱落造成。
- 刘日超乐嘉陵杨顺华郑忠华宋文艳黄渊
- 关键词：雾化
- 三级轴向旋流燃烧室流场结构研究被引量：2
- 2018年
- 为深入了解燃烧室内流场,研究不同来流状态对燃烧室流场结构的影响,基于粒子成像速度仪(PIV)技术,对采用三级轴向旋流器的航空发动机燃烧室进行流场测量,分别在Case 1常温低压(0.49MPa)、Case 2常温中压(0.98MPa)、Case 3常温高压(1.64MPa)、Case 4全状态(高温813K高压2.78MPa)来流条件下进行。研究结果表明,同一燃烧室模型在不同速度、温度和压力来流下有基本相同的流场结构,但在中心回流区尺寸、角落回流区尺寸、主燃孔和掺混孔射流等细节方面仍有明显差异,来流压力较高的流场中心回流区向下游扩展更深入,角落回流区被压缩,主燃孔和掺混孔射流速度增加且进气比例增大。
- 陈柳君乐嘉陵张俊黄渊周瑜
- 关键词：航空发动机燃烧室
- 基于LES的非稳态火焰面/反应进度变量方法模拟部分预混抬举火焰被引量：2
- 2017年
- 为了描述部分预混燃烧的火焰机制,捕捉燃烧现象中的非稳态效应,采用大涡模拟与非稳态火焰面/反应进度变量方法相结合来模拟湍流部分预混燃烧抬举火焰,并将计算结果与实验测量值进行了对比研究。结果表明基于大涡模拟的非稳态火焰面/反应进度变量方法捕捉到了湍流部分预混燃烧中的火焰抬举现象,计算得到的火焰抬举高度大致为35,和实验测量值完全一致。同时采用各热力学参数在混合分数空间的分布来研究部分预混抬燃烧中具有非稳态燃烧特性的局部熄火现象。各不同截面上的燃烧热力学参数与实验测量值吻合较好,说明采用该方法能比较真实地反映湍流部分预混燃烧的火焰机制。
- 熊模友乐嘉陵黄渊宋文艳杨顺华郑忠华
- 关键词：大涡模拟非稳态
- 采用基于火焰面的燃烧模型研究部分预混燃烧被引量：3
- 2017年
- 考虑到火焰面模型的优点,采用基于RANS的稳态火焰面模型,稳态火焰面/反应进度变量模型和非稳态火焰面/反应进度变量模型对部分预混燃烧室进行了数值计算,并将这三种燃烧模型的计算结果和实验值进行对比研究。发现稳态和非稳态火焰面/反应进度变量模型均成功地预测到了部分预混燃烧中的三岔火焰结构和火焰抬举现象,分析了三岔火焰的形成机理及抬举高度。文中稳态火焰面模型计算部分预混燃烧完全失效,无法捕捉到火焰抬举现象,稳态火焰面/反应进度变量方法计算的火焰抬举高度仅仅为10,基于反应进度变量为水、二氧化碳、一氧化碳和氢气的质量分数之和的非稳态火焰面/反应进度变量方法计算的火焰抬举高度仅为20,和实验值之间的误差为42.8%,而基于反应进度变量为一氧化碳和二氧化碳质量分数之和的非稳态火焰面/反应进度变量方法计算的火焰抬举高度大致32,和实验测量值35非常接近,误差仅为8%。对燃烧热力学标量温度和组分的计算,可以发现非稳态火焰面/反应进度变量方法和实验结果吻合最好,其次是稳态火焰面/反应进度变量模型,最差的是稳态火焰面模型。
- 熊模友乐嘉陵黄渊宋文艳杨顺华郑忠华
- 关键词：部分预混燃烧模型
- 同轴射流燃烧室非预混湍流燃烧流场特性大涡模拟研究被引量：1
- 2017年
- 为深入研究带突扩与回流的同轴剪切射流模型燃烧室流场结构与燃烧特性,采用单步快速化学反应、火焰面与反应进度变量三种燃烧模型结合动态亚网格模型在自有CFD平台上对甲烷-空气非预混燃烧流场进行了大涡模拟。计算检验了流场时空尺度及网格和时间步长设置,对流场0~2s内启动并发展进入统计定常状态的非定常流动过程进行了完整模拟。燃烧室内多处监测点湍动能谱统计结果表明,轴向速度在混合剪切层内的脉动存在从约1200Hz开始的多个特征频率序列,且随着测点位置向下游移动,其对应的特征峰从高频到低频逐渐衰减消失。流向多个截面上时间平均的混合分数及温度分布结果表明,反应进度变量模型能够模拟出该流场特有的抬举火焰,与实验数据吻合程度显著优于其余两种模型,即对于该类型强剪切非预混反应流,以燃烧产物质量分数替代采用梯度粘性输运模化的混合分数耗散率进行计算更为合适。
- 周瑜乐嘉陵黄渊
- 关键词：非预混燃烧
