Search Documents

By Xiaowen (Wendy) Ma, Hans Liang, Xiaoying Chen, David Way

"The mineralogical characteristics of the process streams provide key indicators to the potential improvement in the metallurgical performance. Therefore, the process mineralogy and diagnostic metallurgy have been widely used as a tool to assist the metallurgists at the operating plants. The challenge of this study was then how to characterize the diagnostic flowsheet of the minerals of interest in the process streams across the flotation circuit. The investigation from this study indicates that the minerals chematic flowsheet can visibly demonstrate the sources that caused the gold losses during the mineral processing.In this paper, we present studies of the flotation process stream samples, representing the feeds, intermediate products, final concentrate and final tails, which were produced from the flotation concentrators of Mt Carlton mine in Queensland, Australia. This study was focused on characterizing the process streams and learns how the gold in the process streams was associated with other sulphide minerals and host minerals on a size by size basis. The automatic QEMSCAN/MLA Particle Mineral Analysis (PMA) and Trace Mineral Search (TMS) were applied to identify the minerals and to detect the gold in the sized fractions of process streams. The Laser Ablation ICP-MS (LA) spot analysis, QEMSCAN EDS spot analysis and electron microprobe spot analysis were employed to identify and quantify there fractory gold concentrations in the sulphide minerals, like pyrite, and the elemental concentrations of visible gold grains. The diagnostic flowsheet of gold recoveries by the process streams, by mineral association classes and by the particle sizes were generated and discussed."

Jan 1, 2017

By F. H. Wheeler

CONSIDERING the nature of this audience, I perhaps, in addressing you on the subject of air transportation as related to the mining industry, should speak with tongue in check or at least with some reservation. I know full well that many of you, veterans in the great mineral development in Canada, are perhaps more conversant with the overall aviation picture than I, for aviation to you has been the means to an end, and in many instances has been part and ?parcel in the success of your particular projects, in the form of the initial means of transport. I feel, however, that we in this room today have somewhat parallel geographic interests - yours in relation to what is in the ground, mine in the air - and also in the forests and in the streams - aviation, hunting, and fishing - things which have occupied my attention practically all my adult years. In the light of what has gone. before, and in view of present developments, I think you will agree with me when I say that I do not think there could be a more intriguing pursuit than playing a hand in rolling back Canada's tremendous northern frontiers, even if, at this stage, we have made but slight impression on the untold opportunities which lie in our wide-spread, undeveloped areas.

Jan 1, 1956

By L. Don Leet

Abstract Earthborne vibrations from blasting and other sources, including rock bursts, can now be recorded effectively. The nature and magnitude of vibrations from various sources have been investigated and certain applications of the knowledge thus gained have been made. Vibrations from blasting have been measured on special equipment over a range from an electric blasting cap at a few feet, through the largest quarry shot in history to date ( 1,-362,985 lb.), to the first atomic bomb, in uew Mexico on July 16th, 1945. Laws relating amount of explosive to amount of vibration at different distances have been developed and tested. Special studies have been made of the effects of vibrations on buildings, particularly from the standpoint of possible structural damage. Recorded vibrations are used to study the effects of modern short-period delay shooting. It ?has been found that short-period delay firing is extremely effective in reducing vibrations under the right conditions, but with certain limitations. In quarry shooting, there appears to be a minimum number of holes below which vibrations are not appreciably reduced by short-period delay firing. In several instances this minimum has been near 7. Vibrations in the ground from rock bursts and the effects of mining in stresesd ground have been studied both underground and on the surface in a search for a rock burst warning system. Introduction Earthborne vibrations from blasting and other sources, including rock bursts, can now be recorded effectively. Analyses of blasting records have been yielding increasing funds of information of practical as well as theoretical interest. The nature and magnitude of the vibrations have been investigated with special reference to their effects on structures.

Jan 1, 1951

By C. R. Ross

"THE ASSESSMENT of pneumoconiosis-producing dusts in mines and other industries is usually undertaken for one or more of three reasons: (1) To compare dust exposures at particular operations with standards of permissible dustiness in order that hazards may be evaluated and preventive procedures initiated where required.(2) To study particular operations with the objective of improving dust control devices or procedures by locating dust sources. (3) To obtain data which, along with the results of medical investigations, may assist in understanding the causes of dust-induced diseases or aid in the establishment of permissible standards of dustiness. Most dust measurements are made for the first of these reasons -to evaluate dust hazards so that they may be controlled. For this purpose, four quite different assessment techniques have been in widespread use for many years. One of these is based on the konimeter, an instrument familiar to most .Canadian mining men; a second on the impinger, which is also widely used in Canada, and almost exclusively so in the United States; a third on the thermal precipitator, used extensively in Britain and South Africa; and a fourth on the tyndallometer, a light-scattering instrument used in Germany."

Jan 1, 1960

By F. Ajersch, M. Rigaud, M. Jallouli

"Le gonflement des boulettes d'hématite à été étudié à l'aide d'une méthode de mesure standard ISO ainsi qu'à l'aide d 'une méthode de mesure dilatométrique.La première méthode a servi à l'étude des boulettes haute silice (5% SiO2) dont l'indice de basicité (IB = CaO + MgO -i- SiO2 + Al2O3) a été ajustée soit à l'aide d'une addition de dolomie, soit à l'aide d'une addition de pierre à chaux, soit à l'aide d'une addition de ces deux composes. L'indice de basicité variait de 0.05 à 0.12. Les résultats de gonfle ment obtenus à 1000?C dans une atmosphère de 40 CO - 60 N2 ont permis de montrer qu’indépendamment de l'indice de basicité c'est la dolomie qui contribue le plus à diminuer le gonflement par rapport à la pierre à chaux. L'utilisation de l'indice de basicité comme critère de sélection des boulettes est discutée.La méthode dilatométrique a servi à l'étude du gonfle ment de deux types de boulettes préparées à partir de deux concentrés de fer et dont la teneur en silice a été fixée à 2.5% et 0 5%. Ces deux types de boulettes ont reçues les mêmes additions de liant en quantité et en nature. La détermination en continue du gonflement, à 950?C dans 40 CO - 60 N2, a permis d'une part de mieux cerner le niveau de réduction auquel correspond le maximum de gonflement et d'autre part de montrer que l’allure des courbes de gonflement est étroitement liée à la nature du liant. Une augmentation de la teneur en silice des concentres contribue à diminuer le gonflement des boulettes.Des essais de réduction de gonflement non isothermes ont été réalisés à des montées de température à 10 et 20?C/min. Les résultats indiquent que le gonflement déterminé dans les conditions non iso thermes est inferieur à celui obtenu lors des essais isothermes. On discute alors de la validité des essais standards de gonflement qui sont utilisés pour prédire le comportement des boulettes lors de leur réduction dans un haut fourneau."

Jan 1, 1983

By P. Lind, A. Gagnon

"Goldcorp has embarked upon an ambitious, industry-leading, multi-year water use reduction effort. This is a proactive, voluntary effort that reflects our commitment to truly sustainable mining. This effort will focus on increasing water re-use and recycling, and reducing water inputs ultimately towards zero. All of our operating mines are expected to contribute to this effort and are currently in the process of developing plans to meet short and long term objectives. One of the largest needle movers to increase water re-use and recycling is in the area of tailings management. Within five to ten years at Goldcorp, it is envisioned that the need for wet tailings dams will be eliminated and that all new tailings facilities will consist of only filtered or dewatered tailings. Filtered tailings is developing into a mature technology with wide application in mine waste management. This paper describes a conceptual level evaluation on the transition from a slurry to a filtered tailings operation at the Cerro Negro mine in southern Argentina. Key benefits and challenges are highlighted.INTRODUCTION Goldcorp’s Cerro Negro operation is located in the Santa Cruz province in the Patagonian region of Argentina. Cerro Negro has been in operation since 2014 and commercial production was achieved on January 1, 2015. The processing capacity of the plant is 4,000 tonnes of ore per day and in 2016 produced 363,000 ounces of gold and 3.1 million ounces of silver (Goldcorp Inc., 2017). As the mine is located in a semi-arid area, it is important to evaluate technologies that could reduce freshwater consumption and maximize water re-use and recycle rates. The current tailings operation at Cerro Negro incorporates pumping of tailings slurry into a conventionally lined Tailings Storage Facility (TSF). Water is reclaimed from the tailings supernatant pond to the process plant via pumps, however, water losses occur due to evaporation and solution retention within the slurry tailings voids. Converting to a filtered or dewatered tailings operation would allow for increased water recovery prior to tailings disposal, and decreased water losses from evaporation and from retention in the tailings solids void space. This paper describes a conceptual level evaluation on the transition to filtered or dewatered tailings with the primary objective of moving closer Towards Zero Water."

Jan 1, 2017

By R. Elliott, C. A. Pickles, J. Peacey

As the nickel sulfide ore deposits become decreasingly viable as a source of nickel metal, increasing attention is being paid to the exploitation of the nickeliferous laterite ores. Attempts to produce a high grade ferronickel concentrate from these ores by pyrometallurgical methods face several challenges, in particular the selective recovery of the highly disseminated nickel from the reduced ore. The current work explores strategies for maximizing the size of ferronickel particles formed in order to improve the selective recovery of the ferronickel by magnetic separation. The effect of various coal, elemental sulphur, pyrite and sodium sulphate additions on ferronickel particle formation and size at 1000°C have been investigated. It has been observed that increased coal additions leads to a moderate improvement in the size of the ferronickel particles, while the presence of any sulphur bearing reagent leads to a significant increase in particle size. The retention of sulphur in the reduced material has also been studied, showing sodium sulphate has the highest retention, followed by pyrite and then elemental sulphur. An inverse trend was observed in the efficiency of sulphur retention on the basis of residual mass of sulphur per unit mass of reagent added.

Jan 1, 2015

By T. E. W. Nind

"A short description of the chamber lift method for producing oil wells is followed by some comments on compressor requirements and gas utilization. A quantitative analysis of the process is presented and, based on this analysis, criteria are suggested for selecting wells for which this type of lift might be suitable. The analysis leads to the conclusions that chamber lift is an efficient producing mechanism for wells having very low productivity indices (say about 0.05 bbl/d/ p.s.i. or less), but that for such wells optimum utilization of injection gas cannot be attained. It is suggested that the gas injection time per cycle should be kept as low as possible, compatible with a small oil fallback volume, and that the cycle frequency should be in the range of 20-40 cycles/ day. For wells with productivity indices of 0.2 bbl/d/p.s.i. or higher, chamber lift is relatively inefficient; how-ever, if well allowables are low it may offer an attractive, low cost, alternative to gas lift or pumping."

Jan 1, 1960

By M. Dupuis

It has been previously demonstrated [1,2] that using selective collector bar rodding has a big impact on the uniformity of the current density at the top surface of the cathode block and hence on the intensity of the horizontal current in the metal pad, The present study investigates the impact of using selective collector bar rodding on the MHD stability of a 500 kA aluminium electrolysis cell, For this purpose the recently developed MHD package for the aluminium reduction cell time dependent behavior is used.

Jan 1, 2005

By A. R. MacLean

This mine is situated at Point Aconi, Cape Breton , Nova Scotia, and is a replacement for the 100-year-old Princess Collie ry. It is a slope mine which is the continuation of a s trip mine operation. There are fives lopes from the surface driven with Dosco Dintheaders - this marks the first time these machines were used in North America. There are four room -and-pillar sections, each equipped with a Lee Norse HH106 continuous miner, a "Top Dog" roof bolter and two Wagner diesel Teletrams. There are also two Eimco 913 diesel LHD's , which are used for clean-up in the mines. This paper describes the mining operations of this 800,000-tpy mine which, in less than a year and a half , has been developed to the stage where four sections are producing coal at the desired rate.

Jan 1, 1978

By N Robb

This paper gives an account of the laboratory testing, plant testing and final inclusion of the two-stage sulphuric acid leaching process in the Milliken mill. Operational results for September, 1962, are presented, showing a uranium extraction of 95.5 per cent, with a saving of some 16 lbs. of 100 per cent H~S04 and some 9 lbs. of lime per ton milled. Methods for effecting further improvements in the process are indicated. HISTORY P RIOR to the incorporation of two-stage leaching, the process for the treatment of uranium ores at Milliken consisted of a single-stage sulphuric acid leach, followed by a solution recovery step involving washing in thickeners and two-stage filtration. The thickener overflow solution was subsequently concentrated and purified by ion exchange, and the uranium was recovered, by two-stage precipitation, as ammonium diuranate. Figure I illustrates this flowsheet. For satisfactory extraction of uranium from the. Milliken ore, a high residual acidity, of the order of 40 grams per litre of free acid in the discharge solution, must be

Jan 1, 1963

By Tanai Marin

Modelling of pyrometallurgical processes can become extremely challenging when the process involves complex chemical reactions, phase change, multiphase flows, heat and mass transfer and fluid flow. Traditionally modelling involving complex thermodynamic equilibrium and thermochemistry are treated separately from problems in which transport phenomena are dominant, such as CFD and heat transfer studies. The separation is also necessary due to the lack of modelling tools that can effectively couple both solution thermodynamic and transport phenomena or even multiphysics problems. However, there are, quite often, problems in high temperature process metallurgy in which it is not possible to simplify this two modelling approaches. Typical cases of this are for example, the modelling of slag freezing on furnace walls (electric or flash furnaces) to create a protective layer in front of the refractories hot wall, the formation of build up from the matte or metal phases at the hearth of the furnace, calcine smelting from charging piles into slag phase, etc. The problem is even more challenging in nonferrous pyrometallurgy as most, if not all, of the condense phases are highly non ideal solutions involving several components and exhibiting complex crystallization paths and miscibility gaps, such as the Fe-Ni-Cu-S matte, FeO-Fe2O3-SiO2, FeO-MgO-SiO2 slags, etc. In this paper, an approach to model coupled problems involving multiphysics and solution thermodynamic is presented. The approach consists in developing in-house a library containing the Gibbs energy minimization routines and using the local equilibrium concept for process modelling from a multiphysics commercial software (COMSOL). Simple test cases, such as building a non-ideal binary phase diagram on a CFD package is presented, the limitations and of other classic approaches such as equilibrium or Scheil cooling are discussed and compared to this model.

Jan 1, 2015

By R. D. Brown

The Federal Government's White Paper on Tax Reform received further major criticism in the reports of the Committees of the House of Commons and the Senate, which recently completed a lengthy review including public hearings of the government's proposals. These reports both criticized the emphasis placed by the White Paper on its concept of tax equity rather than on a system which would promote economic growth. However, the Finance Committee of the House of Commons was prepared to accept, with modifications, a number of the White Paper's key reform proposals. The Senate Finance Committee was more critical of the basic thrust of the tax reform program and recommended that a modification of the existing tax system would be more appropriate to Canada's needs. Both Committees re-commended a more favourable treatment of mining income than that originally proposed in the White Paper, and to some extent both Committees suggested proposals that were even more favourable to the industry than those set out in Finance Minister Benson's amended proposals made in August of this year. This article out-lines the principal recommendations and comments in the reports of both the House of Commons Committee and the Senate Finance Committee.

Jan 1, 1971

This paper describes the results of a laboratory study on the upgrading of Algoma siderite by the amine flotation of silica and silicate gangue. The experiments indicate that a virtually complete removal of silica is possible in a rougher flotation con-ducted in tap water at natural pH, without prior desliming of the ore. Grinding considerably finer than the apparent liberation size is essential for maximum silica rejection. A cationic collector, Aeromine 3037, gave best results when used in combination with a modifier. Of the several modifiers tested, a Dutch potato dextrin, Dextrin 4356, was found to be the most effective. Both iron recovery and the grade of the final iron product are enhanced by floating the pyrite contained in the ore prior to the amine flotation of the silica, and by scalping the magnetics from the silica concentrate. F OR over twenty years, the Algoma Ore Properties Division of The Algoma Steel

Jan 1, 1963

By J. B. Boisvert, W. Board, C. R. Mooney

"Traditional geostatistical modelling methods are difficult to apply to heavy-tailed deposits because the extreme positive skew of the grade distribution leads to bias. The focus here is on an innovative method of modelling the high-grade population by conditionally simulating discrete volumes (nuggets) of gold with marked Poisson and Cox processes. Capping and other ad-hoc manipulation of the assay data is not required; however, an intensity map for the point processes is required. The method is demonstrated on a dataset from an epithermal deposit and demonstrates an improvement in reducing smearing and controlling the skew of the grade distribution.RÉSUMÉ Les méthodes traditionnelles de modélisation géostatistique sont difficiles à appliquer dans des gisements dont la distribution de la teneur est à queue lourde car l’extrême asymétrie positive de la distribution de la teneur génère des biais. Dans le présent article, l’accent est mis sur une méthode novatrice de modélisation de la population à haute teneur par la simulation conditionnelle des volumes discrets d’or (pépites) au moyen de procédés marqués de Poisson et de Cox. La limite d’utilisation et d’autres manipulations ad hoc des données ne sont pas requises; toutefois, il faut une carte des intensités pour le processus ponctuel. La méthode est démontrée sur un ensemble de données d’essais provenant d’un gisement épithermal et elle montre une amélioration dans le lissage du flou et dans le contrôle de l’asymétrie de la distribution de la teneur."

Jan 1, 2015

By F. Chachula, A. Molinaro

"Sulphide mineral flotation operations continuously strive for a balance between mineralogy, mineral liberation, circuit conditions and reagent solutions to maximize the production rate and recovery of payable metals. The Taseko Gibraltar copper – molybdenum mill operations underwent an extensive project focusing on the fundamentals to identify the effects of process conditions and mineral liberation to improve the bulk circuit copper – molybdenum recovery efficiencies in the bulk flotation circuits. This discussion will provide process data that demonstrates; 1) fundamental relationship between mineral liberation to recovery and grade, 2) flotation circuit condition testing and optimization, 3) reagent combinations resulting in efficient/economical mineral recovery, and 4) the development of a xanthate free reagent solution. The data presented and discussed will also touch on flotation parameters including pH, retention time, kinetics and the effect on bulk circuit recovery performance.INTRODUCTION Taseko's is the owner (75%) and operator of the Gibraltar copper-molybdenum mine, located in south-central British Columbia, the second largest open pit copper mine in Canada and the largest employer in the Cariboo region. Taseko restarted Gibraltar in 2004. From 2006 until 2012, the company invested ~$700 million to expand and modernize the operation which was originally scheduled for demolition when the company acquired it in 1999. In the spring of 2011, construction commenced for Gibraltar Development Plan 3 (GDP3). GDP3 included the construction of a standalone 30,000 ton per day concentrator, built adjacent to the existing 55,000 ton per day facility. The new design capacity of Gibraltar was designed to be 85,000 tons per day. An integral part of the project was a new molybdenum recovery facility. In addition to the new concentrator, the mining fleet was upgraded with a new mining shovel, production drill and ten haul trucks (Taseko, May 2015). In May of 2015, the Company completed an updated mine plan for Gibraltar featuring a 50% decrease in strip ratio. This new mine plan was a result of a detailed, six-month engineering study and forms the basis of an updated NI 43-101 compliant reserve (Taseko Announces Improved Economics At Its Gibraltar Copper Mine, 2015). The new mine plan took into account lower milling costs as a result of the modernized processing facilities (Taseko, May 2015)."

Jan 1, 2018

By J. H. C. Scrimgeour, T. Toong, R. E. Hamilton

"The principles involved in the use and application of mathematical models as an aid to the development of advanced process control computer systems are described briefly. This is followed by a brief survey of known models of this kind for the simulation of mining industry processes.The main part of the paper describes a dynamic model of the flotation process developed by Canadian General Electric. Using flotation rate constants for each species of mineral and gangue, the model provides a numerical solution of the differential equations describing cell behaviour. The cell mode] is used repeatedly as a subroutine in the •over-all-plant model. The features of the cell and plant models are described.Simulation runs from the plant model show stream .assays and other indications of plant behaviour printed typically in 1-to 10-minute steps for several hours of •Operation following typical disturbances such as changes in feed rate or reagent addition rate. The simulation can :Show expected behaviour for both uncontrolled and computer-controlled plants."

Jan 1, 1970

By Celine Y. Zhang, Chloe J. Feng, Charles Q. Jia, Donald W. Kirk

Oil sands petroleum coke (OSPC) is a by-product of coking - a high-temperature process that converts the mined bitumen to the synthetic crude oil (SCO). For every barrel of the SCO produced, about 20 kg of OSPC are generated. Over 70 million metric tons of OSPC have been produced and stockpiled as a result of the development of oil sands industry in the past few decades. Unlike conventional petroleum coke, OSPC has uniquely high impurity contents, particularly sulphur (~ 7 wt%) and metals (e.g. vanadium at ~0.15 wt%), which has hindered its large-scale applications as a solid fuel. Vanadium is a valuable metal that is widely used as an alloy for strengthening steels. The stockpiled OSPC contains about 100,000 tons of vanadium, while the world’s production in 2013 was only 75,600 tons. At the current price of $25US/kg, the OSPC-bound vanadium represents a total value of $2.5 billion US dollars. In this work, we investigate the technical feasibility of extracting vanadium from OSPC. The work focuses on the characterization of vanadium in OSPC and the conditions for effective extraction of vanadium from the OSPC.

Jan 1, 2015

By J. Jerome MacDonald, Stephen B. Flewelling

"In the base metals industry of today, operations required to continuously improve productivity and reduce unit production costs. Consistent with these goals is the improvement of the muck haulage cycle from open stopes with the No. 12 Mine.Existing haulage from open stopes occurs with LHD units experience at Brunswick indicates that LHD units have a higher unit cost per tonne mucked than loader-truck combinations when the haulage exceeds approximately 150 m.The reintroduction of open stoping at Brunswick made existing loader-truck combinations unsuitable. It was recognized that application of recent continuous loading technology could potentially result in significant unit cost savings and productivity increases.A review of the technology resulted in the testing of the Voest Alpine AL60 loader in the open stopes located on the 1000m Le vel. Testing began in January 1988 and finished in late November 1988. The test will provide experience as well as quantified production and maintenenace parameters f or the purpose of quantifying production unit costs and system productivity, ultimately resulting in a decision to further utilize or not utilize the loader as a primarymucking unit at Brunswick .INTRODUCTIONBrunswick Mining and Smelting Corporation Ltd. produces approximately 10500 tonnes per day from the No. 12Mine located 35 km southwest of Bathurst, New Brunswick Mining is the largest (single mine) producer of zinc in the world having mined over 52 million tonnes since 1964. Many years of production still remain with proven and probable reserves in excess of 90 million tonnes. Reserves grade at an average of 3.65% Pb. 8.99070 Zn, 99g/mt Ag and 0.32% Cu . The operation employs 1446 people with approximately 800underground including operations and maintenance staff.In the No. 12 Mine two major mining methods are presently utilized . The primary method is Mechanized Cut and Fill (MCF) accounting for approximately 6500 tonnes per day from four major areas. The remaining 4000 tonnes per day is produced from open stopes with the majority of 2800 tonnes per day being produced by the 1000 level.The 1000Level contains 30 million tonnes segmented into open stope blocks 30 m high, 15m and 30 m on strike , with a transverse width of 10m to 50 m (Fig. I). Stopes are organized on a primary secondary stope arrangement along strike. Primary stopes are undercut and overcut full size with the back fully cablebolted. All stope and lateral development is driven 4.0 m high. All stopes to date are flat bottom (Fig. 2). Primary stopes are filled with cemented rock- fill , secondaries with conventional rockfill, Haulage from the stopes is primarily with 8-yard LHD units ; mucking is approximately 50% conventional and 50% remote. Haul distance to date has averaged approximately 150m with current stope haul distances of approximately 300m. Haulage is presently to a small jaw crusher situated on the 1000 Level. From August 1989on, haulage will be to a set of ore passes feeding to a large gyratory crusher situated on the 1200 Level. Start up of the gyratory crusher corresponds to an increase in tonnage from 1000 Level to 4100 tonnes per day."

Jan 1, 1989

By H. Qi, J. L. Feng, D. J. Li

"In the present paper, the characteristics of the granite burst under loading/unloading conditions was studied by means of both laboratory experiment and finite-discrete element method (FDEM). After gaining the mechanical properties of the granite samples under uniaxial compression conditions, the parameters taken as input arguments in the FDEM modelling were then carefully calibrated so as to provide some basic data including Young’s modulus, Poisson’s ratio, uniaxial strength, Lame constants, cohesive strength, and frictional angle for further being used in simulating granite bursts. All the simulating results by the FDEM were basically in agreement with that obtained from the uniaxial compression test and triaxial granite bursts in laboratory, respectively. This study shows that that such computational mechanics of discontinua can be employed to gain powerful insight into the failure mechanism of granite burst, which could also be a useful tool to clarify the collapse mechanism of rock structures in rock engineering design and rock test scheme optimisation. INTRODUCTIONBased on discontinuum state a priori, discontinuous modeling techniques that are commonly referred to as discrete element methods (DEM) treat the material directly as an assembly of separate blocks or particles, which was originally proposed by Cundall (1971) and Cundall and Strack (1979) from the viewpoint of analogous molecular dynamics simulation to better account for and understand the interaction between the blocks. The distinct feature of these formulations is that a constantly changing problem of geometry and contact conditions can be addressed by solving the equations of motion using an explicit time marching scheme (most often in a ‘leap-frog’ format), while updating contact force histories as a consequence of contacts between different discrete elements and/or resulting from contacts with model boundaries (Bicanic 1998). The most notable implementations of DEMs are arguably represented by the universal distinct element code (UDEC) (Itasca Consulting Group Inc., 2013) and the particle flow code (PFC) (Itasca Consulting Group Inc., 2012). Another formulation is the discrete crack framework that is based upon some criteria for an advancement of crack front, which is coupled with extensive local and global re-meshing while simulating by numerical methods such as finite element technique, commonly using fracture theories. Certainly, it is relatively easy to trace the propagation of single crack front governed by some principles of energy release rates (Ingraffea 1997) but tracing of multiple cracks leads to topological and algorithmic complexities. Otherwise, there is no consideration of interacting cracks in these formulations of finite element or enhanced finite element techniques (e.g., in XFEM and ABAQUS)."

Jan 1, 2015