Computational Study of the Transport Properties of Molten CaO-Si[O.sub.2]-[P.sub.2][O.sub.5]-FeO System.
The transport properties (including diffusion coefficient, viscosity, electrical, and thermal conductivity, etc.) of liquids or melts are important for most physical, chemical and biotechnological processes. In engineering practices, such properties must be gained through measurement, theoretical calculation, or even guessing.
It is well known that CaO-Si[O.sub.2]-[P.sub.2][O.sub.5] system is one of the most important melts in metallurgical, glass and ceramic industries [1-5]. The quaternary CaO-Si[O.sub.2]-[P.sub.2][O.sub.5]-FeO melts have played a crucial role in the steel metallurgical process. The transport properties of CaO-Si[O.sub.2]-[P.sub.2][O.sub.5]-FeO system are critical for understanding the chemical reactions at the interface between molten slag and liquid metal.
Among a large number of transport properties, the diffusion coefficient and viscosity are of the greatest practical and scientific importance and are associated with the transport of mass and momentum, respectively. The diffusion coefficient in molten slag refers to a specific ions rather than to the bulk melt. It describes the rate of chemical or isotopic transport due to thermally activated random walk of atomic species. In the absence of a chemical gradient, we call it self-diffusion coefficient. This parameter cannot be measured without introducing an isotopic gradient. However, in theoretical calculation, such as molecular dynamics (MD) simulations, an isotopie gradient is not required as the motion of individual particles can be traced .
Viscosity of molten slag describes the resistance of fluid against its deformation by a shear stress and can be measured by concentric cylinder viscometer. However, in some metallurgical process, the temperature is always well above the experimental conditions. In this case, viscosity of molten slag also can be computed from MD simulations using the Green-Kubo relation [7, 8].
In the present study, we performed a MD simulation and emphatic studied the effect of slag composition on the transport properties of molten CaO-Si[O.sub.2]-[P.sub.2][O.sub.5]-FeO system. By tracing the mean square displacements of atoms in this melts, MD simulation in the present work not only obtained the self-diffusion coefficient values but also provide an intrinsically link between the viscosity and the structure of the molten system.
Classical molecular dynamics (CMD) simulations were employed to calculate the transport properties and structural information of the quaternary system. All calculations were performed with the parallel package Materials Explorer 5.0 using the Born-Mayer-Huggins (BMH) function [9, 10]:
[mathematical expression not reproducible] (1)
where U([r.sub.ij]) is the interatomic potential; [q.sub.i] and [q.sub.j] are the charges of ion i and ion j; [r.sub.ij] represents the distance between ion i and ion j; [A.sub.ij], [B.sub.ij] and [C.sub.ij] are parameters for BMH potentials.
The charge of Ca, Si, P, Fe and 0 are 2, 4, 5, 2 and -2, respectively. The coulombic potential in eq. (1) (the first term) was calculated by using the Ewald summation. All the noncoulombic parts of the pair potentials were subjected to a short-range cutoff of 10 A. Based on the data from Hirao and Kawamura , the interatomic potential parameters were optimized using Lammps program and shown in Table 1. The slag composition and total number of atoms are detailed in Table 2. In order to ensure the melts is a homogeneous molten slag, the compositions are designed based on the CaO-Si[O.sub.2]-[P.sub.2][O.sub.5]-FeO phase diagram and all the slags are located in the liquid phase region at 1,673 K .
In the calculation, the Newton motion equations were numerically integrated with a step of 1 fs. The initial temperature was set at 4,000 K for 15,000 steps to eliminate the effect of the initial distribution. Then, the temperature was decreased slowly to 2,000 K, 1,873 K and 1,673 K for 10,000 steps, respectively. The calculation was equilibrated at 1,673 K for another 25,000 steps. All simulations were performed in NVT (constant number of particles, volume and temperature) ensemble. The self-diffusion coefficients D of the ions in the molten system can be estimated from the mean squared displacements (MSD) by using following relation ;
[mathematical expression not reproducible] (2)
[mathematical expression not reproducible] (3)
where [r.sub.i](O) and [r.sub.i](t) denote the position of the ion i at time t and time origin, respectively. N is the number of the ion i.
The shear viscosity [eta] of the molten slag system can be calculated by Einstein-Stokes equation .
[eta] = [k.sub.B]T / D[lambda] (4)
where [k.sub.B] is the Boltzmann constant. T is temperature in Kelvin. [lambda] is the diffusion step, which is usually considered to be equal to the diameter of the ion.
Results and discussion
MSD and self-diffusion coefficients
Figure 1 shows the MSD of Ca, Si, P, Fe and 0 ions calculated as a function of time in the No. 2 slag system at 1,673 K. The self-diffusion coefficients of each ion can be obtained from the slope of the linear part of the MSD vs. time curve. The self-diffusion coefficients obtained at 1,673 K of Ca, Si, P, Fe and O ions are 0.056 [A.sup.2]/ps, 0.023 [A.sup.2]/ps, 0.032 [A.sup.2]/ps, 0.053 [A.sup.2]/ps and 0.045 [A.sup.2]/ps, respectively. The results show that the Ca and Fe ions diffuse more rapidly than other ions in the CaO-Si[O.sub.2]-[P.sub.2][O.sub.5]-FeO system. The O ion is somewhat less mobile, and the Si and P ions have the lowest mobility. As Table 3 shows, the diffusivity sequence of Ca>O>Si is also reported in molten CaO-Si[O.sub.2], CaO-MgO-[Al.sub.2][O.sub.3]-Si[O.sub.2] and CaO-Ca[F.sub.2]-Si[O.sub.2] slag systems [15-18]. This suggests the restricted diffusion of the Si ion within a Si[O.sub.4] tetrahedron. To our knowledge, there are few experimental or computational diffusion studies for the CaO-Si[O.sub.2]-[P.sub.2][O.sub.5]-FeO system. Although very limited data are available in the literatures, the Fe ion is more diffusive than O ion, but much less diffusive than Ca ion is confirmed [19, 20]. The self-diffusion coefficient of P is less than 0, but larger than Si, which indicating its weaker bonding in tetrahedral coordination. The self-diffusion coefficients in present work are on the same order and magnitude as those reported molten oxide systems.
Table 4 shows all the self-diffusion coefficients of Ca, Si, P, Fe and O ions in Nos. 1-10 slags. It can be seen that with the slag basicity increases from 0.6 to 1.5, the self-diffusion coefficients of Ca, Si, P, Fe and O ions increase from 0.04 to 0.065 [A.sup.2]/ps, 0.017 to 0.027 [A.sup.2]/ps, 0.023 to 0.038 [A.sup.2]/ps, 0.037 to 0.057 [A.sup.2]/ps and 0.033 to 0.053 [A.sup.2]/ps, respectively. The self-diffusion coefficients of all ions increase with increasing slag basicity. The results represent reasonably well the structural properties of molten CaO-Si[O.sub.2] system, which were restricted to the diffusion of Si and O ions by the formation of silica network structure . With the increase of [P.sub.2][O.sub.5] content from 5 to 14%, the self-diffusion coefficients of Ca, Si, P, Fe and O ions decrease from 0.097 to 0.043 [A.sup.2]/ps, 0.04 to 0.018[A.sup.2]/ps, 0.057 to 0.025 [A.sup.2]/ps, 0.092 to 0.042 [A.sup.2]/ps and 0.08 to 0.035 [A.sup.2]/ps, respectively. As well as Si, increasing the amount of P ions in oxide melts leads to a more close slag structure. Therefore, the polymerization degree of the melts increased and all the self-diffusion coefficients decreased. On the contrary, with the increase of FeO content from 30 to 40%, the self-diffusion coefficients of Ca, Si, P, Fe and O ions increase from 0.056 to 0.097 [A.sup.2]/ps, 0.023 to 0.04 [A.sup.2]/ps, 0.032 to 0.057 [A.sup.2]/ps, 0.053 to 0.092 [A.sup.2]/ps and 0.045 to 0.08 [A.sup.2]/ps, respectively. This is because the Fe ion plays the role of network modifier in the oxide melts. The depolymerization of the slag structure decreases with increasing FeO content, thus, the self-diffusion coefficients increased.
Table 5 shows the viscosities of Nos. 1-10 slags calculated by self-diffusion coefficients. It can be seen that the viscosities decrease from 0.26 to 0.15 Pa-s with the slag basicity increasing from 0.6 to 1.5. Similarly the viscosities decrease from 0.18 to 0.11 Pa-s with the FeO content in slag increasing from 30 to 40%; and increase from 0.14 to 0.28 Pa-s with the [P.sub.2][O.sub.5] content in slag increasing from 5 to 14 %. It is apparent that increasing the amount of network formers (e. g., Si and P ions) in molten CaO-Si[O.sub.2]-[P.sub.2][O.sub.5]-FeO system leads to larger viscosity while adding network modifiers (e. g., Fe and Ca ions) causes viscosity to decrease.
In order to verify the validity of the calculated results, we measured the viscosities of above molten slags by the rotating cylinder method . In the viscosity measurement test, both the crucible and spindle are made of pure iron. The comparisons of the measured and the calculated values for the viscosities are shown in Figure 2. The calculated values are in good agreement with the experimental ones.
Links between transport and structural properties
Each Si and P tetrahedron was classified into five types of structural units, i. e. [Q.sup.0], [Q.sup.1], [Q.sup.2], [Q.sup.3], and [Q.sup.4], where the superscript referred to the number of bridging oxygen atoms in the unit . Figure 3 shows the concentration variation of [Q.sup.n] with changing slag basicity, FeO content and [P.sub.2][O.sub.5] content, respectively. It is important to note that the [Q.sup.n](Si + P) in present work includes both [Q.sup.n](Si) and [Q.sup.n](P). Seen from Figure 3(a), with the increase of slag basicity, the proportion of [Q.sup.0](Si + P) shows a remarkable uptrend while the proportions of [Q.sup.2](Si + P), [Q.sup.3](Si + P) and [Q.sup.4](Si + P) decrease obviously. It suggests that the degree of polymerization of the Si-O and P-O network structures is weakened. This phenomenon is also observed with the increase of FeO content in slag. At the same time, the proportions of [Q.sup.0](Si + P) and [Q.sup.1](Si + P) decrease slightly with the increase of [P.sub.2][O.sub.5] content in slag while the proportions of [Q.sup.2](Si + P) and [Q.sup.3](Si + P) increase slightly. It can be concluded that the network of both P and Si ions are more polymerized and tend to form complex anions.
The ratio of non-bridging oxygen/tetragonal oxygen (NBO/T) is a measure of the de-polymerization of the molten system. Mills et al.  proposed a relational expression between the polymerization of the molten CaO-[Na.sub.2]O-[Al.sub.2][O.sub.3]-Si[O.sub.2] slag and the mole fractions of different oxides.
[mathematical expression not reproducible] (5)
Q = 4 - NBO/T (6)
Similarly, the parameter Q(Si + P), which is a measure of the polymerization of the molten CaO-Si[O.sub.2]-[P.sub.2][O.sub.5]-FeO system and which can be calculated by
[mathematical expression not reproducible] (7)
Figure 4 shows the relationship between the calculated viscosities and the parameter of Q(Si + P) at 1,673 K. It can be intuitively seen that the viscosities vs Q(Si + P) is a linear-related function. However, there is a special point that the deviation between the calculated data and the fitting data is pretty huge. The reason for this may be that the [P.sub.2][O.sub.5] content is extremely high, which going beyond the application of eq. (7). Without thinking about this special point, the coefficient of determination [R.sup.2] increases from 0.65 to 0.92.
FT-IR spectra analysis
The structures of the slags were also detected by FT-IR spectroscopy (Nicolet 5DXC). Figure 5(a) shows the effect of slag basicity on the FT-IR spectra of the slags. In general, the spectra was composed of three regions, i. e. 800-1,200 [cm.sup.-1] with strong intensity, 600-800 [cm.sup.-1] with weak intensity and 400-600 [cm.sup.-1] with middle intensity. The band in the 800-1,200 [cm.sup.-1] range has been generally assigned to the symmetric stretching vibration of Si[O.sub.4] and P[O.sub.4] tetrahedra. The bands shift to lower wavenumber region with increasing slag basicities, which suggests that the proportions of [Q.sup.2](Si + P) and [Q.sup.3](Si + P) decreased. Therefore, the degree of polymerization of the slags decreased . With the increase of FeO content, this phenomenon also showed in Figure 5(b). However, the bands shift from low wavenumber to high wavenumber in Figure 5(c), which indicates that more polymerized structures were formed with increasing [P.sub.2][O.sub.5] content. The weak 600-800 [cm.sup.-1] bands represent the Si-O symmetry stretching [26, 27]. The characteristic peak near 695 [cm.sup.-1] corresponds to the Fe[O.sub.4] tetrahedra stretching bands . The medium intense band in the 400-600 [cm.sup.-1] range represents the bending vibration of T-O-T (T denotes Si or P).
In the present theoretical study, the MSD of Ca, Si, P, Fe and 0 ions in quaternary CaO-Si[O.sub.2]-[P.sub.2][O.sub.5]-FeO melts have been modeled by MD simulations. The self-diffusion coefficients of above ions and the viscosities of the melts were calculated. It was found that the diffusivity of above ions follow the sequence of Ca>Fe>O>P>Si. The self-diffusion coefficients of all ions increase with increasing slag basicity and FeO content, while decrease with increasing [P.sub.2][O.sub.5] content.
The calculated viscosities are in good agreement with the experimental values. It is apparent that increasing the amount of network formers (e. g., Si and P ions) in molten CaO-Si[O.sub.2]-[P.sub.2][O.sub.5]-FeO system leads to larger viscosity while adding network modifiers (e. g., Fe and Ca ions) causes viscosity to decrease.
The links between transport properties and structural properties of the melts were also discussed. The degree of polymerization of the Si-O or P-O network structures weakened with increasing slag basicity and FeO content, while promoted with increasing [P.sub.2][O.sub.5] content. The FT-IR spectra of the slags also show the same structural informations. It is interesting to note that the viscosity vs Q (Si + P) shows a linear-related function.
We expect that this investigation of the transport properties of CaO-Si[O.sub.2]-[P.sub.2][O.sub.5]-FeO system will encourage and motivate further experimental and theoretical investigations, and finally reveal the mechanism of viscosity change with slag composition.
Funding: This work was supported by the National Basic Research Program of China under Grant No. 2013CB632604; Fundamental Research Funds for the Central Universities under Grant CDJZR 14130001; and National Natural Science Foundation of China under Grants Nos. 51090382 and 51604047.
 I. Kansal, A. Goel, D.U. Tulyaganov, R.R. Rajagopal and J.M.F. Ferreira, J. Eur. Ceram. Sco., 32 (2012) 2739-2746.
 U. Hoppe, G. Walter, R. Kranold and D. Stachel, J. Non-Cryst. Solids, 263&264 (2000) 29-47.
 G. Fan, J. Diao, L. Jiang, Z. Zhang and B. Xie, Mater. Trans., 56 (2015) 655-660.
 J. Diao, G.Z. Fan, X. Liu and B. Xie, Metall. Mater. Trans. B, 45 (2014) 1942-1947.
 C. Mercier, C. Follet-Houttemane, A. Pardini and B. Revel, J. Non-Cryst. Solids, 357 (2011) 3901-3909.
 H.W. Ni, H.J. Hui and G.S. Neumann, Rev. Geophys., 53 (2015) 715-744.
 G.B. Martin, F.J. Spera, M.S. Ghiorso and D. Nevins, Am. Mineral, 94 (2009) 693-703.
 N. Sun, L Stixrude, N. Koker and B.B. Karki, Geochim. Cosmochim. Ac, 75 (2011) 3792-3802.
 L. Woodcock, C. Angeli and P. Cheeseman, J. Chem. Phys., 65 (1976) 1565-1577.
 H. Luo, S. Xiao, S. Wang, P. Huai, H. Deng and W. Hu, Comp. Mater. Sci., 111 (2016) 203-208.
 H. Hirao and K. Kawamura, Material Design Using Personal Computer, Shokabo, Tokyo (1994).
 Verein Deutscher Eisenhuttenleute, Slag Atlas, 2nd edition, Verlag Stahleisen mBH, Dusseldorf (1995).
 D.J. Lacks, D.B. Rear and J.A. Van Orman, Geochim. Cosmochim. Ac., 71 (2007) 1312-1323.
 D. Nevins, F.J. Spera and M.S. Ghiorso, Am. Mineral, 94 (2009) 975-980.
 T. Matsumiya, A. Nogami and Y. Fukuda, ISIJ Int., 33 (1993) 210-217.
 W.G. Seo and F. Tsukihashi, ISIJ Int., 44 (2004) 1817-1825.
 K. Shimoda and K. Saito, ISIJ Int., 47 (2007) 1275-1279.
 T. Asada, Y. Yamada and K. Ito, ISIJ Int., 48 (2008) 120-122.
 W.G. Seo and F. Tsukihashi, Mater. Trans., 46 (2005) 1240-1247.
 D.K. Belashchenko and O.I. Ostrovskii, Inorg. Mater., 38 (2002) 958-963.
 D.K. Belashchenko and O.I. Ostrovskii, Inorg. Mater., 38 (2002) 196-203.
 L.J. Zhou and W.N. Wang, Metall. Mater. Trans. B, 47 (2016) 1548-1552.
 Y. Sasaki and M. Iguchi, ISIJ Int., 49 (2009) 602-604.
 K.C. Mills, L. Yuan, Z. Li, G.H. Zhang and K.C. Chou, High Temp. Mater. Proc, 31 (2012) 301-321.
 Y. Sun, K. Zheng, J.L Liao, X.D. Wang and Z.T. Zhang, ISIJ Int., 54 (2014) 1491-1497.
 Y.Q. Sun and Z.T. Zhang, Metall. Mater. Trans. B., 46 (2015) 1549-1554.
 Y.Q. Sun, Z.T. Zhang, L.L. Liu and X.D. Wang, J. Non-Cryst. Solids, 420 (2015) 26-33.
 Z.J. Wang, Q.F. Shu, S. Sridhar, M. Zhang, M. Guo and Z.T. Zhang, Metall. Mater. Trans. B., 46 (2015) 758-765.
Jiang Diao (*), Quan Zhang, Yong Qiao, Lu Jiang and Bing Xie
Quan Zhang, Yong Qiao, Lu Jiang, Bing Xie, College of Materials Science and Engineering, Chongqing University, Chongqing 400044, P. R. China
(*) Corresponding author: Jiang Diao, College of Materials Science and Engineering, Chongqing University, Chongqing 400044, P. R. China, E-mail: firstname.lastname@example.org
Received July 19, 2016; accepted December 30, 2016
Table 1: Potential parameters used in computation. i j [A.sub.ij] (gA/f[s.sup.2]) [B.sub.ij]/ (1/A) [C.sub.ij]/ (g[A.sup.8]/ f[s.sup.2]) 0 0 2.40E-20 5.88E + 00 2.78E-25 Si Si 3.47E-23 6.25E + 00 0 Ca Ca 5.27E-21 6.25E + 00 6.95E-26 Fe Fe 4.70E-23 3.45E + 00 0 P P 4.56E-22 7.06E + 00 0 O Si 1.01E-21 6.06E + 00 0 O Ca 1.15E-20 6.06E + 00 1.39E-25 O Fe 6.41E-22 5.16E + 00 0 O P 3.04E-23 3.45E + 00 0 Si Ca 4.28E-22 6.25E + 00 0 Si Fe 9.22E-23 1.29E + 01 0 Si P 1.73E-23 1.25E + 01 4.49E-25 Ca Fe 3.53E-23 6.25E + 00 0 Ca P 2.64E-21 1.25E + 01 0 Fe P 2.05E-22 6.25E + 00 0 Table 2: Chemical compositions and number of atoms. Nos 1~4 slags concern the parameter series of slag basicity. No. 2 and Nos 5~8 slags concern the parameter series of [P.sub.2][O.sub.5] content. No. 2 and Nos 9~10 slags concern the parameter series of FeO content. Mass fraction/wt% No. CaO Si[O.sub.2] [P.sub.2][O.sub.5] 1 36 24 10 2 32.8 27.2 10 3 28.4 31.6 10 4 22.5 37.5 10 5 35.5 29.5 5 6 33.8 28.2 8 7 31.6 26.4 12 8 30.5 25.5 14 9 30 25 10 10 27.3 22.7 10 Atomic number No. FeO R Ca Si P Fe O 1 30 1.5 675 420 148 437 2,320 2 30 1.2 606 471 146 432 2,345 3 30 0.9 517 537 144 425 2,377 4 30 0.6 401 624 141 416 2,418 5 30 1.2 663 515 74 435 2,313 6 30 1.2 628 489 117 434 2,332 7 30 1.2 583 454 175 430 2,358 8 30 1.2 561 437 203 429 2,370 9 35 1.2 565 440 149 514 2,332 10 40 1.2 524 408 152 598 2,318 Table 3: Diffusivity sequence of ions in different molten systems. References System Matsumiya et al.  CaO-Si[O.sub.2] Seo and Tsukihashi  CaO-Si[O.sub.2] Shimoda and Saito  CaO-MgO-[Al.sub.2][O.sub.3]-Si[O.sub.2] Asada et al.  CaO-Ca[F.sub.2]-Si[O.sub.2] Seo and Tsukihashi  FeO-Si[O.sub.2] Belashchenko and CaO-FeO-Si[O.sub.2] Ostrovskii  Belashchenko and CaO-[P.sub.2][O.sub.5] Ostrovskii  Present work CaO-Si[O.sub.2]-[P.sub.2][O.sub.5]-FeO References T, K Sequence Matsumiya et al.  1,873 Ca>O>Si Seo and Tsukihashi  1,873 Ca>O>Si Shimoda and Saito  1,873 Ca>Mg>Al>O>Si Asada et al.  1,773 F>Ca>O>Si Seo and Tsukihashi  1,600 Fe>O>Si Belashchenko and Ostrovskii  1,873 Ca>Fe>O>Si Belashchenko and Ostrovskii  2,000 Ca>O>P Present work 1,673 Ca>Fe>O>P>Si Table 4: Self-diffusion coefficients of different ions /([A.sup.2]/ps). No. Ca Si Fe P O 1 0.065 0.027 0.057 0.038 0.053 2 0.056 0.023 0.053 0.032 0.045 3 0.045 0.019 0.042 0.026 0.037 4 0.04 0.017 0.037 0.023 0.033 5 0.097 0.04 0.092 0.057 0.08 6 0.083 0.034 0.075 0.049 0.068 7 0.05 0.021 0.049 0.028 0.039 8 0.043 0.018 0.042 0.025 0.035 9 0.069 0.028 0.061 0.04 0.056 10 0.097 0.04 0.092 0.057 0.08 Table 5: Viscosities calculated by self-diffusion coefficients /(Pa-s). No. 1 1 3 4 5 6 7 8 9 10 [eta] 0.15 0.15 0.2 0.26 0.14 0.15 0.19 0.28 0.14 0.11
|Printer friendly Cite/link Email Feedback|
|Author:||Diao, Jiang; Zhang, Quan; Qiao, Yong; Jiang, Lu; Xie, Bing|
|Date:||Feb 1, 2018|
|Previous Article:||Aging of a Low Carbon Heat Resistant Cast Alloy.|
|Next Article:||Melting Features and Viscosity of Ti[O.sub.2]-Containing Primary Slag in a Blast Furnace.|