Search Documents

By PDF

Flotation columns are widely used in mineral processing and coal preparation industry today. Mathematical modeling of flotation columns has conventionally been done using axial dispersion or tank-in-series models. However, for design and scale up purposes these models are not adequate enough because they do not provide detailed information of the flow fields within the column. This paper presents a Two-Dimensional, Two-Phase. Fluid-dynamic Model and its numerical solutions for some cases encountered in column flotation. The MAC (Marker and Cell) method was employed to obtain the numerical solutions, and based on this method a computer code COLUMNFLOW was developed. The simulated flow field reveals the presence of a circulation flow structure in flotation columns, which could be one of the main reasons for mixing in column flotation. The effects of some parameters such as superficial gas and liquid velocities on the characteristics of the flow structure in flotation columns were also tested with the numerical model.

Jan 1, 1994

By Abou Saghafi

In order to investigate the design parameters of a future longwall in a thick coal seam a numerical model of gas emission from the goaf area was developed. A Longwall goaf was modelled as a collection of different size coal blocks. For each coal block the partial differential equation of gas flow was independently solved. In developing the flow equations, both pressure-driven and diffusion driven gas flow regimes were considered. The total gas emission rate was calculated by integrating the gas release from each individual coal block over the whole goaf area. The main output of the computer program are gas flowrate from goaf and residual content of caved coal in the goaf as a function of distance from the face. The model was applied to a longwall face driven in a thick seam where large volumes of coal in the roof were to be left behind the face. The gas contained in the coal was mainly CO, with a high desorption rate. Simulations were under- taken to quantify the effect of size of the caved coal on the flowrate from the goaf. The flow properties of coal, such as diffusivity, permeability and coal swelling due to gas adsorption, were directly measured on coal samples from the seam.

Jan 1, 1996

By R. L. McNearny, K. A. Barker

A study was conducted to compare physical and numerical models of the block-caving mining method. Results from the numerical analysis closely simulated the blocky behavior of the physical model, including bridging and interlocking of the simulated orebody during draw. The draw-down patterns and the rate of draw developed during numerical modeling also closely simulated the physical models. An unexpected result from the numerical modeling was the development of a zone of downward deflection which preceded the development of the zone of material movement.

Jan 1, 1999

By R. L. McNearny

A comparison study was conducted between an actual physical model and a numerical model based on the physical model of the block cave mining method. Results from the numerical analyses closely simulated the blocky behavior of the physical models, including bridging and interlocking of the simulated orebody during draw. The draw-down patterns and the rate of draw developed during numerical modeling also closely simulated the physical models. An unexpected result from the numerical modeling was the development of a zone of downward deflection which preceded the development of the zone of material movement.

Jan 1, 1997

By M. Z. Emad, M. U. Khan

"Coal mines opting for longwall mining or its variations rely on the main and tail gates for production. Stability of these entries or crosscuts is crucial in coal mining operations. Roof or strata failure can be fatal, may lead to property loss and delay in production. Failure of the gates depends upon many factors including the mining method, panel dimensions, pillar width, mine depth, caving rate, dimensions of gates and roof strength properties. In Pakistan coal is mined through a variation of longwall mining method with almost no mechanization while using timber as a primary support. Timber support as chock is very useful in coal mines especially in overhanging face. Similarly gate entries and coal faces are supported by chock supports throughout their length. In this paper the stability of the longwall gate is assessed for a typical coal mine in Pakistan by using finite element code. Numerical modelling considerations for gate entries are discussed. A model parametric study is performed for key factors. The results showed that the existing mining practice is not suitable and there is a need to review the mining approach. INTRODUCTION Main and tail gate entries leading to longwall mining panel play a crucial role in production cycle. Failure of any of these may lead to loss of life, property and production delays. The stability of gate entries depends primarily on the nature of loading and stresses incurred. It is relatively easier to assess the ground response to production for common longwall mining panel. In thin-seam coal mines the gate entries are generally supported throughout their length by timber made chock supports. The tailgates receive huge deformations and stresses during the mining and backfilling cycle. The chock supports are required to take care of rock falls leading to blockage of entries causing production delays and loss of capital. A number of supports are available for supporting the gate entries (Hoek and Brown 1980; Brady and Brown 1985). The support system can be selected based on load-deformation characteristics including the yielding behavior, pre- and post-yielding behavior (Barczak, 2003). The principal of support selection is that the support characteristics must match the ground response in order to maintain the gate entries (Mucho, et al., 1999). It has been observed that thin-seam mining (<1 m) is very challenging and it is necessary to develop an understanding of ground mechanics before the selection of an efficient support system Dolinar (2003). The chock supports have been tested for load-deformation characteristics in the laboratory (Barczak 2003). But the ground and support interaction, response and loads must be assessed for individual site. This can be accomplished by in-situ testing and observations at few locations. The generated data is generally insufficient for complete understanding. Numerical models can combine the data generated in the laboratory and on-site to fulfill the requirements."

Jan 1, 2017

By S. Finsterle, R. C. Trautz

Numerical modeling studies were performed to investigate water seepage into underground openings excavated in unsaturated fractured rock. Water seepage is an important factor affecting the performance of a nuclear waste repository, such as the one proposed at Yucca Mountain, NV. The amount of water dripping into an underground opening depends on the connectivity and permeability of the fracture network, as well as on the capillarity of individual fractures intersecting the opening. A high-resolution numerical model of an unsaturated fracture network was developed, and the appropriateness of seepage predictions using a fracture-continuum model was examined. Model calibrations to liquid-release tests performed at Yucca Mountain are also discussed.

Jan 1, 2002

By S. Finsterle, R. C. Trautz

We performed numerical modeling studies to investigate water seepage into underground openings excavated in unsaturated fractured rock. Water seepage is a key factor affecting the performance of a nuclear waste repository such as the one proposed at Yucca Mountain, Nevada. The amount of water dripping into an underground opening depends on the connectivity and permeability of the fracture network as well as on the capillarity of individual fractures intersecting the opening. We developed a high-resolution numerical model of an unsaturated fracture network and examined the appropriateness of seepage predictions using a fracture-continuum model. Model calibrations to liquid-release tests performed at Yucca Mountain are also discussed.

Jan 1, 2000

By R. L. McNearny, K. A. Barker

Four numerical models written with the "PFC2D " program were compared against four similar physical models. The numerical models were constrained numerically to closely follow the caving behavior of the physical models. Results from the numerical analyses closely simulated both the granular and blocky behaviors of the physical models, including bridging and interlocking of the simulated orebody during draw. The drawdown patterns and the rate of draw developed during numerical modeling also closely simulated the physical models. The same assumptions and governing parameters used in analyzing the first four models were applied to an additional model. The fifth model employed a completely different, uneven, draw sequence than the other four models, allowing the material to flow in a more granular manner with a demonstrated instability of the undercut sections between each drawpoint opening. Drawdown patterns and behavior of the numerical material of the first four models progressed at uneven and inconsistent rates, similar to those seen in the physical models. By inspection, the overall shape and flow lines of both the numerical and physical models were close in the area removed and in the flow characteristics. The fifth model proved to be the most unstable of all the models tested, with poor correlation to the physical models. The uneven draw sequence, coupled with an additional frictionless wall condition and an added drawpoint, caused the material to flow in a very granular manner. Information generated by this study can be used as a first approach in determining the proper draw rate, drawpoint spacing and undercut thickness for a block-caving operation. Numerical models are often difficult to predict and interpret. However, this study should increase the confidence level for use of the PFC2D program in the future.

Jan 1, 2000

By S. Guo

In this paper, a model for the resin anchor bolt system is developed in which the smooth bar is modeled by using one-dimensional truss elements and the resin annulus is modeled by using axisymmetric elements. These two submodels are combined together and then incorporated into the&apos; existing computer code ROOFBOLT. After this development, all of the components in the resin anchor bolt system including smooth bar, bearing plate and resin annulus can be correctly modeled.

Jan 1, 1995

By J. C. Stankus, S. Guo, S. S. Peng

In this paper, a model is developed for a resin anchor bolt system in which the smooth bar is modeled by using one-dimensional truss elements, and the resin annulus is modeled by using axisymmetric elements. These two submodels are combined and then incorporated into an existing computer code called Roofbolt. After this development, all of the components in the resin anchor bolt system, including smooth bar, bearing plate and resin annulus, can be correctly modeled.

Jan 1, 1997

By Matthew E. Fowler, Jean Habimana, Roy F. Cook

The West Side CSO Tunnel, Pump Station & Pipeline Project is currently beginning construction along the downtown waterfront of Portland, Oregon. This4.3 meter diameter, 5,500 meter-long soft ground tunnel will be driven along the westbank of the Willamette River adjacent to buildings and utilities and several major pile-supported bridges. High groundwater pressure, extended lengths of gravel to cobble-size alluvium, and the need to control ground loss and surface settlements were key considerations in the selection of slurry tunnel boring machines for the project This paper presents an overview of the numerical modeling that was performed during design to estimate the strains and additional loads in the pile-supported structures and to assess the efficiency of proposed mitigation measures.

Jan 1, 2003

By D. G. Fredlund, J. A. Wilson

Conventional disposal of waste rock results in the construction of benches with fine and coarse layers dipping at the angle of repose. This interbedded structure influences flow pathways within the waste rock. A numerical modeling program was conducted to study a four-layer system consisting of alternating fine and coarse waste rock. The system was inclined from the vertical to forty-five degrees with heights varying from one to twenty meters. The results show the effects of gravity as preferential flow developed under unsaturated flow conditions.

Jan 1, 2000

By S. S. SHINDE, P. C. JIKAR, N. B. DHOKEY

Mill scale is an oxide waste product of the steel mills, which contains around 70% iron and other allied impurities. In the present context, a dynamic prototype counter current reactor (CCR) has been used for experimentation and simulation of particle bed by using discrete element method (DEM). Mill scale was pulverized to ASTM mesh range 140 (106 μm)/ + 270 (53 μm), and then it was oxidized in CCR to make single oxide Fe2O3 phase by maintaining required thermodynamic conditions such as pure oxygen at 1100 °C. Then this oxidized powder was subjected to reduction reaction in CCR having gas mixture ( H2:N2) ratio of 1:4 at 875 ± 5 °C for filling degree (fd) and reactor revolution (ω) ranging from 5.38 to 16.14% and 2 to 5 rpm, respectively. A creeping fluid flow condition (Re ˂ 1) and natural heat convection (Gr/Re2 ˃˃ 1) has been maintained while undergoing solid–gas reaction. Besides particle bed behavior has been quantified by Froude number (Fr) for the optimum operating window. The simulation was validated by the experimental results of reactor reaction product by X-ray diffraction (XRD) and scanning electron microscope (SEM) images. Prevailing heat transfer mechanism and reduction reaction mechanism has been established for solid–gas counter current conditions.

Nov 22, 2021

By S. S. SHINDE, P. C. JIKAR, N. B. DHOKEY

Mill scale is an oxide waste product of steel mills, containing around 70 percent iron and other allied impurities. In the present context, a dynamic prototype counter-current reactor (CCR) was used for the experimentation and simulation of a particle bed using the discrete element method (DEM). Mill scale was pulverized to ASTM mesh range 140 (106 μm)/+ 270 (53 μm) and oxidized in a CCR to make a single oxide Fe2O3 phase by maintaining the required thermodynamic conditions, such as pure oxygen at 1,100 °C. This oxidized powder was then subjected to a reduction reaction in the CCR, with gas mixture (H2:N2) ratio of 1:4 at 875±5 °C for filling degree, fd , and reactor revolution, ω, ranging from 5.38 to 16.14 percent and 2 to 5 rpm, respectively. A creeping fluid flow condition (Re < 1) and natural heat convection (Gr/Re2 >> 1) were maintained during the solid–gas reaction, and the particle-bed behavior was quantified by Froude number for the optimum operating window. The simulation was validated with the experimental results of reactor reaction product by X-ray diffraction (XRD) and scanning electron microscope (SEM) images.

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 SK SAHOO, UK SINGH, SK SHARMA, GSP SINGH

Surface mining accounts for most coal production in Indian geo-mining conditions compared with underground methods, whose implementations are in deeper coal reserves and for lower stripping ratio in a congested operating space. The bord-and-pillar method is the preferred underground method in India for reasons of flexibility and investment cost. While previous researchers conducted field investigations on intact overburden conditions, this study focuses on the strata control issues in a bord-and-pillar depillaring working under the influence of the dead-load condition of voluminous and fragmented softcover of overburden dump material. The modeling results helped in ascertaining the criticality of strata behavior and rib-pillar stability in actual mining conditions when compared with the conventional intact overburden.

By Juan Andres Jarufe, Carlos Cifuentes

Numerical modelling has been extensively used to estimate plastic displacement in faults and rock mass surrounding mine excavations. Particularly, for fault slip seismicity, modelled plastic displacement is usually utilized to estimate seismic potential, however these methods assume that all the displacement caused during an specific simulated excavation step, is released as a single event at each step. This representation corresponds to a conservative estimation since rock mass displacement does not occur as one single event but as an exponential distribution of many small displacements and a few large ones. When displacement is expressed as seismic magnitude, this exponential distribution is addressed as the Gutenberg Richter (GR) relationship. This paper present a method to decompose the results from numerical models into a distribution of seismic events The reconstruction of the GR relationship from modelling displacement permits not only a more realistic approximation of the expected maximum magnitudes but also provides a tool to develop a probabilistic estimation of the maximum possible magnitudes.

Apr 26, 2022

By A. J. Basu

Geotechnical design plays a key role in mine design. I A proper geotechnical design provides a safe working environment and facilitates efficient mining operation. There are mainly three design approaches I) judgmental, based on observational and experiential facts, 2) empirical, based on empirical design formulae, and 3) analytical, incorporating numerical analysis. Most design applications are based on the design methods one and two as mentioned earlier. However, complexities arise in deep underground mining operations and deep open pit mines. These complexities and difficult geological environment introduce more uncertainties and simple empirical formulae would not be adequate as a design tool. These situations require numerical modelling and analysis of the design problems including interpretation of the results obtained from the analysis. There are several tools available for performing the numerical modelling and analysis. A brief review of the available tools was addressed earlier (Basu and Kuruppu, 1995). The authors selected ANSYS (a general purpose finite element analysis tool) for use in mining geomechanics and geotechnical design. This paper includes several examples demonstrating the applicability of this software in designing structures applicable to mining geomechnics. The authors hope that the design examples will demonstrate the use of the numerical modelling and analysis application in excavation design.

Jan 1, 1996

By W. Keilich

A methodology of subsidence prediction using the Distinct Element code UDEC has been developed as an alternative for subsidence modelling in the Southern Coalfield, New South Wales, Australia. The models have been validated by comparison with empirical results and observed caving behaviour. At this stage, modelling capability is limited to flat lying terrain. It is planned to apply the methodology to areas of high topographical relief to investigate the mechanics of valley closure.

Jan 1, 2006

By Andrzej Walentek

This paper presents the results of numerical modeling of the fractured zone around longwall gateroads. The numerical computations have been based on the results of underground measurements of the fractured zone around longwall gateroads, which were performed for various geological and mining conditions and measured with the use of the endoscopic method with a digital camera. Having results of the underground investigations, trials of duplicating them by the Finite Element Method using software PHASE and the Hoek-Brown failure criterion were performed. In order that the fractured zone predicted by the model-based computations closely simulate that measured underground, some selected parameters in the Hoek-Brown failure criterion were modified based on the earlier determined relations used for prediction of the roof-to-floor convergence in longwall gateroads. As a result, the model-based computations produce results very close to those obtained from the underground measurements. Determination of the fractured-rock zone around the mine roads is of great importance to mining practice, especially in planning and designing supports for longwall gateroads. The influence of the front abutment pressure causes the fractured rock zone around the roads to increase and results in an increase of support load. Moreover, proper determination of the fractured rock zone is one of basic elements in the process of planning roof bolting for mine workings.

Jan 1, 2009