Metallurgical and Materials Transactions B最新文献

The design and optimization of blast furnace (BF) profiles generally rely on limited empirical knowledge and experience. Such exercise can now be improved by means of process modeling and optimization. In this work, the effect of furnace profile on BF performance is numerically investigated across a wide range of furnace volumes (500 to 6000 m³), focusing on the ratio of effective height to belly diameter (H/D ratio). This is done based on the recently developed 3D multi-fluid BF process model. The results indicate that the optimized H/D ratio under different furnace volumes can be determined by minimizing the total energy consumption, namely, the summation of chemical and physical energy consumptions. Comparative analysis indicates that the industrial data on the variation of H/D ratio with furnace volume, collected over years, can be reproduced quantitatively by the current model. The increase in BF size or volume results in high thermal energy efficiency and low coke rate, primarily attributed to the reduced heat dissipation from the top gas and furnace wall. But there exists a size threshold between small and large BFs, approximately 2000 m³ under the conditions considered. Beyond this threshold, the BF performance and in-furnace states do not change significantly. Optimum BF profile needs to consider the coupled effect of H/D ratio and effective furnace volume.

To investigate the qualitative and quantitative effects of flow rate on fluid dynamics within a tundish using channel induction heating, a mathematical model was developed to reflect the macroscopic transport phenomena. The simulation results show that as inlet velocities increased from 0.3 to 0.9, the ratios of plug zone volume to total volume, dead zone volume to total volume, plug zone to dead zone, and plug zone to mixed zone consistently fell within the ranges of 0.257 to 0.263, 0.118 to 0.119, 2.152 to 2.231, and 0.412 to 0.427. Concurrently, the inclusion removal rate in the tundish's dead zone under channel induction heating increased with the inlet velocity. Additionally, analysis of the mean characteristics in the tundish system reveals that the mean residence time linearly decreases with increasing mean flow temperature and linearly increases with the inlet velocity. Furthermore, the inclusion removal rate linearly increases with the inlet velocity. This work enables precise quantitative and qualitative predictions and analyses of physical parameters in the electromagnetic induction heating tundish at various casting speeds, providing an efficient method for evaluating and optimizing the multi-physical fields within the metallurgical container.

In the current study, a three-dimensional mathematical model that integrated the large eddy simulation (LES) turbulent model, volume of fluid (VOF) multiphase model, and the dynamic mesh model, was developed to investigate the influence of the mold oscillation on the steel, slag, and air multiphase flow and the slag entrainment in a slab continuous casting mold. The results indicated that the mold oscillation led to an increase of about 0.35 m/s in the speed of molten steel at the impact point of the upper and lower circulation on the narrow surface. However, it reduced about 0.047 m/s in the maximum speed of the meniscus at 1/4 width of the mold. The Fast Fourier Transform (FFT) analysis revealed that the characteristic frequency of the speed fluctuation in the gravity direction of the meniscus near the narrow surface was 1.27 Hz without the consideration of the mold oscillation. Upon the application of the mold oscillation, a new characteristic frequency of 3 Hz, matching the mold oscillation frequency, emerged with twice the intensity of the original 1.27 Hz frequency. Moreover, the mold oscillation had little effect on the characteristic frequency of the speed fluctuation on the meniscus far away from the narrow surface. A user-defined function (UDF) was employed to quantify the number, size, and spatial distribution of entrained slag droplets. The net slag entrainment rate was reduced from 0.0152 to 0.0113 kg/s after the consideration of the mold oscillation, and the number of entrained slag droplets was also decreased. The effect of mold oscillation on the size distribution of entrained slag droplets and the occurrence location of the slag entrainment on the meniscus were not significant.

The effects of the Al₂O₃ content and MgO/Al₂O₃ ratios on the structure and viscosity of CaO–SiO₂–MgO–Al₂O₃ silicate-based slag were studied by the rotating cylinder method. The evolution of the silicate structure was analyzed by Fourier transform infrared spectroscopy. The results show that when the CaO/SiO₂ ratio is 1.20 and the MgO content is 10 mass pct, the Si–O–Al structure and [AlO₄]⁵⁻ tetrahedron structure of the slag increase with increasing Al₂O₃ content (17 to 22 pct), which increases the viscosity, break-point temperature, and activation energy of the slag. In addition, when the CaO/SiO₂ ratio is 1.20 and the Al₂O₃ content is 18, 20, and 22 mass pct, respectively, with the increase of MgO/Al₂O₃ ratio in the range of 0.4 to 0.7, the viscosity, break-point temperature, and activation energy of slag decrease. The analysis of the Raman spectra in the 800 to 1200 cm⁻¹ range indicates that complex silicate structures ([Si₂O₅]²⁻, [Si₂O₆]⁴⁻) disaggregate into simpler structures ([Si₂O₇]⁶⁻ and [SiO₄]⁴⁻) with the MgO/Al₂O₃ ratio increase. The proportion of structural units Q³ and Q² decreases, the proportion of Q¹ and Q⁰ increases, and the proportion of silicate structural units (Q⁰ + Q¹)/(Q² + Q³) increases, which indicates that the non-bridging oxygen content in the slag increases, resulting in a decrease in slag polymerization.

Production of direct reduced iron (DRI), particularly with green hydrogen, is a key pathway to the decarbonization of the iron and steel industry. However, the sticking tendency during the production of DRI creates serious operational issues and limits production outputs. Coating inert materials on the surface of iron ores can act as a barrier to effectively prevent the bonding between newly formed iron surfaces, and can interfere with the formation of iron whiskers. However, the principle of coating has not been systematically studied. This review covers the mechanism of sticking in both shaft furnaces and fluidized bed-based gaseous DRI production. The factors that influence the reduction kinetics and morphology, including physical and chemical ore properties, pellet induration conditions, and reduction conditions are summarized as well. Understanding the relationship between these factors and morphology change is critical to eliminating the sticking issues of DRI. Findings from this study suggest that coating with inert additives (e.g., metal oxides) can successfully prevent sticking in both shaft furnaces and fluidized bed processes. The types of additives and coating methods, the stage of reduction where the coating is applied, and reduction temperature will dramatically affect the coating performance. The outlook is discussed as well given the need for further work to improve the performance of coating (methods, timing, and cheaper alternatives), to further de-risk DRI technologies.

Understanding the formation of a slag rim on the mold meniscus is crucial for controlling surface defects in the initial shell. However, there is a scarcity of quantitative studies on this matter. This study has developed a comprehensive three-dimensional (3D) numerical model for analyzing the meniscus multi-phase flow, heat transfer, and solidification, considering mold oscillation. The 3D morphology of the solidifying shell and the slag rim in the meniscus region were accurately reproduced. The solidification depth (D_IS), solidification length (L_IS) and oscillation mark depth (D_OM) of the initial shell were used to quantify the morphological characteristics of the initial shell. The results confirmed that the formation position difference of the slag rim along the circumferential direction of the mold significantly affects the initial solidification and uniform growth of the shell. In the corner of the mold, the deeper overflow makes the oscillation mark extend 2.8 to 3.2 mm in the direction of casting. In addition, in order to quantitatively investigate the influence of the slag rim, a two-dimensional (2D) model is established with phenomena and parameters considered the same as those of the 3D model. According to the slag rim morphology obtained by the 3D model, in the 2D model, it is proposed to construct three slag rims with the same maximum thickness of 6 mm at 10, 20 and 35 mm above the meniscus (H_Rim). The simulation of initial shell morphology revealed that a lower formation position of the slag rim led to more severe overflow of molten steel from the meniscus, resulting in non-uniform continuous growth of the initial shell. This increases the likelihood of potential blockage in the liquid slag flow towards the slag channel between solidified shell and mold.

Electric arcs are a necessary heat source in many industrial processes that take place in Submerged Arc Furnaces (SAFs). Arcs exhibit non-linear electrical characteristics and behave in a complex manner. Therefore, an improved understanding of their behavior enables better control of furnace operation. Modeling of industrial arcs is a multiphysics process that involves simultaneously solving several coupled physical phenomena, such as electromagnetics, fluid dynamics, and heat transfer, including a radiative heat transfer from the plasma arc. Coupling fluid dynamics and electromagnetics is known as Magnetohydrodynamics (MHD). For practical applications, however, there are also simpler approaches to arc modeling, either based on simplified physical principles or empirical behavior. In this paper, a combined Cassie–Mayr model (CMM) and a channel arc model (CAM) are implemented and coupled with a submerged arc furnace electrical circuit model. The complete circuit model parameters such as resistances and inductances are estimated using modeling of a full size furnace, and then, actual measurements from a SAF are used to validate the models by comparing current and voltage waveform. Both models are then used to estimate harmonic distortion in a SAF for different arc current ratios, which should help operators to estimate the arc current in real time thus be able to lower and raise the electrode to keep operating conditions constant.

To explore the solutions of saving energy and cost of the rotary hearth furnace (RHF) direct reduction process, this paper constructed an energy consumption model, an economic evaluation model, and a carbon emission calculation model of the RHF direct reduction process. According to the actual production conditions of a steel plant, the influence of combustion air temperature and oxygen enrichment rate on the energy consumption, cost, and carbon emission of the RHF direct reduction process were studied. The calculation results show that for every 50 °C increase in the combustion air temperature, the process energy consumption, comprehensive cost, and carbon emission reduce by about 11 kgce, 42 CHY, and 44 kg, respectively. For every 2 pct increase in the oxygen enrichment rate of the combustion air, the corresponding values are about 10 kgce, 26 CHY, and 37 kg, respectively. In addition, the mathematical model established in this paper can be used to calculate the process energy consumption, cost, and carbon emissions under different raw material and fuel conditions, which is of great theoretical significance for the green and low-carbon transformation of the RHF direct reduction process.

Two ultrasonic modes, i.e., continuous and pulsed ultrasounds, were introduced into the directional solidification process of Cu_68.3Al_27.6Ni_4.1 alloy. A columnar-to-equiaxed structure transition occurred to primary β(Cu₃Al) phase within continuous ultrasonic field, which was accompanied with a grain size reduction by 7.5 times. Under pulsed ultrasound, β phase maintained the fine columnar structures with a similar grain size. In the former case, numerous β phase nucleation sites formed ahead of solid/liquid (S/L) interface because of the large local undercoolings induced by transient cavitation. Meanwhile, intensive acoustic streaming suppressed the liquid temperature gradient from 120 to 85 K/cm, which interrupted the solute transportation along heat flow direction and resulted in equiaxed microstructures. Under the intermittent ultrasonic action in the latter case, fewer nucleation sites were generated near S/L interface but small columnar β grains were split from the original ones under stable cavitation. Since no steady convection was driven, the liquid temperature gradient of 110 K/cm remained almost constant, making those grains grow into refined columnar structures. Under the action of pulsed ultrasound, the yield strength was enhanced by a factor of 1.5 because of grain refinement strengthening, together with 94 pct shape recovery rate due to columnar grain structures.

In the present study, a well-known Iida’s equation for surface tension was modified to improve the predictivity of the surface tension of pure liquid metals. A semi-empirical equation for the surface tensions (({sigma }_{m})) of liquid metal at its melting temperature proposed by Iida et al. uses a generalized (alpha ) value of 0.51 to represent the ratio of the distance required to separate one atomic pair from its equilibrium distance. This study improved the predictability of the equation by refining the (alpha ) value using the equilibrium interatomic distance (({r}_{e})) and atomic radius (({r}_{a})). Assigning an accurate (alpha ) value for each element greatly improves the prediction accuracy of the surface tension for liquid metals. Furthermore, the critical temperature (({T}_{c})) was calculated based on the interatomic distance (({r}_{c})) at ({T}_{c}) and temperature coefficient of density ((d{rho }_{T})/(dT)) and used to predict the temperature dependence coefficient of surface tension ((d{sigma }_{T})/(dT)). As results, more accurate surface tensions of 42 liquid metals were predicted over the entire liquid state temperature.

Graphical Abstract