InTechOpen uses cookies to offer you the best online experience. By continuing to use our site, you agree to our Privacy Policy.

Materials Science » "Solidification", book edited by Alicia Esther Ares, ISBN 978-953-51-3906-5, Print ISBN 978-953-51-3905-8, Published: March 14, 2018 under CC BY 3.0 license. © The Author(s).

Chapter 3

Coupled Model of Precipitates and Microsegregation During Solidification

By Yandong Li, Huamei Duan and Cheng Peng
DOI: 10.5772/intechopen.72597

Article top


The simplified diagram of coupling process (fs—solid fraction).
Figure 1. The simplified diagram of coupling process (fs—solid fraction).
SEM images and EDS analysis of inclusions in the solidification structure.
Figure 2. SEM images and EDS analysis of inclusions in the solidification structure.
The flow chart of technological process of lab scale experiment.
Figure 3. The flow chart of technological process of lab scale experiment.
SEM images and EDS analysis of inclusions in the laboratory sample.
Figure 4. SEM images and EDS analysis of inclusions in the laboratory sample.
The precipitation behavior of inclusions during solidification of 253MA steel.
Figure 5. The precipitation behavior of inclusions during solidification of 253MA steel.
The precipitation behavior of inclusions calculated by FactSage.
Figure 6. The precipitation behavior of inclusions calculated by FactSage.
The inclusions types in heat-resistant steel at 1873 K.
Figure 7. The inclusions types in heat-resistant steel at 1873 K.
The stability areas of inclusions at 1873 K.
Figure 8. The stability areas of inclusions at 1873 K.
The liquid inclusions in heat-resistant.
Figure 9. The liquid inclusions in heat-resistant.

Coupled Model of Precipitates and Microsegregation During Solidification

Yandong Li1, Huamei Duan2 and Cheng Peng1
Show details


To understand the precipitation behavior and solidification process of micro-alloyed steel, abundant thermodynamic data of pure substances were incorporated in the coupled thermodynamic model of inclusions precipitation and solutes micro-segregation during the solidification of heat-resistant steel containing cerium. The liquid inclusions Ce2xAl2ySi1−xyOz where 0 < x < 1, 0 < y < 1−x and z = 1−xy and their generation Gibbs free energy were first introduced to the inclusions reactions according to the Al2O3-SiO2-Ce2O3 phase diagram. Then plant trials, lab experiments and published work in the literature were taken account to valid the established model. Also, the difference of calculated results between FactSage and this model were argued. Finally, the liquid inclusions were found in the samples from experiments in the tube furnace based on the calculations by this model.

Keywords: liquid inclusions, thermodynamic model, rare earth, heat-resistant steel

1. Introduction

The heat resistant steels containing rare earth elements (RE) have absolute advantages in mechanical properties and corrosion resistance [1, 2, 3, 4]. The oxidation resistance of 253MA steel at high temperature is dramatically enhanced as alloyed with cerium element. However, the rare earth elements are used as deoxidizing and desulfurizing agents in the earlier ages owing to the strong oxidation property [5]. As a result, a large number of inclusions such as Ce2O3 and Ce2O2S with high melting temperature and small size are generated in the molten steel. During the continuous casting process these inclusions are easily captured by the submerged nozzle. Then the accumulated inclusions with frozen steel lead to nozzle blocking and the production process is interrupted [6, 7, 8, 9]. As one of the effectual solutions to nozzle blocking calcium treatment can modify the solid alumina clusters to liquid inclusions which have weak affinity force with the nozzle inner wall [10, 11, 12, 13]. Inspired by the above method, the high melting point inclusions containing rare earth elements may be modified to liquid phase by some alloys to solve the nozzle clogging issue of heat resistant steel.

Many prediction models of inclusions formation and evolution are reported in the previous reports. However, the released self-built model mainly involve a few equilibrium reactions of alloy elements with oxygen or sulfur at 1600°C [14, 15]. And the thermodynamic properties of complex liquid inclusions, such as Mn-Si-Al-O, Ca-Al-O, Ca-Al-Mg-O, are calculated by the commercial software [15, 16, 17, 18, 19]. Harada assessed the CaO-Al2O3 liquid inclusions formation setting the fixed ratio of CaO to Al2O3 instead of the whole liquid region [20, 21]. FactSage as the common software in metallurgical field is strong at the equilibrium calculations of multi phases, but the FToxide database is short of solution thermodynamic data containing rare earth and the solidification model is not thoroughly optimized. Liu released one coupling model of micro-segregation and inclusion formation during solidification process in silicon steels [22]. However, the number of reactions was restricted and the solution phase was not considered. Fang proposed a thermodynamic model which includes plenty of formation reactions of rare earth inclusions except the liquid inclusions as a pity [23]. To make up for the deficiency of previous models [24], the more ampler generation reactions of inclusions especially the liquid phase coupling with the solute micro-segregation are introduced in the present model. Then the coupled model was validated by various methods, such as plant trials, lab scale experiments and released experimental data. Finally, the formation regions of liquid phases in the predominance diagrams were given. Moreover, the liquid inclusions with sphere shape and large size were found in the exploitative experiments of lab scale.

2. Modeling process

The coupling model is composed of two parts: the thermodynamic model of inclusions precipitation both in liquid and solid and the dynamic model of solute micro-segregation model in the liquid–solid interface during solidification. The schematic diagram of coupling process is shown in Figure 1. Considering the kinetic factors such as diffusion rate and atomic size, the inclusions with large size are formed in the residual liquid phase above the solidus temperature and the small-size carbonitrides are precipitated in the solidified steel under the liquidus temperature. All the relevant reactions of inclusions generation both in liquid and solidified steel are listed in Table 1 [26, 27, 28, 29, 30, 31, 32, 33] and Table 2 [34, 35, 36, 37], respectively. The activity standard state of solute elements and inclusions was 1 mass% infinite dilute solution and the pure substance, respectively. The activity coefficients of solute elements were arranged by Wagner’s relations expanded to second order interaction coefficient. The solute segregation at liquid–solid interface during solidification was calculated form the Brody-Flemings micro-segregation model modified by Clyne-Kurz [25]. The more details including the mathematic expressions and thermodynamic properties of solutes can be consulted from the previous works [24].


Figure 1.

The simplified diagram of coupling process (fs—solid fraction).

Reaction equations in molten steelΔGθ = A + B × T, J·mol−1
2[Ce] + 3[O] = Ce2O3(s)−1,431,090360.06
[Ce] + 2[O] = CeO2(s)−854274.7249.11
[Ce] + [S] = CeS(s)−422,783120.58
[Ce] + 2[S] = CeS2(s)−131,000121.4
2[Ce] + 3[S] = Ce2S3(s)−1,074,584328.24
3[Ce] + 4[S] = Ce3S4(s)−1,493,010438.9
2[Ce] + 2[O] + [S] = Ce2O2S(s)−1353592.4331.6
[Ce] + [Al] + 3[O] = CeAlO3(s)−1,366,460364
[Ce] + 11[Al] + 3[O] = CeAl11O18(s)−7261120.002282.36
2[Ce] + [Si] + 5[O] = Ce2SiO5(s)−2077830.00632.57
4.67[Ce] + 3[Si] + 13[O] = Ce4.67Si3O13(s)−5299450.001641.95
2[Ce] + 2[Si] + 7[O] = Ce2Si2O7(s)−2800700.00909.66
2x[Ce] + 2y[Al] + (1-x-y)[Si] + z[O] = Ce2xAl2ySi1-x-yOz(l)
[Ce] + [N] = CeN(s)−401,200153
[Ce] + 2[C] = CeC2(s)−202,790125.3
2[Ce] + 3[C] = Ce2C3(s)−224,000205.8
[Al] + [N] = AlN(s)−267,520119.54
2[Al] + 3[O] = Al2O3(s)−1,205,090387.72
[Mn] + [S] = MnS(s)−172676.156.07
[Mn] + [O] = MnO(s)−284,420122.75
[Si] + 2[O] = SiO2(s)−580,550221.03
6[Al] + 2[Si] + 13[O] = Al6Si2O13(s)−4,758,4341579.22
4[Cr] + [C] = Cr4C(s)−195,790218
23[Cr] + 6[C] = Cr23C6(s)−887889.71253.96
7[Cr] + 3[C] = Cr7C3(s)−356120.1417.62
3[Cr] + 2[C] = Cr3C2(s)−182030.1207.39
2[Cr] + 3[O] = Cr2O3(s)−797,190349.72
3[Cr] + 4[O] = Cr3O4(s)−94439.9416.79
[Cr] + [O] = CrO(s)−236,320113.6
2[Cr] + [N] = Cr2N(s)−141,300116.82
[Cr] + [N] = CrN(s)−136,25096.2
[Ni] + [O] = NiO(s)−125,000136.29

Table 1.

The possible reactions and standard Gibbs free energy in liquid.

Reaction equations in solid phase(γ)ΔGθ = A + B × T, J·mol−1
[Al] + [N] = AlN(s)−137536.134.27
[Mn] + [S] = MnS(s)−176,78257.91
[Ti] + [N] = TiN(s)−153158.236.13
[Ti] + [C] = TiC(s)−143585.8452.65
[Nb] + [N] = NbN(s)−162730.6253.6
[Nb] + [C] = NbC(s)−143585.8456.67
[V] + [N] = VN(s)−182030.1207.39
[V] + [C] = VC(s)−797,190349.72

Table 2.

The possible reactions and standard Gibbs free energy in solid.

The generated liquid inclusions whose activity is 1 can be identified as one isolated phase (pure liquid phase) out of the solid inclusions (pure solid phase) and solute elements (solutes in the infinite dilute solution) in the system. As the same as listed in Table 1, the liquid inclusion in the form of Ce2xAl2ySi1−xyOz (l, 0 < x < 1, 0 < y < 1x and z = 1xy) is introduced according to the phase diagram of Al2O3-SiO2-Ce2O3 system. The hard task is the calculation of standard Gibbs free energy of Ce2xAl2ySi1-x-yOz generation as the following reaction:


The chemical reaction of Eq. 1 can be derived from the following reactions:


where the standard Gibbs free energy of reactions 2, 3, 4 are listed in Table 1. And the CALPHD technique was used to optimize the Gibbs free energy change of reaction 8. The Redlich-Kister (R-K) expression was adopted to describe the excess Gibbs free energy of solution phases in Al2O3-SiO2, Al2O3-Ce2O3 and Ce2O3-SiO2 binary systems. The R-K expression is written as Eq. 9.


where GE is the excess Gibbs free energy of binary system; x1, x2 are the mole fraction of end members in the binary system; Lj (j = 0, 1, 2…) are the interaction parameters of components in the binary system. And Lj in Table 3 were obtained by the phase diagram optimization when the restricted experimental points from literature were employed.

Li = a + bT/(J·mol−1)L0L1L2L3

Table 3.

The optimized coefficients of the binary systems.

The excess Gibbs free energy of Ce2O3-Al2O3-SiO2 ternary system was derived by the following geometric Kohler method [38].


where the first term G12Ex1x1+x2x2x1+x2 is the excess Gibbs free energy at x1x1+x2x2x1+x2 point in the 1–2 binary system. The following two terms in Eq. 10 have similar representation as the first term. In order to exactly describe the thermodynamic properties of Ce2O3-Al2O3-SiO2 system, the Gibbs free energy was further modified by the GternaryE term with the value of 825000xCe2O3xAl2O3xSiO22 . After the interaction coefficients in Table 3 are imported, the standard Gibbs free energy of reaction 8 is written as:


When the Gibbs free energy of Eqs. (2)(8) are introduced, the standard Gibbs free energy of Eq. (1) is finally written as follows:


After the standard Gibbs free energy ΔGLθ of liquid inclusions formation was worked out, the inclusions compositions (certain x, y values, z = 1−x−y) can be calculated according to the minimum of Gibbs free energy ΔGL . At the equilibrium status in the molten steel, the Gibbs free energies of chemical reactions in Table 1 are greater than or equal to zero.

3. Modeling validation

The constructed model was validated by plant trial, lab scale experiments and experimental data from others work. Moreover, the calculated results by this model were compared with the popular commercial software FactSage.

The standard chemical compositions of 253MA steel were shown in Table 4. And the industrial samples of 253MA steel were obtained after EAF → AOD → LF → CC process in one stainless steel plant of China. The obtained samples were machined into column of Φ = 4 mm and then polished by emery paper and polishing cloth to make metallographic specimen. Then the field emission scanning electron microscope (FE-SEM) and energy dispersive spectrometer (EDS) were employed to observe the morphology and analyze the compositions of inclusions, as shown in Figure 2. It can be seen that the inclusions with the size of 1-2 μm in the solidification structure of 253 MA steel are mainly Ce2O2S and Ce2O3. CeN maybe exist as the element N was detected as shown in Figure 2(b). The peaks without labels are matrix elements Fe, Cr and Ni.


Table 4.

The chemical compositions of 253MA steel (mass fraction, %).


Figure 2.

SEM images and EDS analysis of inclusions in the solidification structure.

The experiments in laboratory that carefully follows the industrial operations were carried out in the Si-Mo rod heating furnace protected by argon flow at 1873 K, as shown in Figure 3. Firstly, the Fe-Cr-Ni alloys with fixed ratio as in the 253MA steel were loaded in the MgO crucible until completely melted. Then the mass fraction of Si was calculated and added into the crucible after the oxygen content was first determined. After Si was totally melted, the oxygen content was determined for the second time and then the melt was alloyed by Ce. The molten steel was drawn using quartz tube (inner Φ = 4 mm) and immediately quenched into ice water at 30 min after Ce addition. Therefore, the inclusions formed at 1873 K can be preserved since the extremely cooling rate of quenching. The quenched samples were processed and observed as the same as that of industrial samples. It is explicit that the inclusions are Ce2O2S and Ce2O3 with the size smaller than 2 μm as the analysis results shown in Figure 4.


Figure 3.

The flow chart of technological process of lab scale experiment.


Figure 4.

SEM images and EDS analysis of inclusions in the laboratory sample.

The types and amounts of inclusions formed during the solidification of 253MA steel was calculated by the present model and the results are drawn in Figure 5 which (b) is the zoom view of the liquid–solid two-phase region (TL = 1713 K). It can be seen that the inclusions types are mainly Ce2O2S, Ce2O3, CeN and Ce2SiO5. And the calculation results are in agreement with the industrial experiments except Ce2SiO5. Ce2SiO5 is formed as Ce2O3 is reduced during solidification as shown in the magnified image of Figure 5(b). Then Ce2SiO5 must be transformed from Ce2O3 as the increasing Si concentration in the residual liquid phase which owes to the micro-segregation of solutes at the liquid–solid interface. However, it is very difficult for the decomposition of Ce2O3 because the extremely fast cooling rate during the continuous casting process. Therefore, Ce2SiO5 is not found in the industrial samples. What is more, the formed inclusions are only Ce2O3 and Ce2O2S at 1873 K which is completely consistent with the results in the laboratory scale experiments.


Figure 5.

The precipitation behavior of inclusions during solidification of 253MA steel.

The types and amounts of inclusions formed during the solidification of 253MA steel was also calculated by FactSage and the results are drawn in Figure 6. It can be seen that only Ce2O3 and Ce2O2S are formed in the FactSage results where Ce2O2S precipitates at a lower temperature 1848 K comparing with 1873 K in the present model. However, Ce2SiO5 and CeN are not present in the FactSage results since the solutes segregation is not considered in FactSage when the normal equilibrium solidification model is selected. The elements activities were further calculated to understand the lower formation temperature of Ce2O2S in FactSage. The stable pure substances at standard conditions are adopted by FactSage as the activities standard sate when 1mass% infinite dilute solution is chosen in the present work. The activities of solutes in molten steel calculated by the two different standard states can be linked by the following equation.


where aiwH is the activity at the standard state of 1 mass% solution; aiR is the activity when pure substance is standard state; Mi is the relative atomic mass of solutes; MA is the relative atomic mass of solvents, that is Fe in this work; γi0 is the activity coefficient when stable pure substance is standard state in the Henrian solution. The available activity coefficients are listed in Table 5 [39] and the activities transformed from FactSage state to 1 mass% solution comparing with that calculated by this model are present in Table 6.


Figure 6.

The precipitation behavior of inclusions calculated by FactSage.

Standard stateSSLLG

Table 5.

Optimized Henrian activity coefficients lnγ0.

Note: lnγ0 = A/T + B, e: natural constant, S: pure solid, L: pure liquid, G: gas.

This work4.23E-23.105.74E-12.78E-27.15E-41.50E-22.09E-32.82E-4

Table 6.

The activities calculated by this work and FactSage at 1873 K.

As listed in Table 6, the solutes activities calculated from different standard state are basically at the same level except P and N elements. The remarkable differences of P and N activities can owe to the missing γi0 values of P and N. It can be seen that the sulfur activity calculated by FactSage is smaller than that calculated by this model because γi0 of S is not available. It is the direct explanation why the formation temperature of Ce2O2S is lower in the FactSage results.

Besides the industrial trials and lab scale experiments, the types of inclusions reported in the steel containing rare earth elements from literature were compared with the model prediction and the results are present in Table 7. The types of generated inclusions in various grades of steels reported in others work are almost the same as the results calculated by this model.

AuthorsRE addition/%Inclusions foundPredicted by this model
Yue [40]0.016Rare earth sulfidesCe2S3
Liu [41]0.044Rare earth oxysulfidesCe2O2S, CeS
Guo [42]0.034Ce2O2S, CeSCe2O2S, Ce2O3, CeN
Zhai [43]0.08Ce2O2S, Ce2O3, CeO2Ce2O2S, Ce2O3
Factsage0.04Ce2O3, Ce2O2SCe2O3, Ce2O2S, Ce2SiO5, CeN
Present work0.04Ce2O3, Ce2O2S, CeNCe2O3, Ce2O2S, Ce2SiO5, CeN

Table 7.

Comparison of inclusions from literature experiments and predicted by this work.

All the above comparisons send the message that the accomplished coupled thermodynamic in this work is precise and it is competent to exactly assess the precipitation behaviors of inclusions during solidification process of various grades of steels containing rare earth elements.

4. Model application

Aluminum was alloyed in the heat-resistant steel imitating the calcium treatment way to solve the nozzle clogging of 253MA steel. Since the formed inclusions resulting in nozzle clogging are mainly oxides, the inclusions transformation at different initial oxygen contents was studied. And the types of inclusions formed before and after aluminum addition are shown in Figure 7. It can be seen that the stable inclusion at high oxygen content is Ce2Si2O7 in steel without aluminum. By contrast, the inclusions eventually transform to liquid phase as oxygen content increases to 0.018% after 0.01mass% aluminum was added in 253MA steel.


Figure 7.

The inclusions types in heat-resistant steel at 1873 K.

In addition, the stable areas of inclusions (predominance area diagram of inclusions or phase diagram of inclusions) in heat-resistant steel at 1873 K are summarized in Figure 8. The stable regions of liquid inclusions are surrounded by the dash line. There is not liquid inclusions formed when the oxygen content is low (0.005% in this study). And the stable areas of liquid inclusions are growing as the increasing initial oxygen content. There can be no doubt that stable areas of liquid phase will be helpful to solve the nozzle clogging issue of steels alloyed by rare earth elements.


Figure 8.

The stability areas of inclusions at 1873 K.

In 2011, Kojola [6] reported that the clogging rate of 253MA steel dramatically decreased when Al, Si and Ce co-existed in the molten steel. Kojola argued that the small rare earth oxides can be dissolved after Si addition because of the negative interaction of Si to oxygen. However, no convincing proof was given. It is clear that the liquid inclusions are formed according to the prediction by the model in this work. Moreover, based on the calculated stable areas of inclusions the lab scale experiments following the procedure in Figure 3 were carried out when Si-Al alloy was chosen instead of pure Si. The FE-SEM image and EDS analysis of found inclusions were shown in Figure 9(a). The inclusions were complex phases of Ce-Al-Si-O system with the large size of 15 μm. And the complex inclusions were in the liquid region at 1873 K when be marked in the Al2O3-SiO2-Ce2O3 phase diagram.


Figure 9.

The liquid inclusions in heat-resistant.

5. Conclusions

A coupling model of inclusions generation and solutes micro-segregation during the solidification process of steels containing rare earth elements was worked out. The method for calculating the Gibbs free energy of liquid inclusions was first given. Then the reliability of mathematical model was validated by the industrial sampling, experiments in lab scale and inclusions found in literature. The inclusions calculated by FactSage software were different from this model because of the differences in the activity standard state and the solidification model. The inclusions are Ce2O3, Ce2O2S and CeN in the solidification structure of 253MA steel. These small size inclusions containing rare earth elements can be modified to the liquid inclusions after Al adding and the stable liquid regions increase as the initial oxygen content increases. Moreover, the complex liquid inclusions of Ce-Al-Si-O system were found in the preliminary experiments of lab scale. The coupling model is helpful to relieve the nozzle clogging issue of various grades of steel containing rare earth elements.


This work was supported by the Scientific Research Foundation for Scholars in Yangtze Normal University No. 2016KYQD14, No. 2016XJQN34.


1 - Chen L, Ma XC, Wang LM, Ye XN. Effect of rare earth element yttrium addition on microstructures and properties of a 21Cr-11Ni austenitic heat-resistant stainless steel. Materials and Design. 2011;32(4):2206-2212. DOI: 10.1016/j.matdes.2010.11.022
2 - Matway RJ, McGuire MF, Mehta J. Steel alloy having improved creep strength. US. Patent. 1995
3 - Nilsson JO. The influence of nitrogen on high temperature low cycle fatigue behaviour of austenitic stainless steels. Fatigue & Fracture of Engineering Materials & Structures. 1987;7(1):55-64. DOI: 10.1111/j.1460-2695.1984.tb00359.x
4 - Cai GJ, Li CS. Effects of Ce on inclusions and corrosion resistance of low-nickel austenite stainless steel. Materials and Corrosion. 2015;66(12):1445-1455. DOI: 10.1002/maco.201508380
5 - Waudby PE. Rare earth additions to steel. International Metals Reviews. 1978;23(1):74-98. DOI: 10.1179/imtr.1978.23.1.74
6 - Kojola N, Ekerot S, Andersson M, Jönsson P. Pilot plant study of nozzle clogging mechanisms during casting of REM treated stainless steels. Ironmaking and Steelmaking. 2011;38(1):1-11. DOI: 10.1179/030192310X12690127076398
7 - Kojola N, Ekerot S. Pilot plant study of clogging rates in low carbon and stainless steel grades. Ironmaking & Steelmaking. 2011;38(2):81-89. DOI: 10.1179/030192310X12706364542704
8 - Roos E, Karasev A, Jönsson PG. Effect of Si and Ce contents on the nozzle clogging in a REM alloyed stainless steel. Steel Research International. 2015;86(11):1279-1288. DOI: 10.1002/srin.201400344
9 - Yao YK, Zhu MW, Wang DY, Liu CJ, Jiang MF. Investigation of nozzle blockage mechanism during feeding RE into tundish. Chinese Eare Earth. 2004;25(5):17-19. DOI: 10.16533/j.cnki.15-1099/tf.2004.05.006
10 - Lind M, Holappa L. Transformation of alumina inclusions by calcium treatment. Metallurgical and Materials Transactions B. 2010;41(2):359-366. DOI: 10.1007/s11663-009-9337-9
11 - Janke D, Ma Z, Valentin P, Heinen A. Improvement of castability and quality of continuously cast steel. ISIJ International. 2000;40(1):31-39. DOI: 10.2355/isijinternational.40.31
12 - Holappa L, Hämäläinen M, Liukkonen M, Lind M. Thermodynamic examination of inclusion modification and precipitation from calcium treatment to solidified steel. Ironmaking & Steelmaking. 2003;30(2):111-115. DOI: 10.1179/030192303225001748
13 - Yang J, Wang XH, Jiang M, Wang WJ. Effect of calcium treatment on non-metallic inclusions in ultra-low oxygen steel refined by high basicity high Al2O3 slag. Journal of Iron and Steel Research, International. 2011;18(7):8-14. DOI: 10.1016/S1006-706X(11)60083-6
14 - Zhang LF, Ren Y, Duan HJ, Yang W, Sun LY. Stability diagram of Mg-Al-O system inclusions in molten steel. Metallurgical and Materials Transactions B. 2015;46(4):1809-1825. DOI: 10.1007/s11663-015-0361-7
15 - Zhang TS, Min Y, Jiang MF. Effect of magnesium addition on evolution of inclusions in Mn–Si–Al deoxidised molten steels. Canadian Metallurgical Quarterly. 2015;54(2):161-169. DOI: 10.1179/1879139514Y.0000000173
16 - Fujii K, Nagasaka T, Hino M. Activities of the constituents in spinel solid solution and free energies of formation of MgO, MgO·Al2O3. ISIJ International. 2000;40(11):1059-1066. DOI: 10.2355/isijinternational.40.1059
17 - Park JH, Todoroki H. Control of MgO Al2O3 spinel inclusions in stainless steels. ISIJ International. 2010;50(10):1333-1346. DOI: 10.2355/isijinternational.50.1333
18 - Ren Y, Zhang LF, Yang W, Duan HJ. Formation and thermodynamics of Mg-Al-Ti-O complex inclusions in Mg-Al-Ti-deoxidized steel. Metallurgical and Materials Transactions B. 2014;45(6):2057-2071. DOI: 10.1007/s11663-014-0121-0
19 - Zhang TS, Min Y, Liu CJ, Jiang MF. Effect of Mg addition on the evolution of inclusions in Al–Ca deoxidized melts. ISIJ International. 2015;55(8):1541-1548. DOI: 10.2355/isijinternational.ISIJINT-2014-691
20 - Harada A, Maruoka N, Shibata H, Kitamura S. A kinetic model to predict the compositions of metal, slag and inclusions during ladle refining: Part 1. Basic concept and application. ISIJ International. 2013;53(12):2110-2117. DOI: 10.2355/isijinternational.53.2110
21 - Harada A, Maruoka N, Shibata H, Kitamura S. A kinetic model to predict the compositions of metal, slag and inclusions during ladle refining: Part2. Condition to control the inclusion composition. ISIJ International. 2013;53(12):2118-2125. DOI: 10.2355/isijinternational.53.2118
22 - Liu ZZ, Wei J, Cai KK. A coupled mathematical model of microsegregation and inclusion precipitation during solidification of silicon steel. ISIJ International. 2002;42(9):958-963. DOI: 10.2355/isijinternational.42.958
23 - Fang L, Liu HL, Liu CJ, Jiang MF. Model analysis on behavior of cerium in heavy rail steel. Chinese Rare Earth. 2008;29(1):88-91. DOI: 10.3969/j.issn.1004-0277.2008.01.022
24 - Li YD, Liu CJ, Li CL, Jiang MF. A coupled thermodynamic model for prediction of inclusions precipitation during solidification of heat-resistant steel containing cerium. Journal of Iron and Steel Research, International. 2015;22(6):457-463. DOI: 10.1016/S1006-706X(15)30027-3
25 - Clyne TW, Kurz W. Solute redistribution during solidification with rapid solid state diffusion. Metallurgical and Materials Transactions A. 1981;12(6):965-971. DOI: 10.1007/BF02643477
26 - Vahed A, Kay DAR. Thermodynamics of rare earths in steelmaking. Metallurgical Transactions B. 1976;7(3):375-383. DOI: 10.1007/BF02652708
27 - Turkdogan ET. Physical Chemistry of High Temperature Technology. New York: Academic Press; 1980
28 - Alaoua D, Lartigue S, Larere A, Priester L. Precipitation and surface segregation in low carbon steels. Materials Science and Engineering A. 1994;189(1–2):155-163. DOI: 10.1016/0921-5093(94)90411-1
29 - Yamamoto K, Hasegawa T, Takamura JI. Effect of boron on intra-granular ferrite formation in Ti-oxide bearing steels. ISIJ International. 1996;36(1):80-86. DOI: 10.2355/isijinternational.36.80
30 - Taylor KA. Solubility products for titanium-, vanadium-, and niobium-carbide in ferrite. Scripta Metallurgica et Materialia. 1995;32(1):7-12
31 - ParK JY, Park JK, Choo WY. Effect of Ti addition on the potency of MnS for ferrite nucleation in C–Mn–V steels. ISIJ International. 2000;40(12):1253-1259. DOI: 10.2355/isijinternational.40.1253
32 - Chen EB, Dong YC, Guo SZ. Study on thermodynamical properties in Mn-Fe alloy melts. Acta Metallgugica Sinica. 1997;33(8):831-837
33 - Du T. Themodynamics of rare earth elements in iron-base solutions. Journal of Iron and Steel Research. 1994;6(3):6-12. DOI: 10.13228/j.boyuan.issn1001-0963.1994.03.002
34 - Yue EB, Qiu ST, Gan Y. Carbo-nitride precipitation model in austenite of HSLA steel. Journal of Iron and Steel Research. 2006;18(8):49-52. DOI: 10.3321/j.issn:1001-0963.2006.08.012
35 - Tang GB, Wu XY, Yong X, Bai AM, Liu ZD, Yong QL. Numerical simulation of thermodynamic for carbonitride precipitationin austenite of complex microalloyed high strength and low alloy steels. Heat Treatment Metals. 2008;33(8):67-72. DOI: 10.13251/j.issn.0254-6051.2008.08.044
36 - Luo S, Zhu MY, Ji C, Cai ZZ. Solute microsegregation model for continunous casting process of steel. Iron and Steel. 2010;45(6):31-36
37 - Zhang MC, Li SZ, Zhang QY, Zhang JP. Thermodynamic study on carbide and nitride precipitation of X60 pipeline steel in continuous casting. Journal of Anhui University of Technology (Natural Science). 2012;29(3):207-210. DOI: 10.3969/j.issn.1671-7872.2012.03.004
38 - Pelton AD. A general “geometric” thermodynamic model for multicomponent solutions. Calphad. 2001;25(2):319-328. DOI: 10.1016/S0364-5916(01)00052-9
39 - Jung IH, Decterov SA, Pelton AD. A thermodynamic model for deoxidation equilibria in steel. Metallurgical and Materials Transactions B. 2004;35(3):493-507. DOI: 10.1007/s11663-004-0050-4
40 - Yue LJ, Wang LM, Han JS. Effects of rare earth on inclusions and corrosion resistance of 10PCuRE weathering steel. Journal of Rare Earths. 2010;28(6):952-956. DOI: 10.1016/S1002-0721(09)60219-2
41 - Liu X, Yang JC, Gao XZ. Effects of RE on the inclusions and impact toughness of 2Cr13 stainless steel. Journal of University of Science and Technology Beijing. 2010;32(5):605-609
42 - Guo F, Lin Q, Sun XX. Thermodynamic calculation on precipitation and transformation of inclusions in high purity carbon-manganese steels with rare earths and their microstructure analysis. Journal of the Chinese Rare Earth Society. 2004;22(5):624-618. DOI: 10.3321/j.issn:1000-4343.2004.05.006
43 - Zhai QJ. Study on microstructure and formation mechanism of rare earth inclusions during solidification (Chinese edition). Journal of the Chinese Rare Earth Society. 1995;13(1):60-62