## Influence of micropores on structural instability of the combustion wave

**Authors:**Cui H.^{1}, Grigoryevskaya A.A.^{2}, Gulyaev P.Y.^{3}**Affiliations:**- Shandong University of Science and Technology
- Altai State University
- Ugra State University

**Issue:**Vol 15, No 4 (2019)**Pages:**33-40**Section:**New materials and technology**URL:**https://journals.eco-vector.com/byusu/article/view/18941**DOI:**https://doi.org/10.17816/byusu2019433-40- Cite item

#### Abstract

The article shows examples of visualization of the process of heat transfer by radiation in unstable combustion modes of porous powder materials, which are in good agreement with the results of high-speed video recording and micropyrometry. The mathematical model and the results of calculating the structure of the combustion wave in the Ni-Al system are presented. The contribution of radiative heat transfer at an adiabatic combustion temperature in the range of 8–10 % and the effect of its trigger shutdown with decreasing temperature were revealed.

# Introduction

The development of optoelectronic methods for contactless monitoring of combustion and explosion processes [1–3], based on modern television measuring systems of nanosecond resolution [4–6], allows to ensure the stability of technological modes and reproducibility of predetermined properties of materials in the process of obtaining them [7–9]. The high accuracy of micropyrometry of fast moving particles [10; 11; 14], based on image recording in the charge accumulation mode, opens up the possibility of using new methods of panoramic and hyperspectral analysis in the study of explosion welding [12], heat and mass transfer in plasma [8; 13; 14] and detonation spraying [7]. The difference between the technology of self-propagating high-temperature synthesis (SHS) is that it belongs to the methods of out-of-furnace powder metallurgy and the combustion wave propagates along the charge of an exothermically reacting powder mixture [16].

**The aim of the study**, based on the known models of particle packing under conditions of natural bulk gravimetric density [17; 18], is to generate the fractal structure of the mixture of initial powders for which visualization of heat transfer in the combustion wave for the 2-dimensional model taking into account local porosity and radiation heat transfer.

# Description of the problem

The foundations of computer simulation of combustion modes were laid in [16; 18], where the model consists of the heat equation and the kinetics equation:

$c{\rho}_{0}=\frac{\partial T}{\partial t}=\left(\frac{{\partial}^{2}T}{\partial {r}^{2}}+\frac{1}{2}\frac{\partial T}{\partial r}+\frac{1}{{r}^{2}}\frac{{\partial}^{2}T}{\partial {\phi}^{2}}+\frac{{\partial}^{2}T}{\partial {h}^{2}}\right)+{\rho}_{0}Q\frac{\partial \eta}{\partial t},\phantom{\rule{0ex}{0ex}}\frac{\partial \eta}{\partial t}=\left\{\begin{array}{r}{k}_{0}(1-\eta )exp(-\frac{E}{RT})\eta 1,\\ 0\eta \ge 1\end{array}\right\},\phantom{\rule{0ex}{0ex}}$ (1)

with the generally accepted designation of physical quantities, except for cylindrical coordinates: *r* is the radius, *φ* is the polar angle, *h* is the layer height.

The 2-dimensional model described here uses temperature for several layers of porous powder packaging:

$\frac{\partial T}{\partial t}=k\left(\frac{{\partial}^{2}T}{\partial {x}^{2}}+\frac{{\partial}^{2}T}{\partial {y}^{2}}\right)+f;\phantom{\rule{0ex}{0ex}}\frac{\partial A}{\partial t}=L\left(\frac{{\partial}^{2}A}{\partial {x}^{2}}+\frac{{\partial}^{2}A}{\partial {y}^{2}}\right);\phantom{\rule{0ex}{0ex}}\frac{\partial B}{\partial t}=M\left(\frac{{\partial}^{2}B}{\partial {x}^{2}}+\frac{{\partial}^{2}B}{\partial {y}^{2}}\right);\phantom{\rule{0ex}{0ex}}L=\left\{\begin{array}{l}0:T(x,y,t)<0.2\\ 0.01\times T(x,y,t):T(x,y,t)\ge 0.2\end{array}\right.\phantom{\rule{0ex}{0ex}}M=\left\{\begin{array}{l}0:T(x,y,t)<0.2\\ 0.01\times T(x,y,t):T(x,y,t)\ge 0.2\end{array}\right.\phantom{\rule{0ex}{0ex}}k=0.01\phantom{\rule{0ex}{0ex}}$ (2)

Here are indicated: *T* – layer temperature; *A _{Ni}*,

*B*– the corresponding concentration of the components of the mixture in the layer; the coefficients k, L, and

_{Al}*M*can be set in the whole table space from the generated charge structure file [18], while the system is considered closed and adiabatic, and the heat capacity of pores is neglected. Figure 1 shows three versions of the reaction media: a clad powder (

*A*), a mixture of powders of equal size (

*B*), and a model for the Ni-Al system (

*C*). The combustion wave propagates from above, the combustion region is highlighted in red, blue and green color correspond to the initial components of the mixture, yellow to the final synthesis products. So we get the configuration of the combustion front with the allocation of coordinates of hot particles, the temperature of which is dimensionless by the value of the maximum adiabatic temperature.

**Figure 1. Visualization of the front of the SHS reaction wave in 2-dimensional models**

Cell numbers of the graphic screen with the coordinates of the hot particles are listed. After forming the list of boundary points for each of them, a list of the remaining boundary points of the cold particles of the starting products visible from it is compiled. This is checked by the criterion for the intersection of the line connecting the selected points with any other. If the line intersects, the point is considered invisible. When calculating, according to the Stefan-Boltzmann law: S = *σT*^{4}, irradiation heat transfer for all boundary points, the correspondence of particle heating above the threshold temperature is checked. To do this, we believe that all effective heat transfer is carried out at a wavelength corresponding to the maximum density of thermal radiation, which is calculated according to the Wien displacement law: λ_{max} = *b* *T*^{-1}. Given that, the propagation of thermal radiation is limited in pores with a minimum diameter d_{min}, in accordance with the Kirchhoff diffraction limit: d_{min} = λ_{max} / (2*n*), restrictions are imposed on the temperature of the hot powder particles. If it is above the threshold, then it is considered that such a particle participates in radiation heat transfer and this part of the heat is transferred to points that it can irradiate. Figure 2 shows the effect of radiation heat transfer on the burning rate and the geometric shape of the reaction wave front. The results of visualization of three scenarios are presented: L (“low”) – there is no irradiative heat transfer, for example, in an ultrafine powder; M ("medium") – radiation heat transfer of the order of 10 %, for example in ordinary fine powders; H (“high”) – the prevailing radiation heat transfer of more than 50 %, for example, in loose and poorly packed powder mixtures. The numbers indicate the iteration step number when simulating the propagation of a combustion wave. It can be seen that the speed increases significantly due to warming up in the cavities of the embankment. In the complete absence of heat transfer by radiation, combustion propagates along a chain of contacting particles.

**Figure 2. The effect of radiation heat transfer on the speed and stability of the combustion front**

# Pore heat transfer visualization

In the process of calculating radiation heat transfer, an array of transparency of points (obstacles) is used, the output array where the brightness, brightness of the starting point is added. The neighbors (branches) of this point are moving. Each point stores its angle and deviation from it. The sequence of the algorithm is shown in Figure 3. In the process of calculating radiation heat transfer, an array of point transparency (obstacles) is used, the output array where the brightness, brightness of the starting point is added. The neighbors (branches) of this point are moving. Each point stores its angle and deviation from it:

K_{a1} = a1 – b1; (minimum);

K_{b1} = a1 + b1; (maximum).

According to them, for each point, a structure is formed storing the cone (beam) and value (beam intensity). The calculation of this branch with the input parameters of the formed structure is started. To visualize the intensity of thermal radiation, the Bayesian color optimization technique was used [19].

**Figure 3. Visualization of the indicatrix of radiation of a hot particle in pore space**

The following input parameters are used to calculate the branch: point (p_{Top}), cone (fK_{b1}, fK_{a1}) and beam intensity (f_{Val}, initial beam). Zero-weighted points are skipped. The width of its input cone (window) is calculated: f_{Range} = fK_{b1}-fK_{a1}. If the window is zero, then the point is skipped. Next, all branches of this point are sorted. If the cone (window) of the branch (p_{The}) does not intersect with the cone specified by the input parameters (ray), then the point is skipped (the ray does not fall into it). If the cones intersect (beam and window), a new cone is formed from the intersection of the direction areas (the part that passes through the window is cut out of the beam). For the resulting cone (passing through the window), its width is considered: f_{Th} = f_{Max}-f_{Min}; The relation to the width of the initial ray is calculated: f_{K} = f_{Th}/f_{Range}; The value is added to the storage of received radiation at this point: + Δ = f_{Val} * f_{K}, where f_{Val} is the intensity of the initial beam. The calculation of its branches starts, with the parameters: this point (p_{The}), a new cone (f_{Max}, f_{Min}), a new beam intensity (f_{K} * f_{Val} * f_{0pasity}). Under f_{0pasity} is indicated the transparency of this point. If the beam intensity falls below the threshold, its branches are not considered. The result is illustrated in Figure 4, where it can be seen that the space where the study is going is looped in image Fig. 4c. Visualization of the ignition mechanism of the lower layer of particles due to radiation heat transfer is presented in Figure 5.

**Figure 4. The sequence of formation of the shadow cone and the radiation window**

**Figure 5. The mechanism of initialization of the relay mode of combustion by thermal radiation**

# Computer Model Verification

The adequacy of the results of computer modeling by equations (1) and (2) was checked using the porous structure of the initial mixture of Ni and Al powders as an example, as shown in Figure 6.

**Figure 6. Verification of the computer model by the real packing structure of powders**

The computer model was compared in terms of open pycnometric porosity with a mixture of powders of carbonyl nickel, PNK-UT2 grade with a particle size of 3 to 10 microns, and aluminum, PA-4 grade with a particle size of 50 to 100 microns. As a result, the calculated values of porosity for a computer model were obtained in the range from 60 to 65 %. The experimentally measured porosity of a real powder mixture was 65–67 %, which indicates a good adequacy of the computer model.

The successive stages of the computer algorithm when visualizing the process of thermal radiation propagation in open pore space are shown in Figure 7.

**Figure 7. Computer visualization of thermal radiation in the SHS combustion wave**

Based on the results of computer simulation, it can be argued that thermal radiation in the combustion wave freely fills all main pores to a depth of 150-200 microns. The boundaries of the localization of thermal radiation in the region of the combustion front are pores with sizes smaller than the Kirchhoff diffraction limit: d_{min} = λ_{max} / (2n). Obviously, an increase in the effective wavelength of thermal radiation λ_{max} with a decrease in the combustion temperature leads to a collapse of the localization region when the size of the largest pore becomes less than 10λ_{max}. With a characteristic heterogeneity scale Δ = 20 μm, for a real mixture of powders considered in our case, in accordance with the Wien law λ_{max} = *b T ^{-1}*, the value of the “trigger switch-off” of radiation thermal conductivity will be 1500 K [22; 24; 25].

# Conclusions

- The method of computer visualization of radiation heat transfer in the SHS combustion wave is in good agreement with new data and experimental results obtained using high-speed micropyrometry and thermal imaging registration of such processes in [20; 21; 23].
- The effect of radiation heat transfer on the curvature of the combustion front (see Fig. 1C) discovered in the work can be used to develop methods for analyzing the thermal diffusion instability of combustion in the Zel'dovich – Barenblatt model for highly porous powder mixtures.
- The contribution of radiation heat transfer at an adiabatic combustion temperature in the range of 8–10 % and the effect of its trigger switching off due to the Kirchhoff diffraction threshold at a temperature below 1500 K.

### Hong Zhi Cui

Shandong University of Science and Technology
**Author for correspondence.**

Email: ugrasu@ugrasu.ru

China, Qingdao

Distinguished professor of Taishan Scholar, Dean of College of Materials Science and Engineering

### Angelica A. Grigoryevskaya

Altai State University
Email: a.a.grig@bk.ru

Russian Federation, Barnaul

Bachelor of the Faculty of Mathematics and Information Technology

### Pavel Yu. Gulyaev

Ugra State University
Email: p_gulyaev@ugrasu.ru

Russian Federation, Khanty-Mansiysk

Doctor of Technical Sciences, Leading Researcher, Institute of Oil and Gas

- Dolmatov, A. V. Thermal imaging complex with tracking function for joint research of microheterogeneous processes and macrokinetics of SHS phenomenon / A. V. Dolmatov, P. Yu. Gulyaev. – doi: 10.1088/1742-6596/1333/6/062006 // Journal of Physics: Conference Series. – 2019. – Vol. 1333. – P. 062006.
- Гуляев, П. Ю. Физические принципы диагностики в технологиях плазменного напыле-ния / П. Ю. Гуляев, А. В. Долматов. – Текст : непосредственный // Известия Самарского научного центра Российской академии наук. – 2009. – Т. 11. – № 5-2. – С. 382–385.
- Гуляев, П. Ю. Методы оптической диагностики частиц в высокотемпературных потоках / П. Ю. Гуляев, А. В. Долматов, Г. Н. Леонов. – Текст : непосредственный // Ползунов-ский вестник. – 2012. – № 2-1. – С. 4–7.
- Increasing the noise immunity of optical-electronic systems based on video cameras with an optical converter / M. P. Boronenko, P. Yu. Gulyaev, A. E. Seregin, K. G. Poluhina. – doi: 10.1088/1742-6596/643/1/012028 // Journal of Physics: Conference Series. – 2015. – Vol. 643. – P. 012028.
- Dolmatov, A. V. Investigation of structure formation in thin films by means of optical pyrometry / A. V. Dolmatov, I. V. Milyukova, P. Y. Gulyaev. – doi: 10.1088/1742-6596/1281/1/012010 // Journal of Physics: Conference Series. – 2019. – Vol. 1281. – P. 012010.
- Иордан, В. И. Оптико-электронная система диагностики двухфазных потоков динамиче-ским методом счета частиц / В. И. Иордан. – Текст : непосредственный // Известия выс-ших учебных заведений. Физика. – 2008. – Т. 51. – № 9-3. – С. 79–87.
- Гуляев, П. Ю. Автоматизация контроля теплофизических параметров в технологиях де-тонационного напыления / П. Ю. Гуляев, А. В. Долматов. – Текст : непосредственный // Системы управления и информационные технологии. – 2009. – Т. 35. – № 1-2. – С. 230–233.
- Методы контроля температуры и скорости частиц конденсированной фазы в процессе плазменно-дугового напыления / М. П. Бороненко [и др.]. – Текст : непосредственный // Фундаментальные исследования. 2013. – № 10-6. – С. 1194–1199.
- In-situ selfpropagating-hightemperature-synthesis controlled by plasma / H. Z. Cui, P. Yu. Gulyaev, I. P. Gulyaev, I. V. Milyukova // Bulletin of ugra state university. – 2012. – № 2 (25). – P. 28–33.
- Gulyaev, P. Instability of the Ni-Al combustion wave in the Zeldovich-Barenblatt parameters / P. Gulyaev, A. Dolmatov, V. Jordan. – doi: 10.1088/1742-6596/1353/1/012036 // Journal of Physics: Conference Series. – 2019. – Vol. 1353. – № 1. – P. 012036.
- Иордан, В. И. Разрешающая способность виртуальных приборов контроля температуры частиц в плазменных потоках по суммарному спектру / В. И. Иордан. – Текст : непосред-ственный // Ползуновский альманах. – 2008. – № 2. – С. 13–14.
- Saikov, I. V. Shock-wave synthesis in powder mixtures / I. V. Saikov, M. I. Alymov, P. Y. Gulyaev. – doi: 10.1088/1742-6596/1115/4/042012 // Journal of Physics: Conference Series. – 2018. – Vol. 1115. – P. 042012.
- Бороненко, М. П. Измерение скорости и температуры частиц в потоке низкотемператур-ной плазмы / М. П. Бороненко. – Текст : непосредственный // Известия высших учебных заведений. Физика. – 2014. – Т. 57. – № 3-2. – С. 70–73.
- Control of dispersed-phase temperature in plasma flows by the spectral-brightness pyrometry method / A. V. Dolmatov [at al.]. – doi: 10.1088/1757-899X/110/1/012058 // IOP Conference Series : Materials Science and Engineering. – 2016. – Vol. 110. – № 1. – P. 012058.
- Phase formation time evaluation in NiAl combustion systems by the thermal fields visualization method / M. P. Boronenko [at al.] // Scientific Visualization. – 2015. – Vol. 7. – № 5. – P. 102–108.
- Ivleva, T. P. Three-dimensional modes of unsteady solid-flame combustion / T. P. Ivleva, A. G. Merzhanov // Chaos. – 2003. – Vol. 13. – Iss. 1. – P. 80–85.
- Гуляев, П. Ю. Моделирование технологических процессов плазменного напыления по-крытий наноразмерной толщины / П. Ю. Гуляев, И. П. Гуляев // Системы управления и информационные технологии. – 2009. – Т. 35. – № 1-1. – С. 144-148.
- Моделирование фрактальных структур упаковок порошковых СВС-материалов / И. В. Милюкова, А. Л. Трифонов, С. А. Ширяев [и др.]. – Текст : непосредственный // Ползу-новский альманах. – 2007. – № 3. – С. 39–41.
- Гуляев, П. Ю. Байесовское восстановление цвета цифровых изображений / П. Ю. Гуляев, Ю. П. Гуляев, А. В. Долматов. – Текст : непосредственный // Вестник СГУГиТ. – 1997. – № 2. – С. 114–115.
- Boronenko, M. High-speed visualization of combustion synthesis discrete reaction waves: co-herent heat microstructures / M. Boronenko. – doi: 10.1088/1742-6596/1115/4/042021 // Journal of Physics: Conference Series. – 2018. – Vol. 1115. – P. 042021.
- Dolmatov, A. V. Investigation of structure formation in thin films by means of optical pyrometry / A. V. Dolmatov, I. V. Milyukova, P. Y. Gulyaev. – doi: 10.1088/1742-6596/1281/1/012010 // Journal of Physics: Conference Series. – 2019. Vol. 1281. – P. 012010.
- Gulyaev, P. Trace-Analysis of Images of the Differential Chronogram of the Combustion Wave for Recognition of Transitional Modes of SHS / P. Gulyaev, V. Jordan // CEUR Workshop Proceedings : Optoelectronic equipment and devices in systems of pattern recognition, image and symbol information processing. Recognition. – 2017. Vol. 1940. – P. 37–44.
- Borodina, K. Thermal analysis of reaction producing KXTiO2 / K. Borodina, S. Sorokina, N. Blinova. – doi: 10.1007/s10973-017-6840-0 // Journal of Thermal Analysis and Calorimetry. – 2018. – Vol. 131. – № 1. – P. 561–566.
- Гуляев, П. Ю. Trace-анализ дифференциальной хроноскопии волны горения СВС при высокоскоростной телевизионной регистрации / П. Ю. Гуляев, В. И. Иордан, В. В. Му-лярец. – Текст : непосредственный // Доклады ТУСУР. – 2018. – Т. 21. – № 1. – С. 62–67.
- Григорьевская, А. А. Компьютерная визуализация радиационного теплопереноса в волне горения СВС / А. А. Григорьевская, П. Ю. Гуляев. – Текст : непосредственный // Ползуновский альманах. – 2019. – № 4. – С. 5–9.

#### Supplementary files

Supplementary Files | Action | ||

1. | Figure 1. Visualization of the front of the SHS reaction wave in 2-dimensional models | View (244KB) | Indexing metadata |

2. | Figure 2. The effect of radiation heat transfer on the speed and stability of the combustion front | View (411KB) | Indexing metadata |

3. | Figure 3. Visualization of the indicatrix of radiation of a hot particle in pore space | View (189KB) | Indexing metadata |

4. | Figure 4. The sequence of formation of the shadow cone and the radiation window | View (126KB) | Indexing metadata |

5. | Figure 5. The mechanism of initialization of the relay mode of combustion by thermal radiation | View (155KB) | Indexing metadata |

6. | Figure 6. Verification of the computer model by the real packing structure of powders | View (122KB) | Indexing metadata |

7. | Figure 7. Computer visualization of thermal radiation in the SHS combustion wave | View (84KB) | Indexing metadata |