Large Eddy Simulation of Multi-Phase Flow and Slag Entrapment in Molds with Different Widths

: Slag entrapment is a critical problem that affects the quality of steel. In this work, a three-dimensional model is established to simulate the slag entrapment phenomenon, mainly focusing on the slag entrapment phenomenon at the interface between slag and steel in molds with different widths. The large eddy simulation (LES) model and discrete particle model (DPM) are used to simulate the movements of bubbles. The interactions between phases involve two-way coupling. The accuracy of our mathematical model is validated by comparing slag–metal interface ﬂuctuations with practical measurements. The results reveal that the average interface velocity and transverse velocity decrease as the mold width increases, however, they cannot represent the severity of slag entrapment at the interface between slag and steel. Due to the inﬂuence of bubble motion behavior, the maximum interface velocity increases with mold width and causes slag entrapment readily, which can reﬂect the severity of slag entrapment. On this basis, by monitoring the change of impact depths in different molds, a new dimensionless number “C” is found to reveal the severity of slag entrapment at the interface between slag and steel. The results show that the criterion number C increases with mold width, which is consistent with the results of ﬂaw detection. Therefore, criterion number C can be used to reﬂect the severity of slag entrapment in different molds.


Introduction
The mold is an important part of steel metallurgy and can be called the "heart" of a continuous casting machine. During this process, the effect of powder added into the mold can be summarized by three aspects: (1) preventing oxidation of molten steel; (2) absorbing non-metallic inclusions; (3) filling the gap between the mold and the slab surface to improve heat transfer and lubricating the mold surface. However, due to the instability of level fluctuation in the mold, the liquid slag is often incorporated into the molten steel, damaging mechanical properties of steel products, especially for automobile steel. Therefore, slag entrapment has become a critical problem that requires serious concern.
Due to "black box" operation, the slag entrapment cannot be seen during the real casting process. Therefore, visual models have been developed to reproduce the slag entrapment at the interface between slag and steel. Among these studies, Gguta et al. [1] found an apparent asymmetric flow near the nozzle. Afterwards, Li et al. [2] captured the asymmetric vortex distribution of the steel/slag interface by injecting black sesames into water, and developed a mathematical model to predict the asymmetry flow in the mold. The study of Iguchi et al. [3] showed that the slag entrapment caused by shearing strength is the main mechanism for the slag entrapment during high casting speed operation, and that the interfacial tension has a significant influence on it. Watanabe and Yamashita Metals 2021, 11, 374 2 of 13 et al. [4,5] studied slag entrapment by argon blowing, finding that the maximum depth of slag involved would not exceed three times its diameter. In addition, according to Savolainen et al. [6], the effect of the viscosity of slag on the formation and size of slag droplets should be paid attention to take full control of it. Yamada et al. [7] suggested that argon bubbles in a mold become the desirable sites where alumina inclusions are gathered and form large alumina clusters. Despite these studies, some researchers focused on the influence of liquid properties on slag entrapment in the mold, such as liquid metal density and interfacial tension within the slag-metal interface [8][9][10], in which the critical slag entrapment velocity can be determined. However, this velocity can only predict the slag entrapment caused by shearing stress. Lei Hong et al. [11] proposed another theoretical equation for calculating shear entrapment, taking the viscosity of slag into account. However, Chung and Cramb [12] and others believe that the interfacial tension coefficient should be reduced to about three percent of the original value due to the existence of interfacial reactions. Harman [13] considered nine factors and obtained another formula for calculating critical velocity through non-linear fitting. The above studies are of great significance for understanding the slag entrapment phenomenon. However, the physical models cannot meet all the similarity criteria at the same time, which leads to some limitations for slag entrapment results.
With the rapid development of computer science, computational fluid dynamics (CFD) technology has been an important tool in metallurgical process research, and its advantages are increasingly prominent. Saeedipour et al. [14] established a three-phase mathematical model to study the interface wave problem. Liu et al. [15] established a quasi-four-phase model to study the effect of bubbles on the interface fluctuation of the slag-metal interface. Li et al. [16] analyzed three kinds of slag entrapment mechanisms, and gave the transient process of mold slag entrapment in molten steel. Although these studies are of great significance to the study of slag entrapment, the relationship between mold impact depth and velocity, especially the association with slag entrapment, was not explored. There is still a lack of effective evaluation criteria to predict the effect of slag entrapment.
In this work, the movement behavior of mold slag in the mold is studied to reveal the influence of mold structure on slag entrapment and a theoretical foundation for mold width adjustment is laid. The innovations of this paper are composed of three parts: Firstly, the slag entrapment between different molds are elaborated in detail. (2) Secondly, a mathematical model of four-phase (slag-metal-gas-air) flow is established to explain differences in slag entrapment in (1). Thirdly, a new dimensionless value is established to characterize the severity of slag entrapment in molds. The research results can act as a guide for the continuous casting process.

Basic Assumptions
To simplify the calculation, the mathematical model used in this work is based on the following assumptions: (i) Liquid steel is regarded as a Newtonian fluid, and its basic parameters, such as density and viscosity, are considered as constants.
(ii) The heat transfer and solidification process between molten steel and cooling water are not considered, and the thermal characteristics of slag are ignored.
(iii) The discrete phase bubble is assumed to be spherical, and its size change is ignored in the process of floating.
(iv) The taper of the mold, as well as heat transfer between the mold and slab, are all ignored.

Governing Equations
There are four phases existing in the mold: molten steel, liquid slag, air, and argon. In this work, we use the volume of fluid (VOF) method to describe interactions between steel, slag, and air and adopt a discrete particle model (DPM) to track trajectories of argon bubbles. The interactions between these phases involve two-way coupling, following the law of Newton. The continuity equation for these phases can be written as follows: where ρ k is the density of the phase, the foot mark k representing the phases u m is the velocity of the mixture phase, α k is the volume fraction of each phase, which fulfills the equation α l + α s + α g = 1.
Large eddy simulation (LES) is used to solve the Navier-Stokes (N-S) equations of fluid flow: In this equation, fluid motion is solved by a set of equations. The term ρ m represents the density of the mixture phase, P is static pressure. The molecular viscosity µ m and turbulent viscosity µ t in the equation are all weighted average values of the volume fraction of each phase, and the effective viscosity µ effect = µ m + µ t . F is the forces acting on bubbles. The term F γ is the interface tension between these phases.
The expression of subgrid-scale stress τ ij is as follows: where the term τ kk is the isotropic part of the subgrid scale, δ ij is the Kronecker symbol, S ij is the strain rate, turbulent viscosity µ t = ρ m L 2 s S . The calculation formula of mixing length is as follows: where κ is the von Karman constant, taken as 0.4, d is the vertical distance from the fluid to the wall, C s is the Smagorinsky constant, taken as 0.2.

Discrete Particle Model (DPM)
The argon bubbles injected from the nozzle easily float up and escape from the upper surface of the mold. In this process, the momentum exchanges between the argon and steel are treated as two-way coupling, which obeys the second law of Newton: where u p and m p are the velocity and mass of particles, respectively, and F is the resultant force acting on the bubbles, the expression of which can be written as follows: where the terms on the right side of the formula are gravity, buoyancy force, pressure gradient force, drag force, lift force, and virtual mass force, respectively. The calculation equations [17][18][19][20] can be found in Table 1. The bubble size ranges from 0.5-15 mm, following the Rosin-Rammler law: where the variable Y d is the mass fraction of bubbles whose diameters are greater than d.
The averaged bubble diameter d m = 5 mm, and the spread parameter n = 2.

Source Term Formula Annotations
Buoyancy force plus gravity force, The net effect acts on the difference between particle and fluid densities. The variable g is gravity acceleration, d p is particle diameter.
Where drag coefficient ,and u p is particle velocity, Re p is particle Reynolds number. Pressure gradient force F p ρ m /ρ p u p ∇u m Pressure gradient force is significant when ρ/ρ p ≥ 0.1.
Where the virtual mass force coefficient C v−m = 0.5.

Boundary Conditions and Numerical Details
As shown in Figure 1, the whole physical model consists of four parts: nozzle, mold, foot roller zone, and secondary cooling zone. In the production process, the thicknesses of argon and the slag layer are 40 mm and 50 mm, respectively. After entering the mold, the argon gas expands rapidly, resulting in a decrease in the density of argon gas. The bubble density in the molten steel can be calculated through the ideal gas law, as shown in Equation (8). The steel velocity from the nozzle is calculated by the casting speed, and a no-slip condition is used to perform the wall treatment of the mold. In this work, three typical widths of the mold were taken to study the slag entrapment inside the mold, and the cross section are 250 mm × 1100 mm, 250 mm × 1400 mm, 250 mm × 1650 mm, respectively.
where the argon density at 20 • C is 1.78 kg/m 3 . The pressure P is standard atmospheric pressure, which can be assumed to be equal to standard atmosphere pressure near the top of the mold. In order to obtain a fine vortex structure, the mesh is refined near the slag-metal interface and gas-slag interface. The whole region contains 2.1 × 10 6 structural grids with good independence. The calculation step is set to 0.01 s and the total calculation time is 100 s. To save calculation time and cost, the k-ε model is used to calculate the steady-state  In order to obtain a fine vortex structure, the mesh is refined near the slag-metal interface and gas-slag interface. The whole region contains 2.1 × 10 6 structural grids with good independence. The calculation step is set to 0.01 s and the total calculation time is 100 s. To save calculation time and cost, the k-ε model is used to calculate the steady-state field, and then switched to the LES model to simulate the transient flow field. The specific model parameters are shown in Table 2. The movements of argon bubbles consist three key stages, as shown in Figure 2: (1) Entering the mold from the nozzle and being carried downward by the nozzle jet; (2) floating up into the molten steel and passing through the slag-metal interface and entering into the slag; (3) collapsing and subsequently disappearing after floating near the gas/slag interface. The movements of bubbles in the continuous phases and the escape near the gas/slag interface are achieved by coding user-defined functions (UDFs). In addition, the impact point is defined as the position where the vertical shear force equals zero on the center plane of the mold.  Figure 3 shows the surface defects after the rolling process, which is caused by the large inclusions. It is clearly seen that many black lines on the surface of the slabs are stretching across the whole surface of the slab, the lengths of which are more than 90 millimeters. This can degrade the quality of steel significantly. The "black line" problem is a significant problem that constrains the production of automobile steel.  Figure 3 shows the surface defects after the rolling process, which is caused by the large inclusions. It is clearly seen that many black lines on the surface of the slabs are stretching across the whole surface of the slab, the lengths of which are more than 90 millimeters. This can degrade the quality of steel significantly. The "black line" problem is a significant problem that constrains the production of automobile steel. Figure 3 shows the surface defects after the rolling process, which is caused by the large inclusions. It is clearly seen that many black lines on the surface of the slabs are stretching across the whole surface of the slab, the lengths of which are more than 90 millimeters. This can degrade the quality of steel significantly. The "black line" problem is a significant problem that constrains the production of automobile steel.  Figure 4 shows the compositions of inclusions inside the black line characterized by scanning electron microscopy (SEM) and energy dispersive spectroscopy (EDS). The results show that Na, Al, Si, Ca, O, and other elements exist in the black lines, whose compositions are similar to that of the liquid slag in the mold. Therefore, it is certain that the black line defect of hot-rolled sheets is caused by slag entrapment.

Analysis of Inclusions on Slab
(a)  show that Na, Al, Si, Ca, O, and other elements exist in the black lines, whose compositions are similar to that of the liquid slag in the mold. Therefore, it is certain that the black line defect of hot-rolled sheets is caused by slag entrapment. Figure 3 shows the surface defects after the rolling process, which is caused by the large inclusions. It is clearly seen that many black lines on the surface of the slabs are stretching across the whole surface of the slab, the lengths of which are more than 90 millimeters. This can degrade the quality of steel significantly. The "black line" problem is a significant problem that constrains the production of automobile steel.    Figure 5 shows the slag entrapment ratio in slabs with different widths. The slabs are tested after the rolling process. The ratio is defined as the percentage of slabs with slag entrapment problems against totally tested slabs. It can be seen from Figure 5 that the larger the width is, the higher the slag entrapment ratio will be. There must be some special reasons behind this phenomenon, however, few works have reported this problem and tendency. So, in this work, the reason for this phenomenon will be explained in detail, through our mathematical model.   Figure 5 shows the slag entrapment ratio in slabs with different widths. The slabs are tested after the rolling process. The ratio is defined as the percentage of slabs with slag entrapment problems against totally tested slabs. It can be seen from Figure 5 that the larger the width is, the higher the slag entrapment ratio will be. There must be some special reasons behind this phenomenon, however, few works have reported this problem and tendency. So, in this work, the reason for this phenomenon will be explained in detail, through our mathematical model. Figure 5 shows the slag entrapment ratio in slabs with different widths. The slabs are tested after the rolling process. The ratio is defined as the percentage of slabs with slag entrapment problems against totally tested slabs. It can be seen from Figure 5 that the larger the width is, the higher the slag entrapment ratio will be. There must be some special reasons behind this phenomenon, however, few works have reported this problem and tendency. So, in this work, the reason for this phenomenon will be explained in detail, through our mathematical model.

Level Height of Slag-Metal Interface
Firstly, the accuracy of the mathematical model is verified through interface-level detections obtained from an eddy current sensor. The current sensor is installed on the top of the mold (Figure 6a), which can reveal the distance between the mold top and slag-metal interface. The installation position is in the middle of the mold, and the distance between the two probes is 0.4 m, as shown in Figure 6a. It can be seen from Figure  6b that the actual distance between the mold top and slag-metal interface is 83-93 mm. By comparison, the simulation results show that the value is 97-100 mm. It is easily seen that the simulated result is a little higher than the experimental result, the reason for which is the fact that the height of slag-metal interface is enhanced by bulging. Generally, the simulated results agree well with the experimental results, implying that the mathematical model established in this work is reliable (Figure 6b).

Level Height of Slag-Metal Interface
Firstly, the accuracy of the mathematical model is verified through interface-level detections obtained from an eddy current sensor. The current sensor is installed on the top of the mold (Figure 6a), which can reveal the distance between the mold top and slag-metal interface. The installation position is in the middle of the mold, and the distance between the two probes is 0.4 m, as shown in Figure 6a. It can be seen from Figure 6b that the actual distance between the mold top and slag-metal interface is 83-93 mm. By comparison, the simulation results show that the value is 97-100 mm. It is easily seen that the simulated result is a little higher than the experimental result, the reason for which is the fact that the height of slag-metal interface is enhanced by bulging. Generally, the simulated results agree well with the experimental results, implying that the mathematical model established in this work is reliable (Figure 6b).     Figure 7 shows the transient molten steel flow patterns of different cross sections. It can be seen that the flow field in the mold is asymmetric and unstable, with multiple vortices of different scales.  Figure 7 shows the transient molten steel flow patterns of different cross sections. It can be seen that the flow field in the mold is asymmetric and unstable, with multiple vortices of different scales. Generally, the liquid steel that flows out of the nozzle and impinges on the narrow surface of mold may be divided into two streams: upper recirculation flow and lower recirculation flow. Obviously, increasing the mold width leads to smaller impacting velocity and deeper injection of liquid steel. Furthermore, the maximum velocity of the slag-metal interface also increases, as shown in Figure 8.  Generally, the liquid steel that flows out of the nozzle and impinges on the narrow surface of mold may be divided into two streams: upper recirculation flow and lower recirculation flow. Obviously, increasing the mold width leads to smaller impacting velocity and deeper injection of liquid steel. Furthermore, the maximum velocity of the slag-metal interface also increases, as shown in Figure 8.  Figure 7 shows the transient molten steel flow patterns of different cross sections. It can be seen that the flow field in the mold is asymmetric and unstable, with multiple vortices of different scales. Generally, the liquid steel that flows out of the nozzle and impinges on the narrow surface of mold may be divided into two streams: upper recirculation flow and lower recirculation flow. Obviously, increasing the mold width leads to smaller impacting velocity and deeper injection of liquid steel. Furthermore, the maximum velocity of the slag-metal interface also increases, as shown in Figure 8.  The phenomenon can be attributed to the distribution of argon and the velocity field, as shown in Figure 9. The argon bubbles easily float up with a relatively high speed, even higher than the injection velocity of molten steel. When it floats to the gas-slag interface, the bubble collapses and pushes slag firmly into the steel, resulting in an increase in the interface velocity of the slag-metal interface. Therefore, the movements of argon should be strictly controlled in the continuous casting process. Figure 10 shows the impact depth and impact velocity in molds with different widths. Here, the impact depth is defined as the distance between the slag-metal interface and impact point; and the impact velocity is defined as the velocity of the impact point. It can be seen that the flowing strand develops more fully and the impact depth becomes larger with the increasing of mold width. However, due to the flow loss in the flow process, the larger the mold width is, the smaller the impact velocity will be.

Transient Flow Spectrum Analysis of Mold
The phenomenon can be attributed to the distribution of argon and the velocity field, as shown in Figure 9. The argon bubbles easily float up with a relatively high speed, even higher than the injection velocity of molten steel. When it floats to the gas-slag interface, the bubble collapses and pushes slag firmly into the steel, resulting in an increase in the interface velocity of the slag-metal interface. Therefore, the movements of argon should be strictly controlled in the continuous casting process.  Figure 10 shows the impact depth and impact velocity in molds with different widths. Here, the impact depth is defined as the distance between the slag-metal interface and impact point; and the impact velocity is defined as the velocity of the impact point. It can be seen that the flowing strand develops more fully and the impact depth becomes larger with the increasing of mold width. However, due to the flow loss in the flow process, the larger the mold width is, the smaller the impact velocity will be.

Analysis of Velocity Characteristics of Slag-Metal Interface
The average velocity fluctuation of the slag-metal interface is monitored and the results are shown in Figure 11. It can be seen that the fluctuations are quite different in different molds. On the one hand, the average velocity is 0.0732 m·s −1 with a small width of the mold and decreases to 0.0693 m·s −1 with a medium width of the mold. On the other hand, the fluctuation of the velocity with a large width of the mold is more significant, and the velocity continues to decrease with an average value of 0.0647 m·s −1 . This phenomenon is mainly attributed to the loss of molten steel flow. In a word, with the increase in the mold section, the velocity of the slag-metal interface decreases. However, there is little difference in the interface velocity between slag and liquid metal in these three kinds of molds, which cannot explain the why the slag entrapment increases with dif-   Figure 10 shows the impact depth and impact velocity in molds with different widths. Here, the impact depth is defined as the distance between the slag-metal interface and impact point; and the impact velocity is defined as the velocity of the impact point. It can be seen that the flowing strand develops more fully and the impact depth becomes larger with the increasing of mold width. However, due to the flow loss in the flow process, the larger the mold width is, the smaller the impact velocity will be.

Analysis of Velocity Characteristics of Slag-Metal Interface
The average velocity fluctuation of the slag-metal interface is monitored and the results are shown in Figure 11. It can be seen that the fluctuations are quite different in different molds. On the one hand, the average velocity is 0.0732 m·s −1 with a small width of the mold and decreases to 0.0693 m·s −1 with a medium width of the mold. On the other hand, the fluctuation of the velocity with a large width of the mold is more significant, and the velocity continues to decrease with an average value of 0.0647 m·s −1 . This phenomenon is mainly attributed to the loss of molten steel flow. In a word, with the increase in the mold section, the velocity of the slag-metal interface decreases. However, there is little difference in the interface velocity between slag and liquid metal in these three kinds of molds, which cannot explain the why the slag entrapment increases with dif-

Analysis of Velocity Characteristics of Slag-Metal Interface
The average velocity fluctuation of the slag-metal interface is monitored and the results are shown in Figure 11. It can be seen that the fluctuations are quite different in different molds. On the one hand, the average velocity is 0.0732 m·s −1 with a small width of the mold and decreases to 0.0693 m·s −1 with a medium width of the mold. On the other hand, the fluctuation of the velocity with a large width of the mold is more significant, and the velocity continues to decrease with an average value of 0.0647 m·s −1 . This phenomenon is mainly attributed to the loss of molten steel flow. In a word, with the increase in the mold section, the velocity of the slag-metal interface decreases. However, there is little difference in the interface velocity between slag and liquid metal in these three kinds of molds, which cannot explain the why the slag entrapment increases with different widths. Therefore, the average velocity of the slag-metal interface cannot be used to evaluate slag entrapment severity.
Metals 2021, 10, x FOR PEER REVIEW 10 of 13 ferent widths. Therefore, the average velocity of the slag-metal interface cannot be used to evaluate slag entrapment severity. The horizontal velocity fluctuation of the slag-metal interface is shown in Figure 12. The velocity is signed: The negative value of velocity points to the narrow side of the mold, and the positive value points to the direction of the submerged nozzle. In the mold The horizontal velocity fluctuation of the slag-metal interface is shown in Figure 12. The velocity is signed: The negative value of velocity points to the narrow side of the mold, and the positive value points to the direction of the submerged nozzle. In the mold with a small width, most of the velocity value is positive, indicating that the up-flow is strong, corresponding to a flow pattern of double-roll flow (see Figure 7a). However, with the medium width of the mold, the value is mostly negative, indicating that the flow has been transformed into single-roll flow and the flow rate is reduced. When the width of the mold reaches 1650 mm, the flow is completely transformed into single-roll flow, while the flow rate increases. Therefore, the horizontal velocity cannot explain the reasons for the differences of slag entrapment as well. The horizontal velocity fluctuation of the slag-metal interface is shown in Figure 12. The velocity is signed: The negative value of velocity points to the narrow side of the mold, and the positive value points to the direction of the submerged nozzle. In the mold with a small width, most of the velocity value is positive, indicating that the up-flow is strong, corresponding to a flow pattern of double-roll flow (see Figure 7a). However, with the medium width of the mold, the value is mostly negative, indicating that the flow has been transformed into single-roll flow and the flow rate is reduced. When the width of the mold reaches 1650 mm, the flow is completely transformed into single-roll flow, while the flow rate increases. Therefore, the horizontal velocity cannot explain the reasons for the differences of slag entrapment as well.  Figure 13 shows the variation of the maximum interface velocity with time. It can be seen that with the rise in mold width, the maximum velocity of the slag-metal interface also increases, so the probability of slag entrapment increases. This explains the phenomenon that the number of slag droplets entrapped in the mold increases with the increase in the cross section width of the mold, as shown in Figure 3. In addition, it can be found that with the increase in cross section width, the amplitude of velocity fluctuation also increases, that is, a large cross section width easily causes slag entrapment, which is consistent with the changing trend of slag drop number (see Figure 10). Therefore, the maximum interface velocity can reflect the severity of slag entrapment in a different mold, which should be paid more attention to control slag entrapment in steel.  Figure 13 shows the variation of the maximum interface velocity with time. It can be seen that with the rise in mold width, the maximum velocity of the slag-metal interface also increases, so the probability of slag entrapment increases. This explains the phenomenon that the number of slag droplets entrapped in the mold increases with the increase in the cross section width of the mold, as shown in Figure 3. In addition, it can be found that with the increase in cross section width, the amplitude of velocity fluctuation also increases, that is, a large cross section width easily causes slag entrapment, which is consistent with the changing trend of slag drop number (see Figure 10). Therefore, the maximum interface velocity can reflect the severity of slag entrapment in a different mold, which should be paid more attention to control slag entrapment in steel.  Figure 14 reveals the number of slag drops in molds of different widths. It can be seen that the number of slag drops increases with the mold width. The reason for this phenomenon is related to argon blowing, that is, a more considerable amount of argon gas retained in the wide-faced mold and a longer floating time lead to more argon bubbles retained in the slag layer, accompanied by the enhanced emulsification effect. The velocity of bubbles is so high that they may significantly increase the velocity of the steel around them. Therefore, a large mold is more likely to cause slag entrapment.  Figure 14 reveals the number of slag drops in molds of different widths. It can be seen that the number of slag drops increases with the mold width. The reason for this phenomenon is related to argon blowing, that is, a more considerable amount of argon gas retained in the wide-faced mold and a longer floating time lead to more argon bubbles retained in the slag layer, accompanied by the enhanced emulsification effect. The velocity of bubbles is so high that they may significantly increase the velocity of the steel around them. Therefore, a large mold is more likely to cause slag entrapment. Figure 14 reveals the number of slag drops in molds of different widths. It can be seen that the number of slag drops increases with the mold width. The reason for this phenomenon is related to argon blowing, that is, a more considerable amount of argon gas retained in the wide-faced mold and a longer floating time lead to more argon bubbles retained in the slag layer, accompanied by the enhanced emulsification effect. The velocity of bubbles is so high that they may significantly increase the velocity of the steel around them. Therefore, a large mold is more likely to cause slag entrapment.

Slag Entrapment Evaluation Criteria
The floatation of bubbles significantly affects the slag-metal interface, and increases the slag entrapment when mold width gets larger. Simultaneously, the floatation of bubbles can also lift the molten steel that is injected from submerged entry nozzle (SEN). Therefore, a dimensionless criterion number can be defined to reflect the effect of jet rigidity on slag entrapment in steel. As shown in Equation (9)

Slag Entrapment Evaluation Criteria
The floatation of bubbles significantly affects the slag-metal interface, and increases the slag entrapment when mold width gets larger. Simultaneously, the floatation of bubbles can also lift the molten steel that is injected from submerged entry nozzle (SEN). Therefore, a dimensionless criterion number can be defined to reflect the effect of jet rigidity on slag entrapment in steel. As shown in Equation (9): where α is the angle of the jet and θ is the angle of the nozzle (see Figure 2). Therefore, the physical meaning of Equation (9) can be interpreted as the ratio of jet angle to the original nozzle angle. Without bubble injection, the rigid jet is equal to C = 1. Considering the effect of argon, C < 1 when argon floats. Figure 15 shows variations of C value with different mold widths. It can be seen from Figure 15 that the C value increases with mold width. This trend is similar to that of slag drop variations shown in Figure 14, indicating that the C value can be used as a guide to reflect the slag entrapment in molds with different widths. where α is the angle of the jet and θ is the angle of the nozzle (see Figure 2). Therefore, the physical meaning of Equation (9) can be interpreted as the ratio of jet angle to the original nozzle angle. Without bubble injection, the rigid jet is equal to C = 1. Considering the effect of argon, C < 1 when argon floats. Figure 15 shows variations of C value with different mold widths. It can be seen from Figure 15 that the C value increases with mold width. This trend is similar to that of slag drop variations shown in Figure 14, indicating that the C value can be used as a guide to reflect the slag entrapment in molds with different widths. Figure 15. Prediction of slag severity through dimensionless value.

Conclusions
In this work, a three-dimensional model is established to simulate the slag entrapment phenomenon, mainly focusing on the slag entrapment phenomenon in different molds. The large eddy simulation (LES) model is applied for the calculation of the turbulence of molten steel, and the Smagorinsky-Lilly model is used to describe the

Conclusions
In this work, a three-dimensional model is established to simulate the slag entrapment phenomenon, mainly focusing on the slag entrapment phenomenon in different molds. The large eddy simulation (LES) model is applied for the calculation of the turbulence of molten steel, and the Smagorinsky-Lilly model is used to describe the sub-grid scale vortices. Based on the obtained study results, the following conclusions can be drawn: (1) The amount of slag entrapment increases with mold width, which is mainly due to the bubble pulsion. There exists a transformation of flow pattern in the mold when mold width increases. The double-roll flow pattern produces less slag entrapment than a single-roll flow pattern.
(2) The impact depth of molten steel decreases with mold width, while the impact velocity increases with mold width. The larger the mold width, the weaker the rigidity of the jet. Thus, the bubble injection significantly affects the flow field in the mold.
(3) Both of the average velocity and horizontal velocity of slag/metal interface cannot reflect the severity of slag entrapment in the mold. By comparison, the maximum velocity at the interface shows good advantages in predicting the severity of slag entrapment in the mold.
(4) A dimensionless criterion number C with the physical meaning of the ratio of jet angle to nozzle angle is successfully established to predict slag entrapment in different molds.