Search Documents

By K. Baris

Methane (CH4) present in coal seams, generally known as Coalbed Methane (CBM), Coal Mine Methane (CMM) or Coal Seam Gas (CSG), is a major energy source, as it is the primary constituent of natural gas. At the same time, released methane is a potent greenhouse gas (GHG) and a major safety hazard in underground mining operations. Therefore, a considerable effort has been undertaken to develop coal bed methane extraction and utilization technologies that can achieve both environmental and safety benefits. Countries with significant coal reserves and gassy seams, i.e., USA, Australia, Canada, China, India, Poland and Germany have developed mature methodologies and practices to extract, before and during mining, and utilize significant amounts of methane from coal seams. Nevertheless, on a global scale, methane recovery remains untested and many technical and economic challenges require urgent attention. The research presented in this paper, sponsored by the United States Environmental Protection Agency (EPA), is focused on the Zonguldak Basin, the only hard coal region in Turkey. The basin is known to contain very gassy coal seams and consists of three major regions: Armutcuk, Zonguldak and Amasra. Among the three, the Amasra coalfield is considered as the most favorable for CBM recovery, as it comprises of deeper coal seams and where a private company, Hema Enerji A.S. (Hema) is already developing coalbed methane production facilities. This study aims at developing a reservoir model to predict the potential CBM recovery from the CBM-2 well, the only producing CBM well drilled by HEMA Enerji A.S in the Amasra coalfield. History matching was utilized to model the 15-day production data and uncertainty analysis was performed to estimate the optimum values for the reservoir parameters. The produced model yielded an average production rate of 143 m3/day for 10 years. The total amount of CH4 to be produced during this period was 0.53 million m3. Furthermore, the total amount of GHG to be avoided was estimated at 6585 tCO2e.

Aug 1, 2013

By N. R. Green, N. J. Humphreys, D. McBride, T. N. Croft, M. Cross, D. M. Shevchenko

"Centrifugal casting offers a route to high quality products in difficult to cast high temperature low superheat alloys and thin section molds. Under centrifugal forces metal is forced into thin sections and can fill thicknesses of less than a millimeter. However, due to the high liquid metal velocities there is a high risk of surface turbulent flow and air entrainment within the liquid metal. The combination of interacting flow-thermal-solidification phenomena and associated defects is a challenging modeling task which the authors have previously described and validated. Capturing the metal-air interface, on what are inevitably complex three dimensional geometries, results in highly computationally expensive simulations and simulating a single cast can take weeks on a single processor. This contribution reports on modeling a complex centrifugal cast, gas entrainment, bubble transport and solidification, employing meshes of up to a million elements and investigates the scalability of the model on high performance clusters.Introduction The ultimate goal of any casting procedure is to produce a defect-free product. For difficult to cast alloys, such as Titanium, centrifugal casting has a number of advantage over conventional gravity methods [1]. In the casting of thin sections, such as turbine blades, the metal is forced into the thin sections under centrifugal forces. The centrifugal force can overcome the backpressure due to surface tension and the metal can fill mould sections with thicknesses substantially less than a millimeter. However, the materials used in these castings [2] tend to be highly reactive and requires low superheat. High liquid metal velocity gives rise to a high risk of free surface turbulent flow and of associated entrainment through free surface turbulence [3] and trapping of gas present in the flow [1]."

Jan 1, 2012

By C. J. Li, J. L. Feng

In the present study we focuses on developing the poroelasticity model that deals with adsorptioninduced permeability change mostly occurred in porous media constituted several phases such as coal. The underlying pressure induced by adsorption is obtained by Virial expansion form, which is substantially founded on statistical physics. The corresponding thermodynamic pressure applying on the pore surface is then used to detail the pore pressure and further to evaluate the macroscopic volumetric strain and the porosity. In the process of physical adsorption, the adsorbed phase is in thermodynamic equilibrium with the gas phase regardless of pore size, and the corresponding chemical potentials are all equal. The relationship between surface strain and volumetric strain is derived under infinitesimal deformation conditions to gain an explicit formulation of pore surface area with respect to volumetric strain and porosity. After the formulation of surface stress modified by the Langmuir isothermal adsorption is derived, the differential micromechanical constitutive equations of porous media are obtained. By supposing that the ratio of current permeability to initial is equal to the cubic ratio of current porosity to initial porosity, the differential equations for permeability change is consequently derived. The numerical simulations by the present model for permeability change of coal induced by adsorption of CO2 and CH4, respectively, are presented under confined stress conditions. The predictions from the model are good in agreement with the experimental data of sorption-induced permeability change of coals in the literature.

Jan 1, 2015

By Tetsuo Yamazaki

Natural cold seepages are characterized as rapid upward transports of methane from deeper parts of geological structures to the seafloor. Thermogenic and/or biogenic methane generated in deep sediments moves up to seafloor, and major part of methane is consumed in anaerobic oxidation of methane (AOM) and hydrogen sulfide is created in anaerobic sulfate reduction by consortium of microorganisms. Then the hydrogen sulfide is utilized by chemosynthetic communities near the seafloor, such as Acharax, Calyptogena, Bathymodiolus, and the mat-forming bacteria Beggiatoa (Fig. 1). When the supply of methane is large or rapid, remaining methane escapes from seafloor, and it is aerobically oxidized in seawater. A mass balance ecosystem model of chemosynthetic community around natural cold seepages in sediment layer and on seafloor has been created and improved by the authors (Yamazaki et al., 2005; Takeuchi et al., 2007). The model is constructed from 3 main processes as schematically shown in Fig. 2: (1) Methane supply process; (2) Ecosystem process; (3) Water column process. The part of chemosynthetic communities is in the ecosystem process. A numerical model CANDI (Carbon And Nutrient Diagenesis) was developed by Boudreau (1996) to simulate the biogeochemical processes such as organic matter, oxidant, nutrient, bi-product, and pH diagenesis in aquatic surface sediments. Then, the improved version named C. CANDI was developed (Luff and Wallmann, 2003; Luff et al., 2004). The authors created the same type carbon and nutrient diagenesis process model with calcium carbonate formation like C. CANDI (Yamazaki et al., 2005), and have added the sulfur oxidizing and immobilization processes (Takeuchi et al., 2007). The sulfur oxidizing and immobilization processes are connected with the carbon and nutrient diagenesis process model.

Jan 1, 2008

By A. Sand

Within stamp-charged coke making a large volume of coal is compacted to one single coal cake before entering the coke oven chamber. This is done by means of several falling stampers in a stamping machine having a mold nearly of the oven?s dimensions. Producing a high quality coke also from inferior coals requires a cake density of approximately 80% of the coal density. Besides that the industrial stamp-charging process demands a minimum mechanical strength of the coal cake to move it from the stamping mold into the oven chamber without failure. Densification and the build-up of cake strength were investigated earlier in stamping tests using a micro-stamping device. The quantities derived from these tests (as e.g. cake density or porosity) represent average values for the entire cake. Statements describing the local compaction state at different heights are difficult to achieve without destructing the cake. In order to gain better insight into the densification process and the inner structure of the coal cake a computational model based on the discrete element method (DEM) has been set up using 2 and also 3-dimensional simulation software. The stamper?s position and velocity as well as the force acting on the stamper were monitored and the model?s response was compared against measurement data from laboratory stamping tests. It was possible to reproduce the force and displacement pattern of the stamper in response to the viscoelastic properties of the cake using standard DEM bonding and contact models. Furthermore, the rearrangement of particles in response to the compaction by the stamper was tracked by calculating their displacement at the point the stamper hit the coal surface. The latter can also be used as indicators of particle deformation or breakage. By defining control points at different heights the particle displacement, stress and strain rates, porosity could be studied at different heights. Keywords: coke making, coal stamping, compaction, agglomeration, DEM, simulation

Sep 1, 2012

Many investigators believe that column flotation cells offer significant advantages over standard mechanical machines for the flotation of fine particles. However, because of their unique design and operation, conventional techniques for flotation cell scale-up and design cannot be applied to columns. In an attempt to help alleviate this problems a population balance model based on first principles has been developed for fine particle flotation in a column. Two different terms have been considered in the model, i.e., transport and rate. Transport terms, incorporating fluid flow and bouyancy, are used to describe the movement of air bubbles, unattached particles and bubble-particle aggregates along the length of the column. Rate terms, which describe the bubble-particle attachment process, have been derived from first principle considerations. Because the model is based on first principles, it can be useful for the design, control, optimization and scale-up of column flotation cells.

Jan 1, 1987

By G. T. Adel

The growing interest in column flotation has created a need for a fundamental scale-up procedure which is capable of meeting the necessary design and control criteria for industrial columns. Unfortunately, most of the avail able scale-up procedures deal with column diameter only and not overall column volume. Current work focuses on the effect of bubble size distribution and air hold-up on the scale- up of column volume. A series of pressure transducers has been used to measure air hold-up along the length of a bench-scale column. These measurements have been used to validate a mathematical model of the air phase. In addition to scale-up, the model appears to have the potential for use in control of air hold-up and bubble size.

Jan 1, 1990

By Jared Ball, Darren Snashall, Paul Telford, Pengfu Tan, Ben Hogg

"The foamed converter charges were up to 3 8 per month in Xstrata Copper Smelter at Mount Isa. The issue of high converter foamover frequency had remained unsolved for many decades in Mount Isa. A thermodynamic model of copper P-S converter and a viscosity model of converter slag have been developed to simulate the slag blow and copper blow, in order to understand the mechanism of converter foamover. The model predicts the chemistry and viscosity of slag during slag blow and copper blow. The effects of oxygen potential, charge of dirty reverts, slag carryover, under-blowing, over-blowing, and temperature on the foamover and viscosity of slag in copper blow have been discussed. A theory has been developed to understand the mechanism of converter foamover. The practical solutions have also been developed and implemented in Mount Isa Copper Smelter. Some improvements of the industrial operations have been presented. After a series of improvements, the number of foamed charges per month in the converter operations has been decreased significantly from 35 in Aug 2008 to only one per month in 2009/2010.1. IntroductionThe Copper Smelter in Mount Isa receives the concentrate containing about 25% copper and produces anode copper that contains about 99% copper. There are several stages in the process. The first stage carries out in the ISASMELT™ furnace. The ISASMELT™ process produces a molten copper matte containing about 60% copper. An iron-silica slag and a sulfur dioxidebearing off-gas are the by-products of the process. The slag is taken to the slag dams for cooling before it is recycled through the Mount Isa Copper Concentrator.The matte is sent to Peirce-Smith converters. The converters further remove iron and sulphur to produce blister copper which contains about 97% copper. The anode furnaces remove any remaining sulphur and oxygen to produce anode copper that is about 99. 7% pure. This is cast into anodes and sent to the Xstrata Copper Refinery in Townsville for final refining.In the copper pyrometallurgical industrials, slag foaming is an important phenomena which occurs when the gas bubbles have stability and not rupture. Foaming greatly increases the interfacial area between matte, slag and gas and is essential for the operation of high intensity smelting processes. However, foaming or foamover is undesirable in P-S converter since it can damage the equipment and injure the operators. Figure 1 shows that the whole process of a foamover development in a converter."

Jan 1, 2011

By M. M. Botz

Cyanidation tailings which are disposed in a surface impoundment experience a loss of cyanide due to natural degradation, frequently reducing the cyanide concentration to very low levels. Quantifying cyanide losses in terms of pond geometry, local weather conditions and feed solution chemistry has been largely empirical though in many cases mining operations rely on surface impoundments to reduce cyanide to below an internally regulated concentration or an effluent limitation. To permit a more quantitative evaluation of cyanide losses in an impoundment, a computer simulation was developed to estimate the losses of free, WAD and total cyanide due to dissociation, photolysis and volatilization. Results of the model are compared against data collected for a North American tailings impoundment during 1998.

Jan 1, 1999

By M. M. Botz

Cyanidation tailings which are disposed in a surface impoundment experience a loss of cyanide due to natural degradation, frequently reducing the cyanide concentration to very low levels. Quantifying cyanide losses in terms of pond geometry, local weather conditions and feed solution chemistry has been largely empirical though in many cases mining operations rely on surface impoundments to reduce cyanide to below an internally regulated concentration or an effluent limitation. To permit a more quantitative evaluation of cyanide losses in an impoundment, a computer simulation was developed to estimate the losses of free, WAD and total cyanide due to dissociation, photolysis and volatilization. Results of the model are compared against data collected for a North American tailings impoundment during 1998.

Jan 1, 1999

By Yaming Chen, Gerry Papini

This paper presents a case study of three-dimensional groundwater modeling of a proposed open pit at a gold mine project with historic flooded underground (u/g) mine workings in British Columbia, Canada. The modelling approach for dewatering the open pit included passive dewatering by draining water under natural hydraulic gradients and active dewatering by using active dewater wells, together with residual water drains. Target water level drawdowns behind the pit wall were used to determine the number and depths of dewatering wells to meet the slope stability requirements and in estimation of pit inflow rates. Conceptualizing the u/g mine workings and assigning appropriate flow boundary conditions throughout the life of the mine posed some unique challenges.

Sep 30, 2022

By Troy Williams, Chris Preston, Ian Lipchak

"Underground blasting operations are challenging from the standpoint of the distribution of explosivesenergy representative of ring blasting. Energy from both shock and pressure regimes of commercialexplosives may appear concentrated in the collar region of a typical ring blast and diluted at the toes of holes due to the oblique geometries of blastholes. The non-homogenous nature of ore in whichexplosives are distributed via drillholes, adds to the complexities of generating particulate profiles from fragmented material with consistencies that are predictable from blast pattern to blast pattern - well suited for specific underground handling equipment and mill processing. In an ideal world, it would be the blasting operations themselves that represent the primary crushing mechanism, or at least mitigate mechanical crushing that can comprise a large component of the cost in generating suitable muck. This paper presents a dynamic break view in 3D that allows a planner to visualize the potential break zone around a blasthole generated by an explosive load using a Kleine field. Simple as it sounds, this methodology provides information that can be used in conjunction with cavity monitoring surveys (CMS) to potentially judge dilution due to overbreak as well as recovery for a typical blast. As examples, there are two break geometries that are examined regarding circular breaks and elliptical breaks around blastholes. Using a Kleine field to define break, a planner generated isosurface can be generated and compared to CMS data for calibration and prediction, using AEGIS 3D ring design software."

Jan 1, 2016

By James P. Bennett

EAF slag chemistry and its control have been recognized as important in producing quality steel at a low price for a number of decades. EAF steel producers make lime and dolomite additions to furnaces, controlling the CIS ratios and sulfur pickup from the steel. Existing phase diagram data at 1600°C for CaO, MgO, FeO, SiO2, and Al2O3 was used to write a computer model optimizing slag chemistry for slag foaming and reduced refractory wear. The model, its development, and its use in an industrial steelmaking environment will be discussed along with how the model can improve energy efficiency. The importance of EAF mass balance to account for critical materials like Si in scrap and refractory wear using the model will also be emphasized.

Jan 1, 2003

By James P. Bennett

EAF slag chemistry and its control have been recognized as important in producing quality steel at a low price for a number of decades_ EAF steel producers make lime and dolomite additions to furnaces, controlling the CIS ratios and sulfur pickup from the steel. Existing phase diagram data at 1600°C for CaO, MgO, FeO, SiO2, and Al2O3 was used to write a computer model optimizing slag chemistry for slag foaming and reduced refractory wear. The model, its development, and its use in an industrial steelmaking environment will be discussed along with how the model can improve energy efficiency. The importance of EAF mass balance to account for critical materials like Si in scrap and refractory wear using the model will also be emphasized_

Jan 1, 2003

By M Reza Afshar

In this work a mathematical model was developed to study the effect of electromagnetic force on the magnesium melt containing nonmetallic inclusion in an electrical isolated channel. The electromagnetic force exerted on the inclusion induced by DC current field was calculated. A three dimensional turbulent flow equations were solved numerically to compute the velocity field within the channel. On the basis of computed velocity and magnetic fields, the trajectories of inclusions were calculated using the equations of motion for inclusions. Parametric studies were carried out to evaluate the effect of various parameters such as current density, channel diameter and particle density on the removal efficiency of inclusions from magnesium melts.

Jan 1, 2004

By Nagy El-Kaddah, T. T. Natarajan

This paper describes a novel approach for modeling three-dimensional electromagnetically driven flow in inductively stirred systems, and its application to modeling electromagnetic and flow phenomena in sub-mold rotary stirring of billets. It is based on a new T-\jf potential formulation of electromagnetic field, which does not require the calculation of the field in free space. A finite element algorithm has been developed for solving the electromagnetic field and Navier-Stokes equations in the fluid flow domain. The computed results have shown that the electromagnetic force field generates a strong rotational flow within the vertical section covered by the stirrer, and a relatively strong secondary flow beyond the stirrer. It has also been shown that the rotational and secondary flows were driven primarily by the vorticity of the force field at the billet corners. The induced flow in the molten pool was found to be turbulent and the effective mixing region in the molten pool was about three times the length of the stirrer. The mathematical approach presented in this paper provides the theoretical framework for quantifying these effects for any given stirrer configuration and operating parameters.

Jan 1, 1999

By Evgenii Pavlov, Ivan Tronin, Viktor Timofeev, Vladimir Kislov, Sergei Bogovalov

"Modeling of the magnetohydrodynamical (MHD) and thermodynamical processes during the industrial scale Direct Chill (DC) casting of aluminum slabs is performed using a combination of ANSYS/CFX and ANSYS/EMAG software. The electromagnetic field is modeled through ANSYS/EMAG, while hydrodynamic and thermodynamic processes are modeled through ANSYS/CFX. Both packages were modified in order to calculate MHD flows in the liquid pool that forms during aluminum DC casting, and they were unified into a software combination that provides simultaneous simulation of all the above processes. The simulation shows that the electromagnetic stirring of the molten aluminum changes the shape of the solidification front (solid-liquid interface) very dramatically, namely it becomes significantly flatter. This effect should improve the internal structure of aluminum DC slabs and allows us to increase the casting velocity.IntroductionProduction of Aluminum alloys with high mechanical properties and structural uniformity is one of the main goals of industry. The preferred production method of aluminum alloys in most modern plants is Direct Chill (hereafter DC) casting because it offers high productivity and commercial success. One important limitation of this technology is the development of centerline shrinkage in slabs larger than about 1 meter. The improvement of the internal structure of long DC cast slabs can be achieved by electromagnetic stirring (hereafter ES) of the liquid metal above the solidification front. This intensifies heat transport in the zone close to the solidification front and suppresses shrinkage formation."

Jan 1, 2013

By Shawn A. Cefalu, Matthew J. M. Krane, Kent J. VanEvery

"Segregation during electroslag remelting (ESR) of a nickel-chromium-molybdenum alloy (Ni-8 wt% Cr-25 wt% Mo) has been modeled. Features of the model include species transport and AC electromagnetic effects coupled with fluid flow, heat transfer, and solidification. The effect of varying ingot diameter and melt rate, on macrosegregation patterns, local solidification time, and liquid pool profiles obtained from this model are demonstrated. The results show that the macrosegregation in the ingots correlates to the filling velocity.IntroductionElectroslag remelting (ESR) is a secondary melting process used to refine the chemical composition and the microstructure of metals. Macrosegregation of solute elements causes defects in ESR ingots, including freckles and the maldistribution of phases. These defects are not correctable through heat treatment due to the large length scales over which they occur. However, the frequency and distribution of these defects may be minimized, or even prevented, through process control based on an understanding of the conditions that lead to their formation.In the past, studies have been conducted on the effect of the process parameters on the shape of the melt pool and the thickness of the mushy zone [1-7]. In these studies, the shape of the liquid pool is determined by consideration of the heat and fluid flow induced by electromagnetic (Lorentz) force and thermal buoyancy. The transport of solute and the resultant solutal buoyancy effects were not considered in these macroscopic models. Defect formation is tied to formation of secondary phases (e.g., freckling and white spots in nickel superalloys) [5,8]. In cases where secondary phase formation is calculated (e.g., [6]), phase fractions are determined using a micro-scale model given the local solidification time (LST) and local thermal history from the macroscale model.Using the local thermal history from the macroscopic transport model to predict microstructure increases the importance of the accuracy of those results. This investigation models the macrosegregation of molybdenum during the ESR of a Ni-Cr-Mo alloy. The effect of macrosegregation on pool shape and depth is demonstrated, along with its effect on the local solidification time. By comparing the relative effects of the melt rate and ingot radius, a better understanding of the factors that affect macrosegregation, and ultimately secondary phase formation may be obtained."

Jan 1, 2004

By R. Simon

Strain-burst phenomena are generally characterised by sudden releases of energy in a volume of highly stressed rock. Strain bursts are encountered in situations such as crown pillar bursts during cut-and-fill mining, face bursts at development headings cord floor bursts in deep shaft sinking. This paper describes a simple two-dimensional linear elastic FE modeling approach that leads to the calculation of the energy-release rate (ERR) and the total energy storage rate (TESR) in the rock-mass around the openings. The merits of using such energy parameters in the analysis of strain-burst problems for identifying areas of high energy release/storage are discussed through a case stud v. However, much site-specific back analyses are needed to better understand how strain bursts and energy parameters are interrelated.

Jan 1, 1999

Strainburst phenomena, generally characterized by a sudden release of energy in a volume of highly stressed rock, are often encountered in such situations like crown pillar bursts in cut-and-fill mining, face bursts at development headings and floor bursts in deep shaft sinking. This paper describes a simple FE modeling approach which leads to the calculation of the energy release rate (ERR) and the total energy storage rate (TESR) in the rock mass around the openings. The merits of using such energy parameters in the analysis of strainburst problems to identify areas of high energy release/storage are discussed. However, much site-specific back analyses are needed to better understand how strainbursts and energy parameters are inter-related.

Jan 1, 1998