Search Documents

By Amar Prakash, V. M. S. R. Murthy, Ajoy Kumar Singh, Avinash Paul, John Loui Porathur

Weak, layered, and fragile rock mass formation, if not supported properly, is subject to roof failure thereby affecting safety and productivity in underground coal mines. Though Central Mining Research Institute – Indian School of Mines Rock Mass Rating (CMRI-ISM RMR) based well-defined support design guidelines are established, still occurrence of roof failures in underground coal mines is a real matter of concern for mining engineers and researchers. Numerical modeling techniques are successfully used by several researchers by simulating rock mass condition for stability assessment of mine openings. The analysis becomes crucial in case of weak and fragile rock formation. The present research envelops the determination of 31 cases of rock load by CMRI-ISM RMR under different geo-mining conditions followed by the development of modified Rock Mass Rating (RMR), i.e., rock mass rating dynamic ( RMRdyn) by incorporating P-wave velocity as a new parameter. The rock load determined using RMRdyn and numerical models was correlated and found in close agreement. In addition, the deviation in rock load determined by all the three approaches, i.e., CMRI-ISM RMR, numerical modeling, and RMRdyn, was compared with the actual field data. The percentage deviation obtained in RMRdyn and numerical modeling is less compared to CMRI-ISM RMR.

Sep 1, 2023

By J-P Franzidis, E V. Manlap

Froth transportation characteristics have a significant impact on the flotation recovery especially in large flotation cells. The probability of an attached particle being recovered into the concentrate depends strongly on the location where it enters the froth phase. A mathematical model for froth transportation in industrial flotation cells was proposed (Zheng et al, 2004). The model has been tested on the experimental data obtained from testwork conducted at the Western Metals Hellyer lead/zinc concentrator. The results show that the model fits the experimental data well, especially when the froth is well developed (ie for a deep froth). The model fitting can be significantly improved by incorporating an air dispersion factor. However, the models encounter some limitation, ie inconsistent prediction at different froth heights, due to change in froth viscosity which cannot be modelled at the present time.

Jan 1, 2005

By B J. Monaghan, P Zulli, X F. Dong

In the lower zone of blast furnace (BF), slag trickles down through a packed bed of carbonaceous particles in the form of films, rivulets or droplets. During its downward flow, there are significant interactions occurring between slag and other phases, ie gas and packing particles. In particular, the interaction between gas and slag is closely associated with the key flow phenomena such as loading, flooding or channelling in the lower zone of the BF. With gas introduced laterally through tuyeres and slag flowing downward from the cohesive zone to the hearth, this interaction varies throughout the lower zone and in the cohesive zone:  strong cross-flow of gas relative to the slag occurs in front of the raceway  counter-current flow of gas and slag occurs in the upper parts of the lower zone  both strong cross- and counter-current flow of gas and slag around and within the cohesive zone. These interactions can be linked to furnace irregularities, affecting smooth operation and limiting production. Therefore, understanding the gas-slag interaction and its influencing factors remains critical for BF process control. In the current review, the physical behaviour and key influencing factors of gas and slag flow are initially examined. Numerical models describing gas and slag flows in a packed bed at micro-, mesoand macro-scopic levels, are reviewed. Each simulation approach has its advantages and limitations, as well as unique assumptions when used in BF process modelling. Applications of these approaches to gas-slag flow in key identified functional zones will highlight the differences between approaches and modelling outcomes. As no single approach can fully address all aspects of gasslag interactions, an integration of different approaches is considered an effective and beneficial method in providing a more complete understanding of gas-slag flow in the BF process.

Jun 19, 2024

By Du Sichen

The present paper presents the modeling of gas-solid reactions in the case of the decomposition of CaCO3, describing the reaction kinetics on the basis of heat-transfer concepts which could be compared with thermogravimetric results. In the case of solid-solid reactions, the present paper describes the model developed for describing solid-solid reactions from dynamic X-ray diffraction technique, developed by the present team. The model has been successful in explaining the formation of spinel from MgO and Al2O3.

Jan 1, 2006

By O. Konovalov, M. Zhuravkov

In report the results of performance of research works on building up of the connected geoecological and geomechanical monitoring system for modelling and study of a state of mining massif in a region with large-scale underground mining works are discussed. The connected geoecological monitoring requires radical change of methodology of carried out works. To carry out modelling and researches it is impossible without construction of the integrated connected information model of region. Basis of construction of information systems and databases for all sections of Earth Sciences is the set of the approaches, among which, for example, methods of remote sounding, GIS- and GPS-technology etc. Computer technologies are basis on the "supercomputer" technology. The results of performance of some researches for regions of excavation of deposits of potash salts in Russia and Belarus are presented.

Jan 1, 2005

By N. Syna

The selection and preparation of activated carbon suitable for specific application can often involve exhaustive adsorbent preparation and testing. This work examined an emerging technique based on modelling the correlation of observed adsorption capacity to the carbon characteristics as a means of decoupling the roles of carbon textural and surface chemical properties. The insights to these specific roles are then used in tailor-making adsorbents and in establishing the optimal carbon properties for specific application. In this study, a series of activated carbons prepared from bagasse by physical and chemical techniques and commercially available carbons were used in developing the models. These carbons were tested for the adsorption of thiourea-gold (I) chloride complexes. Thiourea is an important alternative gold lixiviant with the potential to address current ecological and safety concerns in the extraction of gold. A four-component model, which correlated the adsorption gold capacity to carbon textural and surface chemical properties, was proposed. The model fit was assessed using multi-linear regression and a partial F-test. This was followed by goodness of fit test achieved through r2, adjusted?r2 and mean square error (MSE) values. The results suggest the adsorption of the gold-thiourea complex is strongly influenced by the internal surface area and pH of the carbon. Tailor-made activated carbon from bagasse showed adsorption gold capacities of up to 198 mgAu/g adsorbent, demonstrating the potential of bagasse as a precursor for gold adsorbents and of this empirical modelling approach in developing adsorbents for metal recovery.

Jan 1, 2014

By Pedro Lobarinhas

In the present work a two dimensional computational model is derived for the description of the time dependent heap leaching process. The flow inside the heap is calculated using the diffusion equations for the flow in porous media. The soil conductivity is also estimated assuming a close packed structure for the grain packing inside the heap. These equations, with appropriate boundary conditions, are solved using a time step splitting scheme. From this calculation the amount of leaching reagent in contact with the ore at any position inside the heap can be calculated and, once coupled with appropriate reaction equations, the amount of metals extracted can be computed and leaching agent concentration updated. The model is tested with available experimental data and the results have shown to be in a very good agreement. The model predicts all the transient features of the data.

Jan 1, 1999

By Roland B. Drugge

The adaptation of pelletizing processes involves the coordination of various factors such as temperatures, gas flows, productions, etc., in order to obtain optimum performance. The consequent problems are considerable due to the large number involved as well as the lack of knowledge of their correlation. It is therefore essential to establish a fundamental theory which takes into account the heat transfer phenomena and the secondary factors which are relevant to the properties of the product. Straight grate and grate kiln pelletizing processes involve heating of a stationary bed on a travelling grate. The heat transfer conditions are therefore of primary importance. This is why this article pays particular attention to conditions which obtain when a stationary bed is heated, in which certain special reactions occur. The conclusions drawn from this penetration of the heat transfer phenomena can be expressed in postulates acting as guide-lines for the establishment of proto-typal processes. Examples of elaborating such processes are also given in the article. Some operational experience at LKAB are also given.

Jan 1, 1975

By D Adhikary, M Khanal, B Poulsen, H Guo

During longwall mining, the ability to reliably estimate the height of longwall mining-induced fractures is critically important in order to effectively manage mining-induced water inflows and aquifer interference. The extent of mining-induced subsurface deformation, fracture, surface subsidence and aquifer interference is largely controlled by local lithological conditions, mining methods and mine layouts. Therefore for the purpose of quantifying mining impacts it is imperative to incorporate the local geology and mine plan parameters in a mine scale model. Two different numerical models based on 3D continuum (Finite Element Method) and 2D discrete (Discrete element method) domains, respectively, are developed to incorporate the mine site parameters of a study site in Australia. Both the models are calibrated and validated with measured mine data. The continuum numerical modelling method is used to investigate the change in permeability of the overburden strata due to underground extraction considering the full 3D nature of geology and mine layout. The discrete element modelling method is used to simulate the initiation and propagation of fractures, allowing the caving of overburden strata to be represented in a realistic manner, and the permeability enhancement of the strata overlying the goaf to be calculated from first principles. Permeability enhancement in the discrete element model is a direct function of the fractures induced in the deformation/caving process; fracture location, extent, aperture and connectedness. From the numerical modelling under the simulated condition for the mine studied here, it is predicted that the permeability of the caving zone from the goaf to five times the mining height increases by almost six orders of magnitude. In the fracture zone to approximately ten times the mining height, the change in permeability is predicted to be around two to three orders of magnitude. At the height of more than ten times the mining height the magnitude of permeability increased is predicted to be one to two orders of magnitude. CITATION:Khanal, M, Poulsen, B, Adhikary, D and Guo, H, 2018. Modelling of height of connective fractures during longwall mining, in Proceedings The Fourth Australasian Ground Control in Mining Conference (AusRock), pp 254–265 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Nov 30, 2018

By B. Csoke

To improve the efficiency of the crushing/screening plant of Bazalt-Középko Kobányák Ltd. (Hungary), the implementation of computer aided quality control is planned. In the current phase, the R & D work aims at the solution of high-priority tasks: 1) Developing the process models, validation of the models by means of industrial measurements. 2) Developing the algorithm for the process control of the crushing/sizing plant. 3) Developing the suitable computer program. To develop the mathematical models of crushing (breakage and probability functions) and to determine the model parameters, the experimental industrial trials were carried out on hydrocone crushers of types Svedala H 4800 EF and H 4800 F, changing the gap size and crushing relatively narrow grain size fractions. Additionally, the parameters of the separation function for the screens operating in the system were determined. It was established that the crushing probability depends on the grain size of the feed (Y) and the gap size (s), so the breakage probability S, for each cone crusher, can be well characterised by the function formed using the ratio of these two parameters: SfYsf==. The distribution of the broken grains can be described reliably using the relative grain size distribution function B(? = x/Y) (usually Schuhmann-Gaudin- or Rosin-Rammler-shaped function). For the screens, the probability function ()2()111/NTx = xd ???---???? is suitable for cha-racterising the separation (where x? grain size, ?-open-area ratio, d ? rectangular screen aperture size, N ? number of appearance of the particle on a screen surface). Using these functions, the comminution circuits were described by means of vector/matrix equations. The calculations and measurements verified the reliability of the developed models and algorithms.

Jan 1, 2014

By M. Morales, J. Yianatos

"In this work, an alternative approach to evaluate the flotation performance along a flotation rougher row is presented. The collection zone as well as the froth zone performances were separately evaluated. For this purpose, sampling for metallurgical balances along with sampling of top-of-froth (TOF) grades in each cell of a flotation row were employed. The modelling of the rougher circuit consisted of using TOF grades for evaluating the collection zone recovery, because they are a good representation of the bubble loading, which identifies the true flotation process, as shown in previous studies for shallow froths. In addition, the feed characterization, size-by-liberation, allowed the identification of the mineral distribution and its evolution along the circuit. The latter result was obtained by fitting the recovery profiles as well as the TOF grades down the flotation bank. The flotation rates and the maximum recoveries in terms of size-by-liberation (size classes: -45 µm, +45-150 µm, +150 µm, and liberation classes: -20%, +20%-80%, +80%) were then estimated. The results were in good agreement with those reported in literature by other authors and showed the new approach allows flotation kinetics to be characterized in industrial circuits.INTRODUCTIONMineral recovery in industrial flotation circuits is a function of the mineral characteristics, such as particle size distribution and liberation (Trahar, 1981, Jameson, 2012), and operating conditions such as air flowrate, pulp flowrate, solids percent and reagent conditioning (Jameson 1977; Klimpel, 1980; Mehrotra and Kapur, 1974). The cell type and size (self-aerated or forced air) and the circuit arrangement (single cells or banks of cells in series) are also relevant factors. Sampling of circuits (e.g. rougher, cleaner or scavenger) are commonly carried out using standard automatic samplers for metallurgical balance assessment. However, the flotation plants are generally not provided with proper accesses for local streams, which hinders the recovery estimations and the flotation rate characterization."

Jan 1, 2018

By A. J. Lynch, E. V. Manlapig, J. S. Hall, G. C. Thorne

Modelling of industrial sulphide flotation based on data collected in several flotation plants is discussed. A mechanistic approach is used in which it is recognised that flotation is a two stage pro- cess with concentration occurring both in the pulp and in the froth. The important features of the model, which is a kinetic model, are – 1. for each particle species there is a range of flotation rates which corresponds to the range of particle sizes, 2. the flotation rates of the liberated gangue particles which have no hydrophobic properties are dependent on the flotation rate of water and are controlled by variables such as the air and frother addition rates, and the pulp level, 3. as an approximation the valuable sulphide particles may be regarded as being present in two groups with high and low flotation rate constants, 4. increase in the addition rate of collector tends to increase the proportion of sulphide particles which are in the high flotation rate constant group. The actual flotation rate of each group of particles is a combination of an "inherent" flotation rate and a flotation rate which is dependent on the flotation rate of water, 5. at particular air and frother addition rates and pulp levels there is a physical limit to the rate at which particles can be transported into the concentrate and this is recognised in the model. The model has been used to study the steady state and dynamic properties of flotation circuits in operating plants. The assessment of the behaviour of composite particles is difficult because of the problem of quantitatively describing the surfaces of mineral particles. A method which is being developed to measure unpolished particles is described briefly.

Jan 1, 1976

By J. Herbst, A. J. Lynch, E. V. Manlapig, G. C. Throne, J. S. Hall, G. Poling

Modelling of industrial sulphide flotation based on data collected in several flotation plants is discussed. A mechanistic approach is used in which it is recognised that flotation is a two stage process with concentration occurring both in the pulp and in the froth. The important features of the model, which is a kinetic model, are – 1. for each particle species there is a range of flotation rates which corresponds to the range of particle sizes, 2. the flotation rates of the liberated gangue particles which have no hydrophobic properties are dependent on the flotation rate of water and are controlled by variables such as the air and frother addition rates, and the pulp level, 3. as an approximation the valuable sulphide particles may be regarded as being present in two groups with high and low flotation rate constants, 4. increase in the addition rate of collector tends to increase the proportion of sulphide particles which are in the high flotation rate constant group. The actual flotation rate of each group of particles is a combination of an "inherent" flotation rate and a flotation rate which is dependent on the flotation rate of water, at particular air and frother addition rates and pulp levels there is a physical limit to the rate at which particles can be transported into the concentrate and this is recognised in the

Jan 1, 1976

By T Krull, T J. Donohue, D Ilic

"Problematic issues in materials handling applications are often amplified in the iron ore industry due to the relatively high density of the ore as well as the abrasive nature of the product. One of the main issues arising from the handling of iron ore is excessive wear, with two applications typically suffering accelerated wear issues being transfer chutes and conveyor belts. These wear issues are problematic due to additional capital costs, but mainly due to the associated downtimes resulting in the accelerated wear. This study will focus on the modelling of iron ore in practical applications where design work was undertaken in order to reduce the wear rates of different materials handling equipment.CITATION:Donohue, T J, Ilic, D and Krull, T, 2013. Modelling of iron ore in materials handling applications with a focus on wear, in Proceedings Iron Ore 2013, pp 355-362 (The Australasian Institute of Mining and Metallurgy: Melbourne)."

Aug 12, 2013

By Y Chung, B Choi, H Lee

Numerous studies have explored refractory corrosion in the presence of slag and/or metal, utilising the Finger Rotating Test (FRT). Recent experimental results revealed that corrosion initially increased with rising rotation velocity but was then suppressed above a critical rotating speed, under conditions involving a two-liquid solution of slag and molten steel. This phenomenon indicated that the liquid/liquid interface movement occurred during spindle rotor. However, understanding of the phenomenon, known as the rod climbing effect, remains inadequate in the case involving slag, molten steel and refractories. To gain deeper insights, a simulation experiment was conducted on the phenomenon using silicon oil/water and the phenomenon was analysed using computerised simulation software (Ansys Fluent). Experimental modelling demonstrated that interfacial movement occurs even at low rotation velocity, as low as 200 rev/min, when the viscosity of silicone oil is at 100 cP. Computational analysis utilising the Volume of Fluid (VOF) method and a realisable k-ε turbulence model elucidated that interface elevation during rotor rotation arises due to viscosity disparities between the two fluids. Furthermore, this movement is more pronounced when the viscosity difference is lower. Additionally, the study findings suggest that higher density differences between two liquids suppress the phenomenon of interface movement.

Jun 19, 2024

By A. A. Hejja

"Liquidus temperature and electrical conductivity data measured on synthetic slags were modelled by the use of neural nets (NN). In this work the applied multilayer feedforward NNs were trained by a conjugate-gradient optimisation for which a three layer formulation was used. Very good fits were obtained for both the liquidus temperature and the conductivity data. The synthetic slags were prepared from pure oxides to represent a wide range of compositions likely to be encountered in ferromanganese and silicomanganese smelting. The slag constituents were in the following range: Mn0;5-30%, Ca0;20-35%, Mg0;5-15%, Si02;27-58%, Al20 3;5%. Liquidus temperatures varied from 1300° to 1380° and increased with increasing basicity ratio. The electrical resistivity of slags decreased with the increase of basicity ratio from 0.55 to 1.1 but above 1.1 basicity ratio the resistivity tended to increase depending upon the MnO content.Neural NetsVarious neural net (NN) topologies such as Hopfied nets, Hamming nets, Carpenter/Grossberg classifiers have been discussed in some detail by Lipmann (1). In this work, the applied multilayer feedforward NNs are not trained by backpropagation (BP) using the generalized delta rule (2) (GDR), but by a conjugate-gradient optimisation [3,4] - a procedure similar to that proposed by Barnard and Cole [5]. A neural net is capable ·of learning the non-linear functional relationship existing between a set of inputs and outputs [ 6]. In this capacity a neural net is used to determine the functional relationship that exists . between N dependent and independent variables, which characterize the system under consideration For a three layer neural net [1] the following function may be formulated, which relates X; inputs to the output Y 0 from node o:"

Jan 1, 1997

By Bo Zhou

During secondary aluminum recovery, complex aluminum scraps are melted and refined often in a rotary melting furnace and the aluminum metal has to pass through a molten salt layer and melt down in a bottom aluminum bath. Previous experimental study and industrial observations have indicated that salt shell formation and re-melting on the metal solids in the early stage plays a critical role in the overall melting process of the scrap. In the present paper, mathematical models are developed to simulate the process of salt shell formation and its re-melting, as well as aluminum melting in the molten melt. Influence of the salt layer properties, the residence time of the metal particles in the salt layer, the particle size, the temperature of bulk melt were investigated. Reasonable agreement with the measurement data was reached. The developed melting model will be used to construct a population balance model of the scrap melting behavior for a rotary melting furnace, taking into account the distributed characteristics of the scrap. Finally, an overall process model based on Computational Fluid Dynamics (CFD) for the rotary furnace will be constructed, coupling with the population balance model for the scrap melting.

Jan 1, 2003

By Bo Zhou

During secondary aluminum recovery, complex aluminum scraps are melted and refined often in a rotary melting furnace and the aluminum metal has to pass through a molten salt layer and melt down in a bottom aluminum bath. Previous experimental study and industrial observations have indicated that salt shell formation and re-melting on the metal solids in the early stage plays a critical role in the overall melting process of the scrap. In the present paper, mathematical models are developed to simulate the process of salt shell formation and its re melting, as well as aluminum melting in the molten melt. Influence of the salt layer properties, the residence time of the metal particles in the salt layer, the particle size, the temperature of bulk melt were investigated. Reasonable agreement with the measurement data was reached. The developed melting model will be used to construct a population balance model of the scrap melting behavior for a rotary melting furnace, taking into account the distributed characteristics of the scrap. Finally, an overall process model based on Computational Fluid Dynamics (CFD) for the rotary furnace will be constructed, coupling with the population balance model for the scrap melting.

Jan 1, 2003

By A Kamkar-Rouhani

The paper describes a mineral potential mapping of gold and other base metals in Mahneshan area, Iran, using geographical information system (GIS), in which, the method of weights of evidence is used to combine spatial data obtained from different sources to make prediction for mineral potential. This is a quantitative tool for mineral resource mapping that can be used to delineate and predict favorable areas for further exploration by using the characteristics of known and similar mineral deposits. In accord with the descriptive model for gold and other base metal deposits, the spatial distribution of 20 known mineral occurrences was used. In addition, digital evidential maps were used in this study. The maps include maps of lithological units, faults and first derivative of total aeromagnetic data, remotely-sensing data of alteration and also geochemical evidential maps produced on the basis of the results of chemical analyses of samples taken from the stream sediments in the area, and analysed for Cu, Pb, Zn, Sb, As, Ag, Au, Ba and Hg elements. The binary maps were produced and then combined with the model of weights of evidence using 20 known mineral occurrences to obtain a final map of the favourable areas for further exploration of gold and other base metals in the study area. The final map is a predictive GIS model of the mineral deposits in the area. An interesting matter about the final predictive GIS model is that the faults are mainly located in favorable areas which possesses the most probability of mineralisation. Therefore, the faults in these areas have priority for further exploration of gold and others base metals.

Jan 1, 2005

By P B. Hills

This paper discusses the computerisation of the drill hole database, some historical records and the creation of a three-dimensional model of the gold bearing structures for resource estimation and mine design at the Beaconsfield Gold Mine in northern Tasmania. The phase of the work reported here commenced in November 1997 and was completed in March 1998. The historical mining on the Tasmania Reef had reached a depth of 450 m below surface when work ceased in 1914. Production up to this time was 854 000 ounces gold. Recent activity has concentrated on testing depth extensions of the previously mined ore structure, to a depth of 1000 m below surface, and the establishment of an underground mining operation to exploit this resource. The lode structure at Beaconsfield is a typical mesothermal vein style deposit. Gold is contained within narrow quartz, carbonate and minor sulphide reefs or veins. Reef width is variable from a few centimetres to several metres with æsplitsÆ or bifurcations being common. The reef averages 2.7 m true width in the economic zones. The key feature of this style of mineralisation is that reef continuity can be confirmed by drilling to a high level of confidence, but the gold grades are somewhat more erratic in nature. The resource estimates were based entirely on drill data. The drill intercepts on the reef were coded for each split directly into a database and these allocations validated graphically. The distribution of grades and widths show an increasing gold grade with reef width for both the drill hole and historical data. The same histogram patterns, mixed populations and overall arithmetic average grade for the historical mine (at 21.52 g/t > 3 g/t Au) and new drill data (at 21.55 g/t > 3 g/t Au) would indicate that there has been no major change in the overall style or grade of the ore from the historically mined areas to the new resource. The reef, host geology, area of historical workings and resource zones at the Beaconsfield Mine were defined as a series of closed wireframes. The Tasmania Reef was modelled as five mutually exclusive wireframes, one each for the Main Reef, Footwall split, Second footwall split, Third footwall split and Hangingwall split. The data for this exercise was extracted directly from the drill hole database, stored in Microsoft Access format. The formation of the solids was largely automated with manual validation and refinements allowing for multiple runs using different interpolation methods and combinations of drill holes. A rotated 3D-block model in the plane of the Tasmania Reef was created to contain the grade model and allow for tonnage, grade and reef width estimates. An inverse distance squared (ID2) interpolation method was used with ranges of 25 m along strike and 125 m down plunge with blank areas infilled with 50 m by 250 m. Each reef split was filled separately using only drill intercepts from that split. The method was independently validated used MIK methods. Longitudinal sections have been prepared showing gold grade, reef width and width +grade. Estimates were also made for arsenic, sulphur and other minor elements. The long sections show three distinct mining zones. A method is presented for the conversion of the resources to reserves by the addition of mining parameters for minimum mining width, minimum pillar, cut-off grade and internal and external dilution. The selection is done on columns within the block model, and allows for all splits to be considered at once. Multiple runs can be made for varying combinations of parameters to test different mining methods and their impact on reserves, mine life and NPV.

Jan 1, 1999