Search Documents

By R. C. Nkhoma, K. A. Annan, S. Ngomane

The effects of welding current, electrode force, and welding time in a resistance spot weld were studied to investigate the effectiveness of welded joints between a thin EN10130: DC04 material and a thicker 817M40 part, through analysis of the microstructural and mechanical properties. All welded specimens were subjected to tensile testing at room temperature (25°C) and sub-zero temperature (–46°C) to test the strength of the welded joints. No full button failure was observed at either room temperature or sub-zero temperature after optimization of the weldng parameters. The fusion zone was observed to consist mainly of martensitic phase, due to rapid quenching, while the HAZ was composed of clusters of martensite in a ferrite and bainite matrix. The base 817M40 metal remained fully ferritic after welding. The hardness was found to increase with increasing welding current. An increase in nugget size, indicating good fusion of the weld, was observed with an increase in the welding current.

Oct 1, 2021

By T. Wajima

Recovery of precious metals from hydrochloric solution using sulfur-impregnated carbonaceous adsorbent was examined. Sulfur-impregnated adsorbent with adsorption abilities of precious metals can be prepared from activated carbon using pyrolysis at 600 oC under 10 % H2S gas flow, and the adsorbed Au and Pd were precipitated as metal phase on the surface of the adsorbent. The adsorbents after precious metals adsorption were burned at 700 oC to obtain the residue mainly containing precious metals. The adsorbent was found to be highly efficient and selective for the uptake of precious metals, especially Au, in the presence of excess base metal ions from simulated spent automotive catalyst eluted solution using a column pack test. Keywords: Gold, platinum, palladium, reduction precipitation

Jan 1, 2020

By J. D. T. Steyl, G. B. Grobbelaar

Two examples are presented to show how first-principles – physical-dynamic, here termed behavioural modelling, is used to build digital twins of unit processes or plant sections. A digital twin equips the metallurgist, engineer, or manager with a powerful tool to conduct model-based design, solve process problems and improve operations. It is particularly relevant in the metallurgical industry where plant behaviour is highly integrated, rates respond to particle dynamics in complicated ways and mineral flows across flowsheet boundaries are variable. The two examples used to illustrate the behavioural modelling methodology are 1) a flotation column unit block and 2) a multi-compartment pressure oxidation (POX) autoclave. It simultaneously accounts for the causality of material and heat flow, takes care of mineral particle dynamics, physical and chemical reaction rates, the plant control system behaviour and operating heuristics. The basic modelling building blocks are briefly reviewed to show how the behavioural modelling attributes constrain the mass and heat flow dynamics to physically realistic values. This is the key aspect that allows project value to be created with minimal input data. Behavioural modelling is especially valuable where the process dynamics are poorly understood or controlled, where data-driven models fail to cope with complex and highly integrated plant behaviour, where quality data is lacking or limited funds are available to generate the data required to fill all the knowledge gaps. The paper uses the flotation column example to show how behavioural modelling and simulation can be used to optimise plant performance, while the multi-compartment POX autoclave model shows how it can be used to diagnose process problems. It also provides a solid basis to add a data-driven component – in other words, it makes a strong case for a hybrid first-principle– data-driven approach to solve process problems, especially where microscopic-level understanding is poor. Keywords: Digital twin, behavioural modelling, simulation, flotation column, pressure oxidation

Jan 1, 2020

By G. Akdogan, C. Dorfling, L. Bumhira

Hydrometallurgical recovery of rare earth elements from phosphors obtained from end-of-life fluorescent lamps involve leaching followed by extraction of dissolved metals from solution. The similar properties of rare earth elements make selective recovery of individual elements challenging. The aim of the paper was to investigate the extraction of primarily Ce and Tb from a hydrochloric acid leach solution also containing other REEs (Y and Eu) as well as impurity metals (Ca and Al). Batch solvent extraction tests were performed with di-(2-ethylhexyl) phosphoric acid (DEHPA) at various extractant concentrations, pH conditions and organic-to-aqueous ratios. Extraction isotherms and separation factors for the range of conditions tested in the project are presented in the paper. Separation of individual REEs was not possible at the conditions investigated, but a mixed REE product can be obtained in a single extraction stage operated at a pH of 0.5 using 1 M extractant at an O/A ratio of 1; these conditions resulted in complete Y and Eu, 95 % Ce and 99 % Tb extraction with Ca and Al coextraction below 7 %. A better quality mixed REE product can be produced by operating at lower pH and extractant concentration conditions, but this would require a larger number of extraction stages. Keywords: Cerium, fluorescent lamps, phosphors, rare earth elements, solvent extraction, terbium

Jan 1, 2020

By S. Banisi, A. Ghasemi, A. Fazeli, M. Rezaei, M. Maleki-Moghaddam

The efficiency of high capacity SAG mills which consume large amounts of energy depends on the grinding and slurry discharge systems. A 9.7 m × 4.8 m SAG mill with radial pulp lifters is used at the Sarcheshmeh copper complex grinding circuit. The process audit indicated that the pulp discharge system was not efficient, judging by the amount of slurry flow back and carry over. These events limited the discharge rate and increased the wear rate of the lifters and grates. With the objective of finding a new design for pulp lifters, a 3D DEM software (KMPCDEM ©), along with the plant operational data and photographic evidence was used. The KMPCDEM © software was used to compare the efficiency of the radial and curved pulp lifter design based on the material flow. The slurry flow was analysed in terms of ease of discharge and lifter wear rate. It was found that changing the radial pulp lifters to curved pulp lifters could improve the discharge efficiency,and decrease wear rate of the pulp lifter in particular at the base of the pulp lifters. The simulation results showed that if curved lifters with a radius of 0.3×mill radius are used, the amount of slurry flow back and carry over will be negligible. It is predicted that by installing the curved pulp lifters, the capacity of the mill increases by 10% and the pulp lifter service life improves by at least 20%. Keywords: SAG mill, radial pulp lifter, DEM,curved pulp lifter, Sarcheshmeh

Jan 1, 2020

By K. Cole, D. Mesa, P. R. Brito-Parada, M. R. van Heerden

Turbulence in stirred tank flotation vessels determines the bulk transport of particles and has an important role in particle-bubble collisions necessary for attachment. Positron emission particle tracking (PEPT) produces quantitative measurements of the kinematic motion of different particle classes inside flotation vessels. This paper presents flow profiles of the vessel, residence time distributions, flow profiles through the mesh and macroturbulent kinetic energy distributions as experienced by valuable (hydrophobic) mineral particles in a laboratory scale flotation vessel. The results highlight the impact of the rotor-stator design and a retrofit mesh modification at the upper part of the vessel which improve recovery by decreasing turbulence in the quiescent zone and at the interface. Keywords: Flotation machines, turbulence, PEPT

Jan 1, 2020

By J. E. Dickinson, R. Campbell, M. J. Cole, K. P. Galvin, M. C. Cole

It is common practice in the beneficiation of iron ore ultrafines to apply reverse flotation, and hence float the gangue minerals, leaving behind a high-grade iron ore product. Given the difficulty in floating particles finer than 20 μm, especially those finer than 10 μm, it is necessary to subject the feed to desliming ahead of flotation. There are two problems with this approach. Firstly, given the known inefficiencies of cyclones, desliming needs to be aggressively applied to particle sizes upwards of 45 μm to ensure all the offending ultrafine gangue is removed. Secondly, desliming also removes a significant portion of the valuable iron ore. Thus, a significant portion of the ultrafines, including the iron ore, is lost to the tailings. This paper examines the direct flotation of iron ore fines using the Reflux Flotation Cell. The RFC operates using a concentrated bubbly zone, in the absence of a conventional froth zone, for delivering strong counter current washing through an inverted fluidised bed of bubbles ahead of the product discharge. A system of parallel inclined channels in the lower part of the cell provides enhanced segregation of the bubbles from the tailings, preventing the loss of air bubbles into the tailings. Experiments conducted under semi-batch conditions produced encouraging results with concentrates of ~65% Fe grade in a single flotation stage. Comparatively, up to five successive batch flotation stages were required when using a laboratory scale batch mechanical flotation cell to achieve similar Fe grades. For a given grade, lower Fe recoveries were typically achieved by the RFC. This was attributed to the effective prevention of ultrafine particles succumbing to hydraulic entrainment under a positive bias flux, as well as an increasing feed volume involved in the semi-batch operation using inverted fluidisation water for cleaning. Overall, similar grade-recovery curves were produced, however, significantly higher grades were obtained by the RFC and with fewer separation stages. Keywords: Desliming flotation, iron ore, modified release analysis, reflux flotation cell

Jan 1, 2020

By X. Cao, R. Wang, H. Han, Y. Hu, Z. Wei, W. Sun, H. Tang, J. Wang

Flotation reagents play a vital role in mineral flotation. In the past century, most studies focused on the modification of anion or cation collectors and the synergistic effect of mixed collectors, which have contributed to research. In this study, we reviewed brand-new reagents, metal-organic/inorganic complexes, with the functional group of metal ion, which are totally different from the traditional anion or cation collector and show absolute specificity for some minerals. For example, Pb-BHA complex collector was developed to realise the atmospheric flotation of scheelite due to the absolute specificity for calcium minerals. Calcium oleate collector was found to be the effective species in the flotation of calcium minerals, and has been employed to selectively separate fluorite from calcite. Selective adsorption of Al-Na2SiO3 polymer on calcite surfaces proved the self-assembly of H2SiO3 polymer using Al3+ had stronger selectively adsorption ability than Na2SiO3. Also that Fe3+-starch complex was a better depressant for iron minerals than caustic starch and contributed to a more efficient reverse flotation. In conclusion, metal ions are good templates for molecular assembly and can be used to control the structure of metal complexes and exhibit special properties, which provide possibilities for the design of new flotation agents. Keywords: Metal-organic/inorganic; flotation; collector; depressant

Jan 1, 2020

By M. E. Rossi, C. Badenhorst

Conditional bias occurs when the expected value of the true grade conditional on the estimated grade is not equal to the estimated grade. It is a quantification of how the average of estimates differs from the true value. Conditional bias is present in all resource estimates due to the unavoidable properties of kriging, and indeed all linear estimation (regression) methods. In a mining resource estimation context, the degree of conditional bias is dependent on the spacing of the conditioning data; more data leads to the expectation of less conditional bias This paper discusses conditional bias from two perspectives: 1. 2. Estimates seeking to eliminate, as much as possible, conditional bias; and Estimates that accept a larger amount of conditional bias in resource models, in order to achieve a specific outcome. The problem of conditional bias in estimation can only be fully appreciated if the objective or purpose of a resource model is understood. There are instances where conditional bias should be avoided because it is not reasonable to make decisions using estimates that are known to be wrong in expected value; for example, when making a final decision in grade control with respect to the destination of the material. In most instances in the mining resource estimation context, however, the estimated values will not be “final estimates”, but “interim estimates”, meaning that in the future more data will be obtained, and the ”interim estimate” will then be updated. In this context, the most useful resource model is an interim model that is able to predict the final model estimates, despite not having all the data available that will be collected in the future. This is the objective, for example, of a recoverable “long-term” model developed prior to operating a mine, for example at the feasibility stage. It attempts to predict the ore that a processing plant will receive, well into the future. In such cases, it can be unequivocally demonstrated that conditional bias is not only unavoidable, it is a necessity.

Sep 1, 2021

By B. M. V. Baloyi, D. Mchabe, D. Chetty

South Africa accounts for the largest global manganese ore production share, primarily from its Kalahari Manganese Field deposit. Despite this, domestic smelting is beset by energy constraints, resulting in most ore being exported for smelting. Mintek is investigating energyefficient pathways for ore reduction, specifically undertaking prereduction tests to reduce energy consumption for manganese smelting. The present contribution considers the effect of methane, anthracite, and holding time, at a fixed temperature, on the extent of reduction of Mamatwan-type ore, primarily from a mineralogical perspective. In particular, the aim was to carburise the monoxide Mn phase to Mn carbide, as this is less energy intensive than producing metallic Mn from the monoxide. Based on thermogravimetric test results, the ore and products from three experiment sets were analysed using X-ray diffraction and scanning electron microscopy to assess the phases and textures present. The findings show that carburisation of manganese monoxide and its equivalents is most advanced under conditions of 8% anthracite, 10% methane, and 60 minutes of holding time, thus confirming with the thermogravimetric data. These results show promise for reducing energy consumption for Mn smelting in ferromanganese production.

Mar 4, 2025

By O. O. Ogunmodimu, D. K. Weatherley, M. S. Powell

Discrete element methods (DEM) have provided the ability to resolve the complex phenomena experienced by ore within comminution devices such as tumbling mills. The new developments in DEM techniques and the corresponding increase in computational power have made it more feasible to study the movement of individual ore particles as they traverse a tumbling mill. A range of information can be extracted on individual particle impact histories, such as kinetic energy at the time of collision, velocities of individual particles, per-particle stress tensors, impulse, and collision energy. DEM codes have been modified to provide required information that can be analysed in associated software to investigate the impact of varied operational parameters on the resultant granular dynamics. This paper explores the reaction of individual ore particles to the breakage environment within a tumbling mill using DEM techniques and studies the stressing conditions that would lead to either impact or abrasive breakage. Keywords: DEM, comminution, tumbling mill, breakage

Jan 1, 2020

By M. Rudolph, A. van Hall, G. Öktem, S. G. Özkan, A. Hassanzadeh

Despite the four-decade study on ultrasound’s (US) impact on mineral’s floatabilities, there is still not clarity regarding its role in mineral surface wettability. For this purpose, the present investigation endeavours ……the hydrophobicity and floatability characteristics of a chalcopyrite-pyrite-quartz (Cp- Py-Qtz) system in the presence and absence of an ultrasonic pre-treatment. The ultrasonic process was carried out by a sonopuls at a constant frequency of 20 kHz with an adjustable power level from 30 to 200 W. Initially, the impact of sonication time (15, 30, 45, 60 and 90 s as well as 10, 20 and 30 min) and power magnitude (30, 60, 90, 120 and 180 W) were evaluated while the dissolved oxygen, temperature, conductivity and pH were also monitored. Collector-less micro-flotation tests were carried out on the non-pre-treated and the US pre-treated samples at 60 W and 15 s. The samples’ hydrophobicities were determined by the drop shape analysis approach. The dissolved-oxygen level was varied using a mini bench pressurised water reactor to study the effect of O2(aq) concentration on the chalcopyrite and pyrite wettability characteristics. The results showed that the minerals’ hydrophilicities were relatively more sensitive to the sonication’s time than its power that resulted in reducing all three minerals’ hydrophobicities. In addition, it was found that the dissolved-oxygen content and creation of sub-micron sized bubbles led to an improvement on chalcopyrite and pyrite hydrophobicities. Finally, we proved that the Cp’s floatability increased and Qtz’s recovery reduced while being subjected to the ultrasonic irradiation (15 s at 60 W); however, Py’s recovery remained constant. This conclusion confirmed the possibility of a selective separation in ultrasound-assisted copper ore flotation. We recommend that further advanced investigations be conducted, using an X-ray photoelectron spectroscopy (XPS) and atomic force microscopy (AFM) to thoroughly understand the surface modifications induced by the ultrasonication

Jan 1, 2020

By M. Stachowiak, M. Nowicki, K. Walas, A. Skoczylas, P. Stefaniak

In this paper the authors present a walking robot dedicated to the inspection of a belt conveyor’s transportation system in an underground mine. The developed robot can independently diagnose all of the key rotating elements of the belt conveyor (pulleys, bearings, drums, gears, engines) based on vibration and acoustic signals as well as typical RGB (red, green, blue) and thermal images. The article presents key diagnostic procedures applied to real industrial data obtained during robot testing.The work was developed as part of the ‘THING–subTerranean Haptic INvestiGator’ project. The authors presented the main assumptions of the research methodology developed for the needs of technical diagnostics of conveyor elements typical of inspection work performed by mining staff. Based on industrial data, the use of algorithms for elements such as rollers, drive stations, and the conveyor belt were shown. The main assumption of the research is to develop a robot adapted to the conditions of a deep mine, capable of carrying out an autonomous inspection of mining infrastructure. The conducted tests and the developed fault detection algorithms were aimed at preparing the robot for underground tests in the operating conditions of the mine.

Sep 1, 2021

By K. C. Corin, M. S. Manono

The mining and mineral concentration of Cu-Ni-PGM sulfide ores in South Africa occurs in semiarid regions. The scarcity and availability of water resources in these regions has become one of the biggest challenges facing mineral concentrators.S Owing to this challenge, concentrators have to find ways through which they can manage and control their water usage. Recycling and re-use of process water in mineral concentrators has therefore become common practice. This practice is beneficial in that it reduces reliance on municipal water, harnesses compliance to stringent environmental regulations on water usage, and it also offers a better response to the Sustainable Development Goals. This could however be somewhat concerning to the operator as recirculated water often has higher concentrations of ions compared to fresh or potable water. The unintended change in the process water quality may potentially affect critical aspects of flotation such as the stability of the froth. This has led to the need to find those process water ions that have the greatest impact on froth stability. Thus, this paper investigates the effect of various ions common in process water of a typical Cu-Ni-PGM concentrator on froth stability. Bubble size, water recoveries, foam height, and dynamic froth stability in both twophase and three-phase flotation systems are used as proxies for froth stability. Findings of this study are deemed crucial towards the development of a process water management protocol in sulfidic Cu-Ni- PGM ore concentrators. Keywords: Bubble diameter, electrolytes, flotation, froth stability, water quality

Jan 1, 2020

By J. Petersen, T. Moyo, D. Maharaj

This study investigates the effect of ammonium thiosulfate on the leaching of gold from the contacts of unpopulated cut electronic printed circuit boards (PCBs). Studies in the literature show ammonium thiosulfate to be a viable option in comparison to the commonly used but hazardous cyanide leaching route. However, the simultaneous dissolution of Cu interferes with the Au leaching reaction, often limiting the extraction to moderate levels. To optimize the system, an ammoniumthiosulfate system containing 0.5M ammonia, 0.5M ammonium thiosulfate, and Cu(II) sulfate at varying concentrations was studied as an alternative hydrometallurgical option for the recovery of gold from PCB. Tests were done with unpopulated standardized model PCBs (cut to 15 x 30mm pieces). Gold extraction of 47.3% was achieved after a 24-hour leaching experiment with a 0.02M copper background concentration. Further experiments focused using different initial background Cu concentrations (0M – 0.1M) to determine its effect on both gold and copper extractions. An initial Cu concentration of 0.045M resulted in the highest gold extraction, with 93.7% extraction after 6 hours of leaching. Further analysis of the data shows that the kinetic rate of gold leaching is correlated in a first-order dependence with the actual copper concentration in the solution, which increases with the concomitant dissolution of copper and gold from the PCBs. A lower threshold of 0.017M and an upper threshold of 0.15M copper in solution were determined, in between which gold extraction shows this first-order rate increase. The upper threshold is potentially linked to the onset of copper precipitation, which blinds the metal surfaces, slowing further leaching. The findings of this study suggest that in an industrial recovery operation of gold from waste PCB using the thiosulfate method, copper concentrations need to be carefully monitored, and the process should preferably operate with uncut or coarsely cut untreated boards to minimize exposure to copper circuits.

Dec 11, 2024

By C. Gatschlhofer, A. Krammer, H. Raupenstrauch, T. Hochsteiner, A. Holzer, L. Wiszniewski

The rising demand for lithium-ion batteries, coupled with international climate targets, emphasises the need for circular economy and zero-waste recycling approaches. However, so far, no industrial-scaled solutions are available that fulfil global recovery targets. Therefore, a novel pyrometallurgical reactor was developed at the Chair of Thermal Processing Technology at Montanuniversitaet Leoben to avoid lithium slagging and recover volatile elements via the gas phase. Nevertheless, previous results show that several impact factors, including interfering elements or pre-treatment steps, strongly limit the recycling efficiency of black mass. To simulate the black mass and provide feedback on possible influences of pre-treatment steps, cathode material from production with additives of interfering elements was analysed regarding melting behaviour. The results show that formation of a further processable alloy was increased by almost 50% when pretreatment steps were excluded, indicating pre-treatment steps as a major influence to be considered within this specific recycling route.

Jan 1, 2023

By N. Snyders, I. Mwandawande, S. M. Bradshaw, M. Karimi

"Flotation columns are normally operated at optimal superficial gas velocities to maintain bubbly flow conditions. However, with increasing superficial gas velocity, loss of bubbly flow may occur with adverse effects on column performance. It is therefore important to identify the maximum superficial gas velocity above which loss of bubbly flow occurs. The maximum superficial gas velocity is usually obtained from a gas holdup versus superficial gas velocity plot in which the linear portion of the graph represents bubbly flow while deviation from the linear relationship indicates a change from the bubbly flow to the churn-turbulent regime. However, this method is difficult to use when the transition from bubbly flow to churn-turbulent flow is gradual, as happens in the presence of frothers. We present two alternative methods in which the flow regime in the column is distinguished by means of radial gas holdup profiles and gas holdup versus time graphs obtained from CFD simulations. Bubbly flow was characterized by saddle-shaped profiles with three distinct peaks, or saddle-shaped profiles with two near-wall peaks and a central minimum, or flat profiles with intermediate features between saddle and parabolic gas holdup profiles. The transition regime was gradual and characterized by flat to parabolic gas holdup profiles that become steeper with increasing superficial gas velocity. The churn-turbulent flow was distinguished by steep parabolic radial gas holdup profiles. Gas holdup versus time graphs were also used to define flow regimes with a constant gas holdup indicating bubbly flow, while wide gas holdup variations indicate churn-turbulent flow. IntroductionTwo types of flow can be distinguished in gasliquid flows in bubble columns: bubbly flow and churn-turbulent flow. The bubbly flow regime is characterized by uniform flow of bubbles of uniform size. Churn-turbulent flow, on the other hand, is characterized by a wide variation in bubble sizes with large bubbles rising rapidly and causing liquid circulation.Flotation columns are normally operated in the bubbly flow regime, which is the optimal condition for column flotation (Finch and Dobby, 1990; Yianatos and Finch, 1990; Ityokumbul, 1992). However, it has been generally observed that flotation column performance deteriorates when the superficial gas velocity is increased beyond a certain limit where the flow regime changes from bubbly flow to churn-turbulent flow (Xu, Finch, and Uribe-Salas, 1991). The identification of this critical or maximum superficial gas velocity is therefore important for optimal operation of flotation columns."

Feb 1, 2019

By D. F. Malan, D. Ile

This study explores the use of backfill in hard rock bord-and-pillar mines to increase the pillar strength and extraction ratio at depth. The use of backfill will also minimize the requirement for tailings storage on surface and the risk of environmental damage. A literature survey indicated that backfill is extensively used in coal mines, but rarely in hard rock bord-and-pillar mines. To simulate the effect of backfill confinement on pillar strength, an extension of the limit equilibrium model is proposed. Numerical modelling of an actual platinum mine layout is used to illustrate the beneficial effect of backfill on pillar stability at greater depths. The magnitude of confinement exerted by the backfill on the pillar sidewalls is unknown, however, and this needs to be quantified using experimental backfill mining sections equipped with suitable instrumentation.

May 3, 2023

By N. Snyders, I. Mwandawande, S. M. Bradshaw, M. Karimi

"Computational fluid dynamics (CFD) was applied to predict the average gas holdup and the axial gas holdup variation in a 13.5 m high cylindrical column 0.91 m diameter. The column was operating in batch mode. A Eulerian-Eulerian multiphase approach with appropriate interphase momentum exchange terms was applied to simulate the gas-liquid flow inside the column. Turbulence in the continuous phase was modelled using the k-_ realizable turbulence model. The predicted average gas holdup values were in good agreement with experimental data. The axial gas holdup prediction was generally good for the middle and top parts of the column, but was over-predicted for the bottom part of the column. Bubble velocity profiles were observed in which the axial velocity of the air bubbles decreased with height in the column. This may be related to the upward increase in gas holdup in the column. Simulations were also conducted to compare the gas holdup predicted with the universal, the Schiller-Naumann, and the Morsi-Alexander drag models. The gas holdup predictions for the three drag models were not significantly different. IntroductionColumn flotation is an important concentration technology that is used in the mineral processing and coal beneficiation industries. The growing interest in the use of column flotation in mineral processing has been attributed to the simpler flotation circuits and improved metallurgical performance compared to conventional flotation cells (Finch and Dobby, 1990). Flotation columns have also found other applications outside mineral processing, such as de-inking of recycled paper (Finch and Hardie, 1999).In column flotation, a rising swarm of air bubbles generated by means of air spargers is employed to collect the valuable mineral particles and separate them from the gangue minerals in a countercurrent process. Wash water, which is continuously fed at the top of the column, is used to eliminate entrained particles and stabilize the froth. The column volume can be divided into two sections – the collection zone in which the bubbles collect the floatable mineral particles, and the cleaning zone (or froth zone) where product upgrading is enhanced through the removal of unwanted particles entrained in the water mixed with bubbles from the collection zone."

Jan 1, 2019