2024 Volume 64 Issue 15 Pages 2238-2244
We measured the density of eight types of molten alkali silicate slag within the range of 1673–1823 K using the Archimedean double-bob method. The melt density of binary alkaline silicate glass linearly increased with a decrease in temperature. The temperature coefficient of the density ranged from −21.5 to −13.2 × 10−5 g·cm−3·K−1, and its order corresponded to the order of the ionic radii of the alkali oxide components. We applied Doolittle’s free volume theory by combining the reported values of melt viscosity with the measured densities and volume expansion coefficients of melts and glass solids to propose an equation for predicting melt density from slag composition and temperature. Melt densities in the range of 1273–1823 K were predicted and compared with experimental values.
In metal and steel production, molten slag plays a role in protecting metal melts from the atmosphere and extracting unwanted impurities.1,2,3) To achieve carbon neutrality by 2050, the development and introduction of high-temperature melting technologies, fuel conversion, and the use of slag by-products4) are being considered. Owing to the aforementioned advanced melting processes, it is important to understand the physical properties of oxide slag melts of various compositions.5,6) This study focuses on the density of high-temperature melts.7) This density is directly related to the weight and volume of molten slag. The coefficient of volume expansion can be calculated from the temperature dependence of the density. Highly reliable melt densities over a wide temperature range are required during the simulation of the melting process of multicomponent silicate glass because it is similar to slag.8) The temperature dependence of melt density is also related to the homogenization of multi-component melts due to thermal convection current and to the phase separation process of slag/metal and slag/matte melts. Additionally, the density is closely concerned with the structure of the melt at the atomic level in microscopic terms. Thus, the density of high-temperature molten slag is an important physical property.
Methods for measuring the density of molten slag at high temperatures mainly include the Archimedean double-bob,9,10) maximum bubble pressure,11) and static sessile drop methods.12) The Archimedean double-bob method is a highly accurate measurement method.9,10,13,14) However, the apparatus used is typically made by researchers. Therefore, the shapes and characteristics of the furnace, heating element, and temperature control device differ. Thus, the error and reproducibility in each measurement must be evaluated. The density of oxide slag melt at high temperatures is highly dependent on the chemical composition and temperature.9,10,13,14) The composition of multicomponent silicate slag melts used in industrial applications varies, and a wide temperature range is typically used during the production process. Thus, a highly accurate prediction equation for the density of the melt from a chemical composition over a wide temperature range is desired. Conventional prediction formulas have been proposed based on reported values and partial molar fractions of slag components.15,16) However, high-accuracy prediction equations have not necessarily been introduced. Highly accurate measured values are required in the development of density prediction equations. In glass solids formed by the rapid cooling of molten slag, lattice oscillations of ion pairs cause thermal expansion.17) Furthermore, for molten slag melt, the contribution of free volume is considered.
This study attempted to create a prediction equation for the density of molten slag based on Doolittle’s free volume theory.18,19) According to Doolittle, the temperature dependence of the viscosity of organomolecules and liquids is associated with free volume. The temperature dependence of density can be converted into the coefficient of volume expansion (temperature dependence of volume). Thus, it is assumed that the density of molten slag can be predicted from the temperature dependence of viscosity, according to Doolittle’s free volume theory. In the free volume theory, the total volume (v) of molten slag is the sum of the occupied volume (v0) and free volume (vf). The occupied volume is the sum of the volume of the molecule and the volume related to vibration; the rest is the free volume. Furthermore, the occupied volume is dependent on temperature. However, the temperature dependence of the viscosity of the melt is governed by changes in the free volume since the temperature dependence of the free volume is considerably greater. Based on these ideas, Doolittle’s viscosity formula is expressed by Eq. (1).
(1) |
where constants A and B are determined by the system. It is known that the viscosities of organic small molecules and polymer liquids are well expressed by Eq. (1). Here, the occupied volume was calculated from the density and thermal expansion coefficient of the glass solid and glass transition temperature. First, we improved the handmade setup of the density measurement of high-temperature melts through the Archimedean double-bob method.14) Thereafter, the measurement conditions of a standard slag sample with relatively high accuracy were established. The temperature dependence of the density of binary alkali silicate glass solids and melts was investigated over a wide temperature range. Alkali silicate glasses are well known for their microscopic structures20,21) in the basic compositional system. Finally, using the measured density, volume expansion coefficient, and glass transition temperature for slag solids and melts, we developed an equation for predicting melt density based on Doolittle’s free volume theory, and the measured and predicted values were compared.
The density of binary alkaline silicate slag melt with a viscosity ranging from 10−2 to 10 Pa·s was measured using the Archimedean double-bob method. Figure 1 shows that a Pt weight of volume V connected to a wire with radius r was immersed in a melt with density ρ. The difference between the present and a previous study14) is that the radiant heat from the electric furnace to the balance was blocked as much as possible with refractory materials and ice water for cooling. The buoyancy (W) received by the weight is expressed by Eq. (2), including the surface tension term on the wire, effect of air convection, and vapor deposition of the substance.
(2) |
where W (10−3 kg) is the buoyancy force that the weight receives from the liquid (melt), and V (10−6 m3) is the volume of the weight in the liquid (melt) (103 kg·m−3). ρ is the density of the melt, r (10−3 m) is the diameter of the wire, and γ (10−3 N·m−1) is the surface tension of the liquid (melt), θ (°) is the contact angle between the hanging thread and liquid surface, and A (10−3 kg) is a constant related to the convection of air and vapor matters on the hanging thread. In Eq. (2), two weights of different volumes are used to offset the effects of the surface tension, vapor deposition, and updraft of the melt on the wire.
(3) |
(4) |
where L and S mean large and small, respectively, and WL is the buoyancy force on the large bob. WS is the buoyancy force on the small bob, VL is the volume of the large bob, and VS is the volume of the small bob. The difference between WL and WS when two bobs were immersed in the melt under the same conditions was obtained. The density of the liquid was obtained from the difference in the buoyancy and the volume difference of the weight using Eq. (5).
(5) |
Considering that the volume of the Pt bob changes because of thermal expansion at high temperatures, the correction of the volume after the thermal expansion of the Pt bob is expressed by Eq. (6).
(6) |
where α is the coefficient of the linear thermal expansion of Pt (1.0102 × 10−5 K−1),22) and ΔT is the temperature difference between the room and density measurements. When determining the density of the melt, the volume difference (VL − VS) of the Pt bob must be determined beforehand. Using the weight of both Pt bobs and the density of Pt,23) the volume difference (ΔV) of the Pt bobs was obtained using Eq. (7).
(7) |
where mL is the weight of a large Pt bob, mS is the weight of a small Pt bob, and ρPt is the density of Pt (21.45 g·cm−3). As shown in Fig. 1, the buoyancy force on the Pt bob was read down to 0.1 mg using an automatic precision balance (CP224S, Sartorius). A wire made of 0.25-mm Pt–Rh13 wt% was used to connect the precision balance and Pt weight. The container used for the molten slag was a Pt–20%Rh crucible, and the weight for immersion was made of Pt. A MoSi2 heating element was used in the setup. The maximum measured temperature with the heat resistance of the aforementioned wire-suspended Pt weight was 1823 K.
To prepare a synthetic slag sample for the measurements, SiO2 (High Purity Chemical Research Institute, 99.9% purity), Li2CO3 (High Purity Chemical Research Institute, 99.9%), Na2CO3 (Kanto Chemical Co., Ltd., 99.8%), and K2CO3 powder (Kanto Chemical Co., Ltd., 99.5%) were used. Table 1 presents the sample composition. The sample names were abbreviated by the content and type of alkali oxides. For example, the 30 mol% Li2O–70 mol% SiO2 sample was denoted by “30LS.” The raw materials were weighed and mixed to obtain the molten slag sample (100 g), which was melted for 2 h in the atmosphere at 1873 K using a Pt–20%Rh crucible. The molten sample was rapidly cooled and vitrified by pouring it onto a stainless-steel mold. The crushed glass sample was remelted in the atmosphere at 1873 K for 2 h, and the density was measured during a temperature-lowering process from 1823 to 1673 K.
Sample | Composition (mol%) | |||
---|---|---|---|---|
Li2O | Na2O | K2O | SiO2 | |
30LS | 30 | – | – | 70 |
35LS | 35 | – | – | 65 |
25NS | – | 25 | – | 75 |
30NS | – | 30 | – | 70 |
35NS | – | 35 | – | 65 |
40NS | – | 40 | – | 60 |
25KS | – | – | 25 | 75 |
30KS | – | – | 30 | 70 |
The buoyancy measurements were started after holding 30 min at each temperature for melt homogenization. Buoyancy readings were continuously performed at each temperature for 30 min, and the average value was used to calculate the density as a buoyancy measurement. The composition and temperature range of the slag for density measurement were determined using the thermodynamic equilibrium calculation software & thermodynamic database, FactSage. A composition of 40Na2O–60SiO2 (mol%), which had a low liquidus temperature and can be measured over a wide temperature range, was selected as the standard. For the molten slag, the reproducibility of the measurement and the comparison between experimental and reported9,24) values were evaluated within the range of 1673–1773 K. The measurement temperature of the density of each slag melt (Table 1) was determined from the liquid phase region in the diagram of the R2O–SiO2 (R = Li, Na, K) binary system up to a maximum of 1823 K.
2.2. Solid Density, Linear Thermal Expansion, and Glass Transition TemperatureStarting materials and sample compositions similar to those used for the determination of the melt density were selected. After weighing and mixing the raw materials (30 g) of the glass samples, they were melted for 2 h in the atmosphere at 1873 K using a Pt crucible. The molten sample was poured onto a carbon mold heated to 673 K and annealed for 1 h at its glass transition temperature plus 20 K and slowly cooled to room temperature at a cooling rate of 1 K·min−1. Finally, a glass sample was obtained. The density of the glass solids at room temperature was measured via the Archimedes method using kerosene as an immersion liquid (accuracy = 0.01 g·cm−3).25) The coefficient of linear expansion of glass, αS, was evaluated using a thermomechanical analyzer (TMA8310, Rigaku). Samples (3.5 × 3.5 × 16 mm3) were heated from room temperature to 873 K at a heating rate of 10 K·min−1. The glass transition temperature (Tg) was calculated from the bend points with different inclinations in the positive thermal expansion curve.25) In addition, αS was calculated from the slope of the expansion curve. The coefficient of volume expansion of the glass solid, βS, was calculated as 3αS ≈ βS.26)
Figure 2 shows the buoyancy measurement results as the raw data of the 40Na2O–60SiO2 (mol%) melt. Figure 3 shows the melt densities calculated from the average of each buoyancy value. Here, the difference between the values of the two measurements was 0.2%. Using the apparatus described in our previous report,14) the difference between two measurements of the slag melt was 0.6%. Here, the reproducibility was improved using the modified apparatus due to the arrangements of ice water cooling and of refractory between the furnace and the electrical balance (Fig. 1). The difference between the density and reported values at each measured temperature was up to 0.23%. These comparisons confirmed that the self-made apparatus (Fig. 1) could measure the melt density with high accuracy and reproducibility.
Figure 4 shows variations in the melt density with temperature, along with the temperature coefficient of density,
Temperature (K) | 30LS | 35LS | 25NS | 30NS | 35NS | 40NS | 25KS | 30KS |
---|---|---|---|---|---|---|---|---|
1673 | 2.143 | 2.124 | 2.216 | 2.216 | 2.220 | 2.201 | 2.180 | 2.167 |
1723 | 2.136 | 2.119 | 2.212 | 2.214 | 2.212 | 2.187 | 2.169 | 2.158 |
1773 | 2.127 | 2.109 | 2.203 | 2.199 | 2.203 | 2.179 | 2.164 | 2.144 |
1823 | 2.124 | 2.105 | 2.193 | 2.191 | 2.185 | 2.165 | 2.155 | 2.135 |
Table 3 lists the values of the measured density (ρS), coefficient of linear thermal expansion (αS), calculated coefficient of body expansion (βS), and glass transition temperature (Tg) of the alkaline silicate glass solids at room temperature. Figure 5 shows the measurement results of the density of the alkaline silicate glass solids and melts, and the molar volume (Vm(T)) is shown as a function of temperature. Vm(T) was calculated using Eq. (8).
(8) |
where Mi and xi are the molecular weight and molar fraction of each component, respectively. ρS and ρL are the densities of the alkaline silicate glass solids and melts, respectively; VS(T) is the molar volume of glass solids, and VL(T) is the molar volume of the melt. Table 4 presents the calculated results of the VS(T) and VL(T). The VS(T) was calculated using Eq. (9), assuming that the volume linearly changed from room temperature to the Tg.
(9) |
Sample | ρs (g·cm−3) | Tg (K) | αs (×10−5 K−1) | βs (×10−5 K−1) |
---|---|---|---|---|
30LS | 2.336 | 748 | 0.977 | 2.932 |
35LS | 2.352 | 741 | 1.125 | 3.375 |
25NS | 2.441 | 752 | 1.199 | 3.598 |
30NS | 2.476 | 753 | 1.397 | 4.191 |
35NS | 2.500 | 729 | 1.544 | 4.632 |
40NS | 2.533 | 724 | 1.663 | 4.990 |
25KS | 2.422 | 781 | 1.382 | 4.146 |
30KS | 2.452 | 793 | 1.702 | 5.107 |
Sample | VS(T) (cm3·mol−1) | VL(T) (cm3·mol−1) |
---|---|---|
30LS | 0.000708T + 21.810 | 0.001754T + 21.336 |
35LS | 0.000784T + 21.012 | 0.002001T + 20.468 |
25NS | 0.000987T + 24.745 | 0.002215T + 24.186 |
30NS | 0.001103T + 24.426 | 0.002542T + 23.758 |
35NS | 0.001243T + 24.248 | 0.002663T + 23.619 |
40NS | 0.001268T + 23.969 | 0.003238T + 23.096 |
25KS | 0.001212T + 28.285 | 0.002808T + 27.522 |
30KS | 0.001620T + 28.580 | 0.003381T + 27.701 |
The VL(T) was calculated using Eq. (10), assuming that the volume change from the melt density, ρL, at a certain temperature to the Tg occurred linearly.
(10) |
In each melt, the molar volume increased with an increase in temperature. This corresponded to a decrease in the density with temperature.
3.4. Prediction of Density of Alkaline Silicate Glass Melts by Multiregression Analysis 3.4.1. Application of Doolittle’s Free Volume TheoryAccording to Doolittle’s free volume theory,18,19) the volume, v(T), of glass is expressed by the sum of the occupied volume, v0(T), and free volume, vf(T) (Eq. (11)).
(11) |
The v0(T) was estimated to be the sum of the volume of the molecule and the volume related to the lattice vibration and equal to the molar volume of the glass solid, VS(T), using Eq. (12).
(12) |
The vf(T) was calculated using Eq. (13).
(13) |
Equations (11) and (12) were substituted into Eq. (1), and a prediction equation expressing the relationship between the viscosity and molar volume was formed.
(14) |
By organizing Eq. (14) with the Arrhenius-type equation, the intercept (A) and slope (B) of the straight line were obtained. Figure 6 shows the relationship between ln η(T) and VS(T)/(VL(T)−VS(T)) for 40NS slag melt as an example. The viscosity η(T) was quoted from the database (INTERGLAD ver. 8).30)
The relationship between the values of ln A and B and the glass composition was investigated, and the regression of the composition of the alkaline silicate glass was obtained via the least squares method. Based on the results, the molar volume of the alkaline silicate glass melts, VL(T), was predicted; the melt density (ρL) was calculated using Eq. (10), and the melt density prediction equation was created. For example, the calculation process of density ρcal for 40NS was described in the following. The obtained values of VS(T) = 25.75 and lnA = −5.940, and B = 0.439 at ln η(T) = 0 were substituted into Eq. (14) and then VL(T) was calculated as 24.40. The density ρcal was calculated as 2.220 from Eq. (10) with the substituted values of 24.40 for VL(T) and 60.8424 for
Using the obtained molar volumes, VS(T) and VL(T), and the viscosity η(T) quoted from the database (INTERGLAD ver. 8),30) the glass composition and calculated values of ln A and B were obtained (Table 5). During the regression analysis, (mol% Li2O, hereinafter referred to as %Li2O), (%Na2O), (%K2O), and (%SiO2) were used as explanatory variables. In Eqs. (15) and (16), the compositions of ln A and B used to perform regression analysis on the alkaline silicate slag composition from the values in Table 5 are shown as a function.
(15) |
(16) |
Sample | ln A | B |
---|---|---|
30LS | −4.597 | 0.244 |
35LS | −4.801 | 0.267 |
25NS | −4.431 | 0.286 |
30NS | −4.937 | 0.329 |
35NS | −5.813 | 0.345 |
40NS | −5.940 | 0.439 |
25KS | −2.590 | 0.279 |
30KS | −4.224 | 0.323 |
30LS | −3.860 | 0.264 |
25KS | −2.443 | 0.319 |
In Eqs. (15) and (16), it is known that coefficients A and B of the alkali oxide, which is a modified oxide,31) were larger than those of the network component, SiO2.
3.4.2. Prediction of Melt DensityThe densities calculated using the prediction equation, ρcal., were compared with the measured densities, ρmea. at 1673 K. Figure 7 shows the difference between the calculated and experimental values. In addition to the measured values used in the regression analysis, the measured values reported in the literature9,24,32) were added. The dotted line in Fig. 7 represents the point where the predicted and measured values coincide. The maximum difference between both values at 1673 K was 0.025 g·cm−3, which was approximately a 1.1% difference from the measured value. Considering that the difference between the predicted and measured values was an average of 0.5%, Eqs. (15) and (16) for the melt density predicted from the composition and temperature were reasonably accurate.
Figure 8 shows the density of the alkaline silicate slag melt measured in this study, ρmea., and a dashed line connecting the density measurement of the glass and the density at the Tg. Furthermore, the melt density predicted over the range of 1273–1823 K in the melt state is shown. ρcal. is indicated by the colored solid lines. Within the range of 1273–1823 K in the melt state, the difference between the measured and predicted values in this study was 0.4% on average.
We evaluated the accuracy of measuring the density of molten alkaline silicate slag with a composition of 40 mol%Na2O–60 mol%SiO2 at high temperatures via the Archimedean double-bob method. In addition, we established a measurement method for a handmade setup based on the comparison and reproducibility between the measured and reported values. The measurements of the melt density of eight binary alkaline silicate glasses linearly increased with a decrease in temperature. The temperature coefficient of density ranged from −21.5 to −13.2 × 10−5 g·cm−3·K−1, and the order of the coefficient corresponded to the order of the ionic radii of the alkali oxide components. The reported values of melt viscosity were combined with the measurement results of the melt density, glass solids, and volume expansion coefficients. Thus, Doolittle’s free volume theory was applied to create a prediction equation for melt density with a good correlation between the predicted and measured values. Although the composition system was limited, melt densities at any temperature within the range of 1273–1823 K were predicted with an error range of 0.4% based on the composition, temperature, molar volume, and viscosity of the alkaline silicate slags.
Authors are requested to declare any conflicts of interest related to the conduct of this research.
No conflicts exist for this research.