Three-Dimensional Numerical Modeling of Macrosegregation in Continuously Cast Billets

Macrosegregation, serving as a major defect in billets, can severely degrade material homogeneity. Better understanding of the physical characteristics of macrosegregation through numerical simulation could significantly contribute to the segregation control. The main purpose of this study was to predict macrosegregation in continuously cast billets with a newly developed three-dimensional macrosegregation model. The fluid flow, solidification, and solute transport in the entire billet region were solved and analyzed. Flow patterns, revealing a typical melt recirculation at the upper region of mold and thermosolutal convection at the secondary cooling zone, significantly affect the solidification and solute distribution. The solute redistribution occurring with thermosolutal convection at the solidification front contributes significantly to continued macrosegregation as solidification proceeds. The results of this study show that the equilibrium partition coefficient is mostly responsible for the magnitude of macrosegregation, while comparison between solute P and S indicated that diffusion coefficients also have some amount of influence. Typical macrosegregation patterns containing a positively segregated peak at the centerline and negatively segregated minima at either side were obtained via the proposed three-dimensional macrosegregation model, which validated by the measured surface temperatures and segregation degree.


Introduction
Macrosegregation is one of the major defects in continuously cast billets, especially for high-carbon steel.It can severely degrade material homogeneity and thus affect the final properties and performance of billet-based products.Macrosegregation is also known as the main cause of cup fracture during wire cord drawing due to the formation of martensite and carbide networks [1].Macrosegregation cannot be eliminated entirely during heat treatment, so preventing it in the casting process is crucial [2].It is of considerable practical significance to analyze the solute distribution in castings to understand the basic physical mechanisms involved in macrosegregation.
Flemings et al. [3] first discovered the importance of convection within the mushy zone during alloy solidification and proposed the Local Solute Redistribution Equation (LSRE) to numerically study the macrosegregation caused by interdendritic flow.The flow of interdendritic liquid was modeled with a fixed dendritic solid network ignoring the effect of mushy zone.Mehrabian et al. [4] applied Darcy's law in Flemings' model by treating the mushy zone as a porous medium, and then established a new model for interdendritic flow behavior and the resulting macrosegregation during ingot solidification.An essential characteristic of solidification is the temperature field; the thermally induced flow exerts a significant effect on macrosegregation, but is simplified as a single where ρ and → u represent mixture density and velocity, respectively, and they are defined using mass fractions f n and volume fractions g n for solid (n = s) and liquid (n = l) phases.
However, as the densities of liquid and solid were assumed to be constant and equal, the volume fractions for solid and liquid are equal to the corresponding mass fractions.

Momentum Equation
where → u s represents the solid velocity, p is the pressure, and the effective viscosity µ e f f describing the sum of laminar and turbulent viscosity was defined to take into account the turbulence effect.The source term S u can be described as follows: where β T and β C,i represent the thermal and solutal expansion coefficients, respectively, T is temperature, C i is concentration of element i, the subscript R identifies the reference value, and ξ is a very small positive number.
The first term on the right-hand of Equation ( 5) accounts for the thermal and solutal buoyancy, which resulted from the temperature gradient and concentration gradient during solidification.The two buoyancy terms were calculated with the Boussinesq approximation.The second term represents the phase interaction force in the mushy zone, which was modeled as a porous medium wherein the Darcy's law can be applied.A m is permeability coefficient, which depends on the morphology of porous media.In this work, A m is calculated using the expression given by Minakawa [15] as A m = 180/λ 2  2 , where λ 2 represents the secondary dendrite arm spacing and the expression [20] used here will be presented in the species equation.

Turbulence Model
To consider the effects of turbulence on the fluid flow and macrosegregation in upper part of mold, the low-Reynolds number k − ε model [21,22] was adopted here, which considered to be preferable than standard two-equation k − ε model or high-Reynolds number k − ε model in continuous casting processes.Governing equations of the model can be written as follows: ) where k and ε are the turbulent kinetic energy and its rate of dissipation, respectively.S k and S ε represent the source terms for k and ε equations, which related to solidification as l +ξ A m k, and The turbulent viscosity can be calculated with the following equation: The coefficients and empirical constants [22] used for the Launder-Sharma low-Reynolds number model are listed in Table 1.
where k T,e f f is the effective heat conductivity, which can be described as follows: where k T,l and k T,s are the thermal conductivity of solid and liquid, respectively, and the turbulent Prandtl number Pr t is listed in Table 1.

of 22
The enthalpy H in the energy conversation equation can be described as a function of temperature: where c ps is the specific heat, H re f is the enthalpy at reference temperature T re f , and L is the latent heat of fusion.The f s − T curve is non-linear in the alloy solidification system.However, for the limitation of the existing models, coupling the non-linear f s − T function with the macrosegregation model is still challenging.Generally, a linear f s − T function [15,17] was assumed to calculate the liquid fraction: where T l and T s are the liquidus and solidus temperature respectively, which can be described as follows: where T f is the fusion temperature of pure iron, m i is the slope of liquidus line, and k p,i is the equilibrium partition coefficient of element i.

Species Equation
where D s,i and D l,i represent the solute diffusion coefficient of solid and liquid, respectively.Sc t is the turbulent Schmidt number, as listed in Table 1.
→ u s is regarded as a constant vector, casting speed here.The mixture species concentration is defined as follows: where α i is the solutal Fourier number, β i is solutal back-diffusion coefficient, α C is the additional term of Fourier number, and the value was set to be 0.1 according to the works of Voller and Beckermann [18].
where t f is local solidification time.X is the length of the microsegregation domain at the end of solidification, which usually is taken to be equal to half the final secondary dendrite arm spacing λ 2 .
Metals 2017, 7, 209 In the present work, the final secondary dendrite arm spacing, which varies with the cooling condition and alloy composition, was calculated using the empirical equations proposed by Won and Thomas [20].
where C R is the cooling rate, and C C is the carbon content.
The liquidus and solidus lines of phase diagram were assumed to be straight, which means that the equilibrium partition coefficient k p,i is a constant.Thus, the solute redistribution at the solid-liquid interface can be described as follows: The solute concentration in liquid phase is equal to that at the solid-liquid interface, due to the complete diffusion was assumed in the liquid.However, the solute concentration in the solid cannot be represented by the interfacial concentration.To obtain the solute concentration in solid phase during remelting, Schneider and Beckermann [23] assumed a constant average solid concentration with the Scheil-type model, and this method is analogous to having a microscopically well mixed solid during solidification.In this study, the following equation was derived according to the Equations ( 16) and (22) to calculate the solid concentration:

Reconstruction of Gravitational Acceleration
Continuous casting is a peculiarly complicated metallurgical process which occurs over a large region.To predict the macrosegregation in continuously cast billets, a vertical 3D geometric model (Figure 1a) was adopted to simulate the continuous casting process here.The sump bottom was determined to be about 9.4 m away from the meniscus according to the factory control software with a casting speed of 1.65 m/min.The size of the straight computational domain used here, to this effect, was designed to be 0.165 × 0.165 × 10 m 3 .The billet caster used in practice is a bow-type (10 m radius), and the billet curvature would influence the fluid flow and solute transport.Moreover, the major reason for liquid flow in the secondary cooling zone is reported to be thermosolutal convection, also called "double diffusive convection", which occurs as a result of the interaction between thermal and solutal gradient.Solute redistribution and solidification in the mushy zone result in thermosolutal buoyancy, which in turn affects the formation of further macrosegregation in the billet.As reflected in Equation ( 5), gravity plays an important role in the calculation of thermosolutal buoyancy.The influence of gravity varies with billet curvature, which can be described with the following variable θ.Thus, a reconstructed gravitational acceleration related to the variable θ was developed containing a x-component and a z-component according to the caster characteristics in this work.The reconstructed acceleration allows the consideration of billet curvature in the calculation of thermosolutal buoyancy with the straight computational domain shown in Figure 1a.
For the bow-type caster, the gravitational acceleration can be divided into radial → R g and tangential → T g components, as shown in Figure 1b.The components and variable θ can be calculated as follows: Considering the feature of straight computational domain and gravitational acceleration, two accelerations → A x and → A z are defined here: In the simulation, the decreasing direction of X was defined as the outer arc of billet caster.Then, the reconstructed gravitational acceleration applied in the calculation of thermosolutal buoyancy can be described as follows: Metals 2017, 7, 209 7 of 22 buoyancy, which in turn affects the formation of further macrosegregation in the billet.As reflected in Equation ( 5), gravity plays an important role in the calculation of thermosolutal buoyancy.The influence of gravity varies with billet curvature, which can be described with the following variable θ .Thus, a reconstructed gravitational acceleration related to the variable θ was developed containing a x-component and a z-component according to the caster characteristics in this work.The reconstructed acceleration allows the consideration of billet curvature in the calculation of thermosolutal buoyancy with the straight computational domain shown in Figure 1a.For the bow-type caster, the gravitational acceleration can be divided into radial  g R and tangential  g T components, as shown in Figure 1b.The components and variable θ can be calculated as follows: Considering the feature of straight computational domain and gravitational acceleration, two accelerations 9 8 0 0 0 0 9 8   .In the simulation, the decreasing direction of X was defined as the outer arc of billet caster.Then, the reconstructed gravitational acceleration applied in the calculation of thermosolutal buoyancy can be described as follows:

Numerical Scheme
The macrosegregation model was solved in the computational domain shown in Figure 1 with the appropriate initial conditions, boundary conditions, and thermophysical parameters in opensource software OpenFOAM; details of the conditions can be found in previous work [24].The highcarbon steel SWRH82B is used in this study, the chemical composition and corresponding diffusion coefficients are listed in Table 2 along with the equilibrium partition coefficients of each element, and the thermophysical parameters of the steel are listed in Table 3.The finite volume method (FVM) was applied to the set of macroscopic transport equations, providing temperature and concentration fields as input parameters for the microsegregation model.The dendritic structure, liquid fraction, and solute distribution at the solidification front calculated by the microscopic model were reapplied to the macrosegregation model to update the flow, temperature, and concentration field values.The large time-step transient solver, PIMPLE (a merged PISO-SIMPLE algorithm), was used to solve the pressure-velocity coupling; this algorithm has shown superior numerical stability in large time-steps

Numerical Scheme
The macrosegregation model was solved in the computational domain shown in Figure 1 with the appropriate initial conditions, boundary conditions, and thermophysical parameters in open-source software OpenFOAM; details of the conditions can be found in previous work [24].The high-carbon steel SWRH82B is used in this study, the chemical composition and corresponding diffusion coefficients are listed in Table 2 along with the equilibrium partition coefficients of each element, and the thermophysical parameters of the steel are listed in Table 3.The finite volume method (FVM) was applied to the set of macroscopic transport equations, providing temperature and concentration fields as input parameters for the microsegregation model.The dendritic structure, liquid fraction, and solute distribution at the solidification front calculated by the microscopic model were reapplied to the macrosegregation model to update the flow, temperature, and concentration field values.The large time-step transient solver, PIMPLE (a merged PISO-SIMPLE algorithm), was used to solve the pressure-velocity coupling; this algorithm has shown superior numerical stability in large time-steps calculations as compared to the Pressure Implicit Splitting Operator (PISO) algorithm [25], which is applied in the PIMPLE algorithm to rectify the pressure-velocity correction, while the Semi-Implicit Method for Pressure-Linked equations (SIMPLE) algorithm allowed us to calculate the pressure from velocity components on a grid by coupling Navier-Stokes equations with an iterative procedure.Both numerical stability and efficiency can be achieved with the PIMPLE algorithm.The initial time step was set at 0.001 s and was adapted gradually (increased to 0.023 s) in the solution.For each time step, the maximum iteration number was set to 600 to reduce the residual values of momentum conservation equation below 10 −6 , as well as the energy and species conservation equations below 10 −7 .About 160 h were required to perform a fully-coupled simulation with 1,587,600 grid cells on a workstation with 3.4 GHz CPU and 16 GB RAM in parallel.Grid independence analysis reported by Ge et al. [13] and Du [26] for the simulation of macrosegregation shows that the variation of grid size has little effect on the predicted segregation pattern.Thus, no grid dependence study was performed in this work.

Boundary Conditions
The boundary conditions considered applicable to the macrosegregation simulation in continuously cast billet are described as follows:

Meniscus
Considering the insulating effect of mold flux, the free surface was set to be adiabatic wall.The normal gradients of other variables were set to be zero.

Inlet and Outlet
Molten steel is poured into the mold through a straight nozzle at a certain flow rate corresponding to the casting speed.The inlet velocity profile was assumed to be flat and determined based on the inlet-outlet mass balance.Thus, boundary conditions of all variables at the inlet can be described as follows: where w b represents the billet width, d in is the SEN inner diameter, and T c is the pouring temperature.The values of k and ε at the SEN inlet were calculated on the basis of the semi-empirical equations [27].At outlet, fully developed condition was adopted for the fluid flow, and normal gradients of other variables were set to be zero.

Moving Wall
Velocity of the moving wall (billet surface) was set to equal casting speed, and the detailed thermal boundary conditions are given as follows: (a) Mold zone: where a is a constant, which value is set to be 2,680,000 W/m 2 , and b is an unknown coefficient: where L m is the mold length, and q is an unknown coefficient: where C W is the specific heat of water, Q W is the flow rate of cooling water, ∆T is the temperature difference of cooling water, and S e f f is effective area of mold.(b) Secondary cooling zone: where T sur f is the surface temperature of billet, T W is the temperature of cooling water, and h is the integrative heat transfer coefficient, which can be calculated using the following equations: Foot-roller zone [28]: where the constants: B = 392.5,c = 0.55, M = 0.0075, and T 0 = 273.15K. Other spray zones: where m, A, and n are constants, m = 116 W/(m 2 • K) , A = 10.44, and n = 0.851.(c) Air cooling zone: where σ is Stefan-Boltzmann constant, e is emissivity of billet surface, and T ambi is ambient temperature.

Results and Discussion
The 3D distribution of temperature, velocity, and solute element can be obtained by solving the governing equations.The flow field directly affects the solute convection, and the solidification determines the dendritic and liquid fraction distributions, then affects the redistribution of solute at solidification front, the results of solidification and flow field would be shown and discussed first.

Fluid Flow and Solidification
The predicted 3D streamline distribution, detailing the flow performance of molten steel in the mold, is illustrated in Figure 2 along with the velocity contour and its corresponding streamlines on the cross section of the Z = 0.4 m plane and diagonal section.The molten steel was poured into the mold through a submerged nozzle with a certain flow rate corresponding to the casting speed.The turbulent flow formed below the nozzle in response to the impinging jet of molten steel.A portion of the molten steel flowed downward into the secondary cooling zone along the caster, while the remainder moved upward along the solidified shell.This created a significant melt recirculation around the impinging jet at the upper region of the caster.The depth of the recirculating region was about 0.45 m, as dependent on the casting speed, mold size, and geometry of the nozzle.Another relatively weak recirculation zone appeared above the large melt recirculation due to the backflow of the molten steel.The cooperation between the backflow and meniscus created the melt recirculation below meniscus, which would cause serious level fluctuation and further increase the likelihood of powder entrapment during continuous casting.The velocity distribution on the cross section (Z = 0.4 m) indicated that larger velocity occurred at the corner compared to the edge near the solidification shell.This can also be reflected in the flow characteristics indicated with streamlines on the cross section shown in Figure 2. The introduced molten steel moved out from the billet center in all directions on the cross section, and the streams focused at the billet corner.The fluid flow along the centerline was transferred to the billet corner due to the solidified shell, enhancing fluid flow at the billet corner.Flow field is the key factor during continuous casting, and it will significantly affect the energy and solute transport.The predicted 3D flow characteristics in the continuously cast billet are crucial for accurate simulation of macrosegregation.
around the impinging jet at the upper region of the caster.The depth of the recirculating region was about 0.45 m, as dependent on the casting speed, mold size, and geometry of the nozzle.Another relatively weak recirculation zone appeared above the large melt recirculation due to the backflow of the molten steel.The cooperation between the backflow and meniscus created the melt recirculation below meniscus, which would cause serious level fluctuation and further increase the likelihood of powder entrapment during continuous casting.The velocity distribution on the cross section (Z = 0.4 m) indicated that larger velocity occurred at the corner compared to the edge near the solidification shell.This can also be reflected in the flow characteristics indicated with streamlines on the cross section shown in Figure 2. The introduced molten steel moved out from the billet center in all directions on the cross section, and the streams focused at the billet corner.The fluid flow along the centerline was transferred to the billet corner due to the solidified shell, enhancing fluid flow at the billet corner.Flow field is the key factor during continuous casting, and it will significantly affect the energy and solute transport.The predicted 3D flow characteristics in the continuously cast billet are crucial for accurate simulation of macrosegregation.Figure 3 shows the predicted solidification profiles, including both the center and corner of the billet surface, along the casting direction.Four surface temperatures that measured with pyrometers are also shown in the figure to validate the solidification model.In the experiment, four infrared thermometers were placed at the 2nd and 3rd zone of the second cooling zone and at two positions in the air-cooling zone.The predicted and experimental surface temperature values are in close agreement, which indicates that the solidification model used here is accurate.The simulated temperature profiles for both the surface center and corner also exhibited similar evolution characteristics: Both decrease sharply in the mold due to the high cooling intensity, then the decreased cooling intensity and release of latent heat create an increase in temperature in the secondary cooling zone.(The temperature either decreased or increased in the secondary cooling zone, but varied slightly in this region.)The difference between the temperature profiles mainly exists in the mold region, wherein the temperature at the surface corner fell below 800 K at the mold exit, while about 1200 K was obtained for the surface center.This can be attributed to the solidification characteristics of billet in the mold, wherein cooling at the corner is two-dimensional while unidimensional solidification is commonly assumed at the billet surface.The liquid fraction profile at the billet center is also shown in the figure as a function of the distance from the meniscus.The sump bottom was predicted to be 9.7 m.
in the mold region, wherein the temperature at the surface corner fell below 800 K at the mold exit, while about 1200 K was obtained for the surface center.This can be attributed to the solidification characteristics of billet in the mold, wherein cooling at the corner is two-dimensional while unidimensional solidification is commonly assumed at the billet surface.The liquid fraction profile at the billet center is also shown in the figure as a function of the distance from the meniscus.The sump bottom was predicted to be 9.7 m. Figure 4 shows the predicted liquid fraction distribution at different sections in the mold including the longitudinal section, diagonal section, and cross section of Z = 0.25 m and 0.4 m planes.As solidification proceeds, the solidified shell with a certain thickness appeared at the billet edge and then developed gradually as the distance to the meniscus increased.The simulated results indicate a thicker billet shell at the diagonal section in comparison to that at the longitudinal section.Distinctions in the liquid fraction distribution were directly revealed by the results at the cross section.About 10 mm in thickness is achieved for the distinction at the mold exit in the present work.The primary cause of that is the cooling characteristic of billet, as shown in Figure 3, which details the evolution of temperature at the billet corner and surface during continuous casting, revealing the distinctions in cooling intensity.There was a sharper variation in temperature at the surface corner compared to the center, which created a relatively larger growth rate of the solidified shell at the surface corner.The flow pattern shown in Figure 2   Figure 4 shows the predicted liquid fraction distribution at different sections in the mold including the longitudinal section, diagonal section, and cross section of Z = 0.25 m and 0.4 m planes.As solidification proceeds, the solidified shell with a certain thickness appeared at the billet edge and then developed gradually as the distance to the meniscus increased.The simulated results indicate a thicker billet shell at the diagonal section in comparison to that at the longitudinal section.Distinctions in the liquid fraction distribution were directly revealed by the results at the cross section.About 10 mm in thickness is achieved for the distinction at the mold exit in the present work.The primary cause of that is the cooling characteristic of billet, as shown in Figure 3, which details the evolution of temperature at the billet corner and surface during continuous casting, revealing the distinctions in cooling intensity.There was a sharper variation in temperature at the surface corner compared to the center, which created a relatively larger growth rate of the solidified shell at the surface corner.The flow pattern shown in Figure 2 has indicated relatively larger velocity at the billet corner, where the fluid flow accelerated the energy transport and enhanced the development of the solidified shell.The development of billet shell at the diagonal section is relatively uniform, while the solidified shell tends to expand outwards at an axial distance about 0.4 m away from the meniscus (Figure 4, part A).The enlarged plot of part A and cross sections evidently indicate a thicker shell at the Z = 0.25 m plane in comparison to the Z = 0.4 m plane, which appears to violate the growth tendency of the solidified shell.To further explore this abnormal finding, the liquid fraction profiles at different cross sections including Z = 0.25 m, 0.40 m, and 0.60 m were analyzed, as shown in Figure 5.The thickest solidified shell appears at the Z = 0.25 m plane among these three billet locations, although with the smallest distance from the meniscus.The abnormality existing in the distribution of the solidified shell is likely related to the flow pattern of molten steel in the mold, so the work zeroed in on the streamline distribution in part A as also shown in Figure 4 to explore the flow performance in this region.As shown, the melt recirculation, which would promote the energy dissipation by increasing the effective heat conductivity and convection, mainly occurs in the region of 0.2 to 0.4 m away from the meniscus.The solidified shell developed extensively and quickly in this region in response to the increased energy dissipation, while a portion of the molten steel with large velocity flowed toward the solidified shell at the end of the impinging jet (i.e., the source of the circulation flow), leading to "washing effect".This caused a decrease and uneven distribution of the solidified shell, as observed in the cross section of Z = 0.4 m.
dissipation by increasing the effective heat conductivity and convection, mainly occurs in the region of 0.2 to 0.4 m away from the meniscus.The solidified shell developed extensively and quickly in this region in response to the increased energy dissipation, while a portion of the molten steel with large velocity flowed toward the solidified shell at the end of the impinging jet (i.e., the source of the circulation flow), leading to "washing effect".This caused a decrease and uneven distribution of the solidified shell, as observed in the cross section of Z = 0.4 m.    performance in this region.As shown, the melt recirculation, which would promote the energy dissipation by increasing the effective heat conductivity and convection, mainly occurs in the region of 0.2 to 0.4 m away from the meniscus.The solidified shell developed extensively and quickly in this region in response to the increased energy dissipation, while a portion of the molten steel with large velocity flowed toward the solidified shell at the end of the impinging jet (i.e., the source of the circulation flow), leading to "washing effect".This caused a decrease and uneven distribution of the solidified shell, as observed in the cross section of Z = 0.4 m.Thereby some measures such as soft-reduction and electromagnetic stirring can be implemented in this region to reduce or suppress the defects.The streamline distribution at each cross section is also shown in Figure 6 to illustrate the melt circulation in flow zone, which created an asymmetrical distribution of the velocity and liquid fraction.This also exerted some influence over the solute distribution, as discussed below.The formed melt circulation shown in Figure 6, which is commonly called "thermosolutal convection" [29], was believed to be induced by thermal and solutal buoyancy.
section decreased from 0.045 m/s to 0.029 m/s as the distance to the meniscus increased; further, it was nearly equal to the casting speed (0.0275 m/s).The maximum liquid fraction decreased from 0.8 to 0.28, and the cross-flow zone gradually shrunk as the solidification proceeds.The limited space and fluid flow shown in the figure corresponding to Z = 7 m plane were responsible for internal flaws such as segregation and porosity.Thereby some measures such as soft-reduction and electromagnetic stirring can be implemented in this region to reduce or suppress the defects.The streamline distribution at each cross section is also shown in Figure 6 to illustrate the melt circulation in flow zone, which created an asymmetrical distribution of the velocity and liquid fraction.This also exerted some influence over the solute distribution, as discussed below.The formed melt circulation shown in Figure 6, which is commonly called "thermosolutal convection" [29], was believed to be induced by thermal and solutal buoyancy.

Solute Distribution and Validation
Involving the effect of thermosolutal buoyancy, solute diffusion, and convection during continuous casting process, the macrosegregation model developed here has predicted the 3D distribution of solute elements in the billet.Figure 7 shows the predicted distribution of different solute elements in mold ((a) C; (b) P; (c) S; (d) Mn), which is qualitatively similar for each element.The various elements encountered similar convection from the fluid flow and were redistributed at the solidification front.The mass fraction of elements in solid and liquid, which is controlled by the partition coefficient, is the only difference between them.For the sake of brevity, this study focuses only on the solute C to analyze the evolution of macrosegregation in the billet.

Solute Distribution and Validation
Involving the effect of thermosolutal buoyancy, solute diffusion, and convection during continuous casting process, the macrosegregation model developed here has predicted the 3D distribution of solute elements in the billet.Figure 7 shows the predicted distribution of different solute elements in mold ((a) C; (b) P; (c) S; (d) Mn), which is qualitatively similar for each element.The various elements encountered similar convection from the fluid flow and were redistributed at the solidification front.The mass fraction of elements in solid and liquid, which is controlled by the partition coefficient, is the only difference between them.For the sake of brevity, this study focuses only on the solute C to analyze the evolution of macrosegregation in the billet.
The distribution pattern of solute C, as shown in Figure 8, confirms that solute redistribution occurred in the mold due to the convection, diffusion, and solidification after the molten steel was poured into mold through the submerged nozzle.Some amount of inhomogeneity in solute distribution was observed, and this is the early stage of macrosegregation in the billet.Solidification plays an important role in solute redistribution; the predicted solidified shell in the mold is also shown in the figure as a red line for comparison.The mass fraction of solute C in the solid zone was lower than the initial value (0.81), i.e., underwent "negative segregation", while aggregation of the solute element in the mushy zone resulted in "positive segregation" (Figure 9).The predicted profiles for the concentration of solute C and the liquid fraction at different distances to the meniscus are also shown in Figure 9.The primary cause of the segregation is believed to be solute redistribution occurring at the solidification front as solidification proceeds.The redistribution of each solute element is mainly dependent on its partition coefficient, which is usually less than 1. (Normally, Metals 2017, 7, 209 14 of 22 there is severer segregation for smaller partition coefficients.)The equilibrium partition coefficient in Equation (22) indicates that solute elements are rejected from the solidified shell into the mushy zone, leading to compositional heterogeneity in the billet.The predicted results also show a relatively high concentration of solute C in the upper region of the mold, which can be attributed to solute redistribution and melt convection in the mold.The flow patterns shown in Figure 2 suggested that a part of the molten steel moved upward along the billet shell, forming a pair of large recirculation zones in response to the impinging jet pouring from the submerged nozzle.Consequently, the solute-enriched liquid adjacent to the billet shell would rise back to the upper zone of the mold, resulting in the high concentration zone, also leading to the negative segregation region near the billet shell.The distribution pattern of solute C, as shown in Figure 8, confirms that solute redistribution occurred in the mold due to the convection, diffusion, and solidification after the molten steel was poured into mold through the submerged nozzle.Some amount of inhomogeneity in solute distribution was observed, and this is the early stage of macrosegregation in the billet.Solidification plays an important role in solute redistribution; the predicted solidified shell in the mold is also shown in the figure as a red line for comparison.The mass fraction of solute C in the solid zone was lower than the initial value (0.81), i.e., underwent "negative segregation", while aggregation of the solute element in the mushy zone resulted in "positive segregation" (Figure 9).The predicted profiles for the concentration of solute C and the liquid fraction at different distances to the meniscus are also shown in Figure 9.The primary cause of the segregation is believed to be solute redistribution occurring at the solidification front as solidification proceeds.The redistribution of each solute element is mainly dependent on its partition coefficient, which is usually less than 1. (Normally, there is severer segregation for smaller partition coefficients.)The equilibrium partition coefficient in Equation (22) indicates that solute elements are rejected from the solidified shell into the mushy zone, leading to compositional heterogeneity in the billet.The predicted results also show a relatively high concentration of solute C in the upper region of the mold, which can be attributed to solute redistribution and melt convection in the mold.The flow patterns shown in Figure 2 suggested that a part of the molten steel moved upward along the billet shell, forming a pair of large recirculation zones in response to the impinging jet pouring from the submerged nozzle.Consequently, the soluteenriched liquid adjacent to the billet shell would rise back to the upper zone of the mold, resulting in the high concentration zone, also leading to the negative segregation region near the billet shell.
The solute distribution changed again below the back-flow region, where positive segregation occurred in the mushy zone adjacent to the solidified shell.The prediction results also indicate a solute-enriched region at the distance of about 0.1 to 0.2 m from the meniscus.To explain this finding, the streamlines distribution in this region is also presented in Figure 8. Local positive segregation was likely a result of recirculation occurring at this region, wherein a relatively high solidification rate was achieved, followed by the mushy zone; the solute-enriched liquid was concentrated at this region as-affected by the recirculation.The "washing effect" [15,17] also created negative segregation in the mushy zone at a distance of 0.3 m from the meniscus, as shown in Figure 9a; positive segregation at the Z = 0.18 m plane and negative segregation at the 0.3 m plane were clearly identified    The solute distribution changed again below the back-flow region, where positive segregation occurred in the mushy zone adjacent to the solidified shell.The prediction results also indicate a solute-enriched region at the distance of about 0.1 to 0.2 m from the meniscus.To explain this finding, the streamlines distribution in this region is also presented in Figure 8. Local positive segregation was likely a result of recirculation occurring at this region, wherein a relatively high solidification rate was achieved, followed by the mushy zone; the solute-enriched liquid was concentrated at this region as-affected by the recirculation.The "washing effect" [15,17] also created negative segregation in the mushy zone at a distance of 0.3 m from the meniscus, as shown in Figure 9a; positive segregation at the Z = 0.18 m plane and negative segregation at the 0.3 m plane were clearly identified with nearly the same liquid fraction distribution.Figure 8    The segregation degree is defined to directly describe the solute distribution: where, i r is segregation degree of element i , and C 0 represents the initial mass fraction.
Figure 10 shows the segregation degree of solute C near the end of solidification as predicted to be 9.7 m. Figure 10a shows contour plots of the segregation distribution at longitudinal and diagonal sections, and Figure 10b shows the segregation profiles at two distances (Z = 9.1 m and 9.9 m) from the meniscus along the X-axis and X-Y direction, respectively.There was negligible variation in the solute distribution when the ongoing solidification nearly complete.The segregation profiles showed no evident distinction between the two plotted locations for both sections, although the central liquid fraction was predicted to be about 0.2 at the Z = 9.1 m plane.This suggests that macrosegregation in the billet formed in the mushy zone with relatively low central liquid fraction, but not the region near the sump bottom and this corresponds with the conclusion reported in previous study [24].Comparison between the segregation profiles of the two sections (longitudinal and diagonal) reveals the similar features in the distribution of solute C. The simulation results indicate that the highest segregation degree (positive segregation) occurred at the central region of billet, while relatively low values of segregation degree (negative segregation) are predicted at either side of the center.This precisely corresponds to the reported typical distribution of centerline segregation with a positively segregated peak at the centerline and negatively segregated minima at the sides [30].The segregation degree is defined to directly describe the solute distribution: where, r i is segregation degree of element i, and C 0,i represents the initial mass fraction.
Figure 10 shows the segregation degree of solute C near the end of solidification as predicted to be 9.7 m. Figure 10a shows contour plots of the segregation distribution at longitudinal and diagonal sections, and Figure 10b shows the segregation profiles at two distances (Z = 9.1 m and 9.9 m) from the meniscus along the X-axis and X-Y direction, respectively.There was negligible variation in the solute distribution when the ongoing solidification nearly complete.The segregation profiles showed no evident distinction between the two plotted locations for both sections, although the central liquid fraction was predicted to be about 0.2 at the Z = 9.1 m plane.This suggests that macrosegregation in the billet formed in the mushy zone with relatively low central liquid fraction, but not the region near the sump bottom and this corresponds with the conclusion reported in previous study [24].Comparison between the segregation profiles of the two sections (longitudinal and diagonal) reveals the similar features in the distribution of solute C. The simulation results indicate that the highest segregation degree (positive segregation) occurred at the central region of billet, while relatively low values of segregation degree (negative segregation) are predicted at either side of the center.This precisely corresponds to the reported typical distribution of centerline segregation with a positively segregated peak at the centerline and negatively segregated minima at the sides [30].As described in Equation ( 12), a linear s f T − function was applied in the calculation of solid fraction in this work.The nonlinear s f T − function corresponding to the Voller-Beckermann microsegregation model should be given as [31]: Obviously, the linear s f T − relation ignores the influence of the partition coefficients and back- diffusion coefficients, which are difficult to determine in the multicomponent alloy system (values for the coefficients vary with different elements).The distinction between the calculated solid fraction and theoretical value would cause some effects on the magnitude of segregation.To discuss this, other two two-dimensional calculations [24] were performed with the linear s f T − function (Equation ( 12)) and nonlinear s f T − relation (Equation (36)) separately.In the calculation, the partition coefficient in Equation ( 36) is determined to be 0.34, and the back-diffusion coefficient is set to be 0.9 according to the calculated values of solute C. As noted already, the coefficients are different for various solute elements.Thus, the discussion only focuses on the predicted results of solute C, which are shown in the Figure 11 with the liquid fraction profiles: (a) Z = 4 m; and (b) Z = 9.9 m.Comparison between the two calculations shows a minor distinction in the liquid fraction at mushy zone, which distinction leads to the difference in the predicted concentration of solute C.However, the influence of solid fraction calculation on the final macrosegregation patterns is shown to be small.Thus, the linear s f T − function applied in this work is acceptable.Even so, the present method of calculating the solid fraction in Voller-Beckermann microsegregation model cannot be considered realistic, and much future work is required to develop a more accurate method.As described in Equation ( 12), a linear f s − T function was applied in the calculation of solid fraction in this work.The nonlinear f s − T function corresponding to the Voller-Beckermann microsegregation model should be given as [31]: Obviously, the linear f s − T relation ignores the influence of the partition coefficients and back-diffusion coefficients, which are difficult to determine in the multicomponent alloy system (values for the coefficients vary with different elements).The distinction between the calculated solid fraction and theoretical value would cause some effects on the magnitude of segregation.To discuss this, other two two-dimensional calculations [24] were performed with the linear f s − T function (Equation ( 12)) and nonlinear f s − T relation (Equation (36)) separately.In the calculation, the partition coefficient in Equation ( 36) is determined to be 0.34, and the back-diffusion coefficient is set to be 0.9 according to the calculated values of solute C. As noted already, the coefficients are different for various solute elements.Thus, the discussion only focuses on the predicted results of solute C, which are shown in the Figure 11 with the liquid fraction profiles: (a) Z = 4 m; and (b) Z = 9.9 m.Comparison between the two calculations shows a minor distinction in the liquid fraction at mushy zone, which distinction leads to the difference in the predicted concentration of solute C.However, the influence of solid fraction calculation on the final macrosegregation patterns is shown to be small.Thus, the linear f s − T function applied in this work is acceptable.Even so, the present method of calculating the solid fraction in Voller-Beckermann microsegregation model cannot be considered realistic, and much future work is required to develop a more accurate method.To validate the 3D macrosegregation model, solute C across the cross section of experimental billet was analyzed with a carbon-sulfur analyzer.Specimens for chemical analysis were extracted by drilling from the billet along the diagonal line; the detailed scheme of the specimens is shown in Figure 12a.The drilling specimens were 5 mm in diameter and 10 mm in depth, with a distance between specimens of 15 mm. Figure 12b gives the comparison in segregation degree profiles between the experiment and simulation, which were in close agreement, confirming the validity of the 3D macrosegregation model developed in this study.

Evolution of Macrosegregation
The computed distribution pattern of solute C on the cross section at different locations of the billet, including 1, 2, 4, 6, 8, 9, and 9.9 m, is shown in Figure 13 to reveal the evolution of macrosegregation during billet continuous casting.The distribution of solute C at the billet edge was relatively uniform and barely changed in comparison to the initial distribution.This can be attributed to the high cooling rate at the billet edge, wherein the solidification was too fast for the solute redistribution.As solidification proceeds, the decreased cooling rate facilitated solute redistribution at the solidification front leading to a low solute concentration in the resulting shell and increased mass fraction in the mushy zone, as shown in the figure corresponding to the Z = 1 m plane.The ringshaped negative segregation zone formed in response to the washing effect in the mold persisted at Z = 9.9 m plane, where accomplished solidification was achieved in the billet.This can be attributed to the solute diffusion coefficient in the solid, as described in Table 2.The limited diffusion of solute elements in the solid did not readily allow variation in the solute distribution at macroscopic scale, and this corresponds exactly to the previous report [30] that macrosegregation in cast products cannot be eliminated by subsequent hot working or annealing.The basic cause of segregation in the resulting To validate the 3D macrosegregation model, solute C across the cross section of experimental billet was analyzed with a carbon-sulfur analyzer.Specimens for chemical analysis were extracted by drilling from the billet along the diagonal line; the detailed scheme of the specimens is shown in Figure 12a.The drilling specimens were 5 mm in diameter and 10 mm in depth, with a distance between specimens of 15 mm. Figure 12b gives the comparison in segregation degree profiles between the experiment and simulation, which were in close agreement, confirming the validity of the 3D macrosegregation model developed in this study.To validate the 3D macrosegregation model, solute C across the cross section of experimental billet was analyzed with a carbon-sulfur analyzer.Specimens for chemical analysis were extracted by drilling from the billet along the diagonal line; the detailed scheme of the specimens is shown in Figure 12a.The drilling specimens were 5 mm in diameter and 10 mm in depth, with a distance between specimens of 15 mm. Figure 12b gives the comparison in segregation degree profiles between the experiment and simulation, which were in close agreement, confirming the validity of the 3D macrosegregation model developed in this study.

Evolution of Macrosegregation
The computed distribution pattern of solute C on the cross section at different locations of the billet, including 1, 2, 4, 6, 8, 9, and 9.9 m, is shown in Figure 13 to reveal the evolution of macrosegregation during billet continuous casting.The distribution of solute C at the billet edge was relatively uniform and barely changed in comparison to the initial distribution.This can be attributed to the high cooling rate at the billet edge, wherein the solidification was too fast for the solute redistribution.As solidification proceeds, the decreased cooling rate facilitated solute redistribution at the solidification front leading to a low solute concentration in the resulting shell and increased mass fraction in the mushy zone, as shown in the figure corresponding to the Z = 1 m plane.The ringshaped negative segregation zone formed in response to the washing effect in the mold persisted at Z = 9.9 m plane, where accomplished solidification was achieved in the billet.This can be attributed to the solute diffusion coefficient in the solid, as described in Table 2.The limited diffusion of solute elements in the solid did not readily allow variation in the solute distribution at macroscopic scale, and this corresponds exactly to the previous report [30] that macrosegregation in cast products cannot be eliminated by subsequent hot working or annealing.The basic cause of segregation in the resulting

Evolution of Macrosegregation
The computed distribution pattern of solute C on the cross section at different locations of the billet, including 1, 2, 4, 6, 8, 9, and 9.9 m, is shown in Figure 13 to reveal the evolution of macrosegregation during billet continuous casting.The distribution of solute C at the billet edge was relatively uniform and barely changed in comparison to the initial distribution.This can be attributed to the high cooling rate at the billet edge, wherein the solidification was too fast for the solute redistribution.As solidification proceeds, the decreased cooling rate facilitated solute redistribution at the solidification front leading to a low solute concentration in the resulting shell and increased mass fraction in the mushy zone, as shown in the figure corresponding to the Z = 1 m plane.The ring-shaped negative segregation zone formed in response to the washing effect in the mold persisted at Z = 9.9 m plane, where accomplished solidification was achieved in the billet.This can be attributed to the solute diffusion coefficient in the solid, as described in Table 2.The limited diffusion of solute elements in the solid did not readily allow variation in the solute distribution at macroscopic scale, and this corresponds exactly to the previous report [30] that macrosegregation in cast products cannot be eliminated by subsequent hot working or annealing.The basic cause of segregation in the resulting billets is that the liquid rejects the solutes during solidification because the solid has less solubility for them as compared to the liquid.The solute rejection was simulated as solute redistribution at the solidification front controlled by the microsegregation model in this work.A solute-enriched liquid formed near the solidified shell as solidification progressed, while the resulting localized solute enrichment caused macrosegregation due to the physical movement of liquid.The thermosolutal buoyancy is believed to play an important role in this movement.The streamline patterns shown in the figure clearly indicate the ongoing thermosolutal convection of molten steel in the billet, which resulted in an uneven, large-scale, and complex solute distribution.The solute C was not equally distributed in the cross sections due to thermosolutal buoyancy and gravity.As solidification progressed, solute elements gradually gathered at the central region of the billet as the mobile molten steel decreased.Macrosegregation formed completely in response to the interaction among solidification, solute redistribution, and thermosolutal convection.It should be noted that remelting and sedimentation of grains will affect the patterns and severity of segregation, especially in continuously cast bloom and slab.Thus, new mathematical model need to be developed considering the sedimentation of grains to more accurately predict the segregation patterns in bloom and slab.

Comparison between Different Elements
Although with similar feature of segregation profile, evident distinctions in the magnitudes of segregation have been obtained through the 3D macrosegregation model for the different solute elements.Among the four typical chemical alloy elements shown in Figure 14, solute P yielded the highest segregation degree at up to 1.920.Solute Mn showed the lowest segregation degree (1.099) and solute C's was 1.322.The equilibrium partition coefficient of various elements is the main contributor to these distinctions.Macrosegregation results from movements of micro-segregated regions over macroscopic distances due to the motion of liquid, while microsegregation arises as a

Comparison between Different Elements
Although with similar feature of segregation profile, evident distinctions in the magnitudes of segregation have been obtained through the 3D macrosegregation model for the different solute elements.Among the four typical chemical alloy elements shown in Figure 14, solute P yielded the highest segregation degree at up to 1.920.Solute Mn showed the lowest segregation degree (1.099) and solute C's was 1.322.The equilibrium partition coefficient of various elements is the main contributor to these distinctions.Macrosegregation results from movements of micro-segregated regions over macroscopic distances due to the motion of liquid, while microsegregation arises as a result of the solute redistribution at solidification front, which is mainly controlled by the equilibrium partition coefficient.The smaller is the partition coefficient, the more disparate is the mass fraction in the liquid and solid and thus the severer is the microsegregation.Although the formation of macrosegregation is primarily affected by the natural convection due to thermal and solutal buoyancy, and partly relevant to the diffusion coefficient of solute elements, the solute redistribution occurring at the solidification front is the determining factor in the magnitudes of the macrosegregation.A smaller equilibrium partition coefficient was adopted for the solute P compared to the solute C and Mn, so severer macrosegregation was predicted for solute P.However, although with smaller equilibrium partition coefficient compared to the solute P, solute S still achieved a lower degree of segregation, due to the distinction in diffusion coefficients between the two elements.Figure 15 shows the variations in diffusion coefficients with temperature for the two elements.The diffusion coefficient of solute S in the solid is much larger than that of solute P, especially at the temperature near the solidus.Increase in diffusion coefficient contributed to the homogenization of the solute elements during solidification, where the resulting microsegregation was relatively weak according to the microsegregation model.Thus, the predicted distinctions in the segregation degree of solutes P and S are reasonable.macrosegregation is primarily affected by the natural convection due to thermal and solutal buoyancy, and partly relevant to the diffusion coefficient of solute elements, the solute redistribution occurring at the solidification front is the determining factor in the magnitudes of the macrosegregation.A smaller equilibrium partition coefficient was adopted for the solute P compared to the solute C and Mn, so severer macrosegregation was predicted for solute P.However, although with smaller equilibrium partition coefficient compared to the solute P, solute S still achieved a lower degree of segregation, due to the distinction in diffusion coefficients between the two elements.
Figure 15 shows the variations in diffusion coefficients with temperature for the two elements.The diffusion coefficient of solute S in the solid is much larger than that of solute P, especially at the temperature near the solidus.Increase in diffusion coefficient contributed to the homogenization of the solute elements during solidification, where the resulting microsegregation was relatively weak according to the microsegregation model.Thus, the predicted distinctions in the segregation degree of solutes P and S are reasonable.macrosegregation is primarily affected by the natural convection due to thermal and solutal buoyancy, and partly relevant to the diffusion coefficient of solute elements, the solute redistribution occurring at the solidification front is the determining factor in the magnitudes of the macrosegregation.A smaller equilibrium partition coefficient was adopted for the solute P compared to the solute C and Mn, so severer macrosegregation was predicted for solute P.However, although with smaller equilibrium partition coefficient compared to the solute P, solute S still achieved a lower degree of segregation, due to the distinction in diffusion coefficients between the two elements.
Figure 15 shows the variations in diffusion coefficients with temperature for the two elements.The diffusion coefficient of solute S in the solid is much larger than that of solute P, especially at the temperature near the solidus.Increase in diffusion coefficient contributed to the homogenization of the solute elements during solidification, where the resulting microsegregation was relatively weak according to the microsegregation model.Thus, the predicted distinctions in the segregation degree of solutes P and S are reasonable.

Conclusions
A 3D macrosegregation model containing both the low Reynolds number k − ε turbulence model and Voller-Beckermann microsegregation model was developed in this study based on the continuum model proposed by Bennon and Incropera.The solidification, fluid flow, and solute transport in the entire billet region between the meniscus and sump bottom were analyzed using this model to reach the following main conclusions.
(1) Flow patterns showed a significant melt recirculation around the impinging jet at the upper region of the mold.Two other relatively weak melt recirculation zones appeared in response to the backflow of molten steel.There was larger velocity at the corner than the edge near the solidification shell.
(2) The solidification model applied in this work was validated according to surface temperature measurements.The relatively high cooling intensity created a thicker solidified shell at the diagonal section compared to the longitudinal section.Melt recirculation and washing effects further created uneven distribution of the solidified shell at the longitudinal section.The simulation results indicated thermosolutal convection at the secondary cooling zone, which resulted in asymmetrical distribution of the velocity and liquid fraction.(3) As convection, diffusion, and solidification progressed, there was measurable inhomogeneity in the solute distribution in the mold mainly occurring as solute enrichment at the upper region and negative segregation in the solidified shell.The distribution of different solute elements was predicted to be qualitatively similar.(4) The results also suggest that the final macrosegregation patterns in the billet formed in the mushy zone with relatively low central liquid fraction, but not at the end of solidification.The segregation profiles of solute C, arising with a positively segregated peak at the centerline and a negatively segregated minima at the sides, were predicted to be similar at the longitudinal and diagonal sections.There was close agreement between the measured and calculated segregation degrees, which supports the validity of the proposed 3D macrosegregation model.(5) The predicted ring-shaped negative segregation zone that formed in the mold still existed at the end of solidification due to the limited solute diffusion coefficient in solid.The solute redistribution occurring at the solidification front created solute-enriched liquid near the solidified shell, which contributed to the evolution of macrosegregation alongside thermosolutal convection.(6) The equilibrium partition coefficient was the main contributor to the distinctions in the magnitude of macrosegregation.Comparison between solutes P and S indicated that diffusion coefficients also have a certain influence on macrosegregation.

Figure 1 .
Figure 1.Vertical mesh adopted here (a); and reconstruction of gravitational acceleration (b).

Figure 1 .
Figure 1.Vertical mesh adopted here (a); and reconstruction of gravitational acceleration (b).

Figure 2 .
Figure 2. Streamline distribution and velocity contour in the mold.Figure 2. Streamline distribution and velocity contour in the mold.

Figure 2 .
Figure 2. Streamline distribution and velocity contour in the mold.Figure 2. Streamline distribution and velocity contour in the mold.

Figure 3 .
Figure 3. Evolution of solidification profiles along the casting direction.
Figure4shows the predicted liquid fraction distribution at different sections in the mold including the longitudinal section, diagonal section, and cross section of Z = 0.25 m and 0.4 m planes.As solidification proceeds, the solidified shell with a certain thickness appeared at the billet edge and then developed gradually as the distance to the meniscus increased.The simulated results indicate a thicker billet shell at the diagonal section in comparison to that at the longitudinal section.Distinctions in the liquid fraction distribution were directly revealed by the results at the cross section.About 10 mm in thickness is achieved for the distinction at the mold exit in the present work.The primary cause of that is the cooling characteristic of billet, as shown in Figure3, which details the evolution of temperature at the billet corner and surface during continuous casting, revealing the distinctions in cooling intensity.There was a sharper variation in temperature at the surface corner compared to the center, which created a relatively larger growth rate of the solidified shell at the surface corner.The flow pattern shown in Figure2has indicated relatively larger velocity at the billet corner, where the fluid flow accelerated the energy transport and enhanced the development of the solidified shell.The development of billet shell at the diagonal section is relatively uniform, while the solidified shell tends to expand outwards at an axial distance about 0.4 m away from the meniscus (Figure 4, part A).The enlarged plot of part A and cross sections evidently indicate a thicker shell at the Z = 0.25 m plane in comparison to the Z = 0.4 m plane, which appears to violate the growth tendency of the solidified shell.To further explore this abnormal finding, the liquid fraction profiles at different cross sections including Z = 0.25 m, 0.40 m, and 0.60 m were analyzed, as shown in Figure 5.The thickest solidified shell appears at the Z = 0.25 m plane among these three billet locations,

Figure 3 .
Figure 3. Evolution of solidification profiles along the casting direction.

Figure 4 .
Figure 4. Contour plots of liquid fraction at different sections.

Figure 5 .
Figure 5. Predicted liquid fraction profiles at different cross sections.

Figure 6
Figure 6 presents the contour plots of velocity and liquid fraction at several cross sections with different distances (Z = 3 m, 5 m, 7 m, and 9 m) to the meniscus.The plots consist of both the velocity and liquid fraction at the same cross section, because the distribution of both variables was symmetrical along the Y-axis.The upper part of each plot shows the velocity distribution and the bottom details the liquid fraction.The simulation results showed consistency in the distribution of

Figure 4 .
Figure 4. Contour plots of liquid fraction at different sections.

Figure 4 .
Figure 4. Contour plots of liquid fraction at different sections.

Figure 5 .
Figure 5. Predicted liquid fraction profiles at different cross sections.

Figure 6
Figure 6 presents the contour plots of velocity and liquid fraction at several cross sections with different distances (Z = 3 m, 5 m, 7 m, and 9 m) to the meniscus.The plots consist of both the velocity and liquid fraction at the same cross section, because the distribution of both variables was symmetrical along the Y-axis.The upper part of each plot shows the velocity distribution and the bottom details the liquid fraction.The simulation results showed consistency in the distribution of

Figure 5 .
Figure 5. Predicted liquid fraction profiles at different cross sections.

Figure 6
Figure6presents the contour plots of velocity and liquid fraction at several cross sections with different distances (Z = 3 m, 5 m, 7 m, and 9 m) to the meniscus.The plots consist of both the velocity and liquid fraction at the same cross section, because the distribution of both variables was symmetrical along the Y-axis.The upper part of each plot shows the velocity distribution and the bottom details the liquid fraction.The simulation results showed consistency in the distribution of velocity and liquid fraction at the same location in the billet.The maximum velocity at the cross section decreased from 0.045 m/s to 0.029 m/s as the distance to the meniscus increased; further, it was nearly equal to the casting speed (0.0275 m/s).The maximum liquid fraction decreased from 0.8 to 0.28, and the cross-flow zone gradually shrunk as the solidification proceeds.The limited space and fluid flow shown in the figure corresponding to Z = 7 m plane were responsible for internal flaws such as segregation and porosity.Thereby some measures such as soft-reduction and electromagnetic stirring can be implemented in this region to reduce or suppress the defects.The streamline distribution at each cross section is also shown in Figure6to illustrate the melt circulation in flow zone, which created an

Figure 6 .
Figure 6.Contour plots of velocity and liquid fraction at different cross sections.

Figure 6 .
Figure 6.Contour plots of velocity and liquid fraction at different cross sections.

Figure 8 .
Figure 8. Detailed distribution of solute C in the mold.

Figure 9 .
Figure 9. Predicted profiles for the concentration of solute C and liquid fraction: (a) Z = 0.18 and 0.30

Figure 8 .
Figure 8. Detailed distribution of solute C in the mold.
also shows the distribution of solute C at different cross sections in the mold.The corresponding predicted concentration and liquid fraction profiles are shown in Figure9b.There were similar distributions of solute C for different cross sections below the turbulent region in the mold.The segregation zone, both positive and negative, gradually moved toward the billet center as solidification progressed.

Figure 8 .
Figure 8. Detailed distribution of solute C in the mold.

Figure 9 .
Figure 9. Predicted profiles for the concentration of solute C and liquid fraction: (a) Z = 0.18 and 0.30 m; and (b) Z = 0.5 and 1.0 m.

Figure 9 .
Figure 9. Predicted profiles for the concentration of solute C and liquid fraction: (a) Z = 0.18 and 0.30 m; and (b) Z = 0.5 and 1.0 m.

Figure 10 .
Figure 10.Segregation degree of solute C near the end of solidification: (a) contour plots; and (b) segregation profiles.

Figure 10 .
Figure 10.Segregation degree of solute C near the end of solidification: (a) contour plots; and (b) segregation profiles.

Figure 11 .
Figure 11.Comparison between the results calculated with linear and nonlinear fs-T functions: (a) Z = 4 m; and (b) Z = 9.9 m.

Figure 12 .
Figure 12.Validation of the model: (a) sampling scheme in experimental billet; and (b) comparison between the calculated and experimental segregation degree.

Figure 11 .
Figure 11.Comparison between the results calculated with linear and nonlinear fs-T functions: (a) Z = 4 m; and (b) Z = 9.9 m.

Figure 12 .
Figure 12.Validation of the model: (a) sampling scheme in experimental billet; and (b) comparison between the calculated and experimental segregation degree.

Figure 12 .
Figure 12.Validation of the model: (a) sampling scheme in experimental billet; and (b) comparison between the calculated and experimental segregation degree.
the liquid rejects the solutes during solidification because the solid has less solubility for them as compared to the liquid.The solute rejection was simulated as solute redistribution at the solidification front controlled by the microsegregation model in this work.A solute-enriched liquid formed near the solidified shell as solidification progressed, while the resulting localized solute enrichment caused macrosegregation due to the physical movement of liquid.The thermosolutal buoyancy is believed to play an important role in this movement.The streamline patterns shown in the figure clearly indicate the ongoing thermosolutal convection of molten steel in the billet, which resulted in an uneven, large-scale, and complex solute distribution.The solute C was not equally distributed in the cross sections due to thermosolutal buoyancy and gravity.As solidification progressed, solute elements gradually gathered at the central region of the billet as the mobile molten steel decreased.Macrosegregation formed completely in response to the interaction among solidification, solute redistribution, and thermosolutal convection.It should be noted that remelting and sedimentation of grains will affect the patterns and severity of segregation, especially in continuously cast bloom and slab.Thus, new mathematical model need to be developed considering the sedimentation of grains to more accurately predict the segregation patterns in bloom and slab.

Figure 13 .
Figure 13.Computed distribution of solute C at different billet cross sections.

Figure 13 .
Figure 13.Computed distribution of solute C at different billet cross sections.

Figure 14 .
Figure 14.Comparison in the segregation profiles between different elements: (a) solutes C and Mn; and (b) solutes P and S.

Figure 15 .
Figure 15.Diffusion coefficients in both the liquid and solid for the solutes P and S: (a) solid zone; and (b) liquid zone.

A
3D macrosegregation model containing both the low Reynolds number − k ε turbulence model and Voller-Beckermann microsegregation model was developed in this study based on the continuum model proposed by Bennon and Incropera.The solidification, fluid flow, and solute transport in the entire billet region between the meniscus and sump bottom were analyzed using this model to reach the following main conclusions.

Figure 14 .
Figure 14.Comparison in the segregation profiles between different elements: (a) solutes C and Mn; and (b) solutes P and S.

Figure 14 .
Figure 14.Comparison in the segregation profiles between different elements: (a) solutes C and Mn; and (b) solutes P and S.

Figure 15 .
Figure 15.Diffusion coefficients in both the liquid and solid for the solutes P and S: (a) solid zone; and (b) liquid zone.

A
3D macrosegregation model containing both the low Reynolds number − k ε turbulence model and Voller-Beckermann microsegregation model was developed in this study based on the continuum model proposed by Bennon and Incropera.The solidification, fluid flow, and solute transport in the entire billet region between the meniscus and sump bottom were analyzed using this

Figure 15 .
Figure 15.Diffusion coefficients in both the liquid and solid for the solutes P and S: (a) solid zone; and (b) liquid zone.

Table 1 .
Coefficients and empirical constants used for the low-Reynolds number model.

Table 2 .
Chemical compositions and corresponding coefficients.

Table 3 .
Operation parameters and thermophysical parameters.