章康康， 丁晟， 陆倩， 王栋
安徽工业大学 建筑工程学院， 安徽 马鞍山 243002
引用格式：
章康康， 丁晟， 陆倩， 等. 固体烟火剂中镁颗粒非均相燃烧及红外辐射特性［J］. 中国粉体技术， 2027， 33（3）： 1-11.
Zhang Kangkang， Ding Sheng， Lu Qian， et al. Analysis of heterogeneous combustion and infrared radiation characteristics of magnesium particles in solid pyrotechnic compositions［J］. China Powder Science and Technology， 2027， 33（3）： 1-11.
DOI：10.13732/j.issn.1008-5548.2027.03.017
收稿日期： 2026-05-08， 修回日期： 2026-06-09， 上线日期： 2026-10-网上出版日期:26。
基金项目： 国家重点研发计划项目，编号：2025YFF0521003；国家科技计划项目雄安新区科技创新专项，编号：2025XAGG0034）；安徽工业大学青年教师科研基金项目，编号：QZ202416。
第一作者简介： 章康康（1994—），男，讲师，博士，硕士生导师，研究方向为金属颗粒燃烧理论与技术技术。E-mail：zkknjust@163.com。
摘要： 【目的】 研究镁基复合含能材料燃烧过程中燃料镁颗粒的红外辐射特性。【方法】 基于欧拉-拉格朗日法建立燃烧流场气固两相耦合特性求解模型，综合考虑燃烧流场中参与性介质发射、吸收和散射作用，采用反向蒙特卡洛法，基于C语言开发非均相燃烧火焰红外辐射传输特性仿真软件，通过燃烧光学诊断实验结果验证数值模型的合理性。【结果】 镁基复合含能材料在红外波长0.8~14.0 μm时呈现灰体辐射与选择性辐射叠加特征。动态扩散燃烧的镁颗粒在红外波长3.0~5.0 μm时贡献了21.32%的总辐射通量，而在红外波长8.0~14.0 μm时贡献率降至4.61%。【结论】 燃烧流场中多相辐射源的中红外波段由颗粒辐射主导，而远红外波段受高温气体辐射支配。
关键词：烟火剂；燃烧特性；红外辐射；数值模拟
Abstract
Objective Magnesium-polytetrafluoroethylene （Mg-PTFE） composite energetic material is an important energy carrier in aerospace and optoelectronic countermeasure applications. This material effectively improves the performance of thermal protection systems and base-bleed devices. Compared with most metal fuel mixtures， the Mg-PTFE composite generates significantly higher energy output and high-intensity spectral signals， making it an ideal payload for infrared decoy flares. However， the analysis of the combustion mechanism is challenging due to the complex multi-physics coupling involved. Consequently， the characterization of infrared radiation properties remains a weak link in numerical simulations. To address this issue， the gas-solid two-phase coupling mechanisms and infrared radiation characteristics are comprehensively analyzed based on computational fluid dynamics and radiation transport theories. The methods and results of this research provide robust theoretical support for the engineering design of related devices.
Methods In this study， experimental samples were first prepared by wet-mixing Mg-PTFE composite powder with dissolved Viton at a specified mass ratio， followed by pressing the mixture into a mold. Second， an optical diagnostic system was established to capture spectral data from the visible to the far-infrared bands. Third， a gas-solid two-phase coupling model was developed based on the Euler-Lagrange approach， wherein the high-temperature gas was treated as the continuous phase and the burning particles as the discrete phase. Source terms for mass， momentum， energy， and species were incorporated into the governing equations. A chemical kinetic mechanism comprising 15 species and 13 elementary reactions was employed to describe the Mg-PTFE combustion system， utilizing Arrhenius-type rate expressions to determine the generation rates of condensed-phase products. Fourth， the influence of participating media on radiative heat transfer was quantitatively described by the radiative transfer equation. Based on the probability density function sampling principle and the reverse Monte Carlo method， a custom simulation code was developed to evaluate the infrared radiation transmission characteristics of the heterogeneous combustion flame. Finally， the multiple radiation sources model was validated by comparing the numerical predictions with the experimental measurements.
Results and Discussion Based on the experimental tests， the combustion process was observed to comprise three stages： deflagration， stable combustion， and flameout. The stable combustion stage was characterized by high average infrared radiance and the formation of a stable effective radiation area. Through spectral data analysis， a continuous spectrum feature was observed， mainly originating from the radiation superposition of high-temperature gas molecules. The radiation peak， calculated by Wien’s displacement law， was located in the wavelength range of 1.2~1.4 µm. The strong selective radiation in the 4.2~4.6 µm wavelength range was attributed to the fundamental frequency vibration transitions of CO and CO2. Approximately 87% of the total radiation energy was distributed within the 1.0~5.0 µm wavelength band， while less than 5% was distributed in the 8.0~14.0 µm far-infrared wavelength band. The numerical simulation results showed high consistency with the experimental data. The simulated axial diffusion height of the high-temperature flame zone agreed well with the experimental measurements. Furthermore， the spatial distribution characteristics of the infrared radiance were evaluated. In the 3.0~5.0 µm wavelength band， 64.2% of the total radiation flux was contributed by the region near the solid burning surface， which was dominated by typical gray-body radiation. The high-temperature gas products contributed 13.65%， whereas the discrete burning particles contributed 21.32%. However， in the 8.0~14.0 µm wavelength band， the contribution rate of the burning particles dropped sharply to 4.61%. This phenomenon was caused by the high temperatures of the discrete burning particles， resulting in the dominance of short-wave radiation.
Conclusion In this study， the infrared radiation energy output characteristics of Mg-PTFE composites are investigated. Based on the experimental and numerical results， the Mg-PTFE composite exhibits a superposition of gray-body and selective radiation in the 0.8~14.0 µm wavelength band， with a peak located near 1.3 µm. Notably， this study innovatively achieves the quantitative decoupling characterization of multiphase radiation sources in the combustion flow field. Dynamically diffusing Mg particles are found to contribute 21.32% of the total radiation flux in the 3.0~5.0 µm wavelength band， whereas this contribution drops sharply to 4.61% in the 8.0~14.0 µm wavelength band. This explicitly demonstrates that the mid-infrared band is dominated by particle radiation， while the long-wave infrared band is governed by high-temperature gas radiation. The proposed approach is a highly effective numerical methodology for analyzing multiphase energetic materials. Furthermore， the numerical model is proven to be accurate. The validated approach can serve as a solid theoretical basis for the engineering design of related devices， such as infrared decoy flares.
Keywords： pyrotechnic compositions； combustion characteristics； infrared radiation； numerical simulation
参考文献（References）
［1］Meiirbekov M， Sadykov M， Kuandyk A， et al. Modern thermal protection systems for solid rocket motors： a comparative review of composite manufacturing technologies［J］. Results in Engineering， 2025， 27： 106550.
［2］Ma Longze， Yu Yonggang. Coupling characteristics of combustion-gas flows generated by two energetic materials in base bleed unit under rapid depressurization［J］. Applied Thermal Engineering， 2019， 148： 502-515.
［3］Elbasuney S， Elsaidy A， Kassem M， et al. Infrared spectra of customized magnesium/Teflon/Viton decoy flares［J］. Combustion， Explosion， and Shock Waves， 2019， 55（5）： 599-605.
［4］Debnath S， Kumar H， Kulkarni P， et al. Investigation of the combustion and radiometric behaviour of magnesium-Teflon-Viton（MTV） decoy flares by partial replacement of magnesium［J］. Central European Journal of Energetic Materials， 2025， 22（2）： 141-165.
［5］Orzechowski A， Nita M， Powała D， et al. Investigation of the ignition of MTV decoy flares［J］. Central European Journal of Energetic Materials， 2021， 18（4）： 529-544.
［6］Gaunekar S， Ambekar A. Effect of hygrothermal aging on the decomposition reactions of magnesium/Teflon/Viton （MTV） pyr-otechnic［J］. Propellants， Explosives， Pyrotechnics， 2025， 50（2）： e202400162.
［7］Lee Y， Oh J， Yoh J J. Effect of hygrothermal aging on the pyrolysis of magnesium and its mixture［J］. Fuel， 2025， 387： 134424.
［8］Ma Shuai， Huang Hongyu， Zou Haoming， et al. Advanced moisture-resistant Mg@PDA@PTFE： proportional optimization and energy release characterization［J］. Fuel， 2025， 381： 133374.
［9］Zhu Chenguang， Xu Chungen， Xue Rui. Study of the spatial distribution of burning particles in a pyrotechnic flame based on particle velocity［J］. Journal of Energetic Materials， 2014， 32（4）： 252-263.
［10］Zhang Kangkang， Lu Qian， Zhou Rongrong， et al. Optical experimental diagnosis and numerical study on the combustion flow field of Mg/Teflon/Viton pyrotechnics［J］. Fuel， 2026， 407： 137434.
［11］Williams M R， Dunne J L. Heat flux measurement from bulk magnesium/Teflon/Viton （MTV） flare composition combustion［J］. Propellants， Explosives， Pyrotechnics， 2013， 38（4）： 563-568.
［12］Zhang Kangkang， Han Yuge， Ren Dengfeng， et al. Experimental and numerical investigations on the condensed combustion products of Mg-based solid fuels［J］. Fuel， 2022， 329： 125423.
［13］林长津， 王浩， 朱晨光. 基于流动与化学反应耦合的镁/聚四氟乙烯烟火剂二维燃烧模型及数值计算［J］. 兵工学报， 2016， 37（2）： 287-292.
Lin Changjin， Wang Hao， Zhu Chenguang. 2D combustion model and numerical simulation of Mg-Teflon pyrolant based on flow-chemistry coupling［J］. Acta Armamentarii， 2016， 37（2）： 287-292.
［14］Xue Xiaochun， Yu Yonggang. Investigation of thermal performance in dual-fuel units［J］. Energy Science & Engineering， 2022， 10（12）： 4782-4803.
［15］Xue Xiaochun， Ma Longze， Yu Yonggang. Effects of Mg/PTFE pyrotechnic compositions on reignition characteristics of base bleed propellants and heating mechanism［J］. Defence Technology， 2022， 18（1）： 94-108.
［16］Zhan Junbiao， Yu Wenhao， Liu Mengying， et al. Numerical and experimental study on combustion instability of a dual-burning-rate solid rocket motor［J］. Acta Astronautica， 2026， 245： 937-953.
［17］Shao Shiyao， Zhang Wenchao， Yue Songchen， et al. Effect of pressure， Mg content and AP size on the microscale flame structure of Mg-based propellant［J］. Aerospace Science and Technology， 2024， 150： 109213.
［18］Liu Yichao， Zhou Bohuai， Shi Chengkuan， et al. Study and numerical analysis of the combustion behavior of Mg-Fe2O3 thermite in sub-atmospheric environment［J］. Fuel， 2026， 407： 137361.
［19］Liu Yichao， Zhou Bohuai， Shi Chengkuan， et al. Numerical analysis of the impact of charge size on the combustion mechanism and radiation-enhanced effect in Mg/PTFE/Viton pyrolants［J］. Infrared Physics & Technology， 2025， 151： 106110.
［20］Guo Haoyang， Yin Mengli， Shi Shengnan， et al. CFD simulation and experimental study on the flow field and typical energetic dust dispersion in the 5 L column vessel［J］. Journal of Loss Prevention in the Process Industries， 2024， 92： 105473.
［21］何金龙， 杨阳， 刘晓星. 多点加载下单颗粒破碎特性的离散单元法模拟［J］. 中国粉体技术， 2025， 31（3）： 89-100.
He Jinlong， Yang Yang， Liu Xiaoxing. DEM simulation of single particle breakage under multi-point loading［J］. China Powder Science and Technology， 2025， 31（3）： 89-100.
［22］史振宇， 刘晓文， 宋来聪， 等. 航空发动机应用领域粉末高温合金的研究进展［J］. 中国粉体技术， 2025， 31（1）： 46-60.
Shi Zhenyu， Liu Xiaowen， Song Laicong， et al. Research progress of powder metallurgy superalloy in aero engine applications［J］. China Powder Science and Technology， 2025， 31（1）： 46-60.
［23］凡凤仙， 许玥. 颗粒声凝并微观动力学离散元模型及数值模拟［J］. 中国粉体技术， 2026， 32（1）： 33-43.
Fan Fengxian， Xu Yue. Microkinetic model and numerical simulations of particulate acoustic agglomeration［J］. China Powder Science and Technology， 2026， 32（1）： 33-43.
［24］余龙， 庞冬冬， 陈双琪， 等. 立式粉体干燥器内二元湿颗粒气-固逆流过程的数值模拟［J］. 中国粉体技术， 2024， 30（2）： 82-95.
Yu Long， Pang Dongdong， Chen Shuangqi， et al. Numerical study on gas-solid countercurrent process of binary wet particles in a vertical powder dryer［J］. China Powder Science and Technology， 2024， 30（2）： 82-95.
［25］顾毅楠， 吴永泽， 俞建峰， 等. 涡流空气分级机补气结构设计及流场仿真［J］. 中国粉体技术， 2024， 30（5）： 158-170.
Gu Yinan， Wu Yongze， Yu Jianfeng， et al. Air supply structure design and flow field simulation of turbo air classifier［J］. China Powder Science and Technology， 2024， 30（5）： 158-170.
［26］章康康， 韩玉阁， 任登凤， 等. Mg/PTFE/Viton富燃烟火剂凝聚相燃烧产物的实验和数值分析［J］. 燃烧科学与技术， 2021， 27（5）： 521-528.
Zhang Kangkang， Han Yuge， Ren Dengfeng， et al. Experimental and numerical analysis of condensed-phase combustion products of Mg/PTFE/Viton fuel-rich pyrolants［J］. Journal of Combustion Science and Technology， 2021， 27（5）： 521-528.
