基于乙烯/空气详细反应机理的旋转爆轰过程数值模拟

为了研究乙烯/空气的旋转爆轰传播特性,采用OpenFOAM7中rhoReactingCentralFoam求解器,以及乙烯/空气21组分36步基元化学反应模型,对旋转爆轰传播过程进行了二维数值模拟研究。通过改变预混气各组分的质量分数,验证了基元反应在旋转爆轰的可适用性,研究了当量比对旋转爆轰波的影响,分析了旋转爆轰燃烧室的流场结构、压力增益和燃料比冲。研究结果表明:在当量比0.7～1.1范围内旋转爆轰波能够自持传播,其传播模态均为单波; 在来流总压0.6 MPa和总温300 K时,压力增益均保持在30%以上,且随当量比的增加旋转爆轰波传播速度、阻塞比、压力增益以及燃料比冲均增大,速度亏损逐渐减小,但随着当量比的增加其增长速率和降低速率呈下降趋势; 爆轰波稳定传播时,中间产物OH集中在爆轰波后方和斜激波下游,反映了化学反应阵面的形状。爆轰波传播过程中,由于爆轰波扫过后压力过高,新鲜燃料无法供应,导致温度曲线图出现波动。     

基于预爆轰点火方式的连续旋转爆轰发动机起爆过程分析

为深入了解预爆轰式连续旋转爆轰发动机点火起爆机理,基于连续旋转爆轰发动机试验平台进行了一系列预爆轰管点火试验和连续旋转爆轰发动机点火起爆试验研究,基于OpenFOAM开展了预爆轰管内初始爆轰波进入发动机环形燃烧室传播过程的数值模拟研究。结果表明:试验研究中,初始爆轰波进入环形燃烧室后波峰压力值迅速降低,爆轰波膨胀解耦并衰减为燃烧波; 初始爆轰波进入环形燃烧室后,并未直接形成旋转爆轰波,而是存在“初始爆轰波解耦—DDT—触发旋转爆轰波”过程; 在总质量流量为380 g/s左右时,随着当量比从0.65提高至0.95,DDT时间迅速从20 ms以上下降至2 ms以下; 当总质量流量下降至280 g/s附近时,出现了旋转爆轰波峰值压力大幅波动等旋转爆轰波不稳定传播现象。数值模拟中,初始爆轰波进入环形燃烧室后逐渐衰减,形成首道激波,该激波在燃烧室内壁面反射后形成反射激波,并伴随首道激波传播; 首道激波在传播过程中在燃烧室进口端面发生反射形成反射激波,该反射激波最终在周向约90°位置处衰减至消失。     

当量比对煤油-空气两相旋转爆轰波的影响

为研究当量比对煤油-空气两相旋转爆轰波的影响,基于欧拉-拉格朗日方法,以常温液态煤油为燃料、高总温空气为氧化剂开展二维数值模拟。通过控制进油量计算不同当量比下的两相旋转爆轰流场,分析当量比对流场结构、传播特性和推力性能的影响。分析结果表明：在一定的当量比范围内,爆轰波波头高度随着当量比的增大而减小;当量比大于0.6时,速度亏损随着当量比的增大总体呈现增加的趋势,变化范围为7.1%~17.2%;当量比大于0.8时产生燃料损失,燃料损失随当量比的增大呈近似线性增加,当量比为2.0时最大为51%;随着当量比的增大,比冲减小,比推力先增大后减小并在当量比为1.5时达到最大值,比冲和比推力的变化范围分别为1 154.3~ 3 912.4 s和1 226.8~1 521.8 N·s/kg。     

环形燃烧室中凹腔对C2H4/Air旋转爆轰流场影响的数值模拟

为研究环形燃烧室中凹腔对旋转爆轰流场的影响,通过在开源计算流体动力学软件OpenFOAM框架内求解Navier-Stokes方程,对C₂H₄/Air旋转爆轰燃烧室（RDC）进行三维数值模拟研究。保持进气总压为0.6 MPa,分别在总温300 K、600 K和800 K条件下对比同轴圆环形和凹腔基环形两种构型RDC的主要流场特征,研究凹腔对旋转爆轰波传播特性的影响,定量分析不同热释放速率下消耗的燃料占比,对比了不同热释放速率下消耗的燃料比例。结果表明：对于凹腔基环形RDC,在凹腔内存在回流区,导致在凹腔上游流速缓慢,但在凹腔收缩段流动明显加速,RDC出口的面平均马赫数大于与之对应的同轴圆环形RDC;受燃料向RDC出口和凹腔内壁方向侧向膨胀的影响,相对于同轴圆环形RDC,爆轰波在凹腔基环形RDC中传播时具有更高的速度亏损;部分新鲜燃料与燃烧产物在凹腔内混合,提高了爆轰波波前反应物的温度;在相同进气条件下,同轴圆环形RDC以爆轰形式消耗的燃料占比更多。所得研究结果对RDC的结构设计和优化具有一定的指导意义。