CC BY
17
Особенности перераспределения электрического и теплового полей на границе раздела фаз при пропускании электрического тока через расплав
В. Н. Цуркин, А. В. Иванов
Институт импульсных процессов и технологий НАН Украины, пр. Богоявленский, 43-А, г. Николаев, 54018, Украина, e-mail: [email protected]
В двумерной постановке методами численного эксперимента проанализированы особенности перераспределения линий однородного электрического поля в объеме металлического расплава на границе раздела наноразмерных фаз, отличающихся проводимостью от матрицы. Показано, что градиенты напряженности электрического и теплового полей определяются формой фазы и соотношением ее проводимости и проводимости матрицы. Анализируемые поверхностные эффекты также определяются расположением сектора поверхности включения относительно линий тока.
Ключевые слова: расплав, обработка электрическим током, наноразмерная частица, кристаллизационная способность, энергетические факторы.
УДК 537.311.33
ВВЕДЕНИЕ
В последние годы в специализированных журналах появляется все больше информации об использовании кондукционной электротоковой обработки (КЭТО) металлических расплавов в технологиях литейного производства. При этом непосредственно через расплав пропускают электрический ток в различных его режимах (постоянный, переменный, импульсный), в отличие от индуцированного в расплаве тока при наложении магнитного поля (индукционная обработка). В этих случаях различаются подходы к постановке и решению электродинамических, электрофизических и электротехнических задач. Результаты использования КЭТО однозначно показывают, что такой способ является эффективным и относительно дешевым, простым в реализации и экологически чистым инструментом для обеспечения высоких показателей качества литого металла независимо от его химического состава, способа выплавки и разливки. При этом обработку расплава проводили в печи, в промежуточной емкости и в форме [1-6].
Авторами многих публикаций отмечаются (на феноменологическом уровне) такие эффекты КЭТО, как увеличение количества зародышей кристаллизации при гомогенном и гетерогенном механизмах зародышеобразования, активизация процессов переноса и растворимости, дегазация и гомогенизация расплава. Но, к сожалению, широкое использование этого метода в промышленности, как, например, электромагнитное перемешивание или магнитогидродинамическая обработка расплава, пока ждет своего часа.
По мнению авторов данной статьи, причины этого обусловлены прежде всего недостаточным уровнем представлений о физических принципах и функциональных возможностях КЭТО изменять структуру и кристаллизационную способность объекта обработки. Набор таких данных позволил бы предложить научно обоснованные рекомендации по проектированию соответствующих высокоэффективных, энергосберегающих технологических модулей, реализующих КЭТО на разных стадиях металлургического передела в системе «шихта-расплав-отливка». В свою очередь это дало бы возможность методу КЭТО активно конкурировать на рынке литейного производства при условии оптимального соотношения «цена-качество».
ХАРАКТЕРИСТИКА ПРОБЛЕМНОЙ СИТУАЦИИ
Ключевая фундаментальная проблема литейного производства, решаемая на всех стадиях металлургического передела в системе «шихта-расплав-отливка», - это повышение показателей качества литого металла или обеспечение требуемых, соответствующих современным запросам [7]. При этом возникают различные проблемные ситуации, определение которых позволяет ставить и решать частные задачи. И уже их совокупность в той или иной мере показывает пути решения вышеупомянутой проблемы. Что касается метода КЭТО, то его, как, впрочем, и другие методы, следует отнести к категории многофакторных, многофункциональных и многоцелевых задач. При этом понимание природы функцио-
© Цуркин В.Н., Иванов А.В., Электронная обработка материалов, 2017, 53(4), 60-66.
нальных возможностей КЭТО требует привлечения модельных представлений металлофизики, электрофизики, электродинамики, гидродинамики и, конечно же, термодинамики. Такая задача в полной постановке решается методами мультифизического моделирования, для которых суть важно выделить сильные и (или) слабые связи анализируемых процессов (strong coupling and (or) weak coupling) [8]. Только такой интегральный подход позволит спрогнозировать внутренние структурно-кинетические и фазовые изменения в расплаве. Но их первопричиной являются внешние структурно-формирующие параметры нагружения металлической системы. Именно их строгий учет позволит перейти к решению мультифизической задачи.
Практически во всех публикациях, рассматривающих пропускание электрического тока через расплав, в качестве такого параметра принимается или амплитудное значение тока, или его плотность. Но следует уточнить, что сферой влияния на расплав при КЭТО является не электрический ток, а электромагнитное поле, которое формируется в объекте обработки. В силу того что ток в расплав подводится электродами конечных размеров и определенной геометрии, с различными вариантами их расположения в расплаве, а также с учетом различных режимов тока, такое поле на макроуровне нельзя назвать однородным [1, 9, 10]. Но расплав не является системой, состоящей из однородной механической смеси атомов как на микро-, суб-микро-, мезо-, так и атомарном уровнях структуры за счет наличия в нем нерастворенных и(или) нерастворимых фаз разных размеров и различного химического состава, а также из-за микронеоднородного строения жидкометаллической системы, определяемой как совокупность относительно упорядоченных кластеров, разделяемых разупорядоченными областями. Формируемые в расплаве зародыши кристаллизации, как правило, также вносят свою лепту в степень неоднородности расплава на соответствующем уровне его структуры. Такие локальные образования при наличии разницы их проводимости с проводимостью матрицы также должны дополнительно локально искривлять потенциалы электрического поля макроуровня. Это, безусловно, может вносить существенный вклад в изменение термодинамической обстановки в расплаве при его кристаллизации и последующем твердении. В ситуации, свободной от какого-либо внешнего воздействия, то есть в условиях естественной кристаллизации, пользуясь понятиями свободной энергии Гиббса, плотности поверхностной энергии на границе раздела фаз, энергии активации процессов переноса, сорбции и переохлаждения, можно
спрогнозировать возможные механизмы и результаты кристаллизации [11-13]. Но при наложении электрического поля, характеристики которого учтены не только на макроуровне, но и на иных уровнях структуры, энергетические параметры такого процесса кристаллизации должны корректироваться.
Не претендуя на полное решение такой задачи, но и не снижая ее общности, определим цель работы: показать методами численного моделирования характер и параметры перераспределения электрического и теплового полей в жидкометаллической системе на частицах (фазах), интерпретируемых по характерному размеру и проводимости как зародыши кристаллизации.
ПОСТАНОВКА ЗАДАЧИ
Для определения особенностей распределения электрического и теплового полей на границе раздела фаз в жидкометаллическом проводнике достаточно в первой постановке ограничиться моделированием ситуации, когда модельная частица находится в сфере влияния электрического поля и имеет проводимость аекл, которая больше или меньше проводимости расплава враспл. Ограничимся размером частицы порядка 10-8 м, что соответствует характерному размеру зародыша при гомогенной кристаллизации, например до- и заэвтектического силумина, а также возможному размеру подложки при гетерогенной кристаллизации [11, 12]. Сопоставляя размер частицы с возможными градиентами плотности тока при разных его режимах [10], ограничимся в первом приближении рассмотрением ситуации, когда частица присутствует в однородном электрическом поле. Рассмотрим единичную частицу, исключая эффекты ее взаимодействия с другими частицами, что существенно усложнило бы вычисления. Достаточно информативным в рамках данной задачи будет ее рассмотрение в двумерной постановке.
Геометрическая модель такой задачи представлена на рис. 1. Здесь 1-4 - границы рассматриваемой области.
Используя приведенные выше допущения, чтобы определить распределение плотности тока, напряженность электрического поля и удельную мощность тепловыделения в режиме постоянного тока, необходимо решить несколько довольно простых уравнений [3, 10]. Для определения скалярного электрического потенциала (Ф):
0 ^Ф (а— 1 = 0,
dx ^ dx) dy ^ dy ^
(1)
где о - удельная электрическая проводимость материала, См/м.
Рис. 1. Геометрическая модель задачи электротокового воздействия на частицу, находящуюся в матрице расплава.
Для определения плотности тока (J):
J = -oV(p.
(2)
Для определения напряженности электрического поля (Е):
J = -E.
(3)
Тогда удельную мощность тепловыделения (q) определяем выражением:
' = J
(4)
где |J| = ^Jx2 + Jу2 - модуль вектора плотности
компоненты
тока; J =- ^ и Jy = - °Ф -
ox dy
вектора плотности тока соответственно.
Поскольку задача стационарная, то начальное условие будет иметь вид: J = J0 \t=0 = const, где
J0 - задаваемая плотность тока.
Граничными условиями для областей 1-4 (рис. 1) будут выражения:
- область 1: n ■ J=0;
- область 2: П ■ J = П ■ J0;
- область 3: П ■ J = П ■ J0;
- область 4: n ■ J = 0.
Уравнения (1-4) решали методом конечных элементов в системе Comsole с соответствующими граничными и начальными условиями.
Зададимся исходными данными для расчета: вариант 1 (проводящее включение): Lp = 1-10-7 м;
d,
= 1-10-'
м;
Jраспл
= 3,8-106 См/м;
= 9,5-10 См/м; вариант 2 (непроводящее включение): Lp = 1-10-7 м; detul = 1-10-'
м;
орасПл = 3,8-106 См/м; овкл = 4-10-4 См/м.
Значение J0 оценим исходя из выражения:
Л = К / 4, (5)
где силу тока (1А) примем равной 10 А.
РЕЗУЛЬТАТЫ РАСЧЕТА И ИХ АНАЛИЗ
Результаты решения уравнений (1-4) для круглой частицы представлены на рис. 2-4 для варианта 1 и на рис. 5-7 - для варианта 2.
Приведенные на рис. 2-13 данные свидетельствуют, что логично, о существенном влиянии формы и проводимости частицы на характер перераспределения электрического и теплового полей в области «частица-расплав». При наличии проводящей частицы наблюдаются ослабление электрического поля внутри ее объема и скачок потенциала на ее границах. Это обусловливает характерное резкое увеличение на 50% тепловой мощности на границах, ориентированных вдоль направления электрического тока, и такое же скачкообразное уменьшение на 60% на границах, перпендикулярных направлению токовых линий. В связи с тем что электрическое поле проникает в проводящую частицу, тепловая мощность (~ 25-30% для данного случая) выделяется и внутри ее объема пропорционально овкл/орасп. Практически не проникая в непроводящую частицу, электрическое поле обусловливает характерный скачок напряженности поля и соответственно удельной мощности тепловыделения на границах, перпендикулярных линиям тока, на 70% для круглой частицы и на 50% в области середины грани кубического включения, но практически уменьшает его до нуля на границах, расположенных вдоль направления протекания тока.
Следует также обратить внимание на пространственное распределение исследуемых характеристик в приграничной с частицей области. Как известно [12, 13], на граничной поверхности с точки зрения термодинамики не наблюдается скачкообразного перехода одной фазы в другую. Приповерхностная область по Гиббсу в математическом представлении имеет объем, равный нулю. По Гугенхайму он учитывается некой толщиной, в которой происходит «формирование» поверхностной энергии, величина этой толщины составляет несколько ангстрем. Эти подходы в том числе отличают математическое описание работы, идущей на образование поверхности раздела фаз. Как видно из представленных на рис. 2-13 графиков, размер переходной зоны, в которой происходит перераспределение анализируемых характеристик, составляет величину, равную примерно 2ёвкл, что на два порядка больше толщины переходного слоя, «формирующего» поверхностную энергию. Таким образом, электрическое поле при КЭТО изменяет термодинамическую обстановку вблизи частицы, а значит, и кристаллизационную способность расплава, внося определенные
(а) (б) (в)
Рис. 2. Распределение напряженности электрического поля, удельной тепловой мощности, плотности тока и электрического потенциала (круглая проводящая частица): (а) - линии электрического потенциала и тока; (б) - распределение напряженности электрического поля; (в) - распределение удельной мощности тепловыделения.
Е, -108 В/м
Е, -108 В/м
3,5 3
2,5
2 1,5
J
2,5
1,5
О 0,5 х, -10"7 м 1 0 0,5 х, -10' м 1
(а) (б)
Рис. 3. Числовое распределение напряженности электрического поля в срединных сечениях, проходящих через области 2-3 (а) и 1-4 (б) (круглая проводящая частица).
д, -1023 Вт/м
д, -1023 Вт/м
0,5 х, -10"' м 1
(а)
0,5 х, -10 м 1
(б)
Рис. 4. Числовое распределение удельной тепловой мощности в срединных сечениях, проходящих через области 2-3 (а), 1-4 (б) (круглая проводящая частица).
(а) (б) (в)
Рис. 5. Распределение напряженности электрического поля, удельной тепловой мощности, плотности тока и электрического потенциала (круглая непроводящая частица): (а) - линии электрического потенциала и тока; (б) - распределение напряженности электрического поля; (в) - распределение удельной мощности тепловыделения.
6 5 4
3 2 1 О
Е, -108 В/м
Е, -108 В/м
1Г
О
х, -10-' м
1
0.? х, -10-' м 1
0,5
(а) (б)
Рис. 6. Числовое распределение напряженности электрического поля в срединных сечениях, проходящих через области 2-3 (а), 1-4 (б) (круглая непроводящая частица).
д, -1023 Вт/м
д, -1023 Вт/м3
10 8 6 4 2 0
о
0,5
(б)
х, -10-7 м
1
Рис. 7. Числовое распределение удельной тепловой мощности в срединных сечениях, проходящих через области 2-3 (а), 1-4 (б) (круглая непроводящая частица).
1 2.
V
\ -
- ....
—
(а) (б) (в)
Рис. 8. Распределение напряженности электрического поля, удельной тепловой мощности, плотности тока и электрического потенциала (кубическая проводящая частица): (а) - линии электрического потенциала и тока; (б) - распределение напряженности электрического поля; (в) - распределение удельной мощности тепловыделения.
3,5 3 2,5 2 1,5 1
Е, -108 В/м
г\
о
0.? х, -10-7 м 1
2,6 2,4 2,2 2 1,8 1,6 1,4
Е, -108 В/м
О
0,5 х, -10-7 м 1
(а)
(б)
Рис. 9. Числовое распределение напряженности электрического поля в срединных сечениях, проходящих через области 2-3 (а), 1-4 (б) (кубическая проводящая частица).
д, -1023 Вт/м
д, -1023 Вт/м
2,8 2,6 2,4 2,2 2 1,8 1,6 1,4 1,2 1
0 5 х, -10" м ]
О
0,5
х, -10-7 м 1
(б)
Рис. 10. Числовое распределение удельной тепловой мощности в срединных сечениях, проходящих через области 2-3 (а), 1-4 (б) (кубическая проводящая частица).
: <
\
V
□
— - —^— _
(в)
(а) (б)
Рис. 11. Распределение напряженности электрического поля, удельной тепловой мощности, плотности тока и электрического потенциала (кубическая непроводящая частица): (а) - линии электрического потенциала и тока; (б) - распределение напряженности электрического поля; (в) - распределение удельной мощности тепловыделения.
Е, -108 В/м
м
о
0,5
(а)
х, -10-7 м
6 5,5 5 4,5 4 3,5 3 2,5 2
Е, -108 В/м
0
0,5 х, -10- м 1
(б)
Рис. 12. Числовое распределение напряженности электрического поля в срединных сечениях, проходящих через области 2-3 (а), 1-4 (б) (кубическая непроводящая частица).
д, -1023 Вт/м
6 5 4 3 2 1 0
д, -1023 Вт/м3
х, -10-7 м ]
0 0,5
(а) (б)
Рис. 13. Числовое распределение удельной тепловой мощности в срединных сечениях, проходящих через области 2-3 (а), 1-4 (б) (кубическая непроводящая частица).
коррективы в ее движущую силу, - химический потенциал, величина и знак которого в разных секторах поверхности частицы будут разными в зависимости от ориентации этого сектора относительно потенциальных линий тока.
ВЫВОДЫ
1. Методами численного моделирования условий кристаллизации расплава, подвергаемого воздействию электрического тока, определены характер пространственного распределения электрического поля и поля тепловой мощности на единичных частицах, интерпретируемых по характерному размеру и проводимости как зародыши кристаллизации. Установлено, что его параметры качественно и количественно зависят от формы частицы и соотношения ее проводимости с проводимостью расплава, а также от ориентации сектора поверхности частицы относительно линий тока.
2. Основываясь на результатах численного моделирования, предложена рабочая гипотеза изменения химического потенциала в области частицы, а значит, и кристаллизационной способности всего расплава.
3. Регулировать величины градиентов напряженности электрического поля и тепловой мощности в приповерхностной зоне частицы и в ее объеме можно изменением характеристик пропускаемого тока.
4. Предложенный алгоритм решения поставленной задачи может быть использован для анализа подобной электрофизической ситуации в случаях наличия ансамбля частиц различной формы, проводимости и размера.
5. Полученные в данной работе результаты могут послужить основой для решения термодинамической задачи, анализирующей образование и рост кристаллов при протекании через расплав электрического тока.
ЛИТЕРАТУРА
1. Цуркин В.Н., Синчук А.В., Иванов А.В. ЭОМ. 2011, (5), 89-98.
2. Цуркин В.Н. Процессы литья. 2015, (5), 3-10.
3. Борисов Г.П., Цуркин В.Н., Синчук А.В., Иванов А.В. Металлургия машиностроения. 2010, (5), 33-39.
4. Chun Yan Ban, Yi Han, Qi Xian Ba, Jian Zhong Cui. Mater Sci Forum. 2007, 546-549, 723-728.
5. Yunhu Zhang, Chanjiang Song, Liang Zhu, Hongxing Zheng et al. Metall Mater Trans B. 2011, 42(3), 604-611.
6. Guangming XU, Jiawei Zheng, Yong Liu, Jianzhong Cu. China foundry. 2005, 2(3), 171-175.
7. Цуркин В.Н. Металл и литье Украины. 2008, (9), 25-28.
8. Подольцев А.Д., Кучерявая И.Н. Техтчна елект-родинамка. 2015, (2), 3-15.
9. Ivanov A.V., Sinchuk A.V., Bogoslavskaya. Surf Eng Appl Electrochem. 2011, 47(1), 28-34.
10. Подольцев А.Д., Кучерявая И.Д. Элементы теории и численного расчета электромагнитных процессов в проводящих средах. Киев: Изд. Института электродинамики НАН Украины. 1999. 363 с.
11. Ефимов В.А. Разливка и кристаллизация стали. М.: Металлургия, 1976. 450 с.
12. Миссол В. Поверхностная энергия раздела фаз в металлах. Пер. с польск., M., 1978. 176 с.
13. Гиббс Дж.В. Термодинамика. Статистическая механика. М.: Наука, 1982. 584 с.
Поступила 30.11.16 После доработки 10.03.17 Summary
The two-dimensional formulation and the methods of numerical experiments were used to analyze the peculiarities of the redistribution of lines of a uniform electric field in the volume of molten metal at the interface of nano-sized phases with conductivity different from that of the matrix. It was shown that the gradients of the electric and thermal fields are determined by the shape of the phase and the ratio of its conductivity and of that of the matrix. The analyzed surface effects are also determined by the location sector of the inclusion of the surface relative to the flow lines.
Keywords: melt, electric current treatment, nano-sized particle, crystallization ability, energy factors.
