Search Documents

By Khaled Morsy

The tragic incidents at Sago and Darby coal mines in 2006 pointed out the need for stronger seals to resist the blast effects of violent gob explosions. Mine Safety and Health Administration (MSHA) has published the final rule on Sealing of Abandoned Areas in April 2008. According to the final rule, the mine seals for unmonitored gob areas should withstand an instantaneous pressure of 120 psi. Reinforced concrete (RC) seal is one of the feasible options that can be used to meet the final rule. Concrete is a heterogeneous, cohesive-frictional material and exhibits complex non-linear inelastic behavior under multi-axial state of stress. The increased use of concrete as primary structural material requires the development of sophisticated material models to have an accurate prediction of the material response to a variety of loading situations. The constitutive modeling of concrete has advanced considerably in the past three decades. Choosing a reasonable concrete material model is important in order to obtain an accurate mine seal design. This paper attempts to address the effect of different concrete material models on the response of RC seal subjected to dynamic loads from methane explosion. Three concrete models, from ABAQUS library, were selected to address the effect of non-linear behavior of concrete and damages due to tensile and compressive stresses. In addition, a typical and verified concrete material model was used for conducting 3D finite element models simulating the RC seal by ABAQUS, EXPLICIT ver. 6.6.

Jan 1, 2008

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 D. P. Roberts

This investigation aimed to simulate the experiments performed by Nic Gay (Gay, N.C. 1976. Fracture growth around openings in large blocks of rock subjected to uniaxial and biaxial compression. International Journal of Rock Mechanics and Minining Sciences and Geomechical Abstracts, vol. 13. pp. 231-243) on fracture growth around openings in blocks of rock subjected to uniaxial and biaxial compression, using the Elfen discrete element code. The results of the physical experiments indicated a number of trends that were successfully replicated in the numerical simulations, including crushing ahead of the face, the formation of a sub-vertical tensile fracture at centre-span, the formation of groups of fractures originating from the crushed face zone, the formation of large corner-to-corner fractures which precede sample failure, and linking of slot corner and sample corner fractures leading to sample failure. In addition, other features were apparent that were not reported by Gay. It is believed that these took the form of conjugate shear fractures (as reported by Gay in a previous paper) and other features that may have been altered or inhibited by the presence of pre-existing structure in the physical samples. It was found that the samples, though nominally made up of the same rock type, must have had significantly differing properties. The presence of pre-existing structures may have contributed to this discrepancy. Strain measurements from the physical experiments were reasonably well matched, except where the gauge lay directly in the path of advancing fractures.

Jan 1, 2012

By D. P. Roberts

This investigation aimed to simulate the experiments performed by Nic Gay (Gay, N.C.1976. Fracture growth around openings in large blocks of rock subjected to uniaxial and biaxial compression. International journal of Rock Mechanics and Minining Sciences and Geomechical Abstracts, vol. 13. pp. 231-243) on fracture growth around openings in blocks of rock subjected to uniaxial and biaxial compression, using the Elfen discrete element code. The results of the physical experiments indicated a number of trends that were successfully replicated in the numerical simulations, including crushing ahead of the face, the formation of a sub-vertical tensile fracture at centre-span, the formation of groups of fractures originating from the crushed face zone, the formation of large corner-to-corner fractures which precede sample failure, and linking of slot corner and sample corner fractures leading to sample failure. In addition, other features were apparent that were not reported by Gay. It is believed that these took the form of conjugate shear fractures (as reported by Gay in a previous paper) and other features that may have been altered or inhibited by the presence of pre-existing structure in the physical samples. It was found that the samples, though nominally made up of the same rock type, must have had significantly differing properties. The presence of pre-existing structures may have contributed to this discrepancy. Strain measurements from the physical experiments were reasonably well matched, except where the gauge lay directly in the path of advancing fractures.

Jan 1, 2012

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

By Kumar Ranjan, Kumari A. V. Asha, G. E. Sreedhar, Banerjee R. Raja, M. Narasimha, R. Shivakumar

In this paper, numerical simulation of magnetite, medium segregation and coal partitioning have been made in a 350 mm diameter cyclone. An attempt is made to assess the nature of impairment in efficiency in relation to the increasing amount of near gravity materials in coal. The CFD model approaches the Algebraic Mixture (ASM) model with the granular option and Reynolds stress model (RSM) to resolve the air core and turbulent mixing of the medium and coal particles. The partition characteristics of the DMC were modelled using both Lagrangian & Eulerian particle tracking for particles ranging in size from 0.5 to 8 mm and in density from 1200 to 2000 kg.m-3. Lagrangian frame work based Discrete Phase model used to simulate the coal particle trajectories inside the DMC and attempted to understand the residence time of different size and density coal particles. It was observed that coal particles having density near to separation-density (having fifty percent of probability to under flow and overflow) are exhibiting increased residence time than any density particles. As expected the smaller size coal particles shows higher residence time than the coarse coal particles. CFD simulations on effect of NGM mass fraction are simulated and analyzed using ASM model including for coal and magnetite. Coal particles with high NGM content show significant effect in misplacement of coal particles towards wrong products at all relative densities. The residence time of the particles increases because of its increased interaction with near dense particles. The effect of NGM fraction on overall cyclone performance as well as detailed flow assessment in terms of residence time distribution, NGM particle concentration profiles with respect to LZVV is presented.

Jan 1, 2018

By Yuting Xue, Brijes Mishra

The relaxation test on sandstone specimens showed typical behavior in pre-failure region and stepwise behavior in the postfailure region. Overall, more significant stress relaxation occurred within the failed specimens than the intact ones. Numerical simulations were conducted with pre-defined failure plane and with Voronoi tessellation to visualize the relaxation behavior. The model with pre-defined failure plane showed the key role of failure plane and asperity in simulating the step-wise post-failure relaxation behavior. The inhomogeneous stress distribution within failed specimen and the viscous deformation of intact rocks induced high stress concentration at the asperities. The observed step-wise relaxation initiated from the failure of asperity. Furthermore, the simulations with Voronoi tessellation showed time-dependent fracture development during relaxation in post-failure region. The presence of fractures completely changed the stress distribution. Stress concentration occurred at the front area of fractures and at the interacting areas between fractures. During relaxation, fractures still developed with time and the sudden significant increase in the “damage” coincided with the acceleration of stress relaxation leading to step-wise relaxation. Finally, the results showed the possibility of using residual strength as the long-term strength of failed rock.

Jun 25, 2020

By Yang Shuxin, Guo Qiliang, He Ping, Lu Lizhi, Yao Rui

"We combined stepwise regression and constrained regression statistical method to obtain the best boundary conditions for the calculation of in-situ stress field in Wushaoling tunnel using finite element method. Our numerical results reveal that: the main geologic factors contributing to the initial background stress are gravity and tectonics; the principal stress values would drop in fault zones; the angle between the maximal horizontal stress and the axis of Wushaoling tunnel is acute; Wushaoling tunnel tends to become stable when this tunnel has 2:1 ratio of the height to its width; the possibility of rock-burst exists under high stress.INTRODUCTIONStress state of engineering region is an important factor affecting the stability of project. Geological disasters in the project area are closely related to their stress environment. Wushaoling tunnel which is 20,050 meters long, lying along Lanzhou-Wuwei segment of Lan-Xin Railway in China, is the key project of Lan-Wu double-track railway. The region surrounding this tunnel is under a high stress environment. Accurate identification of the stress environment of the engineering area and the interaction between geomaterials and geotechnical engineering construction holds significant scientific and practical value for engineering design, construction, and operation.Many scholars explored and discussed the present-day tectonic stress through investigations of tectonic evolution history (Ding, & Lu, 1986), crustal deformation (Zoback, 1992), focal mechanisms (Yang et al., 2003, Xie, Cui, & Zhao, 2004), and so forth. However, the above methods would only provide the approximate direction of the principal stress in a broad region, and could not give the absolute stress values. In-situ stress measurement is one of the primary means to directly gauge the current crustal stress. However, in-situ stress measurement can only be applied in small number limited by measurement conditions."

Jan 1, 2015

By B. Q. Li

A numerical study of the developing fluid flow and thermal aspects during the staning phase of aluminum continuous casting processes is presented. The deforming finite element method, coupled with a quasi-Lagrangian formulation, is used to represent the growth of an ingot and the dynamic fluid flow and thermal behavior that evolves in the growing ingot. The mathematical formulation of the problem and the numerical treatment are discussed. The algorithm solves the free moving nodes as part of total unknowns with an automatic time step control for time integration. This permits an optimal computing efficiency while still keeping the computational errors under control. A selection of the results obtained using the algorithm is presented. The results are concerned with the evolution of fluid flow, temperature and solidification phenomena during a DC round ingot casting process, as it evolves from the beginning of casting operation to a steady state. These results show that the time-varying, complex fluid flow and thermal phenomena during continuous casting processes can be very well represented by the computational algorithm. The results also illustrate that, during the casting process, ~e fluid flow and temperature fields change drastically at the initial stage but evolve slowly afterwards

Jan 1, 1993

In recent years, Computational Fluid Dynamics, CFD, has been commonly utilized in the mining industry to model the fluid flow behavior in underground mine workings. This paper uses CFD modeling to simulate the airflow behavior in underground crosscut regions, where brattice sails are used to direct the airflow into these regions. Brattice sails are cost effective ventilation control devices for temporary or permanent use in underground mining. They can be used to deflect air into the unventilated areas such as crosscut regions. Their design and installation is a fundamental issue for maintaining a sufficient supply of fresh air and achieving effective air circulation and contaminant removal. At the same time, they should have little impact on the mine ventilation system. This paper presents the results of a two-dimensional CFD model, which examines the effects of brattice length on fluid flow behavior in the crosscut regions. The results of this study will assist in understanding the ventilation air behavior and in determining the optimum size of brattice curtain (sails), which provide a highly effective contaminant removal from the unventilated mine areas. This, in turn, helps the mine ventilation designers to meet the mine safety requirements.

Jan 1, 2009

By S. Viswanathan, A. S. Sabau

"Dimensional changes between a pattern die and its corresponding investment casting part occur as a result of complex phenomena such as thermal expansion/contraction and hot deformation (elastic, plastic, and creep) during the processing of the pattern material (wax), mold material (shell) and solidifying alloy. Determining the pattern tooling dimensions is crucial to the dimensional control of the investment casting process. To date, computational methodologies have not been applied to the prediction of dimensional changes in investment castings.This paper deals with the evaluation of wax pattern dimensions, which is one important factor in determining pattern tooling dimensions in investment casting. Cerita 29-51, an industrial wax, is considered in this study. The wax pattern dimensions are determined by considering its thermo physical and thermomechanical properties, restraint of geometrical features by the metal die, and process parameters such as dwell time, platen temperature, injection pressure, and injection temperature.IntroductionThe investment casting process involves making a disposable wax pattern by injecting wax into a metal mold, building a ceramic shell mold around the wax pattern by the application of a series of ceramic coatings to the wax pattern, dewaxing the ceramic shell mold, and casting the alloy into the shell mold. Dimensional changes between the pattern tooling and its corresponding cast part occur as a result of thermal expansion, shrinkage, hot deformation, and creep of the pattern .material (wax), mold material (shell), and solidifying alloy during the processing [1]."

Jan 1, 2001

By C. Sun

In this study, numerical simulations of yield pillars employing the finite element method with time-dependent Drucker-Prager material model are conducted. The stress control mechanism of yield pillars on weak floor strata is studied. The development of pillar yield zone, the distribution of pillar stress and the development of failure zone in the weak floor strata are investigated. The study indicates that a properly selected yield pillar size may significantly reduce the peak stress induced in the weak floor strata and effectively control the development of failure zone in the floor.

Jan 1, 1999

By K. Sasaki

Diffusion phenomena in mine airways are important to trace gas and dust introduced into mine ventilation flows. Authors have carried out tracer gas measurements to check mine ventilation using concentration-time curves measured in mine airways, however dispersion mechanism to enhance effective diffusion have not been cleared. A numerical model using discrete tracer points dosed into tube flow has been presented to simulate their movements in consideration of turbulent velocity profile and fluctuation. Numerical simulations have been carried out to predict distributions of dosed points in tube flows with turbulent velocity profile, turbulent intensity profile and Reynolds stress. The numerical simulation results show that a combination of turbulent velocity profile and fluctuation creates extremely large dispersion of points in longitudinal direction and concentration curve with long tail was observed. This dispersion by velocity profile and fluctuation contributes large effective diffusion coefficient compared with original diffusion coefficient dominating in uniform velocity fields. Furthermore, the effect of Reynolds stress on the concentration-time curves is discussed.

Jan 1, 2009

Numerical simulation for microseismicity due to mining at one of the collieries in Australia was carried out and the results were compared to the observed microscismicity. 3-D Stress analysis showed stress concentration on the edges of the mined areas and on the coal pillars- The stress concentration ahead of the coal face was biased to the tailgate. This is similar to the observed distribution of microseismic events. Stress reduction above and below the mined area was seen. Tensile stress was observed mainly above the mined area. It was realized that compressive failure due to nearly vertical maximum principal stress and nearly horizontal minimum principal stress in retreat direction (NS) might occur ahead of the face. Tensile failure due to nearly vertical minimum principal stress above the mined area was expected. Compressive failure due to nearly horizontal maximum principal stress in EW direction and nearly vertical minimum principal stress above and below the mined area were also expected.

Jan 1, 2001

By Gang Yang, Manchao He, Chengzhang Gao, Jun Zhang, Yayun Li, Jun Yang, Shilin Hou, Ben Liu, Yajun Wang, Jianning Liu

To explore the mechanism of roof directional breaking in gob-side entry formed automatically without coal pillars, the numerical simulation was used to study the deformation movement of roof, the deformation evolution of roadway, and the stress distribution of roadway under roof directional breaking and non-roof directional breaking. Based on the results, the deformation movement of roof is divided into three stages: the goaf roof collapse stage, the basic roof deformation stage, and the basic roof stability stage. The goaf is filled with gangue due to the broken expansion characteristics of the collapsed rock under the roof directional breaking. Due to the supporting force of gangue, the deformation of the basic roof under the roof directional breaking is smaller than that under the non-roof directional breaking, which reduces the influence of the deformation of the basic roof on the deformation of the roadway roof and improves the integrity of the roadway. At the same time, roof directional breaking cuts off the stress transfer path of partial roof, which weakens the peak stress convergence on the arc side, transfers the peak stress convergence to deeper, and improves the stability of roadway. Finally, a field test was carried out to further verify the accuracy of the numerical simulation. The field test results are in good agreement with the numerical simulation results. The roadway deformation is small and meets the requirements of safety production. It is proved that roof directional breaking plays a positive role in improving the stability of the roadway.

Dec 22, 2022

By Mark S. Rayson, Jessy Lee, David O. Johnson, Matthew R. Rawls

The ability to manipulate the structure of liquid chemicals in an explosive is important in mining applications, where precision blasting is crucial for the safe, reliable, economical and eco-friendly extraction of ores. Static mixers allow explosives, emulsion-based in particular, to be custom-made as they offer precise and reproducible mixing without the need for moving mechanical parts. In this paper, we perform numerical simulations on the laminar mixing of two nitrate-based fluids in a Kenics static mixer unit for various mixer parameters. Two numerical methods were investigated in the study of the mixing pattern: solving the mass transfer equation and using particle tracking. The former introduces artificial numerical diffusion which can lead to spurious smearing of the component boundaries, a problem that can be avoided with the latter technique. We have also demonstrated that the particle tracing method is capable of qualitatively predicting the blending patterns under the conditions studied.

Jan 1, 2014

By M. Kamyab

The Deep Exploration Technology project, as part of the Australian Cooperative Research Centre, is developing technologies to drill deep boreholes for mineral exploration purposes. The aim of the initiative is to transfer drilling technologies from oil and gas into mining applications. Several modifications have to be pondered due to this technology exchange: smaller cuttings sizes because of using impregnated diamond bits; narrower annulus space; higher drilling fluid velocity and likelihood of turbulent flow regime in the annulus; and higher pressure loss. In particular the results of our lab experiments have shown that the increase in concentration of small size cuttings can significantly change the rheology of the mud when drilling in hard rocks. This, in return, could have a noticeable increase in the pressure loss encountered along the annulus space. As a continuation of the above study, numerical simulations were conducted to investigate the cuttings transport behaviour in small size annulus space. Fluid velocity, annulus size, cuttings concentration, cuttings density, and cuttings size were some of the parameters that were varied to see their effects on cuttings transportation and the pressure loss along the annulus space. The results indicate that, due to the differences between oil and gas and mining drilling, conventional oil well planning is not applicable for mineral exploration drilling and numerical simulation would help in this matter.

Aug 1, 2013

By Valerio Viti

In High Pressure Die Casting a molten metal is injected through a thin gate into the cast cavity. High injection pressures and high gate velocities create atomization phenomena which can negatively affect the final quality of the cast. In order to control atomization process it is essential to understand the atomization patterns and the two-phase flowfield that exist at the gate exit. In the present work, 2D numerical simulations of the flow of molten Magnesium through a high aspect-ratio rectangular gate are performed using a water analogue for an open and closed cavity. The numerical simulations made use of VOF-type (Eulerian) physical models and Lagrangian models. Further, a sub-grid scale model was implemented that, in conjunction with VOF-type equations can efficiently predict the general atomization pattern without the need for high-resolution grids. The numerical results were compared to experimental data for validation.

Jan 1, 2006

By Dermot M. Ross-Brown

The objectives of this generic study were to use numerical modeling techniques to determine under what conditions seismic waves generated by an earthquake might cause instability to an underground opening, or cause fracturing and joint movement that would lead to an increase in the permeability of the rock mass. In general, the study showed that moderate earthquakes (up to 0.41 g) did not cause instability of the tunnel or major fracturing of the rock mass. An artificial rock-burst tremor with high accelerations, however, was found to be amplified around the tunnel, and fracturing occurred as a result of the seismic loading in salt and granite. In shale, even moderate seismic loading resulted in tunnel collapse.

Jan 1, 1981

By S. D. Anitha Kumari, T. G. Sitharam

"The huge urbanisation of the world cities have invariably resulted in enormous requirements of development in transport and other facilities. As the surface space available for the construction of these infrastructural facilities is found to be inadequate, underground space is widely used recently for this purpose. However, to understand the behaviour of these underground structures through experiments under various loading conditions is not practically feasible. Hence numerical methods have widely been used to model these structures and predict their behaviour. The underground structures generally rest on soil or rock. Depending upon the properties of the material surrounding the tunnel, the deformation pattern and the stress contours vary. If the tunnels are placed on soft ground at a shallow depth, the presence of existing structures will add to the stresses and strains developed. In the case of densely populated cities, tunnelling operation may result in huge damages if the settlements are not within the desirable limits. This may in turn lead to the collapse of the tunnel along with adjoining structures. So the determination of stress levels associated with the yield of wall due to excavation is very important and necessary supports and reinforcement should be provided. Numerical simulations using discrete element method (DEM) and finite element method (FEM) have been done to understand the stress variation around single and twin tunnels. Studies are also done to quantify the improvement of stability on the introduction of lining around the tunnels in a weak cemented soil mass like sandstone.INTRODUCTIONMost of the life line facilities which include pipelines for electricity, water, communication, sewage, gas, roads, metro, utility tunnels etc in modern urban areas are underground structures. Infrastructure development is crucial for any city and sustainable development urges for underground constructions. Accordingly, there is already a great demand for tunnels to be constructed in difficult and crowded urban settings. This results in the fact that a large number of risk factors are associated with tunnels in these densely populated urban areas. These constraints pose threats to tunnel construction and also demands for minimal disturbance to the public and to the surrounding utilities. One of the most important factor to be considered while designing an underground structure is its stability with respect to the surrounding medium and adjoining existing structures. In this study, a comprehensive presentation of the numerical modeling of single and twin tunnels adopting two different methods is presented."

Jan 1, 2015