Математическая модель кинетики движения двухфазной зоны при нагреве цилиндрического тела из двухкомпонентного сплава Текст научной статьи по специальности «Физика»

УДК 544.3 + 519.63 + 517.956.4

DOI: 10.14529/теМ60407

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ КИНЕТИКИ ДВИЖЕНИЯ ДВУХФАЗНОЙ ЗОНЫ ПРИ НАГРЕВЕ ЦИЛИНДРИЧЕСКОГО ТЕЛА ИЗ ДВУХКОМПОНЕНТНОГО СПЛАВА

АД. Дрозин

Южно-Уральский государственный университет, г. Челябинск

В случае если в теле при его нагреве или охлаждении происходит фазовый переход, процесс распространения в нем тепла претерпевает сильные изменения. В отличие от одноком-понентного металла, при охлаждении многокомпонентного сплава, при переходе через температуру ликвидус, энтальпия, плотность и коэффициент теплопроводности остаются непрерывными, но их производные и, в частности, теплоемкость терпит разрыв. При нагреве поверхности твердого цилиндрического тела из двухкомпонентного металлического сплава выше температуры ликвидус у поверхности возникает двухфазная зона, перемещающаяся к оси цилиндра. В работе предложен метод введения «эффективной» теплоемкости, позволяющая рассчитывать скорость движения этой зоны, а также температуру сплава в любой его точке и в любой момент времени. Использовалось обычное уравнение теплопроводности с переменными коэффициентами, которые, в данном случае, зависят от температуры. Жидкая, твердая и переходная зоны образца рассчитывались по одному и тому же уравнению, коэффициенты его рассчитывались по одному и тому же методу. Для каждой точки образца для каждой соответствующей температуры, рассчитывались доли жидкой и твердой фаз. В данной работе, для простоты, раствор считался совершенным. Удельный объем, коэффициент теплопроводности и удельная энтальпия рассчитывались как средневзвешенные от соответствующих величин для жидкого и твердого состояний. Удельная теплоемкость рассчитывалась как производная по температуре от энтальпии. Полученную систему дифференциальных уравнений сводили к конечно-разностным уравнениям. Для решения полученной системы разностных уравнений была разработана компьютерная программа. В статье приведены результаты одного из таких расчетов. Результаты работы могут быть интересны не только металлургам, но и метрологам при разработке теории самотестирующегося датчика температуры.

Ключевые слова: теплопроводность; фазовый переход; сплав; модель.

Введение

Расчет скорости прогрева образца (детали при термообработке, кусков шихты в начальный период ее плавления) значительно осложняется наличием фазового перехода. На поверхности раздела фаз терпят разрыв все термические коэффициенты, выделяется теплота фазового перехода, граница раздела фаз движется. Однако положение еще более осложняется, если мы имеем дело не с однокомпонентным металлом, а со сплавом, состоящим из нескольких компонентов. В этом случае границей раздела твердой и жидкой фаз служит не поверхность, а двухфазный слой, в котором сосуществуют жидкость и твердые кристаллы [1-3].

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

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

В данной работе приведен метод расчета кинетики прогрева первоначально твердого двух-компонентного цилиндрического образца.

Физическая химия и физика металлургических систем

За базу было взято обычное уравнение теплопроводности [4-6] (ф)с(ф)5Т = 1 г—"

рч 'ср ' — =--1 А4 'г

дх г дг ^ дг

(ф)

в котором теплоемкость ср ' рассчитывалась как производная по температуре от энтальпии всей

сложной системы, плотность р(ф) и коэффициент теплопроводность рассчитывались как

средне взвешенные величины от соответствующих значений жидких и твердых составляющих сплава при данной температуре.

Для этого вначале рассчитаем необходимые термодинамические характеристики двухкомпо-нентного сплава, для простоты считая все фазы совершенными растворами.

1. Термодинамика двухкомпонентного сплава

Рассмотрим сплав с приведенной на рис. 1 диаграммой состояния.

1.1. Удельная энтальпия и удельный объем сплава

Будем снабжать величины, относящиеся к твердой фазе, верхним индексом 5", а к жидкой - L . Пусть Т5 - температура солидус для данного сплава, то есть температура, при превышении которой в твердом сплаве начинает появляться жидкая фаза, Ть - температура ликвидус для данного сплава, то есть температура, при снижении ниже которой в жидком сплаве начинает появляться твердая фаза. Обозначим символом X7 массовую долю 7-го компонента в исходном сплаве.

Рис. 1. Диаграмма состояния двухкомпонентного сплава

1.1.1. Энтальпии отдельных компонентов сплава

Принимаем, что для температур, соответствующих твердому состоянию, удельная энтальпия (¡7(5)) и удельный объем ( у7(5)) 7-го компонента сплава равны

¡(5)(Т) = ¡7 + ср )(Т - Т0 ), vг(5)=-1- (7 = 1,2) , (1.1)

Р7

где = Ь (т° ) - удельная энтальпия чистого 7 -го компонента при некоторой температуре Т7° ;

(5)

Ср7 - удельная теплоемкость 7 -го компонента; р7 удельная плотность 7 -го компонента сплава в твердом состоянии.

Для температур, соответствующих жидкому состоянию, удельную энтальпию (^ Ь)) и удельный объем (V-Ь ^) /-го компонента будем считать равными соответственно

= $ + ср\гз - Т-0 ) + с^(Т - Т3 ), (' = 1,2) . (1.2)

Р-

Здесь Т5 - температура плавления чистого - -го компонента; АИ-5 Ь) - удельная энтальпия его плавления.

1.1.2. Энтальпия и объем сплава

1.1.2.1. При температуре ниже Т5 весь сплав твердый. Его энтальпия и объем равны

НА (Т ) = £ т(5)h(S )(Т ), К = £ т(5 )v(5) . -=1 /=1

(5) 5Иа (5) 5УА -

Здесь п/ ' =—р-, V/ ' =—- парциальные удельные энтальпия и объем - -го компонента

дт] ) дт) )

сплава. За неимением лучших данных будем считать их равными соответствующим удельным величинам чистых компонентов и вычислять по формулам (1.1).

Так как весь сплав твердый, то т(5) = т/ = тАХ/, где т/ - масса - -го компонента во всем

т

сплаве, а X, = —'--его валовая доля. Тогда

МА

НА (Т) = тА £h'(S) (Т)X,, V = тА £ V'(5)Х-. (1.3)

-=1 -=1

1.1.2.2. При температуре Т5 < Т < Ть из твердой фазы начинает выделяться жидкая фаза. Процесс может протекать в два этапа. При снижении температуры ниже Ть вначале выделяется твердая фаза (5) - нужно рассматривать равновесие между жидкой фазой (Ь) и твердой

фазой (5) . Состав любой фазы Ф сплава характеризуется двумя долями его компонентов:

(ф) (ф)

„(ф) т\ '

Ф)= х2Ф)= ^ ( т = т 1Ф) + т2Ф)) . т т V '

Так как х[Ф) + х2Ф) = 1, то из двух концентраций вторая определяется первой, так что состав

фазы можно задавать долей первого компонента (х|Ф)) .

Пусть т( ) (Т) и т( ' (Т) - массы твердой и жидкой фаз сплава при данной температуре Т . Найдем их:

т(Ь) + т(5)= тА, (5) х1Ь) - X,

< —ч т ' = т а----—

| (£-)„(£), (5) „(5) „ А х(Ь) х(51) •

т 'Х\ '+ т 'Х\ ' = тАХ1 л} - л}

Введя обозначения

(^) Х1Ь)- Х1 (Ь) л (5) п ,,

а ил-, а 1 , (1.4)

Х|Ь) - Х|51)

получим

т(5^ = тАа(5), т(Ь) = тАа(Ь^. (1.5)

Величины (х1 Ь),х1 5)) для каждой температуры должны быть известны из диаграммы состояния (см. рис. 1):

Х1(ф)=Ф((ф)(Т), (ф = Ь, 5) .

Массы отдельных компонентов каждой фазы определяются по формуле

(ф) (ф)чл(ф)

т ' = тАау ' у к '

Х}Ф), (Ф = L, S; i _ 1,2).

Энтальпия и объем всего сплава при наших допущениях равны: НА (Т) = тА £ Рф)£(¡(ф)Хг(Ф))1 Ка (Т) = тА £ 1а(ф)Х (^ф)Х<ф))

ф=Ь,5 V 7=1 ) ф=Ь,5 V 7=1

(1.6)

(1.7)

(1.8)

(ф) дНА (ф) дКа л

Здесь ¡7 ' =—-ф-, V7 ' =—(ф) - парциальные удельные энтальпия и объем 7 -го компонента

дт( ) дт7( )

сплава. Так же, за неимением точных данных, будем считать их равными соответствующим удельным величинам чистых компонентов и вычислять по формулам (1.1), (1.2).

1.1.2.3. При температуре выше Ть весь сплав жидкий. Его энтальпия и объем равны

На (Т) = £т(ь¡Ь)(Т), V = £т^).

i=1

i=1

(ь) дНА (ь) дКА -

Здесь ¡7 ' =—р), V 7 =—- парциальные удельные энтальпия и объем 7 -го компонента

дт) ) дтг( )

сплава. Так же будем считать их равными соответствующим удельным величинам чистых компонентов и вычислять по формулам (1.2).

Так как весь сплав жидкий, то т(ь) = т7 = тАХ7. Тогда

HA (T) = mA £ h(L) (T)X,, FA (T) = mA £ v(L)X.

(1.9)

i=1

i=1

1.2. Определение эффективной теплоемкости

Эффективную теплоемкость сплава, которую мы будем использовать в уравнении теплопро-

е дН

водности, определим как сР = -

' A

dT dm,

1.2.1. Определение ¡А, v*4' для температур Т < Т5

Для температур ниже Т5, в соответствии с формулами (1.3), получим

сР=1 ¡;( 5) х,=£ $) х, .

(1.10)

i=1

i=1

1.2.2. Определение ИА', vJ для температур Т5 < Т < Ть Применяя формулы (1.8), получим

сР = I а(ф);£ (¡(ф) Х(ф)) + £ а(ф)£(

¡(ф) х(ф)')

ф=5 ,Ь V 7=1 ) ф=5 ,Ь V 7=1 )

Вычислим производные, входящие в это выражение. В соответствии с (1.4), (1.6),

а

(S >_.

xiS)'(x{L)- X1)+x{ L)'( X1 - x{ S))

а( L)'_-a(S)'.

( х< ь)- X 5 ))2

Хг(ф)'=ф(Ф)'(Т)(ф = 5,Ь; 7 = 1,2).

1.2.3. Определение НА V для температур Т > Ть

Для температур выше Тц, в соответствии с формулами (1.9), получим

сР =1 ь); X, =£ сьг> X,.

(1.11)

(1.12)

(1.13)

(1.14)

i=1

i=1

1.3. Алгоритм расчета эффективной теплоемкости

А1. Известны данные по сплаву: Состав сплава (массовые доли компонентов) - (Х1, Х2); температуры солидус, ликвидус и начала выпадения двухфазной эвтектики сплава Т5, Ть ; функции зависимости состава твердой и жидкой фаз сплава от температуры Т5 < Т < Ть; удельные те-

(5) (Ь) (5) (Ь)

плоемкости компонентов сплава в твердом и жидком состояниях - ср1', ср1', ср2, ср2 , теплоты плавления чистых компонентов сплава - ДЛ1(5 Ь), ДН25 Ь^; удельные плотности компонентов сплава в твердом и жидком состояниях - р(5),р|Ь^,р25), р2Ь) ; коэффициенты теплопроводности сплава в твердом и жидком состояниях - 5), Ь). А2. Расчет сплава: А2.1. Если Т < Т5 . А2.1.1. Для - = 1,2

н(5 )(Т ) = Н0 + ср/ )(Т - Т0).

А2.1.2. НА (Т ) = £ ^ )(Т) X, ср = Н*А = £ с? >Хг, vА = 0; vА = £Х-; V' = 0.

( 5 )'

х,

-=1

-=1

'=1 р/

( 5 )

А2.2. Если Т5 < Т < ТЬ . А2.2.1. Для ' = 1,2

Н5 )(т )=н?+ср: )(т - Т-°), V;5 ^

р/

(5 )

н(Ь)(т) = н? + ср5)(т -Т-0) + ДйГЬ)+сръг)(Т-т), v)ь) = *.

р/;

Для Ф = 5,Ь : х(Ф) = Ф(Ф) (Т), Х'(Ф)' = Ф(Ф> (Т).

Л")")") (5) Х1(Ь)- Х1 (Ь) , (51)

А2.2.2. а ;=—ттт-—; а ; = 1 -а ;;

(Ф)' _ „(Ф)',

Х1( Ь)- х<5 У

а

(51)'=.

151)'( х<Ь)-х1) + х1Ь)' (х1 - х151)) (х<Ь)- х<51))2 '

а( Ь)' = -а(51)'.

( 2

( 2

А2.2.3. Н*А = £ а(Ф) £h'(Ф)X'(Ф) , vА = £ а(Ф) £v)Ф)X'(Ф)

Л

Ф=Ь,5 V'=1 2

Ф=5 ,Ь

V '=1

2

сер = НА' = £ а(Ф)|£(срф)Х'(Ф)+ НФ)Х'(Ф)') + а(Ф)'£(н(Ф)х(Ф)) Ф=5 ,Ь V '=1 '=1

^ ДФ)^ у(Ф) Х (Ф)'+а(Ф)'^ V(Ф) Х (Ф)'^

А= £ а(Ф)£ ^Ф)Х}'

+ а-'£ v]ф) ХГ /=1

Ф=51,Ь V '=1

А2.3. Если Т > Ть . А2.3.1. Для / = 1,2

н(Ь (Т) = Н0 + с£) Т - Т0) + дН5-Ь} + ср, > (Т - Т ).

А2.3.2. НА = £Н'(Ь)(Т)X, , ср = НА'= £ср)Х1. vА = Ха£-|у, vА' = 0.

/=1

/=1

'=1 р/

V

1.4. Пример расчета

По изложенной методике рассчитаем необходимые термодинамические характеристики сплава Sn-Bi с содержанием В^ равным 20 % (ат.). По приведенной на рис. 2 диаграмме [7-9] определили уравнения соответствующих линий:

АпА1Е: хВ = -5,666962 • 10-5Т2 +1,488313 • 10-2 Т - 0,4038416; АВА2Е: хВ = 1,743208-10-5Т2 -3,90639•Ю-3Т + 0,776183; АпА3Е^ : хВ = 2,250585•Ю-6Т2 -3,181351 • 10-3Т + 0,6167243; АВЕВ : хВ = 1,28399•Ю-4Т + 0,965152568.

£п. № % £п

а0,6%В|) 58

Ав, 200

100

0 10 20 30 40 50 60 70 ВО 90 100 Эп В!,% (ат) В|

Рис. 2. Диаграмма состояния системы Sn-Bi

Из справочных данных [10-11] подобрали необходимые для расчетов данные (см. таблицу).

Параметры расчета

Параметр Sn Bi

Удельная теплота плавления (Дж/кг) hS ~L =59 000 hS ~L =11 300

Уд. теплоемкость (Дж/(кг К) тв. 41 =243 cp} =127

Уд. теплоемкость (Дж/(кг К) ж. 4? =240 4? =140

Плотность (кг/м3) тв. p1S > =7280 p2S> =9800

Плотность (кг/м3) ж. p1L) =6850 pL =10 300

Коэф. теплопроводностьи (Вт/(м К)) тв. l{S} =65,8 ) =6,37

Коэф. теплопроводности (Вт/(м К)) ж. a,il) =31,6 l!2L) =14,65

Атомная масса 118,69 208,98

Результаты расчетов приведены на рис. 3.

Рис. 3. Зависимости удельной энтальпии (а), эффективной теплоемкости (б), удельного объема (в) и коэффициента теплопроводности сплава Sn-Bi с 30,6 % (мас.) Bi от температуры

2. Расчет прогрева цилиндра из двухкомпонентного сплава 2.1. Постановка задачи

Рис. 4. Схема прогрева цилиндрического образца из двухкомпонентного сплава

Имеется цилиндр радиуса R, поверхность которого нагревается по закону Т R ) = ) (рис. 4). В результате цилиндр расплавляется (цилиндр находится в некоторой оболочке, которая

не дает ему растечься). В каждой точке сплава температура, в какой-то момент времени, поднимается до температуры солидус и начинает появляться жидкая фаза. Через какое-то время достигается температура ликвидус, оставшиеся кристаллы растворяются в жидкости, все система в этом месте переходит в жидкое состояние и продолжает нагреваться. Требуется найти функцию T (t, r) - температуру цилиндра в любой точке, находящейся на расстоянии r от оси цилиндра в момент времени t .

Подобная постановка задачи имеет отношение к термообработке деталей, к расплавлению кусков шихты. Она также интересна метрологам для обоснования теории самотестирующегося датчика [12-15].

2.2. Система уравнений

дт 1 д(. дт Л

pcp — = "4Xr— 1, 0 < r < R (2Л)

dt r dr ^ dr )

dT (t ,0)

-^-—t- = 0, (2.2)

dr

T(t,R) = |(t), t > 0, (2.3)

T(0,r) = ф(r), 0 < r < R. (2.4)

Cp = cP (T),p = pe (T), X = Xe (T). (2.5)

2.3. Метод решения

Уравнения системы записывались в конечно-разностном виде. Полученная система разностных уравнений решалась методом прогонки. Величины cp, pe, Xе рассчитывались по методике, описанной в части 1. Была составлена соответствующая компьютерная программа, которая и осуществляла соответствующие расчеты.

2.4. Пример расчета

В качестве примера приведены результаты расчета прогрева цилиндра диаметром 5 мм, первоначально имеющего температуру 20 °С, температура поверхности которого с момента t = 0 стала поддерживаться равной 300 °С (рис. 5, 6). Расчеты приведены для сплава Sn-Bi: с 69,4 % олова (температура ликвидус 200 °С, температура солидус 139 °С).

'ШВп=69,4 Centre

300 250 -200 ^ 150 -

100 -50 -

0-'-1-'-'-'-

0 0.02 0.04 0.06: 0.08 0.1 0.1-2

time, s

Рис. 5. Изменение температуры в центре образца со временем

Рис. 6. Распределение температуры по сечению образца через 0,001 с (а); 0,01 с (б); 0,05 с (в);

0,1 с (г) после начала прогрева

Заключение

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

Литература

1. Любое, Б.Я. Теория кристаллизации в больших объемах /Б.Я. Любое. - М.: Наука, 1975. -256 с.

2. Скрипов, В.П. Спонтанная кристаллизация переохлажденных жидкостей / В.П. Скрипов, В.П. Коеерда. - М.: Наука, 1984. - 232 с.

3. Авдонин, Н.А. Математическое описание процессов кристаллизации / Н.А. Авдонин. -Рига: Зинатне, 1980. - 180 с.

4. Тихонов, А.Н. Уравнения математической физики / А.Н. Тихонов, А.А. Самарский. - М. : Наука, 1972. - 736 с.

5. Самарский, А.А. Теория разностных схем /А.А. Самарский. - М. : Наука, 1977. - 656 с.

6. Самарский, А.А. Методы решения сеточных уравнений / А.А. Самарский, Е.С. Николаев. -М. : Наука, 1978. - 592 с.

7. Диаграммы состояния двойных металлических систем: справ. / под. ред. Н.П. Лякишева. -М.: Машиностроение, 1996. - 992 с.

8. Захаров, А.М. Диаграммы состояния двойных и тройных систем / А.М. Захаров. - М: Металлургия, 1990. - 250 с.

9. Справочник по пайке / под ред. И. Е. Петрунина. - М.: Машиностроение, 2003. - 480 с.

10. Физические величины. Справочник / под ред. И.С. Григорьева, Е.З. Мейлихова. - М: Энер-гоатомиздат, 1991. - 1250 с.

11. Рабинович, В.А. Краткий химический справочник / В.А. Рабинович, З.Я. Хавин. - Л: Химия, 1978. - 392 с.

12. Шестаков, А.Л. Модель самодиагностирующегося датчика параметра с нелинейной функцией преобразования /А.Л. Шестаков, А.С. Семенов //Измерение. Мониторинг. Управление. Контроль. - 2015.- № 1 (11). - С. 17-22.

13. Шестаков, А.Л. Методы теории автоматического управления в динамических измерениях /А.Л. Шестаков. - Челябинск: ЮУрГУ, 2013.

14. Самокалибрующийся термометр на основе точек плавления. Конструкция и алгоритмы работы / М.Д.Белоусов, В.В. Дьячук, Д.А. Мирзаев, А.Л. Шестаков // Вестник ЮУрГУ. Серия «Компьютерные технологии, управление, радиоэлектроника». - 2013. - Т. 13, № 1. - С. 26-33.

15. Белоусов, М.Д. Оценка собственного состояния термометров сопротивлений / М.Д. Белоусов, А.Л. Шестаков, Н.М. Япарова // Вестник ЮУрГУ. Серия «Компьютерные технологии, управление, радиоэлектроника». - 2012. - № 35. - С. 105-109.

Дрозин Александр Дмитриевич, д-р техн. наук, профессор кафедры материаловедения и физико-химии материалов, Южно-Уральский государственный университет, г. Челябинск; drozi-nad@susu.ru.

Поступила в редакцию 7 октября 2016 г

DOI: 10.14529/met160407

MATHEMATICAL MODEL OF THE TWO-PHASE ZONE MOTION DURING HEATING OF THE TWO-COMPONENT CYLINDRICAL OBJECT

A.D. Drozin, drozinad@susu.ru

South Ural State University, Chelyabinsk, Russian Federation

If due to heating or cooling of an object, a phase transition occurs, the propagation of heat significantly changes. Unlike one-component metal, upon cooling of the multicomponent alloy, when passing through the liquidus temperature, enthalpy, density and thermal conductivity are continuous, but their derivatives and, in particular, the heat capacity become discontinuous. By heating the surface of the solid metal cylindrical object with two-component composition above the liquidus temperature, the two-phase zone occurs at the surface, and moves to the cylinder axis. We propose a method of introducing "effective" specific heat that allows calculating the speed of this zone movement, as well as the temperature of the object at any point and at any time. The ordinary equation of heat conductivity with variable coefficients which, in this case, depend on temperature was used. Liquid, solid and transitional zones of a sample were counted by the same equation, its coefficients were counted by the same method. For each point of the sample for each respective temperature the proportion of solid and liquid phases was calculated. In this paper, for simplicity, the solution was considered as ideal. Specific volume, thermal conductivity and specific enthalpy were calculated as the weighted average of the corresponding values for the liquid and solid states. The specific heat was calculated as the derivative of the enthalpy with respect to temperature. The resulting system of

differential equations was reduced to finite-difference equations. A computer program was developed to solve the resulting system of difference equations. Results of one of such calculations are given in the article. The developed method allows to calculate the velocity of the phase boundary, as well as the object temperature at any point and at any time. The results may be interesting not only for metallurgists, but for metrology to develop self-testing temperature sensor theory.

Keywords: heat conduction; phase transition; alloy; model.

References

1. Lyubov B.Ya. Teoriya kristallizatsii v bol'shikh ob"emakh [Theory of Crystallization in Large Volumes]. Moscow, Nauka Publ., 1975. 256 p.

2. Skripov V.P., Koverda V.P. Spontannaya kristallizatsiya pereokhlazhdennykh zhidkostey [Spontaneous Crystallization of Undercooled Liquids]. Moscow, Nauka Publ., 1984. 232 p.

3. Avdonin N.A. Matematicheskoe opisanie protsessov kristallizatsii [Mathematical Description of Crystallization Processes]. Riga, Zinatne Publ., 1980. 180 p.

4. Tikhonov A.N., Samarskiy A.A. Uravneniya matematicheskoy fiziki [Equations of Mathematical Physics]. Moscow, Nauka Publ., 1972. 736 p.

5. Samarskiy A.A. Teoriya raznostnykh skhem [Theory of Difference Schemes]. Moscow, Nauka Publ., 1977. 656 p.

6. Samarskiy A.A., Nikolaev E.S. Metody resheniya setochnykh uravneniy [Methods of Solving Grid Equations]. Moscow, Nauka Publ., 1978. -592 p.

7. Diagrammy sostoyaniya dvoynykh metallicheskikh sistem: Spravochnik [Phase Diagrams of Binary Metal Systems. Reference Book]. Vol. 1. Lyakishev N.P. (Ed.). Moscow, Mashinostroenie Publ., 1996.992 p.

8. Zakharov A.M. Diagrammy sostoyaniya dvoynykh i troynykh sistem [Phase Diagrams of Binary and Ternary Systems]. Moscow, Metallurgiya Publ., 1990. 250 p.

9. Spravochnikpo payke [Soldering Reference Book]. Petrunin I.E. (Ed.). Moscow, Mashinostroenie Publ., 2003. 480 p.

10. Fizicheskie velichiny. Spravochnik [Physical Values. Reference Book]. Grigor'ev I.S. and Mey-likhov E.Z. (Eds). Moscow, Energoatomizdat Publ., 1991. 1250 p.

11. Rabinovich V.A., Khavin Z.Ya. Kratkiy khimicheskiy spravochnik [Brief Chemical Reference Book]. Leningrad, Khimiya Publ., 1978. 392 p.

12. Shestakov A.L., Semenov A.S. [Model of a Self-Diagnosing Parameter Sensor with Nonlinear Transformation Function]. Izmerenie. Monitoring. Upravlenie. Kontrol', 2015, no. 1 (11), pp. 17-22. (in Russ.)

13. Shestakov A.L Metody teorii avtomaticheskogo upravleniya v dinamicheskikh izmereniyakh [Methods of Automatic Control Theory in Dynamic Measurements]. Chelyabinsk, South Ural St. Univ. Publ., 2013.

14. Belousov M.D., D'yachuk V.V., Mirzaev D.A., Shestakov A.L. [Self-Calibrating Thermometer Based on Melting Points. Construction and Operation Algorythms]. Bulletin of the South Ural State University. Ser. Computer Technologies, Automatic Control & Radioelectronics, 2013, vol. 13, no. 1, pp. 26-33. (in Russ.)

15. Belousov M.D., Shestakov A.L., Yaparova N.M. [Estimation of Eigenstate of Resistive Thermometers] . Bulletin of the South Ural State University. Ser. Computer Technologies, Automatic Control & Radioelectronics, 2012, no. 35 (294), issue 17, pp. 105-109. (in Russ.)

Received 7 October 2016

ОБРАЗЕЦ ЦИТИРОВАНИЯ

Дрозин, А.Д. Математическая модель кинетики движения двухфазной зоны при нагреве цилиндрического тела из двухкомпонентного сплава / А.Д. Дрозин // Вестник ЮУрГУ. Серия «Металлургия». - 2016. -Т. 16, № 4. - С. 67-77. Ш1: 10.14529/теШ0407

FOR CITATION

Drozin A.D. Mathematical Model of the Two-Phase Zone Motion During Heating of the Two-Component Cylindrical Object. Bulletin of the South Ural State University. Ser. Metallurgy, 2016, vol. 16, no. 4, pp. 67-77. (in Russ.) DOI: 10.14529/met160407