Search Documents

A digital computer program has been developed to mathematically simulate the effect of limited changes in feed rate on the grinding circuit of an operating mill. The data on which the program is based was supplied by the staff of the mill in question. Samples taken around the complete grinding circuit were analyzed for solids content and particle size distribution for two feed rates. Matrix models were used to represent the grinding units: rod mill and pebble mills. Parameters for these models were computed using a method developed by Mular(l). Classifier parameters were computed by the Powell search method(2) for the Plitt model (3) . The simulation program, which permits the computing of the size distributions in the circuit resulting from the changes in feed rate between the two measured values, assumes that a linear relationship exists between the corresponding parameter coefficients. The program uses an iteration method to bring the simulated grinding circuit to a steady-state material balance for any specified' feed rate between those measured. The program is, in addition, extremely fast and economic in computing time. The results show that the simulations are not everywhere consistent with practice. For example, higher feed rates on simulation predict finer products. These errors are shown to be attributable to anomalies in data gathering techniques. An improved sampling procedure for grinding-circuit data gathering is suggested for future simulations of this type. It is concluded that the

Jan 1, 1971

By R. F. Pilgrim

A digital computer program to mathematically simulate the effect of limited changes in feed rate in the range of 80 to 95 tph on a plant-scale wet grinding circuit has been developed and tested. The circuit studied included an 8-ft-diam X 12-ft rod mill operating in open circuit and feeding two parallel 12-ft-diam x 10-ft pebble mills, each one operating in closed circuit with an 8 X 32-ft rake classifier. The simulation is an empirical one based on two sets of data collected manually from the operating circuit under what was said to be steady-state feed rates of 80 and 95 tph with all other operating variables also at steady state. The procedure is a combination of several published approaches to the modeling of grinding mills and classifiers and it is applicable to a wide range of circuits. The program permits the computing of the size distributions of all feeds and discharges of all grinding units and classifiers in the circuit for a specified feed rate. The data used in this report for demonstrating the program are unsatisfactory since the circuit was not operating under steady-state conditions when it was sampled. Consequently the circuit model parameters could not be accurately estimated. The predictive power would be greatly improved if new data were obtained under true steady-state conditions.

Jan 1, 1973

By J. Ramirez-Castro, J. A. Finch

"Galena is the finest component of the product of a closed-circuit grind treating a lead-zinc ore, largely due to preferential recycling of this high-density mineral by the cyclone. For example, at the Pine Point Mines concentrator, for ~ 55% -200 mesh over-all, the galena will be close to 95% -200 mesh. This sliming of the lead can limit recovery, selectivity and dewatering. A proposed circuit modification to reduce the sliming is to incorporate a unit cell. A simulator was constructed to study the No. 3 grinding circuit at Pine Point modified by a unit cell on ball mill discharge (BMD).The simulator used a cumulative size basis rate-of-breakage with plug flow ball mill model, the parameters being obtained by laboratory grinding tests. Flotation and the cyclone models used a cumulative recovery/size matrix. The behaviour of each mineral, galena, sphalerite, pyrite and calcite/dolomite, was simulated. Simulation was adapted to a programmable calculator. A close fit to the existing circuit performance at two tonnages and mineral mixes was established.The simulator showed that the modified circuit could significantly coarsen the Pb grind, from 95% to —50% -200 mesh, depending on recovery, without substantially altering the size of the other minerals. The other principal findings were: compared to laboratory flotation tests on existing circuit BMD, the unit cell concentrate (CONC) will be coarser and markedly lower in grade; relatively low unit cell stage recoveries are desirable, because the cyclone scavenges the tail to give a circuit recovery which is higher than the unit cell stage recovery; unit cell stage recovery should not exceed the cyclone stage recovery, as the latter is the maximum circuit recovery; low unit cell stage recovery helps offset the disadvantage of the high flow at BMD; other high-density minerals will also be scavenged into the CONC, as illustrated by pyrite, making CONC quality difficult to maintain; density control of BMD flotation does not greatly alter feed density to the cyclone or ball mill.Possible plant metallurgy is forecast for a CONC diverted to final lead concentrate."

Jan 1, 1980

By Sukumar Bandopadhyay, Arunapuram Sundararajan

"For the last three decades, the application of network planning has provided a basis for a formal and general approach to the discipline of project management. For fixed activity duration, a very simple algorithm commonly known as critical path method (CPM), gives the length of time required for the total project as well as indicates the activities which are vital to the completion of a project. Project Evaluation Review Techniques (PER T) on the other hand, generalize the above approach to recognize uncertainty in the activity durations by allowing them to be random variables. One of the misleading aspects of current "" PER T"" solution method is, however, the implication that there is a unique critical path. In general, any number of paths could be critical, depending on the particular realization of the random activity durations that actually occur. In this case, simulation can be utilized to estimate completion time of a project.In this paper, the concept of criticality indexing is discussed, and a simulation model application which defines the critical path based on the criticality indexing concept, to a longwall development - extraction network, is presented. The model also provides suitable scheduling of resources based on the resource constraints and smooth s out period-t o-period assignments. IntroductionThe successful management of large-scale projects require careful planning, scheduling , and coordinating numerous interrelated activities. To aid in these task s, formal procedures based on the networks and network techniques (variously called PERT or CPM) were developed. In its original form the completion time for individual activities of a project was assumed fixed and known in advance: this in turn led to the familiar algorithms for computing the critical path, earliest start time, slack, and so forth. The unreality of such an assumption is, of course , apparent in many contexts. Consequently, attempts were made to introduce probalistic activity completion times, allowing representation of the stochastic nature of most projects."

Jan 1, 1987

By A. Haghighat, B. Lattimer

"Mine face ignition due to a high levels of methane is one of the most common fire incidents in underground coal mines. This work exams ignition scenarios in continuous miner headings. Depending on the magnitude of the ignition or resulting fire, the ventilation curtain in the continuous miner (CM) heading may be damaged or dislodged, affecting the ventilation into the area. A computational study was conducted to investigate the effects of different levels of damage to the ventilation curtain on the ventilation flow into the CM area. The computational fluid dynamics software Fire Dynamics Simulator (FDS), Version 6.0 was used to predict the conditions that develop due to a 200 kW fire in the coal mine face following a methane ignition with different levels of curtain damage. Smoke layer depth, temperatures and mass flow into the CM area were used to evaluate the impact of no damage to the curtain, partially dislodged curtain, and fully removed of the curtain. Partial removal of the curtain allowed more smoke to accumulate in the CM area and a reduction in ventilation into the CM area. With full removal of the curtain, the smoke accumulated through the entire depth of the CM area and flowed out the entrance to the CM area and the ventilation completely bypassed the region. INTRODUCTION Active working faces in underground coal mines are one the most hazardous locations in terms of fire events. According to Mine Safety and Health Administration (MSHA) data analysis, the highest number of non-fatal fire incidents was reported at mining faces, intersections and crosscut areas in comparison with other locations from 2000 to 2012. About 464 fires happened at active working faces from 2000 to 2012, most consisting of a methane ignition, in some cases followed by localized fire. These are high risk incidents due to the potential to raise float dust and damage ventilation systems, particularly auxiliary systems such as curtain or tubing. Depending on the magnitude of the ignition or fire, the ventilation curtain in the continuous miner (CM) area may be damaged or dislodged, affecting the ventilation into the area. Five computational fluid dynamics (CFD) scenarios were considered to investigate the smoke layer depth, temperature, visibility, mass flow rate, and the behavior of airflow in an entry with a continuous miner (CM) located in it. The state of the auxiliary ventilation was varied based on the assumption that as an ignition increases it magnitude it will affect ventilation curtain proportionally. Visibility and temperature along the height of the entry at different sections were investigated and the interface height at each section was determined. Also, the fresh air height was studied for tenability analysis in future. FDS, Version 6.0 was used to predict the conditions that develop due to a 200 kW fire on the coal mine face following a methane ignition with different levels of curtain damage. Smoke view software was used as the post processor to visualize and analyze methane fire behavior at the face. Most entry fires are reported as fuel-controlled fires as they are rarely restricted to air access (Ingason et al, 2014). The working face fire in the following models was assumed to be a fuel-controlled fire."

Jan 1, 2016

By K. Luxbacher, A. Haghighat, B. Lattimer

Mine face ignition due to high levels of methane is one of the most common fire incidents in underground coal mines. In continuous miner headings, depending on the magnitude of the ignition or resulting fire, auxiliary ventilation controls may be damaged or dislodged, affecting the ventilation into the area. A computational study was conducted to investigate the effects of different levels of damage to an exhausting ventilation curtain in a continuous miner heading, and the impact on ventilation to the face. The computational fluid dynamics software Fire Dynamics Simulator (FDS) Version 6.0 was used to predict the conditions that develop due to a 200-kW fire in the face following a methane ignition with different levels of assumed ventilation curtain damage. Smoke layer depth, temperature and mass flow inby the last open crosscut were used to evaluate the impact of no damage to the curtain, partially dislodged curtain, and fully removed curtain. As more of the ventilation curtain was damaged, smoke depth and temperature outby the fire increased. With the full removal of the curtain, the smoke accumulated from the last open crosscut to the face, with the ventilation completely bypassing the continuous miner region.

Jan 1, 2016

By J. M. Cabal, M. Chica Olmo, J. P. Laille

One of the main goals of a simulation of a coal deposit mining is to predict a series of parameters considered by the miner to be useful, e.g. ash content, sulphur content, seam thickness etc). The idea is to simulate different ways of mining a multi-seam brown coal deposit (the Mine Puentes de Garcia Rodriguez deposit, ENDESA, Spain ), which is exploited by several bucket-wheel-excavators (BWE) . Two mining models are used: a geostatistical model which provides a numerical way of re- presenting the three-dimensional ore-body, then a deterministic model, consisting of simulating the technical aspects of the mining. The main purpose of this simulation is to show how the mining operations can be reproduced using these two techniques.

Jan 1, 1983

By Sukumar Bandopadhyay, Arunapuram Sundararajan

"Unique production and environmental problems are associated with multiple seam mining where dragfines are used for overburden and interburden removal. On the production side, the problems are decreased dragfine efficiency in chopdown operations, the general need for spoil rehandle, complications in sequencing stripping and coal handling operations, and potential stability problems in highly stacked spoil. Many variations of dragline methods have been reported in the literature to encounter these problems. Because of the physical properties and peculiar hydrology of sub-arctic region mining of multiple seams in such areas probably could not be done with the methods used in temperate regions.The properties of thawed material would make spoil hard to place in a stable manner and thus would influence the mini sequence. For instance, multiple seam mining which requires some rehandling of spoil and prolonged exposure of overburden and interburden material might be precluded by the instability of thawed spoil and therefore mining techniques have to be modified. Since the deployment and operation of a dragline in modes different from conventional schemes will affect productivity, it will be necessary to study the effects on production of changing the dragfine operational schemes. This paper presents a heuristic simulation model of multi-seam dragline operation of a sub-arctic mine that can be used to analyze the handling of overburden and detailed study of the dragfine performance."

Jan 1, 1986

By Majid Ataeepour, Ernest Y. Baafi

High level programming languages such ass C, FORTRAN and BASIC have been used in the mining industry to simulate truck-shovel systems. In addition, some special languages for simulation such as GASP, GPSS, SLAM, PROOF and SIMAN are commonly used. The most recent simulation system ARENA uses graphical modelling approach as well as animation notion. In this paper, the main elements of ARENA required to perform a truck- shovel simulation and to view the simulation results by means of the model animation are presented. The layout of the haulage system including shovels, trucks and dump sites was generated using draw tool in ARENA. All the components of the haulage system were represented by graphical modules. Animation tool was used to depict the performance of the haulage system as the simulation progressed.

Jan 1, 1996

By A. Exikis

Simulating a ventilation network of an underground mine can help in solving complex ventilation systems with many features, selecting proper fan specifications, location and settings. Early computer systems running simple network simulations were able to provide solutions to the most common and immediate problems regarding underground mine airflow. The more advanced systems available today can perform more sophisticated modeling and simulation of a mine?s ventilation network and present their results in a more understandable and easy to communicate form. This paper presents a ventilation network case study from an underground mine. All aspects of mine ventilation network modeling and simulation are discussed as well as the results from the simulation process. The algorithms and techniques used are thoroughly explained. Some useful guidelines for the application of ventilation network software are also provided in the text.

Jan 1, 2006

By G. G. Li

"A dynamic ventilation-modeling package, VREX, was used to simulate the mine activities through tracking the equipment movements in underground workings and determine the airflow demand for a new production scenario in a ventilation-on-demand (VOD) system. It incorporates the mine development data from a discrete event simulation and balances the ventilation network using the CANVENT ventilation solver. A case study was performed to demonstrate the capabilities of the software, evaluate the performance of the network and determine the potential energy savings for implementing VOD in a deep mine operation.INTRODUCTIONThe ventilation cost of an underground mine can be significant—up to 50% of the mine’s energy consumption could be used to run fans and provide heat/refrigeration to meet the specified production environment for workers and equipment. The VOD system has proven to be more efficient in operating the fans and supplying the amount of air corresponding to the mine activities taking place in underground workings (Hardcastle, Kocsis, & O’Connor, 2006). In order to properly solve and evaluate the VOD system, a dynamic ventilation modeling software VREX—the ventilation rules engine has been developed to simulate the mine activities and determine the airflow demand for each production scenario within an operation period (Dasys & Cameron, 2011).VREX is a decision support tool, backed by the CANVENT mine ventilation solver (Li, Kocsis, & Hardcastle, 2009), for determining the outcomes of network configuration changes over time. The software incorporates the mine development data with ventilation modeling through the discrete event simulation (DES) performed by SimMineTM. As it simulates the air requirements for different vehicles and activities, VREX can evaluate both the performance of the ventilation system and the potential economic benefits of making operational or technological changes (Pritchard, 2010). For example, the model can be used to assess the economic and environmental impacts arising from modifying regulations to support the use of electric or lower emission vehicles. The combination of DES and ventilation modeling provides a new technology to generate the production scenarios according to mine schedules and activities and to determine the potential benefits for various levels of ventilation control, such as time interval and work shift change."

Jan 1, 2019

By Pierre Roy, Claude Gosselin, Maryam Sadeghi, Michel Bourassa, Christian Rochefort, Claude Bazin, Gautier Mahieu, François Lavoie, Marilène Renaud, David Cataford

"Spirals are gravity concentrators used for the concentration of iron oxide minerals from ores. This paper presents the mineral size recovery curves observed for spirals operating in two Canadian iron ore processing plants. The size recovery curves are presented for rougher, cleaner and re-cleaner spirals. The partition curves observed for the spirals are significantly different than those observed for typical gravity concentration devices such as hydrocyclones or hydraulic classifiers. The size recovery curves show the complex behaviour of the minerals as they flow down a spiral and are analysed with a view toward process modelling and simulation. Initial simulation results show that gains in the iron recovery can be achieved by increasing the recovery of coarse iron oxide particles in the rougher spirals.INTRODUCTION Spiral concentrators are gravity based separation devices used for the preparation of coal, iron and different heavy minerals ores (Wills, 1992). The classification of the particles in spirals is based on the combination of gravity, centrifuge, drag, lift and Bagnold forces (Atasoy and Spottiswood, 1995; Burt, 1984; Bouchard, 2001). Several papers are dedicated to the analysis of the forces acting on the particles flowing down a spiral (Das and al., 2007; Sivamohan and Forssberg, 1985; Holland-Batt and Holtham, 1991; Holtham , 1992). Most of these papers aim at the development of a general model for the simulation of the classification of the particles for coal ore processing. Few papers deal with the operation of spirals for processing iron ores (Hyma and Meech, 1989; Srivastava, 2001; Dallaire et al., 1978) and fewer papers present results for spirals used as cleaners following a rougher stage. Since the classification in spiral classifiers is dependent on the particle size and specific gravity a partition or performance curve commonly applied to classifiers (Mai Manga et al., 2012) is used here to analyse spiral operation. Few papers were found to describe the operation of spirals in terms of partition curves (Atasoy and Spottiswood, 1995; Holland-Batt and Holtham, 1991; Richards et al., 2000). Most of the investigations presented the application of a model to simulate the position of coal particles of various sizes and specific gravities along the trough of a spiral (Kapur et al., 1996; Holland-Batt and Holtham, 1991). In addition most of the available literature on spirals deals with one stage of concentration by spirals (Dallaire et al., 1978; Das and al., 2007; Sivamohan and Forssberg, 1985) while industrial circuits usually consist of rougher and cleaner spirals."

Jan 1, 2014

By Hossein Aminian, Daniel Hodouin, Claude Bazin, Claude Jacob

"Solvent extraction followed by electrowinning (SXEW) is an economical option for the processing of low grade and oxidized copper ore. Phenomenological models were developed to simulate the copper solvent extraction and the electrowinning processes. The two models were linked together to make a SX/EW simulator. The simulator is used to predict the operation of the mines Gaspe pilot plant put in operation in 1996. Results confirm the potential of the simulation that would subsequently be used for student training, process optimization and to assess the performances of control strategies.INTRODUCTIONSolvent extraction followed by electrowinning (SXEW) has become in the last two decades a key process for the recovery of copper from solution obtained by leaching low grade copper ore. Several operations have been commissioned in North and South Americas in the last decade (Arbiter, 1994). In Canada, Gibraltar Mine (Scott, 1991) operates a SX plant followed by electrowinning (EW) to produce copper cathode from leached low-grade ore dumps. From 1996 to 1997 Mines Gaspe operated a pilot plant to test the economical feasibility of leaching an oxidized copper ore, followed by solvent extraction and electrowinning of copper. Despite copper SX/EW processes have the reputation of being low cost operations, plant engineers should always look at way to make the plant more efficient. Process models provide an economic platform to learn, control and ultimately optimize the operation of a plant.Several models have been independently proposed to simulate the SX process (Cadman and Hsu, 1970; Medina Gomez, 1980; Wilkinson and Ingham, 1983; Garcia et al., 1988; Aminian et al., 1997) and the EW process (Bryson, 1981; Barton and Scott, 1992; Barton and Scott, 1993). However few papers, if any, deal with the integration of the two process models into a single simulator of a SX-EW plant. This paper describes the development of an integrated simulator for the SX-EW process and shows some steady-state simulation runs of an actual SX-EW pilot plant. The paper consists of four sections. The first section shows the general layout and principles of operation of a SX-EW plant.The second section describes the SX model while the third section deals with the process model used for the copper electrowinning process. The last section shows how the two process models are linked to simulate the operation of the Gaspe pilot plant."

Jan 1, 1999

By Arthur B. Steiker

Digital simulation was used to model a proposed underground' train system. The trains are to haul ore from numerous loading points beneath the production area of the mine to a central underground dump location, and then return to the mine. The purpose of the simulation was to aid mine engineers in studying various alternatives and configurations regarding track layout, crusher (dump point) specifications, dispatch procedures, fleet size, and train length, among others. The simulation' model can also be used to study the same aspects of existing haulage systems. To accommodate a wide and initially unknown range of system designs, a modular model structure was developed with major components devoted to track routing, train movement, loading process, dumping process, and dispatch. Key statistics tabulated were number of trains dumped during a simulated week, queueing delays in the mine and at the crusher, crusher utilization, cycle time, and percent of time for each: activity (loading, moving, dumping, etc.). The simulation results were instrumental in identifying a train system design which would reduce queueing delays and meet production goals with a minimum of trains and track mileage.

Jan 1, 1982

By Chen Ziming, A. Burkert, A. L. Mular

A mechanism has been proposed for autogenous grinding that involves three distinct types of breakage: self impact breakage, media breakage and abrasion. Based on this analysis a mathematical model and an autogenous grinding simulator have been developed. Using pilot plant data fitted to the model, the simulator was tested and employed to examine both the scale-up problem and the effect of feed size distribution on performance.

Jan 1, 1995

By B. K. Mishra

The motion of the charge inside a grinding mill is of considerable practical importance. The dynamics of grinding mills have been the subject of study by a number of workers because it can accurately predict: i) the power drawn by the mill, ii) the distribution of forces inside the mill, and hence the wear rate of balls as well as the liners, and iii) mass transport of milled product through the porous ball charge. Some of the earlier analyses of the motion of a charge of particulate material in a tumbling mill were reported by White (1904), Davis (1919), and Dyer (1945) and have been reviewed by Rose and Sullivan (1958). In these works, frictional force at the wall of the mill and between the colliding balls, and the interaction of different ball trajectories, was not considered. Moreover, it was assumed that the profile of the tumbling charge is independent of the degree of filling, which is generally not true. However, all this work led to the identification of so-called "equilibrium surface" which is the plane that separates the rotating charge into two regions. The inner region is static; here there is no motion relative to the mill shell. The other region is a shear region where either the balls flow down the surface of the charge (cascading motion), or are projected out of the surface and follow approximately parabolic paths (cataracting motion).

Jan 1, 1990

By D. Vavadakis

Abandoned bauxite surface mines in the Parnassos area may be rehabilitated using bauxite residues as filler. The objective of this study was to estimate the percentage of precipitation that may infiltrate through the bauxite residues to the lower geological formations when they are disposed off in mine open pits. Two scenarios were investigated simulating: a) available field tests, and b) a demonstration test that was planned to be performed in an abandoned open pit close to the Aluminium of Greece plant. It was indicated that the annual infiltration from bauxite residues disposed off in the demonstration test for the most probable scenario is very low, equal to 13.17mm corresponding to 3.1% of the annual precipitation. Furthermore, it was estimated that the annual percolation from the field tests is approximately 4.4%, very close to the infiltration values measured in the field tests for the same period.

Jan 1, 2006

By Navid Bahrani

The grain-based model (GBM) has resolved the difficulty of matching both the directtension (sc) and unconfined-compressive (st strengths of typical brittle rock by providing a synthetic material that mimics deformable, breakable polygonal grains cemented along their adjoining sides with deformable, breakable interfaces modeled as smooth-joint contacts. In this paper, a GBM utilizing a parallel-bonded base material is calibrated to match both st and sc of Lac du Bonnet granite and then used to simulate a Brazilian test. The Brazilian strength (sB) is slightly overestimated and the failure mode includes ductile post-peak response that is not representative of what occurs in laboratory Brazilian tests. By reducing the grain tensile strength, grain failure occurs during the Brazilian test and a more realistic failure mode that includes brittle post-peak response is obtained. However, the unconfinedcompressive strength of the modified material is too low. Sensitivity analysis on the grain and grainboundary strength properties demonstrates that the grain tensile strength controls the sB/sc ratio. Because the grain strength is controlled by the parallel-bond properties and the lowest achievable st/sc ratio in a parallel-bonded material is ~0.25, it is concluded thatsc, st and sB for a typical brittle rock cannot be matched with a GBM that utilizes a parallel-bonded base material.

Jun 1, 2012

By Tomohide Takeyama, Yuki Yoshida, Masaki Kitazume

"Abstract A series of centrifuge model tests was conducted to investigate the failure pattern and stability of cement treated clay ground having an unsupported vertical slope. The model ground was subjected to monotonically increasing centrifugal acceleration until it failed due to the enhanced self-weight. The test results were compared with the numerical analyses in which the treated soil was simulated by the Smoothed Particle Hydrodynamics (SPH) method. In the simulation, the Dracker- Prager model is employed as the constitutive model. In order to simulate the brittle behavior of the cement treated soil, it is assumed that the unconfined strength is gradually reduced to a certain value which corresponds to the residual strength due to accumulated plastic strain. It was found that the simulation by using SPH method can express the brittle failure observed in the centrifuge model tests.IntroductionBecause of economical and environmental concerns it has become more difficult in recent years to obtain enough soil to construct reclaimed land for airports, electric power plants, manufacture plants and so on. On the other hand, it has become difficult to find and construct disposal sites for waste subsoil or dredged soil. These circumstances have encourage geotechnical engineers at present to use these poor quality materials as filling materials more than they were used one or two decades ago. Several soil admixture methods have been developed to improve their properties and to manufacture high quality fill materials.The deformation and strength characteristics of cement treated soils have been investigated in Japan, which have revealed that the treated soil is characterized by a brittle behavior with quite high compressive strength and quite small strain at failure, and the treated soil with relatively high strength shows tensile failure instead of shear failure. Regarding to evaluate the stability of structures adjacent to or buried in the treated ground, Kitazume (1998) and Kitazume et al. (2003) investigated the earth pressure of cement treated ground by the centrifuge model tests. And they emphasized that the failure pattern of the cement treated ground was much different from that of ordinary ground. After the Hyogoken- Nambu earthquake in 1995, many design codes for the seismic stability were revised in Japan, so that seismic stability must be evaluated more precisely. The dynamic earth pressure of ordinary grounds is usually calculated by a pseudo-static limit equilibrium method based on the Mononobe- Okabe’s method (Okabe, 1924; Mononobe and Matsuo, 1929). This calculation method is also applied to the cement treated ground, although its applicability has not been thoroughly clarified."

Jan 1, 2014

By Wilhelm Brandstätter, Johannes Steiner, Alois Triessnig, Christian Redl

"The results of numerical simulations of bubble detachment from an orifice in a horizontal plate submerged in liquid are presented. The bubbles are created by forcing gas through this orifice into stagnant liquid. The simulations are performed with a two-dimensional incompressible Lattice Boltzmann (LB) scheme. In previous works the same problem was investigated by two-phase simulations in which surface tension effects were considered by an interaction potential formulation between the phases. In contrast, in the present work a volume of fluid (VOF) model consisting of a carrier fluid and a passively advected scalar function for the calculation of gravitational and surface tension forces, respectively, is used. The aims of this study are to examine the applicability of the one-phase LB method to the numerical analysis of bubble formation and to verify the theoretical and experimental predictions of the dependencies of the departure diameter on gravitational forces, surface forces and wettability effects. To this purpose the simulations were performed at different values of gravitational acceleration g, surface tension coefficient s and contact angle ß. Nonlinear fits to the data yield functional dependencies which are in good agreement with the predictions. IntroductionBubble formation behavior is an extensively studied subject, especially in chemical engineering, power industries and metallurgical processing. Knowledge of the underlying physics of bubble motion is, for example, required for the design of chemical reaction towers and agitation tanks. Another important application is the argon gas injection through the wall of the tundish nozzle into liquid steel stream to deter clogging with solid inclusions. Argon bubble motion and its effect on flow depends greatly on different parameters, e.g. gas flow rate and size of the bubbles. Thus, if these bubbles are not managed properly, they can become encapsulated upon solidification and result in undesirable impurities in the material [1, 2]."

Jan 1, 2004