Search Documents

By S. K. Sahoo, S. K. Sharma, U. K. Singh, G. S. P. Singh

This research focuses on the strata control issues in a Bord and Pillar depillaring working under the influence of dead load condition of voluminous and fragmented softcover of overburden dump material. It reports the modeling of a mine working through the FLAC 2D software under the Indian geo-mining conditions. The modeling results of an actual mine working under softcover were compared with that under the intact overburden condition to assess the influence of the soft overburden over the depillaring working. The modeling work involved progressive caving of the strata and the cyclic goaf filling following the modelsimulated occurrence of main fall and periodic caving. Mohr-Coulomb Failure criteria along with a simplified strain-softening (MCSS) model was used to study the failure and caving mechanism. Salamon’s compaction model was used to simulate compaction and resultant stress recovery in the goaf material. Comparison of the model findings was done in terms of load on supports, abutment stress, and convergence, apart from the stress redistribution and failure mechanism. The load transfer mechanism of the soft overburden was also studied for an overall understanding of the strata behavior in such workings.

Jun 16, 2020

By Andrew Theobald, Peter Dare-Bryan, Trevor Byers

The use of advanced modelling techniques and the timing flexibility of electronic initiation systems allows for unprecedented control over blast outcomes as shown in several complex blasting problems experienced by the open pit operation Xstrata Ernest Henry Mine (EHM) in north-west Queensland, Australia.

Jan 1, 2010

By Ashok Kumar, Rajendra Singh, Petr Waclawik

High-grade coal is locked in the shaft protective pillars at 900 m below surface in the Czech part of Upper Silesian Coal Basin. To optimise recovery of coal from the larger size shaft protective pillar, room and pillar (called bord and pillar in India) is trialled to further develop these on smaller pillars. This development followed a constant gallery width of 5.2 m and 3.5 m height. The monitored pillar size in the first developed panel varied from 860 m2 to 1225 m2 of area. The stability of these pillars (i.e., natural supports) at this depth of cover is observed to be a vital issue to control ground movement and safety of the surface structures. A simple calculation of safety factor revealed that application of an empirical method for estimating size and shape of the pillar is inappropriate due to complex geo-mining conditions of the site. Accordingly, geo-technical instruments are installed in the panel for monitoring the stability of the diamond shaped pillar at different stages of the development. Extensive monitoring of the pillar behaviour by these instruments provided an important set of data for such optimum recovery. An established process of simulation-based pillar strength estimation is often used by the mining industry, but calibration is generally done with the help of empirical formulation. However, for the deep-seated pillar, the strength suggested by the existing empirical formulations (except the CMRI formula) is found to be much different from the instrument observations. Using a study conducted at Czech hard coal mines, an attempt is made to develop a strain-softening-based numerical modelling approach to assess pillar performance at deeper cover. The recorded value and nature of the pillar spalling in the field along with the CMRI formula are used to calibrate the developed simulation approach. An assessment of the stability of the developed pillars using simulated models did not indicate any collapse of the natural supports. This paper briefly presents difficulties in pillar design at great depth of cover and field and simulation results of pillar strength estimation in Czech hard coal mines.

Jan 1, 2018

By O. K. Mahabadi

The scope of this study is to simulate the behaviour of a layered rock sample under standard laboratory Brazilian Disc test using an innovative combined finite-discrete element method (FEM/DEM) research code. The influence of layering and the direction of loading on the samples behaviour are studied through various layering and loading configurations. A parametric study including the effect of loading rate on the tensile strength of samples has been carried out. This paper demonstrates the suitability of the numerical approach to explicitly model rock deformation and failure.

May 1, 2009

By Ian Masters, Saul B. Davies, David T. Gethin

The purpose of this model is to simulate and analyze the energy flows within a rotary aluminum recycling furnace, using the global properties of its main components, the gas, metal and refractory walls. The model will simulate a series of operations such as loading, heating and stirring which will constitute the processing of a batch. Each main component of the furnace is treated as a one-dimensional heat transfer medium, the thickness and area of each medium being calculated from the actual geometry of the furnace. The physical properties are taken as the weighted average of those for the various layers forming the medium. Incorporated into the model are elements describing the radiative exchange within the furnace, melting of the charge, and forced convection in the melt due to the rotation of the furnace. This approach has produced a model with a reasonable computing time that can be used to run a whole series of simulations for the prediction and comparison of the furnace performance under different operating conditions.

Jan 1, 2000

By K Koskinen, S Ronkainen

The flow field induced in an Outokumpu TankCell and by an OK rotor-stator mixing device in a laboratory and a process size flotation cell was modelled numerically using computational fluid dynamics (CFD). The primary target was to discover and reproduce the various interactions between the flow field and the flotation cell by numerical means. As power consumption depends sensitively on flow field, it can be regarded as a top-level derived quantity (or measure/metric) indicating the validity of the numerically obtained flow field. Having properly validated numerical means to produce the flow field the influence of various design parameters can be captured with certain accuracy. In addition, separate gas and solid phase models for bubbles and solids, respectively, were implemented. The comparisons between velocities and turbulence quantities measured using the laser doppler velocimetry (LDV) technique and CFD modelling were performed. The grid generation was performed with the commercial Gambit 1.3 and the subsequent integration and solving of the resulting algebraic equations with the commercial CFD-code Fluent 6.22. The computationally economic multiple reference frame method was included to model the moving rotor blades. Gas and solid phase models were modified mixture models.

Jan 1, 2005

A mathematical model for spontaneous combustion of coal with moisture included is presented in this paper. The one-dimensional unsteady state model consists of conservation equations for oxygen, water vapour and inherent moisture of coal and energy for both gaseous and solid phases. A first order Arrhenius reaction rate for oxidation under both pore diffusion and chemically controlled reaction regime is considered. The rate of evaporation or condensation is also considered as a function of temperature of coal, water content of coal and gas streams. The equation of the model is solved numerically by the finite difference technique. Influences of different parameters on the process of spontaneous combustion can be examined by using this model. Also the model can be used to simulate full-scale storage conditions.

Jan 1, 2005

By Mohamed Sassi, I. Kolo, R. L. Sousa, R. K. Abu Al-Rub, Manhal Sirat

"Fractures induced by folds are pivotal in hydrocarbon exploration, ground water transport and harnessing of geothermal energy. This is because of the need to predict and understand fracture propagation in reservoirs which have folded rock formations. Despite its prominence, there is a lack of reliable numerical models for simulating fractures resulting from rock folding. This is partially due to the difficulty involved in fracture modelling. As a contribution to fold-fracture modelling, this work applies a coupled plasticity-damage model to rock fracturing and anticlinal folds. Parametric studies on a single folded layer give insight into the fold formation and behaviour. The anisotropic continuum damage model which considers both tension and compression is formulated using the power damage evolution law. For ease in formulation, strain equivalence hypothesis is adopted whereby the strain is the same for both damaged and undamaged configurations. Lubliner plasticity yield criterion is adopted for plastic deformation. The model is coded in Abaqus user subroutine UMAT and is applicable to quasi-brittle materials.INTRODUCTIONThe understanding of fractures in the earth’s crust helps in the prediction, evaluation and characterization of fractured reservoirs, whether they are petroleum, geothermal or underground water reservoirs. This is partly because reservoir permeability and stability are dependent on fracture properties (Nelson, 2001). To achieve this, various approaches have been developed to analyze reservoir rock fractures. These could be grouped into (a) discrete approaches – where the rock mass is represented by a finite number of well-defined components and (b) continuum approaches – where the mathematical assumption of an infinitesimal element is made (Jing, 2003). Ivanova (Ivanova, 1998) comprehensively studied fracture systems in nature and highlighted folds as one of the major geologic settings that produce fractures. Naturally occurring folds in reservoirs sometimes compound fracture analysis. It might not be clear whether the fractures formed before, during or after rock folding. Fractures formed during the folding process are expected to show stress orientations related to folding and are thus of primary concern in rock folding simulation. It could be argued that there are no numerical models that as"

Jan 1, 2015

By C. V. Deutsch, O. Babak, J. G. Manchuk

High levels of geological heterogeneity occurring over relatively short distances pose a challenge to the unconventional open-pit techniques used in oil sands mining. There is strong motivation to mine selectively and minimize dilution and loss; however, the large equipment and high mining rate used to meet production targets tend to reduce selectivity. A numerical model is developed to account for geological complexity, mining equipment, regulatory criteria, and other operational parameters. This framework is used for more accurate forecasts, optimization of equipment, and other operational decisions. An example demonstrates the approach in light of Alberta Energy Regulator criteria.

Jan 1, 2014

By J. A. Jankowski

The situation in the raw material market has brought into consideration the commercial exploitation of the mineral deposits on the deep sea bed and in particular the manganese nodules. They occur generally in depths of 4000-5000 m, partially covered or buried in bottom sediments. The deposits and their potential recovery have been intensively explored and researched during the past twenty years (Bischoff 1979, Halbach 1988, Padan 1990). Although the present market situation limits the rentability of deep sea mining on an industrial scale, the mining technology is being further developed in many countries. The stagnation phase in the deep mining issue allows sufficient time for conducting research on the precautionary protection of the deep sea environment and for consultation in the development of technologies with the least possible impact on it. The subject has gained in importance since autumn 1994, when the United Nations International Law of the Sea was signed. This also regulates human activities in international waters from the environmental point of view. As fields tests have shown in the past, even with the most modern technologies, mining of the manganese nodules will inevitably be accompanied by a resuspension of some bottom sediments. The generated sediment plumes will be further transported with the ocean currents and in the worst case form stable turbidity layers in different depths (Ozturgut 1981a, Baturin 1991). One of the problems of the environmental impact assessment is the estimation of the plume persistence before it eventually dilutes to approximately the ambient concentration level or settles at the bottom (Lavelle 1981a, Lavelle 1981b, Lavelle 1982, Thiel 1991a, Thiel 1991b). The aim of the research is to develop a numerical sediment transport model in connection with potential manganese nodule mining in a reference area located in the southeastern equatorial Pacific Ocean off Peru (Peru Basin, Discol Experimental Area). In this region, a German mining claim is registered and interdisciplinary field works of the TUSCH (German: TiefseeUmweltSCHutz, deep sea environmental protection) research group are carried out in order to obtain the required scientific background for evaluating of the deep sea mining impacts and to provide the necessary supporting data for the modelling (Thiel 1991a, Thiel 1991h).

Jan 1, 1995

By B BussiFre, J Molson, M Aubertin, G W. Wils

Waste rock piles containing reactive sulfide minerals can serve as long-term sources of acid mine drainage (AMD). Within these piles, many physical, geochemical and biological processes can contribute to the production of AMD. The flow of water, for example, depends on internal pile structure, which depends on the dumping methods used to construct the pile. In this paper, the internal structure and grain size distribution of typical waste rock piles are first reviewed with respect to their effect on internal water flow. A numerical model is then applied to simulate unsaturated flow within a waste pile assuming two primary waste materials, different atmospheric conditions and several geometric configurations. The simulations suggest that water movement can be controlled by introducing fine-grained material layers within the coarser host material. However, these layers should not be horizontal because they can lead to localised high-saturation regions as well as preferential flow favouring the production of AMD. A shallow slope towards the pile boundary may be introduced to force drainage to the outside while maintaining low saturation within the interior, thus reducing the production of AMD.

Jan 1, 2003

By Henry Hu, Alfred Yu

"A numerical model has been developed to simulate heat and mass transfer phenomena occurring during a novel magnesium reduction process-the vertical retort technology. The model was based on the control - volume finite difference approach. Simulations were run to determine the effect of various processing parameters on magnesium reduction cycle time. The model predicted temperature distributions, heating curves, recovery ratio of magnesium, and total process times. The predictions were used to optimize the magnesium reduction process including dimensions of the retort, geometries of the reactant compounds, and productivities. The results indicate that the implementation of the optimized processing parameters significantly decreases the reduction time, and consequently increases the production rate as well as lowers the energy consumption.IntroductionMagnesium consumption has been growing considerably in the past decade [1]. Of the total magnesium production, 30% is produced by the thermal reduction (Pidgeon) process, and it is anticipated to be still the main technology for producing magnesium in future [2]. One of the advantages of the Pidgeon process is that it is capable of producing high purity magnesium demanded by both the automotive and aerospace industries. However, the inefficiency of heat and mass transfer taking place in the horizontal Pidgeon reduction process results in low productivity, high energy consumption, labor-intensive operation, large floor space and high capital investment. This limits the production scale of a horizontal Pidgeon reduction process [3].Vertical reduction processTo overcome the shortcomings of the horizontal magnesium reduction process, a novel vertical magnesium reduction process [3] for large-scale production has been developed, in which heat and mass transfer in the retort is enhanced. The new process shortens reduction time, reduces energy consumption and increases productivity. The vertical reduction system is comprised of a vertically-placed retort, a heating furnace and a vacuum system. A section of the vertical retort is illustrated in Figure 1. Since the retort is set vertically, instead of horizontally, the entire procedure of loading the reactants material in the new process is different from the one for the horizontal reduction process."

Jan 1, 2001

By Hongbo Yang, Tongbo Zhang, Ting Huang, Jingshe Li, Fangfang Song

"Liquid core length is a very important parameter of the process of continuous casting, which directly reflects the cooling status in secondary cooling zone and solidification rate of the slab. In this paper, the billet caster in some steel plant in China was chosen as the prototype of simulating. A numerical simulating model for calculating the liquid core length of continuous casting billet was built by ProCAST software, which was verified reasonable and accurate by pin-shooting technique. Based on the model, the relationships between casting speed, casting temperature and liquid core length of casting billet were found, which was very practically useful for the workers on sit to improve casting billet quality. IntroductionContinuous casting process developed very rapidly because of the high metal yield, energy saving, and low-cost 111• And the quality requirements of continuous casting billet for customers became higher and higher with the fast development of economic121• Internal defect usually is the major quality problem for continuous casting billet[3•61 , which is influenced by the solidification process of casting billet, especially the solidification process of liquid core17•81 • So it is very important to study on the relationship between casting parameters and the liquid core length which means the distance between the crystallizer meniscus and the casting billet solidification end-point. In practical production process, the workers usually control the quality of casting billet according to the liquid core length."

Jan 1, 2013

The paper describes a simulating method.for ground water contamination originating from industrial and mining activrty. The finite difference model simulates heat, mass and solute transport. The applied stochastical method for dispersion provides more accurate solution of mass transport modelling. This method gives a more realistic approach of chemical processes in an aquifer depending on temperature.

Jan 1, 1988

By T. Miura, Y. Mori, H. Aoki

In a flash smelting furnace, easily combustible concentrate particles excessively oxidized and melt. The excessively oxidized molten particles then collide with less oxidized solid ones to create uniform oxidized molten particles. This collision phenomenon plays a very important role in achieving high smelting performance. The authors developed a mathematical model describing above combustion phenomena. This model incorporated fluid flow, heat and mass transfer, chemical reactions and collisions of concentrate particles. Both gas flow and particle motion were calculated using the Eulerian method. Copper concentrate was assumed to consist of chalcopyrite(CuFeS2). Considered in the concentrate reactions were the decomposition of CuFeS2, and the oxidation of the resulting sulfur(S) and pyrrhotite(FeS) to produce magnetite(Fe30 4) and sulfur dioxide(S02). Moreover a reaction due to concentrate particle collisions, in which Fe304 in the excessively oxidized molten particle reacts with FeS in the less oxidized particle, was also considered. The particle collisions were assumed to occur according to a collision probability rule obtained by computational experiments. Particle growth in diameter was described as a change in the volume fractions of particle phases. When applied to a commercial flash smelting furnace, this model reproduced the particle growth and component transfer as observed in the actual furnace.

Jan 1, 2007

By Seyed Ahmad Jenabali-Jahromi

In this paper a meshless method called Finite Point Method (FPM) is developed to simulate the solidification process of continuously cast bloom steel in both primary and secondary cooling region. The method is based on the use of a weighted least-square interpolation procedure. A transverse slice of bloom as it moves with casting speed is considered as computational domain. The two dimensional heat transfer equation together with temperature dependent thermophysical properties is solved in the computational domain. The present work is verified in the mold region by the comparison of the surface temperature simulated by FPM and finite volume method (FVM) and also comparison of solidified shell thickness simulated by FPM and measured on a breakout bloom. For secondary cooling region the validation is done by comparison of the surface temperature simulated by FPM, FVM and thermovision measurements.

Jan 1, 2009

By S. H. Seyedein

A mathematical model was developed to study fluid flow and solidification in a single-roll casting system. Fluid flow coupled with energy equations were solved numerically using a finite volume method. A body-fitted coordinate transformation technique was adopted to handle the irregularly shaped domain. In order to couple the continuity and momentum equations, a SIMPLEC algorithm was used. A fixed-grid enthalpy method was considered to model the solidification in the system. The turbulence effect on flow field was modeled using a low-Re k-e model The model developed in this work was further run under various conditions to optimize the operating parameters of the process.

Jan 1, 2004

By Murali M. Gadde, Syd S. Peng, C. T. Holland

Cutter roof failure is one of the most common ground control problems affecting the safety and economy of an underground coal mine operation. Several attempts were made in the past to understand the processes responsible for cutter development and to devise mitigating measures. Yet, many cutter failures still happen and their advance prediction remains as elusive as before. One of the reasons for the impasse in predicting cutters lies in the limitations of the current numerical modeling technologies. The development and propagation of cutters is a progressive phenomenon and continuum models can only approximate this process. Since the discontinuum modeling is severely restricted by the lack of reliable input data, alternatives must be sought to better understand cutter failures using the continuum models. This paper proposes one approach using strain-softening material model to simulate the cutter initiation and propagation process. An attempt has also been made to provide some insights into the reasons for the change in cutter location at different parts of a development layout.

Jan 1, 2005

By K. Krishnapisharody

A mathematical model has been developed to predict the interface dynamics and droplet generation of the buoyant phase in a two-fluid system. The model employs the well-known SIMPLE algorithm for obtaining flow fields and the Volume of Fluid (VOF) method, for interface tracking. The implementation of the model in water-oil systems, where the oil layer is placed above a recirculating water bath, is discussed. The evolution of the water-oil interface leading to the formation of oil droplets is computed realistically.

Jan 1, 2008

By Guo Chen, Yifeng Shi, Xiaobiao Shang, Nanshang Shen, Bangqi Zhang, Jinhui Peng, Junruo Chen

"The energy is introduced directly into the volume of the honeycomb ceramic regenerator due to the complex interaction of microwaves with minerals in microwave gas heating system. A three-dimensional numerical model is developed to simulate the electromagnetic field by solving Maxwell's equations using finite-difference time-domain method. Computational result is presented for the microwave power absorption of the regenerator. Effects of the dielectric constant and the effective loss factor on the microwave power absorption in the regenerator are studied as well. Numerical results show that the heating pattern strongly depends on the dielectric property of the regenerator. The demonstration of the model of electromagnetic field distribution can be applied effectively and efficiently to optimize the fundamental parameters in the microwave heating processing.IntroductionIn recent years, microwave energy has been utilized in many fields, such as food processing, chemical engineering and mineral processing [l]. In contrast to conventional heating, where heat is transferred from a surface to the interior, microwave energy is introduced directly into the volume of minerals due to the complex interaction of the microwave with the mineral. Based on volumetric heating effects, microwave heating has several advantages such as pollution free, rapid heating technique, internal heating, high efficiency over conventional heating [ 1-4]."

Jan 1, 2013