Research Article  Open Access
Numerical Simulation and Experimental Investigation on InclusionArgonLiquid Steel Phenomenon in RH Refining of HighManganese and HighAluminum Steel
Abstract
In this paper, factors influencing the inclusion removal in highmanganese and highaluminum steel in RH refining process were studied by numerical simulations, production practice, and metallographic experiments. A mathematical model for inclusion removal was established, and the phenomenon of inclusions mixing in RH upleg region was verified due to fluid circulation. Removal efficiency of RH circulation time 120 s is much better than 600 s, and it was the lowest efficiency after 600 s. After 600 s circulation time, it shall not apply in production practice. The mass concentration of inclusions in practical steel was 11.64% probability error than values obtained by numerical simulation, because the numerical simulation did not consider the problem of inclusions adsorbing to the walls of refractory materials and corrosion of refractories. This work lays the foundation for the optimization and upgrading of process technology and establishes big data for full automation of RH out of furnace refining.
1. Introduction
With the optimization of quality and process, significant progress has been made in the applications of highmanganese and highaluminum steel. They possess excellent strength, low density, high quality, and good elongation, thereby showing utility in automobile inner plates, electrical motor rotors, marine engineering, military applications, and refrigeration materials. However, the problem of nonmetallic inclusions in metallurgical industry seriously affects the quality of highmanganese and highaluminum steel, which leads to lower service life. With the rapid development in the steel industry, there is increasing demand for highquality steel. Researchers have found that oxide nonmetallic inclusions play an important role in the performance of clean steel [1–3]. The removal of inclusions and the control of small inclusions are important factors that directly affect the quality of steel. RH refining in steelmaking process has prominent technical advantages for smelting clean steel [4]. It is conducive to the removal of inclusions and the improvement of steel quality. Therefore, removal of inclusions in steel in RH refining process has been the focus of metallurgical research for several years [5–12].
Chen et al. analyzed the process of inclusion collision and growth removal, and considered the influence of gas bubbles in the formation model. However, the influence factors of inclusions adsorbing RH refractory material walls were neglected [13]. The model also did not meet the overall mass conservation of inclusions, and did not consider the distribution of inclusions in RH flow field [14]. Taking ladle as the research object, the inclusion removal in the refining process was studied using the equation model of particle number concentration. The research focused on the formation mechanism of inclusions and the distribution of gas stirring in clean steel. Wang et al. studied the crystal structure and cluster mechanism of inclusions, especially alumina [15]. However, there was no reasonable demonstration on the distribution of inclusions or a mathematical model. Since there were no detailed descriptions and explanation of inclusions, it is difficult to use their experimental conclusions to guide the process production practice [16, 17]. Describing the RH flow field of inclusions and establishing the physical model, the mathematical model did not consider the problem of inclusions growth and collision. It also did not consider the probability error between the model establishment and production practice. Due to the shortcomings of existing models, this study uses actual production data to modify the RH inclusions removal model and analyzes the influence factors of RH circulation flow rate and circulation time on inclusions removal. Predicting the accuracy and reliability of RH inclusion removal provides theoretical guidance for process production practice.
In this paper, a mathematical model of RH inclusionargonliquid steel is established in order to provide the foundation for process production by combining mathematical simulation with production practice. By determining the influence of RH circulation time and circulation flow rate on inclusionargonmolten steel, the RH inclusionsargonmolten steel are analyzed by using the methods of numerical simulation and production practice, and the feasibility of the mathematical model is demonstrated on the basis of experimental data. The proposed method overcomes the probability error between the mathematical model and experimental data, making it more valuable for production practice. Two technical routes of numerical simulation and experimental research are used to analyze and predict the accuracy and reliability of RH inclusion removal, which provide a theoretical basis for process production practice and lay a basis for fully automated steelmaking.
2. Model Building
2.1. Basic Hypothesis
The hypothesis of the model is as follows.
The liquid steel is an incompressible Newtonian fluid, and its multiphase flow is unsteady. The fluctuation of liquid level in ladle top slag and vacuum chamber and the influence of temperature on other physical quantities are neglected. Inclusions nucleate into spheres instantaneously in molten steel, ignoring the effect of slag on other physical quantities. Ladle refractories and RH refractories set wall function. The effect of inclusions on the physical quantity of molten steel is neglected in fluid flow. When inclusions reach the molten steel level, they can be considered to be removed. Argon bubbles are assumed to have a circular [18, 19].
2.2. Control Equation
In the simulation, continuity equations were used to model the reference materials, as follows:where refers to the average volume density, kg/m^{3}; is the average mass velocity, m/s. A homogeneous flow model is used to simulate the fluid field, and the average density of gas and liquid are and , respectively, kg/m^{3}. is the gas fraction, %. RH under vacuum conditions, the gas fraction is calculated as follows:where H refers to the liquid level height of the vacuum chamber, m; Z_{0} corresponds to nozzle height, m; Z refers to the height of the twophase region, m; refers to vacuum chamber pressure, Pa; refers to gravitation acceleration, m/s^{2}; and refers to the formula for determining the distribution of gas fraction hold up, %.
Momentum equation:where refers to turbulent viscosity, Pa·s; refers to mixing viscosity, Pa·s; refers to gas fraction, %; and refers to gas velocity vector, m/s.
Energy equation:
The turbulence model constants are and . refers to the mixing velocity vector, m/s. refers to the turbulent energy dissipation rate, m^{2}/s^{3}.
is expressed as follows:
Turbulent flow energy dissipation equation:
In the formulas (4) and (6), the values of the constants in the above turbulence models are as follows: , , , , . Inclusion collision growth model of density distribution function and radius of inclusions R are related by the expression . The number density of inclusions is . A and B refer to inclusion distribution factor. The mass concentration can be expressed as [2, 11, 20, 21]:
Here, there is no top slag on the liquid surface of the RH vacuum chamber. Inclusions adsorbed by bubbles return to the molten steel. Inclusions adsorbed by bubbles float to the liquid level of the ladle for removal. The liquid level in RH vacuum chamber has no ability to remove inclusions directly. The mass conservation equation of inclusion collision removal can be expressed as follows [22]:where refers to the decrease of the quantity density caused by the collision polymerization of inclusions, m^{−3}/s; refers to effective diameter, m; and refers to inclusion of bulk density, ppm. RH circulating fluid flow is mainly comprises of turbulent collision and Stokes collision.
In the formula (9), is the kinematic viscosity, m^{2}/s; is the constant; μ is the dynamic viscosity, Pa·s; and is the average radius of inclusions, m. Boundary conditions ignore the surface shear stress of RH steel and vacuum chamber. The refractory materials area set wall function, the wall function node is a nonslip boundary condition, and the gradient of wall velocity component is zero. Table 1 lists the inclusionsargonmolten steel physical parameters, where we set the initial RH vacuum to 100 Pa, the gas is argon, and the initial gas circulation flow is 120 Nm^{3}/h, 110 Nm^{3}/h and 100 Nm^{3}/h. The initial mass concentration of inclusions is 187 ppm. The VOF (volume of fluid) model is used for RH fluid circulation, and the DPM (discrete phase model) is used for RH inclusions movement. Figure 1 shows the RH ladle physical model and grid model. SolidWorks is used for physical model, ICEMCFD (the integrated computer engineering and manufacturing code for computational fluid dynamics) is used for grid model, ANSYS FLUENT is imported for numerical simulation calculation, and CFDPOST is used for posttreatment analysis.

(a)
(b)
3. Calculation Results and Analysis
3.1. Analysis of the Impact of Inclusions
As shown in Figure 2, the mass concentration of inclusions is calculated by numerical simulation every 120 s, with the total numerical simulation time being 600 s. The probability error between simulated and experimental mass concentration curve of inclusions is low. The average mass concentration probability error of inclusions as determined by the mathematical model differs from the experimental data by 11.64%. The numerical simulation calculates the mass concentration of inclusions every 120 s and compares the probability error between the corrected mass concentration curves of inclusions. The mass concentration comparison between numerical simulation and experimental data for different circulation times is listed in Table 2. For the mass concentration in the experimental and numerical simulation process, there is error for numerical simulation result. This paper analyzes the root of the probability error from the numerical simulation and process practice. This paper is based on the mathematical model of highmanganese and highaluminum steel, and the mathematical model lays the theoretical foundation for RH fullautomatic refining.

The removal of RH inclusions is analyzed by numerical simulation, the circulation time is set to 1200 s. Figure 3 shows removal of RH inclusions with different circulation flow rates. The removal of inclusions in RH with different circulation time is compared with a circulating flow rate of 120 Nm^{3}/h, where the removal efficiency of inclusion mass concentration is the highest. With the increase of circulation time, the removal efficiency of inclusions is the highest in time zone A (0 s–120 s), time zone B (120 s–240 s), time zone C (240 s–360 s), time zone D (360 s–480 s), time zone E (480 s–600 s), and F (600 s–1200 s) within 600 s of circulation time. The removal efficiency of inclusions in time zone F (600 s–1200 s) decreases and achieves the lowest value. Therefore, it is not recommended to utilize the F (600 s–1200 s) circulation time in production practice. This is further corroborated by the following three points: (a) increased RH production time leads to higher inclusions; (b) the inclusion removal efficiency is lower; and (c) since the numerical simulation does not consider the problem of inclusions adsorbing to the refractory material walls and corrosion of refractories, it exists in the process of practice.
Meanwhile, the residual amount of inclusions is influenced by many factors, among which the turbulence and buoyancy plays a dominant role. Figure 4 shows the relationship between mass concentration and circulation time with inclusions of different sizes during circulation. It can be seen that the inclusion decreases gradually with the increase of time. At the same time, small size inclusion removal efficiency is lower. In comparison, the buoyancy of large inclusions effect is larger, so it is more easy for large inclusions to enter slag. When the circulation time reaches 600 s, buoyancy plays a dominant role in the removal of inclusions, and the diameter of inclusions also has a great impact on the removal efficiency.
RH circulation time affects the efficiency of inclusion removal and the process production. Optimization of RH circulation time is one of the effective methods to remove inclusions. For RH circulation times of 120 s, 360 s, and 600 s, mass concentration of 165 ppm, 136 ppm, and 114 ppm are observed, respectively. Figure 5 shows the inclusions removal at different circulation times and the distribution of inclusions in the fluid domain at different circulation times. From the inclusions in RH particle trajectory in Figure 5, it can be seen that the number of inclusions is the largest when the RH circulation time is 120 s. This is because the circulation time is insufficient for complete removal and only a few inclusions float upwards to the molten steel level. Since there are many inclusions in the steel in the initial 120 s, the probability of inclusions colliding, aggregating, and floating to the liquid steel level is the greatest, so the removal rate of inclusions in the initial 120 s is the fastest. When the RH circulation time reaches 600 s, the amount of inclusions remaining in the steel is very small because the inclusions have enough time to float upward, such that they can be removed. The number of inclusions in molten steel decreases at the later stage of the circulation, which reduces the probability of inclusion collision and aggregation, which significantly reduces the removal rate of inclusions. Compared with the inclusion removal in RH circulation times of 120 s, 360 s, and 600 s, the efficiency of inclusion removal in RH circulation time of 600 s is low. At this time, the distribution of inclusions in steel is the lowest. At the later stage of RH circulation, small size inclusions become a limiting factor for RH refining.
(a)
(b)
(c)
Figure 6 shows the inclusions in ladle particle trajectory from the fluid dynamics analysis. Argon is injected into the upleg, which drives the molten steel due to floatation. The molten steel flows to the vacuum chamber at the upleg and then into the ladle from the downleg due to gravity in the vacuum chamber, thus completing one circulation flow. Due to the interaction of the inclusionsargonmolten steel, the inclusions are adsorbed on the argon bubbles. The fluid flows and the argon bubbles float upward, and the inclusions are brought to the steel slag liquid surface by the argon bubbles. Differential pressure ruptures argon bubbles, and a portion of the inclusions is removed after reaching the steel liquid surface. The rest of the inclusions are rolled and mixed into the molten steel. The trajectory of inclusions in the ladle is shown in Figure 6, where the RH vacuum degree is 100 Pa and circulation flow rate is 120 Nm^{3}/h. As shown in Figure 6, the inclusions are mostly distributed at the ladle liquid level and the upleg area. This is because the inclusions follow the flow of argon and molten steel. Due to floatation, the inclusions from the upleg flow into the vacuum chamber. Under the gravity of molten steel in the vacuum chamber, the inclusions rush into the ladle, while the RH vacuum chamber does not have the ability to remove the inclusions. In addition, due to the large turbulent kinetic energy of the fluid, there is entrainment in the upleg area of RH suction nozzle. The entrainment area of the upleg drives the inclusions on the steel surface to return to the molten steel, which is not conducive to the removal of inclusions. Reasonable control of turbulence intensity is the key to reduce entrainment of inclusions on the steel surface. Reducing the entrainment intensity of the upleg is one of the methods to improve the removal efficiency of inclusions.
3.2. Production Practice Data and Numerical Simulation Analysis
Based on the production of highmanganese and highaluminum steel, RH vacuum is set to 100 Pa and the circulation flow is 120 Nm^{3}/h. Samples are taken every 120 seconds: A0, A120, A240, A360, A480, and A600, until a time of 600 s is achieved. As listed in Table 2, the difference between mass concentration of inclusions measured by numerical simulation and sample electrolytic extraction without aqueous solution is 11.64%. This mass concentration probability error of inclusions occurs because the numerical simulation does not consider the problem of inclusions adsorbing to the refractory material walls and corrosion of refractories. At the later stage of RH circulation, small size inclusions have become the limiting link of RH refining.
Samples of 5 mm × 10 mm are polished by inclusion A0, A120, A240, A360, A480, and A600 and then map scanning Al_{2}O_{3}. The morphology of inclusions is characterized by strips that are irregular, long, and angular, as shown in Figure 7. Most of the samples analyzed by map scanning of inclusions are nonmetallic inclusions, with Al_{2}O_{3} being the main component. The size of A0 inclusion is larger in production practice, whereas the size of A600 inclusion decreases significantly in the later stage of RH circulation. This is because large sized inclusions are removed with the increase of RH circulation time. At the later stage of RH circulation, small size inclusions become the limiting factor of RH refining, which is consistent with the conclusion of numerical simulation.
(a)
(b)
(c)
(d)
(e)
(f)
4. Conclusions
(1)In the production of highmanganese and highaluminum steel by RH, a mathematical model of inclusions removal is established, which lays a theoretical foundation for predicting removal of inclusions from liquid steel.(2)Optimizing RH circulation flow rate and circulation time are effective for removing inclusions by numerical simulations. Compared with the RH circulation times of 120 s and 600 s, the RH circulation time of 120 s has the highest removal efficiency. A 600 s circulation time makes the removal efficiency the lowest, and thus it is not recommended for application in production practice.(3)The turbulence intensity of molten steel in RH fluid field is large, and there exists mixing phenomenon in the upleg area of RH suction nozzle, which leads to “mixing” of inclusions into vacuum chamber or ladle, which affects the removal efficiency of inclusions. In the later stage of RH circulation, removal of small size inclusions becomes a limiting step in RH refining.(4)Samples of highmanganese and highaluminum steel produced by RH are taken for metallographic analysis every 120 seconds after the beginning of circulation time. The size of inclusion A0 is larger in the early stage of RH circulation, and its content decreases in the later stage of circulation.(5)The mass concentration probability error of the mathematical model is 11.64%, because the numerical simulation does not consider the problem of inclusions adsorbing to the refractory material walls and corrosion of refractories. This meets the technological requirements of producing highmanganese and highaluminum steel by RH treatment.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The authors are grateful for the financial support by the National Key R&D Project of China (2017YFC0805100), the Science and Technology Projects of Liaoning Province (2018307003), the Excellent Talent Project of University of Science Technology of Liaoning (2017RC01), and the project (SKLMEAUSTIN201903) supported by State Key Laboratory of Metallic Materials for Offshore Equipment and Applications.
References
 L. C. Trindade, J. J. M. Peixoto, C. A. da Silva et al., “Influence of obstruction at gasinjection nozzles (number and position) in RH degasser process,” Metallurgical and Materials Transactions B, vol. 50, no. 1, pp. 578–584, 2019. View at: Publisher Site  Google Scholar
 W. Liu, S. Yang, J. Li, and H. Yang, “Numerical investigation of inclusion motion at molten steelliquid slag interface during Ruhrstahl Heraeus (RH) process,” JOM, vol. 70, no. 12, pp. 2877–2885, 2018. View at: Publisher Site  Google Scholar
 B. Han, R. Zhu, Y. Zhu et al., “Research on selective oxidation of carbon and aluminum with introduction of CO_{2} in RH refining of lowcarbon steel process,” Metallurgical and Materials Transactions B, vol. 49, no. 6, pp. 3544–3551, 2018. View at: Publisher Site  Google Scholar
 H. Ling and L. Zhang, “Investigation on the fluid flow and decarburization process in the RH process,” Metallurgical and Materials Transactions B, vol. 49, no. 5, pp. 2709–2721, 2018. View at: Publisher Site  Google Scholar
 Y. H. Li, Y. P. Bao, R. Wang, L. F. Ma, and J. S. Liu, “Modeling study on the flow patterns of gasliquid flow for fast decarburization during the RH process,” International Journal of Minerals Metallurgy and Materials, vol. 25, no. 2, pp. 153–163, 2018. View at: Publisher Site  Google Scholar
 G. Chen, S. He, and Y. Li, “Investigation of the airargonsteelslag flow in an industrial RH reactor with VOFDPM coupled model,” Metallurgical and Materials Transactions B, vol. 48, no. 4, pp. 2176–2186, 2017. View at: Publisher Site  Google Scholar
 D.Q. Geng, J.X. Zheng, K. Wang et al., “Simulation on decarburization and inclusion removal process in the RuhrstahlHeraeus (RH) process with ladle bottom blowing,” Metallurgical and Materials Transactions B, vol. 46, no. 3, pp. 1484–1493, 2015. View at: Publisher Site  Google Scholar
 D.Q. Geng, H. Lei, and J.C. He, “Numerical simulation of the multiphase flow in the RheinsahlHeraeus (RH) system,” Metallurgical and Materials Transactions B, vol. 41, no. 1, pp. 234–247, 2010. View at: Publisher Site  Google Scholar
 M. A. van Ende, Y. M. Kim, M. K. Cho, J. Choi, and I. H. Jung, “A kinetic model for the RuhrstahlHeraeus (RH) degassing process,” Metallurgical and Materials Transactions B, vol. 42, no. 3, pp. 477–489, 2010. View at: Publisher Site  Google Scholar
 D. You, S. K. Michelic, C. Bernhard, D. Loder, and G. Wieser, “Modeling of inclusion formation during the solidification of steel,” ISIJ International, vol. 56, no. 10, pp. 1770–1778, 2016. View at: Publisher Site  Google Scholar
 J.H. Wei, N.W. Yu, Y.Y. Fan, S.L. Yang, J.C. Ma, and D.P. Zhu, “Study on flow and mixing characteristics of molten steel in RH and RHKTB refining processes,” Journal of Shanghai University (English Edition), vol. 6, no. 2, pp. 167–175, 2002. View at: Publisher Site  Google Scholar
 C. Zhang, Y. Bao, D. Tang et al., “Effect of size of ferrotitanium on the melting behavior and the yield in the refining of interstitial free Steel by RH degasser,” ISIJ International, vol. 58, no. 12, pp. 2262–2270, 2018. View at: Publisher Site  Google Scholar
 X. Chen, G. Cheng, J. Li, Y. Hou, J. Pan, and Q. Ruan, “Characteristics and formation mechanism of inclusions in 304L stainless steel during the VOD refining process,” Metals, vol. 8, no. 12, pp. 1024–1031, 2018. View at: Publisher Site  Google Scholar
 S.I. Chung, Y.H. Shin, and J.K. Yoon, “Flow characteristics by induction and gas stirring in ASEASKF ladle,” ISIJ International, vol. 32, no. 12, pp. 1287–1296, 1992. View at: Publisher Site  Google Scholar
 G. Wang, Y. Xiao, C. Zhao, J. Li, and D. Shang, “Atomic cluster aggregates in nucleation of solid alumina inclusion in the aluminum deoxidation for liquid iron,” Metallurgical and Materials Transactions B, vol. 49, no. 1, pp. 282–290, 2018. View at: Publisher Site  Google Scholar
 D.Q. Geng, H. Lei, and J.C. He, “Effect of traveling magnetic field on flow, mixing, decarburization and inclusion removal during RH refining process,” ISIJ International, vol. 52, no. 6, pp. 1036–1044, 2012. View at: Publisher Site  Google Scholar
 S.G. Zheng and M.Y. Zhu, “Modelling effect of circulation flow rate on inclusion removal in RH degasser,” Journal of Iron and Steel Research International, vol. 23, no. 12, pp. 1243–1248, 2016. View at: Publisher Site  Google Scholar
 A. Huang, C. Ji, H. Z. Gu, C. H. Dai, and L. Zhong, “Numerical simulation of inclusion removal in tundish with babbling curtain and electromagnetic field,” Journal of Iron and Steel Research, vol. 19, no. 12, pp. 162–166, 2012. View at: Google Scholar
 J. Guo, S.S. Cheng, H.J. Guo, and Y.G. Mei, “Novel mechanism for the modification of Al_{2}O_{3}based inclusions in ultralow carbon Alkilled steel considering the effects of magnesium and calcium,” International Journal of Minerals, Metallurgy, and Materials, vol. 25, no. 3, pp. 280–287, 2018. View at: Publisher Site  Google Scholar
 M. I. Epov, E. P. Shurina, and M. K. Artem’ev, “Numerical homogenization of electric characteristics of media with contrast smallscale inclusions,” Doklady Earth Sciences, vol. 442, no. 1, pp. 97–99, 2012. View at: Publisher Site  Google Scholar
 L.F. Guo, H. Li, H.T. Ling, Y. Wang, and W.C. Song, “Numerical simulation of fractal agglomeratinggrowth of Al_{2}O_{3} inclusions in tundish,” Journal of Iron and Steel Research International, vol. 20, no. 11, pp. 7–12, 2013. View at: Publisher Site  Google Scholar
 R. Tsujino, J. Nakashima, M. Hirai, and I. Sawada, “Numerical analysis of molten steel flow in ladle of RH process,” ISIJ International, vol. 29, no. 7, pp. 589–595, 1989. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Lifeng Chen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.