Next Article in Journal
Expired Glucosamine Drugs as Green Corrosion Inhibitors for Carbon Steel in H2SO4 Solution and Synergistic Effect of Glucosamine Molecules with Iodide Ions: Combined Experimental and Theoretical Investigations
Previous Article in Journal
Li-Rich Layered Oxides: Structure and Doping Strategies to Enable Co-Poor/Co-Free Cathodes for Li-Ion Batteries
Previous Article in Special Issue
Influence of Fly Ash on the Fluidity of Blast Furnace Slag for the Preparation of Slag Wool
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessment of Inclusion Removal Ability in Refining Slags Containing Ce2O3

1
College of Materials Science and Engineering, Taiyuan University of Technology, Taiyuan 030024, China
2
State Key Laboratory of Advanced Metallurgy, University of Science and Technology Beijing, Beijing 100083, China
3
Department of New Energy and Materials Engineering, Shanxi Electronic Science and Technology Institute, Linfen 041075, China
4
The Third Geological Center Laboratory, China Metallurgical Geology Administration, Taiyuan 030027, China
5
Shanxi Taigang Stainless Steel Co., Ltd., Taiyuan 030030, China
6
Department of Chemistry & Chemical Engineering, Lvliang University, Lvliang 033001, China
7
College of Materials Science and Engineering, Chongqing University, Chongqing 400044, China
*
Authors to whom correspondence should be addressed.
Crystals 2023, 13(2), 202; https://doi.org/10.3390/cryst13020202
Submission received: 21 November 2022 / Revised: 17 January 2023 / Accepted: 17 January 2023 / Published: 23 January 2023
(This article belongs to the Special Issue Metallurgical Slag (Volume II))

Abstract

:
The elimination of inclusions in steelmaking processes has been widely studied. The removal of inclusions by slags containing the rare earth oxide Ce2O3 are studied using an integrated numerical model. The integrated model involves the inclusion motion model, interfacial tension calculation model, surface tension calculation model of slag, and the mass action concentration model, based on ion and molecule coexistence theory. The motion behaviors of both solid Al2O3 inclusions and 50%wtAl2O3–50%wtCaO liquid inclusions of varied sizes at CaO-Ce2O3-SiO2-Al2O3(-MgO) slag systems are evaluated. The results show that it is more difficult to remove the inclusions with smaller sizes and in slag with a higher viscosity. Liquid inclusions are more difficult to remove than solid inclusions. It is found that the CaO-Ce2O3-SiO2-Al2O3-MgO refining slag shows a better ability to remove Al2O3 inclusions than that of the CaO-SiO2-Al2O3-MgO slag. The reason for this is that the addition of the rare earth oxide Ce2O3 can decrease the viscosity of slags, as well as improving the wetting effects of slags on Al2O3 inclusions. For two slags systems, the CaO-Ce2O3-SiO2-Al2O3-MgO slag system shows a better ability to remove Al2O3 inclusions than the CaO-Ce2O3-SiO2-Al2O3 slag system. The addition of 5% to 8% Ce2O3 in a CaO-SiO2-Al2O3-MgO slag is an optimized case for industrial applications.

1. Introduction

Non-metallic inclusions are harmful to steel product quality [1,2,3], and they also influence the steelmaking process via clogging of the nozzles [4,5,6]. A great number of studies have focused on the analysis of inclusions during the whole steelmaking process [7,8,9,10,11,12,13,14,15]. Meanwhile, the inclusion formation and control in a steelmaking unit, for example, a ladle furnace [16,17,18], vacuum degassing [19,20,21,22], Ruhrstahl–Heraeus [23,24,25,26,27], and tundish [28,29,30,31] have also been widely studied.
Absorption by slag [32,33,34] and adhesion to refractories [35] are two common ways to remove inclusions. Adherence to rising bubbles is a good method for the floatation of inclusions. The process of inclusion absorption by the top slag can be divided into four steps [36]: (1) inclusions are transported to the turbulent boundary layer of the steel–slag interface; (2) inclusions transport through the boundary layer to the steel–slag interface; (3) inclusions separate to the slag; and (4) inclusions are dissolved in the slag phase.
For step (1), a great number of studies have focused on the nucleation, motion, collision, and transport of inclusions in secondary refining units by using fundamental water models [37,38,39,40,41,42], mathematical models [43,44,45,46,47], and Computational Fluid Dynamics (CFD) models [48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66]. Both the Eulerian approach (the statistical evaluation of inclusion species) and the Lagrangian approach (the evaluation of particles) were used in the models. However, as reviewed by C. Chen et al. [66], steps (2)–(4) of the removal of inclusions were neglected in most CFD studies. This is challenging as it is a combination of all the steps. The set-up of inclusion removal in slag is according to the Stokes rising velocity in most Eulerian approaches, following the work of Tacke and Ludwig [67]. In Lagrangian approaches, the inclusion particles were commonly treated as being fully absorbed by the slag once the particle touched the interface. A novel approach was developed [68] by considering a critical velocity condition to evaluate the trap or to reflect the conditions of inclusions. A more complicated approach was developed recently [69].
As mentioned above, step (1) can be modeled elaborately in the CFD models. One attempt was made to combine steps (1)–(2) in a CFD model [66,70]. It was found that the interface slip velocity has a greater influence on the removal of small-sized inclusions. Step (2), i.e., the deposition of inclusions at the boundary layer of steel–slag interface, is dominant and cannot be neglected for small inclusions. In contrast, the influence of the boundary layer on large-sized inclusions is negligible since the inclusions are likely to penetrate the boundary layers. Thus, steps (3) and (4) are considered to be key steps for the removal of inclusions in general descriptions [32].
The studies on step (3) could be dated back to a pioneer work by Nakajima and Okamura [71] in 1992. They developed a model for single solid spherical inclusion motion at the steel–slag interface. This model could be categorized into a Lagrangian approach, which was based on Newton’s second law. The model has been continuously developed by introducing the turbulent burst theory to calculate the possibility of the re-entrainment of inclusions [72], by applying it to solid inclusions of other compositions [15,33,73] and liquid inclusions [74,75], by extending the model to different shapes of inclusions, such as octahedral and plate-like, etc. [76,77,78], by introducing the slag dissolution kinetics [79], by improving the theory of a slag film capture [80], and by applying it to a study on the water model of particles [81,82]. Recently, CFD models coupled flow–solid interaction [83,84] or the use of the dynamic overset grid technique [85] were developed to study the inclusion motion and separations. Furthermore, the motions of inclusions under a supergravity field were studied from a new perspective [86].
For step (4), inclusions dissolution in slags has been widely studied using experiments [34,87,88,89,90,91,92,93,94] and numerical models [79,95,96,97,98]. It worth mentioning that the oxide dissolution in slags have been studied using a novel in situ observation method [92,98]. In addition, phase field models have been developed to study the dissolution of inclusions [95,96,97,98]. Recently, a concept of the inclusion capacity of slag is put forward [94].
Typically, slag is formed in many processes [99,100,101,102,103,104,105,106,107]; for example, ironmaking, convertor steelmaking, refining, continuous casting, the electroslag remelting process, and ferroalloys production. The thermodynamic properties and physicochemical properties of many types of slags have been widely evaluated and modeled [108,109,110,111,112,113,114]. Typical refining slags are CaO-Al2O3-SiO2 slag and CaO-SiO2-Al2O3-MgO slag systems. The inclusions dissolution in two slag systems have been studied [90,91,92,115].
For rare earth (RE) oxides, it has been added in liquid steel, as mentioned, and reviewed by Hasegawa and Sakuma [116] in 1956. It is noted that the rare earth metal addition in liquid steel in ingot casting processes have been widely investigated in 1940–1970s [116,117,118]. The clogging of the submerged entry nozzles of RE oxide inclusions have hindered the application of RE element addition in steels in continuous casting routes [119]. Later, in the 1970s, the phase diagram of Al2O3-REOX was systematically measured by Mizuno et al. [120,121]. An attempt on the activity measurement in steelmaking slags CaO-Al2O3-Ce2O3 was first performed by Ueda et al. [122]. The thermodynamic properties, including phase diagram [123,124,125], structure [126,127,128], the activity and mass action concentration [129,130,131]; and physicochemical properties, including the viscosity [127,128,131], melting point [132,133], and surface tension [134] of the slag were studied. According to the research results, the addition of rare earth oxides in refining slag can reduce the activity of Al2O3 in slag and increase the viscosity and the melting point of slag. Attempts were also made to evaluate the effects of rare earth oxides-containing slags on steel cleanliness [135], as well as the performance of the microalloying of rare earth elements in steels by reducing the rare earth oxides in slag [136,137]. The microalloying method may be an alternative way to add rare earth elements into steel. Specifically, the balance between the Ce2O3-containing refining slag and the molten steel-containing aluminum will lead to a small amount of cerium being dissolved into the molten steel. This rare earth element can perform a vital role in purifying, modifying, and micro-alloying steels [138], modifying the inclusions [139,140,141,142], refining the solidification structure [143], and improving the performance [144,145,146] of RE steels. The effects of rare earth elements on the modification of inclusions have been widely studied, both in lab experiments and within industrial scale production. The rare earth can react with O and S to form rare earth oxygen sulfide [139]. The sizes of inclusions are also decreased. However, the mechanism of inclusion absorbed by the rare earth oxide-containing slag is still unclear.
In this paper, the separation of inclusions to the rare earth oxide-containing slag are studied. The motion behaviors of both solid Al2O3 inclusions and the liquid inclusions of varied sizes at CaO-Ce2O3-SiO2-Al2O3 and CaO-Ce2O3-SiO2-Al2O3-MgO slag systems are evaluated using an integrated numerical model. The integrated model involves the inclusion motion model, the interfacial tension calculation model, the surface tension calculation model of slag, and the mass action concentration model, based on the ion and molecule coexistence theory. This paper is organized as follows: all the model descriptions, validations, and calculation process parameters are presented in Section 2. The results and discussions are presented in Section 3. The results can be used to predict the movement of inclusions at the steel–slag interface, with different overall wettabilities of the steel–slag-inclusion system. It will provide hints for the design of the refining slag composition, and an evaluation method of inclusion removal ability in the slags.

2. Modeling Methodology and Model Details

The studied CaO-Ce2O3-SiO2-Al2O3 and CaO-Ce2O3-SiO2-Al2O3-MgO slag systems are denoted as slag A and slag C, respectively. The scope and flow chart of the models are shown in Figure 1.
The inclusion motion model requires the input of the properties of slags and inclusions; for example, the density and viscosity of slags, and the overall wettability of the steel–slag-inclusion system. The density and viscosity data of the two slag systems are referred from [87,133]. The overall wettabilities of steel–slag inclusions require interfacial tension between each of the two components, i.e., the interfacial tension of steel–slag, steel–inclusion, and slag–inclusion. The estimation of the interfacial tension via the Girifalco-Good equation may further require the surface tension of steel, inclusion, and slag, respectively. Thereafter, the data of the inclusion and steel have been widely studied. The surface tension of rare earth oxide-containing slag is a challenging work. The estimation model using the Butler equation and the individual surface tension of the oxide component can be utilized. In the model, the molar fraction and the activity of the oxide component are required data. Alternatively, the mass action concentration (activity) of the oxide could be calculated using the well-known ion and molecule coexistence theory.
The previous analysis followed a sequence from the modeling strategy perspective. In this section, the inverse order of model descriptions will be presented. Specifically, the surface tension of slag will be described in Section 2.1. The ion and molecule coexistence theory, and the Butler equation-based model will be briefly introduced in Section 2.1.1 and Section 2.1.2, respectively. The specific surface tension of the CaO-Ce2O3-SiO2-Al2O3 and CaO-Ce2O3-SiO2-Al2O3-MgO slag systems will be illustrated in Section 2.1.3 and Section 2.1.4, respectively. Secondly, the interfacial tension of steel–inclusion, slag–inclusion, and steel–slag will be calculated in Section 2.2.1, Section 2.2.2 and Section 2.2.3, under Section 2.2. The yellow dotted line box in Figure 1 shows the calculation flow in Section 2.2, and the arrow indicates that the parameter value obtained in the previous step is substituted into the formula for the next step. For example, the yellow arrow indicates substituting the value of the surface tension of solid Al2O3 or liquid 50%wtAl2O3–50%wtCaO inclusion and the surface tension of liquid steel into the Girifalco-Good equation. The interfacial tension between liquid steel and inclusion can be calculated by substituting the determined value of the interaction coefficient at the same time. After this, the overall wettability of the steel–slag-inclusions is obtained. The inclusion motion model will be described in Section 2.3. The results of the motion behavior of the solid and liquid inclusions will be presented in the Results and Discussion sections. It should be noted that the overall studied 11 slag compositions are presented in Table 1, Table 2, Table 3 and Table 4.
The studied composition and viscosity of the CaO-Ce2O3-SiO2-Al2O3 slag system are listed in Table 1; the viscosity data are referred to from [133].
The composition of slag is converted to mole fraction and shown in Table 2.
The studied composition and viscosity of the CaO-Ce2O3-SiO2-Al2O3-MgO slag system are listed in Table 3; the viscosity data are referred from [87,88].
The composition of slag is converted to mole fraction and shown in Table 4.

2.1. Calculation Model of the Surface Tension of Slag

2.1.1. The Ion and Molecule Coexistence Theory

The ion and molecule coexistence theory of molten slag is based on the analysis of the structural units existing in molten slag. According to the equilibrium reaction equation and the equilibrium constant between the components, and from calculating the mass action concentration of each component, to determine the activity of each component in molten slag.
The ion and molecule coexistence theory describes the slag structure as [147]:
  • Slag is composed of various ions, compounds (molecules), and electrons. For the slag system studied in this paper, slag contains Ca2+, O2−, Ce2+, Mg2+, Al2O3, SiO2, and aluminate and silicate generated by their mutual reactions.
  • The coexistence of ions and molecules is continuous in the slag system of any composition.
  • Dynamic equilibrium reactions are achieved between ions and molecules.
  • The chemical reaction in slag obeys the law of mass action.
It is considered that this theory can be extended from inside of the slag to the surface of slag, so that the surface tension of the slag can be calculated [148]. Recently, many scholars have used this theory to carry out theoretical calculation research, and the results are in good agreement with the experimental results [149,150,151,152].

2.1.2. The Butler Equation

The Butler equation proposes that the surface of the solution is regarded as the surface phase, and that the relationship between the surface tension and the thermodynamic parameters is deduced based on the assumption that the components in the surface phase and the components in the bulk phase are in thermodynamic equilibrium [153,154]. This model calculates the surface tension of slag based on the Butler equation, which can be expressed as:
σ = σ i 0 + R T A i l n N i S u r f N i B u l k
where σi0 is the surface tension of pure component i, R is a gas constant, T is the temperature, and NiSurf and NiBulk are the mass action concentrations of component i in the surface phase and bulk phase, respectively. Ai is the surface area of the monolayer atoms of component i, which can be expressed as:
A i = L N A 1 / 3 V i 2 / 3
where L is the correction factor caused by the surface structure, which is related to the arrangement of monolayer atoms in the lattice. For a molten salt and ionic oxide mixture, it is usually set to 1.091. NA is Avogadro’s constant, NA = 6.02 × 1023; Vi is the molar volume of component i.
Taking a multi-components slag system as an example, the calculation steps of the model are:
  • The mass action concentration NiBulk of each component in a multi-components slag system are calculated by using the ion and molecule coexistence theory;
  • According to Equation (1), the expression of the mass action concentration NiSurf of each component on the surface of a multi-components slag system can be written;
  • For a combination of the surface tension and the molar volume data of each component, the equation ∑NiSurf = 1 and the expansion of the Butler equation, the mass action concentrations NiSurf of each component on the surface of a multi-components slag system are obtained;
  • The surface tension is acquired by solving the Butler equation.

2.1.3. Surface Tension Calculation of the CaO-Ce2O3-SiO2-Al2O3 Slag System

(1)
Mass Action Concentration Calculation of Each Component
Based on the ion and molecule coexistence theory and phase diagram, the CaO-Ce2O3-SiO2-Al2O3 slag system was analyzed. At a refining temperature of 1773~1873 K, the slag contains three kinds of ions Ca2+, O2−, Ce2+, two kinds of simple molecules, Al2O3 and SiO2, and 13 kinds of complex molecules CaO·Al2O3, CaO·2Al2O3, CaO·6Al2O3, 3CaO·Al2O3, 12CaO·7Al2O3, CaSiO3, Ca2SiO4, Ca3SiO5, Ce2O3·Al2O3, Ce2O3·11Al2O3, 2CaO·Al2O3·SiO2, CaO·Al2O3·2SiO2, and 3Al2O3·2SiO2 [130].
Defining b1 = ∑nCaO, b2 = ∑nCe2O3, a1 = ∑nSiO2 = 0.1, a2 = ∑nAl2O3, these are the mass fractions of CaO, Ce2O3, SiO2, and Al2O3, respectively. The mass action concentrations of each component are defined as follows: N1 = NCaO, N2 = NCe2O3, N3 = NSiO2, N4 = NAl2O3, N5 = NCaO·Al2O3, N6 = NCaO·2Al2O3, N7 = NCaO·6Al2O3, N8 = N3CaO·Al2O3, N9 = N12CaO·7Al2O3, N10 = NCaSiO3, N11 = NCa2SiO4, N12 = NCa3SiO5, N13 = NCe2O3·Al2O3, N14 = NCe2O3·11Al2O3, N15 = NCaO·Al2O3·2SiO2, N16 = N2CaO·Al2O3·SiO2, and N17 = N3Al2O3·2SiO2.
The chemical equilibriums in the temperature range of 1773~1873 K can be defined as follows [130]:
( Ca 2 + + O 2 ) + Al 2 O 3 = CaOAl 2 O 3
Δ G Θ = 18120 19.62 T , N 5 = K 1 N 1 N 4
( Ca 2 + + O 2 ) + 2 Al 2 O 3 = CaO 2 Al 2 O 3
Δ G Θ = 16400 26.8 T , N 6 = K 2 N 1 N 4 2
( Ca 2 + + O 2 - ) + 6 Al 2 O 3 = CaO 6 Al 2 O 3
Δ G Θ = 17430 37.2 T , N 7 = K 3 N 1 N 4 6
3 ( Ca 2 + + O 2 - ) + Al 2 O 3 = 3 CaO Al 2 O 3
Δ G Θ = 17000 32.0 T , N 8 = K 4 N 1 3 N 4
12 ( Ca 2 + + O 2 - ) + 7 Al 2 O 3 = 12 CaO 7 Al 2 O 3
Δ G Θ = 86100 205.1 T , N 9 = K 5 N 1 12 N 4 7
( Ca 2 + + O 2 - ) + SiO 2 = CaSiO 3
Δ G Θ = 81416 10.498 T , N 10 = K 6 N 1 N 3
2 ( Ca 2 + + O 2 - ) + SiO 2 = Ca 2 SiO 4
Δ G Θ = 160431 + 4.016 T , N 11 = K 7 N 1 2 N 3
3 ( Ca 2 + + O 2 - ) + SiO 2 = Ca 3 SiO 5
Δ G Θ = 93366 23.03 T , N 12 = K 8 N 1 3 N 3
( 2 Ce 3 + + 3 O 2 - ) + Al 2 O 3 = Ce 2 O 3 Al 2 O 3
Δ G Θ = 58555.5 15.04 T , N 13 = K 9 N 2 N 4
( 2 Ce 3 + + 3 O 2 - ) + 11 Al 2 O 3 = Ce 2 O 3 11 Al 2 O 3
Δ G Θ = 49323 80.56 T , N 14 = K 10 N 2 N 4 11
( Ca 2 + + O 2 - ) + Al 2 O 3 + 2 SiO 2 = CaO × Al 2 O 3 2 SiO 2
Δ G Θ = 13816.44 55.266 T , N 15 = K 11 N 1 N 4 N 3 2
2 ( Ca 2 + + O 2 - ) + Al 2 O 3 + SiO 2 = 2 CaO Al 2 O 3 SiO 2
Δ G Θ = 61964.64 60.29 T , N 16 = K 12 N 1 2 N 4 N 3
3 Al 2 O 3 + 2 SiO 2 = 3 Al 2 O 3 2 SiO 2
Δ G Θ = 4351 10.46 T , N 17 = K 13 N 4 3 N 3 2
The equilibrium constants Ki of all reactions satisfy the expression ΔGΘ = −R·T·lnKi; Ki can be calculated using the following formula:
K i = exp ( Δ G Θ R T )
where ΔGΘ is the change of the reaction-free energy of each reaction.
The total equilibrium mole number of all structural units in the CaO-Ce2O3-SiO2-Al2O3 slag system can be expressed as:
n i = 2 n 1 + 5 n 2 + n 3 + n 4 + n 5 + n 6 + n 7 + n 8 + n 9 + n 10 + n 11 + n 12 + n 13 + n 14 + n 15 + n 16 + n 17
where, ni corresponds to the mole number of 17 substances existing in the slag, Ni = ni/∑ni, but for the ion-pairs (2Ce3+ + 3O2−) and (Ca2+ + O2−) existing in the slag, NCe2O3 = NCe3+ + NO2− = (2nCe3+ + 3nO2−)/(∑ni ) = 5nCe2O3/∑ni, (Ca2+ + O2−) is calculated in the same way as (2Ce3+ + 3O2−).
According to mass equilibrium, the following formulas could be obtained:
b 1 = n i ( 0.5 N 1 + N 5 + N 6 + N 7 + 3 N 9 + N 10 + 2 N 11 + 3 N 12 + N 15 + 12 N 16 )
b 2 = n i ( 0.2 N 2 + N 13 + N 14 )
a 1 = n i ( N 3 + N 10 + N 11 + N 12 + 2 N 15 + N 16 + 2 N 17 )
a 2 = n i ( N 4 + N 5 + 2 N 6 + 6 N 7 + N 8 + 7 N 9 + N 13 + 11 N 14 + N 15 + N 16 + 3 N 17 )
N 1 + N 2 + N 3 + N 4 + N 5 + N 6 + N 7 + N 8 + N 9 + N 10 + N 11 + N 12 + N 13 + N 14 + N 15 + N 16 + N 17 = 1
Therefore, Equations (3)–(15) and (18)–(22) are the governing equations of the developed thermodynamic model for calculating the mass action concentrations Ni of the structural units or ion couples in the CaO-Ce2O3-SiO2-Al2O3 slag. Specifically, the values of N1, N2, N3, and N4 can be obtained by solving the Equation set (18)–(22), by taking Equations (3)–(15) and the specific values of T, b1, b2, and a2 as the initial conditions. Then, N5N17 can be solved by substituting N1N4 into Equations (3)–(15). Newton’s iteration method is used in the process of solving the equation set.
The calculated results of the mass action concentration (activity) and the surface tension are identical to the calculation results of Wu et al. [130], which have been validated with experimental results. The total average deviation of the model was 4–9% [154,155,156].
(2)
Surface tension calculation
Extending the ion and molecule coexistence theory from inside of the slag to the surface, the surface mass action concentrations are defined as follows: N1S = NSCaO, N2S = NSCe2O3, N3S = NSSiO2, N4S = NSAl2O3, N5S = NSCaO·Al2O3, N6S = NSCaO·2Al2O3, N7S = NSCaO·6Al2O3, N8S = NS3CaO·Al2O3, N9S = NS12CaO·7Al2O3, N10S = NSCaSiO3, N11S = NSCa2SiO4, N12S = NSCa3SiO5, N13S = NSCe2O3·Al2O3, N14S = NSCe2O3·11Al2O3, N15S = NSCaO·Al2O3·2SiO2, N16S = NS2CaO·Al2O3·SiO2, and N17S = NS3Al2O3·2SiO2.
According to the mass conservation:
N 1 S + N 2 S + N 3 S + N 4 S + N 5 S + N 6 S + N 7 S + N 8 S + N 9 S + N 10 S + N 11 S + N 12 S + N 13 S + N 14 S + N 15 S + N 16 S + N 17 S = 1
A new equation containing only N1S, N2S, N3S, and N4S is obtained by substituting the surface mass action concentration into the equilibrium constant expression in Equations (3)–(15), and then substituting the results into Equation (23).
According to the Equation (1), the Butler equation corresponding to the CaO-Ce2O3-SiO2-Al2O3 slag system can be written in the following form:
σ C a O C e 2 O 3 S i O 2 A l 2 O 3 = σ C a O 0 + R T A C a O l n N 1 s N 1
= σ C e 2 O 3 0 + R T A C e 2 O 3 l n N 2 s N 2
= σ S i O 2 0 + R T A S i O 2 l n N 3 s N 3
= σ A l 2 O 3 0 + R T A A l 2 O 3 l n N 4 s N 4
Three new equations are obtained by subtracting (25)–(27) from (24), respectively, and the equation simplified by (23) is formed into a set of equations. The values of N1S, N2S, N3S, and N4S can be gained by solving the equations. The surface tension of the CaO-Ce2O3-SiO2-Al2O3 slag system can be acquired by substituting the values of N1S, N2S, N3S, and N4S into one of the Equations of (24)–(27).
To calculate the surface tension of this slag system, the surface tension and molar volume data of each component are also required, which can be found in references [109,148], and are listed in Table 5.
The surface tension and molar volume of Ce2O3 cannot be found in the available literatures. However, the data of rare earth lanthanide elements can be found in references [134,147], and they are listed in Table 6.
According to the data in Table 6, the relationship between the molar volume (m3/mol) of Ce2O3 and temperature (K) is reasonably estimated to be 49.0 [1 + 1·10−4·(T-1773)]·10−6, and the surface tension of Ce2O3 at 1873 K is estimated to be 680 × 10−3 N/m. The surface tension at 1773 K is 690 × 10−3 N/m.
After calculating via MATLAB software, the surface tension data of the CaO-Ce2O3-SiO2-Al2O3 slag system are obtained and listed in Table 7. It should be noted that the surface tension varies within a limited range.

2.1.4. Surface Tension Calculation of the CaO-Ce2O3-SiO2-Al2O3-MgO Slag System

Compared with the CaO-Ce2O3-SiO2-Al2O3 slag system, the CaO-Ce2O3-SiO2-Al2O3-MgO slag system increases by one more ion Mg2+, one more simple molecule MgO, and five more complex molecules MgO·Al2O3, MgO·SiO2, Mg2SiO4, CaO·MgO·2SiO2, and 2CaO·MgO·2SiO2. The mass action concentration of each component was defined as follows: N18 = NMgO, N19 = NMgO·Al2O3, N20 = NMgO·SiO2, N21 = NMg2SiO4, N22 = NCaO·MgO·2SiO2, and N23 = N2CaO·MgO·2SiO2. The following five equilibrium reactions are required:
( Mg 2 + + O 2 - ) + Al 2 O 3 = MgO Al 2 O 3
Δ G Θ = 35530 17.39 T , N 19 = K 14 N 18 N 4
( Mg 2 + + O 2 - ) + SiO 2 = MgO SiO 2
Δ G Θ = 41089.4 6.10 T , N 20 = K 15 N 18 N 3
2 ( Mg 2 + + O 2 - ) + SiO 2 = 2 MgO SiO 2
Δ G Θ = 67130.8 2.508 T , N 21 = K 16 N 18 2 N 3
( Ca 2 + + O 2 - ) + ( Mg 2 + + O 2 - ) + 2 SiO 2 = CaO MgO 2 SiO 2
Δ G Θ = 162602 + 18.81 T , N 22 = K 17 N 1 N 18 N 3 2
2 ( Ca 2 + + O 2 - ) + ( Mg 2 + + O 2 - ) + 2 SiO 2 = 2 CaO MgO 2 SiO 2
Δ G Θ = 73688 63.639 T , N 23 = N 1 2 N 18 N 3 2
The surface tension of the CaO-Ce2O3-SiO2-Al2O3-MgO slag system is calculated by using the existing data and calculation methods described in Section 2.1.3. The results are shown in Table 8.

2.2. Calculation Model of Interfacial Tension

According to Girifalco-Good equation, the interfacial tension between inclusions and molten steel, between inclusions and slag, and between slag and steel are calculated.
Girifalco-Good equation [157]:
σ a b = σ a + σ b 2 φ a b σ a σ b
where the subscripts a and b represent the two phases, σab is the interfacial tension between phases a and b, σa and σb are the surface tensions of each phase, and φab represents the interaction coefficient between a and b.

2.2.1. Interfacial Tension between Inclusion and Molten Steel

This section mainly refers to the research results of Nakajima [158]. According to the Girifalco-Good equation, the interfacial tension between inclusions and molten steel can be expressed as:
σ I M = σ I + σ M 2 φ I M σ I σ M
The surface tension σI of Al2O3 inclusions is 0.75 N/m. To calculate the surface tension σM of molten steel, as a simplification, only a certain content of O and S are considered in steel. Then, the surface tension of liquid steel is expressed as follows:
σ M = [ 1910 825 l o g 10 ( 1 + 210 [ O ] ) 540 l o g 10 ( 1 + 185 [ S ] ) ] 10 3 ( N / m )
( [ S ] 0.18 w t % , [ O ] 0.015 w t % )
In this study, [S] = 0.005%, [O] = 0.0025%, and σM = 1.6052 N/m.
The interaction coefficient between solid inclusions and molten steel:
φ I M = 0.351 0.507 ( R 1 )
where R is the surface roughness of the inclusion, and its value ranges over 1.0~1.36.
Ideally, the surface roughness of the inclusion is 1, so that φIM = 0.351. Substituting all the data into Equation (35), we can obtain σIM = 1.585 N/m.

2.2.2. Interfacial Tension between Inclusion and Slag

The interfacial tension between inclusion and slag can be expressed as follows:
σ I S = σ I + σ S 2 φ I S σ I σ S
The surface tension between the inclusion and slag has been obtained. The interaction coefficient φIS between the inclusion and slag is required.
The interaction coefficient between the Al2O3 inclusion and the CaO-SiO2-Al2O3-MgO-CaF2-Na2O slag system was studied in reference [158]; φIS is expressed as:
φ I S = 0.947 × X A l 2 O 3 + 0 × X S i O 2 + 0.602 × X C a F 2 + 0.601 × X M g O + 0.545 × X C a O + 0.475 × X N a 2 O + 0.705 × X A l 2 O 3 × X C a O 0.609 × X A l 2 O 3 × X S i O 2 + 1.564 × X S i O 2 × X C a O + 0 × X C a F 2 × X C a O + 0.755 × X C a F 2 × X S i O 2 + 1.292 × X S i O 2 × X M g O + 1.156 × X S i O 2 × X N a 2 O
where XCaO, XSiO2, XAl2O3, XMgO, XCaF2, and XNa2O are the molar fractions of each component.
The interaction coefficient between the slag system and the Al2O3 inclusion was studied in the research results of Nakajima [158], as shown in Table 9.
According to the data in the literature and theoretical analysis, φIS = 0.8 is estimated for the CaO-Ce2O3-SiO2-Al2O3 slag system.
For the CaO-Ce2O3-SiO2-Al2O3-MgO slag system, the interaction coefficient is calculated by Equation (38), ignoring the effect of Ce2O3. The calculation results are shown in Table 10.
Based on the calculated data, considering the stronger interaction between slag containing Ce2O3 and Al2O3, φIS = 0.75 is used for the inclusion–slag interfacial tension of slag C0, and φIS = 0.78 is used for the inclusion–slag interfacial tension of slag C1–C5.
By substituting σI, σS, and φIS of slags with different compositions into Equation (37), σIS is calculated.

2.2.3. Interfacial Tension between Slag and Steel

The interfacial tension between slag and steel can be expressed as:
σ M S = σ M + σ S 2 φ M S σ M σ S
For the CaO-SiO2-Al2O3-MgO-CaF2-Na2O slag system studied in reference [158], the interaction coefficient between slag and steel is expressed as:
φ M S = 0.576 × X S i O 2 + 0.399 × X A l 2 O 3 + 0.326 × X C a F 2 + 0.687 × X M g O + 0.607 × X C a O 0.675 × X S i O 2 × X C a O + 0.338 × X S i O 2 × X A l 2 O 3 0.391 × X A l 2 O 3 × X C a O 0.892 × X A l 2 O 3 × X C a F 2 1.701 × X M g O × X S i O 2 0.27 × X M g O × X A l 2 O 3
The relevant data of Ce2O3 cannot be found in the literature, so that the interaction coefficient between the CaO-Ce2O3-SiO2-Al2O3 slag system and molten steel needs to be reasonably estimated. The CaO-SiO2-Al2O3-MgO slag system is referred, and the interaction coefficient between slag and steel in reference [158] is listed in Table 11.
It is observed that the interaction coefficient between the CaO-Al2O3-SiO2-MgO slag system, which is similar to the CaO-Ce2O3-SiO2-Al2O3 slag system studied in this paper, and steel, ranges from 0.428 to 0.454. The average value of the interaction coefficient between CaO-Ce2O3-SiO2-Al2O3 slag and liquid steel is about 0.43. In order to improve the accuracy of the interaction coefficient, it is modified according to the following analysis.
In Section 2.3, the motion model of inclusion at the slag–steel interface, the overall wettability cosθIMS of the steel–slag inclusion system is required:
c o s θ I M S = σ I M σ I S σ M S
The relationship between σIM, σIS, σMS, and θIMS is shown in Figure 2. There are two wetting states between slag and inclusion, as shown in Figure 2a; cosθIMS > 0, which means the inclusion shows a good wettability for slag. If cosθIMS <0, as shown in Figure 2b, this means that the inclusion shows a poor wettability for slag.
It can be seen that cosθIMS increases with the decrease in σIS or σMS, so that the slag–steel interaction coefficient should be reduced to a smaller value to compensate for the deficiency where σIS may be too large, which is beneficial for improving the accuracy of the calculation model. Therefore, for CaO-Ce2O3-SiO2-Al2O3 slag system, φMS = 0.4 is used in the calculation.
By substituting σM, σS, and φMS into Equation (39), the interfacial tension σMS between the CaO-Ce2O3-SiO2-Al2O3 slag system and molten steel can be obtained.
For the CaO-Ce2O3-SiO2-Al2O3-MgO slag system, the slag composition is close to the slags in Table 11. Thus, the interaction coefficient between slag and steel is referred to in Table 11. Therefore, φMS = 0.43 is used in the calculation of the steel–slag interfacial tension of slag C0 without Ce2O3. Considering that the interaction between slag containing Ce2O3 and molten steel will be strengthened, φMS = 0.45 is estimated in the calculation of slag C1–C5. In the same way, the interfacial tension σMS between the CaO-Ce2O3-SiO2-Al2O3-MgO slag system and molten steel can be obtained.

2.2.4. Interfacial Tension Summary

All of the calculation results of the two slag systems are summarized in Table 12 and Table 13.

2.3. Motion Model of Inclusion at the Steel–Slag Interface

The interfacial tension data calculated in the above sections will be used in the motion model of inclusions that will be described in this section.

2.3.1. Assumptions of This Model

The basic assumptions of this mathematical motion model of inclusions at the steel–slag interface are listed as follows:
(1) The inclusions are spherical in constant volume;
(2) No chemical reaction occurs between the phase interfaces;
(3) All fluids are incompressible and isothermal;
(4) The slag phase exists in liquid form;
(5) The steel–slag interface is smooth;
(6) The movement process of inclusions depends on the force balance between the buoyancy force, rebound force, drag force, and the fluid added mass force;
(7) The surface tension is uniform at the interface.
As proposed by Nakajima and Okamura [71], when the inclusion approaches the slag–steel interface, a steel film forms when Re ≥ 1, otherwise, the inclusion contacts with the slag directly without steel film formation. The Reynolds number of inclusions is defined as:
R e = 2 ρ M u R I μ M
where RI is the radius of inclusion, ρM is the density of molten steel, μM is the viscosity of the molten steel, and u is the velocity of inclusion.
In the studied small size of the inclusion cases, for example, if the diameter is less than 150 μm, the Reynolds number is less than 1. Thus, the model of non-steel film could be used. In this paper, the motion behavior of the inclusion is studied at the steel–slag interface. Specifically, the inclusion is at an initial position below the slag; that is, when t = 0, the displacement from the center of inclusion to the interface Z = 0, as shown in Figure 3. Then, there are three kinds of interfacial motion behaviors of the inclusions: Pass (P), Oscillate (O), and Remain (R). “Pass” means that the displacement of the inclusion is higher than two times the radius, the motion path of inclusions is 1→2→3, as shown in Figure 3a. “Oscillate” means that the inclusions oscillate but do not penetrate the slag–steel interface, and the motion path of the inclusions is 1→2→3→4, as shown in Figure 3c. “Remain” means that the displacement of inclusion is less than two times the radius, and the motion path of the inclusions is 1→2→3, as shown in Figure 3b. When it moves up by two times their radius 2RI, as shown in Figure 3a, it is considered to be completely separated from the molten steel into the slag.
The initial floating velocity of inclusions is expressed by the Stokes velocity, and the Stokes velocity of liquid inclusions (such as 50%wtAl2O3–50%wtCaO) can be expressed as:
u = 2 3 R I 2 ( ρ M ρ I ) g μ M ( 1 + κ ) ( 2 + 3 κ )
where ρI is the density of inclusions; g is the acceleration of gravity; κ = μIM, μI is the viscosity of inclusion, and μM is the viscosity of molten steel.
In this study, ρM = 7000 kg/m3, μM = 0.006 Pa·s. ρI = 2814 kg/m3, μI = 0.118 Pa·s, and κ = μIM =19.67 for liquid inclusion 50%wtAl2O3–50%wtCaO.
For solid inclusions (such as Al2O3), κ→∞, ρI = 3990 kg/m3, the Stokes velocity of this inclusion can be expressed as:
u = 2 9 R I 2 ( ρ M ρ I ) g μ M

2.3.2. Motion Equation of Inclusion

When the inclusions move across the slag–steel interface, four forces act on the inclusions, i.e., the buoyancy force Fb, rebound force Fr, drag force Fd, and fluid added mass force Ff.
The buoyancy force is in an upward direction, while the rebound force, fluid added mass force, and drag force are either upward or downward, depending on the behavior of the inclusions at the interface.
According to Newton’s second law:
4 3 π R I 3 ρ I d 2 Z d t 2 = F b F d F r F f
where Z is the displacement of the inclusion and t is the movement time of the inclusion.
The rebound force is defined as:
F r = 2 π R I σ M S H ( Z * )
where Z* is the dimensionless displacement of the inclusion and H is a function of Z*.
H ( Z ) = Z 1 cos θ I M S
The buoyancy force is defined as:
F b = 4 3 π R I 3 g ( ρ S J ( Z * ) ρ I )
where J(Z*) is the density-varying term, which describes the degree of the inclusion’s entering into the slag.
J ( Z * ) = 1 4 ( ρ M ρ S 1 ) Z * 3 3 4 ( ρ M ρ S 1 ) Z * 2 + ρ M ρ S
The drag force is defined as:
F d = 4 π g R I 3 μ S A I ( Z * ) d Z * d t *
where t* is the dimensionless movement time of the inclusion, A is a function of κ, which is related to the viscosity of inclusions and molten steel; I(Z*) is the viscosity-varying term, which is related to the position of the inclusion relative to the interface.
A = 2 + 3 κ 2 ( 1 + κ )
I ( Z * ) = ( μ M μ S 1 ) Z * 2 2 ( μ M μ S 1 ) Z * + μ M μ S
The fluid added mass force is defined as:
F f = 2 3 π R I 3 g ρ S J ( Z * ) d 2 Z * d t * 2
Substituting Equations (46), (48), (50) and (53) into (45), the dimensionless motion equation of inclusion can be written as follows:
d 2 Z * d t * 2 = 2 ( ρ S J ( Z * ) ρ I ) ( ρ S J ( Z * ) + 2 ρ I ) 3 * D ( Z * ) * H ( Z * ) 6 A E ( Z * ) I ( Z * ) * d Z * d t *
where both D and E are functions of Z*.
D ( Z * ) = σ M S g R I 2 ( ρ S J ( Z * ) + 2 ρ I )
E ( Z * ) = g R I 3 ( ρ S J ( Z * ) + 2 ρ I ) μ S
The dimensionless form of displacement, time, velocity, and acceleration are defined as:
Z = Z R I
t * = t g R I
d Z * d t * = 1 g R I d Z d t
d 2 Z * d t * 2 = 1 g d 2 Z d t 2
The densities of slag ρs are 2700 kg/m3 and 2750 kg/m3 for the CaO-Ce2O3-SiO2-Al2O3 slag system and the CaO-Ce2O3-SiO2-Al2O3-MgO slag system, respectively.
By solving Equation (54), the displacement, velocity, and acceleration changes of the inclusions with time can be solved, and the force changes of the inclusions can be calculated by using the values of these parameters.
The methods for solving the force balance and the motion equation of inclusions are as follows: the Equations (54) and (46), (48), (50) and (53) are the governing equations of the developed mathematical model for calculating the dimensionless displacement Z*(t) and the force on inclusions. Specifically, by initializing the inclusion size, the governing Equation (54) of Z*(t) is solved using the ode45 function (Runge kutta algorithm) in MATLAB software with input of the calculated physical properties, as mentioned in previous sections. After the dimensionless displacement is obtained, the force can be calculated using Equations (46), (48), (50) and (53) with specific sizes of inclusions. This model is originally validated using experimental data by Nakajima and Okamura [71]. Their modeling work on inclusions movement in slags were reproduced, and the results are in good agreement with the original paper [71]. This ensures the reliability of the present simulation.

2.3.3. Model Validation

Based on the similarity theory, a water model experiment was established to verify the mathematical model. A schematic diagram of the experimental equipment is shown in Figure 4. In the experiment, silicon oil was used to simulate slag, water was used to simulate molten steel, and hollow alumina balls with 2.1 mm diameter were used to simulate inclusions. A cylindrical container with a particle release port at the bottom was filled with water, and the height of the water was 40 cm. In order to eliminate the refraction generated by the cylindrical camber, a square container was added around the cylindrical container, and the water level was higher than 40 cm. The water was covered by silicon oil with a thickness of more than 5 mm. Hollow alumina balls were released from the bottom release port stationary. A high-speed camera (FuHuang AgileDevice Revealer 5KF20S) was used to record the frame of the alumina balls touching the water–oil interface until the motion stopped.
The dimensionless displacement curve of the alumina ball moving at the water–oil interface could be obtained by processing the video shot using the high-speed camera. Meanwhile, the initial conditions of the water model in Table 14 were substituted into the mathematical model for calculation [159]. The comparison between the calculated results and the experimental results is shown in Figure 5. It can be seen from the figure that the overall trend of the calculated value is consistent with the experimental value. After the particle reaches the water–oil interface, it floats up for a certain distance and is bounced back. The final dimensionless displacements for the experiment and model are 0.594 and 0.655, respectively. It is in good agreement, which verifies the feasibility and accuracy of the model.

3. Results and Discussion

In this section, the motion behavior of solid Al2O3 inclusions and liquid inclusions in two slag systems are presented. The evolution of force on inclusions, and the displacement of inclusions as a function of movement time are studied. By identifying the displacement of inclusion, the motion behavior of the inclusions can be categorized as Remain, Oscillate, and Pass. Using a large number of studied cases, the force and displacement results of two typical slag systems, A5 and C3, with respect to different sizes of inclusions, are given as examples. Finally, the inclusion removal abilities of the two slag systems are compared.

3.1. Motion Behavior of Solid Al2O3 Inclusions in the CaO-Ce2O3-SiO2-Al2O3 Slag System

The force of the solid Al2O3 inclusions with different sizes in molten steel are shown in Figure 6. The maximum value of the t axis is the time taken for the inclusions to reach force equilibrium. The force equilibrium for the inclusions is a function of moving time with different sizes. As can be seen from the figure, with the increase in the inclusion size, the time for inclusions to reach the force equilibrium in the process of motion is increasing. For all of the inclusions, the rebound force and drag force play a major role during the movement of the inclusion. With the increase in the size of inclusion, the rebound force, drag force, and fluid added mass force are increased.
The displacements of solid Al2O3 inclusions with different sizes at the steel–slag interface in slag A5 are shown in Figure 7.
According to the dimensionless displacement shown in Figure 7 and the classifications of Pass, Oscillate, and Remain, the inclusions with radiuses of 50 μm, 80 μm, and 140 μm can be categorized as the “Pass” group, and the time of inclusions passing through the slag–steel interface is prolonged with the increase in inclusion size. The inclusions with a radius of 10 μm and 20 μm can be categorized as the “Oscillate” group, and the inclusion with a radius of 5 μm can be categorized as the “Remain” group.
The motion behaviors of inclusions with different sizes in slags A1–A5 were calculated, and the results are shown in Table 15.
From Table 15, the smaller the size of the inclusion is, the more difficult it is to remove. The inclusions with a radius that is smaller than 20 μm cannot be removed in all the A series slag systems, and the inclusions with a radius of larger than 50 μm can be removed perfectly. The inclusion removal abilities of the five slags are ranked as A5 > A2 > A3 > A1 > A4. It should be noted that the viscosity of A5 is the smallest among the slags A1 to A5, indicating that the reduction in viscosity is more conducive to improving the inclusion removal abilities of the slags.

3.2. Motion Behaviors of Liquid 50%wtAl2O3–50%wtCaO Inclusions in the CaO-Ce2O3-SiO2-Al2O3 Slag System

The force of liquid 50%wtAl2O3-50%wtCaO inclusion in molten steel is basically similar to that of solid Al2O3, and it is not shown here. The displacements of the liquid 50%wtAl2O3–50%wtCaO inclusion at the steel–slag interface with different sizes in slag A5 are shown in Figure 8. The motion behaviors of liquid 50%wtAl2O3–50%wtCaO inclusions are similar to those of solid Al2O3 inclusions. Specifically, the inclusions with radiuses of 5 μm and 10 μm show Remain behaviors at the steel–slag interface. The 20 μm-radius inclusion exhibits oscillatory behavior. The inclusions with radiuses of 50 μm, 80 μm, and 130 μm can pass the interface and can be removed, but the time for liquid inclusions passing through the slag–steel interface is longer than for solid inclusions under the same slag condition.
The motion behaviors of liquid 50%wtAl2O3–50%wtCaO inclusions with different sizes in slags A1–A5 were calculated, and the results are shown in Table 16.
For liquid inclusions, the removal abilities in slags A1 to A5 are not as good as those of solid inclusions. The inclusions with a radius smaller than 50 μm are difficult to remove, and the inclusions with a radius larger than 80 μm can be removed in some of the slags. Similar to the results of solid inclusions, the inclusion removal abilities of the five slags are ranked as A5 > A2 > A3 > A1 > A4; the lower the viscosity of the refining slag is, the better the removal of inclusions.
Compared with Table 15 and Table 16, from the perspective of the separation behavior of inclusions at the steel–slag interface, the solid Al2O3 inclusion with a radius of 20 μm showed Remain behavior in slag A4, whose viscosity is the highest of the five slags, and it showed Oscillate behavior in other slags. For the liquid inclusion with the same radius of 20 μm, it showed Oscillate behavior merely in slag A5, which posed the lowest viscosity among the five slags. In addition, the solid Al2O3 inclusion with a radius of 80 μm can pass the interface of all five slags, while the liquid inclusion with the same radius can only pass the interfaces of slags A2 and A5. Obviously, the size range in which the solid inclusions can pass the interface is larger than that of liquid inclusions. Compared with Figure 7 and Figure 8, it is easy to find that the removal time of liquid inclusions is longer than that of solid inclusions. All of these are indicating that liquid 50%wtAl2O3–50%wtCaO inclusions are more difficult to remove than solid Al2O3 inclusions.

3.3. Motion Behaviors of Solid Al2O3 Inclusions in the CaO-Ce2O3-SiO2-Al2O3-MgO Slag System

The displacement of solid Al2O3 inclusions with different sizes at the steel–slag interface in slag C3 is shown in Figure 9. As can be seen from Figure 9, the inclusions with a radius of 20 μm, 50 μm, 80 μm, and 140 μm can be removed from the slag. The inclusions with a radius of 5 μm and 10 μm showed oscillatory behavior. Compared with Figure 7 and Figure 9, the solid inclusions removal ability of slag C3 is improved compared to that of slag A5. Furthermore, the removal time is shorter.
Similarly, the motion behaviors of inclusions in slag C0–C5 were calculated, and the results are shown in Table 17.
The inclusions with a radius of smaller than 10 μm cannot be removed in all of the C series slag systems, and the inclusions with a radius of larger than 20 μm can be removed perfectly in slags C3 to C5. This C series slags showed a better performance than the A series slags for solid inclusion removal.
By comparing the removal of the Al2O3 inclusions between slag C0 and slags C1–C5, it can be seen that the refining slag containing rare earth oxide Ce2O3 shows a better ability to remove inclusions than ordinary refining slag. The inclusion removal ability of the six slags is ranked as C3 > C5 > C4 > C1 > C2 > C0. By comparing the physical properties of slag C0 and slags C1–C5, it can be found that slags C1–C5 have a higher wettability, as well as a better performance of the wetting effect on the Al2O3 inclusion than slag C0. However, the overall wettability varies within a limited range. As shown in Table 3, the viscosity of the six slags is ranked as C3 < C5 < C4 < C1 < C2 < C0. This is the same order as the removal ability of the slags. Thus, the viscosity may be a crucial factor. This is the reason for why the slag containing rare earth oxide Ce2O3 shows a better ability to remove Al2O3 inclusions.
Compared with Table 15 and Table 17, the removal ability of solid Al2O3 inclusions in the C series slag system is significantly higher than that in the A series slag system. The limiting size of the inclusions that can be removed is reduced in the C series slag system.

3.4. Discussion

According to the analysis above, the ability of slags containing Ce2O3 to remove inclusions is better than the slag without Ce2O3, and the viscosity plays a vital role in the ability of the two slag systems to remove inclusions. In the C series slags, the addition of MgO can effectively reduce the viscosity. The sizes of the inclusions that can effectually be removed is related to the statement of the inclusion. In this paper, for solid inclusions, it can be removed with a size of more than 20 μm in a special refining slag with a specific composition. For liquid inclusions, the size of the inclusion needs to be larger than 80 μm. That is to say, the liquid inclusions more easily remain in the molten steel. After the analysis, the ability of the refining slag to remove inclusions of specific compositions can be evaluated.
Nevertheless, the model is used to predict the motion of single particle inclusions at the slag–steel interface in a relatively ideal state, but there are some shortcomings. First, the liquid steel and slag are assumed to be static, without considering the real flow state. Second, the inclusion is a single particle floating at the interface, without considering the collision aggregation between the inclusions. Third, the inclusion may be accelerated and the rising velocity may be higher than the Stokes velocity. Although the results cannot be compared with industrial operations, the integrated model can predict a tendency for the different compositions of the slags. For a specific slag, the composition of the slag could influence the viscosity and surface tension, and further, the overall wettability, and many other issues. In this study, the C series slag with a lower viscosity shows the best performance of inclusion removal. The 45%wt CaO–35.5%wt Al2O3–4.5%wt SiO2–10%wt MgO–5%wt Ce2O3 refining slag is optimized in this study. It is noted that a systematic evaluation of the inclusion removal ability of selected slag systems, as well as the effect of the tailored slag on slag eye in ladles [160], are ongoing work.

4. Conclusions

After summarizing and analyzing the motion behaviors of solid and liquid inclusions with different sizes at the steel–slag interface of the two slag systems, the following conclusions can be obtained:
(1) In the refining process, whether solid or liquid inclusions, large inclusions are easier to remove from the steel–slag interface than small inclusions. Reducing the viscosity of the slag is more conducive to the removal of the inclusions. For the same CaO-Ce2O3-SiO2-Al2O3 slag system, liquid 50%wtAl2O3-50%wtCaO inclusions are more difficult to remove than solid Al2O3 inclusions.
(2) The CaO-Ce2O3-SiO2-Al2O3-MgO refining slag shows a better ability to remove Al2O3 inclusions than that of CaO-SiO2-Al2O3-MgO slag. The reason for this is that the addition of rare earth oxide Ce2O3 can decrease the viscosity of slags, as well as improving the wetting effect of slags on Al2O3 inclusions.
(3) For two slags systems, the CaO-Ce2O3-SiO2-Al2O3-MgO slag system shows a better ability to remove Al2O3 inclusions than that of the CaO-Ce2O3-SiO2-Al2O3 slag system. The addition of 5% to 8% Ce2O3 in the CaO-SiO2-Al2O3-MgO slag is an optimized case for industrial applications.
(4) The integrated model mainly involves an inclusion motion model and the slag properties models, including the interfacial tension, surface tension, and mass action concentration model, based on ion and molecule coexistence theory. However, the predictions are limited to the force balance of inclusions. A combined model considering the transport of inclusions by fluid mechanics, the interaction of inclusions with a turbulent boundary layer, and the movement of inclusions at the steel–slag interface will be a future work.

Author Contributions

Conceptualization, C.C.; methodology, W.L., Y.T., J.D. and C.C.; software, J.C. (Jianqi Cao), Y.L. and C.C.; investigation J.C. (Jianqi Cao), Y.L., F.Z. and C.C.; resources, W.L., Y.T., J.D. and C.C.; data curation, J.C. (Jianqi Cao) and Y.L.; writing—original draft preparation, J.C. (Jianqi Cao) and Y.L.; writing—review and editing, W.L., J.C. (Julong Che), F.Z., Y.T., D.L., J.D. and C.C.; supervision, W.L., J.C. (Julong Che) and C.C.; funding acquisition, C.C. and D.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 51904204; Research Project Supported by Shanxi Scholarship Council of China, grant number 2022-040; Lyuliang High level Talents Key R&D Project, grant number 2019112; and Scientific and Technological Innovation Programs of Higher Education Institutions in Shanxi, grant numbers 2019L0971, 2019L0938.

Data Availability Statement

Not applicable.

Acknowledgments

The discussions with Professor Keiji Nakajima of KTH-Royal Institute of Technology and Professor Emeritus Shebin Wang of Taiyuan University of Technology are acknowledged. The anonymous reviewers are acknowledged for their valuable comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kiessling, R. Non-Metallic Inclusions in Steel: Part III: The Origin and Behavior of Inclusions and Their Influence on the Properties of Steels; The Metals Society: London, UK, 1968. [Google Scholar]
  2. Zhang, L.F. Non-Metallic Inclusions in Steels: Fundamentals; Metallurgical Industry Press: Beijing, China, 2019. [Google Scholar]
  3. Park, J.H.; Zhang, L.F. Kinetic Modeling of Nonmetallic Inclusions Behavior in Molten Steel: A Review. Metall. Mater. Trans. B 2020, 51, 2453–2482. [Google Scholar] [CrossRef]
  4. Nadif, M.; Lehmann, J.; Burty, M.; Domgin, J.F. Control of Steel Reoxidation and CC Nozzle Clogging: An overview. In Proceedings of the International Conference on Clean Steel 7, Balatonfüred, Hungary, 4–6 June 2007. [Google Scholar]
  5. Huang, R.K.; Zhang, L.F.; Jiang, R.B.; Guo, Y.T.; Zhou, Q.Y.; Duan, H.J. Effect of FeOx in an ultra-low carbon Al-killed steel on nozzle clogging. Iron Steel 2021, 56, 43–50. [Google Scholar] [CrossRef]
  6. Hua, C.J.; Wang, M.; Bao, Y.P. Effect of Nozzle Clogging on the Fluid Flow Pattern in a Billet Mold with Particle Image Velocimetry Technology. Metall. Mater. Trans. B 2020, 51, 2871–2881. [Google Scholar] [CrossRef]
  7. Fu, B.H.; Chen, C.; Cheng, G.G.; Pan, J.X.; Li, Y.; Pan, W. Inclusions in 430 Stainless Steelmaking During AOD-LF-CC Process. Iron Steel 2012, 47, 40–43. [Google Scholar] [CrossRef]
  8. Chen, C.; Cheng, G.G. Delta-ferrite distribution in a continuous casting slab of Fe-Cr-Mn austenitic stainless steel. Metall. Mater. Trans. B 2017, 48, 2324–2333. [Google Scholar] [CrossRef]
  9. Li, Y.H.; He, Y.B.; Ren, Z.F.; Bao, Y.P.; Wang, R. Comparative Study of the Cleanliness of Interstitial-Free Steel with Low and High Phosphorus Contents. Steel Res. Int. 2021, 92, 2000581. [Google Scholar] [CrossRef]
  10. Wang, Y.Z.; Jiang, M.; Wang, X.H. Formation and evolution of inclusions in Q345D steel during secondary reefing process. Iron Steel 2022, 57, 63–72. [Google Scholar] [CrossRef]
  11. Wang, K.P.; Wang, Y.; Xu, J.F.; Chen, T.J.; Xie, W.; Jiang, M. Investigation on evolution of inclusions in bearing steel during secondary refining. Iron Steel 2022, 57, 42–49. [Google Scholar] [CrossRef]
  12. Cao, J.Q.; Chen, C.; Xue, L.Q.; Li, Z.; Zhao, J.; Lin, W.M. Analysis of the inclusions in the whole smelting process of non-oriented silicon steel DG47A. Iron Steel 2022. [Google Scholar] [CrossRef]
  13. Yuan, B.H.; Liu, J.H.; Zeng, J.H.; Zhang, M.; Huang, J.H.; Yang, X.D. Evolution of Inclusions and Cleanliness in Ti-Bearing IF Steel Produced via the BOF–LF–RH–CC Process. Metals 2022, 12, 434. [Google Scholar] [CrossRef]
  14. Yuan, X.B.; Zhong, M.; Wu, Y.W.; Wang, C. Characterizing Inclusions in the Weld Metal of Eh36 Shipbuilding Steel Processed by CaF2-30 Wt Pct TiO2 Flux. Metall. Mater. Trans. B 2022, 53, 656–661. [Google Scholar] [CrossRef]
  15. Qiao, T.; Cheng, G.G.; Huang, Y.; Li, Y.; Zhang, Y.L.; Li, Z.C. Formation and Removal Mechanism of Nonmetallic Inclusions in 42CrMo4 Steel during the Steelmaking Process. Metals 2022, 12, 1505. [Google Scholar] [CrossRef]
  16. Jiang, M.; Liu, J.C.; Li, K.L.; Wang, R.G.; Wang, X.H. Formation Mechanism of Large CaO-SiO2-Al2O3 Inclusions in Si-Deoxidized Spring Steel Refined by Low Basicity Slag. Metall. Mater. Trans. B 2021, 52, 1950–1954. [Google Scholar] [CrossRef]
  17. Podder, A.; Coley, K.S.; Phillion, A.B. Modeling Study of Steel–Slag–Inclusion Reactions During the Refining of Si–Mn Killed Steel. Steel Res. Int. 2022, 2100831. [Google Scholar] [CrossRef]
  18. Meng, Y.Q.; Li, J.L.; Zhu, H.Y.; Wang, K.P. Effect of LF soft bubbling time on oxide inclusions in Si-killed spring steel. Iron Steel 2022, 57, 48–54. [Google Scholar] [CrossRef]
  19. Zhang, Y.X.; Chen, C.; Lin, W.M.; Yu, Y.C.; E, D.Y.; Wang, S.B. Numerical Simulation of Tracers Transport Process in Water Model of a Vacuum Refining Unit: Single Snorkel Refining Furnace. Steel Res. Int. 2020, 91, 2000022. [Google Scholar] [CrossRef]
  20. Ouyang, X.; Lin, W.M.; Luo, Y.; Zhang, Y.Z.; Fan, J.P.; Chen, C.; Cheng, G.G. Effect of Salt Tracer Dosages on the Mixing Process in the Water Model of a Single Snorkel Refining Furnace. Metals 2022, 12, 1948. [Google Scholar] [CrossRef]
  21. Cheng, G.; Zhang, L.F.; Ren, Y.; Yang, W.; Zhao, X.L.; Wang, G.C. Evolution of Nonmetallic Inclusions with Varied Argon Stirring Condition during Vacuum Degassing Refining of a Bearing Steel. Steel Res. Int. 2021, 92, 2000364. [Google Scholar] [CrossRef]
  22. Yang, R.G.; Yan, Z.H.; Cao, C.G.; Shi, S.D.; Wang, G.H.; Liu, J.G. Practice of VD furnace oxygen decarburization process. Iron Steel 2021, 56, 39–43. [Google Scholar] [CrossRef]
  23. Li, Y.H.; Zhu, H.W.; Wang, R.; Ren, Z.F.; Lin, L. Prediction of two phase flow behavior and mixing degree of liquid steel under reduced pressure. Vacuum 2021, 192, 110480. [Google Scholar] [CrossRef]
  24. Wei, G.S.; Dong, J.F.; Zhu, R.; Han, B.C. Effect of ladle bottom blowing on RH dehydrogenation and inclusions. Iron Steel 2021, 56, 63–68. [Google Scholar] [CrossRef]
  25. Li, Y.H.; Zhu, H.W.; Wang, R.; Ren, Z.F.; He, Y.B. Bubble behavior and evolution characteristics in the RH riser tube-vacuum chamber. Int. J. Chem. React. Eng. 2022, 20, 1053–1064. [Google Scholar] [CrossRef]
  26. Zhu, G.S.; Deng, X.X.; Ji, C.X. Effect of vacuum degree on inclusion removal in ultra low carbon steel during RH refining process. Iron Steel 2022. [Google Scholar] [CrossRef]
  27. Wang, K.P.; Wang, Y.; Xu, J.F.; Xie, W.; Chen, T.J.; Jiang, M. Influence of Temperature Drop from 1773 K (1500 °C) to 1743 K (1470 °C) on Inclusion of High Carbon Chromium Bearing Steel Melts Treated by RH Vacuum Degassing. Metall. Mater. Trans. B 2022. [Google Scholar] [CrossRef]
  28. Wang, Q.; Tan, C.; Huang, A.; Yan, W.; Gu, H.Z.; He, Z.; Li, G.Q. Numerical simulation on refractory wear and inclusion formation in continuous casting tundish. Metall. Mater. Trans. B 2021, 52, 1344–1356. [Google Scholar] [CrossRef]
  29. Chang, S.; Zou, Z.S.; Li, B.K.; ISAC, M.; Guthrie, R.I.L. Modeling Inclusion Removal when Using Microbubble Swarm in a Full-Scale Tundish with an Impact Pad. Metall. Mater. Trans. B 2021, 53, 526–536. [Google Scholar] [CrossRef]
  30. Xin, Z.Y.; Cui, H.N.; Li, T.; Tang, G.Z.; Zhu, Y.L.; Yan, J.C.; Li, J.G. Numerical Simulation of Inclusion Removal by Bubble Injection in the Submerged Nozzle With Swirling Flow. Metall. Mater. Trans. B 2022, 53, 2570–2586. [Google Scholar] [CrossRef]
  31. Fan, J.P.; Li, Y.Q.; Chen, C.; Ouyang, X.; Wang, T.Y.; Lin, W.M. Effect of Uniform and Non-Uniform Increasing Casting Flow Rate on Dispersion and Outflow Percentage of Tracers in Four Strand Tundishes under Strand Blockage Conditions. Metals 2022, 12, 1016. [Google Scholar] [CrossRef]
  32. Reis, B.H.; Bielefeldt, W.V.; Vilela, A.C.F. Absorption of non-metallic inclusions by steelmaking slags—A review. J. Mater. Res. Technol. 2014, 3, 179–185. [Google Scholar] [CrossRef] [Green Version]
  33. Wang, R.; Li, Y.H.; Li, D.Z.; Kang, Y.; Bao, Y.P.; Yan, Z.J. Inclusions Absorbed by Slags in Interstitial-Free Steel Production. Steel Res. Int. 2020, 91, 1900440. [Google Scholar] [CrossRef]
  34. Wen, X.; Ren, Y.; Zhang, L.F. Effect of CaF2 Contents in Slag on Inclusion Absorption in a Bearing Steel. Steel Res. Int. 2022, 2200218. [Google Scholar] [CrossRef]
  35. Salgado, U.D.; Weiß, C.; Michelic, S.K.; Bernhard, C. Fluid force-induced detachment criteria for nonmetallic inclusions adhered to a refractory/molten steel interface. Metall. Mater. Trans. B 2018, 49, 1632–1643. [Google Scholar] [CrossRef]
  36. Wikström, J.; Nakajima, K.; Jonsson, L.; Jönsson, P. Application of a Model for Liquid Inclusion Separation at a Steel-Slag Interface for Laboratory and Industrial Situations. Steel Res. Int. 2008, 79, 826–834. [Google Scholar] [CrossRef]
  37. Chen, C.; Cheng, G.G.; Yang, H.K.; Hou, Z.B. Physical modeling of Fluid flow Characteristics in a Delta Shaped, Four-Strand Continuous Casting Tundish with Different Flow Control Devices. Adv. Mater. Res. 2011, 284–286, 1071–1079. [Google Scholar] [CrossRef]
  38. Chen, C.; Cheng, G.G.; Sun, H.B.; Hou, Z.B.; Wang, X.C.; Zhang, J.Q. Effects of salt tracer amount, concentration and kind on the fluid flow behavior in a hydrodynamic model of continuous casting tundish. Steel Res. Int. 2012, 83, 1141–1151. [Google Scholar] [CrossRef]
  39. Chen, C.; Rui, Q.X.; Cheng, G.G. Effect of salt tracer amount on the mixing time measurement in a hydrodynamic model of gas-stirred ladle system. Steel Res. Int. 2013, 84, 900–907. [Google Scholar] [CrossRef]
  40. Conejo, A.N.; Zhao, G.A.; Zhang, L.F. On the Limits of the Geometric Scale Ratio Using Water Modeling in Ladles. Metall. Mater. Trans. B 2021, 52, 2263–2274. [Google Scholar] [CrossRef]
  41. Yang, R.W.; Chen, C.; Lin, Y.C.; Zhao, Y.; Zhao, J.; Zhu, J.J.; Yang, S.S. Water model experiment on motion and melting of scarp in gas stirred reactors. Chin. J. Process Eng. 2022, 22, 954–962. [Google Scholar] [CrossRef]
  42. Pei, K.H.; Chen, C.; Zhao, Y.; Lin, Y.C.; Yang, R.W.; Zhu, J.J.; Wang, T.; Yang, K.; Lin, W.M. Water model experiment on the motion, melting, and mixing of scrap in bottom stirred reactors. Chin. J. Process Eng. 2022, 22, 1601–1612. [Google Scholar] [CrossRef]
  43. Lindborg, U.; Torssell, K. A Collision Model for the Growth and Separation of Deoxidation Products. Trans. Metall. Soc. AIME 1968, 242, 94–102. [Google Scholar]
  44. Linder, S. Hydrodynamics and Collisions of Small Particles in a Turbulent Metallic Melt with Special Reference to Deoxidation of Steel. Scand. J. Metall. 1974, 3, 137–150. [Google Scholar]
  45. Engh, T.A.; Lindskog, N. A Fluid Mechanical Model of Inclusion Removal. Scand. J. Metall. 1975, 4, 49–58. [Google Scholar]
  46. Söder, M.; Jönsson, P.; Alexis, J. Most relevant mechanisms of inclusion growth in an induction-stirred ladle. Scand. J. Metall. 2002, 31, 210–220. [Google Scholar] [CrossRef]
  47. Lei, H.; Nakajima, K.; He, J.C. Mathematical Model for Nucleation, Ostwald Ripening and Growth of Inclusion in Molten Steel. ISIJ Int. 2010, 50, 1735–1745. [Google Scholar] [CrossRef] [Green Version]
  48. Nakanishi, K.; Szekely, J. Deoxidation Kinetics in a Turbulent Flow Field. Trans. ISIJ 1975, 15, 522–530. [Google Scholar] [CrossRef]
  49. Johansen, S.T.; Boysan, F.; Engh, T.A. Numerical calculations of removal of inclusions and dissolution of refractory in bubble stirred ladles. In Proceedings of the 4th Japan-Nordic Countries Joint Symposiums on Science and Technology of Process Metallurgy, Tokyo, Japan; 1986; pp. 182–215. [Google Scholar]
  50. Miki, Y.; Shimada, Y.; Thomas, B.G.; Denissov, A. Model of Inclusion Removal during RH Degassing of Steel. Iron Steelmak. 1997, 28, 31–38. Available online: https://experts.illinois.edu/en/publications/model-of-inclusion-removal-during-rh-degassing-of-steel (accessed on 19 January 2023).
  51. Wakoh, M.; Fuchigama, K.; Endoh, K.; Imamura, N.; Kiyose, A.; Sawada, I. Analysis of Inclusion Behavior in a ladle refining process by using newly developed coagulation model. In Proceedings of the Scanmet I, Luleå, Sweden, 7–8 June 1999; pp. 267–274. [Google Scholar]
  52. Sheng, D.Y.; Söder, M.; Jönsson, P.; Jonsson, L. Modeling micro-inclusion growth and separation in gas-stirred ladles. Scand. J. Metall. 2002, 31, 134–147. [Google Scholar] [CrossRef]
  53. Söder, M.; Jönsson, P.; Jonsson, L. Inclusion growth and removal in gas-stirred ladles. Steel Res. Int. 2004, 75, 128–138. [Google Scholar] [CrossRef]
  54. Hallberg, M.; Jönsson, P.G.; Jonsson, L.T.I.; Eriksson, R. Process model of inclusion separation in a stirred steel ladle. Scand. J. Metall. 2005, 34, 41–56. [Google Scholar] [CrossRef]
  55. Kwon, Y.J.; Zhang, J.; Lee, H.G. A CFD-based Nucleation-growth-removal Model for Inclusion Behavior in a Gas-agitated Ladle during Molten Steel Deoxidation. ISIJ Int. 2008, 48, 891–900. [Google Scholar] [CrossRef] [Green Version]
  56. Geng, D.Q.; Lei, H.; He, J.C. Numerical Simulation for Collision and Growth of Inclusions in Ladles Stirred with Different Porous Plug Configurations. ISIJ Int. 2010, 50, 1597–1605. [Google Scholar] [CrossRef] [Green Version]
  57. Claudotte, L.; Rimbert, N.; Gardin, P.; Simonnet, M.; Lehmann, L.; Oesterlé, B. Behaviour of Oxide Inclusions in Liquid Steel Multi-QMOM Simulation. Steel Res. Int. 2010, 81, 630–636. [Google Scholar] [CrossRef]
  58. Daoud, I.L.A.; Rimbert, N.; Jardy, A.; Oesterlé, B.; Hans, S.; Bellot, J.P. 3D modeling of the aggregation of oxide inclusions in a liquid steel ladle: Two numerical approaches. Adv. Eng. Mater. 2011, 13, 543–549. [Google Scholar] [CrossRef]
  59. Gardin, P.; Gauthier, S.; Simonnet, M.; Lehmann, J. Multiphase Model for Predicting the Elimination of Inclusions inside a Liquid-Steel Ladle. Adv. Eng. Mater. 2011, 13, 538–542. [Google Scholar] [CrossRef]
  60. Lou, W.T.; Zhu, M.Y. Numerical simulations of inclusion behavior in gas-stirred ladles. Metall. Mater. Trans. B 2013, 44, 762–782. [Google Scholar] [CrossRef]
  61. Chen, C.; Jonsson, L.T.I.; Tilliander, A.; Cheng, G.G.; Jönsson, P.G. A mathematical modeling study of tracer mixing in a continuous casting tundish. Metall. Mater. Trans. B 2015, 46, 169–190. [Google Scholar] [CrossRef]
  62. Chen, C.; Jonsson, L.T.I.; Tilliander, A.; Cheng, G.G.; Jönsson, P.G. A mathematical modeling study of the influence of small amounts of KCl solution tracers on mixing in water and residence time distribution of tracers in a continuous flow reactor-metallurgical tundish. Chem. Eng. Sci. 2015, 137, 914–937. [Google Scholar] [CrossRef]
  63. Guo, X.P.; Godinez, J.; Walla, N.J.; Silaen, A.K.; Oltmann, H.; Thapliyal, V.; Bhansali, A.; Pretorius, E.; Zhou, C.Q. Computational Investigation of Inclusion Removal in the Steel-Refining Ladle Process. Processes 2021, 9, 1048. [Google Scholar] [CrossRef]
  64. Wang, F.; Sun, B.Y.; Liu, Z.Q.; Li, B.K.; Huang, S.; Zhang, B.J. Numerical Simulation on Motion Behavior of Inclusions in the Lab-Scale Electroslag Remelting Process with a Vibrating Electrode. Metals 2021, 11, 1784. [Google Scholar] [CrossRef]
  65. Siddiqui, M.I.H.; Albaqami, A.; Arifudin, L.; Alluhydan, K.; Alnaser, I.A. Simulation of Inclusion Particle Motion Behavior under Interfacial Tension in Continuous Casting Mold. Materials 2022, 15, 7458. [Google Scholar] [CrossRef]
  66. Chen, C.; Ni, P.Y.; Jonsson, L.T.I.; Tilliander, A.; Cheng, G.G.; Jönsson, P.G. A Model Study of Inclusions Deposition, Macroscopic Transport, and Dynamic Removal at Steel–Slag Interface for Different Tundish Designs. Metall. Mater. Trans. B 2016, 47, 1916–1932. [Google Scholar] [CrossRef]
  67. Tacke, K.H.; Ludwig, J.C. Steel flow and inclusion separation in continuous casting tundishes. Steel Res. 1987, 58, 262–270. [Google Scholar] [CrossRef]
  68. Rückert, A.; Warzecha, M.; Koitzsch, R.; Pawlik, M.; Pfeifer, H. Particle Distribution and Separation in Continuous Casting Tundish. Steel Res. Int. 2009, 80, 568–574. [Google Scholar] [CrossRef]
  69. Gutiérrez, E.; Garcia-Hernandez, S.; Barreto, J.J. Mathematical Analysis of the Touching Inclusions Parameters at the Tundish Free Surface to Predict More Realistic Inclusion Removal Rates. Steel Res. Int. 2019, 90, 1900328. [Google Scholar] [CrossRef]
  70. Chen, C. Some Aspects on Macroscopic Mixing in a Tundish. Ph.D. Thesis, KTH Royal Institute of Technology, Stockholm, Sweden, 2015. [Google Scholar]
  71. Nakajima, K.; Okamura, K. Inclusion Transfer Behavior across Molten Steel-Slag Interface. In Proceedings of the 4th International Conference on Molten Slags and Fluxes, Sendai, Japan, 8–11 June 1992; pp. 505–510. [Google Scholar]
  72. Bouris, D.; Bergeles, G. Investigation of Inclusion Re-Entrainment from the Steel-Slag Interface. Metall. Mater. Trans. B 1998, 29, 641–649. [Google Scholar] [CrossRef]
  73. Strandh, J.; Nakajima, K.; Eriksson, R.; Jönsson, P. Solid inclusion transfer at a steel-slag interface with focus on tundish conditions. ISIJ Int. 2005, 45, 1597–1606. [Google Scholar] [CrossRef] [Green Version]
  74. Strandh, J.; Nakajima, K.; Eriksson, R.; Jönsson, P. A mathematical model to study liquid inclusion behavior at the steel-slag interface. ISIJ Int. 2005, 45, 1838–1847. [Google Scholar] [CrossRef] [Green Version]
  75. Xuan, C.; Persson, E.S.; Sevastoplev, R.; Nzotta, M. Motion and detachment behaviors of liquid inclusion at molten steel-slag interfaces. Metall. Mater. Trans. B 2019, 50, 1957–1973. [Google Scholar] [CrossRef]
  76. Shannon, G.N.; Sridhar, S. Modeling Al2O3 inclusion separation across steel-slag interfaces. Scand. J. Metall. 2005, 34, 353–362. [Google Scholar] [CrossRef]
  77. Shannon, G.N.; Sridhar, S. Film-Drainage, Separation and Dissolution of A12O3 Inclusions across Interfaces between Molten Steel and Ladle-, Tundish- and Mold-Slags. High Temp. Mater. Proc. 2005, 24, 111–124. [Google Scholar] [CrossRef]
  78. Zhou, Y.L.; Deng, Z.Y.; Zhu, M.Y. Consideration on removal of non-spherical solid inclusions at slag-steel interface. Chin. J. Process Eng. 2022, 22, 222–231. [Google Scholar] [CrossRef]
  79. Valdez, M.; Shannon, G.S.; Sridhar, S. The Ability of Slags to Absorb Solid Oxide Inclusions. ISIJ Int. 2006, 46, 450–457. [Google Scholar] [CrossRef] [Green Version]
  80. Shannon, G.; White, L.; Sridhar, S. Modeling inclusion approach to the steel/slag interface. Mater. Sci. Eng. A 2008, 495, 310–315. [Google Scholar] [CrossRef]
  81. Liu, C.; Yang, S.F.; Li, J.S.; Zhu, L.B.; Li, X.G. Motion Behavior of Nonmetallic Inclusions at the Interface of Steel and Slag. Part I: Model Development, Validation, and Preliminary Analysis. Metall. Mater. Trans. B 2016, 47, 1882–1892. [Google Scholar] [CrossRef] [Green Version]
  82. Yang, S.F.; Li, J.S.; Liu, C.; Sun, L.Y.; Yang, H.B. Motion Behavior of Nonmetallic Inclusions at the Interface of Steel and Slag. Part II: Model Application and Discussion. Metall. Mater. Trans. B 2014, 45, 2453–2463. [Google Scholar] [CrossRef]
  83. Liu, W.; Yang, S.F.; Li, J.S. Calculation of Static Suspension Depth and Meniscus Shape of a Solid Spherical Inclusion at the Steel–Slag Interface. Metall. Mater. Trans. B 2020, 51, 422–425. [Google Scholar] [CrossRef]
  84. Liu, W.; Liu, J.; Zhao, H.X.; Yang, S.F.; Li, J.S. CFD Modeling of Solid Inclusion Motion and Separation from Liquid Steel to Molten Slag. Metall. Mater. Trans. B 2021, 52, 2430–2440. [Google Scholar] [CrossRef]
  85. Zhang, X.M.; Pirker, S.; Saeedipour, M. Numerical investigation of particle motion at the steel–slag interface in continuous casting using VOF method and dynamic overset grids. Exp. Comput. Multiph. Flow 2023, 5, 178–191. [Google Scholar] [CrossRef]
  86. Du, Y.N.; Guo, L.; Yang, Y.; Zhang, S.; Yu, H.Z.; Guo, Z.C. Numerical simulation of inclusions floating behavior under supergravity field in liquid steel. Iron Steel 2022, 57, 54–62. [Google Scholar] [CrossRef]
  87. Liu, Y.Q.; Wang, L.J.; Chou, K.C. Dissolution Behavior of Al2O3 in Refining Slags Containing Ce2O3. ISIJ Int. 2014, 54, 728–733. [Google Scholar] [CrossRef]
  88. Wang, L.J.; Wang, Q.; Li, J.M.; Chou, K.C. Dissolution mechanism of Al2O3 in refining slags contained Ce2O3. J. Min. Metall. Sect. B-Metall. 2016, 52, 35–40. [Google Scholar] [CrossRef] [Green Version]
  89. Shi, G.Y.; Zhang, T.A.; Dou, Z.H.; Niu, L.P. Dissolution Behavior of Al2O3 Inclusions in CaO-Al2O3 Based Slag Representing Aluminothermic Reduction Slag. Crystals 2020, 10, 1061. [Google Scholar] [CrossRef]
  90. Zhao, S.; Li, Z.S.; Xu, R.Z.; Khasraw, D.; Song, G.Y.; Xu, D. Dissolution Behavior of Different Inclusions in High Al Steel Reacted with Refining Slags. Metals 2021, 11, 1801. [Google Scholar] [CrossRef]
  91. Yeo, S.; Um, H.; Chung, Y. The Effect of Alumina Activity on Dissolution Behavior of Alumina Particles in CaO-Al2O3-SiO2 Slags. Metall. Mater. Trans. B 2021, 52, 3938–3945. [Google Scholar] [CrossRef]
  92. Ren, C.Y.; Zhang, L.F.; Zhang, J.; Wu, S.J.; Zhu, P.; Ren, Y. In Situ Observation of the Dissolution of Al2O3 Particles in CaO-Al2O3-SiO2 Slags. Metall. Mater. Trans. B 2021, 52, 3288–3301. [Google Scholar] [CrossRef]
  93. Song, S.Y.; Li, J.; Yan, W. Dissolution mechanism of alumina inclusions into CaO–SiO2–Al2O3–FetO–MgO slag via a converter tapping process. Int. J. Appl. Ceram. Technol. 2022, 19, 2734–2748. [Google Scholar] [CrossRef]
  94. Zhang, L.F.; Ren, Y. Concept of inclusion capacity of slag and its application. Iron Steel 2022. [Google Scholar] [CrossRef]
  95. Xuan, C.J.; Mu, W.Z. A phase-field model for the study of isothermal dissolution behavior of alumina particles into molten silicates. J. Am. Ceram. Soc. 2019, 102, 6480–6497. [Google Scholar] [CrossRef]
  96. Xuan, C.J.; Mu, W.Z. A mechanism theory of dissolution profile of oxide particles in oxide melt. J. Am. Ceram. Soc. 2020, 104, 57–75. [Google Scholar] [CrossRef]
  97. Mu, W.Z.; Xuan, C.J. Phase-field study of dissolution behaviors of different oxide particles into oxide melts. Ceram. Int. 2020, 46, 14949–14956. [Google Scholar] [CrossRef]
  98. Xuan, C.J.; Mu, W.Z. Dissolution kinetics of arbitrarily-shaped alumina in oxide melt: An integration of phase-field modelling and real-time observation study. J. Alloys Compd. 2020, 834, 155168. [Google Scholar] [CrossRef]
  99. Pan, X.Y.; Shen, F.M.; Gao, Q.J.; Jiang, X.; Zheng, H.Y. Effect of Erosion Behavior of FeO-CaO-SiO2-MgO-Al2O3 Blast Furnace Primary Slag on Al2O3 Substrate. Crystals 2021, 11, 957. [Google Scholar] [CrossRef]
  100. Zheng, H.Y.; Wang, S.; Guo, Y.C.; Jiang, X.; Gao, Q.J.; Shen, F.M. Thermodynamic model for calculating sulfur distribution ratio between CaO-SiO2-MgO-Al2O3 slag and hot metal. Iron Steel 2021, 56, 16–23. [Google Scholar] [CrossRef]
  101. Hu, X.G.; Shen, F.M.; Zheng, H.Y.; Guo, Y.C.; Jiang, X.; Gao, Q.J. Effect of temperature on Al2O3 activity in CaO-SiO2-Al2O3-MgO blast furnace slag system. Iron Steel 2022, 57, 25–33. [Google Scholar] [CrossRef]
  102. Liu, K.; Han, Y.H.; Yang, F.; Zhu, L.G. Effect of MgO on microstructure of CaO-Al2O3 based quaternary system mold fluxes. Iron Steel 2020, 55, 52–58. [Google Scholar] [CrossRef]
  103. Zhao, B.B.; Zhang, J.; Yan, B.J. Interfacial Phenomena and Reaction Kinetics between High Al Molten Steel and CaO-SiO2-Type Flux. Metals 2022, 12, 391. [Google Scholar] [CrossRef]
  104. Wang, X.F.; Wang, Q.Q.; Zhang, X.B.; Wang, Q.; He, S.P. Effect of AlN on properties of CaO-SiO2 based mold fluxes for high-Al steel. Iron Steel 2022, 57, 64–71. [Google Scholar] [CrossRef]
  105. Silva, A.M.B.; Oliveira, M.A.; Peixoto, J.J.M.; Silva, C.A.D. Slag-Steel Emulsification on a Modified RH Degasser. Metall. Mater. Trans. B 2021, 52, 2111–2126. [Google Scholar] [CrossRef]
  106. Li, S.J.; Cheng, G.G.; Yang, L.; Chen, L.; Yan, Q.Z.; Li, C.W. A Thermodynamic Model to Design the Equilibrium Slag Compositions during Electroslag Remelting Process: Description and Verification. ISIJ Int. 2017, 4, 713–722. [Google Scholar] [CrossRef] [Green Version]
  107. Bublik, S.; Tangstad, M.; Einarsrud, K.R.E. Interfacial Behaviour in Ferroalloys: The Influence of FeMn Slag Composition. Metall. Mater. Trans. B 2022, 53, 3276–3291. [Google Scholar] [CrossRef]
  108. Bai, C.G.; Yan, Z.M.; Pang, Z.D.; Jiang, Y.Y.; Ling, J.W.; Lü, X.W. Advances of measurement and calculation model of slag viscosity. Iron Steel 2020, 55, 27–37. [Google Scholar] [CrossRef]
  109. Lü, X.W.; Yan, Z.M.; Pang, Z.D.; Bai, C.G.; Liang, D.; Xie, H. Effect of Al2O3 on physicochemical properties and structure of blast furnace slag: Review. Iron Steel 2020, 55, 1–10. [Google Scholar] [CrossRef]
  110. Fang, J.L.; Pang, Z.G.; Xing, X.D.; Xu, R.S. Thermodynamic Properties, Viscosity, and Structure of CaO–SiO2–MgO–Al2O3–TiO2–Based Slag. Materials 2021, 14, 124. [Google Scholar] [CrossRef]
  111. Haghdani, S.; Tangstad, M.; Einarsrud, K.R.E. A Raman-Structure Model for the Viscosity of SiO2-CaO-Al2O3 System. Metall. Mater. Trans. B 2022, 53, 1733–1746. [Google Scholar] [CrossRef]
  112. Zhang, R.; Meng, Y.F.; Wang, Z.; Jiao, S.Y.; Jia, J.X.; Min, Y.; Liu, C.J. Surface Tension Prediction of Multicomponent Slags of the CaO–SiO2–Al2O3–FeO–Fe2O3 System Based on Microstructure Analysis. Metall. Mater. Trans. B 2022, 53, 571–583. [Google Scholar] [CrossRef]
  113. Wang, L.; Cheng, S.S.; Liu, P.B.; Chen, Y.B. Effect of K2O and Na2O on activity of components and MgO content of high Al2O3 slag. Iron Steel 2022, 57, 48–56. [Google Scholar] [CrossRef]
  114. Zhao, M.G.; Li, G.; Li, Z.; Wang, Q.Q.; He, S.P. Study of Thermodynamic for Low-Reactive CaO-BaO-Al2O3-SiO2-CaF2-Li2O Mold Flux Based on the Model of Ion and Molecular Coexistence Theory. Metals 2022, 12, 1099. [Google Scholar] [CrossRef]
  115. He, S.P.; Chen, G.J.; Guo, Y.T.; Shen, B.Y.; Wang, Q. Morphology control for Al2O3 inclusions without Ca treatment in high aluminum steel. Metall. Mater. Trans. B 2015, 46, 585–594. [Google Scholar] [CrossRef] [Green Version]
  116. Hasegawa, M.; Sakuma, K. Use of Rare Earth Elements in Steel Making. Tetsu-to-Hagané 1956, 42, 503–517. [Google Scholar] [CrossRef] [Green Version]
  117. Malm, S. On the Precipitation of Slag Inclusions during Solidification of High-carbon Steels Deoxidized with Aluminium and Misch Metal. Scand. J. Metall. 1976, 5, 248–257. [Google Scholar]
  118. Emi, T.; Haida, O.; Sakuraya, T.; Sanbongi, K. Mechanism of sulfides precipitation during solidification of calcium- and rare earth treated ingots. In Proceedings of the Steelmaking Proceedings, Detroit, MI, USA; 1979; Volume 62, pp. 574–584. [Google Scholar]
  119. Zhao, S.; Han, Q.; Shen, F.; Shao, G.; Xie, Z. Study on the Nozzle Blockage Mechanism of Steel Containing Rare Earth Metals during Casting. Scand. J. Metall. 1985, 14, 175–182. [Google Scholar]
  120. Mizuno, M.; Berjoan, R.; Coutures, J.P.; Foex, M. Phase Diagram of the System Al2O3-La2O3 at Elevated Temperature. Yogyo-Kyokai-Shi 1974, 82, 631–636. [Google Scholar] [CrossRef]
  121. Mizuno, M.; Berjoan, R.; Coutures, J.P.; Foex, M. Phase Diagram of the System Al2O3-CeO2 at Liquidus Temperature. Yogyo-Kyokai-Shi 1975, 83, 90–96. [Google Scholar] [CrossRef]
  122. Ueda, S.; Morita, L.; Sano, N. Activity of AlO1.5 for the CaO–AlO1.5–CeO1.5 System at 1773 K. ISIJ Int. 1998, 38, 1292–1296. [Google Scholar] [CrossRef]
  123. Wu, P.; Pelton, A.D. Coupled thermodynamic-phase diagram assessment of the rare earth oxide-aluminium oxide binary systems. J. Alloys Compd. 1992, 179, 259–287. [Google Scholar] [CrossRef]
  124. Navrotsky, A.; Lee, W.; Mielewczyk-Gryn, A.; Ushakov, S.V.; Anderko, A.; Wu, H.H.; Riman, R.E. Thermodynamics of solid phases containing rare earth oxides. J. Chem. Thermodyn. 2015, 88, 126–141. [Google Scholar] [CrossRef] [Green Version]
  125. Sun, L.F.; Liu, Z.Y.; Jiang, M.F. Phase Equilibria of CaO-SiO2-La2O3-Nb2O5 System in Reducing Atmosphere. Metals 2022, 12, 768. [Google Scholar] [CrossRef]
  126. Jeong, S.J.; Kim, T.S.; Park, J.H. Relationship Between Sulfide Capacity and Structure of MnO-SiO2-Al2O3-Ce2O3 System. Metall. Mater. Trans. B 2017, 48, 545–553. [Google Scholar] [CrossRef]
  127. Qi, J.; Liu, C.J.; Zhang, C.; Jiang, M.F. Effect of Ce2O3 on Structure, Viscosity, and Crystalline Phase of CaO-Al2O3-Li2O-Ce2O3 Slags. Metall. Mater. Trans. B 2017, 48, 11–16. [Google Scholar] [CrossRef]
  128. Ma, Z.; Zhao, Z.W.; Guo, W.T.; Wang, Z. Influence of Nb2O5 and basicity on the viscosity and structure of CaO-SiO2-Nb2O5-CeO2-CaF2 slag system. Metall. Res. Technol. 2020, 117, 307. [Google Scholar] [CrossRef]
  129. Wu, C.C.; Cheng, G.G.; Long, H. Calculating model of action concentration of slag systems containing Ce2O3. Chin. J. Nonferrous Met. 2013, 23, 2999–3005. [Google Scholar] [CrossRef]
  130. Wu, C.C.; Cheng, G.G.; Long, H.; Yang, X.H. A Thermodynamic Model for Evaluation of Mass Action Concentrations of Ce2O3-contained Slag Systems Based on the Ion and Molecule Coexistence Theory. High Temp. Mater. Proc. 2013, 32, 207–214. [Google Scholar] [CrossRef]
  131. Wu, C.C.; Cheng, G.G.; Tian, J. A Thermodynamic Model for Evaluation of Mass Action Concentrations of La2O3-Al2O3-CaF2-CaO-MgO Slags for Electroslag Remelting Based on the Ion and Molecule Coexistence Theory. High Temp. Mater. Proc. 2013, 32, 541–550. [Google Scholar] [CrossRef]
  132. Wu, C.C.; Cheng, G.G.; Long, H. Effect of Ce2O3 and CaO/Al2O3 on the Phase, Melting Temperature and Viscosity of CaO-Al2O3-10 Mass% SiO2 Based Slags. High Temp. Mater. Proc. 2014, 33, 77–84. [Google Scholar] [CrossRef]
  133. Long, H.; Cheng, G.G.; Wu, B.; Wu, Y. Melting and Fluidity Properties of Refining Slag Containing Ce2O3 for Steelmaking. J. Chin. Rare Earth Soc. 2010, 28, 721–727. Available online: http://www.re-journal.com/WKD3/WebPublication/paperDigest.aspx?paperID=28ff3254-3d8b-4c9d-83fc-f0ff6e2f874b (accessed on 20 November 2022).
  134. Wu, C.C.; Cheng, G.G. Calculating models on surface tension of RE2O3-MgO-SiO2 (RE = La, Nd, Sm, Gd and Y) melts. Trans. Nonferrous Met. Soc. China 2014, 24, 3696–3701. [Google Scholar] [CrossRef]
  135. Yang, X.H.; Long, H.; Cheng, G.G.; Wu, C.C.; Wu, B. Effect of refining slag containing Ce2O3 on steel cleanliness. J. Rare Earths 2011, 29, 1079–1083. [Google Scholar] [CrossRef]
  136. Wu, C.C. Research on the Rare Earth Addition Technology and Theory in Refining Process for Special Steel. Ph.D. Thesis, University of Science and Technology Beijing, Beijing, China, 2015. [Google Scholar]
  137. Babenko, A.A.; Smirnov, L.A.; Upolovnikova, A.G.; Shartdinov, R.R. Study of possibility of cerium reduction from slags of CaO-SiO2-Ce2O3-15%Al2O3-8%MgO system. IOP Conf. Ser. Mater. Sci. Eng. 2020, 966, 012010. [Google Scholar] [CrossRef]
  138. Li, D.Z.; Wang, P.; Chen, X.Q.; Fu, P.X.; Luan, Y.K.; Hu, X.Q.; Liu, H.W.; Sun, M.Y.; Chen, Y.; Cao, Y.F.; et al. Low-oxygen rare earth steels. Nat. Mater. 2022, 21, 1137–1143. [Google Scholar] [CrossRef]
  139. Kang, J.; Yu, Y.C.; Zhang, J.L.; Chen, C.; Wang, S.B. Effect of Rare Earth on Inclusion Evolution in Industrial Production of HRB500E Steel. Metall. Res. Technol. 2021, 118, 220. [Google Scholar] [CrossRef]
  140. Wang, H.; Bao, Y.P.; Zhi, J.G.; Duan, C.Y.; Gao, S.; Wang, M. Effect of Rare Earth Ce on the Morphology and Distribution of Al2O3 Inclusions in High Strength IF Steel Containing Phosphorus during Continuous Casting and Rolling Process. ISIJ Int. 2021, 61, 657–666. [Google Scholar] [CrossRef]
  141. Zhang, J.; Ren, Y.; Ren, Q.; Zhang, L.F. Transient Evolution of Nonmetallic Inclusions in a Si–Mn-Killed Stainless Steel with Cerium Addition. Steel Res. Int. 2022, 93, 2100773. [Google Scholar] [CrossRef]
  142. Yang, H.X.; Ren, Y.; Ji, S.; Zhang, L.F. Modification of Sulfides in a High Sulfur Steel by Cerium Addition. Metall. Mater. Trans. B 2022. [Google Scholar] [CrossRef]
  143. Bao, D.H.; Cheng, G.G.; Huang, Y.; Qiao, T.; Dai, W.X. Refinement of Solidification Structure of H13 Steel by Rare Earth Sulfide. Steel Res. Int. 2022, 93, 2100304. [Google Scholar] [CrossRef]
  144. Liu, D.Q.; Liu, M.; Jiang, X.; Yu, Y.H.; Zhang, Z.M.; Feng, X.M.; Lai, C.B. Effect of Yttrium-Based Rare Earth on Inclusions and Cryogenic Temperature Impact Properties of Offshore Engineering Steel. Crystals 2022, 12, 305. [Google Scholar] [CrossRef]
  145. Jiang, Y.Y.; Li, Y.C.; Dai, Q.W.; Chai, S.S.; Li, K.J. Evolution of inclusions with cerium addition and effects of C-containing rare earth inclusions on the toughness of ultra-high-strength structural steel. Ironmak. Steelmak. 2022, 2099696. [Google Scholar] [CrossRef]
  146. Huang, F.; Li, J.; Geng, R.M.; Wang, C.C. Effect of rare earth on inclusion evolution and corrosion resistance of HRB400E steel. Mater. Corros. 2022, 13104. [Google Scholar] [CrossRef]
  147. Zhang, J. Computational Thermodynamics of Metallurgical Melt and Solutions, 1st ed.; Metallurgical Industry Publisher: Beijing, China, 2007; p. 7. [Google Scholar]
  148. Cheng, G.G.; Liao, N.B. Calculation Model for Surface Tension of Slag Melt. J. Iron Steel Res. Int. 1999, 6, 17–20. Available online: https://kns.cnki.net/KCMS/detail/detail.aspx?dbcode=CJFD&filename=YING199902004 (accessed on 20 November 2022).
  149. Zhou, Y.; Zhu, R.; Wang, H.Y.; Zhang, H.J. Effect of various components on the distribution of phosphorus in CaO–FeO–MgO–SiO2–MnO–TiO2–V2O5–P2O5 slag based on IMCT. Ironmak. Steelmak. 2021, 48, 570–578. [Google Scholar] [CrossRef]
  150. Li, B.; Li, L.; Guo, H.J.; Guo, J.; Duan, S.C.; Sun, W.X. A phosphorus distribution prediction model for CaO–SiO2–MgO–FeO–Fe2O3–Al2O3–P2O5 slags based on the IMCT. Ironmak. Steelmak. 2020, 47, 771–780. [Google Scholar] [CrossRef]
  151. Zhang, G.L.; Cheng, G.G.; Wang, Y.P.; Zhang, X.; Tan, B.; Yuan, C.W.; Miao, H.S. A thermodynamic model of the titanium distribution behaviour between slag and molten steel during LF refining process: Description and application. Ironmak. Steelmak. 2022, 49, 679–692. [Google Scholar] [CrossRef]
  152. Sun, H.; Yang, J.; Yang, W.K.; Zhang, R.H. Evaluation of Phosphorus Enrichment Capacity of CaO–SiO2–FeO–MgO–MnO–P2O5–Al2O3 Dephosphorization Slag Based on Ion-Molecule Coexistence Theory. Steel Res. Int. 2022, 2200662. [Google Scholar] [CrossRef]
  153. Butler, J.A.V. The Thermodynamics of the Surfaces of Solutions. Roy. Soc. 1932, 135, 348–375. [Google Scholar] [CrossRef]
  154. Wu, C.C.; Cheng, G.G.; Ma, Q.Q. Calculation model of Surface Tension of CaF2-CaO-Al2O3-MgO-SiO2 Slag System. Shanghai Metals 2014, 36, 33–36, 41. [Google Scholar] [CrossRef]
  155. Wu, C.C.; Cheng, G.G. Calculating Models on the Surface Tension of CaO-Al2O3-TiO2 Molten Slag. Shanghai Metals 2014, 36, 35–39. [Google Scholar] [CrossRef]
  156. Wu, C.C.; Cheng, G.G. Computational Model on Surface Tension of CaO-MnO-SiO2 Slag System. Hot Work. Technol. 2014, 43, 58–62. [Google Scholar] [CrossRef]
  157. Girifalco, L.A.; Good, R.J. A theory for the estimation of surface and interfacial energies. Ⅰ. Derivation and application to interfacial tension. J. Phys. Chem. 1957, 61, 904–909. [Google Scholar] [CrossRef]
  158. Nakajima, K. Estimation of Interfacial Tensions between Phases in the Molten Iron-Slag-Inclusion (Alumina) System. Tetsu-to-Hagané 1994, 80, 383–388. [Google Scholar] [CrossRef] [Green Version]
  159. Yang, H.B. Study on Moving Behavior of Inclusion during Process of Passing Steel-Slag Interface. Ph.D. Thesis, University of Science and Technology Beijing, Beijing, China, 2015. [Google Scholar]
  160. Conejo, A.N.; Feng, W.H. Ladle Eye Formation Due to Bottom Gas Injection: A Reassessment of Experimental Data. Metall. Mater. Trans. B 2022, 53, 999–1017. [Google Scholar] [CrossRef]
Figure 1. Flow chart of integrated calculation model of inclusion motion.
Figure 1. Flow chart of integrated calculation model of inclusion motion.
Crystals 13 00202 g001
Figure 2. The relationship between σIM, σIS, σMS, and θIMS. (a) cosθIMS > 0, (b) cosθIMS < 0.
Figure 2. The relationship between σIM, σIS, σMS, and θIMS. (a) cosθIMS > 0, (b) cosθIMS < 0.
Crystals 13 00202 g002
Figure 3. Three kinds of interfacial motion behaviors of inclusions (a) Pass, (b) Remain, (c) Oscillate.
Figure 3. Three kinds of interfacial motion behaviors of inclusions (a) Pass, (b) Remain, (c) Oscillate.
Crystals 13 00202 g003
Figure 4. Schematic diagram of experimental equipment.
Figure 4. Schematic diagram of experimental equipment.
Crystals 13 00202 g004
Figure 5. Comparison of experimental result and calculated result.
Figure 5. Comparison of experimental result and calculated result.
Crystals 13 00202 g005
Figure 6. Force analysis of solid Al2O3 inclusions with different sizes, (a) 5 μm, (b) 10 μm, (c) 20 μm, (d) 50 μm.
Figure 6. Force analysis of solid Al2O3 inclusions with different sizes, (a) 5 μm, (b) 10 μm, (c) 20 μm, (d) 50 μm.
Crystals 13 00202 g006aCrystals 13 00202 g006b
Figure 7. The displacements of solid Al2O3 inclusions with different radiuses of inclusions in slag A5.
Figure 7. The displacements of solid Al2O3 inclusions with different radiuses of inclusions in slag A5.
Crystals 13 00202 g007
Figure 8. The displacements of liquid 50%wtAl2O3–50%wtCaO inclusions with different radiuses of inclusions in slag A5.
Figure 8. The displacements of liquid 50%wtAl2O3–50%wtCaO inclusions with different radiuses of inclusions in slag A5.
Crystals 13 00202 g008
Figure 9. The displacement of solid Al2O3 inclusions with different radiuses of inclusions in slag C3.
Figure 9. The displacement of solid Al2O3 inclusions with different radiuses of inclusions in slag C3.
Crystals 13 00202 g009
Table 1. The viscosities of CaO-Ce2O3-SiO2-Al2O3 of different compositions at 1773 K (mass fraction).
Table 1. The viscosities of CaO-Ce2O3-SiO2-Al2O3 of different compositions at 1773 K (mass fraction).
SerialCaO
%wt
Al2O3
%wt
SiO2
%wt
Ce2O3
%wt
Basicity
C/A
Viscosity (Pa·s)
A154.6430.361051.80.416
A251.0034.001051.50.363
A346.3638.641051.20.398
A451.4328.5710101.80.497
A548.2126.7910151.80.289
Table 2. The mole fractions of studied CaO-Ce2O3-SiO2-Al2O3 slag cases.
Table 2. The mole fractions of studied CaO-Ce2O3-SiO2-Al2O3 slag cases.
SerialXCaOXAl2O3XSiO2XCe2O3
A10.67050.20450.11450.0105
A20.63870.23380.11690.0107
A30.59620.27280.12000.0110
A40.65800.20070.11940.0218
A50.64440.19660.12480.0342
Table 3. Mass fraction composition and viscosity of CaO-Ce2O3-SiO2-Al2O3-MgO slag system at 1823 K.
Table 3. Mass fraction composition and viscosity of CaO-Ce2O3-SiO2-Al2O3-MgO slag system at 1823 K.
SerialCaO
%wt
Al2O3
%wt
SiO2
%wt
MgO
%wt
Ce2O3
%wt
Basicity
C/A
Viscosity (Pa·s)
C04837.54.51001.270.495
C14537.54.51031.20.307
C246.236.34.51031.270.374
C34535.54.51051.270.225
C443.434.14.51081.270.265
C54532.54.51081.3850.264
Table 4. The mole fraction of studied CaO-Ce2O3-SiO2-Al2O3-MgO slag cases.
Table 4. The mole fraction of studied CaO-Ce2O3-SiO2-Al2O3-MgO slag cases.
SerialXCaOXAl2O3XSiO2XMgOXCe2O3
C00.55310.23710.04840.16130
C10.53380.24420.04980.16610.0061
C20.54450.23490.04950.1650.006
C30.53860.23330.05030.16760.0102
C40.53130.22920.05140.17140.0167
C50.54610.21650.0510.16990.0166
Table 5. The relationships between surface tension and molar volume of pure components and temperature.
Table 5. The relationships between surface tension and molar volume of pure components and temperature.
ComponentRelationship between Surface Tension (×10−3 N/m) and Temperature (K)Relationship between Molar Volume (m3/mol) and Temperature (K)
CaO791 − 0.0935 T20.7 [1 + 1·10−4·(T-1773)]·10−6
Al2O31024 − 0.177 T28.3 [1 + 1·10−4·(T-1773)]·10−6
SiO2243.2 + 0.031 T27.516 [1 + 1·10−4·(T-1773)]·10−6
MgO1770 − 0.636 T16.1 [1 + 1·10−4·(T-1773)]·10−6
Table 6. Surface tension and molar volume data for several lanthanide oxides.
Table 6. Surface tension and molar volume data for several lanthanide oxides.
ElementOxideSurface Tension
(×10−3 N/m, 1873 K)
Surface Tension (×10−3 N/m, 2573 K)Relationship between Molar Volume (m3/mol) and Temperature (K)
LaLa2O3686 (Solid)560 (Liquid)50.1 [1 + 1·10−4·(T-1773)]·10−6
CeCe2O3To be estimatedUnknownTo be estimated
NdNd2O3677 (Solid)Unknown46.62 [1 + 1·10−4·(T-1773)]·10−6
SmSm2O3Unknown815 (Liquid)47.93 [1 + 1·10−4·(T-1773)]·10−6
GdGd2O3664 (solid)Unknown49.09 [1 + 1·10−4·(T-1773)]·10−6
Table 7. Surface tension of CaO-Ce2O3-SiO2-Al2O3 slag system at 1773 K.
Table 7. Surface tension of CaO-Ce2O3-SiO2-Al2O3 slag system at 1773 K.
SerialCaO
%wt
Al2O3
%wt
SiO2
%wt
Ce2O3
%wt
Basicity
C/A
Surface Tension
(N/m)
A154.6430.361051.80.5693
A251.0034.001051.50.5690
A346.3638.641051.20.5693
A451.4328.5710101.80.5682
A548.2126.7910151.80.5670
Table 8. Surface tension of CaO-Ce2O3-SiO2-Al2O3-MgO slag system at 1823 K.
Table 8. Surface tension of CaO-Ce2O3-SiO2-Al2O3-MgO slag system at 1823 K.
SerialCaO
%wt
Al2O3
%wt
SiO2
%wt
MgO
%wt
Ce2O3
%wt
Surface Tension (N/m)
C04837.54.51000.5992
C14537.54.51030.6725
C246.236.34.51030.6735
C34535.54.51050.6736
C443.434.14.51080.6740
C54532.54.51080.6754
Table 9. Measured data of interaction coefficient between slag system and Al2O3 inclusion.
Table 9. Measured data of interaction coefficient between slag system and Al2O3 inclusion.
XCaOXAl2O3XSiO2XMgOφIS
0.4530.1250.42300.630
0.4860.0590.45400.701
0.3420.1250.53200.564
0.6730.2870.0350.0050.820
0.6730.2870.0350.0050.821
Table 10. Interaction coefficient between inclusion and slag (ignoring the effect of XCe2O3).
Table 10. Interaction coefficient between inclusion and slag (ignoring the effect of XCe2O3).
SerialXCaOXAl2O3XSiO2XMgOXCe2O3φIS
C00.55310.23710.04840.161300.7603
C10.53380.24420.04980.16610.00610.7588
C20.54450.23490.04950.1650.0060.7542
C30.53860.23330.05030.16760.01020.7499
C40.53130.22920.05140.17140.01670.7424
C50.54610.21650.0510.16990.01660.7361
Table 11. Measured data of interaction coefficient between CaO-SiO2-Al2O3-MgO slag system and steel.
Table 11. Measured data of interaction coefficient between CaO-SiO2-Al2O3-MgO slag system and steel.
XCaOXAl2O3XSiO2XMgOφMS
0.5360.2460.1680.0500.439
0.6040.2770.0680.0510.433
0.4810.2870.1650.0670.432
0.5140.3060.1090.0720.453
0.5440.3240.0560.0760.454
0.5050.2330.1640.0980.438
0.5720.2620.0670.0990.428
0.4660.2150.1600.1590.435
0.5310.2430.0650.1610.429
Table 12. Interfacial properties of the CaO-Ce2O3-SiO2-Al2O3 slag system at 1773 K.
Table 12. Interfacial properties of the CaO-Ce2O3-SiO2-Al2O3 slag system at 1773 K.
SerialσIM (N/m)σIS (N/m)σMS (N/m)Overall Wettability
cosθIMS
A11.5850.27381.40970.9301
A20.27381.40960.9302
A30.27381.40970.9301
A40.27371.40940.9304
A50.27361.40900.9307
Table 13. Interfacial properties of the CaO-Ce2O3-SiO2-Al2O3-MgO slag system at 1823 K.
Table 13. Interfacial properties of the CaO-Ce2O3-SiO2-Al2O3-MgO slag system at 1823 K.
SerialσIM (N/m)σIS (N/m)σMS (N/m)Overall Wettability
cosθIMS
C01.5850.34361.36100.9121
C10.31461.34260.9462
C20.31481.34290.9459
C30.31481.34300.9458
C40.31491.34310.9457
C50.31511.34350.9452
Table 14. Water model parameters required in calculation.
Table 14. Water model parameters required in calculation.
SubstanceDensity
(kg·m−3)
Viscosity
(Pa·s)
Interfacial Tension (N·m)cosθIMS
Water-Oil
σMS
Al2O3-Water
σIM
Al2O3-Oil
σIS
Water9970.0010.0550.6140.635−0.382
Silicon oil9630.096
Al2O3710-
Table 15. Motion behaviors of solid Al2O3 inclusions with different radiuses in slags A1–A5.
Table 15. Motion behaviors of solid Al2O3 inclusions with different radiuses in slags A1–A5.
RadiusSlags with Different Compositions
A4A1A3A2A5
51%CaO
29%Al2O3
10%SiO2
10%Ce2O3
55%CaO
30%Al2O3
10%SiO2
5%Ce2O3
46%CaO
39%Al2O3
10%SiO2
5%Ce2O3
51%CaO
34%Al2O3
10%SiO2
5%Ce2O3
48%CaO
27%Al2O3
10%SiO2
15%Ce2O3
5 μmRRRRR
10 μmRRRRO
20 μmROOOO
50 μmOOOPP
80 μmPPPPP
140 μmPPPPP
Note: R = Remain, O = Oscillate, and P = Pass.
Table 16. Motion behaviors of liquid 50%wtAl2O3–50%wtCaO inclusions with different radiuses in slags A1–A5.
Table 16. Motion behaviors of liquid 50%wtAl2O3–50%wtCaO inclusions with different radiuses in slags A1–A5.
RadiusSlags of Different Compositions
A4A1A3A2A5
51%CaO
29%Al2O3
10%SiO2
10%Ce2O3
55%CaO
30%Al2O3
10%SiO2
5%Ce2O3
46%CaO
39%Al2O3
10%SiO2
5%Ce2O3
51%CaO
34%Al2O3
10%SiO2
5%Ce2O3
48%CaO
27%Al2O3
10%SiO2
15%Ce2O3
5 μmRRRRR
10 μmRRRRR
20 μmRRRRO
40 μmROOOO
50 μmOOOOO
80 μmOOOPP
130 μmOPPPP
Note: R = Remain, O = Oscillate, and P = Pass.
Table 17. Motion behaviors of Al2O3 inclusions with different radiuses in slags C0–C5.
Table 17. Motion behaviors of Al2O3 inclusions with different radiuses in slags C0–C5.
RadiusWithout Ce2O3Slags containing Ce2O3 in Different Compositions
C0C2C1C4C5C3
48%CaO
37.5%Al2O3
4.5%SiO2
10%MgO
46.2%CaO
36.3%Al2O3
4.5%SiO2
10%MgO
3%Ce2O3
45%CaO
37.5%Al2O3
4.5%SiO2
10%MgO
3%Ce2O3
43.4%CaO
34.1%Al2O3
4.5%SiO2
10%MgO
8%Ce2O3
45%CaO
32.5%Al2O3
4.5%SiO2
10%MgO
8%Ce2O3
45%CaO
35.5%Al2O3
4.5%SiO2
10%MgO
5%Ce2O3
5 μmRRROOO
10 μmROOOOO
20 μmROOPPP
30 μmROPPPP
50 μmOPPPPP
80 μmOPPPPP
140 μmPPPPPP
Note: R = Remain, O = Oscillate, and P = Pass.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cao, J.; Li, Y.; Lin, W.; Che, J.; Zhou, F.; Tan, Y.; Li, D.; Dang, J.; Chen, C. Assessment of Inclusion Removal Ability in Refining Slags Containing Ce2O3. Crystals 2023, 13, 202. https://doi.org/10.3390/cryst13020202

AMA Style

Cao J, Li Y, Lin W, Che J, Zhou F, Tan Y, Li D, Dang J, Chen C. Assessment of Inclusion Removal Ability in Refining Slags Containing Ce2O3. Crystals. 2023; 13(2):202. https://doi.org/10.3390/cryst13020202

Chicago/Turabian Style

Cao, Jianqi, Yao Li, Wanming Lin, Julong Che, Feng Zhou, Yunfang Tan, Dongliang Li, Jie Dang, and Chao Chen. 2023. "Assessment of Inclusion Removal Ability in Refining Slags Containing Ce2O3" Crystals 13, no. 2: 202. https://doi.org/10.3390/cryst13020202

APA Style

Cao, J., Li, Y., Lin, W., Che, J., Zhou, F., Tan, Y., Li, D., Dang, J., & Chen, C. (2023). Assessment of Inclusion Removal Ability in Refining Slags Containing Ce2O3. Crystals, 13(2), 202. https://doi.org/10.3390/cryst13020202

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop