SL - PAP - G.Urbain - 1987 - Viscosity Estimation of Slags
SL - PAP - G.Urbain - 1987 - Viscosity Estimation of Slags
The viscosity of silicate melts and slags is a property difficult to measure accurately using sophisticated apparatus, it is both
time consuming and expensive. The use of models to estimate the viscosity is the first step to be carried out, mainly with
metallurgical slags involving numerous components. The proposed model is an experimental one based on accurate data for
pure binary and ternary melts. From this data, using a simple hypothesis on the behaviour of unusual compounds, it is possible
to estimate the viscosity of a complex melt. This estimation is used as a first approximation and can be amended, if necessary,
by only a few direct measurements.
Berechnung der Viskositaten von Schlacken. Die Viskositiitsmessungen von Schlacken und flussiqen Silikatschmelzen sind
zeitaufwendig und schwierig. Urn genaue Messungen zu erhalten, benotiqt man komplizierte Geriite, da man bei hohen Tem-
peraturen arbeiten rnuls, Fur viele Anwendungen und ganz speziell fur die metallurgischen Schlacken mit zahlreichen Kompo-
nenten konnen berechnete Werte genugen. Das Versuchsmodell fur diese Werte beruht auf Viskositiitsmessungen, die an
biniiren und terniiren Schmelzen vorgenommen wurden. Die berechneten Werte beziehen sich auf die Eigenschaften der Kom-
ponenten, fur die es keine direkten Messungen gibt. 1m allgemeinen ist der berechnete Wert ein erster Niiherungswert, der
allenfalls durch eine bestimmte Anzahl von direkten Messungen verbessert und vervollstiindigt werden muB.
For the moment, a coherent theory of the liquid state does The probabilities P; and P; are related to the variables of
not exist and theoretical models are not within our reach, state and to the structure of the liquid (degree of polymeri-
especially for associated liquids such as silicates or phos- zation). The "jumping" probability P; is deduced from the
phates. A pragmatic way is to use empirical models based on balance between the kT and the energy level of the potential
measured viscosities. Two models are favoured by geolog- well E w • Following a statistical approach proposed by H. D.
ists: those proposed in 1972 by Y. Bottinga I) and by H. Weymanrr') the equation of viscosity is:
Shaw"). Both used nearly the same published viscosity data 1 1
(mostly geological melts) and empirical rules to calculate the Jl = (kTI E w)2. (21tkT)2. V -2/3. P; - I. exp(EwlkT), (I)
viscosities. The heterogeneity of the accuracies of these data
m and v are mass and volume of the structural unit, E; en-
present some difficulties for application to slags.
ergy of the well, k and T Boltzmann's constant and tempera-
The proposed model originated in a collective work carried ture in K, P; is the "hole" probability connected with the
out in the sixties in Irsid's Physical Chemistry Department structural model of the liquid.
under the direction of Professor P. Kozakevitch. The probability P; is proportional to the concentration N; of
"holes" given at T by the equilibrium:
(2)
Background for estimations
!i.Hv and !i.Sv are the partial molar enthalpy and entropy as-
Viscous flow is a momentum transport process in the fluid sociated to the "hole" formation. Using equations (I) and
phase, just as thermal conductance is a transport of energy. (2) it is possible to obtain a two parameter expression for the
On a macroscopic scale, shear viscosity is the ratio between viscosity:
shear stress and the rate of flow. Units for dynamic viscosity
Jl=A· T·exp(1000·BIT). (3)
are:
The pre-exponential parameter A is a function of m, E w , v
- the S.1. unit Poiseuille = Pas = N s/rn? or and !i.Sv. Parameter B is a function of E; and !i.Hv.
- the CGS unit poise = 0.1 Pas.
In this publication the unit "poise" is used as the purpose is Relation between A and B. Parameters A and B expressed in
to give useful information needed by practicians. molar quantities are from (1):
A =Ao·exp( -!i.SJR) (4)
Expressions for the viscous flow. On a molecular scale, vis-
I
cous flow involves relative movements of the structural ele-
with Ao=R ·(2MIE)2 ·(11 V)2/3 .(11 N)I/3,
ments of the liquid. Two conditions must be realised for
such a displacement: B=(E+!i.Hv)IR. (5)
- an energetic condition: the "jumping" probability Pe, As the "hole" equilibrium is governed by the free energy
- a geometric condition: the "hole" or "free volume" prob- !i.Gv and the equation:
ability r..
!i.GvCT)=!i.Hv- T·!i.S"
Simultaneous "satisfaction" of these two conditions is nec-
it is always possible for T= T; to have !i.Gv=O, then:
essary to move some structural units when external stress is
applied. The fluidity, the reciprocal of viscosity, could be !i.Sv=!i.HvITc.
written:
Combination with equations (4) and (5) gives:
InA =lnAo+EI(RTc)-IOOOBITc
Dr-Inq. Georges Urbain, Oir. Recherche C.N.R.S., Universite de Per- as A o, E and T; are constant for a given liquid, a simple for-
pignan, France mulation is:
'"
....
. ......
' .
07 • ". ". '" ". 1D. .
, : ", ". ", " . . '38 ,'
.: \ :.. ...•... ' .)6 .
0,6 .: ~:. .... :.. .. ,. '. . 31. .'
.......
0.. 0. .32
.... ' ... '. . 30
0
: : :. °•• ••
, 28
: ": : .... 0 ••
"
. 26 .'
'.
-,
-, 2S
21.
.:
"
.'
'.
'. . \16 ':~6 ".20 ···.22 ' •. 23 ... 22,: 20.:
•.12 ".14 ". ".
1,0 0,9 0,8 0;1 D,6 o.S 0,4
N(CoOI
equation is convenient. For alumina-oxide melts the small Table 2. Parameters a, band c of the equation B=a+b·a+c·a 2 for
number of experiments limits the equation to a parabolic ap- binary melts alumina-oxide. a is equivalent to the mole fraction of the
oxide
proximation. Pure oxides, such as MgO, CaO, MnO, PbO,
Tafel 2. Parameter a, b und c der Gleichung B=a+b·a+c·a 2 wobei
are not accessible to measurement. The extrapolated value a=Molenbruch des Oxid fur binare Schmelzen ist
of B obtained from binaries is in agreement with the theore-
tical value given by J. O. M. Bockris") using the "hole" mod- binary melt ref. a b c iv
el. All information for each system is given in sequence:
Alumina + MgO 3 13.2 16.9 -18.6 O<a<1
I number of measured melts, Alumina + CaO 5 14 13.2 40 -42.5 O<a<1
ii origin of data, Alumina + MnO 4 15 13.2 20 -25.6 O<a<1
III equation B - X (mole fraction),
iv validity limits,
v B values for pure oxides.
Binary melts silica-oxide. For all these composrtions Table 3. Numerical values of coefficients a(i), b(i) and c(i) for the four
equations (8) with i=O, 1,2 and 3 for the three oxides MgO, CaO and
X = N(SiO z). For pure silica B = 64.1 5 ) , see table 1.
MnO
Tafel 3. Parameter a, b und c der 4 Gleichungen (8), wobei i=O, 1,2
Binary melts of type alumina-oxide. Melted pure alumina has und 3 der Schmelzen der 3 ternaren Systeme SI0 2-AI20 3-MO. Ka-
a measured B= I3.2±0.5 and a calculated B=8.6 4 ) , for tion M ~ Mg, Ca oder Mn
Tm = 2054 °C, this difference is discussed in 6). For all sys-
tems alumina-oxide the composition parameter a defined i a (i) b(i) c(i)
all Mg Ca Mn Mg Ca Mn
before a = MO/(MO + A z0 3 ) is equivalent to the mole frac-
tion MO of the oxide: a = N(MO), see table 2. 0 13.2 15.9 41.5 20.0 -18.6 -45 -25.6
1 30.5 -54.1 -117.2 26 33 130 -56
Ternary melts. Only three ternary systems are studied in the 2 -40.4 138 232.1 - 110.3 -112 -298.6 186.2
3 60.8 -99.8 - 156.4 64.3 97.6 213.6 -104.6
whole composition range: silica-alumina-MO oxide with
Table 4. Example of B calculations for three melts M = Mg, Ca or Mn. The value of parameter B is obtained
Tafel 4. Beispiele der Berechnung von B fur drei Schmelzen from equations (7) and (8) for the composition given by X
and a. Table 3 recorded the parameters a(i), b(i) and c(i) of
sample no. 29 30 31
equation (8) for the three cations.
composition:
silica 44.9 42 33.14 Figure 2 gives iso-B curves for melts silica-alumina-CaO is-
alumina 29.9 19.8 38.36 sued from 10).
CaO 25.2 38.2 29.5
For the following numerical applications the use of a simple
X 0.5016 0.4439 0.3721 basic program is a good solution. A copy on floppy disc for
a 0.6051 0.7781 0.5830
Apple computer is available on request.
B 28.7 24.23 23.38
Numerical applications
modifier
1) all iron (Fe + + +) is equivalent to AI + + +,
iron behaviour amphoteric
tetrahedral cations T =0.6469 T =0.6469 2) all iron is Fe + + .
amphoteric cations A =0.1335 A =0.113
modifier cations M=0.2328 M=0.2922 For a summary of the calculations see table 9.
total: 1.0133 1.0528
The calculation of X and a is obvious following their defini-
calculated: X 0.6385 0.6145 tion and paying attention to the fact that
a 0.6355 0.7198 N(TiO z)=2·N(MO) and that (Fez03)=3·N(MO). For the
B(CaO) 33.75 31.12 two ternaries silica-alumina-CaO and silica-alumina- MgO
B(MgO) 32.24 29.8 the B values are obtained with equations (7) and (8) using a
B (mean) 33.09 30.56
micro computer.
Table 10. Measured and estimated viscosities x 10' in poise Table 12. Summary of calculations
Tafel 10. Gemessene und geschiitzte Viskositiiten x 10' in Poise Tafel 12. Zusammenfassung der Ergebnisse
Again, the calculations are summarized, see table 15. Data under well defined conditions. The different possibilities of
published by the manufacturer are interpolated from viscos- this model are:
ity measurements made in Irsid laboratories. The good
agreement obtained here is a control of the hypothesis made - estimations used as a first approximation,
for the calcium fluoride behaviour in this slag. Some other - an aid for planning experiments and avoidance of numer-
examples can be found in 12). ous preliminary trials,
- a test for structural hypothesis, and
- viscosity-composition-temperature equations compatible
Conclusions with micro computers.
The scope of such an experimental model has to be as sim- (A 00065; received: 8. August 1986
ple as possible in order to give estimations of slag viscosities in revised form: 17. November 1986)
References
') Bottinqa, Y.; Weill. D. F.: Am. J. Sci. 272 (1972), p. 438/75. '0) Urbain. G.; Cambier, F.; Deletter, M.; Anseau, M. R.: Trans. J. B. Cer-
') Shaw. H.: Am. J. Sci. 272 (1972), p. 870/93. am. Soc. 80 (1981), p. 139/41.
J) Weymann, H. D.: Kolloid Z. Polymer 181 (1962), p. 131/37. ") Ikeda. K.; Tamura. A.; Shiraishi, Y.; Saito. T.: Bul. Res. Inst. Min.
4) Bockris, J. O. M.; Reddy. A. K. N.: Modern Electrochemistry, v. 1, Tohoku 29 (1973), p. 24/36.
chapter 6. Plenum Press, New York, 1970. 12) Urbain. G.: J. Mat. Education 1985 Vol. 7 No.6, p. 1007/78.
5) Urbain. G.; Bottinqa, Y.; Richet, P.: Geochim. et Cosmo. Acta 46 13) Urbain. G.: Rev. Int. H. Temp. Refrat. 21 (1984), p. 107/11.
(1982), p. 1061/72. 14) Urbain. G.: Rev. Int. H. Temp. Refract. 20 (1983), p. 135/39.
6) Urbain. G.: Rev. Int. H. Temp. Refract. 19 (1982), p. 55/57. IS) Bernard. G.; Riboud, P. v.; Urbain. G.: Rev. Metallurgic 78 (1981)
7) Urbain, G.; Millon. F.; Cariset, S.: C. R. Acad. Sc. Paris 290 (1980), p. No.5, p. 421134.
137/40. 16) oa« I. N.; Cambier, F.; Deletter, M.; Urbain. G.: B. Ceram. Trans. J.
') Bockris, J. O. M.; Mackensie, J. D.; Kitchener, J. A.: Trans. Farad. 83 (1984), p. 105/12.
Soc. 51 (1955), p. 1734/48. ") Piwinski, A. J.; Weed. H. c.: Thermochimica Acta 37 (1980), p. 189/95.
9) Panek. z.. Silicaty 22 (1979), p. 97/102. IS) Murase. T.; Mcliirney, A. R.: Science 167 (1970), p. 1491/93.
Letter to the Editor from Professor Dr.-Ing. H. Berns, Ruhr- tool steel is of a hot working type and in service carefully
Universitiit Bochum, FRG preheated to > 150°C. This enlarges the plastic zone size
several times, thus enhancing void formation and dimple
Concerning: Pacyna, J.: Effect of non-metallic inclusions fracture in CT specimens. Under these conditions the frac-
(NMI) on fracture toughness of tool steels. steel res. 57 ture toughness, in transverse specimens at 150°C, is for in-
(1986) No. 11, p. 586/92. stance raised by about 50% if the sulphur content is low-
ered from 0.0070 to 0.0027% via a Ca-treatment*). As the
The author's last conclusion is that "application of argon- author tested longitudinal specimens at room temperature
ing, electroslag remelting or vacuum which radically de- only, his results are misleading as far as the effect of refining
crease the NMI fraction is useless as regards fracture tough- processes on the fracture toughness is concerned. For ESR
ness; there is only an increase in the cost of steel produc- steel with a hardness of - 460 HY a 40% increase was mea-
tion". sured at 20°C in longitudinal and transverse specimens*).
(ALOOOl)
In my opinion this is a strong statement based on insuffi-
cient evidence. The three melting procedures are compared
*) Wendl, F.: Influence of production procedures on the microstructure
at a hardness (500 to 575 HY) and a testing temperature and toughness of hot working too) steels with 5% Cr. Dr.-lng. thesis,
(presumably 20°C) where quasicleavage prevails. Here, of Ruhr-Universitat Bochum, 1985; see also: Fortschr. Ber. VOl Reihe 5,
course, little influence of NMI is to be expected. But this Nr. 91, VOl-Verlag Diisseldorf 1985.