Теплофизическое моделирование процесса нанесения внутритрубного покрытия методом центробежного СВС
Автор: Садыков Руслан Артурович, Потапов Виктор Иванович, Ермоленко Андрей Александрович, Трофимов Евгений Алексеевич
Журнал: Вестник Южно-Уральского государственного университета. Серия: Металлургия @vestnik-susu-metallurgy
Статья в выпуске: 2 т.13, 2013 года.
Бесплатный доступ
Разработана и реализована в форме программы для ЭВМ модель теплообмена между смесью для СВС и стенкой металлической трубы, реализующегося в ходе получения слоистых металлокерамических труб посредством центробежного СВС.
Теплообмен, математическое моделирование, свс-металлургия, центробежное свс-литье
Короткий адрес: https://sciup.org/147156841
IDR: 147156841
Текст научной статьи Теплофизическое моделирование процесса нанесения внутритрубного покрытия методом центробежного СВС
Способ получения композиционных слоистых металлокерамических труб с внешним металлическим и внутренним керамическим слоями с помощью реакций самораспространяющегося высокотемпературного синтеза (СВС) разработан в 70-х годах прошлого века. Основоположниками этого направления стали классики применения метода СВС: Osamu Odawara (Япония) и А.Г. Мержанов с сотрудниками (СССР).
Стальные композиционные трубы, покрытые внутренним слоем металлокерамики, используются для транспортировки абразивных материалов в угольной и металлургической промышленности, теплоэнергетике и т. п., поскольку обладают устойчивостью к эрозии, коррозии, механическому износу, нагреву [1–3].
Чтобы с помощью метода центробежного СВС получить покрытие, удовлетворяющее эксплуатационным требованиям, необходимо обеспечить оптимальные условия его синтеза. Поэтому возникает необходимость создания математической модели процесса нанесения внутритрубного покрытия методом центробежного СВС, которая бы учитывала все множество теплофизических и физико-химических параметров, формирующих это покрытие.
В работе [4] при создании математической модели было принято, что теплообмен происходит только радиально, и не учитывалась скорость фронта волны горения, от которого во многом зависят условия протекания процесса. Задачей настоящего исследования стала разработка математической модели центробежного СВС, учитывающей теплофизические факторы процесса, в том числе скорость распространения фронта горения.
В процессе формирования покрытия методом центробежного СВС по длине трубы образуется три области: сухой смеси, жидкой и твердой фаз покрытия. В общем виде математическая модель распределения температур в ходе центробежного СВС должна учитывать теплообмен между среда- ми: порошкообразной смесью для СВС; твердой фазой, образующей покрытие; жидкой фазой, металлом стенки трубы, а также газовой фазой внутри трубы и снаружи.
В настоящей работе представлена методика, использовавшаяся для моделирования теплообмена между стенкой трубы и смесью для СВС в процессе прохождения фронта горения вдоль трубы [5, 6].
При составлении математического описания процесса теплообмена сделаны следующие допущения: толщина слоя смеси СВС постоянна по радиусу и по длине трубы, теплофизические пара- метры сред постоянны и не зависят от температуры, фронт волны горения распространяется в режиме устойчивого, стационарного горения (в поршневом режиме).
С учетом принятых допущений уравнения энергии для стенки трубы и смеси, начальные и граничные условия примут вид:
<
^^ = 7 -T ) + dt XCT,S 1 ( S 1 ст )
+ Y (T -T } + a д Т'1 х ст,см V см ст / a ст _ 2 ;
dX ат ат см см д + Ur Д = Хсм,ст (Тст Тсм) + д t дx,
+ V А +а д Тсм +__Q■
+ хсм, s 2 (Ts 2 тсм ) + a см _ 2 '
’ d X 2 С см Р см
Т ст ( x , 0) = Ф ст ( x ), Т см ( x , 0) = Ф см ( x );
Т ст (0, t ) = Ф ст ( t ), Т см (0, t ) =ф см ( t);
^ Хст етст^0,!) = аст (Тст(0,t)-Тсм(0,t));(3)
a x
^ см д Т см(02 ^ )- = « см,ст ( Т см (0, t ) - Т ст (0, t )),
I дx’ где Tст , Tсм , Ts1 , Ts2 – температуры стенки, смеси, окружающей среды снаружи и внутри трубы соответственно;
_ а ст, s 1 p ст, s 1 _ а ст,см p ст,см
Х ст, s 1 = Z ; Х ст,см = Z ;
р ст С ст S ст р ст c ст S ст
_ а см, s 2 p см, s 2 _ _ а см,ст p см,ст
Х см, s 2 = Z ; Х см,ст = Z ;
рсм c см S см рсм Ссм S см aст, aсм – температуропроводности стенки, смеси;
U r - скорость продвижения фронта горения;
а ст, s 1 , а ст,см , а см, s 2 , а см,ст - коэффициенты теплоотдачи от стенки к внешней окружающей среде, от стенки к смеси, от смеси к внутренней окружающей среде, от смеси к стенке;
p ст, s 1 , p ст,см , p см, s 2 , p см,ст – периметры раздела между стенкой и внешней окружающей средой, между стенкой и смесью, между смесью и внутренней окружающей средой, между смесью и стенкой;
р ст, р см - плотности стенки, смеси;
c ст, c см – удельные теплоемкости стенки, смеси;
S ст, S см – площади поперечного сечения стенки трубы и смеси;
Q – мощность источника энергии;
Х ст, Х см - теплопроводность стенки трубы, смеси.
Точный расчет количества теплоты, выделяющейся в процессе СВС, достаточно сложен, так как количество выделяющейся теплоты будет существенно зависеть от температуры протекающего процесса. Представление о характере и количественных характеристиках измерений, происходящих с величиной теплового эффекта в ходе повышения температуры процесса, дает информация, представленная на рис. 1, которая получена с помощью баз данных FactSage, представленных в открытом доступе .
Упрощенно можно рассчитать тепловой эффект реакции как сумму средневзвешенных значений тепловых эффектов реакций алюминия с оксидами в составе смеси.
Для численного решения краевой задачи (1)–(3) использовали метод конечных разностей. Систему (1), условия (2), (3) представили в дискретном блочно-матричном виде:
Т П + 1 = MT m + NTS + Q s + 1 , (4)
где T m + 1, M , T m , N , TS , Q , I - блочные матрицы вида:
Г0 |
n + 1 " |
F |
G |
0 " |
|||||||||||
0 |
n + 1 2 |
A |
F |
C |
0 |
||||||||||
rp n + 1 Tm |
= |
, M |
= |
||||||||||||
0 |
n + 1 m - 1 |
0 |
A |
F |
C |
||||||||||
.0 |
n + 1 m |
_ |
L 0 |
G |
F |
||||||||||
Г |
0 1 n |
Г E |
0 |
" |
|||||||||||
0 2 |
0 |
E |
0 |
||||||||||||
n T m = |
, |
N = |
= |
, |
|||||||||||
0 n v m |
1 |
0 |
E |
0 |
|||||||||||
0 n L m |
J |
. 0 |
E |
. |
|||||||||||
Г0 |
S " |
Гф" |
Г 1 1 " |
||||||||||||
0 |
S |
0 |
0 |
||||||||||||
T s = |
, |
Q s = |
, |
I |
= |
: |
. |
||||||||
0 |
S |
0 |
0 |
||||||||||||
.0 |
S _ |
L |
0 |
. 1 2 _ |
Элементы блочных матриц – матрицы коэффициентов 0 m , 0 m +1, A , G , E , F , C , V , 0 s , 1 1 , I 2 имеют следующий вид:
n
m
n стm n смm
, 0
n + 1
m
n + 1 ст m
n + 1 см m
, 0 S =
T s 1
L T s 2 J

1000 1200 1400 1600 1800 2000 2200 2400 2600 2800 T, К
Рис. 1. Зависимость количества теплоты, выделяющейся в ходе реакции алюмотерми-ческого восстановления из одного килограмма стехиометрических смесей на основе различных оксидов, от температуры процесса
A =
Т a ст h 2
Т a см +TU r h 2 2 h
' Т Х ст, 5 1 ТХ ст,СМ
E =
ТХ ст, 5 1
ТХ см, 5 2
I 2 =
Тaст «ст,™ / _ n _ n стm смm h ^ст V
h а см,ст f Т а см Ти г )
I+"'
n стm
F =
2т а ст h 2
ТХ ст,см
¥ =
C =
1 1 =
ТХ см,ст
Q _,
Т a ст h 2
G =
h 2
1 Чем, 5 2 Чем, ст
Т a см Тиг h2 2h
2 т а см h 2 _
2 т a см h 2
ст ст,см (@ ” @ n
(^ ст1 ^ с м| )
h а см,ст f т а см ТЦ ) n
1,2 1 1 ^ см1
^ см I h 2 2 h Г
- ®стП )
Для адаптации математической модели использовали интервальный метод для определения параметра c см .
Построенная модель теплообмена между смесью для СВС и стенкой трубы реализована в виде программы для ЭВМ, написанной на языке программирования C# с использованием возможностей параллельного программирования. Программа позволяет рассчитывать распределение температур сред, контактирующих в ходе процесса СВС по длине трубы, в зависимости от времени, прошедшего с момента начала процесса (примеры результатов таких расчетов, проведенных для условий, представленных в таблице, показаны на рис. 2–4).
В процессе дальнейших работ по созданию полной модели процесса планируется разработать модели других частных процессов теплообмена, реализующихся в процессе СВС, а затем объединить их.
Данные, использованные в процессе расчета
Параметры |
Стенка |
Смесь СВС |
||
Теплофизические |
Плотность, кг/м3 |
7800 |
4000 |
|
Теплоемкость, Дж/(кг·К) |
500 |
1000 |
||
Теплопроводность, Вт/(м·К) |
47 |
1,00 |
||
Теплопередача, Вт/(м2·К) |
Внешняя среда |
1000 |
10000 |
|
Труба/Смесь |
30 000 |
30 000 |
||
Геометрические |
Внешний радиус трубы, м |
0,15 |
||
Внутренний радиус трубы/внешний радиус смеси, м |
0,14 |
0,14 |
||
Внутренний радиус смеси, м |
– |
0,09481 |
||
Периметр раздела сред, м |
Внешняя среда |
0,94248 |
0,59573 |
|
Труба/Смесь |
0,87964 |
0,87964 |
||
Площадь сечения, м2 |
0,00911 |
0,03333 |
||
Механические |
Скорость горения, мм/с |
– |
20,0 |
|
Скорость вращения, мм/с |
1500 |
1500 |

Рис. 2. Пример результатов расчета распределения температур взаимодействующих сред в сечении трубы ( x = 0,15 м) в зависимости от времени, прошедшего от начала процесса

Рис. 3. Пример результатов расчета распределения температур взаимодействующих сред вдоль трубы спустя 10 с от начала процесса

Рис. 4. Пример результатов расчета распределения температур взаимодействующих сред вдоль трубы спустя 150 с от начала процесса
Выводы
В рамках работ по построению теплофизической математической модели процесса получения слоистых металлокерамических труб с внешним металлическим и внутренним керамическим слоями с помощью СВС разработана и реализована в форме программы для ЭВМ модель теплообмена между смесью для СВС и стенкой металлической трубы.
Часть работы, связанная с термохимическими расчетами процесса СВС, выполнена при поддержке РФФИ, грант № 13-08-00545.
Список литературы Теплофизическое моделирование процесса нанесения внутритрубного покрытия методом центробежного СВС
- Концепция развития СВС как области научно-технического прогресса/под ред. А.Г. Мержанова. -Черноголовка: Территория, 2003. -368 с.
- Мержанов, А.Г. Твердопламенное горение/А.Г. Мержанов, А.С. Мукасьян. -М.: Торус Пресс, 2007. -336 с.
- Санин, В.Н. СВС-металлургия труб с износостойким защитным покрытием с использовани ем техногенных отходов металлургических производств/В.Н. Санин, Д.Е. Андреев, В.И. Юхвид//Изв. высш. учеб. заведений. Порошковая металлургия и функциональные покрытия. -2011. -№ 2. -С. 37-43.
- Wang, Y.-F. Finite element analysis of residual thermal stress in ceramic-lined composite pipe prepared by centrifugal-SHS/Yu-Fei Wang, Zhen-Guo Yang//Materials Science and Engineering: A. -2007. -Vol. 460-461. -P. 130-134.
- Потапов, В.И. Математические модели теплофизических процессов в объектах многослойной структуры: моногр./В.И. Потапов. -Челябинск: Изд-во ЮУрГУ, 2004. -270 с.
- Ермоленко, А.А. Математическая модель теплообмена при нанесении внутритрубного покрытия методом центробежного СВС/А.А. Ермоленко, Р.А. Садыков, В.И. Потапов//Материалы IX Междунар. науч.-практ. конф. «Современные научные достижения -2013». Т. 68. Секция: Математика/Издат. дом «Образование и наука». -Прага, 2013. -72 с.
- Kontseptsiya razvitiya SVS kak oblasti nauchno-tekhnicheskogo progressa [The Concept of SHS as a Field of Scientific and Technical Progress] Ed. by A.G. Merzhanova. Chernogolovka, Territoriya, 2003. 368 p.
- Merzhanov A.G., Mukas'yan A.S. Tverdoplamennoe gorenie [Solid-Flame Combustion]. Moscow, Torus Press, 2007. 336 p.
- Sanin V.N., Andreev D.E., Yukhvid V.I. SHS Metallurgy Tubes with Wear-resistant Protective Coating with the Use of Industrial Waste Iron and Steel Industries [SVS-metallurgiya trub s iznosostoykim zashchitnym pokrytiem s ispol'zovaniem tekhnogennykh otkhodov metallurgicheskikh proizvodstv]. Izvestiya vuzov. Poroshkovaya metallurgiya i funktsional'nyepokrytiya, 2011, no. 2, pp. 37-43.
- Wang Y.-F., Yang Z.-G. Finite Element Analysis of Residual Thermal Stress in Ceramic-Lined Composite Pipe Prepared by Centrifugal-SHS. Materials Science and Engineering A, 2007, vol. 460-461, pp. 130-134 DOI: 10.1016/j.msea.2007.01.017
- Potapov V.I. Matematicheskie modeli teplofizicheskikh protsessov v ob"ektakh mnogosloynoy struktury [Mathematical Models of Thermal Processes in the Objects of the Multilayer Structure]. Chelyabinsk, SUSU Publ., 2004. 270 p.
- Ermolenko A.A., Sadykov R.A., Potapov V.I. Mathematical Model of Heat Transfer for Application-Tube Coating by Centrifugal SHS [Matematicheskaya model' teploobmena pri nanesenii vnutritrubnogo pokrytiya meto-dom tsentrobezhnogo SVS]. Materialy IX Mezhdunarodnoy nauchno-prakticheskoy konferentsii “Sovremennye nauchnye dostizheniya -2013 ”. Tom 68. Sektsiya: Matematika [Proceedings of the IX International Scientific-Practical Conference “Modern Scientific Achievements -2013”. Vol. 68. Section: Mathematics]. Praga, Education and Science, 2013, p. 72.