Search Documents

By B. K. Mishra

This paper focuses on the use of a numerical tool known as the discrete element method (DEM) to study the motion of ball charge in ball mills (14-18 ft). DEM is employed to simulate the motion of individual balls. While doing so, the effect of liner design on the grinding performance of the mill is studied. Two important factors such as the configuration of the liner and the geometry of the lifting surfaces are correlated with the overall mill performance. An attempt is made to show the effectiveness of the lifter bars in transmitting power from the mill shell to the interior of the mill with particular emphasis to large diameter mills.

Jan 1, 1992

The paper describes a simulating method.for ground water contamination originating from industrial and mining activrty. The finite difference model simulates heat, mass and solute transport. The applied stochastical method for dispersion provides more accurate solution of mass transport modelling. This method gives a more realistic approach of chemical processes in an aquifer depending on temperature.

Jan 1, 1988

By A. T. Iannacchione

Numerical simulation of coal pi1lar loading has traditionally been a difficult task due to the unique and highly variable properties of coal and the inability of numerical procedures to duplicate these properties. This simulation was based upon material properties developed from laboratory tests of coal cubes and examined against coal pillar strength data developed from actual underground measurements. The coal cube measurements indicated the coal followed a nonlinear Mohr- Coulomb failure criterion. The in situ measurements showed that a pillar yield zone developed similar to Wilson's hypothesis and, at high confinements, the pillar core seemed to followed a pseudo-ductile behavior. The strain-softening material model, based on a 1inear Mohr- Coulomb failure criterion, simulated coal cube and pillar behavior at moderate load conditions. However, at high loads the models were inaccurate, due to an inadequate failure mechanism and/or to changes in the pillar constraints at the coalbed interface.

Jan 1, 1989

By T. Miura, Y. Mori, H. Aoki

In a flash smelting furnace, easily combustible concentrate particles excessively oxidized and melt. The excessively oxidized molten particles then collide with less oxidized solid ones to create uniform oxidized molten particles. This collision phenomenon plays a very important role in achieving high smelting performance. The authors developed a mathematical model describing above combustion phenomena. This model incorporated fluid flow, heat and mass transfer, chemical reactions and collisions of concentrate particles. Both gas flow and particle motion were calculated using the Eulerian method. Copper concentrate was assumed to consist of chalcopyrite(CuFeS2). Considered in the concentrate reactions were the decomposition of CuFeS2, and the oxidation of the resulting sulfur(S) and pyrrhotite(FeS) to produce magnetite(Fe30 4) and sulfur dioxide(S02). Moreover a reaction due to concentrate particle collisions, in which Fe304 in the excessively oxidized molten particle reacts with FeS in the less oxidized particle, was also considered. The particle collisions were assumed to occur according to a collision probability rule obtained by computational experiments. Particle growth in diameter was described as a change in the volume fractions of particle phases. When applied to a commercial flash smelting furnace, this model reproduced the particle growth and component transfer as observed in the actual furnace.

Jan 1, 2007

By Seyed Ahmad Jenabali-Jahromi

In this paper a meshless method called Finite Point Method (FPM) is developed to simulate the solidification process of continuously cast bloom steel in both primary and secondary cooling region. The method is based on the use of a weighted least-square interpolation procedure. A transverse slice of bloom as it moves with casting speed is considered as computational domain. The two dimensional heat transfer equation together with temperature dependent thermophysical properties is solved in the computational domain. The present work is verified in the mold region by the comparison of the surface temperature simulated by FPM and finite volume method (FVM) and also comparison of solidified shell thickness simulated by FPM and measured on a breakout bloom. For secondary cooling region the validation is done by comparison of the surface temperature simulated by FPM, FVM and thermovision measurements.

Jan 1, 2009

By S. H. Seyedein

A mathematical model was developed to study fluid flow and solidification in a single-roll casting system. Fluid flow coupled with energy equations were solved numerically using a finite volume method. A body-fitted coordinate transformation technique was adopted to handle the irregularly shaped domain. In order to couple the continuity and momentum equations, a SIMPLEC algorithm was used. A fixed-grid enthalpy method was considered to model the solidification in the system. The turbulence effect on flow field was modeled using a low-Re k-e model The model developed in this work was further run under various conditions to optimize the operating parameters of the process.

Jan 1, 2004

By Murali M. Gadde, Syd S. Peng, C. T. Holland

Cutter roof failure is one of the most common ground control problems affecting the safety and economy of an underground coal mine operation. Several attempts were made in the past to understand the processes responsible for cutter development and to devise mitigating measures. Yet, many cutter failures still happen and their advance prediction remains as elusive as before. One of the reasons for the impasse in predicting cutters lies in the limitations of the current numerical modeling technologies. The development and propagation of cutters is a progressive phenomenon and continuum models can only approximate this process. Since the discontinuum modeling is severely restricted by the lack of reliable input data, alternatives must be sought to better understand cutter failures using the continuum models. This paper proposes one approach using strain-softening material model to simulate the cutter initiation and propagation process. An attempt has also been made to provide some insights into the reasons for the change in cutter location at different parts of a development layout.

Jan 1, 2005

By Sang-Joon Kim, Hae-Geon Lee, Hyun-Jin Cho

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

Jan 1, 2011

By K. Krishnapisharody

A mathematical model has been developed to predict the interface dynamics and droplet generation of the buoyant phase in a two-fluid system. The model employs the well-known SIMPLE algorithm for obtaining flow fields and the Volume of Fluid (VOF) method, for interface tracking. The implementation of the model in water-oil systems, where the oil layer is placed above a recirculating water bath, is discussed. The evolution of the water-oil interface leading to the formation of oil droplets is computed realistically.

Jan 1, 2008

By Guo Chen, Yifeng Shi, Xiaobiao Shang, Nanshang Shen, Bangqi Zhang, Jinhui Peng, Junruo Chen

"The energy is introduced directly into the volume of the honeycomb ceramic regenerator due to the complex interaction of microwaves with minerals in microwave gas heating system. A three-dimensional numerical model is developed to simulate the electromagnetic field by solving Maxwell's equations using finite-difference time-domain method. Computational result is presented for the microwave power absorption of the regenerator. Effects of the dielectric constant and the effective loss factor on the microwave power absorption in the regenerator are studied as well. Numerical results show that the heating pattern strongly depends on the dielectric property of the regenerator. The demonstration of the model of electromagnetic field distribution can be applied effectively and efficiently to optimize the fundamental parameters in the microwave heating processing.IntroductionIn recent years, microwave energy has been utilized in many fields, such as food processing, chemical engineering and mineral processing [l]. In contrast to conventional heating, where heat is transferred from a surface to the interior, microwave energy is introduced directly into the volume of minerals due to the complex interaction of the microwave with the mineral. Based on volumetric heating effects, microwave heating has several advantages such as pollution free, rapid heating technique, internal heating, high efficiency over conventional heating [ 1-4]."

Jan 1, 2013

By H. Singh, Anupam Agrawal, H. S. Grewal

"Erosion in industrial appliances causes huge loss of valuable resources especially when solid particles are present in the working medium. A number of erosion models that are available in the literature claimed to predict the erosion response. In this work, five erosion models have been evaluated using computational fluid dynamics (CFD) methodology. Three most commonly used experimental test rigs (solid particle and unsubmerged slurry jet type and pot type) have been simulated and comparisons with the experimental results have been performed. The study helps in understanding the effect of fluid medium and flow pattern on the erosion process. The knowledge acquired shows that phenomenological based erosion models give better results considering the wide range of operating parameters. The in-depth study of the flow field shows the possibility that erosion mechanism, being parameter dependent, could change with location on the same specimen.1.0 IntroductionErosion is a surface degradation phenomenon which causes the removal of material due to impacts either of solid particles and/or bubbles or gas entrained in a fluid medium. Erosion has various categories depending upon the erosion medium such as, solid particles erosion (solid particles entrained in gaseous medium), slurry erosion (solid particles entrained in liquid medium) or cavitation erosion (erosion caused by exploding gas or vapor bubbles). Humphrey [1] has listed 21 different factors that decide the erosion response of the components under a given working condition. The number of factors involved in the erosion process gives us an idea about its complex nature."

Jan 1, 2012

Numerical simulation of the fluid flow and heat transfer has been carried out for a waste-heat boiler in an Outokumpu type copper flash smelter. In this boiler, SO2-containing off-gas flows into the waste-heat boiler radiation section at about 1350°C. A fraction of cooler off-gas from the electrostatic precipitator at about 200°C was recirculated into the radiation section, and a small amount of ambient air was leaked into the boiler from the bottom. The computational results show that the flow is strongly recirculating and distinctly three dimensional in most of the boiler, particularly in the radiation section. The hot off-gas cooled down rapidly after entering the radiation section due to radiation. The predicted flow pattern and temperature distribution were in a good agreement with laboratory models and industrial measurements.

Jan 1, 1995

By Nan Wang

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

Jan 1, 2009

By Bryan J. Travis, Thomas L. Cook

INTRODUCTION Our growing concern for adequate and secure sources of energy and minerals has stimulated vigorous exploration for new sources, research toward a better understanding of geological processes, and development of new extraction technologies. The need for control, or at least prediction, of subsurface fluid flow is important for many of these technologies for example: primary, secondary and tertiary recovery of oil; ground water and waste management; in-situ fossil energy extraction (oil shale, coal, tar sands); and solution mining of uranium, copper and other minerals. These technologies, especially the last two, are characterized by highly complex systems. A partial list of the physical processes occurring would include flow in porous/fractured media, multi-phase and multi-component flow with heat transport, chemically active fluids and soils, tracers, diffusion and dispersion, fracturing and dissolution. A great deal of understanding of how these processes behave and interact can be obtained from models. MODELS Theoretical models are valuable because they: 1) provide a frame of reference for interpreting results of laboratory and field experiments; 2) once validated by experiments, allow a variety of geometries, injection/production strategies, etc., to be examined (at relatively low cost) for efficiency and stability; 3) can provide guidance to the design of experiments and field operations. Most models are based on the fundamental principles of mass, momentum and energy balance. But from this starting point, many paths can be taken. For example, there is the theory of Payatakes (1973) which concentrates on the microscale dynamics. In this approach, the rock or soil matrix is represented by a complex of characteristic channels (such as periodically constricted smooth tubes). Detail of flow within the characteristic channel is calculated very accurately. A difficulty with this model is that description of a representative channel can require several parameters which may not be easily measured. Also, it is not clear how the channels can be combined practically to model a large scale flow. Another type of model is the "global" one described by Bear (1972). Here, the continuum equations for conservation of mass, momentum and energy are averaged over a distribution of pore channels, resulting in a set of conservation equations in which the small scale structure of the medium is replaced by quantities such as porosity, permeability, dispersion coefficients and tortuosity. This approach allows computation of large scale flows. However, additional constitutive equations are needed which relate averaged quantities such as permeability to observables such as local saturation, particle size distributions, and others. An important difference between models is the way they handle the momentum equation. Payatakes' model solves the full equation. In others, it is replaced by a simpler relation such as Darcy's law (valid for slow flow rates) or by Forchheimer's equation which extends Darcy's law to higher flow rates. These simple relations can nevertheless match a great deal of experimental data (Dullien 1975). The permeability term which appears in Darcy's law and in Forchheimer's expression has been related to other quantities such as porosity and particle surface area. The Kozeny-Carmen equation is a well-known example, valid for some ranges of porosity and particle sizes and for some materials. Several other semi-theoretical, semi-empirical formulations have been devised, none of which are entirely satisfactory. One goal of researchers has been the ability to predict permeability of a porous material from basic measurable quantities such as grain size distributions without recourse to adjustable parameters. This effort has proceeded from consideration of distributions of idealized, non-intersecting channels, to intersecting channels, to consideration of both distributed pore radii and pore neck radii. This last approach has been used successfully to predict permeabilities (ranging over several orders of magnitude) in sandstones (Dullien 1975). Additionally, studies on explicit networks of channels (e.g., Fatt 1956) and percolation theory (e.g., Larson et. al. 1981) have been used to examine interconnectivity effects in porous media with the hope of eventually being able to predict permeability accurately. To be of more than theoretical interest, a model must be compared s controlled laboratory conditions where boundary and initial conditions and relevant material properties can be accurately determined. Also, extraneous forces can be eliminated so that the interactions between processes of interest can be clearly seen and flow models can be vigorously tested. In contrast to the well-defined environment and

Jan 1, 1982

By T. W. Doe, J. Noorishad

The problem of fluid flow from a well into a horizontal fracture has been studied using a nonlinear coupled hydromechanical finite element model. Such an approach is required due to the strong nonlinearity and anisotropy of deformation and fluid flow moduli introduced by the presence of the fracture. The resulting nonuniform total stress field invalidates the application of the conventional fluid flow techniques to the problem. Stable convergent solution from the above analysis are verified by a coupled steady state algorithm.

Jan 1, 1982

By Menglong Bian, Xiaotong Dong

Thermal damage is an urgent problem faced by deep mining. High-temperature rock mass, hot water and electromechanical equipment are the sources of the thermal damage. The increase in air temperature in mine tunnel is the result of fluid–solid coupling heat transfer from various heat sources. This paper takes a roadway with hot water and a mine car as the research object. Based on the Fluent, Static Structural and Steady-State Thermal modules in ANSYS, the thermal-fluid–structure coupling analysis in the roadway is carried out. It is found that the hot water has a significant effect on the thermal environment of the mine. Under the experimental parameters, the maximum heat flux density and equivalent elastic strain at the interface between the hot water and the rock are 9789.3 W/m2 and 0.00054532 m/m, respectively. The heat flux distribution in the rock area near the mine car does not change much, and the equivalent elastic strain value in this area is high. After the air is heated in the hot water area, the temperature curve changes gradually and then rises steadily. By changing the inlet air temperature parameter, it is found that the lower the inlet air temperature, the better the cooling effect, and the change of the inlet air temperature can change the air temperature field in the roadway as a whole. By changing the parameter of inlet air speed, it is found that higher wind speed brings better cooling effect, but with the increase in wind speed, the cooling effect of roadway thermal environment is very limited. When the wind speed increases from 5 m/s to 10 m/s, the air temperature of most measuring points drops by more than 2℃; when the wind speed increases from 15 m/s to 20 m/s, the air temperature of most measuring points drops below 0.7℃. After the wind speed reaches a certain value, the cooling effect of ventilation is very limited.

May 5, 2022

By T. F. Adams, L. G. Margolin

The Bedded Crack Model (BCM) is a constitutive model for brittle materials. It is based on effective modulus theory and makes use of a generalized Griffith criterion for crack growth. It is used in a solid dynamic computer code to simulate stress wave propagation and fracture in rock. A general description of the model is given and then the theoretical basis for it is presented. Some effects of finite cell size in numerical simulations are discussed. The use of the BCM is illustrated in simulations of explosive fracture of oil shale. There is generally good agreement between the calculations and data from field experiments.

Jan 1, 1982

By C H. Ryu

This paper presents a simplified two-dimensional numerical model of block fragmentation in a Jointed rock mass during the late or throw stage of a blast after wave effects have subsided. The model is based on the distinct element method. Examples illustrate the approach and allow for a comparison of effects of Joint patterns on fragmentation and resulting block motion.

Jan 1, 1986

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

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

Jan 1, 2010

By Y. Jin, H. You, H. Pang, B. Lin

Successful application of steam assisted gravity drainage (SAGD) on the recovery of land facies ultra-heavy oil in Xinjiang oilfield depends largely on shortening the preheating period with the aid of the micro-fracturing technology, which is conducted through controlled water injection to a pair of horizontal wells to create homogenously distributed microcracks in the zone between the two wells. The numerical simulation of this mechanical process is yet to be investigated. Up to date, there has been research work on rock mechanics during the SAGD production period of marine deposited ultra-heavy oilsands, but the counterpart research on land facies ultra-heavy oilsands of Xinjiang oilfield is left undone. In this regard, this study is attempted to perform numerical simulation of the fracturing process on land facies unconsolidated oilsands. With this objective, the study first investigates geological and physical identification of the field collected cores, then carries out mechanical analysis on laboratory triaixal test results, based on which numerical simulation of the fractured zone is performed. The simulation is applied through a finite element platform adopting extended Drucker-Prager constitutive model, with the parameters determined by experimental work. It is verified that hydraulic fracturing of unconsolidated oilsands generates a microcrack feature in terms of a dilative zone instead of linear apertures. The extent to which the fractured zone propagates is visually demonstrated based on the magnitude distribution of equivalent plastic strain or volumetric dilation. Evaluation of the connection of injection and production wells can possibly be determined by the expansion of the saturated region around each well. Finally, the implications on the efficiency of hydraulic fracturing in stimulation of the pay zone are investigated.

Jan 1, 2015