CC BY
129
МАШИНОСТРОЕНИЕ MECHANICAL ENGINEERING
https://doi.org/10.21122/2227-1031-2018-17-3-177-182 UDC 621.791
Calculation of Particles Flow Temperature during Plasma Spraying of Mixture Consisting of Self-Fluxing Powder and Ceramics
A. S. Kalinichenko4, Yu. K. Krivosheyev1), V. V. Meshkova1), O. G. Devoino1)
1)Belarusian National Technical University (Minsk, Republic of Belarus)
© Белорусский национальный технический университет, 2018 Belarusian National Technical University, 2018
Abstract. Plasma spraying is one of the most effective methods allowing both to restore worn surfaces of parts and create wear-resistant coatings on new parts aiming the increase of their service life. Properties of the produced coatings depend on number of parameters, such as a plasma temperature, a chemical and fractional composition of the sprayed mixture, a distance from a plasma torch to the surface of a part, etc. Mathematical modeling of the process can significantly reduce the cost of processing of technological modes and is widely used at present for a calculation of technological parameters. The paper is devoted to mathematical simulation aiming to determine an effect of the injected ceramics content on the change in a temperature of a particles flow, as well as finding the modes in which the particles of high-temperature ceramics will be in the liquid state when they are deposited on the surface of a product. A mathematical model of particles heating in plasma has been formulated and a system of equations has been compiled. The system of equations has been solved numerically in Mathcad by a standard procedure using the Rkadapt function. Calculations have been carried out for a volume concentration of Al2O3 ceramics in a mixture from 5 to 50 % and for a plasma temperature at the exit from the plasma torch in the range from 6000 to 10000 K. Calculations have shown that the concentration of ceramics does not significantly affect the temperature of a mixture. The temperature of the particles depends to a large extent on the temperature of the plasma and the diameter of particles. It has been determined that for the entire range of calculated values the temperature of the self-fluxing powder in contact with the substrate exceeds a melting point. Fractional particle size has a strong effect on the temperature of particles at the moment of contact with the substrate. The dependences of a temperature of the ceramic phase on the particle size at different concentrations and plasma temperature have been determined. Analysis of the coatings microstructures has shown a good correlation with the results of the calculation.
Keywords: plasma spraying, ceramic phase, mathematical modeling, temperature of particles, microstructure For citation: Kalinichenko A. S., Krivosheyev Yu. K., Meshkova V. V., Devoino O. G. (2018) Calculation of Particles Flow Temperature during Plasma Spraying of Mixture Consisting of Self-Fluxing Powder and Ceramics. Science and Technique. 17 (3), 177-182. https://doi.org/10.21122/2227-1031-2018-17-3-177-182
Расчет температуры потока частиц при плазменном напылении смеси самофлюсующегося порошка и керамики
Докт. техн. наук, проф. А. С. Калиниченко1), канд. физ.-мат. наук Ю. К. Кривошеев1), инж. В. В. Мешкова4, докт. техн. наук, проф. О. Г. Девойно1)
^Белорусский национальный технический университет (Минск, Республика Беларусь)
Реферат. Плазменное напыление является одним из эффективных методов, позволяющих как восстанавливать изношенные поверхности деталей, так и создавать износостойкие покрытия на новых деталях с целью увеличения срока
Адрес для переписки
Калиниченко Александр Сергеевич
Белорусский национальный технический университет
просп. Независимости, 65,
220013, г. Минск, Республика Беларусь
Тел.: +375 17 296-66-86
Н Наука
итехника. Т. 17, № 3 (2018)
Address for correspondence
Kalinichenko Alexander S.
Belarusian National Technical University
65 Nezavisimosty Ave.,
220013, Minsk, Republic of Belarus
Tel.: +375 17 296-66-86
их службы. Свойства создаваемых покрытий зависят от ряда параметров, таких как температура плазмы, химический и фракционный составы напыляемой смеси, расстояние от плазмотрона до поверхности и др. Математическое моделирование процесса позволяет значительно снизить стоимость отработки технологических режимов и широко используется в настоящее время для расчета технологических параметров. В настоящей работе была поставлена цель проведения математического моделирования для определения влияния содержания вводимой керамики на изменение температуры потока частиц, а также нахождения режима, при котором частицы высокотемпературной керамики будут в жидком состоянии при осаждении на поверхность изделия. Сформулирована математическая модель нагрева частиц в плазме и составлена система уравнений, которая решалась численно в пакете МаШСа<1 стандартной процедурой с использованием функции Кка(!ар1 Расчеты проводились для объемной концентрации керамики А1203 в смеси от 5 до 50 % и для температуры плазмы на выходе из плазмотрона в интервале от 6000 до 10000 К. Вычисления показали, что концентрация керамики не влияет значительно на температуру смеси. Температура частиц в большей мере зависит от температуры плазмы. Определено, что для всего диапазона расчетных величин температура самофлюсующегося порошка при контакте с подложкой превышает температуру плавления. Фракционный размер частиц оказывает сильное влияние на температуру частиц в момент соприкосновения с подложкой. Определены зависимости температуры керамической фазы от размера частиц при различных концентрациях и температуре плазмы. Анализ микроструктур покрытий показал хорошую корреляцию с результатами расчета.
Ключевые слова: плазменное напыление, керамическая фаза, математическое моделирование, температура частиц, микроструктура
Для цитирования: Расчет температуры потока частиц при плазменном напылении смеси самофлюсующегося порошка и керамики / А. С. Калиниченко [и др.] // Наука и техника. 2018. Т. 17. № 3. С. 177-182. https://doi.org/10. 21122/2227-1031 -2018-17-3-177-182
Increasing the wear resistance of working surfaces is an important scientific and practical task for surface engineering. One of the methods of the wear resistance increase of steel surfaces is plasma spraying of self-fluxing Nickel-chromium-based powder [1]. From the point of view of reducing of the coating cost while maintaining its high wear resistance, the practical interest is the introduction into the sprayed mixture of ceramic powder based on aluminum oxide [2].
The properties of the coatings created depend on several parameters such as plasma temperature, chemical and fractional composition of the sprayed mixture, the distance from the plasma torch to the part's surface etc. Experimental determination of optimal regimes of plasma spraying process requires big expenditures of time and money, so mathematical modeling is widely used to save time and money.
Mathematical simulation and computer experiments can significantly reduce the cost of the processing technological modes and determine optimal spraying parameters. For example, the simulation of the main parameters of plasma arc burner: non-equilibrium thermodynamic and chemical models of plasma formation, energy transfer during plasma motion, boundary conditions, etc. was performed [3].
It is known from practice that the deformation of particles and their cooling significantly affects the adhesion of the deposited layer to the base and
the occurrence of stresses. The conducted computer simulation allowed determining the optimal temperature conditions for the formation of high-quality coatings [4].
Various programs are developed to simulate the plasma spraying process. Thus, in the work [5] the efficiency of the IPS Virtual Point program developed in Fraunhofer-Chalmers Centre for calculating the thickness of the plasma coating on parts for the automotive industry is considered. The possibility of application of this program for the analysis of the thickness of the layer deposited is shown.
A model for assessing the thickness non-uniformity of the layer deposited during plasma spraying is presented in [6]. On the basis of the developed model a theoretical study of the influence of process parameters on the cross-sectional character of the layer is carried out.
To estimate the influence of the modulation mode of the plasma torch electrical parameters on the plasma spraying process, a mathematical model was developed. This model allowed optimizing parameters of the plasma torch and the deposition process. Using the model it was possible to predict process parameters when a strong and continuous coating on the part surface with good adhesion was formed [7]. The problems of mathematical modeling of different stages of the plasma spraying process are considered in [8-10], first of all, the pro-
HayKa
wTexHMKa. T. 17, № 3 (2018)
cess of formation of the structure and porosity of the coating.
Despite the large number of publications, mathematical simulations of plasma spraying process are of considerable interest, especially in relation to specific conditions, because there are no models that completely describe this process. For example, an important aspect is the analysis of changes in the temperature of the particles' flow in the presence of a ceramic phase, which is characterized by a high melting point.
So, to determine the influence of the volume fraction of ceramic powders in sprayed mixture on changes in the temperature of the particle flow it is advisable to carry out mathematical simulate-on that allows you to choose the optimal modes of plasma spraying without conducting a large number of experiments.
In this paper, the purpose of mathematical simulation was set to determine the effect of the content of the injected ceramics on the change in the temperature of the particles flow, as well as finding the modes in which the particles of high-temperature ceramics will be in the liquid state when deposited on the surface of the product. Taking into account the complexity of the mathematical description of the real plasma spraying process, we introduce a number of simplifications.
1. The particles are uniformly distributed over the volume of the flow and there is no density stratification during the flight.
2. A mixture of powders is injected into the plasma jet at the plasma torch slice.
3. There is no change in the trajectory of the plasma jet and powder material.
4. Plasma properties do not change when the temperature decreases.
The mixture of the sprayed material consists of a self-fluxing nickel-chromium powder and aluminum oxide powder Al2O3. The diameter of the self-fluxing powders in the calculations adopted was equal to 100 ^m and ceramic particles were in range from 25 to 100 ^m. The distance from the plasma gun to the surface of part is 100 mm.
In the flow of plasma-forming gas (nitrogen), the speed of which is 47 m/s and the volume flow rate is 40 l/min, particles of self-fluxing Nickel-chromium powder are introduced with a mass flow rate of 5 kg/h. To improve the wear resistance, a ceramic powder based on aluminum oxide was
injected into the self-fluxing powder. Calculations were carried out for the volume concentration of Al2O3 ceramics in a mixture of 5 to 50 % and for the plasma temperature at the exit from the plasma torch in the range of 6000 to 10000 K.
Equations for the change of particles' temperature and gas one taking into account convective heat exchange have the following form:
mlcludT- = aF(T3 -T1 ); dy
m2c2u —2 = aF(T3 - T2 ); dy
(1)
(2)
irjl
Mc3u-^ = aF [n1(T1 - T3) + n2(T2 - T3)], (3)
where m1, m2 - masses of particles; c1, c2 - specific heat capacity of materials; T1, T2 - temperature of particles; F - surface area of the particle; M, c3, T3 - mass flow, specific heat and temperature of the gas, correspondingly; a - coefficient of convective heat transfer; u - flow speed of gas and particles; y - distance from the nozzle exit of the plasma torch; n1, n2 - number of particles of self-fluxing alloy and ceramics that interact with the gas flow; subscripts 1, 2, 3 - devoted to self-fluxing powder, ceramic powder and gas, correspondingly.
The form of equations (1) and (2) is close to the equations given in [1]. Calculations for the particle velocity using the expressions given in [2] showed that the particle velocity does not differ from the gas flow velocity. The convective heat exchange coefficient a can be written as follows
a = — Nu, D
where A - thermal conductivity of the gas; D -diameter of particles; Nu - Nusselt number that determines the convective heat exchange of the particle surface with the gas flow around it.
In the case of a spherical particle flow by a uniform gas stream with constant properties it is possible, as in [2], to use known dependence of the Ranz - Marshall
Nu = 2 + 0,6 Re0,5 Pr0,33.
■ HayKa
urexHMKa. T. 17, № 3 (2018)
Equations (1)-(3) can be simplified by substituting the values of mass and surface area of particles:
dTi _ 6a (T -T);
dy c1 uD1
dT2 - t2 );
-t * ^ v 3 2 / '
dy c2uD2
dT3 6am
dy c3uM p1
1 x(Tj -T3)+-X-(t2 -T)
D
D,
The resulting system of equations was solved numerically in the MathCad package by a standard procedure using Rkadapt function - the solution of a system of ordinary differential equations by Runge - Kutt method with automatic step selection.
In the case when the replacement of self-fluxing powder with ceramic one is carried out by volume fractions, the mass of the material introduced into the flow is decreased, since the density of the main component (self-fluxing alloy) is more than twice the density of ceramics. It is worth to note that the volume specific heat capacities (product of density and specific mass heat capacity) differ slightly, as:
PC = c;= 4,8332 kJ/(m3 • K); p2c2 = c* = 4,7697 kJ/(m3 • K).
That is, the difference in the products of density and heat capacity does not exceed 1.5 %. Therefore, the content of ceramics in the mixture, as shown by calculations, does not significantly affect the nature of the temperature distribution and its values (fig. 1). By increasing the volume content of ceramics in the mixture from 5 to 50 %, the temperature of aluminum oxide particles decreases only by 50-60 °C, that is, by 2.3-3.0 %. More significantly the temperature of the particles depends on the temperature of the plasma.
Analysis of calculation results (fig. 1) indicates that when the plasma temperature increases from 7000 to 10000 K, the temperature of the ceramic particles increases by about 800 degrees. At the size of ceramic particles of 100 microns, they reach a melting point for plasma temperature
not lower than 7400 K at a concentration of 5 % and 7700 K - for a concentration of 50 %.
Г.К 3100
1500
6000
7000
8000
9000
10 000 Tplasma, К
Fig. 1. Effect of plasma temperature on temperature of Al2O3 particles at the time of their contact with a sample surface for different volume content of ceramics in the mixture: 1 - 5 % of ceramics; 2 - 25 % of ceramics; 3 - 50 % of ceramics; particle diameter - 100 |m
Analysis of temperature change calculations at different distances from the plasma gun shows that a noticeable difference in the temperature of metallic and ceramic particles is observed at a distance of more than 5 cm from the outlet, when the temperature of the heated particles is higher by 10 degrees than the less heated ones. Temperature dependence of ceramic and self-flux particles is exponential (fig. 2). Moreover, for the entire range of calculated values, the temperature of the self-fluxing powder in contact with the substrate exceeds the melting point.
2000 T, K 1600
1400 1200 1000 800 600 400
■A
/
✓
1 2 3 4 5 6 7 81, cm 10
Fig. 2. Calculated temperature change in mixture components from the moment of entering into plasma to the contact with the substrate for particles of 100 microns in size and plasma temperature of 6500 K
Analysis of the calculation results shows that the average volumetric flow temperature decreases
Наука
итехника. Т. 17, № 3 (2018)
by about 15 degrees with an increase in the proportion of ceramic particles from 5 to 25 vol. %. The consequence of the decrease in the flow temperature can be deterioration in the adhesion of the developed layer to the substrate, as well as a weaker relationship between the particles of Nickel-chromium alloy and aluminum oxide.
Particle size has a strong effect on the temperature of the particles at the moment of contact with the substrate (fig. 3). For rplasma = 7000 K and the content of ceramics in the mixture is equal to 25 % the melting temperature is reached by only particles with size less than 40 microns. When the plasma temperature rises to 9000 K, the melting point is reached by particles with a diameter equal to or less than 63 microns. Moreover, when the plasma temperature rises, the dependence curve has a flatter character (fig. 3).
4000
Т, К 3000
2000
1000
25
50
75
d, |m 120
Fig. 3. Influence of ceramic particles size on their temperature at the moment of contact with the surface for volumetric
content of ceramics 25 %; 1 - Tplasma = 9000 K; 2 - Tplasma = 7000 K (dashed line shows the temperature of ceramics melting point)
To check the adequacy of the calculation results some experiments were carried out on the coating formation made of Nickel-based alloy (system Ni-Fe-Cr-Si-B-C) containing oxide ceramics using the plasma spray installation UPU-3D with the plasma torch PP-25. The plasma temperature was 7500 K.
Optical and electron microscopy were used for microstructure analysis. The microstructure of the coating with the volume content of ceramics 20 and 33 % is shown in fig. 4.
Analysis of microstructures shows that at such plasma temperature not all ceramic particles have melted (fig. 4a), that confirms the calculation results. From fig. 3 it follows that at such plasma temperature the size of the ceramic particles for complete melting should not exceed 50 microns,
■ Наука
итехника. Т. 17, № 3 (2018)
and in experiments the ceramics particles had diameter up to 63 microns. The increase in the volume content of ceramics leads to a decrease in the average temperature of the particle flow and, as a consequence, the presence of larger number of un-melted particles can be seen (fig. 4b).
а b
m
Fig. 4. Microstructure of plasma coatings with different content of ceramics: а - 20 %; b - 33 % (*500)
This is confirmed by SEM-microscopy as well. With the increase in ceramics concentration one can observe not only more non-uniform structure but also marks of tungsten particles (light particles in fig. 5b) confirmed by spectral analysis.
а b
, . «
- 'Wft, „^W
Fig. 5. SEM-microstructure of plasma coatings with different ceramics content: а - 15 %; b - 33 %
So, microstructural studies of deposited coatings are in good correlation with predicted temperature of particles establishing relations between some process parameters and particles temperature during deposition. Based on these results it is possible to estimate a formed coating structure which depends on temperature.
CONCLUSION
In the present work the mathematical simulations of temperature of the particles flow consisting of self-fluxing powder and oxide ceramics in plasma spraying process have been considered. It is determined that the plasma temperature and granulometric composition of ceramics are main factors for the formation of a high-quality coating.
0
REFERENCES
1. Nejat Y. S., Muharrem Y. (2008) Improvement of Wear Resistance of Wire Drawing Rolls with Cr-Ni-B-Si+WC Thermal Spraying Powders. Surface and Coatings Technology, 202 (13), 3136-3141. https://doi.org/10.1016/). surfcoat.2007.11.022.
2. Kalinichenko A. S., Devoino O. G., Meshkova V. V. (2015) Influence of Oxide Ceramics Content on Structure and Properties of Nickel-Chrome Plasma Coatings. Sov-remennye Metody i Tekhnologii Sozdaniya i Obrabotki Materialov: Sb. Nauch. Trudov. Kn. 2. Tekhnologii i Obo-rudovanie Mekhanicheskoi i Fiziko-Tekhnicheskoi Obra-botki [Modern Methods and Technologies for Development and Treatment of Materials: Collection of Research Papers. Book 2. Technologies and Equipment for Mechanical and Physical-Technical Treatment]. Minsk, Physical-Technical Institute of the National Academy of Sciences of Belarus, 171-174 (in Russian).
3. Trelles J. P., Chazelas C., Vardelle A., Heberlein J. V. R. (2009) Arc Plasma Torch Modeling. Journal of Thermal Spray Technology, 18 (5-6), 728-752. https://doi.org/10. 1007/s11666-009-9342-14.
4. Fukanuma H., Huang R., Tanaka Y., Uesugi Y. (2009) Mathematical Modeling and Numerical Simulation of Splat Cooling in Plasma Spray Coatings. Journal of Thermal Spray Technology, 18 (5-6), 965-974. https://doi.org/10. 1007/s11666-009-9366-6.
5. Berce A. (2011) Simulation of Thermal Spraying in IPS Virtual Paint. Master's Thesis in Solid and Fluid Mechanics. Sweden, Goteborg. Available at: http://publications. lib.chalmers.se/records/fulltext/152889.pdf. (Accessed 6 September 2017).
6. Sadovoy A. (2014) Modeling and Offline Simulation of Thermal Spray Coating Process for Gas Turbine Applications. Available at: http://tuprints.ulb.tu-darmstadt.de/4042/ 1/Modeling%20and%20offline%20simulation%20of%20t hermal%20spray%20coating%20process%20for%20gas% 20turbine % 20 applications% 20-%201400617d%20final.pdf (Accessed 6 September 2017).
7. Kadyrmetov A. M., Smolentsev E. V., Mal'tsev A. F., Sukhochev G. A. (2012) Modeling of Plasma Spraying Coating Process in the Mode of Plasma Arch Power Modulation Applied on Details of Transport Machines. Nauch-nyi Zhurnal KubGAU = Scientific Journal of KubSAU [Kuban State Agrarian University], 84 (10), 1-10.
8. Mostaghimi J. (2007) Understanding Plasma Spray Coating: a Modeling Approach. 18th International Symposium on Plasma Chemistry. Kyoto, Japan, August 26-31, 2007. Available at: https://plas.ep2.rub.de/ispcdocs/ispc18/ispc 18/content/slide00237.pdf. (Accessed 6 September 2017).
9. Ivanov E. M. (1983) Engineering Calculation of Thermal-Physical Processes during Plasma Spraying. Saratov, Publishing House of Saratov University. 138 (in Russian).
10. Boronenko M. P., Gulyayev I. P., Seregin A. E. (2012) Model of Motion and Heating in Plasma Jet. Vestnik Yugorskogo Gosudarstvennogo Universiteta = Yugra State University Bulletin, 25 (2), 7-15 (in Russian).
Received: 29.01.2018 Accepted: 05.03.2018 Published online: 29.05.2018
ЛИТЕРАТУРА
1. Nejat Y. S. Improvement of Wear Resistance of Wire Drawing Rolls with Cr-Ni-B-Si+WC Thermal Spraying Powders / Y. S. Nejat, Y. Muharrem // Surface & Coatings Technology. 2008. Vol. 202, No 13. P. 3136-3141.
2. Калиниченко, А. С. Влияние содержания оксидной керамики на структуру и свойства никельхромовых плазменных покрытий / А. С. Калиниченко, О. Г. Де-войно, В. В. Мешкова // Современные методы и технологии создания и обработки материалов: сб. науч. тр.: в 3 кн. / редкол. С. А. Астапчик (гл. ред.) [и др.]. Минск: ФТИ НАН Беларуси, 2015. Кн. 2: Технологии и оборудование механической и физико-технической обработки. С. 171-174.
3. Arc Plasma Torch Modeling / J. P. Trelles [et al.] // Modeling. Journal of Thermal Spray Technology. 2009. Vol. 18, № 5-6. P. 728-752. https://doi.org/10.1007/s11666-009-9342-14.
4. Mathematical Modeling and Numerical Simulation of Splat Cooling in Plasma Spray Coatings / H. Fukanuma [et al.] // Journal of Thermal Spray Technology. 2009. Vol. 18, No 5-6. P. 965-974. https://doi.org/10.1007/ s11666-009-9366-6.
5. Berce, A. Simulation of Thermal Spraying in IPS Virtual Paint. Master's Thesis. Chalmers Reproservice [Electronic Resource] / А. Berce. Sweden, Goteborg, 2011. Mode of access: http://publications.lib.chalmers.se/records/fulltext/ 152889.pdf. Date of Access: 06.09.2017.
6. Sadovoy, A. Modeling and Offline Simulation of Thermal Spray Coating Process for Gas Turbine Applications [Electronic Resource] / А. Sadovoy. Darmstadt, 2014. Mode of access: http://tuprints.ulb.tu-darmstadt.de/4042/ 1/Modeling%20and%20offline%20simulation%20of%20t hermal %20spray%20coating%20process%20for%20gas% 20turbine%20applications%20-%201400617d%20final.pdf. Date of Access: 06.09.2017.
7. Моделирование процесса плазменного напыления покрытий на детали транспортных машин в режиме модуляции мощности дуги плазмотрона / А. М. Кадыр-метов [и др.] // Научный журнал КубГАУ. 2012. Т. 84, № 10. С. 1-10.
8. Mostaghimi, J. Understanding Plasma Spray Coating: a Modeling Approach [Electronic Resource] / J. Mostaghimi // 18th International Symposium on Plasma Chemistry. Kyoto, Japan, August 26-31, 2007. Japan: Kyoto, 2007. Mode of Аccess: https://plas.ep2.rub.de/ispcdocs/ispc18/ispc18/ content/slide00237.pdf. Date of Access: 06.09.2017.
9. Иванов, Е. М. Инженерный расчет теплофизических процессов при плазменном напылении / Е. М. Иванов. Саратов: Изд-во Сарат. ун-та, 1983. 138 с.
10. Бороненко, М. П. Модель движения и нагрева в плазменной струе / М. П. Бороненко, И. П. Гуляев, А. Е. Серегин // Вестник Югорского гос. ун-та. 2012. Т. 25, вып. 2. С. 7-15.
Поступила 29.01.2018 Подписана в печать 05.03.2018 Опубликована онлайн 29.05.2018
Наука
итехника. Т. 17, № 3 (2018)
