Search Documents

By M. Gonzalez Fernandez

"If the final goal of metallurgy is the extraction of metals from a given ore, everyone knows that this achievement relies on many separation steps involving much more than a biphasic material made of metal and gangue. Nevertheless, still today, most of the design and follow-up of industrial processes only consider metal grades; in best case metal grades per size class. Working with such heterogeneous material, made of several minerals, with different grain size and texture, implies to consider finer ore models to be able to predict the behavior of material and the subsequent operation performances. Analytical tools are now available to provide such relevant and accurate data which are applied in geometallurgy. To achieve a true geometallurgical approach for process design, but also for optimization, accurate process simulation absolutely needs to combine these fine ore models with the corresponding unit operation models. Based on generic, but configurable, criteria of description and a hierarchical approach, the material modelling tool embedded in the USIM PAC process simulator provides the necessary flexibility to describe any type of ore with almost no limit in the level of details and associated physical properties that can be met in the mining industry. The simulation engine and the mathematical models of unit operations are designed to manage such a diversity of material modelling. The paper present various ore models and some of their applications in flowsheet simulation. It proves the interest of this detailed and advanced approach to accurately predict plant performances. INTRODUCTION Process simulation used for process design and optimization requires accurate unit operation modeling. During the last decades, main efforts of the different teams developing software simulators have been dedicated to those models for comminution and physical separations such as flotation or gravity concentration (Kojovic et al., 2012, Soni and Yoon, 2014). In most case, the ore is modelled with a size distribution and metal or mineral grades. A water phase is considered for wet processes. The association of unit operations and such material models generally allows, after some parameters adjustment, to achieve good agreement between experimental and simulated data. The use of the obtained simulator to test alternative configuration of circuits or settings of equipment devices can be accurate enough as far as the ore processed is close to the ore used for simulator adjustment. However, when the ore gets very different grades, coming from different mineral proportions or grain sizes, the simulation results can be of poor accuracy. Because any physical separations such as gravity, magnetic or flotation, do not sort metal and gangue but particles, made of minerals, containing metals, the only way to be accurate enough for allowing extrapolations, is to model particles. The phase model tools embedded in the USIM PAC process simulator gives this possibility. However, as such level of details are not always available for analyses, the tool has been designed with a high level of flexibility, to adapt the model to the available data. In this paper, the principles of this conceptual modelling tool are presented, and some application cases illustrate the material model flexibility and its interest."

Jan 1, 2018

By Masakazu Niinae, Kenji Aoki

"Electrokinetic soil remediation by electroosmosis and ionic migration is restricted to soluble substances. At a higher pH, the soil particles sorb more heavy metals than at a lower pH and the precipitation of heavy metals are promoted. Both mechanisms make heavy metals immobile, rendering cleanup more difficult and electrokinetic extraction inefficient.In this study, the effectiveness of introducing a complexing agent to enhance the removal of cadmium from kaolinite in the region of neutral pH by ionic migration is examined. The complexing agent selected was ethylenediaminetetraacetic acid (EDTA). The investigation was carried out in a laboratory-scale column containing kaolinite contaminated with cadmium. As the results of the study, cadmium migrates toward the anode without accumulation of cadmium and is successfully removed at the anode reservoir."

Jan 1, 2003

By Graeme Lane, Yuqing Feng, Kaixi Jiang, Junwu Zhou, Zhengchang Shen, Tao Song

"Airlift reactors are used in a wide range of industries, such as hydrometallurgy, biochemical processes, chemical process industry and waste water treatment. With increasing superficial gas velocity three circulation regimes can be defined: no gas bubbles in the downcomer, gas bubbles remain stationary in the downcomer and gas bubbles in the circulating regime. In this paper a numerical modelling method was presented for gas-liquid flow in airlift reactors under different bubble recirculation regimes. Gas-liquid flow was modelled using the Eulerian two-fluid equations, and extra user defined subroutines were incorporated to consider the complex physics, such as bubble induced turbulence and turbulent dispersion force. Some alternative correlations of drag coefficient were used to capture the bubble distributions in the riser and downcomer of the airlift reactor, with consideration of the interaction between bubbles. As using constant bubble size is hard to describe global characteristics of the flow in airlift reactors for three circulation regimes, a model including multiple bubble sizes was presented to get more accurate simulation results of gas holdup and water velocity. Also, the inhomogeneous MUSIG model was used, as an exploration, to compute the complex flow regimes. The MUSIG model is a population balance method to predict the bubble sizes, taking break-up and coalescence into account. The modelling method was applied to a laboratory internal loop airlift reactor, and the simulation results were compared with the experimental measurements of Simcik et al. (2011) for gas holdup and water velocity. The three bubble circulation regimes were observed from simulation figures. Reasonable agreement was obtained over a range of operating conditions, and an improvement was demonstrated using the proposed methods."

Jan 1, 2014

By C. Y. Sun

Flotation of chalcopyrite and unactivated malachite with xanthates (ethylxanthate, butylxanthate, octylxanthate and dodecylxanthate) is studied in this paper. Minerals interacted with xanthates are detected by FTIR, and simulated by molecular mechanics. It shows that chalcopyrite can be floated with 4 kinds of xanthates, and unactivated malachite can?t be floated with ethylxanthate but can be floated directly with later 3 kinds of xanthates while characteristic absorption peaks of those four kinds of xanthates are all detected by FTIR on chalcopyrite and malachite surface. After simulating the interaction between xanthate and mineral, it shows that H2O and OH-can interact with chalcopyrite and malachite, and their interaction energy may be so called energy bulwark which must be overcome by collector when collector interacts on minerals; the two sulfur atoms in xanthate group interact with copper and neighbor iron atoms on chalcopyrite. It also shows that molecular mechanics results with logP of collectors can interpret the mechanism of xanthate-mineral on atom scale.

Jan 1, 2014

By Áureo Del Vecchio, Lídia Dely Alves de Sousa, José Pereira

The potassium demand for agriculture purposes has increased in Brazil, which has great deficit of internal demands. Potassium-rich rocks have been used as fertilizers as technological alternatives for agriculture aimed to reduce the amount of chemicals and the sector dependence for such products. Vermiculite tailings from União Brasileira de Mineração - UBM, at Santa Luzia - PB, was studied aiming their application in agriculture by potassium release. The sample was analyzed by XRF and XRD; studies on release kinetic with particles smaller than 2 mm added of organic acids (citric and oxalic) at 10-2 mol/L concentration; samples were analyzed by XRD after kinetic assays. The work was conducted using a shaker incubator at 300-rpm agitation and time variation to better understand the potassium availability in the ground. Results show the presence of K2O at 3.46% concentration. Citric and oxalic acids increased the release speed to almost the same level with speed coefficients (k) of 2.61x10-2 h-1/2 and 2.55x10-2 h-1/2 respectively. In statistical terms for both acids, the linear equation represents more than 99.0% of the experimental assays. Stonemeal was applied to banana seedlings in experiment conducted in a greenhouse. The use of vermiculite tailings in banana seedlings promoted potassium absorption, which influenced its growth. The amount of potassium liberated from vermiculite tailings could be used as a long-term alternative to fertilizers in agriculture cultures.

Jan 1, 2014

By Danuta Szyszka

This paper presents results of tests performed to determine minimum concentration of a surfactant used at which coalescence of air bubbles during flotation has low probability to occur. The measurements of the critical coalescence concentration, also known as CCC, were carried out in a device consisting of glass cell into which air was pumped by peristaltic pump through capillaries. The flotation foam was recorded with a high resolution digital camera. Basing on the video records, average diameters of the air bubbles were measured and the values obtained were subsequently used to determine the CCC values. The author of paper determined the CCC values of surfactants tested using the linear regression method, and subsequently analyzed the results obtained. To define initially the capability of a tested surfactant to flotation, efficiency of multiparticle model flotation and flotation of individual size fractions were determined in a flotation apparatus for laboratory use in the presence of only foam making agent at CCC value pre-determined during the tests previously carried out. The tests of coalescence phenomenon were carried out with use of selected flotation frothing agents.

Jan 1, 2014

By Chao Li, Kym Runge, Saeed Farrokhpay, Dee Bradshaw

"It is well known that froth phase has a controlling influence on the overall flotation performance. Froth recovery is expected to be affected by froth transportation behaviour. The latter has been represented by different volume balance and energy balance based models. In froth transportation models based on a volume balance, the froth phase is divided into different zones. The total froth residence time is the sum of the time spent in each zone and is dependent on the cell geometry, bubble burst rate and the position that bubbles enter into the froth phase. Alternatively, the froth transportation model based on an energy balance treats the froth phase as a block. In this model, froth height above the concentrate lip is being predicted by minimising the energy of the system which is dependent on the froth geometry, froth stability and froth rheology. The froth height can be used to calculate the froth efficiency of delivering mineralized bubbles from the pulp to the concentrate lip. This paper critically analyses the current models of froth transportation in flotation, and examines the limitations of these models. The importance of froth rheology, which is the resistance to the froth flow towards the discharging weir, in determining froth transportation, is also highlighted."

Jan 1, 2014

By P. Yu, L. X. Liu

"A mechanistic model for tumbling mills was developed during PhD studies at the JKMRC. The concept was presented at the IMPC of 2014, followed by progress in a sub-process of the model presented at the IMPC of 2016. Additionally, a number of papers on the sub-models and breakage function have been published. This paper provides a consolidated summary of the outcomes and status of the model. The overall model structure is presented along with the sub-models such as appearance functions, breakage rate functions, energy distribution, transport, and multicomponent grinding interaction model. The strengths and capabilities of the model structure as achieved to date are presented. The gaps in knowledge and sub-processes are also presented, providing a guide to the areas that require further research and development. It is proposed that the approach developed in this paper provides a structure that can accommodate future mechanistic, dynamic mill model evolution that can incorporate the latest research outcomes and provide a platform to test modelling improvements. This can assist in substantially improving the predictive mill modelling capability in the industry, the importance of which should not be underestimated for this workhorse of mineral processing.INTRODUCTION Grinding is almost the largest single energy-consuming process in the mineral processing industry, which was estimated that about 0.2% of global, and 1.3% of Australia’s electricity is consumed by comminution of gold and copper ores (Ballantyne and Powell, 2014). Although there are some relatively new grinding mills, such as stirred mills, vibrating mills and HPGRs (High-Pressure Grinding Rolls), etc., tumbling mills, as the traditional grinding facilities, are more extensively used."

Jan 1, 2018

By Michael Young, Rodrigo Araya, Cesar Gomez, Le Huynh, James Finch

The performance of flotation machines is closely related to the gas dispersion characteristics achieved in their operation. It is known that frothers play a major role not only in the size distribution of the generated bubble population, but also on the swarm rise velocity. A program to characterize frother effects on two-phase hydrodynamics (air-water) of Jameson cells was undertaken in a fully instrumented pilot unit installed at McGill University. Gas dispersion parameters (gas velocity, gas holdup and bubble size distributions) of the gas-liquid mixture generated in the downcomer were measured after its release into the separation chamber. The unit was operated under different conditions including gas velocity, liquid jet velocity and frother concentration (MIBC). Bubble size distributions for concentrations above the critical coalescence concentration (CCC) were narrow with Sauter mean diameter between 0.6 to 1.3 mm. The bubble size followed the general trend as encountered in other types of flotation machines, i.e., rapid decrease with increasing frother concentrations below 10 µL/L, and almost no change in size for concentrations above 30 µL/L. The CCC was between 10 and 15 µL/L. Strong correlations were found between bubble size and air-to-liquid flow ratio (APR), and between the bulk density of the aerated solution in the separation tank and the vacuum inside the downcomer. An equation is proposed to estimate gas holdup in the separation tank from vacuum and bulk density readings.

Jan 1, 2014

By Francesco Vegliò, Luigi Toro, Francesca Pagnanelli

Different hydro-metallurgical processes have been studied in order to recover manganese(II) from oxidised minerals such as MnO2. An innovative acid reductive leaching carried out in mild experimental conditions was here considered using agricultural and industrial wastes containing carbohydrates as reducing agents. The leaching process was investigated from a kinetic and equilibrium point of view in order to optimise the experimental conditions, develop a kinetic model and obtain high purity concentrated MnSO4 solutions. To study the leaching kinetic a Shrinking Core Model (SCM) was modified by a variable activation energy term to take into account both the variable reagent concentrations and the reaction rate change with respect to the manganese dioxide conversion. The influence of the particle size in the kinetic model is discussed with the support of SEM analysis before and after leaching. Manganese extraction yield was maximised (170-190 g/L) using appropriate carbohydrates and acid excesses, while impurity concentrations were kept low (Fe<10ppm and Ca<300ppm) pre-treating the ores, using hydro-alcoholic leaching medium and different combinations of down-stream processes. In particular precipitation, solvent extraction and adsorption onto active carbon were tested to eliminate Fe and Ca from MnSO4 solutions.

Jan 1, 2003

By Samira Rashidi, Robert Dunne, Raj Rajamani

"High Pressure Grinding Rolls (HPGR) are known for more energy saving, better particle weakening and higher tonnage compared to semi-autogenous/ball mill circuits. These machines are scaled up based on the roll diameter. Therefore, studies in lab-scale HPGRs are always preferred because less weight of feed sample is required to run batch tests. The 50-hp laboratory-scale machine at the University of Utah is the first-ever lab-scale machine designed by late Prof Schönert, the inventor of HPGR. This 200x100 mm [DxL] HPGR is equipped with a robust data acquisition system, which enables running steady state tests using a feed sample as small as 6 kg to 12 kg. This paper presents the HPGR grinding study on two different categories of ores, one hard ore and the other soft. A test series has been performed to evaluate the energy consumption and the relation with size reduction on open multi-pass and locked-cycle HPGR circuits. The effect of feed size distribution and operating conditions has been investigated. Locked-cycle tests result in higher reduction ratio and higher rate of fines (100 mesh, 48 mesh and 28 mesh) generation when compared to open multi-pass tests. Data shows that scale-up trend lines stay consistent even when the circuit configuration and operating conditions are varied."

Jan 1, 2014

By Shankar Narasimhan

Data Reconciliation is a technique that has been developed more than fifty years ago and has been refined over the past few decades and is now being deployed routinely in several process industries such as refining, petrochemicals, mineral processing and hydro-metallurgy. The main aim of data reconciliation is to obtain accurate estimates of all process variables and parameters using plant measurements that are invariably corrupted by random errors and may also contain biases and other gross errors. These reconciled estimates are subsequently used either in offline exercises such as simulation, retrofitting, and de-bottlenecking or in online optimization and control applications. In this presentation, a brief history of data reconciliation, its assumptions, and uses are reviewed. New developments in this area are highlighted such as robust data reconciliation to handle model uncertainties, and reconciliation using identified models derived purely from operating data. Several case studies in the mineral processing industries reported in the literature are also highlighted to indicate the potential benefits of this technique.

Sep 1, 2012

By A. R. Videla, P. R. Brito-Parada, I. Mackay

"Linking results from laboratory scale experiments to industrial flotation behaviour is challenging. Typically, such experiments involve batch tests in which the system does not operate at steady-state, making it difficult to infer the effects that operating conditions have on flotation performance. In order to overcome this limitation a 4-litre recirculating tank was previously developed at Imperial College London. This tank is capable of reaching, and operating at, steady-state by recycling overflowing concentrate back into the feed. As well as instruments to control operating conditions, it is fitted with a system of sensors that allow the surface of the froth to be dynamically monitored. From this information, it is possible to measure the air recovery—a proxy for froth stability. Thus, this bench-scale tank can be used to understand the effect of differing operating conditions on flotation performance at steady state. However, so far, this cell has only been used to investigate idealised systems with only one or two species.Reprocessing of tailings dams is not only environmentally desirable but also increasingly economically feasible due to the declining head grades of primary deposits. There is also the added benefit of no further milling being required prior to flotation. However, the effects of fine and ultrafine particles on froth stability are not yet fully understood. In this work, the bench-scale continuous tank has been used for the first time to determine the flotation response of a complex feed, consisting of samples from a copper tailings dam, to changes in operating conditions. It was shown that the froth stability in the system is comparable to that of previous work and industrial tests, with a peak in air recovery being found at a superficial gas velocity of 1.13 cm/s. There is scope to optimize the froth stability of tailings flotation for enhanced metallurgical performance."

Jan 1, 2018

By Alan Bye, Nick Beaton, Wayne Stange, Emmy Manlapig, James Groutsch

"For over 50 years the JKMRC has applied a unique education-based research strategy to develop modelling and simulation tools that are "fit for purpose" for industrial application. These tools are widely used for the analysis and optimisation of mining and mineral processing systems. As the mining industry context evolves to deal with declining head grades, falling commodity prices, professional skill shortages, and an exponential growth in the volume of data, a re-evaluation of what "fit for purpose" means in the future is required. It is contended that future industry needs for simulation-based decision support requires tools and systems that incorporate: • Multi-component modelling, • Models of complex sensor behaviour, • Integration across the physical value-chain e.g. resource to market, • Financial analysis including risk and uncertainty assessment, • Enabling multi-criterion decision making (technical, financial and environmental), • Ore-body knowledge at multiple scales, • The use of fundamental numerical models to provide insight for the development of practically useful models, and • Decision support delivered through remote operations centres. JKMRC and CRC-ORE are collaborating to develop a next-generation simulation-based decision support system that assists in addressing these challenges. The Integrated Extraction Simulator (IES) combines many years of pragmatic and successful work integrated with modern technology trends - in particular the rapid move to cloud-based computing."

Jan 1, 2014

By C. Eswaraiah

A drop weight tester (DWT) is designed and fabricated for the purpose of analyzing single particle impact breakage characteristics of different ores. The test results were evaluated through the breakage distributions of different size fractions at five different input energy levels. The breakage parameter, t10 (Narayanan, 1986) is used to represent the degree of size reduction, which is assumed to be representative of the breakage product size distribution obtained from single particle breakage. Representation of breakage distribution with this method consist one more breakage dependent parameter (a). To evaluate the impact breakage characteristic of different ores, a mathematical relationship was established to analyze the dependency of breakage parameters, t10 and a on specific input energy (EIS) or fracture energy (Tavares et al, 1998, Fandrich, 1998), specimen size (?) and Bond?s operating work index (WI0). Mathematical results developed gives adequate model for all ball mill modeling. It can be concluded that the drop weight test method is a useful and practical way of evaluating the impact strengths of various ores and results of which can be used in modeling mills.

Sep 1, 2012

By Kirsten Corin, Margreth Tadie, Michael Nicol, Cyril O’Connor, Jenny Wiese

"The recovery of base metal sulfides has been emphasized in the study of the flotation of platinum group minerals (PGMs). However, in some ore bodies such as the Platreef of the Bushveld Complex, sulfide minerals do not make up the majority of the more significant PGM mineralization. It is therefore important to develop a fundamental understanding of how non-sulfide and sulfide PGMs respond respectively to the current reagents being used in the flotation of these minerals. The interactions between platinum and palladium sulfides and tellurides with the thiol collectors, ethyl xanthate and diethyl dithiophosphate, have been studied using the electrochemical techniques of mineral rest potential and voltammetry measurements. The differences in the interaction between the two types of minerals with these reagents were of key interest. Rest potential measurements indicated stronger electrochemical interactions of the minerals with ethyl xanthate than diethyl dithiophosphate. These measurements also indicated that there are different forms of oxidized collector present on the mineral surfaces after interaction with ethyl xanthates. Voltammetry measurements support these observations and suggest differences of adsorption on Pt and Pd sulfides than on Pt and Pd tellurides. Voltammetry also showed that there was minimal electrochemical interaction between the Pt and Pd minerals and diethyl dithiophosphate."

Jan 1, 2014

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

By Rodrigo Carvalho, Marcelo Tavares

"The accurate prediction of the outcomes in SAG grinding has been a great challenge to researchers and process operators. The use of Discrete Element Methods for modelling charge motion in tumbling mills allowed extraction of not only visual data, such as the position of the charge, but also the energy transfer in terms of the energy dissipated in the collisions. One of the advantages of DEM technique is its scalability to different mill sizes and some operational conditions such as speed and mill filling. Recently the Authors proposed the coupling between DEM outputs and modified Population Balance Models that with proper dedicated ore characterization tests allows the prediction of the outcomes of an energy transfer example. This modelling framework is called the Multi-scale Mechanistic Modelling. This paper shows the application of the Multi-scale Mechanistic Model to a set of pilot plant tests ran at different operational conditions. The prediction of the grinding outcomes and even the stabilization of the charge was achieved by the use of the Multi-scale Mechanistic Modelling framework. This approach has a great potential to become a powerful tool to design and optimize AG and SAG mills."

Jan 1, 2014

By Robert Peretti, Antonio Zucca, Andrea Manni, Silvia Serranti, Paolo Massacci, Marcello Ghiani

"This research study aims to identify the best way of processing colemanite concentrates, which is sometimes characterised by a high content of As sulphides, finely disseminated in the matrix of the boron mineral.Scanning Electron Microscopy analyses have been carried out to identify all the mineralogical constituents and to evaluate their distribution, frequency and mutual relationships.Colemanite concentrates are mainly used by chemical plants producing boric acid and gypsum. A high As content is unacceptable for further use of tailings coming from the leaching process and is very dangerous if the tailings are released into the environment.Grinding tests were conducted under dry and wet conditions in order to assess how As sulphide minerals are released from colemanite minerals.Some of the difficulties involved in obtaining selective liberation of As sulphides from colemanite minerals were overcome, although a certain decrease in As grade in the higher size classes was registered.All comminution products have been analysed by SEM, monitoring the degree of liberation of the mineralogical constituents.Flotation tests have been carried out, following different comminution flow sheets and adopting different modulating conditions in order to evaluate the influence of the following parameters:- collector concentration versus contaminants concentration (As minerals) in the froths;- kind and quantity of the frothing agent to maximise mineral recovery in the froths;- kind and quantity of the modulating agents to maximise the selectivity of the process;- conditioning modalities in order to obtain a suitable dispersion of fines and avoid segregation of Asbearing particles together with colemanite ones;- solid/liquid ratio of the pulp to minimise the water flow into the plant.Information obtained from SEM analysis has been adopted to evaluate segregation phenomena in pulp and in flotation products."

Jan 1, 2003

By Erico Tabosa, Kym Runge, Peter Holtham

"The distribution of turbulence in an industrial scale mechanical flotation cell and its effect on flotation kinetics was studied using a fully instrumented 3m³ cell. The cell was operated at a copper concentrator in Australia with different combinations of airflow rates, impeller speeds and sizes and cell aspect ratio providing a wide range of hydrodynamic conditions. An analysis of the overall flotation cell performance showed that the overall copper recoveries were very similar for the conditions tested. However, by decoupling pulp effects from froth effects it was possible to evaluate whether the changes made affected the pulp zone and/or froth zone responses. The analysis showed that the overall recovery has the potential to be up to 7% higher, if not limited by froth recovery. The cell hydrodynamics were characterised by the overall energy dissipation rate and, to understand the distribution of turbulence within the cell, localised measurements were made using a piezoelectric vibration sensor. By comparing the metallurgical performance of the cell with the different hydrodynamic conditions it was found that the collection zone flotation rate constant was not directly related to overall energy dissipation, as is commonly observed at laboratory scale, where energy is usually only changed by varying impeller speed. In this study, energy dissipation rate was varied using not only impeller speed, but also the size of the impeller and the cell aspect ratio. In these circumstances, the best correlation was found to be between the collection zone flotation rate constant and the distribution of turbulence, i.e. the size of the turbulent zone."

Jan 1, 2014