ISIJ International
Online ISSN : 1347-5460
Print ISSN : 0915-1559
ISSN-L : 0915-1559
Regular Article
Calculation of Physical Properties for Use in Models of Continuous Casting Process-Part 2: Steels
Kenneth C. Mills Shyamprasad KaragaddePeter David LeeLang YuanFatemeh Shahbazian
Author information
JOURNAL OPEN ACCESS FULL-TEXT HTML

2016 Volume 56 Issue 2 Pages 274-281

Details
Abstract

The objective of the present study was to calculate physical property values for steels from their chemical compositions for subsequent use in mathematical models of the fluid flow, heat transfer and shell solidification in the continuous casting mould. Values of the following properties of steels are calculated for temperatures between 298 K and 2000 K; Heat Capacity (Cp) Density (ρ) Thermal conductivity (k) and diffusivity (a) Electrical resistivity (R) Viscosity (η) Surface (γm) and Interfacial tension (γmsl). In addition temperatures of transitions (Liquidus Tliq, Solidus Tsol) and various solid state transitions were also calculated. Ferritic and austenitic phases of Carbon - and stainless steels are both covered. The associated software is available on the following websites (i) http://www.mxif.manchester.ac.uk/resources/software (ii) https://sites.google.com/site/shyamkaragadde/software/thermophysical-properties.

1. Introduction

Recent mathematical models, coupling fluid flow, heat transfer and shell solidification in the continuous casting mould, have been shown to be capable of providing (i) reliable values for the horizontal heat flux and powder consumption and (ii) insight into the mechanisms responsible for defect formation.1,2) However, these models require input data for the physical properties of both the steel and the mould powders used in the casting. Physical property determinations are both time-consuming and expensive. Thus it is unreasonable to use measured property data for all the steels being continuously cast. Consequently, in this study relations between property values for various types of steel and their chemical compositions are explored to provide property data for a wide range of steels since chemical composition data are routinely available. A similar study has been carried out for the mould powders in Part 1 of this communication.3) The predicted physical property values can then be used as input data in the mathematical model to then calculate the fluid flow, microstructure, etc.

Thus the objectives of the study were:

(i) To provide values of the physical properties of various steels calculated from their chemical compositions, and

(ii) To continually improve the methods used for property estimations.

2. Calculation of Steel Properties

The physical properties of steels in the solid state are affected by the nature of the phase formed in the steel (i.e. ferrite or austenite). Thus differences in elemental concentrations (e.g. % C) tend to have a much larger effect on the property values when they result in a phase change from ferrite to austenite or visa-versa than when there is no phase change involved. Thus it important to determine whether a steel is ferritic or austenitic, or a mixture of the two phases. In order to demonstrate these property differences, calculations have been carried out for two steels (a peritectic, low alloy steel and an austenitic, stainless steel); the compositions are given in Table 1.

Table 1. Chemical composition of two steels (mass%). LA=low-alloy and SS=stainless steel.
FeCSiMnCrNiMoTiVAlSPN
LA98.30.100.21.00.150.050.050.050.050.050.010.010.007
SS50.50.172.241.7424.420.60.240.050.050.050.0050.0150.023

2.1. Calculation of Steel Composition

The chemical composition of the steel is first normalised to 100% and then mole fractions (X) are calculated. The mean molecular weight (M) of the steel is also calculated.

2.2. Classification of Steels

A directory has been developed to interrogate the composition to determine whether:

(i) a steel was to be classified as a low alloy or a stainless steel, and

(ii) it forms austenite or ferrite or a mixture of the two phases (i.e. peritectic).

In this study a steel was classified as a stainless steel if it contained >8% Cr.

The scheme used to differentiate between the steels is that devised by Wolf4,5) and involves the calculation of the parameters, ferrite potential (FP) and carbon potential (CP) which are defined in Eqs. (1) and (2).

Wolf4,5) introduced the term “ferrite potential” as a measure of the amount of peritectic reaction occurring in the steel. Although the ferrite potential is affected by the carbon content, it is also affected by other alloying elements; some of these elements stabilise the ferrite (e.g. Cr, W, Mo, Al and Si) whilst others stabilise the austenite phase (e.g. Ni, Mn, Co, N and Cu). For low-alloy steels, the ferrite potential (FP) can be calculated by Eq. (1) where CP is the carbon potential which is determined using Eq. (2).   

Low alloy:FP=2.5( 0.5-CP ) (1)
  
CP=( wt%C ) +0.04( %Mn ) +0.1( %Ni ) +0.7( %N ) -0.14( %Si ) -0.04( %Cr ) -0.1( %Mo ) -0.24( %Ti ) (2)

The phase is considered to be ferrite when FP>1.05 and to be austenite when FP=<0.8.The peritectic reaction is at a maximum at FP=1.0 and the range, FP=(0.8 to 1.05), corresponds to the crack-sensitive range for steels. The peritectic range occurs between CP values of 0.06 and 0.17 and approximate values of fraction of ferrite (denoted fα) formed are calculated by Eq. (3). For accurate values of fractions ferrite and austenite the reader is advised to consult thermodynamic software packages.   

f α =1- f γ =( 0.17-CP ) /( 0.17-0.06 ) (3)

For stainless steels the FP is calculated from the (“Ni”/ “Cr”) ratio using Eq. (4)4,5) where “Ni” and “Cr” are defined in Eqs. (5) and (6), respectively.   

Stainless steel:FP=5.26( 0.74-Ni/Cr ) (4)
  
Ni=( wt% Ni ) +0.31( %Mn ) +22( %C ) +17.5( %N ) (5)
  
Cr=( wt% Cr ) +1.5( %Si ) +1.65( %Mo ) +2( %Nb ) 3( %Ti ) (6)

An approximate value for the fraction ferrite in stainless steels is calculated4,5) from Eq. (7).   

f α =1- f γ =5.28{ 0.74-( Ni/Cr ) } (7)

2.3. Calculation of Transition Temperatures

The thermodynamic temperature scale (K) is used throughout this paper except where certain authors have used Celsius in their equations; these were changed to K in the calculations. Routines for calculating approximate values of the transition temperatures are described below since the properties are affected by the nature of the phase formed in the steel; for accurate values of these transition temperatures the reader should refer to thermodynamic software predictions.

2.3.1. Liquidus Temperature (Tliq)

Several studies have reported equations relating liquidus temperatures as a function of the chemical composition of the steel.6,7,8,9,10,11,12) In addition, Jablonka et al.12) reported equations to calculate Tliq values for steels based on the Fe–C system for ferrite (Tliq(K)=1811−79%C) and austenite i.e. for steels with C%>0.16 (56.02.%C=(1811 K−Tliq) −1.13.10−3(1811 K−Tliq)2). The Tliq values for the steel calculated here for both low alloy and stainless steels were derived using Eqs. (8), (9), (10), (11) due to Miettinen and Howe.11) The steels are classified according to the phase formed on solidification and interactions between elements are taken into account eg. interaction between Cr and C on Tliq is covered in Eq. (8) by {0.033.%Cr. (%C)}. Miettinen and Howe cite an average uncertainty of ± 4 K from a comparison of calculated with experimental Tliq values. The equations used to calculate Tliq values for the various categories of steel11) are given in Eqs. (8), (9), (10), (11).   

Low alloy Ferritic steels: T liq ,K=1   811-76.2%C-10.35. ( %C ) 2 -11.66%Si-4.35%Si.( %C ) -5.62%Mn -0.223%Mn.( %C ) -1.95.%Cr-0.033.%Cr.( %C ) -2.2%Mo-0.845%Mo.( %C ) -3.58%Ni-0.836%Ni( %C ) -24.78%P-12.94%P( %C ) -32.8%S-17.72%S( %C ) (8)
  
Austenitic T liq ,K=1   801-60.1%C-6.14. ( %C ) 2 -11.49%Si -5.61%Si.( %C ) -4.26*%Mn+0.453*%Mn.( %C ) -2.47.%Cr+1.33.%Cr.( %C ) -4.36%Mo+0.07%Mo.( %C ) -1.97%Ni-0.589%Ni( %C ) -30.92%P-5.35%P( %C ) -33.2%S-10.09%%S( %C ) (9)

Stainless steels   

Ferritic T liq ,K=1   811-76.28%C-1.399%Cr.%C-0.0837%Ni%C -9.43%Si-0.2128%Cr.%Si+0.2391%Si.%Ni -7.55%Mn+0.1779%Mn.%Cr+0.01%Mn.%Ni -1.56%Cr-0.0236 ( %Cr ) 2 -2.87%Mo+0.044%Cr%Mo -0.0205%Ni%Mo-4.29%Ni+0.0095%Cr%Ni0.0428 ( %Ni ) 2 12%Nb+0.0789%Cr.%Nb+1.77`%Ni.%Nb -18.55%Ti+0.0538%Cr.%Ti+0.3514%Ni.%Ti -42.64%N-1.2721%Cr%N-1.0988%Ni%N +0.002084 ( %Cr ) 3 +0.000695 ( %Ni ) 3 +238.96%Ti%N (10)
  
Austenitic T liq ,K=1   801-62.62%C-0.5175%Cr.%C -0.20267%Ni%C-8.29%Si-0.4646%Cr.%Si +0.1755%Si.%Ni-4.22%Mn-0.143%Mn.%Cr -0.0377%Mn.%Ni-2.59%Cr-0.034 ( %Cr ) 2 -0.26%Mo+0.0494%Cr%Mo-0.0244%Ni%Mo -1.69%Ni+0.0546%Cr%Ni-0.0175 ( %Ni ) 2 -10.17%Nb -0.0225%Cr.%Nb+0.0781`%Ni.%Nb-16.91%Ti -0.0907%Cr.%Ti+0.1422%Ni.%Ti-21.6%N +0.1153%Cr%N-0.2366%Ni%N+0.000392 ( %Cr ) 3 -0.000356 ( %Ni ) 3 +154.98%Ti%N (11)

2.3.2. Solidus Temperature (Tsol)

Jablonka et al.12) reported the following relations for the solidus temperature (Tsol) of ferritic (Tsol(K)= 1811−460% C) and austenitic steels (Tsol(K)= 1811−185% C) based on the Fe–C system. Since the alloying elements have a significant effect on both Tsol and Tliq (see Eqs. (8), (9), (10), (11)) it is not possible to use these relations to estimate Tsol for the wide compositional range covered in this study. Consequently, in the calculation of physical properties it has been assumed that fusion takes place abruptly at Tliq.

2.3.3. Solid State Transition Temperatures (T)

There are several solid state transitions in the Fe–C system. Jablonka et al.12) reported the following relations to calculate transition temperatures for the Fe–C system; at the present time no corrections have been made for the effect of alloying elements on the various transition temperatures. Eqs. (12), (13) and (14) have been adopted to calculate the effects of C content on the transition temperatures for (γδ) and (α+γγ), respectively.   

Transition( γδ ) T γδ ( K ) =1   665+624.4%C (12)
  
C<0.76% T α+γγ ( K ) = A 3 =1   171-584%C+593% C 2 -190% C 3 (13)
  
C>0.76% T α+γFe3Cγ ( K ) = A cm =865+24.1( %C ) +159.2 ( %C ) 2 (14)
  
 Transition ( αγ ) T αγ ( K ) = A 1 =990   K (15)

2.4. Heat Capacity (Cp) and Enthalpy (HT-H298)

The Cp and enthalpy are calculated using partial molar values (shown in Eq. (16) where P is the property being studied, Xi is the mole fraction of components of steel such as Fe, Cr, Mn, Ni etc.).   

P steel = i=1 N ( X i P i ) (16)

The temperature dependence of heat capacity (Cp) for the solid phase is usually expressed in the form of Eq. (17) where a*, b* and c* are constants and the enthalpy (HT-H298) is given by Eq. (18).   

C p =a*T+b*T-c*/ T 2 (17)
  
( H T - H 298 ) = 298 T C p dT =a*( T-298 ) +0.5b*( T 2 - 298 2 ) +( c*/ T ) -( c*/298 ) (18)

However, Cp-T relations for Fe alloys contain a sharp peak in Cp around 1020 K due to the magnetic transition but the magnitude of the enthalpy associated with this transition tends to decrease with increasing alloy addition in steels, This makes it difficult to model. Consequently, the coefficients for a*, b* and c* for Fe were obtained in the form of Eq. (18) from a best fit of enthalpy values for Fe. Thus the resultant Cp values do not show the enhanced peak in Cp associated with the magnetic transition but this has little effect on the enthalpy values.

The enthalpy of transition (ΔHtr ) of ferritic steels for the various transitions are given below   

Δ H tr αγ =18.1   kJ k g -1 . (19)
  
Δ H tr γδ =14.8 kJ k g -1 . (20)
  
Δ H fus = T liq Δ S fus kJ k g -1 (21)

The values of Cp and (HT-H298) are given in JK−1kg−1 and kJ kg−1, respectively; values were converted from mol−1 to kg−1 by multiplying by (1000/M) where M is the average molecular weight of the steel (M = ΣXi Mi). The calculated values of Cp and (HT-H298) for low alloy and austenitic stainless steels are shown in Figs. 1(a) and 1(b), respectively. It can be seen that Cp and enthalpy are little affected by composition.

Fig. 1.

Calculated values for (a) Cp and (b) Enthalpy (HT-H298) as functions of temperature for a peritectic low alloy steel (Δ) and an austenitic stainless steel (■); Transition temperatures shown as vertical lines; dashed=low alloy; dotted= stainless steel.

2.5. Density (ρ)

Low alloy steels

The molar volume (V) of the steel at 298 K was calculated using partial molar volumes of the various steel components (Eq. (22)).   

V 298 steel = i=1 N ( X i V i ) (22)

The density (ρ) can be calculated from the molar volume by using Eq. (23) where M is the molecular weight of the steel   

ρ steel =M/ V steel (23)

Ferritic steels

The density of steel is also affected by the carbon content, Consequently, a correction was applied for the effect of carbon. The relations reported by Jablonka et al.,12) shown in Eq. (24), were used where g=2.62×10−2 for the α and δ phases and 1.46×10−2 for the γ phase.   

ρ T, steel = ρ T,Fe (1-g.   C%) (24)

Jablonka et al.12) expressed the temperature dependence of ρT,Fe by Eq. (25) where T*={T−273}. The temperature dependence of ρT, steel was obtained by substituting ρ298, steel for 7876 in Eq. (25).   

α andδ-phases: ρ T,Fe (kg m -3 )=7   876-0.297T*- 5.62.10 -5 T * 2 (25)

The transition (αγ) is accompanied by an increase in density (Δραγ = +80 kgm−3.12)

Austenitic steels

Austenite has a higher density than ferrite. Jablonka et al.12) reported the following relation for austenite, where T*={T(in K)−273}. This equation was extrapolated to lower temperatures to calculate density values for the range, 298 K to Tαγ. However, corrections were made for the effect of C content and the alloying content on ρ298 γ,Fe.   

γ-phase: ρ T γ ,Fe (kg m -3 )=8   100-0.506T* (26)

The transition of austenite to δ –phase is accompanied by a decrease in density, Δρ=−42 kg m−3.12) The transition of δ –phase to liquid is also accompanied by a decrease in density, Δρfus=−250 kgm−3.12)

Liquid steel

The density of the liquid at Tliq (denoted by the superscript m) was calculated using Eq. (27). Jablonka et al.12) reported Eq. (28) for the effect of carbon on ρT,Fe for the liquid phase. Thus the density of the liquid was calculated by combining ρmliq with the dependencies of density on the C content and temperature12)]; as shown in Eq. (28).   

ρ m liq (kg m -3 )= ρ m sol -250 (27)
  
ρ T liq (kg m -3 )= ρ m liq -0.835( T- T liq ) (28)

It should also be noted that ρmliq can be calculated from the relation (Vm steel=Σ(X.Vm)Fe+(X.Vm)Cr+(XVm)Ni+(XVm)Mn+ (XVm)Si+ . . .) and (ρmliq= M/Vm) where the superscript, m indicates at Tliq. Values for (dV/dT) can be also calculated on a partial molar basis; this was used in the calculation of the temperature dependency of the liquid alloy. The effect of C content can be introduced by multiplying ρmliq by (1−0.01%C).12) These values were found to be in good agreement with the adopted values.

The calculated density values for low alloy and austenitic stainless steels are shown in Fig. 2.

Fig. 2.

Calculated values of the density of low-alloy steel (Δ) and austenitic, stainless steels (■) as a function of temperature; Transition temperatures are shown as vertical lines; dotted=stainless and dashed=low-alloy steels.

Fig. 3.

Calculated values of (a) Thermal conductivity and (b) Thermal diffusivity as function of temperature for low alloy (Δ) and austenitic stainless (■) steels.

Stainless steels

The densities of stainless steels were calculated using the empirical relation shown in Eq. (29).13) It should be noted that although the ferrite to austenite transition is accompanied by an increase in density (Δρ) the values of ρ298 for the two phases are similar because Δρ is largely offset by the higher thermal expansion coefficient for austenite.   

ρ 298 steel ( kg m -3 ) =79.6( %Fe ) +78.3( %Cr ) +85.4( %Ni ) +76.9( %Mn ) +60.2( %Mo ) +47.1( %Si ) (29)

The linear thermal expansion coefficient (α) for austenitic stainless steels is given in Eq. (30).13)   

10 6 α( K -1 )=15.8+0.6× 10 -2 ( T-298 ) (30)

The expansion coefficient for ferrite is significantly lower than that of austenite for temperatures between 298 and 1000 K as can be seen from Eq. (31)13) but increases abruptly for temperatures above 1000 K with the formation of austenite (Eq. (32))   

( 298-1   000   K )                  10 6 α( K -1 ) =10.2+0.6× 10 -2 ( T-298 ) (31)
  
(1   000 K- T liq ) 10 6 α( K -1 ) =14.2+2.4× 10 -2 ( T-1   000 ) (32)

Densities at various temperatures were calculated using Eq. (33) where Tref is 298 K or 1000 K in the case of ferrite at temperatures above 1000 K.   

ρ T steel ( kg m -3 ) = ρ Tref steel ( 1-  α T ( T- T ref ) ) 3 (33)

Liquid phase

The density of liquid stainless steel can be calculated using Eq. (34).   

ρ m liq ( kg m -3 ) =69.4( %Fe ) +66.3( %Cr ) +71.4( %Ni ) +57.2( %Mn ) +51.5( %Mo ) +49.3( %Si ) (34)

However, it can also be calculated by using the calculated value of ρm sol and the following relation, ρm liq=ρm sol/1.04 (derived from Fe14)). The calculated density values obtained with Eqs. (29) and (34) are subject to uncertainties of ca. ± 2%. Consequently, in order to avoid the situation where ρm liq>ρm sol, a mean of the ρm liq values obtained by the two methods was used to avoid this anomaly.

2.6. Thermal Conductivity (k) and Diffusivity (a) and Electrical Resistivity (R)

2.6.1. Thermal Conductivity (k)

Experimental measurements of the thermal conductivity (k) of molten metals are prone to errors arising from contributions from convection; these errors are accentuated by the low viscosities of metals. Values for the solid state are affected by the microstructure, which is, in turn, affected by the thermal and mechanical histories of the sample. Consequently, values cited here refer to annealed alloys i.e. with minimum resistivity. At high temperatures electronic conduction in metals is dominant15) Consequently, it has been proposed15) that the most reliable method of estimating thermal conductivities of liquid metals is via the Wiedemann-Franz-Lorenz (WFL) rule (Eq. (35) where R is the electrical resistivity).   

k el T =2.45 10 -8 ( T/R ) (35)

Low alloy steels

The WFL rule does not apply to solids at lower temperatures where there are two contributions, namely,

(i) electronic conduction, kelT (which can be calculated by combining Eqs. (36)16) and (35).

(ii) lattice thermal conductivity kTlat which is calculated by an empirical relationship (Eq. (37))16)   

R 298 ( 10 -8 Ω  m -1 )=100{ 0.11%Fe+1.92%Cr+3.53%Mn +0.97%Nb+2.1%Ni+- 9.9%Si+3.34%Ti+7.9%Al +2%Co+8.28%C+2%S+1%Ca+0.11%Nb 0.6%W +4.5%P } (36)
  
k 298 lat ( W m -1 K -1 ) =8.1%Fe+0.35%Cr+0.11%Mn+0.39%Nb -0.28%Ni+-96.8%Si-1%Ti   100%Al-0.28%Co+246%C +1%S+1%Ca+0.81%Nb+4.6%W+1%P  (37)

The effective conductivity (keff) at 298 K. can then be calculated using Eq. (38).   

k eff T = k el T + k lat T (38)

The thermal conductivity of ferrite (kα) at 298 K is significantly higher than that of austenite (kγ). However, α- ferrite transforms to γ- austenite around. 1180 K in most low alloy steels and kγ has a value of 25.4 Wm−1K−1 at this temperature.16) The k-T relations for temperatures between 298 K and Tαγ can be expressed in the form of Eq. (39) for both austenitic and ferritic steels. Note that values for the temperature coefficient (dk/dT) in this range will be negative for ferritic steels and positive for austenitic low-alloy steels   

( 298-  T α γ )             k T ( W m -1 K -1 ) = k 298 eff +{ ( k Tαγ - k 298 ) ( T-298 ) /( T αγ -298 ) } (39)

The keff values for austenite in the temperature range (Ttrαγ to Ttrγδ ) were calculated using the empirical relation shown in Eq. (40); the values for the δ-phase were derived by extrapolating the relation for the γ- phase to Tliq. Note the calculated values are independent of the alloy composition in this range; any uncertainties due to composition were found to be within experimental uncertainties   

:( T tr αγ - T liq γδ ) k eff T ( W m -1 K -1 ) =25.4+0.013( T- T tr αγ ) (40)

Values of k for the liquid phase were obtained by assuming that (i) km(l)/km(s)=1.07=Rm(s)/Rm(l) (ii) the WFL Rule applied and (iii) the temperature coefficient (dk/dT)=0.015 Wm−1K−2.14) The calculated values of keffT for low alloy and austenitic stainless steels are shown in Fig. 4(a) and the equivalent thermal diffusivity values are given in Fig. 4(b). It can be seen from Fig. 4(a) that (i) kα> kγ for temperatures between 298 K and 1000 K (ii) the temperature coefficient (dk/dT) is negative in ferrite but is positive in both austenite and the liquid (iii) there is a small decrease in thermal conductivity and diffusivity on fusion.

Fig. 4.

The electrical resistivity (in 10−8 Ω m) as a function of temperature for a low alloy steel (●) and an austenitic stainless steel (■).

Stainless steels

Values of keff 298 can be calculated by deriving values for kel using a combination of Eqs. (35), (36), (37), (38). The values of kT for austenitic steels are little affected by changes in chemical composition, so values of kT above 1100 K can be calculated using Eq. (37). Alternatively, keff can be calculated for austenitic stainless steels by using Eq. (41) for temperatures between 298 K and Tliq.13)   

Austenitic: k eff T ( W m -1 K -1 ) =9.2+0.0175 T-2× 10 -6 T 2 (41)

Values of the thermal conductivities of ferritic stainless steels are also little affected by differences in concentrations of the alloying elements. However, ferritic stainless steels undergo a magnetic transition at ca. 1000 K so two equations are needed to describe the temperature dependence of k. Thus Eqs. (42) and (43), respectively, describe k values in the temperature ranges (298–1100 K) and (1100−Tliq K). Uncertainties of ± 5% and ± 10% in the k values were reported for austenitic and ferritic stainless steels, respectively.   

Ferritic( 298-1   100   K ) : k eff T ( W m -1 K -1 ) =23.5+0.016( T-298 ) (42)
  
Ferritic( 1   100- T liq ) : k eff T ( W m -1 K -1 ) =25.4+0.013( T-1   100 ) (43)

2.6.2. Thermal Diffusivity (a)

The thermal diffusivity values (aeff) are derived from the calculated values of Cp, ρ and k and the relation a=(k/Cp ρ). The calculated values for the thermal diffusivity are given in Fig. 4(b). It should be noted that experimental values of thermal diffusivity show a sharp valley around 1000 K due to the peak in Cp caused by the magnetic transition. The calculated Cp values do not contain this Cp peak; thus there is no pronounced valley in the calculated thermal diffusivity values.

Electrical resistivity (R) and conductivity (κ)

There is little demand for electrical resistivity values in the modelling of the continuous casting process. However, electrical resistivity values are calculated in the derivation of thermal conductivity values and electrical resistivity data for steels are needed for other steelmaking processes. Electrical resistivity values are cited in units of 10−8 Ω m (i.e. 70×10−8 Ω m=700 nΩ m).

Low alloy steels

The electrical resistivity of a steel at 298 K is calculated using Eqs. (36), (37), (38). Austenitic steels have a higher resistivity than ferritic steels. At high temperatures the values of electrical resistivity are little affected by chemical composition so constant values for specific temperatures can be adopted with little loss in accuracy. Values for the liquid can be calculated with an uncertainty of ±10% using Eq. (44) where Rim is the electrical resistivity at Tliq for the alloying constituents (e.g. Fe, Mn, Si etc.). However, Rm(l) values can be calculated from the calculated Rm(s) value and the relation Rm(s)/Rm(l)=1.07 (i.e. equal to that for pure Fe14)); these latter values are considered to be more reliable and have been adopted.   

R m steel ( 10 -8 Ω  m -1 ) = i=1 N ( X i R i m ) (44)

Stainless steels

The electrical resistivities are calculated for austenitic and ferritic steels using Eqs. (45) and (46), respectively.13)   

Austenitic: R( 10 -8 Ω m -1 )=51.9+8.6× 10 -2 T-2.35× 10 -5 T 2 (45)
  
Ferritic: R( 10 -8 Ω m -1 )=23.8+12.9× 10 -2 T-4.17× 10 -5 T 2 (46)

2.7. Viscosity (η)

Low alloy steels

The experimental uncertainties for measurements on pure iron are about ± 0.5 m Pas or about ± 10%. Steels are less pure and contain non-metallic inclusions which increase the viscosity; thus experimental uncertainties associated with viscosity measurements on steels are probably ca. ± 20%.With this level of experimental uncertainty it is difficult to determine the effect of alloy concentrations (e.g. Si%, C%) on the viscosity. For this reason, the viscosity of steel was taken to be that of pure iron (given by Eq. (47).14) However, this relation was amended (see Eq. (48)) to allow for the effects of (i) alloying elements and (ii) non-metallic inclusions on the viscosity. The calculated results are shown in Fig. 5.   

  η Fe ( mPas ) =0.1896 exp ( 6   206.5/T ) (47)
  
η steel =1.1    η Fe (48)
Fig. 5.

Calculated values for the viscosities of a low alloy steel (Δ) and stainless steel (■) as a function of temperature.

Stainless steels

The following relation was derived from measurements of the viscosity of a stainless steel using oscillation viscometry.13,17)   

ln η ( mPas ) =-2.396+( 7   950/T ) (49)

It can be seen from Fig. 5 that the viscosity values calculated for the two steels using the different methods lie within ± 10% of one another.

2.8. Surface Tension (γm) and Interfacial Tension (γms)

Surface tension of steel (γm)

The surface tension of steel (γ) and its temperature dependence (dγ/dT) are both very sensitive to the concentrations (ppm levels) of soluble O and S contents of the steel. However, most steels contain Al contents of > 50 ppm, and these Al levels are sufficient to hold the soluble O content to ca. 5 ppm. Most alloying elements have little effect on the soluble S levels (the exceptions being Ca and the rare-earths). For this reason, the soluble S content approximates to the total S content. Thus, the surface tensions of steels are sensitive to the S content of the steel. Consequently, Eq. (50) due to Su18) has been used to calculate the surface tensions of liquid iron. The effects of alloying elements can be taken into account by (i) using density data for the elements and the method of Butler19) or (ii) by using a partial molar approach for the alloy (excluding the S,O i.e. (γalloy= i=1 N ( X i γ i ) .20) However, the densities and surface tensions of the principal alloying elements tend to be similar to those of iron. Consequently, the simplified relation shown in Eq. (50) can be used for steels without much loss of accuracy.   

γ(mN m -1 )=1   880-0.41(T-1   803   (K)) -0.009T   ln( 1+ e 19   411 T -4.6849 ×(wt%S) ) (50)

The calculated surface tension values for low alloy and stainless steels are shown in Fig. 6. It can be seen that:

Fig. 6.

Calculated values of surface tension as a function of temperature for (a) a low alloy containing 100 ppm S (Δ) and austenitic stainless steel containing 50 ppm S (■) and (b) for steels with identical composition but containing different S contents: 10 ppm=▲; 20 ppm=O; 50 ppm=●: and 100 ppm=X.

(i) The surface tension decreases with increasing S content.

(ii) The temperature coefficient (dγ/dT) becomes increasingly positive with increasing c S content; (dγ/dT) is negative for 10 and 20 ppm S but is slightly positive for 50 ppm S but becomes progressively more positive with increasing S contents.

Stainless steels

A similar relation (Eq. (51)) has been reported which allows the surface tension of stainless steels to be calculated.13,21) The temperature coefficient (dγ/dT) can be calculated using Eq. (52).   

γ( mN m -1 ) =1   840-0.4( T-1   823 K ) -0.056 ln {1+ e [( 28   798/T ) -8.5647] ( 0.68   wt%S ) (51)
  
(dγ/dT)=-0.4-0.056{ [ 28   798( 1-B ) /BT ]+ln B } (52)
Where B=1+e [(28798/T)−8.5647] (0.68 wt%S)

Interfacial energy or surface tension (γmsl)

The interfacial energy or interfacial tension (γmsl) is calculated from Eq. (53) where γsl =slag surface tension, γm can be derived from Eq. (34) and φ is an interaction coefficient, the magnitude of which was calculated in Part 1 (Eqs. (44), (45), (46)). The value of γmsl is largely determined by γm since γm≈ (4 to 5) γsl and the value of γm, in turn, is largely determined by the S content of the steel.   

γ msl = γ m + γ sl -2φ ( γ m γ sl ) 0.5 . (53)

2.9. Emissivity (ε)

Low alloy steel

Emissivity is also a surface property. In practice, it is affected by the nature of the oxides formed in the surface layers with ε values varying from (0.1–0.3) for high vacuum conditions to 0.7 for oxidised conditions. Consequently, the emissivity values are affected by the partial pressure of oxygen (pO2). Since the atmosphere in the continuous casting mould is very reducing, the emissivity values reported here refer to low pO2 values. The values of the total normal emissivity (εTN) of pure Fe14) were adopted for the steel; the following relations were used: The calculated values for a low alloy steel are shown in Fig. 7.   

298   K- T αγ :   α-Fe: ε TN =0.075+0.000113( T-298 ) (54)
  
T αγ T γδ :γ-Fe( austenite ) ε TN =0.075+0.000208( T- T αγ ) (55)
  
T γδ T δliq :δ-Fe ε TN =0.235+0.000136( T- T γδ ) (56)
  
T> T liq :liquid ε TN =0.275 (57)
Fig. 7.

Calculated values of the total normal emissivity as a function of temperature for a low alloy steel (Δ), an austenitic stainless steel (■)22,23) and liquid steel (●).

Stainless steel

Total normal emissivity values22,23) reported for an austenitic, stainless steel are higher than those cited for low alloy steels; these data can be represented by Eq. (58). It can be seen that the values are in good agreement with the values for austenite in the low alloy steel.   

ε TN =0.122+0.0001( T-298 ) (58)

It can be seen from Fig. 7 that εTN values for the austenitic stainless steel are higher than those for the low alloy, peritectic steel.

Acknowledgements

This work was made possible by the facilities and support provided by the EU (RFSR-PR-10005 DDT), the Research Complex at Harwell, and the EPSRC (EP/I02249X/1).

Symbols, Abbreviations, Units

a= Thermal diffusivity (m2s−1)

Cp= Heat capacity (JK−1kg−1)

H= Enthalpy (J kg−1)

k= Thermal conductivity (Wm−1K−1)

R= Electrical resistivity (10−8Ω m−1)

T= Temperature (K)

X= mole fraction

α= Thermal expansion coefficient (K−1)

γ= Surface tension (mNm−1)

εTN= Total normal emissivity

η= Viscosity (Pas)

κ= Electrical conductivity (10−8 Ωm=10nΩm)

ρ= Density (kgm−3)

Appendix

User to Insert

a) Composition: %Fe; %Cr; %Ni; %Mo; %Ni; %Si; %Ti; %V; %Al; %Co; %C; %S; %O; %Ca; % Nb; %W. Note 5 ppm S=0.0005%.

b) Tliq and other Ttr values: Default values will be calculated (±5 K) if no insert is made by the user.

Output

All calculated data as a function of temperature given on “Collected” Worksheet.

References
 
© 2016 by The Iron and Steel Institute of Japan
feedback
Top