A Combined Hybrid 3-D / 2-D Model for Flow and Solidification Prediction during Slab Continuous Casting

A combined hybrid 3-D/2-D simulation model was developed to investigate the flow and solidification phenomena in turbulent flow and laminar flow regions during slab continuous casting (CC). The 3-D coupling model and 2-D slicing model were applied to the turbulent flow and laminar flow regions, respectively. In the simulation model, the uneven distribution of cooling water in the width direction of the strand was taken into account according to the nozzle collocation of secondary cooling zones. The results from the 3-D turbulent flow region show that the impact effect of the molten steel jet on the formation of a solidification shell is significant. The impact point is 457 mm below the meniscus, and the plug flow is formed 2442 mm below the meniscus. In the laminar flow region, grid independence tests indicate that the grids with a cell size of 10 × 10 mm2 are sufficient in simulations to attain the precise temperature distribution and solidification profile. The liquid core of the strand is not entirely uniform, and the solidification profile agrees well with the integrated distribution of cooling water in secondary cooling zones. The final solidification points are at a position of 400–500 mm in the width direction and are 17.66 m away from the meniscus.


Introduction
Continuous casting (CC) technology has become the primary method of producing steel strands in the steelmaking industry.During the CC process, the molten steel is continuously fed into the water-cooled mold through a submerged entry nozzle (SEN) and a solidified shell of sufficient thickness is formed.Subsequently, the strand is pulled into the secondary cooling zones and cooled by water spray or air-mist spray in order to solidify completely.The strand quality, particularly regarding surface and inner cracks, is closely related to the turbulent flow and the heat transfer during the solidification involved in a CC process [1][2][3][4][5].This is particularly true given that, in the slab CC process, the solidification profile of the slab transverse section-which closely relates to the integrated distribution of cooling water in secondary cooling zones-is not entirely uniform [6][7][8].This effect on centerline segregation is significant [7,9,10].
To date, both a two-dimensional (2-D) slicing model and three-dimensional (3-D) coupling model are widely used to predict the flow and solidification phenomena during the CC process.Owing to the high calculation efficiency, the 2-D slicing model is widely used to predict the temperature distribution and solidification profile during the CC process [7,[10][11][12].The calculation efficiency of the 2-D slicing model is improved by using an effective thermal conductivity concept to indirectly take the flow effect into account.However, the 2-D slicing model assumes that the solidification shell is a heat-transfer slice which moves from the meniscus to the solidification end.Evidently, the effect of turbulent flow on temperature distribution and the solidification shell is not directly considered in the 2-D slicing model.In fact, especially for slab the CC process, the molten steel jet which impacts the narrow face has a significant effect on the temperature distribution and solidification shell [13].Hence, to consider the effect of turbulent flow, the 3-D coupling model has been applied during the CC process [13][14][15][16][17]. Compared with the 2-D slicing model, the computational domain is much larger.Thereby, the amount of calculation is very large.Furthermore, the integrated distribution of cooling water in the width directions of secondary cooling zones is not considered in the previous 3-D calculations [13,17].This simple treatment of the heat transfer condition would evidently affect the accurate prediction of the solidification profile.Recently, Sun et al. [18] proposed a method to divide the computational domain of the CC bloom to improve the calculation efficiency.However, in that simulation model, the grid independence tests were not carried out and the heat transfer coefficient was assumed to be constant around the strand transverse section.Moreover, compared with the bloom, the heat transfer boundary of the CC slab is more complex.In the width direction, the distribution of the cooling water is not entirely uniform.Thus far, the combined hybrid 3-D/2-D numerical model has not been used in the simulation of a slab CC process.
In the present work, to simultaneously consider the effect of turbulent flow and improve calculation efficiency, a combined hybrid 3-D/2-D numerical model was established and used to explore the transport phenomena during the slab CC process.The 3-D coupling model and 2-D slicing model were adopted in the turbulent flow region and laminar flow region, respectively.The effects of flow in the turbulent region on the temperature distribution and the location of the final solidification point are considered by the data of the interface in the 3-D turbulent flow region.These are transmitted to the 2-D slice as the initial conditions.The grid independence tests were carried out to find a suitable mesh cell size.This was then adopted in the simulation of the laminar flow region.In the present model, the distribution of cooling water in the width direction was taken into account.

Mathematical Formulation
The continuity equation is ∂ρ ∂t The momentum equation is and where ρ, t, u i , u j , µ eff , P, g i , f L , A mushy , v i,p , T, T Solidus , and T Liquidus are molten steel density (kg/m 3 ), time (s), i-component of velocity (m/s), j-component of velocity (m/s), effective viscosity (kg•m −1 •s −1 ), pressure (N/m 2 ), i-component of acceleration due to gravity (m/s 2 ), liquid fraction, morphology constant, i-component of casting speed (m/s), temperature (K), solidus temperature (K), and liquidus temperature (K), respectively.The value of A mushy is usually between 10 5 and 10 8 in the numerical modelling of CC processes.In addition, the higher the value of A mushy , the steeper the transition of the velocity of the material to zero as it solidifies.Values between 10 4 and 10 7 are recommended for most computations (based on the user's guide of Ansys Fluent 14.0).Based on these two things, the value of A mushy was set to 10 7 in the present work.Standard k-ε equations are ∂(ρε) ∂t where where G, k, ε, µ, and µ t are the generation of turbulence kinetic energy (kg•m −1 •s −3 ), turbulent kinetic energy (m 2 /s 2 ), dissipation rate of turbulence energy (m 2 /s 3 ), molecular viscosity (kg•m −1 •s −1 ), and turbulent viscosity (kg•m −1 •s −1 ), respectively.The standard values of σ k , σ ε , c 1 , c 2 , and c µ recommended by Launder and Spalding are 1.0, 1.3, 1.44, 1.92, and 0.09, respectively.The energy equation is where where H, Pr t , h, L, λ, and c p , are enthalpy (J/kg), turbulent Prandtl number, sensible enthalpy (J/kg), pure solvent melting heat (J/kg), thermal conductivity (W•m −1 •K −1 ), and specific heat (J•kg −1 •K −1 ), respectively.The value of Pr t was set to 0.85 (based on the theory guide of Ansys Fluent 14.0).More details about the mathematical formulation are available in reference [14].

Computational Domain
According to the geometric symmetry, the strand was assumed to be ideally symmetrical and a quarter of the strand was selected to be calculated.The geometry and working parameters of the CC process are listed in Table 1.The whole 3-D computational domain (1530 × 190 × 20,362 mm 3 ) is very large.If the coupled model is adopted in it directly, the cost of computation will be very large.To cut down the great amount of calculation and improve the calculation efficiency, the computational domain is suitably divided into two parts, the 3-D turbulence region and the 2-D laminar flow region.The division method is illustrated in Figure 1.When the z-component velocity of molten steel is equal to the casting speed, the plug flow is formed.This means that the local molten steel and shell are in a relatively static status.That is, the first transverse section of the strand-where z-component velocities equal to the casting speed-is the interface between the turbulent flow region and laminar flow region.The data of the interface in the 3-D turbulent flow region are transmitted to the initial 2-D slice in the laminar flow region using a coordinate interpolation algorithm.

Boundary Conditions and Physical Properties
The inlet velocity of SEN was calculated based on mass conservation.The k and ε values were estimated using the semi-empirical relations.The values of these parameters are v in = 1.502 m/s, k = 0.003 m 2 /s 2 , ε = 0.0134 m 2 /s 3 , and T = 1811 K.The boundary condition of the outlet in the 3-D part was set as 'outflow'.The wall condition was employed on the meniscus without heat transfer.
The wall condition was employed on the strand surface.The heat flux of the mold surface was calculated using Equation ( 13) [19].
where L z and v c are the distance from the meniscus (m) and the casting speed (m/min), respectively.
In the secondary cooling zone, the cooling types include four aspects: water spray cooling, radiation cooling, water evaporation cooling, and roll contact cooling.In the present work, four aspects of cooling types have been considered.First, we determined the heat conditions for different cooling types with the same method described in our publications [7,9].The heat flux for the water spray cooling region was calculated using (14).More detailed descriptions about the other three cooling types in the secondary cooling zone are available in our publications [7,9].Then, to simplify this heat transfer process, we use the integrated heat transfer coefficient (in which the four cooling types have been taken into account) to describe the heat transfer process.
where w, T w , and T 0 are the spray water impingement density (L•m −2 •s −1 ), slab surface, and water temperature, respectively.Due to the existence of spray water overlap, the distribution of spray water for each row of nozzles along the width direction should be identified based on the nozzle collocation, as shown in Figure 3a.More details about this technique are available in our publication [7].According to the nozzle collocation of secondary cooling zones, the integrated distribution of cooling water from Zone 1 to Zone 8 is obtained, as shown in Figure 3b.In order to deal with the heat transfer condition conveniently, the strand surface in the width direction is divided into 26 pieces (sections No. 1 to No. 25 being 30 mm wide and section No. 26 being 15 mm wide) for heat transfer.The cooling water sprayed onto each piece is determined based on the nozzle collocation of each cooling zone.The cooling water on each piece is supposed to be uniformly distributed.
The physical properties of steel Q345 used in the calculation are listed in Table 2.

Boundary Conditions and Physical Properties
The inlet velocity of SEN was calculated based on mass conservation.The k and ε values were estimated using the semi-empirical relations.The values of these parameters are v in = 1.502 m/s, k = 0.003 m 2 /s 2 , ε = 0.0134 m 2 /s 3 , and T = 1811 K .The boundary condition of the outlet in the 3-D part was set as 'outflow'.The wall condition was employed on the meniscus without heat transfer.
The wall condition was employed on the strand surface.The heat flux of the mold surface was calculated using Equation ( 13) [19].
where L z and v c are the distance from the meniscus (m) and the casting speed (m/min), respectively.
In the secondary cooling zone, the cooling types include four aspects: water spray cooling, radiation cooling, water evaporation cooling, and roll contact cooling.In the present work, four aspects of cooling types have been considered.First, we determined the heat conditions for different cooling types with the same method described in our publications [7,9].The heat flux for the water spray cooling region was calculated using Equation (14).More detailed descriptions about the other three cooling types in the secondary cooling zone are available in our publications [7,9].Then, to simplify this heat transfer process, we use the integrated heat transfer coefficient (in which the four cooling types have been taken into account) to describe the heat transfer process.
where w, T w , and T 0 are the spray water impingement density (L•m 2 •s 1 ), slab surface, and water temperature, respectively.Due to the existence of spray water overlap, the distribution of spray water for each row of nozzles along the width direction should be identified based on the nozzle collocation, as shown in Figure 3a.More details about this technique are available in our publication [7].According to the nozzle collocation of secondary cooling zones, the integrated distribution of cooling water from Zone 1 to Zone 8 is obtained, as shown in Figure 3b.In order to deal with the heat transfer condition conveniently, the strand surface in the width direction is divided into 26 pieces (sections No. 1 to No. 25 being 30 mm wide and section No. 26 being 15 mm wide) for heat transfer.The cooling water sprayed onto each piece is determined based on the nozzle collocation of each cooling zone.The cooling water on each piece is supposed to be uniformly distributed.
The physical properties of steel Q345 used in the calculation are listed in Table 2.

Solution Procedure
The coupled calculation model was solved using the SIMPLE method with ANSYS Fluent 14.0.For both the 3-D turbulent flow region and the 2-D laminar flow region, the time step size was set to 0.05 s.Based on the relations between the length of model and casting speed, the total number of time steps for 3-D turbulent flow region and 2-D laminar flow region were 3500 and 17,802, respectively.The calculation was performed on a computer with a 3.50-GHz Intel Core i7-3770k processor and 16.0-GB RAM.For each time step, the tolerances of continuity, x-velocity, y-velocity, z-velocity, k, ε, and energy were set as 0.001, 0.001, 0.001, 0.001, 0.001, 0.001, and 1 × 10 −6 , respectively.

Flow and Solidification Phenomena in the Turbulent Flow Region
Figure 4a shows the melt flow pattern in the turbulent flow region.It is seen that the molten steel jet flows through the port of SEN and impacts on the narrow face of the mold.The molten steel jet is divided into two parts at the impact point and forms upper and lower recirculation vortexes.The fresh molten steel with higher temperature is brought to the impact point, where the solidification shell is remelting.The remelting phenomenon of the solidification shell is shown in Figure 4b.It is obvious that the impact effect of the molten steel jet on the formation of the solidification shell near the impact point is remarkable.

Solution Procedure
The coupled calculation model was solved using the SIMPLE method with ANSYS Fluent 14.0.For both the 3-D turbulent flow region and the 2-D laminar flow region, the time step size was set to 0.05 s.Based on the relations between the length of model and casting speed, the total number of time steps for 3-D turbulent flow region and 2-D laminar flow region were 3500 and 17,802, respectively.The calculation was performed on a computer with a 3.50-GHz Intel Core i7-3770k processor and 16.0-GB RAM.For each time step, the tolerances of continuity, x-velocity, y-velocity, z-velocity, k, ε, and energy were set as 0.001, 0.001, 0.001, 0.001, 0.001, 0.001, and 1 × 10 −6 , respectively.

Flow and Solidification Phenomena in the Turbulent Flow Region
Figure 4a shows the melt flow pattern in the turbulent flow region.It is seen that the molten steel jet flows through the port of SEN and impacts on the narrow face of the mold.The molten steel jet is divided into two parts at the impact point and forms upper and lower recirculation vortexes.The fresh molten steel with higher temperature is brought to the impact point, where the solidification shell is remelting.The remelting phenomenon of the solidification shell is shown in Figure 4b.It is obvious that the impact effect of the molten steel jet on the formation of the solidification shell near the impact point is remarkable.The impact point of the jet on the narrow face is obtained from Figure 5a,b.At the impact point, the z-component of velocity is equal to zero and the velocity magnitude reaches a valley value.In the present case, the impact point is at the position 457 mm below the meniscus.In order to find the interface between the turbulent flow region and the laminar flow region, the z-component of velocities of characteristic lines at the central symmetry plane of the strand narrow face were investigated.
From Figure 5c, it can be seen that the z-component of the molten steel velocity is equal to the casting speed when arriving at the location 2442 mm below the meniscus, where it is roughly placed at the end of secondary cooling Zone 3.
To deal with the heat transfer condition for each secondary cooling zone separately, the transverse section at the end of secondary cooling Zone 3 (2560 mm below meniscus) is selected as the interface between the turbulent flow region and laminar flow region.Therefore, the 2-D slicing model is adopted in the laminar flow region from Zone 4 to Zone 8.
Metals 2018, 8, x FOR PEER REVIEW 7 of 10 present case, the impact point is at the position 457 mm below the meniscus.In order to find the interface between the turbulent flow region and the laminar flow region, the z-component of velocities of characteristic lines at the central symmetry plane of the strand narrow face were investigated.From Figure 5c, it can be seen that the z-component of the molten steel velocity is equal to the casting speed when arriving at the location 2442 mm below the meniscus, where it is roughly placed at the end of secondary cooling Zone 3.
To deal with the heat transfer condition for each secondary cooling zone separately, the transverse section at the end of secondary cooling Zone 3 (2560 mm below meniscus) is selected as the interface between the turbulent flow region and laminar flow region.Therefore, the 2-D slicing model is adopted in the laminar flow region from Zone 4 to Zone 8.

Grid Independence Tests for the Laminar Flow Region
To begin with, grid sizes for the laminar flow region have been carefully chosen to ensure the grid independence solution.As seen in the results of the grid independence tests shown in Figure 6, relatively larger differences exist in the solidification profile and temperature extreme in the transverse section (17.56 m from the meniscus) with different cell sizes.It is obvious that the cell size 25 × 25 mm 2 , as shown in Figure 6a, is too large for this case to precisely predict the solidification profile and temperature distribution.The solidification profile does not correspond with the integrated distribution of cooling water in secondary cooling zones (shown in Figure 3).Also, the lowest temperature is higher than the other cases.In comparing Figure 6b,c, where the cell size decreases from 15 × 15 mm 2 to 10 × 10 mm 2 , there is a similar solidification profile which is corresponding with the integrated distribution of cooling water in the width direction of strand.The maximum and minimum temperatures in the case where the cell size is 15 × 15 mm 2 are very close to those in the case where the cell size is 10 × 10 mm 2 .When the cell size reduces from 10 × 10 mm 2 to 5 × 5 mm 2 , as shown in Figure 6d, the temperature extreme of 1181 K to 1733 K shows little change compared with that using grids with a cell size of 10 × 10 mm 2 .Thereby, in the present simulation, the grids with cell size 10 × 10 mm 2 are sufficient to obtain the precise temperature distribution and solidification profile.

Grid Independence Tests for the Laminar Flow Region
To begin with, grid sizes for the laminar flow region have been carefully chosen to ensure the grid independence solution.As seen in the results of the grid independence tests shown in Figure 6, relatively larger differences exist in the solidification profile and temperature extreme in the transverse section (17.56 m from the meniscus) with different cell sizes.It is obvious that the cell size 25 × 25 mm 2 , as shown in Figure 6a, is too large for this case to precisely predict the solidification profile and temperature distribution.The solidification profile does not correspond with the integrated distribution of cooling water in secondary cooling zones (shown in Figure 3).Also, the lowest temperature is higher than the other cases.In comparing Figure 6b,c, where the cell size decreases from 15 × 15 mm 2 to 10 × 10 mm 2 , there is a similar solidification profile which is corresponding with the integrated distribution of cooling water in the width direction of strand.The maximum and minimum temperatures in the case where the cell size is 15 × 15 mm 2 are very close to those in the case where the cell size is 10 × 10 mm 2 .When the cell size reduces from 10 × 10 mm 2 to 5 × 5 mm 2 , as shown in Figure 6d, the temperature extreme of 1181 K to 1733 K shows little change compared with that using grids with a cell size of 10 × 10 mm 2 .Thereby, in the present simulation, the grids with cell size 10 × 10 mm 2 are sufficient to obtain the precise temperature distribution and solidification profile.

Temperature Field and Solidification Profile
The slab center and wide surface center temperature along the casting direction predicted by the present model are compared to that predicted by our previous model, which is proven and widely employed in temperature prediction during the CC process.As the results show in Figure 7, the temperature predicted by the present model agrees well with that of the models by Long et al. [7] and Long and Chen [9].The slab center temperature, as shown in Figure 7a, decreases slowly before complete solidification.However, after complete solidification, the slab center temperature decreases sharply.This is because there is no latent heat to maintain the temperature when the local region solidifies completely.The wide surface center temperature is shown in Figure 7b.Due to the strong cooling of the mold, the surface temperature decreases radically, and then recalesces when the strand is pulled out of the mold.In the secondary cooling zone, the surface temperature decreases gradually.

Temperature Field and Solidification Profile
The slab center and wide surface center temperature along the casting direction predicted by the present model are compared to that predicted by our previous model, which is proven and widely employed in temperature prediction during the CC process.As the results show in Figure 7, the temperature predicted by the present model agrees well with that of the models by Long et al. [7] and Long and Chen [9].The slab center temperature, as shown in Figure 7a, decreases slowly before complete solidification.However, after complete solidification, the slab center temperature decreases sharply.This is because there is no latent heat to maintain the temperature when the local region solidifies completely.The wide surface center temperature is shown in Figure 7b.Due to the strong cooling of the mold, the surface temperature decreases radically, and then recalesces when the strand is pulled out of the mold.In the secondary cooling zone, the surface temperature decreases gradually.

Temperature Field and Solidification Profile
The slab center and wide surface center temperature along the casting direction predicted by the present model are compared to that predicted by our previous model, which is proven and widely employed in temperature prediction during the CC process.As the results show in Figure 7, the temperature predicted by the present model agrees well with that of the models by Long et al. [7] and Long and Chen [9].The slab center temperature, as shown in Figure 7a, decreases slowly before complete solidification.However, after complete solidification, the slab center temperature decreases sharply.This is because there is no latent heat to maintain the temperature when the local region solidifies completely.The wide surface center temperature is shown in Figure 7b.Due to the strong cooling of the mold, the surface temperature decreases radically, and then recalesces when the strand is pulled out of the mold.In the secondary cooling zone, the surface temperature decreases gradually.The temperature distribution in transverse sections with different distances from the meniscus (as shown in Figure 8a) have a similar profile.The temperature decreases as the distance from the meniscus increases.Figure 8b shows the solidification profile evolution process during continuous casting.It is obvious that the liquid core of the strand is not entirely uniform.At the local region of the strand in the width direction, from roughly 400 to 500 mm, the molten steel solidifies completely owing to the minimum amount of the cooling water sprayed on it in the secondary cooling zone.The final solidification points (another one in the other half of the geometry) are 17.66 m away from the meniscus.Moreover, at the position about 150 mm in width direction, the molten steel solidifies completely a little earlier than that at the region from 400 to 500 mm.The reason is that the amount of cooling water at the former position is a little more than that at the region from 400 to 500 mm.This results from the valley of the distribution curve being very narrow at the position of about 150 mm, although the proportion of cooling water is slightly lower than that at the region from 400 to 500 mm.The solidification profile agrees well with the integrated distribution of cooling water in the secondary cooling zones.The temperature distribution in transverse sections with different distances from the meniscus (as shown in Figure 8a) have a similar profile.The temperature decreases as the distance from the meniscus increases.Figure 8b shows the solidification profile evolution process during continuous casting.It is obvious that the liquid core of the strand is not entirely uniform.At the local region of the strand in the width direction, from roughly 400 to 500 mm, the molten steel solidifies completely owing to the minimum amount of the cooling water sprayed on it in the secondary cooling zone.The final solidification points (another one in the other half of the geometry) are 17.66 m away from the meniscus.Moreover, at the position about 150 mm in width direction, the molten steel solidifies completely a little earlier than that at the region from 400 to 500 mm.The reason is that the amount of cooling water at the former position is a little more than that at the region from 400 to 500 mm.This results from the valley of the distribution curve being very narrow at the position of about 150 mm, although the proportion of cooling water is slightly lower than that at the region from 400 to 500 mm.The solidification profile agrees well with the integrated distribution of cooling water in the secondary cooling zones.

Conclusions
A combined hybrid 3-D/2-D simulation model with a high calculation efficiency has been developed to investigate the flow and solidification behaviors during the slab CC process.The main conclusions are summarized as follows: (1) The impact effect caused by the molten steel jet on the formation of a solidification shell is significant.The impact point is at the position 457 mm below the meniscus, and the plug flow is formed 2442 mm below the meniscus.(2) For the simulation of the laminar flow region, the grids with a cell size of 10 × 10 mm 2 are sufficient to attain a precise temperature distribution and solidification profile.(3) The solidification profile of the strand is not entirely uniform.The final solidification points, roughly being at the position from 400-500 mm in the width direction, are 17.66 m away from the meniscus.

Figure 1 .
Figure 1.Illustration of 3-D turbulent flow and 2-D laminar flow regions during the CC (continuous casting) process.In order to find out the interface between the 3-D turbulent flow region and the 2-D laminar region, a computational domain with a size of 1530 × 190 × 3250 mm 3 was built and calculated.The position of this interface can be obtained by analyzing the z-component velocity variation of molten steel at the central symmetry plane of strand narrow face.The meshes created by ANSYS ICEM CFD 14.0 and adopted in 3-D turbulent flow region and 2-D laminar flow region are shown in Figure 2. The number of elements for the meshes of the 3-D part and 2-D part are 517,793 and 770, respectively.

Figure 2 .
Figure 2. Meshes adopted in the (a) 3-D turbulent flow region and (b) 2-D laminar flow region.

Figure 1 .
Figure 1.Illustration of 3-D turbulent flow and 2-D laminar flow regions during the CC (continuous casting) process.

Figure 1 .
Figure 1.Illustration of 3-D turbulent flow and 2-D laminar flow regions during the CC (continuous casting) process.In order to find out the interface between the 3-D turbulent flow region and the 2-D laminar region, a computational domain with a size of 1530 × 190 × 3250 mm 3 was built and calculated.The position of this interface can be obtained by analyzing the z-component velocity variation of molten steel at the central symmetry plane of strand narrow face.The meshes created by ANSYS ICEM CFD 14.0 and adopted in 3-D turbulent flow region and 2-D laminar flow region are shown in Figure 2. The number of elements for the meshes of the 3-D part and 2-D part are 517,793 and 770, respectively.

Figure 2 .
Figure 2. Meshes adopted in the (a) 3-D turbulent flow region and (b) 2-D laminar flow region.Figure 2. Meshes adopted in the (a) 3-D turbulent flow region and (b) 2-D laminar flow region.

Figure 2 .
Figure 2. Meshes adopted in the (a) 3-D turbulent flow region and (b) 2-D laminar flow region.Figure 2. Meshes adopted in the (a) 3-D turbulent flow region and (b) 2-D laminar flow region.

Figure 3 .
Figure 3. (a) Schematic of nozzle collocation; (b) integrated distribution of cooling water from Zone 1 to Zone 8.

Figure 3 .
Figure 3. (a) Schematic of nozzle collocation; (b) integrated distribution of cooling water from Zone 1 to Zone 8.

Figure 4 .
Figure 4. (a) Flow field in the turbulent flow region; (b) liquid fraction on the different transverse sections.The impact point of the jet on the narrow face is obtained from Figure5a,b.At the impact point, the z-component of velocity is equal to zero and the velocity magnitude reaches a valley value.In the

Figure 4 .
Figure 4. (a) Flow field in the turbulent flow region; (b) liquid fraction on the different transverse sections.

Figure 5 .
Figure 5. (a) z-component of velocity and (b) velocity magnitude variation of the line with 10 mm from the narrow face at the central symmetry plane.(c) z-component of the velocity variation of different lines at the central symmetry plane of the strand narrow face.

Figure 5 .
Figure 5. (a) z-component of velocity and (b) velocity magnitude variation of the line with 10 mm from the narrow face at the central symmetry plane.(c) z-component of the velocity variation of different lines at the central symmetry plane of the strand narrow face.

Figure 7 .
Figure 7. (a) Slab center and (b) wide surface center temperature along the casting direction.

Figure 7 .
Figure 7. (a) Slab center and (b) wide surface center temperature along the casting direction.Figure 7. (a) Slab center and (b) wide surface center temperature along the casting direction.

Figure 7 .
Figure 7. (a) Slab center and (b) wide surface center temperature along the casting direction.Figure 7. (a) Slab center and (b) wide surface center temperature along the casting direction.

Figure 8 .
Figure 8.(a) Temperature field and (b) solidification profile in transverse sections with different distances from the meniscus.

Table 1 .
The geometry and working parameters of the CC (continuous casting) process.

Table 1 .
The geometry and working parameters of the CC (continuous casting) process.

Table 1 .
The geometry and working parameters of the CC (continuous casting) process.

Table 2 .
Physical properties of steel Q345.

Table 2 .
Physical properties of steel Q345.