Flow Field in Slab Continuous Casting Mold with Large Width Optimized with High Temperature Quantitative Measurement and Numerical Calculation

: In this paper, the rod deﬂection method was applied to quantitatively measure velocity near the mold surface at high temperatures and the k- ε model coupled with a discrete phase model (DPM) was adopted to simulate the ﬂow ﬁeld in the mold. The calculated results match very well with the measured results under all the present conditions. Under the conditions of the large mold width of 1800 mm, 1.1 m/min casting speed and 140 mm submerged entry nozzle (SEN) immersion depth, the velocity near the mold surface decreases with increasing the argon gas ﬂow rate. When the argon gas ﬂow rate is 6 L/min, the ﬂow pattern is the double roll ﬂow (DRF). When the argon gas ﬂow rate is increased to 10 L/min and 14 L/min, the ﬂow pattern is the single roll ﬂow (SRF), and the risk of slag entrainment increases. With an argon gas ﬂow rate of 10 L/min, and an immersion depth of 160 mm, the velocity near the mold surface sensitively increases with increasing the casting speed. When the casting speed is 1.1 m/min, an intermediate ﬂow (IF) is formed with the intensiﬁed mold surface ﬂuctuation, which can easily result in slag entrainment defects. When the casting speed is only increased to 1.2 m/min, the velocity near the mold surface changes drastically and is close to the upper limit velocity of 0.4 m/s. When the casting speed is 1.1 m/min, and the argon gas ﬂow rate is 10 L/min, the velocity near the mold surface is obviously increased with increasing the immersion depth. When the immersion depth of the nozzle increases from 140 mm and 160 mm to 180 mm, the ﬂow pattern changes from SRF or IF to DRF. When the bottom shape of the SEN changes from mountain to well, the velocity near the mold surface decreases. We suggest adopting the well-bottom nozzle to reduce the risk of slag entrainment.


Introduction
With the development of the automobile industry, the demand for an exposed panel with high surface quality is on the rise. In order to reduce the incidence of linear defects on the surface of the automobile exposed panel, studies aiming to find improvements in steelmaking and continuous casting (CC) technologies have been continuously carried out by many integrated steel companies in the world. Generally, the width of the mold usually ranges from 880 to 2300 mm and the defect ratio of the slab with a width larger than 1600 mm shares a high proportion, according to the practical industrial experiments. Since there is a close relationship between the defects caused by the CC process and the flow field in the mold, the flow field in a CC mold with a width larger than 1600 mm requires urgent study.
The effects of the parameters of the CC process and the structure of the nozzle on the flow field in the CC mold with a wide cross-section have been studied extensively.
As regards the investigations with water model and numerical simulations, Bessho [1] studied the flow field in a mold with a width of 1250-2000 mm by means of a numerical simulation and water model, and found that the two factors determining the flow pattern of molten steel in the mold are the kinetic energy of molten steel at the side port of submerged entry nozzle (SEN) and the buoyancy force generated by the argon gas bubbles. Chen [2][3][4], He [5] and Liu [6] used a water model to study the flow field in a mold with four different widths ranging from 1800 to 3250 mm and concluded that increasing the air blowing volume easily leads to an increase in fluctuation in the liquid level at 1/4 the width of the mold. Lin [7] reported that the surface velocity of a slab CC mold changed significantly with the casting speed when using a large slab width of 1650 mm and the water model. Wang [8] investigated the flow field of liquid steel in the mold, including large widths of 1850 and 2150 mm, and indicated that the optimal argon gas flow rate is positively correlated with the casting speed. Hernandez [9] used the k-ε model coupled with the volume of fluid model (VOF) to find that the rectangular side port of SEN can effectively reduce the occurrence of backflow compared with the circular side port. As regards the flow field in the mold with a width of 1736 mm, by using a Large Eddy Simulation (LES), Liu [10] found that the lower circulation area in the mold is usually asymmetric, while the upper circulation area has two typical flow patterns, namely clockwise flow and counter clockwise flow. For the flow field in a mold with a width of 1880 mm, Ramos [11] pointed out that the flow pattern in the mold depends on the kinetic energy of the steel stream at the side port of the SEN and the peak velocity is near the side port. Xiong [12] used a water model to study the influence of the side port shape of SEN on the level fluctuation on the surface of the mold with a width of 1800 mm, and demonstrated that the amplitude of the liquid level fluctuation is in the order of elliptic nozzle < rectangular nozzle < square nozzle. Ren [13] found that the average liquid level contour and flow field tend to be symmetrical when the casting speed is low for slab continuous casting in a mold with a width of 2000 mm.
As regards the investigations including industrial experiments, Kato [14] utilized the water model and industrial experiments to study the flow field in a mold with a width of 1600 mm, and proposed that two optimized conditions for the flow field are that the liquid steel must fill the nozzle and that the flow must be a piston flow in the nozzle. Again, using the water model and industrial experiments to study the flow field in molds with widths of 1650 mm and 1850 mm, Tsukaguchi [15] indicated that the SEN with swirling flow can increase the CC speed and improve the slab surface quality. Deng [16][17][18] conducted a large number of tests with the water model and industrial experiments for a wide slab CC mold, and concluded that the ratio of casting speed to mold width and argon gas flow rate is the decisive factor affecting the flow pattern in the mold, and the critical argon gas flow rate from double roll flow (DRF) to single roll flow (SRF) increases significantly with increasing casting speed. For the 1930 mm mold width, Chen [19] demonstrated that the advantage of the concave bottom nozzle is able to slow down the appearance of the asymmetric flow on the two side ports of the SEN. Cho [20] made use of LES and VOF models, in addition to industrial experiments to study the flow field in a mold with a width of 1600 mm, and found that an SEN with a negative inclination angle is helpful for the formation of DRF.
In order to study the flow field in the mold more accurately and directly, a large number of researchers have developed various velocity measurement methods to measure the velocity of molten steel in the mold. Iguchi [21] put forward a measurement method to measure molten metal with a low melting point and spherical velocimeter. The velocity of the molten metal can be calculated by measuring the resistance of a metal ball in the molten metal. Iguchi [22] proposed another method for measuring the flow rate of molten steel by measuring the shedding frequency of a Karman vortex street. The resistance measurement method could also obtain the liquid steel flow rate by measuring the resistance of a nondeflectable rod inserted into the liquid steel. Thomas [23] used the nail dipping method to obtain the velocity of liquid steel and free surface profile of the mold by measuring the slope and direction of the solidified liquid steel on the nail plate.
However, the quantitative high temperature measurements of velocities near the mold surface under the different CC parameters, such as the different casting speeds, argon gas flow rates, immersion depths of SEN and different nozzle bottom shapes were rarely carried out together with a comparison of the numerical simulation results. These quantitative results are very important for steel plants to optimize the flow field in the mold to minimize the slag entrainment defects for automobile exposed panel production.
In our previous research studies, the rod deflection method [24] was proposed as a technique to obtain the velocity near the mold surface by measuring the deflection angle of a detecting rod inserted into the liquid steel. With the combination of the measurement of velocity near the mold surface and mathematical modelling, the flow fields in the mold with the narrow width [25] and the medium width [26] were also investigated under the different CC parameters. In this paper, under the different CC parameters of the casting speed, argon gas flow rate, SEN immersion depth and bottom shape, the rod deflection method was used to quantitatively measure the velocity near the surface of the mold at a high temperature. Furthermore, the standard k-ε coupled with DPM model was simultaneously used to quantitatively study the effects of the above CC parameters on the velocity near the liquid surface, flow pattern, and liquid level fluctuation in the mold with a large width of 1800 and 2000 mm. Here. the influence of various CC parameters on the flow field is clarified, which is helpful in order to reduce the incidence of slag entrainment defects in the CC process with molds with large widths for producing automobile exposed panels.

High Temperature Velocity Measurement Method
For the ultra low carbon steel used for producing automobile exposed panels, a high temperature velocity measurement was conducted. In order to measure the velocity near the mold surface more accurately and simply, the rod deflection method was put forward in our previous research [24][25][26][27]. In Figure 1, the velocity measuring device of the rod deflection method is composed of a balance block, deflection bearing, deflection angle indicator and detecting rod. The function of the balance block is to make the center of gravity of the whole velocity measuring device close to the center of the deflection shaft, which means that a slight change in the liquid steel flow rate can be reflected in the deflection angle of the detecting rod. When measuring the velocity near the surface of the mold, it is only necessary to insert the detecting rod into the position at the 1/4 width of the mold. In this process, the detecting rod should be perpendicular to the mold surface. The impact of the molten steel flow on the rod will cause its deflection of a certain angle. The stainless steel is used for the detecting rod material because it has no magnetism and is not affected by the electromagnetic field. One stainless steel detecting rod with a diameter of 10 mm can stay in the molten steel for about 30 s before it melts and 10 deflection angles can be collected after that the detecting rod tend to be stable in the molten steel. Under each process condition, three detecting rods were used to obtain more than 30 deflection angles.
When the three forces of gravity, buoyancy force and impact force acting on the detecting rod reach a balanced state, the detecting rod leans to a certain deflection angle (θ), the relationship of three forces can be described by the following equation: where G is the gravity of the detecting rod (N). L 1 is the distance between the barycenter of detecting rod and the rotational pivot (m). F f is the buoyancy force on detecting rod (N). F D is the impact force on detecting rod (N). L 2 is the distance between the acting point of the impact force and the rotational pivot (m). The impact force can be obtained by the Equation (2), because the impact force of liquid steel on the detecting rod is equal to the drag force. where U 0 is the velocity near the mold surface (m/s). C D is the resistance coefficient of the flow. ρ is the density of the molten steel (kg/m 3 ). A is the projected area of the immersion part of the detecting rod in the direction vertical to the flow of molten steel (m 2 ). So, the velocity near the mold surface can be calculated from the deflection angle (θ) of the detecting rod from Equation (3).
The velocities near the mold surface can be quantitatively measured at high temperatures by using this method. For a detailed description of the measuring method of the rod deflection method refer to our previous papers [24][25][26][27]. When the three forces of gravity, buoyancy force and impact force acting on the detecting rod reach a balanced state, the detecting rod leans to a certain deflection angle (θ), the relationship of three forces can be described by the following equation: where G is the gravity of the detecting rod (N). L1 is the distance between the barycenter of detecting rod and the rotational pivot (m). Ff is the buoyancy force on detecting rod (N). FD is the impact force on detecting rod (N). L2 is the distance between the acting point of the impact force and the rotational pivot (m). The impact force can be obtained by the Equation (2), because the impact force of liquid steel on the detecting rod is equal to the drag force.
where U0 is the velocity near the mold surface (m/s). CD is the resistance coefficient of the flow. ρ is the density of the molten steel (kg/m 3 ). A is the projected area of the immersion part of the detecting rod in the direction vertical to the flow of molten steel (m 2 ). So, the velocity near the mold surface can be calculated from the deflection angle (θ) of the detecting rod from Equation (3).

Numerical Simulation
In this paper, a standard k-ε model coupled with the DPM model aiming to build a connection between the continuous phase and the discrete phase, was applied to simulate the turbulent flow in the SEN and CC mold. In order to simplify the complex solution process, some reasonable assumptions were applied in the numerical simulation in this paper. These include: that molten steel is an incompressible fluid; the argon gas bubbles are rigid spheres with a uniform diameter, which are not affected by pressure and temperature; the fusion between bubbles and bubble rupture are ignored; the discrete phases such as argon gas bubbles do not occupy the volume in the computational domain and are regarded as point particles; the influence of the solidified shell of slab on the flow field is ignored; the influence of mold powder on the flow field is ignored; the mold oscillations are not considered.

Governing Equation
In the calculation of the liquid steel and argon two-phase flow in the mold, the Euler method was used for the liquid steel continuous phase and the Lagrange method was used for the argon gas discrete phase. The discrete phase and continuous phase are coupled by source term F in governing Equation (5). The continuous phase and the discrete phase are a bidirectional coupling. The movement of the particles is affected by the continuous phase, and conversely, the particles also affect the movement of the continuous phase.
The mass conservation (Equation (4)) and the momentum conservation (Equation (5)) equations are used to calculate the motion of continuous phase.

∂ρ ∂t
where u is the flow velocity (m/s); ρ is the density of molten steel (kg/m 3 ); p is the pressure; µ l + µ t is the effective viscosity (Pa·s); g is the gravitational acceleration (m/s 2 ); F is the force of the argon gas bubbles acting on the molten steel (N/m 3 ). In this paper, the turbulent flow is simulated with the standard k-ε model. The turbulent viscosity µ t is a function of the turbulent kinetic energy k and the turbulent energy dissipation rate ε.
The standard k-ε model ignores the effect of molecular viscosity, and establishes two transport equations of k and ε.
where the constants in the model are C µ = 0.09, σ k = 1.00, σ ε = 1.30, C 1 = 1.44, C 2 = 1.92. G k is the turbulent kinetic energy generation, which is expressed as: To calculate F for Equation (5), the DPM model solves a force balance on each argon gas bubble governed by Newton's second law and can be described as following: where m b is the mass of bubble (kg); µ b is the velocity of bubble (m/s); F d is the drag force (N); F g is the gravitational force (N); F b is the buoyancy force (N); F v is the virtual mass force (N); F p is the pressure gradient force (N); F L is the lift force (N). All the detailed information about the forces mentioned above has been carefully explained in our previous work [24]. Figure 2 shows the geometry and computational meshes of the mold. The computational domain consists of 680,000 cells. In particular, the mesh size in the SEN is refined to calculate the fluid flow more accurately. The orthogonal quality is 0.47 (minimum) and the skewness is 0.52 (maximum). The parameters of the experiment and numerical simulation for each part are presented in Table 1.

Computational Domain and Boundary Conditions
For continuous phase molten steel, the inlet velocity can be regarded as a fixed value due to the mass conservation between the inlet and outlet. The outlet boundary condition is set to be the pressure outlet. The stable sliding surface is applied to be the free surface boundary condition. Additionally, the other wall boundary conditions are taken as the stable non sliding surface.
For the discrete phase of the argon gas bubbles, the argon gas bubbles are set as rigid spheres with a diameter of 1 mm. The boundary conditions of the top free surfaces and the outlet of the mold are set to escape, which means that the discrete phase particles will escape from the computational domain when touching these regions. The boundary condition of the SEN is taken as a reflection. Namely, the particles will rebound away from the wall with the recovery coefficient specified by the user. All the boundary conditions of the remaining walls are set to be a trap, and particles will be adsorbed on the walls when contacting with the walls. The convergence criterion is set to be that the residual value is smaller than 0.001. The unsteady simulations are chosen to describe the flow of the molten steel in the mold because of the transient flow field in the actual mold. For the present work, the time step for the simulation is 0.005 s, and the total time of calculation for the continuous casting process is 30 s; when the flow field in the mold has reached a relatively stable state. Consequently, the simulated results after calculation for 30 s are presented in this paper. Figure 2 shows the geometry and computational meshes of the mold. The computational domain consists of 680,000 cells. In particular, the mesh size in the SEN is refined to calculate the fluid flow more accurately. The orthogonal quality is 0.47 (minimum) and the skewness is 0.52 (maximum). The parameters of the experiment and numerical simulation for each part are presented in Table 1.
For continuous phase molten steel, the inlet velocity can be regarded as a fixed value due to the mass conservation between the inlet and outlet. The outlet boundary condition is set to be the pressure outlet. The stable sliding surface is applied to be the free surface boundary condition. Additionally, the other wall boundary conditions are taken as the stable non sliding surface.
For the discrete phase of the argon gas bubbles, the argon gas bubbles are set as rigid spheres with a diameter of 1 mm. The boundary conditions of the top free surfaces and the outlet of the mold are set to escape, which means that the discrete phase particles will escape from the computational domain when touching these regions. The boundary condition of the SEN is taken as a reflection. Namely, the particles will rebound away from the wall with the recovery coefficient specified by the user. All the boundary conditions of the remaining walls are set to be a trap, and particles will be adsorbed on the walls when contacting with the walls. The convergence criterion is set to be that the residual value is smaller than 0.001. The unsteady simulations are chosen to describe the flow of the molten steel in the mold because of the transient flow field in the actual mold. For the present work, the time step for the simulation is 0.005 s, and the total time of calculation for the continuous casting process is 30 s; when the flow field in the mold has reached a relatively stable state. Consequently, the simulated results after calculation for 30 s are presented in this paper.

Results and Discussion
In this paper, the rod deflection method was used to measure the velocity near the mold surface under the different CC parameters (casting speed, argon gas flow rate, SEN immersion depth and bottom shape). Moreover, the calculated and measured results under these conditions are compared and analyzed. The effects of the different CC parameters on the flow pattern and liquid level fluctuation in the mold are further studied.

Effect of Argon Gas Flow Rate
In Figure 3, the calculated velocities near the mold surface are compared with the measured ones under the different argon gas flow rates. When the cross-section width is 1800 mm, the casting speed is 1.1 m/min, and the nozzle immersion depth is 140 mm. The calculated velocity is the velocity 10 mm below the top surface of the mold. It is clear that the calculated results match well with the measured results. With increasing argon gas flow rates, both the measured and calculated values show a decreasing trend. When the argon gas flow rate is 6 L/min, both the measured and calculated values are about 0.25 m/s, which indicates that the molten steel near the mold surface flows from the narrow wall to the SEN, so the flow pattern in the mold is the DRF. When the argon gas flow rates are 10 L/min and 14 L/min, the velocity near the mold surface is smaller than −0.2 m/s, indicating the flow pattern in the mold may be the intermediate flow (IF) or SRF. Therefore, the argon gas flow rate has a great influence on the flow pattern in the mold with a large width.

Results and Discussion
In this paper, the rod deflection method was used to measure the velocity near the mold surface under the different CC parameters (casting speed, argon gas flow rate, SEN immersion depth and bottom shape). Moreover, the calculated and measured results under these conditions are compared and analyzed. The effects of the different CC parameters on the flow pattern and liquid level fluctuation in the mold are further studied.

Effect of Argon Gas Flow Rate
In Figure 3, the calculated velocities near the mold surface are compared with the measured ones under the different argon gas flow rates. When the cross-section width is 1800 mm, the casting speed is 1.1 m/min, and the nozzle immersion depth is 140 mm. The calculated velocity is the velocity 10 mm below the top surface of the mold. It is clear that the calculated results match well with the measured results. With increasing argon gas flow rates, both the measured and calculated values show a decreasing trend. When the argon gas flow rate is 6 L/min, both the measured and calculated values are about 0.25 m/s, which indicates that the molten steel near the mold surface flows from the narrow wall to the SEN, so the flow pattern in the mold is the DRF. When the argon gas flow rates are 10 L/min and 14 L/min, the velocity near the mold surface is smaller than −0.2 m/s, indicating the flow pattern in the mold may be the intermediate flow (IF) or SRF. Therefore, the argon gas flow rate has a great influence on the flow pattern in the mold with a large width.  Figure 4 shows the streamlines on center plane Z = 0 of the mold, when the argon gas flow rates are 6 L/min, 10 L/min and 14 L/min, respectively, under the same conditions mentioned above. The figures below are the enlarged views corresponding to the marked area of the figures above. When the argon gas flow rate is 6 L/min, a part of the molten steel named stream A will flow up directly to the mold surface near the nozzle wall [28].  Figure 4 shows the streamlines on center plane Z = 0 of the mold, when the argon gas flow rates are 6 L/min, 10 L/min and 14 L/min, respectively, under the same conditions mentioned above. The figures below are the enlarged views corresponding to the marked area of the figures above. When the argon gas flow rate is 6 L/min, a part of the molten steel named stream A will flow up directly to the mold surface near the nozzle wall [28]. The other part molten steel stream flows to the narrow wall of the mold and then is divided into two streams when approaching the narrow wall. Stream B moves up along the narrow wall to reach the mold surface, and then flows to the nozzle. When meeting with stream A on the mold surface, it drives stream A and the SEN together, due to the fact that stream B is much stronger than stream A. Stream C moves down along the narrow wall to form a clockwise circulation flow in the lower part of the mold. Therefore, the flow pattern is DRF at the argon gas flow rate of 6 L/min, which has low risk of slag entrainment. Just as Kunstreich concluded, maintaining a stable DRF is key to eliminating slab defects [29].
When the argon gas flow rate is 10 or 14 L/min, the enlarger view below shows that stream A has occupied the whole upper circulation area, and stream B nearly disappears under the effect of the floating of the argon gas bubbles. The flow stream flowing out of the side port of the SEN no longer moves directly to the narrow wall, but instead moves directly up to the mold surface near the nozzle, thus forming an SRF. In this case, it is easy for slag entrainment to be caused due to the downward flow near the narrow wall of the mold. This results in the incidence of slag entrainment defects being greatly increased, which is consistent with the previous results [30,31]. stream A has occupied the whole upper circulation area, and stream B nearly disappears under the effect of the floating of the argon gas bubbles. The flow stream flowing out of the side port of the SEN no longer moves directly to the narrow wall, but instead moves directly up to the mold surface near the nozzle, thus forming an SRF. In this case, it is easy for slag entrainment to be caused due to the downward flow near the narrow wall of the mold. This results in the incidence of slag entrainment defects being greatly increased, which is consistent with the previous results [30,31]. The top-surface levels obtained from the surface pressure are related to the stability of the top surface. The top-surface liquid displacement (hi) was estimated from a simple potential energy balance [32]: where the P i is the pressure of one position on the top surface, P mean is the mean value of pressure across the whole top surface; ρ steel is the density of the molten steel; ρ slag is the density of the top slag and g is the gravitational acceleration rate. Figure 5 is the top surface level profile of the mold under the same conditions and it can reflect the fluctuation of the top free surface of the mold. When the argon gas flow rate is 6 L/min, the fluctuation of the whole liquid level is small, as shown in Figure 5a. As the stream moves out of the nozzle, it impacts the narrow wall and consumes part of the kinetic energy, and then enters the top free surface at a certain angle. So, a small fluctuation of the free surface will be formed between the narrow wall and the SEN. When the argon gas flow rate is 10 L/min, the fluctuation of the liquid level near the nozzle increases, as shown in Figure 5b. The high fluctuation area increases obviously, resulting in the increased risk of slag entrainment. When the argon gas flow rate is 14 L/min-because a large amount of molten steel rushes to the top free surface near the nozzle-the fluctuation of the whole free surface will inevitably deteriorate sharply, and the risk of slag entrainment will further increase, as presented in Figure 5c. Therefore, when the mold width is 1800 mm, the immersion depth is 140 mm, and the casting speed is 1.1 m/min, the argon gas flow rate should be controlled at or below 6 L/min. pressure across the whole top surface; ρsteel is the density of the molten steel; ρslag is the density of the top slag and g is the gravitational acceleration rate. Figure 5 is the top surface level profile of the mold under the same conditions and it can reflect the fluctuation of the top free surface of the mold. When the argon gas flow rate is 6 L/min, the fluctuation of the whole liquid level is small, as shown in Figure 5a. As the stream moves out of the nozzle, it impacts the narrow wall and consumes part of the kinetic energy, and then enters the top free surface at a certain angle. So, a small fluctuation of the free surface will be formed between the narrow wall and the SEN. When the argon gas flow rate is 10 L/min, the fluctuation of the liquid level near the nozzle increases, as shown in Figure 5b. The high fluctuation area increases obviously, resulting in the increased risk of slag entrainment. When the argon gas flow rate is 14 L/min-because a large amount of molten steel rushes to the top free surface near the nozzle-the fluctuation of the whole free surface will inevitably deteriorate sharply, and the risk of slag entrainment will further increase, as presented in Figure 5c. Therefore, when the mold width is 1800 mm, the immersion depth is 140 mm, and the casting speed is 1.1 m/min, the argon gas flow rate should be controlled at or below 6 L/min.

Effect of Casting Speed
The calculated velocity near the mold surface is compared with that measured at the different casting speeds in Figure 6. The experimental conditions are a mold width of 1800 mm, argon gas flow rate of 10 L/min, immersion depth of 160 mm, and casting speed of 1.1 and 1.2 m/min, respectively. The calculated results match well with the measured re-

Effect of Casting Speed
The calculated velocity near the mold surface is compared with that measured at the different casting speeds in Figure 6. The experimental conditions are a mold width of 1800 mm, argon gas flow rate of 10 L/min, immersion depth of 160 mm, and casting speed of 1.1 and 1.2 m/min, respectively. The calculated results match well with the measured results. When the casting speed is 1.1 m/min, both the calculated and measured liquid steel flow velocities are negative, which indicates that the flow direction of the molten steel near the mold surface is from the SEN to the mold narrow wall and the flow pattern may be IF or SRF. However, when the casting speed is increased to 1.2 m/min, the measured and calculated values are both positive and the flow velocity is close to 0.4 m/s. Namely, the flow direction of the molten steel near the mold surface is from the mold narrow wall to the SEN and the flow pattern is DRF. On the other hand, the velocity of 0.4 m/s is generally regarded as the upper limit velocity near the mold surface. When the velocity is larger than 0.4 m/s, the risk of shear entrainment of the mold powder will be greatly increased [33]. Therefore, the casting speed should be controlled to be smaller than 1.2 m/min. It can be seen, that the change of casting speed is only 0.1 m/min, but the flow velocity near the mold top free surface changes drastically and the flow direction is reversed sharply. Therefore, the flow field in the CC mold with a large width is very sensitive to the casting speed. The main reason is that the change in steel throughput is very large when the casting speed is changed for the CC mold with a large width. versed sharply. Therefore, the flow field in the CC mold with a large width is very sensitive to the casting speed. The main reason is that the change in steel throughput is very large when the casting speed is changed for the CC mold with a large width.   Figure 7a, it can be clearly seen that when the casting speed is 1.1 m/min, stream A-which flows directly to the mold surface near the nozzle-is stronger than stream B. Therefore, the molten steel near the mold surface at 1/4 the width of the mold flows from SEN to the narrow wall. The reason is that the jet stream from the side port is relatively weak and is lifted and dissipated by the flotation of the argon gas bubbles. Moreover, owing to the shallow immersion depth of the SEN, the upward molten steel near the port of the SEN can reach the mold surface with a small kinetic energy loss and flow to the narrow wall, thus forming an IF flow pattern between the SRF and DRF. This was discussed in Kunstreich and Liu's research, where it was stated that this flow state is unstable and the meeting point changes with time [29,34]. The risk of slag entrainment in this flow state is also high.
In contrast, when the casting speed is 1.2 m/min, the flow of molten steel from the side port is relatively strong and most of the molten steel flows directly to the narrow wall and is divided into two streams. The upward stream B moves up along the narrow wall to the mold top surface, and then flows towards the SEN. Stream C flows downwards along the narrow wall to form a counterclockwise circulation flow in the left lower part of the mold. Under the influence of the floating of the argon gas bubbles, the small stream A with a small amount of kinetic energy will move up near the SEN wall, but the flow is very weak. At the high casting speed of 1.2 m/min, stream B is much stronger than steam A and so plays the dominant role, thus forming the DRF.   Figure 7a, it can be clearly seen that when the casting speed is 1.1 m/min, stream A-which flows directly to the mold surface near the nozzle-is stronger than stream B. Therefore, the molten steel near the mold surface at 1/4 the width of the mold flows from SEN to the narrow wall. The reason is that the jet stream from the side port is relatively weak and is lifted and dissipated by the flotation of the argon gas bubbles. Moreover, owing to the shallow immersion depth of the SEN, the upward molten steel near the port of the SEN can reach the mold surface with a small kinetic energy loss and flow to the narrow wall, thus forming an IF flow pattern between the SRF and DRF. This was discussed in Kunstreich and Liu's research, where it was stated that this flow state is unstable and the meeting point changes with time [29,34]. The risk of slag entrainment in this flow state is also high.   Figure  8a. When the casting speed is 1.2 m/min-owing to the strong stream B-the molten steel can still flow to the top free surface with a large amount of kinetic energy after impacting the narrow wall, and then combines with the weak stream A to flow to the SEN together on the mold surface to form a DRF. It has also been reported that the fluctuation in the mold liquid level is relatively gentle under the DRF [35]. Therefore, when the mold width is 1800 mm, the argon gas flow rate is 10 L/min, and the immersion depth is 160 mm, the casting speed should be greater than 1.1 m/min and smaller than 1.2 m/min. In contrast, when the casting speed is 1.2 m/min, the flow of molten steel from the side port is relatively strong and most of the molten steel flows directly to the narrow wall and is divided into two streams. The upward stream B moves up along the narrow wall to the mold top surface, and then flows towards the SEN. Stream C flows downwards along the narrow wall to form a counterclockwise circulation flow in the left lower part of the mold. Under the influence of the floating of the argon gas bubbles, the small stream A with a small amount of kinetic energy will move up near the SEN wall, but the flow is very weak. At the high casting speed of 1.2 m/min, stream B is much stronger than steam A and so plays the dominant role, thus forming the DRF. Figure 8 shows the top surface level profile of the mold under the same conditions mentioned above. When the casting speed is 1.1 m/min, most of the molten steel flows upwards to the top free surface near the nozzle under the effect of the floating of the argon gas bubbles, so that the liquid level near the SEN fluctuates violently, as shown in Figure  8a. When the casting speed is 1.2 m/min-owing to the strong stream B-the molten steel can still flow to the top free surface with a large amount of kinetic energy after impacting the narrow wall, and then combines with the weak stream A to flow to the SEN together on the mold surface to form a DRF. It has also been reported that the fluctuation in the mold liquid level is relatively gentle under the DRF [35]. Therefore, when the mold width is 1800 mm, the argon gas flow rate is 10 L/min, and the immersion depth is 160 mm, the casting speed should be greater than 1.1 m/min and smaller than 1.

Effect of Immersion Depth
Under the conditions of a mold width of 1800 mm, casting speed of 1.1 m/min, and argon gas flow rate of 10 L/min, the numerical simulation results of the surface velocities are compared with the results obtained by the rod deflection method with the different SEN immersion depths in Figure 9. When the immersion depth of the SEN is 140 and 160 mm, the surface velocity of the mold is negative from the calculation and measurement results, which means that the flow direction of the molten steel near the mold surface is from the SEN to the narrow wall of the mold. So, it can be judged that the flow pattern is the SRF or IF. When the immersion depth is 180 mm, the velocity near the mold surface is positive. The flow direction of the molten steel near the mold surface is from the mold narrow wall to the SEN. Owing to the small surface velocity, the flow pattern can be judged as a weak DRF.

Effect of Immersion Depth
Under the conditions of a mold width of 1800 mm, casting speed of 1.1 m/min, and argon gas flow rate of 10 L/min, the numerical simulation results of the surface velocities are compared with the results obtained by the rod deflection method with the different SEN immersion depths in Figure 9. When the immersion depth of the SEN is 140 and 160 mm, the surface velocity of the mold is negative from the calculation and measurement results, which means that the flow direction of the molten steel near the mold surface is from the SEN to the narrow wall of the mold. So, it can be judged that the flow pattern is the SRF or IF. When the immersion depth is 180 mm, the velocity near the mold surface is positive. The flow direction of the molten steel near the mold surface is from the mold narrow wall to the SEN. Owing to the small surface velocity, the flow pattern can be judged as a weak DRF.  Figure 10 shows the streamlines on the center plane of the mold when the immersion depth is 140, 160 and 180 mm, under the above conditions. It can be seen from the streamlines that with increasing the immersion depth, the strengths of stream A and stream B change obviously. When the immersion depth is 140 mm, stream A is dominant, so the flow pattern in the mold is close to the SRF. When the immersion depth is 160 mm, the strength of stream A is almost equivalent to that of stream B, and the flow pattern in the mold should be IF. When the immersion depth is further increased to 180 mm, stream B dominates the upper circulation flow in the mold, and the flow pattern of the mold is weak DRF. The reason for this is that increasing the SEN immersion depth leads to an increase in the upwards distance of stream A, which makes it more difficult to directly impact the top free surface. Therefore, with an increasing SEN immersion depth, it is easier to form DRF, which can reduce the risk of slag entrainment.  Figure 11 shows the top surface level profiles of the mold under the same above conditions. When the casting speed is 1.1 m/min and the argon gas flow rate is 10 L/min, the kinetic energy of molten steel flowing out from side port is relatively small and the buoyancy force of the argon gas bubbles is relatively larger. When the immersion depth of the SEN is 140 mm, the climbing distance of the molten steel to the top surface is short, and  Figure 10 shows the streamlines on the center plane of the mold when the immersion depth is 140, 160 and 180 mm, under the above conditions. It can be seen from the streamlines that with increasing the immersion depth, the strengths of stream A and stream B change obviously. When the immersion depth is 140 mm, stream A is dominant, so the flow pattern in the mold is close to the SRF. When the immersion depth is 160 mm, the strength of stream A is almost equivalent to that of stream B, and the flow pattern in the mold should be IF. When the immersion depth is further increased to 180 mm, stream B dominates the upper circulation flow in the mold, and the flow pattern of the mold is weak DRF. The reason for this is that increasing the SEN immersion depth leads to an increase in the upwards distance of stream A, which makes it more difficult to directly impact the top free surface. Therefore, with an increasing SEN immersion depth, it is easier to form DRF, which can reduce the risk of slag entrainment.  Figure 10 shows the streamlines on the center plane of the mold when the immersion depth is 140, 160 and 180 mm, under the above conditions. It can be seen from the streamlines that with increasing the immersion depth, the strengths of stream A and stream B change obviously. When the immersion depth is 140 mm, stream A is dominant, so the flow pattern in the mold is close to the SRF. When the immersion depth is 160 mm, the strength of stream A is almost equivalent to that of stream B, and the flow pattern in the mold should be IF. When the immersion depth is further increased to 180 mm, stream B dominates the upper circulation flow in the mold, and the flow pattern of the mold is weak DRF. The reason for this is that increasing the SEN immersion depth leads to an increase in the upwards distance of stream A, which makes it more difficult to directly impact the top free surface. Therefore, with an increasing SEN immersion depth, it is easier to form DRF, which can reduce the risk of slag entrainment.  Figure 11 shows the top surface level profiles of the mold under the same above conditions. When the casting speed is 1.1 m/min and the argon gas flow rate is 10 L/min, the kinetic energy of molten steel flowing out from side port is relatively small and the buoyancy force of the argon gas bubbles is relatively larger. When the immersion depth of the SEN is 140 mm, the climbing distance of the molten steel to the top surface is short, and  Figure 11 shows the top surface level profiles of the mold under the same above conditions. When the casting speed is 1.1 m/min and the argon gas flow rate is 10 L/min, the kinetic energy of molten steel flowing out from side port is relatively small and the buoyancy force of the argon gas bubbles is relatively larger. When the immersion depth of the SEN is 140 mm, the climbing distance of the molten steel to the top surface is short, and the strong energy stream A impinging on the top free surface leads to a sharp fluc-tuation in the top free surface near the SEN and 1/4 width of the mold. With increasing the SEN immersion depth to 160 mm, stream A impinging on the free surface becomes weak due to large distance to the top surface. Stream B becomes stronger and tends to form IF. It can be seen in the liquid level fluctuation diagram that the change in liquid level gradually tends to be gentle, but the fluctuation near the nozzle is still intense. This is because the liquid steel can impact the top free surface under the buoyancy force of the argon gas bubbles, resulting in a sharp fluctuation in the liquid level near the nozzle. When the immersion depth of the SEN is further increased to 180 mm, more of molten steel of stream A will combine into stream B because of insufficient energy floating up. It is easier to form a DRF. Therefore, under the present process conditions, the mold flow field is the optimized when the immersion depth is 180 mm. It can be seen in the liquid level fluctuation diagram that the change in liquid level gradually tends to be gentle, but the fluctuation near the nozzle is still intense. This is because the liquid steel can impact the top free surface under the buoyancy force of the argon gas bubbles, resulting in a sharp fluctuation in the liquid level near the nozzle. When the immersion depth of the SEN is further increased to 180 mm, more of molten steel of stream A will combine into stream B because of insufficient energy floating up. It is easier to form a DRF. Therefore, under the present process conditions, the mold flow field is the optimized when the immersion depth is 180 mm.   Under the same conditions mentioned above, Figure 13 shows the streamlines on the center plane of the mold when the bottom shape is a mountain and a well. When the bottom shape is a mountain, the molten steel impinges on the bottom of the SEN without much consumption of kinetic energy due to the mountain bottom shape and ejects from the ports on both sides, and then impacts the narrow wall to split into two strong streams. A part of the flow stream moves intensively upwards along the narrow wall, enters the top free surface at a certain angle, then moves towards the SEN. Therefore, a strong DRF will be formed with the mountain bottom shaped SEN.

Effect of SEN Bottom Shape
On the other hand, when the bottom shape of the SEN is a well, the molten steel will whirl and rebound in the groove when it impacts the bottom of the nozzle, resulting in a   Figure 12 shows the comparison between the surface velocities measured by rod deflection method and numerical simulation under the conditions of the mold width of 2000 mm, the casting speed of 1.0 m/min, the argon gas flow rate of 8 L/min and the immersion depth of 150 mm. When the bottom shape of the SEN is a mountain, the measured and calculated velocities near the mold surface are about 0.4 m/s. When the bottom shape of the SEN is a well, the measured and calculated velocities are smaller than 0.25 m/s. Although both are DRF, the former velocities are much larger than the latter ones. Under the same conditions mentioned above, Figure 13 shows the streamlines on the center plane of the mold when the bottom shape is a mountain and a well. When the bottom shape is a mountain, the molten steel impinges on the bottom of the SEN without much consumption of kinetic energy due to the mountain bottom shape and ejects from the ports on both sides, and then impacts the narrow wall to split into two strong streams. A part of the flow stream moves intensively upwards along the narrow wall, enters the top free surface at a certain angle, then moves towards the SEN. Therefore, a strong DRF will be formed with the mountain bottom shaped SEN.

Effect of SEN Bottom Shape
On the other hand, when the bottom shape of the SEN is a well, the molten steel will whirl and rebound in the groove when it impacts the bottom of the nozzle, resulting in a Under the same conditions mentioned above, Figure 13 shows the streamlines on the center plane of the mold when the bottom shape is a mountain and a well. When the bottom shape is a mountain, the molten steel impinges on the bottom of the SEN without much consumption of kinetic energy due to the mountain bottom shape and ejects from the ports on both sides, and then impacts the narrow wall to split into two strong streams. A part of the flow stream moves intensively upwards along the narrow wall, enters the top free surface at a certain angle, then moves towards the SEN. Therefore, a strong DRF will be formed with the mountain bottom shaped SEN. great loss of kinetic energy of the jet stream. The strength of the flow ejecting from the ports is much weaker than that of the mountain bottom shaped SEN, so that the DRF formed is obviously weaker. Under the same conditions, Figure 14 shows the top surface level profiles of the mold with different bottom shapes of SEN. When the bottom shape is a mountain, the fluctuation of the free surface near the 1/4 width of the mold is quite intense, which is caused by the strong upward circulation flow [36]. The strong streams will eject from the ports on both sides without much consumption of kinetic energy due to mountain bottom shape. A part of flow stream moves intensively upward along the narrow wall to form the strong DRF. Under this condition, the surface velocity is as high as 0.4 m/s, so the liquid level fluctuation is large and the shear stress near the top surface is intensified, and the risk of slag entrainment is greatly increased. Lu also concluded that the strong upper circulation flow easily results in slag entrainment by the strong shear stress [31]. When the bottom shape is well, the free surface fluctuation of the mold is obviously alleviated. Under the present conditions, the SEN with a well bottom shape is much better than that with mountain shape bottom. On the other hand, when the bottom shape of the SEN is a well, the molten steel will whirl and rebound in the groove when it impacts the bottom of the nozzle, resulting in a great loss of kinetic energy of the jet stream. The strength of the flow ejecting from the ports is much weaker than that of the mountain bottom shaped SEN, so that the DRF formed is obviously weaker.
Under the same conditions, Figure 14 shows the top surface level profiles of the mold with different bottom shapes of SEN. When the bottom shape is a mountain, the fluctuation of the free surface near the 1/4 width of the mold is quite intense, which is caused by the strong upward circulation flow [36]. The strong streams will eject from the ports on both sides without much consumption of kinetic energy due to mountain bottom shape. A part of flow stream moves intensively upward along the narrow wall to form the strong DRF. Under this condition, the surface velocity is as high as 0.4 m/s, so the liquid level fluctuation is large and the shear stress near the top surface is intensified, and the risk of slag entrainment is greatly increased. Lu also concluded that the strong upper circulation flow easily results in slag entrainment by the strong shear stress [31]. When the bottom shape is well, the free surface fluctuation of the mold is obviously alleviated. Under the present conditions, the SEN with a well bottom shape is much better than that with mountain shape bottom.

Conclusions
In this paper, the rod deflection method was applied to quantitatively measure the velocity near the mold surface at high temperatures under the different casting speeds: argon gas flow rates; the submerged entry nozzle (SEN) immersion depths; and bottom shapes of SEN. The k-ε model coupled with discrete phase model was utilized to simulate the flow field in the mold. The main conclusions are as follows.
(1) Under the conditions of the large mold width of 1800 mm, 1.1 m/min casting speed and 140 mm SEN immersion depth, the velocity near the mold surface decreases with increasing the argon gas flow rate. The calculated results match very well with the measured results. When the argon gas flow rate is 6 L/min, the flow pattern is the double roll flow (DRF), and the liquid level fluctuation is small. When the argon gas flow rate is increased to 10 and 14 L/min, the flow pattern is the single roll flow (SRF), the liquid level fluctuation is large, and the risk of slag entrainment increases. (2) When the mold width is 1800 mm, the argon gas flow rate is 10 L/min, and the immersion depth is 160 mm, the velocity near the mold surface sensitively increases with increasing casting speed. The calculated results are consistent with the measured ones. When the casting speed is 1.1 m/min, the direction of the velocity near the mold surface is from SEN to narrow wall, forming an intermediate flow (IF). The

Conclusions
In this paper, the rod deflection method was applied to quantitatively measure the velocity near the mold surface at high temperatures under the different casting speeds: argon gas flow rates; the submerged entry nozzle (SEN) immersion depths; and bottom shapes of SEN. The k-ε model coupled with discrete phase model was utilized to simulate the flow field in the mold. The main conclusions are as follows.
(1) Under the conditions of the large mold width of 1800 mm, 1.1 m/min casting speed and 140 mm SEN immersion depth, the velocity near the mold surface decreases with increasing the argon gas flow rate. The calculated results match very well with the measured results. When the argon gas flow rate is 6 L/min, the flow pattern is the double roll flow (DRF), and the liquid level fluctuation is small. When the argon gas flow rate is increased to 10 and 14 L/min, the flow pattern is the single roll flow (SRF), the liquid level fluctuation is large, and the risk of slag entrainment increases. (2) When the mold width is 1800 mm, the argon gas flow rate is 10 L/min, and the immersion depth is 160 mm, the velocity near the mold surface sensitively increases with increasing casting speed. The calculated results are consistent with the measured ones. When the casting speed is 1.1 m/min, the direction of the velocity near the mold surface is from SEN to narrow wall, forming an intermediate flow (IF). The mold surface fluctuation is obviously intensified, which is easy to cause slag entrainment defects. When the casting speed is only increased to 1.2 m/min, the velocity near the mold surface changes drastically and is close to 0.4 m/s, which is the upper limit velocity near the mold surface. Therefore, the optimized casting speed should be larger than 1.1 m/min and smaller than 1.2 m/min. (3) When the mold width is 1800 mm, the casting speed is 1.1 m/min, and the argon gas flow rate is 10 L/min, the velocity near the mold surface is obviously increased with the increasing immersion depth. The calculated results are also in accord with the measured ones. When the immersion depth of the nozzle is 140 and 160 mm, the velocity near the mold surface is from SEN to narrow wall, the flow pattern is SRF or IF, and the liquid level fluctuates greatly. When the immersion depth of the SEN is 180 mm, the velocity near the mold surface is from narrow wall to the SEN, the flow pattern is DRF, and the fluctuation of the liquid level is small, which is conducive to control the risk of slag entrainment.