马杰, 王晓冰, 王红丽, 牛青林, 董士奎. 发动机尾喷焰复燃化学反应模型评价与重构[J]. 红外与激光工程, 2024, 53(7): 20240112. DOI: 10.3788/IRLA20240112
引用本文: 马杰, 王晓冰, 王红丽, 牛青林, 董士奎. 发动机尾喷焰复燃化学反应模型评价与重构[J]. 红外与激光工程, 2024, 53(7): 20240112. DOI: 10.3788/IRLA20240112
MA Jie, WANG Xiaobing, WANG Hongli, NIU Qinglin, DONG Shikui. Evaluation and reconstruction of afterburning reaction kinetics of rocket exhaust plume[J]. Infrared and Laser Engineering, 2024, 53(7): 20240112. DOI: 10.3788/IRLA20240112
Citation: MA Jie, WANG Xiaobing, WANG Hongli, NIU Qinglin, DONG Shikui. Evaluation and reconstruction of afterburning reaction kinetics of rocket exhaust plume[J]. Infrared and Laser Engineering, 2024, 53(7): 20240112. DOI: 10.3788/IRLA20240112

发动机尾喷焰复燃化学反应模型评价与重构

Evaluation and reconstruction of afterburning reaction kinetics of rocket exhaust plume

  • 摘要: 复燃效应的准确预估对于精细描述尾喷焰反应流场参数和提高尾喷焰红外辐射计算精度至关重要。文中以固体火箭发动机为研究对象,建立尾喷焰复燃有限速率化学反应模型,结合流体计算动力学(Computational Fluid Dynamics,CFD)方法和尾喷焰红外辐射计算模型,评估不同化学反应动力模型在尾喷焰流场参数和红外光谱辐射计算方面的精度,基于各化学反应速率曲线与试验数据重构适用于尾喷焰CO/H2反应体系的10步气相化学反应动力模型,并验证和校核复燃化学反应模型可靠性。结果表明:不同化学反应模型计算所得的尾喷焰流场结构差异微弱,轴向温度峰值最高相差200 K左右,差异主要发生在复燃区域;化学反应动力模型对不稳定产物CO影响最为显著,CO2分布差异主要发生在高含量区域,最大差异达到近50%,且低含量组分的差异高达两三个量级;在2.7 μm和4.3 μm典型波段内,不同化学反应工况下的尾喷焰光谱辐射峰值强度差异达到近40%;基于反应速率试验数据构建的9组分10步反应的CO/H2反应体系的尾喷焰辐射计算值与BEM-II试验数据的差异低于6%。该研究可为准确预测火箭发动机尾喷焰反应流场的红外辐射特性提供高保真化学反应动力模型。

     

    Abstract:
    Objective Diatomic/polyatomic molecules in rocket exhaust plumes emit specific bands of infrared radiation during high-temperature vibrational transitions, making them crucial radiation sources of concern in the measurement field. Usually, afterburning occurs when rocket exhaust plumes mix with air, releasing a large amount of heat and significantly raising the temperature level and infrared radiation of the plumes. Therefore, afterburning is a crucial step in accurately calculating the reacting flow field parameters and infrared radiation of rocket exhaust plumes. Using computational fluid dynamics (CFD) methods to predict the reaction flow field of rocket exhaust plumes and evaluate the degree of afterburning has become a feasible technical approach. This underscores the importance of constructing an accurate chemical reaction kinetics model to predict afterburning in rocket exhaust plumes. A highly accurate flow field structure is essential, as different chemical reaction dynamics models can lead to significant differences in the composition, content, and distribution of the flow field of rocket exhaust plumes. However, the infrared radiation characteristics of rocket exhaust plumes are extremely sensitive to the flow field temperature, component content and distribution. To improve the accuracy of infrared radiation calculation for rocket exhaust plumes, higher requirements are placed on the chemical reaction dynamics model for accurate flow field parameters of rocket exhaust plumes.
    Methods With solid rocket engines as the research focus, the central difference scheme method is employed to solve the three-dimensional Navier-Stokes (N-S) equations with chemical reaction sources. Based on the finite rate chemical reaction model expressed in the Arrhenius formula, a 10-step gas-phase chemical reaction kinetics model is developed for the CO/H2 reaction system. The gas radiation physical properties parameters are computed using the statistical narrow spectral band model, and the radiation transport equation is solved using the apparent light method. Through fitting of experimental data using the three-parameter Arrhenius formula, a chemical reaction kinetics model is structured to closely match the experimental data.
    Results and Discussions  There is little difference in the plume structure calculated by different chemistry models (Fig.4), and a maximum temperature difference of 200 K exists in the area where afterburning occurs (Fig.6). The influence of chemical reaction kinetics on the unstable product CO is most significant (Fig.9(b)). The maximum difference in CO2 reaches nearly 50% (Fig.9(a)), mainly occurring in high-content regions, and the impact on low mole fraction components varies by 2-3 orders of magnitude (Fig.10). The impact of different chemical reactions on the peak intensity of spectral radiation varies by nearly 40% in the 2.7 μm and 4.3 μm bands (Fig.11), and the difference in integrated spectral band intensity within different bands reaches about 40% (Fig.13). Based on the proposed CO/H2 reaction system with 9 components and 10 steps, the difference between the calculated spectral intensity of BEM-II plumes' infrared radiation and the measured data is less than 6% (Tab.4).
    Conclusions  Different chemical reaction kinetics models have a relatively small impact on the structure of the rocket exhaust plume flow field but have a significant effect on temperature, component generation, and infrared radiation characteristics. For different chemical reaction models, the collision frequency, temperature correlation index and activation energy parameters corresponding to each chemical reaction kinetic equation are different. Under the same operating conditions, the reaction rate and heat release (absorption) between components are different, which affects the temperature, component content, and distribution of the reaction flow field. A chemical reaction kinetics model with small errors was constructed by combining the trend of positive reaction rate curves corresponding to each reaction kinetics equation and experimental data. This study can provide high-fidelity chemical reaction models for accurately predicting the reaction flow field and infrared radiation of rocket exhaust plumes.

     

/

返回文章
返回