УДК 519.6:621.793+532.517+544.015.4
DOI: 10.18698/1812-3368-2019-6-18-39
ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ КРИСТАЛЛИЗАЦИИ МОДИФИЦИРОВАННОЙ МЕТАЛЛИЧЕСКОЙ КАПЛИ ПРИ РАСТЕКАНИИ НА ПОДЛОЖКЕ
В.Н. Попов А.Н. Черепанов
popov@itam.nsc.ru ancher@itam.nsc.ru
ИТПМ СО РАН, Новосибирск, Российская Федерация
Aннотация
Ключевые слова
Проведено численное моделирование процессов Численное моделирование, при падении капель расплавленного металла на металлическая капля, подлож-подложку. Рассмотрено затвердевание на металли- ка, модифицирование, сплав ческой поверхности бинарного сплава алюминия, алюминия, наноразмерные модифицированного активированными тугоплав- тугоплавкие частицы, кими наноразмерными частицами, которые являют- кристаллизация ся центрами зарождения кристаллической фазы. Сформулирована математическая модель, описывающая термо- и гидродинамические явления в капле при взаимодействии с твердой подложкой, гетерогенное зародышеобразование в процессе охлаждения расплава и последующую кристаллизацию. Течение в жидкости описано уравнениями Навье — Стокса в приближении Буссинеска. Положение свободной границы расплава зафиксировано частицами-маркерами, перемещающимися с локальной скоростью жидкости. На поверхности контакта расплав-подложка учтено термическое сопротивление. Гидродинамическая задача рассмотрена в условиях кристаллизации расплавленного металла. Описаны температурные режимы, кинетика роста твердой фазы в затвердевающем сплаве алюминия при различных размерах формирующихся сплэтов. Установлено удовлетворительное совпадение формы сплэта, полученного по результатам численного моделирования, с имеющимися экспериментальными данными. Подтверждена адекватность модели кристаллизации при наличии ультрадисперсных тугоплавких частиц в бинарном сплаве. Определено, что независимо от размера капли имеет место объемная кристаллизация металла. Установлено, что
при большой скорости соударения капли с подложкой за время растекания расплава возможно обра- Поступила 19.02.2019 зование небольшой доли твердой фазы © Автор(ы), 2019
Работа выполнена в рамках Программы фундаментальных научных исследований государственных академий наук на 2013-2020 годы (проект
№ 0323-2018-0012)
Введение. Одним из простых и доступных способов нанесения стойких защитных покрытий на поверхность деталей либо восстановления геометрии их изношенных поверхностей является газопламенное напыление расплавленных частиц различных (медных, алюминиевых, цинковых и др.) сплавов [1]. При газопламенном напылении происходит последовательная укладка отдельных сплэтов — растекшихся и затвердевших на обрабатываемой поверхности капель расплава, и очевидно, что обеспечение высоких механических и прочностных свойств покрытий во многом определяется их высокой дисперсностью и однородностью кристаллического зерна [2, 3].
Сплавы алюминия (А1 + Б1) широко используют в автомобильной и аэрокосмической отраслях. В настоящее время одну из главных проблем представляет повышение прочности алюминия, которая может быть увеличена механизмом измельчения зерна при добавлении в расплав ультрадисперсных тугоплавких частиц. В качестве модификаторов могут использоваться предварительно специально подготовленные тугоплавкие наноразмерные частицы (БЮ, Т1С, Т1Ы и др.), которые даже при относительно малой массовой доле в расплаве, не более 0,05 %, позволяют увеличить число центров кристаллизации [4, 5]. В настоящее время существует несколько публикаций, посвященных экспериментальным исследованиям влияния ультрадисперсных тугоплавких частиц в сплавах алюминия на измельчения зерна в слитке и повышение прочности затвердевшего металла [6-9]. Экспериментально подтверждено уменьшение среднего размера зерна в затвердевшем слитке, способствовавшее улучшению прочностных свойств металла [6]. В ходе исследований измерена температура в металле во время его кристаллизации [7, 10].
Обеспечение высоких механических и прочностных свойств алюминиевых покрытий на поверхностях деталей и конструкций во многом должно обусловливаться структурой нанесенных сплэтов. Для получения мелкого зерна в затвердевшем металле необходимо возникновение как можно большего числа центров кристаллизации в расплаве, что достигается вве-
дением в напыляемый материал модифицирующих компонентов [2]. Применение такой технологии дает возможность повысить эксплуатационные свойства восстанавливаемых деталей и их поверхностного слоя. При совершенствовании методов нанесения стойких защитных покрытий актуально детальное исследование происходящих процессов. В связи с этим использование численного моделирования часто является единственным способом, позволяющим получить представление о рассматриваемых процессах при сопоставлении получаемых результатов с имеющимися экспериментальными данными.
Результаты численных исследований неравновесной объемной кристаллизации модифицированной металлической капли на холодной подложке приведены в [11]. Предложена математическая модель и рассмотрена кинетика кристаллизации однокомпонентного расплава. Выявлено существование переохлаждения в начальной стадии кристаллизации, величина которого зависит от размера частиц и интенсивности охлаждения.
3Б-Моделирование процессов при послойном наложении металлических капель на поверхность подложки проведено в [12]. Получена оценка термического сопротивления на поверхности контакта расплава с подложкой. Результаты численного моделирования взаимодействия капель расплавленного алюминиевого сплава (А1 + Б1) и стальной пластины сравнивались с данными, полученными при экспериментах, однако кинетика кристаллизации, во многом влияющая на структуру затвердевшего металла, не рассмотрена.
Математическая модель кристаллизации двухкомпонентного сплава предложена в [13]. Однако по результатам анализа представленных в этой работе результатов определено, что модель не описывает эвтектическую кристаллизацию в малых объемах металла и имеет узкую область применения. Поэтому адекватное и точное описания процессов при затвердевании многокомпонентных сплавов остается актуальной задачей.
В настоящей работе рассмотрен процесс растекания и затвердевания на подложке капли двухкомпонентного сплава алюминия системы А1 + Б1, модифицированного активированными тугоплавкими наноразмерными частицами, которые являются центрами зарождения и роста кристаллической фазы. Сформулирована математическая модель, описывающая термо- и гидродинамические явления в капле при ее соударении с твердой подложкой и кинетику кристаллизации в процессе охлаждения расплава. Использовано полученное по результатам экспериментов значение термического сопротивления на поверхности контакта расплав-подложка [12]. Гидродинамические процессы рассмотрены в условиях затвердевания ме-
талла. Определены кинетика роста твердой фазы и температурные режимы в затвердевающем расплаве при различных размерах формирующихся сплэтов. Проведено сравнение формы сплэтов, полученных по результатам численного моделирования, с имеющимися данными [12], а также подтверждена адекватность предложенной модели кристаллизации при наличии ультрадисперсных тугоплавких частиц в сплаве алюминия имеющимися в литературе результатами экспериментов [10].
Математическая модель и алгоритм реализации. Предполагается, что сферическая частица расплава диаметром йо соударяется с подложкой под прямым углом к ее поверхности со скоростью Уо и растекается по поверхности (рис. 1). Расплав модифицирован тугоплавкими нанораз-мерными частицами, которые служат активными центрами кристаллизации и их массовое содержание тр не превышает 0,05 %. Диаметр частиц йр много меньше размера капли, поэтому они не влияют на теплофизиче-ские характеристики расплава.
Рис. 1. Схема деформирующейся капли после соударения с подложкой: 1 — капля; 2 — свободная поверхность расплава; 3 — подложка
Для численного моделирования поведения расплава после его соударения с поверхностью подложки используются уравнения Навье — Стокса для вязкой несжимаемой жидкости. Значения теплофизических характеристик материала капли в жидком и твердом состояниях полагаются постоянными и равными средним значениям в рассматриваемых интервалах значений температуры. Поверхность подложки плоская и абсолютно смачиваемая. Динамический угол контакта в перемещающейся точке соприкосновения поверхности капли и твердой поверхности определяется из решения задачи. Температура подложки ниже температуры
кристаллизации металла. Радиационным теплообменом на поверхности подложки пренебрегаем.
С учетом принятых допущений уравнения движения жидкости и неразрывности в цилиндрической системе координат (г, z), согласно схеме, приведенной на рис. 1, имеют вид
du du du 1 dp 1
— + u — + v— =---+ —
dt dr d z p dr p
dv dv dv 1 dp 1 — + u — + v— =---- + —
dt dr dz p dz p
1 d ir du"j^öf du^ u r dr ^ dr ) dz v dz ) r2
1 d ( dv ^ d ( dv --1 r^— 1 + — I д —
r dr v dr) dz v dz
g;
(1)
1 dru dv
= 0.
г дг дг
Здесь и, V — компоненты вектора скорости вдоль координат г и г; t — время; р — давление; д — динамическая вязкость, значение которой для жидкой и дисперсной сред определяется, согласно [14], по формуле
Ц0
при
fs - 0;
Цю exp [5 fs /3 (1 - fs)] при 0 < fs < 0,4,
где fs — доля твердой фазы в расплаве при кристаллизации (если fs > 0,4, то и = 0, V = 0).
Начальные условия ^ = 0) для уравнений (1) и = 0, V = р = = 4а / й0 + ра, где а — поверхностное натяжение расплавленного металла; ра — атмосферное давление.
Граничные условия для гидродинамической задачи при г = 0 опреде-
лены как u = 0, r жидкой частицы —
3v
— | = 0, r dr
dp
— \ = 0, а на свободной поверхности dr
„ „ du ^ f du dv . p-2ок-2^—1nr—+— Inz =0;
„ ^ „ dv Л (du dv .
p - 2<JK - 2^— I nz —+ — I nr = 0.
dz) \dz dr
Здесь пг, пг — компоненты единичного вектора нормали п к свободной поверхности капли по осям г и г; К — кривизна поверхности капли,
K _ r2(r'zff-z'r") + [(r')2 + (z')2]rz' ~ 2r2 [(r ')2 + (z ')2]3/2 '
где штрих обозначает производную относительно длины дуги s (см. рис. 1).
До начала затвердевания материала капли на поверхности подложки (z = 0) и далее на границе расплав-твердая фаза выполняются условия прилипания u = 0, v = 0.
Уравнение, описывающее перенос теплоты в капле:
dT dT dT 1 б ( dT Л б ( 6T Л к dfs — + u — + v — =--1 rae— 1 + — I ae— 1+--—, (2)
dt dr dz r dr v dr J dz v dz J c2 dt
где ae — локальное значение коэффициента температуропроводности, ae = ai при fs = 0, ae = a2 приf = 1, ae = ai(1 -fs) + a2fs в случае 0 < f < 1 и u = 0, v = 0 при fs > 0,4. Здесь ai = Xi /(ci pi) — температуропроводность жидкой (i = 1) и твердой (i = 2) фаз материала капли, материала подложки (i = 3); X, c, р — теплопроводность, теплоемкость и плотность; к — удельная теплота плавления.
Перенос теплоты в подложке описывается уравнением
dT
— = a3 dt
(1 d dT d2T)
r dr dr dr2
(3)
Начальные значения (^ = 0) температуры и доли твердой фазы в расплаве Т = Гю и = 0, температура подложки Т = Т30.
Граничное условие на поверхности соприкосновения капли и подложки г = 0, г < Гс (^ согласно [12]:
Х- — 1 dz
. dT
= A3 —
0+ dz
= R"(Tl0+-Tl0-), i = 1 2'
0- Rres
где Яе — тепловое сопротивление; Гс (^ — радиус пятна контакта (см. рис. 1).
Теплообмен между поверхностью капли и окружающей средой
än
h — = sa0(T340 - T4),
где в — степень черноты материала; а0 — постоянная Стефана — Больц-мана.
Условия симметрии в капле и подложке (г = 0): г^Т / dг) = 0.
Граничные условия на нижней (г = -hs, 0 < г < гс), верхней несмочен-ной (г = 0, г > гс) и боковой (г > гз, ^ < г < 0) поверхностях подложки имеют вид
*3 ^ = 0.
дг
Здесь ^, г3 — толщина и радиус подложки (2г3 » й0).
Предполагается, что рост кристаллической фазы подчиняется нормальному механизму и определяется линейной зависимостью скорости роста от переохлаждения Я = KvАТ [4, 15], где Я — радиус кристалла, точка означает производную по времени; Kv — физическая константа; АТ = Т/ - Т, Т/ — температура ликвидуса.
Полагая, что рост твердой фазы происходит одновременно на всех наноразмерных частицах в области переохлажденного расплава начиная со времени t = ^0, когда температура достигла исходной температуры ликвидуса Т/о, согласно [16], определяем fs соотношением
/5 = 1 - ехр {33 кЫрЯ31, (4)
t
где Я(г, г, t) = гр + | ^АТйт, гр = йр/2; Ыр — число наночастиц в единице объема. Переохлаждение АТ = Т/ - Т имеет вид
АТ = ТА--- Т.
(1 - и )1_к
Здесь температура ликвидуса Т/ аппроксимирована линейной зависимостью от концентрации С растворенного компонента (81); Та — температура плавления чистого металла-растворителя (А1); Р — модуль коэффициента наклона линии ликвидуса на диаграмме состояния А1 + 81. Концентрацию легирующего компонента определяем из уравнения неравновесного рычага [15] С = С0 / (1 - fs)1-k, где С0 — исходная концентрация; к — коэффициент распределения растворенного компонента (81). Рост твердой фазы а-компоненты (алюминия) сплава А1 + 81 происходит в температурном интервале Т/0 > Т > Те, где Т/0 = Та - РС0, Те — температура эвтектики. Величина дfs / дt определяется как
и = 4пЫр^(1 - и,)Я2 (Т/ - Т).
При расчетах кристаллизации в уравнении (2) для к используется значение удельной теплоты плавления а-компоненты (алюминия) сплава. Полагаем, что при T = Te доля твердой фазы равна fsa.
После охлаждения металла до температуры эвтектики происходит затвердевание Р-компоненты сплава. Ввиду малой взаимной растворимости алюминия и кремния полагаем, что при дальнейшем охлаждении расплава рост твердой фазы подчиняется нормальному механизму, характеризуемому константой роста Ke. Доля твердой фазы fse в затвердевающем эвтектическом расплаве определяется соотношением
fse = 1 - exp j3 nNp (R3 - яЦ, (5)
t
где R(r, z, t) = Ra + J Ke (TE - T)dт — радиус зародыша, растущего после
te
момента времени t = tE, когда температура расплава достигла температу-
tE
ры TE; Ra (r, z, tE) = rp + J KvATdx. При расчетах кристаллизации эвтек-
tio
тики в уравнении (2) вместо к используется k$i — удельная теплота плавления для Si и формула
% = 4KNpKe(1 - fse)R2(Te - T). ot
Эвтектическая кристаллизация проходит в диапазоне Te > T > Tend, где Tend — температура полного затвердевания расплава. Доля твердой фазы fs в процессе полного затвердевания сплава определяется как fs = fs а + fse.
Алгоритм решения основан на использовании конечно-разностной аппроксимации уравнений Навье — Стокса и теплопереноса. В области перемещения капли строится пространственная сетка: n = ihr, i = 0, ..., I, Zj = jhz, j = 0, ..., /. Стороны ячеек обозначаются целыми значениями индексов i, j, а центры — полуцелыми значениями i - 1/2 и j - 1/2. Область неустановившихся течений расплава ограничивается системой последовательно пронумерованных частиц-маркеров, расположенных по контуру капли, которые перемещаются с локальной скоростью жидкости. Определено шесть вариантов пересечения ячейки движущейся границей жидкости, различное сочетание которых используется при расчетах. В подложке применяется сетка вида ri = ihr, i = 0, ..., I, Zj = -jhzp, j = 0, ..., Jp. Распределение температуры описывается значениями в узлах сетки,
составляющие скоростей и, V определяются в серединах боковых граней ячеек, а давление р — в центрах ячеек [18-20].
При аппроксимации уравнений (1) используется неявная конечно-разностная схема. Аппроксимация вязкостных членов осуществляется со вторым порядком точности. Для выполнения требования монотонности схемы пространственные производные в левой части уравнений аппроксимируются с учетом знака скорости в окрестности узла сетки. Основное критическое замечание может сводиться к тому, что погрешность аппроксимации имеет порядок 0(т, h), однако окончательная проверка численного метода сводится к сравнению полученного результата с данными соответствующего физического эксперимента. Обоснованность применения таких аппроксимаций обсуждалась многими авторами (см. [21] и имеющуюся там библиографию).
Давление и поле скоростей, соответствующие моменту времени ^, определяются из поэтапного решения системы (1). Первоначально, с использованием значений скоростей с предшествующего временного шага в уравнениях количества движения находились составляющие поля скоростей. Далее, с использованием метода искусственной сжимаемости [20] рассчитывалось давление. Проводилось несколько итераций для согласования найденных скоростей и давления внутри капли и в граничных ячейках до выполнения неравенства тах |У-и| < у, где и — вектор скорости в расплавленном металле с компонентами и, V в радиальном и осевом направлениях (в твердом металле и = 0, V = 0); у — заданное малое число (у = 10-4). Решение алгебраических систем, получаемых при аппроксимации уравнений движения, осуществляется итерационным методом сопряженных градиентов.
После определения скоростей вычисляется температурное поле в двух сопряженных областях в капле и подложке по уравнениям (2)-(5). Для аппроксимации уравнений используется неявная конечно-разностная схема. Пространственные производные в левой части уравнения (2) в области расплавленного металла аппроксимируются с учетом знака скорости в окрестности каждого узла сетки для выполнения требования о монотонности схемы, имеющей порядок О(т, К). В области затвердевшего расплава и подложке исходная задача аппроксимируется с порядком О(т, к2). Решение алгебраических систем, получаемых при аппроксимации уравнений теплопереноса, осуществляется итерационным методом. Для согласования температурных полей в расчетных областях применяется метод простой итерации.
Соотношения
Гт = Гт + ит^н+1; гт = гт + чт^п+1; т = 1, •••, М(^),
используются для вычисления новых координат маркеров гц , на
(п+1)-м временном шаге. Здесь ит, чЦ — скорости перемещения частиц, определенные в соответствии со скоростями жидкости в ближайшей окрестности каждой из них; М(^) — число маркеров в момент времени Ъ. На каждом временном шаге проводится удаление или добавление частиц для их равномерного расположения по контуру капли и вводится новая нумерация (М(и) ~ 2500).
Временной шаг тп+1 для проведения расчетов в момент определяется по соотношению
тт{ЫГ, Нг}
^П+1 — I I I >
К0 тах{иЦ,}
т
К > 1 (при расчетах К0 = 20).
Адекватность модели и алгоритма ее реализации подтверждена качественным и количественным совпадением результатов расчетов с данными физических экспериментов [10, 12, 17].
Результаты численных экспериментов. Численное исследование динамики растекания и кристаллизации модифицированной металлической капли на подложке проводилось при следующих параметрах. Исходные данные для капли сплава А1 + 6,5 % [10, 12, 22]: = 3,92-103, 6-105 м; у0 = 3, 50 м/с; |Д0 = 0,00109 кг/(м-с); а = 1,07 Н/м; С1 = 1000 Дж/(кг-К);
= 70 Вт/(м-К); рх = 2,57-103 кг/м3; к = 3,89-105 Дж/кг; С2 = 1000 Дж/(кг-К); ^2 = 152,5 Вт/(м-К); р2 = 2,57-103 кг/м3; кя = 1,64-105 Дж/кг; Т0 = 903 К, Та = 933 К, Т/0 = 886,85 К; Те = 850 К; р = 7,1 К/%; СЪ = 6,5 % масс., к = 0,14, К = 0,005 м/(с-К); Ке = 0,005 м/(с-К); в = 0,4; ра = 105 Н/м2; = = 5,7-10-8 Вт/(м2-К4). Размер модифицирующих наночастиц йр = 10-7 м, плотность рр = 5440 кг/м3, концентрация тр = 0,03 % масс. Для подложки (сталь) [12]: Ы = 2-10-3 м, Г3 = 1,5-10-2 м, 1,540"3 м, Т30 = 473 К, С3 = = 550 Дж/(кг-К), р3 = 7,8-103 кг/м3; ^3 = 25,1 Вт/(м-К), = 1,64-10-6 (м2-К)/Вт. Параметры сетки: I = 440, ] = 200, ¡р = 100.
Для проверки модели кристаллизации рассмотрено объемное затвердевание алюминиевого сплава А1 + 6,5 % в цилиндрической форме с инокулированными тугоплавкими наноразмерными частицами Т1Ы
в количестве 0,03 % от массы (объемная концентрация наночастиц Ыр = = 2,7-1017 1/м3). При расчетах для определения начальных и граничных условий использовались условия эксперимента, в котором фиксировалось изменение температуры при остывании и кристаллизации в центре цилиндрической отливки диаметром 34 мм и высотой 51 мм [10]. Результаты, полученные при численном моделировании и в процессе замеров на опытной установке, приведены на рис. 2. Результаты расчетов хорошо коррелируют с данными физического эксперимента и свидетельствуют о том, что предложенная модель удовлетворительно описывает реальный физический процесс кристаллизации модифицированного сплава при принятых значениях кинетических констант роста.
Г, К
900
800
0 100 200 и с
Рис. 2. Изменение температуры при охлаждении и кристаллизации
сплава А1 + 6,5 % 81: точки — эксперимент; сплошная линия — расчет
Полученная по результатам расчетов динамика деформации капли диаметром 3,92-10 м при v0 = 3 м/с приведена на рис. 3. Расплав модифицирован наноразмерными частицами. Рассчитанные контуры меняющейся свободной поверхности капли для различных моментов времени сравниваются с результатами, представленными в работе [12] и полученными при формировании сплэта в сходных условиях. Определено, что при растекании капли толщина сформировавшегося сплэта в центральной части составляет 0,23 мм, а кристаллизация металла происходит в течение 6,45 мс. После того как расплав, примыкающий к подложке, затвердевает, расходящееся течение жидкости в верхнем слое образует валик на периферии формирующегося сплэта. Движение расплава происходит до тех пор, пока растущая твердая фаза не образует кристаллический каркас.
в
Рис. 3. Динамика поверхности капли после соприкосновения с подложкой при 0,3 (а), 3 (б) и 7 мс (б)
Изменение температуры, переохлаждения расплава, доли твердой фазы в металле за первые 15 мкс в точке соприкосновения с подложкой (r = 0) показано на рис. 4. Переохлаждение определяется как AT = Ti - T при Tio > T > Te и AT = TE - T в случае Te > T > Tend. После контакта расплава с подложкой перегрев снимается в течение 0,15 мкс и происходит переохлаждение металла ниже температуры ликвидуса Ti более чем
на 4 К (рис. 4, а, б). Время метастабильного состояния, когда доля кристаллической фазы практически равна нулю, составляет около 0,1 мкс. Кристаллизация а-компоненты (А1) сплава А1 + Б1 вблизи подложки происходит в течение 4 мкс, достигает значения /¡а = 0,494, и после небольшого переохлаждения расплава (« 0,1 К) начинается эвтектическая кристаллизация Р-компоненты (Б1). Полное затвердевание завершается к моменту времени 11 мкс (см. рис. 4).
t, мкс
Рис. 4. Изменение температуры (а), переохлаждения (б) и доли твердой
фазы (б) в расплаве в точке соприкосновения капли с подложкой
0 5 10 и мкс
в
Изменение доли твердой фазы и температуры в слое металла толщиной 0,23 мм вблизи подложки при г = 0 в различные моменты времени показано на рис. 5. Согласно полученным результатам, за период времени 15 мкс у поверхности соприкосновения с подложкой образуется двухфазная зона толщиной 0,02 мм (рис. 5, а). Перегрев в расплавленном металле на расстоянии 0,1 мм от подложки сохраняется, так как происходит приток теплоты в ходе движения расплава (рис. 5, б). К моменту времени 0,5 мс в расплаве у поверхности подложки образуется слой затвердевшего металла толщиной 30 мкм, а ширина двухфазной зоны составляет 70 мкм (рис. 5, б).
Изменение температуры и доли твердой фазы в металле до полного затвердевания расплава в точке соприкосновения с подложкой г = 0 мм и в плоскости г = 0,23 мм при г = 0 представлено на рис. 6. По мере того как
Численное моделирование кристаллизации модифицированной металлической капли. fs
0,5
0,2 z, мм
0,2 z, мм
0,2 z, мм
0,2 z, мм
Рис. 5. Распределение доли твердой фазы (а, б) и температуры (б, г) в слое металла вблизи плоскости соприкосновения с подложкой в течение 15 мкс (а, б)
и 0,5 мс (б, г)
Г, К 900
850
800
fs
=\_^ 7/0 1
^^ ТЕ
0,5 2 ^--
/
^—-----
0
6 t, мс 0
6 t, мс
Рис. 6. Изменение температуры (а) и доли твердой фазы (б) в расплаве у подложки (1) и в плоскости свободной поверхности (2) формирующегося сплэта
движение жидкости в вертикальном направлении замедляется, а слои металла становится тоньше, к моменту времени 1 мс в точке г = 0,23 мм перегрев снимается. Рост твердой фазы начинается при переохлаждении металла ниже температуры ликвидуса не более чем на 0,01 К, время мета-стабильного состояния кристаллической фазы составляет около 0,4 мс
(см. рис. 6), кристаллизация а-компоненты сплава длится около 4,3 мс, а эвтектическая кристаллизация Р-компонены — 0,75 мс (рис. 6, б). С момента снятия перегрева в расплаве кристаллизация происходит во всем объеме незатвердевшего металла.
Переохлаждение, темп кристаллизации и время затвердевания существенно меняются с удалением от подложки (см. рис. 4, 6). Особенностью кинетики кристаллизации является наличие различных значений переохлаждений на температурных кривых. При этом температура рекалесценции (температура, при которой происходит локальное выравнивание интенсивностей выделения скрытой теплоты кристаллизации и ее отвод в затвердевшую корку металла и подложку) растет в направлении к свободной поверхности. В результате этого возникает осцилляция температуры и возможен даже нагрев нижних затвердевших слоев капли (рис. 6, а).
На форму сформировавшегося сплэта в основном влияют скорость падения капли на поверхность и ее размер. От температуры подложки существенно зависят переохлаждение, темп кристаллизации и время затвердевания расплава [11]. Полученная при расчетах динамика деформации капли диаметром ¿0 = 6-105 м при У0 = 50 м/с и совпадении остальных условий с рассмотренными выше показана на рис. 7. Растекание металла происходит в течение 2 мкс, а толщина сформировавшегося сплэта в центральной части составляет 4,5 мкм.
z/d0
Рис. 7. Форма капли после соприкосновения с подложкой при 0,5 (1),
1,5 (2), 2,0 мкс (3)
Изменение температуры, переохлаждения расплава и доли твердой фазы в затвердевающем металле на поверхности соприкосновения с подложкой и в плоскости г = 4,5 мкм формирующегося сплэта при г = 0 показано на рис. 8. У поверхности подложки перегрев снимается в течение 0,04 мкс, время метастабильного состояния около 0,07 мкс, кристаллизация а-компоненты сплава происходит за 2,8 мкс, а полное затвердевание — за 3,7 мкс, максимальное переохлаждение 7,5 К, переохлаждение для начала эвтектической кристаллизации 0,15 К. В центральной точке на поверхности
Г, К
Рис. 8. Изменение температуры (а), переохлаждения (б) и доли твердой фазы (б) в расплаве в плоскостях пятна контакта (1) и свободной поверхности (2) формирующегося сплэта
t, мкс
формирующегося сплэта продолжительность метастабильного состояния составляет 0,42 мкс, кристаллизация а-компоненты происходит за 5 мкс, а эвтектическая кристаллизация Р-компоненты — за 2,25 мкс, максимальное переохлаждение 1,2 К.
Во время деформации капли температура расплава на ее свободной поверхности фактически не меняется, а движение жидкости способствует переносу теплоты в направлении подложки, в которую оно отводится. Когда слой металла на подложке становится тонким, скорость в вертикальном направлении снижается, перегрев в расплаве снимается и с 1,5 мкс кристаллизация происходит во всем объеме металла. К моменту времени, когда движение жидкости прекращается (2 мкс), во всем ее объеме существует твердая фаза с долей, не превышающей 0,4 (рис. 9).
Особенность кинетики кристаллизации — наличие двух максимумов переохлаждений на температурных кривых, которые могут принимать значения от нескольких градусов до десятых долей градуса. Особенности изменения переохлаждения расплава и существование второго максимума при эвтектической кристаллизации, который особенно заметен в области интенсивного отвода теплоты вблизи поверхности соприкосновения с подложкой, продемонстрированы на рис. 8, б. Более высокие значения скорости охлаждения способствуют большим значениям переохлаждений.
0 2 t, мкс 0 2 At, мкс
а б
Рис. 9. Изменения температуры (а) и доли твердой фазы (б) в металле растекшейся капли от плоскости соприкосновения с подложкой
до свободной поверхности
Ввиду того, что переохлаждения имеют относительно небольшие значения, на основании данных работы [11] можно сделать следующий вывод: гомогенного зарождения кристаллической фазы во время формирования сплэта не происходит.
Заключение. Предложена математическая модель, описывающая процессы деформации и неравновесной кристаллизации жидкой капли бинарного сплава с модифицирующими тугоплавкими наноразмерными частицами при ее соударении с подложкой. Рассчитанные размеры затвердевшего на подложке сплэта и температурные режимы кристаллизации сплава удовлетворительно согласуются с известными экспериментальными данными.
Проведено численное моделирование кинетики затвердевания жидкой капли бинарного сплава (Al + Si), модифицированного наноразмерными частицами, после ее соударения с твердой подложкой. Описаны особенности кристаллизации расплава в зависимости от размеров капли и, как следствие, режимов теплообмена с подложкой. Установлено, что при большой скорости соударения капли с подложкой за время растекания расплава возможно образование лишь небольшой доли твердой фазы. Определено, что независимо от размера капли имеет место объемно-последовательная кристаллизация металла.
ЛИТЕРАТУРА
[1] Коренев В.Н., Родичев А.Ю., Семенов А.В. Методы восстановления и упрочнения деталей газопламенным напылением. Сварочное производство, 2014, № 2, с. 40-43.
[2] Solonenko O.P., Kudinov V.V., Smirnov A.V., et al. Micro-metallurgy of splats: theory, computer simulation and experiment. JSME Int. J. Ser. B, 2005, vol. 48, iss. 3, pp. 366-380. DOI: https://doi.org/10.1299/jsmeb.48.366
[3] Ильиных С.А., Кирнос И.В., Крашанинин В.А. и др. Физико-химические свойства покрытий, получаемых дозвуковым и сверхзвуковым плазменным напылением порошков металлов и их композиций. Известия вузов. Порошковая металлургия и функциональные покрытия, 2015, № 1, с. 49-54.
DOI: https://doi.org/10.17073/1997-308X-2015-1-49-54
[4] Фомин В.М., Черепанов А.Н., ред. Плазмохимический синтез ультрадисперсных порошков и их применение для модифицирования металлов и сплавов. Новосибирск, Наука, 1995.
[5] Cherepanov A.N., Poluboyarov V.A., Solonenko O.P., et а1 Impact of mechanical activation and modification of initial powder with refractory nanoparticles on the properties of plasma-sprayed coatings. In: Progress in plasma processing of materials 2003. Begell House Inc., 2003, pp. 507-513.
[6] Borodianskiy K., Zinigrad M., Gedanken A. Aluminum A356 reinforcement by carbide nanoparticles. JNanoR, 2011, vol. 13, pp. 41-46.
DOI: https://doi.org/10.4028/www.scientific.net/JNanoR.13.41
[7] Kuzmanov P.M., Popov S.I., Yovkov L.V., et al. Investigation the effect of modification with nanopowders on crystallization process and microstructure of some alloys. AIP Conf. Proc., 2017, vol. 1893, iss. 1, art. 030104. DOI: https://doi.org/10.1063/L5007562
[8] Bozhanova N.M., Panov I.T., Manolov V.K., et al. Modification of properties of aluminum protective a nodes by nanopowder materials. Thermophys. Aeromech., 2018, vol. 25, iss. 5, pp. 759-764. DOI: https://doi.org/10.1134/S0869864318050128
[9] Lazarova R., Bojanova N., Dimitrova R., et al. Influence of nanoparticles introducing in the melt of aluminum alloys on castings microstructure and properties. Inter. Metal-cast., 2016, vol. 10, iss. 4, pp. 466-476. DOI: https://doi.org/10.1007/s40962-016-0033-7
[10] Popov S., Manolov V., Cherepanov A. Mathematical modelling of crystallization of metal alloys. Comptes rendus de l'Academie Bulgare des Science, 2012, vol. 65, no. 4, pp. 441-448.
[11] Cherepanov A.N., Popov V.N., Solonenko O.P. Numerical analysis of the spreading and crystallization dynamics of the modified metal droplet on the substrate. Thermophys. Aeromech., 2008, vol. 15, iss. 3, pp. 483-489.
DOI: https://doi.org/10.1134/S0869864308030141
[12] Xue M., Heichal Y., Chandra S., et al. Modeling the impact of a molten metal droplet on a solid surface using variable interfacial thermal contact resistance. Mater. Sci., 2007, vol. 42, iss. 1, pp. 9-18. DOI: https://doi.org/10.1007/s10853-006-1129-x
[13] Popov S., Manolov V., Kuzmanov P., et al. Mathematical model of crystallization of multicomponent alloy at presence of nanoparticles. J. Mater. Sci. Technol., 2014, vol. 22, no. 3, pp. 167-174.
[14] Мастрюков Б.С. Теплофизика металлургических процессов. М., МИСиС, 1996.
[15] Flemings M.C. Solidification processing. McGraw-Hill, 1974.
[16] Баландин Г.Ф. Основы теории формирования отливки. Ч. 2. М., Машиностроение, 1979.
[17] Попов В.Н. Моделирование затвердевания металлической капли на холодной подложке. Математическое моделирование, 2001, т. 13, № 9, с. 119-127.
[18] Harlow F.H., Welch J.E. Numerical calculation of time-depend viscous incompressible flow of fluid with free surface. Phys. Fluids., 1965, vol. 8, iss. 12, pp. 2182-2189. DOI: https://doi.org/10.1063/L1761178
[19] Patankar S.V., Spalding D.B. A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows. Int. J. Heat Mass Trans., 1972, vol. 15, iss. 10, pp. 1787-1806. DOI: https://doi.org/10.1016/0017-9310(72)90054-3
[20] Chorin A.J. A numerical method for solving incompressible viscous flow problems. J. Comput. Phys. 1967, vol. 137, iss. 2, pp. 118-125.
DOI: https://doi.org/10.1006/jcph.1997.5716
[21] Roache P. Computational fluid dynamics. Hermosa Publ., 1976.
[22] Зиновьев В.Е. Теплофизические свойства металлов при высоких температурах. М., Металлургия, 1989.
Попов Владимир Николаевич — д-р физ.-мат. наук, старший научный сотрудник, главный научный сотрудник ИТПМ СО РАН (Российская Федерация, 630090, Новосибирск, ул. Институтская, д. 4/1).
Черепанов Анатолий Николаевич — д-р физ.-мат. наук, профессор, главный научный сотрудник ИТПМ СО РАН (Российская Федерация, 630090, Новосибирск, ул. Институтская, д. 4/1).
Просьба ссылаться на эту статью следующим образом:
Попов В.Н., Черепанов А.Н. Численное моделирование кристаллизации модифицированной металлической капли при растекании на подложке. Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки, 2019, № 6, с. 18-39. DOI: 10.18698/1812-3368-2019-6-18-39
NUMERICAL SIMULATION OF CRYSTALLIZATION OF A MODIFIED METAL DROP DURING METAL SPREADING ON A SUBSTRATE
Siberian branch of Russian Academy of Sciences, Khristianovich Institute of Theoretical and Applied Mechanics, Novosibirsk, Russian Federation
The purpose of the research was to numerically simulate Numerical simulation, metal
the processes when melting drops fall on a substrate. drop, substrate, modification,
The paper deals with the solidification on the metal aluminum alloy, nanosized surface of a binary aluminum alloy modified by refractory particles, crystallization
V.N. Popov A.N. Cherepanov
popov@itam.nsc.ru ancher@itam.nsc.ru
Abstract
Keywords
activated refractory nanosized particles, which are the
centers of crystalline phase nucleation. We formulated a
mathematical model which describes the thermo- and
hydrodynamic phenomena in the drop upon interaction
with a solid substrate, heterogeneous nucleation during
melt cooling, and subsequent crystallization. The flow in
a liquid is described by the Navier — Stokes equations in
the Boussinesq approximation. The position of the free
boundary of the melt is fixed by marker particles moving
with the local liquid velocity. On the melt — substrate
contact surface, thermal resistance is taken into account.
The hydrodynamic problem is considered under
conditions of crystallization of molten metal. The
temperature conditions and the kinetics of the growth of
the solid phase in the solidifying aluminum alloy are
described for various sizes of formed splats. Satisfactory
agreement was found between the shape of the splat
obtained by the results of numerical simulation and the
available experimental data. The adequacy of the
crystallization model in the presence of ultradisperse
refractory particles in a binary alloy is confirmed. It was
determined that, regardless of the size of the drop, bulk
crystallization of the metal takes place. It was found that
at a high rate of collision of a drop with a substrate
during the melt spreading, a small fraction of the solid Received 19.02.2019
phase can be formed © Author(s), 2019
The work was carried out as part of the Program of Fundamental Scientific Research of State Academies of Sciences for 2013-2020 (project no. 0323-2018-0012)
REFERENCES
[1] Korenev V.N., Rodichev A.Yu., Semenov A.V. Rework and strengthening techniques of parts by gas-flame sputtering. Svarochnoe proizvodstvo, 2014, no. 2, pp. 40-43 (in Russ.).
[2] Solonenko O.P., Kudinov V.V., Smirnov A.V., et al. Micro-metallurgy of splats: theory, computer simulation and experiment. JSME Int. J. Ser. B, 2005, vol. 48, iss. 3, pp. 366-380. DOI: https://doi.org/10.1299/jsmeb.48.366
[3] Il'inykh S.A., Kirnos I.V., Krashaninin V.A., et al. Physicochemical properties of coatings fabricated by subsonic and supersonic plasma deposition of powders of metals and their compositions. Izvestiya vuzov. Poroshkovaya metallurgiya i funktsional'nye pokrytiya [Universities' Proceedings. Powder Metallurgy and Functional Coatings], 2015, no. 1, pp. 49-54 (in Russ.). DOI: https://doi.org/10.17073/1997-308X-2015-1-49-54
[4] Fomin V.M., Cherepanov A.N., eds. Plazmokhimicheskiy sintez ultradispersnykh po-roshkov i ikh primenenie dlya modifitsirovaniya metallov i splavov [Plasma chemical synthesis of ultradisperse powders and their application for modified metals and alloys]. Novosibirsk, Nauka Publ., 1995.
[5] Cherepanov A.N., Poluboyarov V.A., Solonenko O.P., et al. Impact of mechanical activation and modification of initial powder with refractory nanoparticles on the properties of plasma-sprayed coatings. In: Progress in plasma processing of materials 2003. Begell House Inc., 2003, pp. 507-513.
[6] Borodianskiy K., Zinigrad M., Gedanken A. Aluminum A356 reinforcement by carbide nanoparticles. JNanoR, 2011, vol. 13, pp. 41-46.
DOI: https://doi.org/10.4028/www.scientific.net/JNanoR.13.41
[7] Kuzmanov P.M., Popov S.I., Yovkov L.V., et al. Investigation the effect of modification with nanopowders on crystallization process and microstructure of some alloys. AIP Conf. Proc., 2017, vol. 1893, iss. 1, art. 030104. DOI: https://doi.org/10.1063/L5007562
[8] Bozhanova N.M., Panov I.T., Manolov V.K., et al. Modification of properties of aluminum protective a nodes by nanopowder materials. Thermophys. Aeromech., 2018, vol. 25, iss. 5, pp. 759-764. DOI: https://doi.org/10.1134/S0869864318050128
[9] Lazarova R., Bojanova N., Dimitrova R., et al. Influence of nanoparticles introducing in the melt of aluminum alloys on castings microstructure and properties. Inter. Metal-cast., 2016, vol. 10, iss. 4, pp. 466-476. DOI: https://doi.org/10.1007/s40962-016-0033-7
[10] Popov S., Manolov V., Cherepanov A. Mathematical modelling of crystallization of metal alloys. Comptes rendus de l'Academie Bulgare des Science, 2012, vol. 65, no. 4, pp. 441-448.
[11] Cherepanov A.N., Popov V.N., Solonenko O.P. Numerical analysis of the spreading and crystallization dynamics of the modified metal droplet on the substrate. Thermophys. Aeromech., 2008, vol. 15, iss. 3, pp. 483-489.
DOI: https://doi.org/10.1134/S0869864308030141
[12] Xue M., Heichal Y., Chandra S., et al. Modeling the impact of a molten metal droplet on a solid surface using variable interfacial thermal contact resistance. Mater. Sci., 2007, vol. 42, iss. 1, pp. 9-18. DOI: https://doi.org/10.1007/s10853-006-1129-x
[13] Popov S., Manolov V., Kuzmanov P., et al. Mathematical model of crystallization of multicomponent alloy at presence of nanoparticles. J. Mater. Sci. Technol., 2014, vol. 22, no. 3, pp. 167-174.
[14] Mastryukov B.S. Teplofizika metallurgicheskikh protsessov [Thermal physics of met-allurgic processes]. Moscow, MISIS Publ., 1996.
[15] Flemings M.C. Solidification processing. McGraw-Hill, 1974.
[16] Balandin G.F. Osnovy teorii formirovaniya otlivki. Ch. 2 [Foundation of cast forming theory. P. 2]. Moscow, Mashinostroenie Publ., 1979.
[17] Popov V.N. Modelling of the solidification of a metal drop on a cold substrate. Ma-tematicheskoe modelirovanie, 2001, vol. 13, no. 9, pp. 119-127 (in Russ.).
[18] Harlow F.H., Welch J.E. Numerical calculation of time-depend viscous incompressible flow of fluid with free surface. Phys. Fluids., 1965, vol. 8, iss. 12, pp. 2182-2189.
DOI: https://doi.org/10.1063/L1761178
[19] Patankar S.V., Spalding D.B. A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows. Int. J. Heat Mass Trans., 1972, vol. 15, iss. 10, pp. 1787-1806. DOI: https://doi.org/10.1016/0017-9310(72)90054-3
[20] Chorin A.J. A numerical method for solving incompressible viscous flow problems. J. Comput. Phys. 1967, vol. 137, iss. 2, pp. 118-125.
DOI: https://doi.org/10.1006/jcph.1997.5716
[21] Roache P. Computational fluid dynamics. Hermosa Publ., 1976.
[22] Zinovyev V.E. Teplofizicheskie svoystva metallov pri vysokikh temperaturakh [Thermal physical properties of metals at high temperatures]. Moscow, Metallurgiya Publ., 1989.
Popov V.N. — Dr. Sc. (Phys.-Math.), Senior Researcher, Chief Researcher, Siberian branch of Russian Academy of Sciences, Khristianovich Institute of Theoretical and Applied Mechanics (Institutskaya ul. 4/1, Novosibirsk, 630090 Russian Federation).
Cherepanov A.N. — Dr. Sc. (Phys.-Math.), Professor, Chief Researcher, Siberian branch of Russian Academy of Sciences, Khristianovich Institute of Theoretical and Applied Mechanics (Institutskaya ul. 4/1, Novosibirsk, 630090 Russian Federation).
Please cite this article in English as:
Popov V.N., Cherepanov A.N. Numerical simulation of crystallization of a modified metal drop during metal spreading on a substrate. Herald of the Bauman Moscow State Technical University, Series Natural Sciences, 2019, no. 6, pp. 18-39 (in Russ.). DOI: 10.18698/1812-3368-2019-6-18-39