Search Documents

By M. E. Duncan Fama

The aim of the paper was to model numerically some aspects of the extraction of coal pillars by split and fendering. A two-dimensional elasto-plastic plane strain model with strain softening and interface modelling was used. The model was able to simulate the improvement of conditions experienced in lifting from the fender, when the roof adjacent to the fender has caved. In contrast, when a stook was left, providing support to a span of uncaved roof, the model predicted the fender would suffer considerably more deformation and strain.

Jan 1, 1992

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

By R. J. Pine, J. M. Rance, J. S. Coggan

ABSTRACT: The safe design of excavations created in a fractured rock mass is of paramount importance in rock engineering. Depending on the nature of the stress regime acting on an excavation, and the associated deformational behaviour of the rock mass surrounding the excavation, a number of potential failure mechanisms exist. These involve a combination of shear failure on existing joints, extension of critically oriented joints and propagation of new fractures through previously intact rock. Conventional numerical codes usually adopt either a continuum or discontinuum approach; this often results in significant over simplification and kinematic limitations particularly when modelling complex rock mass behaviour in the close vicinity of excavations. This paper describes the results from several case examples using a combined finite element-discrete element code, ELFEN that incorporates a crack propagation mode. Through the use of selected rock engineering examples, the authors demonstrate the potential of the code in the realistic simulation of complex rock failure both in underground and surface excavation applications. RÉSUMÉ: La conception sûre des excavations créées dans une masse de roche fracturée est d'importance primordiale dans la géotechnique. Selon la nature du régime d'effort autour d’une excavation, et le déformation associé de la masse de roche autour de l'excavation, il exsite un certain nombre de mécanismes potentiels d'échec. Ceux-ci comportent une combinaison d'échec de cisaillement sur les joints existants, la prolongation des joints d’orientation critique et la propagation de nouvelles ruptures dans la roche précédemment intacte. Les codes numériques conventionnels adoptent habituellement l’approche, soit d’un continuum, soit de discontinuum; ceci a souvent comme conséquence une sursimplification significative et les limitations cinématiques en particulier en modelant le comportement complexe de la masse de roche très proche des excavations. Cet article décrit les résultats de plusieurs exemples de cas en utilisant un code qui s’agit d’une combination d’élément discret et d'élément fini, ELFEN qui incorpore un mode de propagation de fissure. Par l'utilisation des exemples choisis géotechniques, les auteurs démontrent le potentiel du code dans la simulation réaliste de l'échec complexe de roche pour les applications aux excavations soit au fond soit à l’extérieure. ZUSAMMENFASSUNG: Das sichere Design der Aushöhlungen, die in einer zerbrochenen Felsenmasse verursacht werden, ist vom paramount Wert in der Felsentechnik. Abhängig von der Natur des Druckregimes, das auf einer Aushöhlung und dem dazugehörigen deformationalen Verhalten der Felsenmasse umgibt die Aushöhlung fungiert, bestehen eine Anzahl von möglichen Ausfalleinheiten. Diese beziehen eine Kombination des Scherausfalls auf vorhandenen Verbindungen, Verlängerung der kritisch orientierten Verbindungen und Ausbreitung der neuen Brüche durch vorher intakten Felsen. Herkömmliche numerische Codes nehmen normalerweise entweder ein kontinuum oder discontinuum Annäherung an; dieses ergibt häufig bedeutende Übervereinfachung und kinematische Beschränkungen, besonders wenn es kompliziertes Felsenmasseverhalten in der Nähe der Aushöhlungen modelliert. Dieses Papier beschreibt die Resultate von einigen Fallbeispielen mit einem kombinierten begrenzten Element - getrennten Elementcode, ELFEN, das einen Sprungausbreitungmodus enthält. Durch den Gebrauch von vorgewählten Felsentechnikbeispielen, demonstrieren die Autoren das Potential des Codes in der realistischen Simulation des komplizierten Felsenausfalls bei Anwendungen für beide Tagebaue und Untergrundsaushöhlungen.

Jan 1, 2003

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 D. J. Reddish

The paper describes a case study where numerical modelling has been applied to predict potential risks associated with inflows of water in active mine workings from an overlying sandstone aquifer. The objective of the modelling was to predict the development of strata fractures and strains as the longwall panels were extracted and investigate their relationship to an sandstone aquifer. The numerical modelling study also included the investigation of the effect of mining on a fault that existed close to the start position of one of the panels. The potential effects of mining on the fault were undertaken by parametric studies adjusting properties for the fault and the relative position of extractions, which allowed the mine planner to position the panel layout accordingly.

May 1, 2003

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 J W. Lee, E Baafi, N Aziz

Ground subsidence due to mining is a common dilemma confronting the underground coal mining industry. The effects of underground longwall mining on river valleys have come under scrutiny, especially when mining is located underneath catchment areas. As public awareness of mining impacts on the environment increases, there is the need to develop damage mitigation strategies. Results of numerical modelling show that river valley response to underground longwall mining can be simulated and evaluated.

Jan 1, 2005

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 R. I. Jeffery, J. W. Summers

An account is given of the verification and calibration of a numerical modelling technique that has been developed for the prediction of subsurface strata deformation following the extraction of coal by longwall mining methods.

Jan 4, 1992

By S. P. Girrens

A coupled finite element geomechanical- hydrology code is currently under development for application to the problem of predicting groundwater disturbances associated with mine subsidence. The structural-fluid coupling is addressed by calculating the subsided mine geometry, with emphasis placed on determining the strata disturbance and locating damaged regions, for input into a hydrology code, which determines localized volume flow rates and aquifer fluctuations. Benefits from coupling will be best realized when field measurements, an additional aspect of the study concurrent with analytical investigations, indicating the relationship between increasing rock strain and increasing permeability are incorporated into hydraulic material descriptions. Hydrologic and struc- tural calculations are presented to demonstrate computational capabilities applicable to mine subsidence

Jan 1, 1982

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

By F. Liu, H. A. Meier, D. M. Valiveti, P. A. Gordon

Hydraulic fracture is a frequently used technique in the Oil and Gas industry to create high permeability pathways for extraction of hydrocarbons. Numerical simulation of this method poses challenges such as representation of 3D fractures, propagation of randomly oriented fractures in 3D media, enabling fluid leak-off from fracture to poro-elastic medium. We have developed a numerical model using the Generalized Finite Element Method (GFEM) to simulate fluid driven fracture in a poro-elastic medium that addresses these challenges. In the GFEM framework, special enrichment functions are used to augment the solution space in the displacement and pore-pressure fields in order to accurately capture the response in the vicinity of the fracture. Adaptive refinement of the finite element mesh is also employed to facilitate accurate and efficient solution of the model. Our model fully couples fluid flow in the fracture, represented with the lubrication approximation, to the geo-mechanical behavior of the bulk poroelastic rock, comprised of an elastic rock skeleton saturated with fluid. Fracturing fluid exerts traction on the fracture surface and is also allowed to leak-off into the poroelastic medium using a transfer function approach. Stress intensity factors (SIFs) are evaluated based on the Contour Integral Method, accounting for the fluid pressure on the fracture surface. The SIFs are used to govern a full 3-D evolution of the fracture. Verification of the model is performed for known behaviors for the poroelastic response and fluid flow in the fracture.

Jan 1, 2015

By Kathy S. Kalenchuk

Damages have been observed in shafts, backfill infrastructure and the slick line at the Leeville Mine. Shaft damage (shearing and spalling of the concrete liner) is attributed to mine-induced rock mass deformations concentrated where regional faults intersect the shafts. Aggregate bin, mixer station and conveyor drive damage (shotcrete shearing, heave of concrete floor slabs and bending of a vertical steel beams) is attributed to rock mass squeezing induced by elevated magnitudes of the vertically orientated major principal stresses. Separation of the slick line has occurred do to extension of the stressshadowed rock mass immediately above the mined out areas. Three-dimensional numerical models have been developed and calibrated to simulated rock mass deformations with magnitudes similar to those believed to have produced the observed damages. By stepping through the mine sequence, models reproduce deformations at the same mining stage as when damage has occurred. Calibrated models have been used to forward project further damage over the remaining mine life.

Jun 1, 2012

By B. K. Mishra

This paper focuses on the use of a numerical tool known as the discrete element method (DEM) to study the motion of ball charge in ball mills (14-18 ft). DEM is employed to simulate the motion of individual balls. While doing so, the effect of liner design on the grinding performance of the mill is studied. Two important factors such as the configuration of the liner and the geometry of the lifting surfaces are correlated with the overall mill performance. An attempt is made to show the effectiveness of the lifter bars in transmitting power from the mill shell to the interior of the mill with particular emphasis to large diameter mills.

Jan 1, 1992

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 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