## Services on Demand

## Article

## Indicators

## Related links

- Cited by Google
- Similars in Google

## Bookmark

## Journal of the Southern African Institute of Mining and Metallurgy

*Print version* ISSN
0038-223X

### J. S. Afr. Inst. Min. Metall. vol.111 no.10 Johannesburg Oct. 2011

**TRANSACTION PAPER**

**Estimating the physical properties of slags**

** K.C. Mills ^{I}; L. Yuan^{I}; R.T. Jones^{II}**

^{I}Department of Materials, Imperial College, London

^{II}Pyrometallurgy Division, Mintek, Randburg, South Africa

**SYNOPSIS**

The objective of this work was to provide process engineers with values of the physical properties of various slag systems involved in high-temperature processes. Software that calculates the thermophysical properties of slags from chemical composition is available on the www.pyrometallurgy.co.za website. This paper outlines the principles underlying the various models available in the program. The software calculates the following properties of crystalline, glassy, and liquid slags (where appropriate) as a function of temperature: heat capacity, enthalpy, density, viscosity, thermal conductivity, electrical conductivity, and surface tension. We hope, in the future, to update the program to (i) add new models as they become available, (ii) remove any 'bugs' discovered in existing programs, and (iii) provide guidance on the limitations of individual programs.

**Keywords:** pyrometallurgy, molten, slag, physical properties, heat capacity, enthalpy, density, viscosity, thermal conductivity, electrical conductivity, surface tension.

** Introduction**

Slags play an important part in many processes in metal production and refining, coal gasification, continuous casting, etc. The properties of these slags have a significant effect on the performance of the process, e.g. the selection of viscosity and solidification temperature are key to minimizing defects in the continuous casting of steel^{1}, and the electrical conductivity of the slag determines the power supply requirements of many electric smelting processes^{2}. The physical properties of these slags are needed to solve process problems and to improve product quality. However, in recent years there has been an exponential increase in the use of mathematical models to solve process and quality problems, and these models require reliable input data for the physical properties of the slags used in the process. As there are a large number of processes, a wide range of slag compositions used in any one process, and a large number of properties involved, there is a need for a vast amount of property data. The measurement of thermo-physical properties of molten slags is both difficult and expensive. Consequently, there has been a drive to develop empirical rules and models to estimate the various properties of slags from their chemical compositions, as this information is available on a routine basis. Furthermore, the properties of glasses and enamels are also related to those of slags, and so the models developed can, in some cases, be applied to these materials too.

A significant number of models have been reported to estimate the thermo-physical properties of slags. For some properties (e.g. viscosity) there are a large number of reported models available, whereas, in other cases, there are few reported models (e.g. electrical and thermal conductivities) and these are frequently available only for either simple slag systems (e.g. ternaries) or for a limited compositional range.

Silicates are the basis of most metallurgical slags, and their physical properties are very dependent on the silicate structure developed in the molten slag^{3}. The dependence of property on structure is in the hierarchy, viscosity (*η*) > electrical conductivity (κ) and thermal conductivity (*k*) > thermal expansion coefficient (α) > density (ρ) > surface tension (γ) > *C _{p} *and enthalpy. In fact, some workers have proposed that the viscosity is a convenient form to represent slag structure

^{4,5}. In recent years, the capabilities of commercial thermodynamic software packages have expanded from the estimation of thermo-dynamic properties to the estimation of parameters to represent the slag structure and which lead to the calculation of slag densities and viscosities

^{6-10}. The reliability of the property values estimated by the various thermodynamic packages will be evaluated in the future. Most process engineers do not have either the time or the specific expertise required to implement the reported models. Consequently, we have developed here a simple program that will enable engineers to estimate thermo-physical properties from the chemical composition of the slag in mass per cent.

** Objectives of this work **

The objectives of this work are:

(i) To develop a program to calculate various physical properties of slags from their chemical compositions and to make this openly available by placing it on the www.pyrometallurgy.co.za website

(ii) To write a set of notes to accompany the program to aid its use, point out the limitations of predictions of individual models, range of applicability, etc., and provide guidance on which model to use for the slag system in question

(iii) To identify any 'bugs' within the program and correct them continually

(iv) To add new models as they become available.

** Effect of structure and other factors on properties **

A full list of the symbols and abbreviations used is given at the end of this paper.

*Structure of silicate and alumino-silicate slags *

Slags are formed of ions, and slag reactions are electrochemical in nature, involving the exchange of ions. Slags contain two forms of bonds: (i) covalent Si-O bonds that form into chains, rings etc., and (ii) ionic bonds involving cations such as Na^{+} or Ca^{2+} that break the silicate chains, etc., to form Na^{+}-O-bonds^{11}. The silicate structural unit (SiO_{4}^{4-}) consists of one Si^{4+} ion surrounded tetrahedrally by four O_{2}ions (Figure 1)^{11}. Each of these O ions has a negative charge, and so can connect to either another O^{-} ion and thereby add to the network (denoted as a bridging oxygen (BO or Oº)) or to a cation, thereby breaking the chain (denoted as a nonbridging oxygen (NBO or O^{-}). Oxygen ions that are not associated with Si ions are referred to as free oxygen (denoted O^{2-})^{11}.

The structure of a slag can thus be represented by the mole fractions (*X *) of Oº, O^{-}, and O^{2-}present. Consequently, it is customary to divide various constituents into either network formers (e.g. *SiO*_{2}) or network breakers (*CaO, Na*_{2}*O*, etc.). However, when *Al*_{2}*O*_{3} is added to a silicate slag, the *Al *^{3+} ion can be absorbed into the *Si*^{4+} chain but requires charge-balancing cations (e.g. *Na*+ or 0.5*Ca*^{2+}) to create a (*NaAl *)^{4+} ion which must be located near the *Al *^{3+} as shown in Figure 2. Cations on charge-balancing duties cannot act as network-breakers.

*Parameters used to represent structure *

Several parameters have been used to represent the structure of the slag. The earliest models used the basicity (%CaO/ %SiO_{2}) or basicity indices where different weights were given to different basic oxides (e.g. CaO, MgO, or FeO). The most widely used parameter is the ratio of NBO to tetragonal ions (NBO/T), which is calculated using Equation [1]. NBO/T is a measure of the depolymerization of the slag.

where *X *= mole fraction, *f *= fraction of *M*_{2}*O*_{3} acting as a network breaker *MO=CaO, FeO *etc. *M*_{2}*O=Na*_{2}*O *etc. and *M*_{2}*O*_{3}=*Fe*_{2}*O _{3}, Cr*

_{2}

*O*

_{3}etc.

Some workers prefer using the parameter *Q *(defined in Equation [2]) since this provides a measure of the polymerization of the slag that is easier to visualize; *Q *will be used here. Table I provides some examples of the physical reality of various *Q *and (NBO/T ) values.

Other models^{10}calculate the mole fractions of Oº, O^{-}, and O^{2}-(which are usually calculated using the Gaye modification^{12}of the Kapoor-Frohberg cell model^{13}) and relate these to the property values.

Structural data for slags can be determined using Raman spectroscopy; in these cases the structural parameters can be expressed as the mole fractions of *Q*1 (equivalent to 2*CaO.SiO*_{2} or *O*^{2}-units), *Q*_{2} (3*CaO. *2*SiO*_{2} units), *Q*_{3} (*CaO.SiO*_{2} units), and *Q*_{4}(*SiO*_{2} or O^{0 }units).

One difficulty when using both (*NBO/T*) and *Q *is that they do not differentiate between the effects of the different cations on the structure or property (e.g. CaO, MgO, and *Na*_{2}*O*). For this reason, the optical basicity (Λ) has been used to represent the structure as it does differentiate between cations but this parameter (Λcorr) is usually corrected for the cations used in charge-balancing duties^{14}.

Du and Seetharaman^{15,16 } took a different approach and represented the activation energy for viscous flow as a function of the excess free energy (Δ*G*^{xs}) of the molten slag.

*Effect of different cations on properties *

Property values are affected by the nature of the cation as can be seen from Table II. The cations affect the various properties in different ways:

(i) Divalent ions (e.g. Ca

^{2+}) will break two different chains and the resulting ionic bonds hold these two chains together (O^{-}-Ca^{2+}-O^{-}), whereas monovalent ions (O^{-}-Na^{+}) will break only one chain; consequently, for the viscosities of slags with equivalent compositions,ηCaO>ηNa_{2}O.(ii)For equivalent slag compositions (identical

Q) there are twice as many Na^{+}ions as Ca^{2+}ions; since the electrical conductivity is dependent upon the number of ions available, it follows that the electrical conductivity of κNa_{2}O> κCaO.(iii) The electrical conductivity is also dependent upon the mobility of the cations; thus it would be expected that the conductivity would be greater for smaller cations; this is valid for Li, Na, and K ions, but the reverse trend can be seen with Ba, Ca, and Mg. This is probably due to the effects of increased polarization with decreasing ionic radius that could also affect the cationic mobility.

Inspection of Table II indicates that, although property values are affected by different cations, the effect differs according to the property.

*Effect of temperature on physical properties *

Increasing temperature tends to loosen the network structure. Thus properties such as viscosity and thermal conductivity tend to decrease with increasing temperature^{3}. For other properties where the network offers a resistance to movement (e.g. electrical conductivity and thermal expansion), property values tend to increase with increasing temperature.

When looking at the effect of temperature on solid slags it is necessary to divide slags into (i) glassy slags and (ii) crystalline slags. It is customary to express structure-related properties (e.g. viscosity, electrical conductivity) by either the Arrhenius or Weymann equations shown in Equations [3] and [4] respectively, where *A *and *B *are constants and (*B =E/R*) where *E *= activation energy for viscous flow and *R *= gas constant.

When heating a glassy slag from room temperature, there is an abrupt change in property (e.g. *C _{p}*, α) at the glass transition temperature (

*T*) where the slag transforms into a supercooled liquid (

_{g}*scl*). This transition at

*T*is accompanied by a step-like increase in

_{g}*C*and a 3-fold increase in α above

_{p}*T*. The viscosity of the

_{g}*scl*decreases smoothly with increasing temperature above

*T*and is usually expressed in the form of the Vogel-Fulcher relation (Equation [5]) with the viscosity value at

_{g}*T*,

_{g}*η*

*T*(

_{g}*dPas*) = 1013.4 (Figure 3a). Thus for supercooled liquids there are no abrupt changes in property at

*T*(e.g. no enthalpy of fusion). In contrast, when a liquid slag is cooled to the point where crystals are precipitated (known as the break temperature,

_{liq}*T*) there is a sudden increase in viscosity (Figure 3b).

_{br}

Heating a crystalline slag results in no abrupt changes in property at *T _{g}*, and both

*C*and α values tend to be lower than those of glassy slags for the region where

_{p}*T*<

_{g}*T*<

*T*. However, for crystalline slags, there is an abrupt change in both enthalpy (Δ

_{liq}*H*) and volume (Δ

^{fus}*V*) when the crystalline solid transforms to liquid in the fusion region. Thus it is necessary to differentiate between glassy and crystalline slags when estimating property values in the solid state; as a rough rule of thumb, glassy slags occur when Q>2.5 and crystalline slags when Q=<2.

^{fus}

** Modelling thermo-physical properties of slags **

The fundamentals underlying the various methods used to calculate individual properties in the program are described below. All temperatures in this paper refer to the thermodynamic temperature, *K*.

*Liquidus temperature (Tliq) *

It is essential to have a value for the liquidus temperature when modelling the thermo-physical properties of both solid and liquid slags. In the commercial thermodynamic models *T _{liq} *can be calculated for the temperature where the activity of the liquid is unity. This facility is not available for the present software and a value should be provided by the user. If the user does not provide a

*T*value, the model will insert a default value (Equation [6]). This was obtained by carrying out regression analysis on 124

_{liq}*T*values taken from several sources

_{liq}**20,21**. This approach provides general trends for the effect of individual oxides on

*T*. However, the liquidus surface contains 'peaks and valleys' that frequently occur at compositions corresponding to compound formation, e.g. '2CaO.SiO

_{liq}_{2}'. The database contains a number of

*T*values corresponding to such compounds. The mean uncertainty of the estimated

_{liq}*T*value is ± 130K. This is due mainly to large deviations (± 500K) that can occur for compound 'peaks and valleys' which constitute about a quarter of the database; for the remainder of the database, the uncertainty is < ± 100K. For these reasons, users are recommended to provide a value if they have one.

_{liq}*Glass transition temperature (Tg) *

It is also essential to have a value of *T _{g}*for the estimation of the thermo-physical properties of solid, glassy materials. The following equation was derived by regression analysis of experimental data for

*T*and chemical compositions of the slags.

_{g}It should be noted that the constants associated with MgO, MnO, and FeO are 'unrealistic' since they are based on values for slags with ca. 1% of these oxides, which leads to unrealistic estimates of *T _{g}*when applied to slags with much higher concentrations of these oxides.

*Heat capacity (C _{p}) and enthalpy (HT-H298)*

* Crystalline slags *

The heat capacities of slags are little affected by the structure of the slag; thus it is possible to obtain reasonable estimates of *C _{p} *from partial molar

*C*values for individual components (Equation [8] where 1, 2 = CaO,

_{p}*SiO*

_{2}, etc.).

It is customary to express the temperature dependence of *C _{p} *of crystalline slags in the form:

Thus it is possible to derive individual values of the parameters, *a, b*, and *c *by:^{19}

The enthalpy (*H _{T}-H*

_{298}) is given by

Consequently, (*H _{T}-H*

_{298}) for the crystalline state can be calculated from the

*a, b*, and

*c*values.

It takes energy for a crystalline solid to transform into a liquid slag, and this energy (enthalpy of fusion, Δ*H *^{fus}) can be calculated from the entropy of fusion (Δ*S *^{fus}) that represents the structural changes accompanying this transition. The Δ*H *^{fus }can be calculated in the following way^{22}:

*Glassy slags *

The *C _{p} *values for glassy slags are very similar to those of crystalline slags in the range (298K to

*T*), but glassy slags show a stepwise increase of ca. 0.2 kJ K

_{g}^{-1}kg

^{-1}at

*T*(

_{g}*C*usually occurs around 1100 JK

_{p}Tg^{-1}kg

^{-1}(=1.1 kJ K

^{-1}kg

^{-1}). The enthalpy (

*H*

_{T}-H_{298}) value at

*T*(=

_{liq}*Tm*) for the liquid phase must be identical for both glassy and crystalline phases of the same composition. Thus we can calculate

*C*and (

_{p}*H*) values by Equation [14] for the (

_{T}-HT_{g}*T*) range by assuming that

_{g}-T_{liq}*C*= a'+ b'T in this range.

_{p}where Δ = (*H _{T}m -H*

_{298})

*crys*- (

*H*

_{Tg-H}_{298})glass and

*a'=C*. There is no enthalpy of fusion term at

_{p}scl T_{g}*T*for the supercooled liquid.

_{liq}*Liquid slags *

Values for the liquid can be calculated from Equation [8] using the values *C _{p}*(l) for the various slag constituents (e.g.

*CaO, SiO*

_{2}, etc.)

^{22}

*Density (***ρ***) and Thermal expansion (***α,β***) *

*Liquid slags *

The densities of most of the components of the slag tend to be similar, so Keene^{23 }proposed a simple relation for calculating the densities of molten steelmaking slags at 1673K with an uncertainty of ca. ± 5%.

The densities of molten slags (ρ) can be modelled^{19}using the partial molar volumes (*V*, Equations [16] and [17], where *M*= molecular weight) of the various slag constituents. The effect of the *SiO*_{2} and *Al*_{2}*O*_{3} on the molar volume (i.e. structure) of the slag can be represented by simple relations (Equations [18] and [19] respectively) derived from an analysis of density measurements^{22}. The densities are calculated for a reference temperature of 1773K and then adjusted to other temperatures by applying a temperature coefficient of -0.01% K^{-1}. The calculated densities usually lie within ± 2% of measured values.

*Solid slags *

Thermal expansion coefficients (α) for crystalline and glassy slags are very similar for temperatures in the range (298*K *to *T _{g}*). However, for glassy slags, the transition into a super-cooled liquid (

*scl*) is accompanied by a three-fold increase in α, whereas, in contrast, for crystalline slags, α maintains its gradual increase with increasing temperature. Consequently, for any specific temperature in the range (

*T*to

_{g}*T*),

_{liq}*VT*

^{scl}> VT^{crys}. However, at

*T*, there is a sudden increase in volume (Δ

_{liq}*V*

^{fus}) as the crystalline solid transforms into a liquid; in contrast, there is no volume change associated with the transition of a

*scl*into a liquid (i.e. Δ

*V*

^{fus}=0). Thus the enhanced volume change associated with α

^{scl }in the (

*T*to

_{g}*T*) range is offset by Δ

_{liq}*V*

^{fus }.

Crystalline slags

Slags with crystalline phases have a higher density than glassy, amorphous slags, and will show a marked change in density at *T _{liq} *(Δρ

^{fus}.). The slag densities are estimated at 298

*K*by using a very similar method to that adopted for the liquid phase. The values at

*T*>298

*K*were calculated by using an average linear thermal expansion coefficient (α) of 9

**x**10

^{-6}K

^{-1}.

Glassy slags

Molar volumes (*V _{T}*) are calculated in an identical manner to that for crystalline phases up to

*T*. For temperatures above

_{g}*T*, the molar volume of

_{g}*scl*(

*V*) is calculated using Equation [20] and the liquid molar volume at

_{r}*T*(

_{liq}*V*

^{m}) and the molar volume of the

*scl*at

*T*(

_{g}*VT*):

_{g}*Viscosity (**η**) *

Viscosity is a measure of the resistance encountered when moving one set of atoms over a lower layer of atoms. Thus, as the network structure becomes more polymerized (i.e. *Q *increases) the resistance to viscous flow will increase and thus the viscosity (*η*) will increase. Since increased temperature loosens up the structure, the viscosity will decrease with increasing temperature. The viscosity has been used as a measure of the slag's structure^{5}.

The Riboud^{24}, Urbain^{25}, and Iida^{26 }models have been implemented in the software.

*Riboud model*^{24}

This model is simple and is applicable to a wide range of slags. The model divides the slag constituents into five different categories (see below). Additional slag constituents are covered in the software and these have been allocated to the various groups and are denoted by {}. The model applies to the following ranges: *SiO*_{2} (28-48%), *CaO *(13-52%), *Al*_{2}*O*_{3} (0-17%), *CaF*_{2} (0-21%), *Na*_{2}*O*(0-27%).

The temperature dependence is expressed via the Weymann equation {*η*(*dPas *) = *A _{W} T *exp (

*B*)}and the viscosity is calculated from this relation where

_{W}/T*A*and

_{W}*B*are calculated using Equations [25] and [26].

_{W}*Urbain model*^{25 }

This model divides the various slag constituents into the following groups:

Network modifiers:

Here it has been assumed that *Fe*_{2}*O*_{3} and *Cr*_{2}*O*_{3} behave both as network breakers and as amphoterics, where f is the fraction (*f *) behaving as network modifiers and a value *f*=0.6 is assumed.

The Urbain model works predominantly on a basis of *MxO *so this creates extra ions and it is necessary to normalize *X _{G}, X_{M}*, and

*X*by dividing by the term (1+ 0.5

_{A}*XFeO*

_{1.5}+

*XTiO*

_{2}+

*X*

_{ZrO2}+

*X*

_{CaF2}) to give

*X** and

_{G}*X_{M}*X**.

_{A}The model assumes the Weymann relation

Urbain^{25 }found that *A *and *B *were linked through the equation

The *B _{W} *value must be calculated via the equations

where subscript *i *can be 0, 1, 2, or 3 and *a, b*, and *c *are given constants for each case e.g. 0, 1, 2, 3.

Different values for *a, b*, and *c *are given for 0,1,2,3 to calculate *B *values for the *CaO-Al*_{2}*O*^{3+}*SiO*_{2}, *MgO-Al*_{2}*O*^{3+}*SiO*_{2}, and *MnO-Al*_{2}*O*^{3+}*SiO*_{2} systems. In this study, we have modified *X*_{'MnO}' to represent *X _{MnO} +XFeO +X_{NiO} +X_{CrO} *+0.6(

*X*

_{Fe2}

*O*

_{3}+

*X*

_{Cr2}

*O*

_{3}). The global

*B*(

*B*) is given by

_{globa}lThe values calculated with this software are in good agreement with the values cited by Urbain^{24 }but the model does not have *Bi *values covering monovalent oxides (e.g. Na_{2}O).

*Iida model*^{26}

The model^{26} makes use of the basicity index (*B _{i}*) to represent structure:

where *A *= pre-exponential term, *E *= activation energy; *η*_{0}= hypothetical viscosity for each slag constituent (*i*); the parameters *A, E*, and *η*_{0} are all given as functions of temperature.

where *V _{m} *= molar volume for each constituent, and

*H*= 5.1 (

_{i }*T*

^{m })

*and*

_{i }*R** = gas constant

The various constituents are divided into the following categories: (1) acidic (*SiO*_{2}, *ZrO*_{2}, *TiO*_{2}) denoted by subscript A, (2) basic (*CaO, MgO, Na*_{2}*O, K*_{2}*O, Li*_{2}*O, FeO, MnO, CrO, CaF*_{2}, etc.) denoted by subscript B, and (3) amphoteric (*Al*_{2}*O*_{3}, *B*2*O*_{3}, *Fe*_{2}*O*_{3}, *Cr*_{2}*O*_{3}). The basicity index *Bi *is calculated by

where α* _{i}*.=constant for each constituent expressing its relative basicity and %i = mass %.

This was the original Iida model. Subsequently, the model was modified^{26}to account for the amphoterics where their basicity changed according to the temperature. This was done through the modified basicity index (*B _{i}*

^{j}), particularly for

*Al*

_{2}

*O*

_{3}, which was done through back-calculation from experimental viscosity data (Equation [42]). It was concluded that

*Fe*

_{2}

*O*

_{3}and

*Cr*

_{2}

*O*

_{3}worked basically as basic oxides so they appear on the top line. The model gets exceedingly complicated because the modified α (denoted α*) for

*Al*

_{2}

*O*

_{3}was obtained from experimental viscosity data for certain systems, e.g.

*CaO +MgO +Al*

_{2}

*O*

^{3+}

*SiO*

_{2}, and then expressed by Equation [42].

It is difficult to apply the Iida model to systems where there is no experimental data since α* _{i}*.*values are determined for each system/family and there is no general overall value for α

**. For example, one can calculate α*

_{i}** from a, b, and c values for*

_{i}*CaO +MgO +Al*

_{2}

*O*

_{3}+

*SiO*

_{2}but other a, b, c values (leading to a different α

*i.*Al*

_{2}

*O*

_{3}) are given for

*Li*2

*O +Al*

_{2}

*O*

_{3}+

*SiO*

_{2}. The high accuracy claimed with this model comes from its calibration with experimental data for each family of slags.

*Thermal conductivity (k) *

There are no extant models for estimating thermal conductivities but Mills^{3}noted that the thermal conductivity of liquid silicate slags at *T _{liq} *(

*k*

^{m}) increased linearly as (NBO/T) decreased (i.e.

*Q*increased). Phonon conduction is considered to occur as lattice waves. Slags exhibit covalent bonding (in chains, etc.) and ionic bonding. Since

*k*

^{m }increases with increasing

*Q*, it is obvious that thermal conduction is greater along the covalent chain than across the cationic bonds. Consequently, a relationship between thermal conductivity and viscosity might be expected, and recent work

^{4,27 }has shown that such a link exists. Further evidence for this link between thermal conductivity and viscosity lies in the temperature dependence of the thermal conductivities of molten slags that can be satisfactorily represented by the Arrhenius relation, which is widely used for the temperature dependence of viscosity

^{27,28}.

Thermal conductivity (*kR*) measurements at temperatures on glassy and liquid slags contain a significant and unknown contribution from radiation conduction**29**. These *kR *contributions are much smaller in transient hot wire (THW) measurements than for those obtained with the laser pulse (LP) method, because the emitting surface area in LP experiments is 10x that in THW studies. Therefore, only thermal conductivity data from THW studies were accepted here. Typical examples of the temperature dependence of glassy, crystalline, and liquid slags are shown in Figure 5. The maximum in *k *occurs at *T _{g}*and it is our contention that the thermal conductivity in the range (

*T*to

_{g}*T*) can be represented by a Vogel-Fulcher relation (see Equation. [5]), which is widely used for viscosities in this range.

_{liq}

*Liquid slags *

Method 1-relation with viscosity

This uses the relation of thermal conductivity (*k*) of the liquid at *T _{liq} *(

*k*

^{m}) with viscosity. Reported line source data for ln

*k*

^{m}

^{4,27, 30-32 }were correlated with Riboud calculations of ln

*η*

^{m}

^{24}. The data showed some scatter because of (i) uncertainties in the calculated ln

*η*

^{m }values, (ii) experimental uncertainties associated with ln

*k*

^{m }, and (iii) sensitivities of both

*η*

^{m }and

*k*

^{m }in the region around

*T*and the uncertainties in

_{liq}*T*itself. The following relation was obtained.

_{liq}Method 2-relation to both *Q *and ln *η*^{m }

This method is similar to Method 1 in that it uses a relation between ln *η*^{m }and ln *k*^{m}. The data is shown in Figure 6, and has been expressed as an exponential relation.

Then ln *k*^{m }can be calculated using the equation

The values of *k*^{m }obtained are valid only for the range *Q *= 2 to 3.2. Care should be taken particularly for slags with *Q *> 3.2 and for slags with high *Al*_{2}*O*_{3} contents (e.g. 45% CaO + 10% *SiO*^{2+} 45% *Al*_{2}*O*_{3} which would show a high *Q *value but, in actual fact, is largely made up of calcium aluminates and will show relatively low *η*^{m }and *k*^{m }values).

*Method 3-relation with Q *

The relation between ln *k*^{m }and *Q *is shown in Figure 6b and it can be seen that it has similarities to the equation between ln *η*^{m }and *Q *(Figure 6a, Equation [47]). There is significant dispersion in the scattered data for *k*^{m }but the upward trend with increasing *Q *is obvious. Unfortunately, the experimental values for *k*^{m }lie within the range *Q *= 2 to 3.2. The following relation was derived:

This equation is valid only for the range Q = 2 to 3.2.

The temperature dependence would be expected to follow an Arrhenius-type relation because of the correlation between *η* and *k*.

There are insufficient data to determine B values as a

function of Q, so an approximate relation for *dk/dT *was obtained but should be used with caution.

*Solid slags *

Method 1

The experimental data for ln *k*_{298} and ln *k _{Tg} *show appreciable scatter when plotted vs

*Q*, but they do show k increasing with increasing

*Q*. The data were assumed to follow an exponential relation similar to that used for ln

*η*

^{m }(Equation [47], Figure 6b). It was noted that the experimental data for solid

*Li*2

*O*- containing slags were significantly higher than the curve values and so a correction for

*Li*

_{2}

*O*was added to the exponential relation:

This equation should not be used out of range i.e. *Q *> 3.3 or *Q *< 2.

Values can be calculated at other temperatures by linear interpolation between 298 *K *and *T _{g}*, and values have also been calculated in the same manner for

*T*to

_{g}*T*(but this should in fact follow a Vogel-Fulcher relation, Equation [5], but there was not enough data to establish the necessary constants). This procedure should be accepted only for slags in the range

_{liq}*Q*= 2 to 3.3.

Method 2

A second method is also shown. It was found that for slags with *Q *around 2.5 the thermal diffusivity between 298 K and *T _{g}*tends to have a constant value

*a*= 4 x 10-7 m2s-1. Values calculated from

*k*= a.

*C*.ρ using the calculated values of

_{p}*C*and ρ are also shown. These values apply only when

_{p}*Q*= 2 to 3.

*Crystalline slags *

Nishioka *et al*.^{31 }report that *k *values for crystalline slags are between 1.5 and 2.5 times higher than that for the glassy slag and are independent of temperature (Figure 5). Thus it has been assumed that *k _{cryst} *= 2

*k*and is independent of temperature from 298 to 1000K.

_{glass }** Electrical conductivity (**κ

*) of liquid slags*The electrical conductivity involves the movement of cations under the influence of an applied field. In this case, the polymeric silicate network hinders the movement of the cations. Thus the factors affecting the electrical conductivity (κ) are (i) the concentration, charge, and size of the cations and (ii) the polymeric network (as represented by *Q*) which hinders the movement of cations. Smaller cations are more mobile, but the increased mobility of smaller cations is offset by the increased polarization in the M-O bond. For instance, for *M*_{2}*O-SiO*_{2} slags, the hierarchy is κ* _{Li}*> κ

*> κ*

_{Na}*i.e. decreasing with increasing cation radius. However, the trend is reversed for the*

_{K }*M*O-

*SiO*

_{2}slags where κ

*> κ*

_{Ba}*> κ*

_{S}r*Ca*> κ

_{Mg}^{5,16 }i.e. κ decreases as (

*z/r*2), possibly because of the higher charge (

*z*) in this slag series which would be expected to increase polarization. As increasing temperature gradually loosens up the structure, the hindrance to the movement of cations by the silicate network is reduced, and hence the conductivity increases with increasing temperature. It should also be noted that factors causing an increase in electrical conductivity would cause a decrease in thermal conductivity. The electrical conductivities of slags containing Li

^{+}, Na

^{+}, and

*K*

^{+ }are significantly higher than those containing

*Ca*

^{2+},

*Mg*

^{2+}, etc., because, for an equivalent slag structure (

*Q*),

*n*=

_{Na+}

_{nCa }^{2+}where n = number of cations present.

The electrical conductivities have been estimated using recent papers by Zhang and Chou^{5, 36}, both of which involve using the relation between conductivity and viscosity.

Method 1

Chou and Zhang^{35, 5 }proposed a relationship for slags containing *CaO, MgO, Al*_{2}*O*_{3}, and *SiO*_{2}.

Zhang calculated values of ln *k *using this equation and parameters to calculate the viscosity. These parameters are not generally available for all slag systems, so ln *η* values were calculated here using the Riboud, Iida, and Urbain models.

Method 2

Zhang *et al. *^{5,36 }reported that the following equations apply for *M*_{2}*O-SiO*_{2} (where *M*_{2}*O =Na*_{2}*O *etc.) and *MO-SiO*_{2} (e.g. CaO) systems.

The model assumes (i) that cations on charge balancing duties are not available for electrical conduction and (ii) the charge balancing is directly related to the concentrations of the various cations present (i.e. a statistical distribution). The differences between Equations [53] and [54] were attributed to the fact that there are twice as many Na^{+ }ions as Ca^{2+}for equivalent compositions. Here the model calculates the ratio (*r**) of *M*^{+ }ions present in the slag (Equation [55]) and then calculates the conductivity using Equation [56] (which was derived from the product of *r** multiplied by the differences between Equations [53] and [54]). It is recommended that this method be used when the slag contains *M*^{+ }ions e.g. *Na*^{+}, *K*^{+}, or *Li*^{+}.

*Surface tension (***γ***)*

Method 1

Surface tension is not a bulk property but is a surface property. The free surface of the molten slag contains higher concentrations of the constituents with lower surface tension. Values of surface tension and the trends in surface tension can be calculated using a partial molar approach (Equation [57]), and are shown in Figure 7^{37}. The model divides slag components into two classes: (i) oxides with higher surface tension values where values of *X*_{2} γ_{2} tend to be similar to that shown in Figure 7a and (ii) components with lower γ values (surfactants *B*_{2}*O*_{3}, *CaF*_{2}, *Na*_{2}*O, K*_{2}*O, Fe*2*O*_{3}, and *Cr*_{2}*O*_{3}) where values of *X*_{2} γ_{2} tend to be similar to those shown in Figure 7a.

These surface active components cause a rapid decrease in surface tension and, in these cases, *X*_{2} γ_{2} and can be represented by two curves determined by a critical point N (corresponding to the minimum in Figure 7b) (i) for <N as a polynomial by *X*_{2}. γ_{2} = *a +bX +cX*^{2 }and (ii) >*N *by *X*_{2} γ_{2} = *a*' +*b'X *(which are shown by the dashed line in Figure 7b). The method has the advantage that it can easily be applied to multicomponent industrial slags, but uncertainties in the estimated values are ca. ± 10%. The biggest problem with this method lies in its inability to deal with two, or more, surface-active components simultaneously (e.g. *CaF*_{2} and *B*_{2}*O*_{3}) and, for these conditions, Method 1 tends to exaggerate the decrease in γ in these cases (i.e. _{γ}_{meas > }_{γ}_{calc}).

**Method 2 **

Method 2 was devised to address this problem. It was assumed that (i) occupation of the surface occurs in the hierarchy (of lowest surface tension) *B*_{2}*O*_{3}> *K*_{2}*O*>*Na*_{2}*O*>*CaF*_{2}> *Fe*_{2}*O*_{3}>*Cr*_{2}*O*_{3} and (ii) surface saturation occurs at N=0.12. The various surfactant contributions were calculated from the various oxides until the value N=0.12 was attained. When N>0.12 all contributions from these surfactants were taken as positive contributions (=X_{i}_{γ}* _{i}*) i.e. they were not considered to be present in the surface layer and thus contributed normally to the bulk surface tension. The values calculated using Method 2 are preferred.

Most slags have negative temperature coefficients (*d*γ/*dT *)^{20}. However, it is known that slags with high *SiO*_{2} contents exhibit positive (*d*γ/*dT *) values^{38}. This is related to the fact that *SiO*_{2} has a relatively low surface tension value. Ideally, (*d*γ/*dT *) could be calculated on a partial molar basis but we do not have (*d*γ/*dT *) data for all the slag constituents. Consequently, a constant value of (*d*γ/*dT *) has been assumed. However, (*d*γ/*dT *) values were calculated on a partial molar basis for those slag constituents where data were available (e.g. Equation [58]).

**Using the software **

The intention behind this work is to provide a practical set of models that can be used by process engineers as working tools in everyday pyrometallurgical work. The software is available from http://www.pyrometallurgy.co.za/KenMills/, initially in the form of a spreadsheet, but other interactive calculation methods may be available in due course. The attractive feature of having the software available on a website is that the latest version can always be available to all users, and the burden of sending out software updates is relieved. Specific detailed instructions for the use of the software are available on the same website, so are not repeated here. The software requires the specification of the composition and liquidus temperature of the slag. If the glass transition temperature for glassy slags is known, this can also be specified; otherwise this will be estimated by the software.

** Acknowledgements**

The authors would like to express their thanks to the following for their valuable comments during the preparation of this paper: Professor Seshadri Seetharaman (KTH, Stockholm), Professors M. Susa and M. Hiyashi (Tokyo Institute of Technology), Professor K. Morita and Dr Y. Kang (University Tokyo), and Professor K.C. Chou and Dr G.H. Zhang (University of Science and Technology Beijing).

** References**

1. MILLS, K.C. and FOX, A.B. The role of mould fluxes in continuous casting- so simple yet so complex. *ISIJ International*, vol. 43, 2003. pp. 1479-1486. [ Links ]

2. JONES, R.T., REYNOLDS, Q.G., CURR, T.R., and SAGER, D. Some myths about DC arc furnaces. Southern African Pyrometallurgy 2011. Jones, R.T. and Den Hoed, P. (eds.). *Southern African Institute of Mining and Metallurgy*, Johannesburg, 6-9 March 2011. pp.15-32. [ Links ]

3. MILLS, K.C. The influence of structure on the physico-chemical properties of slags. *ISIJ International*, vol. 33, 1993. pp. 148-156. [ Links ]

4. HAYASHI, M., ISHII, H., SUSA, M., FUKUYAMA, H., and NAGATA, K. Effect of ionicity of nonbridging oxygen ions on thermal conductivity of molten alkali silicates. Physics and Chemistry of Glasses -*European Journal of Glass Science and Technology Part B*, vol. 42, no. 1, 2001. pp. 6-11. [ Links ]

5. ZHANG, G.H. and CHOU, K.C. Simple method for estimating the electrical conductivity of oxide melts with optical basicity. *Metallurgical and Materials Transactions B*, vol. 41, no. 1, February 2010. pp. 131-136. [ Links ]

6. NATIONAL PHYSICAL LABORATORY. MTDATA. www.npl.co.uk/mtdata/ [ Links ]

7. FactSage www.factsage.com [ Links ]

8. THERMO-CALC SOFTWARE www.thermocalc.com [ Links ]

9. TASKINEN, P., GISBY, J., PIHLASALO, J., and KOLHINENT, T. Validation of a new viscosity database with industrial smelting slags. *Proceedings, European Metallurgical Conference EMC2009*, Innsbruck, Austria, 28 June-1 July 2009. GDMB, Clausthal-Zellerfeld, Germany. vol. 3, pp. 915-930. [ Links ]

10. LING ZHANG, SUN, S., and JAHANSHAHI, S. Review and modeling of viscosity of silicate melts: Part I. Viscosity of binary and ternary silicates containing CaO, MgO, and MnO. Part II. Viscosity of melts containing iron oxide in the CaO-MgO-MnO-FeO-Fe_{2}O_{3}-SiO_{2} system. *Metallurgical and Materials Transactions B*, vol.29, 1998. pp. 177-186 and 187-195. [ Links ]

11. MYSEN, B.O. Structure and Properties of Silicate Melts. Elsevier, Amsterdam, 1988. [ Links ]

12. GAYE, H and WELFRINGER, J. Modelling of the thermodynamic properties of complex metallurgical slags. Proceedings of the. 2nd International Symposium on Metallurgical Slags and Fluxes. Fine, H.A. and Gaskell, D.R. (eds.). *The Metallurgical Society of AIME*, Warrendale, PA, 1984. pp. 357-373. [ Links ]

13. KAPOOR, M.L. and FROHBERG, M.G. Chemical Metallurgy of Iron and Steel. *Iron and Steel Institute*, London, 1971. pp. 17-21. [ Links ]

14. MILLS, K.C. and SRIDHAR, S. Viscosities of ironmaking and steelmaking slags. *Ironmaking & Steelmaking*, vol.26, no. 4, 1999. pp. 262-268. [ Links ]

15. DU SICHEN, BYGDEN, J., and SEETHARAMAN, S. A model for estimation of viscosities of complex metallic and ionic melts. *Metallurgical and Materials Transactions B*, vol. 25, 1994. pp. 519-525. [ Links ]

16. DU SICHEN, BYGDEN, J., and SEETHARAMAN, S. Estimation of the viscosities of binary metallic melts using Gibbs energies of mixing. *Metallurgical and Materials Transactions B*, vol. 25, 1994. pp. 589-595. [ Links ]

17. MORINAGA, K., SUGINOHARA, T., and YANAGASE, T. The electrical conductivity of CaO-SiO_{2}-Fe_{2}O_{3} and Na_{2}O-SiO_{2}-Fe_{2}O_{3} melts. *Journal of the Japan Institute of Metals*, vol. 39, 1975. pp. 1312-1317. [ Links ]

18. BOCKRIS, J. O'M., MACKENZIE, J.D., and KITCHENER, J.A. Viscous flow in silica and binary liquid silicates. *Transactions of the Faraday Society*, vol. 57, 1955. pp. 1734-1748. [ Links ]

19. URBAIN, G. *Revue Internationale des Hautes Temperatures et des Refractaires*, vol. 11, 1974. p. 133. [ Links ]

20. Slag Atlas. VDEh , Dusseldorf, 1995. [ Links ]

21. FORSBACKA, L. Experimental study and modelling of viscosity of chromium containing slags. PhD thesis, Helsinki University of Technology, Finland. 2007. (TKK-MT-196). [ Links ]

22. MILLS, K.C. and KEENE, B.J. Physical properties of BOS slags. *International Materials Reviews*, vol. 32, 1987. pp. 1-120. [ Links ]

23. KEENE, B.J. Ref 20 page 44. [ Links ]

24. RIBOUD, P.V., ROUX, Y., LUCAS, L.D., and GAYE, H. Improvements of continuous casting powders. Fachber. *Huttenpraxis Metallweiterverarb*., vol. 19, 1981. pp. 859-869. [ Links ]

25. URBAIN, G. Viscosity estimation of slags. *Steel Research*, vol. 58, 1987. pp. 111-116. [ Links ]

26. IIDA, T. Proceedings of the Mills Symposium, London, 2002. Aune, R.E. and Sridhar, S. (eds.). *National Physical Laboratory*, Teddington, UK. see also CAMP-ISIJ 2002 15(4). [ Links ]

27. MILLS, K.C., YUAN, L., and LI, Z. Unpublished results. *Imperial Colleg*e, London, 2011. [ Links ]

28. KANG, Y. and MORITA, K. Thermal conductivity of the CaO-Al_{2}O_{3}-SiO_{2} system. *ISIJ International*, vol. 46, 2006. pp. 420-426. [ Links ]

29. GARDON, R. *Proceedings of the 2nd International Thermal Conductivity Conference*, Ottawa, Canada, 1962. pp. 167. [ Links ]

30. NAGATA, K. *Proceedings Of the 2nd International Conference on Metallurgical Slags and Fluxes*, Lake Tahoe, NV, USA (1984) pp. 875. [ Links ]

31. NISHIOKA, K., MAEDA, T., and SHIMIZU, M. Application of square-wave pulse heat method to thermal properties measurement of CaO-SiO_{2}-Al_{2}O_{3} system fluxes. *ISIJ International*. vol.46, 2006. pp. 427-433. [ Links ]

32. OZAWA, S., ENDO, R., and SUSA, M. Thermal conductivity measurements and prediction for R_{2}O-CaO-SiO_{2} (R=Li, Na, K) slags Tetsu-to-Hagane, vol. 93, no. 6, 2007. pp. 416-423. [ Links ]

33. OZAWA, S. and SUSA, M. Effect of Na_{2}O additions on thermal conductivities of CaO-SiO_{2} slags. *Ironmaking and Steelmaking*, vol. 32, no. 6, 2005. pp. 487-493. [ Links ]

34. SUSA, M., KUBOTA, S., HAYASHI, M., and MILLS, K.C. Thermal conductivity and structure of alkali silicate melts containing fluorides. *Ironmaking and Steelmaking*, vol. 28, no.5, 2001. pp. 390-395. [ Links ]

35. SUSA, M., WATANABE, M., OZAWA, S., and ENDO, R. Thermal conductivity of CaO-SiO_{2}-Al_{2}O_{3} glassy slags: its dependence on molar ratios of Al_{2}O_{3}/CaO and SiO_{2}/Al_{2}O_{3}. *Ironmaking and Steelmaking*, vol. 34, no. 2, 2007. pp. 124-130. [ Links ]

36. ZHANG, G.H., YAN, B.J., CHOU, K.C., and LI, F.S. Relation between viscosity and electrical conductivity of silicate melts. *Metallurgical and Materials Transactions B*, vol. 42, no.2, 2011. pp. 261-264. [ Links ]

37. MILLS, K.C. Estimation of physicochemical properties of coal slags. ACS Symposium Series 301, Mineral Matter and Ash in Coal. Vorres, KS (ed.). *American Chemical Society*, Washington DC, 1986. pp. 195-214. [ Links ]

38. MUKAI, K. and ISHIKAWA, T. *Nippon Kinzoku Gakkaishi*, vol. 45, no. 2, 1981. pp. 147. [ Links ]

This paper was first presented at the Southern African Pyrometallurgy Conference, 6-9 March 2011, Misty Hills, Muldersdrift.

*© The Southern African Institute of Mining and Metallurgy, 2011. SA ISSN 0038-223X/3.00 + 0.00. *

**Symbols, abbreviations, units **

α | = Thermal diffusivity (m2s-1) |

C_{p} | = Heat capacity (JK^{-1}/ mol^{-1 }or JK^{-1}/kg^{-1}) |

E | = Activation energy (kJ mol^{-1}) |

f | = Fraction of M^{3+} ions acting as network breaker |

(H_{T}-H_{298}) | = Enthalpy (Jmol^{-1 }or Jkg^{-1}) |

ΔH^{fus} | = Enthalpy of fusion (Jmol^{-1 }or Jkg^{-1}) |

k | = Thermal conductivity (Wm-1K^{-1}) |

(l) | = Liquid phase |

M | = Molecular weight (g mol^{-1}) |

p | = Partial pressure (atm, bar) |

Q | = 4-(NBO/t)= measure of polymerization |

R | = Gas Constant= 8.314 (J mol^{-1} K^{-1}) |

r | = Cation radius (10^{-10}m) |

(s ) | = Solid phase |

scl | = Supercooled liquid |

T | = Temperature (K or ºC) |

T_{g} | = Glass transition temperature (K or ºC) |

T_{liq} | = Liquidus temperature (K or ºC) |

V | = Molar volume (m3 mol^{-1}) |

Z | = Charge on cation (=2 for Fe^{2+}) |

α | = Linear thermal expansion coefficient (K^{-1}) |

β | = Volume thermal expansion coefficient (K^{-1}) |

γ | = Surface tension (mNm^{-1}) |

γms | = Slag/metal interfacial tension (mNm^{-1}) |

κ | = Electrical conductivity (Ω^{-1}m^{-1}) |

Λ | = Optical basicity |

η | = Viscosity (Pas or dPas) |

θ | = Contact angle (º) |

ρ | = Density (kgm-3) |

BO | = Bridging oxygen (Oº) |

NBO | = Non-bridging oxygen ((O-)) |

**Subscripts **

br | = Break (temperature) |

liq | = Liquidus |

m^{(subscript)} | = Metal |

s | = Slag |

sol | = Solidus |

**Superscripts **

(superscript) ^{m} | = Value at T _{liq} |