Search Documents

By Pengfu Tan

In order to understand and control the formation of dioxins and furans in iron ore sinter plants and to possibly avoid expensive end-of-pipe removal stages, thermodynamic calculations on the stability of PCDD/Fs have been combined with sinter process modeling by means of kinetic simulations and computational fluid dynamics. Three thermodynamic databases of PCDD/Fs, derived using the Group Additivity approach and two computational molecular modeling methods, MNDO and PM3, respectively, combined with the SGTE database have been used to model the PCDDIF formation and the PCDD/F isomer distributions in iron ore sintering. The results show very similar conditions for the PCDD/F formation using the different databases. Concerning the isomer distributions of toxic PCDD/Fs, the calculated data derived from the MNDO method showed optimum agreement with measured values from industrial furnaces. The thermodynamic calculations also proved that the isomer distributions found in combustion processes are equilibrium values, while the concentration ratio of PCDDs and PCDFs in the wind boxes does not represent an equilibrium state. A CFD model of iron ore sintering was developed including the kinetics of 15 relevant chemical reactions with the aim of describing the temperature fields and the gas flow patterns in the sintering bed. On the basis of these information and experimental observations reported in the literature possible pathways for the PCDD/F formation in sinter plants have been identified: PCDD/Fs are formed in the critical temperature region of 250 to 450°C below the combustion zone, are then carried downwards with the gas and condensed close to the bottom of the sinter bed. Transported to the discharge end with the solid mixture, they are again released into the gas phase when the flame front approaches the bottom. PCDD/Fs are also formed in the last wind boxes when the hot off-gases cool down and reach the critical temperature range. Consequently, methods to reduce the dioxin formation reactions are proposed including inhibitor addition to the sinter feed and off-gas quenching.

Jan 1, 2003

By Michael L. Free, Shukla. Abhijeet

The short-circuiting tendency of a copper electrowinning cell (EW) based on the roughness of cathode deposits was investigated with varying parameters such as Cu concentration, guar concentration, current density, and temperature. A bench scale EW cell was set up to produce copper under conditions similar to industrial Cu EW. In addition, a statistical analysis and finite element analysis were performed to analyze the deposited thickness, surface roughness and maximum feature height. Interaction data for the cathode deposits suggest that current density and Cu concentration together affect the roughness of deposit more than the combined effect of guar and temperature. Model predictions of roughness were similar to experimental values. The predicted time to short-circuiting was similar to the measured time. This study provides an industrial relevant baseline against which new operating parameters can be compared.

Jan 1, 2014

By MOHAMED WAEL, ANASTASIA GOEVA, ABDULAZIZ TAWFEEQI, IHTISHAM TARIQ, S. M. ZAKIR HOSSAIN, NIAZ ALI KHAN, MUHAMMAD FAISAL IRFAN

Carbon capture and storage (CCS) is an attractive area of research in such fields as CO2 mineral carbonation, global warming and sustainable energy systems. In this study, carbonation efficiency for aqueous mineral carbonation (MC) was achieved through two steps,which include leaching of calciumfrom cement kiln dust (CKD) followed by the reaction of pure CO2 with the calcium hydroxide precipitates formed by the hydroxylation using NaOH. Response surface methodology (RSM) with a Box-Behnken design (BBD) was applied to optimize the calciumleaching yield,while the carbonation efficiency from CKD was assessed using RSMwith central composite design (CCD). Optimization of calcium leaching is highly important, as it is a rate-limiting reaction step in MC and also influences and enhances carbonation efficiency. Different parameters including acid concentration (HNO3), leaching temperature, leaching time, and dose of CKD sample were considered in order to optimize the maximum yield of Ca leaching from the CKD sample. In addition, different CO2 flow rates and temperatures were used as parameters for optimizing carbonation efficiency. Two quadratic regression models were developed for each process, i.e. calciumleaching and carbonation. For calcium leaching, a maximum of 98.55% calcium was extracted under the optimal set of acid concentration 4.13 M, 90 °C, 28 min leaching time, and 13.8 g of CKD sample. For carbonation, the maximum carbonation efficiency of 89.2% was achieved for a CO2 flow rate of 1163 cm3/min at 90 °C. Calcium leaching results indicate that the leaching yield was significantly affected by all the input parameters except leaching time. For carbonation, both factors affected the carbonation efficiency, with the effect temperature shown to be greater than that of the CO2 flow rate. Additionally, the predicted results agreed well with the experimental values for both calcium leaching and carbonation processes, with errors of less than 1% and 5%, respectively.

Jun 5, 2020

By MOHAMED WAEL, ANASTASIA GOEVA, ABDULAZIZ TAWFEEQI, IHTISHAM TARIQ, S. M. ZAKIR HOSSAIN, NIAZ ALI KHAN, MUHAMMAD FAISAL IRFAN

Carbon dioxide (CO2) mineral carbonation is a carbon capture and storage (CCS) technique that can sequester CO2 permanently in the form of durable carbonate minerals. Extendable permanent storage capacity — onshore and offshore storage reservoirs — and the exothermic nature of the mineralization process make this technique attractive [1]. In the present study, maximum carbonation efficiency for aqueous mineral carbonation was achieved through two steps: optimization of calcium leaching from cement kiln dust followed by the reaction of pure CO2 with the calcium hydroxide precipitates formed by hydroxylation using sodium hydroxide (NaOH). Optimization of calcium leaching is highly important, as it is a rate-limiting reaction step in mineral carbonation and also influences and enhances the carbonation efficiency. To maximize the calcium leaching as well as carbonation efficiency, response surface methodology — with either Box-Behnken design or central composite design — was applied. It is the most well-known statistical, fully randomized and bias-free method for determining optimized conditions using quadratic regression models.

By H. Chaouki

The aim of this work is to simulate the forming process of a green anode. For this purpose, a nonlinear compressible viscoelastic constitutive law is presented. The concept of natural reference configuration is considered. Within an isothermal thermodynamic framework, a Helmholtz free energy is proposed to take into account the nonlinear compressible deformation process occurring between natural reference configuration and current configuration. A dissipation potential is introduced in order to characterize the irreversible aspect of compaction process. Thus, the constitutive law is formulated through two equations: 1) an expression of Cauchy stress tensor and 2) a differential equation characterizing the evolution of the natural reference configuration. Material parameters are assumed to be a function of the bulk density. This allows capturing the micro-mechanical phenomena within the chosen macroscopic approach. An experimental study based on compaction tests is carried out in order to characterize the mechanical behaviour of the anode paste. A user's material VUMAT subroutine for finite element dynamic explicit analysis has been developed and implemented in the ABAQUS commercial software. To evaluate the model predictive capability, numerical simulations of the compaction forming process of anode paste were performed. Simulation results show that the constitutive law predicts the experimental trends and gives interesting physical responses. This constitutes a first step toward making a benchmark with experimental results on the forming process of anode paste.

Jan 1, 2011

By D. Minson, G. Lindner, S. Miskovic

Thickening is an essential process in the mining industry and currently one of the most costeffective technologies for recovering water from tailings streams. Improved performance of thickeners can lead to significant energy savings, and increased recovery of water close to the mill. Its performance is directly affected by the rheological properties of slurry at the raking zone. To further understand the dynamics of slurry inside thickeners, a multiphase particle-in-cell method (MP-PIC), is used to investigate its behaviour on a sliding channel. Parameters like particle-size distribution (PSD), rheological properties, and domain discretisation were benchmarked to study its influence on highly dense flows. The slurry concentration by weight of the samples ranged between 60 to 70% of solids by mass, under which conditions it is known to behave like a non-newtonian fluid. An experimental flume, 0.2 m wide and 2 m long, was used for validation under different angles and slurry composition. Visual analysis of the slurry height and shape, as well as velocity measurements, were compared with numerical data. The preliminary results show the ability of MP-PIC method to represent different types of highly dense flows, while still being computationally affordable and practical for industrial use. Keywords: High density thickener, CFD, MP-PIC, multiphase-flow

Jan 1, 2020

MODELLING FOR DESIGN AND CONTROL AT KALGOORLIE NICKEL SMELTER Modelling at KNS has been driven by a need to have better tools to answer questions about control, design and optimisation of the plant. This paper describes a number of modelling tools developed for use at KNS and how the development of modelling tools fits with a design and control strategy. In the area of control, a thermodynamic model has been developed to allow process parameters to be calculated. The model allows for the automatic selection of data from a laboratory data base, desired operating set points are input and the operating parameters are calculated. The model has the ability to adjust the burner efficiency and the furnace heat loss depending upon conditions. To assist with the testing of proposed furnace design changes, various mathematical models of the burner and shaft have been created using PHOENICS. These models are particularly well suited to combustion heat transfer and flow aspects of both the burners and shaft. Validation of these models have been based on novel sampling and measuring techniques in the shaft and also by using scale and cold models. Overall plant energy and cost structures can be analysed using a flowsheet simulation, developed to evaluate potential alterations to the plant flowsheet. This model system allows the global evaluation of process changes.

Jan 1, 1993

By Raghavan Srinivasan, Y. V. R. K. Prasad, William G. Frazier, Steve C. Medeiros

"Aerospace gas turbine disks operate in an environment of relatively high stresses caused by centrifugal forces and elevated temperatures. These severe conditions necessitate the need for materials with good high temperature strength and low cycle fatigue resistance. One class of alloys used for this task is the nickel base superalloys, out of which, lN 718 is the most widely used in the aerospace industry [I]. The properties of IN 718 are attributed to the combined effects of the chemistry, heat treatment, and microstructure. The chemistry is tailored not only for solid solution strengthening but also for precipitation hardening developed during heat treatment, which combined with a fine grained microstructure lead to excellent mechanical properties such as low cycle fatigue resistance and elevated temperature strength. The properties of a gas turbine disk are sensitive to the microstructure, in particular the grain size, which is dependent on the processing history. The ability to precisely control the microstructural development during forging is dependent on controlling the process so that the workpiece is deformed within a ""safe"" region where no microstructural damage or flow instabilities occur. The microstructural mechanisms during deformation may themselves vary within the ""safe"" region and it is desirable to determine them within the range of parameters that are commonly used in industrial processing. The objective of this work is to establish a relationship between the grain size and the process control parameters i.e., temperature and strain rate, in the hot working of IN 178.BackgroundIN 718 belongs to a class of nickel-iron base superalloys with a face-centered cubic (FCC) austenitic matrix strengthened by precipitation of ordered intermetallic or carbide precipitates and is a nickel-rich alloy strengthened by ordered body centered tetragonal (BCT) precipitates, Y"". The alloying of nickel and iron results in the formation of the austenitic matrix, y, while the addition of chromium and molybdenum cause solid solution strengthening. Also, alloying with niobium causes precipitation of the metastable hardening constituent Y"" (Ni3(Nb,Ti)). Extended aging treatments in the temperature range 650 °C - 980 °C allow the more stable orthorhombic 8 phase (Ni3Nb) to form. Trace amounts of boron are also added to serve as grain boundary strengtheners for high temperature creep resistance. Hot deformation of IN 718 is performed between the Y"" and 8 solvus temperatures to facilitate recrystallization and to control grain growth. While deformation below the Y"" solvus (""""962°C) produces no change in grain size, deformation above the Y"" solvus but below the 8 solvus will lead to a dynamically recrystallized microstructure. The presence of the 8 precipitates will pin the grain boundaries after recrystallization and restrict grain growth so that a fine grain microstructure is obtained. Deformation above the 8 solvus ("""" 1038°C) will lead to a microstructure which has undergone grain growth during deformation."

Jan 1, 2000

By N. El-Kaddah, V. L. Acoff, M. Arenas

"Control of the weld penetration, weld pool shape and temperature field in arc welding is critical to achieve the desired microstructure and mechanical properties in the resultant welds. This is particularly important in welding gamma titanium aluminides (y-TiAl) where the weld pool is relatively shallow and the temperature gradient in the heat-affected zone is quite steep. This paper investigates theoretically and experimentally the evolution of the shape of the weld pool for stationary gas tungsten arc welds in y-TiAl alloys. A mathematical model has been developed for predicting the electric current flow, the resultant electromagnetic forces in the metal, and the transient development of the temperature and velocity fields in the molten pool generated by a spatially variable heat flux at the arc-striking zone. The governing equations were solved numerically using the control volume technique. Welding experiments were carried out to determine the shape of the weld pool and the extent of the heat affected zone at different welding times for model validation. The numerical predictions were found to be in reasonable agreement with measurements. It is suggested that the model presented in this paper is likely to provide the information needed for predicting the microstructure and mechanical properties of y-TiAl welds. IntroductionIn recent years, there has been considerable interest in the welding of gamma titanium aluminides as a means of manufacturing high-temperature, high-performance components for turbine and automobile engines. These alloys are highly reactive, and must be welded in vacuum or inert gas atmosphere, rendering electron beam, and inert gas tungsten arc welding (GTA) to be the only suitable welding techniques for joining these materials [1-3]. The GTA welding is especially attractive because of its flexibility and relative cost-effectiveness. However, previous studies on the weldability and weld properties of these alloys using this process revealed a decrease in the mechanical properties in the fusion zone [4]. This is generally attributed to the change of the original duplex structure to a dendritic structure amplified by the shallowness of weld spot. Quantifying the effect of arc operating parameters on the geometry and temperature profile of the welding pool prior cooling is the first step for controlling the fusion zone microstructure and mechanical properties."

Jan 1, 1999

By Tsvetan Damyanov, Nikolay Hristov, Evgeni Kostadinov, Lubomir Kuzev

"The process of grinding in drum mills is directly dependent on the dynamic forces occuring in the free movement of the grinding bodies in the volume. The subject of discussion are the potential impacts between grinding bodies having different shape and equal weight in a drum mill that has a diameter of 1200 mm. The study utilised ball grinding bodies having a diameter of 120 mm and tetraballpebs having a conditional diameter of 120 mm. The study was carried out using the Solid Works Simulation software. The following cases were examined: • Impact between two identical grinding bodies; • Impact between a grinding body and the shell surface. The study provided numerical values for impact force, kinetic energy of impact, linear speed of impact, linear displacment of the center of gravity and duration of impact deformation. The impact force of the balls is approximately 179000 N. The impact force of the tetraballpebs with typical vertices, edges and spherical surfaces depending on the area of impact is 177000 to 185000 N. Upon impact between two spherical surfaces of tetraballpebs, the contact area is 263.55 mm2, which is more than twice the contact area upon impact between two balls, which is 114.39 mm2. This determines better conditions for the transmission of kinetic energy by reducing opportunities for overgrinding."

Jan 1, 2014

By K. D. Luxbacher

Understanding methane emissions in underground coal mines is critical for a safe and productive mine. In addition to reasonable estimation of initial coalbed reservoir parameters, it is also crucial that changes in effective stress due to mining and pore pressure reduction are taken into account due to their effects on porosity and permeability. Primary parameters for estimation of emissions or modeling of the mining environment for this purpose are porosity and permeability which can change dramatically as a result of stress redistribution associated with mining and gas desorption from a large coal volume. These parameters affect the emission rates and ventilation requirements, as well as water inflow into the working environment. Stopping leakage, on the other hand, is a secondary stress dependent factor in estimation of emissions, as convergence of the roof and floor strata, compromising the integrity of the stopping, may result in leakage, making prediction of ventilation requirements difficult. This paper aims to examine the effects of porosity and permeability changes of the coal seam on methane emissions in an underground continuous miner section. The models were developed and executed in a dynamic fashion to simulate an advancing section. Through this process, the changes of effective stress in coal, particularly their change paths, on porosity and permeability were incorporated into the models and methane emissions, concentrations, air requirements, water inflow and possible leakage from the stoppings were investigated using a conventional coalbed methane reservoir model.

Jan 1, 2009

By Vaughan R. Voller

"Over the last ten years an increasing number of heat and mass transfer models of solidification processes have been presented in the literature. A key component in these models is the coupling between thermal and solute fields. This is achieved by developing relationships based on the behavior at the microscopic scale in the casting (typically a secondary dendrite arm space). An important class of models in this respect are models that describe the local scale redistribution of solute alloy components-collectively referred to as microsegregation models. This paper explores this class of models.IntroductionModeling takes a key role in processing materials for properties. This is particularly time in solidification processes where current computer techniques allow for a detailed treatment of the critical macroscopic phenomena. Typically these macro scale models consist ·of three main coupled components.1. A model of flow in the melt and solid + liquid mushy region.2. A model of heat transfer.3. A model of solute transport.The focus of this paper is the solute transport component, in particular the manner in which process scale transport is coupled to transport ·at the local scale of the solid-liquid interface. This step requires a modeling of the redistribution of solutes at the scale of the secondary arm spaces in the dendritic mushy region-a process referred to as microsegregation. Basic microsegregation models use the limit assumptions of no mass diffusion in the solid (Gulliver-Scheil [l]) or complete diffusion in the solid (equilibrium lever rule) in a fixed and space. Recently, however, more sophisticated microsegregation models that account for finite mass diffusion and coarsening of the arm space have been developed [2-19]. In this paper governing equations for a general microsegregation model for a multi-component alloy in a coarsening and spacing _are introduced. From these general equations various models for a binary system are developed and tested"

Jan 1, 2000

By C. V. Deutsch, D. A. Silva

The geologic modeling of multiple domains is required for the prediction of mineral resources and reserves. In general, explicit methods have wide acceptance among modelers due to the ability to manually add deemed important trends and geological features into the models. As an alternative, implicit geologic modeling techniques are increasingly adopted to reduce professional time requirements. These methods aim for objective, repeatable and flexible models by replacing the manual time-consuming interpretation of sections with the interpolation of signed distance values. These techniques, however, have been developed for the binary case of two domains at a time. The current paper presents an algorithm that extends the binary case to multiple domains. The method is able to transfer large-scale geological features from data to geologic models. A case study illustrates the benefits of the algorithm and practical implementation considerations.

Jan 1, 2016

An appropriate scale-up methodology is very useful for ball mill design and optimization of the grinding circuit. To date, the scale-up is based on the selection and breakage function method using lab mill data. However, a major drawback of that approach is that the charge motion in a lab mill is vastly different from that in a large industrial mill. Therefore, in the methodology proposed here, first the energy of individual collision in the mill is computed employing the discrete element method (DEM). Then the breakage of particles as a function of collision energy is studied in a separate experiment. Together, this information constitutes what is known as collision energy based population balance method. The appropriateness of this method is shown for 38-cm ball mill.

Jan 1, 1998

By N. R. Green, N. J. Humphreys, D. McBride, T. N. Croft, M. Cross, D. M. Shevchenko

"Centrifugal casting offers a route to high quality products in difficult to cast high temperature low superheat alloys and thin section molds. Under centrifugal forces metal is forced into thin sections and can fill thicknesses of less than a millimeter. However, due to the high liquid metal velocities there is a high risk of surface turbulent flow and air entrainment within the liquid metal. The combination of interacting flow-thermal-solidification phenomena and associated defects is a challenging modeling task which the authors have previously described and validated. Capturing the metal-air interface, on what are inevitably complex three dimensional geometries, results in highly computationally expensive simulations and simulating a single cast can take weeks on a single processor. This contribution reports on modeling a complex centrifugal cast, gas entrainment, bubble transport and solidification, employing meshes of up to a million elements and investigates the scalability of the model on high performance clusters.Introduction The ultimate goal of any casting procedure is to produce a defect-free product. For difficult to cast alloys, such as Titanium, centrifugal casting has a number of advantage over conventional gravity methods [1]. In the casting of thin sections, such as turbine blades, the metal is forced into the thin sections under centrifugal forces. The centrifugal force can overcome the backpressure due to surface tension and the metal can fill mould sections with thicknesses substantially less than a millimeter. However, the materials used in these castings [2] tend to be highly reactive and requires low superheat. High liquid metal velocity gives rise to a high risk of free surface turbulent flow and of associated entrainment through free surface turbulence [3] and trapping of gas present in the flow [1]."

Jan 1, 2012

By M. M. Botz

Cyanidation tailings which are disposed in a surface impoundment experience a loss of cyanide due to natural degradation, frequently reducing the cyanide concentration to very low levels. Quantifying cyanide losses in terms of pond geometry, local weather conditions and feed solution chemistry has been largely empirical though in many cases mining operations rely on surface impoundments to reduce cyanide to below an internally regulated concentration or an effluent limitation. To permit a more quantitative evaluation of cyanide losses in an impoundment, a computer simulation was developed to estimate the losses of free, WAD and total cyanide due to dissociation, photolysis and volatilization. Results of the model are compared against data collected for a North American tailings impoundment during 1998.

Jan 1, 1999

By M. M. Botz

Cyanidation tailings which are disposed in a surface impoundment experience a loss of cyanide due to natural degradation, frequently reducing the cyanide concentration to very low levels. Quantifying cyanide losses in terms of pond geometry, local weather conditions and feed solution chemistry has been largely empirical though in many cases mining operations rely on surface impoundments to reduce cyanide to below an internally regulated concentration or an effluent limitation. To permit a more quantitative evaluation of cyanide losses in an impoundment, a computer simulation was developed to estimate the losses of free, WAD and total cyanide due to dissociation, photolysis and volatilization. Results of the model are compared against data collected for a North American tailings impoundment during 1998.

Jan 1, 1999

By Troy Williams, Chris Preston, Ian Lipchak

"Underground blasting operations are challenging from the standpoint of the distribution of explosivesenergy representative of ring blasting. Energy from both shock and pressure regimes of commercialexplosives may appear concentrated in the collar region of a typical ring blast and diluted at the toes of holes due to the oblique geometries of blastholes. The non-homogenous nature of ore in whichexplosives are distributed via drillholes, adds to the complexities of generating particulate profiles from fragmented material with consistencies that are predictable from blast pattern to blast pattern - well suited for specific underground handling equipment and mill processing. In an ideal world, it would be the blasting operations themselves that represent the primary crushing mechanism, or at least mitigate mechanical crushing that can comprise a large component of the cost in generating suitable muck. This paper presents a dynamic break view in 3D that allows a planner to visualize the potential break zone around a blasthole generated by an explosive load using a Kleine field. Simple as it sounds, this methodology provides information that can be used in conjunction with cavity monitoring surveys (CMS) to potentially judge dilution due to overbreak as well as recovery for a typical blast. As examples, there are two break geometries that are examined regarding circular breaks and elliptical breaks around blastholes. Using a Kleine field to define break, a planner generated isosurface can be generated and compared to CMS data for calibration and prediction, using AEGIS 3D ring design software."

Jan 1, 2016

By M Reza Afshar

In this work a mathematical model was developed to study the effect of electromagnetic force on the magnesium melt containing nonmetallic inclusion in an electrical isolated channel. The electromagnetic force exerted on the inclusion induced by DC current field was calculated. A three dimensional turbulent flow equations were solved numerically to compute the velocity field within the channel. On the basis of computed velocity and magnetic fields, the trajectories of inclusions were calculated using the equations of motion for inclusions. Parametric studies were carried out to evaluate the effect of various parameters such as current density, channel diameter and particle density on the removal efficiency of inclusions from magnesium melts.

Jan 1, 2004

By Nagy El-Kaddah, T. T. Natarajan

This paper describes a novel approach for modeling three-dimensional electromagnetically driven flow in inductively stirred systems, and its application to modeling electromagnetic and flow phenomena in sub-mold rotary stirring of billets. It is based on a new T-\jf potential formulation of electromagnetic field, which does not require the calculation of the field in free space. A finite element algorithm has been developed for solving the electromagnetic field and Navier-Stokes equations in the fluid flow domain. The computed results have shown that the electromagnetic force field generates a strong rotational flow within the vertical section covered by the stirrer, and a relatively strong secondary flow beyond the stirrer. It has also been shown that the rotational and secondary flows were driven primarily by the vorticity of the force field at the billet corners. The induced flow in the molten pool was found to be turbulent and the effective mixing region in the molten pool was about three times the length of the stirrer. The mathematical approach presented in this paper provides the theoretical framework for quantifying these effects for any given stirrer configuration and operating parameters.

Jan 1, 1999