<<

Hindawi Publishing Corporation Mathematical Problems in Volume 2016, Article ID 4170371, 8 pages http://dx.doi.org/10.1155/2016/4170371

Research Article Dimensional Analysis of Average Diameter of Bubbles for Bottom Blown Furnace

Dongxing Wang, Yan Liu, Zimu Zhang, Shao, and Ting’an Zhang

School of of Northeastern University, Key Laboratory of Ecological Utilization of Multi- Intergrown of Ministry, Shenyang 110819, China

Correspondence should be addressed to Yan Liu; [email protected]

Received 8 March 2016; Accepted 21 July 2016

Copyright © 2016 Dongxing Wang 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.

Average diameter of bubbles is important in copper furnace. Based on the principle of similarity, a slice model of a furnace with bottom-blown oxygen in matte- process was established. A high-speed camera was used to record images continuously and clearer pictures were selected for treatment. Finally, image processing software was used for obtaining the average diameter of the bubbles. The effects of different injection conditions and equipment factors such as the diameter of nozzle, the nozzle installing angle, and gas velocity on the average diameter of bubbles were studied with cold water model experiment, exploring the dispersion and disintegration rules of bubbles. According to experimental data and Buckingham’s theorem, by using dimensional analysis 0.29374 −0.46572 −0.16725 method, an empirical formula on average diameter of bubbles was established (𝑑𝐵 = 0.41666𝑑 𝜃 V ). It can be seen from the formula that nozzle installing angle and diameter of nozzle make the most impact on the average diameter of bubbles in bottom blown oxygen copper furnace.

1. Introduction of bottom blowing is only 10 to 20 seconds, and top and bottom combined blowing is 20∼50 seconds. From this, we About 80% of copper is produced by pyrometallurgical can see a significant advantage of bottom blowing. process. The other 20% is obtained by . Pyro- The bottom blown technique has been widely used in metallurgical process includes matte smelting, converting, the , such as oxygen bottom blown converter; , and electrorefining of blister copper [1]. The objective of the smelting process is to melt concentrate and oxygen is blown into the molten pool through the bottom oxidizeSandFeintheconcentrateandproducematte(aCu- of the converter and makes molten smelting into steel. enriched molten phase) containing 50∼75% copper. Thetechniqueofblowingargonfromthebottomofladle The molten matte is processed in a converter resulting in takes advantage of the stirring effect of bottom blown argon, crudecopperof98.5∼99.5% copper content [2]. Converting is making the composition and temperature of liquid steel mostly carried out in the cylindrical Peirce-Smith converter. uniform and the inclusions and bubbles in steel rise quickly, Oxygen-enriched air is often used as the oxidant. thus increasing the purity of liquid steel [3, 4]. In the field Bathsmeltingisoneofthemostimportantmattesmelting of nonferrous metallurgy industry, the application of bottom techniques with oxygen-enriched air blown into the matte, blown technique is very limited, which is mainly used for including top blown such as Ausmelt and Isa smelt, side smelting and copper smelting [5–8]. blown such as Noranda and Teniente process, and bottom In recent years, oxygen bottom blown copper smelting blown such as SKS process. (SKS) has been used more and more widely in Injection technology is widely used in hot metal pretreat- China. SKS includes two generations of technology; the first ment,strengtheningthereactioninfurnace,secondaryrefin- is “bottom blown smelting-PS furnace converting” and the ing of liquid steel, and other fields. In converter , second is “bottom blown smelting-bottom blown continuous for example, mixing time of top blowing is 90 to 120 seconds, converting.” 2 Mathematical Problems in Engineering

SKSprocesshasmanyadvantagessuchashighoxygen coalesce to form a horizontal unstable gas envelope. enrichment, low temperature, and autothermal operation. It If the submergence is shallower than 300 mm, bubbles has a very strong adaptability of raw materials; the prepara- may break through the bath surface and form a vertical tion of is very simple. It can directly deal with channel. Then the oxygen utilization is poor. the wet material, block material, and garbage material. Its Valencia et al. [15] had numerically studied the fluid thermal efficiency is very high, without additional fuel. This dynamics in a model of Teniente type copper converter new technology has great potential for expanding production with combined lateral and bottom injection. The volume capabilityandenergysaving[8–10].Sincethefirstbottom fraction of water distribution, velocity vectors, agitation of blownoxygencopperfurnacewasputintoproductionin bathsurface,andgaspenetrationinthreedifferentmodels 2008, there will be more than 8 bottom blown furnaces was studied. With the help of bottom injection, the model running in China in 2015. produced more jet interaction, bath inner agitation, and Bottom blown oxygen copper smelting process is a new mixing.ZhaoandIrons[16]usedK-HandR-Tinstability technology whose significant feature is the application of high theory to explain the breakup of a bubble forming at a speed bottom jet. From the viewpoint of metallurgical reac- submergednozzle.Theypresentthefirsttheorytopredict tion engineering, the bath smelting method can be divided the transition from bubbling to jetting for submerged gas into two kinds, namely, the bubble bath smelting and jetting injection. bath smelting. When the oxygen-enriched air is injected from Lou et al. [17] studied the formation and motion rules of tuyere into the bath at a low velocity, the air is actually sprayed bubbles in Isasmelt furnace by water model and numerical into the molten pool at a pulse type, belonging to the bubble simulation method. Their conclusion is that, with the increase generation system, namely, the bubble bath smelting. When of gas flow rate 𝑄, the bubble frequency decreases and the bubble diameter increases. When the gas flow rate 𝑄 is equal the oxygen-enriched air is blown from tuyere into the bath 3 at a high velocity, approaching or exceeding sonic velocity, to 24 m /h, the equivalent diameter 𝑑 of the bubble is about the air is blown into the molten pool at stability state of 25 mm, and residence time of bubble ranges from 0.14 to continuous injection, belonging to the jet system, namely, the 0.16 s. jetting bath smelting. Bottom blown oxygen copper smelting process is still in a For jetting blown, the air separates into many small development stage with many shortcomings and deficiencies bubbles. The smaller bubble diameter means the larger gas- [8–10]. It is high speed jetting blowing and different from liquid boundary area at the same quantity of gas. With small side-blown Peirce-Smith copper converter and other tradi- bubbles, molten pool has no fierce nuisance spitting, so small tional bubble bath smelting. The mechanisms of high speed bubblesarehelpfultoextendtheservicelifeoftherefractory bottom jet still need to be studied which is very important material. for the improvement of the process. The air injected into the moltenpoolwillbreakupintomanysmallbubbles;thereis In bath smelting process with gas injection, the size, very little research about the process, such as the dispersion rising velocity, and dispersion of the bubbles, and the gas and size of bubbles, the factors which influence bubbles size holdup have enormous influence on the interface reaction and their relationship. or multiphase reaction, directly affecting the reaction rate In this paper, a water model of the bottom blown furnace and oxygen utilization. Besides, the residence time of smaller was produced to study the mean bubble diameter in the SKS bubbles is longer. process. Some effects of variables on the bubbles dispersions Alargenumberofstudieshavebeencarriedoutonthe and gas holdup were researched in our former paper [18]. investigation of the flow regimes in submerged gas injection. This paper mainly analyzes the relationship between the Davidson and coworkers [11, 12] formulated a theoretical parameters and derived the empirical equation of mean model of continuous bubble formation at submerged orifices bubble diameter. From empirical equation, the effects of in water model. They resulted in the relationship between the diameter of nozzle, the nozzle installing angle, and gas 𝑉 𝐺 𝑉 𝐺1.2 −0.6 bubble volume, , and gas flow rate, : =1.138 g . velocity on the average diameter of bubbles were studied. A camera was used to take the photographs of the bubble formation. The mean volume of bubbles measured from the photographs. Hoefele and Brimacombe [13] investigated 2. Experimental the behavior of gas injected into liquids by high-speed cinematography and pressure measurement in the tuyere. 2.1. Principles. For geometric similarity, a water model with They found that matte converters are in the bubbling regime. a scale down of 1 : 5 was built. The main geometrical sizes of Testsonanoperatingnickelconverterconfirmedthatwith the model were proportional to those of the prototype. Dynamic similarity relies on the fact that the modified low pressure blowing practice air discharges into the bath 󸀠 in the form of large bubbles whose corresponding diameters Froud number Fr between the model and prototype should are 400 mm at room temperature. The bubble frequency is be equal. Using (1), the quantity of gas injected into the model −1 between 10 and 12 s . Bustos et al. [14] thought the blockage can be calculated: of tuyeres and the life of the lining in a side-blown 4 1/2 air copper converter are related to the injection of air. Bubble 𝜌𝑔𝑝 𝜌 𝑑 𝐻 𝑄 =( ⋅ 𝑙𝑚 ⋅ 𝑚 ⋅ 𝑚 ) ⋅𝑄 , (1) volume and equivalent radius must be considered. In the 𝑚 𝜌 𝜌 𝑑4 𝐻 𝑝 Peirce-Smith copper converter, bubbles growing at adjacent 𝑔𝑚 𝑙𝑝 𝑝 𝑝 Mathematical Problems in Engineering 3

Table 1: Attributes of model and materials. Radius of the Width of the Height of liquid Density of Density of System −3 −3 pool/mm pool/mm level/mm liquid/(kg⋅m ) gas/(kg⋅m ) Water model 342 300 342 1000 1.25 Prototype 1710 15000 1710 4600 1.37

Main Feed inlets Gas outlet oil jet Optional oil jet

Figure 1: Bottom blown oxygen copper smelting furnace.

Figure 2: Sketch of experimental setup. where 𝜌 stands for the density; 𝑑 is the diameter of nozzle; 𝐻 is the depth of the pool; 𝑄 stands for the quantity of injected gas; the subscripts 𝑝, 𝑚, 𝑙,and𝑔 standfortheprototype, model liquid, and gas, respectively. must be dimensionally homogeneous, which means both sides of an equation must have the same units [20]. In engi- neering, empirical formula between the parameters can be 2.2. Apparatus. The water model used in our experiment derived from experiment data by the method of dimensional references a company’s bottom blown oxygen copper smelt- analysis which also enabled weighing the magnitude of each ingfurnaceasshowninFigure1.Theexternaldiameter operating parameter [21]. For example, by the means of of the furnace is 4.4 m which lined with 380 mm thickness dimensional analysis, Yan et al. [22] obtained a penetration refractory bricks and the length is 16.5 m. Nine nozzles are depth criterion equation of 1/2 anode structure electrolytic installed on the bottom of the furnace in two rows. For the cell. convenience of observation, a cross-sectional model which contains only one nozzle was built up. The specifications areshowninTable1.Thematerialofthemodelisorganic 3. Dimensional Analysis of Average Diameter . Inside diameter of the model is 684 mm and as a cross- of Bubbles sectional model the length is 300 mm, so that the experiment phenomena can be easily observed from the cross section. Through water model experiments, the average diameter of Water model is shown in Figure 2. There are different bubbles was obtained. By the method of dimensional analysis, positions on the bottom of the model for nozzles to be the influence of different experimental conditions on the installed. averagediameterofbubblescanbestudiedandthefitting equation was obtained. Considering the main factors that affect bubbles in 2.3. Experimental Procedure and Methods. The model was 𝜋 half full with water and air was injected into the water from the experiments, according to Buckingham’s theorem ( - principle), the independent variables were found out. the nozzles on the bottom of the model. Combined with high- 𝑑 speed camera, image processing technology is an important The average diameter of bubbles 𝐵 was mainly affected method to study the characteristics of bubbles [19]. A high- by the following factors: speed camera was used to record the experimental phenom- (1) The diameter of the nozzle 𝑑. ena continuously; single frame images with bubbles could be 𝐻 obtained for the water model experiment. Clearer pictures (2) The height of liquid level . were selected, then identifiable bubbles were marked, and the (3) The angle between the nozzle and horizontal plane 𝜃. edges were reinforced treatment. Finally, image processing (4) The diameter of the model 𝐷. software was used for obtaining the mean diameter of the 𝑃 bubbles. (5) The gas pressure . (6) The viscosity of water 𝜇𝑙.

2.4. Dimensional Analysis. Dimensional analysis is an impor- (7) The viscosity of gas 𝜇𝑔. tant method to research the relationships of variables that 𝜎 describe the physical phenomena. The foundation of dimen- (8) The gas-liquid interfacial tension 𝑔-𝑙. sional analysis is that the empirical and rational equations (9) The density of water 𝜌𝑙. 4 Mathematical Problems in Engineering

𝛼8 𝛽8 𝛾8 Table 2: Dimension of variable. Π8 =𝜌𝑙 𝐻 𝜇𝑙 𝑔, 𝑑 𝑑𝐻𝜃𝐷 𝑃𝜇 𝜇 𝜎 𝜌 𝜌 𝑔 V 𝛾 𝐵 𝑙 𝑔 𝑔-𝑙 𝑙 𝑔 𝛼9 𝛽9 9 Π9 =𝜌𝑙 𝐻 𝜇𝑙 V. 𝑀 000001 1 1 1 1 1 0 0 𝐿 11101−1 −1 −10−3 −31 1 (6) 𝑇 00000−2 −1 −1 −200−2 −1 For Π4, plug in the dimensions of every variable; equation can be obtained as [𝑀0𝐿0𝑇0] (10) The density of gas 𝜌𝑔. (7) 𝑔 −3 𝛼4 𝛽 −1 −1 𝛾4 −1 −2 (11) The acceleration of gravity . =[𝑀𝐿 ] [𝐿] 4 [𝑀𝐿 𝑇 ] [𝑀𝐿 𝑇 ]. (12) The gas velocity V. The exponential equation set can be obtained as By the above analysis, using dimensional analysis we can draw general functional form as follows: 𝑀: 0=𝛼4 +𝛾4 +1, 𝑑 = 𝑓 (𝑑, 𝐻, 𝜃, 𝐷,𝑃,𝜇 ,𝜇 ,𝜎 ,𝜌,𝜌 ,𝑔,V). 𝐿 0=−3𝛼 +𝛽 −𝛾 −1, 𝐵 𝑙 𝑔 𝑔-𝑙 𝑙 𝑔 (2) : 4 4 4 (8)

Or 𝑇: 0=−𝛾4 −2.

𝑓(𝑑𝐵,𝑑,𝐻,𝜃,𝐷,𝑃,𝜇𝑙,𝜇𝑔,𝜎𝑔 𝑙,𝜌𝑙,𝜌𝑔,𝑔,V)=0. (3) Solve for 𝛼4 =1, 𝛽4 =2,and𝛾4 =−2.SoΠ4 = - 2 −2 2 2 𝜌𝑙𝐻 𝜇𝑙 𝑃=𝜌𝑙𝐻 𝑃/𝜇𝑙 . The variable dimension is listed in Table 2. 𝜋 Similarly, By the analysis of -principle, the total number of 𝜇 variables is 13 and the number of independent variables is 𝑔 Π5 = , 𝑘=3,sothat𝑛−𝑘= 10dimensionless quantity can be set up 𝜇𝑙 𝜌 𝐻 𝜇 [22–24]. 𝑙, ,and 𝑙 were selected as independent variables: 𝜌 𝐻𝜎 −2 𝑙 𝑔-𝑙 𝛼 𝛽 𝛾 1 −3 0 Π =𝜌𝐻𝜇 𝜎 = , 𝜌 =𝑀 1 𝐿 1 𝑇 1 =𝑀𝐿 𝑇 , 6 𝑙 𝑙 𝑔-𝑙 2 dim 𝑙 𝜇𝑙

𝛼2 𝛽2 𝛾2 0 1 0 𝜌 dim 𝐻=𝑀 𝐿 𝑇 =𝑀𝐿 𝑇 , (4) 𝑔 Π7 = , (9) 𝜌𝑙 𝛼3 𝛽3 𝛾3 1 −1 −1 dim 𝜇𝑙 =𝑀 𝐿 𝑇 =𝑀𝐿 𝑇 , 𝜌2𝐻3𝑔 󵄨 󵄨 󵄨 󵄨 Π =𝜌2𝐻3𝜇−2𝑔= 𝑙 , 󵄨𝛼 𝛽 𝛾 󵄨 󵄨1−30󵄨 8 𝑙 𝑙 2 󵄨 1 1 1󵄨 󵄨 󵄨 𝜇 󵄨 󵄨 󵄨 󵄨 𝑙 󵄨𝛼 𝛽 𝛾 󵄨 = 󵄨01 0󵄨 =−1=0.̸ 󵄨 2 2 2󵄨 󵄨 󵄨 (5) 𝜌 𝐻V 󵄨 󵄨 󵄨 󵄨 Π =𝜌𝐻𝜇−1V = 𝑙 . 󵄨𝛼3 𝛽3 𝛾3󵄨 󵄨1−1−1󵄨 9 𝑙 𝑙 𝜇𝑙 𝜌 𝐻 𝜇 Because (5) is not equal to 0, 𝑙, ,and 𝑙 are independent So of each other and it is reasonable that they were selected as 2 𝜇 𝜌 𝐻𝜎 𝜌 2 3 𝑑𝐵 𝑑 𝐷 𝜌𝑙𝐻 𝑃 𝑔 𝑙 𝑔-𝑙 𝑔 𝜌𝑙 𝐻 𝑔 independent variables. For the variables 𝑑𝐵, 𝑑,and𝐷,asthey 𝑓( , , ,𝜃, , , , , , 𝐻 𝐻 𝐻 𝜇2 𝜇 𝜇2 𝜌 𝜇2 only contain length dimension, therefore in the construction 𝑙 𝑙 𝑙 𝑙 𝑙 of dimensionless Π, they can be directly expressed with the (10) 𝜌 𝐻V dimensionless independent variable 𝐻. Each Π is expressed 𝑙 )=0. as follows: 𝜇𝑙 𝑑 𝐵 𝑑𝐵 Π0 = , In order to get the expression of ,theaveragediameter 𝐻 of bubbles can be expressed as an explicit function: 𝑑 2 𝜇 𝜌 𝐻𝜎 𝜌 2 3 𝑑 𝑑 𝐷 𝜌𝑙𝐻 𝑃 𝑔 𝑙 𝑔 𝑙 𝑔 𝜌 𝐻 𝑔 Π1 = , 𝐵 =𝑓 ( , ,𝜃, , , - , , 𝑙 , 𝐻 𝐻 1 𝐻 𝐻 2 𝜇 2 𝜌 2 𝜇𝑙 𝑙 𝜇𝑙 𝑙 𝜇𝑙 𝐷 (11) Π2 = , 𝜌 𝐻V 𝐻 𝑙 ). 𝜇𝑙 Π3 =𝜃, In addition, for the specific circumstances of our exper- 𝛼 𝛽 𝛾 4 4 4 𝐷 𝐻 𝑃 𝜌 𝜇 𝜌 𝜇 𝜎 𝑔 Π4 =𝜌𝑙 𝐻 𝜇𝑙 𝑃, imental study, where , , , 𝑙, 𝑙, 𝑔, 𝑔, 𝑔-𝑙,and are quantitative, (11) can be written as 𝛼5 𝛽5 𝛾5 Π5 =𝜌 𝐻 𝜇 𝜇𝑔, 𝑙 𝑙 𝑑𝐵 𝑑 𝜌𝑙𝐻V =𝑓2 ( ,𝜃, ). (12) 𝛼 𝛽 𝛾 𝐻 𝐻 𝜇 Π =𝜌 6 𝐻 6 𝜇 6 𝜎 , 𝑙 6 𝑙 𝑙 𝑔-𝑙 The above formula needs to be determined by experi- 𝛼7 𝛽7 𝛾7 Π7 =𝜌𝑙 𝐻 𝜇𝑙 𝜌𝑔, ments. Mathematical Problems in Engineering 5

4. The Derivation of Empirical Formula 4.4 4.2 According to the established equations (11) and (12), com- bined with the specific conditions, the specific empirical 4.0 formula can be derived [22–24]. 3.8 In various phenomena, criterion equation can be expressed in the form of power function with the independ- 3.6 ent variable within a certain range. By the above analysis, 3.4 empirical formula can be fitted:

Calculated results Calculated 3.2 𝑏1 𝑏3 𝑑𝐵 𝑑 𝑏 𝜌𝑙𝐻V =𝐴( ) (𝜃) 2 ( ) . (13) 3.0 𝐻 𝐻 𝜇𝑙 2.8 Among them 𝐴, 𝑏1, 𝑏2,and𝑏3 are fitting coefficient. Do logarithmic operations on the above formula; then 2.6 2.6 2.8 3.0 3.2 3.4 3.6 3.8 4.0 4.2 4.4 𝑑 𝑑 𝜌 𝐻V 𝐵 = 𝐴+𝑏 +𝑏 𝜃+𝑏 𝑙 . Experimental results ln 𝐻 ln 1 ln 𝐻 2 ln 3 ln 𝜇 (14) 𝑙 Figure 3: Relationship between experimental results and calculated Based on the linear relationship of the above form, fitting results. was carried out on the data of our injection experiment process [18]. The above fitting coefficients were obtained: ln 𝐴 = 5. Discussions 1.98584, 𝐴 =7.28519,𝑏1 = 0.29374, 𝑏2 = −0.46572,and𝑏3 = −0.16725.So 5.1. Effects of Nozzle Diameter on the Average Diameters of Bubbles. The gas velocity was 340 m/s, and the angle between 𝑑 𝑑 ∘ 𝐵 = 7.28519 + 0.29374 − 0.46572 𝜃 the nozzle and horizontal plane 𝜃 is 90 . With different ln 𝐻 ln 𝐻 ln diameters of the nozzles, the mean diameters of bubbles were (15) 𝜌 𝐻V measured in the experiment. From (18) and experimental − 0.16725 𝑙 . ln conditions, (19) was obtained: 𝜇𝑙 0.29374 So 𝑑𝐵 = 0.019331𝑑 . (19) 0.29374 −0.16725 𝑑 𝑑 −0.46572 𝜌𝑙𝐻V 𝐵 = 7.28519 ( ) (𝜃) ( ) . (16) The predicted values by (19) and experimental results 𝐻 𝐻 𝜇𝑙 wereshowninFigure4.Whenthediameterofthenozzle Then the empirical formula on average diameter of bub- increases, the mean diameter of bubbles has a trend of bles was established as increase.

𝑑𝐵 5.2. Effects of Angle between the Nozzle and Horizontal Plane 0.29374 −0.16725 (17) 𝑑 −0.46572 𝜌𝑙𝐻V ontheAverageDiametersofBubbles.The gas velocity was = 7.28519𝐻 ( ) (𝜃) ( ) . 340 m/s, and the diameter of the nozzle was 3.0 mm. With 𝐻 𝜇𝑙 different angles between the nozzles and horizontal plane, the In this formula, 𝐻 is a constant value which is equal mean diameters of bubbles were measured in the experiment. 𝜌 to 0.36 m. 𝑙 is the density of water which is equal to From (18) and experimental conditions, (20) was obtained: ⋅ −3 𝜇 × 1000 kg m . 𝑙 is the viscosity of water which is equal to 1.0 −0.46572 −3 −2 𝑑 = 0.028531𝜃 . 10 N⋅s⋅m .So 𝐵 (20) 0.29374 −0.46572 −0.16725 𝑑𝐵 = 0.41666𝑑 𝜃 V . (18) The predicted values by (20) and experimental results wereshowninFigure5.AsillustratedinFigure5,withthe Itcanbeseenfromthedimensionlesscriterionequation increase of the angles, the mean diameter of bubbles has a that,amongthefactorsdiscussed,nozzleinstallingangleand decreasing tendency. diameter of nozzle make the most impact on the average diameter of bubbles in bottom blown oxygen copper furnace. 5.3. Effects of Gas Velocity on the Average Diameters of Bubbles. There is a negative correlation between the nozzle installing 𝜃 ∘ angle and average diameter of bubbles. There is a positive The angle between the nozzle and horizontal plane was 90 correlationbetweenthediameterofnozzleandtheaverage and the diameter of the nozzle was 3.0 mm. With different gas diameter of bubbles. velocity, the mean diameters of bubbles were measured in the Use the experimental results and the calculated results as experiment. From (18) and experimental conditions, (21) was theabscissaandordinate,respectively,toplotagraphand obtained: compare with 𝑦=𝑥.TheresultsareshowninFigure3.The −0.16725 𝑑𝐵 = 0.0093021V . (21) experimental results accord with the calculated results well. So the model can predict the mean bubbles diameter in the The predicted values by (20) and experimental results water model of bottom blown furnace. were shown in Figure 5. It can be seen from Figure 6 that 6 Mathematical Problems in Engineering

4.2

4.0 3.6

3.8

3.4 (mm)

(mm) 3.6 B B d d 3.4 3.2 3.2

3.0 300 350 400 450 500 550 600 2.0 2.5 3.0 3.5 4.0 4.5 5.0 (m/s) d (mm) Measured Measured Predicted Predicted Figure 4: Effects of nozzle diameter on the average diameters of Figure 6: Effects of nozzle diameter on the mean diameters of bubbles. bubbles.

4.4 is difference between cold model results and the industrial situation, cold model is an interim solution to research high 4.2 temperature metallurgical reactor. The cold water model we used follows the principle 4.0 of similarity, including geometric similarity and dynamic similarity.KeepingthesameFroudenumberallowsthe 3.8 (mm) fluid flow patterns obtained in cold model similar with the B

bubble,” Journal of Fluid Mechanics, vol. 12, no. 3, pp. 408–416, [28] L.Li,Z.Liu,B.Li,H.Matsuura,andF.Tsukihashi,“Watermodel 1962. and CFD-PBM coupled model of gas-liquid-slag three-phase [12] J. K. Walters and J. F. Davidson, “The initial motion of a flow in ladle metallurgy,” ISIJ International,vol.55,no.7,pp. gasbubbleformedinaninviscidliquidPart2.Thethree- 1337–1346, 2015. dimensional bubble and the toroidal bubble,” Journal of Fluid [29] Y.Liu, F.-Y.Su, Z. Wen, Z. Li, H.-Q. Yong, and X.-H. Feng, “CFD Mechanics,vol.17,no.3,pp.321–336,1963. modeling of flow, temperature, and concentration fields in a [13] E. O. Hoefele and J. K. Brimacombe, “Flow regimes in sub- pilot-scale rotary hearth furnace,” Metallurgical and Materials merged gas injection,” Metallurgical Transactions B,vol.10,no. Transactions B,vol.45,no.1,pp.251–261,2014. 4, pp. 631–648, 1979. [14] A. A. Bustos, G. G. Richards, N. B. Gray, and J. K. Brimacombe, “Injection phenomena in nonferrous processes,” Metallurgical Transactions B,vol.15,no.1,pp.77–89,1984. [15] A. Valencia, M. Rosales-Vera, and C. Orellana, “Fluid dynamics in a teniente type copper converter model with one and two tuyeres,” Advances in Mechanical Engineering,no.1,pp.323–335, 2013. [16] Y.-F. Zhao and G. A. , “The breakup of bubbles into jets during submerged gas injection,” Metallurgical Transactions B, vol. 21, no. 6, pp. 997–1003, 1990. [17]W.-T.Lou,B.-Q.Zhang,andZ.Shi,“Simulationofbubble dynamics in water model of ISA vessel,” China Nonferrous Metallurgy,no.1,pp.48–53,2010. [18] D.Wang,Y.Liu,Z.Zhang,P.Shao,andT.Zhang,“Experimental study of bottom blown oxygen copper smelting process for water model,” in Proceedings of the 7th International Conference on Micromechanics of Granular Media (Powders and Grains), vol.1542,pp.1304–1307,,Australia,July2013. [19] L.-L. Wang, X.-D. Yong, R. Li, and Q. Ma, “Image process- ing technology in the application of bubble characteristics research,” Journal of Sichuan University (Engineering Science Edition), vol. 44, no. 2, pp. 188–192, 2012. [20] E.-R. Avram, E. Mamut, L. Oancea, and C. Bacu,ˆ “The method of dimensional analysis using similitude theory applied for modeling of fluid flow in microchannels,” in Proceedings of the 7thInternationalSymposiumonAdvancedTopicsinElectrical Engineering (ATEE ’11), pp. 1–6, Bucharest, Romania, May 2011. [21] G. Mary, S. Mezdour, G. Delaplace, R. Lauhon, G. Cuvelier, and F. Ducept, “Modelling of the continuous foaming operation by dimensional analysis,” Chemical Engineering Research and Design, vol. 91, no. 12, pp. 2579–2586, 2013. [22]Y.Liu,Y.Li,T.Zhang,andN.Feng,“Researchonthepenetra- tion depth in aluminum reduction cell with new type of anode and cathode structures,” JOM,vol.66,no.7,pp.1202–1209,2014. [23] Z. Zhang, The Principle of Dimensional Analysis,People’sRail- way Press, Beijing, China, 1979. [24] Y. Liu, T. Zhang, Q. Wang, W. Yao, M. Sano, and J. He, “Dimen- sion analysis theory in the application of fluid mechanics,” Industrial Furnace,vol.29,no.6,pp.9–12,2007(Chinese). [25] G. Solorio-D´ıaz, R. D. Morales, and A. Ramos-Banderas, “Effect of a swirling ladle shroud on fluid flow and mass transfer in a water model of a tundish,” International Journal of Heat and Mass Transfer, vol. 48, no. 17, pp. 3574–3590, 2005. [26] M. Iguchi, H. Kawabata, Y. Ito, K. Nakajima, and Z.-I. Morita, “Continuous measurements of bubble characteristics in a molten iron bath with Ar gas injection,” ISIJ International,vol. 34, no. 12, pp. 980–985, 1994. [27] Y.Liu,M.Sano,T.A.Zhang,Q.Wang,andJ.He,“Intensification of bubble disintegration and dispersion by mechanical stirring in gas injection refining,” ISIJ International,vol.49,no.1,pp. 17–23, 2009. Advances in Advances in Journal of Journal of Operations Research Decision Sciences Applied Mathematics Algebra Probability and Statistics Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014

The Scientific International Journal of World Journal Differential Equations Hindawi Publishing Corporation Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014