Search Documents

By Lin-fang Shen, Shao-jun Li, Zhi-liang Wang

"The meso-structure of real soil consists of numerous particles with different sizes, while the traditional reconstruction method of quartet structure generation set (QSGS) produce relatively uniform soil particles resulting in significant discrepancy between the actual structure and reconfiguration model. In order to remedy the deficiency, the QSGS is improved by considering the influence of porosity and autocorrelative function of soil, which make the models closer to the real soil. Based on the reconfiguration model, Lattice Boltzmann method is combined in the algorithm to conduct further simulation of fluid flow in soil pores. The D2Q9 model is applied to construct the two dimensional model for simulating mesoscopic seepage field of saturated soil by setting the non-equilibrium extrapolation scheme at the inlet and outlet boundary, and the bounce-back scheme at the soil particles boundary as well as the left and right boundary. According to the case study, the proposed program is realized to simulate mesoscopic seepage field of reconstructed soil with constant inlet seepage velocity. The results show that the seepage direction of soil is along the passage with good connectivity, therefore the permeability is controlled by these passages. The seepage velocity in such channels is faster than in bad connectivity channels. Even if the channel has larger pores, the seepage velocity is still small. INTRODUCTIONUnderstanding the seepage mechanism in porous materials is significant in both theoretical and engineering problems. As a typical porous media, soil consists of solid particles and voids, which form connected channels allowing water to flow through. Such seepage in soil exists widely among geotechnical project of slope, foundation pit, and tunnel, etc. Meanwhile, a number of engineering accidents are caused by limited knowledge of seepage characters in soil (Huang, 2006; Liu, 2013). The conventional method in describing seepage features mostly based on the parameter of permeability, which is a macro statistic mean evaluating the whole connectivity of soil. Obviously, this parameter cannot exactly reflect the details of how water flows through the inner soil structure, especially the flow velocity in local mesoscopic area, hence is unable to represent the actual seepage properties. Therefore, it is necessary to develop mesoscopic approaches to fully study the seepage mechanisms."

Jan 1, 2015

By Hang Zhao, Xinping Li, Pingan Wu, Yi Luo

The deep-buried rock surrounding the Jinping II Hydro-power Station was taken as the engineering background. Using model specimens of classified rock materials to simulate fractured rock mass under high in-situ stress conditions, the test about acoustic propagation law on single fracture model of specimen under different pressures and different angles were carried out. The variation of p-wave velocity in fractured rock mass and axial stress-strain relationship in the whole test procedure was also studied. Based on damage mechanics principle and the acoustic wave propagation in cracked rock mass, comparative study on damage constitutive relation of fractured rock mass under high in-situ stress in consideration of equivalent strain and equivalent energy was performed. As the results showed, with the change of axial stress level under constant confining pressure, the change in wave velocity can be divided into three stages: slow increase (small variation in amplitude), rapid reduce (large variation in amplitude ) and accelerate reduction. When the confining pressure is constant, the maximum p-wave velocity and the peak strength of specimens increases with the decrease of the crack angle. When the crack angle is constant, the maximum p-wave velocity and the peak strength of specimens increases with the increase of confining pressure.It is identified that the theoretical calculation results of constitutive relation about elasto-plastic damage model containing variable wave velocity in the fractured rock mass is consistent with the model test results.Therefore,the constitutive relation can provide reasonable suggestions for engineering practice.

Jan 1, 2015

By B. K. Gorain

"Mining of low grade deposits along with the necessity to design and operate process plants requiring high throughputs presents unique metallurgical environmental and sustainability challenges that need to be addressed to reduce project risks, improve project economics and to sustain the profitability of operations over the life-of-mine.Some of the key challenges include poor recovery of valuables and concentrate quality issues due to complex mineralogy and high ore variability, water scarcity and quality issues, high energy costs, complex environmental issues including tailings management, high capital and operating costs, worldwide constraints on resource availability, complicated regulations in emerging economies, changing demographics of mining operations with ever-changing cultural or community issues and the urgent need for a better assessment of project viability to reflect the dynamics of present mining scenario. Developing breakthrough solutions to these problems is the key to success of any new project and maintaining profitability of a mining operation. An integrated multi-disciplinary approach to problem solving is crucial along with the need for a more lateral thinking, strategic partnerships and collaboration involving the right players. This paper presents some of the key findings based on experiences from a large number of challenging projects and operations focused on solving complex metallurgical problems.IntroductionThere is a plethora of challenges in developing new ore bodies. The metallurgical complexity of these deposits has increased significantly due to the need for mining lower grade ores along with the presence of deleterious elements that have environmental and processing implications. In addition, many of the mining projects are capital intensive due to high throughput requirements to be economical along with the need for building infrastructure due to remote locations, while operating costs steadily increase due to high reagent, energy and labour costs, increased state royalties and community obligations, to name a few of the key factors. There is also scarcity of water in many regions requiring the use of sea water or brackish water, which has both environmental and metallurgical implications. The shortage of qualified technical manpower adds significantly to the woes of many mining companies. Considering these multiplicity of challenges, it has become imperative for many mining companies to be innovative in addressing these issues."

Jan 1, 2015

By J-F. Béland, N. Parson

Aluminum extrusions are widely used for automotive crash management systems due to their superior specific energy absorption capabilities, as compared to steel. The Al-Mg-Si alloy series can be readily extruded into complex multi-hollow profiles, as required in crash systems, and can be used to target a range of desired strength levels. This paper focuses on the development of an axial crush test protocol and evaluation of the crush performance of common 6XXX series alloys for a variety of tempers. The observed crush behaviour is then compared to tensile test results.

Jan 1, 2015

By R. O. Suzuki, S. Natsui, R. Nashimoto, T. Kumagai, T. Kikuchi, H. Takai

A newly developed computational fluid dynamics model based on a multi-phase particle simulation is presented for predicting the entrainment behavior of liquid metal into slag due to rising single bubble. By comparing results calculated using this model against experimental data, it was found that the transient behavior of bubbles and the two immiscible liquids can be accurately estimated by this method. The rupturing of the thin water film surrounding the bubble was less reliably predicted, but this 3-dimensional unsteady numerical model still nevertheless provides valuable new information for directly predicting the change in the liquid-liquid interface over time. This model predicted the Young-Laplace pressure, and the interface shape was simulated with volume conservation using this model as opposed to the conventional particle method. This simulation model can calculate both of dispersed phase and continuous phase seamlessly.

Jan 1, 2015

By Chris Twigge-Molecey, John Peacey

To meet the business objectives of a metallurgical project, there is often a need for new and untried technology or equipment, whose very novelty can present serious risks to meeting those objectives. The key risk management approach is a rigorous stepwise progression from theoretical evaluations and laboratory tests to pilot plant trials. These must be followed by demonstration-scale trials before final construction of a commercial scale plant, paying particular attention to the magnitude of the scale-up between phases. Success depends on following a well-defined program, not eliminating any steps, and ensuring adequate resources are available to the development team. Chapter 6 discusses the nature of technology risks to metallurgical plant projects, and the rigorous management practices needed to identify, analyze, evaluate, and avoid or mitigate them.

Jan 1, 2015

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 A. Makinen

There is an increasing trend between increasing consumption of water and the production of metal concentrate. Along with this trend we see that ore grades are declining and getting more complex to concentrate around the world. This will not only complicate processing but also generate new requirements for the management of water and waste streams on mining sites. The mining industry is a large user of water and its water consumption is growing rapidly. The processing of lower quality ores requires greater volumes of water. For these reasons, the need for efficient use of water and water management at mine sites is becoming more urgent. The perception of water in the mining industry has changed from thinking about it only in terms of "water treatment" to actually seeing it as a valuable resource, a resource that needs proper management. This resource management needs to involve various aspects, such as optimizing water usage of the entire mine or concentrator site, minimizing fresh water usage, minimizing or eliminating waste generated in water treatment and implementing practices for water reuse and water recycling. This paper presents trends and new developments relating to waste and water management in the mining industry.

Jan 1, 2015

By Chunfeng Zhao, Chong Bao, Cheng Zhao, Jiashen Tian, Matsuda Hiroshi

"Based on digital image correlation (DIC) method and self-developed code, the global strain field changing and failure properties of rock-like materials with pre-existing double flaws are experimentally studied under uniaxial compression. By theoretical analysis with linear elastic fracture mechanics (LEFM), the strain approach proves to be practicable to investigate cracking process. Thus, two types of process zone are defined to discuss the specimen’s coalescence evolution at different loading stages. Crack initiation, propagation, and coalescence is a process of zone development and nucleation. DIC strain field results are studied on a meso-level using a strain approach. In the shear coalescence mode, the SPZ (shear process zone) in a bridge area coalesces with each other, while the TPZ (tensile process zone) grows independently. In general, the paper tries to establish the link between microscopic mechanical mechanisms and macroscopic mechanical responses in flawed rock. Further research should be done to study the cracking process and its effect on rock strength in detail by the DIC method.INTRODUCTIONNatural rock contains discontinuities including pores, fractures, inclusions or other defects. The existence of these discontinuities in the rock can decrease the strength and stiffness of the rock and they are a source of initiation of new discontinuities which may in turn propagate and link with other cracks and further decrease the strength and the stiffness of the rock (Sagong & Bobet, 2002). Lots of theoretical, experimental and numerical researches have been done in studying crack initiation, propagation and coalescence in rocks (Hoek & Bieniawski, 1965; Kranz, 1983; Reyes & Einstein, 1991; Park & Bobet, 2009; Gonçalves & Einstein, 2013). But when it comes to multiple pre-flaws, the issue that how the preexisting flaws affect each other becomes complicated and uncertain. The stress field around the tips of flaw can either enhance the impact of single flaw or weaken it. According to the observations have been done before, it is found that the way and function in which these cracks propagate and coalescence depends on the geometric arrangements (Bobet, 1997; Bobet & Einstein, 1998; Wong & Chau, 2001)."

Jan 1, 2015

By Joseph Hirschi, Behrooz Abbasi, Harrold Gurley, Yoginder P. Chugh

"The research goal was to assess if analytical tool/s can be developed to quantify displacements around a geological fault as longwall face moves toward it and past it. A company was faced with crossing a 2-3 m (6-10 ft.) normal up-thrown fault in the development entries and the face while mining a 375 m (1,200 ft.) wide longwall face. A field monitoring and numerical modeling study was performed to develop data. Field measurements consisted of roof-to-floor convergence and rectangular rosettes in the immediate roof. Numerical modeling consisted of constructing a 3-D structural model using the FLAC3D software with the fault plane incorporated as ubiquitous joint. Based on analyses and field observations data, the coal company was informed of impending movements along the fault plane and likely impacts on the development entries. Additional stiff supports were installed in head gate entries around the stage loader and equipment about 60 m (200 ft.) outby of the longwall face (solid). As the longwall face approached the fault, predicted displacements were documented and manifested by severe loading on the installed supports.BACKGROUND AND STUDY OVERVIEWThis paper summarizes results of a cooperative study aimed at improving longwall mining operations in Illinois. More specifically, the study focused on: 1) In-mine measurements in development entries while mining ahead of and behind a faulted zone, 2) Making observations on stability of pillars and entries while mining through the fault, and 3) Developing simplified analytical structural models using FLAC3D to assess if predictive techniques can be developed to plan mining operations through fault zone. Available data on geometry of the fault with estimates of the cohesion and frictional parameters were used for modeling."

Jan 1, 2015

By Gheorghe Marin, Franco Chiesa, David Levasseur, Bernard Duchesne

Over the past 20 years, pressure die casting has evolved to a point where structural parts can now be produced by this process. This has been made possible by the development of vacuum die casting and by the development of die-castable, low iron content alloys. This resulted in the production of low porosity, weldable die-cast parts on which heat treatment may be applied. In the present paper, the influence porosities on the tensile properties of as-cast A360 tensile test bars was evaluated. This was done by casting 24 clusters, comprising three 6mm diameter test bars each, in a 280T cold chamber pressure die casting machine after a steady dynamic equilibrium had been reached. Modeling of filling allowed to explain the differences in porosity observed in the 3 samples of a same cluster while the results of solidification modeling were related to the microstructure. Tensile samples were radiographed on a real time X-cube apparatus and separated into 2 lots: samples where no porosity defect were visible and samples with visible defects. By measuring the tensile properties of the 2 lots it was possible to quantify the drop in tensile properties brought about by porosity, which is the most common defect observed in pressure die-cast parts.

Jan 1, 2015

By C. S. Ha, H. S. Choi, D. Y. Shin

POSCO installed a new 100 ton oxygen bottom-blown converter in 2013 for the development of new steelmaking process named as PS-BOP which is reading POSCO-Saarstahl Basic Oxygen Process. This process was developed by the collaboration research between POSCO and Saarstahl AG in Germany that is based on substantial increase of CO post combustion in a converter by blowing hot air from the top of converter. The collaboration research has been finished successfully in November 2013 with 17% decrease in hot metal ratio (HMR) when hot air was blown from the top of the oxygen bottom-blown converter. After the collaboration research the converter has been tried to be operated as oxygen top and bottom-blown converter to reduce bottom wear and to decrease production cost with increased use of low cost Fe source by the help of increased CO post combustion. The presentation includes experience on the start up of oxygen bottom-blown converter, HMR reduction trials with hot air top blowing and transition from oxygen bottom-blown to top and bottom-blown converter operation.

Jan 1, 2015

By E. Johansson, T. Siren, K. Kemppainen

"As part of the preparations for the design of an underground repository for radioactive waste at Posiva’s Olkiluoto site on the west coast of Finland, there has been considerable site investigation work. Except the extensive surface based investigations, Posiva started in 2004 the construction of an underground characterisation facility termed the ‘ONKALO’. The ONKALO has now been completely excavated to the anticipated repository depth i.e. -430–450 m. The research conducted in the ONKALO gives further information on the bedrock and groundwater conditions of the final disposal site, as well as on the impact of the construction. The ONKALO provides an excellent opportunity to investigate the rock at tunnel scale, to conduct in-situ testing in rock, to develop excavation and final disposal techniques in realistic conditions. The ONKALO has aided in collecting the data needed, supported by a Preliminary Safety Assessment, for the application of the construction license that was submitted in the end of 2012. The target is to begin disposal operations in 2022. According to current plans, the final disposal would end in 2112 and the repository would be sealed up by 2120. The testing in the ONKALO still continues several years. This paper presents the current characterisation and research activities in the ONKALO. INTRODUCTIONAfter the Finnish Government's favourable policy decision in 2001, Posiva Oy, which is an expert organisation responsible for the final disposal of spent nuclear fuel of the owners in Finland, has focused further investigations at the Olkiluoto site, western Finland (Figure 1) and has started in 2004 the construction of an underground facility termed the ‘ONKALO’. The ONKALO is not a nuclear waste repository, but more than just an underground research laboratory URL – it is an underground rock characterisation facility. The ONKALO has now been completely excavated to the anticipated repository depth i.e. -430–450 m. In addition to facilitating bedrock research, the ONKALO also provides an opportunity to develop excavation techniques and final disposal techniques in realistic conditions. Later, the ONKALO facilities can be put into use when building and using the repository. The construction methods and materials have been selected in order not to jeopardise the long-term safety of the repository site. The findings of the ONKALO and other investigations will provide the knowledge needed for an application, supported by a Preliminary Safety Assessment, to construct a final repository for spent nuclear fuel at the site. Posiva has submitted an application for the construction licence for a disposal facility in the end of 2012. The target is to begin disposal operations in 2022 and the repository would be sealed up by 2120."

Jan 1, 2015

By Xiang Li, Xibing Li

It is paramount that the stability of rock slopes be correctly accounted for under diverse operating and environmental conditions. Owing to different sources of uncertainty entailed in such stability problems, conventional deterministic approach through the calculation of a factor of safety (FOS) cannot reflect quantified uncertainties explicitly and sufficiently. Thus it is highly desirable to invoke the reliability approach to consider uncertainties in a systematic manner. Nowadays, the current design standards have an increasing tendency to ensure rock slope stability by pre-defining a target reliability index (or target probability of failure), viz., the reliability level of a rock slope is pre-set as a target to be satisfied. In such a situation and considering the fact that the FOS used as an indicator of slope performance still has a very wide acceptance, we focus on proposing an inverse-reliability-based methodology which aims at developing a two-fold indicator strategy, i.e., using the FOS and reliability together to deal with slope stability problems. The basic idea behind it is to utilize a so-called inverse first-order reliability method (IFORM) to directly calibrate the FOS provided the target reliability index has been given. In this way, the calibration scheme relating the FOS can be rationally done by not only incorporating uncertainties, but also satisfying a desired level of the prescribed reliability. Based on this, to make the IFORM computationally feasible and also reduce its computational cost, we introduce the response surface method (RSM) to generate a surrogate for the FOS-based limit state function (LSF), considering the lengthy and complicated form of the FOS in some recognized methods (e.g., within the framework of limit equilibrium). Distinguished from the existing RSMs that potential inadequacies may arise, we present an RSM based on a hybrid approach combining support vector machine (SVM; used as a regression technique) and uniform design (UD; pertaining to experimental design). Using a representative example of rock slope stability problem, we demonstrate the feasibility and validity of the developed methodology, during which the associated sensitivity study and further discussion are also performed. By dint of such a methodology, we can establish a link between the conventional deterministic approach and the reliability approach, where a logical and realistic framework is provided, for cases achieving the specified reliability, to assist in selecting the FOS which remains adequate for rock slope stability even in the face of uncertainty.

Jan 1, 2015

By Randy McMeekin, Chris Twigge-Molecey, J. Peter Blake

The business objectives of a metallurgical project determine targets and constraints, such as the range of products, production capacities, in-service dates, and capital and operating cost parameters. Permitting is usually left to the project team, as are decisions regarding site location, technology, and construction approach. Success requires an orderly methodology, represented in this book by front end loading (FEL) phases: FEL1 examines several alternative approaches; FEL2 ranks and eliminates less economically or technically attractive options; and FEL3 develops the preferred option to the point of major permit completion, with defined scope, schedule, and cost to submit for approval to proceed with implementation in FEL4. Working through FEL1, FEL2, and FEL3, the project definition improves, estimates become more accurate, and project risks are reduced. The project design basis at the end of FEL3 is the control document for project implementation in FEL4.

Jan 1, 2015

By J. P. Harrison, K. Gao

It is well known that in situ stress measurements at a given site will generally display significant variability. However, it appears that a robust analysis of such measurements – i.e., one that is faithful to the tensorial nature of the data – is not available. Despite the fact that stress is a tensor, stress magnitude and orientation have long been analyzed separately using classical statistics. In earlier work we demonstrated the applicability of tensor statistics to the analysis of stress variability. Here, we investigate further the difference between tensor and classical statistics with respect to the stress tensor correlation matrix. Our analysis commences with the selection of a series of valid correlation matrices. Then, assuming that standardized tensor components follow the multivariate standard normal distribution, these correlation matrices are used to generate marginal probability distributions of the eigen-parameters (i.e., principal stress magnitude and rotation angle). From these, random stress tensors are generated, their eigenparameters determined, and the probability density distributions of these obtained using both classical and tensor statistics. Our results indicate that, for all correlation matrices, both methods give practically identical distributions for principal stress magnitude. However, the discrepancy between the two methods with regard to probability density distribution of rotation angle is highly dependent on stress component correlation. This discrepancy is usefully quantified in terms of the integral, with respect to orientation, of the difference between the two probability density distributions. This parameter, which we term Area Difference, is closely linked to the magnitude of the correlation matrix determinant, which itself is related to the normal stress correlation coefficient (i.e., ?x ,?y ??). Assuming that Area Difference should not exceed 0.1, we identify the critical determinant for any particular value of normal stress component correlation coefficient. Thus, using published in situ stress data and the obtained critical determinants, we show that for a specific normal stress component coefficient, when the correlation matrix determinant is less than its critical determinant, the more precise tensor statistics should be used to calculate the distribution of rotation angle. These results demonstrate that the customary analysis of stress measurement data may lead to erroneous results, and that tensor statistics appears to be preferred approach.

Jan 1, 2015

By Rong. G., T. Q. Lin, W. Du

Based on analysis of the evolutionary history of the valley in Baihetan, the mechanism of asymmetrical evolution, gradual migration-incision evolved from left bank to the right, is revealed in Baihetan river valleys. The process of in-situ stress evolution is carried out by using the Universal Distinct Element Code (UDEC). Consequently, the numerical results which indicate the main regulations of that stress field are as follow: (1) The phenomenon termed "high in-situ stress fold" appears in the bottom of the valley. With the maximum principal stress being 28-34 MPa, the stress increases significantly at a vertical depth of approximately 130 m. (2) The in-situ stress presents zonal distribution along the bank slope. It is apparent that the in-situ stress field consists of stress relaxation zone (SRZ), stress concentration zone (SCZ) and original stress zone (OSZ). The horizontal depth of stress relaxation zone varies from the left bank (roughly 200 m) to the right (about 120 m). Furthermore, the stress evolution of the left bank slope is studied. As a result, it is found that asymmetrical erosion contributes to the difference of the in-situ stress distribution as a critical factor, between the left bank and the right. And for an asymmetrical valley, the evolutionary process of in-situ stress field is a circular procedure, which includes the incision of valley, the accumulation of landslides and the migration of valley.

Jan 1, 2015

By F. S. Su, P. Z. Pan, X. T. Feng, F. Yan, W. W. Ji

"This paper aims to simulate the spalling failure behavior in hard rock based on a self-developed numerical model, i.e. a rock discontinuous cellular automaton (RDCA). A crack initiation and propagation criterion is proposed based on the plastic strain, which is used to link the continuity and discontinuity during the failure process. The discontinuities of rock, such as cracks and other internal boundaries, are represented by discontinuous displacement functions, i.e. tip enrichment and step enrichment. The partition of unity is used to improve integral precision of the elements intersected internal boundaries, which is traced by a level set method. Arbitrary crack geometries can be conveniently represented and it does not require any remeshing in the modeling of crack initiation, propagation and coalescence. Numerical modeling of excavation in hard rock using RDCA is conducted. The spalling failure process of hard rock induced by excavation is well simulated. The depth and shape of hard rock spalling failure are good agreement with in situ observation.INTRODUCTIONThe excavation in rock mass is more and more popular. At the same time, the disasters, such as rock burst, rock spalling and falling etc, induced by excavation are quite prominent. Numerical modeling is an important way to help us understand the failure mechanism of rocks. Many numerical methods have been used to simulate the excavation of rock mass. For example, the finite element method (FEM) (Karakus and Fowell, 2003), the discrete element method (DEM)(Cundall and Strack, 1979, Labra et al., 2008), Finite difference method (FDM)(Zhu et al., 2009), Element Free Galerkin Method (Belytschko et al., 1994), Discontinuous displacement analysis (DDA) (Shi, 1988) and Boundary element method (BEM) (Katsikadelis, 2002, Chen and Hong, 1992) are commonly found in the stability analysis of excavation in rocks. In the modeling of crack propagation or excavation, most of these methods are tedious since the internal boundaries have to be consistent with the element boundaries. And in many cases, the remesh has to be conducted. In recent years, people are interested in particular types of FEMs, such as the extended finite-element method (XFEM), the particle discretization scheme finite-element method(PDS-FEM(Hori et al., 2005, Oguni et al., 2009)), the combined finite-discrete element method (FEM/DEM) (Munjiza and Wiley, 2004, Cai, 2008) and spectral-element method(Gharti et al., 2012) etc."

Jan 1, 2015

By J. -F. Blais, M. -O. Simonnot, J. Mocellin, G. Mercier, J. -L. Morel

Large quantities of ferromanganese sludge are generated as waste material by blast furnace during the manufacturing of ferromanganese. These slags are very rich in manganese, zinc and lead (5 to 40 %). These residues have been deposited in a pond in periphery of the steel industry. Considering the market value of the metals, these brownfield sites can now be considered as secondary resources. We have developed a hydrometallurgical process for selectively leaching Zn, Mn and Pb in order to produce compounds of Mn and Zn pure enough to be economically recoverable and a residue rich in Pb. Batch laboratory experiments were carried out to determine appropriate leaching conditions to maximize zinc extraction (100 %) and manganese extraction (90.0 %) with sulfuric acid and to generate a residue containing 15 to 30% of lead. This Pb residue is recyclable by another process.

Jan 1, 2015

By B. Chareyre, H. Pourpak, L. Scholtès, E. Papachristos

"Numerical methods used to study hydraulic fracturing in low permeability rock masses need to take into account the interaction between the progressive failure mechanism of the intact rock matrix and the pre-existing natural fractures. A fully coupled hydro-mechanical 3D discrete elements model dealing with such interaction is presented here. The ultimate objective is to grasp the role of the rock-matrix failure on the propagation of the main fluid driven fracture under different stress fields and pre-existing fracture orientations. The proposed method provides the spatiotemporal distribution of the induced crack events which can be tracked during the simulations, making possible the identification of the local mechanisms leading to macro-failure in the rock matrix or reactivation of the discontinuity planes. Fluid flow inside a fracture, pressurized crack and hydro-fracking experimental tests were simulated as part of the validation process of the model and, are presented in this paper.INTRODUCTIONWith the expansion of shale gas and geothermal reservoirs exploitation, the technique of Hydraulic Fracturing (HF) has been widely used over the past years (McClure, 2012; Britt, 2012). Despite recent advances in micro-seismicity analysis or numerical modeling (Nagel et al., 2013), the way progressive failure of rocks affects the fluid flow depending on the in situ stress conditions or the rock heterogeneities is not fully controlled yet (Zoback, 2010).In this work we focused on setting up a three-dimensional discrete elements model to study fluid driven fractures propagation. The intact rock is modeled by a discrete element method (DEM) able to reproduce fracture initiation and propagation (Scholtès and Donzé, 2013, Garcia et al., 2013). The pre-existing fractures are explicitly modeled and can be individually defined or be part of a discrete fracture network (DFN) (Harthong et al., 2012). Their mechanical response is fully coupled with the fluid within a HM-DEM framework. The fluid flow and it’s interaction with the intact matrix and the pre-existing discontinuities is modeled through a Pore-scale Finite Volume (PFV) scheme specially enhanced for fractured rock modeling (Chareyre et al., 2012, Catalano et al., 2014) ."

Jan 1, 2015