The Effects of a Submerged Entry Nozzle on Flow and Initial Solidification in a Continuous Casting Bloom Mold with Electromagnetic Stirring

The melt flow, level fluctuation, temperature field, and solidification behavior coupled with electromagnetic stirring (EMS) effects in the continuous casting mold region of U71Mn steel bloom were numerically analyzed by commercial computational fluid dynamics (CFD) software named ANSYS FLUENT. The effects of submerged entry nozzle (SEN) structures and the installation methods for optimized four-port SEN on the flow pattern, level fluctuation, heat transfer and initial solidification behavior in a bloom mold loaded with EMS were investigated. The aim is to propose a better SEN condition for the big bloom casting of high railway steel. The water simulation experiments were conducted to show the flow characteristics under different SEN conditions and verify the numerical model of flow pattern. The experimental and numerical simulation results showed that the optimized four-port SEN with diagonal installation cannot only improve the flow pattern of the molten steel by alleviating the level fluctuation and reducing the impact pressure to the wall. It is also beneficial for temperature variation at both bloom surface and corner, as well as the local solidified shell thinning phenomena due to the elimination of impingement effect.


Introduction
The big bloom continuous caster is widely used to produce high-end steel products, such as heavy rail steel, bearing steel and free-cutting steel.Most defects affecting steel quality in the continuous casting (CC) process are associated with metallurgical behavior in the mold, which is largely determined by the submerged entry nozzle (SEN).The SEN can affect the flow pattern, solidified shell growth, free surface status and inclusion removal in the CC mold [1,2].Previous works have developed both experimental and numerical models for analyzing the flow [3][4][5] and initial solidification [6][7][8] in the CC mold under different SEN conditions.These studies found that an appropriate SEN structure cannot only enhance the flow field, but also improve the heat transfer process in the mold of both slabs and blooms significantly.To ensure good mold metallurgical behavior during a big bloom casting process, a proper type of SEN should be proposed and investigated.
At present, the SEN with straight single outlet is unfavorable for flotation removal of non-metallic inclusions and may lead to inactive meniscus.The two-port SEN may cause the solidified shell thinning at the narrow face, but can obtain a high removal rate of inclusions.Both types of SENs are widely applied for bloom castings.Recently, quad-furcated SENs with outlets perpendicular to the mold wall has become a concern with regard to which can alleviate level fluctuation and non-uniform solidification, and improve the removal rate of inclusions.Wu et al. [9] analyzed the metallurgical differences among the three kinds of SEN mentioned above, and found that four-port SEN is beneficial to both initial solidification and inclusion floating.Dang et al. [10] proposed a five-port SEN and compared with the four-port SEN.The new SEN can eliminate the non-uniform growth of solidifying shell and ensure flow stability.Sun et al. [11,12] designed a quad-furcated SEN with outlets at tangential directions to enhance the flow pattern by generating swirling flow and eliminating the non-uniform growth of solidified shell.This type of SEN has a complex structure, and may easily trigger nozzle clogging.The proper SEN structure for bloom casting is becoming multi-port, however it should not be too complicated.Therefore, the four-port SEN needs deeper investigation to obtain maximum metallurgical benefit.Meanwhile, as a widely applied technology for big bloom casting, in-mold electromagnetic stirring (M-EMS) affects the flow and solidification behavior in the bloom mold significantly.Ren et al. [13,14] investigated the flow and solidification process in the bloom mold with EMS under normal one-port SEN and found EMS can eliminate the superheat and weaken the immersion depth.To investigate the flow, heat transfer, and initial solidification process in the bloom mold under different SEN conditions more precisely, the effects of M-EMS should not be ignored.
The two-port SENs were formerly applied for the bloom casting of high rail steel U71Mn produced by the 3# con-caster in No. 1 Steelmaking Plant of Wuhan Iron & Steel Corporation (WISCO), where the slag entrapment and disqualification of non-metallic inclusions frequently occurred during the casting process.To solve this problem, an optimized four-port SEN with certain divergent and downward angle was designed to reduce the impingement effect of jet flow to the mold wall.Furthermore, researchers proposed a new idea to change the installation method of the optimized SEN from conventional orthogonal with outlets in the directions of mold wall centers to diagonal with outlets pointing at the mold corners.To prove the exact metallurgical advantage of the diagonally installed four-port SEN over others, the metallurgical differences in the bloom mold under the SENs of two-port, four-port with orthogonal and four-port with diagonal were discussed and compared in this paper.Water model experiments without EMS were conducted to discuss the basic difference of the flow pattern under the different SEN conditions to verify the applicability of the numerical turbulence model.Furthermore, the flow pattern, level fluctuation, heat transfer and initial solidification under different SENs were investigated and compared by a mathematical model considering the effect of M-EMS.

Basic Assumptions
In order to simulate the thermal-physical phenomena in continuous casting bloom mold, the following assumptions are made in the formulating of models: (1) The turbulent flow and solidification phenomena of the steel in the mold are assumed to be at steady state.(2) The influence of mold oscillation and mold taper is ignored.
(3) The molten steel is assumed to be homogeneous incompressible Newtonian fluid.(4) Only the liquid slag layer is considered, other state slag layers are ignored.(5) The curvatures of both mold and strand are not taken into account.

Model of Fluid Flow
The flow pattern in the mold is numerically simulated by solving the continuity and momentum equations under turbulent conditions.The equations are as follows: where P is the pressure and u i is the flow velocity in direction i; ρ is the volume fraction average density; µ e f f is the effective viscosity coefficient, which is the sum of the laminar (µ l ) and turbulent (µ t ) viscosity.µ t can be defined as The low-Reynolds number k-ε turbulent model [15] is applied to describe the turbulence flow and solidification phenomena in the computational domain.The governing equations can be written as follows: where S k = µ l k K P and S ε = µ l ε K P are the solidification induced source terms; the empirical constants and functions in the low-Reynolds number model are presented in Table 1.

Heat Transfer Model
To obtain a precise prediction of the temperature field and solidification behavior of steel, the enthalpy equation is applied [16]: where λ and ρ 1 are the thermal conductivity and density of the steel, respectively, and u is the fluid velocity.The total enthalpy H can be described as where the c p is specific heat, L is latent heat of molten steel solidification, and the liquid fraction f l is updated by following formulations:

Electromagnetism Model
The source term F E in Equation (2) represents the electromagnetic force, which contains radial force (F r ) and tangential force (F θ ) [17]: where B 0 is the magnetic induction density at the boundary of molten steel; σ is the electrical conductivity of liquid steel; f is the current frequency of the EMS; V θ is the tangential velocity; µ 0 is the magneto-conductivity; and r is the value of radial displacement.The radial force (F r ) is ignored for F r , and is much smaller than F θ in a rotating magnetic field.

Effect of Phase Transfer
The source term of moment in Equation (2), S P , which is called Darcy's source term, is utilized to account for the effect of phase change on the convection and turbulence in the mushy zone, which can be considered as porous medium.It can be described by the Darcy's law, the permeability K p is presented as follows [18]: Permeability parallel to primary dendrite arms: Permeability perpendicular to primary dendrite arms: where λ 1 , λ 2 are the primary arm space and secondary arm space, respectively; and u s, i is the velocity of solid in the direction i.
In the calculations, the Darcy's source terms in three directions are calculated by Equation ( 13) and added into FLUENT calculations by User Defined Functions (UDF), in which K par p is applied in x and y direction, and K per p is used in z direction.

VOF Model
In this paper, the steel/slag interface fluctuation in the mold is tracked by volume of fluid (VOF) method which can model two or more immiscible fluids by solving a single set of momentum equations and tracking the volume fraction of each of the fluids throughout the domain.In the VOF model, for two kinds of fluids, a set of momentum equations are adopted in common; for incompressible fluid, it is assumed that the densities of liquid steel and slag are constant, so the volume fraction of liquid steel should conform to the following equation [16,19]: The primary-phase volume fraction will be computed based on the following constraint: The volume fraction average density should be: Metals 2017, 7, 146 5 of 16 where α q means the phase fraction of q and ρ q is the density of phase q.In this paper, the behavior of steel and slag interface is discussed, so the value of n is set to be 2.All other properties (for example, viscosity) are computed in this manner.The surface tension source term in Equation (2), F T , which is calculated by following equation: where σ t is the surface tension coefficient, and k i is the curvature.
In the FLUENT algorithm, the geometric reconstruction scheme is used to describe the morphology of the phase interface [20].

Simulation Models and Parameters
A quarter of three-dimensional (3D) schematics of the top part of the computational domain for the two-port SEN is given in Figure 1, in which the immersion depth of SEN is 210 mm, the thickness of slag film is 30 mm, and the distribution of the two phases can be obtained.where means the phase fraction of and is the density of phase .In this paper, the behavior of steel and slag interface is discussed, so the value of n is set to be 2.All other properties (for example, viscosity) are computed in this manner.
The surface tension source term in Equation ( 2), , which is calculated by following equation: where is the surface tension coefficient, and is the curvature.In the FLUENT algorithm, the geometric reconstruction scheme is used to describe the morphology of the phase interface [20].

Simulation Models and Parameters
A quarter of three-dimensional (3D) schematics of the top part of the computational domain for the two-port SEN is given in Figure 1, in which the immersion depth of SEN is 210 mm, the thickness of slag film is 30 mm, and the distribution of the two phases can be obtained.Full geometrical models are taken for all cases.For example, the meshed computational model equipped with the two-port SEN is shown in Figure 2, where the technology of local grid refinement is applied to simulate the behavior of the initial solidified shell more accurately.The meshes of FLUENT computational domain include non-uniform grids with approximately 1,700,000 cells.To better understand the effect of SEN structure and SEN installation on the metallurgical behavior in the CC bloom mold, water model experiments and numerical simulation cases under different SEN conditions are conducted in this paper.The SEN structures and install methods discussed in this paper are presented in Figure 3.In all SENs, the immersion depth, and nozzle outlet angle are fixed at 210 mm and 15°, respectively.Parallel calculations are performed on an Intel cluster using one node.Each calculation for different SEN conditions lasts about one day.Full geometrical models are taken for all cases.For example, the meshed computational model equipped with the two-port SEN is shown in Figure 2, where the technology of local grid refinement is applied to simulate the behavior of the initial solidified shell more accurately.The meshes of FLUENT computational domain include non-uniform grids with approximately 1,700,000 cells.To better understand the effect of SEN structure and SEN installation on the metallurgical behavior in the CC bloom mold, water model experiments and numerical simulation cases under different SEN conditions are conducted in this paper.The SEN structures and install methods discussed in this paper are presented in Figure 3.In all SENs, the immersion depth, and nozzle outlet angle are fixed at 210 mm and 15 • , respectively.Parallel calculations are performed on an Intel cluster using one node.Each calculation for different SEN conditions lasts about one day.The geometrical parameters and operation conditions of caster are summarized in Table 2. Table 3 lists the thermophysical properties of U71Mn steel and mold flux.As the slag in the upper part of the bloom mold does not participate in the solidification process, and the adiabatic condition is employed at the free surface, some of the slag thermophysical parameters are not necessary to this calculation.Viscosity and thermal conductivity are very important thermophysical parameters to predicting flow and solidification behavior of the steel, especially in the mushy zone.The variation of viscosity and thermal conductivity with local temperature are calculated by software code JMatPro7.0(Sente Software Ltd., Surrey, UK) and applied in the calculation by piecewise linear form using the value points presented in Figure 4.The geometrical parameters and operation conditions of caster are summarized in Table 2. Table 3 lists the thermophysical properties of U71Mn steel and mold flux.As the slag in the upper part of the bloom mold does not participate in the solidification process, and the adiabatic condition is employed at the free surface, some of the slag thermophysical parameters are not necessary to this calculation.Viscosity and thermal conductivity are very important thermophysical parameters to predicting flow and solidification behavior of the steel, especially in the mushy zone.The variation of viscosity and thermal conductivity with local temperature are calculated by software code JMatPro7.0(Sente Software Ltd., Surrey, UK) and applied in the calculation by piecewise linear form using the value points presented in Figure 4.The geometrical parameters and operation conditions of caster are summarized in Table 2. Table 3 lists the thermophysical properties of U71Mn steel and mold flux.As the slag in the upper part of the bloom mold does not participate in the solidification process, and the adiabatic condition is employed at the free surface, some of the slag thermophysical parameters are not necessary to this calculation.Viscosity and thermal conductivity are very important thermophysical parameters to predicting flow and solidification behavior of the steel, especially in the mushy zone.The variation of viscosity and thermal conductivity with local temperature are calculated by software code JMatPro7.0(Sente Software Ltd., Surrey, UK) and applied in the calculation by piecewise linear form using the value points presented in Figure 4.

Boundary and Initial Conditions
The initial and boundary conditions for velocity, turbulent kinetic energy and temperature are set according to set procedure.For inlet, = = × = 0.5679 m/s, = = 0, = 0.01 , = ./ , = 1765.15K.It is noted that the kin and are calculated by the semi-empirical equations [21].For the computational outlet, the fully developed flow condition is applied, where the normal gradients of all variables are set to zero.For the free surface, the normal derivative of all variables is set to zero, and the adiabatic condition is employed at the free surface.
For the bloom wall, the detailed boundary conditions for heat transfer at the bloom wall are presented.In the mold region, the heat exchange effect in the domain is described by the average heat flux [22]: where is the average heat flux density, W/m 2 ; is the density of cooling water, kg/m 3 ; is specific heat capacity of cooling water, J/(kg • K); is the flow rate of cooling water, m 3 /s; ∆ is the water temperature difference in and out of cooling wall, K; and is the effective contact area of liquid steel and mold wall, m 2 .As the initial solidified shell is formed, the air gap between the mold copper and shell would lead to the uneven heat transfer of the steel bloom, reducing the heat flux of the corner [23].In this work, a correction factor ( = 0.5) of heat flux density on the corner ( ) is necessarily adopted while calculating the heat transfer behavior at the bloom corner: when the bloom is pulled out of the mold, the heat transfer coefficient of zone I is calculated by following equation, in which water cooling is applied [24].

Boundary and Initial Conditions
The initial and boundary conditions for velocity, turbulent kinetic energy and temperature are set according to set procedure.For inlet, in /D, T in = 1765.15K.It is noted that the k in and ε are calculated by the semi-empirical equations [21].For the computational outlet, the fully developed flow condition is applied, where the normal gradients of all variables are set to zero.For the free surface, the normal derivative of all variables is set to zero, and the adiabatic condition is employed at the free surface.
For the bloom wall, the detailed boundary conditions for heat transfer at the bloom wall are presented.In the mold region, the heat exchange effect in the domain is described by the average heat flux [22]: where q is the average heat flux density, W/m 2 ; ρ w is the density of cooling water, kg/m 3 ; c w is specific heat capacity of cooling water, J/(kg•K); W is the flow rate of cooling water, m 3 /s; ∆T is the water temperature difference in and out of cooling wall, K; and S is the effective contact area of liquid steel and mold wall, m 2 .As the initial solidified shell is formed, the air gap between the mold copper and shell would lead to the uneven heat transfer of the steel bloom, reducing the heat flux of the corner [23].
In this work, a correction factor (α = 0.5) of heat flux density on the corner (q g ) is necessarily adopted while calculating the heat transfer behavior at the bloom corner: Metals 2017, 7, 146 when the bloom is pulled out of the mold, the heat transfer coefficient of zone I is calculated by following equation, in which water cooling is applied [24].
where h is the heat transfer coefficient at the bloom surface, W/ m 2 •K ; w is the spraying water density, L/ m 2 •s ; T w represents the ambient temperature, K; and β is the correction factor of heat transfer coefficient in the secondary cooling zone.

Water Experiment and Flow Validation
There is no proper experimental method to validate EMS modeling.An indirect way could be the verification of its effect, for example, the amount of solidified shell in function of meniscus distance.In this paper, researchers conducted several water model experiments without EMS to obtain a brief understanding of the flow pattern in the mold under different SEN conditions and verify the numerical model without EMS.Subsequently, the electromagnetic force is calculated by Equation ( 10) and added into FLUENT by user defined functions as a momentum source to further investigate the mold metallurgical behavior coupled with EMS.The distribution of electromagnetic intensity in the mold is measured by Gauss meter, and the result is presented in Figure 5.
where h is the heat transfer coefficient at the bloom surface, W/(m • K); is the spraying water density, L/(m • s); represents the ambient temperature, K; and is the correction factor of heat transfer coefficient in the secondary cooling zone.

Water Experiment and Flow Validation
There is no proper experimental method to validate EMS modeling.An indirect way could be the verification of its effect, for example, the amount of solidified shell in function of meniscus distance.In this paper, researchers conducted several water model experiments without EMS to obtain a brief understanding of the flow pattern in the mold under different SEN conditions and verify the numerical model without EMS.Subsequently, the electromagnetic force is calculated by Equation ( 10) and added into FLUENT by user defined functions as a momentum source to further investigate the mold metallurgical behavior coupled with EMS.The distribution of electromagnetic intensity in the mold is measured by Gauss meter, and the result is presented in Figure 5.

Water Model Experiment
A Water model experiment is useful to investigate fluid flow phenomena and helpful for the guidance of quality improvement and process parameter optimization throughout the continuous casting process.The detailed information about the experiment theory, procedure and device can be found elsewhere [25].In this paper, the water model experiment (ratio 1:1) under three SEN conditions without EMS are conducted, and the experiment results of turbulent flow under different SEN conditions are present in Figure 6.In the water experiment, the molds are filled with transparent water in the initial stage, then the liquid with blue tracer is poured into the mold through different SEN conditions.The flow pattern at the early stage of pouring process now can be easily tracked.The level heights of four key positions on the free surface during the water experiment are detected to characterize the level fluctuation under different SEN conditions.The key positions are at center of narrow face (1#), nearby the nozzle (2#), mold corner (3#), and center of wide face (4#), respectively.The results of average level fluctuations at these positions under different SEN conditions are listed in Table 4. Compared with A-type SEN, the impact effect of jet flow on the narrow face is reduced through the B-type SEN with orthogonal installation, but a significant impact on wide face is formed.This means B-type SEN with orthogonal installation may alleviate the uneven growth of solidifying shell around the narrow side, but cause the non-uniform shell at the wide side near the impingement zone.The level fluctuations under B-type SEN with orthogonal installation apparently decrease at 1#, 2# and 3# positions, but increase to 0.5 mm at 4# position.When changing the four-port SEN to diagonal installation, the level fluctuations at 1#, 2# and 4# positions are small, but would lead to a higher fluctuation at the mold corner (0.7 mm).The

Water Model Experiment
A Water model experiment is useful to investigate fluid flow phenomena and helpful for the guidance of quality improvement and process parameter optimization throughout the continuous casting process.The detailed information about the experiment theory, procedure and device can be found elsewhere [25].In this paper, the water model experiment (ratio 1:1) under three SEN conditions without EMS are conducted, and the experiment results of turbulent flow under different SEN conditions are present in Figure 6.In the water experiment, the molds are filled with transparent water in the initial stage, then the liquid with blue tracer is poured into the mold through different SEN conditions.The flow pattern at the early stage of pouring process now can be easily tracked.The level heights of four key positions on the free surface during the water experiment are detected to characterize the level fluctuation under different SEN conditions.The key positions are at center of narrow face (1#), nearby the nozzle (2#), mold corner (3#), and center of wide face (4#), respectively.The results of average level fluctuations at these positions under different SEN conditions are listed in Table 4. Compared with A-type SEN, the impact effect of jet flow on the narrow face is reduced through the B-type SEN with orthogonal installation, but a significant impact on wide face is formed.This means B-type SEN with orthogonal installation may alleviate the uneven growth of solidifying shell around the narrow side, but cause the non-uniform shell at the wide side near the impingement Metals 2017, 7, 146 9 of 16 zone.The level fluctuations under B-type SEN with orthogonal installation apparently decrease at 1#, 2# and 3# positions, but increase to 0.5 mm at 4# position.When changing the four-port SEN to diagonal installation, the level fluctuations at 1#, 2# and 4# positions are small, but would lead to a higher fluctuation at the mold corner (0.7 mm).The level behavior is more acceptable than the other two conditions.The impingement effect of jet flow to both narrow and wide faces can be limited, and the uniform growth of solidifying shell at both wide and narrow side may be obtained by this SEN condition.However, the detailed behavior of solidifying shell under different SEN conditions should be further investigated by numerical method.
Metals 2017, 7, 146 9 of 16 level behavior is more acceptable than the other two conditions.The impingement effect of jet flow to both narrow and wide faces can be limited, and the uniform growth of solidifying shell at both wide and narrow side may be obtained by this SEN condition.However, the detailed behavior of solidifying shell under different SEN conditions should be further investigated by numerical method.

Flow Validation
The simulation result of flow pattern under the four-port SEN with orthogonal installation is verified by water model experiment under the same SEN condition.Figure 7 shows the comparison of water simulation and numerical modeling results at an early stage of the pouring process.It can be seen that the simulated flow pattern is basically consistent with the flow path and distribution of water with blue tracer, which means that it is feasible to analyze the influence of different nozzle structures and installation methods on the metallurgical behavior in the mold by numerical simulation.

Flow Validation
The simulation result of flow pattern under the four-port SEN with orthogonal installation is verified by water model experiment under the same SEN condition.Figure 7 shows the comparison of water simulation and numerical modeling results at an early stage of the pouring process.It can be seen that the simulated flow pattern is basically consistent with the flow path and distribution of water with blue tracer, which means that it is feasible to analyze the influence of different nozzle structures and installation methods on the metallurgical behavior in the mold by numerical simulation.
Metals 2017, 7, 146 9 of 16 level behavior is more acceptable than the other two conditions.The impingement effect of jet flow to both narrow and wide faces can be limited, and the uniform growth of solidifying shell at both wide and narrow side may be obtained by this SEN condition.However, the detailed behavior of solidifying shell under different SEN conditions should be further investigated by numerical method.

Flow Validation
The simulation result of flow pattern under the four-port SEN with orthogonal installation is verified by water model experiment under the same SEN condition.Figure 7 shows the comparison of water simulation and numerical modeling results at an early stage of the pouring process.It can be seen that the simulated flow pattern is basically consistent with the flow path and distribution of water with blue tracer, which means that it is feasible to analyze the influence of different nozzle structures and installation methods on the metallurgical behavior in the mold by numerical simulation.

Simulation Results and Discussion
The effects of two kinds of SENs and installation methods for four-port SEN on the level fluctuation, turbulent flow, heat transfer, and initial solidification in the CC mold zone are also investigated by numerical simulation.All the simulation cases are performed with mold EMS (600 A, 2 Hz), by which the rotational flow field is formed in the mold zone and the shock pressure on the mold wall can be lightened.This is beneficial to inclusion floatation and improvement of heat transfer, reducing the impacting depth of the poured steel [26].

Turbulent Flow
Figure 8 shows the velocity contour and vector at the symmetric x-z plane (y = 0) for A-type SEN (a), B-type SEN with orthogonal installation (b), and at the diagonal plane for B-type SEN with diagonal installation (c).The velocity contour and vector at four different x-y planes (bloom cross sections) under the three cases are also presented in Figure 8 to visually display the different methods of steel injection, which are z = 0 m at free surface, z = 0.234 m at the nozzle outlets center, z = 0.42 m at the EMS center, and z = 0.7 m at mold outlet, respectively.To conveniently describe the simulation conditions and results, the three cases are named as Case A for two-port SEN, Case B for four-port SEN with orthogonal installation and Case C for four-port SEN with diagonal installation.Figure 9 shows the three-dimensional level fluctuations under the different cases, in which the plane of steel volume fraction value 0.5 is chosen to express the status of level fluctuation.

Simulation Results and Discussion
The effects of two kinds of SENs and installation methods for four-port SEN on the level fluctuation, turbulent flow, heat transfer, and initial solidification in the CC mold zone are also investigated by numerical simulation.All the simulation cases are performed with mold EMS (600 A, 2 Hz), by which the rotational flow field is formed in the mold zone and the shock pressure on the mold wall can be lightened.This is beneficial to inclusion floatation and improvement of heat transfer, reducing the impacting depth of the poured steel [26].

Turbulent Flow
Figure 8 shows the velocity contour and vector at the symmetric x-z plane (y = 0) for A-type SEN (a), B-type SEN with orthogonal installation (b), and at the diagonal plane for B-type SEN with diagonal installation (c).The velocity contour and vector at four different x-y planes (bloom cross sections) under the three cases are also presented in Figure 8 to visually display the different methods of steel injection, which are z = 0 m at free surface, z = 0.234 m at the nozzle outlets center, z = 0.42 m at the EMS center, and z = 0.7 m at mold outlet, respectively.To conveniently describe the simulation conditions and results, the three cases are named as Case A for two-port SEN, Case B for four-port SEN with orthogonal installation and Case C for four-port SEN with diagonal installation.Figure 9 shows the three-dimensional level fluctuations under the different cases, in which the plane of steel volume fraction value 0.5 is chosen to express the status of level fluctuation.It can be intuitively seen from the z = 0.234 m planes under different cases that the impingement effect of jet flow on the mold wall differs from each other.The strong velocity at the center of narrow sides but without impingement on the wide sides for Case A, the remarkable impingement effect on both wide and narrow sides for Case B, and the feeble impingement and velocity around the center of both wide and narrow sides for Case C, respectively.The velocity field difference is gradually reduced along the casting direction and turned to swirling flow by the effect of EMS.The velocity of molten steel feeding into the mold through A-type SEN is about 0.75 m/s and 0.15 m/s larger than B-type SEN.This may cause a higher shock pressure of the narrow face in the mold, then forms a larger and stronger circular flow upon the impacting point, which can easily induce the severe fluctuation of the free surface.This also increases the chance of slag entrapment.The center of circulations upon the impacting point in Case A and Case B is 0.185 m and 0.23 m, and the big range and quantity of level fluctuation in Case A is more obvious than in Case B. For the Case B and Case C having the same structure of SEN, but with different installation methods, the flow velocity of molten steel injected from the nozzle outlets are the same.An installation method such as Case C can increase the distance between the nozzle outlets and mold wall, which means low shock velocities on mold corner walls, and less shock pressure of liquid steel to the mold wall.For Case C, the flow intensity of upper circulation with the center at approximately 0.21 m  It can be intuitively seen from the z = 0.234 m planes under different cases that the impingement effect of jet flow on the mold wall differs from each other.The strong velocity at the center of narrow sides but without impingement on the wide sides for Case A, the remarkable impingement effect on both wide and narrow sides for Case B, and the feeble impingement and velocity around the center of both wide and narrow sides for Case C, respectively.The velocity field difference is gradually reduced along the casting direction and turned to swirling flow by the effect of EMS.The velocity of molten steel feeding into the mold through A-type SEN is about 0.75 m/s and 0.15 m/s larger than B-type SEN.This may cause a higher shock pressure of the narrow face in the mold, then forms a larger and stronger circular flow upon the impacting point, which can easily induce the severe fluctuation of the free surface.This also increases the chance of slag entrapment.The center of circulations upon the impacting point in Case A and Case B is 0.185 m and 0.23 m, and the big range and quantity of level fluctuation in Case A is more obvious than in Case B. For the Case B and Case C having the same structure of SEN, but with different installation methods, the flow velocity of molten steel injected from the nozzle outlets are the same.An installation method such as Case C can increase the distance between the nozzle outlets and mold wall, which means low shock velocities on mold corner walls, and less shock pressure of liquid steel to the mold wall.For Case C, the flow intensity of upper circulation with the center at approximately 0.21 m It can be intuitively seen from the z = 0.234 m planes under different cases that the impingement effect of jet flow on the mold wall differs from each other.The strong velocity at the center of narrow sides but without impingement on the wide sides for Case A, the remarkable impingement effect on both wide and narrow sides for Case B, and the feeble impingement and velocity around the center of both wide and narrow sides for Case C, respectively.The velocity field difference is gradually reduced along the casting direction and turned to swirling flow by the effect of EMS.The velocity of molten steel feeding into the mold through A-type SEN is about 0.75 m/s and 0.15 m/s larger than B-type SEN.This may cause a higher shock pressure of the narrow face in the mold, then forms a larger and stronger circular flow upon the impacting point, which can easily induce the severe fluctuation of the free surface.This also increases the chance of slag entrapment.The center of circulations upon the impacting point in Case A and Case B is 0.185 m and 0.23 m, and the big range and quantity of level fluctuation in Case A is more obvious than in Case B. For the Case B and Case C having the same structure of SEN, but with different installation methods, the flow velocity of molten steel injected from the nozzle outlets are the same.An installation method such as Case C can increase the distance between the nozzle outlets and mold wall, which means low shock velocities on mold corner walls, and less shock pressure of liquid steel to the mold wall.For Case C, the flow intensity of upper circulation with the center at approximately 0.21 m below the meniscus is weak, which may reduce the level fluctuation.The highest level fluctuations for the three cases are 11.5 mm, 9.5 mm and 5.5 mm, respectively.In local regions, the maximum height of level fluctuation in Case A and Case B is much higher than Case C, which may easily lead to damage of the steel/slag interface and slag entrainment.The largest value of level fluctuation under Case C is acceptable for the movement of slag, which happens at four corners of the mold free surface due to the circular flow.Therefore, the range of level fluctuation and chance of slag entrapment can be reduced under Case C.

Temperature Field and Solidification
Figure 10 presents the liquid fraction distribution at various cross sections of x-y planes and x-z plane (y = 0) for the three cases.Figure 11 shows the growth of shell thickness at the x-z planes (a) and y-z planes (b) along the casting direction under different SEN conditions.The results of solidified shell behavior show that the solidifying shell increases non-uniformly at the cross section of the x-z plane but uniformly at the y-z plane along the casting direction under A-type SEN due to a lack of flow impingement on the shell of mold wide face.It is worth noting impingement on the narrow face, which can stunt the growth of the solidifying shell near the impacting zone.The thickness of solidified shell at narrow side decreases from 2.78 mm at around 0.2 m below the meniscus to 1.67 mm at 0.35 m below meniscus, while the shell at wide side grows uniformly.The similar phenomenon of solidifying shell growth can also be observed in Case B, in which the four outlets of the nozzle are perpendicular to the mold walls, the uneven phenomenon of the solidifying shell growth at x-z plane is mostly eliminated, because of the lower outlet velocities of molten steel and lower shock pressure on the mold narrow face compared with Case A. But the solidifying shell at the y-z plane has an uneven growth in the flow impacting zone because the distance from nozzle outlets to the mold wide face is 50 mm shorter than to the narrow face.The thickness of solidified shell at wide side increases from 4.92 mm at 0.2 m below the meniscus to 5.82 mm at 0.32 m below the meniscus, which is thicker than other two cases.The thickness of solidified shell at the width direction of mold outlet (z = 0.7 m) and computational outlet (z = 1.23 m) under Case B are 14.28 mm and 26.13 mm respectively, which is about 3.3 mm and 1.51 mm thicker than Case A, respectively.While the thickness of solidified shell at the thickness direction of mold outlet and computational outlet under Case B are 13.98 mm and 26.05 mm, which is about 1.92 mm and 0.81 mm thinner than Case A. Compared with Case B, Case C increases the distance of jet flow to mold wall by changing the installation method of four-port SEN from orthogonal to diagonal, which has a great effect on decreasing the shock pressure on the wall, lowering the level fluctuation range and reducing the chance of slag entrapment.The thickness of solidified shell at mold outlet (z = 0.7 m) is 15.43 mm at width direction and 14.02 mm at thickness direction, respectively.It increases to 26.65 mm at width direction and 26.16 mm at thickness direction for the computational outlet (z = 1.23 m).The results of Figures 10 and 11 show that the solidifying shell increase gradually and uniformly along the casting direction under Case C at both x-z plane and y-z plane due to a lack of flow impingement on the four surfaces of mold on the corner wall.Less shock pressure to the mold wall leads to small effect on the uniformity of solidifying shell.Therefore, the growth behavior of solidifying shell in the bloom mold under the four-port SEN with diagonal installation is more acceptable than the other two cases, and can reduce the risk of surface cracks and liquid breakout.Figure 12 shows the temperature variation at the centerline of the mold corner (dashed line) and the centerline of the bloom wide face (solid line) along the casting direction for three different cases.The results show that the surface and corner temperature of the bloom upon the injection impacting zone has the same variation under different cases, but the downward trend of surface temperature under Case B has been changed when the bloom reaches the impacting zone.There is a smaller increase, then the surface temperature continues to decrease after the bloom is pulled out of the region.The wide surface temperature at mold outlet (z = 700 mm) under Case A, Case B and Case C are 1310 K, 1344 K and 1326 K, respectively.Then, the temperature in all cases rises to 1345 K at the computational outlet (z = 1230 mm) because of the change of heat transfer condition from mold region to secondary cooling zone.The reason for the mold corner temperature variation around the impacting zone under Case C is similar to the wild wall surface temperature variation under Case B. The impacting location of the Case C is on the mold corner.The corner temperature at the mold outlet for the three cases are 1150 K, 1146 K and 1183 K, which means that Case C is less   Figure 12 shows the temperature variation at the centerline of the mold corner (dashed line) and the centerline of the bloom wide face (solid line) along the casting direction for three different cases.The results show that the surface and corner temperature of the bloom upon the injection impacting zone has the same variation under different cases, but the downward trend of surface temperature under Case B has been changed when the bloom reaches the impacting zone.There is a smaller increase, then the surface temperature continues to decrease after the bloom is pulled out of the region.The wide surface temperature at mold outlet (z = 700 mm) under Case A, Case B and Case C are 1310 K, 1344 K and 1326 K, respectively.Then, the temperature in all cases rises to 1345 K at the computational outlet (z = 1230 mm) because of the change of heat transfer condition from mold region to secondary cooling zone.The reason for the mold corner temperature variation around the impacting zone under Case C is similar to the wild wall surface temperature variation under Case B. The impacting location of the Case C is on the mold corner.The corner temperature at the mold outlet for the three cases are 1150 K, 1146 K and 1183 K, which means that Case C is less Figure 12 the temperature variation at the centerline of the mold corner (dashed line) and the centerline of the bloom wide face (solid line) along the casting direction for three different cases.The results show that the surface and corner temperature of the bloom upon the injection impacting zone has the same variation under different cases, but the downward trend of surface temperature under Case B has been changed when the bloom reaches the impacting zone.There is a smaller increase, then the surface temperature continues to decrease after the bloom is pulled out of the region.The wide surface temperature at mold outlet (z = 700 mm) under Case A, Case B and Case C are 1310 K, 1344 K and 1326 K, respectively.Then, the temperature in all cases rises to 1345 K at the computational outlet (z = 1230 mm) because of the change of heat transfer condition from mold region to secondary cooling zone.The reason for the mold corner temperature variation around the impacting zone under Case C is similar to the wild wall surface temperature variation under Case B. The impacting location of the Case C is on the mold corner.The corner temperature at the mold outlet for the three cases are Metals 2017, 7, 146 14 of 16 1150 K, 1146 K and 1183 K, which means that Case C is less likely to cause the corner crack on the continuous casting bloom.The corner temperature at computational outlet rises to approximately 67 K, 81 K and 21 K to 1217 K, 1227 K and 1202 K at the secondary cooling zone (z = 700 mm-1230 mm).A lower temperature increase means less chance to cause the longitudinal corner crack of the continuous casting bloom.The temperature difference of corner and wide surface at mold outlet under these three cases are 160 K, 198 K and 143 K respectively, a low temperature difference between the bloom surface and corner along the casting direction can be observed under the Case C.This means less chance to form transverse cracks of the steel.To avoid repetition, the temperature variation at narrow face under three cases along the casting direction is not described in this paper, because the variation is similar to the wide surface, where the temperature at impingement zone of narrow surface also has a sudden rise under Case A due to the strong impact to the wall.
likely to cause the corner crack on the continuous casting bloom.The corner temperature at computational outlet rises to approximately 67 K, 81 K and 21 K to 1217 K, 1227 K and 1202 K at the secondary cooling zone (z = 700 mm-1230 mm).A lower temperature increase means less chance to cause the longitudinal corner crack of the continuous casting bloom.The temperature difference of corner and wide surface at mold outlet under these three cases are 160 K, 198 K and 143 K respectively, a low temperature difference between the bloom surface and corner along the casting direction can be observed under the Case C.This means less chance to form transverse cracks of the steel.To avoid repetition, the temperature variation at narrow face under three cases along the casting direction is not described in this paper, because the variation is similar to the wide surface, where the temperature at impingement zone of narrow surface also has a sudden rise under Case A due to the strong impact to the wall.
By comprehensively considering and comparing the results of Case A, Case B and Case C, the results of flow pattern, level fluctuation and solidification behavior in the computational zone under Case C are more likely to fulfill the requirements for the quality of mold metallurgical behavior for bloom casting.

Conclusions
The water model experiments with different SENs were conducted, and the effects of SENs on the flow pattern, temperature field, and solidification behavior in the mold of heavy railway steel were calculated and compared by FLUENT 14.5 (FLUENT Inc., New York, NY, USA).The main conclusions of this paper are summarized as follows: (1) Water model simulation results illustrate that four-port SEN with diagonal installation can alleviate the impingement effect of jet flow both on the narrow and wide face, which may be beneficial to the uniform growth of solidifying shell.(2) When using the four-port SEN with diagonal installation, the maximum level fluctuation can be reduced to 5.5 mm, compared with two-port SEN and four-port SEN with orthogonal installation of which is 11.5 mm and 9.5 mm respectively, and the shock pressure of the injecting flow to the mold wall can be lightened because of low injecting velocity and long impacting distance to the wall.(3) When using the four-port SEN with diagonal installation, the temperature trend at both wide and narrow surfaces is more acceptable compared with two-port SEN or four-port SEN with orthogonal installation.The corner temperature in the mold zone is about 30 K higher because of the liquid steel impacting effect on the mold corner.By comprehensively considering and comparing the results of Case A, Case B and Case C, the results of flow pattern, level fluctuation and solidification behavior in the computational zone under Case C are more likely to fulfill the requirements for the quality of mold metallurgical behavior for bloom casting.

Conclusions
The water model experiments with different SENs were conducted, and the effects of SENs on the flow pattern, temperature field, and solidification behavior in the mold of heavy railway steel were calculated and compared by FLUENT 14.5 (FLUENT Inc., New York, NY, USA).The main conclusions of this paper are summarized as follows: (1) Water model simulation results illustrate that four-port SEN with diagonal installation can alleviate the impingement effect of jet flow both on the narrow and wide face, which may be beneficial to the uniform growth of solidifying shell.(2) When using the four-port SEN with diagonal installation, the maximum level fluctuation can be reduced to 5.5 mm, compared with two-port SEN and four-port SEN with orthogonal installation of which is 11.5 mm and 9.5 mm respectively, and the shock pressure of the injecting flow to the mold wall can be lightened because of low injecting velocity and long impacting distance to the wall.
(3) When using the four-port SEN with diagonal installation, the temperature trend at both wide and narrow surfaces is more acceptable compared with two-port SEN or four-port SEN with orthogonal installation.The corner temperature in the mold zone is about 30 K higher because of the liquid steel impacting effect on the mold corner.(4) Compared with two-port SEN or four-port SEN with orthogonal installation, local solidified shell thinning phenomenon on both wide and narrow sides is eliminated under four-port SEN with diagonal installation due to a lack of impingement of jet flow.

Figure 1 .
Figure 1.Top part of computational domain.

Figure 1 .
Figure 1.Top part of computational domain.

Figure 2 .
Figure 2. The meshed computational model equipped with two-port SEN and mold chamfer.

Figure 3 .
Figure 3. (a,b) The geometric schematic of three types of the SENs and their installation methods.

Figure 2 . 16 Figure 2 .
Figure 2. The meshed computational model equipped with two-port SEN and mold chamfer.

Figure 3 .
Figure 3. (a,b) The geometric schematic of three types of the SENs and their installation methods.

Figure 3 .
Figure 3. (a,b) The geometric schematic of three types of the SENs and their installation methods.

4 Figure 4 .
Figure 4. Variations of viscosity (a) and thermal conductivity (b) with local temperature.

Figure 4 .
Figure 4. Variations of viscosity (a) and thermal conductivity (b) with local temperature.

Figure 5 .
Figure 5.The distribution of measured electromagnetic intensity in the mold.

Figure 5 .
Figure 5.The distribution of measured electromagnetic intensity in the mold.

Figure 7 .
Figure 7.Comparison of water simulation and numerical modeling results.

Figure 6 .
Figure 6.Water model experiment results under different SEN conditions.

Table 4 .
The average level fluctuations at four positions under different SEN conditions.

Figure 7 .
Figure 7.Comparison of water simulation and numerical modeling results.

Figure 7 .
Figure 7.Comparison of water simulation and numerical modeling results.

Figure 8 .
Figure 8. Velocity contour and vector at the symmetric x-z plane (y = 0) for Case A with A-type SEN (a), Case B with B-type SEN with orthogonal installation (b), and at diagonal plane for Case C with B-type SEN with diagonal installation (c).

Figure 9 .
Figure 9. (a-c) The comparison of three-dimensional level fluctuations in the mold for different cases.

Figure 8 . 16 Figure 8 .
Figure 8. Velocity contour and vector at the symmetric x-z plane (y = 0) for Case A with A-type SEN (a), Case B with B-type SEN with orthogonal installation (b), and at diagonal plane for Case C with B-type SEN with diagonal installation (c).

Figure 9 .
Figure 9. (a-c) The comparison of three-dimensional level fluctuations in the mold for different cases.

Figure 9 .
Figure 9. (a-c) The comparison of three-dimensional level fluctuations in the mold for different cases.

Figure 10 .
Figure 10.(a-c) Liquid fraction distribution at various cross sections and x-z plane for different cases.

Figure 11 .
Figure 11.Shell thickness distributions at the x-z plane (a) and y-z plane (b) along the casting direction for the three cases.

Figure 10 .
Figure 10.(a-c) Liquid fraction distribution at various cross sections and x-z plane for different cases.

Metals 2017, 7 , 146 13 of 16 Figure 10 .
Figure 10.(a-c) Liquid fraction distribution at various cross sections and x-z plane for different cases.

Figure 11 .
Figure 11.Shell thickness distributions at the x-z plane (a) and y-z plane (b) along the casting direction for the three cases.

Figure 11 .
Figure 11.Shell thickness distributions at the x-z plane (a) and y-z plane (b) along the casting direction for the three cases.

Figure 12 .
Figure 12.Temperature variation at centerline of mold corner (dashed line) and bloom wide face (solid line) along the casting direction for three different cases.

Figure 12 .
Figure 12.Temperature variation at centerline of mold corner (dashed line) and bloom wide face (solid line) along the casting direction for three different cases.

Table 1 .
Empirical constants and functions used in the Low-Reynolds number model.

Table 2 .
Parameters and operating conditions.

Table 2 .
Parameters and operating conditions.

Table 2 .
Parameters and operating conditions.

Table 3 .
Thermal physical properties of U71Mn steel.

Table 3 .
Thermal physical properties of U71Mn steel.