Effect of Side-blowing Arrangement on Flow Field and Vanadium Extraction Rate in Converter Steelmaking Process

Similar documents
Simulation and Application of Bottom-Blowing in Electrical Arc Furnace Steelmaking Process

Computational Analysis of Blast Furnace Pulverized Coal Injection For Iron Making

CFD modelling of argon stirred ladle. FIMECC Advanced Melt Metallurgy (AMMe) -project

Control of Gas Bubbles and Slag Layer in a Ladle Furnace by Electromagnetic

LECTURE 13. Fundamentals of Converter Steelmaking Technology. Key words: BOF steelmaking, jet penetration, refining reaction, soft and hard jet

Numerical Simulation on Effects of Electromagnetic Force on the Centrifugal Casting Process of High Speed Steel Roll

Lecture 14 Modern trends in BOF steelmaking

Investigating Two Configurations of a Heat Exchanger in an Indirect Heating Integrated Collector Storage Solar Water Heating System

CFD simulation on the turbulent mixing flow performance of the liquid-liquid ejector

NUMERICAL STUDY ON FILM COOLING AND CONVECTIVE HEAT TRANSFER CHARACTERISTICS IN THE CUTBACK REGION OF TURBINE BLADE TRAILING EDGE

Experimental Research on Reducing the Dust of BOF in CO 2 and O 2 Mixed Blowing Steelmaking Process

Metallurgy and lining life in basic oxygen converters

Study of Solid Accretion Formation inside Pyrometallurgical Vessels by a Wax Model and Similarity Conversion for Gas Bottom-Blown Process

NUMERICAL SIMULATION OF SLAG FOAMING IN OXYGEN STEELMAKING

Numerical Simulation of the Flue Gas Desulfurization in the Sieve Tray Tower

Computational Fluid Dynamics (CFD) Simulations of Dispersed Flow of Viscous Oil in Water in a Horizontal Pipe

Study of water falling film over horizontal drop shaped and inverted drop shaped tubes

Investigating two configurations of a heat exchanger in an Indirect Heating Integrated Collector Storage Solar Water Heating System (IHICSSWHS)

EFFECT OF COMBINED HOLE CONFIGURATION ON FILM COOLING WITH AND WITHOUT MIST INJECTION

Numerical Analysis of Two-phase Flow and Heat Transfer in External Cyclic Purification Hot-dip Galvanizing Bath

Numerical analysis of eccentric orifice plate using ANSYS Fluent software

Lecture 12 Converter Steelmaking Practice & combined blowing. Key words: Top blown steelmaking, combined blowing, bottom stirring, hot metal refining

Modelling of the off-gas exit temperature and slag foam depth. of an Electric Arc Furnace

Influence of TiC on the Viscosity of CaO MgO Al 2 O 3 SiO 2 TiC Suspension System

ADVANCES in NATURAL and APPLIED SCIENCES

Modeling of Fluid Flow, Mixing and Inclusion Motion in Bottom Gas- Stirred Molten Steel Vessel

Dr. J. Wolters. FZJ-ZAT-379 January Forschungszentrum Jülich GmbH, FZJ

Simulazione. Numerical Simulation of Magnetic Field and Flow Control in Mold with Vertical Combination Electromagnetic Brake

NUMERICAL SIMULATION OF AIR NATURAL CIRCULATION AND THERMAL RADIATION IN PASSIVE CONTAINMENT COOLING SYSTEM

Analysis of Hot Metal Desiliconization Behavior in Converter Experiments by Coupled Reaction Model

Physical modelling of slag-metal dispersion

Development of Dynamic Models for Oxygen Steelmaking

CFD-study of a 230 MWe coal fired boiler to predict the influence of secondary fuels on slagging, fouling, CO corrosion and NOx formation

JJMIE Jordan Journal of Mechanical and Industrial Engineering

Numerical modelling of shielding gas flow and heat transfer in laser welding process

OPERATING PARAMETER EFFECTS ON THE SOLID CIRCULATION RATE IN A DUAL FLUIDIZED-BED GASIFICATION SYSTEM

COOLING EFFECT ENHANCEMENT IN MAGNETRON SPUTTERING SYSTEM

Available online at ScienceDirect. Procedia Engineering 105 (2015 )

Modeling of Electric Arc Furnaces (EAF) with electromagnetic stirring

IMPROVEMENT OF AIR BLAST SUPPLY OF UPRIGHT CONVERTER

Analysis of bloated droplet theory using steelmaking process models

International Journal of Scientific and Research Publications, Volume 8, Issue 8, August ISSN

Thermal Management of Densely-packed EV Battery Set

Numerical Modeling of Buoyancy-driven Natural Ventilation in a Simple Three Storey Atrium Building

The effect of foamy slag in electric arc furnaces on electric energy consumption

Haiqi YU and Miaoyong ZHU

Effects of electrode immersion depth and remelting rate on electroslag remelting process

Influence of Hot Metal Flow State to the Hearth Flow Field during Blast Furnace Tapping

FLUID STRUCTURE INTERACTION MODELLING OF WIND TURBINE BLADES BASED ON COMPUTATIONAL FLUID DYNAMICS AND FINITE ELEMENT METHOD

Numerical simulation of multiphase flow in a Vanyukov furnace

Tananop Piemsinlapakunchon, and Manosh C. Paul

A study on the effect of the airflow rate of the ceiling type air-conditioner on the ventilation performance

Recent Advances in the Fluid Dynamics of Ladle Metallurgy

Analysis of Droplet Generation in Oxygen Steelmaking

Lecture 25: Principles of degassing

Heat transfer enhancement in fire tube boiler using hellically ribbed tubes

Research Article Influence of the Slag Density on the Splashing Process in a Steelmaking Converter

Modelling of plasma heating of the substrate

Simulation of Atmospheric Air Micro Plasma Jet for Biomedical Applications

Theory Comparison between Propane and Methane Combustion inside the Furnace

Effect of Separation Angle and Nozzle Radial Position on Mixing Time in Ladles with Two Nozzles

CCC Annual Report Quantitative Measurement of Molten Steel Surface Velocity with Nail Boards and SVC Sensor in CC Molds

ASSESSMENT OF BATH MOVEMENT INTENSITY ON PHYSICAL MODEL LADISLAV KOVÁŘ

Reduction of FeO in Molten Slags by Solid Carbon in. the Electric Arc Furnace Operation*

Project Background (long term goals)

Numerical Investigation of the Flow Dynamics of a Supersonic Fluid Ejector

PROBLEM 2: PREDICTING GLASS RIBBON SHAPE IN

Design optimization of a bioreactor for ethanol production using CFD simulation and genetic algorithms

APPLICATION OF THE COMBINED REACTORS METHOD FOR ANALYSIS OF STEELMAKING PROCESS

GAS MIXING AND CHEMICAL HOMOGENIZATION OF STEEL IN 100 T LADLE FURNACE

Numerical modelling of the filling of formworks with self-compacting concrete

Energetica e Macchine, University of Udine, Udine, Italy 3 Dipartimento di Ingegneria Chimica, chimica industriale e scienza dei materiali,

CFD simulations of RTD of a strawberry pulp in a continuous ohmic heater

CFD Modelling and Analysis of Different Plate Heat Exchangers

Mathematical Model of Scrap Melting in EAF using Electromagnetic Stirring

Part III: Slag Practices and Oxygen/Carbon Injection when Melting Direct Reduced Iron

Lecture 23: Injection ladle metallurgy

Study on the permeability index of COREX melter-gisifer and its influencing factors

Numerical analysis of CO concentration distribution on extraction ventilation type tunnel Tao Guo 1,a, Jun Zhou 2,b, Xiaonan Liu 3,c and Bo Shang 4

CFD investigation of coal gasification: Effect of particle size

Lecture 17 Alternative Charge Materials in EAF

PYREJET A MULTI-FUNCTION COMBUSTION/INJECTION SYSTEM FOR EAF STEELMAKING

DRUG DELIVERY PARTICLES FABRICATED USING THE ELECTROHYDRODYNAMIC ATOMIZATION (EHDA) METHOD

A Novel Bottom Stirring Scheme to Improve BOF Performance through Mixing and Mass Transfer Modelling

Simulation of Binary Mixture Droplet Evaporation using VOF Methods

Research on the cavitation characteristic of Kaplan turbine under sediment flow condition

MODELLING COMBUSTION AND THERMAL NO X FORMATION IN ELECTRIC ARC FURNACES FOR THE PRODUCTION OF FERRO-SILICON AND SILICON-METAL

A simulation of the coupled problem of magnetohydrodynamics and a free surface for liquid metals

Metallurgy and Materials Engineering Department, Golpayegan University of Technology, Golpayegan , Iran

DEVELOPMENT OF ULTRA-CLEAN AND LOW-COST PROCESS FOR GALVANIZED STEEL BY LADLE FURNACE CHI-WEI TAI*, CHUNG-YI LI**, TONG-FA YEH***

Characterization of Metal Droplets in Slag after Desulfurization of Hot Metal

Study on reduction behavior of molybdenum trioxide in molten steel

Effects of Reducing Time on Metallization Degree of Carbothermic Reduction of Tall Pellets Bed

A Numerical Investigation on bubble formation in Coolant channel with various Ethylene Glycol and water compositions and change in contact angle

ADAPTING THE FORCHHEIMER EQUATION FOR THE FLOW OF AIR THROUGH WHEAT STRAW BEDS

Comparative Analysis of Different Orifice Geometries for Pressure Drop

Computer modelling of a convective steam superheater

Numerical simulation of vitrification processes: Modeling of air bubbling in molten glass

SIMULATION OF SOLIDIFICATION PROCESS USED IN TWIN ROLL CASTING

Analysis of Fuel Injection System in Blast Furnace with the Help of CFD Software Approach

Transcription:

ISIJ International, Vol. 58 (2018), ISIJ International, No. 5 Vol. 58 (2018), No. 5, pp. 852 859 Effect of Side-blowing Arrangement on Flow Field and Vanadium Extraction Rate in Converter Steelmaking Process Fuhai LIU, 1,2) Dongbai SUN, 1) Rong ZHU, 2) * Kai DONG 2) and Ruiguo BAI 3) 1) National Center for Materials Service Safety, University of Science and Technology Beijing, Beijing, 100083 China. 2) Beijing Key Laboratory of Research Center of Special Melting and Preparation of High-end Metal Materials, University of Science and Technology Beijing, Beijing, 100083 China. 3) Chengde Iron and Steel Company, Chengde, Hebei, 067002 China. (Received on August 3, 2017; accepted on January 11, 2018; J-STAGE Advance published date: March 22, 2018) To improve the vanadium extraction rate, the iron ore powder was used to be injected into the molten bath by side-blowing. Both simulation numerical and water experiment model for analyzing the behaviors of molten bath flow field had been developed to investigate the mixing time, velocity profile and dead-zone volume of molten steel under various injection arrangements and gas flow rates. The result showed that the mixing time would reduce, with increasing the distance between side-blowing nozzle and molten bathe surface. Moreover, although the kinetic energy would be removed by the mutual influenced between bottom-blow bubbles and the side-blowing jet at a certain content, the stirring effect of molten bath is still improved by the side-blowing jet. Based on the result of simulation, an injection arrangement was used in the 150 t vanadium extraction converter. Based on the results, the vanadium content of semi-steel and T. Fe in the slag with side-blowing arrangement is 0.033% and 34.1%, respectively. And the vanadium content of semi-steel and T. Fe in the slag with side-blowing arrangement is 0.044% and 32.1%, respectively. KEY WORDS: side-blowing; powder injection; vanadium extraction; simulation. 1. Introduction As a valued element, the vanadium could improve the plasticity and corrosion resistant of steel. For the traditional vanadium extraction method, the vanadium would be oxidized to V 2 O 5 and transmitted into liquid slag in converter steelmaking process. 1 3) However, there would be an optional oxidation process between vanadium and carbon. Therefore, both the thermodynamic and dynamic conditions should be analyzed to suppress the carbon oxygen reaction for increasing the vanadium extraction rate. Based on the literatures, 4 7) with an appropriate FeO content in liquid slag, lowering the smelting temperature and improving the stirring effect could achieve bigger vanadium extraction rate. At present, optimizing the stirring effect of molten bath is a proper method to enhance the vanadium extraction rate, and the smelting temperature is controlled by adding the slag material during steelmaking process. Both the oxygen lance structure and bottom-blowing arrangement have been researched to improve the dynamic condition of molten bath for a long time, and the powder injection technology is also carried out in recent years. Kruskopf 8) researched the effects of bottom-blowing structure nozzle and flow rate on the flow field of molten bath by * Corresponding author: E-mail: zhurong1206@126.com DOI: http://dx.doi.org/10.2355/isijinternational.isijint-2017-463 numerical simulation with the scrap melting and chemical reaction model. Higuchi and Tago 9) proposed that twisted nozzle oxygen lance could suppress the spitting rate of liquid slag, and increase the stirring effect of molten bath in steelmaking process. Tang and Wang et al. 10) adopted the water experiment to analyze the powder injection effect in converter steelmaking process, and shown the optimum operation was lance height being 1 550 mm, powder size being 1 5 mm and bottom blowing rate being 450 Nm 3 h 1. Yamanoglu and Guthrie et al. 11) found that the diffusion effect of powder would be improved with increasing carrier gas flow rate, and the appropriate solid-gas ratio could enhance the impaction area and stirring effect. When the slag material with the size of powder is injected into molten steel by the side-blowing, the reactions in molten bath could be accelerated, compared to conventional injection method. At the same time, the powder with large kinetic energy could also improve the stirring effect of molten bath, during powder injection process. As a result, both metallurgical and operation benefits could be achieved by the powder injection technology, as shown in the literatures. Hence, Chengde Iron and Steel Company wants to apply the powder injection technology in a 150 t vanadium extraction converter. At present, many research achievements have been obtained mainly on the studies of lime-powder injection method, but the lime is banned in the vanadium extraction process. Since, CaO addition would decrease the level 2018 ISIJ 852

of vanadium content in slag, leading to disadvantage for economic recovery. Hence, it was necessary to carry out the research of water experiment and numerical simulation for the powder injection method before industrial application. In the present study, behavior of flow field and stirring ability of molten bath with various kinds of powder injection arrangement had been analyzed, and the results would be compared with conventional vanadium extraction process. Mixing time, average velocity and dead-zone volume were measured by water experiment and numerical simulation. Based on the results, an appropriate arrangement was applied in a 150 t vanadium extraction converter to analyze smelting time, V 2 O 5 and T. Fe content in slag. 2. Water Experiment Model 2.1. Experimental Principle Based on the similarity of kinetics and geometry, the water experiment model of converter with one powder injection nozzle has been built. The effect of the temperature on the dynamic pressure of the jet is relatively small, according to the research. 12 14) Therefore, the temperature effects were not considered in this paper. As shown in Eq. (1), the modified Froude number was adopted to make the model dynamically similar to the 150 t vanadium extraction converter, and the inertial force and gravity was took into account. 2 u g u1 2 g1 Fr Fr1,that was... (1) gd l g gd1 l1 g1 The flow rate of the water model could be calculated by the Eq. (2): Q Q d l g d l g 1 1 5 1 1 g1 5 l g g1 M g l1 g1 g... (2) where, M was geometric similarity ratio of the model; u, u 1 were gas velocities (m s 1 ); Q, Q 1 presented gas flows (Nm 3 h 1 ); d, d 1 standed for nozzle diameters (mm); ρ l, ρ l1 were liquid density (kg m 3 ); ρ g, ρ g1 presented gas density (kg m 3 ); g acted as acceleration of gravity (m s 2 ); the subscripts l, g were short for liquid and gas. The nomenclature with subscript (1) represent the parameters of model, and the nomenclature without subscript represent the parameters of prototype. The prototype and model parameters were shown in Table 1. 2.2. Experiment Instrument The Fig. 1 showed the main experiment instruments, and the converter was made from plexiglass with the geometric ratio 1:4 in this research. The supersonic oxygen lance nozzle and powder injection nozzle were fabricated from brass material In the experiment, tap water, corn oil, and compressed air were used to represent molten metal, liquid slag and blowing gas, respectively. Eight porous plugs were adopted to simulate the bottom-blowing process. The bottom-blowing arrangement of the model was the same as that of the prototype. Before the water experiment, both water and oil would be added into the plexiglass model according to the Table 1, and then the compressed air would be blown into. After the compressed air passed to run 5 minutes, the Table 1. Geometrical parameters of model and prototype. Laval nozzle Parameters Prototype Model Mach number 2.01 2.01 Diameter of nozzle s throat (mm) 32.6 8.2 Diameter of nozzle s exit (mm) 42.0 10.5 Diameter of powder injection nozzle(mm) 25 6.25 Molten bath depth (mm) 1 478 369.5 Slag layer thickness (mm) 100 25 Furnace diameter (mm) 5 444 1 361 Bottom diameter (mm) 5 008 1 252 Total height (mm) 7 790 1 947.5 Density of molten steel (kg/m 3 ) 7 200 1 000 Density of top-blowing gas (kg/m 3 ) 1.43 1.29 Density of side-blowing gas (kg/m 3 ) 1.36 1.29 Density of bottom-blowing gas (kg/m 3 ) 1.36 1.29 Fig. 1. Experimental instruments of water experiment. tracer element concentration and mixing time would be measured, and other details of experiment procedure could be found in the literature 11. For measuring the tracer element concentration in the bath, three conductivity electrodes (A, B and C) had been installed at different locations, as shown in Fig. 1 During each water experiment, fifty milliliters KCl solution was used as tracer to change the initial conductivity of bath, and the injection location of KCl was on the front end of the funnel, which delivered compressible air into the powder injection nozzle. In this paper, the mixing time was defined as that when the tracer concentration at the monitor points reached 99% of the average tracer concentration in the bath. The T A, T B, and T C presented the mixing times measured by three conductivity electrodes, respectively. The mean mixing time (T M ) was achieved as an arithmetical mean value by Eq. (3): TA TB TC TM... (3) 3 As depicted in the Fig. 2, the injection nozzle were installed 150 mm or 300 mm below the surface of molten bath (red line), based on the angle between nozzle and midline being 0, 22.5, or 45. That means there were six kinds of powder injection arrangement. In this paper, the injection arrangement would be described as A (injection angle)-b (injection depth). For example, the 22.5 150 arrangement means the injection nozzle would be installed 150 mm below the surface of molten bath with the angle being 22.5. 853 2018 ISIJ

domain and sum of all volume fractions must be one in every cell. Moreover, both effective density ρ (kg m 3 ) and effective viscosity μ (Pa s) were solved by the weighted averaging method, according to the volume fraction of each phase. The continuity, momentum, and energy equations would be expressed as: t ui x 0... (4) Fig. 2. The powder injection arrangement in the water model. Table 2. Gas flow rates and lance heights of prototype and model. Side-blowing rate (Nm 3 h 1 ) The oxygen lance height, top-blowing and bottom-blowing flow rate were constants in all experiments. The gas flow rates and lance heights of prototype and model were listed in Table 2. 3. Numerical Model Top-blowing rate (Nm 3 h 1 ) Lance height (mm) Bottomblowing rate (Nm 3 h 1 ) Prototype 80/100/120/140/160 18 000 1 500 300 Model 0.98/1.23/1.47/1.72/1.96 220.71 375 3.68 To analyze the effects of different powder injection arrangements on flow field characteristics of molten bath, a 3-D model was developed with multi-phase flow distribution. In numerical simulation process, the oxygen lance height, side-blowing top-blowing and bottom-blowing rate was 1 500 mm, 120 Nm 3 /h, 18 000 Nm 3 /h and 300 Nm 3 /h, respectively, and the assumptions about numerical model were presented as follows: (1) All chemical reactions or powder injection process were not taken into consideration; (2) The gas phase was regarded as compressible Newtonian fluid, while the liquid slag and molten steel were addressed as incompressible Newtonian fluid; (3) A no-slip condition was applied to all of the walls, and the standard wall function was used to solve the averaged velocity near the wall. 3.1. Governing Equations In the numerical simulation process, the VOF model had been applied to calculate the variation of free interface sharp of multi phases. The fluid phase was defined as a single continuum (α), which did not interpenetrate with other phases. Therefore, the average value of α phase in a cell presented the volume fraction of the cell occupied by the fluid. The i th fluid s volume fraction in the cell was denoted as α i, using to control each phase volume within the domain. 15,16) As a result, three conditions could be defined as follows: α i =0 The cell was empty (of the i th fluid). α i =1 The cell was full (of the i th fluid). 0<α i <1 The cell contained the interface between the i th fluid and one or more other fluids. In the VOF model, the fluids phases in the computational u i uu i j p ij gi fa... (5) t x x x j 2 ct P uct j P p ui T uj ij keff 2... (6) x j x j x j x j For the transport equation, the average effective density (ρ) effective viscosity (μ) and specific heat (c p ) were determined as: gg ll mm... (7) i gg ll mm... (8) cp gcpg, lcp, l mcpm,... (9) where, u i and u j were the velocity components in the direction of i and j, respectively (m s 1 ); t was the physical time (s); p was the pressure (Pa); g was the gravity acceleration (m s 1 ); T was the temperature (C); c p was the specific heat (J kg 1 K 1 ); τ ij and k eff were the viscous stress tensor (N m 2 ), and effective thermal conductivity (W m 1 K 1 ), respectively; The subscripts g, l, and m represented the gasphase, liquid slag and molten steel, respectively. Based on the VOF transport equation, 17,18) the standard k-ε model was adopted to describe the turbulence kinetic energy (k) and its dissipation rate (ε), which could be calculated by following equation: k kui xi... (10) t k Gk Gb YM x j x j ui xi 2 t C1 Gk C3Gb C 2... (11) xj xj In the Eqs. (10) and (11), the G k and G b were the generation of turbulence kinetic energy (kg m 1 s 3 ), due to the mean velocity gradient and the buoyancy, respectively. Y M represented the contribution of the fluctuating dilatation (kg m 1 s 3 ) in compressible turbulence to the overall dissipation rate. The σ κ and σ ε were the turbulent Prandtl numbers for k and ε, respectively. The turbulent viscosity (μ t ) was achieved by combining k and ε as follows: 2 j t C ;... (12) k 2018 ISIJ 854

where C 1ε, C 2ε, C μ, σ κ and σ ε were constants and had the following default values: 18) C 1ε =1.44, C 2ε =1.92. C μ =0.09, σ κ =1.0, σ ε =1.3. 3.2. Computation Methodology The geometry ratio of converter model was 1:1 in the numerical simulation, due to the 150 t vanadium extract converter. The calculation domain of the model included five nozzles oxygen lance, molten bath, side-blowing and bottom-blowing nozzles, and the space from oxygen nozzle tip to molten steel surface. The side-blowing nozzle had the same parameters for all side-blowing arrangements. The thermo-physical properties and specific parameters of multiphase model were depicted in Tables 1 and 3, respectively. The inlet boundary inlet of top-blowing, side-blowing and bottom-blowing were used mass flow inlet, and pressureoutlet boundary condition was adopted for the outlet position of converter (Gauge pressure =101 325 Pa). Initially, the simulation model was calculated with no oxygen or nitrogen blowing through oxygen lance, side-blowing or bottom-blowing nozzles. Because of the large gradient of velocity at the Laval nozzle, the initial time step size was 10 5 s. And then an adaptive method was chosen for the time steps with the global Courant number being 1. During numerical simulation process, the pressure velocity was coupled in PISO scheme with the unsteady solution mode. The Body Force Weighted discretization and Geometric reconstruction scheme were used for calculating the pressure based on the gravity of molten bath and the interface interpolation method, respectively. The other equations were calculated by second-order upwind schemes. Moreover, for each numerical simulation case, the residuals of dependent variables would be addressed as convergent with the residual errors of energy and other variables being less than 10 6 and 10 3, respectively. 3.3. Mesh Sensitivity Test To analyze the mesh sensitivity of the simulation model, the mesh of vanadium converter was addressed as following levels: fine mesh (180 342 cells), medium mesh (143 171 cells), and coarse mesh (100 547 cells). The average flow velocity profile of molten bath for three kinds of mesh levels was depicted in Fig. 3. In order to achieve an accurate result, there were 100 monitor points had been selected randomly to measure and calculate the average flow velocity of molten bath, during numerical simulation process. When the variation of flow velocity was relatively small or fluctuated within a stable section, the quasi-steady-state flow field was developed. The average percentage of variation of the average flow velocity calculated by the coarse and medium mesh level was 9.2%. For the medium and fine mesh levels, the variation was negligible (less than 1%). However, the computational time required for the fine mesh was approximately 1.7 times that for the medium mesh. Therefore, the result calculated by the medium mesh was chosen to research and discuss in this paper. 4. Results and Discussions 4.1. Analysis of Water Experiment The mixing time could represent the stirring effect of molten bath, and the stirring effect improves with the mixing time reducing. The relation between the mixing time and the side-blowing rate has been shown in the Fig. 4, under the same lance height, top-blowing and bottom-blowing rate Fig. 3. Flow velocity of molten bath for three mesh levels. Fig. 4. The mixing time of water experiment with various sideblowing rate and arrangement. Table 3. Thermo-physical properties of three-phase. Molten steel Liquid slag Oxygen Nitrogen Density/(kg m 3 ) 7 000 3 000 Ideal gas 1.36 10 5 Viscosity/(kg m 1 s 1 ) 0.0065 0.35 1.919 10 5 1.663 10 5 Thermal conductivity/(w m 1 K 1 ) 15 1.2 0.0246 0.0242 C P/(J kg 1 K 1 ) 670 1 200 Piecewise-polynomial 1 040.67 Temperature/(K) 1 470 1 470 298 298 855 2018 ISIJ

condition. The mixing time of six kinds of bottom-blowing arrangements are different. Under the combined-blowing condition, the mixing time decreases with rising side-blowing rate. Compared with the average mixing time with different arrangements, the 45 300 arrangement obtains the shortest average mixing time being 38.8 s, and the 0 150 arrangement achieves the longest average mixing time being 47.4 s, at the tested conditions. When the side-blowing flow rate is range from 0.98 to 1.96 Nm 3 /h, the 45 300 arrangement also has the shortest mixing time, which proves the stirring effect of molten bath would be the best with 45 300 arrangement. To research the relationship between the mixing time and the side-blowing rate, the average mixing time of all side-blowing arrangement has been analyzed. Based on the result, the average slopes of the linear regression are 15.2, 9.5, 6.1 and 2.1, with flow rate in the range of 0.98 to 1.23 Nm 3 /h, 1.23 to 1.47 Nm 3 /h, 1.47 to 1.72 Nm 3 /h and 1.72 to 1.96 Nm 3 /h, respectively. That means the downtrend of mixing time would be suppressed with flow rate increasing, at the tested condition. Figure 5(a) depicts the average mixing time of different side-blowing angles. The average mixing time decreases with increasing the injection angel. Considering the sideblowing arrangement and changing rule of mixing time shown in the Figs. 2(a) and 5(a), it seems that kinetic energy of side-blowing jet would be suppressed by the bottomblowing bubbles, making the stirring ability of side-blowing jet can not be fully utilized. As a result, when the distance between the side-blowing jet and bottom-blow bubbles reduces, there is more side-blowing kinetic energy would be removed, which reduces flow velocity and extends mixing time of molten steel. In order to take full advantage of side-blow stirring ability, the side-blowing nozzle should be inserted far from bottom-blowing nozzle as much as possible. Hence, stirring ability of side-blowing jet could be maintained to a certain extent before reaching the bottomblowing bubbles. Figure 5(b) presents the average mixing time of 150 and 300 arrangements. Compared with 300 arrangement, the side-blowing nozzle is more closed to the surface of molten steel with 150 arrangement, according to Fig. 2(b). It seems that the molten bath can be stirred more intensively, when the distance between side-blowing nozzle and bath surface increases. Hence, the mixing time reduces with the 300 arrangement. The shorter mixing time represents the faster reaction rate of molten bath would be improved in an vanadium extraction converter. Comparing with the average mixing time of various side-blowing arrangements being 43.0 s, 45 300 arrangement obtains the shortest average mixing time being 38.8 s, a decreasing ratio of 10.4%. Hence, the 45 300 arrangement is helpful to stir molten bath in steelmaking process, at the tested condition. 4.2. Analysis of Numerical Simulation Figure 6 shows the flow velocity distributions on the vanadium converter cross section at 54 s with different sideblowing arrangements, and the section is 225 mm below the molten bath surface. Different colors represent different flow velocities, with maximum velocity being addressed as red and minimum velocity being blue. The detail values are shown in the figure. As shown in Fig. 6, the flow velocity near the converter wall is quite low with none side-blowing arrangement. Based on the result, the flow velocity, which region is near the side-blowing nozzle, would increase obviously As the arrangement changing, this flow velocity pattern is moreor-less preserved. The average flow velocity of molten bath with 0 150, 0 300, 22.5 150, 22.5 300, 45 150 and 45 300 arrangement is 0.071 m/s, 0.073 m/s, 0.074 m/s 0.076 m/s, and 0.080 m/s, respectively. Compared with the six kinds of side-blowing arrangement, the average flow velocity of molten bath with the 45 300 arrangement is the biggest, so is the area surrounded by the 0.03 m/s isovelocity. The bigger average flow velocity is, the larger area is surrounded by 0.03 m/s isovelocity, and the smaller area of the 0.01 m/s isovelocity is. At the tested condition, when the injection angle of side-blowing nozzle is 0, 22.5 and 45, the average flow velocity is 0.072 m/s, 0.075 m/s and 0.079 m/s, respectively. Moreover, the average flow velocity is 0.074 and 0.076 with the nozzle depth being 150 mm and 300 m, respectively. Based on the results of water experiment and numerical simulation, it seems that the flow velocity of molten bath would improve when the side-blowing nozzle is inserted far from bottom-blowing nozzle and bath surface. Fig. 5. Comparison of mixing time. (a) Average mixing time of different side-blowing angles. (b) Average mixing time of different side-blowing depths. 2018 ISIJ 856

Fig. 6. Flow velocity distribution on cross section at various molten bath depth. Fig. 7. Flow velocity distribution on longitudinal section at different injection angle. As analyzed in water experiment, the velocity interference between side-blowing jet and bottom-blowing bubbles would suppress the stirring ability of side-blowing. Although velocity interference is hard to detect without appropriate tracer material in water experiment, the flow field of molten steel calculated by the numerical simulation could easily show this phenomena. The Fig. 7 present the flow velocity distributions on longitudinal section at different injection angle in the molten bath. As shown in the Fig. 7, there are two main low-velocity regions (the flow velocity is less than 0.01 m/s) in molten bath, which is addressed as region B and region D. The region B is the outcome of combined action of top-blowing and bottom-blowing, and the side-blowing seems show a little effect on this region. It is hard to analyze the effect of side-blow injection angle on the volume of region D, because the bottom-blowing is also a great influence factor on the molten bath. Besides, there is low-velocity region has been formed in the molten bath with the 0 150 arrangement, which is addressed as region C. The Fig. 8 presents the stream of molten flow with the 0 150 arrangement. Based Fig. 8. Flow field of molten bath in region C with 0 150 arrangement. on the result, the movement direction of molten flow, which is formed by the bottom-blowing and side-blowing, is in the opposite direction. Hence, the kinetic energy of molten bath has been removed, and the region C would be formed. As a result, the average velocity of region A with 0 side-blowing injection angle seems lower than other injection angle in the Fig. 6. For the 45 150 arrangement, the velocity of the region D would be improved, since the side-blowing nozzle is inserted far away from the bottom of converter. Therefore, it would take more time to let the bubbles of side-blowing rise to the bath surface, and the interaction of molten bath and side-blowing bubbles should be strengthened. As a 857 2018 ISIJ

result, the stirring effect of molten bath is improved, and the volume of the region D is reduced. In this paper, the zone, where flow velocity of molten steel is lower than 0.001 m/s, is defined as the velocity deadzone, and the average flow velocity and volume of deadzone has been shown in the Fig. 9. Based on the results, the volume of dead-zone varies oppositely with average flow velocity of molten bath. The minimum volume is 3.01 m 3 generated by 45 300 arrangement with the biggest average flow velocity being 0.080 m/s, and the maximum volume is 3.52 m 3 generated by 0 150 arrangement with the minimum average flow velocity being 0.071 m/s. The result shows that appropriate side-blowing arrangement could improve the stirring ability of side-blowing jet, and suppress the velocity dead-zone volume. Comparing with six kinds of side-blowing arrangement, the 45 300 arrangement could achieve the biggest average flow velocity of molten bath being 0.080 m/s, and the minimum volume of dead-zone being 3.01 m 3. 4.3. Industrial Application Research Based on the results of water experiment and numerical simulation, the 45 300 arrangement showed the best stirring effect on molten bath in vanadium extraction process. In order to analyze the metallurgical and operation benefits of the powder injection technology, both the 45 300 arrangement and none powder injection arrangement were adopted in a 150 t vanadium extraction converter, and the powder was the iron ore with diameter being 0.75 mm. Hereafter, none powder injection arrangement would be addressed as none arrangement. There were 160 heats collected in the industrial application research, including 80 heats with 45 300 arrangement and 80 heats with none arrangement, and the operation conditions were same with both the 45 300 and none arrangements in melting process. Semisteel component, vanadium extraction rate and T. Fe content in the slag were studied in this research. The semi-steel, which is before and after melting process, would be defined as semi-steel A and semi-steel B, respectively. The average components and temperatures of semi-steels are shown in Table 4. The initial conditions of semi-steel A is stable, which means the effect of semi-steel A conditions on the vanadium extraction rate could be negligible. Moreover, the smelting process for the vanadium extraction is same with both kinds of powder injection arrangement. Comparing with none arrangement, the content of vanadium in semi-steel B is dropped by 0.011 mass% with the 45 300 arrangement, and the vanadium extraction rate in the converter is improved by 4.8% from 81.0% to 85.8% with the temperature of semi-steel B increasing little. Figure 10 illustrates the distribution of vanadium in semi-steel B using two kinds of powder injection arrangements for vanadium extraction. When the none arrangement is used, the content of vanadium in semi-steel B distributes from 0.027 mass% to 0.061 mass%. And the content of vanadium in semi-steel B is from 0.021 mass% to 0.046 mass% with the 45 300 arrangement. Therefore, the 45 300 arrangement would decrease the content and fluctuation range of vanadium content in semi-steel B, which means the powder injection technology could improve the stablity of melting process. V 2 O 5, T. Fe, M. Fe, and FeO contents in slag are shown in Fig. 11. Compared with the none arrangement, stirring effect of molten bath would be enhanced with 45 300 arrangement. The content of V 2 O 5 is increased by 1.4 mass%, the content of FeO is decreased by 2.1 mass%, and the content of T. Fe loss is reduced by 2.0 mass% in slag. However, the content of M. Fe in slag is increased by 0.7 mass%. During the vanadium extraction process, fluidity of slag would be reduced with FeO content dropping. As a result, the M. Fe is hard to be transmitted into molten bath, because the slag becomes sticky due to the reduction of FeO. Table 4. Average values analysis of semi-steels with different powder injection arrangement. Semi-steel A Semi-steel B Label None arrangement 45 300 arrangement C (mass%) 4.46 4.46 V (mass%) 0.231 0.232 Si (mass%) 0.21 0.22 Ti (mass%) 0.15 0.15 Temperature (K) 1 586 1 585 C (mass%) 3.61 3.60 V (mass%) 0.044 0.033 Temperature (K) 1 659 1 650 Fig. 9. The relationship between the average flow velocity and velocity dead-zone. Fig. 10. Distribution of vanadium contents in semi-steel B. 2018 ISIJ 858

Fig. 11. Based on the results, it seems that with the iron ore powder injected into the molten steel, the temperature of molten bath could be control better, which is good for the thermodynamics condition of vanadium extraction. At the same time, the side-blow method would improve the stirring effect of molten bath and the mass transfer of vanadium, which agrees with the results of water experiment and numerical simulation. Above all, the powder injection technology could inhibit the condition that the carbon oxygen reaction well, and achieve both metallurgical and operation benefits. 5. Conclusions Distribution of V 2O 5, FeO, T. Fe and M. Fe contents in slag. In this paper, the effect of side-blowing arrangement on flow field of molten bath has been researched, and different flow rates, injection angles and injection depths of sideblowing nozzle have been adopt in the water experiment and numerical simulation. The mixing time, average flow velocity and velocity dead-zone is used to analyze the suitable side-blowing arrangement. With the industrial application research, an optimum side-blowing arrangement is determined for 150 t vanadium extraction converter. The main results of this study are summarized as follows: (1) In order to take full advantage of side-blow stirring ability, the side-blowing nozzle should be inserted far from bottom-blowing nozzle as much as possible, and the molten bath can be stirred more intensively, when the distance between side-blowing nozzle and bath surface increases, at the tested conditons. (2) Compared with the six kinds of side-blowing arrangement, the 45 300 arrangement could achieve the biggest average flow velocity of molten bath being 0.0329 m/s, and the minimum volume of dead-zone being 1.90 m 3. (3) The 45 300 arrangement can accelerate the smelting process in 150 t vanadium extraction converter. Compared with none side-blowing arrangement, the content of vanadium in semi-steel B with 45 300 arrangement is 0.033 mass% and vanadium extraction rate is increased by 4.8%. The content of V 2 O 5 is increased by 1.4 mass%, the content of FeO is reduced by 2.1 mass% and the content of T. Fe loss is dropped by 2.0 mass% in slag. However, slag becomes sticky due to the reduction of FeO, which makes the content of M. Fe in slag be increased by 0.7 mass%. Acknowledgements The authors would like to express their thanks for the support by the National Nature Science Foundation of China (NSFC 51334001), Fundamental Research Funds for the Central Universities (FRF-AS-17-005) and the National Key Technology R&D Program of the 12 th Five-Year Plan (12FYP 2012BAC27B01). REFERENCES 1) B. V. Slobodin, T. I. Krasnenko, E. M. Gorbunova and T. P. Sirina: Metallurgist, 48 (2004), 85. 2) L. N. Monakhov, S. V. Khromov and P. I. Chemousov: Metallurgist, 48 (2004), 381. 3) P. R. Taylor, S. A. Shuey, E. E. Vidal and J. C. Gomez: Miner. Metall. Process., 23 (2006), 80. 4) L. A. Smirnov, V. A. Rovnushkin, A. L. Smirnov and S. P. Burmasov: Steel Transl., 40 (2010), 472. 5) Y. Han, M. Zhang, Q. He, P. Wang and D. Zhang: Special Steel, 38 (2017), 27. 6) L. Liu, W. Sun, L. Zhang and X. Li: Chin. J. Eng., 38 (2016), 1590. 7) G. Wang, Y. Wang, B. Xie, X. Li, X. Guo and X. Li: J. Iron Steel Res. Int., 19 (2012), 36. 8) A. Kruskopf: Metall. Mater. Trans. B, 48 (2017), 619. 9) Y. Tago and Y. Higuchi: ISIJ Int., 43 (2003), 209. 10) B. Tang, X. Wang, Z. Zou and A. Yu: Can. Metall. Q., 55 (2016), 124. 11) F. Qian, R. Mutharasan and B. Farouk: Metall. Mater. Trans. B, 27 (1996), 911. 12) C. S. Kumar and A. S. Kumar: ISIJ Int., 46 (2006), 1171. 13) X. Zhou, M. Ersson, L. Zhong, J. Yu and P. Jönsson: Steel Res. Int., 85 (2014), 273. 14) I. Sumi, Y. Kishimoto, Y. Kikuchi and H. Igarashi: ISIJ Int., 46 (2006), 1312. 15) H. J. Odenthal, W. H. Emling, J. Kempken and J. Schlüter: Iron Steel Technol., 11 (2007), 71. 16) F. Liu, R. Zhu, Q. Wang and R. Bai: ISIJ Int., 55 (2015), 1633. 17) W. Malalasekera and H. K. Versteeg: An Introduction to Computational Fluid Dynamics, Pearson, Harlow, (2007), 1. 18) B. E. Launder and D. B. Spalding: Lectures in Mathematical Model of Turbulence, Academic Press, London, (1972), 124. 859 2018 ISIJ