Nowadays, around 95% of the world raw steel is cast by means of the continuous casting process  . The quality of the steel slabs is greatly influenced by the fluid flow in the continuous casting mold  , and the geometry of the Submerged Entry Nozzle (SEN) is the most important parameter   . In the continuous casting mold, mold powder is added in order to prevent heat losses, provide good lubrication, and protect the molten steel from air oxidation. However, sometimes the mold powder is entrapped into the molten steel, and this entrapment deteriorates the quality of the solid steel  .
Formation of Karman’s vortexes on the meniscus of the mold is promoted by improper designs of the SEN. These vortexes cause powder entrapment and macro inclusions in the steel products  -  . The design of the SEN aims to minimize instabilities of the free surface to prevent the Karman’s vortexes. An elliptic SEN is proposed in  . In accordance to their authors, this design reduces the pressure difference around the nozzle, avoiding powder entrapment. In  the transient flow inside a bifurcated SEN is physically and numerically studied. It is reported that the flow pattern inside the SEN is periodic in nature. The multiple frequencies are determined from the power spectrum of the numerically generated time series. In  , the effect of the configuration parameters of the SEN in the flow field in a slab mold is analyzed by means of water modeling. The design of a swirling-flow SEN is proposed in  . Its performance is evaluated through water modeling. This SEN supposedly allows an increase of the casting speed without affecting the surface quality of solidified steel. In  it is analyzed the optimization of the molten steel flow in the mold by changes in the port design of the SEN.
An anchor-shaped design for a SEN for the slab continuous casting of steel is presented here. Its performance is evaluated by means of mathematical modeling. Transient 3D numerical simulations are carried out using the Computational Fluid Dynamics (CFD) technique. The performance of the proposed SEN is compared with that of a conventional cylindrical SEN.
2. Design of the Anchor-Shaped SEN
The geometry of the considered two-port anchor-shaped SEN is shown in Figure 1. The SEN dimensions are as follows: height = 0.5 m, arm-span = 0.40 m, plate thickness at arm tips = 0.02 m, inlet diameter = 0.06 m, exit port diameter = 0.04 m, wall thickness =
Figure 1. The conventional cylindrical SEN (left) and the anchor-shaped SEN (right).
0.01 m, pool height = 0.05 m, discharge angle = 12 degrees to the horizontal.
3. Mathematical Model
The flow of an isothermal incompressible Newtonian fluid and the mass conservation are represented by the Navier-Stokes equations and the continuity equation  . Both equations come from force balances and mass conservation, respectively:
where ρ is the fluid density; ui is the ith component of the fluid velocity u; t is time; xj is j spatial coordinate; p is pressure; and μeff is the effective fluid viscosity. To maintain the mass balance in the system, the continuity equation must be solved. Turbulence is simulated by means of the classical two equations K-ε model  :
The effective viscosity μeff of Equation (1) is determined from the expression
, where μ0 is the laminar viscosity and μt is the turbulent viscosity. On the other hand, μt is obtained from the expression, where K and ε are calculated from Equation (2) and Equation (3). Besides, the values of σK, σε, C1, C2 and Cμ are 1.0, 1.3, 1.44, 1.92 and 0.09, respectively  . Boundary conditions for K and ε at the inlet nozzle is calculated from and  where Uin and Dn are the inlet nominal velocity and the nozzle diameter, respectively. The Pressure-Implicit with Splitting of Operators (PISO) algorithm was employed for the pressure-velocity coupling.
The Volume of Fluid (VOF) model is employed to issue the multiphase flow. VOF is based on the assumption that two or more phases are not interpenetrating  . For each additional phase q its volume fraction αq is introduced as a variable. In each control volume the volume fractions of all phases sum to unity. The tracking of the interface between the phases is accomplished by solving the continuity equation for each phase   :
4. Numerical Solution
To analyze the performance of the conventional and anchor-shaped SEN, transient 3-D two-phase (air, molten steel) isothermal computer simulations were carried out using a time step of 0.001 s. This time step value was determined from numerical experiments, and corresponds to the maximum value which still provides numerical stability. Computer runs of 90 s of integration time were considered to obtain a well developed molten steel flow in the mold. The mold dimensions were as follows: height, 2 m; width, 0.6 m; thickness, 0.2 m. The initial level of molten steel was 1.82 m, the remaining volume of the mold was air. The physical properties of molten steel were as follows: density, 7100 kg/m3; viscosity, 0.0067 kg/(m.s). Boundary conditions were as follows: velocity inlet 2.83 m/s, which corresponds to a casting speed of 1.5 m/min; turbulent kinetic energy 0.08 m2/s2, turbulent dissipation rate 0.755 m2/s3. The computational mesh employed in the computer simulations is shown in Figure 2. This mesh contains 464,000 tetrahedral/hybrid cells.
5. Results and Discussion
Figure 3 shows the phase distribution and the velocity vectors for the conventional cylindrical SEN (left) and the anchor-shaped SEN (right) for the vertical plane with Y = 0. An elapsed time of 90 s is chosen. The cylindrical SEN exhibits a deeper depression of the free surface and larger magnitude of the velocity vectors in the upper section of the mold than the anchor-shaped SEN. This means that the cylindrical SEN is prone to generate Karman’s vortexes on the meniscus of the mold. Therefore, the chance of powder entrapment and capturing of macro inclusions is greater for the conventional SEN than for the anchor-shaped one.
Figure 2. The continuous casting mold and the computational mesh employed.
Figure 3. Phase distribution and velocity vectors for the conventional cylindrical SEN (left) and the anchor-shaped SEN (right) in the vertical plane Y = 0. Red phase is molten steel, blue phase is air. Elapsed time = 90 s.
of the velocity vectors around the anchor-shaped SEN is significantly smaller than the case of the cylindrical SEN. This phenomenon is due to the presence of the solid arms of the anchor-shaped SEN which act as physical barriers to the flow of molten steel.
Figure 5 depicts the velocity vectors in the horizontal plane for Z = 0.8. The plane Z = 0.8 is important given that is very close to the free surface of the mold, where the Karman's vortexes are formed. Two recirculation zones are observed on the left side of the cylindrical SEN. Besides, a downward flow of molten steel is appreciated on the right side of this SEN. The anchor-shaped SEN shows a homogeneous flow on both sides of the plane. This indicates that the probability of Karman’s vortexes formation and powder entrapment is lower than for the cylindrical SEN.
Finally, Figure 6 shows the evolution of the average velocity magnitude in the mold as time proceeds from the beginning of the cast where the molten steel is motionless. This figure was obtained through averaging the molten steel velocity in the mold for different elapsed time intervals. Until 60 s of time has elapsed, the behavior of both SEN is analogous. However, from this time the anchor-shaped SEN presents a lower average
Figure 4. Magnification of Figure 3 in the vicinity of each SEN.
Figure 5. Velocity vectors in the horizontal plane Z = 0.8, close to the free surface. Cylindrical SEN (left), anchor-shaped SEN (right).
Figure 6. Evolution of the average velocity magnitude in the mold for the cylindrical SEN (solid line) and anchor SEN (dotted line).
velocity magnitude than the cylindrical SEN. This means that the anchor-shaped SEN causes less turbulence in the mold.
An anchor-shaped Submerged Entry Nozzle design for the continuous casting of steel was proposed in this work. The performance of this nozzle was numerically compared with the performance of a conventional cylindrical nozzle. From the computer results, the following conclusions arise:
1) The cylindrical SEN exhibits the deepest depression of the free surface and the largest velocity vectors in the upper section of the mold.
2) The velocity vectors near the free surface in the vicinity of the SEN are greater for the cylindrical SEN than that corresponding to the anchor-shaped SEN.
3) The chance for the formation of Karman’s vortexes and powder entrapment becomes small for the anchor-shaped SEN.
 Slezak, W., Korolczuk, M. and Migas, P. (2015) High Temperature Rheometric Measurements of Mold Powder. Archives of Metallurgy and Materials, 60, 289-294. http://dx.doi.org/10.1515/amm-2015-0046
 Thomas, B.G. (2005) Modeling of Continuous-Casting Defects Related to Mold Fluid Flow. The 3rd International Congress on Science & Technology of Steelmaking, AIST, Charlotte, 9-12 May 2005, 847-861.
 Garcia-Hernandez, S., Morales, R.D., Barreto, J.J. and Morales-Higa, K. (2013) Numerical Optimization of Nozzle Ports to Improve the Flu Dynamics by Controlling Backflow in a Continuous Casting Slab Mold. ISIJ International, 53, 1794-1802. http://dx.doi.org/10.2355/isijinternational.53.1794
 Iguchi, M., Yoshida, J., Shimizu, T. and Mizuno, Y. (2000) Model Study on the Entrapment of Mold Powder into Molten Steel. ISIJ International, 40, 685-691. http://dx.doi.org/10.2355/isijinternational.40.685
 Iguchi, M. and Ilegbusi, O.J. (2011) Modeling Multiphase Materials Processes: Gas Liquid Systems. Springer, New York. http://dx.doi.org/10.1007/978-1-4419-7479-2
 Tsutsumi, K., Watanabe, K., Suzuki, M., Nakada, M. and Shiomi, T. (2004) Effect of Properties of Mold Powder Entrapped in Molten Steel in a Continuous Casting Process. 7th International Conference on Molten Slags Fluxes and Salts, The South African Institute of Mining and Metallurgy, Johannesburg, 25-28 January 2004, 803-806.
 Ueda, Y., Kida, T. and Iguchi, M. (2004) Unsteady Pressure Coefficient around an Elliptic Immersion Nozzle. ISIJ International, 44, 1403-1409. http://dx.doi.org/10.2355/isijinternational.44.1403
 Real, C., Miranda, R., Vilchis, C., Barron, M.A., Hoyos, L. and Gonzalez, J. (2006) Transient Internal Flow Characterization of a Bifurcated Submerged Entry Nozzle. ISIJ International, 46, 1183-1191.
 Liang, Z., Wang, N., Zou, Z. and Yu, A. (2009) Optimization of Submerged Entry Nozzle of Continuous Casting Slab. Proceedings of the 2009 TMS Annual Meeting and Exhibition, San Francisco, 16-19 February, 569-574.
 Tsukaguchi, Y., Hayashi, H., Kurimoto, H., Yokoya, S., Marukawa, K. and Tanaka, T. (2010) Development of Swirling-Flow Submerged Entry Nozzle for Slab Casting. ISIJ International, 50, 721-729.
 Kuroda, T., Mizobe, A. and Kurisu, J. (2011) Flow Optimization in the Mould by Port Design Improvement of Submerged Entry Nozzle. Proceedings UNITECR 2011 Congress: 12th Biennial Worldwide Conference on Refractories—Refractories-Technology to Sustain the Global Environment, The American Ceramic Society, 30 October-2 November, Kyoto, 434-437.
 Thomas, B., Yuan, Q., Sivaramakhrisnan, S., Shi, T., Vanka, S.P. and Assar, M.B. (2001) Comparison of Four Methods to Evaluate Fluid Velocities in a Continuous Slab Casting Mold. ISIJ International, 41, 1262-1271. http://dx.doi.org/10.2355/isijinternational.41.1262
 Solorio-Díaz, G., Morales, R.D., Palafox-Ramos, J., García-Demedices, L. and Ramos-Ban- deras, A. (2004) Analysis of Fluid Flow Turbulence in Tundishes Fed by a Swirling Ladle Shroud. ISIJ International, 44, 1024-1032. http://dx.doi.org/10.2355/isijinternational.44.1024
 Hirt, C.W. and Nichols, B.D. (1981) Volume of Fluid (VOF) Method for the Dynamics of Free Boundaries. Journal of Computational Physics, 39, 201-225. http://dx.doi.org/10.1016/0021-9991(81)90145-5