Show Summary Details
More options …

# High Temperature Materials and Processes

Editor-in-Chief: Fukuyama, Hiroyuki

Editorial Board: Waseda, Yoshio / Fecht, Hans-Jörg / Reddy, Ramana G. / Manna, Indranil / Nakajima, Hideo / Nakamura, Takashi / Okabe, Toru / Ostrovski, Oleg / Pericleous, Koulis / Seetharaman, Seshadri / Straumal, Boris / Suzuki, Shigeru / Tanaka, Toshihiro / Terzieff, Peter / Uda, Satoshi / Urban, Knut / Baron, Michel / Besterci, Michael / Byakova, Alexandra V. / Gao, Wei / Glaeser, Andreas / Gzesik, Z. / Hosson, Jeff / Masanori, Iwase / Jacob, Kallarackel Thomas / Kipouros, Georges / Kuznezov, Fedor

IMPACT FACTOR 2018: 0.427
5-year IMPACT FACTOR: 0.471

CiteScore 2018: 0.58

SCImago Journal Rank (SJR) 2018: 0.231
Source Normalized Impact per Paper (SNIP) 2018: 0.377

Open Access
Online
ISSN
2191-0324
See all formats and pricing
More options …
Volume 34, Issue 8

# Kinetics Study on Reduction of CaWO4 by Si from 1423 K to 1523 K

Qifeng Shu
• Corresponding author
• School of Metallurgical and Ecological Engineering, University of Science and Technology Beijing, Beijing 100083, Beijing, China
• State Key Laboratory of Advanced Metallurgy, University of Science and Technology Beijing, Beijing 100083, Beijing, China
• Email
• Other articles by this author:
/ Jing Wu
• School of Metallurgical and Ecological Engineering, University of Science and Technology Beijing, Beijing 100083, Beijing, China
• State Key Laboratory of Advanced Metallurgy, University of Science and Technology Beijing, Beijing 100083, Beijing, China
• Email
• Other articles by this author:
/ Kuochih Chou
• School of Metallurgical and Ecological Engineering, University of Science and Technology Beijing, Beijing 100083, Beijing, China
• State Key Laboratory of Advanced Metallurgy, University of Science and Technology Beijing, Beijing 100083, Beijing, China
• Email
• Other articles by this author:
Published Online: 2015-01-10 | DOI: https://doi.org/10.1515/htmp-2014-0161

## Abstract

Investigation on reduction kinetics of scheelite(CaWO4) provides important fundamental knowledge to control and optimize process of Ferrotungsten manufacturing and direct reduction of scheelite in steelmaking. In this work, the isothermal reduction kinetics of CaWO4 powder by Si powder at 1,423 K, 1,473 K and 1,523 K, where all reactants and products are in solid state, was investigated by using X-ray diffraction quantitative analysis. Scanning electronic microscopy (SEM) attached with energy dispersive spectra (EDS) was also employed to investigate the microstructure of reaction interface. Both Jander model and Ginstling–Brounshtein model could describe kinetics of reduction of CaWO4 well, whereas 3D interfacial reaction model could not describe the kinetics of reduction of CaWO4. The rate determining step for reduction of CaWO4 could be diffusion of Si across product layer. The values of activation energies obtained by fitting kinetic data using Jander and Ginstling–Brounshtein model are as great as 379.93 kJ/mol and 387.16 kJ/mol respectively. The oxygen partial pressure during reduction has impact on the kinetics of reduction. The reduction was retarded by increase of oxygen partial pressure.

PACS.: 82.20.Pm

## Introduction

Tungsten steel is one of the most important high speed steels. Addition of tungsten in steel will bring great hardness and good stability at high temperature. Usually, ferro-tungsten alloy was added into steel during refining process of steel. The ferro-tungsten alloy is manufactured by reduction of scheelite (CaWO4) or Wolframite with Si-Fe alloy or aluminum alloy [1]. Since high melting point of tungsten, the manufacturing of ferro-tungsten requires high temperature and leads to large energy consumption, serious environmental pollution and other problems.

There have been many attempts to direct use of oxides instead of ferroalloys as alloying agents in steelmaking [26]. Chychko [24] performed plant trial (3 ton scale) using three different alloying mixtures (MoO3+C; MoO3+C+FeOx; and MoO3+C+CaO) as Mo-alloying agents in steelmaking. Mo yield up to 93% had been achieved in plant trial. Alternative method using scheelite directly for productions of tungsten steel has also been proposed [7, 8]. According to this method, scheelite ore or scheelite ore mixed with reductants (Si-Fe alloy, Aluminum alloy) was directly added into molten steel, and reduced by Si, Al in steel or reductants. Successful plant trail had been carried out in a 20 ton Electrical Arc Furnace to produce high W die steel using scheelite [8]. The overall yield of W could be up to 92%.

Investigation on reduction kinetic of scheelite provides important fundamental knowledge to control and optimization of process of ferrotungsten manufacture or direct reduction of scheelite in steelmaking. However, regarding reduction of scheelite, only carbothermic reaction was studied aiming at producing WC from scheelite [9, 10]. There are very few reports on reduction kinetics of CaWO4 by Si. In this work, the isothermal reduction kinetics of CaWO4 by Si from 1,423 K to 1,523 K, where all reactants and products are in solid state, was investigated by using X-ray diffraction quantitative analysis. Scanning electronic microscopy (SEM) attached with energy dispersive spectral (EDS) was employed to investigate the microstructure of reaction interface.

## Experimental

CaWO4 for reduction experiments was synthesized through solid-state reaction between CaCO3 and WO3. CaCO3 (purity>99%) and WO3(purity>99%) powder was taken as raw materials. Raw powders were weighed and then mixed thoroughly in an agate mortar. Powder mixtures were pressed into a pellet with diameter of 8 mm and then placed in a corundum crucible. The pellets were heated at 1,343 K for 24 h, then crushed and ground into powder. As shown in Figure 1, XRD analysis on synthesized powder has verified that CaWO4 with high purity has been synthesized.

Figure 1

XRD pattern of synthesized CaWO4

Analytical grade Si powder (200 seize) was adopted as reductant. Synthesized CaWO4 and Si powder with molar ratio of 3:2 were mixed in an agate mortar with ethanol as mixing media. The mixed powder was pressed into pellet with mass of about 0.55 g and diameter of 8 mm under pressure of 2.6 MPa.

The reduction of CaWO4 by silicon was carried out in a vertical corundum tube furnace with silicon carbide as heating elements. An S-type (Pt-Pt/10 Rh) thermocouple was mounted underneath the even temperature zone of furnace to control and measure temperature inside furnace. Argon with high purity (>99.999%) was introduced from bottom of corundum tube and flow rate of argon was maintained at 200 ml/min (STP) using a rotameter. After temperature in the furnace reached right temperature, pellets were put in a graphite crucible (ID: 9 mm; OD: 12 mm; H: 35 mm), and then crucible was placed in the even temperature zone of furnace. The crucible were quickly pulled out from furnace and quenched in air when heating was completed. To investigate the effect of different atmosphere, reduction experiments were also carried out in air or in argon with corundum crucible.

The samples were ground into fine powder prior to XRD measurements. XRD examinations were performed using a 18 KW X-ray diffractometer (RIGAKU TTR III) in air with Cu Ka radiation (λ=1.5418 Å, target voltage: 40 kV, target current: 300 mA). Patterns were collected by scanning from 2θ=10° to 90° with step size of 0.01° and a scanning rate of 10°/min. Before XRD quantitative analysis, the background and Kα2 radiation contribution of acquired XRD patterns were removed by employing Crystallographica Search-March (Oxford Cryosystems) software.

Before SEM examination, samples were embedded with resin, ground, polished and coated with gold. SEM examinations were performed using SEM (FEI MLA 250). The working voltage was 25 kV.

## Isothermal reduction kinetics of CaWO4 by Si from 1,423 K to 1,523 K

The reduction of CaWO4 by Si could be expressed as following equation: $3\mathrm{S}\mathrm{i}+2\mathrm{C}\mathrm{a}\mathrm{W}{\mathrm{O}}_{4}=2\mathrm{C}\mathrm{a}\mathrm{O}+2\mathrm{W}+3\mathrm{S}\mathrm{i}{\mathrm{O}}_{2}$(1)Figures 24 showed XRD patterns for samples treated for various periods at 1,423 K, 1,473 K and 1,523 K respectively. It could be seen that diffraction peaks of CaWO4, Si and W could be distinguished in the patterns. However, there is not any product phase except W which could be distinguished in term of XRD patterns. This could be due to the formation of some non-crystalline phases. Nevertheless, the fraction conversion could be determined by only employing diffraction peaks of CaWO4 and W.

Figure 2

XRD patterns for samples heated at 1,523 K for various periods

Figure 3

XRD patterns for samples heated at 1,473 K for various periods

Figure 4

XRD patterns for samples heated at 1,423 K for various periods

The fractional conversions were determined using an X-ray quantitative analysis method based on RIR (Relative Intensity Ratio) values [11]. The method of X-ray quantitative analysis adopted in the present work was described briefly as follows.

The ratio of mass percentage of CaWO4 to W could be calculated from intensity of the strongest peaks for CaWO4 (peak(1,1,2)) and W(peak (1,1,0)) according to following equation. $\frac{\mathrm{\omega }\left(\mathrm{C}\mathrm{a}\mathrm{W}{\mathrm{O}}_{4}\right)}{\mathrm{\omega }\left(\mathrm{W}\right)}=\frac{\mathrm{R}\mathrm{I}\mathrm{R}\left(\mathrm{W}\right)}{\mathrm{R}\mathrm{I}\mathrm{R}\left(\mathrm{C}\mathrm{a}\mathrm{W}{\mathrm{O}}_{4}\right)}×\frac{\mathrm{I}\left(\mathrm{C}\mathrm{a}\mathrm{W}{\mathrm{O}}_{4}\right)}{\mathrm{I}\left(\mathrm{W}\right)}$(2)where ω(CaWO4) and ω(W) are mass percentages of CaWO4 and W respectively. RIR(W) and RIR(CaWO4) are relative intensity ratio values for W and CaWO4 respectively. I(CaWO4) and I(W) are intensities of the strongest peaks for CaWO4 (peak(1,1,2)) and W(peak (1,1,0)) respectively.

The fractional conversion could be further calculated from ω(CaWO4)/ω(W) according to the following equations. $\alpha =\frac{1}{1+\frac{\omega \left({\text{CaWO}}_{\text{4}}\right)}{\omega \left(\text{W}\right)}×\frac{183.84}{287.84}}$(3)

RIR(W)/RIR(CaWO4) could be calculated from relative intensity ratios of W and CaWO4 obtained from ICDD cards. Alternatively, it could be also determined from XRD measurements on mixture of W and CaWO4 with ω(CaWO4)/ω(W) = 1. Figure 5 showed XRD patterns for CaWO4+W mixture (ω(CaWO4)/ω(W) = 1). The XRD measurements were repeated twice. Mean value for RIR(W)/RIR(CaWO4) is 1.467.

Figure 5

XRD patterns for CaWO4+W mixture (ω(CaWO4)/ω(W)=1)

The fractional conversions for various periods at 1,423 K, 1,473 K and 1,523 K were determined from XRD patterns in Figures 24 and shown in Figure 6. According to phase equilibrium data for the present system [12], all reactants and products are in solid state at 1,423–1,523 K. Reaction between solid Si and solid CaWO4 particles takes places during reduction process. After contact between Si and CaWO4, a product layers could be generated around reactant. With progress of reaction, this layer could be gradually accumulated. This product layer isolates the reactant from each other, therefore hinders the overall reaction. Continuous reaction could only persist through (1) diffusion of reactant (in the form of atoms or ions) through the product layer to the reaction interface; (2) chemical reactions among reactants at the reaction interface. In most cases, the diffusion of reactant through the product layer is slower than interfacial chemical reaction and often becomes the rate-controlling step of the overall reaction [13, 14]. Jander model [15] and Ginstling–Brounshtein model [16] are most frequently used to describe kinetics of three-dimensional diffusion controlling reaction. The expressions of Jander model and Ginstling–Brounshtein model are as follows: $\text{Jander}:{\left(1-{\left(1-\alpha \right)}^{1/3}\right)}^{2}=kt$(4) $\text{Ginstling}-\text{Brounshtein}:\text{\hspace{0.17em}}1-\frac{2}{3\alpha }-{\left(1-\alpha \right)}^{2/3}=kt$(5)

Figure 6

Fractional conversions as functions of reaction time at various temperatures

where α is fractional conversions. t is time, k is rate constant and follows Arrhenius equation $k=A\mathrm{exp}\left(-E/RT\right)$ where E and A designate the apparent activation energy and pre-exponential factor respectively.

In some cases, the interfacial chemical reaction is slower than diffusion through the product layer and should be considered as rate-controlling step. The equation for three-dimensional chemical reaction could be as follow. $1-\left(1-\mathrm{\alpha }{\right)}^{1/3}=kt$(6)Experimental kinetics data were fitted using Jander, Ginstling–Brounshtein and 3D chemical reaction model. Fitting results are shown in Figures 79. It could be seen from Figures 7 and 8 that both Jander and Ginstling–Brounshtein model fit the kinetic data with fairly accuracy. The correlation coefficients attained by the least square fitting are as large as 0.9945 and 0.9925 for Jander and Ginstling–Brounshtein model respectively. The good agreements between fitted and experimental kinetic data using Jander and Ginstling–Brounstein indicate that 3D diffusion could be rate-determining step of overall reaction. In comparison, as shown in Figure 9, there are large differences between fitted kinetic data using 3D interfacial reaction model and experimental kinetic data, and the correlation coefficient is only. The bad fitting with 3D interfacial reaction model indicates that interfacial chemical reaction could not be rate-determining step of overall reaction. Accordingly, it could be inferred from above fittings that the reduction of CaWO4 by Si from 1,423 K to 1,523 K could be controlled by three-dimensional diffusion through the product layer.

Figure 7

Result of fitting kinetic data using Ginstling–Brounshtein model

Figure 8

Result of fitting kinetic data using Jander model

Figure 9

Result of fitting kinetic data using 3D interfacial chemical reaction model

The slopes of straight line in Figures 7 and 8 directly provided rate constants k at different temperatures. As shown in Figures 10 and 11, the natural logarithms of rate constants k were plotted versus reciprocal of absolute temperature to yield the apparent activation energy and pre-exponential values. The temperature dependences of rate constant obtained from the two fitting lines are expressed as follows: $\begin{array}{rl}k=\phantom{\rule{thickmathspace}{0ex}}& 3.80×{10}^{7}\left({\mathrm{s}}^{-1}\right)exp\left(\frac{-379.93\left(\mathrm{k}\mathrm{J}/\mathrm{m}\mathrm{o}\mathrm{l}\right)}{RT}\right)\\ & \text{\hspace{0.17em}}\left(\mathrm{J}\mathrm{a}\mathrm{n}\mathrm{d}\mathrm{e}\mathrm{r}\text{\hspace{0.17em}}\mathrm{m}\mathrm{o}\mathrm{d}\mathrm{e}\mathrm{l}\right)\end{array}$(7) $\begin{array}{rl}k=\phantom{\rule{thickmathspace}{0ex}}& 2.72×{10}^{2}\left({\mathrm{s}}^{-1}\right)exp\left(\frac{-387.16\left(\mathrm{k}\mathrm{J}/\mathrm{m}\mathrm{o}\mathrm{l}\right)}{RT}\right)\\ & \text{\hspace{0.17em}}\left(\mathrm{G}\mathrm{i}\mathrm{n}\mathrm{s}\mathrm{t}\mathrm{l}\mathrm{i}\mathrm{n}\mathrm{g}\phantom{\rule{negativethinmathspace}{0ex}}\phantom{\rule{negativethinmathspace}{0ex}}-\phantom{\rule{negativethinmathspace}{0ex}}\phantom{\rule{negativethinmathspace}{0ex}}\mathrm{B}\mathrm{r}\mathrm{o}\mathrm{u}\mathrm{n}\mathrm{s}\mathrm{h}\mathrm{t}\mathrm{e}\mathrm{i}\mathrm{n}\text{\hspace{0.17em}}\mathrm{m}\mathrm{o}\mathrm{d}\mathrm{e}\mathrm{l}\right)\end{array}$(8)Only few reports could be found in literature on kinetics of silicothermic reduction. These studies mainly focused on kinetics of silicothermic reduction of calcined dolomite with regard to Pidgeon process [17, 18]. These studies all suggested that the reduction of dolomite by silicon was controlled by solid-state diffusion of reactants, which is similar to the reaction mechanism revealed in the present work. Wulandari et al. [17] reported an activation energy value of 299~322 kJ/mol for silicothermic reduction of calcined dolomite in flowing argon, which is close to our activation energy value.

Figure 10

Natural logarithm of rate constant k obtained by fitting Ginstling–Brounshtein model as a function of reciprocal of temperature

Figure 11

Natural logarithm of rate constant k obtained by fitting Jander model as a function of reciprocal of temperature

## Microstructure of samples

To further elucidate the mechanism for reaction, microstructures for samples with various reaction time and temperature were investigated using SEM-EDS. Figure 12 showed some examples of SEM micrographs for samples reduced for various periods at various temperatures. The light grey phase and dark grey phase were distinguished as CaWO4 and Si phase respectively. The white small particles found between light grey and dark grey phases were determined as W. The composition of black phase was determined as mixture of CaO, SiO2 and WO3, which could be formed by reaction between product (CaO, SiO2) and CaWO4. This phase (denoted by CaxSiyWzOt, x, y, z, t could be variable) was not found in XRD analysis (see Figures 24), which may be due to its amorphous nature.

It could be seen from the SEM micrographs that produced tungsten particle has a mean particle size less than 1 µm. The nature of tungsten particle determined that tungsten particle should not be seen as an obstacle for diffusion of reactant to interface. Si particle was separated from CaWO4 by CaxSiyWzOt phase. The continuous reduction requires transfer of Si or CaWO4 across CaxSiyWzOt phase to contact the other reactant phase. As seen in Figure 12(a-d), CaWO4 particle was surrounded by small tungsten particles. It is postulated that Si diffuses across CaxSiyWzOt phase to react with CaWO4 at interface between CaWO4 and CaxSiyWzOt phase. The main resistance for overall reduction could come from the diffusion of Si in CaxSiyWzOt phase.

Figure 12

Microstructures of samples heated at different temperature and for various periods (a): 1,423 K, 60 min (b) 1,473 K, 10 min (c) 1,473 K, 90 min (d) 1,523 K, 30 min

## Effects of different atmosphere on reduction of CaWO4 by Si

The kinetics of CaWO4 reduction by Si was investigated in argon with graphite crucible. To investigate the effect of different atmosphere on reduction, reduction experiments were also conducted in two different atmospheres: in air and in argon with corundum crucible. Since reaction between graphite and trace oxygen in argon would deplete the oxygen in argon, the oxygen pressure around graphite crucible is relatively lower than that around corundum crucible. Therefore, oxygen pressure in graphite crucible is lowest. The XRD patterns for samples in different atmosphere with reaction time of 90 min at 1,473 K were shown in Figure 13. It could be seen from the figure that the amount of produced W increases with decrease of oxygen pressure. The retard of CaWO4 reduction at high oxygen pressure could be attributed to the surface oxidation of Si. Si powder tends to be oxidized at high oxygen pressure and a thin oxidation layer could form on the surface of powder. The oxides “coating” is an obstacle for diffusion of Si to reaction interface, therefore overall reaction slowed down. It is proposed that weak reducing atmosphere would be beneficial to the reduction of Si.

Figure 13

XRD pattern for samples heated at 1,473 K for 90 min in (1) air with a corundum crucible; (2) argon with a corundum crucible; (3) argon with a graphite crucible

## Conclusions

The kinetics of reduction of CaWO4 by Si was investigated by using quantitative XRD analyses at 1,423 K, 1,473 K and 1,523 K. Following conclusions could be reached:

• 1.

Both Jander model and Ginstling–Brounshtein model could describe well the kinetics of reduction of CaWO4, whereas 3D interfacial reaction model could not describe the kinetics of reduction of CaWO4.

• 2.

The rate-determining step for reduction of CaWO4 could be diffusion across product layer.

• 3.

The values of activation energies obtained by fitting kinetic data using Jander and Ginstling–Brounshtein model are as great as 379.93 kJ/mol and 387.16 kJ/mol, respectively

• 4.

The oxygen partial pressure during reduction has impact on the kinetics of reduction. The reduction was retarded by increase of oxygen partial pressure.

## Acknowledgment

Financial supports from NSFC (grants no. 51174018), Beijing Higher Education Young Elite Teacher Project (No. YETP0346) and State key laboratory of advanced metallurgy are gratefully acknowledged.

## References

• 1.

Yücel O, Özçelebi MA. Reduction smelting of bursa-uludağ tungsten concentrates by the aluminothermic process. Scand J Metall 2000;29:108–13.

• 2.

Chychko A, Teng L, Nzotta M, Seetharaman S. Fe2MoO4 as a precursor material for Mo alloying in steel (Part II): up-scaling test. Steel Res Int 2011;82:886–97.

• 3.

Chychko A, Teng L, Nzotta M, Seetharaman S. Synthesis and characterization of Fe2MoO4 as a precursor material for Mo alloying in steel. Steel Res Int 2011;82:269–76.

• 4.

Chychko A, Teng L, Seetharaman S. MoO3 evaporation Studies from binary systems towards choice of Mo precursors in EAF. Steel Res Int 2010;81:784–91.

• 5.

Song SQ, Xue ZL, Zhang J, Dong J. Silicothermic self-reducing MoO3 briquettes for direct alloying of Mo in molten steel. Ironmaking Steelmaking 2014;41:628–32.

• 6.

Song SQ, Xue ZL, Zhang B, Wang WX, Liu RN, Wang GL. Smelting vanadium steel by direct alloying with self-reduction briquette of V2O5. Ironmaking Steelmaking 2014;41:7–11.

• 7.

Chen ZX, Li JR. Investigation of scheelite instead of ferrotungsten in steelmaking. Iron Steel 1992;27:15–18 (in Chinese). Google Scholar

• 8.

Li ZB, Guo PM, Zhang HS. Thermodynamic analysis of smelting alloy steels with scheelite, molybdenum oxide and vanadium slag. J Iron Steel Res 1999;11:15 (in Chinese). Google Scholar

• 9.

Johnston RF, Nguyen HT. Direct reduction carburisation of scheelite with carbon. Miner Eng 1996;9:765–73.

• 10.

Singh H, Pandey OP. Direct synthesis of nanocrystalline tungsten carbide from scheelite ore by solid state reaction method. Ceram Int 2013;39:785–90.

• 11.

Chung BF. Quantitative interpretation of X-ray diffraction patterns of mixtures. I. matrix-flushing method for quantitative multicomponent analysis. J Appl Crystallogr 1974;7:519.

• 12.

Chang LLY, Scroger MG, Phillips B. Alkaline-earth tungstates: equilibrium and stability in the M-W-O systems. J Am Ceram Soc 1966;49:385–90.

• 13.

Yan BJ, Liu JH, Zhang JY. Kinetic studies of LaCoO3-δ preparation from solid state reactions. High Temp Mater Processes 2004;23:163–76. Google Scholar

• 14.

Shu QF, Liu JH, Yan BJ, Zhang JY. Phase formation mechanism and kinetics in solid-state synthesis of undoped and calcium-doped lanthanum manganite. Mater Res Bull 2009;44:649–53.

• 15.

Jander WZ. Reaktionen im festen Zustande bei höheren Temperaturen. Reaktionsgeschwindigkeiten endotherm verlaufender umsetzungen. Z Anorg Chem 1927;163:1.

• 16.

Ginstling M, Brounshtein BI. Concerning the diffusion kinetics of reactions in spherical particles. J Appl Chem USSR 1950;23:1327. Google Scholar

• 17.

Wulandari W, Brooks GA, Rhamdhani MA, Monaghan BJ. Kinetic analysis of silicothermic process under flowing argon atmosphere. Can Metall Q 2014;53:17–25.

• 18.

Morsi LM, Ali HH. Kinetics and mechanism of silicothermic reduction process of calcined dolomite in magnetherm reactor. Int J Min Proc 2014;127:37–43.

Accepted: 2014-11-23

Published Online: 2015-01-10

Published in Print: 2015-12-01

Citation Information: High Temperature Materials and Processes, Volume 34, Issue 8, Pages 805–811, ISSN (Online) 2191-0324, ISSN (Print) 0334-6455,

Export Citation