Search Documents

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 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 W. M. Ma, Daniel W. H. Su, K. Centofanti, Yi Luo, W. L. Zhong, Syd S. Peng

NTRODUCTION Surface subsidence monitoring is the most effective means for investigating the nature of, and developing the techniques for pre¬dicting the surface movement and deformation process. However, because of imperfection in survey equipment and field practices, the surface movement data collected from field survey are not error-free. Also due to lack of proper numerical techniques for computing surface deformation from the measured surface move¬ment data, a relatively large portion of the subsidence data col¬lected could not provide information other than surface movement. In a paper by Luo and Peng (1992), the sources and magnitude of errors in the measured surface movements using an electronic total station in subsidence survey are discussed. If the errors in the measured surface movement components are left untreated before calculation for deformation indices (i.e., slope, strain and curva¬ture), the errors will be carried over to the deformation indices and enlarged if improper numerical procedures are used in the calcu¬lation. It should be noticed that in protecting structures affected by ground subsidence, accurate surface deformation indices are much more important than surface movement. In order to improve the accuracy of the collected surface movement data and the calcu¬lated deformation indices, a number of numerical procedures are also proposed in the same paper. The proposed techniques include: a data smoothing technique for data pre-processing; and a number of numerical methods and/or procedures for calculating surface deformation indices from the measured surface movement data. DATA SMOOTHING Because of various factors, the measured surface movement data always contain errors, small or large. In order to better utilize the collected data, a data smoothing method is proposed to correct the errors in the measured movement data. The method uses the same concept as the influence function method (e.g., Knothe, 1957) on which many existing subsidence prediction methods are based. The data smoothing method makes necessary corrections on the measured movement for the point of interest according to not only the measured movement at this point but also those at the surrounding points. The contribution of a surrounding point in making the correction for the point of interest is inversely propor¬tional, in exponential form, to the distance between them. The new magnitude of a movement component at the point of interest (Zi*) is calculated by

Jan 1, 1992

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 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 Sang-Joon Kim, Hae-Geon Lee, Hyun-Jin Cho

"Decarburization phenomena on the surface of liquid iron have been simulated by numerical method and the results were compared to the experimental data. During the DeC process in a converter, a large amount of splashed droplets are generated from the top surface of the bath and consequently the reaction surface area is increasing tremendously. Since the decarburization mainly occurs at the interface, it is important to understand the reactions on the liquid iron droplets. The flow and species concentration fields around liquid droplets taking into account the reaction of carbon with oxygen at the interface were calculated by a commercial CFD code FLUENT which solving multiphase Navier-Stokes equation with Volume Of Fluid method, species transport equations and energy transport equation. The initial carbon concentration was fixed to 4 wt% and the concentrations of oxygen in the atmosphere were varied from 5 to 50 wt%. As comparing the numerical calculation with experimental data written in the reference, the results have a good agreement for decarburization rate.IntroductionIn a primary steelmaking process, one of the main functions is to decarburize from hot metal and the 4% of carbon is decreased to about 0.5% of carbon. Even though there are some possible reactions leading to the removal of carbon from hot metal, the main reaction occurs on the free surface by direct oxidation of carbon at the gas-liquid interface with producing the carbon monoxide gas. In many metallurgical processes involving a jet impinging onto a steel bath surface, especially in a Basic Oxygen Furnace (BOF) converter, it is important to understand the behavior of molten iron droplets splashing out of the steel. The iron droplets are generated by an oxygen or air jet blowing and it increases the total surface area of liquid iron tremendously. Since the decarburization reaction is mainly affected by the surface area, it is necessary to verify the decarburization phenomena around the splashed droplets from liquid iron in order to allow an analysis ofBOF process possible."

Jan 1, 2011

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

By Nan Wang

Based on the commercial computational fluid dynamics software, Fluent, the fluid flow and temperature field in an iron-bath reactor were simulated numerically, and the effects of top-blown and side-blown operating parameters on the mixing time of iron bath, distribution range of oxygen in top space, and post-combustion ratio were investigated. The simulated results showed that the mixing time of iron bath increased initially and then decreased with the increase of inclined degree and immersion depth of side-blown lance. The distribution range of oxygen became wider as pulling the top lance away from the iron bath. Furthermore, post-combustion ratio was reduced given the top lance was too close the bath. The simulation results provide theoretical base for the process optimization of iron-bath reactor and further development of new type iron-bath smelting reduction process.

Jan 1, 2009

By M. K. Ameen, R. El-Hacha

This paper presents the analytical results obtained from numerical simulation of reinforced concrete slabs strengthened with various FRP systems and subjected to blast loading. A three-dimensional finite element model using AUTODYN software was used to predict the overall behavior of the slabs. The numerical results were validated by comparing with experimental results from a program carried out to investigate the efficiency of different strengthening systems on the blast resistance improvement of small-scale reinforced concrete (RC) slabs. The experimental program consisted of testing nine RC slabs strengthened with three different externally bonded (EB) FRP systems attached at the bottom side of the slabs and one RC slab tested without strengthening for comparison purposes. The EB FRP strengthening systems investigated included two layers of Carbon Fiber Reinforced Polymer (CFRP) sheets, two layers of Glass Fiber Reinforced Polymer (GFRP) sheets, and one layer of Glass-Fiber Resin Polymer (GRP) sheet. Different epoxy resin, patching mortars and impregnating epoxy adhesives were used for the slab surface preparation and FRP system application of the strengthening layers. The aim was to evaluate the residual strength of the test specimens after exposure to blast loads. The blast loading parameters were evaluated for each test using AUTODYN simulation and compared with the computer program CONWEP. All specimens were subjected to same spherical TNT charge weight and stand-off distance with the respect to the slab top surface. Extensive comparison is made between the experimental results and the AUTODYN simulations for the control unstrengthened slab and the strengthened slabs including time history records of the bottom steel strain and mid-span acceleration and the peak dynamic deflection of the slabs.

Jan 1, 2010

By Yuting Xue, Danqing Gao, Brijes Mishra

Heterogeneity and discontinuity significantly affect rockmass strength. Past studies have largely used arbitrary scaling approaches to produce rockmass strength from laboratory rock tests. This paper proposes a stochastic approach to produce rockmass strength from laboratory rock tests. Based on the laboratory data, MATLAB with an extreme value stochastic model generated a database for each physico-mechanical property of the rock. Then finite-difference software FLAC2D simulated laboratory-sized rock specimens. Random material properties in MATLAB were generated using the grids developed in the numerical model. The random properties were then applied to the final FLAC2D model. Model runs simulated the approach performed in the laboratory. The results from the model indicate that a stochastic approach produced compressive strengths that were lower than those from the deterministic approach. Failure modes for each specimen were different from those observed in the laboratory. In addition, random density also influenced the failure mode, highlighting the importance of stochastic analysis in rocks.

By Yuting Xue, Danqing Gao, Brijes Mishra

Heterogeneity and discontinuity significantly affect rock strength. For accurate stability prediction, intact rock behavior is imperatively included in rock mass behavior. However, past research largely used arbitrary scaling approaches to produce rock strength. This paper considers a stochastic approach in order to determine the strength of a rock. Based on the laboratory data, a MATLAB® with extreme value stochastic model generates a database for each physico-mechanical property. Then FLAC® simulates laboratory-sized rock specimens. The grids developed in the numerical model can in turn develop random material properties in MATLAB®, which researchers then applied to the final FLAC® model. Model runs simulate the approach performed in the laboratory. The results from the model indicate that a stochastic approach produces strengths that are lower than a deterministic approach. Failure modes for each specimen are different, similar to observations in the laboratory. In addition, random density also influences the failure mode, highlighting the importance of stochastic analysis in rocks.

By Toshio Suzuki, Machiko Ode, Won Tae Kim, Seong Gyoon Kim

"The Ostwald ripening and growth and coarsening of the secondary dendrite arm are studied using a phase-field model. The phase-field parameters are obtained at a thin interface limit proposed by Kim et al. In Ostwald ripening analysis, the difference of the concentration due ·to· the curvature effect is shown. The third power of the average particle radius Rav is proportional to time and the ripening constant shows a good agreement with the value by Deppers et al. The particle radii distribution shows the transition behavior to attain a self-similar stationary shape. The growth and coarsening behavior of the secondary dendrite arm is also analyzed using a phase-field model. The estimated arm spacing and its dependence on time are in good agreements with experimental data.IntroductionThe ripening or coarsening process controls the microstructure during the late stage of solidification. The ripening phenomena are recognized firstly by the Ostwald [1]. The first quantitative analysis for the Ostwald ripening in the limit of vanishing solid fraction was provided by Lifshitz, Slyozov and Wagner (LSW) [2,3]. In their study, that the mean radius grows proportional to the cubic root of time and the size distribution of precipitations shows a ""self-similar"" shape are shown. Most of the studies for the Ostwald ripening confirm or predict these two general rules above. In case of finite solid fraction, the general rules are observed but the ripening constant is larger and the self-similar shape is wider than the LSW prediction respectively [4-6]. To analyze· the high solid fraction system, the theories that consider the solute field correlation are proposed [4,7-8]. In the theory, however, since the particle is treated as the point, the effect of the concentration field deviation from the monopole field should be taken into account. Recently the numerical simulation, in which not only diffusive transportation of solute but also convection effect is considered, has been done using a phase-field model [9]."

Jan 1, 2000

By T. G. Sitharam, Resmi Sebastian

"Wave propagation studies across jointed rocks have significance in the fields of rock mechanics and geophysics. The passage of a seismic wave induces strains in rocks and often, these strains are of low magnitudes except at the vicinity of source. These low magnitude strains are generally in the elastic range and hence elastic constants necessary for the design and analysis of strength of rock mass, can be derived from wave velocities. Measurement of travel time velocities and measurement of resonance frequencies are the two methods usually adopted for obtaining wave propagation velocities.Resonant column apparatus measures wave velocities from the resonant frequencies obtained during testing. Depending on the type vibration applied, torsional or flexural, shear and compression wave velocities are secured. By varying the amplitude of vibrations applied, various strain levels can be introduced and corresponding resonant frequencies and wave velocities can be calculated. The resonant column experiments were performed on a model material, hardened Plaster of Paris (POP) samples, in laboratory. The paper explains the preparation of samples and the working mechanism of resonant column apparatus along with the estimation of wave velocity from resonant frequency. The numerical simulation of resonant column apparatus using distinct element method and its validation with the experimental results are presented in the paper.A distinct element code is usually used to simulate the movement of materials with intersecting discontinuities. Finite rotations and displacements of individual blocks representing the material are permitted and new contacts between the blocks are identified as the calculation progresses. Three dimensional Distinct Element Code (3DEC), the software intended for the numerical simulation of rock engineering projects, can well simulate the behavior of jointed rocks under static and dynamic loading conditions. The numerical simulation of resonant column apparatus has been performed with 3DEC considering its efficiency in modeling jointed rocks.The numerical modeling of the resonant column apparatus with application of torsional and flexural vibrations to the sample is elaborated in the paper. The response of samples to a range of frequencies is obtained and the frequency that induces maximum displacement in the sample is measured as the resonant frequency, for any particular strain level. The displacements measured at various points of the sample during application of vibrations are compared and presented. The numerical modeling of resonant column apparatus with an intact sample (without joint) and a sample with a joint incorporated are explained at a particular strain level. The resonant frequencies obtained from laboratory experiments and numerical simulations are compared and illustrated. It has been found that resonant frequencies obtained using numerical and laboratory experiments are in good comparison and can be extended for various strain levels and sample conditions."

Jan 1, 2015

"A numerical procedure based on Smoothed Particle Hydrodynamics (SPH) is developed for rock blasting phenomena involving multiple intersecting joints. In this procedure, an elasto-plastic damage constitutive model has been applied in the framework of SPH to simulate yielding response of the surrounding rock medium. In this paper, a square shape rock medium with centrally located emulsion explosive in a borehole is numerically blasted with and without the presence of two diagonal rock joints. The model is simulated by assuming all four sides of the rock medium as free surfaces. The paper elaborates the propagation of stress wave, expansion dynamics of gaseous products, and development of fracture, fragmentation of the medium and dissipation of energy due to damage.INTRODUCTIONDamage prediction of rock mass and its mean fragmentation size are of important to design blast parameters for a given explosive. Current blasting practice solely depends on empirical relationships between fragmentation sizes to those blast design parameters while mathematical descriptions of the physical laws governing the process of rock blasting are limited. A constitutive relationship based on properties of explosives and rock mass is of special interest because of its capabilities in describing the key features of rock fracture and fragmentation, resulting in understanding of the physical process involved in rock blasting. It is well known that blast-induced damage and fragmentation of rock mass occur due to detonation induced stress wave and product gas driven fracture propagation. Nevertheless, the fundamental mechanism of the failure process and the behavior of rock characteristics under these two different loading are not well understood in many cases, which are closely related to growth of fractures and interaction of cracks under different loading conditions.Continuum damage mechanics is considered to be appropriate to describe the failure process under such high-amplitude stress wave. Based on this approach, several theoretical models have been developed to analyze the dynamic damage of brittle material such as rock and concrete with microstructure (Chen, 1999; Taylor et al., 1986; Wu et al., 2004; Yazdchi et al., 1996). These models are mostly based on the theory of damage mechanics or fracture mechanics with some statistical treatment to incorporate random distribution of microcracks. In the application of fracture mechanics theory, sometimes it is difficult to estimate the micromechanical parameters such as crack density and fracture toughness in different modes. In this study, damage in rock material is evaluated according to the earlier model developed by Grady and Kipp (1980), which is a fracture model coupled with a material description for stress wave propagation to predict quantitatively fracture and fragmentation under explosive loading conditions."

Jan 1, 2015

By B. Panton, J. Carvalho, E. T. Brown

"A major consideration in the design of high capacity tiedown anchors for dam, bridge and tower foundations is the tensile resistance of the rock mass to pullout, typically as a result of overturning moments or hydrostatic uplift. Rock mass pullout capacity for installed anchors is developed from the tensile strength and fracture propagation properties of intact rock, the orientation and physical properties of the discontinuities and anchor confinement at depth. The typically assumed, but generally conservative, design approach is to calculate anchor pullout capacity using the dead weight of a uniformly shaped inverted “cone” with an assumed initiation point and breakout angle. In some cases, some tensile or shear strength is considered across the area of the assumed pullout cone, however, these strength estimates are based primarily on based on limited historic testing data, empirical relationships for rock mass strength developed for other applications, and engineering judgement.As an alternative to the current foundation anchor design method, this paper presents a new approach to the rock cone pullout problem using Discrete Fracture Networks (DFN) combined with numerical simulations. The simulations presented in the paper investigate the influence of fractures in a synthetic rock mass on ultimate anchor strength. In addition, the global influence on fracture network statistics on rock mass pullout failure is also discussed and it is contended that some of the aleatory and epistemic uncertainty inherent in the design of rock anchors can be reduced by considering fracture statistics in anchor design."

Jan 1, 2015

By A. J. Basu

Rock mass properties testing for determining strength parameters is an integral part of geomechanics. These tests are performed either in a laboratory environment or in a field (in-situ). The authors share their experience in using numerical modelling in simulating the laboratory test procedures. This approach helps studying the rock breakage and failure mechanism more efficiently and provides a high degree of visualization. Moreover, imperfections and discontinuities can be easily incorporated in these models. This capability of providing arbitrary fractures and discontinuities facilitates in building a large knowledgebase of rock strength and failure mechanisms. This knowledgebase may then be used to generate heuristics for determining rock strengths for rock mass with fractures and discontinuities.

Jan 1, 1997

By Bernhard Hiebl

Global trends that are constantly challenging the iron and steel producers are related to worsening of the iron ore grades, higher substitution rates of main iron ores to cheaper raw materials due to economic pressure and increased environmental standards. Changing the raw materials in an operating plant can significantly influence the sintering or pelletizing process. Adjusting all process parameters to regain stable and high productive operation takes considerable time and causes production losses. Furthermore changes in raw materials lead to variations in the process gas flows and emission concentrations.Therefore SVAI developed two models for numerical simulation of sintering process and pelletizing process. Amongst others the models are based on Ergun equation and include chemical reactions of the related process. For already operating plants the simulations shall predict changes of the process parameters when raw materials are exchanged or when operational parameters are varied. Furthermore the simulations give advantageous information for the design of green field plants and allow the simulation of different plant configurations.The objective of this paper is to demonstrate what the latest developments are in these numerical simulations and what Siemens VAI is doing to overcome the remaining challenges in this field.

Aug 2, 2018