Search Documents

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

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 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 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

By K. S. Rao, R. P. Tiwari, V. K. Shrivastava

ABSRACT: In this paper numerical simulation of jointed block mass behaviour obtained in an experimental study has been presented. Singh et al. (2002) conducted extensive tests on the jointed model materials with several joint configurations. The results of this experimental study were used to validate the strength behaviour of numerical block mass model developed using professionally available software. This software basically contains different kinds of material models. From the results of numerical simulation the existence of strength anisotropy was observed. The strength anisotropy so obtained has been compared with the experimental values reported by Singh et al. (2002). The anisotropic curves in both cases are U-shaped with flat base. They show good agreement in sliding mode of failure zone in which inclination ß (= 90° - ?) of joint set-I lie between 20° to 40° for all cases of interlocking. ZUSAMMENFASSUNG: In dieser numerischen Papiersimulation des verbundenen Blockmasse Verhaltens, das in einer experimentellen Studie erhalten wird, ist dargestellt worden. Singh et al. (2002) leiteten umfangreiche Tests auf den verbundenen vorbildlichen Materialien mit einigen gemeinsamen Konfigurationen. Die Resultate dieser experimentellen Studie wurden verwendet, um das Stärke Verhalten des numerischen Blockmasse Modells zu validieren, das mit professionell vorhandener Software entwickelt wurde. Diese Software enthält im Allgemeinen unterschiedliche Arten der materiellen Modelle. Von den Resultaten der numerischen Simulation wurde das Bestehen von Stärke Anisotrophie beobachtet. Die Stärke Anisotrophie also erreicht ist mit den experimentellen Werten verglichen worden, die von Singh et al. (2002) berichtet werden. Die anisotropen Kurven in beiden Fällen sind mit flacher Unterseite u-förmig. Sie zeigen gute Vereinbarung in gleitender Erscheinungsform des Ausfalls Zone, in der Neigung ß (= 90° - ?) von der Verbindung stellenSie Lüge zwischen 20° auf 40° für alle Fälle vom Ineinandergreifen ein. RESUME: Dans cette simulation numérique de papier du comportement joint de la masse de bloc obtenu en étude expérimentale a été présenté. Singh et autres (2002) ont effectué les essais étendus sur les matériaux modèles joints avec plusieurs configurations communes. Les résultats de cette étude expérimentale ont été employés pour valider le comportement de résistance du modèle numérique de la masse de bloc développé en utilisant le logiciel professionnellement disponible. Ce logiciel contient fondamentalement différents genres de modèles matériels. Des résultats de la simulation numérique on a observé l'existence de l'anisotropie de force. L'anisotropie de force ainsi obtenu a été comparée aux valeurs expérimentales rapportées par Singh et autres (2002). Les courbes anisotropes dans les deux cas sont en U avec la base plate. Elles montrent la bonne concordance dans le mode coulissant de la zone d'échec en lequel la inclination ß (= 90° - ?) du joint placez-Je le mensonge entre 20° à 40° pour tous les cas d'enclencher.

Jan 1, 2003

By Yuting Xue, Danqing Gao, Brijes Mishra

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

By D. F. Malan, J. A. L. Napier

A number of shallow coal or hard rock mines employ pillar mining systems as a strategy for roof failure control. In certain platinum mine layouts, pillars are designed to ‘crush’ in a stable manner as they become loaded in the panel back area. The correct sizing of pillars demands some knowledge of the pillar strength and the overall layout stress distribution. It is particularly important to understand the impact of the layout geometry on the effective regional ‘stiffness’ of the rock mass around each pillar. An important design strategy is to model relatively detailed layout configurations which include a precise representation of the local pillar layout geometry and to analyse multiple mining scenarios and extraction sequences to select optimal pillar sizes and barrier pillar spacing. Although computational solution techniques are now impressive in terms of run time efficiency, a major difficulty is often encountered in assigning suitable material properties to the pillars and in devising an effective material description of the layered rock strata overlying the mine excavations. This paper outlines an efficient numerical strategy that can be used to assess large-scale pillar layout performance while retaining the ability to modify individual pillar constitutive behaviour. The proposed method is applied to selected layouts to compare estimated average pillar stress values against values determined by detailed modelling and against observed behaviour.

May 3, 2023

By Cheng Li, Yanlei Wang, Lingrui Kong, Jianxin Tang, Fang Yuan

Coal mined-out areas have a significant impact on mining safety and affect the pressure distributions of neighboring coal faces. The PFC3D software is used to study the goaf mechanical characteristics, and the core loading test verifies the rationality of the numerical simulation method. The results from the numerical modeling showed that the void ratio, particle size ratio, and particle size distribution of the rock particle sample (RPS) had a significantly remarkable influence on the compression deformation characteristics of the rock block. The parallel-bond tensile strength and parallel-bond cohesion had a notable impact on the compression deformation traits of the rock block only when the axial stress level was lower than 12 MPa. By contrast, when the axial stress exceeded 5 MPa, the effective modulus significantly influenced the compression deformation features. However, the stiffness ratio, loading rate, and parallel-bond friction angle had almost no effect on the compression deformation characteristics. Furthermore, when the stress level was less than 10 MPa, the stress–strain curves from the rock core samples’ numerical simulation and loading tests were approximately linear. In addition, the Pearson correlation coefficients of the stress–strain curves suggest that the goaf zone’s core numerical simulation results had a robust correlation with the loading test results. The findings of this study can serve as a reference for the field application of PFC3D software and safety assessments in goaves.

Jan 18, 2022

By Hwasoo Park

Keywords: FC mold, Argon, Magnetic, Solidification Fluid flow plays a significant role for steel surface quality in continuous casting. Although water models are able to model fluid flow, the flow pattern is not of practical interest since the moving solidifying shell cannot be accounted. Computer modeling has become a popular and dependable tool to analyze the complex phenomena occurring in continuous casting. In this numerical simulation we coupled momentum and energy equations considering the solidification of steel along with turbulence ?-e model. The vertical bending mold is equipped with bifurcated nozzle with 15 degrees of port angle. FC mold is supplied with DC induced magnetic fields. Influence of argon and magnetic fields on molten steel flow and heat transfer is investigated by solving momentum equations adding source terms due to multiphase ? magneto hydrodynamics. The results show that magnetic fields reduce the jet velocity from port by lower magnetic core induced field but the influence on meniscus velocity is complex to describe. Magnetic fields increase surface temperature and homogenize its distribution on surface and inside a mold.

Jan 1, 2005

By Y. M. Chen, H. L. Zhang, W. U. Bing, C. Q. Zhou

Multiphase flow in the widely used Vanyukov furnace was numerically studied. An unsteady three-dimensional and three-phase flow model was firstly built using the computational fluid dynamics (CFD) software ANSYS FLUENT®, and then solved with the volume of fluid (VOF) and k - ??model. The results showed that the proposed model could be used to predict the multiphase movement, the slag/air fluctuation, the vortex formation, and effects of structural and operational parameters. By fast Fourier transform (FFT), the dominant frequency of density with time signal was calculated as 0.29 Hz. The analysis of different injection flow rates of enriched air indicated that this variable has a major effect on the mean slag velocity. The peak mean velocity increased from 2.17 to 4.99 m/s while the flow rate of enriched air varied from 70 to 160 m/s. The proposed model provides a method to optimize the furnace structure and operating conditions for the best furnace performance and lowest energy consumption.

Jan 1, 2015