Effects of EMS Induced Flow on Solidification and Solute Transport in Bloom Mold

The flow, temperature, solidification, and solute concentration field in a continuous casting bloom mold were solved simultaneously by a multiphysics numerical model by considering the effect of in-mold electromagnetic stirring (M-EMS). The mold metallurgical differences between cases with and without EMS are discussed first, and then the solute transport model verified. Moreover, the effects of EMS current intensity on the metallurgical behavior in the bloom mold were also investigated. The simulated solute distributions were basically consistent with the test results. The simulations showed that M-EMS can apparently homogenize the initial solidified shell, liquid steel temperature, and solute element in the EMS effective zone. Meanwhile, the impingement effect of jet flow and molten steel superheat can be reduced, and the degree of negative segregation in the solidified shell at the mold corner alleviated from 0.74 to 0.78. However, the level fluctuation and segregation degree in the shell around the center of the wide and narrow sides were aggravated from 4.5 mm to 6.2 mm and from 0.84 to 0.738, respectively. With the rise of current intensity the bloom surface temperature, level fluctuation, stirring intensity, uniformity of molten steel temperature, and solute distribution also increased, while the growth velocity of the solidifying shell in the EMS effective zone declined and the solute mass fraction at the center of the computational outlet (z = 1.5 m) decreased. M-EMS with a current intensity of 600 A is more suitable for big bloom castings.


Introduction
Most defects affecting steel quality in the continuous casting (CC) process are associated with metallurgical behavior in the mold including flow pattern, heat transfer, initial solidification, and solute transport etc.At present, in-mold electromagnetic stirring (M-EMS) used to optimize molten steel flow is a widely accepted technique in big bloom casting, which has been proved to be beneficial to enhance the columnar-to-equiaxed transition during solidification and reduce the surface and subsurface defects [1][2][3].Geng et al. [4] and Yu et al. [5] indicated that the EMS parameters affects the metallurgical behavior and the steel quality significantly, an optimum EMS parameter for a certain bloom should be proposed for a higher quality.Recently, numerous mathematical and experimental studies on the CC process involving M-EMS were conducted to investigate the flow field distribution of molten steel [6][7][8][9], and the heat transfer and initial solidification behaviors influenced by EMS induced fluid flow were also studied in past years [10][11][12][13].Yang et al. [10] and Ren et al. [13] recently discussed the effects of EMS parameters on the flow pattern and initial solidification in the big bloom mold fed by a normal nozzle with a single outlet, and stated that M-EMS can promote superheat dissipation and enhance the percentage of the equiaxed zone, while the mutual effect between the EMS induced flow and the growth of the solidifying shell lacked discussion.Sun et al. [14] who designed a swirling flow nozzle for round bloom, concluded that a horizontal swirling flow in the mold can eliminate the superheat and depress local shell thinning.However, the solute transport phenomena associated with the EMS induced turbulent flow during the initial solidification process in the bloom mold has been rarely reported, although this deeply affects the macrosegregation profile of the final products.Aboutalebi et al. [15] and Yang et al. [16] developed a fully coupled flow, heat and solute transport model which predicted a turbulent flow induced solidification and segregation profile for CC round bloom and slab, respectively, while the relation of flow and segregation profile was not clearly explained and the effect of M-EMS was ignored.Lei et al. [17] and KwanGu et al. [18] preliminarily investigated the flow, solidification, and solute transport in CC slab mold with an electromagnetic brake, but without detailed explanation of the effect of the electromagnetic field on the solidification process and solute transport, while the verification of the solute transport model was neglected.Sun et al. [19] studied the effect of a swirling flow nozzle (SFN) on the flow, solidification, and macrosegregation behavior of casting blooms, and concluded that the formation of positive and negative segregation in the initial solidified shell can contribute to the flotation of solute-richer molten steel and the "solute washing effect" [20], which could offer the reader useful enlightenment for analyzing the swirling flow induced macrosegregation in the CC mold.However, the effects of M-EMS parameters on the macro-transport phenomena in CC bloom mold, especially the solute transport behavior, need detailed investigation.
In order to gain a deep insight into the M-EMS effects, the flow, temperature field, solidification behavior, and solute concentration field in the CC bloom mold of U71Mn steel were solved simultaneously by a three-dimensional multiphysics numerical model by considering the electromagnetic stirring force.In this paper, the metallurgical effects with EMS on the flow, level fluctuation, heat transfer, initial solidification, and solute transport phenomena in the CC mold were discussed first, and compared to that without EMS.The simulated segregation profiles of solute element carbon (C) in the region of initial solidified shell with a thickness of 30 mm at both the wide and the narrow sides were examined by the infrared carbon-sulfur determinator.Furthermore, the influences of EMS current intensity on the mold metallurgical behavior were also investigated and an optimized current intensity was suggested for big bloom casting.

Basic Assumptions
A three dimensional coupled model was established by combining electromagnetism, fluid flow, heat, and solute transport based on the following assumptions: (1) The transport phenomena in the mold are assumed to be at steady state and the influence of inclusions on the fluid flow, heat transfer, and species transport is neglected to simplify the simulation.(2) The impact of fluid flow on the internal heat transfer of molten steel is ignored in this model, and the liquid steel is assumed to be an incompressible Newtonian fluid.The influence of mold oscillation and mold taper on the fluid flow is also ignored.(3) The effect of thermal contraction on the fluid flow and temperature field in the bloom is neglected.
(4) The mold arc is neglected, and the computational zone is assumed to be a vertical model.(5) The effect of the melt flow on the electromagnetic field is ignored due to the small magnetic Reynolds number (about 0.01) in the stirring process.

Turbulent Flow
The melt flow pattern in the mold is numerically simulated by solving the continuity and momentum equations under turbulent conditions: where µ eff is the effective viscosity coefficient, which is the sum of the laminar (µ l ) and turbulent (µ t ) viscosity.µ t can be defined as: where K and ε are the turbulent kinetic energy and dissipation rate of this energy, respectively.The low-Reynolds number K-ε turbulent model [21,22] is applied to describe the turbulence flow and solidification phenomena in the computational domain.The empirical constants are set as C µ = 0.09 and f µ = exp [ −3.4 The last term in Equation( 2) accounts for the phase interaction force within the mushy zone, and it can be described by Darcy's law, the permeability K p is presented as follows [23]: 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; u s, i is the velocity of solid in the direction of i; f l is the local liquid fraction.The source term F B in the momentum equation represents the full buoyancy force which includes the solutal buoyancy and the thermal buoyancy effect.The detailed information is described by the following equation.
where T l is liquidus temperature, C l,i is locally averaged concentration of solute element i in liquid phase, and C i l, 0 is solute concentration at the liquidus temperature.The density ρ is assumed to be a constant except in the buoyancy terms according to the Boussinesq approximation [24].

Heat Transfer Model
To obtain a precise prediction of the temperature field and solidification behavior, the enthalpy equation is employed.
where u i represents fluid flow velocity in direction of i, the turbulent Prandtl number Pr t is set to the be 0.9 [25], the total enthalpy H can be obtained by the following equation: where c p is specific heat, L is latent heat, and liquid fraction f l is updated by: Metals 2017, 7, 72 where T l and T s are the liquidus and solidus temperature, respectively; T pure is the fusion temperature of pure iron, m i , k i are the slope of liquidus line and equilibrium partition coefficient of solute element i, respectively.C * s, i and C * l, i are the interface concentration of solute element i in the solid and liquid phase, respectively.

Solute Transport Model
The conservation equation for multicomponents in the CC system is expressed as: where S c l = µ l ρD l and S c t = 1 [15] represent the laminar and turbulent Schmidt numbers, respectively.D l is the liquid diffusive coefficient.S c is the source term: where D s is the solid diffusive coefficient, c s and c l are the species mass fraction in the solid and liquid steel, respectively.

Electromagnetism Model
For the source term in Equation ( 2), F E represents the electromagnetic force, which contains radial force (F r ) and tangential force (F θ ) [26]: 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 magnetoconductivity; r is the value of radial displacement.The radial force (F r ) is ignored for F r is much smaller than F θ in a rotating magnetic field.The electromagnetic force is added into FLUENT by user defined functions as a momentum source.The distribution of electromagnetic intensity in the mold at the axis under different current intensity is measured by the Gauss meter, the result is presented in Figure 1.

VOF Model
The steel/slag interface fluctuation in the mold is tracked by the volume fraction of fluid (VOF) method, which can track the shape of the fluctuating steel/slag interface.The detailed information and equations on the VOF model in representing the movement of the mold slag layer can be found

VOF Model
The steel/slag interface fluctuation in the mold is tracked by the volume fraction of fluid (VOF) method, which can track the shape of the fluctuating steel/slag interface.The detailed information and equations on the VOF model in representing the movement of the mold slag layer can be found elsewhere [27].

Simulation Procedure and Verification
The fluid flow, heat, initial solidification behavior, and solute transport in the CC bloom mold were calculated using FLUENT 15.0.The segregation profiles of C in the region of the initial solidification shell at both the wide and the narrow face were measured by an infrared carbon-sulfur determinator and compared with the calculation results.

Operating Condition and Parameters
The target chemical composition of U71Mn steel in the plant is listed in Table 1.The construction and installation method of submerged entry nozzle (SEN) for bloom casting is presented in Figure 2, where the immersion depth is 210 mm and the port downward angle is 15 • , respectively.SEN has been proved to be a practical design in optimizing the fluid and heat transfer in the mold while producing U71Mn steel [28].Table 2 lists the industrial conditions for producing the steel, and the thermo-solutal and physical parameters of U71Mn steel are presented in Table 3.

VOF Model
The steel/slag interface fluctuation in the mold is tracked by the volume fraction of fluid (VOF) method, which can track the shape of the fluctuating steel/slag interface.The detailed information and equations on the VOF model in representing the movement of the mold slag layer can be found elsewhere [27].

Simulation Procedure and Verification
The fluid flow, heat, initial solidification behavior, and solute transport in the CC bloom mold were calculated using FLUENT 15.0.The segregation profiles of C in the region of the initial solidification shell at both the wide and the narrow face were measured by an infrared carbon-sulfur determinator and compared with the calculation results.

Operating Condition and Parameters
The target chemical composition of U71Mn steel in the plant is listed in Table 1.The construction and installation method of submerged entry nozzle (SEN) for bloom casting is presented in Figure 2, where the immersion depth is 210 mm and the port downward angle is 15°, respectively.SEN has been proved to be a practical design in optimizing the fluid and heat transfer in the mold while producing U71Mn steel [28].Table 2 lists the industrial conditions for producing the steel, and the thermo-solutal and physical parameters of U71Mn steel are presented in Table 3.    Viscosity and thermal conductivity are very important thermo-physical parameters for predicting the flow and solidification behavior of the steel, especially the mushy zone.The variations of viscosity and thermal conductivity with local temperature for U71Mn steel shown in Figure 3, which are calculated by the code JMatPro7.0,were applied in the calculation process.

Nozzle adaption
Figure Viscosity and thermal conductivity are very important thermo-physical parameters for predicting the flow and solidification behavior of the steel, especially the mushy zone.The variations of viscosity and thermal conductivity with local temperature for U71Mn steel shown in Figure 3, which are calculated by the code JMatPro7.0,were applied in the calculation process.

Model Building
A quarter of the three-dimensional (3D) schematic sketch of the top part of a real mold is given in Figure 4a.To simulate the behavior of the initial solidified shell more accurately, the technology of local grid refinement is applied, the meshes of FLUENT computational domain include non-uniform grids with cells about 3,100,000.The meshed geometrical model for the computational

Model Building
A quarter of the three-dimensional (3D) schematic sketch of the top part of a real mold is given in Figure 4a.To simulate the behavior of the initial solidified shell more accurately, the technology of local grid refinement is applied, the meshes of FLUENT computational domain include non-uniform grids with cells about 3,100,000.The meshed geometrical model for the computational domain is presented in Figure 4b.The meshed computational model has the same geometrical parameters as the real caster.

Initial and Boundary Conditions
The initial and boundary condition for velocity, temperature, turbulent kinetic energy, dissipation rate, and species concentration are set as follows: For

Initial and Boundary Conditions
The initial and boundary condition for velocity, temperature, turbulent kinetic energy, dissipation rate, and species concentration are set as follows: For computational inlet, in /D, T in = 1765.15K, and C in = C 0 = 0.7246%.k in and ε are calculated by the semiempirical equations [29].For outlet, the fully developed flow condition is adopted, where the normal gradients of all variables are set to zero.For the free surface, the normal derivative of all variables is set to be zero, and the adiabatic condition is adopted for the free surface.
For bloom wall, the detailed boundary conditions for heat transfer at the bloom wall are presented as: In the mold region, the average heat flux is applied to describe the heat exchange effect in the actual domain [30]: where q is the average heat flux density, W/m 2 ; ρ w is the water density, kg/m 3 ; c w is the specific heat capacity of the cooling water, J/(kg•K); W is the water flow rate, m 3 /s; ∆T is the water temperature difference between in and out, K; S is the effective contact area of the steel and the mold wall, m 2 .A correction factor (α) is necessary while calculating the solidification behavior at the bloom mold corner.
The heat transfer coefficient in the secondary cooling zone I is calculated by the following equation, in which water cooling is applied.
where h is the heat transfer coefficient, W/(m 2 •K); w is the water density, L/(m 2 •s); and T w represents the ambient temperature; β is the correction factor.

Verification of Solute Transport
To validate the mathematical coupled model, the mass fraction of the chemical element C in the initial solidified shell thickness of 30 mm is evaluated by an infrared carbon-sulfur determinator.The specific locations of the drilling samples are presented in Figure 5.Where the samples are drilled with a 4.2 mm diameter drill up to a depth about 20 mm and the weight is around 2 g.The test and simulated results for the verification are without EMS.Figure 6 shows the comparison between the simulation results and the result from the infrared carbon-sulfur determinator.A basically consistent variation tendency can be observed between the simulation results and testing results along both width direction and thickness direction.Figure 6 shows the comparison between the simulation results and the result from the infrared carbon-sulfur determinator.A basically consistent variation tendency can be observed between the simulation results and testing results along both width direction and thickness direction.Figure 6 shows the comparison between the simulation results and the result from the infrared carbon-sulfur determinator.A basically consistent variation tendency can be observed between the simulation results and testing results along both width direction and thickness direction.

Metallurgical Effects of M-EMS
Two simulation cases in this paper are discussed first and compared to find out the exact effect of mold EMS on flow, heat transfer, initial solidification, and solute transport in the mold region.One is the case without EMS, the other is the case loaded with EMS, in which the current intensity is 600 A and the electromagnetic frequency is 2.0 Hz.

Flow Field
Figure 7 shows the contour and vector of flow patterns in the upper part of the bloom without EMS (a) and with EMS (b).Also, the velocity contour of several key cross sections in the mold zone for the two cases is presented in Figure 7, which are z = 0.234 m for the center of the SEN outlets, z = 0.42 m for the EMS center and z = 0.7 m for the mold outlet, respectively.In the bloom mold without EMS, the bulk flow leaving the SEN outlets impinges on the mold corner wall and is split into two opposing directions.Upper circulation is confined by the walls and meniscus and the lower circulation is formed by the mold corner and the return flow near the mold center.The center of the upper circulation is located at around 0.09 m bloom from the meniscus, and the lower circulation center is locating at approximately 0.45 m.When the mold is loaded with EMS, four obvious circulations are formed upon and under the EMS center，the two lower circulations can

Metallurgical Effects of M-EMS
Two simulation cases in this paper are discussed first and compared to find out the exact effect of mold EMS on flow, heat transfer, initial solidification, and solute transport in the mold region.One is the case without EMS, the other is the case loaded with EMS, in which the current intensity is 600 A and the electromagnetic frequency is 2.0 Hz.

Flow Field
Figure 7 shows the contour and vector of flow patterns in the upper part of the bloom without EMS (a) and with EMS (b).Also, the velocity contour of several key cross sections in the mold zone for the two cases is presented in Figure 7, which are z = 0.234 m for the center of the SEN outlets, z = 0.42 m for the EMS center and z = 0.7 m for the mold outlet, respectively.In the bloom mold without EMS, the bulk flow leaving the SEN outlets impinges on the mold corner wall and is split into two opposing directions.Upper circulation is confined by the walls and meniscus and the lower circulation is formed by the mold corner and the return flow near the mold center.The center of the upper circulation is located at around 0.09 m bloom from the meniscus, and the lower circulation center is locating at approximately 0.45 m.When the mold is loaded with EMS, four obvious circulations are formed upon and under the EMS center, the two lower circulations can improve the heat transfer mechanism of molten steel in the EMS effect area, eliminate steel superheat, and promote uniform growth of the solidifying shell.The location of the upper circulation center in the diagonal plane changes to 0.21 m, with a higher velocity than the case without EMS, but the circulation under the nozzle outlets is eliminated and the flow is changed to the casting direction because of the increasing stirring effect of the EMS.In the center cross section of the SEN outlets, the four jet flows impinge on the mold corner, rebound to the bloom center and form four horizontal circulations near the narrow sides with a high velocity around the mold corner without EMS.The horizontal circulations remain for a while, then are eliminated at the mold outlet.The horizontal circulations created by the mold wall and molten steel in the bloom center can be changed into a strong horizontal swirling flow gradually by EMS, and the strongest flow occurs in the cross section of the EMS center.In the horizontal swirling flow, the largest tangential velocity locates at the solidification front near the center of the wide and narrow sides.The strong tangential velocity can affect the solidification behavior in the mold, effectively prevent the growth of columnar crystal, promote the growth of equiaxed grain [31], and be beneficial to the homogeneity of temperature and solute element in the molten steel.
changed into a strong horizontal swirling flow gradually by EMS, and the strongest flow occurs in the cross section of the EMS center.In the horizontal swirling flow, the largest tangential velocity locates at the solidification front near the center of the wide and narrow sides.The strong tangential velocity can affect the solidification behavior in the mold, effectively prevent the growth of columnar crystal, promote the growth of equiaxed grain [31], and be beneficial to the homogeneity of temperature and solute element in the molten steel.The comparison of the three-dimensional time-averaged level fluctuation with and without EMS is presented in Figure 8.The height difference between the zenith and the nadir of the free surface level is defined as the largest level fluctuation value of the mold.The largest level fluctuation value without EMS is about 4.5 mm, which happens at very few positions, for most of the other places the level fluctuation is less than 1 mm, the free surface is relatively stable.After the EMS is loaded, the level fluctuation at the center and edge of the free surface in the mold moves up and down respectively, which is affected by the electromagnetic force at the upper mold.The largest value of level fluctuation is 6.2 mm in the case with EMS.For bloom castings, the value of level fluctuation should be controlled as low as possible while the M-EMS can aggravate the movement of mold flux, which will increase the chance of slag entrapment.A favorable electromagnetic stirring effect along with a relatively stable free surface should be obtained to guarantee proper metallurgical behavior.and down respectively, which is affected by the electromagnetic force at the upper mold.The largest value of level fluctuation is 6.2 mm in the case with EMS.For bloom castings, the value of level fluctuation should be controlled as low as possible while the M-EMS can aggravate the movement of mold flux, which will increase the chance of slag entrapment.A favorable electromagnetic stirring effect along with a relatively stable free surface should be obtained to guarantee proper metallurgical behavior.

Heat Transfer and Solidification
Figure 9 presents the temperature contour at the diagonal plane ((a) without EMS, (b) 600 A, 2.0 Hz) and the temperature distribution of molten steel at the computational centerline for the two cases (c).The temperature at the molten steel core without EMS rapidly decreases from 1765 K to 1735.7 K at 0.3 m below the meniscus, and with little change in the rest of the computational zone, at which the centerline temperature decreases from 1735.7 K at 0.3 m below the meniscus to 1728 K at the computational outlet.While the temperature at the liquid core at 0.3 m below the meniscus is 1737 K after the EMS is loaded, the average temperature in the molten steel between 0.3 m and 0.48 m below the meniscus is higher than the case without EMS, and the temperature decreases from 1735.2 K at 0.48 m below the meniscus to 1725.4K at the computational outlet.The temperature of molten steel with EMS is uniformly distributed and remains at around 1725 K from 0.8 m to 1.5 m below the meniscus, while the molten steel temperature without EMS is non-uniformly distributed and remains in remarkable decline at about 2.6 K higher at the center of the computational outlet.Therefore, the mold EMS can effectively reduce and even eliminate the superheat of the molten steel by retaining the molten steel at a high temperature in the upper zone of the mold.
Figure 10 shows the temperature variation at the centerline of the bloom wide face and chamfered corner in the two cases, in which the solid and dashed lines represent the corner and surface temperature, respectively.The variation rule of temperature at both the bloom surface and the corner with and without can be easily understood by combining the flow results mentioned above.The surface and corner temperature at 0-0.2 m below the meniscus is almost the same in the two cases because of the small EMS effect in this region.In the EMS effect area, the descent velocity of the mold corner temperature in the case loaded with EMS is larger than the other one, because the electromagnetic force can apparently lighten the scouring effect of the jet flow on the corner wall and accelerate the heat transfer process of the corner.The corner temperatures at the mold outlet in these two cases are 1286 K and 1180 K, respectively, and the temperature difference decreases gradually as the casting process continues, which is changed to 1241 K and 1195 K at the computational outlet, respectively.The corner temperature changes about ±30 K till the bloom reaches to the straightening area, which effectively avoids the high-temperature brittle zone of the

Heat Transfer and Solidification
Figure 9 presents the temperature contour at the diagonal plane ((a) without EMS, (b) 600 A, 2.0 Hz) and the temperature distribution of molten steel at the computational centerline for the two cases (c).The temperature at the molten steel core without EMS rapidly decreases from 1765 K to 1735.7 K at 0.3 m below the meniscus, and with little change in the rest of the computational zone, at which the centerline temperature decreases from 1735.7 K at 0.3 m below the meniscus to 1728 K at the computational outlet.While the temperature at the liquid core at 0.3 m below the meniscus is 1737 K after the EMS is loaded, the average temperature in the molten steel between 0.3 m and 0.48 m below the meniscus is higher than the case without EMS, and the temperature decreases from 1735.2 K at 0.48 m below the meniscus to 1725.4K at the computational outlet.The temperature of molten steel with EMS is uniformly distributed and remains at around 1725 K from 0.8 m to 1.5 m below the meniscus, while the molten steel temperature without EMS is non-uniformly distributed and remains in remarkable decline at about 2.6 K higher at the center of the computational outlet.Therefore, the mold EMS can effectively reduce and even eliminate the superheat of the molten steel by retaining the molten steel at a high temperature in the upper zone of the mold.Figure 10 shows the temperature variation at the centerline of the bloom wide face and chamfered corner in the two cases, in which the solid and dashed lines represent the corner and surface temperature, respectively.The variation rule of temperature at both the bloom surface and the corner with and without can be easily understood by combining the flow results mentioned above.
The surface and corner temperature at 0-0.2 m below the meniscus is almost the same in the two cases because of the small EMS effect in this region.In the EMS effect area, the descent velocity of the mold corner temperature in the case loaded with EMS is larger than the other one, because the electromagnetic force can apparently lighten the scouring effect of the jet flow on the corner wall and accelerate the heat transfer process of the corner.The corner temperatures at the mold outlet in these two cases are 1286 K and 1180 K, respectively, and the temperature difference decreases gradually as the casting process continues, which is changed to 1241 K and 1195 K at the computational outlet, respectively.The corner temperature changes about ±30 K till the bloom reaches to the straightening area, which effectively avoids the high-temperature brittle zone of the steel and prevents the chance of transverse cracks.The surface temperature in the EMS effective zone is higher than the case without EMS, because the horizontal swirling flow generated by EMS can strongly move the solidification front into the bloom center, then the solidification front is replaced by high temperature molten steel.The bloom surface temperature increases with the change of the cooling condition, which increases from 1272 K at the mold outlet to 1411 K at the computational outlet without EMS and from 1355 K to 1440 K with EMS, respectively.The values of temperature rise are 136 K without EMS and 85 K with EMS, which indicates that EMS in the bloom mold can reduce the possibility of surface cracks.
temperature rise are 136 K without EMS and 85 K with EMS, which indicates that EMS in the bloom mold can reduce the possibility of surface cracks.Figure 11 shows the liquid fraction distributions at the diagonal plane and cross section of the EMS center (z = 0.42 m) for the two different cases.The growth velocity of the solidifying shell at the mold corner in the impact area in the case without EMS is lower than the other because of a remarkable impingement of the jet flow on the corner wall, and the solidified shell thickness at the corner of the computational outlet is apparently thinner.Compared with the case without EMS, the solidified shell in the cross section z = 0.42 m in the case loaded with EMS is thicker at the corner, but thinner at both the wide and narrow sides, and the solidified shell distribution is more uniform.Figure 11 shows the liquid fraction distributions at the diagonal plane and cross section of the EMS center (z = 0.42 m) for the two different cases.The growth velocity of the solidifying shell at the mold corner in the impact area in the case without EMS is lower than the other because of a remarkable impingement of the jet flow on the corner wall, and the solidified shell thickness at the corner of the computational outlet is apparently thinner.Compared with the case without EMS, the solidified shell in the cross section z = 0.42 m in the case loaded with EMS is thicker at the corner, but thinner at both the wide and narrow sides, and the solidified shell distribution is more uniform.
Figure 12 presents the shell thickness and liquid/solid distributions at the y-z plane (x = 0 m) along the casting direction for the two cases, where the solid line refers to the case without EMS, while the dashed line represents the case loaded with EMS.The growth rhythm of the solidifying shell in the width direction is associated with the results shown in Figure 11, the growth velocity of the solidifying shell at the narrow side in the case loaded with EMS is slower in the EMS effective zone and the shell thickness at the mold outlet (z = 0.7 m) is 14.2 mm, which is about 3 mm thinner than the case without EMS.With the decrease of electromagnetic force, the shell thickness difference between the two cases decreases as well, the solidified shell at the computational outlet for the two cases is 26.86 mm with EMS and 28.3 mm without EMS, respectively.The difference in mushy zone length variation between the case with and without EMS can also be observed in Figure 12, in which the length of the mushy zone under the EMS effective zone is significantly smaller, especially in the EMS center, where the mushy zone lengths are 3.25 mm and 5.4 mm respectively.The main reason for the difference in mushy zone length is the remarkable tangential velocity at the solidification front of the narrow side, which is generated by the electromagnetic force.M-EMS can reduce the growth velocity of the solidifying shell and the length of the mushy zone in the EMS effective zone, which makes the shell thickness at the computational outlet 1.44 mm thinner than the case without EMS.
than the case without EMS.With the decrease of electromagnetic force, the shell thickness difference between the two cases decreases as well, the solidified shell at the computational outlet for the two cases is 26.86 mm with EMS and 28.3 mm without EMS, respectively.The difference in mushy zone length variation between the case with and without EMS can also be observed in Figure 12, in which the length of the mushy zone under the EMS effective zone is significantly smaller, especially in the EMS center, where the mushy zone lengths are 3.25 mm and 5.4 mm respectively.The main reason for the difference in mushy zone length is the remarkable tangential velocity at the solidification front of the narrow side, which is generated by the electromagnetic force.M-EMS can reduce the growth velocity of the solidifying shell and the length of the mushy zone in the EMS effective zone, which makes the shell thickness at the computational outlet 1.44 mm thinner than the case without EMS.intuitively the effect of EMS on the distribution of solute element C at the bloom cross section.For the case without EMS, the original liquid steel with initial mass fraction of solute element C (0.725%) is sprayed into the mold through the four outlets of SEN which are pointed directly at the mold corners, the primitive steel channels are formed by the effect of the inertia force, and the liquid steel with a high solute content floats upward to the meniscus dead zone because of the density change of liquid steel with temperature, which leads to the appearance of positive segregation in the 4 mm initial solidified shell.When arriving at the flow impact area, the impingement of liquid steel on the mold corner can dispel the enrichment of the solute at the solidification front and as there is not enough time to supply a rich solute content while the solidification continues, a negative position will be developed under this condition, and the maximum negative segregation in the diagonal plane will occur at the impact point.The high level of negative segregation is gradually reduced after the bloom is pulled out of the EMS effective zone, where the enriched solute steel will again gather at the solidification front.Compared with the result of solute element distribution at the bloom corner, a much weaker impingement effect happens around the wide and narrow sides of the bloom, the degree of negative segregation in the solidified shell of the narrow side is obviously lower.The maximum negative segregation is 0.74 when the corner solidified shell thickness is 10 mm and 0.84 when the shell thickness at the width direction is 15 mm.For the case loaded with EMS, the jet trajectory is gradually changed by the electromagnetic force, the impact pressure is lightened, and the maximum negative segregation at the bloom corner improves to 0.78 where the shell thickness is 12.5 mm compared with the case without EMS.The degree of negative segregation in the solidified shell around the wide and narrow sides in the EMS effective zone is more serious than the case without EMS, because the high tangential velocity near the solidification front around the center of the bloom narrow and wide sides can wash out the solute enriched mushy steel and move it into the bloom center, thereby enhancing the average mass fraction of the solute element C in the molten steel.The largest negative segregation at the width centerline of the computational outlet is 0.738, where the thickness is about 9.2 mm.Compared with the case without EMS, the degree and area of negative segregation in the initial solidified shell in the case loaded with EMS is larger, so the average mass fraction of solute element C in the bloom center is relatively higher due to the mass conservation law.The mass fractions of element C at the center of the computational outlet with and without EMS are 0.7904% and 0.7743%, respectively.
For big bloom castings, the main functions of M-EMS are in eliminating the superheat of the molten steel, making the temperature and solute in the molten steel uniform, improving the non-metallic inclusion floatation, promoting the ability of fluxes for adsorption of inclusion, and improving the solidification structure of the bloom.There is no obvious improvement of M-EMS on centerline segregation of big bloom casting, the effective technologies to reduce the centerline segregation for big bloom casting are soft reduction and final electromagnetic stirring [32].In many steel plants, the combination of M-EMS with soft reduction or final EMS is widely applied to improve the internal quality of the casting bloom.
Metals 2017, 7, 72 14 of 20 fraction of solute element C in the bloom center is relatively higher due to the mass conservation law.The mass fractions of element C at the center of the computational outlet with and without EMS are 0.7904% and 0.7743%, respectively.For big bloom castings, the main functions of M-EMS are in eliminating the superheat of the molten steel, making the temperature and solute in the molten steel uniform, improving the non-metallic inclusion floatation, promoting the ability of fluxes for adsorption of inclusion, and improving the solidification structure of the bloom.There is no obvious improvement of M-EMS on centerline segregation of big bloom casting, the effective technologies to reduce the centerline segregation for big bloom casting are soft reduction and final electromagnetic stirring [32].In many steel plants, the combination of M-EMS with soft reduction or final EMS is widely applied to improve the internal quality of the casting bloom.

Effect of Current Intensity
After analyzing and comparing the results of the flow pattern, heat transfer, initial solidification, and solute transport in the CC bloom mold with and without EMS (600 A, 2.0 Hz), the exact effect of EMS current intensity on the metallurgical behavior in the mold will also now be discussed in this paper.The EMS parameters for the four cases are 450 A, 500 A, 550 A, and 600 A respectively, of which the frequencies are fixed at 2.0 Hz.To avoid repeating, the influence of different EMS frequencies is ignored because of the similar changing rule to current intensity.

Effect of Current Intensity
After analyzing and comparing the results of the flow pattern, heat transfer, initial solidification, and solute transport in the CC bloom mold with and without EMS (600 A, 2.0 Hz), the exact effect of EMS current intensity on the metallurgical behavior in the mold will also now be discussed in this paper.The EMS parameters for the four cases are 450 A, 500 A, 550 A, and 600 A respectively, of which the frequencies are fixed at 2.0 Hz.To avoid repeating, the influence of different EMS frequencies is ignored because of the similar changing rule to current intensity.The velocity vector (a) and distribution of element C (b) at the cross section of EMS center (z = 0.42 m) are presented in Figure 16.With the increase of current intensity, the tangential velocity at the solidification front is obviously increased, especially at the places near the center of the wide and narrow sides; the largest velocity at this moment for these four cases is 0.197 m/s, 0.219 m/s, 0.24 m/s, and 0.271 m/s respectively.The stirring effect is not satisfactory when the current intensity is less than 500 A, at which the distribution of solute element C in the cross section of the EMS center still remains part of the initial injection trajectory of the molten steel, the stirring force is not strong enough to improve the uniformity of solute and temperature in the molten steel, and the ability of inclusion flotation is not significantly enhanced.When the current intensity is larger than 500 A, the stirring effect increases with the rise of current intensity, and the distribution of solute is more and more uniform in the bloom cross section of the EMS effective zone.Figure 17 shows the distribution of solute element C at computational outlets (a) and the centerline along the casting direction (b) under different current intensities.It can be seen that the distribution law of the solute element in the solidified shell under different EMS current intensities is very similar, the maximum mass fraction of solute element C is presented at the solidification front, where the positive segregation will be formed in the mushy zone about to solidify.Compared with the other three cases, the area of solute enriched steel at the solidification front is smaller and the mass fraction of the solute element C at the bloom center is lower in the case with 600 A EMS.The mass fractions of C at the center of the computational outlet are 0.79275% for 450 A, 0.79389% for 500 A, 0.79214% for 550 A, and 0.79042% for 600 A, respectively.The velocity vector (a) and distribution of element C (b) at the cross section of EMS center (z = 0.42 m) are presented in Figure 16.With the increase of current intensity, the tangential velocity at the solidification front is obviously increased, especially at the places near the center of the wide and narrow sides; the largest velocity at this moment for these four cases is 0.197 m/s, 0.219 m/s, 0.24 m/s, and 0.271 m/s respectively.The stirring effect is not satisfactory when the current intensity is less than 500 A, at which the distribution of solute element C in the cross section of the EMS center still remains part of the initial injection trajectory of the molten steel, the stirring force is not strong enough to improve the uniformity of solute and temperature in the molten steel, and the ability of inclusion flotation is not significantly enhanced.When the current intensity is larger than 500 A, the stirring effect increases with the rise of current intensity, and the distribution of solute is more and more uniform in the bloom cross section of the EMS effective zone.The velocity vector (a) and distribution of element C (b) at the cross section of EMS center (z = 0.42 m) are presented in Figure 16.With the increase of current intensity, the tangential velocity at the solidification front is obviously increased, especially at the places near the center of the wide and narrow sides; the largest velocity at this moment for these four cases is 0.197 m/s, 0.219 m/s, 0.24 m/s, and 0.271 m/s respectively.The stirring effect is not satisfactory when the current intensity is less than 500 A, at which the distribution of solute element C in the cross section of the EMS center still remains part of the initial injection trajectory of the molten steel, the stirring force is not strong enough to improve the uniformity of solute and temperature in the molten steel, and the ability of inclusion flotation is not significantly enhanced.When the current intensity is larger than 500 A, the stirring effect increases with the rise of current intensity, and the distribution of solute is more and more uniform in the bloom cross section of the EMS effective zone.Figure 17 shows the distribution of solute element C at computational outlets (a) and the centerline along the casting direction (b) under different current intensities.It can be seen that the distribution law of the solute element in the solidified shell under different EMS current intensities is very similar, the maximum mass fraction of solute element C is presented at the solidification front, where the positive segregation will be formed in the mushy zone about to solidify.Compared with the other three cases, the area of solute enriched steel at the solidification front is smaller and the mass fraction of the solute element C at the bloom center is lower in the case with 600 A EMS.The mass fractions of C at the center of the computational outlet are 0.79275% for 450 A, 0.79389% for 500 A, 0.79214% for 550 A, and 0.79042% for 600 A, respectively.Figure 17 shows the distribution of solute element C at computational outlets (a) and the centerline along the casting direction (b) under different current intensities.It can be seen that the distribution law of the solute element in the solidified shell under different EMS current intensities is very similar, the maximum mass fraction of solute element C is presented at the solidification front, where the positive segregation will be formed in the mushy zone about to solidify.Compared with the other three cases, the area of solute enriched steel at the solidification front is smaller and the mass fraction of the solute element C at the bloom center is lower in the case with 600 A EMS.The mass fractions of C at the center of the computational outlet are 0.79275% for 450 A, 0.79389% for 500 A, 0.79214% for 550 A, and 0.79042% for 600 A, respectively.By comprehensively considering the mold metallurgical behavior at different current intensities, the M-EMS with a current intensity of 600 A was found to be more suitable for big bloom castings than the other cases.

Conclusions
The multi-physical metallurgical behavior in the bloom mold loaded with EMS was investigated by a three-dimensional numerical model.The effects of EMS induced flow on the heat  Figure 18 displays the variation of surface temperature and shell thickness along with casting direction at different current intensities.With the increase of current intensity (from 450 A to 550 A), the surface temperature increases as well, when the current intensity is over 550 A, the surface temperature does not result in any obvious change.The surface temperatures at the mold outlet (z = 700 mm) for each case are 1325.7 K for 450 A, 1343.3K for 500 A, 1359.4K for 550 A, and 1355 K for 600 A, respectively.Also, 1391.5 K, 1400.9K, 1405.8K, and 1403.1 K at the outlet of the secondary cooling zone I (z = 1.23 m), which then increase to 1444.2 K, 1450.7 K, 1447.2K and 1445 K respectively at the computational outlet due to the change of heat transfer conditions, respectively.With the increase of current intensity, the electromagnetic force rises, while the growth velocity of the solidifying shell in the EMS effect zone deceases.The solidified shell in the thickness direction at the mold outlet at four different current intensities is about 15.3 mm for 450 A, 15.2 mm for 500 A, 14.8 mm for 550 A, and 14.2 mm for 600 A, respectively, and increases to 27.55 mm, 27.2 mm, 27.57mm, and 26.86 mm at the outlet of computational zone (z = 1500 mm), respectively.By comprehensively considering the mold metallurgical behavior at different current intensities, the M-EMS with a current intensity of 600 A was found to be more suitable for big bloom castings than the other cases.

Conclusions
The multi-physical metallurgical behavior in the bloom mold loaded with EMS was investigated by a three-dimensional numerical model.The effects of EMS induced flow on the heat By comprehensively considering the mold metallurgical behavior at different current intensities, the M-EMS with a current intensity of 600 A was found to be more suitable for big bloom castings than the other cases.

Metals 2017, 7 , 72 5 of 20 Figure 1 .
Figure 1.Measured distribution of electromagnetic intensity in the mold under different current intensity.

Figure 1 .
Figure 1.Measured distribution of electromagnetic intensity in the mold under different current intensity.

Figure 2 .
Figure 2. The three-dimensional view (a) and installation method (b) of submerged entry nozzle (SEN).Figure 2. The three-dimensional view (a) and installation method (b) of submerged entry nozzle (SEN).

Figure 2 .
Figure 2. The three-dimensional view (a) and installation method (b) of submerged entry nozzle (SEN).Figure 2. The three-dimensional view (a) and installation method (b) of submerged entry nozzle (SEN).

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

Figure 3 .
Figure 3. Variations of viscosity (a) and thermal conductivity (b) with local temperature.
in Figure4b.The meshed computational model has the same geometrical parameters as the real caster.

Figure 4 .
Figure 4.The schematic sketch of mold and meshed model of computational domain.(a) mold top part; (b) meshed model.

Figure 4 .
Figure 4.The schematic sketch of mold and meshed model of computational domain.(a) mold top part; (b) meshed model.

Figure 5 .
Figure 5. Distributions of drilling samples in the initial solidified shell.

Figure 5 .
Figure 5. Distributions of drilling samples in the initial solidified shell.

Figure 5 .
Figure 5. Distributions of drilling samples in the initial solidified shell.

Figure 6 .
Figure 6.Comparison between measured and simulated segregation profiles of C in the region of initial solidification shell.(a) width direction; (b) thickness direction.

Figure 6 .
Figure 6.Comparison between measured and simulated segregation profiles of C in the region of initial solidification shell.(a) width direction; (b) thickness direction.

Figure 7 .
Figure 7. Velocity contour and vector in the upper part of the mold for the cases without (a) and with (b) electromagnetic stirring (EMS).

Figure 7 .
Figure 7. Velocity contour and vector in the upper part of the mold for the cases without (a) and with (b) electromagnetic stirring (EMS).

Figure 8 .
Figure 8.Comparison of three dimensional level fluctuations in the mold without (a) and with (b) EMS.

Figure 8 .
Figure 8.Comparison of three dimensional level fluctuations in the mold without (a) and with (b) EMS.
Metals 2017, 7, 72 11 of 20 steel and prevents the chance of transverse cracks.The surface temperature in the EMS effective zone is higher than the case without EMS, because the horizontal swirling flow generated by EMS can strongly move the solidification front into the bloom center, then the solidification front is replaced by high temperature molten steel.The bloom surface temperature increases with the change of the cooling condition, which increases from 1272 K at the mold outlet to 1411 K at the computational outlet without EMS and from 1355 K to 1440 K with EMS, respectively.The values of temperature rise are 136 K without EMS and 85 K with EMS, which indicates that EMS in the bloom mold can reduce the possibility of surface cracks.

Figure 9 .
Figure 9. Temperature contour at the diagonal planes ((a) without EMS; (b) 600 A, 2.0 Hz) and molten steel temperature at computational centerline (c) for the two cases.

Figure 9 .
Figure 9. Temperature contour at the diagonal planes ((a) without EMS; (b) 600 A, 2.0 Hz) and molten steel temperature at computational centerline (c) for the two cases.

Figure 9 .
Figure 9. Temperature contour at the diagonal planes ((a) without EMS; (b) 600 A, 2.0 Hz) and molten steel temperature at computational centerline (c) for the two cases.

Figure 10 .
Figure 10.Temperature variations at the centerline of the bloom wide face and chamfered corner.

Figure 10 .
Figure 10.Temperature variations at the centerline of the bloom wide face and chamfered corner.

Figure 11 .
Figure 11.Liquid fraction distributions at the diagonal plane and the EMS center for different cases.Figure 11.Liquid fraction distributions at the diagonal plane and the EMS center for different cases.

Figure 11 .
Figure 11.Liquid fraction distributions at the diagonal plane and the EMS center for different cases.Figure 11.Liquid fraction distributions at the diagonal plane and the EMS center for different cases.Metals 2017, 7, 72 13 of 20

Figure 12 .
Figure 12.Liquid and solid distributions at the y-z plane along the casting direction for different cases.

Figure 12 .
Figure 12.Liquid and solid distributions at the y-z plane along the casting direction for different cases.

4. 1
Figure 13 presents the distributions of solute element C on the diagonal plane ((a) the case without, (b) the case loaded with EMS) and at the diagonal line of the computational outlet(c) in the two different cases.Figure 14 shows the distributions of solute element C on the y-z plane, various cross sections for the bloom ((a) without EMS, (b) 600 A, 2.0 Hz) and at the centerline of the computational outlets (c) in the two different casting cases, in which the cross sections of z = 0.234 m (center cross section of SEN outlets), z = 0.42 m( EMS center), z = 0.7 m (mold outlet) and z = 1 m are chosen to display more

Figure 13 .
Figure 13.Distributions of solute element C on the diagonal plane ((a) without EMS; (b) 600 A, 2 Hz) and at the diagonal line of the computational outlets (c).

Figure 13 .
Figure 13.Distributions of solute element C on the diagonal plane ((a) without EMS; (b) 600 A, 2 Hz) and at the diagonal line of the computational outlets (c).

Figure 14 .
Figure 14.Distributions of solute element C on the y-z plane, at various cross sections of the bloom (a,b), and the centerline of the computational outlets (c).

Figure 15 presents
the three dimensional level fluctuation in the mold under different current intensities.The fluctuation range of the liquid level in the mold rises with the increase of current intensity and the maximum site of the level fluctuation occurs near the four corners of the bloom, and the largest fluctuating values are 5.3 mm, 5.7 mm, 5.9 mm, and 6.2 mm, respectively.The fluctuation degree increases with the rise of current intensity.For big bloom castings, the value of level fluctuation should be controlled within around ±3 mm, and the maximum fluctuation of the four current intensities are close to this range.

Figure 14 .
Figure 14.Distributions of solute element C on the y-z plane, at various cross sections of the bloom (a,b), and the centerline of the computational outlets (c).

Figure 15 presents
the three dimensional level fluctuation in the mold under different current intensities.The fluctuation range of the liquid level in the mold rises with the increase of current intensity and the maximum site of the level fluctuation occurs near the four corners of the bloom, and the largest fluctuating values are 5.3 mm, 5.7 mm, 5.9 mm, and 6.2 mm, respectively.The fluctuation degree increases with the rise of current intensity.For big bloom castings, the value of level fluctuation should be controlled within around ±3 mm, and the maximum fluctuation of the four current intensities are close to this range.

Figure 15 .
Figure 15.Comparison of three dimensional level fluctuations under different current intensity.

Figure 16 .
Figure 16.Velocity vector (a) and distribution of element C (b) at EMS center.

Figure 15 .
Figure 15.Comparison of three dimensional level fluctuations under different current intensity.

Figure 15 .
Figure 15.Comparison of three dimensional level fluctuations under different current intensity.

Figure 16 .
Figure 16.Velocity vector (a) and distribution of element C (b) at EMS center.

Figure 16 .
Figure 16.Velocity vector (a) and distribution of element C (b) at EMS center.

Figure 17 .
Figure 17.Distribution of element C at computational outlet (a) and at centerline (b).

Figure 18
Figure18displays the variation of surface temperature and shell thickness along with casting direction at different current intensities.With the increase of current intensity (from 450 A to 550 A), the surface temperature increases as well, when the current intensity is over 550 A, the surface temperature does not result in any obvious change.The surface temperatures at the mold outlet (z = 700 mm) for each case are 1325.7 K for 450 A, 1343.3K for 500 A, 1359.4K for 550 A, and 1355 K for 600 A, respectively.Also, 1391.5 K, 1400.9K, 1405.8K, and 1403.1 K at the outlet of the secondary cooling zone I (z = 1.23 m), which then increase to 1444.2 K, 1450.7 K, 1447.2K and 1445 K respectively at the computational outlet due to the change of heat transfer conditions, respectively.With the increase of current intensity, the electromagnetic force rises, while the growth velocity of the solidifying shell in the EMS effect zone deceases.The solidified shell in the thickness direction at the mold outlet at four different current intensities is about 15.3 mm for 450 A, 15.2 mm for 500 A, 14.8 mm for 550 A, and 14.2 mm for 600 A, respectively, and increases to 27.55 mm, 27.2 mm, 27.57mm, and 26.86 mm at the outlet of computational zone (z = 1500 mm), respectively.

Figure 18 .
Figure 18.Variation of surface temperature and shell thickness along with casting direction at different current intensities.

Figure 17 .
Figure 17.Distribution of element C at computational outlet (a) and at centerline (b).

Figure 18 20 Figure 17 .
Figure18displays the variation of surface temperature and shell thickness along with casting direction at different current intensities.With the increase of current intensity (from 450 A to 550 A), the surface temperature increases as well, when the current intensity is over 550 A, the surface temperature does not result in any obvious change.The surface temperatures at the mold outlet (z = 700 mm) for each case are 1325.7 K for 450 A, 1343.3K for 500 A, 1359.4K for 550 A, and 1355 K for 600 A, respectively.Also, 1391.5 K, 1400.9K, 1405.8K, and 1403.1 K at the outlet of the secondary cooling zone I (z = 1.23 m), which then increase to 1444.2 K, 1450.7 K, 1447.2K and 1445 K respectively at the computational outlet due to the change of heat transfer conditions, respectively.With the increase of current intensity, the electromagnetic force rises, while the growth velocity of the solidifying shell in the EMS effect zone deceases.The solidified shell in the thickness direction at the mold outlet at four different current intensities is about 15.3 mm for 450 A, 15.2 mm for 500 A, 14.8 mm for 550 A, and 14.2 mm for 600 A, respectively, and increases to 27.55 mm, 27.2 mm, 27.57mm, and 26.86 mm at the outlet of computational zone (z = 1500 mm), respectively.

Figure 18 .
Figure 18.Variation of surface temperature and shell thickness along with casting direction at different current intensities.

Figure 18 .
Figure 18.Variation of surface temperature and shell thickness along with casting direction at different current intensities.

Table 1 .
Target chemical composition of U71Mn steel.Measured distribution of electromagnetic intensity in the mold under different current intensity.

Table 1 .
Target chemical composition of U71Mn steel.