Search Documents

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 M. Kilfoil

Many minerals processing operations result in crushed rock being mixed with a liquid. Sometimes, the rock particles have to be dried before the next process begins and there may be a considerable distance between the equipment carrying out these processes. Bateman Engineered Technologies (BET) developed, patented and commercialized a method of simultaneous pneumatic conveying and drying. A joint project to assist with optimizing the process is under-way. Part of it is a numerical simulation of the process. The simulation is done using a dedicated program written in MATLAB, (commercial software which is widely used in engineering). It accounts for heat and mass transfer as well as particle transport. The main part of the program consists of a loop containing a series of equations which are evaluated sequentially. During any particular iteration, the properties are changed one at a time. Since the properties are changing relatively slowly, this method gives good results. The outputs of the simulation are consistent with the operational experience of BET’s engineers. The output is very comprehensive and all numerical. Since the program is currently being used for research purposes, flexibility is very important. Thus the output is converted into a spreadsheet file. Variations in the output quantities are displayed graphically using MS Excel. Some simulation could have been done using computational fluid dynamics (CFD) software. The advantages and disadvantages (technical, financial and skills requirements) of the two approaches for this particular type of simulation are compared briefly. Keywords: Pneumatic conveying, Heat and mass transfer, Computer model

Jan 1, 2003

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 Li. Jiarong, Zuojian Liang, Qingyan Xu, Baicheng Liu, Shizhong Liu, Hailong Yuan

Single crystal solidification of turbine blades is a key technology for the production of reliable and high efficiency gas turbines. The solidification process of single crystal investment castings causes a number of problems that have to be solved such as the increase of casting defects and complexity of production process. Computer modeling and simulation is helpful to predict the defects and improve the process design. In the paper, a mathematical model for three dimensional simulation of solidification and single crystal growth of investment castings was developed based on CA, mechanism of crystal growth and basic transfer equations such as heat, mass transfer equations. Many factors including constitutional undercooling, curvature undercooling and anisotropy, which had vital influences on microstructure evolution, were considered in the model. Temperature, crystal growth interface and structure defects of investment castings of single crystal superalloy were investigated at withdrawal rates of 4.5mm/min and 7mm/min. Simulated results showed that microstructure evolution and structure defects were predicted reasonably and calculated results were coincident with experimental ones. The study indicated that three dimensional solidification simulation technique is a powerful tool for understanding the fundamental of solidification and the formation of structure defects in single crystal investment castings and finally for optimizing process design.

Jan 1, 2004

By Miaoyong Zhu, Sen Luo

"A two-dimensional CAFD model was developed for simulating the solute diffusion-controlled dentritic growth in the low Peclet number regime. Mass conservation and local equilibrium at the solid/liquid interface were solved to determine the growth velocity of interface with consideration of constitutional and curvature undercooling. The modified decentered square growth proposed by Wang et al was adopted to capture the first eight neighboring cells. The comparison between the classic Lipton-Glicksman-Kurz (LGK) analytical model for steady state tip variables and the simulated values showed that the present model was capable of quantitatively determining the dendritic growth, and the simulation tests demonstrated that the effects of the mesh anisotropy in growth kinetics and crystallographic orientation were eliminated. The numerical simulations were performed to the dendritic growth ofFe-0.6%C binary alloy in undercooled melt, and the results showed that the present model could predict the dendrite morphology with artificial growth orientation in acceptable level.IntroductionDuring the solidification of metals and alloys, the dendrite grows from the substance or the undercooling melt and solute redistributes between the dendrite and the interdendritic liquid. The solidification microstructure is close related to the phase and solute evenness and has a close linkage with the properties of final production.In the past few decades, various kinds of deterministic and stochastic models have been developed to predict the evolution of microstructures in solidification, with the fast development of computer power. Among these methods, Phase Field (PF) model [1, 2] has been known as one of the most adequate deterministic models for directly simulating the dendritic growth. However, there is a vital drawback that the PF method requires considerable computation time and is limited to simulate very small domains with a few dendrites [3]. On the other hand, stochastic models, such as Monte Carlo (MC) method [4, 5] and Cellular Automaton (CA) method [6-8], have been usually applied to predict the solidification microstructure evolution. Comparing with the MC method, the dendritic growth kinetics is explicitly integrated in the CA model, while the MC model doesn't. Therefore, the CA method is more reasonable to predict the dendritic growth during the solidification of melt."

Jan 1, 2013

By Min Chen

A numerical simulation of the stress field of the solidified shell in the mold during continuous casting of peritectic steel slab was calculated with commercial software. The results showed that the maximal stress on wide face of the solidified shell was located at the corner and near middle of the slab, and it increased with increasing the width of slab. In addition, the stress was also influenced by operating parameters such as superheat degree, drawing speed and cooling intensity of the mold, and the maximal stress value and its position changed with these parameters. For the slab with section size of 3200mm×150mm, the proper drawing speed was around 1.2m/min, with cooling intensity of 5500L/min while the superheat degree was 15°C to 25°C in order to prevent longitudinal crack occurring during peritectic steel continuous casting.

Jan 1, 2009

By Huimin Lu, Lian Zhou, Qiang Li

"A three-dimensional mathematical model is developed for the regenerative rotary hearth furnace with 40 meters medium diameter. The influence of combustion switchover and outlet pressure on the fluctuation and distribution of temperature and pressure is studied with numerical simulation method. The results showed that the temperature field and flow field is suitable for requirement of industrial production; the fluctuation of the temperature field and pressure is acceptable; the influence of outlet pressure on the furnace pressure is great.IntroduceRotary hearth furnace which is originally utilized as heat treating and reheating industrial furnace, now is also used for carbon based DRl production. With the need of market and the concern of resource recycling and environmental issues, the DRl production from ore or from waste material by rotary hearth furnace attracts more and more attentions from steel industry [1,2].When rotary hearth furnace works, the feed pellets are charged on the hearth and rotate within an annular furnace chamber with the process of reduction. For the sake of rapid reduction, usually the depth of the pellet is one or two layers, accounting for less than 90% of the volume of furnace. As result, the exhaust gases carry away a lot of heat, so that energy utilization efficient is low [3]. Currently, in order to reduce the consumption of energy and improve the production efficiency, many methods have been adopted. The regenerative combustion technology is one of the most promising."

Jan 1, 2013

By Dong-dong Liu, Jian-ping Niu, Qi Liu, Shu-sen Cheng

A numerical model was developed to simulate the three-dimensional heat transfer in the belly and lower shaft region of the blast furnace wall consisting of copper steel composite stave, packing and shell. The temperature distribution of copper steel composite stave without the slag layer or refractory lining was investigated and then compared with the copper stave. The results showed that copper steel composite stave can satisfy the requirement of overheatingfree and had an analogous heat transfer behavior as copper stave but more stable cooling intensity. This model had been validated by the good agreement between the computed and measured temperature in the thermal test.

Jan 1, 2015

By Dengfu Chen, Huabiao Chen, Yunwei Huang, Mujun Long, Youguang Ma, Lintao Gui

In the current work, a three-dimensional model coupling turbulent flow, heat and solute transport was developed to investigate the solute transport and redistribution in the mold of slab continuous casting where the Reynolds number is very high and the molten steel flows strongly. The conservation equations of momentum, energy and species for a multicomponent system which includes solute element C, Si, Mn, P, S were solved with the commercial software ANSYS Fluent. The parameters used in the model are based on the actual continuous casting. The fluid flow, temperature distribution and solute elements distribution were analyzed. Results showed that negative segregation occurred near the strand surface, and the species concentration reached a peak value at the solidification front during solidification. The segregation degree of solute element S is higher than the other solute elements because of its lower equilibrium partition coefficient. The species concentration in the liquid pool is homogeneous because of the turbulent flow effect.

Jan 1, 2015

By S Servant

In the mine of this study, cemented hydraulic backfill is placed under the form of a slurry at a percentage of solid by weight varying from 50 to 65 per cent. This slurry settles and cures in the stope. During this process an important amount of water must be drained from the stope. Field observations in the case studied here show that the geometry of the drains system has a significant impact on the drainage of the stopes. The mode of filling produces a backfill showing an evident stratification typical of a sedimentation process. In situ observations on this material show that this stratification translates into a very strong heterogeneity in the physical properties of the cemented backfill. This heterogeneity should affect the hydraulic conductivity of the backfill in the stope. In order to study this aspect, finite element simulations were performed. This paper presents the results of these simulations, where the impact of this stratification on hydraulic conductivity is related to the performance of the backfill drainage system for various geometry. This preliminary study is aimed at providing a better understanding of the role of the drains in the dewatering process of the backfilled stope in the case of this mine.

Jan 1, 1998

By Servant S

In the mine of this study, cemented hydraulic backfill is placed under the form of a slurry at a percentage of solid by weight varying from 50 to 65 per cent. This slurry settles and cures in the stope. During this process an important amount of water must be drained from the scope. Field observations in the case studied here show that the geometry of the drains system has a significant impact on the drainage of the stopes. The mode of filling produces a backfill showing an evident stratification typical of a sedimentation process. In situ observations on this material show that this stratification translates into a very strong heterogeneity in the physical properties of the cemented backfill. This heterogeneity should affect the hydraulic conductivity of the backfill in the scope. In order to study this aspect, finite element simulations were performed. This paper presents the results of these simulations, where the impact of this stratification on hydraulic conductivity is related to the performance of the backfill drainage system for various geometry. This preliminary study is aimed at providing a better understanding of the role of the drains in the dewatering process of the backfilled stope in the case of this Mine.

Jan 1, 1998

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 Y. Ren, B. Yang

The magnetic field characteristics due to electromagnetic stirring in the secondary cooling zone of an extra-thick slab during continuous casting are numerically determined using Maxwell's equations, with particular emphasis on the distribution of magnetic induction intensity and electromagnetic force, as well as the influence of process parameters on the magnetic field. The orientation of the electromagnetic force is the same when the electromagnetic stirrers at the 30° and 60° posistions of the arc-shaped slab are powered in the same direction; otherwise, the orientation is opposite. The magnetic induction at the centre of the section of the stirrer's central axis at the 30° positio, as well as the current, frequency, and fitting formulae, are as follows: B = {0.075I.(1/A)-0.3} mT and B = {-1.75f.(1/Hz)+22} mT. The electromagnetic force and current fitting formula at the central point of the B-B section is F = {0.03168I2.(1/A2)-03.762I.(1/A)+326.7} N/m3, and the fitting formula of electromagnetic force and frequency at the centre point is divided into two segments. The first takes the form of a parabola: = {-21((f-4.3).(1/Hz)2)+970} N/m3, and the second segment is negative linear: F = {-1.65f.(1/Hz)+51.5} N/m3.

Apr 16, 2024

By Tingyao Liu, Lihui Han, Qing Liu, Bao Wang, Teng Yang, Yong Sheng

This paper mainly investigates the mass transfer phenomena of rare-earth concentrate from leaching process in Rushton stirred tank reactor with rotate speed of 100-600 rpm. In the research, a certain amount of sphere particles with 50-450 μm in diameter is treated as the mineral particles. Comparison of laboratory experiment data with numerical simulation is performed. The results show that the critical impeller speed for complete dispersion is 300 rpm. And the maximum value of turbulence kinetic energy appeared 15 mm away from paddle, where the mixing time is decreased 55% and 80%, respectively, when mineral particles are injected into stirred tank through one or four filling ports. Moreover, with mixing process going on, the numbers of particles near the paddle region accounts for as much as 58% of the total amount. Meanwhile, following the particle diameter increasing, the particle concentration of the upper part of stirred tank reactor is increased.

Jan 1, 2015

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 Lala Zhao, Yuemin Zhao, Guangyuan Xie, Chusheng Liu, Chenlong Duan

"Screening is an important means of size classification used on granular or particulate materials. The vibrating modes of screen separator have a great effect on the screening process. Elliptical vibration mode is a complex motion mode which gets the common characteristics of both linear and circular vibration modes. However, the fundamental motion and screen penetration of the particles on the elliptically vibrating screen deck are still poorly understood. A numerical study of the motion of particulates flow along an elliptically vibrating screen deck is presented by means of a three dimensional Discrete Element Method (3D-DEM). The motion characteristic of the particle flow is analyzed at a particle scale. The effects of vibration amplitude, throwing index, vibration direction angle and screen deck inclination angle on the screening process are discussed. The results show that the particles average velocity increases linearly with the increase of vibration amplitude, throwing index and screen deck inclination angle, but decreases linearly with the increase of vibration direction angle. The screening efficiency is highest when the vibration amplitude, throwing index, vibration direction angle and screen deck inclination angle are 4.5 mm, 3.3, 40° and 6°, respectively. This work is helpful for developing a deep understanding of particle motion on the elliptically vibrating screen and for optimizing screen separator designs."

Jan 1, 2014

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