Search Documents

By Peter N. Holtham, Geoffrey J. Lyman, Matthew S. Brennan, Vijay Subramanian, Rui X. Rong, Tim Napier-Munn

Studies of the dense medium and classifying cyclone have been conducted using the Fluent computational fluid dynamics code. Experimental measurements of the segregation of medium in a DM cyclone by gamma-ray tomography and the performance of a novel classifying cyclone design have been used to test and interpret the 3D CFD predictions. The predictions of the CFD models were generally plausible, though more work is required to incorporate turbulent mixing of the slurry in the cyclone. The benefits of the novel cyclone design were explained by the model in terms of conditions favouring reduction of short-circuiting at the vortex finder, and the removal of the air core from the wall at the apex. Work is continuing.

Jan 1, 2003

By Güven Önal, Zeki Dogan, Neset Acarkan

In this study, it was aimed to produce glass sand with less than 0.15% Fe2O3 for window glass manufacture and glass ware grade with less than 0.05% Fe2O3 from Yalõköy sand dunes having 0.48% Fe2O3. In order to reduce the iron content of the Yalõköy sand dune, magnetic separation, flotation, scrubbing, H2SO4 leaching and drying were used. The optimum iron removal was achieved by H2SO4 baking at 150oC and 8 kg/t acid consumption. In this way, sand concentrate was obtained with 0.04% Fe2O3 and sand recovery 95%.

Jan 1, 2003

By L. X. Zhang

It is very important to improve the A/S of the bauxite to keep sustaining development of Chinese alumina industry. This article summarized the issues of several desilication technologies, which are named as normal flotation, reverse flotation and bi-flotation. Normal flotation has been got industrialization because of technologic advantages. The campaign experiment of reverse flotation which has principle advantage was carried out in 2004, and bi-flotation which is combined the advantages of both the normal flotation & reverse flotation has the most attractive future. The development of bauxite flotation will be focused on multi-utilization of bauxite by mineral process technology and corresponding equipment which satisfied with the flotation technology and reagent.

Jan 1, 2014

By Jussi Liipo, Pertti Lamberg

Image analysis of process samples gives precise data on particles and their composition. The results presented as liberation tables and maximum grade-recovery curves are informative and good for comparing different feed materials, but basing on them it is impossible to predict the results of the mineral separation process (flotation, magnetic separation, gravity separation) in details. The exact and extensive particle data has remained large unused, because there has not been a model nor a modeling program to include particle information in mineral separation processes. This study introduces a method, which incorporates the particle data in the kinetic model of flotation. Similar approach can be used in modeling other separation processes, too.

Jan 1, 2003

By Christopher Greet, Alfred Van den Bosch

"A range of media alloys with different chromium content has been developed for different milling applications. When selecting the grinding media alloy for a specific application, consideration must be given to the mineralogy of the ore as well as the water chemistry and mill operating conditions. The paper describes the influence of the microstructure of the alloy on the wear resistance and the influence of the chrome content on the pulp chemistry and downstream processing. It also provides information about the in-house tools which Magotteaux have developed to predict the media wear rate with different chromium content alloys and predict their flotation efficiencies that can be expected when the right alloy is selected for a particular application."

Jan 1, 2014

By Lev O. Filippov, Inna V. Filippova

"Flotation of metal-organic precipitate in column is sensible to the hydrodynamic conditions determined by the diameter and dispersion of rising bubbles. Thus, the dissipation energy in the system resulting of ratio between superficial fluid velocities (gas and feed) in the column is main parameter for precipitate flotation in pilot or industrial scales. The knowledge of aggregate structure can lead to better understanding and control of many solid-liquid separation processes (sedimentation, filtration, thickening and flotation). Many authors highlight the importance of structural properties (i.e. floc density) for different industrial applications.This approach was applied to achieve the flotation of molybdenum precipitate obtained with organic collector a- benzoin oxime for molar ratio 1:2. Flotation tests were carried out at lab and pilot scale with rigorous control of the size distribution of bubbles and the values of dissipation energy in the systems.Structural evolution of Mo precipitate and the destruction phenomenon were described by using of the various fractal dimensions, which were deduced from experimental data of laser sizer devices and optic microscopy pictures of precipitate. The boundary fractal dimensions Df of large aggregates allows the detection of a structural scaled organisation of the precipitate flocs (elementary germs-microflocs-flocs). The fractal dimension of the particle mass distribution Dm=1.94+/-0.15 suggests, for the floc structure, a cluster-cluster aggregation process either with structuration or breakage phases.The destruction process characterised by a fractal dimension Dg controls the final size distribution. The floc destruction occurs for the lowest values of the dissipation energy (0.02 – 0.05 W/kg). The relation between Dg and the system dissipation energy Et (for various energy E dissipated in the agitation cell) involves the different breakage mechanisms of aggregates in function of the aggregation bonds (Fig. 1).Besides, the destruction of floc produces very fine particles (elementary germs or/and microfloc according to E), and thus induces the difficulties of bubble-particle attachment during the precipitate flotation in column. Yet, very small bubbles, which increase the collision probability, would make the flotation feed rate decrease. The optima ratio between superficial gas and feed velocities are determined for bubble size and flotation results control at pilot scale. Thus, it is possible to providing the catastrophically changes of precipitate column flotation by using the following established experimental relation: the mean diameter d of the flocs under the action of dissipated energy Et follows a power like law (d~Et -0.5). This allows estimating the size distribution of flocs for the E values which can appear in column."

Jan 1, 2003

By Peter Ay, Elena Geldt, Satyanarayana Narra, Claudia Glaser

"As part of an extensive research and development project, the technical and technological implementation of the patents DE 10 358 245 B4 and WO/2006/056172 was developed with the goal to present a system for the innovative roll- agglomeration method exhibiting a) a modular setup, b) high productivity and c) suitability for the use of fibres of different applications and product qualities (high technical availability). The processability of the rarely adapted raw fibre materials into high-value intermediates brings forward their application in various processing industries dealing with the production of fibre reinforced plastic composites, heat and sound insulation materials as well as with the gas and liquid cleaning systems. The innovative fibre agglomeration presented in this work is performed by a combination of rolling and pressing the fibres inside a geometrically pre-defined roller table, in which the fibres are forced to bound through entanglement / interlocking and develop into agglomerate. The agglomeration process is carried out in a circular roller path with a rectangular cross-section. For improving the adhesive properties, and thus the quality of the agglomerates, the boundary surfaces can be covered / patterned / lined / coated with various materials. The resulting roll-agglomerates are characterized by different quality parameters such as mass, shape, porosity and stability. The properties of the agglomerates can be modified based on the end-product requirements by slightly modifying the process parameters such as filling ratio, residence time, roll-geometry, procurement of the boundary surfaces, etc. as well as the input material properties such as fibre type, fibre fineness, fibre length, etc. This work carries extensive work on the roll-agglomeration system to assess the effectiveness and the variability based on the product specifications. Agglomerates were produced with wide variations of system parameters with varying roller table liners and based on evaluation of the produced agglomerates with respect to their mass, the geometrical dimensions, their porosity and stability."

Jan 1, 2014

By E. D. Gálvez, L. A. Cisternas, R. Acosta-Flores, F. A. Lucay

The optimal design of flotation circuits is a complex task as the mathematical model is nonconvex and includes binary variables, making it difficult to obtain the global optimum solution. Furthermore, there are epistemic uncertainties in the evaluation of the performance of flotation and grinding units. This epistemic uncertainty is generated because semi-empirical models are used with parameters that are a function of the circuit design. Therefore, it is not possible to adjust these models based on experimental data as the experimental conditions are not known at the design stage. This paper analyzed the effect of the uncertainties in flotation units, grinding units, and both units simultaneously. As a result, this paper demonstrated that, for a specific mineral, few structures exist that are optimal for a wide range of values in flotation stage recoveries and liberation fractions in the grinding units. A remarkable aspect of this work is that it introduces grinding, one of the operations that most influences energy consumption in mining and has not been studied in depth in the design of flotation circuits. To demonstrate that there are few optimal structures for a specific mineral when there is high uncertainty in the process units, an a priori approach or first principles was used. Two methods of proof were used: proof by construction and proof by exhaustion. Proof by construction is the construction of a concrete example with a property to demonstrate that something having that property exists. In this work, two examples were considered in the proof by construction. In proof by exhaustion, the conclusion was established by dividing it into a finite number of cases and proving each one separately. Here, 30,000 cases were included for each example. This knowledge can be useful in circuit design including grinding, e.g., after selecting the circuit (or the most promising circuits), and the equipment design parameters and operating conditions can be defined based on simulation and laboratory tests.

Jan 1, 2018

By Barbara Peszko, Tadeusz Tumidajski

"Various methods are used to identify the grain composition of fine-grained materials because so far the best notion of the grain size has not found and there is no single method to define grain size. In practice the most important are the following grain sizes: sieve diameter dt, effective diameter de, substitute diameter d2 and statistical diameter dst. The mentioned sizes are connected with the methods of analyses. Methods are correct but measure of different diameters of grain and their results are different.The most important methods of granulometric composition analyses are the following: the sieve analysis, the sedimentation method and the conductometric (Coulter-Counter) methods. In each of these methods various physical features connected with the grain size are measured.Sizes of grains defined by different formulas or measurement methods are associated mutually by dependencies, which are results of differences in their determination or calculation ways.In order to determine principles of results comparison the classification analysis results on linearly or non-linearly connected was made. Analyses which results demonstrate the same type of grain size distribution are treated as linearly connected. It is based on the general formulae which determines the probability density function of random variable which is the function of another random variable having the known density function.In order to show the difference in the grain composition of the same material qualified by different methods, three methods of investigations were distinguished. They are: the sieve analysis performed on precise sieves, the sedimentation analysis performed by means of a sedimentation balance “Sartorius” and the microscopic analysis performed after preparing a microscopic slide where the grain size was measured by the Martin’s method.The materials selected for the analysis were quartzite and glass, which are the material of high purity.The obtained results of mentioned above three types of analyses were approximated by distribution function of four distributions: Weibull, log-normal, cut Weibull and Gaudin-Schuhman."

Jan 1, 2003

By Wieslaw Wszolek, Kazimierz Trybalski, Ryszard Tadeusiewicz

"The present work is an attempt of neural networks application to the description of complex technological processes of useful minerals enrichment. Neural networks are among the most convenient tools that have been added recently to the professional toolbox of modelling specialists, and they are particularly useful in realisation of multiple practical tasks, characterised by great complexity and relatively low amount of available a priori knowledge about the nature of the modelled processes. They are successfully applied in a wide variety of problems related to data processing and analysis, their prediction, classification or process control. The object of the study was the technological site for copper ore enrichment by floatation, being a part of technological line in the Ore Enrichment Plant ZG „Polkowice”, KGHM „Polska Miedz” S.A.The main floatation is described by parameters measured by the supervision and control system, which can be grouped with respect to the functions performed for the object:a) input function: aFG - copper contents in the feed material for the main floatation [% Cu]b) control: LFG I - level in the first section of the floatation unit [%], LFG II - level in the second section of the flotation unit [%], LFG III - level in the third section of the floatation unit [%], LFG IV - level in the fourth section of the floatation machine [%]. These quantities compose the five elements long input vector.c) output function: ? - copper contents in the waste from main flotation (final waste) [% Cu]The measurement data, containing the above-mentioned parameters, have been collected during the normal operation of the plant. The data contain momentary values of the parameters, registered every 10 minutes, covering the period of 196 hours of the plant’s operation (1176 measurements for each quantity).For construction of the model used in the study classical multi-layer perceptron-type neural network has been used, with one hidden neuron layer, trained by the error back propagation method."

Jan 1, 2003

By Andrei Kolesnikov

"The results of CFD modelling of complex heat transfer and fluid flow in DC electric arc furnace are presented. The model includes energy and momentum conservation equations for gas phase (freeboard), liquid phase (molten bath) and solid phase (walls and graphite electrode). Simulation revealed the temperature distribution in all elements of the furnace. The influence of the arc and melt resistance on the temperature distribution is discussed.Computational Fluid Dynamics (CFD) becomes a powerful tool for better understanding of the heta and mass transfer processes in industrial equipment. In current work commercial CFD software STAR-CD was used to analyse complex heat transfer in DC EAF, including forced and natural convection, conduction and radiation. The finite volume solver using body-fitted grid was used. The grids are non-staggered hexahedral grids and all variables are evaluated at cell centers. Steady-state calculations were performed with PISO algorithm. For the convective terms in equations the hybrid differencing scheme was used. An upwind differencing scheme was used for other terms. Conjugative heat transfer was taken into account in heat transfer calculations on gas-liquid and gas-solid interfaces. Radiation from the arc and inner surfaces of the furnace was calculated by a “discrete beam” approach. Source term in energy conservation equation was used to represent heat release in the zone occupied by the DC electric arc. The volume of this zone was calculated from semi-empirical equations. The Joule heat and heat due to chemical reaction in the melt were calculated through the source term expression in the energy balance equations for the slag and metal in the molten bath. The gas acceleration in the electric arc zone and stirring of the slag due to gas bubbles motion were accounted for by momentum source terms in the momentum conservation equations.All simulations were carried out using axisymmetric two-dimensional grid. The grid resolution was 20 cm per cell in most part of the furnace, except the arc region, where resolution was 5 mm per cell."

Jan 1, 2003

By Grigori I. Voiloshnikov

"The long-term studies on development of the technology for precious metals recovery using activated carbons executed by Irkutsk Research Institute of Precious and Rare Metals and Diamonds (IRGIREDMET) were finished by successful introduction of ""Carbon-In-Pulp"" process in Russia (1984).The review of some studies, accompanied by development, testing, introduction and subsequent improvement of considered technology, carried out by IRGIREDMET is submitted in this paper.The basic regularities of gold and silver adsorption onto activated carbons from cyanide, thiourea, chloride and other solutions are investigated using various phisical-and-chemical methods ( X-ray photoelectron spectroscopy, Voltammetry, X-ray microanalysis etc.). Using voltammetry it was found that gold cyanide is adsorbed onto activated carbon without chemical changes and partial reduction of gold is possible only during drying of carbon samples before analysis. Metallic gold and silver may be dissolved in boiling cyanide solution at argon atmosphere when they are in contact with activated carbon surface. Gold adsorption mechanism for activated carbons has great differences from that for the ion exchange resins. So the most probable mechanism for gold adsorption from cyanide solution is adsorption of neutral molecules Me n+[Au(CN)2-]n.X-ray quantitative microanalysis technique was developed to investigate gold adsorption kinetics onto spherical activated carbon. Gold profiles in granules of activated carbon for various adsorption conditions were tested (carbon particle size, mixing efficiency, adsorption time etc.). Mathematical methods were used to determine gold diffusion coefficient within activated carbon granule.The properties of various types domestic and foreign activated carbons were tested and the demands to activated carbons suitable for gold extraction were determined. As a result foreign and Russian activated carbons AG-90, AG-95 were recommended for use at Russian gold plants. Domestic adsorbents are cheaper in comparison with foreign activated carbons. The methods to restore carbon adsorption activity including acid treatment and thermal reactivation were tested and implemented."

Jan 1, 2003

By Olgac Kangal, Ali Guney, Gulay Bulut, Guven Onal, Ayhan A. Sirkeci

"This paper presents the results of pilot scale column flotation applications in the two different copper concentrators in Turkey. The Kure Copper Plant is located in the Black Sea Region of Turkey. There are four types of copper ore formations in the Küre region. The so called disseminated ore contains 60-70% siliceous minerals of which Cu content is around 0.5%. The remaining three ore types in the region are generally named “massive ores”. The massive ores may assay up to 15% Cu and the major gangue mineral is pyrite associated with relatively little amount of silica. The main copper mineral in all formations is chalcopyrite that is generally finely disseminated in pyrite therefore satisfactory liberation requires ultrafine grinding especially in the case of the massive ores.A blend of these ores, containing 1.7% Cu on average, is fed to the Küre Copper Concentrator whose capacity is 1.2 million tons of ore annually. Unfortunately, today the plant operates below its designed capacity treating around 550,000 tons of ore per year. The design of Küre Copper Plant was undertaken by the Outokumpu Company of Finland. The operation of the plant is almost automatically controlled and the feed, the concentrate and the tailings are assayed for copper by on-line XRF analysers. The particle size analysers are also connected to the grinding circuit. The feed material is ground almost below 50 µm and the flotation is conducted around pH 11.2 using Hostaflot X-231 (a trade mark of Hoechst) as the collector. Apart from the rougher cells, there are four cleaning stages after which the final concentrate is produced. A very fine grinding requirement for an adequate level of liberation makes successful flotation difficult and consequently the final concentrate can hardly assay above 16% Cu leading to copper losses around 20% in the tailings that contains 0.3% copper."

Jan 1, 2003

By E. V. Manlapig

"Some researchers postulate that the floatability of a particle in the streams of an industrial flotation circuit is a function of only its size and liberation properties. They do not discount the importance of reagent or surface contamination on particle floatability but assume that the reagent coverage or surface contamination of the particle surfaces is essentially uniform for particles of similar size and mineralogical properties when subjected to a particular chemical regime.In this paper, samples collected from laboratory batch flotation tests performed on different streams of Cominco Teck’s Red Dog lead cleaning circuit have been characterised with respect to their size, mineral association and mineral surface exposure. An analysis of this data set indicates that the flotation rate of a particle in an industrial flotation environment is a strong function of its size and also of the mineral composition of its surface (Figure 1). It is a function not only of the proportion of the most floatable mineral on its surface but of the floatability of the other minerals on its surface as well.Nodal analysis is a technique developed by the authors that can be used to determine if the ore floatability of a particle class is changing within an industrial flotation circuit. This technique is used within the paper to show that the floatability of all the different types of particle classes studied within the Red Dog circuit did not change significantly during the lead cleaning stage. However, it was also concluded that particles with the same size and mineral composition did not all float with the same rate in the individual streams of the flotation circuit, but with a distribution of rates. It is suspected that the particles within a size mineralogical class in this study have varying degrees of surface oxidation and that this results in varying particle flotation rates. Thus basing particle floatability on size and mineralogical effects, alone, in a flotation model is not always valid."

Jan 1, 2003

By Milind Chikhale, Yogesh Reja, Keshav Tyagi

"Using computer models for process design is a well-established practice in the mineral industry. In recent years, there has been emphasis on modeling of processing plant as a whole. This paper covers development of a model for iron ore processing plant. Basic purpose of this model is to use simulation for efficient plant design and help in minimizing risk at the very initial stage of project. Feed to the plant is Hematite iron ore from a deposit located in Eastern India. Major processes in plant include multi-stage crushing, screening, grinding, classification, gravity and magnetic separation. Models for different process equipment were developed using METSIM and Bruno simulation software. Customized black box models and known mathematical models were used for predicting behaviour of processes. Data needed was identified based on the model type. The required data was arranged as per project requirement and processed for development and validation of model. Equipment models were developed and validated on the basis of process design basis, theoretical and empirical data from technical references, and bench scale and pilot scale testwork results. Validated process circuits were assembled to simulate and represent the entire plant process. The result shows expected variation in plant output as a result of variation in the feed. Study concludes by discussing potential applications of the model in plant design."

Jan 1, 2014

By Cristian Melipichun, Leopoldo Gutierrez

The effect of the presence of fine particles of kaolinite and montmorillonite on the induction time for attachment of bubbles to particles of a copper concentrate rich in chalcopyrite was evaluated in this study. The results obtained from induction time measurements were correlated with data of concentrate recovery achieved from micro-flotation experiments. The effect of the presence of fine clays on the stability of air-water foams was also studied through measurements of the foam volume generated in a glass column under different conditions of air flowrate and MIBC concentrations. The experimental results showed that induction time values increased under the presence of fine clays. These data agreed with results obtained from micro-flotation experiments that showed that copper concentrate recovery decreased when clay particles were added. Although a complete explanation for these results is not available, the fact that experimental observations indicated that fine clay particles tend to adsorb on gas bubbles, suggests that the interactions between fine phyllosilicate particles and air bubbles have to be taken into account when studying the effect of clays on froth flotation. Foamability tests showed that the presence of fine clay particles increased the foam volumes only at high additions of clays and low dosages of MIBC. The results obtained at higher dosages of MIBC showed that the foam volumes decreased with the addition of clays under all the experimental conditions tested.

Jan 1, 2014

By Peter Ay, Matthias Schmidt, Christian Stollberg

"The aim of the study is the development of a seamless reprocessing procedure for plastics from automotive parts which are to be recycled (figure 1). Before automatic sorting can be performed, the materials first have to be identified by a near and middle infrared (NIR/MIR) system.Materials reprocessing basically includes research on size reduction, decoating, and roughening of the materials. Variants of the necessary technologies such as different shredders and mills and physico-chemical or mechanical decoating and roughening procedures together with the fundamental approaches underlying these methods are investigated and benchmarked. This allows the theoretical determination of the impact of each method not only on the overall process but also on the result achieved within the identification and sorting procedures."

Jan 1, 2003

By V. Panayotov

Flotation practice usually uses the Sheridan-Grizvold selection method, where zinc minerals are activated by addition of copper sulphate and depressed by addition of sodium cyanide, ferrous sulphate, zinc sulphate or various mixtures of these reagents. The method is suitable for the flotation of sulphide ores, which, do not bear precious metals, such as silver and gold. When classical technology with cyanide is applied to depress zinc minerals, the precious metals, which are contained mainly into lead minerals, are dissolved into the solution, which represents a technological loss. In order to avoid this process, we have made use from the knowledge on the semiconducting properties of lead and zinc sulphide minerals. We have changed the surface charge of minerals by means of an electrochemical impact realized by direct current application and specially designed set of insoluble or semi-soluble electrodes. Thus, flotation separation of lead and zinc sulphide minerals was achieved, at increased recovery of lead, zinc and silver. The technology proposed was implemented in a Bulgarian flotation factory; it worked over 5 years providing significant, positive economic, technological and environmental impact. Silver recovery was increased by 0.8-1.0 %, and lead and zinc recovery was increased by 0.8-1.0 %. The technology is easily controllable and adaptable to raw materials and parameters of the technological process. It can be paid back for 6-8 months after the implementation. Keywords: lead-zinc flotation separation, electrochemical impact

Sep 1, 2012

By Takashi Furuyama, Sukeyuki Mori, Tetsuo Ikejiri, Abel Bissombolo, Yoshitsugu Koga, Masamichi Hata

In Japan, the increasing number of automobiles newly put into circulation is followed by a larger amount of end-of-life vehicles from which valuable metals are being recovered in shredder and related plants. Nevertheless, the processing of these scrapped automobiles which has reach a recycling rate of about 75%, generates a huge volume of waste, locally known as .shredder dust., which is commonly dumped on landfills. The saturation and depletion of valuable disposal sites in the near future, as well as the social and environmental problems related to the accumulation of these residues disposals are a major challenge for the automobile shredding industry. To promote an alternative way for the management of these disposals, the authors have initiated and developed a system for their processing, on the basis that more than 90% of these residues are combustible materials that could be used as refuses derived fuel. The developed processing system points out the mechanical/physical separation of shredded automobiles residues into two streams by means of an ascending airflow passing through a vertical zigzag separator and, the compaction of the light-density fraction of materials (that is about 60% of auto shredder residue recovered as float products and mainly composed of cotton, sponges and light plastics) to form a refuse derived fuel (RDF), the density, calorific value and ignition loss of which are 1.19 g/cm3, 4470 kcal/kg and 60.7%, respectively.

Jan 1, 2003

By I. Zimbovski, T. Ivanova, E. Chanturiya

"The mechanism of selective action of new reagent-collector AMD of phenylpyrazolone class with sulphide minerals of copper-zinc ores was theoretically and experimentally investigated. The validation was carried out using spectrophotometry, Fourier-transformed spectroscopy, laser scanning microscopy, analytical scanning electronic microscopy and flotation. The reagent AMD in alkaline medium on chalcopyrite and sphalerite surfaces is chemically sorbed, forming strong hydrophobic complexes with copper and zinc. On the pyrite surface, AMD acidifies ferric iron ions with formation of hydrophilic oxidates. The selectivity of collector AMD with separation of sphalerite and pyrite in alkaline condition (??=10) is confirmed by flotation experiments. With AMD, the obtained recovery of sphalerite was 34.4 % at extraction of pyrite of 3 and 6 %. The low recovery of sphalerite is caused by weak collective ability of the reagent. The use of copper sulfate increases the flotability of sphalerite up to 82.0 % at extraction of pyrite of 5 and 2 %. It has been experimentally determined that during the flotation of sphalerite and pyrite by AMD reagent with addition of butyl potassium xanthate, and in the presence of copper sulfate and regulator of formation of ammonium thiocyanate complexes, the recovery of sphalerite is increased to 91 %, while the recovery of pyrite is decreased to 3.8 %. The chalcopyrite is floated actively (90 % recovery). The results of all theoretical research and experiments were confirmed by flotation tests on ores and products of dressing in industrial conditions."

Jan 1, 2014