A New Approach to Calculate the Velocity of Interdendritic Fluid Flow during Solidification Using Etched Surface Height of Actual Metal Ingot

Macrosegregation remains one of main defects affecting metal materials properties, which is mainly caused by interdendritic fluid flow during solidifying. However, as for controlling actual specific segregation, it is still difficult to effectively measure or simulate this kind flow instead of pure fluid flow, especially in complex casting processes of high-grade materials. Herein, a new method for obtaining velocity magnitude and direction of interdendritic fluid flow during metal solidifying is proposed from boundary layer and standard deviation obtained by measuring etched surface heights of the actual ingot and using statistical principles. Taking continuous casting bloom of GCr15 bearing steel as an example, it is indicated that the calculated velocity magnitudes under different sides and superheats can be explained by process features and, hence, solidification mechanism. The velocity magnitude and fluctuation are higher on the inner curve side and under low superheat. Meanwhile, it is found that the fluctuation extent of secondary arm spacing is more relevant with interdendritic fluid flow, although its magnitude is mainly determined by the cooling rate. Moreover, on the basis of the calculated velocity directions and magnitudes, there is a positive correlation between segregation area ratio and the effective ratio between interdendritic flow velocity and growth velocity especially in the equiaxed grain zone, which corresponds with classic macrosegregation formation theory. The above findings and comparison with other results demonstrate the validity of the new approach, which can obtain the magnitude and the direction of interdendritic fluid velocity for two or three-dimensional multiscale velocity distribution by tailoring measuring length and numbers.


Introduction
Macrosegregations are very typical defects in metal materials now, which can cause performance failures due to nonuniform composition distribution and structure [1][2][3]. Regarding macrosegregation, the fundamental formation theory was pioneered by Fleimings and co-workers about 50 years ago and is now extensively accepted [4]. In the classical theory, microsegregation and interdendritic fluid flow are two essential aspects for macrosegregation formation, and the latter usually plays a decisive role during the process [5,6]. It can safely be said that all types of macrosegregation form within the liquid-solid zone, and it is the result of slow interdendritic flow in most cases, but there is much left to study to fully understand and control macrosegregation in practice [4,6]. Hence, numerous studies have attempted to measure, simulate, and control different characteristics of interdendritic liquid flow for specific processes [7]. For example, because the velocity magnitude can probably be calculated by Darcy's Law if the permeability is obtained, many researchers are interested in the experimental measurement of permeability and mathematical models of how permeability varies with fraction solid, columnar or equiaxed structures, different processing, etc. [8][9][10][11][12][13][14]. However, because the structures are somewhat conceptual and the pressure gradients are assumed or approximate, the calculated velocities always deviate from the real. At the same time, flow velocity in the pure liquid zone is always simulated by many researchers, but it is still very difficult to correctly simulate the flow velocity in the mushy zone, especially for complex fluid flow or large-size casting or with different actual dendritic grains of complex morphology [2,6,7,15,16]. Additionally, if based on the actual result (like the segregation content ratio) and local solute redistribution equation [5,6], the velocity magnitude can be calculated by at least two contents, but it is only a relative value with respect to a reference content that may not be available to compare (e.g., original liquid content that ignores element diffusion or liquid transfer) and is always a cubic average value with a large scale (e.g., centimeter) for the shortcomings of content measuring methods. Therefore, as to the complex solidifying process, it is always difficult to investigate precisely the characteristics of interdendritic liquid flow for the formation mechanism of different actual segregation morphology.
Solute redistribution is the result of different kinds of segregation, especially for some elements with low distribution coefficients, such as carbon in steel. So it should be useful to obtain firstly the actual element distribution in order to study the specific segregation formation mechanism. However, because some complicated phenomena like eutectic reaction occur in the final solidification, it is always not easy to measure directly the element change from the casting dendrite grain center to its margin (almost always at least 0.5 mm long) by existing technologies (e.g., wet chemistry, optical emission spectrometry, X-ray spectroscopy or fluorescence) [17], not to mention many casting grains. Therefore, on the basis of the nonuniform distribution, hot pickling is always used to etch the samples in order to reveal the structure and related element distribution [18][19][20][21]. By different grayscales, different structures (e.g., columnar, equiaxed, macrosegregation) can be distinguished from the etched structure picture [15,[20][21][22]. Meanwhile, the height of the sample surface becomes uneven after being etched because the electrochemical reaction of hot pickling appears at different levels or kinds in different zones of smooth surface. For example, different surface heights can induce different grayscales, and different grayscales can be indicative of segregated carbon element content in steel alloy [15,22,23], so the height on the sample surface after etching can also probably be indicative of carbon element content. Moreover, the grayscale value can be more easily influenced by measuring processes such as light source features, and the large area and multiscale distribution of element content is very difficult to measure directly by existing methods, so the etched surface height could be a possible method to study the element redistribution phenomenon in a large area. Although an entirely accurate quantitative correlation equation between the etched surface height and the element content at a certain position is now difficult to establish, plenty of etched surface heights (abbreviated as h es in this article) may be partly suitable for investigating some segregation formation processes by simultaneously using some statistical methods.
In this study, the h es values at different surface places are measured firstly in GCr15 bearing steel bloom of continuous casting for the reason that its high-content carbon element (mass fraction, 1.02%) is easy to segregate during solidifying and the fluid flow is strong enough during the continuous casting process [15,19]. Secondly, the magnitude and the direction of velocity in interdendritic fluid flow (abbreviated as V i f f in this article) are calculated by establishing corresponding models on the basis of the h es . Finally, the results and the verification are discussed. In all, this study demonstrates a new approach for measuring the velocity of interdendritic fluid flow during metal solidification, which may provide a unique alternative method for understanding the interdendritic fluid flow, macrosegregation formation, and hence, material properties. After solidifying, the sample of GCr15 steel was obtained by cutting along the transverse section of the bloom. Then, the sample surface was polished and smooth as a mirror.

Sample Preparation
In order to get the es h , which is caused by different carbon contents, the commonly used hot pickling experiment with 50% hydrochloric acid-water solution (the bath temperature is 60-80 °C and the etching time is 25 min) [19,21] was carried out. The surface was cleaned at once by swiping in warm water and blow-dried quickly after removing it from the pickling solution. The sample etched casting macrostructures are shown in Figure 2a,c, which are mainly comprised of casting columnar and equiaxed dendritic grains (white part). The bloom casting processes are shown in Table 1; the only difference between 1# sample and 2# sample is the casting superheat, which is measured in the tundish. The superheat of 1# sample (20 °C) was lower than that of 2# sample (35 °C). After solidifying, the sample of GCr15 steel was obtained by cutting along the transverse section of the bloom. Then, the sample surface was polished and smooth as a mirror. In order to get the h es , which is caused by different carbon contents, the commonly used hot pickling experiment with 50% hydrochloric acid-water solution (the bath temperature is 60-80 • C and the etching time is 25 min) [19,21] was carried out. The surface was cleaned at once by swiping in warm water and blow-dried quickly after removing it from the pickling solution. The sample etched casting macrostructures are shown in Figure 2a,c, which are mainly comprised of casting columnar and equiaxed dendritic grains (white part). The bloom casting processes are shown in Table 1; the only difference between 1# sample and 2# sample is the casting superheat, which is measured in the tundish. The superheat of 1# sample (20 • C) was lower than that of 2# sample (35 • C). The es h was measured by a 3D Measuring Laser Microscope produced by Ol company (LEXT OLS4000, Tokyo, Japan). Due to the microscope limiting, the recta areas (12.5 mm × 220 mm) from the inner curve side to the outer curve side we chosen to implement measuring, as shown in Figure 2b,d.
Taking the inner curve side as an example, line sweeping measurement w formed at 13 positions (A to M in Figure 3) from down to up, as shown by the v dashed arrow line in Figure 3. The h es was measured by a 3D Measuring Laser Microscope produced by Olympus company (LEXT OLS4000, Tokyo, Japan). Due to the microscope limiting, the rectangular areas (12.5 mm × 220 mm) from the inner curve side to the outer curve side were first chosen to implement measuring, as shown in Figure 2b,d.
Taking the inner curve side as an example, line sweeping measurement was performed at 13 positions (A to M in Figure 3) from down to up, as shown by the vertical dashed arrow line in Figure 3.  The es h was measured by a 3D Measuring Laser Microscope produced by Olympus company (LEXT OLS4000, Tokyo, Japan). Due to the microscope limiting, the rectangular areas (12.5 mm × 220 mm) from the inner curve side to the outer curve side were first chosen to implement measuring, as shown in Figure 2b,d.
Taking the inner curve side as an example, line sweeping measurement was performed at 13 positions (A to M in Figure 3) from down to up, as shown by the vertical dashed arrow line in Figure 3.    The distances to the bloom edges at A-M were 1, 3, 5, 10, 20, 30, 40, 50, 60, 70, 80, 90 and 100 mm, respectively. Thus, as to each bloom, there were 26 positions of line sweeping for both inner curve and out curve sides. Because the microscope can measure a height every 0.25 µm, the total measuring points at a position (the length is 12.5 mm) were 50,001, which was enough for statistical calculation. The measurement in Figure 3 is for velocity magnitude calculation. For velocity direction, 1 mm long measuring was conducted at the left side and the right side of each position, as shown by the horizontal dashed arrow lines (α and β) from bloom edge toward center in Figure 4 where the vertical dashed arrow corresponds to the same arrow in Figure 3.  The distances to the bloom edges at A-M were 1, 3, 5, 10, 20, 30, 40, 50, 60, 70, 80, 90  and 100 mm, respectively. Thus, as to each bloom, there were 26 positions of line sweeping  for both inner curve and out curve sides. Because the microscope can measure a height  every 0.25 μm, the total measuring points at a position (the length is 12.5 mm) were 50,001, which was enough for statistical calculation. The measurement in Figure 3 is for velocity magnitude calculation. For velocity direction, 1 mm long measuring was conducted at the left side and the right side of each position, as shown by the horizontal dashed arrow lines (α and β) from bloom edge toward center in Figure 4 where the vertical dashed arrow corresponds to the same arrow in Figure 3.

Etched Surface Height
After measuring, the surface heights at A position in the outer curve side under 35 °C superheat were obtained as shown in Figure 5a. In addition to the height fluctuation, we found that the heights decreased linearly with increasing the distance as a whole (marked schematically by the inclined arrow line in Figure 5a). In fact, this distinct whole decreasing trend could have been caused by the original uneven surface before hot pickling, and it is very difficult to remove completely with the existing polishing technology. In order to eliminate the unexpected influence of the original unevenness, the cutoff length λc = 2.5 mm was chosen to filter the original data on the basis of the surface roughness theory [24,25]. The corresponding filtered heights at A position are shown in Figure  5b, which is characteristic of the almost flat basic surface marked schematically by a horizontal arrow line. Therefore, all heights at different positions were filtered before model construction.   Figure 3).

Etched Surface Height
After measuring, the surface heights at A position in the outer curve side under 35 • C superheat were obtained as shown in Figure 5a. In addition to the height fluctuation, we found that the heights decreased linearly with increasing the distance as a whole (marked schematically by the inclined arrow line in Figure 5a). In fact, this distinct whole decreasing trend could have been caused by the original uneven surface before hot pickling, and it is very difficult to remove completely with the existing polishing technology. In order to eliminate the unexpected influence of the original unevenness, the cutoff length λc = 2.5 mm was chosen to filter the original data on the basis of the surface roughness theory [24,25]. The corresponding filtered heights at A position are shown in Figure 5b, which is characteristic of the almost flat basic surface marked schematically by a horizontal arrow line. Therefore, all heights at different positions were filtered before model construction.            Figure 9a shows the macrostructure of the red rectangular area from Figure 4. The white part displays the dendrite grain morphology, and the black part displays the interdendritic area, which is characteristic of high element content. The horizontal dashed arrow corresponds to the vertical dashed arrow in Figure 3. Figure Figure 9a shows the macrostructure of the red rectangular area from Figure 4. The white part displays the dendrite grain morphology, and the black part displays the interdendritic area, which is characteristic of high element content. The horizontal dashed arrow corresponds to the vertical dashed arrow in Figure 3. Figure Figure 9a shows the macrostructure of the red rectangular area from Figure 4. The white part displays the dendrite grain morphology, and the black part displays the interdendritic area, which is characteristic of high element content. The horizontal dashed arrow corresponds to the vertical dashed arrow in Figure 3.   In addition, the heights may seem not to be larger or smaller when the point is closer to or farther from the grain center. This is because the surface macrostructure is photographed by a SELP1650 digital camera, which cannot display three-dimensional morphology (i.e., the height difference) and some dendrites are cut off by cross section. Notwithstanding, the above discussion indicates that there is some correlation between the etched surface height and segregating element content. Carbon is the element segregated most in GCr15 bearing steel solidification [17,19]. Therefore, based on the characteristics of h es distribution, a quantitative correlation equation between h es and the dimensionless quantity of carbon content is put forward firstly, given by Equation (1). C d is the dimensionless quantity of carbon content. At a certain position, 50,001 effective carbon contents are rearranged in descending order. Then, if the 12.5 mm measuring line at each position is treated as a solidification unit for the reason that every point on a measuring line has almost the same solidifying condition, the rearranged contents can be displayed by liquid fraction as a horizontal axis from a statistical average angle, as shown in Figure 10.

Model Description
In addition, the heights may seem not to be larger or smaller when the poin to or farther from the grain center. This is because the surface macrostructure graphed by a SELP1650 digital camera, which cannot display three-dimensional ogy (i.e., the height difference) and some dendrites are cut off by cross section. standing, the above discussion indicates that there is some correlation between t surface height and segregating element content. Carbon is the element segregate GCr15 bearing steel solidification [17,19]. Therefore, based on the characterist distribution, a quantitative correlation equation between es h and the dim quantity of carbon content is put forward firstly, given by Equation (1). d C is t sionless quantity of carbon content. At a certain position, 50,001 effective carbon are rearranged in descending order. Then, if the 12.5 mm measuring line at eac is treated as a solidification unit for the reason that every point on a measurin almost the same solidifying condition, the rearranged contents can be displayed fraction as a horizontal axis from a statistical average angle, as shown in Figure   In order to analyze the interdendritic fluid flow of macrosegregation that m pens in the middle and late period, the period with 0 to 0.7 liquid fraction in Figu be the focus of research (e.g., Figure 11). According to the near-equilibrium sol model (NESM) where no diffusion occurs in solids and finite diffusion occurs [7,26,27], it is assumed that the distribution in Figure 11 should be expressed wel tion (3)  That is to say, the highest C d can be regarded the final solidifying point, and the lowest is the original solidifying point. The content in solidus increases when liquid fraction is decreasing.
where C d is the dimensionless quantity of carbon content at a point. h i es is the etched surface height at a point. h a es is the average height of the etched surface at each position. In order to analyze the interdendritic fluid flow of macrosegregation that mainly happens in the middle and late period, the period with 0 to 0.7 liquid fraction in Figure 10 will be the focus of research (e.g., Figure 11). According to the near-equilibrium solidification model (NESM) where no diffusion occurs in solids and finite diffusion occurs in liquids [7,26,27], it is assumed that the distribution in Figure 11 should be expressed well by Equation (3) at all positions. Equation (2) is obtained by using the Scheil equation form to effectively describe element distribution in the near-equilibrium solidification model [28,29].  δ γ + .
The difference of the distribution characteristics between the 0 to and the 0.7 to 1.0 liquid fraction may be attributed to rare interdendriti period with nucleation and original growth. In fact, it is not easy to ob ment distribution of a dendrite grain, especially for the distribution i solidification, because some more complicate reactions (e.g., eutectic r As C d can be considered as solid content C * S on the basis of Equation (1), Equation (3) can be obtained on the basis of Equation (2) and be used to fit the content distribution like that in Figure 10 using the Levenberg-Marquardt algorithm [30]. These fitting calculations were carried out using Microsoft Excel software.
where C * S is the solute content at the solid phase of the solidifying interface. k E is the effective solute redistribution coefficient at the solidifying interface in the NESM. C 0 is the original solute content before solidification. f L is the liquid fraction.
where µ = k E C d 0 , γ = k E − 1, and C d o is the original dimensional quantity of carbon content. By Equation (1), C d o should be equal to 1 under ideal conditions and is the calculated result of δ/(γ + 1).
The difference of the distribution characteristics between the 0 to 0.7 liquid fraction and the 0.7 to 1.0 liquid fraction may be attributed to rare interdendritic flow in the initial period with nucleation and original growth. In fact, it is not easy to obtain the whole element distribution of a dendrite grain, especially for the distribution in the final stage of solidification, because some more complicate reactions (e.g., eutectic reaction) will occur for multicomponent alloys such as steel. Moreover, only in the two-dimensional surface of the actual ingot, the morphology appearing in the surface may not be representative of the three-dimensional grains structure. Owing to the above reasons, 12.5 mm long is measured to obtain enough data, which is more than 10 times a dendrite's width. Thus, even though the surface morphology information displays incompletely for the whole grain structures at some points, the statistically average features of a rearranged distribution of a great many heights could represent effectively the segregating element content distribution of the position by eliminating the deviation of original data as far as possible.
All fitted results at each positions are shown in Tables 2-5. The determination coefficient R 2 is calculated as Equation (4) [30]. When R 2 is closer to 1, the fitted formula can express the original data more effectively. Because all R 2 are almost bigger than 0.9 (more than half of them bigger than 0.95), it means that Equation (3) is a good formula to quantitatively describe the content distribution with 0 to 0.7 liquid fraction at all positions.
where y i is the original value. y is the average of the original value.ŷ i is the corresponding regression value.  (3), and R 2 is the corresponding determination coefficient. k E is the effective solute redistribution coefficient. M y is growth velocity magnitude in the solidification unit along the measuring line at each position. δ Ca is the average thickness of the concentration boundary layer).  (3), and R 2 is the corresponding determination coefficient. k E is the effective solute redistribution coefficient. M y is growth velocity magnitude in the solidification unit along the measuring line at each position. δ Ca is the average thickness of the concentration boundary layer).   (3), and R 2 is the corresponding determination coefficient. k E is the effective solute redistribution coefficient. M y is growth velocity magnitude in the solidification unit along the measuring line at each position. δ Ca is the average thickness of the concentration boundary layer). Based on fitted parameters and Equation (3), k E is obtained on the basis of γ value. Meanwhile, the fitted C d o is almost equal to 1 in Tables 2-5, which indicates that Equation (2) and the NESM are suitable to describe the distribution because the ideal C d o is equal to 1, further indicating that the above method is suitable to calculate the k E value.

Distance to
Meanwhile, according to the theory suggested by Burton and his co-workers [27,31], k E at the solidifying interface in the NESM can also be calculated theoretically by Equation (5). Then, the average thickness of concentration boundary layer δ Ca can be obtained by Equation (6). Then, where δ Ca is the average thickness of the concentration boundary layer. D is the mass diffusivity of the carbon element, 0.005 m 2 /s [32]. k 0 is the ideal solute redistribution coefficient and is chosen as 0.1 for the carbon element [33]. M y is the corresponding growth velocity magnitude in the solidification unit along the measuring line at each position, which can be calculated by Equation (7).
where L is the measuring length at each position, 12.5 mm. T l and T s are the liquidus temperature (1466 • C) and solidus temperature (1328 • C) of GCr15. R is the cooling rate ( • C/s), which is calculated by thermal simulation of the finite element model, and the governing equations are Equations (8) and (9) [33]. Axial heat conduction in the casting direction was ignored. ∂T ∂t where T is the temperature. t is time. λ is the thermal conductivity. ρ is the density. S h is the specific heat. The evolution of the latent heat during the solidification was incorporated into the calculation by using the effective specific heat method, as shown in Equation (9) [34].
where S h is the effective specific heat. H L is the latent heat and f s is the solid fraction.
Combining measured k E in Tables 2-5 with the aforementioned equations, δ Ca can be calculated as shown in Tables 2-5. Regardless of whether it is laminar flow or turbulent flow, there is a boundary layer during the initial stage that can be assumed to feature laminar flow characteristics, especially in the limited length on a flat plate [35]. Thus, in this research, the interdendritic flow on dendritic grain structure in each position can be considered as the flat plate flow model because there is only enough space to flow during the solidifying process. According to the dimensionless equation in the transport phenomena of the boundary layer theory [32,35], the average thickness of the velocity boundary layer δ Va can be calculated by Equation (10) after getting δ Ca .
where Sc is Schmidt number. µ is the dynamic viscosity, 0.0067 Pa·s, and D is the mass diffusivity of the carbon element, 0.005 m 2 /s, and ρ is the density, 7600 kg/m 3 [32]. Based on the boundary layer theory [34,35], the change of thickness of the velocity boundary layer δ V with velocity V and flow length φ can be described by Equation (11), the average thickness δ Va in φ long zone can be expressed by Equation (12). Then, From Equation (12), the equivalent flow velocity V δ Va can be calculated by Equation (13): where δ Va is obtained by Equation (10). φ is chosen as 1 mm based on the measuring size. Meanwhile, there is also the original thickness of the concentration boundary layer without fluid flow, which can be calculated by Equation (14), only considering molecular diffusion [36].
Based on Equation (14), when C L /C 0 = 1.01, the original thickness of the concentration boundary layer δ Ca 0 can be obtained by Equation (15). Then the original thickness of the velocity boundary layer δ Va 0 can also be calculated by Equation (10).
On the basis of δ Ca o and δ Ca , δ Va o and δ Va can be obtained. Then, the corresponding velocity can be calculated. Finally, the difference between V δ Va and V δ Vao is the magnitude of interdendritic fluid flow velocity V i f f , as shown in Equation (16). The corresponding flowchart is also shown in Figure 12. Following is the calculation method for the direction of interdendritic flow v with respect to the above iff V . In Figure 4, the measuring heights of 1 mm length α Following is the calculation method for the direction of interdendritic flow velocity with respect to the above V i f f . In Figure 4, the measuring heights of 1 mm length α and β should have the same features if there is no fluid flow because the left and right sides are almost under the same solidifying conditions for their limited distance. Then the standard deviation S d as shown in Equation (17) [30], is firstly introduced to measure the fluctuation extent of left or right side heights.
where h i es is the etched surface height at i point. h a es is the average value of the etched surface height. n is the number of the point.
where E VD is the difference value. S L d is the standard deviation of the left side, and S R d is the standard deviation of the right side.
As shown in Equation (18), the standard deviation difference between left side and right side is introduced to measure the velocity direction. As a matter of fact, because interdendritic fluid flow can induce more stochastic solidification behavior [7,36], it was assumed in this research that the main fluid flow direction is from right side to left side if E VD is positive, and the main fluid flow direction is from left side to right side if E VD is negative. The absolute value of E VD should indicate the extent along the growth direction. Actually, interdendritic fluid flow exists in three-dimensional space, and the calculated velocity magnitude above is also a three-dimensional value. Considering the surface from the inner curve side to outer curve side as x-positive axis, the grain growth direction from bloom left edge to center as y-positive axis, and the casting direction as z-positive axis, as shown in Figure 13, only the velocity direction component θ on xz plane can be obtained because only left and right side heights on the cross section are measured here.
where VD E is the difference value. L d S is the standard deviation of the left side, and R d S is the standard deviation of the right side.
As shown in Equation (18), the standard deviation difference between left side and right side is introduced to measure the velocity direction. As a matter of fact, because interdendritic fluid flow can induce more stochastic solidification behavior [7,36], it was assumed in this research that the main fluid flow direction is from right side to left side if  As to the degree θ , its value can be calculated by Equation (19), assuming the position with the largest absolute value of difference has the same or symmetric degree (0° or 180°) to grain growth direction. The degree calculation flowchart is also shown in Figure 14.  As to the degree θ, its value can be calculated by Equation (19), assuming the position with the largest absolute value of difference has the same or symmetric degree (0 • or 180 • ) to grain growth direction. The degree calculation flowchart is also shown in Figure 14.
where θ is the degree of fluid flow direction on the xz plane with respect to the x-positive axis in Figure 13, and arccos is the anticosine sign. E max VD is the standard deviation difference, which has the maximum absolute value among all points. The denominator in Equation (19) is the negative value of the absolute of E max VD .
Metals 2021, 11, x FOR PEER REVIEW 18 of where θ is the degree of fluid flow direction on the xz plane with respect to the x-positiv axis in Figure 13, and arccos is the anticosine sign.
max VD E is the standard deviation diffe ence, which has the maximum absolute value among all points. The denominator in Equa tion (19) is the negative value of the absolute of max VD E .

Results and Discussion
Based the above model, the magnitude and direction of velocity of interdendrit fluid during metal solidifying can be calculated. We will discuss these results under di ferent conditions in this section.  Table 6 and Figure 1 We found that the magnitude decreased from the bloom edge to center as a whole. Th high velocity magnitude near bloom edge could be attributed to the high flow velocit when the fluid moved from tundish to mold. Meanwhile, there was also electromagnet stirring in the initial solidification in the mold.

Results and Discussion
Based the above model, the magnitude and direction of velocity of interdendritic fluid during metal solidifying can be calculated. We will discuss these results under different conditions in this section.

Magnitude of Interdendritic Fluid Flow Velocity
The V i f f magnitudes at different positions shown in Figure 3 in inner and outer curve sides under different superheats (20 and 35 • C) are shown in Table 6 and Figure 15. We found that the magnitude decreased from the bloom edge to center as a whole. The high velocity magnitude near bloom edge could be attributed to the high flow velocity when the fluid moved from tundish to mold. Meanwhile, there was also electromagnetic stirring in the initial solidification in the mold.
Notwithstanding, in order to compare more clearly, each side was divided into three zones, i.e., chill zone, middle zone, and central zone in Figure 16. The chill zone contained three positions whose distances to the edge were 1, 3, and 5 mm, respectively. The middle zone contained five positions whose distances to the edge were 10, 20, 30, 40, and 50 mm, respectively. The central zone contained five positions whose distances to the edge were 60, 70, 80, 90, and 100 mm, respectively. Then we found that, a little different from the other zones, the velocity magnitude in the central zone decreased firstly, then showed an increasing trend and decreased again finally, as shown by arrow line in Figure 16c. This phenomenon was in correspondence with the extent of point and centerline segregation morphology of the central zone in Figure 2. The increasing trend could be attributed to the suction force to compensate the centerline porosity, i.e., the V segregation formation in continuous casting [6,7,37]. The reason for decreasing again could be M position being too close to the bloom center where solidifying completes quickly.  Notwithstanding, in order to compare more clearly, each side was divided into three zones, i.e., chill zone, middle zone, and central zone in Figure 16. The chill zone contained three positions whose distances to the edge were 1, 3, and 5 mm, respectively. The middle zone contained five positions whose distances to the edge were 10, 20, 30, 40, and 50 mm, respectively. The central zone contained five positions whose distances to the edge were 60, 70, 80, 90, and 100 mm, respectively. Then we found that, a little different from the other zones, the velocity magnitude in the central zone decreased firstly, then showed an increasing trend and decreased again finally, as shown by arrow line in Figure 16c. This phenomenon was in correspondence with the extent of point and centerline segregation morphology of the central zone in Figure 2. The increasing trend could be attributed to the suction force to compensate the centerline porosity, i.e., the V segregation formation in continuous casting [6,7,37]. The reason for decreasing again could be M position being too close to the bloom center where solidifying completes quickly. There was no apparent difference between different sides, but the corresponding average value displayed some regularity, if only for the middle and central zones. As shown in Figure 17a, not only under low superheat, but also under high superheat, both of the velocity magnitudes in the inner curve sides had a greater value than those in the outer curve sides. Meanwhile, as to the standard deviation calculated by Equation (16) Figure 15c).
There was no apparent difference between different sides, but the corresponding average value displayed some regularity, if only for the middle and central zones. As shown in Figure 17a, not only under low superheat, but also under high superheat, both of the velocity magnitudes in the inner curve sides had a greater value than those in the outer curve sides. Meanwhile, as to the standard deviation calculated by Equation (16) in Figure 17b, the values in the inner curve side also had a greater value than those in the outer curve side. These phenomena indicated that there was more active interdendritic fluid flow in the inner curve side by comparison with the outer curve side, which corresponded with the gravity features of the curved segment of the continuous casting machine in Figure 1. However, as to the chill zone, as shown in Figure 18, there was no same trend under different superheats, and the extent of changing proportion was relatively minor. This is because the chill zone almost solidified completely in the mold instead of the curved segment. Meanwhile, Figure 19a,b display related average value and standard deviation, respectively. We found that the velocity magnitude and related fluctuation were higher under 20 • C. There was no apparent difference between different sides, but the corresponding average value displayed some regularity, if only for the middle and central zones. As shown in Figure 17a, not only under low superheat, but also under high superheat, both of the velocity magnitudes in the inner curve sides had a greater value than those in the outer curve sides. Meanwhile, as to the standard deviation calculated by Equation (16) in Figure  17b, the values in the inner curve side also had a greater value than those in the outer curve side. These phenomena indicated that there was more active interdendritic fluid flow in the inner curve side by comparison with the outer curve side, which corresponded with the gravity features of the curved segment of the continuous casting machine in Figure 1. However, as to the chill zone, as shown in Figure 18, there was no same trend under different superheats, and the extent of changing proportion was relatively minor. This is because the chill zone almost solidified completely in the mold instead of the curved segment. Meanwhile, Figure 19a,b display related average value and standard deviation, respectively. We found that the velocity magnitude and related fluctuation were higher under 20 °C.    (20) where K is the permeability. ε is the volume fraction of interdendritic liquid. μ is the viscosity. P Δ is the pressure loss and L is the related distance.   (20) [38], besides the SDAS, which positively influences the permeability, the velocity magnitude is also determined by the pressure gradient ∆P/L in the mushy zone.
where K is the permeability. ε is the volume fraction of interdendritic liquid. µ is the viscosity. ∆P is the pressure loss and L is the related distance. Thus, the pressure magnitude should be greater under lower superheat, which has a narrower mushy zone owing to a shorter liquid core length in continuous casting bloom. Meanwhile, there was a positive relationship between the fluctuation of velocity magnitude and SDAS standard deviation, i.e., SDAS showed more fluctuant under 20 °C. Therefore, it was indicated that although the SDAS magnitude was mainly determined by the cooling rate, the fluctuation extent of SDAS was more relevant to interdendritic fluid flow by changing solute and thermal distribution, breaking dendrite, etc. Thus, the pressure magnitude should be greater under lower superheat, which has a narrower mushy zone owing to a shorter liquid core length in continuous casting bloom. Meanwhile, there was a positive relationship between the fluctuation of velocity magnitude and SDAS standard deviation, i.e., SDAS showed more fluctuant under 20 • C. Therefore, it was indicated that although the SDAS magnitude was mainly determined by the cooling rate, the fluctuation extent of SDAS was more relevant to interdendritic fluid flow by changing solute and thermal distribution, breaking dendrite, etc.

Direction of Interdendritic Fluid Flow Velocity
The differences of standard deviation between left side height and right side height E VD at different positions are shown in Figure 21. G position has the maximum absolute value of E VD among all points. Thus, the pressure magnitude should be greater under lower superheat, which has a narrower mushy zone owing to a shorter liquid core length in continuous casting bloom. Meanwhile, there was a positive relationship between the fluctuation of velocity magnitude and SDAS standard deviation, i.e., SDAS showed more fluctuant under 20 °C. Therefore, it was indicated that although the SDAS magnitude was mainly determined by the cooling rate, the fluctuation extent of SDAS was more relevant to interdendritic fluid flow by changing solute and thermal distribution, breaking dendrite, etc.

Direction of Interdendritic Fluid Flow Velocity
The differences of standard deviation between left side height and right side height VD E at different positions are shown in Figure 21. G position has the maximum absolute value of VD E among all points. Because the grain growth direction was from the edge to the center, which is the more negative, the fluid flow was closer to the grain grow direction based on the above model.
In Figure 22, the degree θ of fluid flow direction on the xz plane, as shown in Figure 13, and Equation (19) were obtained on the basis of assuming the position with the largest Because the grain growth direction was from the edge to the center, which is the more negative, the fluid flow was closer to the grain grow direction based on the above model. In Figure 22, the degree θ of fluid flow direction on the xz plane, as shown in Figure 13, and Equation (19) were obtained on the basis of assuming the position with the largest negative value of difference in Figure 21 (G position) had the same degree to grain growth direction, i.e., the degree was zero.  Because interdendritic fluid flow exists in three-dimensional space and the calculated velocity magnitude above is also a three-dimensional value, in order to compare them Because interdendritic fluid flow exists in three-dimensional space and the calculated velocity magnitude above is also a three-dimensional value, in order to compare them equally, the effective ratio between interdendritic flow velocity and growth velocity (V i f f /M x )e was calculated as Equation (21) which can contain the influence of velocity direction. ( where M x is the growth velocity whose direction is from the bloom edge toward bloom center, and the magnitude is equal to 2000 mm/R, R is the cooling rate ( • C/s), which was calculated by the above thermal simulation of finite element model. θ is the degree of fluid flow direction on the xz plane with respect to the x-positive axis. At the same time, the segregation area ratios at different positions were calculated by Equation (22) for the reason that it can reflect in part the related two-dimensional segregation extent caused by interdendritic fluid flow [21,39]. With the SELP1650 digital camera produced by Sony Corporation (Tokyo, Japan), the macrostructures whose locations are marked by a red rectangle in Figure 4 at 10 positions (B, C, D, E, F, G, H, I, J, K) in the inner curve side under 20 • C were obtained as shown in Figure 23. These macrostructures were used to calculate the corresponding segregation area ratios that display in Figure 24. B, C, D, J, and K mainly contained equiaxed grains, but E, F, G, H, and I mainly contained columnar grains.
where, R seg is the segregation area ratio, %, A t is the total analyzing area, 25 mm 2 , A seg is the area of segregation points (black part), mm 2 .   the macrosegregation formation theory proposed by Fleimings and co-workers [4][5][6]. As to the columnar grain positions marked by the dot-dash line, they also show a minor increasing trend if excluding E position marked by the full line. This may be because, compared with the dendritic columnar grain, the morphology in E position is a cell columnar grain, which grows under a higher temperature gradient. Meanwhile, the non-apparent increasing trend in the columnar grain zone may be attributed to its deviation from the assumption of the same solidifying condition for left side and right side at a certain position in the columnar grain zone of approximate directional solidification.  Figure 25 shows the relationship between R seg and (V i f f /M x )e. Firstly, the whole points do not display an apparent trend. However, if only analyzing the equiaxed grain positions (i.e., B, C, D, J, and K, which are marked by the dashed line.), we found that there was a positive correlation between R seg and (V i f f /M x )e which corresponded to the macrosegregation formation theory proposed by Fleimings and co-workers [4][5][6]. As to the columnar grain positions marked by the dot-dash line, they also show a minor increasing trend if excluding E position marked by the full line. This may be because, compared with the dendritic columnar grain, the morphology in E position is a cell columnar grain, which grows under a higher temperature gradient. Meanwhile, the non-apparent increasing trend in the columnar grain zone may be attributed to its deviation from the assumption of the same solidifying condition for left side and right side at a certain position in the columnar grain zone of approximate directional solidification.

Model Verification Discussion and Further Work
From the aforementioned research, by establishing the model, the magnitude and direction of velocity were obtained. As for the verification of the model, the following

Model Verification Discussion and Further Work
From the aforementioned research, by establishing the model, the magnitude and direction of velocity were obtained. As for the verification of the model, the following results should first be provided: (1) The δ/(γ + 1) values are almost equal to 1 in Tables 2-5, which indicates that the above method is suitable to calculate k E . (2) The comparison results of velocity magnitudes of different sides and under different superheats can be explained by continuous casting features and solidification theory. (3) There is a positive correlation between segregation area ratio and the effective ratio between interdendritic flow velocity and growth velocity, especially in the equiaxed grain zone.
Meanwhile, Takahashi et al. studied the relationship between interdendritic flow velocity and k E using a complex hot model experimental of steel ingot [40]. The fitted empirical formula is shown as Equation (23) [37,40]. By this equation, the k E is obtained on the basis of the calculated interdendritic flow velocity in this research. Figure 26 shows the empirical k E by Equation (23)  The growth velocity at the same distance to the bloom edge was considered as the approximate average value of different superheats. By comparing E k in Figure 26, we The growth velocity at the same distance to the bloom edge was considered as the approximate average value of different superheats. By comparing k E in Figure 26, we found that the measured value and the empirical value had the same change trend with increasing interdendritic flow velocity, and the values were very close, especially in high flow velocity (i.e., Figure 26a). Meanwhile, the differences between the measured value and the empirical value were more distinct when the velocity was decreasing. The reason for this could be that the empirical formula Equation (23) was obtained on the basis of experimental high flow velocity. Additionally, the above velocity magnitude in Table 6 was from 0.001 to 1 m/s, which was approximately half of pure fluid velocity in some other simulated works by rough contrast [7,[41][42][43]. In all, from the aforementioned discussion, what can be concluded is that the established method is available to obtain the magnitude and direction of interdendritic flow velocity.
where, Sh = (V i f f /M y )/(7500 + V i f f /M y ).
This work provides a feasible method to obtain simultaneously the magnitude and the direction of V i f f for two or three-dimensional multiscale velocity distribution by adjusting measuring length and numbers. Meanwhile, it should be noted that the calculated velocity in this paper is only a two-dimensional value, and it will need more related heights to calculate the whole three-dimensional velocity. However, for some specific segregation, the corresponding two-dimensional velocity distribution should be enough to resolve the problems. The above method was carried out under the assumption of the same solidifying condition of 12.5 mm length in the position line for magnitude and 1.0 mm length for direction. Thus enough attention should be paid to the scale condition according to specific solidification characteristics. In the future, more delicate simulation or experiment is necessary to modify the model, including getting a more accurate correlation between the etched surface height and element content at different points.

1.
We have put forward a full model that can calculate the magnitude and direction of interdendritic fluid flow during solidifying at different positions by measuring the corresponding etched surface heights of casting metal bloom.

2.
Calculated results show that the velocity magnitude decreases from continuous casting bloom margin toward the center as a whole, and it decreases in the central zone first, then shows an increasing trend, and decreases again finally, which is in correspondence with segregation morphology. Besides the chilled zone, the velocity magnitude of the interdendritic fluid flow in the inner curve side is higher than that in the outer curve side. The velocity magnitude and related fluctuation is higher under 20 • C. Moreover, it is indicated that the fluctuation extent of SDAS is more relevant with interdendritic fluid flow, although the SDAS magnitude is mainly determined by the cooling rate.

3.
Meanwhile, on the basis of the magnitude and direction of velocity, there is a positive correlation between segregation area ratio and the effective ratio between interdendritic flow velocity and growth velocity, especially in equiaxed grain zones.

4.
On the whole, the model is demonstrated as a useful approach for determining the velocity of interdendritic fluid flow. The velocity distribution and associated defects formation in two or three-dimensional multiscale space can be better understand based on this model.