1. Introduction
In continuous casting, molten steel is poured from a tundish into a water-cooled mould and a partially solidified billet or slab is withdrawn from the bottom of the mould. The billet or slab is then cooled by a water spray (this is the secondary cooling process) so that the solidified billet or slab is produced constantly and continuously. Continuous casting is a bridge between steelmaking and rolling. Secondary cooling is an essential part of continuous casting and can strongly influence the quality of billets or slabs. The cooling must be controlled relative to the casting speed to avoid the formation of internal and surface cracks.
Steel solidification behavior is influenced by heat transportation under the specific cooling conditions. In the secondary cooling zone, the heat transfer behavior of the billet or slab surface is directly linked to the characteristics of the spray. These can be manipulated to control the solidification process, and in turn the billet or slab quality and the casting productivity. In an attempt to meet high billet or slab quality requirements, heat transfer in the mould and the secondary cooling have received much research attention. Due to requirements for high quality, the secondary cooling technology of continuous casting has been developed with specially designed nozzles, finer waterway control, and a more effective water injection strategy.
In recent years, numerical simulation has been widely applied to design and optimize the secondary cooling process. An accurate simulation, however, is strongly dependent on the level of understanding of the physical mechanism of the casting process, the high-temperature properties of the material involved, and thermal boundary conditions. Laboratory experiments can provide a database that can be used to specify the boundary conditions in mathematical models of secondary cooling. Thermal boundary conditions that characterize the boiling heat transfer in secondary cooling are applicable for industrial processes after validation with reliable experiments.
Important work on spray cooling can be found in the literature [
1,
2,
3,
4,
5,
6]. The main cooling effect derives from evaporation. Heat transfer is influenced by water flow, droplet size, and velocity. Most studies either focused on the macroscopic effect of the sprays on the heat transfer rate or attempted to characterize the heat transfer processes during the impact of the spray on a surface. Many efforts have been made to understand the effect of changing the magnitude of the heat flux or the heat transfer coefficient (HTC) in different boiling heat transfer regimes. Thomas et al. [
7,
8] carried out laboratory measurements of water flow and heat transfer during spray cooling. Their research focused on the conditions of the surface of the steel strand in the secondary spray cooling zones with water jet–air mist cooling. The steel surface temperature range was 1200 °C to 200 °C. Ramstorfer et al. [
9] developed a dynamic spray cooling experimental platform where they measured the HTCs due to spray cooling using an experimental setup, allowing spray cooling up to a surface temperature of 1250 °C. Horský et al. [
10] developed experimental methods and numerical models for spray-cooled surface heat transfer. Those papers discussed heat transfer during spray cooling and optimization of the cooling process. Ito et al. [
11] also investigated the effects of the hydraulic pressure and water flow rate of a cooling water spray on cooling intensity and developed a more efficient secondary cooling system with a high-pressure water spray. Tsutsumi et al. [
12] did laboratory experiments on the cooling capacities of hydraulic and mist spray cooling by several kinds of spray nozzles. They proposed an equation for HTC that considered the spray thickness and collision pressure. El-Bealy [
13] studied the degree to which homogeneity of cooling conditions with air mist nozzles improved slab quality. The Brno University of Technology Heat Transfer and Fluid Flow Laboratory [
14] established reliable techniques to measure the effects of individual nozzles and combinations of nozzles to determine the HTCs. They also discussed methods for determining the HTCs and using them in a solidification model.
Most industrial spray nozzles are polydispersed, which allows them to generate water droplets with a wide range of sizes and velocities. This, in addition to the complex interaction of spray droplets, makes it difficult to predict (by single or multiple sprays) the performance of an industrial spray nozzle. Extensive laboratory studies have been done on the heat transfer of drops and sprays at low mass flux or for a single spray nozzle. However, data on industrial spray nozzles with high mass flux or on multiple spray nozzles are very scarce. Meanwhile, thinner samples with short cooling times are used in laboratory studies, which have led to a measured data shortage in high-temperature ranges. Due to fewer measurement data, the measurement error has become larger. It has become essential to use thicker samples and closer to the actual billet or slab. Experiments closely emulating the actual process of secondary cooling in continuous casting are particularly important.
In the paper, we developed a multifunctional experimental setup, enabling a quantitative understanding of the secondary cooling of a continuous casting. The objective was to explore and obtain data that characterized heat transfer at high temperatures using multiple air mist spray nozzles. The calculated temperature profile, heat flux and HTCs provided insight and useful data for the development of cooling strategies for continuous casting.
2. Establishing the Experimental Setup
Compared with industrial secondary cooling of continuous casting, a laboratory-scale simulator is cost effective because it is not necessary to interrupt actual steel production. The laboratory simulator provides more flexibility to change the range of operational parameters required for optimizing the secondary cooling process and developing new continuous casting technologies. By separating the secondary cooling from the casting in the experimental setup, we avoid the operational risk associated with using liquid metal. This enables us to simulate air-mist spray cooling in the secondary cooling process.
An overview of the system configuration is shown in
Figure 1a,b. The experimental setup consisted of an air–water spray, slab heating, slab cooling and data acquisition and analysis systems. To simulate the features of an industrial air mist spray, a laboratory air–water spray system was designed, as shown in
Figure 2. The system consisted of water tanks, pumps, an air compressor, an air pressure tank, electromagnetic flowmeters, temperature transmitters, pressure transmitters, electric control valves, metal hoses and spray racks.
Cooling water was supplied to the setup from the water tank. A valve was used to adjust the flow to the desired flowrate and pressure. The flow parameters were measured with the flowmeters and pressure gauges, which were controlled by the control software.
Nozzle alignments and cooling parameters are the main factors influencing the secondary cooling of continuous casting. The experimental setup was designed to meet the requirements of different nozzles with an extensive range of mass fluxes. The system could be adjusted to test the cooling effects of different spray configurations. The slab heating system comprised a feeding car, a heating furnace and a discharging electric actuator. After the desired temperature was reached, the slab was transported into the cooling station by a rack and pinion structure actuated by a motor. The slab cooling system included the cooling station, a slab depressing device, a reciprocating motion device, a drainage pipe, water pumps, electromagnetic flowmeters and the steam drainage system. It enabled spray rack replacement, slab movement simulation, adjustment of the slab cooling angle, and water and vapor discharge. Depending on the dynamic adjustment of the slab position by hydraulic pressure, it was possible to continuously simulate the position of the slab throughout the entire secondary cooling zone. The cooling station was placed in the positioner with 0° to 90° tilting to meet the slab arc change from initial mould position to final level position.
Figure 3 shows the top side cooling device that used three rows and four columns of nozzles with a 120° impingement angle. According to the nozzle angle, the overlap between the two nozzles was designed with nozzle manufacturer’s suggestion. The jets coalesced into a mist curtain on the slab surface.
In this way, the slab contact conditions (e.g., surface structure and clamping roller) could be manipulated to be similar to those used in industrial practice. There were four clamping rollers with a 120-mm diameter and a 160-mm separation distance that could clamp a slab with a 100- to 350-mm thickness.
Figure 4 shows the clamping rollers with the top spray configuration. During the tests, the slab surface was exposed to the air mist nozzles through the gap between the rollers. The heat was conducted to the slab surface and removed by the cooling air and water. The cooling process consumed much water and air, and heat was released during the process. Vapor and cooling water were collected and directed into a drain while the apparatus operated.
A data acquisition system was used to monitor the temperature, pressure, flow rate, and so forth during heating and cooling. It comprised a computer, a data acquisition box, K-type thermocouples and a software package, which were used for process parameter acquisition, calculation and data display. The thermocouples were inserted into the subsurfaces of the slabs. The slabs were then preheated to a desired initial temperature and transported to the cooling system. The measured temperature was used as an input for calculating the inverse heat conduction. By using the calculated results, we could evaluate the local heat flux at the boundaries. The heat flux was expressed via a boiling curve with different variables. The experimental setup provided the flexibility to adjust the operational parameters (mass flux, heating temperature, casting speed, etc.), enabling study of the heat transfer in secondary cooling. The technical parameters of the setup are shown in
Table 1.
3. Inverse Heat Conduction Problem
Normally, it is difficult to measure surface heat flux and temperature of a slab undergoing air mist spray cooling. When thermocouples are set at a certain distance from the surface of the slab to measure temperatures at different positions, a mathematical model can be used to calculate the surface heat flux and surface temperature [
15,
16]. The heat boundary condition can be calculated by recording the temperature as a function of time, which is in the form of an inverse heat transfer problem [
17].
We used thermocouples embedded in the slab at certain depths to measure cooling curves at different positions. The surface temperature and heat flux of the slab were calculated using the Beck’s sequential function specification method [
18,
19]. As a result of the calculation, the relation between the HTC and the slab surface temperature could be obtained. In this study, we regarded the internal heat transfer in the slab as one-dimensional unsteady heat conduction along the direction of the spray cooling, ignoring the heat loss from the side of the slab.
Figure 5 shows the physical model of the inverse heat conduction problem during the cooling.
The slab was heated to
T0 and cooled by the heat flux
q(
t). A thermocouple was embedded at
x = δ to record the temperature history at a time interval of Δ
t. The aim was to calculate
q(
t) using the measured temperature data with known initial conditions and the thermophysical properties of the materials involved in the system. The governing equation and boundary conditions were
where
ρ is the slab density,
δ is the distance between the outer surface of the slab and the thermocouple,
k is the conductivity,
T is the slab temperature,
T0 is the initial slab temperature,
Cp is heat capacity,
x is the thickness direction coordinate, and
Y(
t) is the measured temperature. In the model,
T0 and
Y(
t) were known; only the surface heat flux
q(
t) had to be solved. Since air mist cooling is a transient process, the temperature field and the heat flux
qM at a certain time
tM can be solved for if the temperature field
TM−1(
x) and the heat flux
qM−1 at
tM−1 have been determined. Assuming
q(
t) =
qM, which is constant when
tM−1 <
t <
tM, Equation (1) can be modified as Equation (2):
In this work, a sensitivity coefficient
was introduced to evaluate the sensitivity of the measured temperature point error. In Beck’s method, the heat flux
q(
t), is discretized into a series of
qi over a measurement interval Δ
t. The heat flux guess is kept constant within a period of time,
rΔ
t, where
r is the number of future time steps. Solved as a forward heat conduction problem, the predicted temperature by using the applied heat flux guess
qi can be obtained at
t =
i +
rΔ
t. In this way, the value of
qi in each time interval can be found to minimize the difference between the measured and calculated temperatures. The surface heat flux at different times is therefore calculated by defining the least squares error function. A program was written in Fortran to solve the one-dimensional transient heat conduction problem. In the process of calculation, the measurement interval Δ
t was 0.25 s, and the number of future time steps
r was 10. Thermocouples were embedded at
x = 20 mm to record the temperature history. The thermal properties of the materials used in the calculation were functions of temperature and are shown in
Table A1. In the course of calculation, it was considered that the steel density 7850 Kg/m
3 was constant. After programming, one of the parameters about the number of mesh points needed to be tested to investigate the grid independence of the results. When the number of mesh points was greater than 400, the results do not change. The number of mesh points used in this paper was 1000.
4. Experimental Results and Discussion
One of the goals of the setup was to obtain experimental data under industrial heat flux boundary conditions for describing the boiling heat transfer in the secondary cooling of continuous casting. Precise temperature measurement was a key for a reliable heat transfer measurement [
20]. Several 8-mm diameter holes were drilled 180 mm deep from the side of the 200-mm-thick slab perpendicular to the casting direction. The holes were drilled with flat bottoms to ensure that the tips of the thermocouples could touch the bottoms of the holes. A short distance from the cooling surface to the holes was required to reduce the time delay as the thermocouples responded to the changes in surface heat flux. In this experiment, the thermocouples were inserted 20 mm from the cooling surface. The thermocouples’ fastening devices were welded to the surface of the slab to push the thermocouples’ tips into contact with the bottoms of the holes. The gaps in the holes were filled with refractory cotton.
Figure 6 shows the armored K-type thermocouples that were inserted into the holes and “locked” inside by fastening devices. The locations of thermocouples numbered 1 to 5 are also shown in the figure. The thermocouples were far away from the edge of the slab to ensure a reduced heat loss effect. The thermocouple NO.1 and NO.3 were arranged inside the slab under the clamping roller respectively. The NO.2 was in the middle of the two clamping rollers. The slab was tilted 15° from the horizontal position to simulate the slab position.
The air mist spray characteristics, such as water distribution profiles, droplet velocity and the Sauer mean diameter, had an effect on heat transfer. These parameters depended on atomization and gas-water parameters of the nozzle. The main experiment parameters are shown in
Table 2. Due to the swing of the spray cooling racks and the overlapping between the nozzles, the uniformity of the water distribution in the vertical and horizontal casting direction was guaranteed. It could be considered that the slab surface cooling was homogeneous. The paper focused on the average convection heat transfer coefficient in secondary cooling. The average water flux was defined as a ratio of water flow to the area of the cooling surface of the slab. The water flux was 0.84 to 3.0 L/(m
2∙s) in the experiment. The nozzle arrangements with three rows and four columns were derived from the industrial continuous casting in this paper. The experimental nozzle model used was HPZ5.0-120B2.
Figure 7. shows a typical slab cooling in the experiment setup for a top spray configuration at a water temperature of 15.2 °C. Clearly, much vapor was generated during cooling, and the slab was cooled from the surface.
Typical temperature readings relative to time curves during the spray cooling are shown in
Figure 8. Regardless of the thermocouple location, the time–temperature curves were similar. A high cooling rate was observed, and the slab temperature dropped drastically. Due to the strong relation between the heat flux and the surface temperature, heat extraction rates changed rapidly with time.
The inverse heat transfer algorithm was applied to calculate heat fluxes by using the measured temperature profile data. The calculated heat fluxes were plotted as functions of surface temperatures at each thermocouple location to determine the boiling curve.
Figure 9 shows the relation between the surface temperature and the heat flux at different measuring positions. It can be seen that the heat flux was not linearly related to the surface temperature. The five thermocouples’ average heat flux of the slab first increased to a maximum value when the slab was cooled to approximately 574.2 °C, and then decreased from the peak value with further cooling. The maximum average heat flux was 5.0 × 10
5 W/m
2 at 572.4 °C. At 1000 °C, the average heat flux was 2.2 × 10
5 W/m
2. In the temperature range from 1100 to 800 °C, the average heat flux was 2.9 × 10
5 W/m
2. Comparing the heat flux at the different thermocouple positions, better vapor discharge conditions resulted in faster cooling and a higher heat flux.
As the simulated casting speed was 1.8 m/min, all thermocouples areas were cooled by sprays and support rolls. Each thermocouple reflected the average heat transfer on the surface. The NO.5 thermocouple was located in the lower part of the sample. Good steam exhaust conditions resulted in the greater heat flux.
The surface temperature history experienced two distinct regimes. The boiling curve allowed a better understanding of the physical state of the surface [
21]. The first turn in the curve coincided with the starting time of the spray cooling, from radiation to transition boiling. The second turn coincided with the critical heat flux (CHF), from transition boiling to nucleate boiling. At high surface temperatures, where film boiling was dominant, a vapor film near the solid hot surface minimized direct droplet contact time with the surface, resulting in a low heat-transfer rate. As the surface temperature decreased, the droplets began to penetrate the vapor film, and a sharp increase in the heat transfer rate was observed. After reaching the CHF, the heat flux decreased, going through the nucleate boiling regime and finally to convection or one-phase cooling.
The air mist spray heat transfer curve was similar to that of pool boiling in all the boiling regimes. However, we found that the cooling curve did not experience the Leidenfrost point due to the absence of stable film boiling. This can be explained by the following four circumstances. First, the distance between the thermocouple and the cooled surface was 20 mm, the thermocouple’s measurement delay resulted in the experiment not showing any Leidenfrost temperature. Second, the air and water flow rate are high in industrial conditions, but the residence time was short in the high-temperature region under experiment conditions. Therefore, the temperature point of Leidenfrost was not obvious. Third, the liquid droplets on the slab did not form a layer under the air mist spray cooling because the spray carried substantial momentum, which pushed the residual droplets away [
22]. Fourth, the cooling steel surface oxidized, forming an oxide layer that increased the slab roughness. The oxidation of the surface had the effect of raising the gasification nucleation number.
Figure 10 shows the relation between surface temperature and HTC. It could be seen that the HTC linearly varies with the temperature from 1100 to 500 °C. As temperatures decreased from 1100 to 800 °C, the average HTC linearly increased from 120.0 W/m
2∙K to 542.6 W/m
2∙K. At a CHF temperature of 574.2 °C, the average HTC was 888.9 W/m
2∙K.
Using the method developed by this paper, we measured the heat transfer boiling curve of another condition with water fluxes of 0.84 L/(m2∙s) and 3.0 L/(m2∙s). The effect of the cooling water flux was assessed by comparing the heat transfer characteristics of the tests.
Similar slopes for the transition boiling and nucleate boiling regime were observed for different water fluxes during cooling and shown in
Figure 11a. The HTC curve showed a typical linear relation, as shown in
Figure 11b. The heat flux with a water flux of 3.0 L/(m
2∙s) was higher, especially the CHF. There was a clear trend of increased CHF with water flux. Higher water flux led to more droplets breaking the vapor film and accelerating the steam discharge, strengthening the heat transfer on the slab surface.
The HTC as a function of the temperature difference between the steel surface and the water flux was found to be particularly useful to apply to the solidification model, as shown in
Figure 12. We concluded that the HTC was a close linear function of water mass flux for surface temperatures of 850 to 950 °C. The integrated HTC for the section can be described by the generic equation HTC =
aGb, where a and b are constants determined experimentally. The following relation for heat transfer in the transition film boiling region for surface temperatures of 850 to 950 °C is suggested:
HTC indicated the heat transfer efficiency between the surface of slab and cooling water, and the effect of heat transfer was high when HTC was large. In general, it needed to be measured by experiment and statistics, and expressed by empirical formula. Significant work regarding spray cooling could be found in the literature. Most researchers studied the macroscopic heat transfer rate. Different researchers had given different empirical formulas according to the test conditions and these formulas had different forms. Gan [
23] summarized some of the formula as follows:
HTC = 360 × G0.556, where 0.8 < G <2.5 (L/m2∙s), 727 °C < Ts < 1027 °C
HTC = 423 × G0.556, where 1 < G < 7 (L/m2∙s), 627 °C < Ts < 927 °C
HTC = 581G0.451 (1 − 0.0075Tw), where Tw is the water temperature
HTC = α(708G0.75Ts−1.2 + 0.116), where α is a calibration factor, Ts is the slab temperature
HTC = 157G0.55(1 − 0.0075Tw), where Tw is the water temperature
HTC = 130 + 350G
Compared with the formula in the paper, there were some differences between the researches. In some cases, the difference was striking. HTC was related to the factors of water flux, spray pressure, spray distance, nozzle structure, surface temperature of slab, water temperature etc. All of these factors had an impact on the relation.
An optimized transition boiling can help to improve heat transfer for specific spray characteristics. The test model was for carbon steel, but the results equally apply to other materials [
24]. However, deviations should be expected due to the surface characteristics of the spray-cooled surface [
25].
In summary, heat flux and HTC are critical for the theoretical heat design of an industrial process. The experimental setup provides a powerful tool for a quantitative understanding of the heat transportation phenomena of not only continuous casting but also other industries.