ISIJ International
Online ISSN : 1347-5460
Print ISSN : 0915-1559
ISSN-L : 0915-1559
Fundamentals of High Temperature Processes
Kinetics on Formation, Growth, and Removal of Alumina Inclusions in Molten Steel
Katsuhiro Sasai
Author information
JOURNAL OPEN ACCESS FULL-TEXT HTML

2022 Volume 62 Issue 3 Pages 426-436

Details
Abstract

A series of Al deoxidation mechanisms from the nucleation and growth of Al2O3 nuclei immediately after the addition of Al to the growth, agglomeration, and removal of Al2O3 inclusions after the deoxidation equilibrium has been analyzed in light of the kinetics taking into consideration the influences of the interfacial properties on the basis of Al deoxidation experiments of molten steel. The nucleation number density of Al2O3 is (0.72 to 1.62) × 1014 m−3 and increases as the degree of supersaturation increases and the interfacial tension between the nuclei and molten steel decreases. These tendencies can be explained by the homogeneous nucleation theory, and the average interfacial tension, frequency factor, nucleation time, and average nucleation rate are respectively estimated to be 1.43 N·m−1, 4.27 × 1035 m−3·s−1, 0.01 s, and 1.96 × 1016 m−3·s−1 for the nucleation of Al2O3. Al2O3 nuclei rapidly grow to Al2O3 single inclusions having diameters of 2.0 to 2.6 µm through diffusion growth of supersaturated O in molten steel within 2.2 to 3.7 s after the addition of Al, and the molten steel reaches the deoxidation equilibrium. In the subsequent deoxidation equilibrium, the growth rate of Al2O3 single inclusions increases as the O concentration in molten steel increases, and their growth mechanism can be explained by Ostwald ripening. Meanwhile, Al2O3 cluster inclusions grow with the increase in the agglomeration force while agglomerating not only with single inclusions dispersed in molten steel but also with other cluster inclusions existing in the floating paths.

1. Introduction

The deoxidation control implemented in the final process of steelmaking plays an extremely important role in determining the quality and properties of steel products. Thus, there have been many studies to reveal the principle of deoxidation, and these research outcomes were reviewed in detail by, for example, Sakao et al.1) An overview of these studies shows that the thermodynamic data on the deoxidation equilibrium of almost all types of industrially available deoxidizing agents have been measured and evaluated, and have been effectively used for operation analysis in steelmaking based on the equilibrium theory. In contrast, according to research on kinetics, the agglomeration and floating separation of inclusions are considered to be the rate-determining step for deoxidation rates. However, because of the difficulty of conducting dissolution experiments for each elementary step and uncertainty about the influences of the interfacial properties, it seems that there is no common understanding of the mechanisms of the formation and growth of inclusion nuclei immediately after deoxidation as well as the mechanisms of the growth and agglomeration of inclusions after the deoxidation equilibrium in each elementary process.

The author has conducted fundamental studies on the agglomeration mechanism of Al2O3 inclusions in molten steel from the viewpoint of interface chemistry. In the course of these studies, the agglomeration force acting between Al2O3 particles in molten steel was successfully measured and it was verified that Al2O3 inclusions are subject to a strong mutual attracting force generated by the cavity bridge due to their low wettability with molten steel.2,3,4) Also, through Al deoxidation experiments on molten steel, the author quantitatively showed that the particle diameters of Al2O3 cluster inclusions increased in proportion to the agglomeration force generated by the cavity bridge, and that Al2O3 cluster inclusion removal rates also increased as particle diameters increased.5) Based on these fundamental studies, it can be concluded that the Al2O3 inclusions produced in molten steel through Al deoxidation agglomerate and coalesce by the agglomeration force due to the cavity bridge force and are removed as Al2O3 cluster inclusions. However, in order to improve the quality and properties of steel products by controlling the diameters and quantity of Al2O3 inclusions in molten steel, it is still necessary to kinetically ascertain the overall mechanisms of Al deoxidation, including the formation, growth, agglomeration, and removal of Al2O3 inclusions.

In the present study, the formation and growth of Al2O3 nuclei immediately after the addition of Al and the growth, agglomeration, and removal of Al2O3 inclusions in the subsequent deoxidation equilibrium were analyzed through mathematical models derived by assuming appropriate mechanisms on the basis of Al deoxidation experiments that controlled the concentrations of oxygen and sulfur in molten steel, which are interfacial active elements. Through these examinations, the series of Al deoxidation mechanisms were clarified from a perspective of the kinetics in consideration of the influences of the interfacial properties.

2. Experimental Method

2.1. Al Deoxidation Experiments of Molten Steel

The Al deoxidation experiments were conducted using a graphite resistance heating furnace with a capacity of 30 kW.2,3,4,5) In the experiments, 500 g of electrolytic iron (with C concentration of 0.001 mass%, S concentration of 0.0001 mass%, O concentration of 0.005 mass%, and Mn concentration of 0.0001 mass%) was melted in an alumina crucible with an inner diameter of 40 mm and a height of 150 mm under an Ar gas atmosphere. The molten steel temperature was constantly set at 1600°C. After adjusting the O concentration to 0.018 to 0.044 mass% by adding Fe2O3, Al of 0.02 to 0.06 mass% was added to the molten steel based on the stoichiometric proportion of Eq. (1) (reaction between O of 0.018 mass% and Al of 0.02 mass%) so as to adjust the oxygen concentration of total Al2O3 inclusions immediately after the addition of Al to approximately 0.017 mass%, and so as to vary the O concentrations after the deoxidation equilibrium in the range of 0.0006 to 0.0261 mass%. In this study, X represents the dissolved state of the element X.   

2 Al _ +3 O _ =A l 2 O 3 (1)
The experiments also included cases of Al deoxidation with the O concentration and S concentration in molten steel adjusted to 0.0009 mass% or less and 0.018 to 0.073 mass%, respectively. These experiments were started 60 s after the addition of Al and continued for 600 s. The samples of molten steel were taken with transparent quartz tubes with an inner diameter of 6 mm immediately before the addition of Al, after the start of the experiments, and during the experiments; the total oxygen, Al, and S concentrations of the molten steel were analyzed using the samples.

The O concentrations during the experiments were obtained from the analytical values of the Al concentrations using the thermodynamic reevaluation value for the Al deoxidation equilibrium proposed by Itoh et al.6) except for the case of high O concentrations, with Al concentrations less than the analytical limit, in which the O concentrations were calculated using the thermodynamic reevaluation value for the Al deoxidation equilibrium and mass balances based on the O concentrations before the addition of Al and the additive amounts of Al. The oxygen concentrations of total Al2O3 inclusions [I.O]T were obtained by subtracting the O concentrations during the experiments from the analytical values of the total oxygen concentrations at the respective sampling times.

2.2. Observation of Inclusions by the Shapes with Optical Microscope

10 mm long specimens subject to microscopic observation of inclusions were taken from the central portions of rod-like molten steel samples, and the circular cut faces of the specimens were mirror polished. Using an optical microscope, the particle size distribution was obtained for the following: the Al2O3 cluster inclusions having diameters of 10 μm or more existing in the cross-sectional areas with 5 mm diameter excluding the outermost peripheries at 100x magnification; and the Al2O3 single inclusions having diameters of 0.5 μm or more existing in the areas of 1 to 4 mm2 at 1000x magnification. Based on the particle size distribution of the Al2O3 single and cluster inclusions, the average particle diameters and the volume number density of the inclusions in the molten steel samples were calculated using DeHoff’s equation.7) Also, using an EPMA (Electron Probe Micro Analyzer), composition analyses were conducted for typical molten steel samples at each experimental level, and it was confirmed that all inclusions were Al2O3.

2.3. Evaluation of the Formation Number Density of Total Al2O3 Single Inclusions

Because the Al2O3 cluster inclusions are the aggregates of Al2O3 single inclusions, the number density of total Al2O3 single inclusions is considered to be the sum of the number density of the single inclusions suspended in molten steel and that of the single inclusions that form clusters. Because the number density of total Al2O3 single inclusions NV,N (m−3) formed immediately after the addition of Al could not be obtained directly by experiment, it was evaluated by Eq. (2) which corrects the number density of total Al2O3 single inclusions at the start time of the experiments (60 s after the addition of Al) with the ratios of the oxygen concentrations of total Al2O3 inclusions immediately after the addition of Al to those at the start time of the experiments.   

N V,N = [ N V,S(60) + ( d CI(60) / d SI(60) ) 3 ε N V,C(60) ] [I.O] T(0) / [I.O] T(60) (2)
where NV,C is the volume number density of Al2O3 cluster inclusions (m−3), NV,S is the volume number density of Al2O3 single inclusions (m−3), ε is the filling rate of Al2O3 in Al2O3 cluster inclusions with a value of 0.244,5) dCI is the average particle diameter of Al2O3 cluster inclusions (m), dSI is the average particle diameter of Al2O3 single inclusions (m), and the subscripts of (0) and (60) represent the time immediately after the addition of Al and the start of the experiment 60 s after that, respectively.

3. Experimental Results

3.1. Volume Number Density of Al2O3 Nuclei

The degree of supersaturation SS of the Al deoxidation reaction in Eq. (1) is defined by Eq. (3).   

S S = K S / K E (3)
where KS is the solubility product in Eq. (1) expressed by [Al]2·[O]3 and KE is the equilibrium constant of Eq. (1) which uses the reevaluation value of the Al deoxidation equilibrium proposed by Itoh et al.6) [X] represents the X concentration. Figure 1 shows the relation between NV,N and SS. As can be seen, NV,N increased as SS increased. Figure 2 shows the interfacial tension between molten steel and Al2O3, σ A l 2 O 3 -Fe (N·m−1), obtained by substituting into Young’s Eq. (4) the surface tension of molten steel, σFe (N·m−1), and the contact angle between molten steel and Al2O3, θ A l 2 O 3 -Fe (°), noted in previous studies.4,5)   
σ A l 2 O 3 -Fe = σ A l 2 O 3 - σ Fe cos θ A l 2 O 3 -Fe (4)
where σ A l 2 O 3 is the surface tension of Al2O3 with a value of 0.75 N·m−1.8) Then, σ A l 2 O 3 -Fe after the deoxidation equilibrium in each experiment was estimated from Fig. 2, and the relation between σ A l 2 O 3 -Fe and NV,N was obtained as shown in Fig. 3. NV,N demonstrated an increasing trend as σ A l 2 O 3 -Fe decreased, but the influence of σ A l 2 O 3 -Fe was smaller than that of SS in Fig. 1. According to the homogeneous nucleation theory, NV,N is predicted to increase when nucleation is enhanced since the critical nucleus radii become smaller as the supersaturation degree increases and the interfacial tension decreases. Thus, the relation of Figs. 1 and 3 can be explained qualitatively.
Fig. 1.

Relation between degree of supersaturation SS and formation number density NV,N of total Al2O3 single inclusions in Al deoxidation reaction.

Fig. 2.

Influences of the concentrations of O and S in molten steel on interfacial tension σ A l 2 O 3 -Fe between molten steel and Al2O3.

Fig. 3.

Relation between molten steel-Al2O3 interfacial tension σ A l 2 O 3 -Fe and formation number density NV,N of total Al2O3 single inclusions.

When critical nuclei Al2O3 uniformly distributed in molten steel with a volume number density of NV,N independently grow by promptly consuming supersaturated Al and O, Al2O3 single inclusions with a particle diameter d are produced according to Eq. (5).   

d=2 [ 3( C B(0) - C Eq )/(4π N V,N C P ) ] 1/3 (5)
where CB is the O molar concentration in molten steel (mol·m−3), CEq is the equilibrium O molar concentration in molten steel (mol·m−3), and CP is the O molar concentration in Al2O3 calculated by 3 ρ A l 2 O 3 / M A l 2 O 3 = 1.17 × 105 mol·m−3. In the Al deoxidation experiments, because the production amount of Al2O3 single inclusions corresponds to an oxygen mass percentage of 0.017 mass%O, CB(0) − CEq can be calculated by 0.017/100·ρFe/MO = 74.4 mol·m−3, where ρFe is the density of molten steel with a value of 7000 kg·m−3, ρ A l 2 O 3 is the density of Al2O3 with a value of 3970 kg·m−3, MO is the atomic mass of O (kg·mol−1), and M A l 2 O 3 is the molecular mass of Al2O3 (kg·mol−1). Figure 4 shows a comparison of the relation between dSI and NV,N with the calculation results of Eq. (5). At a relatively early stage 60 s after the addition of Al, dSI gradually decreased with the increase in NV,N and roughly corresponded to the particle diameter predicted by Eq. (5). Thus, NV,N can be regarded to be the volume number density of critical nuclei Al2O3. Also, the decrease in dSI with the increase in NV,N was barely observed after 60 s elapsed since the addition of Al; contrarily, dSI showed an increase from the values at 60 s after the addition of Al with the increase in NV,N. Accordingly, it was considered that the particle diameters of Al2O3 single inclusions were dominated by nucleation until approximately 60 s after the addition of Al, but the influence of the nucleation became diluted as time passed thereafter and another particle growth mechanism contributed to the increase in the particle diameters of the Al2O3 single inclusions.
Fig. 4.

Relation between formation number density NV,N of total Al2O3 single inclusions and average particle diameter dSI of Al2O3 single inclusions.

3.2. Particle Diameter Changes and Decrease Rates of Al2O3 Inclusions in Molten Steel

Figure 5 shows the changes in dCI and dSI with time. dCI gradually increased as time advanced before decreasing after peaking at 240 s. The maximum value of dCI was reduced with increases in the concentrations of O and S in the molten steel. This is because the agglomeration force acting among Al2O3 inclusions in the molten steel is reduced with increases in the concentrations of O and S. The details of this matter are discussed later in 4.2.4. Also, dSI showed an increasing tendency as the O concentration increased. This growth mechanism of Al2O3 single inclusions also observed in Fig. 4 is discussed in detail in 4.2.3. The time change in the logarithms of [I.O]T shown in Fig. 6 was a linear decreasing with time according to the first-order kinetics expressed by Eq. (6).   

-d [I.O] T /dt= k T [I.O] T (6)
where t is time (s) and kT is the decrease rate constant for the oxygen concentration of total Al2O3 inclusions (s−1). The decrease rate of [I.O]T was reduced with the increases in the concentrations of O and S in the molten steel. This is because dCI decreased as the concentrations of O and S increased as shown in Fig. 5, thereby reducing the floating separation rates of the Al2O3 cluster inclusions.
Fig. 5.

Time changes in average particle diameter dCI of Al2O3 cluster inclusions and average particle diameter dSI of Al2O3 single inclusions.

Fig. 6.

Time change in oxygen concentration [I.O]T of total Al2O3 inclusions.

4. Discussion

Because when a deoxidizing agent is added to quiescent molten steel the dissolution of the deoxidation agent and its chemical reaction with oxygen are quite fast,9) it is considered that nucleation and the growth of nuclei occur immediately with a rapid reduction in dissolved oxygen to an equilibrium state, and then deoxidation products grow relatively gradually while floating and separating from molten steel.10,11) Also, according to previous Al deoxidation experiments, it is known that when the O concentration increases after reaching the deoxidation equilibrium, Al2O3 single inclusions gradually grow concurrently with the agglomeration growth of Al2O3 cluster inclusions. Thus, in what follows, a series of Al deoxidation mechanisms are examined in a manner that classifies the Al deoxidation process into the following elementary steps: ① nucleation and ② growth of Al2O3 nuclei immediately after the addition of Al, and ③ growth of Al2O3 single inclusions and ④ agglomeration growth and removal of Al2O3 cluster inclusions during floating after reaching the deoxidation equilibrium; and derives rate equations capable of appropriately expressing the time changes in particle number density and particle diameters in each elementary step.

4.1. Al2O3 Nucleation Mechanism ①

4.1.1. Al2O3 Nucleation Theory

Because even a slight decline in supersaturation significantly reduces nucleation rates, Al2O3 nucleation is analyzed by applying the assumption that Bogdandy12) used to calculate the nucleation of ferric oxide vapor to Al deoxidation. That is, Eq. (7) holds when assuming that the nucleation of Al2O3 occurs at a constant nucleation rate I0 (nucleus·m−3·s−1) during a period tk (s) when the initial number density of molecules of supersaturated Al2O3 N0 (molecule·m−3) is reduced to 0.9N0 due to nucleation and, after that, no additional nucleation occurs.   

0.1 N 0 =1/2 I 0 0 t k τ t k wdtdτ (7)
where τ is time (s), w is the growth rate of Al2O3 nuclei (molecule·s−1·nucleus−1), and I0 is the nucleation rate when the number density of molecules of dissolved Al2O3 in supersaturation becomes 0.95N0 (with the degree of supersaturation SS0) that can be obtained by Eqs. (8), (9), (10), (11), (12) according to the homogeneous nucleation theory.13,14)   
I 0 = K V exp(-Δ G N0 * /( k B T)) (8)
  
Δ G N0 * =16/3π σ A l 2 O 3 -Fe 3 [ V A l 2 O 3 /(RTln( S S0 )) ] 2 (9)
  
K V = n Al 2 O 3 * [ σ A l 2 O 3 -Fe /( k B T) ] 1/2 [ 2 v A l 2 O 3 /(9π) ] 1/3 n Al 2 O 3 ( k B T/h)exp(-Δ G A /( k B T)) (10)
  
n Al 2 O 3 * =4π r *2 / { 2 [ 3 v A l 2 O 3 /(4π) ] 1/3 } 2 (11)
  
r * =2 σ A l 2 O 3 -Fe V A l 2 O 3 /(RTln( S S0 )) (12)
where KV is the frequency factor (m−3·s−1), kB is the Bolzmann constant (J·K−1), Δ G N0 * is the change in the free energy of formation of Al2O3 critical nuclei at the degree of supersaturation SS0 (J), T is the absolute temperature of 1873 K, V A l 2 O 3 is the molar volume of Al2O3 (m3·mol−1), R is the gas constant (J·K−1·mol−1), v A l 2 O 3 is the volume of a single molecule of Al2O3 (m3·molecule−1), h is the Planck constant (J·s), ΔGA is the diffusion activation energy with a value of 1.36 × 10−19 J,15) r* is the radius of the critical nucleus of Al2O3 (m), n Al 2 O 3 * is the number of molecules on the surface of the critical nucleus of Al2O3 approximated by Eq. (11), and n A l 2 O 3 is the number of molecules of Al2O3 in the unit volume of molten steel with a value of 1% of N0 ([I.O]T = 0.017 mass%) on the assumption that almost all supersaturated Al2O3 was in a dissolved state at the time of nucleation (with the degree of supersaturation SS0). Given that Al2O3 nuclei grow through steady diffusion of O as with the analysis of Si deoxidation by Sano et al.,11) w in Eq. (7) can be expressed by Eqs. (13) and (14).   
w=4π/(3 v A l 2 O 3 )d ( r *2 +2 k D t) 3/2 /dt (13)
  
k D = D O ( C B0 - C I ) V A l 2 O 3 /3 (14)
where DO is the diffusion coefficient of O in molten steel with a value of 2.3 × 10−9 m2·s−1,16) CB0 is the molar concentration of O in molten steel at the degree of supersaturation SS0 (mol·m−3), and CI is the molar concentration of O on the interface between molten steel and Al2O3 at the molten steel side (mol·m−3), which can be considered equal to CEq on the basis of the deoxidation equilibrium. Equation (15) can be obtained by substituting Eq. (13) into Eq. (7) and integrating it. With Eq. (15) nucleation time tk can be calculated.   
0.1 N 0 =2π I 0 [ ( r *2 +2 k D t k ) 3/2 t k - ( r *2 +2 k D t k ) 5/2 /(5 k D )+ r *5 /(5 k D ) ]/ (3 v A l 2 O 3 ) (15)

4.1.2. Al2O3 Nucleation Rates and Nucleation Time

The volume number density of Al2O3 nuclei can be expressed by the equality NV,N = I0·tk. Equation (16) can be obtained by substituting Eq. (8) into I0 in the equality and taking the logarithm of both sides.   

ln( N V,N )=ln( K V t k )-Δ G N0 * /( k B T) (16)

Assuming that KV and tk are constant under Al deoxidation, a linear relation with an inclination of −1 can be expected between ln(NV,N) and Δ G N0 * /(kB·T).10) Figure 7 shows the relation between ln(NV,N) and Δ G N0 * /(kB·T) obtained for several cases. In the figure, ☐ represents ln(NV,N) corresponding to Δ G N0 * /(kB·T) calculated by Eq. (9) using σ A l 2 O 3 -Fe of each experiment obtained from Fig. 2 in accordance with the concentrations of O and S in molten steel after deoxidation, and ○ represents ln(NV,N) corresponding to Δ G N0 * /(kB·T) calculated by Eq. (9) using σ A l 2 O 3 -Fe = 2.40 N·m−1 that Turpin and Elliott17) applied to the nucleation of Al2O3. The dispersion of ○ with a constant value of σ A l 2 O 3 -Fe was smaller than ☐, but no linear relation could be found between ln(NV,N) and Δ G N0 * /(kB·T). In contrast, ■, ●, and ▲ in the figure represent ln(NV,N) corresponding to Δ G N0 * /(kB·T) calculated by Eq. (9) using the average of σ A l 2 O 3 -Fe obtained from Fig. 2 for each experiment assuming that the deoxidation ratio fDI on the interfaces between molten steel and Al2O3 nuclei were 1.0, 0.95, and 0.9, respectively. Table 1 shows the average σ A l 2 O 3 -Fe for each fDI used in the calculation. fDI is the deoxidation ratio on the interfaces of Al2O3 nuclei defined in the same way as in Eq. (19), and an fDI value of 1 indicates the deoxidation equilibrium. The linear relations with an inclination of −1 with respect to each fDI were calculated through the least square method and drawn as solid lines in the figure. Every case of using the average σ A l 2 O 3 -Fe for each fDI in the calculation shows a linear relation with an inclination of −1 between ln(NV,N) and Δ G N0 * /(kB·T) with less dispersion than the cases of □ and ○. It is difficult to accurately determine the interfacial tension between molten steel and Al2O3 nuclei while taking into consideration the influences of the interfacial active elements in the nucleation process, where ultrafine Al2O3 nuclei form in an extremely short period of time. Using different σ A l 2 O 3 -Fe for each experiment causes the dispersion of the experimental results to increase because, as can be seen in Eq. (9), Δ G N0 * /(kB·T) varies as the cube of σ A l 2 O 3 -Fe . Because the influence of σ A l 2 O 3 -Fe on NV,N is small according to the experimental results shown in Fig. 3, in what follows, a discussion on nucleation is carried out using a constant value for σ A l 2 O 3 -Fe , namely the average, so as to easily obtain the linear relation and alleviate the influence of the differences in σ A l 2 O 3 -Fe .

Fig. 7.

Relation between volume number density NV,N of Al2O3 nuclei and changes in nucleation free energy Δ G N0 * .

Table 1. Average interfacial tension σ A l 2 O 3 -Fe , average critical nucleus radius r*, frequency factor KV, and average nucleation rate I0 with respect to interfacial deoxidation ratio fDI.
Interfacial deoxidation ratio fDIAverage interfacial tension σ A l 2 O 3 -Fe (N·m−1)Average nucleus radius r* (m)Frequency factor KV (m−3·s−1)Average nucleation rate I0 (m−3·s−1)
11.594.87×10−105.56×10351.46×109
0.951.434.38×10−104.27×10351.96×1016
0.91.263.86×10−103.11×10351.84×1022
0.851.183.62×10−102.64×10353.60×1024
0.81.113.40×10−102.26×10352.05×1026

The interceptions of the solid lines in Fig. 7 are equal to ln(KV·tk) as is evident in Eq. (16). Table 1 shows the average r*, KV and the average I0 calculated by Eqs. (8), (9), (10), (11), (12) using the average σ A l 2 O 3 -Fe and the average SS0 of 47000 corresponding to each fDI. The values of tk calculated on the basis of ln(KV·tk) read from Fig. 7 and KV in Table 1 for each fDI are indicated by ○ marks in Fig. 8. Also, the dashed-dotted line in Fig. 8 represents the tk theoretically predicted by Eq. (15) using the averages of r* and I0 shown in Table 1. When fDI was 0.95, the experimental and theoretical values of tk was approximately 0.01 s and both values was coincident. Therefore, the nucleation time of Al2O3 can be considered to be on the order of 0.01 s. On the interfaces between Al2O3 nuclei and molten steel, the O concentration was higher than the equilibrium value by approximately 5% (fDI = 0.95), and σ A l 2 O 3 -Fe was 1.43 N·m−1. Also, as can be seen in Table 1, I0 can be expressed by Eq. (8) when KV is 4.27 × 1035 m−3·s−1, and the average I0 in these experiments was 1.96 × 1016 m−3·s−1. Wakou and Sano18) claimed an I0 of 8.3 × 1015 m−3·s−1 on the basis of the evaluated values of the inclusion number density 1 s after Al deoxidation. This value is close to the average I0 obtained in this study. Based on the above explanation, the discussion about the nucleation of Al2O3 is considered to be appropriate, and also the quantitative results for the interfacial tension, frequency factor, nucleation time, and nucleation rate are considered to be valid.

Fig. 8.

Relation between nucleation time tk and interfacial deoxidation ratio fDI.

4.2. Mechanisms for Growth and Removal of Al2O3 Nuclei and Al2O3 Inclusions

4.2.1. Mathematical Models for Growth and Removal of Al2O3 Nuclei and Al2O3 Inclusions

The following mechanisms are considered to contribute to the growth and removal of Al2O3 inclusion particles in elementary steps ② to ④: (1) growth through diffusion, (2) diffusion growth based on the difference in solubility due to the sizes of inclusion particles (Ostwald ripening), (3) collision agglomeration due to Brownian motion, and (4) collision agglomeration due to the difference in floating rates and floating separation. Below, rate equations based on the respective mechanisms are derived in order to discuss the mechanisms that contribute to the growth and removal of Al2O3 inclusion particles.

(1) Growth through diffusion

The mechanism by which Al2O3 inclusion particles individually grow by consuming supersaturated Al and O with Al2O3 inclusion particles uniformly dispersed in molten steel is discussed in the following. When dividing molten steel into spherical elements with Al2O3 inclusion particles individually positioned at the center of the respective spherical elements, the radii RN (m) of the spherical elements can be expressed by Eq. (17) using the volume number density NV (m−3) of the Al2O3 inclusion particles.   

R N = [ 3/(4π N V ) ] 1/3 (17)
Assuming that deoxidation equilibrium is reached on the interface between Al2O3 and molten steel while the Al2O3 inclusion particles grow through the diffusion of O in the molten steel, the growth of Al2O3 inclusion particles can be calculated by Eqs. (18), (19), (20) derived by Turkdogan19) while taking into consideration the reduction in solute concentrations associated with the growth of deoxidation products. Here, because CP of the Al2O3 inclusion particles is significantly larger than CEq, (CP − CEq) in Eq. (18) is approximated by CP.   
D O t/ R N 2 [ ( C B(S) - C Eq )/( C P - C Eq ) ] 1/3 = 1/6ln[ ( u 2 +u+1)/ (u-1) 2 ]-1/ 3 ta n -1 ( (2u+1)/ 3 ) +1/ 3 ta n -1 ( 1/ 3 ) (18)
  
u 3 = f DB =( C B(S) - C B )/( C B(S) - C Eq ) (19)
  
d=2 R N [ ( C B(S) - C B )/ C P ] 1/3 (20)
where the subscript of (S) represents the origin of time, and u3 and fDB are the deoxidation ratios of molten steel defined by Eq. (19) along with the deoxidation ratios of interfaces. In these equations, the diameters of the Al2O3 inclusion particles at early stages are considered to be negligibly small.

(2) Diffusion growth based on the difference in solubility due to the sizes of inclusion particles (Ostwald ripening)

Regarding Ostwald ripening of precipitated compounds in steel, Hasegawa et al.20) have already analytically derived the growth model of MnS. Using this growth model, the growth rate equations of Al2O3 inclusion particles based on O diffusion rate determination can be obtained as shown in Eqs. (21) and (22).   

d r ˜ /dt= K OW (1/ r ˜ )(1/ r Cr -1/ r ˜ ) (21)
  
K OW =2 σ A l 2 O 3 -Fe D O V A l 2 O 3 /(ηRT)[ C B /( C P - C B ) ] (22)
where r ˜ is the radius of Al2O3 inclusion particles (m), η is the stoichiometric factor of compounds with a value of 3/5 in the case of Al2O3, and rCr is the critical radius at which a particle neither grows nor dissolves and corresponds to the average particle radius. Solving Eq. (21) by applying the Lifshitz and Slyozov theory21) considering the particle diameter distribution gives Eq. (23), which expresses the time changes in average particle diameters due to Ostwald ripening.   
(d/2) 3 - ( d (S) /2) 3 =4/9 K OW t (23)
When η is equal to 1, Eq. (23) is consistent with the equation, capable of being applied to the growth of single precipitates, suggested by Lindborg and Torssell.22)

(3) Collision agglomeration due to Brownian motion

The mechanism assumed in this section is that fine Al2O3 inclusion particles (primary particles) that have initially dispersed uniformly in molten steel agglomerate through two particles collision with each other due to Brownian motion. By applying the agglomeration rate of two particles, which is obtained on the assumption that the Brownian motion is equivalent to the diffusion process of particles, to the balance equation of particle number density, and solving this equation under the initial agglomeration conditions (i.e., the conditions in which the majority of particles are primary particles with nearly identical particle diameters), Eq. (24) is obtained as the rate equation of the number density of Al2O3 inclusion particles.23)   

d N V /dt=- K BR N V 2 (24)
where KBR is the rate constant of Brownian agglomeration (m3·s−1) and expressed by Eq. (25) according to Smoluchowski.24)   
K BR =4 k B T/(3μ) (25)
where μ is the coefficient of viscosity of molten steel with a value of 0.005 Pa·s. Assuming that NV(S) is the initial number density of Al2O3 inclusion particles, Eq. (24) can be solved to obtain Eq. (26).   
N V = N V(S) /(1+ K BR N V(S) t) (26)
Because the total volume of fine particles is constant as they are not removed by flotation, Eq. (26) can be used to obtain Eq. (27), which describes the time changes in particle diameters due to Brownian agglomeration.   
(d/2) 3 - ( d (S) /2) 3 = K BR N V(S) ( d (S) /2) 3 t (27)

(4) Collision agglomeration due to the difference in floating rates and floating separation

i) Collision agglomeration of floating Al2O3 cluster inclusions with Al2O3 single inclusions

The mechanism assumed in this section is that Al2O3 cluster inclusions floating in molten steel grow mainly through agglomeration with fine Al2O3 single inclusions suspended in the molten steel. Given molten steel having a spherical element of volume VE (m3) with one large Al2O3 cluster inclusion located at the center and fine Al2O3 single inclusions suspended around it, because the agglomeration rate of Al2O3 single inclusions onto the Al2O3 cluster inclusion is equal to the decrease rate of Al2O3 single inclusions, the time changes in the oxygen concentration of an Al2O3 single inclusion can be expressed by Eqs. (28) and (29).   

d C SI /dt=- k S C SI (28)
  
k S =( A C / V E ) k m (29)
where CSI is the molar concentration of O included in Al2O3 single inclusions in the molten steel (mol·m−3), AC is the surface area of Al2O3 cluster inclusions (m2), km is the mass transfer coefficient of Al2O3 single inclusions to Al2O3 cluster inclusions at the molten steel side (m·s−1), and kS is the agglomeration rate constant of Al2O3 single inclusions onto floating Al2O3 cluster inclusions and, at the same time, the decrease rate constant of Al2O3 single inclusions (s−1).5) In addition, Eq. (30) holds considering that the growth rate of an Al2O3 cluster inclusion is equal to the agglomeration rate of Al2O3 single inclusions onto the cluster inclusion.   
d( d CI )/dt=2 M A l 2 O 3 k S C SI /[ 3ε ρ A l 2 O 3 ( A C / V E ) ] (30)
By integrating Eq. (30) after substituting CSI obtained from Eq. (28) into it on the grounds that AC/VE during the Al deoxidation experiments are nearly constant,5) Eq. (31) is obtained as an expression of the diameter of Al2O3 cluster inclusion.   
d CI = d CI(S) +2 M A l 2 O 3 C SI(S) [ 1-exp(- k S t) ]/[ 3( A C / V E )ε ρ A l 2 O 3 ] (31)

ii) Collision agglomeration of floating Al2O3 cluster inclusions with Al2O3 single inclusions and other Al2O3 cluster inclusions

The mechanism assumed in this section is that Al2O3 cluster inclusions in molten steel grow through agglomeration with other Al2O3 cluster inclusions during floating in addition to the agglomeration with Al2O3 single inclusions. Given that the floating rate vC of Al2O3 cluster inclusions in molten steel (m·s−1) follows Stokes’ law expressed by Eq. (32), the growth rate of Al2O3 cluster inclusions through agglomeration with other Al2O3 cluster inclusions in the volume of molten steel in the floating path is given by Eq. (33).11)   

v C =2g( ρ Fe - ρ C ) ( d CI /2) 2 /(9μ) (32)
  
d( d CI )/dt= α C v C /2= α C g( ρ Fe - ρ C ) d CI 2 /(36μ) (33)
where g is gravitational acceleration (m·s−2) and ρC is the density of Al2O3 cluster inclusions (kg·m−3), which can be expressed by Eq. (34) using the filling rate of Al2O3 cluster inclusions.   
ρ C =ε ρ A l 2 O 3 +(1-ε) ρ Fe (34)
where αC is the volume of Al2O3 cluster inclusions in the unit volume of molten steel that can be expressed by Eq. (35) using the initial value of the oxygen concentration [I.O]C in Al2O3 cluster inclusions in the molten steel.   
α C = ρ Fe M A l 2 O 3 [I.O] C(S) /(300 M O ε ρ A l 2 O 3 ) (35)
Because the growth of Al2O3 cluster inclusions occurs by the mechanisms of both agglomeration with Al2O3 single inclusions (Eq. (30)) and agglomeration with other Al2O3 cluster inclusions (Eq. (33)), Eq. (36) holds.   
d( d CI )/dt=[ α C g( ρ Fe - ρ C ) /(36μ) ] d CI 2 +2 M A l 2 O 3 k S C SI(S) /[ 3ε ρ A l 2 O 3 ( A C / V E ) ] = A 4 d CI 2 + B 4 (36)
Here, CSI in Eq. (30) is approximated by CSI(S) for simplification. Integrating Eq. (36) obtains Eq. (37) as an expression of the diameter of Al2O3 cluster inclusion.   
d CI = tan( ( A 4 B 4 ) 0.5 t+ta n -1 ( ( A 4 / B 4 ) 0.5 d CI(S) ) ) / ( A 4 / B 4 ) 0.5 (37)

iii) Floating separation rates of total Al2O3 single inclusions

In Al deoxidation molten steel, Al2O3 single inclusions having diameters of a few μm coexist with coarse Al2O3 cluster inclusions having diameters of a few dozen μm formed through the agglomeration of Al2O3 single inclusions. Thus, [I.O]T can be expressed by Eq. (38) using the volume number density NV,T (m−3) of total Al2O3 single inclusions combining the single and cluster inclusions.   

[I.O] T =100(4π/3) ( d SI /2) 3 N V,T 3 M O ρ A l 2 O 3 /( M A l 2 O 3 ρ Fe ) (38)
By substituting Eq. (38) into Eq. (6), the removal rate of total Al2O3 single inclusions is given by Eq. (39).   
-d N V,T /dt= k T N V,T (39)

4.2.2. Growth Mechanism of Al2O3 Nuclei ②

Inclusion nuclei grow very quickly. It is reported that the diameters of inclusion nuclei increase to 0.3 to 0.8 μm in about 1 second after addition of a deoxidizing agent in the case of Al deoxidation.18) In the case of Si deoxidation, the diameters reportedly reach 1.4 μm10) concurrently with a rapid decrease in dissolved oxygen to equilibrium values.9,10) As is clear from the derivation of the growth rate equations explained above, the process of rapid growth of Al2O3 nuclei while consuming supersaturated Al and O can be explained by the diffusion growth model expressed by Eqs. (17), (18), (19), (20). Figure 9 shows the time changes in fDB and dSI calculated by the diffusion growth model. The solid and dotted lines in the figure respectively represent the results calculated using the maximum Al2O3 nucleus volume number density (NV,N = 1.62 × 1014 m−3) and the minimum Al2O3 nucleus volume number density (NV,N = 7.15 × 1013 m−3) obtained from experiments in terms of NV. Also, these time changes are based on CP − CEq ≈ CP = 1.17 × 105 mol·m−3 and CB(S) − CEq = 74.4 mol·m−3. As can be seen in the figure, dSI increased along with the rapid decrease in supersaturated O immediately after the addition of Al, and grew to 2.0 to 2.6 μm after 2.2 to 3.7 s had elapsed (fDB = 0.9999) since the addition of Al depending on the number density of Al2O3 nuclei. These particle diameters roughly correspond to the experimental values 60 s after the addition of Al. Accordingly, the Al2O3 nuclei that formed immediately after the addition of Al are considered to have grown rapidly due to the diffusion of O in molten steel.

Fig. 9.

Changes in deoxidation ratio fDB and particle diameters dSI of Al2O3 single inclusions associated with diffusion growth of Al2O3 nuclei.

4.2.3. Growth Mechanism of Al2O3 Single Inclusions ③

As can be seen in Fig. 5, Al2O3 single inclusions grew as the O concentration in the molten steel increased. Considering that such growth occurred after reaching the deoxidation equilibrium 60 s after the addition of Al, the growth mechanism cannot be (1) growth due to diffusion. Also, because Al2O3 single inclusions have particle diameters of up to approximately 3 μm and are scarcely able to float according to Stokes’ law, (4) collision agglomeration due to the differences in floating rates can be excluded from the possible growth mechanisms of Al2O3 single inclusions. So, either (2) Ostwald ripening or (3) collision agglomeration due to Brownian motion can be considered to be the growth mechanism of Al2O3 single inclusions. Figure 10 shows the particle diameters of Al2O3 single inclusions based on the growth rate equations obtained as Eqs. (23) and (27). Despite some dispersion, linear relations can be found between (dSI/2)3 − (dSI(S)/2)3 and time as indicated by the solid lines in the figure. Figure 11 shows the influences of the concentrations of O and S on the inclinations of the lines (= [(dSI/2)3 − (dSI(S)/2)3]·t−1) obtained from Fig. 10. The dotted line represents the inclination (=KBR·NV,S(S)·(dSI(S)/2)3) of the Brownian agglomeration obtained on the basis of Eqs. (25) and (27) using the average volume number density of 6.42 × 1013 m−3 and the average particle diameter of 2.10 μm at the start of the experiments. The solid lines are the inclinations of Ostwald ripening (=4/9KOW) calculated from Eqs. (22) and (23) using σ A l 2 O 3 -Fe of each experiment in Fig. 2 in accordance with the concentrations of O and S in the molten steel. The inclination of Brownian agglomeration was close to the experimental values for the cases with an O concentration of 0.0009 mass% or less, but smaller than the other experimental values for the cases with a higher O concentration because Brownian agglomeration does not depend on O concentration. Thus, although it cannot be completely excluded from the possible growth mechanisms of Al2O3 single inclusions, (3) the collision agglomeration due to Brownian motion cannot be considered to be the primary one. When including the dependency on the concentrations of O and S, the inclination of Ostwald ripening roughly corresponds to the experimental results. Considering that the influences of (3) the collision agglomeration due to Brownian motion on the inclination are small, it is considered that the growth mechanism of Al2O3 single inclusions can be explained mainly by (2) diffusion growth based on the difference in solubility due to the sizes of inclusion particles.

Fig. 10.

Relation between (dSI/2)3 − (dSI(S)/2)3 and time.

Fig. 11.

Influences of concentrations of O and S on [(dSI/2)3 − (dSI(S)/2)3]·t−1.

4.2.4. Agglomeration Growth and Removal Mechanisms of Al2O3 Cluster Inclusions ④

The author has experimentally and theoretically verified that Al2O3 inclusions in molten steel are subject to the action of the agglomeration force due to the cavity bridge force, which is far stronger than the van der Waals force and the capillary force on the surfaces of molten steel,2,3,4) and that agglomeration rates increase as the agglomeration force increases, thereby increasing the particle diameters of Al2O3 cluster inclusions.5) Figure 12 shows the relation between the maximum particle diameters dCI,Max (corresponding to dCI after 240 s in Fig. 5) of Al2O3 cluster inclusions of each experiment and the agglomeration force FA,S (acting between spherical Al2O3 inclusions with dSI = 2 μm) due to the cavity bridge force reported in the previous paper.5) Also, for comparison purposes, the particle diameters of Al2O3 cluster inclusions calculated by Eqs. (31) and (37) are plotted in Fig. 12 as □ and ◇, respectively. In the calculation of the agglomeration model, the ratio of the oxygen concentration in Al2O3 cluster inclusions to the oxygen concentration of total Al2O3 inclusions ([I.O]T = [I.O]C + [I.O]S, [I.O]S: oxygen concentration in Al2O3 single inclusions in molten steel) is 0.533,5) and the average [I.O]T(S) is 0.0152 mass%. Therefore, CSI(S) is 0.0152·(1 – 0.533)/100·ρFe/MO = 31.1 mol·m−3, and αC is 0.0152·0.533·ρFe· M A l 2 O 3 /(300MO·ε· ρ A l 2 O 3 ) = 1.24 × 10−3 according to Eq. (35). In the previous paper,5) the experimental values of kS were obtained by applying the first-order kinetics, Eq. (28), to the changes in the oxygen concentration of Al2O3 single inclusions in molten steel, and the appropriateness of the experimental values has been verified with a mass transfer model using the penetration theory. Also, AC/VE during experiments can be considered to have a constant value of approximately 195 m−1 5) due to mutual compensation of the influences of particle coarsening and floating separation, and dCI(S) was determined to be 24.9 μm, which was the average of the actual measurements. According to the favorable linear relation between dCI,Max and FA,S, it is obvious that the agglomeration growth of Al2O3 cluster inclusions is strongly affected by the agglomeration force due to the cavity bridge force acting between Al2O3 inclusions. Also, dCI,Max is consistent with the agglomeration model involving single inclusions only when FA,S is relatively small. However, when FA,S becomes larger, dCI,Max increases and is more consistent with the agglomeration model that involves both cluster and single inclusions than the model that only involves single inclusions. Accordingly, with the increase in the agglomeration force between inclusions, it is considered that Al2O3 cluster inclusions were removed by flotation while growing such that they not only agglomerated with suspended fine single inclusions but also strongly attracted and agglomerated with other cluster inclusions existing in the floating paths.

Fig. 12.

Comparison of maximum Al2O3 cluster inclusion diameters dCI,Max and calculated Al2O3 cluster inclusion diameters using the agglomeration model.

4.3. Influences of Deoxidation Mechanisms in Each Elementary Step on the Formation, Growth, and Removal of Al2O3 Inclusion Particles

In 4.1 and 4.2, the Al deoxidation process in molten steel was divided into several elementary steps, and the mechanisms of formation, growth, and removal of Al2O3 inclusion particles in each elementary step were clarified through proposals of mathematical models based on said mechanisms. In this section, to coherently verify a series of mechanisms, the changes in the diameters and volume number density of Al2O3 inclusion particles are calculated from immediately after addition of Al to after the deoxidation equilibrium by using the aforementioned mathematical models while taking into consideration the continuity of the elementary steps. Figures 13 and 14 show comparisons of the calculation results and experimental values of the diameters and volume number density of Al2O3 inclusion particles, respectively. Because the growth, agglomeration, and removal of Al2O3 inclusion particles after reaching deoxidation equilibrium change in accordance with the O concentration in the molten steel, calculated values with a low O concentration (0.0006 mass%) and high O concentration (0.0261 mass%) are shown as solid and dashed-dotted lines in Figs. 13 and 14, respectively. In these calculations, it was assumed in the calculation that the rapid agglomeration and removal of Al2O3 inclusions due to the stirring of the molten steel by the addition of Al occurred instantaneously immediately after reaching the Al deoxidation equilibrium, and almost no further changes occurred until 60 s elapsed thereafter. Thus, immediately after reaching the deoxidation equilibrium, the diameter of Al2O3 cluster inclusions was the average diameter at the start of the experiments, and the number density of total Al2O3 single inclusions was 89% of that immediately after the addition of Al (= [I.O]T(60)/[I.O]T(0) = 0.0152/0.017). Also, the experimental values for the number density of total Al2O3 single inclusions were obtained by converting the actually measured [I.O]T with Eq. (38), which uses the average particle diameter of Al2O3 single inclusions of 2.10 μm as dSI instead of using microscopic observation results with large dispersion.

Fig. 13.

Time changes in diameters of Al2O3 inclusion particles.

Fig. 14.

Time changes in total volume number density of Al2O3 inclusion particles.

During the period from immediately after the addition of Al until tk = 0.01 s elapsed once nucleation was completed with Al2O3 nuclei having a diameter of 8.8 × 10−4 μm formed at I0 = 2.0 × 1016 m−3·s−1, the volume number density of Al2O3 nuclei reached 2.0 × 1014 m−3, which is close to the experimental value range of NV,N = 0.72 to 1.62 × 1014 m−3. The formed Al2O3 nuclei then consumed all supersaturated O and, after 1.9 s, equilibrium was reached through the diffusion growth of the Al2O3 nuclei into Al2O3 single inclusions having a diameter of 1.83 μm. During this period, considering the very slow floating rates of fine Al2O3 single inclusions, it was assumed that Al2O3 single inclusions did not undergo any change in volume number density. In the equilibrium state after 1.9 s from the addition of Al, Al2O3 single inclusions and Al2O3 cluster inclusions were removed while concurrently growing through different mechanisms. The Al2O3 single inclusions with a diameter of 1.83 μm grew through Ostwald ripening with a growth rate determined by O diffusion, and the particle diameter increased as the O concentration increased, which corresponded well with the experimental values. Meanwhile, the Al2O3 cluster inclusions with a diameter of 24.9 μm which rapidly agglomerated due to stirring of the molten steel in association with the addition of Al were removed by flotation while agglomerating single inclusions when the O concentration was high, and both single inclusions and other cluster inclusions when the O concentration was low. Thus, the diameters of Al2O3 cluster inclusions increased as the O concentration decreased and with the acceleration of the decrease in the number density of total Al2O3 single inclusions accordingly. All these calculation results roughly reproduced the experimental values.

As explained above, the Al deoxidation mechanisms proposed in this research were verified to be appropriate because the mathematical models of the formation, growth, and removal of Al2O3 inclusion particles derived by appropriately assuming mechanisms for each elementary step could coherently reproduce the Al deoxidation behavior from immediately after the addition of Al to after the deoxidation equilibrium. In addition, by the Al deoxidation mechanisms elucidated in this study, it is possible to optimize the rate controlling factors such as the degree of supersaturation and O, which is the interfacial active element, as a whole in consideration of the effects of the factors on each elementary step. As a result, the diameter and quantity of Al2O3 inclusion can be controlled at a high lebel, and it is considered that this can greatly contribute to the improvement of the quality and properties of steel products.

5. Conclusions

The following conclusions were obtained as a result of analyses of a series of Al deoxidation behavior from nucleation and growth of Al2O3 nuclei immediately after the addition of Al to the growth, agglomeration, and removal of Al2O3 inclusions after reaching deoxidation equilibrium in light of the kinetics on the basis of Al deoxidation experiments in which the interfacial tension between molten steel and Al2O3 inclusions was controlled by changing the concentrations of O and S in the molten steel.

(1) The nucleation number density of Al2O3 was (0.72 to 1.62) × 1014 m−3 and increased as the degree of supersaturation increased and the interfacial tension between the nuclei and molten steel decreased. However, the influence of the interfacial tension was relatively small. These tendencies can be explained by the homogeneous nucleation theory and the interfacial tension, frequency factor, nucleation time, and average nucleation rate in the Al deoxidation experiments were estimated to be 1.43 N·m−1, 4.27 × 1035 m−3·s−1, 0.01 s, and 1.96 × 1016 m−3·s−1, respectively.

(2) The growth mechanism of Al2O3 nuclei can be explained by the diffusion growth of supersaturated O in molten steel, and the critical nuclei rapidly grew to Al2O3 single inclusions having diameters of 2.0 to 2.6 μm within 2.2 to 3.7 s after addition of Al in accordance with the number density of Al2O3 nuclei and the deoxidation equilibrium was reached.

(3) In the deoxidation equilibrium, Al2O3 single inclusions and Al2O3 cluster inclusions concurrently grew and were removed through different growth mechanisms.

(4) The growth rate of Al2O3 single inclusions increased as the O concentration in the molten steel increased, and their growth mechanism can be explained by diffusion growth (Ostwald ripening) based on the difference in solubility due to the sizes of the inclusion particles.

(5) Al2O3 cluster inclusions floated and separated while growing along with the increase in the agglomeration force originating from the cavity bridge force such that they not only agglomerated with fine single inclusions dispersed in the molten steel but also attracted and agglomerated with other cluster inclusions existing in the floating paths.

(6) The Al deoxidation mechanisms proposed in this research were verified to be appropriate because the mathematical models of the formation, growth, agglomeration, and removal of Al2O3 inclusions derived by appropriately assuming mechanisms for each elementary step could coherently reproduce the Al deoxidation behavior from immediately after the addition of Al to after reaching the deoxidation equilibrium.

References
 
© 2022 The Iron and Steel Institute of Japan.

This is an open access article under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs license.
https://creativecommons.org/licenses/by-nc-nd/4.0/
feedback
Top