Effects of Impeller Rotational Speed and Immersion Depth on Flow Pattern, Mixing and Interface Characteristics for Kanbara Reactors Using VOF-SMM Simulations

: The Kanbara Reactor (KR) is a primary desulfurization technology in the hot metal pretreatment reﬁning process that is widely employed in the modern steelmaking industry. The operating parameters of KR impeller immersion depth (IID) and rotation speed (IRS) have a crucial impact on the process performance and the desulfurization effect. Still, their inﬂuences have not been fully understood. This study systematically investigated the effects of IID and IRS on the ﬂow pattern, mixing behavior, vortex core depth, and free surface characteristics for KR processes based on a 3D Volume of Fluid (VOF) model coupled with the sliding mesh method (SMM). The model was validated via scale-down water model experiments and then applied to the KR process, and simulations found that IID and IRS have different impacts on the ﬂow pattern. Speciﬁcally, the discharge ﬂow location moves downward with IID increasing, but the discharge strength and mean velocity hardly changes. Comparatively, the rise of IRS signiﬁcantly increases the mean velocity, but few changes occur to the discharge ﬂow position. Increasing IRS improves bath hydrodynamics, strengthens recirculation, and efﬁciently shortens mixing time, but IID has a neglectable effect on these features. The minimum mixing time is 55 s at a maximum IRS of 260 rpm. Moreover, the vortex core depth and free surface velocity visibly increase with the increase of IRS. Comparatively, IID has a limited effect on the ﬂow and mixing behavior but directly impacts the distribution of recirculation regions at the axial direction and the velocity gradient on the free surface at the radial direction. Furthermore, the correlation equations of these critical parameters as a function of the operating parameters were obtained. The results from this study may provide references for operating optimizations and industrial practices of KRs.


Introduction
With the increasing demand for high-quality steel, cleanliness is of paramount importance in each stage of steel refining and production.Except for special applications, the sulfur existing in products is usually a detrimental impurity that adversely affects the temperature toughness, weldability, and hydrogen-induced cracking resistance during the service life of products [1,2].Thus, it is essential to reduce the sulfur content to as low as possible to satisfy the specification of the product.Desulfurization is usually carried out before the hot metal is mixed into the converter because of favorable thermodynamic conditions, and this stage is thus called the pretreatment procedure for desulfurization.Approaches used here can be classified into two categories: mechanically stirred and powder injection (PI).Among them, the Kanbara Reactor (KR) [3], as a typical mechanically stirred method (MSM) [4][5][6], has significant advantages with high efficiency, short treatment time, and lower cost, providing the intensely hydrodynamic conditions through the revolving motion of a submerged impeller.Practically, the rotational speed and immersion depth of the impeller, as two vital operating parameters, remarkably impact the entrainment of the reagent powder, mixing characteristics, and refining efficiency, as well as the service life of an impeller and the ladle because of higher refractory wear induced by the washing of revolving, high-temperature metal.
Some studies have been conducted on different aspects of MSM or KRs, including fluid flows, bath mixing behaviors, reagent dispersion, and impeller design based on scale-down water models, numerical simulations, and laboratory-scale high-temperature experiments.For example, Visuri et al. [1] have given a recent review in modeling hot metal desulfurization and summarized the previous contributions and work on fluid flow phenomena, mixing patterns, and particle/bubble dispersions in KR processes.Yamamoto et al. [7,8] investigated the flow characteristic of aluminum melt by mechanically stirring and found that the impeller depths, rotation speeds, and impeller twisting angles impact the mixing time.Kato et al. [9] evaluated the flow pattern and mass transfer in an aluminum melting furnace using a numerical simulation.They found that the region with high kinetic energy was enlarged, and the turbulent kinetic energy at the free surface increased significantly with impeller rotation speed.This result illustrated that increasing impeller rotation speed is an effective method to enhance the mass transfer behavior of reactors.Their more recent research [10] found that increasing rotation speed would increase vortex depth and the intensity of turbulent fluctuations on the free surface.He et al. [11] investigated the coalescence and sedimentation behavior of falling droplets in a vessel of MSM.They found that increasing rotation speed is conducive to the dispersion of droplets in the mechanical stirring vessel.Later, they [12] also investigated the dispersed behavior of desulfurization particles at different rotation speeds and found that rotation speed affects the liquid dispersion in the mechanical stirring vessel and has a noticeable effect on the distribution of particles in the hot metal.Li et al. [13,14] used a Volume of Fluid with a dispersed particle model (VOF-DPM) to analyze the drawdown mechanism and the dispersion characteristic of light particles at different rotation speeds.They found that a higher impeller rotation speed can provide more power consumption and increase particle velocity, which means the mixing efficiency would be better.Nakai et al. [15,16] described the variation of the gas-liquid interface at different stages of the mixing process in a 1/8 scale water model.They found that the vortex height and operation parameters, such as impeller immersion depths, rotation speeds, and the geometric size of KR, have a relationship.Ji et al. [17] investigated the mixing behavior of KR vessels using computational fluid dynamics (CFD) simulations.By comparing fixed-and variable-velocity stirring, they found the latter stirring mode to be more favorable in the mixing of a bath.Sahu et al. [18] used CFD simulation, which optimized the operating parameters in the synthesis of hybrid aluminum matrix composites by stir casting.They found that blade angle, impeller size, and stirring speed are the typical factors that affect the flow pattern.Liu et al. [19] investigated impeller immersion depths, rotation speeds, and eccentric mechanical stirring modes.They found that increasing rotation speed can effectively strengthen the shear effect around the impeller, which intensifies bubble disintegration and dispersion.With an increasing impeller immersion depth, the distance between the vortex and impeller rises gradually.All the related research above shows that the impeller immersion depth and rotation speed significantly influence the fluid flow and mixing efficiency.However, to date, only a few reports focus on this topic, and it is not fully understood or well documented due to the limits of the measurement method and high costs of in-situ experimentation at a high ambient temperature and within a hazardous environment.
This study aims to provide a detailed analysis of the effect of impeller operating variables on fluid flow pattern, mixing characteristics, vortex core depth, and free surface velocity of KRs based on the developed Volume of Fluid model with the sliding mesh method (VOF-SMM).This study extends our previous work [20] in which the effect of the Metals 2021, 11, 1596 3 of 18 impeller dimension was discussed.Here, our focus is the rotational speed and immersion depth of impellers.

Mathematical Model
Our model has been reported in recent studies on the effect of impeller dimension for the KR process [20] and on a new dephosphorization technology for the converter steelmaking process [21].In those studies, the grid sensitivity, calibration of the turbulence model, and model validation have been carefully examined, but considering some improvements and the completion of the present work, the outline is given as follows.

Governing Equations of VOF Model
To track the free surface of gas-liquid two-phase flow in a KR, some assumptions were specified, and details were given in our previous work [20].As a result, the hydrodynamics of this multiphase system is governed by a set of single Navier-Stokes equations and a continuity equation.
∂ρ ∂t Here, u i and u j are the velocity components in the direction of i and j, respectively (m•s −1 ); t is the time (s); p is the pressure (Pa); g is the gravity acceleration (m•s −2 ); and τ ij is the Reynold stress tensor (N•m −2 ), which is calculated as follows: where δ ij is the Kronecker delta (δ ij = 1 if I = j and δ ij = 0 if i = j); µ eff is the effective viscosity (Pa•s) calculated by: where µ is the molecular viscosity (Pa•s) and the turbulent viscosity µ t (Pa•s) is calculated by the turbulence model.
In the present study, a modified approach was used based on the multi-fluid formulation proposed by Ubbink and Isssa [22].The evolution of the interface/free surface [23] is described by an additional transport equation for the indicator function representing the volume fraction of one phase that needs to be solved together with the continuity and momentum equations: where, u i,r = u i,q − u i,p is the relative velocity (m•s −1 ) between phase q and phase p.
Here, the third term on the left hand of Eqn. ( 5) is added to bring in a counter-gradient artificial interface compression term in which the u i,r ensures compression with a suitable compression velocity, while the ∂/∂x i guarantees conservation and α q 1 − α q guarantees boundedness with the compression only in the interface region.
In each cell, the volume fractions of all phases sum to unity, i.e., The one-fluid formulation relies on the fact that multiple fluids (or phases) are not interpenetrating.Because the same immiscible fluids are considered as an effective fluid throughout the domain, the momentum equations shared by all phases are solved with an effective density ρ (kg•m −3 ) and an effective viscosity µ (Pa•s).The fluid properties are calculated by a weighted averaging method based on the volume fraction of each phase.The volume-fraction-averaged density and viscosity are calculated, respectively, by: In Equations ( 1)-( 8), the subscript q denotes the liquid and gas phases.
To consider the effect of surface tension, an additional term is included in the momentum equation as f σ on the interface S(t), calculated per unit volume using the continuum surface force model with the following expression: where σ is surface tension (N•m −1 ) and κ is the mean curvature of the free surface determined by:

Turbulence Model
We have tested four kinds of classic turbulence models, i.e., the standard k-ε, standard k-ε with swirl flow correction term, RNG k-ε, and SST k-ω, and found that the RNG k-ε model can give an excellent solution compared with the experiments [20].As a result, this study employs the RNG k-ε model to close the turbulence.Thus, the turbulence kinetic energy k equation is given by where G k represents the generation of turbulent kinetic energy due to the mean velocity gradients.G b represents the generation of turbulent kinetic energy due to buoyancy.Turbulent energy dissipation rate ε equation is given by where µ t is the turbulence viscosity, µ t0 is the value of turbulent viscosity without swirl modification, Ω is a characteristic swirl number, and αs is a swirl constant set as 0.07 in the current study.The constants used here are C 1ε = 1.44 and C 2ε = 1.92.

Sliding Mesh Method for Impeller Motion
The simulation of KRs involves the impeller motion, which requires specified numerical tactics.In the current study, the sliding mesh method (SMM) [24][25][26][27][28][29] is adapted, and full transient simulations are carried out using two grid zones (see Figure 1).Compared with other technologies used to numerically deal with the region motion, such as the multiple reference frame (MRF) method, there are no additional source terms inserted into the momentum equations, e.g., centrifugal, Coriolis, and Euler forces, but the mesh will be moving.In SMM, one grid zone (a rotational zone) is attached to the impeller to model the rotating motion during the mechanically stirred process.In contrast, the other grid zone is attached to the vessel wall (a fixed region), as shown in Figure 1.Since the nodes in the rotating region move rigidly in a given rotational zone for SMM, it is imperative that the interface between the zones consists of two overlapping faces, one of which is attached to the rotating region while the other is attached to the fixed region.Specifically, the sliding interface consists of a set of identical surface elements accessible from both sides of the interface.In a single movement step, the mesh in the moving region slides with a predefined velocity across the mesh in the fixed region.After each step, the interface vertices in the moving and static region will be reattached according to the initially computed vertex map list.The interaction between the two regions is modeled by interpolating the information across these interfaces.The SMM algorithm is used to allow for velocity and pressure data transfer from one domain to the other and vice versa.
full transient simulations are carried out using two grid zones (see Figure 1).Compared with other technologies used to numerically deal with the region motion, such as the multiple reference frame (MRF) method, there are no additional source terms inserted into the momentum equations, e.g., centrifugal, Coriolis, and Euler forces, but the mesh will be moving.In SMM, one grid zone (a rotational zone) is attached to the impeller to model the rotating motion during the mechanically stirred process.In contrast, the other grid zone is attached to the vessel wall (a fixed region), as shown in Figure 1.Since the nodes in the rotating region move rigidly in a given rotational zone for SMM, it is imperative that the interface between the zones consists of two overlapping faces, one of which is attached to the rotating region while the other is attached to the fixed region.Specifically, the sliding interface consists of a set of identical surface elements accessible from both sides of the interface.In a single movement step, the mesh in the moving region slides with a predefined velocity across the mesh in the fixed region.After each step, the interface vertices in the moving and static region will be reattached according to the initially computed vertex map list.The interaction between the two regions is modeled by interpolating the information across these interfaces.The SMM algorithm is used to allow for velocity and pressure data transfer from one domain to the other and vice versa.

Tracer Transport Equation
To measure the mixing performance of KRs, the mixing time as a significant indicator has been adapted.Consequently, a passive scalar transport equation of the tracer, which is essentially governed by a convective-diffusion equation, has to be solved at each time step.

𝜌 𝜕(𝛼
where C is the local tracer concentration, αL is the liquid phase volume fraction, and Deff is the effective diffusivity.The effective diffusion coefficient may be calculated as the sum of the molecular Dlam and turbulent diffusivities Dt.Here, Dlam is set as 10 −9 m 2 /s, a typical value for a solute in liquids, and Dt is calculated from the turbulent kinetic viscosity νt as

Tracer Transport Equation
To measure the mixing performance of KRs, the mixing time as a significant indicator has been adapted.Consequently, a passive scalar transport equation of the tracer, which is essentially governed by a convective-diffusion equation, has to be solved at each time step.
where C is the local tracer concentration, α L is the liquid phase volume fraction, and D eff is the effective diffusivity.The effective diffusion coefficient may be calculated as the sum of the molecular D lam and turbulent diffusivities D t .Here, D lam is set as 10 −9 m 2 /s, a typical value for a solute in liquids, and D t is calculated from the turbulent kinetic viscosity ν t as D t = ν t /Sc t , where Sc t is the turbulent Schmidt number and its value 0.7 is used.The velocity field used in the above tracer equation was obtained from a steady solution of VOF simulation and remained as a constant value during the calculation of this passive scalar.
The initial tracer concentrations are designated as zero everywhere in the vessel except at the blue Region A, as shown in Figure 2. To monitor the variation of tracer with time when measuring the mixing time, seven probe points were placed in several typical locations.Their coordinates are also listed in Figure 2. Dt = νt/Sct, where Sct is the turbulent Schimidt number and its value 0.7 is used.The velocity field used in the above tracer equation was obtained from a steady solution of VOF simulation and remained as a constant value during the calculation of this passive scalar.
The initial tracer concentrations are designated as zero everywhere in the vessel except at the blue Region A, as shown in Figure 2. To monitor the variation of tracer with time when measuring the mixing time, seven probe points were placed in several typical locations.Their coordinates are also listed in Figure 2.

Computational Details and Boundary Conditions
A mesh with 450,000 cells, which is a reasonable balance between the computational cost (630,000 cells) and accuracy (260,000 cells), was used after the grid independence test in the previous research, which showed that the change of simulated results was <4% when further refining mesh.The typical mesh is shown in Figure 3. Figures 3a,b show the mesh details on the central vertical sections and the impeller.

Computational Details and Boundary Conditions
A mesh with 450,000 cells, which is a reasonable balance between the computational cost (630,000 cells) and accuracy (260,000 cells), was used after the grid independence test in the previous research, which showed that the change of simulated results was <4% when further refining mesh.The typical mesh is shown in Figure 3. Figure 3a,b show the mesh details on the central vertical sections and the impeller.ity field used in the above tracer equation was obtained from a steady solution of VOF simulation and remained as a constant value during the calculation of this passive scalar.
The initial tracer concentrations are designated as zero everywhere in the vessel except at the blue Region A, as shown in Figure 2. To monitor the variation of tracer with time when measuring the mixing time, seven probe points were placed in several typical locations.Their coordinates are also listed in Figure 2.

Computational Details and Boundary Conditions
A mesh with 450,000 cells, which is a reasonable balance between the computational cost (630,000 cells) and accuracy (260,000 cells), was used after the grid independence test in the previous research, which showed that the change of simulated results was <4% when further refining mesh.The typical mesh is shown in Figure 3. Figures 3a,b show the mesh details on the central vertical sections and the impeller.As for boundary conditions in the current calculation, the top surface of the vessel is treated as a pressure outlet boundary, and the ladle wall, bottom, shaft, and impeller are treated as the no-slip wall boundary condition.The PISO algorithm is used to decouple the pressure-velocity.The discretization of the convective terms in the momentum and turbulence model equations employed the QUICK scheme, and the time step was set to 0.001 s.The calculations were carried out using the commercial computational fluid dynamic software FLUENT 14.5.The main convergence criterion is that the residuals for torque, velocity, and pressure achieve a steady state, i.e., they do not change significantly with further increasing time steps or iterations.The second criterion for checking the residuals was committed for all governing equations to less than 1 × 10 −5 .The detailed geometric and physical parameters in this study are shown in Table 1.

Model Validation
Model validations have been reported [20].Because the availability of the industrial KR process data is limited by the measurement technology and the high costs of in-situ experimentation due to a high temperature under hazardous conditions, the water model experiments were resorted to in order to validate the model.Here, similar conditions of the water model for KR have been adapted.In earlier studies, the effect of the impeller dimension on the fluid flow and interface behavior was investigated based on numerical simulation.It was found that the flow field and interface behavior of the process of KR hot metal pretreatment can be represented by a 3D Volume of Fluid (VOF) model coupling the sliding mesh method and RNG k-ε model.The water model experiment was designed to validate the present mathematical model.The simulated results based on the current model are in good agreement with the measured results in the water model, which corroborates the applicability of the proposed model.
Furthermore, concerning the rotational speed and immersion depth, the comparison of the vortex profile between a water model and numerical simulation at the rotation velocity of 180 rpm and different impeller immersion depths is shown in Figure 4. Figure 5 shows the vortex depth of the water model and the numerical simulation at the impeller immersion depth of 110 mm with different rotation speeds.It is observed that the shape and vortex depth can be well captured, and specific agreement between the simulated and experimental results can reasonably arrive.

Fluid Flow Pattern and Qualifying Inactive Zone
The bath flow in KRs is driven by the high-speed rotating impeller, which directly transfers kinetic energy to the liquid and mixes it with reagent powders for desulfurization.Consequently, it is worth investigating the effect of impeller immersion depth (IID) and impeller rotation speed (IRS) on the stirring and mixing characteristic of KRs. Figure 6 shows the central vertical cross-section streamlines and a close-view of the local region around the impeller at an immersion depth of 100 mm, 110 mm, 120 mm, and 130 mm.The simulation conditions can be found in Table 1, and here the rotation speed of the impeller is fixed to 180 rpm.Similar to our previous study, two main recirculation flow patterns, where the flow field is separated into the upper and lower zones by the discharge flow driven by the impeller rotation, were also observed in Figure 6.With the increase of the impeller immersion depth, the root location of discharge flow moves downward, but the discharge strength, i.e., the maximum velocity, appears to remain invariant for all cases; see the corresponding close-view of the local region around the impeller in Figure 6a-d

Fluid Flow Pattern and Qualifying Inactive Zone
The bath flow in KRs is driven by the high-speed rotating impeller, which directly transfers kinetic energy to the liquid and mixes it with reagent powders for desulfurization.Consequently, it is worth investigating the effect of impeller immersion depth (IID) and impeller rotation speed (IRS) on the stirring and mixing characteristic of KRs. Figure 6 shows the central vertical cross-section streamlines and a close-view of the local region around the impeller at an immersion depth of 100 mm, 110 mm, 120 mm, and 130 mm.The simulation conditions can be found in Table 1, and here the rotation speed of the impeller is fixed to 180 rpm.Similar to our previous study, two main recirculation flow patterns, where the flow field is separated into the upper and lower zones by the discharge flow driven by the impeller rotation, were also observed in Figure 6.With the increase of the impeller immersion depth, the root location of discharge flow moves downward, but the discharge strength, i.e., the maximum velocity, appears to remain invariant for all cases; see the corresponding close-view of the local region around the impeller in Figure 6a-d   In contrast, the flow pattern under different IRSs is also shown in Figure 7a-e, w gives the vertical cross-section flow field and streamlines at IRSs of 180 rpm, 200 rpm rpm, 240 rpm, and 260 rpm.Flow velocity in the bath, both local and global, has increased significantly with the increase of IRS.However, the root location of disch flow with increasing IRS shows almost no change.These results reveal that the effec IRS and IID on the flow pattern are different.These results also make it possible to ationally control the flow pattern in a KR via the complementarity of the two param Further, some quantities comparisons were also conducted.Figures 8(a  Further, some quantities comparisons were also conducted.Figure 8a,b show the axial and radial velocity variations along Lines A and B. Line A is a vertical line of 0.075 m from the axial centerline of the vessel, and Line B is a horizontal line of 0.85 m from the bottom surface, which is illustrated in the subgraph of Figure 8a(1).Through comparison, it can be seen that the changing trend of velocities on Line A is the same.However, the level of maximum velocity decreases with the increase of IID, which can also be observed from the velocity contour of the central vertical section at different IIDs. Figure 8a(2) shows that the radial velocity distribution is almost axis-symmetrical.The velocity near the outside edge of the impeller is greater than at the center and near the wall.It can also be seen in Figure 8a(2) that the velocity along Line B below the impeller slightly increases when increasing IID, and the closer it is to the edge of the impeller, the more pronounced this increase is.Besides, the velocity below the impellers is always minimal, reaching almost zero.Finally, it should be stressed that the velocity change caused by the IID change is finite compared with the influence of IRS.
Figure 8b(1,2) show the axial and radial velocity variations at different IRSs.It can be observed that the axial and radial velocity increases significantly, and the diameter of the lower-velocity region decreases as the impeller rotation speed increases.However, for all the impeller rotation velocities, the velocity below the impeller is still miniscule, making it particularly difficult to achieve global mixing of the species in the reactor.Furthermore, Figure 9 shows the velocity distribution contour of horizontal sections with different heights under the different IRSs.It can be seen that the liquid velocity of different horizontal sections increases noticeably with increasing IRS.However, the velocity below the impeller is still very small no matter how fast the IRS is.These phenomena were also observed by T. Yamamoto et al. [7][8][9][10] in the flow characteristic of aluminum melt by mechanically stirring.miniscule, making it particularly difficult to achieve global mixing of the species in the reactor.Furthermore, Figure 9 shows the velocity distribution contour horizontal sections with different heights under the different IRSs.It can be seen that the liquid velocity of different horizontal sections increases noticeably with increasing IRS.However, the velocity below the impeller is still very small no matter how fast the IRS is.These phenomena were also observed by T. Yamamoto et al. [7][8][9][10] in the flow characteristic of aluminum melt by mechanically stirring.To further study the effect of operation parameters on the bath flow of KRs, the mean velocity of the fluid in the bath under different IIDs was also investigated by indirectly measuring the dynamic energy of the bath obtained from the mechanical stirring.Figure 10a shows the comparison of the mean velocities for different IIDs.As expected, when the immersion depth increases, the mean velocity is basically unchanged.The comparison of the mean velocity of the bath under different IRSs is shown in Figure 10b.The mean ve- To further study the effect of operation parameters on the bath flow of KRs, the mean velocity of the fluid in the bath under different IIDs was also investigated by indirectly measuring the dynamic energy of the bath obtained from the mechanical stirring.Figure 10a shows the comparison of the mean velocities for different IIDs.As expected, when the immersion depth increases, the mean velocity is basically unchanged.The comparison of the mean velocity of the bath under different IRSs is shown in Figure 10b.The mean velocity increases linearly as the rotation speed increases but remains invariant when IID changes.The mean velocity of the liquid rises 0.022 m/s as the impeller rotation velocity increases by 10 rpm.Unlike related research, the relationship between the mean velocity of the bath and IRS has been further obtained.A fitting regression, where the mean velocity v is correlated with the impeller rotation speed ω, is calculated as follows: To further study the effect of operation parameters on the bath flow of KRs, the mean velocity of the fluid in the bath under different IIDs was also investigated by indirectly measuring the dynamic energy of the bath obtained from the mechanical stirring.Figure 10a shows the comparison of the mean velocities for different IIDs.As expected, when the immersion depth increases, the mean velocity is basically unchanged.The comparison of the mean velocity of the bath under different IRSs is shown in Figure 10b.The mean velocity increases linearly as the rotation speed increases but remains invariant when IID changes.The mean velocity of the liquid rises 0.022 m/s as the impeller rotation velocity increases by 10 rpm.Unlike related research, the relationship between the mean velocity of the bath and IRS has been further obtained.A fitting regression, where the mean velocity ̅ is correlated with the impeller rotation speed ω, is calculated as follows:  Additionally, to further examine the effect of IID and IRS on the volume fraction of the inactive zone, an indicator γ is defined, which is the ratio of the region size of flow velocity when it is less than 0.1 m/s to the fluid volume in the bath.The γ at different IIDs is shown in Figure 11.With the increase of IID, the γ decreased from 7.63% to 4.92%.Therefore, it can be inferred that the impeller immersion depth does significantly affect local mixing.A fitting regression of the volume fraction of the inactive zone γ as a function of impeller immersion depth I is obtained as follows.The volume fraction of the inactive zone at different rotation velocities is shown in Figure 12.The γ of the stirred molten pool decreases from 6.49% to 4.32% as the impeller rotation velocity increases from 180 rpm to 260 rpm.The γ decreases by 0.92% as the impeller rotation velocity increases from 180 rpm to 200 rpm, which is the maximum decline The volume fraction of the inactive zone at different rotation velocities is shown in Figure 12.The γ of the stirred molten pool decreases from 6.49% to 4.32% as the impeller rotation velocity increases from 180 rpm to 260 rpm.The γ decreases by 0.92% as the impeller rotation velocity increases from 180 rpm to 200 rpm, which is the maximum decline in the range of 180 rpm to 260 rpm.The γ drops linearly and monotonically in this range, i.e., the γ decreases by 0.2% as the impeller rotation velocity increases by 10 rpm.The volume fraction of the inactive zone at different rotation velocities is shown in Figure 12.The γ of the stirred molten pool decreases from 6.49% to 4.32% as the impeller rotation velocity increases from 180 rpm to 260 rpm.The γ decreases by 0.92% as the impeller rotation velocity increases from 180 rpm to 200 rpm, which is the maximum decline in the range of 180 rpm to 260 rpm.The γ drops linearly and monotonically in this range, i.e., the γ decreases by 0.2% as the impeller rotation velocity increases by 10 rpm.Further, a fitting regression of γ as a function of impeller rotation speed ω is obtained as follows:

Mixing Time
Mixing time is an essential index to measure mixing characteristics.In order to obtain the mixing time of the reactor, the stimulus-response method was used here.When simulations/experiments are fully developed, a certain amount of tracer, specifically, 50 mL of saturated KCl solution, is added to the reactor.The change of tracer concentration is recorded, and the mixing time defined as the time span from the moment the tracer is supplied to that when the concentration curves of all monitoring points approximate a constant value with a fluctuation range restricted to ±5 pct.Some related computational details are given in Section 2.4. Figure 13 plots the specific tracer concentration change at those monitoring points defined in Figure 2. Here, the response curves of tracer with time under different IRSs are shown in Figure 13a-e, which correspond to the IRSs of 180 rpm, 200 rpm, 220 rpm, 240 rpm, and 260 rpm, respectively.In Figure 13, the red vertical line indicates the mixing time based on the above definition.The curve of each point shows how to respond to the change of tracer concentrations and reveals the flow pattern to a certain extent.For example, monitoring point 7, the curves converge gradually on the equilibrium concentration for IRSs of 180, 200 and 220 rpm.Next, the curve has one peak for the IRS at 240 rpm.As the IRS increases to 260 rpm, the curves quickly converge on the equilibrium concentration; moreover, the period shortens with the rise of the IRS.Since the rotating impeller drives the fluid flows, the rotating velocity directly impacts the kinetic energy transferred to the fluid.Thus, with the increase of IRS, the response of the tracer to the fluid hydrodynamic becomes more and more intense, and the tracer concentration at each monitoring point arrives more quickly at the theoretical equilibrium concentration, which suggests the mixing should be reinforced significantly with increasing IRS.Furthermore, the mixing time at different IRSs is shown in Figure 14.It decreases significantly with the increase of IRS.As can be seen from the above, with the increase of IRS, the mean velocity of the bath increases significantly.The minimum mixing time is 55 s for the IRS at 260 rpm-a maximum value in this study's range.Meanwhile, the volume fraction of the inactive zone decreases visibly, which produces a decrease in mixing time.Consequently, the rise of the IRS is an effective means to shorten the mixing time.Furthermore, the mixing time at different IRSs is shown in Figure 14.It decreases significantly with the increase of IRS.As can be seen from the above, with the increase of IRS, the mean velocity of the bath increases significantly.The minimum mixing time is 55 s for the IRS at 260 rpm-a maximum value in this study's range.Meanwhile, the volume fraction of the inactive zone decreases visibly, which produces a decrease in mixing time.Consequently, the rise of the IRS is an effective means to shorten the mixing time.The mixing time under different IIDs is shown in Figure 15.It can be seen that the variation trend of mixing time with IID is first reduced and then remains almost unchanged.The minimum mixing time is 71.5 s at an IID of 130 mm.Practically, when the IID range is between 110 mm and 130 mm, it always takes about 72 s to achieve global mixing.Obviously, IID has a limited impact on mixing time in this IID range.As is indicated above, with the increase of IID, the mean velocity of the bath has almost no change, but the volume fraction of the inactive zone decreases visibly.The volume fraction of the inactive zone at an IID of 100 mm is significantly greater than the other IIDs, making global mixing more difficult and increasing the mixing time.Consequently, compared with IID, IRS has a more significant effect on the mixing performance of the bath in KRs.

Vortex Core Depth
Vortex characteristics have an important influence on the dispersion behavior of particles and the refining efficiency, as reported in some studies [20].It can be seen that there is almost no change to the profile of the gas-liquid interface and vortex core depth at different impeller immersion depths, as shown in Figure 16.The distance between the bottom of the vortex core and the top edge of the impeller increases continuously as IID increases.However, this doesn't mean that the immersion depth of the impeller does not affect the behavior of the gas-liquid interface; this is discussed further in Section 3.4.The mixing time under different IIDs is shown in Figure 15.It can be seen that the variation trend of mixing time with IID is first reduced and then remains almost unchanged.The minimum mixing time is 71.5 s at an IID of 130 mm.Practically, when the IID range is between 110 mm and 130 mm, it always takes about 72 s to achieve global mixing.Obviously, IID has a limited impact on mixing time in this IID range.As is indicated above, with the increase of IID, the mean velocity of the bath has almost no change, but the volume fraction of the inactive zone decreases visibly.The volume fraction of the inactive zone at an IID of 100 mm is significantly greater than the other IIDs, making global mixing more difficult and increasing the mixing time.Consequently, compared with IID, IRS has a more significant effect on the mixing performance of the bath in KRs.The mixing time under different IIDs is shown in Figure 15.It can be seen that the variation trend of mixing time with IID is first reduced and then remains almost unchanged.The minimum mixing time is 71.5 s at an IID of 130 mm.Practically, when the IID range is between 110 mm and 130 mm, it always takes about 72 s to achieve global mixing.Obviously, IID has a limited impact on mixing time in this IID range.As is indicated above, with the increase of IID, the mean velocity of the bath has almost no change, but the volume fraction of the inactive zone decreases visibly.The volume fraction of the inactive zone at an IID of 100 mm is significantly greater than the other IIDs, making global mixing more difficult and increasing the mixing time.Consequently, compared with IID, IRS has a more significant effect on the mixing performance of the bath in KRs.

Vortex Core Depth
Vortex characteristics have an important influence on the dispersion behavior of particles and the refining efficiency, as reported in some studies [20].It can be seen that there is almost no change to the profile of the gas-liquid interface and vortex core depth at different impeller immersion depths, as shown in Figure 16.The distance between the bottom of the vortex core and the top edge of the impeller increases continuously as IID increases.However, this doesn't mean that the immersion depth of the impeller does not affect the behavior of the gas-liquid interface; this is discussed further in Section 3.4.

Vortex Core Depth
Vortex characteristics have an important influence on the dispersion behavior of particles and the refining efficiency, as reported in some studies [20].It can be seen that there is almost no change to the profile of the gas-liquid interface and vortex core depth at different impeller immersion depths, as shown in Figure 16.The distance between the bottom of the vortex core and the top edge of the impeller increases continuously as IID increases.However, this doesn't mean that the immersion depth of the impeller does not affect the behavior of the gas-liquid interface; this is discussed further in Section 3.4.Increasing rotation speed is an effective way to increase the vertex core depth (VCH), favoring the entrainment and dispersion of the reagent by an impeller.The vortex depth and gas-liquid interface profile under different rotation speeds are shown in Figure 17.This figure vividly depicts how impeller rotation speed strongly impacts the vortex core depth and gas-liquid interface profile.The vortex core depth gradually increases as the IRS increases.Until reaching a rotation speed of 260 rpm, the bottom of the vortex passes the upper interface of the impeller, which is advantageous to the dispersion behavior of particles and the refining efficiency.The results of the vortex core depth with respect to the numerical simulation and water model experiment under different rotation speeds were also plotted in Figure 17.T. Yamamoto et al. [7,8] have reported the morphological changes to a free surface with different IRSs.The present work further analyses the relationship between the vortex depth and impeller rotation speed.To further investigate the effect of IID and IRS on the interfacial behavior of KRs, the free surface velocity is reported in Section 3.4.A fitting regression of the vortex core depth HVCD as a function of impeller rotation speed ω is obtained as follows.Increasing rotation speed is an effective way to increase the vertex core depth (VCH), favoring the entrainment and dispersion of the reagent by an impeller.The vortex depth and gas-liquid interface profile under different rotation speeds are shown in Figure 17.This figure vividly depicts how impeller rotation speed strongly impacts the vortex core depth and gas-liquid interface profile.The vortex core depth gradually increases as the IRS increases.Until reaching a rotation speed of 260 rpm, the bottom of the vortex passes the upper interface of the impeller, which is advantageous to the dispersion behavior of particles and the refining efficiency.The results of the vortex core depth with respect to the numerical simulation and water model experiment under different rotation speeds were also plotted in Figure 17.T. Yamamoto et al. [7,8] have reported the morphological changes to a free surface with different IRSs.The present work further analyses the relationship between the vortex depth and impeller rotation speed.To further investigate the effect of IID and IRS on the interfacial behavior of KRs, the free surface velocity is reported in Section 3.4.A fitting regression of the vortex core depth H VCD as a function of impeller rotation speed ω is obtained as follows.Increasing rotation speed is an effective way to increase the vertex core depth (VCH), favoring the entrainment and dispersion of the reagent by an impeller.The vortex depth and gas-liquid interface profile under different rotation speeds are shown in Figure 17.This figure vividly depicts how impeller rotation speed strongly impacts the vortex core depth and gas-liquid interface profile.The vortex core depth gradually increases as the IRS increases.Until reaching a rotation speed of 260 rpm, the bottom of the vortex passes the upper interface of the impeller, which is advantageous to the dispersion behavior of particles and the refining efficiency.The results of the vortex core depth with respect to the numerical simulation and water model experiment under different rotation speeds were also plotted in Figure 17.T. Yamamoto et al. [7,8] have reported the morphological changes to a free surface with different IRSs.The present work further analyses the relationship between the vortex depth and impeller rotation speed.To further investigate the effect of IID and IRS on the interfacial behavior of KRs, the free surface velocity is reported in Section 3.4.A fitting regression of the vortex core depth HVCD as a function of impeller rotation speed ω is obtained as follows.

Free Surface Velocity
The velocity on the gas-liquid interface strongly influences the entrainment and entrapment of the desulfurization particles.The interface velocity at different IIDs is shown in Figure 18.It clearly indicates that the velocity around the vortex core is much bigger than that near the walls, i.e., there is a gradient of velocity along the radial direction.It also suggests that the regents around the vortex core are much more active than those near the walls.Furthermore, with the increase of IID, the distance between the bottom of the vortex core and impeller increases, resulting in the rotational impeller's effect on the interface velocity around the bottom of the vortex core being less pronounced.Consequently, the velocity gradient between the axis-centerline and walls becomes insignificant when the impeller is taken down.For example, the velocity on the gas-liquid interface around the vortex core is obviously greater than those near the walls for immersion depths of 100 mm, i.e., the change of interface velocity along with radial direction is relatively steep for the small IID (see Figure 18a).However, for the case of IID at 130 mm, the change of velocity, along with radial direction, is relatively even (see Figure 18d).As a result, the interface velocity near the vortex core decreases significantly with the increasing IID, which reveals that the entrainment and dispersion effect of the vortex is affected by the arrangement of IID.

Free Surface Velocity
The velocity on the gas-liquid interface strongly influences the entrainment and entrapment of the desulfurization particles.The interface velocity at different IIDs is shown in Figure 18.It clearly indicates that the velocity around the vortex core is much bigger than that near the walls, i.e., there is a gradient of velocity along the radial direction.It also suggests that the regents around the vortex core are much more active than those near the walls.Furthermore, with the increase of IID, the distance between the bottom of the vortex core and impeller increases, resulting in the rotational impeller's effect on the interface velocity around the bottom of the vortex core being less pronounced.Consequently, the velocity gradient between the axis-centerline and walls becomes insignificant when the impeller is taken down.For example, the velocity on the gas-liquid interface around the vortex core is obviously greater than those near the walls for immersion depths of 100 mm, i.e., the change of interface velocity along with radial direction is relatively steep for the small IID (see Figure 18a).However, for the case of IID at 130 mm, the change of velocity, along with radial direction, is relatively even (see Figure 18d).As a result, the interface velocity near the vortex core decreases significantly with the increasing IID, which reveals that the entrainment and dispersion effect of the vortex is affected by the arrangement of IID.In addition, the velocity distribution on the gas-liquid interface at different IRSs was investigated, and the simulation results are shown in Figure 19.The velocity on the gasliquid interface increases significantly with the rise of IRS; this phenomenon can be explained in terms of a high-speed rotating impeller that intensely powers the velocity of the gas-liquid interface.Notably, for the case of IRS at 260 rpm, the bottom of the vortex core sunk below the top edge of the impeller.Thus, it makes the interface velocity near the impeller sharply increase and dramatically promotes the dispersion behavior of the reagents and desulfurization efficiency.In addition, the velocity distribution on the gas-liquid interface at different IRSs was investigated, and the simulation results are shown in Figure 19.The velocity on the gas-liquid interface increases significantly with the rise of IRS; this phenomenon can be explained in terms of a high-speed rotating impeller that intensely powers the velocity of the gas-liquid interface.Notably, for the case of IRS at 260 rpm, the bottom of the vortex core sunk below the top edge of the impeller.Thus, it makes the interface velocity near the impeller sharply increase and dramatically promotes the dispersion behavior of the reagents and desulfurization efficiency.

Free Surface Velocity
The velocity on the gas-liquid interface strongly influences the entrainment and entrapment of the desulfurization particles.The interface velocity at different IIDs is shown in Figure 18.It clearly indicates that the velocity around the vortex core is much bigger than that near the walls, i.e., there is a gradient of velocity along the radial direction.It also suggests that the regents around the vortex core are much more active than those near the walls.Furthermore, with the increase of IID, the distance between the bottom of the vortex core and impeller increases, resulting in the rotational impeller's effect on the interface velocity around the bottom of the vortex core being less pronounced.Consequently, the velocity gradient between the axis-centerline and walls becomes insignificant when the impeller is taken down.For example, the velocity on the gas-liquid interface around the vortex core is obviously greater than those near the walls for immersion depths of 100 mm, i.e., the change of interface velocity along with radial direction is relatively steep for the small IID (see Figure 18a).However, for the case of IID at 130 mm, the change of velocity, along with radial direction, is relatively even (see Figure 18d).As a result, the interface velocity near the vortex core decreases significantly with the increasing IID, which reveals that the entrainment and dispersion effect of the vortex is affected by the arrangement of IID.In addition, the velocity distribution on the gas-liquid interface at different IRSs was investigated, and the simulation results are shown in Figure 19.The velocity on the gasliquid interface increases significantly with the rise of IRS; this phenomenon can be explained in terms of a high-speed rotating impeller that intensely powers the velocity of the gas-liquid interface.Notably, for the case of IRS at 260 rpm, the bottom of the vortex core sunk below the top edge of the impeller.Thus, it makes the interface velocity near the impeller sharply increase and dramatically promotes the dispersion behavior of the reagents and desulfurization efficiency.

•
Impeller immersion depth and rotation speed have different effects on the fluid flow pattern of the bath under the present study range.The root location of the discharge flow moves downward with the impeller immersion depth increasing, but the discharge strength and the mean velocity of the bath show hardly any change.
Comparatively, the increase in the impeller rotation speed significantly improves the mean velocity, but there is little change in the position of discharge flow.Furthermore, increasing impeller immersion depth or rotation speed can effectively reduce the volume fraction of the inactive zone, but it cannot eliminate it.As a result, the correlation equations of γ as a function of ω or I are formulated under the range of 180 ≤ ω ≤ 260 and 100 ≤ I ≤ 130, respectively.

•
Increasing impeller rotation speed is the most direct and effective way to shorten the mixing time.However, the impeller immersion depth has a limited impact on the mixing time by comparison.In the present study, a minimum mixing time of 55 s is achieved at the maximum impeller rotation speed of 260 rpm.

•
The vortex core depth and the velocity at the gas-liquid interface increase significantly with the increasing impeller rotation speed, and a linear fitting regression has been proposed.However, while the impeller immersion depth has little effect on the vortex core depth, it has a visible influence on the velocity distribution of the free surface.
The velocity gradient on the gas-liquid interface between the axis to walls becomes steep with the decreasing impeller immersion depth.

Figure 1 .
Figure 1.The computational domain consisted of the fixed and rotational regions and separated the sliding mesh interface.

Figure 1 .
Figure 1.The computational domain consisted of the fixed and rotational regions and separated the sliding mesh interface.

Figure 2 .
Figure 2. Schematic illustration of tracer addition region and the position of monitoring points with detailed coordinates.

Figure 3 .
Figure 3. Numerical mesh of (a) central vertical profiles and (b) impeller.

Figure 2 .
Figure 2. Schematic illustration of tracer addition region and the position of monitoring points with detailed coordinates.

Figure 2 .
Figure 2. Schematic illustration of tracer addition region and the position of monitoring points with detailed coordinates.

Figure 3 .
Figure 3. Numerical mesh of (a) central vertical profiles and (b) impeller.Figure 3. Numerical mesh of (a) central vertical profiles and (b) impeller.

Figure 3 .
Figure 3. Numerical mesh of (a) central vertical profiles and (b) impeller.Figure 3. Numerical mesh of (a) central vertical profiles and (b) impeller.

Figure 4 .
Figure 4. Comparison of vortex depth of numerical simulation and water model at impeller immersion depths of 100 mm and 120 mm.Figure 4. Comparison of vortex depth of numerical simulation and water model at impeller immersion depths of 100 mm and 120 mm.

Figure 4 .
Figure 4. Comparison of vortex depth of numerical simulation and water model at impeller immersion depths of 100 mm and 120 mm.Figure 4. Comparison of vortex depth of numerical simulation and water model at impeller immersion depths of 100 mm and 120 mm.

Figure 4 .
Figure 4. Comparison of vortex depth of numerical simulation and water model at impeller immersion depths of 100 mm and 120 mm.

Figure 5 .
Figure 5.Comparison of vortex depth of numerical simulation and water model at impeller rotation velocities of 220 rpm and 260 rpm.
. As a result, the upper recirculation zone was expanded.This flow pattern change would influence the local and global mixing characteristics.

Figure 5 .
Figure 5.Comparison of vortex depth of numerical simulation and water model at impeller rotation velocities of 220 rpm and 260 rpm.

19 Figure 6 .
Figure 6.The streamlines of the central vertical cross-section at an immersion depth of (a) 100 mm, (b) 110 mm, (c) 120 mm, and (d) 130 mm.In contrast, the flow pattern under different IRSs is also shown in Figure 7a-e, which gives the vertical cross-section flow field and streamlines at IRSs of 180 rpm, 200 rpm, 220 rpm, 240 rpm, and 260 rpm.Flow velocity in the bath, both local and global, has been increased significantly with the increase of IRS.However, the root location of discharge

Figure 6 .
Figure 6.The streamlines of the central vertical cross-section at an immersion depth of (a) 100 mm, (b) 110 mm, (c) 120 mm, and (d) 130 mm.

Figure 6 .
Figure 6.The streamlines of the central vertical cross-section at an immersion depth of (a) 100 (b) 110 mm, (c) 120 mm, and (d) 130 mm.
), (b) show axial and radial velocity variations along Lines A and B. Line A is a vertical line of m from the axial centerline of the vessel, and Line B is a horizontal line of 0.85 m from bottom surface, which is illustrated in the subgraph of Figure 8 (a1).Through compar

Figure 8 .
Figure 8.Comparison of velocity distributions along with Line A for the axial direction in subgraph (1) and Line B for the radial direction in sub-graph (2) at different (a) IIDs and (b) IRSs.

Figure 8 . 19 Figure 9 .
Figure 8.Comparison of velocity distributions along with Line A for the axial direction in sub-graph (1) and Line B for the radial direction in sub-graph (2) at different (a) IIDs and (b) IRSs.Metals 2021, 11, x FOR PEER REVIEW 11 of 19

Figure 10 .
Figure 10.Comparison of the mean velocities of the fluid in the bath at different (a) impeller immersion depths (IIDs) and (b) IRSs.

Figure 10 .
Figure 10.Comparison of the mean velocities of the fluid in the bath at different (a) impeller immersion depths (IIDs) and (b) IRSs.

Figure 11 .
Figure 11.Comparison of the volume fraction of inactive zones at different IIDs.

17 )Figure 11 .
Figure 11.Comparison of the volume fraction of inactive zones at different IIDs.

Figure 14 .
Figure 14.Comparison of mixing times at different IRSs.

Figure 15 .
Figure 15.Comparison of mixing times at different IIDs.

Figure 14 .
Figure 14.Comparison of mixing times at different IRSs.

Figure 14 .
Figure 14.Comparison of mixing times at different IRSs.

Figure 15 .
Figure 15.Comparison of mixing times at different IIDs.

Figure 15 .
Figure 15.Comparison of mixing times at different IIDs.

Figure 16 .
Figure 16.Vortex depth and gas-liquid interface profile at different IIDs with IRSs of 180 rpm.

) 19 Figure 16 .
Figure 16.Vortex depth and gas-liquid interface profile at different IIDs with IRSs of 180 rpm.

Table 1 .
The geometrical, physical, and operating parameters used in this study.