UDK 621.74.047:519.

68 ISSN 1580-2949
Original scientific article/Izvirni znanstveni ~lanek MTAEC9, 47(6)815(2013)

M. VELICKA et al.: RESEARCH OF THERMAL PROCESSES FOR THE CONTINUOUS CASTING OF STEEL

RESEARCH OF THERMAL PROCESSES FOR THE

CONTINUOUS CASTING OF STEEL
RAZISKAVE TERMI^NIH DOGAJANJ PRI KONTINUIRNEM
ULIVANJU JEKLA

Marek Velicka, David Dittel, Rene Pyszko, Miroslav Prihoda, Miroslav Vaculik,
Pavel Fojtik, Jiri Burda
VSB – Technical University of Ostrava, Faculty of Metallurgy and Materials Engineering, 17. listopadu 15, 708 33 Ostrava, Czech Republic
marek.velicka@vsb.cz

Prejem rokopisa – received: 2013-04-03; sprejem za objavo – accepted for publication: 2013-04-08

The article deals with the determination of the basic indicators of heat transfer in the continuous casting of steel, which can be
described as an unsteady process with complicated boundary conditions for the solution. An analytical solution of this problem
is practically impossible and, therefore, mathematical modelling is applied with a certain simplification of the real conditions
and with a description of those criteria that influence the most the process of solidification and cooling. Using a simulation
program and the knowledge of input parameters, it was possible to predict the distribution of the thermal field of a continuously
cast blank in the course of its casting. Simulations also allowed us to deal with the issues of the inner structure, surface quality,
mechanical properties of a continuously cast blank, metallurgical length, change in the thickness of a strand shell and over-
heating of steel. Some results obtained with numerical simulations are documented for concrete examples.
Keywords: continuous casting of steel, modelling, heat transfer, shell, mould

^lanek obravnava dolo~anje osnovnih pokazateljev prenosa toplote pri kontinuirnem ulivanju jekla, ki se lahko opi{ejo kot
nestabilen proces s kompliciranimi robnimi pogoji za re{itev. Analitska re{itev tega problema je prakti~no nemogo~a, zato je
bilo uporabljeno matemati~no modeliranje z nekaterimi poenostavitvami realnih pogojev in z opisom tistih meril, ki najbolj
vplivajo na proces strjevanja in ohlajanja. S programom za simulacijo in poznanjem vhodnih parametrov je bilo mogo~e pred-
videti razporeditev temperaturnega polja kontinuirno ulite gredice med njenim ulivanjem. Uporaba simulacije je omogo~ila opis
notranje strukture, kvalitete povr{ine, mehanskih lastnosti kontinuirno ulite gredice, metalur{ke dol`ine, spreminjanja debeline
strjene skorje in pregretja jekla. Nekateri rezultati, dobljeni z numeri~no simulacijo, so prikazani za konkretne primere.
Klju~ne besede: kontinuirno ulivanje jekla, modeliranje, prenos toplote, skorja, kokila

calculations of stress conditions, and also to predict

1 INTRODUCTION segregation during the continuous casting of steel. The
Mathematical precise description of thermal pro- finite-element method or finite-difference method are
cesses during the continuous casting of steel is very used most frequently as the calculation algorithm in this
difficult, since the process of cooling and solidification software.
of a continuously cast blank is influenced by many This paper deals with the solution of thermal pro-
factors. cesses during the continuous casting of round steel
For this reason it is necessary to find, with the help of blanks with a diameter of 410 mm, using the simulation
mathematical and physical modelling, the criteria that software based on the explicit finite-element method.
have the biggest influence on the solidification and cool-
ing of a continuously cast blank. Understanding the 2 THERMAL FIELD OF A CONTINUOUSLY
thermal processes taking place during the continuous CAST BLANK
casting of steel is of crucial importance because it ena-
bles a prediction of a formation of defects, an enhance- The kinetics of an unsteady thermal field is described
ment of the thermal processes during continuous casting, with the Fourier partial differential equation.2 If we want
the optimum locations for cooling nozzles or a minimisa- to describe the thermal field of a moving, continuously
tion of breakout risks, etc.1 cast blank, it is necessary to take into account, in the
It is evident, that it is impossible to optimise the pro- transverse direction, also the rate of casting, by trans-
cess of continuous casting of steel only by modelling. A forming the classical Fourier equation to the Fourier-
very close interaction with the results of experimental Kirchhoff equation:
measurements is always necessary since these results ∂t qv
= a ⋅ ∇ 2t + (1)
introduce, into the system, characteristic features of a ∂t cp ⋅ r
concrete continuous-casting machine. At present, it is
possible to perform, with the use of sophisticated where a/(m2/s) is the temperature-conductivity coeffi-
software, not only the thermal calculations but also the cient, ∇ 2 /(m–2) is the Laplace's operator, cp/(J/kg K) is

Materiali in tehnologije / Materials and technology 47 (2013) 6, 815–818 815

M. VELICKA et al.: RESEARCH OF THERMAL PROCESSES FOR THE CONTINUOUS CASTING OF STEEL

the specific heat capacity, r/(kg/m3) is the density and

qv/(W/m3) is the inner heat source.
This equation can only be solved with the conditions
of monovalency of the solution, which, in the classical
concept, are divided into geometric, physical, surface
and initial conditions. Physical conditions are generally
created with the chemical composition of the cast steel,
geometric conditions are created with the shape of the
cast products, and the initial condition is closely linked
to the temperature of the liquid steel in a tundish.
Numerous researches concerning the continuous casting
of steel3–6 combined with modelling physical models use
the Neumann surface condition for the mould and the
Fourier condition for the secondary and tertiary cooling.

3 SIMULATION MODELLING
For numerical modelling, the ProCAST program tool
was used, which is based on the implicit finite element Figure 2: Calculation mesh of a blank
method. Alternatively, an original program code Tefis Slika 2: Ra~unska mre`a gredice
based on the explicit finite difference method, specially
designed for numerical modelling of the continuously An application of program ProCAST consists of five
cast blank temperature field, was used. The program modules, which are functionally mutually inter-
Tefis enables real-time simulations too. The second connected. The first module is used for defining the body
method requires compliance with a numeric stability geometry and for selecting the fineness and the shape of
condition, which describes mutual dependence between the calculation mesh. It is then necessary to determine
fineness of the calculation mesh and calculation time the physical, initial and surface conditions that may be
step. modified in line with the requirements of the next
module. After correctly entering the initial calculation
parameters, the main calculation module of the software
launches the program.
The program allows a graphical visualisation of the
calculated results and an export of the values, graphs and
images for further processing. The program can launch a
calculation in two modes – the thermal and flow modes.
The results of the thermal mode give the thermal fields
of a round, continuously cast blank without the velocity
vectors of liquid steel. If it is also required to visualise
the velocity vectors, then it is necessary to activate the
flow mode (Figure 1).
Distribution of the temperatures during the solidifi-
cation or cooling is calculated in the nodal points of the
whole volume of a continuously cast blank (Figure 2).
The program also comprises a vast database of the
information about cast steel grades, including their ther-
mophysical properties (density, specific heat capacity,
heat conductivity, viscosity, etc.).

4 RESULTS OF NUMERICAL SIMULATIONS

Within the research two steel grades of different
chemical compositions used for the production of con-
tinuously cast blanks with a diameter of 410 mm were
subjected to a simulation (Table 1). The following
influences were analysed: the influence of the chemical
Figure 1: Velocity vectors of a blank composition, the casting rate, overheating of the steel
Slika 1: Vektor hitrosti v gredici above the liquidus temperature, the influence of the mol-

816 Materiali in tehnologije / Materials and technology 47 (2013) 6, 815–818

 M. VELICKA et al.: RESEARCH OF THERMAL PROCESSES FOR THE CONTINUOUS CASTING OF STEEL

ten-steel level in the mould on the resulting technolo-

gical parameters, generally represented with the metal-
lurgical length, the length of the liquid phase, the surface
temperatures of a continuously cast blank at the outlets
of the primary and secondary cooling zones, and the
thickness of the strand shell.
Table 1: Chemical composition of the steel in mass fractions, w/%
Tabela 1: Kemijska sestava jekla v masnih dele`ih, w/%

Brand of Chemical composition, w/%

the steel C Mn Si P S
Steel A 0.168 1.350 0.381 0.013 0.007
Steel B 0.913 0.342 0.246 0.009 0.004
Figure 4: Increase in the thickness of the casting shell depending on
The chemical compositions, especially the carbon the casting speed in the case of steel A
contents, have a principal influence on the heat removal Slika 4: Pove~anje debeline skorje, v odvisnosti od hitrosti litja v pri-
in all the zones of the continuous-casting machine. In the meru jekla A
case of low-carbon steels the biggest shrinkage of the
strand shell occurs as a result of a peritectic reaction as bolic law may be used for the solidification of an ingot to
well as a considerable deceleration of its growth, which a distance of approximately 75 % of the diameter of a
is manifested by a reduction of the heat-flow density continuously cast blank from the ingot mould surface.
through the mould wall to its minimum. Another result When this critical value is exceeded, the differences are
of this reaction is an increased occurrence of surface more pronounced, as it is known that the solidification of
defects of a continuously cast blank, particularly of round, continuously cast blanks accelerates towards their
surface cracks. centre, in comparison to the ingot mould.
The casting rate is related to the dimensions of a con- The thickness of the shell can determine the results of
tinuously cast blank, the type of steel and the type of the logarithmic equation for steel A:
mould. Higher values of the casting rate cause a shorter x = (−1.490 ⋅10 −2 ⋅ ln v z +1683
. ⋅10 −1 ) ⋅
stay of a continuously cast blank in the mould, which (2)
( −6. 378 ⋅10 −3 ⋅ln v z + 5. 874 ⋅10 −1 )
increases the surface temperatures and, simultaneously, ⋅t
also the local values of the thermal-flow density.
And for steel B, the following applies:
It may be stated on the basis of the performed simu-
lations that, with respect to investigating the influence of x = (−13690
. ⋅10 −2 ⋅ ln v z +1869
. ⋅10 −1 ) ⋅
(3)
the casting rate on the metallurgical length and the length ( −9. 369 ⋅10 −3 ⋅ln v z + 5. 456 ⋅10 −1 )
⋅t
of the liquid phase, their very strong linear dependence is
evident and is best characterised with a high value of the where vz/(m/s) is the specific casting rate and t/s is the
slope in the regression equations (Figure 3). casting time.
The growth of the strand shell with the changing Finally, an analysis of the influence of overheating of
casting rate (Figure 4) confirms the previous scientific steel on the investigated technological parameters was
research works,7,8 which document the fact that the para- made. Overheating of steel can be defined as the diffe-
rence between the casting temperature and the liquidus

Figure 3: Influence of the metallurgical length depending on the cast-

ing speed in the case of steel A Figure 5: Blank-temperature change depending on the overheating of
Slika 3: Vpliv metalur{ke dol`ine v odvisnosti od hitrosti litja v pri- steel B
meru jekla A Slika 5: Sprememba temperature gredice glede na pregrevanje jekla B

Materiali in tehnologije / Materials and technology 47 (2013) 6, 815–818 817

M. VELICKA et al.: RESEARCH OF THERMAL PROCESSES FOR THE CONTINUOUS CASTING OF STEEL

temperature. This temperature difference should be, from steel. The casting rate, which influences practically all
the viewpoint of the operation, as small as possible, since the casting parameters, appeared to be most important
it reduces the thermal stress and the energy intensity of quantity.
the continuous casting of steel, as well as the scrap factor
and the formation of surface cracks (Figure 5). Acknowledgements
The results of the simulations also show that the
greater is the overheating of steel, the later is the solidi- The work was realised within the frame of the grant
fication of a continuously cast blank and the greater is project No. SP 2013/53, under the financial support of
the distance of this solidification from the molten-steel The Ministry of Education, Youth and Sports.
level. The overheating of steel in the interval from 0 °C to
45 °C increases the metallurgical length, on average, by 6 REFERENCES
0.5 m, and a similar trend is characteristic of all steel
1
grades. Even a greater influence of this parameter on the Y. Xie, H. Yu, X. Ruan, B. Wang, Y. Ma, Mathematical modelling of
increase in the length of the liquid phase was confirmed. mould temperature field during continuous casting of steel, Journal
of Materials Processing Technology, 199 (2008) 1–3, 49–55
An increase in the overheating temperature by 40 °C 2
I. L. Ferreira, J. E. Spinelli, J. C. Pires, A. Garcia, The effect of melt
above the liquidus temperature caused an elongation of temperature profile on the transient metal/mold heat transfer coeffi-
the metallurgical length, on average by 3 %; however, cient during solidification, Materials Science and Engineering A, 408
the length of the liquid phase was increased even by (2005) 1–2, 317–325
3
20 %, while the surface temperature of a continuously J. X. Fu, J. S. Li, H. Zhang, X. Z. Zhang, Industrial Research on
cast blank at the output from the primary and secondary Broadening of Slab in Continuous Casting Process, Journal of Iron
and Steel Research International, 17 (2010) 8, 20–24
cooling zones increased by up to 3 %. 4
M. Janik, H. Dyja, Modelling of three-dimensional temperature field
inside the mould during continuous casting of steel, Journal of Mate-
rials Processing Technology, 157 (2004) 1, 177–182
5 CONCLUSION 5
S. Sanyal, J. K. Saha, S. Chandra, A Mathematical Model to Opti-
mise Aluminium Wire Injection in Steel Melts, Steel Research Inter-
The optimum configuration of the thermal mode of national, 80 (2009) 8, 559–567
the continuous casting machine substantially influences 6
J. Stetina, F. Kavicka, The influence of the chemical composition of
the quality of the cast products. The processes partici- steels on the numerical simulation of a continuously cast slab, Mater.
pating in the cooling and solidification of a continuously Tehnol., 45 (2011) 4, 363–367
7
cast blank show their physical characters in the compli- K. Stransky, J. Dobrovska, F. Kavicka, V. Gontarev, B. Sekanina, J.
Stetina, Two numerical models of the solidification structure of
cated transfer phenomena of the heat and mass transfer. massive ductile cast-iron casting, Mater. Tehnol., 44 (2010) 2, 93–98
In order to allow a general solvability of these tasks, 8
J. X. Fu, W. S. Hwang, J. S. Li, S. F. Yang, Z. Hui, Effect of Casting
computer simulations are used making it possible to defi- Speed on Slab Broadening in Continuous Casting, Steel Research
ne the optimum parameters of the continuous casting of International, 82 (2011) 11, 1266–1272

818 Materiali in tehnologije / Materials and technology 47 (2013) 6, 815–818