Search Documents

By L. Hekkala, J. Ollila, J. Keihäs, P. Mäkelä

The Outokumpu Steel Belt Sintering (SBS) Technology is used for manufacturing chromite pellets that are charged into a smelting furnace for ferrochromium production. One of the main elements in successful furnace operation is the control of gas flows inside the furnace, giving improved temperature distribution in the bed during the process. Computational Fluid Dynamics (CFD) gives excellent opportunities for studying these phenomena. An investigation using CFD can be performed relatively fast. For example, implications of different configurations can be studied within a day and an optimum design chosen. Therefore, Outokumpu Technology has developed a Steel Belt Sintering Model (SBSM), which takes advantage of a commercial CFD code. Each main physicochemical event occurring in a steel belt sintering furnace is incorporated independently in the SBS model as a sub-model, enabling every event to be studied separately. The commercial CFD code solves mass, momentum, species and energy equations for determination of the gas flow, composition and temperature. In addition, the Outokumpu SBSM takes into account the moving pellet bed, drying of pellets interms of water content, coke combustion, oxidation and reduction of chromite and the decomposition of carbonates and silicates. During the development of the model, numerous laboratory-scale tests were performed. The sub-models were validated against the laboratory tests. The reliability of the Outokumpu Steel Belt Sintering Model was validated against a pilot-scale batch-sintering furnace by comparing the modelling results against measurements. The aim of this paper is to give a brief presentation of the Outokumpu steel belt sintering model by showing some general modelling examples. Emphasis has been placed on demonstrating some of the main features of the model.

Jan 1, 2007

By S. Eckert, M. Reuter, D. Obiso, S. Kriebitzsch, M. Akashi

The present paper focuses the application of CFD techniques to investigate the Top-Submerged- Lance (TSL) gas injection in liquid metal. Previous works of the authors have shown that up- and down- scaling procedures based on the modified Froude number have some shortcomings, as this approach does not take into account the interfacial and viscous forces. Indeed, surface tension and dynamic viscosity of the smelting slags (σ = 0.4-0.5 N/m, μ = 0.2 Pa·s) are higher than the operating fluids that have been used in literature (water, paraffin oil) to study TSL injection in down-scaled furnaces. In order to get closer to real systems, the authors study the TSL injection of Argon in a liquid metal. An experimental campaign was carried out in the X-ray lab at Helmholtz-Zentrum Dresden-Rossendorf (HZDR), where the eutectic alloy GaInSn was used as liquid phase. The alloy is liquid at room temperature, and X-Ray imaging is used to picture the multiphase flow in a quasi-2D vessel. The aim of the present work is to demonstrate the applicability of CFD techniques to model multiphase flows involving liquid metals, and validate the model using the data produced at HZDR. The commercial software ANSYS Fluent® was used together with the Volume of Fluid model to directly resolve the gas-liquid interphase. Some features of the flow, such as the void fraction distribution and bubble detachment frequency are tracked with CFD and compared to the experimental data. INTRODUCTION In the last thirty years the TSL smelting technology has increasingly been used for industrial non- ferrous smelting processes and is now one of the significant pyrometallurgical processes. Developed in the early 1970’s at CSIRO, the furnace is nowadays commercialized in different market sectors, including copper and lead smelting, zinc fuming, waste treatment and slag recovery (Floyd, 2005). Despite the wide application of these reactors in the pyrometallurgical market, some aspects of the process are still not fully clarified, since in-situ measurements and optical accesses are prohibitive because of the high temperatures and aggressive conditions. Phenomena such as the bubbling regime, the slag solidification and the submerged combustion still need a deeper scientific investigation to obtain a better understanding of the process.

Jan 1, 2019

By A. Ochieng

Solids suspension influences the quality of mixing and energy requirement in a solid-liquid system, both of which determine the efficiency of industrial processes such as nickel precipitation. Nickel solids concentration distribution in a stirred tank was investigated using computational fluid dynamics (CFD) and experimental methods. The concentration distribution of the nickel solids was compared with that of sand and glass. The laser Doppler velocimetry (LDV) method was used to measure the velocity field for the liquid-only system and an optical technique was employed to determine the axial solids concentration distribution. Regions of in homogeneity in the tank were identified. It was found that, for a given solids loading, the solids concentration distribution depended on both particle size and particle size distribution. High solids loadings were investigated and a difference in the concentration distribution pattern was obtained with nickel, flint glass and sand particles. The CFD simulation results highlighted problems that could be associated with some conventional experimental methods of determining solids concentration distribution in a stirred tank. Keywords: nickel; CFD; simulation; mixing; solids suspension.

Jan 1, 2010

By Yuqing Feng, Ping Fang, Dezhou Wei, Zichuan Guan, Shuling Gao, Baoyu Cui

"Air core, as an inevitable phenomenon of rotational flow field inside hydrocyclone, has a major impact on hydrocyclone’s separation characteristics. The air core formation process, its structural features and its influencing factors were investigated systematically in this paper. Numerical study on the air core of a f50 mm hydrocyclone was carried out using a commercial computational fluid dynamics (CFD) modelling software package ANSYS/Fluent with the selection of Reynolds Stress Model (RSM) for the turbulence model and Volume of Fluid (VOF) for capturing the air-liquid interface. Physical experiments were conducted under the same geometry and operation conditions and the particle image velocimetry (PIV) technology was used to measure the liquid flow. The CFD results are compared with the experimental measurement with good agreement demonstrated. The air core was formed generally at around 0.7 s, air outsides of the air core was eliminated completely after 1.7 s. The centreline axial velocities and the water blocking hydrocyclone simulations showed that air of the air core entered the hydrocyclone from the underflow area and exhausted from the overflow area, but there was also some air flowing into the air core from the overflow area. The largest swing of the air core occurred on the top of the underflow pipe, and a protuberance similar to the shape of Adam’s apple appeared when the air entered the vortex finder. CFD simulation results also showed that the air core diameter was decided by axial position and parameters of hydrocyclone. The structural parameters have more effects on the air core diameter than the operational parameters."

Jan 1, 2014

By Jesus Gonzalez

Peirce-Smith converters (PSC) are long cylindrical chemical reactors where copper matte reacts with air. Traditionally, air is injected laterally through submerged tuyeres at different subsonic velocities. However, the accumulated experience in steelmaking, particularly in the basic oxygen furnace (BOF) process, shows that the gas bottom injection may improve the mixing efficiency of a liquid pool. The aim of bottom injection is to maximize the copper kinetic energy meanwhile avoiding excessive reactor splashing. When gas is injected from a single orifice at the bottom of a liquid pool, two main gas flow regimes can be identified, bubbling and jetting. In this work, a PSC with bottom air injection is studied by means of multiphase 3D CFD numerical simulations considering multiple air velocities. Special attention is paid to the bubbles and jet formation mechanism. The mixing efficiency is measured by calculating the turbulent kinetic energy of the copper matte.

Jan 1, 2008

By Graeme Lane, Yuqing Feng, Kaixi Jiang, Junwu Zhou, Zhengchang Shen, Tao Song

"Airlift reactors are used in a wide range of industries, such as hydrometallurgy, biochemical processes, chemical process industry and waste water treatment. With increasing superficial gas velocity three circulation regimes can be defined: no gas bubbles in the downcomer, gas bubbles remain stationary in the downcomer and gas bubbles in the circulating regime. In this paper a numerical modelling method was presented for gas-liquid flow in airlift reactors under different bubble recirculation regimes. Gas-liquid flow was modelled using the Eulerian two-fluid equations, and extra user defined subroutines were incorporated to consider the complex physics, such as bubble induced turbulence and turbulent dispersion force. Some alternative correlations of drag coefficient were used to capture the bubble distributions in the riser and downcomer of the airlift reactor, with consideration of the interaction between bubbles. As using constant bubble size is hard to describe global characteristics of the flow in airlift reactors for three circulation regimes, a model including multiple bubble sizes was presented to get more accurate simulation results of gas holdup and water velocity. Also, the inhomogeneous MUSIG model was used, as an exploration, to compute the complex flow regimes. The MUSIG model is a population balance method to predict the bubble sizes, taking break-up and coalescence into account. The modelling method was applied to a laboratory internal loop airlift reactor, and the simulation results were compared with the experimental measurements of Simcik et al. (2011) for gas holdup and water velocity. The three bubble circulation regimes were observed from simulation figures. Reasonable agreement was obtained over a range of operating conditions, and an improvement was demonstrated using the proposed methods."

Jan 1, 2014

By Oszkar Bíró, Peter Mittag, Bernhard Aigner, Stefan Pirkeri, Philipp Glttleri

"This paper considers CFD modeling of processing and melting of small metallic particles by means of electrical heating. The particles fall continuously onto a liquid metal bath which is heated by an electric arc. After melting of the particles liquid metal is tapped.The charging behavior of the particles in the supply unit as well as in the reactor freeboard is studied by means of Euler-Euler granular simulations. The flow situation in the metal bath due to gas injection and magnetic fields is calculated by combining Navier-Stokes and Maxwell solvers. The macroscopic melting process occurring during the continuous charging of the particles is studied by kinetic laws for melting reactions. The temperature field is evaluated by balancing heat sources due to Joule's heating and heat losses due to latent heat of melting as well as convection and radiation.As a result of these simulations the fully three-dimensional flow fields of particles and gas in the atmosphere is obtained. Furthermore the flow field as well as the magnetic field in the metal bath can be studied. As a main result the three-dimensional concentration field of the still unsolved particles in the liquid metal can be evaluated.IntroductionIn this paper a metallurgical process for melting small iron particles is investigated by means of CFD-simulations. In a first process step the preheated and already mostly reduced iron particles fall continuously into the reactor through a vertical tube mounted at the lid of the furnace. In a second step the particles are melted in the liquid metal bath. Meanwhile the melt is heated by electric arcs to provide the heat needed for melting the incoming particles.In Figure I a sketch of the geometry setup considered in this paper is provided. The particles fall from a height of approximately 12 meters onto the bath surface in the center of the furnace. The melt is heated by electric arcs that are placed symmetrically around the center axis."

Jan 1, 2001

By Andrei Kolesnikov

"The results of CFD modelling of complex heat transfer and fluid flow in DC electric arc furnace are presented. The model includes energy and momentum conservation equations for gas phase (freeboard), liquid phase (molten bath) and solid phase (walls and graphite electrode). Simulation revealed the temperature distribution in all elements of the furnace. The influence of the arc and melt resistance on the temperature distribution is discussed.Computational Fluid Dynamics (CFD) becomes a powerful tool for better understanding of the heta and mass transfer processes in industrial equipment. In current work commercial CFD software STAR-CD was used to analyse complex heat transfer in DC EAF, including forced and natural convection, conduction and radiation. The finite volume solver using body-fitted grid was used. The grids are non-staggered hexahedral grids and all variables are evaluated at cell centers. Steady-state calculations were performed with PISO algorithm. For the convective terms in equations the hybrid differencing scheme was used. An upwind differencing scheme was used for other terms. Conjugative heat transfer was taken into account in heat transfer calculations on gas-liquid and gas-solid interfaces. Radiation from the arc and inner surfaces of the furnace was calculated by a “discrete beam” approach. Source term in energy conservation equation was used to represent heat release in the zone occupied by the DC electric arc. The volume of this zone was calculated from semi-empirical equations. The Joule heat and heat due to chemical reaction in the melt were calculated through the source term expression in the energy balance equations for the slag and metal in the molten bath. The gas acceleration in the electric arc zone and stirring of the slag due to gas bubbles motion were accounted for by momentum source terms in the momentum conservation equations.All simulations were carried out using axisymmetric two-dimensional grid. The grid resolution was 20 cm per cell in most part of the furnace, except the arc region, where resolution was 5 mm per cell."

Jan 1, 2003

By B. Wan

The process of synthesizing silica nanopowder by the gas phase thermal oxidation of tetraethylorthosilicate (TEOS) in a diffusion flame reactor was simulated using a commercial computational fluid dynamic (CFD) code. The fuel combustion process and silica particle formation and growth in the flame are modeled. The temperature, velocity and particle size distribution (PSD) fields inside the reactor are computed. Chemical reaction rate and population balance model (PBM) were used to calculate the particle formation and growth and PSD.

Jan 1, 2006

By Swati Mohanty

Spiral concentrator, which works on the principle of gravity separation, is widely used in the mineral industries for separation of solid particles based on the size and density. In the present study, CFD simulation of the flow profile on a laboratory scale spiral concentrator of height 2.29 m, pitch 0.435 m and width 0.17m was carried out. A single computational domain covering the entire spiral was considered. Volume-of-Fluid algorithm was used to track the interface and predict the flow profile. In the present study the Donor-Acceptor method was used for the calculation of volume fraction in the cell. Simulation results show that the velocity increases as the water flows down due to the convective force and gravitational force, but the viscous force acting in the opposite direction retards the flow. Thus the maximum velocity is attained at turn three at which the viscous force overpowers the other forces and further down the velocity decreases. Measured water depths along the radius of the spiral at various heights were compared with simulated depths, which were in good agreement. The simulated liquid phase velocity profile was used to carry out preliminary studies on the flow profile of solid particles of size 500 micron and densities 2600kg/m3 and 5100kg/m3, based on the Lagrangian model taking into consideration fluid-particle interaction.

Sep 1, 2012

By Oluwafemi Salami

Conducting full-scale fire experiments in underground mines is difficult and expensive. Sometimes, it is impossible to carry out fire tests in mines due to the danger they pose to mine workers and the underground facilities. However, understanding the potential occurrence of fires in a mine and their dynamics is critical to developing an appropriate emergency evacuation plan, fire suppression, and firefighting techniques. In this study, a numerical fire model using computational fluid dynamics (CFD) with fire dynamics simulator (FDS) 6.7.6 was developed to examine the influence of tunnel bifurcation on maximum temperature beneath a tunnel ceiling. The model was validated using full-scale experimental data. The maximum ceiling gas temperature predicted using FDS and measured maximum smoke temperature during the experiment shows good agreement. The study shows that bifurcation could lower the maximum ceiling temperature of fire-induced smoke in a subsurface environment. The findings from this study apply to the broader field of ventilation and fire life safety in tunnels and other underground environments.

Jun 23, 2024

By X. S. Fan, T. R. Feng, D. Chen

To meet the flotation process requirements, the large-scale flotation cell must be designed for sufficient pulp mixing, reasonable air dispersion and solid suspension, which is very difficult to design for super large-scale flotation cell. Hence, the flotation cell scale-up and optimization for a short period is a challenge for super large-scale flotation cell design. Computational fluid dynamics simulation has been applied in flotation cell study, which is effective, timesaving and low cost in the study of the flow field distribution for super large-scale flotation cell. This paper demonstrates the flow field characteristics in the 600m3 flotation cell by CFD. The flotation cell model is meshed with hexahedral grid and the high mesh quality assures the computational precision. MRF model, namely a steady-state approximation, ensures that individual cell zones can be assigned different rotational speed. Euler-Euler approach is chosen for multiphase flow simulation to treat different phases as interpenetrating continua. Numerical simulation of single phase, two-phase and three-phase are investigated respectively. Results suggested that the flotation cell has the apparent upper and lower flow circulation. In condition of the impeller rotation, the air is dispersed reasonably. In addition, the solid phase volume concentration distribution in the cell has a little difference.

Jan 1, 2015

By A. Juganda, Jr. Bogin, S. Saki, R. Gilmore, J. Brune, J. Grubb, S. Lolon

"In longwall mines, atmospheric or barometric pressure fluctuations can disturb the pressure balance between the gob and the ventilated working area of the mine, resulting in a phenomenon known as “gob breathing”. Gob breathing triggers a gas flow across the gob and the working areas and may result in a condition where a methane accumulation in the gob flows into the face area forming an explosive mixtures. This paper discusses results of Computational Fluid Dynamics (CFD) modeling carried out to analyze this phenomenon and its impact on the explosive mixture development under a bleeder-ventilated longwall gob panel scheme. Modeling results indicate that the gas inflow and outflow across the gob and the formation of Explosive Gas Zones (EGZs) are directly affected by the barometric pressure changes. Methane gas and EGZs in the gob expand out toward the face and bleeder entries during the falling barometric pressure. Where methane zones interface with mine air, EGZ fringes may form along the face and in the bleeder entries. When the atmospheric pressure increases, an ingression of oxygen into the gob is observed that can also increase EGZs in volume. The findings from this study help assess the methane ignition and explosion risks associated with fluctuating atmospheric pressures. INTRODUCTION Methane explosions continue to be a daunting risk for underground coal miners, although the number of related fatalities and injuries in the U.S. coal mining industry has steadily decreased since the establishment of the U.S. Bureau of Mines in 1910. Still, the consequences of a methane explosion are often disastrous with multiple fatalities and property damage often leading to permanent shutdown of the mine. Methane is formed during the coalification process and is released from the coal seam to the mine atmosphere when the coalbed is disturbed by mining or natural causes such as earthquakes. If not properly diluted by appropriate mine ventilation, this methane may accumulate in the active mine workings and gob areas. Researchers at the Colorado School of Mines have developed numerical models showing how and where explosive methane may accumulate in gob areas. Historical mine explosions appear to show a connection between mine explosions and fluctuating barometric pressure as a result of stormy weather. McIntosh (1957), Boyer (1964), and Kissell et al. (1973) studied the influence of barometric drops on major coal mine disasters in the United States prior to 1970. Their statistical analyses found that a majority of these disasters occurred in the fall and winter months when the barometric pressures were influenced by unstable weather conditions and noted increased methane content in the mine workings during times of falling pressure. Studies conducted by Fauconnier (1992) and Hemp (1994) found similar connections between methane explosions and barometric pressure fluctuations in a majority of gas explosions in South African mines. Ten out of twelve major mine explosions with five or more fatalities in the U.S. after 1970 were found to have occurred during the months of November through April when barometric pressure swings were more abrupt and intense (Lolon et al., 2015)."

Jan 1, 2016

By J. A. Marts, R. C. Gilmore, Jr. Bogin, S. A. Saki, J. F. Brune, J. W. Grubb

The main purpose of mine ventilation design is to provide sufficient quantity and quality of air to the workers and to dilute methane and other contaminants. It is generally perceived that the additional air along the longwall face will improve methane dilution on the face and in the tailgate. However, computational fluid dynamics (CFD) modeling efforts at the Colorado School of Mines (CSM) under a National Institute for Occupational Safety and Health (NIOSH) funded research project have found that higher flow velocities along the longwall face are shown to increase the pressure differential between the gob and longwall face allowing more methane to be entrained in the active face and tailgate (TG) area, thereby negating the dilution effect. The increased face ventilation causes an increase in pressure along the headgate side which allows more oxygen to ingress into the gob area, thereby increasing the amount of oxygen available to form explosive methane air mixtures in the gob and to support spontaneous combustion of the coal. Two things are happening at two different locations along the longwall face, increase in pressure as the flow enters the active area cause an increase in pressure which forces more oxygen into the gob, but as the flow is developed along the longwall face, the higher velocities causes greater pressure drop allowing methane to enter the active face further downstream from the headgate (HG) towards the tailgate side. In this paper, a parametric study will be presented to discuss the effect that varying the face air quantity has on methane concentrations in the tailgate and formation of explosive gas zones (EGZs) in the gob. Counter to conventional wisdom, it appears that increased longwall face air quantities may increase the explosion hazard as they result in increased EGZ volumes in the gob, along with increased methane quantities in the tailgate return.

Jan 1, 2015

By S. Lolon

Spontaneous combustion in underground coal mines has become a serious problem particularly in the caved area (gob). Recent statistic has shown that approximately 17% of total 87 underground coal mine fires in the United States are attributed to spontaneous combustion (De Rosa, 2004). Spontaneous combustion results from a self-heating process in exothermic conditions. The accumulated heat, if not removed, is conducive to the rapid increase of temperature and may result in mine fires or explosions. It is well accepted that the interaction between oxygen and coal substances is the main cause for spontaneous combustion, while other factors such as pyrite, moisture, and bacteria play a secondary role to the self-heating of coal. Therefore, only coal oxidation is considered in this study. Coal oxidation occurs as coal comes into contact with air. This process involves complex phenomena in terms of heat transfer, chemical surface absorption, and energy balance related to inherent properties of coal (Wang et al., 2003). For simulation purpose, overall reactions can be simplified as suggested by Mitchell (1990): C + O2 CO2 + heat [65°- 94°C] (1) CO2 + C 2CO + heat [100°- 150°C] (2) Stoichiometric reactions (1) and (2) show that 2.66 grams of oxygen are required to oxidize 1.0 gram of carbon. Once the temperature exceeds 100°C, the reaction has little chance to stop. The crucial step in reducing spontaneous combustion risk is locating the ignition point of spontaneous combustion (hot spot). This information is useful in the effort of developing a preventive method effectively. In this study, this was conducted based on the best gathered information: ventilation surveys conducted in an existing longwall mine located in the western US; the laboratory experiments performed on a physical gob model; and the CFD models to investigate the flow behavior in the gob, the oxidation of coal, and heat transfer phenomena. The hot spot location is determined as a function of oxygen concentration and gob temperature. The critical values are the following: 5% (by volume) for oxygen and 100°C for gob temperature. CHARACTERISTICS OF GOB MATERIAL AND MODELS Permeability Concept Permeability is one of the key parameters in any study dealing with porous medium. This is determined by particle size and the ratio of void volume to the total volume (porosity) of porous medium. This relationship is known as the Carman-Kozeny equation (Scheidegger, 1957) and is given mathematically as [ ] where k* = theoretical specific permeability, m2 dm = the mean particle size, m n = porosity However, permeability tests with various particle sizes were carried out at the University of Utah?s ventilation model (Figure 1) to verify this relationship. These tests indicated a necessity to modify this theoretical relationship. A factor of 0.898 was obtained to make up the different between theoretical permeability and the experimental one. The modified version of (3) is given by: [ ] This equation is used to determine permeability of the simulated mine gob throughout the study. [ ] Mine Gob Material The gob is represented by zones filled with material of given size distribution. For physical measurements and CFD simulations, the material used were crushed rock and coal. Due to compaction, three different characteristics of gob material were used in the model: unconsolidated, semiconsolidated, and consolidated. In the real condition, the largest coal-rock particles are more likely to be located in the area behind the shields. This material is freshly broken and unconsolidated. The size of these broken particles is based on a study carried out by analyzing images taken from the area behind the shields in three coal mines in the United States (Pappas and Mark, 1993). The results have shown that the mean size of materials behind the shields is about 1.22 m. The permeability associated with this size was then obtained from (4) and permeability tests. This was 4.68 x10-7 m2. For semiconsolidated and consolidated zones, due to the lack of experimental data, these were determined through CFD simulations. These simulations were designed by assigning a permeability of 3.15 x 10-8 m2 for consolidated zone and 7.98 x 10-9 m2 to the other zone. The airflow pattern in the gob assigned with given permeabilities was expected to follow the experimental distribution (Brunner, 1985). The mean particle sizes for the semi-consolidated and consolidated zones were 0.02 and 0.006 m, respectively; smaller than those of the unconsolidated zone. Figure 2 shows the 3D view of permeability changes in the gob.

Jan 1, 2009

By S. Ragab

The objective of flotation models is to predict the recovery rate of minerals from a flotation cell. In this paper a CFD?based flotation models is presented. The model uses first order rate equation, and the collision frequency is estimated using Abrahamson [1] model with input from CFD two-phase simulations. Spatial distributions of dissipation rate and air volume fraction (also called void fraction) are determined by the two-phase hydrodynamic simulations. Knowing these parameters throughout the machine is essential in understanding the effectiveness of different components (rotor, stator or disperser, jets) on the flotation efficiency. Not only the average pulp recovery rate can be estimated but also the regions of high/low recovery rate can be identified. The CFD-based flotation model presented here can be used to determine the dependence of recovery rate constant at any locality within the pulp on particle diameter, particle specific gravity, contact angle, and surface tension.

Jan 1, 2012

By Tom Liu

This paper is a fundamental study of bubble-particle collision efficiency in turbulent flows. A CFD-based microscale scheme for the prediction of turbulent bubble-particle collision efficiency in a 3D isotropic turbulent flow is developed, and in this way, turbulence effects on particle-bubble collision efficiency are taken into account for the first time. By comparing the computed results with experimental data from the literature, the CFD simulation is shown to correctly predict the bubble-particle collision efficiency in a quiescent environment. The CFD method was then applied to the investigation of bubble-particle collision efficiency in a turbulent flow, and the simulation results are shown to be qualitatively consistent with the real observations. The predicted collision efficiency can be used in a multiscale modelling scheme, which would account for all control variables in a flotation cell, such as impeller structure and tip speed, air nozzle structure and air flow rate This study indicates how microscale simulations of bubble and particle behaviours, coupled to macroscale CFD simulations of flotation cells, can provide a powerful tool predicting the flotation rate in a real flotation cell.

Jan 1, 2014

By Yiming Rong, Gang Wang

"Water quenching is one of common heat treatment processes but with complex physical phenomenon, in which vaporizing, bubbling and convection are involved. The general simulation of liquid quenching processes is conducted by using temperature dependent value of heat transfer coefficient (HTC) obtained from experimental measurement under a strong assumption. A preliminary multiphase flow model developed by using Ansys-CFX for water quenching has been presented in this paper, which shows a distributed HTC with a time dependent feature. The conjugate heat transfer is used to describe thermal conduction in the simplified workpiece, and generation, growth and flow of bubbles are modeled by using Eulerian-Eulerian method for the gas-liquid multiphase flow. The simulation result shows the intense variation of heat exchange through the interface between workpiece and liquid, and the unsteady distribution around the solid surface.IntroductionQuenching is a rapid cooling process which makes the desired solid phase transformations happen for the certain material properties. For aluminum alloy, rapid quenching will prevent all phase transformation and generate a supersaturated single phase."

Jan 1, 2012

By J. P. Jylhä, A. Jokilaakso

The Outotec Flash Smelting process is widely used in primary copper production. A layer of slag on top of matte (copper–iron sulphide melt) creates a barrier through which matte droplets descend. Droplets of different sizes settle through the layer at different rates, which increases possibility of copper losses as smaller matte droplets may be removed with the slag. In this work CFD-DEM (computational fluid dynamics – discrete element method) coupling was used to simulate matte droplets settling through the slag during slag tapping. Discrete droplets and slag flow were simulated with DEM and CFD, respectively. This new approach enables tracking and solving movement and interactions of individual droplets. User-defined models were created to take coalescence and reaction kinetics into account. Coalescence created larger droplets, while the reactions decreased droplet size and increased its density. Drag caused funnelling effect, which pulled droplets towards the centre line of the descending cluster. This research is the only study known by the present authors in modelling liquid-liquid systems with the CFD-DEM coupling. Simulating coalescence and chemical reactions for discrete droplets enables more realistic simulations. However, CFD- DEM simulations are computationally resource intensive already without additional models, so a full scale industrial settler modelling still waits for further increase in computer performance. INTRODUCTION Flash smelting is a method used in copper and nickel production. In the flash smelting process, oxidation of sulphides in a concentrate is used to release energy needed to smelt the feed. In the process, copper matte is formed. Matte has copper content in the range of 55 – 70% (Sridhar, Toguri, & Simeonov, 1997), the rest is mainly sulphur and iron but has still many other impurities (Taskinen, Seppälä, Laulumaa, & Poijärvi, 2001). Oxidized iron forms iron-silicate slag while SO2 is collected for acid production (Davenport & Partelpoeg, 2015). The Outotec flash smelting furnace (FSF) can be divided in four major component: concentrate burner, reaction shaft, settler and uptake. Concentrate, flux and oxygen enriched air are injected into the reactor through the concentrate burner. The mixture then reacts in the reaction shaft forming matte and slag, which are collected in the settler, while off-gases flow out of the reactor through uptake for heat and dust recovery, and SO2 collection for H2SO4 production. The matte and slag form separate layers due to immiscibility and density difference. Raw material feed is continuous, but the matte and the slag are tapped periodically. A flash smelting furnace is depicted in Figure 1.

Jan 1, 2019

By Sanja Miskovic, Manjunath Basavarajappa

In this study, a coupled CFD-DEM (Computational Fluid Dynamics-Discrete Element Method) approach is used to simulate two-phase, solid-liquid flow within the flotation cell. Turbulence plays a key role in flotation cells and level of turbulence is affected locally by the presence of particles. Turbulence in the flow is modelled using Reynolds Averaged Navier-Stokes (RANS) based k-e model. Using particle transport and spatial distribution, mixing mechanism for differed sized particles are analysed. By including the effects of particle-particle collision, particle rotation, and particle translation, as well as the volume occupied by the particles in the fluid, CFD-DEM approach can be used to accurately predict the motion of particles in the turbulent solid-liquid flows with strong coupling. Electrical resistance tomography (ERT) is used to investigate the effect of impeller size, shape, off-bottom clearance and impeller speed on mixing and suspension behaviour inside flotation cell. The effect of flotation impeller design, impeller size, and off-bottom clearance on flow hydrodynamics is investigated. In the case of moderate and high solids loading (> 10 %w/w), energy is dissipated through inter-particle interactions and four-way coupling becomes important (Derksen 2003). The main goal of this study is to understand the effect of particle size and concentration on flow dynamics at high Re in highly turbulent systems using advanced modelling and measurement techniques. Experimental measurements obtained by ERT) system allowed detailed mapping of temporal hold-up profiles of the dispersed phase. Impeller size and off-bottom clearance are found have significant effect on the particle suspension characteristics.

Jan 1, 2014