Анализ устойчивости ячеистого фронта кристаллизации и формирование структуры затвердевающего расплава

Качество твердых металлов, получаемых путем направленного затвердевания расплавов, определяется в значительной мере рельефом поверхности раздела фаз.

В данной работе будет рассматриваться бинарный сплав, то есть сплав, состоящий из основного компонента и растворенного. Растворенный компонент имеет более низкую температуру плавления но сравнению с основным (в дальнейшем его будем называть примесью).

Если на межфазной поверхности имеет место равновесие и отсутствует переохлаждение перед зарождающимся кристаллом, а также переохлаждение, связанное с кривизной поверхности раздела, то затвердевание происходит плоским фронтом. Реально же при кристаллизации перед плоским фронтом образуется граничный слой, обогащенный растворенным компонентом.

Образование граничного слоя происходит вследствие зависимости температуры ликвидуса T_L от состава расплава (рис 1а): чем меньше примеси, тем выше Ть поэтому сначала кристаллизуется сплав с меньшей долей примеси (eg <Со), чем исходный состав, при этом избыток растворенного компонента оттесняется в жидкую фазу. Тогда непосредственно пред фронтом кристаллизации оказывается расплав, обогащенный примесью, концентрация которой убывает по мере удаления от фронта до Со (рис. 16). Накопление примеси вблизи фронта приводит к тому, что в этой области Tl ниже, чем в окружающем расплаве (рис.1в). Плоский фронт будет устойчив, если все точки перед фронтом находятся при температуре выше Ть, так как любой выступ, образовавшийся на гладкой поверхности раздела, оказывается окруженным перегретой относительно T_L жидкостью и снова расплавится, иначе плоский фронт теряет устойчивость [8].

По мере своего роста возмущения формы фронта стабилизируются, при этом граница раздела приобретает форму регулярных гексагональных ячеек [3].

Кроме того, при описании причин неустойчивости плоского фронта следует учесть различные условия теплоотвода от выступающих участков фронта и от плоской межфазной поверхности, а также роль межфазного поверхностного натяжения, которое затрудняет процесс роста выступа, смещая температуру равновесия фаз в область более низких температур (выступ менее устойчив по отношению к

пульсациям температуры, чем плоский фроит)[3,7].

Расстояние до фронта

а                                 б                                  в

Рис. 1 Схемы к анализу граничного слоя: а - участок равновесной фазовой диаграммы; б - распределение примеси вблизи фронта, в - распределение температуры вблизи фронта

Итак, при нарушении устойчивости плоского фронта происходит образование правильных ячеек. Появление же дендритов связано с неустойчивостью ячеистых структур, которая обусловлена, по видимому, теми же причинами, что и неустойчивость плоского фронта кристаллизации. Таким образом, устойчивость ячеистого фронта определяется, с одной стороны, концентрационным переохлаждением, с другой, - кривизной фронта и межфазным поверхностным натяжением Анализ устойчивости позволяет найти зависимости критических характеристик неустойчивых возмущений от критических значений параметров процесса, это позволит выявить те параметры процесса (параметры порядка), которые управляют поведением системы в момент перехода от устойчивости к неустойчивости. Характеристики неустойчивых возмущений можно связать с параметрами структуры кристаллизующегося расплава; в частност и, длина волны возмущений плоского фронта может рассматриваться как приближенная оценка характерного расстояния между дендритами, а длина волны возмущений ячеистого фронта - между вторичными ветвями дендрита

В работах [2,7] рассматривается устойчивость плоского фронта кристаллизации, в первом случае - с учетом концентрационного переохлаждения, во втором - без учета. В представленной работе делается попытка рассмотреть устойчивость ячеистого фронта кристаллизации с учетом как концентрационного переохлаждения, так и межфазного поверхностного натяжения и кривизны поверхности раздела. Эго позволит выявить основные параметры, влияющие на устойчивость фронта, а значит, и на формирование структуры кристаллизующегося расплава. Согласно [3] первичная структура закристаллизовавшегося расплава определяет его служебные свойства, поэтому анализ устойчивости дает ключ к возможности управления этими свойствами

Постановка задачи

Будем рассматривать устойчивость ячеистого фронта кристаллизации, считая, что ячейки имеют форму параболоида вращения (z- ось вращения, р0 - радиус кривизны вершины параболы), эго соответствует росту в переохлажденном расплаве иглообразного кристалла. Такая форма поверхности и температурное поле вокруг нее были получены при некоторых предположениях ГП. Иванцовым |4]

Процесс кристаллизации металлического расплава с учетом перераспределения растворенного компонента описывается следующей системой уравнений, дополненной условиями Стефана на поверхности раздела, описанными ниже:

30:             до

—  _ а Д0     —- = од0  1 - 12,

(О

3t             at где ©, - температура в жидкой и 02 - в твердой фазах, о-распределение концентрации примеси в жидкой фазе. Диффузией в твердой фазе и движением расплава пренебрегаем. Зависимость физических характеристик материала от температуры не учитывается. Плотность твердой и жидкой фаз считается одинаковой На фронте кристаллизации выполняются условия баланса тепла и массы:

302     30,

^2 "Т" - ^1     ~ YL" ■ v ’ dn      dn

do

D—- =-(1 - k)on • v, reU dn где U = (г / E(t, г) = о}, E(t, г) - 0 задает границу раздела фаз; у , L- плотность, удельная теплота плавления; к - коэффициент распределения примеси (равен отношению концентрации в момент времени t к начальной концентрации Со);

= п • V() - производная по направлению единичного вектора нормали к поверхности фронта, v - скорость перемещения поверхности раздела, поверхность раздела пока будем считать изотермической:

где 0L - температура фазового перехода без учета перераспределения примеси.

Вдали от фронта задаем условия:

• •    в расплаве

• •    в твердой фазе

т-^с0, ©2 ->©

vo >

В начальный момент времени:

02 •

о(о,г) = с0,

0,   ®2(°,»')=©02.

Считается, что теплоотвод осуществляется только через твердую фазу.

Изотермическая поверхность раздела, имеющая форму параболоида вращения, является устойчивой в квазистационарном режиме, то есть поля температур и концентрации перемещаются без изменения конфигурации, как жесткое целое, со скоростью V. Так как поверхность раздела фаз изотермическая, вектор ее нормали параллелен градиенту температуры (и концентрации)[1], тогда на поверхности

d©2(r,t) а©2

d t 3©2 ИЛИ —-- 3t

at з©7

+ V • VO2 = О

—— = 0 Зп

V = -

Перемещение поверхности кристалла происходит в направлении, противоположном тепловому потоку, поэтому v • n = v и условие (2,) на поверхности преобразуется следующим образом:

а©2    а©, а©2     а©2 з©2

yL-™ = X,                    reU.

St       dn dn    ² dn dn

Введем безразмерные параметры:

о

• •    безразмерная концентрация примеси: с = — - 1; ^со

безразмерная теплоемкость: Л -=

Х2

V

безразмерная теплота кристаллизации: R =

ЬГУРо

2X20l ’

безразмерная величина наклона линии ликвидуса на фазовой диаграмме: М = тсо ©L

• •    безразмерная кривизна поверхности раздела (для невозмущенного параболоида вращения она равна 1, но отлична от 1, когда на поверхность раздела накладываются возмущения): К] -Кр₀;

Г

• •    безразмерный коэффициент поверхностного натяжения: Tj = —— ;

1 к безразмерный коэффициент распределения примеси: Р = —

где v - модуль вектора скорости перемещения поверхности раздела вдоль оси z, m -тангенс угла наклона линии ликвидуса на фазовой диаграмме, К - кривизна поверхностности раздела, Г - коэффициент межфазного поверхностного натяжения

Учитывая осевую симметрию поверхности раздела, можно ограничиться рассмотрением сечения параболоида плоскостью yOz, то есть будет решаться плоская задача. Пусть параболический фронт движется таким образом, что вершина параболы перемещается вдоль оси z со скоростью v (рис.2), ось у направлена перпендикулярно оси z, начало координат находится в фокусе параболы, тогда уравнения (1) запишутся в виде

Рис.2. Расположение системы координат в начальный момент времени

ат, Га2^ а2т7 at ay2 + az2 J ас (а2с а2с) at -Dlay2 +az2)'

Решения уравнений (8) будем искать в виде функции сложного безразмерного аргумента [2]:

^(^■’ИЛ =Ц%.)-«с* т

р =z - vt+77 ~ ^vo² + у²

Для каждого фиксированного момента времени аргумент и меняется от 0 - бесконечно удаленные точки в твердой фазе, до <х> - бесконечно удаленные точки в расплаве, и=1 соответствует поверхности раздела фаз. Подставляя (9) в систему (8) и используя правила дифференцирования сложной функции (Т; = Т; (u(t, у, z)) ), получим d2T dT fl рЛ du2 du k2u 2) d2c def 1 qf —г + — — + — = 0, du2 duk2u 27

vpo i = 1,2, где p: =---- ai

vp0 где q = —

условия на границе раздела фаз:

dT2

dT, de    q

__L _ n ___ „ JL p du   ’du   2 ’

T] = T2 = 1, u = 1;

вдали от границы раздела фаз:

'1 г0 ’

Tj ~> Г()2>

Решения уравнений (10) с учетом соответствующих граничных условий имею!

вид:

Т.(и) = Т0 +

c(u) = ^ ехр(%) Р Erfc^^j;

2 z 2

где Erf(z) = -т= J e"1 dt, Erfc(z) = 1 - Erf(z) [8]; соотношения (13) определяют 0

квазистационарное решение задачи. Условия теплового баланса на границе раздела фаз примут вид

G =МР, то есть G=M— и = 1, с                        с du

L

Q = —, где cv - теплоемкость.

cv

Область концентрационного переохлаждения перед фронтом появляется, если градиент фактической температуры в жидкости пересекает равновесную кривую температуры ликвидуса (рис. 1в) [3]. Чтобы учесть влияние концентрационного переохлаждения и поверхностного натяжения, нужно изменить температурные условия на границе раздела фаз (11):

Т] =Т2 =1-М(с + 1) + Г]К], и = 1.                        (15)

Заметим, что кривизна параболы различна в вершине и вдали от нее, это приводит к невыполнению условия изотермичности поверхности раздела фаз и приведенные выкладки становятся неправомерными. Однако, если рассматривать область непосредственно вблизи вершины параболы, то условие изотермичное™ поверхности раздела будет приближенно выполнено, и тогда можно в соотношениях (13,14) изменить константы интегрирования с учетом (15).

Ограничение рассмотрения только области в окрестности вершины параболы оправдано тем, что именно на вершине дендрита при потере устойчивости появляются вторичные ветви [3]. Исследуем устойчивость вершины параболы, внося возмущения на фронте в виде функции z - vt = Z,(t,y), где Z](t,y) - функция, описывающая возмущения поверхности раздела (будем считать их возмущениями плоского фронта z-vt = Ро/2 , но при этом будем учитывать характер распределения температуры и концентрации вблизи вершины параболического фронта ); эти возмущения сопровождаются возмущениями Т; = 1) - Tjs , с =c-cjs ,1=1,2 квазистационарных полей. Величины возмущений считаем малыми, чтобы можно было привести граничные условия при u' = z-vt = Z,(t,y) к условиям на. u = z-vt = p0/2, используя разложение в ряд Тейлора

Компоненты вектора нормали и кривизна вычисляются по формулам

а затем

Записывая систему (1)-(4) для возмущенного и невозмущенного решения, вычитая вторые из первых, получим уравнения эволюции пульсаций:

ат  f а2т'   а2т)   ас'  f а2с   а2с') .

a  Lay2   az2 /  at  lay2   az2;

граничные условия на Zi приводятся к условиям на u = z-vt = p0/2 по методике, изложенной в [3], используя разложение в ряд по Zi и сохраняя только линейные члены. Учтем концентрационное переохлаждение на вершине параболы dT      dT. dT. de de

—-  <—-    T—-  = -M— du u=,    du U=1 de du U=1 du U=1

в виде Tj(u')-T1(1)<-m(c(u')-c(1)), где и’ обозначена возмущенная поверхность, после разложения в ряд по Z], сохраняя только линейные члены:

или P^^GcZ,. Тогда граничные условия (2)-(3) на поверхности раздела для пульсаций представляются в виде ряда по Zi (с учетом разложения выражения для кривизны (16) в биномиальный ряд с отрицательным показателем [1]).

, ат' ,      5с' a2Z,

Т +-— Z, + Мс + М— -Z, -Г,—7 + P.Z, - G„Z,= 0 . 1 = 1,2,

‘ dz                  dz ^{1    1} ay*

8г' аг az. ат₂ az, ат ат₂ ат az, ^Л "аГ "' аГ ” ^Л Lay""^ ⁺ "а_у^{4 Л21} "az"^{4 Z1} az" ~ ^R"aT

ас                        i - k az,                 р0

-—+ (l-k)c +(l-k)Z, +——57 ^{= 0}’ ^пр^{и и}'⁼у

Граничные условия на бесконечном удалении от поверхности раздела(4):

с' —> О Т<' < <ю z —><ю

1 1(19)

T2'->0,z->-oo.

Начальные условия следуют из вида возмущений, которые задаются гармониками [ti.c Л^АДгХВДгХСДгХО^1^(20)

Тогда решения системы с учетом граничных условий (19):

(т2, Т], с’ , Z,,} = {А,e“2Z, В.е0'7, С, D, }ept+i“’y,(21)

Г 2 р             2 р [ 2 р где а^ = Jro +—, a2=Jco + —, Р =        , Аь Вь Ci - константы интегрирования, определяемые из (18). Итак, подставляя (21) в (18), получим систему линейных алгебраических уравнений, из которой найдем константы интегрирования:

О-A! ч-е^^В, тМе^»''2^ Ц^го2 +Pj Gc)d! =0, е«2Ро/2д +0 В, тМе^^С, + (Гро2 + Р, -Gc)d, = 0, J        1               111        1     V J 1     1

' Xa₂e^azPo/IA._v + a₁e-“-^p<,/2B₁ + 0-RpD, =0,                             ^(22j

O' А1 + 0B, + (p + l-k)e"Ppo/2C1 +^l + k + -^-p]D1 =0.

Система (22) разрешима, если ее определитель равен 0, то есть

M^l - k + —- (р + 1 - к)(Г]® 2

+ Р] - Gc) • (Ла 2 + a 1) - Rp(p +1 - к) = 0.

Соотношение (23) определяет в неявном виде зависимость инкремента нарастания возмущений р от физических свойств расплава и параметров процесса; оно задает некоторую гиперповерхность р(со, Pi, R, Gc,...) в пространстве параметров [3]. Нарушению устойчивости соответствует изменение знака действительной части р с отрицательного на положительный. Полагая в (23) р=0, получаем уравнение границ устойчивости и неустойчивости в пространстве параметров и волнового числа со . Критическое значение некоторого параметра, например Gc= Gc* определяется из условия, что кривая ц=ц(со, Gc) (при фиксированных остальных параметрах),представляющая собой сечение гиперповерхности плоскостью Gc= Gc* , лежит ниже плоскости р-0 везде, за исключением точки со=со., в которой она касается плоскости. Тогда при Gc< Gc* спонтанно возникающие возмущения формы фронта затухают; при Gc= Gc* появляются незатухающие возмущения с волновым числом со.; при Gc> Gc* существует непрерывный интервал волновых чисел, соответствующих нарастающим возмущениям. Таким образом, можно установить зависимость параметров дендритной структуры от какого-либо параметра процесса, зафиксировав остальные параметры. Например, длину волны возмущения (измеряемую в метрах) 1, --2к/®«р0 на ячеистом фронте можно рассматривать в первом приближении как расстояние между вторичными ветвями дендрита, так как появление вторичных ветвей на главной оси дендрита связано с потерей устойчивости вершины дендрита; ветви, появившиеся на вершине дендрита, в дальнейшем или продолжают расти, или расплавляются. Линейный анализ устойчивости дает лишь критерий их появления; для того, чтобы проследить эволюцию ветвей, необходима нелинейная постановка задачи.

Используя физические константы, приведенные в [9], были получены некоторые результаты. На рис.3-5 приведены зависимости длины волны возмущения от параметров М, Gc ,Fi , характеризующих наклон линии ликвидуса на фазовой диаграмме, влияние коэффициента диффузии в жидкой фазе и кривизну поверхности раздела; с увеличение М и Gc длина волны устойчивого возмущения уменьшается, с увеличением Гi - увеличивается.

М0Л-4 М0л-4

2 к

1,8

1’6| \

• 1,4 '            '

1,2

• 1                                                              .

м

Рис.З. Длина волны критического возмущения в зависимости от параметра М

5 р.

_ ;   100000 200000 300000 400000 500009е

Рис.4. Длина волны критического возмущения в зависимости от параметра Gc

МО-М

10000 20000 30000 40000 50000

Рис. 5. Длина волны критического возмущения в зависимости от параметра Fj

Итак, в рамках линейной теории устойчивости рассматривается устойчивость возмущений на ячеистом фронте кристаллизации с учетом концентрационного переохлаждения и межфазного поверхностного натяжения. Получен критерий устойчивости фронта, который позволяет установить зависимости некоторых параметров дендритной структуры от параметров процесса кристаллизации и физических констант.